def METHOD_NAME(self, ED1D1Stack_conc_recirc):
    m = ED1D1Stack_conc_recirc
    solver = get_solver()
    init_arg = {
        ("flow_vol_phase", ("Liq")): 5.2e-4,
        ("conc_mol_phase_comp", ("Liq", "Na_+")): 34.188,
        ("conc_mol_phase_comp", ("Liq", "Cl_-")): 34.188,
    }  # Corresponding to C_feed = 2g/L
    m.fs.feed.properties.calculate_state(
        init_arg,
        hold_state=True,
    )
    m.fs.EDstack.voltage_applied[0].fix(10)
    m.fs.recovery_vol_H2O.fix(0.7)
    edfs.condition_base(m)
    check_dof(m)
    # Initialize and solve the model
    edfs.initialize_system(m, solver=solver)
    edfs.solve(m, solver=solver)
    assert value(m.fs.feed.properties[0].flow_vol_phase["Liq"]) == pytest.approx(
        5.2e-4, rel=1e-3
    )
    assert value(m.fs.prod.properties[0].flow_vol_phase["Liq"]) == pytest.approx(
        3.64e-4, rel=1e-3
    )
    assert value(m.fs.disp.properties[0].flow_vol_phase["Liq"]) == pytest.approx(
        1.56e-4, rel=1e-3
    )
    assert value(
        sum(
            m.fs.feed.properties[0].conc_mass_phase_comp["Liq", j]
            for j in m.fs.properties.ion_set
        )
    ) == pytest.approx(2.00, rel=1e-3)
    assert value(m.fs.product_salinity) == pytest.approx(1.1333, rel=1e-3)
    assert value(m.fs.disposal_salinity) == pytest.approx(4.0223, rel=1e-3)
    assert value(m.fs.mem_area) == pytest.approx(18.5338, rel=1e-3)
    assert value(m.fs.costing.specific_energy_consumption) == pytest.approx(
        0.1192, abs=0.001
    )
    assert value(m.fs.costing.LCOW) == pytest.approx(0.397526, rel=1e-3)
    assert value(m.fs.EDstack.inlet_concentrate.pressure[0]) == pytest.approx(
        169278.127, rel=1e-3
    )
    assert value(m.fs.EDstack.outlet_concentrate.pressure[0]) == pytest.approx(
        101325.00, rel=1e-3
    )
    assert value(m.fs.EDstack.inlet_diluate.pressure[0]) == pytest.approx(
        169278.127, rel=1e-3
    )
    assert value(m.fs.EDstack.outlet_diluate.pressure[0]) == pytest.approx(
        101325.00, rel=1e-3
    )
    assert value(m.fs.EDstack.inlet_concentrate.temperature[0]) == pytest.approx(
        298.15, rel=1e-3
    )
    assert value(m.fs.EDstack.outlet_concentrate.temperature[0]) == pytest.approx(
        298.15, rel=1e-3
    )
    assert value(m.fs.EDstack.inlet_diluate.temperature[0]) == pytest.approx(
        298.15, rel=1e-3
    )
    assert value(m.fs.EDstack.outlet_diluate.temperature[0]) == pytest.approx(
        298.15, rel=1e-3
    )