.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "examples/00-fluent/conjugate_heat_transfer.py"
.. LINE NUMBERS ARE GIVEN BELOW.
.. only:: html
.. note::
:class: sphx-glr-download-link-note
:ref:`Go to the end <sphx_glr_download_examples_00-fluent_conjugate_heat_transfer.py>`
to download the full example code.
.. rst-class:: sphx-glr-example-title
.. _sphx_glr_examples_00-fluent_conjugate_heat_transfer.py:
.. _conjugate_heat_transfer:
Conjugate Heat Transfer
-----------------------
.. GENERATED FROM PYTHON SOURCE LINES 30-43
Objective
=====================================================================================
This tutorial demonstrates how to model forced convection in a louvered fin heat
exchanger. This case solves equations for both Fluid and Solid domain.
As a result, temperature field evolved together.
This tutorial demonstrates how to perform the following tasks:
* Calculate the fin heat transfer rate.
* Use periodic boundaries to reduce the size of the computational domain.
* Use a convective thermal boundary condition to represent heat transfer.
* Examine and understand the relationship between flow and temperature.
.. GENERATED FROM PYTHON SOURCE LINES 43-45
.. code-block:: Python
:dedent: 1
.. GENERATED FROM PYTHON SOURCE LINES 47-49
Import required libraries/modules
=================================
.. GENERATED FROM PYTHON SOURCE LINES 49-68
.. code-block:: Python
import csv
import os
from pathlib import Path
import matplotlib.pyplot as plt
import pyvista as pv
import ansys.fluent.core as pyfluent
from ansys.fluent.core import examples
from ansys.fluent.visualization.matplotlib import Plots
from ansys.fluent.visualization.pyvista import Graphics, pyvista_windows_manager
from ansys.fluent.visualization.pyvista.pyvista_windows_manager import PyVistaWindow
geom_filename = examples.download_file(
"cht_fin_htc_new.scdoc",
"pyfluent/examples/CHT",
)
.. GENERATED FROM PYTHON SOURCE LINES 69-71
Fluent Solution Setup
=====================
.. GENERATED FROM PYTHON SOURCE LINES 73-75
Launch Fluent session with meshing mode and print Fluent version
================================================================
.. GENERATED FROM PYTHON SOURCE LINES 75-85
.. code-block:: Python
meshing = pyfluent.launch_fluent(
product_version="25.1.0",
mode="meshing",
dimension=3,
precision="double",
processor_count=4,
)
print(meshing.get_fluent_version())
.. GENERATED FROM PYTHON SOURCE LINES 86-88
Start Watertight Geometry Meshing Workflow
==========================================
.. GENERATED FROM PYTHON SOURCE LINES 88-120
.. code-block:: Python
meshing.workflow.InitializeWorkflow(WorkflowType=r"Watertight Geometry")
meshing.workflow.TaskObject["Import Geometry"].Arguments = dict(FileName=geom_filename)
meshing.workflow.TaskObject["Import Geometry"].Execute()
meshing.workflow.TaskObject["Add Local Sizing"].Execute()
meshing.workflow.TaskObject["Generate the Surface Mesh"].Arguments = dict(
{
"CFDSurfaceMeshControls": {
"MinSize": 0.3,
"MaxSize": 1,
"ScopeProximityTo": "faces",
},
}
)
meshing.workflow.TaskObject["Generate the Surface Mesh"].Execute()
meshing.workflow.TaskObject["Describe Geometry"].UpdateChildTasks(SetupTypeChanged=True)
meshing.workflow.TaskObject["Describe Geometry"].Arguments.setState(
{
r"CappingRequired": r"No",
r"InvokeShareTopology": r"No",
r"NonConformal": r"Yes",
r"SetupType": r"The geometry consists of both fluid and solid regions and/or voids",
}
)
meshing.workflow.TaskObject["Describe Geometry"].Execute()
.. GENERATED FROM PYTHON SOURCE LINES 121-123
Update Interface Boundaries; Create Region
==========================================
.. GENERATED FROM PYTHON SOURCE LINES 123-335
.. code-block:: Python
meshing.workflow.TaskObject["Update Boundaries"].Arguments.setState(
{
r"BoundaryLabelList": [
r"interface-out-solid-a",
r"interface-out-high-a",
r"interface-out-low-a",
r"interface-4-solid-sweep",
r"interface-4-high-sweep",
r"interface-4-low-sweep",
r"interface-3-solid-sweep",
r"interface-3-high-sweep",
r"interface-3-low-sweep",
r"interface-2-solid-sweep",
r"interface-2-high-sweep",
r"interface-2-low-sweep",
r"interface-1-solid-sweep",
r"interface-1-high-sweep",
r"interface-1-low-sweep",
r"interface-solid-in-a",
r"interface-in-high-a",
r"interface-in-low-a",
r"interface-tube-2-solid-a",
r"interface-tube-2-high-a",
r"interface-tube-2-low-a",
r"interface-tube-1-solid-a",
r"interface-tube-1-high-a",
r"interface-tube-1-low-a",
r"interface-4-fluid-high-tet",
r"interface-4-fluid-low-tet",
r"interface-3-fluid-low-tet",
r"interface-3-fluid-high-tet",
r"interface-2-fluid-high-tet",
r"interface-2-fluid-low-tet",
r"interface-1-fluid-high-tet",
r"interface-1-fluid-low-tet",
r"interface-1-solid-tet-4",
r"interface-1-solid-tet-3",
r"interface-1-solid-tet-2",
r"interface-1-solid-tet-1",
],
r"BoundaryLabelTypeList": [
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
r"interface",
],
r"OldBoundaryLabelList": [
r"interface-out-solid-a",
r"interface-out-high-a",
r"interface-out-low-a",
r"interface-4-solid-sweep",
r"interface-4-high-sweep",
r"interface-4-low-sweep",
r"interface-3-solid-sweep",
r"interface-3-high-sweep",
r"interface-3-low-sweep",
r"interface-2-solid-sweep",
r"interface-2-high-sweep",
r"interface-2-low-sweep",
r"interface-1-solid-sweep",
r"interface-1-high-sweep",
r"interface-1-low-sweep",
r"interface-solid-in-a",
r"interface-in-high-a",
r"interface-in-low-a",
r"interface-tube-2-solid-a",
r"interface-tube-2-high-a",
r"interface-tube-2-low-a",
r"interface-tube-1-solid-a",
r"interface-tube-1-high-a",
r"interface-tube-1-low-a",
r"interface-4-fluid-high-tet",
r"interface-4-fluid-low-tet",
r"interface-3-fluid-low-tet",
r"interface-3-fluid-high-tet",
r"interface-2-fluid-high-tet",
r"interface-2-fluid-low-tet",
r"interface-1-fluid-high-tet",
r"interface-1-fluid-low-tet",
r"interface-1-solid-tet-4",
r"interface-1-solid-tet-3",
r"interface-1-solid-tet-2",
r"interface-1-solid-tet-1",
],
r"OldBoundaryLabelTypeList": [
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
r"wall",
],
r"OldLabelZoneList": [
r"interface-out-solid-a",
r"interface-out-high-a",
r"interface-out-low-a",
r"interface-4-solid-sweep",
r"interface-4-high-sweep",
r"interface-4-low-sweep",
r"interface-3-solid-sweep",
r"interface-3-high-sweep",
r"interface-3-low-sweep",
r"interface-2-solid-sweep",
r"interface-2-high-sweep",
r"interface-2-low-sweep",
r"interface-1-solid-sweep",
r"interface-1-high-sweep",
r"interface-1-low-sweep",
r"interface-solid-in-a",
r"interface-in-high-a",
r"interface-in-low-a",
r"interface-tube-2-solid-a.2",
r"interface-tube-2-solid-a.1",
r"interface-tube-2-solid-a",
r"interface-tube-2-high-a.2",
r"interface-tube-2-high-a.1",
r"interface-tube-2-high-a",
r"interface-tube-2-low-a.2",
r"interface-tube-2-low-a.1",
r"interface-tube-2-low-a",
r"interface-tube-1-solid-a.2",
r"interface-tube-1-solid-a.1",
r"interface-tube-1-solid-a",
r"interface-tube-1-high-a.2",
r"interface-tube-1-high-a.1",
r"interface-tube-1-high-a",
r"interface-tube-1-low-a.2",
r"interface-tube-1-low-a.1",
r"interface-tube-1-low-a",
r"interface-4-fluid-high-tet",
r"interface-4-fluid-low-tet",
r"interface-3-fluid-low-tet",
r"interface-3-fluid-high-tet",
r"interface-2-fluid-high-tet",
r"interface-2-fluid-low-tet",
r"interface-1-fluid-high-tet",
r"interface-1-fluid-low-tet",
r"interface-1-solid-tet-4",
r"interface-1-solid-tet-3",
r"interface-1-solid-tet-2",
r"interface-1-solid-tet-1",
],
}
)
meshing.workflow.TaskObject["Update Boundaries"].Execute()
meshing.workflow.TaskObject["Create Regions"].Execute()
.. GENERATED FROM PYTHON SOURCE LINES 336-338
Custom Journal for Creating Periodicity due to Non-Conformal Objects
====================================================================
.. GENERATED FROM PYTHON SOURCE LINES 338-355
.. code-block:: Python
meshing.workflow.TaskObject["Describe Geometry"].InsertNextTask(
CommandName=r"RunCustomJournal"
)
meshing.workflow.TaskObject["Run Custom Journal"].Rename(NewName=r"set-periodicity")
meshing.workflow.TaskObject["set-periodicity"].Arguments = dict(
{
r"JournalString": r"""/bo rps translational semi-auto periodic-1-high periodic-2-high periodic-3-high periodic-4-high , 0 0 -2.3
/bo rps translational semi-auto periodic-5* , 0 0 -2.3
/bo rps translational semi-auto periodic-6-high , 0 0 -2.3
/bo rps translational semi-auto periodic-7-high , 0 0 -2.3
""",
}
)
meshing.workflow.TaskObject["set-periodicity"].Execute()
.. GENERATED FROM PYTHON SOURCE LINES 356-358
Update Boundary Layer Task
==========================
.. GENERATED FROM PYTHON SOURCE LINES 358-411
.. code-block:: Python
meshing.workflow.TaskObject["Update Regions"].Execute()
meshing.workflow.TaskObject["Add Boundary Layers"].AddChildToTask()
meshing.workflow.TaskObject["Add Boundary Layers"].InsertCompoundChildTask()
meshing.workflow.TaskObject["smooth-transition_1"].Rename(NewName=r"aspect-ratio_1")
meshing.workflow.TaskObject["aspect-ratio_1"].Arguments.setState(
{
"BLControlName": r"aspect-ratio_1",
"BLRegionList": [
r"fluid-tet-4",
r"fluid-tet-3",
r"fluid-tet-2",
r"fluid-tet-1",
r"fluid-sweep-fin2",
r"fluid-sweep-fin1",
r"fluid-sweep-fin5",
r"fluid-sweep-fin3",
r"fluid-sweep-fin6",
r"fluid-sweep-fin4",
r"fluid-in",
r"fluid-out",
],
r"BLZoneList": [
r"wall-fluid-tet-4-solid-tet-4",
r"wall-fluid-tet-3-solid-tet-3",
r"wall-fluid-tet-2-solid-tet-2",
r"wall-fluid-tet-2-solid-tet-2-wall-fluid-tet-3-solid-tet-3-fluid-tet-2-solid-tet-2",
r"wall-fluid-tet-1-solid-tet-1",
r"wall-fluid-sweep-fin-solid-sweep-fin.1",
r"wall-fluid-sweep-fin-solid-sweep-fin",
r"wall-fluid-sweep-fin-solid-sweep-fin.5",
r"wall-fluid-sweep-fin-solid-sweep-fin.4",
r"wall-fluid-sweep-fin-solid-sweep-fin.3",
r"wall-fluid-sweep-fin-solid-sweep-fin.2",
],
r"BlLabelList": r"wall*",
r"CompleteBlLabelList": [
r"wall-fluid-sweep-fin-solid-sweep-fin",
r"wall-fluid-tet-1-solid-tet-1",
r"wall-fluid-tet-2-solid-tet-2",
r"wall-fluid-tet-3-solid-tet-3",
r"wall-fluid-tet-4-solid-tet-4",
],
r"FaceScope": {
r"GrowOn": r"selected-labels",
},
r"OffsetMethodType": r"aspect-ratio",
}
)
meshing.workflow.TaskObject["aspect-ratio_1"].Execute()
.. GENERATED FROM PYTHON SOURCE LINES 412-414
Generate Mesh
=============
.. GENERATED FROM PYTHON SOURCE LINES 414-417
.. code-block:: Python
meshing.workflow.TaskObject["Generate the Volume Mesh"].Execute()
.. GENERATED FROM PYTHON SOURCE LINES 418-420
Improve Volume Mesh
===================
.. GENERATED FROM PYTHON SOURCE LINES 420-439
.. code-block:: Python
meshing.workflow.TaskObject["Generate the Volume Mesh"].InsertNextTask(
CommandName=r"ImproveVolumeMesh"
)
meshing.workflow.TaskObject["Improve Volume Mesh"].Arguments.setState(
{
r"CellQualityLimit": 0.05,
r"VMImprovePreferences": {
r"ShowVMImprovePreferences": False,
r"VIQualityIterations": 5,
r"VIQualityMinAngle": 0,
r"VIgnoreFeature": r"yes",
},
}
)
meshing.workflow.TaskObject["Improve Volume Mesh"].Execute()
.. GENERATED FROM PYTHON SOURCE LINES 440-442
Save Mesh File
==============
.. GENERATED FROM PYTHON SOURCE LINES 442-446
.. code-block:: Python
save_mesh_as = str(Path(pyfluent.EXAMPLES_PATH) / "hx-fin-2mm.msh.h5")
meshing.tui.file.write_mesh(save_mesh_as)
.. GENERATED FROM PYTHON SOURCE LINES 447-449
Switch to Solution Mode
=======================
.. GENERATED FROM PYTHON SOURCE LINES 449-452
.. code-block:: Python
solver = meshing.switch_to_solver()
.. GENERATED FROM PYTHON SOURCE LINES 453-455
Auto-create Mesh Interfaces
===========================
.. GENERATED FROM PYTHON SOURCE LINES 455-458
.. code-block:: Python
solver.tui.define.mesh_interfaces.create("int", "yes", "no")
.. GENERATED FROM PYTHON SOURCE LINES 459-461
Mesh Check; Review Fluent transcript for errors
===============================================
.. GENERATED FROM PYTHON SOURCE LINES 461-464
.. code-block:: Python
solver.mesh.check()
.. GENERATED FROM PYTHON SOURCE LINES 465-467
Create a few boundary list for display and post-processing
==========================================================
.. GENERATED FROM PYTHON SOURCE LINES 467-483
.. code-block:: Python
graphics_session1 = Graphics(solver)
mesh1 = graphics_session1.Meshes["mesh-1"]
wall_list = []
periodic_list = []
symmetry_list = []
for item in mesh1.surfaces_list.allowed_values:
if len(item.split("wall")) > 1:
wall_list.append(item)
if len(item.split("periodic")) > 1:
periodic_list.append(item)
if len(item.split("symmetry")) > 1:
symmetry_list.append(item)
.. GENERATED FROM PYTHON SOURCE LINES 484-486
Display Mesh
============
.. GENERATED FROM PYTHON SOURCE LINES 486-497
.. code-block:: Python
mesh1.show_edges = True
mesh1.surfaces_list = wall_list
mesh1.display("window-1")
p = pyvista_windows_manager.get_plotter("window-1")
p.view_isometric()
p.add_axes()
p.add_floor(offset=1, show_edges=False)
light = pv.Light(light_type="headlight")
p.add_light(light)
.. GENERATED FROM PYTHON SOURCE LINES 498-503
Temperature, Energy, Laminar Viscous Model
==========================================
* Set Temperature Unit
* Enable Energy Equation
* Enable Laminar Viscous Model
.. GENERATED FROM PYTHON SOURCE LINES 503-510
.. code-block:: Python
solver.setup.general.units.set_units(
quantity="temperature", units_name="c", scale_factor=1.0, offset=0.0
)
solver.setup.models.energy.enabled = True
solver.setup.models.viscous.model.set_state("laminar")
.. GENERATED FROM PYTHON SOURCE LINES 511-513
Change a few material properties of default Air
===============================================
.. GENERATED FROM PYTHON SOURCE LINES 513-521
.. code-block:: Python
air_dict = solver.setup.materials.fluid["air"].get_state()
air_dict["density"]["value"] = 1.2
air_dict["viscosity"]["value"] = 1.5e-5
air_dict["thermal_conductivity"]["value"] = 0.026
air_dict["specific_heat"]["value"] = 1006.0
solver.setup.materials.fluid["air"].set_state(air_dict)
.. GENERATED FROM PYTHON SOURCE LINES 522-524
Change a few material properties of default Aluminum
====================================================
.. GENERATED FROM PYTHON SOURCE LINES 524-531
.. code-block:: Python
al_dict = solver.setup.materials.solid["aluminum"].get_state()
al_dict["density"]["value"] = 2719.0
al_dict["thermal_conductivity"]["value"] = 200.0
al_dict["specific_heat"]["value"] = 871.0
solver.setup.materials.solid["aluminum"].set_state(al_dict)
.. GENERATED FROM PYTHON SOURCE LINES 532-534
Copy Copper and change a few material properties of default Copper
==================================================================
.. GENERATED FROM PYTHON SOURCE LINES 534-542
.. code-block:: Python
solver.setup.materials.database.copy_by_name(type="solid", name="copper")
cu_dict = solver.setup.materials.solid["copper"].get_state()
cu_dict["density"]["value"] = 8978.0
cu_dict["thermal_conductivity"]["value"] = 340.0
cu_dict["specific_heat"]["value"] = 381.0
solver.setup.materials.solid["copper"].set_state(cu_dict)
.. GENERATED FROM PYTHON SOURCE LINES 543-545
Set Tube Cell Zone Material as Copper
=====================================
.. GENERATED FROM PYTHON SOURCE LINES 545-554
.. code-block:: Python
tube_dict = solver.setup.cell_zone_conditions.solid["solid-tube-1"].get_state()
tube_dict["material"] = "copper"
solver.setup.cell_zone_conditions.solid["solid-tube-1"].set_state(tube_dict)
tube_dict = solver.setup.cell_zone_conditions.solid["solid-tube-2"].get_state()
tube_dict["material"] = "copper"
solver.setup.cell_zone_conditions.solid["solid-tube-2"].set_state(tube_dict)
.. GENERATED FROM PYTHON SOURCE LINES 555-557
Set Boundary Condition for Inlet and Outlet
===========================================
.. GENERATED FROM PYTHON SOURCE LINES 557-567
.. code-block:: Python
solver.setup.boundary_conditions.velocity_inlet["inlet"].momentum.velocity = 4.0
solver.setup.boundary_conditions.velocity_inlet["inlet"].thermal.temperature = (
293.15 # Need to specify in Kelvin
)
solver.setup.boundary_conditions.pressure_outlet[
"outlet"
].thermal.backflow_total_temperature = 293.15
.. GENERATED FROM PYTHON SOURCE LINES 568-570
Set Thermal Boundary Condition for Wall Inner Tube
==================================================
.. GENERATED FROM PYTHON SOURCE LINES 570-583
.. code-block:: Python
solver.setup.boundary_conditions.wall["wall-inner-tube-1"].thermal.thermal_condition = (
"Convection"
)
solver.setup.boundary_conditions.wall[
"wall-inner-tube-1"
].thermal.heat_transfer_coeff = 1050.0
solver.setup.boundary_conditions.wall["wall-inner-tube-1"].thermal.free_stream_temp = (
353.15
)
solver.setup.boundary_conditions.copy(from_="wall-inner-tube-1", to="wall-inner-tube-2")
.. GENERATED FROM PYTHON SOURCE LINES 584-586
Enable HOTR
===========
.. GENERATED FROM PYTHON SOURCE LINES 586-589
.. code-block:: Python
solver.solution.methods.high_order_term_relaxation.enable = True
.. GENERATED FROM PYTHON SOURCE LINES 590-592
Define Report Definitions
=========================
.. GENERATED FROM PYTHON SOURCE LINES 592-677
.. code-block:: Python
solver.solution.report_definitions.surface["outlet-enthalpy-flow"] = {}
solver.solution.report_definitions.surface["outlet-enthalpy-flow"].report_type = (
"surface-flowrate"
)
solver.solution.report_definitions.surface["outlet-enthalpy-flow"].field = "enthalpy"
solver.solution.report_definitions.surface["outlet-enthalpy-flow"].surface_names = [
"outlet"
]
solver.solution.report_definitions.surface["avg-pressure-inlet"] = {}
solver.solution.report_definitions.surface["avg-pressure-inlet"].report_type = (
"surface-areaavg"
)
solver.solution.report_definitions.surface["avg-pressure-inlet"].field = "pressure"
solver.solution.report_definitions.surface["avg-pressure-inlet"].surface_names = [
"inlet"
]
solver.solution.report_definitions.volume["max-vel-louvers4"] = {}
solver.solution.report_definitions.volume["max-vel-louvers4"].report_type = "volume-max"
solver.solution.report_definitions.volume["max-vel-louvers4"].field = (
"velocity-magnitude"
)
solver.solution.report_definitions.volume["max-vel-louvers4"].cell_zones = [
"fluid-tet-4"
]
solver.solution.report_definitions.surface["wall-shear-int"] = {}
solver.solution.report_definitions.surface["wall-shear-int"].report_type = (
"surface-integral"
)
solver.solution.report_definitions.surface["wall-shear-int"].field = "wall-shear"
solver.solution.report_definitions.surface["wall-shear-int"].surface_names = [
"wall-fluid-sweep-fin-solid-sweep-fin-shadow",
"wall-fluid-tet-1-solid-tet-1",
"wall-fluid-tet-2-solid-tet-2",
"wall-fluid-tet-3-solid-tet-3",
"wall-fluid-tet-4-solid-tet-4",
]
solver.solution.monitor.report_plots.create(name="outlet-enthalpy-flow-plot")
solver.solution.monitor.report_plots["outlet-enthalpy-flow-plot"].report_defs = (
"outlet-enthalpy-flow"
)
solver.solution.monitor.report_files["outlet-enthalpy-flow-file"] = {}
solver.solution.monitor.report_files["outlet-enthalpy-flow-file"] = {
"report_defs": ["outlet-enthalpy-flow"],
"file_name": r"outlet-enthalpy-flow.out",
}
solver.solution.monitor.report_plots["avg-pressure-inlet-plot"] = {}
solver.solution.monitor.report_plots["avg-pressure-inlet-plot"] = {
"report_defs": ["avg-pressure-inlet"]
}
solver.solution.monitor.report_files["avg-pressure-inlet-file"] = {}
solver.solution.monitor.report_files["avg-pressure-inlet-file"] = {
"report_defs": ["avg-pressure-inlet"],
"file_name": r"avg-pressure-inlet.out",
}
solver.solution.monitor.report_plots["max-vel-louvers4-plot"] = {}
solver.solution.monitor.report_plots["max-vel-louvers4-plot"] = {
"report_defs": ["max-vel-louvers4"]
}
solver.solution.monitor.report_files["max-vel-louvers4-file"] = {}
solver.solution.monitor.report_files["max-vel-louvers4-file"] = {
"report_defs": ["max-vel-louvers4"],
"file_name": r"max-vel-louvers4.out",
}
solver.solution.monitor.report_plots["wall-shear-int-plot"] = {}
solver.solution.monitor.report_plots["wall-shear-int-plot"] = {
"report_defs": ["wall-shear-int"]
}
solver.solution.monitor.report_files["wall-shear-int-file"] = {}
solver.solution.monitor.report_files["wall-shear-int-file"] = {
"report_defs": ["wall-shear-int"],
"file_name": r"wall-shear-int.out",
}
.. GENERATED FROM PYTHON SOURCE LINES 678-680
Hybrid Initialization; Slit Interior between Solid Zones; Save Case
===================================================================
.. GENERATED FROM PYTHON SOURCE LINES 680-687
.. code-block:: Python
solver.solution.initialization.initialization_type = "hybrid"
solver.solution.initialization.hybrid_initialize()
solver.setup.boundary_conditions.slit_interior_between_diff_solids()
solver.file.write(file_type="case", file_name="hx-fin-2mm.cas.h5")
.. GENERATED FROM PYTHON SOURCE LINES 688-690
Set Aggressive Length Scale Method; Run Calculation & Save Data
===============================================================
.. GENERATED FROM PYTHON SOURCE LINES 690-702
.. code-block:: Python
solver.solution.run_calculation.pseudo_time_settings.time_step_method.time_step_method = (
"automatic"
)
solver.solution.run_calculation.pseudo_time_settings.time_step_method.length_scale_methods = (
"aggressive"
)
solver.solution.run_calculation.iterate(iter_count=250)
solver.file.write(file_type="case-data", file_name="hx-fin-2mm.dat.h5")
.. GENERATED FROM PYTHON SOURCE LINES 703-705
Post-Processing Mass Balance Report
===================================
.. GENERATED FROM PYTHON SOURCE LINES 705-720
.. code-block:: Python
inlet_mfr = solver.scheme_eval.exec(
('(ti-menu-load-string "/report/fluxes/mass-flow no inlet () no")',)
).split(" ")[-1]
outlet_mfr = solver.scheme_eval.exec(
('(ti-menu-load-string "/report/fluxes/mass-flow no outlet () no")',)
).split(" ")[-1]
net_mfr = solver.scheme_eval.exec(
('(ti-menu-load-string "/report/fluxes/mass-flow no inlet outlet () no")',)
).split(" ")[-1]
print("Mass Balance Report\n")
print("Inlet (kg/s): ", inlet_mfr)
print("Outlet (kg/s): ", outlet_mfr)
print("Net (kg/s): ", net_mfr)
.. GENERATED FROM PYTHON SOURCE LINES 721-723
Heat Balance Report
===================
.. GENERATED FROM PYTHON SOURCE LINES 723-730
.. code-block:: Python
htr = solver.scheme_eval.exec(
('(ti-menu-load-string "/report/fluxes/heat-transfer yes no")',)
).split(" ")[-1]
print("Heat Balance Report\n")
print("Net Imbalance (Watt): ", htr)
.. GENERATED FROM PYTHON SOURCE LINES 731-733
Plot Monitors
=============
.. GENERATED FROM PYTHON SOURCE LINES 733-766
.. code-block:: Python
fig, axs = plt.subplots(2, 2, figsize=(10, 8))
fig.suptitle("Monitor Plots")
outFilesList = []
fileList = os.listdir(os.getcwd())
for tempFile in fileList:
fName, ext = os.path.splitext(tempFile)
if ext == ".out":
outFilesList.append(tempFile)
outFilesList.sort()
index = 0
for ax in axs.flat:
X = []
Y = []
i = -1
with open(outFilesList[index], "r") as datafile:
plotting = csv.reader(datafile, delimiter=" ")
for rows in plotting:
i += 1
if i == 1:
var = rows[1]
if i > 2:
X.append(int(rows[0]))
Y.append(float(rows[1]))
ax.plot(X, Y)
ax.set(xlabel="Iteration", ylabel=var, title=var)
index += 1
plt.tight_layout()
.. GENERATED FROM PYTHON SOURCE LINES 767-769
Show graph
==========
.. GENERATED FROM PYTHON SOURCE LINES 769-772
.. code-block:: Python
plt.show()
.. GENERATED FROM PYTHON SOURCE LINES 773-776
.. image:: ../../_static/cht_avg_pressure.png
:align: center
:alt: Average Pressure
.. GENERATED FROM PYTHON SOURCE LINES 778-779
Average Pressure
.. GENERATED FROM PYTHON SOURCE LINES 781-783
Contour Plot
============
.. GENERATED FROM PYTHON SOURCE LINES 783-812
.. code-block:: Python
graphics_session1 = Graphics(solver)
contour1 = graphics_session1.Contours["contour-1"]
contour1.field = "temperature"
contour1.surfaces_list = wall_list
contour1.display("window-2")
p = pyvista_windows_manager.get_plotter("window-2")
p.view_isometric()
p.add_axes()
p.add_floor(offset=1, show_edges=False)
p.add_title(
"Contour of Temperature on Walls", font="courier", color="grey", font_size=10
)
light = pv.Light(light_type="headlight")
p.add_light(light)
p.remove_scalar_bar()
p.add_scalar_bar(
"Temperature [K]",
interactive=True,
vertical=False,
title_font_size=20,
label_font_size=15,
outline=False,
position_x=0.5,
fmt="%10.1f",
)
.. GENERATED FROM PYTHON SOURCE LINES 813-815
Create Iso-Surface of X=0.012826 m
==================================
.. GENERATED FROM PYTHON SOURCE LINES 815-820
.. code-block:: Python
solver.results.surfaces.iso_surface["x=0.012826"] = {}
solver.results.surfaces.iso_surface["x=0.012826"].field = "x-coordinate"
solver.results.surfaces.iso_surface["x=0.012826"] = {"iso_values": [0.012826]}
.. GENERATED FROM PYTHON SOURCE LINES 821-823
Vecotor Plot
============
.. GENERATED FROM PYTHON SOURCE LINES 823-855
.. code-block:: Python
graphics_session1 = Graphics(solver)
vector1 = graphics_session1.Vectors["vector-1"]
vector1.surfaces_list = ["x=0.012826"]
vector1.scale = 2.0
vector1.skip = 5
vector1.display("window-3")
p = pyvista_windows_manager.get_plotter("window-3")
p.view_isometric()
p.add_axes()
# p.add_floor( offset=1, show_edges=False)
p.add_title("Vector Plot", font="courier", color="grey", font_size=10)
light = pv.Light(light_type="headlight")
p.add_light(light)
p.remove_scalar_bar()
p.add_scalar_bar(
"Velocity [m/s]",
interactive=True,
vertical=False,
title_font_size=20,
label_font_size=15,
outline=False,
position_x=0.5,
fmt="%10.1f",
)
o = PyVistaWindow(None, None)
o._fetch_mesh(mesh1)
o._display_mesh(mesh1, p)
.. GENERATED FROM PYTHON SOURCE LINES 856-858
XY Plot of Pressure
===================
.. GENERATED FROM PYTHON SOURCE LINES 858-866
.. code-block:: Python
plots_session1 = Plots(solver)
p1 = plots_session1.XYPlots["p1"]
p1.surfaces_list = ["x=0.012826"]
p1.y_axis_function = "pressure"
p1.x_axis_function = "direction-vector"
p1.direction_vector.set_state([0, 1, 0])
.. GENERATED FROM PYTHON SOURCE LINES 867-869
Show graph
==========
.. GENERATED FROM PYTHON SOURCE LINES 869-872
.. code-block:: Python
p1.plot("p1")
.. GENERATED FROM PYTHON SOURCE LINES 873-876
.. image:: ../../_static/cht_xy_pressure.png
:align: center
:alt: XY Plot of Pressure
.. GENERATED FROM PYTHON SOURCE LINES 878-879
XY Plot of Pressure
.. GENERATED FROM PYTHON SOURCE LINES 881-883
Exit Fluent Session
===================
.. GENERATED FROM PYTHON SOURCE LINES 883-884
.. code-block:: Python
solver.exit()
.. _sphx_glr_download_examples_00-fluent_conjugate_heat_transfer.py:
.. only:: html
.. container:: sphx-glr-footer sphx-glr-footer-example
.. container:: sphx-glr-download sphx-glr-download-jupyter
:download:`Download Jupyter notebook: conjugate_heat_transfer.ipynb <conjugate_heat_transfer.ipynb>`
.. container:: sphx-glr-download sphx-glr-download-python
:download:`Download Python source code: conjugate_heat_transfer.py <conjugate_heat_transfer.py>`
.. container:: sphx-glr-download sphx-glr-download-zip
:download:`Download zipped: conjugate_heat_transfer.zip <conjugate_heat_transfer.zip>`
.. only:: html
.. rst-class:: sphx-glr-signature
`Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_