PyFluent入门之旅（4）算例求解
在网格划分完成或已有网格的情况下，可以进行算例的求解。
1. 切换/打开求解器
一般启动求解器前有两种情况：
- 已启动FluentMeshing并生成了网格，需要在不退出FluentMeshing的情况下直接切换至Fluent求解器。
- 已经有现成的网格文件，只需要启动Fluent求解器线程并直接读取。
针对需要直接切换的情况，在“PyFluent入门之旅（3）网格划分”的代码基础上继续编写即可，使用下述代码直接切换至求解器，并默认读取之前生成的网格
solver = meshing.switch_to_solver()
对于需要重新编写代码直接启动求解器并读取现有网格的情况。需要重新导入pyfluent库，启动求解器并使用tui读取网格文件。
import ansys.fluent.core as pyfluent
solver = pyfluent.launch_fluent(precision="double", processor_count=4, mode="solver")
solver.tui.file.read_case(case_file_name = 'mixing_elbow.msh.h5')
读取网格后可以顺便检查一下网格参数，其结果会直接输出至python界面
solver.tui.mesh.check()
2. 基本求解设置
进行求解前需要进行一些必要的设置。
首先是单位的修改，可以根据自身需求决定是否修改。例如将长度单位修改为‘in’
solver.tui.define.units("length", "in")
接着根据问题需要开启所需功能，例如该案例需要开启能量方程
solver.tui.define.models.energy("yes", ", ", ", ", ", ", ", ")
创建材料功能也十分常用，该案例我们需要创建液态水
solver.tui.define.materials.copy("fluid", "water-liquid")
将流体域材料改为液态水
solver.tui.define.boundary_conditions.fluid(
"elbow-fluid",
"yes",
"water-liquid",
"no",
"no",
"no",
"no",
"0",
"no",
"0",
"no",
"0",
"no",
"0",
"no",
"0",
"no",
"1",
"no",
"no",
"no",
"no",
"no",
)
修改完成后会输出以下结果
边界条件设置。
solver.tui.define.boundary_conditions.set.velocity_inlet(
"cold-inlet", [], "vmag", "no", 0.4, "quit"
)
solver.tui.define.boundary_conditions.set.velocity_inlet(
"cold-inlet", [], "ke-spec", "no", "no", "no", "yes", "quit"
)
solver.tui.define.boundary_conditions.set.velocity_inlet(
"cold-inlet", [], "turb-intensity", 5, "quit"
)
solver.tui.define.boundary_conditions.set.velocity_inlet(
"cold-inlet", [], "turb-hydraulic-diam", 4, "quit"
)
solver.tui.define.boundary_conditions.set.velocity_inlet(
"cold-inlet", [], "temperature", "no", 293.15, "quit"
)
# hot inlet (hot-inlet), Setting: Value:
# Velocity Specification Method: Magnitude, Normal to Boundary
solver.tui.define.boundary_conditions.set.velocity_inlet(
"hot-inlet", [], "vmag", "no", 1.2, "quit"
)
solver.tui.define.boundary_conditions.set.velocity_inlet(
"hot-inlet", [], "ke-spec", "no", "no", "no", "yes", "quit"
)
solver.tui.define.boundary_conditions.set.velocity_inlet(
"hot-inlet", [], "turb-intensity", 5, "quit"
)
solver.tui.define.boundary_conditions.set.velocity_inlet(
"hot-inlet", [], "turb-hydraulic-diam", 1, "quit"
)
solver.tui.define.boundary_conditions.set.velocity_inlet(
"hot-inlet", [], "temperature", "no", 313.15, "quit"
)
# pressure outlet (outlet), Setting: Value:
# Backflow Turbulent Intensity: 5 [%]
# Backflow Turbulent Viscosity Ratio: 4
solver.tui.define.boundary_conditions.set.pressure_outlet(
"outlet", [], "turb-intensity", 5, "quit"
)
solver.tui.define.boundary_conditions.set.pressure_outlet(
"outlet", [], "turb-viscosity-ratio", 4, "quit"
)
可以设置是否在计算中启用残差图显示，这里选择显示
solver.tui.solve.monitors.residual.plot("yes")
为方便后处理获取数据可以在所需边界创建参数报告
solver.tui.solve.report_definitions.add(
"outlet-temp-avg",
"surface-massavg",
"field",
"temperature",
"surface-names",
"outlet",
"()",
"quit",
)
3. 进行求解
完成必要的设置后，我们可以进行求解了。
首先进行初始化，使用默认的混合初始化。
solver.tui.solve.initialize.hyb_initialization()
接着保存以下算例文件
solver.tui.file.write_case("mixing_elbow1.cas.h5")
先求解100步试试看
solver.tui.solve.iterate(100)
求解完成后保存数据文件，方便后续的查看。
solver.tui.file.write_data("mixing_elbow1.dat.h5")
可以看到计算至93步时已经收敛。
原文地址：https://blog.csdn.net/laughzjh/article/details/140250329
免责声明：本站文章内容转载自网络资源，如本站内容侵犯了原著者的合法权益，可联系本站删除。更多内容请关注自学内容网（zxcms.com）！