Comparison of Nested Cylinder Test

Comparison of Nested Cylinder Test#

V1 Client (old)

 1with fl_v1.flow360_unit_system:
 2    params = fl_v1.Flow360Params(
 3        volume_output=fl_v1.VolumeOutput(
 4            output_format="paraview",
 5            output_fields=["primitiveVars", "Mach", "VelocityRelative"],
 6        ),
 7        surface_output=fl_v1.SurfaceOutput(
 8            output_format="paraview",
 9            output_fields=["primitiveVars", "Cp", "Cf", "yPlus"],
10        ),
11        navier_stokes_solver=fl_v1.NavierStokesSolver(
12            absolute_tolerance=1e-9,
13            linear_solver=fl_v1.LinearSolver(max_iterations=35),
14        ),
15        turbulence_model_solver=fl_v1.NoneSolver(),
16        freestream=fl_v1.FreestreamFromMachReynolds(
17            Mach=0.05, Reynolds=10, temperature=288.15, alpha=-90, beta=0
18        ),
19        time_stepping=fl_v1.SteadyTimeStepping(
20            CFL=fl_v1.AdaptiveCFL(), max_pseudo_steps=2000
21        ),
22        boundaries={
23            "blk-1/InletOutlet": fl_v1.FreestreamBoundary(),
24            "blk-1/OuterWall": fl_v1.NoSlipWall(velocity=[0, 0, 0]),
25            "blk-1/Cylinder": fl_v1.NoSlipWall(),
26        },
27        volume_zones={
28            "blk-1": fl_v1.FluidDynamicsVolumeZone(
29                reference_frame=ReferenceFrameOmegaRadians(
30                    center=(0, 0, 0), axis=(0, 0, -1), omega_radians=0.25
31                )
32            )
33        },
34    )

V2 Client (new ✅)

 1volume_mesh["blk-1"].center=(0, 0, -1.32392) * fl.u.m
 2volume_mesh["blk-1"].axis=[0,0,-1]
 3
 4with fl.SI_unit_system:
 5    params = fl.SimulationParams(
 6        operating_condition=fl.operating_condition_from_mach_reynolds(
 7            mach=0.05, 
 8            reynolds=10, 
 9            temperature=288.15, 
10            alpha=90 * fl.u.deg, 
11            beta=0 * fl.u.deg,
12            project_length_unit = 1 * u.m,
13        ),
14        models=[
15            fl.Fluid(
16                navier_stokes_solver=fl.NavierStokesSolver(
17                    absolute_tolerance=1e-9,
18                    linear_solver=fl.LinearSolver(max_iterations=35),
19                ),
20                turbulence_model_solver=fl.NoneSolver(),
21            ),
22            fl.Wall(surfaces=volume_mesh["blk-1/Cylinder"]),
23            fl.Wall(surfaces=volume_mesh["blk-1/OuterWall"], velocity=[0.0, 0.0, 0.0]),
24            fl.Freestream(entities=volume_mesh["blk-1/InletOutlet"]),
25            fl.Rotation(entities=[volume_mesh["blk-1"]], spec=fl.AngularVelocity(812.31 * fl.u.rpm)),
26        ],
27        time_stepping=fl.Steady(max_steps=2000),
28        outputs=[
29            fl.VolumeOutput(
30                output_format="paraview",
31                output_fields=["primitiveVars", "Mach", "VelocityRelative"],
32            ),
33            fl.SurfaceOutput(
34                entities=[volume_mesh["blk-1/Cylinder"], volume_mesh["blk-1/OuterWall"]],
35                output_format="paraview",
36                output_fields=["primitiveVars", "Cp", "Cf", "yPlus"],
37            ),
38        ],
39    )