.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "examples/00-fluent/brake.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_brake.py>`
to download the full example code.
.. rst-class:: sphx-glr-example-title
.. _sphx_glr_examples_00-fluent_brake.py:
#################################
Automotive Brake Thermal Analysis
#################################
Objective:
==========
Braking surfaces get heated due to frictional heating during braking.
High temperature affects the braking performance and life of the braking system.
This example demonstrates:
* Fluent setup and simulation using PyFluent
* Post processing using PyVista (3D Viewer) and Matplotlib (2D graphs)
.. GENERATED FROM PYTHON SOURCE LINES 41-43
Import required libraries/modules
==================================================================================
.. GENERATED FROM PYTHON SOURCE LINES 43-57
.. code-block:: Python
import csv
import matplotlib.pyplot as plt
import ansys.fluent.core as pyfluent
from ansys.fluent.core import examples
import ansys.fluent.visualization.pyvista as pv
import_filename = examples.download_file(
"brake.msh.h5",
"pyfluent/examples/Brake-Thermal-PyVista-Matplotlib",
)
.. GENERATED FROM PYTHON SOURCE LINES 58-60
Fluent Solution Setup
==================================================================================
.. GENERATED FROM PYTHON SOURCE LINES 62-64
Launch Fluent session with solver mode and print Fluent version
---------------------------------------------------------------
.. GENERATED FROM PYTHON SOURCE LINES 64-68
.. code-block:: Python
session = pyfluent.launch_fluent(precision="double", processor_count=2, version="3d")
print(session.get_fluent_version())
.. GENERATED FROM PYTHON SOURCE LINES 69-71
Import mesh
------------
.. GENERATED FROM PYTHON SOURCE LINES 71-74
.. code-block:: Python
session.settings.file.read_case(file_name=import_filename)
.. GENERATED FROM PYTHON SOURCE LINES 75-77
Define models and material
--------------------------
.. GENERATED FROM PYTHON SOURCE LINES 77-81
.. code-block:: Python
session.settings.setup.models.energy = {"enabled": True}
session.settings.setup.general.solver.time = "unsteady-2nd-order-bounded"
session.tui.define.materials.copy("solid", "steel")
.. GENERATED FROM PYTHON SOURCE LINES 82-84
Solve only energy equation (conduction)
---------------------------------------
.. GENERATED FROM PYTHON SOURCE LINES 84-86
.. code-block:: Python
session.tui.solve.set.equations("flow", "no", "kw", "no")
.. GENERATED FROM PYTHON SOURCE LINES 87-91
Define disc rotation
--------------------
(15.79 rps corresponds to 100 km/h car speed
with 0.28 m of axis height from ground)
.. GENERATED FROM PYTHON SOURCE LINES 91-114
.. code-block:: Python
session.settings.setup.cell_zone_conditions.solid["disc2"] = {
"solid_motion": {
"solid_motion_zone_motion_function": "none",
"solid_motion_axis_direction": [0, 1, 0],
"solid_motion_axis_origin": [-0.035, -0.821, 0.045],
"solid_motion_velocity": [0, 0, 0],
"solid_omega": -15.79,
"solid_relative_to_thread": "absolute",
"enable": True,
}
}
session.settings.setup.cell_zone_conditions.solid["disc1"] = {
"solid_motion": {
"solid_motion_zone_motion_function": "none",
"solid_motion_axis_direction": [0, 1, 0],
"solid_motion_axis_origin": [-0.035, -0.821, 0.045],
"solid_motion_velocity": [0, 0, 0],
"solid_omega": -15.79,
"solid_relative_to_thread": "absolute",
"enable": True,
}
}
.. GENERATED FROM PYTHON SOURCE LINES 115-119
Apply frictional heating on pad-disc surfaces
----------------------------------------------
Wall thickness 0f 2 mm has been assumed and 2e9 w/m3 is the heat generation which
has been calculated from kinetic energy change due to braking.
.. GENERATED FROM PYTHON SOURCE LINES 119-127
.. code-block:: Python
session.settings.setup.boundary_conditions.wall["wall-pad-disc2"] = {
"thermal": {"q_dot": {"value": 2000000000}, "wall_thickness": {"value": 0.002}}
}
session.settings.setup.boundary_conditions.wall["wall_pad-disc1"] = {
"thermal": {"q_dot": {"value": 2000000000}, "wall_thickness": {"value": 0.002}}
}
.. GENERATED FROM PYTHON SOURCE LINES 128-132
Apply convection cooling on outer surfaces due to air flow
-----------------------------------------------------------
Outer surfaces are applied a constant htc of 100 W/(m2 K)
and 300 K free stream temperature
.. GENERATED FROM PYTHON SOURCE LINES 132-146
.. code-block:: Python
session.tui.define.boundary_conditions.set.wall(
"wall-disc*",
"wall-geom*",
"()",
"thermal-bc",
"yes",
"convection",
"convective-heat-transfer-coefficient",
"no",
100,
"q",
)
.. GENERATED FROM PYTHON SOURCE LINES 147-150
Initialize
----------
Initialize with 300 K temperature
.. GENERATED FROM PYTHON SOURCE LINES 150-154
.. code-block:: Python
session.settings.solution.initialization.initialization_type = "standard"
session.settings.solution.initialization.standard_initialize()
.. GENERATED FROM PYTHON SOURCE LINES 155-161
Post processing setup
---------------------
* Report definitions and monitor plots
* Set contour plot properties
* Set views and camera
* Set animation object
.. GENERATED FROM PYTHON SOURCE LINES 161-270
.. code-block:: Python
session.settings.solution.report_definitions.volume["max-pad-temperature"] = {}
session.settings.solution.report_definitions.volume[
"max-pad-temperature"
].report_type = "volume-max"
session.settings.solution.report_definitions.volume["max-pad-temperature"] = {
"field": "temperature",
"cell_zones": ["geom-1-innerpad", "geom-1-outerpad"],
}
session.settings.solution.report_definitions.volume["max-disc-temperature"] = {}
session.settings.solution.report_definitions.volume[
"max-disc-temperature"
].report_type = "volume-max"
session.settings.solution.report_definitions.volume["max-disc-temperature"] = {
"field": "temperature",
"cell_zones": ["disc1", "disc2"],
}
session.settings.solution.monitor.report_plots.create(name="max-temperature")
session.settings.solution.monitor.report_plots["max-temperature"] = {
"report_defs": ["max-pad-temperature", "max-disc-temperature"]
}
session.settings.solution.monitor.report_files.create(name="max-temperature")
session.settings.solution.monitor.report_files["max-temperature"] = {
"report_defs": ["max-pad-temperature", "max-disc-temperature"],
"file_name": "max-temperature.out",
}
session.settings.solution.monitor.report_files["max-temperature"].report_defs = [
"max-pad-temperature",
"max-disc-temperature",
"flow-time",
]
session.settings.results.graphics.contour.create(name="contour-1")
session.settings.results.graphics.contour["contour-1"] = {
"boundary_values": True,
"range_option": {"auto_range_on": {"global_range": True}},
"field": "temperature",
"draw_mesh": False,
"coloring": {"smooth": False},
"color_map": {
"user_skip": 9,
"log_scale": False,
"visible": True,
"width": 6,
"show_all": True,
"font_name": "Helvetica",
"font_size": 0.032,
"font_automatic": True,
"length": 0.54,
"size": 100,
"format": "%0.2e",
"position": 1,
"color": "field-velocity",
},
"mesh_object": "",
"node_values": True,
"contour_lines": False,
"display_state_name": "None",
"filled": True,
}
session.settings.results.graphics.contour["temperature"] = {
"field": "temperature",
"surfaces_list": "wall*",
"color_map": {
"visible": True,
"size": 100,
"color": "field-velocity",
"log_scale": False,
"format": "%0.1f",
"user_skip": 9,
"show_all": True,
"position": 1,
"font_name": "Helvetica",
"font_automatic": True,
"font_size": 0.032,
"length": 0.54,
"width": 6,
"bground_transparent": True,
"bground_color": "#CCD3E2",
"title_elements": "Variable and Object Name",
},
}
session.settings.results.graphics.contour["temperature"].range_option.option = (
"auto-range-off"
)
session.settings.results.graphics.contour["temperature"].range_option.set_state(
{
"auto_range_off": {"maximum": 400.0, "minimum": 300, "clip_to_range": False},
}
)
session.settings.results.graphics.views.restore_view(view_name="top")
session.settings.results.graphics.views.camera.zoom(factor=2)
session.settings.results.graphics.views.save_view(view_name="animation-view")
session.settings.solution.calculation_activity.solution_animations[
"animate-temperature"
] = {
"animate_on": "temperature",
"frequency_of": "flow-time",
"flow_time_frequency": 0.05,
"view": "animation-view",
}
.. GENERATED FROM PYTHON SOURCE LINES 271-276
Run simulation
---------------
* Run simulation for 2 seconds flow time
* Set time step size
* Set number of time steps and maximum number of iterations per time step
.. GENERATED FROM PYTHON SOURCE LINES 276-282
.. code-block:: Python
session.settings.solution.run_calculation.transient_controls.time_step_size = 0.01
session.settings.solution.run_calculation.dual_time_iterate(
time_step_count=200, max_iter_per_step=5
)
.. GENERATED FROM PYTHON SOURCE LINES 283-286
Save simulation data
--------------------
Write case and data files
.. GENERATED FROM PYTHON SOURCE LINES 286-288
.. code-block:: Python
session.settings.file.write(file_type="case-data", file_name="brake-final.cas.h5")
.. GENERATED FROM PYTHON SOURCE LINES 289-291
Post processing with PyVista (3D visualization)
===============================================
.. GENERATED FROM PYTHON SOURCE LINES 293-295
Create a graphics session
-------------------------
.. GENERATED FROM PYTHON SOURCE LINES 295-297
.. code-block:: Python
graphics_session1 = pv.Graphics(session)
.. GENERATED FROM PYTHON SOURCE LINES 298-300
Temperature contour object
--------------------------
.. GENERATED FROM PYTHON SOURCE LINES 300-302
.. code-block:: Python
contour1 = graphics_session1.Contours["temperature"]
.. GENERATED FROM PYTHON SOURCE LINES 303-305
Check available options for contour object
-------------------------------------------
.. GENERATED FROM PYTHON SOURCE LINES 305-308
.. code-block:: Python
contour1()
.. GENERATED FROM PYTHON SOURCE LINES 309-311
Set contour properties
----------------------
.. GENERATED FROM PYTHON SOURCE LINES 311-328
.. code-block:: Python
contour1.field = "temperature"
contour1.surfaces_list = [
"wall-disc1",
"wall-disc2",
"wall-pad-disc2",
"wall_pad-disc1",
"wall-geom-1-bp_inner",
"wall-geom-1-bp_outer",
"wall-geom-1-innerpad",
"wall-geom-1-outerpad",
]
contour1.range.option = "auto-range-off"
contour1()
contour1.range.auto_range_off.minimum = 300
contour1.range.auto_range_off.maximum = 400
.. GENERATED FROM PYTHON SOURCE LINES 329-331
Display contour
---------------
.. GENERATED FROM PYTHON SOURCE LINES 331-334
.. code-block:: Python
contour1.display()
.. GENERATED FROM PYTHON SOURCE LINES 335-338
.. image:: ../../_static/brake_surface_temperature.png
:align: center
:alt: Brake Surface Temperature Contour
.. GENERATED FROM PYTHON SOURCE LINES 340-341
Brake Surface Temperature
.. GENERATED FROM PYTHON SOURCE LINES 343-345
Post processing with Matplotlib (2D graph)
===============================================
.. GENERATED FROM PYTHON SOURCE LINES 347-349
Read monitor file
-----------------
.. GENERATED FROM PYTHON SOURCE LINES 349-363
.. code-block:: Python
X = []
Y = []
Z = []
i = -1
with open("max-temperature.out", "r") as datafile:
plotting = csv.reader(datafile, delimiter=" ")
for rows in plotting:
i = i + 1
if i > 2:
X.append(float(rows[3]))
Y.append(float(rows[2]))
Z.append(float(rows[1]))
.. GENERATED FROM PYTHON SOURCE LINES 364-366
Plot graph
----------
.. GENERATED FROM PYTHON SOURCE LINES 366-375
.. code-block:: Python
plt.title("Maximum Temperature", fontdict={"color": "darkred", "size": 20})
plt.plot(X, Z, label="Max. Pad Temperature", color="red")
plt.plot(X, Y, label="Max. Disc Temperature", color="blue")
plt.xlabel("Time (sec)")
plt.ylabel("Max Temperature (K)")
plt.legend(loc="lower right", shadow=True, fontsize="x-large")
.. GENERATED FROM PYTHON SOURCE LINES 376-378
Show graph
----------
.. GENERATED FROM PYTHON SOURCE LINES 378-381
.. code-block:: Python
plt.show()
.. GENERATED FROM PYTHON SOURCE LINES 382-385
.. image:: ../../_static/brake_maximum_temperature.png
:align: center
:alt: Brake Maximum Temperature
.. GENERATED FROM PYTHON SOURCE LINES 387-388
Brake Maximum Temperature
.. GENERATED FROM PYTHON SOURCE LINES 390-392
Close the session
==================================================================================
.. GENERATED FROM PYTHON SOURCE LINES 392-395
.. code-block:: Python
session.exit()
.. _sphx_glr_download_examples_00-fluent_brake.py:
.. only:: html
.. container:: sphx-glr-footer sphx-glr-footer-example
.. container:: sphx-glr-download sphx-glr-download-jupyter
:download:`Download Jupyter notebook: brake.ipynb <brake.ipynb>`
.. container:: sphx-glr-download sphx-glr-download-python
:download:`Download Python source code: brake.py <brake.py>`
.. container:: sphx-glr-download sphx-glr-download-zip
:download:`Download zipped: brake.zip <brake.zip>`
.. only:: html
.. rst-class:: sphx-glr-signature
`Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_