From 8b9e919705f879e73071fd64b410216a1be783a9 Mon Sep 17 00:00:00 2001
From: Anurag Kulshrestha <anurag.kulshrestha@ansys.com>
Date: Sun, 28 Aug 2022 23:34:18 +0530
Subject: [PATCH] New tests added

---
 tests/test_solvermode/test_controls.py        | 133 ++++++++++++++++++
 tests/test_solvermode/test_general.py         |  49 +++++++
 tests/test_solvermode/test_initialization.py  |  71 ++++++++++
 tests/test_solvermode/test_materials.py       |  37 +++++
 tests/test_solvermode/test_methods.py         |  60 ++++++++
 tests/test_solvermode/test_models.py          | 127 +++++++++++++++++
 .../test_solvermode/test_named_expressions.py |  44 ++++++
 tests/util/fixture_fluent.py                  |   2 +-
 8 files changed, 522 insertions(+), 1 deletion(-)
 create mode 100644 tests/test_solvermode/test_controls.py
 create mode 100644 tests/test_solvermode/test_initialization.py
 create mode 100644 tests/test_solvermode/test_methods.py
 create mode 100644 tests/test_solvermode/test_named_expressions.py

diff --git a/tests/test_solvermode/test_controls.py b/tests/test_solvermode/test_controls.py
new file mode 100644
index 00000000000..892caadb37b
--- /dev/null
+++ b/tests/test_solvermode/test_controls.py
@@ -0,0 +1,133 @@
+import pytest
+
+
+@pytest.mark.integration
+@pytest.mark.quick
+@pytest.mark.setup
+def test_controls(load_mixing_elbow_mesh):
+    solver = load_mixing_elbow_mesh
+    solver.setup.models.multiphase.models = "vof"
+    assert solver.setup.models.multiphase.models() == "vof"
+    solver.setup.general.gravity = {"enable": True, "components": [0.0, 0.0, -9.81]}
+    assert solver.setup.general.gravity.components() == [0, 0, -9.81]
+    solver.setup.general.solver.time = "steady"
+    assert solver.setup.general.solver.time() == "steady"
+    solver.solution.controls.advanced.multi_grid.algebraic_mg_controls.coupled_parameters.coarsening_parameters.laplace_coarsening = (
+        True
+    )
+    assert solver.solution.controls.advanced.multi_grid.algebraic_mg_controls.coupled_parameters.coarsening_parameters() == {
+        "max_coarse_levels": 40,
+        "coarsen_by_interval": 8,
+        "conservative_coarsening": False,
+        "aggressive_coarsening": False,
+        "laplace_coarsening": True,
+    }
+    solver.solution.controls.advanced.multi_grid.algebraic_mg_controls.coupled_parameters.coarsening_parameters.max_coarse_levels = (
+        45
+    )
+    assert (
+        solver.solution.controls.advanced.multi_grid.algebraic_mg_controls.coupled_parameters.coarsening_parameters.max_coarse_levels()
+        == 45
+    )
+    solver.solution.controls.advanced.multi_grid.algebraic_mg_controls.coupled_parameters.coarsening_parameters = {
+        "max_coarse_levels": 48,
+        "coarsen_by_interval": 9,
+        "conservative_coarsening": True,
+        "aggressive_coarsening": True,
+        "laplace_coarsening": True,
+    }
+    assert (
+        solver.solution.controls.advanced.multi_grid.algebraic_mg_controls.coupled_parameters.coarsening_parameters.max_coarse_levels()
+        == 48
+    )
+    assert (
+        solver.solution.controls.advanced.multi_grid.algebraic_mg_controls.coupled_parameters.coarsening_parameters.coarsen_by_interval()
+        == 9
+    )
+    assert (
+        solver.solution.controls.advanced.multi_grid.algebraic_mg_controls.coupled_parameters.coarsening_parameters.conservative_coarsening()
+        == True
+    )
+    assert (
+        solver.solution.controls.advanced.multi_grid.algebraic_mg_controls.coupled_parameters.coarsening_parameters.aggressive_coarsening()
+        == True
+    )
+    solver.solution.controls.advanced.multi_grid.algebraic_mg_controls.scalar_parameters.fixed_cycle_parameters.max_cycle = (
+        300
+    )
+    assert (
+        solver.solution.controls.advanced.multi_grid.algebraic_mg_controls.scalar_parameters.fixed_cycle_parameters.max_cycle()
+        == 300
+    )
+    solver.solution.controls.advanced.multi_grid.algebraic_mg_controls.scalar_parameters.fixed_cycle_parameters = {
+        "pre_sweeps": 1,
+        "post_sweeps": 2,
+        "max_cycle": 350,
+    }
+    assert solver.solution.controls.advanced.multi_grid.algebraic_mg_controls.scalar_parameters.fixed_cycle_parameters() == {
+        "pre_sweeps": 1,
+        "post_sweeps": 2,
+        "max_cycle": 350,
+    }
+    solver.solution.methods.p_v_coupling.flow_scheme = "Coupled"
+    assert solver.solution.methods.p_v_coupling.flow_scheme() == "Coupled"
+    solver.solution.methods.p_v_coupling.coupled_form = True
+    assert solver.solution.methods.p_v_coupling.coupled_form() == True
+    solver.solution.controls.advanced.multi_grid.algebraic_mg_controls.scalar_parameters.smoother_type = (
+        "gauss-seidel"
+    )
+    assert (
+        solver.solution.controls.advanced.multi_grid.algebraic_mg_controls.scalar_parameters.smoother_type()
+        == "gauss-seidel"
+    )
+    assert solver.solution.controls.pseudo_time_explicit_relaxation_factor() == {
+        "global_dt_pseudo_relax": {
+            "turb-viscosity": 1.0,
+            "body-force": 1.0,
+            "k": 0.75,
+            "omega": 0.75,
+            "mp": 0.5,
+            "density": 1.0,
+        }
+    }
+    solver.solution.controls.pseudo_time_explicit_relaxation_factor = {
+        "global_dt_pseudo_relax": {"turb-viscosity": 0.9, "body-force": 0.8}
+    }
+    assert solver.solution.controls.pseudo_time_explicit_relaxation_factor() == {
+        "global_dt_pseudo_relax": {
+            "turb-viscosity": 0.9,
+            "body-force": 0.8,
+            "k": 0.75,
+            "omega": 0.75,
+            "mp": 0.5,
+            "density": 1.0,
+        }
+    }
+    solver.solution.methods.p_v_coupling.flow_scheme = "SIMPLE"
+    solver.solution.controls.under_relaxation = {"pressure": 0.9}
+    solver.solution.controls.under_relaxation = {"density": 0.9}
+    assert solver.solution.controls.under_relaxation() == {
+        "mom": 0.7,
+        "turb-viscosity": 1,
+        "density": 0.9,
+        "omega": 0.8,
+        "mp": 0.5,
+        "body-force": 1.0,
+        "pressure": 0.9,
+        "k": 0.8,
+    }
+    solver.solution.controls.under_relaxation = {
+        "body-force": 0.7,
+        "density": 0.75,
+        "mom": 0.8,
+    }
+    assert solver.solution.controls.under_relaxation() == {
+        "mom": 0.8,
+        "turb-viscosity": 1,
+        "density": 0.75,
+        "omega": 0.8,
+        "mp": 0.5,
+        "body-force": 0.7,
+        "pressure": 0.9,
+        "k": 0.8,
+    }
diff --git a/tests/test_solvermode/test_general.py b/tests/test_solvermode/test_general.py
index 2141dec0701..c47f4b8cb23 100644
--- a/tests/test_solvermode/test_general.py
+++ b/tests/test_solvermode/test_general.py
@@ -17,3 +17,52 @@ def test_solver_import_mixingelbow(load_mixing_elbow_mesh):
         solver_session.scheme_eval.scheme_eval('(units/quantity-info "length")')[-1]
         == "in"
     )
+    solver_session.setup.general.solver.time.get_attr("allowed-values")
+    solver_session.setup.general.solver.time = "unsteady-2nd-order"
+    solver_session.setup.general.solver.time = "unsteady-1st-order"
+    solver_session.setup.general.solver.time = "unsteady-2nd-order-bounded"
+    solver_session.setup.general.solver.time = "steady"
+
+    # solver.setup.general.gravity = {"gravity": True, "y_component": -9.81}
+    # solver.mesh.scale(x_scale=0.001, y_scale=0.001, z_scale=0.001)
+    solver_session.setup.general.solver.type.get_attr("allowed-values")
+    solver_session.setup.general.solver.type = "density-based-implicit"
+    assert solver_session.setup.general.solver.type() == "density-based-implicit"
+    solver_session.setup.general.solver.type = "density-based-explicit"
+    assert solver_session.setup.general.solver.type() == "density-based-explicit"
+    solver_session.setup.general.solver.type = "pressure-based"
+    assert solver_session.setup.general.solver.type() == "pressure-based"
+
+    solver_session.file.auto_save.data_frequency = 10
+    assert solver_session.file.auto_save.data_frequency() == 10
+    solver_session.file.auto_save.case_frequency = "each-time"
+    assert solver_session.file.auto_save.case_frequency() == "each-time"
+    solver_session.file.auto_save.root_name = "file_auto_save"
+    assert solver_session.file.auto_save.root_name() == "file_auto_save"
+    solver_session.setup.reference_values.compute(from_zone_name="outlet")
+
+
+@pytest.mark.quick
+@pytest.mark.setup
+def test_disk_2d_setup(load_disk_mesh):
+    session = load_disk_mesh
+    assert session._root.get_attr("active?")
+    assert session.check_health() == "SERVING"
+    ###
+    assert not session.setup.models.energy.enabled()
+    assert session.scheme_eval.scheme_eval("(case-valid?)")
+    session.tui.mesh.check()
+    assert session.setup.general.solver.two_dim_space.get_attr("allowed-values") == [
+        "swirl",
+        "axisymmetric",
+        "planar",
+    ]
+    assert session.setup.general.solver.two_dim_space() == "planar"
+    session.setup.general.solver.two_dim_space = "axisymmetric"
+    assert session.setup.general.solver.two_dim_space() == "axisymmetric"
+    session.setup.general.solver.two_dim_space = "swirl"
+    assert session.setup.general.solver.two_dim_space() == "swirl"
+    session.setup.general.solver.two_dim_space = "planar"
+    assert session.setup.general.solver.two_dim_space() == "planar"
+    # Bug 682773
+    # session.setup.general.gravity = {"gravity": True, "x_component": -9.81}
diff --git a/tests/test_solvermode/test_initialization.py b/tests/test_solvermode/test_initialization.py
new file mode 100644
index 00000000000..6b0fc51c7a2
--- /dev/null
+++ b/tests/test_solvermode/test_initialization.py
@@ -0,0 +1,71 @@
+import pytest
+from util.fixture_fluent import download_input_file
+
+
+@pytest.mark.quick
+@pytest.mark.setup
+def test_initialize(launch_fluent_solver_3ddp_t2):
+    solver = launch_fluent_solver_3ddp_t2
+    input_type, input_name = download_input_file("pyfluent/wigley_hull", "wigley.msh")
+    solver.file.read(file_type=input_type, file_name=input_name)
+    solver.parallel.partition.set.use_laplace_smoothing = True
+    solver.parallel.partition.method(partition_method="metis", count=2)
+    solver.setup.materials.database.copy_by_name(type="fluid", name="air")
+    solver.setup.materials.database.copy_by_name(type="fluid", name="water-liquid")
+    solver.setup.models.multiphase.models = "vof"
+    solver.setup.general.gravity = {"enable": True, "components": [0.0, 0.0, -9.81]}
+    solver.setup.general.solver.time = "steady"
+    solver.execute_tui("/define/models/multiphase/vof-sub-models yes no ")
+    solver.execute_tui(
+        "/define/phases/set-domain-properties/change-phases-names? water air "
+    )
+    solver.setup.boundary_conditions.pressure_inlet["inflow"].phase["mixture"] = {
+        "open_channel": True,
+        "direction_spec": "Direction Vector",
+        "vmag": 1.452,
+        "ht_bottom": -0.941875,
+        "turb_intensity": 0.01,
+        "turb_viscosity_ratio": 1,
+    }
+    solver.setup.boundary_conditions.pressure_outlet["outflow"].phase["mixture"] = {
+        "open_channel": True,
+        "ht_bottom": -0.941875,
+        "den_spec": "From Free Surface Level",
+        "direction_spec": "Normal to Boundary",
+        "turb_intensity": 0.01,
+        "turb_viscosity_ratio": 1,
+        "p_backflow_spec_gen": "Static Pressure",
+    }
+
+    solver.solution.methods.p_v_coupling.flow_scheme = "Coupled"
+    solver.solution.methods.p_v_coupling.coupled_form = True
+    solver.solution.controls.advanced.multi_grid.algebraic_mg_controls.coupled_parameters.coarsening_parameters.laplace_coarsening = (
+        True
+    )
+    solver.solution.initialization.open_channel_auto_init = {
+        "boundary_thread": 3,
+        "flat_init": True,
+    }
+    assert solver.solution.initialization.open_channel_auto_init() == {
+        "boundary_thread": 3,
+        "flat_init": True,
+    }
+    # solver.solution.initialization.hybrid_initialize()
+    # solver.exit()
+
+
+@pytest.mark.quick
+@pytest.mark.setup
+def test_fmg_initialize(launch_fluent_solver_3ddp_t2):
+    solver = launch_fluent_solver_3ddp_t2
+    input_type, input_name = download_input_file(
+        "pyfluent/vki_turbine", "vki_turbine.cas.gz"
+    )
+    solver.file.read(file_type=input_type, file_name=input_name)
+    solver.execute_tui("(rpsetvar 'fmg-init/reorder? #f)  ")
+    solver.mesh.check()
+    solver.solution.initialization.standard_initialize()
+    solver.solution.initialization.fmg_initialize = True
+    # assert solver.solution.initialization.fmg_initialize() == True
+    solver.execute_tui("(benchmark '(iterate 2))  ")
+    # solver.solution.initialization.hybrid_initialize()
diff --git a/tests/test_solvermode/test_materials.py b/tests/test_solvermode/test_materials.py
index 94fe0ec990f..55d83df9665 100644
--- a/tests/test_solvermode/test_materials.py
+++ b/tests/test_solvermode/test_materials.py
@@ -21,3 +21,40 @@ def test_solver_material(load_mixing_elbow_mesh):
         "water-liquid"
         in solver_session.setup.cell_zone_conditions.fluid["elbow-fluid"].material()
     )
+
+    solver_session.setup.materials.database.copy_by_name(type="fluid", name="air")
+    solver_session.setup.cell_zone_conditions.fluid["elbow-fluid"].material = "air"
+    assert (
+        "air"
+        in solver_session.setup.cell_zone_conditions.fluid["elbow-fluid"].material()
+    )
+    assert (
+        "air"
+        in solver_session.setup.cell_zone_conditions.fluid["elbow-fluid"].material()
+    )
+    assert solver_session.setup.materials.child_names == [
+        "database",
+        "fluid",
+        "solid",
+        "mixture",
+        "inert_particle",
+        "droplet_particle",
+        "combusting_particle",
+        "particle_mixture",
+    ]
+    assert solver_session.setup.materials.database.get_active_command_names() == [
+        "copy_by_formula",
+        "copy_by_name",
+        "list_materials",
+        "list_properties",
+    ]
+
+    solver_session.setup.materials.database.copy_by_formula(
+        type="fluid", formula="c2h6"
+    )
+
+    solver_session.setup.cell_zone_conditions.fluid["elbow-fluid"].material = "ethane"
+    assert (
+        "ethane"
+        in solver_session.setup.cell_zone_conditions.fluid["elbow-fluid"].material()
+    )
diff --git a/tests/test_solvermode/test_methods.py b/tests/test_solvermode/test_methods.py
new file mode 100644
index 00000000000..f89bee4174e
--- /dev/null
+++ b/tests/test_solvermode/test_methods.py
@@ -0,0 +1,60 @@
+import pytest
+
+
+@pytest.mark.quick
+@pytest.mark.setup
+def test_methods(load_mixing_elbow_mesh):
+    solver = load_mixing_elbow_mesh
+    solver.setup.models.multiphase.models = "vof"
+    solver.setup.general.gravity = {"enable": True, "components": [0.0, 0.0, -9.81]}
+    solver.setup.general.solver.time = "steady"
+    solver.solution.methods.p_v_coupling.flow_scheme = "Coupled"
+    solver.solution.methods.p_v_coupling.coupled_form = False
+    assert solver.solution.methods.p_v_coupling() == {
+        "flow_scheme": "Coupled",
+        "coupled_form": False,
+    }
+    solver.solution.methods.discretization_scheme = {"pressure": "presto!"}
+    assert solver.solution.methods.discretization_scheme() == {
+        "mom": "second-order-upwind",
+        "omega": "second-order-upwind",
+        "mp": "compressive",
+        "pressure": "presto!",
+        "k": "second-order-upwind",
+    }
+    solver.solution.methods.gradient_scheme = "least-square-cell-based"
+    assert solver.solution.methods.gradient_scheme() == "least-square-cell-based"
+    solver.solution.methods.warped_face_gradient_correction.enable(
+        enable=True, gradient_correction_mode="fast-mode"
+    )
+    solver.solution.methods.warped_face_gradient_correction.enable(
+        enable=False, gradient_correction_mode="fast-mode"
+    )
+    solver.solution.methods.expert.numerics_pbns.velocity_formulation = "relative"
+    assert (
+        solver.solution.methods.expert.numerics_pbns.velocity_formulation()
+        == "relative"
+    )
+    solver.solution.methods.expert.numerics_pbns = {
+        "implicit_bodyforce_treatment": True,
+        "velocity_formulation": "absolute",
+        "physical_velocity_formulation": True,
+        "disable_rhie_chow_flux": True,
+        "presto_pressure_scheme": False,
+        "first_to_second_order_blending": 2.0,
+    }
+    assert solver.solution.methods.expert.numerics_pbns() == {
+        "implicit_bodyforce_treatment": True,
+        "velocity_formulation": "absolute",
+        "physical_velocity_formulation": True,
+        "disable_rhie_chow_flux": True,
+        "presto_pressure_scheme": False,
+        "first_to_second_order_blending": 2.0,
+    }
+    solver.solution.methods.expert.numerics_pbns.presto_pressure_scheme = True
+    assert solver.solution.methods.expert.numerics_pbns.presto_pressure_scheme() == True
+    solver.solution.methods.gradient_scheme = "green-gauss-node-based"
+    assert solver.solution.methods.gradient_scheme() == "green-gauss-node-based"
+    solver.solution.methods.warped_face_gradient_correction.enable(
+        enable=True, gradient_correction_mode="memory-saving-mode"
+    )
diff --git a/tests/test_solvermode/test_models.py b/tests/test_solvermode/test_models.py
index 16be8cb5b94..83ae7d062dd 100644
--- a/tests/test_solvermode/test_models.py
+++ b/tests/test_solvermode/test_models.py
@@ -9,3 +9,130 @@ def test_solver_models(load_mixing_elbow_mesh):
     assert not solver_session.setup.models.energy.enabled()
     solver_session.setup.models.energy.enabled = True
     assert solver_session.setup.models.energy.enabled()
+    solver_session.setup.models.multiphase.models = "vof"
+    assert solver_session.setup.models.multiphase.models() == "vof"
+    solver_session.setup.models.viscous.model = "laminar"
+    assert solver_session.setup.models.viscous.model() == "laminar"
+    solver_session.setup.models.viscous.model = "k-epsilon"
+    assert solver_session.setup.models.viscous.model() == "k-epsilon"
+    solver_session.setup.models.viscous.near_wall_treatment.wall_function = (
+        "enhanced-wall-treatment"
+    )
+    solver_session.setup.models.multiphase.models = "eulerian"
+    assert solver_session.setup.models.multiphase.models() == "eulerian"
+
+
+@pytest.mark.quick
+@pytest.mark.setup
+def test_disk_2d_models(load_disk_mesh):
+    solver_session = load_disk_mesh
+    solver_session.setup.general.solver.two_dim_space = "axisymmetric"
+    solver_session.setup.general.solver.two_dim_space = "swirl"
+    solver_session.setup.models.viscous.model = "k-epsilon"
+    assert solver_session.setup.models.viscous.model() == "k-epsilon"
+    solver_session.setup.models.viscous.near_wall_treatment.wall_function = (
+        "enhanced-wall-treatment"
+    )
+    assert (
+        solver_session.setup.models.viscous.near_wall_treatment.wall_function()
+        == "enhanced-wall-treatment"
+    )
+    solver_session.setup.models.viscous.near_wall_treatment.wall_function.get_attr(
+        "allowed-values"
+    )
+    solver_session.setup.models.viscous.near_wall_treatment.wall_function = (
+        "standard-wall-fn"
+    )
+    assert (
+        solver_session.setup.models.viscous.near_wall_treatment.wall_function()
+        == "standard-wall-fn"
+    )
+    solver_session.setup.models.viscous.near_wall_treatment.wall_function = (
+        "non-equilibrium-wall-fn"
+    )
+    assert (
+        solver_session.setup.models.viscous.near_wall_treatment.wall_function()
+        == "non-equilibrium-wall-fn"
+    )
+    solver_session.setup.models.viscous.near_wall_treatment.wall_function = (
+        "menter-lechner"
+    )
+    assert (
+        solver_session.setup.models.viscous.near_wall_treatment.wall_function()
+        == "menter-lechner"
+    )
+    solver_session.setup.models.viscous.near_wall_treatment.wall_function = (
+        "scalable-wall-functions"
+    )
+    assert (
+        solver_session.setup.models.viscous.near_wall_treatment.wall_function()
+        == "scalable-wall-functions"
+    )
+    solver_session.setup.models.viscous.near_wall_treatment.wall_function = (
+        "user-defined-wall-functions"
+    )
+    # As udf is not available, so it will change automatically to standard wall function
+    assert (
+        solver_session.setup.models.viscous.near_wall_treatment.wall_function()
+        == "standard-wall-fn"
+    )
+
+    solver_session.setup.models.viscous.model = "k-omega"
+    assert solver_session.setup.models.viscous.model() == "k-omega"
+    assert solver_session.setup.models.viscous.k_omega_model.get_attr(
+        "allowed-values"
+    ) == [
+        "standard",
+        "geko",
+        "bsl",
+        "sst",
+    ]
+    assert solver_session.setup.models.viscous.k_omega_model() == "standard"
+    solver_session.setup.models.viscous.k_omega_model = "geko"
+    assert solver_session.setup.models.viscous.k_omega_model() == "geko"
+    solver_session.setup.models.viscous.k_omega_model = "bsl"
+    assert solver_session.setup.models.viscous.k_omega_model() == "bsl"
+    solver_session.setup.models.viscous.k_omega_model = "sst"
+    assert solver_session.setup.models.viscous.k_omega_model() == "sst"
+    assert solver_session.setup.models.viscous.model.get_attr("allowed-values") == [
+        "inviscid",
+        "laminar",
+        "k-epsilon",
+        "k-omega",
+        "mixing-length",
+        "spalart-allmaras",
+        "k-kl-w",
+        "transition-sst",
+        "reynolds-stress",
+        "scale-adaptive-simulation",
+        "detached-eddy-simulation",
+    ]
+    solver_session.setup.models.viscous.k_omega_options.kw_low_re_correction = True
+    assert (
+        solver_session.setup.models.viscous.k_omega_options.kw_low_re_correction()
+        == True
+    )
+    solver_session.setup.models.viscous.turbulence_expert.kato_launder_model = True
+    assert (
+        solver_session.setup.models.viscous.turbulence_expert.kato_launder_model()
+        == True
+    )
+    solver_session.setup.models.viscous.turbulence_expert.production_limiter.clip_factor = (
+        9
+    )
+    assert (
+        solver_session.setup.models.viscous.turbulence_expert.production_limiter.clip_factor()
+        == 9
+    )
+    solver_session.execute_tui(
+        "/define/models/viscous/turbulence-expert/thermal-p-function? yes "
+    )
+    solver_session.setup.models.viscous.turbulence_expert.turb_non_newtonian = True
+    assert (
+        solver_session.setup.models.viscous.turbulence_expert.turb_non_newtonian()
+        == True
+    )
+    solver_session.setup.models.viscous.model = "laminar"
+    assert solver_session.setup.models.viscous.model() == "laminar"
+    solver_session.setup.models.viscous.model = "spalart-allmaras"
+    assert solver_session.setup.models.viscous.model() == "spalart-allmaras"
diff --git a/tests/test_solvermode/test_named_expressions.py b/tests/test_solvermode/test_named_expressions.py
new file mode 100644
index 00000000000..a06144f3680
--- /dev/null
+++ b/tests/test_solvermode/test_named_expressions.py
@@ -0,0 +1,44 @@
+import pytest
+
+
+@pytest.mark.quick
+@pytest.mark.setup
+def test_expression(load_mixing_elbow_mesh):
+    solver_session = load_mixing_elbow_mesh
+    solver_session.setup.models.energy.enabled = True
+    solver_session.setup.named_expressions["r"] = {}
+    solver_session.setup.named_expressions["r"] = {
+        "definition": "(Position.z**2.0 +Position.x**2.0)**0.5"
+    }
+    solver_session.setup.named_expressions["r1"] = {}
+    solver_session.setup.named_expressions["r1"] = {"definition": "1-(r/0.014[m])"}
+    solver_session.setup.named_expressions["v1"] = {}
+    solver_session.setup.named_expressions["v1"] = {"definition": "r1**(1.0/6.0)"}
+    solver_session.setup.named_expressions["vel_cold"] = {}
+    solver_session.setup.named_expressions["vel_cold"] = {
+        "definition": "1.264 * 1.43 [m s^-1] * max(0,v1)"
+    }
+    assert (
+        solver_session.setup.named_expressions["r"].definition()
+        == "(Position.z**2.0 +Position.x**2.0)**0.5"
+    )
+    assert solver_session.setup.named_expressions["r1"].definition() == "1-(r/0.014[m])"
+    assert solver_session.setup.named_expressions["v1"].definition() == "r1**(1.0/6.0)"
+    assert (
+        solver_session.setup.named_expressions["vel_cold"].definition()
+        == "1.264 * 1.43 [m s^-1] * max(0,v1)"
+    )
+    solver_session.setup.boundary_conditions.velocity_inlet["cold-inlet"] = {
+        "vmag": "vel_cold",
+        "turb_intensity": 0.0999999,
+        "turb_hydraulic_diam": 1.0,
+    }
+    solver_session.setup.boundary_conditions.velocity_inlet["hot-inlet"] = {
+        "vmag": "max(vel_cold, 1.5 [m/s])"
+    }
+    assert solver_session.setup.boundary_conditions.velocity_inlet[
+        "cold-inlet"
+    ].vmag() == {"option": "value", "value": "vel_cold"}
+    assert solver_session.setup.boundary_conditions.velocity_inlet[
+        "hot-inlet"
+    ].vmag() == {"option": "value", "value": "max(vel_cold, 1.5 [m/s])"}
diff --git a/tests/util/fixture_fluent.py b/tests/util/fixture_fluent.py
index 3f8afa22d40..71d7e4df2ce 100644
--- a/tests/util/fixture_fluent.py
+++ b/tests/util/fixture_fluent.py
@@ -33,7 +33,7 @@ def download_input_file(directory_name, full_file_name, data_file_name=None):
         globals()[dat_name] = None
         if not globals()[dat_name]:
             globals()[dat_name] = download_file(
-                filename=full_file_name,
+                filename=data_file_name,
                 directory=directory_name,
             )
         file_type = "case-data"