CST Studio Suite - Thermal and Mechanical Simulation
CST Studio Suite - Thermal and Mechanical Simulation
Trademarks
3DS_Document_2019
CST, the CST logo, Cable Studio, CST BOARDCHECK, CST EM
STUDIO, CST EMC STUDIO, CST MICROWAVE STUDIO, CST
2024
PARTICLE STUDIO, CST Studio Suite, EM Studio, EMC Studio,
Microstripes, Microwave Studio, MPHYSICS, MWS, Particle Studio,
Suite
| CST Studioref.:
PCB Studio, PERFECT BOUNDARY APPROXIMATION (PBA),
SystèmesInformation|
Studio Suite, IdEM, Spark3D, Fest3D, 3DEXPERIENCE, the
Compass icon, the 3DS logo, CATIA, BIOVIA, GEOVIA,
| Confidential
SOLIDWORKS, 3DVIA, ENOVIA, NETVIBES, MEDIDATA, CENTRIC
PLM, 3DEXCITE, SIMULIA, DELMIA and IFWE are commercial
© Dassault
trademarks or registered trademarks of Dassault Systèmes, a French
© Dassault Systèmes
"société européenne" (Versailles Commercial Register # B 322 306
440), or its subsidiaries in the United States and/or other countries. All
other trademarks are owned by their respective owners. Use of any
Dassault Systèmes or its subsidiaries trademarks is subject to their
express written approval.
3DS.com/SIMULIA
2
Table of contents
Chapter 1 – Introduction ......................................................................................................... 5
Welcome ............................................................................................................................. 5
How to Get Started Quickly ............................................................................................. 5
What is CST MPHYSICS STUDIO? ................................................................................ 5
Who Uses CST MPHYSICS STUDIO? ............................................................................ 5
CST MPHYSICS STUDIO Key Features ............................................................................ 6
General ............................................................................................................................ 6
Structure Modeling ........................................................................................................... 6
Mechanics Solver ............................................................................................................ 6
Thermal Steady State Solver ........................................................................................... 7
Thermal Transient Solver ................................................................................................ 7
Conjugate Heat Transfer Solver ...................................................................................... 7
SAM (System and Assembly Modeling) .......................................................................... 8
Visualization and Secondary Result Calculation ............................................................. 8
Result Export ................................................................................................................... 8
3
Simulation Results ......................................................................................................... 37
Stationary Bi-Directional EM-CHT Link Set-up .............................................................. 39
Chapter 3 – Solver Overview ................................................................................................ 42
Solvers and Sources ......................................................................................................... 42
Mechanical Solver ............................................................................................................. 42
Thermal and Conjugate Heat Transfer Solvers ................................................................. 43
Background Material ...................................................................................................... 43
Material Properties ......................................................................................................... 44
Boundary Conditions ..................................................................................................... 47
Imports/Sources and Loads ........................................................................................... 50
Monitors at Points .......................................................................................................... 62
Monitors on Faces ......................................................................................................... 63
3D Field Monitors ........................................................................................................... 64
Steady State Thermal Solver Parameters ..................................................................... 64
Steady State and Transient Conjugate Heat Transfer Solver Parameters .................... 65
Excitation Signal Settings .............................................................................................. 66
Transient Thermal Solver Settings ................................................................................ 67
4
Chapter 1 – Introduction
Welcome
Welcome to CST Studio Suite®, the powerful and easy-to-use electromagnetic field
simulation software. This program combines a user-friendly interface with unsurpassed
simulation performance. CST Studio Suite contains a variety of solvers for carrying out
Thermal and Mechanical Simulation. They are all grouped as a specific Thermal and
Mechanical Module, also known as CST MPHYSICS® STUDIO.
Please refer to the CST Studio Suite - Getting Started manual first. The following
explanations assume that you have already installed the software and familiarized
yourself with the basic concepts of the user interface.
A key feature of CST MPHYSICS STUDIO is its tight integration with the other CST
Studio products. This allows an easy to use workflow for coupled EM-Multiphysics
simulations.
A further outstanding feature is the full parameterization of the structure modeler, which
enables the use of variables in the definition of your component. In combination with the
built-in optimizer and parameter sweep tools, CST MPHYSICS STUDIO is capable of
analyzing and designing thermal and mechanical aspects of devices.
5
CST MPHYSICS STUDIO Key Features
The following list gives you an overview of the CST MPHYSICS STUDIO main features.
Note that not all of these features may be available to you because of license restrictions.
Contact a sales office for more information.
General
Graphical user interface for Windows 10, Windows Server 2016/2019, Windows
11 and Windows Server 2022.
The structure can be viewed either as a 3D model or as a schematic. The latter
allows for easy coupling of thermal simulation parameters with circuit simulation.
Various independent types of solver strategies (based on hexahedral as well as
tetrahedral meshes) allow accurate simulations with a high level of performance
for a wide range of multi-physical applications.
For specific solvers highly advanced numerical techniques offer features like
Perfect Boundary Approximation® (PBA) for hexahedral grids and curved and
higher order elements for tetrahedral meshes.
Structure Modeling
Advanced ACIS-based, parametric solid modeling front end with excellent
structure visualization
Feature-based hybrid modeler allows quick structural changes
Import of 3D CAD data from ACIS® SAT/SAB, CATIA®, SOLIDWORKS®,
Autodesk Inventor, IGES, VDA-FS, STEP, PTC Creo, Siemens NX, Parasolid,
Mechanics Solver
Temperature dependent Young’s modulus
Displacement boundary condition
Traction boundary condition
Thermal expansion
Neo-Hookean material model for simulation of large deformations
Various stress plots: von Mises, hydrostatic and tensor components
Strain plots including visualization of the volumetric strain
Nonlinear solver computes the Green-Lagrange and the Almansi-strain as well
as the 2nd Piola-Kirchhoff and Cauchy stress tensors
Displacement plot including visualization of deformed mesh
Import of force densities from EM-solvers
Export of deformed structure to CST MICROWAVE STUDIO
6
Thermal Steady State Solver
Isotropic and anisotropic material properties
Bioheat material properties
Nonlinear material properties (Bioheat properties and thermal conductivity)
Thermal contact resistance
Moving media
Convection for human voxel models
Heat transfer by conduction in volumes
Heat transfer by convection and radiation through surfaces
Sources: fixed and floating temperatures, heat sources, eddy current and
stationary current loss fields, volume/surface power loss distributions in dielectric
or lossy metal materials imported from CST MICROWAVE STUDIO, CST EM
STUDIO or CST PCB Studio, crashed particle loss distribution from CST
PARTICLE STUDIO
Adiabatic / fixed or floating temperature / open boundary conditions
Automatic parameter studies using built-in parameter sweep tool
Automatic structure optimization for arbitrary goals using built-in optimizer
Network distributed computing for optimizations, parameter sweeps and remote
calculations
Thermal conductance matrix calculation
Equivalent Circuit EMS/MPS/DS Co-Simulation for linear problems
7
External heat sources imported from CST MICROWAVE STUDIO or CST EM
STUDIO
Axial/Centrifugal fan model support
Planar and volume flow resistance model support
Two-resistor/Delphi component model support
Thermal contact properties: resistance, capacitance
Thermal surface properties: surface emissivity and heat transfer coefficient
Heat pipes support
Thermoelectric coolers (TEC) based on the Peltier effect
ECXML file format import
Full GPU acceleration support
Import of 2D losses computed with the LF-TD planar solver
Temperature dependent fan curve for time variant simulations
Solar radiation: Solar radiation calculator based on ASHRAE Earth model.
Direct, diffused and reflected components supported.
Various 2D and 3D field visualization options for thermal fields, heat flow
densities, displacement fields, stress fields, etc.
Animation of field distributions
Result Export
Export of result data such as fields, curves, etc. as ASCII files
Export screen shots of result field plots
Automation
Powerful VBA (Visual Basic for Applications) compatible macro language with
editor and macro debugger
OLE automation for seamless integration into the Windows environment
(Microsoft Office®, MATLAB®, AutoCAD®, MathCAD®, Windows Scripting
Host, etc.)
8
About This Manual
This manual is primarily designed to enable a quick start of CST MPHYSICS STUDIO.
It is not intended to be a complete reference guide to all the available features but will
give you an overview of key concepts. Understanding these concepts will allow you to
learn how to use the software efficiently with the help of the online documentation.
Document Conventions
Buttons that should be pressed within dialog boxes are always written in italics,
e.g. OK.
Key combinations are always joined with a plus (+) sign. Ctrl+S means that you
should hold down the Ctrl key while pressing the S key.
The program’s features can be accessed through a Ribbon command bar at the
top of the main window. The commands are organized in a series of tabs within
the Ribbon. In this document a command is printed as follows: Tab name: Group
name Button name Command name. This means that you should activate
the proper tab first and then press the button Command name, which belongs to
the group Group name. If a keyboard shortcut exists, it is shown in brackets after
the command.
Example: View: Change View Reset View (Space)
The project data is accessible through the navigation tree on the left side of the
Your Feedback
We are constantly striving to improve the quality of our software documentation. If you
have any comments regarding the documentation, please send them to your support
center: 3DS.com/Support.
9
Chapter 2 – Simulation Workflows
This chapter contains two workflow examples demonstrating the basic features of CST
MPHYSICS STUDIO. In the first example, a very simple structural mechanics model of
an accelerometer is created. This workflow describes in detail, how to generate a model
geometry, assign material properties and sources, generate a mesh and run the
simulation. Besides, the visualization and interpretation of structural mechanics results
are discussed.
The second example describes the detailed workflow for setting up uni- and bi-
directional EM-Thermal coupled simulations. In the uni-directionally coupled simulation,
a high frequency electromagnetic solver calculates the ohmic losses in the walls of an
HF-filter which are imported by the Conjugate Heat Transfer (CHT) solver. Then the
CHT solver conducts a thermal analysis to get the temperature distributions in the filter
and its surroundings. In the bi-directionally coupled simulation, the EM solver and the
CHT solver exchange not only ohmic losses but also the temperature field, i.e., the high
frequency EM solver can import the temperature field calculated by the CHT solver to
update its temperature-dependent EM properties of its materials.
Studying these examples carefully will help to become familiar with many standard
operations that are important when performing simulations with CST MPHYSICS
STUDIO.
In the subsequent chapters some remarks concerning the extended features of the
The following explanations describe the “long” way to open a particular dialog box or to
launch a particular command. Whenever available, the corresponding Ribbon item will
be displayed next to the command description. Because of the limited space in this
manual, the shortest way to activate a particular command (i.e. by either pressing a
shortcut key or by activating the command from the context menu) is omitted. You should
regularly open the context menu to check available commands for the currently active
mode.
The Structure
The following picture demonstrates the spatial structure of a simple accelerometer. It
consists of two fixed flat conductors with a potential difference applied, and a movable
conductor between them.
10
1
3
2
If the system moves with acceleration, the inertial force pushes the movable conductor
towards one of the fixed ones. The potential difference, e.g., between the conductors 2
and 3 changes proportionally.
The following dialog box should now be positioned in the upper right corner of the main view:
11
The blue arrow always indicates the next step necessary for your problem definition. Usually,
a project template will already set the problem type and initialize some basic settings like
units and background properties. You do not need to process the steps in this order, but we
recommend that you follow this guide at the beginning in order to ensure all necessary steps
have been completed.
Look at the dialog box as you follow the various steps in this example. You may close the
assistant at any time. Even if you re-open the window later, it will always indicate the next
required step.
If you are unsure of how to access a certain operation, click on the corresponding line. The
Quick Start Guide will then either run an animation showing the location of the related menu
entry or open the corresponding help page.
1. Select the brick creation tool from the main menu: Modeling: Shapes Brick
2. Press the Escape key in order to open the dialog box.
3. Fill up the brick size fields as it is shown in the table below.
12
4. In order to select the material, click on the corresponding combo box and select
Copper (annealed). This material is predefined for CST MPHYSICS STUDIO
projects.
The dialog box Material Parameters: Copper (annealed) appears where various
properties of copper can be modified. Select the tab Mechanics in this dialog box.
13
© Dassault Systèmes | CST Studio Suite 2024
In this tab you can change the mechanical properties of the selected material. These
are the three most important mechanical properties:
Young's modulus defines the stiffness of the isotropic elastic material. It is normally
measured in GPa, or kN/mm2. The typical values vary between 0.01 GPa (rubber)
and over 1000 GPa (diamond). It is important to know the value of this material
parameter very well, since it has a large influence on the accuracy of the solution.
Poisson's ratio defines the scale of the transverse contraction of a longitudinally
stretched body. This parameter can vary between -1 and 0.5, whereas most of the
materials are characterized by a positive Poisson's ratio.
Thermal expansion coefficient is the strain of a body if its temperature changes by
1 K. This value is utilized to compute strain induced by an external temperature
field.
2. Now press Cancel and start creating a new brick (Modeling: Shapes Brick ) with
the following size (please do not press OK yet):
Xmin -6 Xmax 6
Ymin -1 Ymax 1
Zmin 0.05 Zmax 0.7
3. In order to change the material for the new solid, select [Load from Material Library…]
in the Material combo box. The dialog box Load from Material Library appears. Select
the material Steel-1010 and press the button Load.
14
© Dassault Systèmes | CST Studio Suite 2024
4. Now press the button OK in the Brick dialog box. A new brick consisting of Steel-
1010 is created.
5. By selecting Modeling: Picks Picks (S) activate the general pick tool to pick two
edges of the second brick, as shown in the picture below:
7. Again, open the Brick dialog and enter the following values:
15
For the new brick a new material should be created. Select [New Material…] in the
Material combo box. The New Material Parameters dialog is shown. In the General
tab, enter Plastic for the Material name.
16
After that, switch to the Mechanics tab in this dialog, select Normal for material
Type, set the Young’s modulus to 2 GPa and the Poisson’s ratio to 0.4.
8. Pick and chamfer one of the upper edges with the chamfer width of 0.7, and keep the
default angle of 45° (Modeling: Tools Blend Chamfer Edges ) in order to obtain
the following structure:
17
9. Create the following bricks:
One of Plastic with the following size:
The last one made of Copper (annealed) with the following dimensions:
18
11. In the dialog Transform Selected Object select the operation Mirror, check the
boxes Copy and Unite and set the mirror plane normal to 0, 0, 1, as shown in the
following picture. If necessary, uncheck Shape Center for the Mirror plane normal:
In the present example let us fix the both sides of the model and apply a pressure to the
middle electrode, which would mimic the influence of inertial forces during acceleration.
The following steps must be performed:
19
1. Press the toolbar button Simulation: Boundaries Displacement Boundary .
2. Select the side faces of the model, as shown in the picture below (you have to select
five faces at x-min and 3 faces at x-max):
After pressing the Return key, the dialog box Define Displacement Boundary will appear:
After pressing the Return key, the dialog box Define Traction Boundary will appear:
20
6. Put the value of -2e-6 GPa as the Z-coordinate of the traction vector. This means that
the pressure of 2 kPa is applied towards the negative direction of the Z-axis. This
pressure would roughly correspond to the acceleration of 17*g, or 170 m/s 2, into the
positive Z-direction.
Mesh Settings
The structural mechanics solver is quite sensitive to the quality of discretization. In order
to obtain reliable results, the default mesh density needs to be increased. To do this,
press the toolbar button Simulation: Mesh Global Properties .
In the Mesh Properties – Tetrahedral dialog, change the Cells per max model box edge
setting for Model to 20:
21
Press the OK button to accept the changes and close the window.
You can click the Help button in order to learn more about the controls in this dialog box.
For now, the default settings are good enough, so just click the Start button. After the
calculation has been started, you can control the execution of the solver in the Progress
and Messages windows.
22
The directory NT: 1D Results Adaptive Meshing contains information on the adaptive
mesh refinement performed by the solver. Here you can inspect the number of cells in
the mesh for each iteration step, time used by the solver to generate the solution, as
well as the relative error of the solution. For example, in the picture below you can see
the number of degrees of freedom in the solution for each step of mesh refinement.
Please note that the exact values may differ slightly on different systems.
A click on the item NT: 2D/3D Results Displacement displays a deformation plot of
the body deformation.
23
Here the original shape of the model is shown semi-transparently whereas the scalar
plot of absolute displacement is shown on the solid deformed shape.
Select Arrows from the plot type pull-down menu in 2D/3D Plot: Plot Properties to display
a vector plot of the body deformation, as shown in the following picture.
Selecting Contour from the plot type pull-down menu displays a scalar plot and enables
the vector component pull-down menu in the 2D/3D Plot ribbon, which contains the
following items: X, Y, Z, Abs, Normal, Tangential.
The item Abs demonstrates the distribution of the absolute value of displacement
within the solution domain.
The items Normal and Tangential demonstrate the length of the corresponding
projection of displacement vector onto each body surface.
Navigation Tree item NT: 2D/3D Results Strain contains the following sub-items:
Directory Components contains the components XX, YY, ZZ, XY, XZ and YZ of the
strain tensor.
Sub-item Volumetric displays the distribution of the volumetric strain in the model,
which means the relative volume change in each node of the solution domain. The
negative values mean contraction, whereas the positive values mean expansion.
24
Finally, Navigation Tree item NT: 2D/3D Results Stress contains the following entries:
Directory Components contains the components XX, YY, ZZ, XY, XZ and YZ of the
stress tensor.
The tree-entry Von Mises displays the distribution of von Mises stress within the
solution domain. If this stress at some location is higher than the yield strength of the
corresponding material, plastic deformation takes place in this location. Von Mises
stress is always positive.
The tree-entry Hydrostatic displays the hydrostatic stress distribution, reproducing
the change of the volume in the stressed body. The negative values mean contraction
forces.
The tree-entry First Principal Stress displays the distribution of the largest eigenvalue
In the following picture the distribution of the absolute value of displacement vector is
shown on the cutting plane perpendicular to the Y-axis.
25
Vector fields can be visualized on a cutting plane in the same manner. Just select the
Arrow plot type in the 2D/3D Plot ribbon (of course, for a plot with vector data like
Displacement). In this case the Fields on Plane mode stays activated.
Summary
This example should have given you an overview of the key concepts of CST
MPHYSICS STUDIO. Now you should have a basic idea of how to do the following:
If you are familiar with all these topics, you have a very good starting point for further
improving your usage of CST MPHYSICS STUDIO.
For more information on a particular topic, we recommend you browse through the online
help system which can be opened via the Help button in the upper right corner. If you
have any further questions or remarks, do not hesitate to contact your technical support
team. We also strongly recommend that you participate in one of our special training
classes held regularly at a location near you. Ask your support center for details.
Two types of EM-CHT couplings are supported, i.e., uni-directional and stationary bi-
directional couplings. In a uni-directional EM-CHT coupling the EM solver first solves the
electromagnetic fields and the resultant thermal losses. The CHT solver then imports
those losses as heat sources and performs a thermal analysis to obtain the temperature
field in the computational domain. Currently, the EM solvers which can be uni-
directionally coupled with the CHT solver include the HF frequency domain solver, HF
transient solver, LF frequency domain solvers (Fullwave, Eletroquasistatic and
26
Magnetoquasistatic), LF time domain solvers (Eletroquasistatic and
Magnetoquasistatic), and Stationary Current solver.
A stationary bi-directional EM-CHT coupling not only allows the CHT solver to import the
thermal losses obtained from an EM solver but also allows the EM solver to import the
temperature field calculated by the CHT solver. The bi-directional coupling can find its
application in problems involving materials whose EM properties are dependent on
temperature. Currently, the EM solvers which can be coupled with the CHT solver in
stationary bi-directions include the HF frequency domain solver, HF transient solver, LF
frequency domain solver (Magnetoquasistatic), and Stationary Current solver.
The typical workflow for setting up a uni-directional EM-CHT coupled project is first
demonstrated. The simulated device consists of a filter placed on a horizontal support
and surrounded by air. The HF frequency domain solver is first used to perform a
frequency domain analysis of the filter. The ohmic losses from the filter will be obtained
along with the electromagnetic fields. Physically, the ohmic losses will be transformed
into heat and result in the increase of the filter’s temperature. The corresponding thermal
analysis is conducted by using the CHT solver which imports the ohmic losses as heat
sources and computes the temperature distribution in the filter with the cooling effect of
the surrounding air being considered.
The workflow for setting up a stationary bi-directional EM-CHT coupled project shares a
lot with its counterpart for the uni-directional coupling and, therefore, only the difference
To open the project you have to download a copy first, by clicking on the Download
symbol. Once this is done you are ready to open the example by clicking on the Open
Project symbol.
27
Start the automatic creation of a coupled electromagnetic/thermal computation by
selecting Home: Simulation Simulation Project EM-Thermal Coupling Uni-
directional.
A dialog box appears to create the first part of the EM-Thermal link:
After you click the OK button, the dialog box below appears to create the second part of
the EM-Thermal link:
28
In the dialog box, the project type Thermal & Mechanics is already chosen. So you only
need to select the thermal solver type Conjugate Heat Transfer and rename the project
to CHT. Once you are done, press the OK button. The EM and thermal simulation tasks
whose names are EM1 and CHT, respectively, will be created and added to the Tasks
folder in the Navigation Tree in the Schematic view sequentially.
Click OK. The corresponding monitors and field imports are configured automatically in
both simulation projects. Now you may switch to the thermal project (select the CHT
project tab in the main view) and configure additional material properties, necessary
thermal sources, boundary conditions and calculation parameters.
In the navigation tree a field source called EM1 (named after the name of the EM
project of the EM-Thermal link) has been automatically added and configured. Edit it to
29
reconfigure it if necessary. The exclamation mark indicates that losses are missing
because the EM1 simulation has not yet been performed.
Select the Thermal tab and check the properties before closing the dialog with OK.
30
© Dassault Systèmes | CST Studio Suite 2024
Note: solids made of material PEC will be automatically replaced in the solver by the
material Copper (annealed).
Make sure that the material named “MWSSCHEM1/dielectric” has the following thermal
properties (the density value can also be entered in the Density tab):
31
© Dassault Systèmes | CST Studio Suite 2024
Please check the material named “MWSSCHEM1/brass” in the same manner:
32
© Dassault Systèmes | CST Studio Suite 2024
Open Simulation: Settings Boundaries to exploit the symmetry of the model with
respect to the YZ plane and set-up the YZ symmetry boundary condition:
33
Mesh Settings
Please change the mesh type to CFD (Simulation: Mesh Global Properties CFD)
and open the mesh properties dialog. Adjust the Minimum cell setting to a fraction of 40
for the maximum cell in the model.
Solver Parameters
Please change the default temperature unit to Celsius in the project Units dialog (Home:
Settings Units ) to have the results presented in Celsius later.
34
Then open the solver parameter dialog (Simulation: Solver Setup Solver ). Activate
the Gravity check box to simulate the effect of natural convection. Switch off the
Turbulence model check box. Adjust the Ambient temperature unit to Celcius and
specify the ambient conditions of 20°C.
Open the CHT Solver Special Settings dialog by pressing the button Specials… and limit
the Number of iterations to a maximum of 80 in order to shorten the otherwise lengthy
simulation time. Please note that the results obtained after 80 iterations are not fully
converged. If more accurate results are desired change the number of iterations to
automatic calculation (i.e. switch off the Maximum checkbox). With current settings, the
overall simulation time should be less than 20 minutes.
35
© Dassault Systèmes | CST Studio Suite 2024
Apply all settings with the OK or Apply buttons before closing the CHT Solver Special
Settings and Conjugate Heat Transfer Solver Parameters dialog box.
Coupled Run
Switch back to schematic of the master project (first tab) and therein to the Schematic
view (select the appropriate tab at the bottom of the main view). Press the button Home:
Simulation Update . At first, the EM calculation will be started. Next, the losses will
be computed. Finally, these losses will be imported into the thermal project, and the
thermal calculation will be performed.
Alternatively, right click on NT: Tasks Coupled EM-Thermal1 and Update the task.
A progress bar will appear in the progress window which will update you on the task
progress. You can activate this window by selecting View: Window Windows
Progress Window. Information text regarding the simulation will appear above the
progress bar. The most important stages are listed below for the CHT solver:
1. Updating tasks: 1of 1: the selected task includes the previously created EM and
CHT simulation.
2. …. the EM simulation is performed….
36
3. CHT solver: Surface mesh generation: the solid surfaces are triangulated.
4. CHT solver: Octree grid generation: the CFD mesh is constructed by using the
solid surface triangulations.
5. CHT solver: Importing surface/volume losses: the losses from the EM simulation
are imported and mapped into the CFD mesh.
6. CHT solver: Upgrade grid: inactive cells are removed from the CFD mesh.
7. CHT solver: Iterations: the simulation is performed.
Simulation Results
Once the EM simulation has been completed please leave the schematic and return to
the CHT simulation. Follow the progresses of the CHT simulation by looking at the
convergence monitors in the NT: 1D Results Convergence monitors Equation
residuals and NT: 1D Results Convergence monitors Equation balances.
The simulation completes 80 iterations before stopping. The following dialog box pops
up because the simulation has not fully converged (i.e. several convergence criteria
have not been met due to the low maximum number of iterations):
37
To visualize the loss imported from the EM simulation, select NT: 2D/3D Results Heat
source densities and a cut plane, for instance X=0. Please note that the losses can only
be visualized on cut planes (check ribbon 2D/3D Plot: Sectional View Fields on
Plane). Observe that the losses are the highest on the walls of the coaxial feeds and of
the cylinders.
Once the simulation has stopped, visualize in the same cut plane the temperature by
selecting NT: 2D/3D Results Temperature.
The temperature increases are the highest where the losses are also the highest. The
air in contact with the walls of the filter heats up and carries the heat away which cools
down the filter. One can observe that the simulation has not converged to a steady-state
solution in the whole domain because the heat carried by the air flow has not yet reached
the top boundary of the computational domain. Still, the simulation has enough
progressed to show a correct temperature distribution inside the filter.
The CHT solver takes into account the air cooling effect by simultaneously calculating
the heat transfer in the fluid and solid domains and the air flow caused by the
temperature gradients and gravity. This key feature differentiates the CHT solver from
the thermal solvers which do not solve for the air flow and thus can simulate neither
natural nor forced convection.
38
The air flow can be visualized by selecting NT: 2D/3D Results Velocity
The velocity vector plot shows the air circulating inside the filter as well as the heated
air ascending and being replaced by air at ambient temperature.
The dialog boxes for creating the EM and CHT simulation projects are the same as those
in setting up the uni-directional coupling:
39
© Dassault Systèmes | CST Studio Suite 2024
However, after both the EM and thermal simulation projects are created and before the
loss import dialog appears, a dialog box which is unique to the bi-directional coupling
will pop up. It asks for the number of iterations between the EM and the thermal tasks:
After you input your desired number and press OK, the above dialog box will disappear
and a loss import dialog will pop up. From now on, the set-up process is identical to the
previous workflow for uni-directional couplings.
Note that the bi-directional EM-Thermal coupling is generally used for the cases where
the EM properties of material are temperature dependent. Project “Combline Filter
Draft” is adopted here just for demonstration of the workflow. The properties of its
materials are not dependent on temperature. Special attention also needs to be paid to
40
the bi-directional coupling between the HF frequency domain solver and the thermal
solvers. After both the HF frequency domain solver and the CHT sub-projects are
created, “Rebuild simulation projects” in the property of “Sequence” in Schematic view
should be checked to ensure all frequency samples are re-calculated by the HF
frequency domain solver in every iteration.
41
Chapter 3 – Solver Overview
Solvers and Sources
Various simulation types differ in the definition of materials, boundary conditions and
sources. The way to define materials in CST MPHYSICS STUDIO is quite similar for all
solvers, whereas there are larger differences in the definition of sources and boundary
conditions. For this reason, an overview of the sources, loads and boundaries for each
solver are explained below.
Mechanical Solver:
Displacement boundary: Simulation: Boundaries Displacement Boundary
Traction boundary: Simulation: Boundaries Traction Boundary
External temperature and/or force distribution:
Simulation: Sources Field Import
Mechanical Solver
The mechanical solver is a tetrahedral based solver for structural mechanic problems.
Its main application is computing deformations driven by thermal expansion and external
forces. The deformation results can be used for a subsequent High Frequency
Electromagnetic analysis with the tetrahedral based frequency domain solvers from CST
MICROWAVE STUDIO.
Refer to the chapter Simulation Workflow for a description of the basic features. The
import of temperature and force density distributions is described in the section Workflow
for Coupled Simulations.
42
Thermal and Conjugate Heat Transfer Solvers
CST MPHYSICS STUDIO includes a thermal and a conjugate heat transfer (CHT)
solver. The thermal solver is optimized to simulate thermal conduction in the steady state
and transient regime and supports hexahedral and tetrahedral grids. The CHT solver is
a CFD based heat transfer solver capable of solving thermal conduction, convection and
radiation simultaneously in the steady state and transient regime. The effect of solar
radiation is also supported. The main applications of these solvers include solving
steady state or transient temperature problems resulting from various types of losses.
Both solvers are also well suited to compute standalone thermal problems. The following
sub-sections will demonstrate the most important aspects of a thermal simulation with
CST MPHYSICS STUDIO.
Background Material
The first step for setting up a thermal simulation is to define the units for temperature
and dimension, like it has been described in the chapter Simulation Workflow.
Afterwards an appropriate background material should be selected. Open the material
background properties dialog box by selecting Modeling: Materials Background :
43
© Dassault Systèmes | CST Studio Suite 2024
The easiest way to assign the necessary values is to copy the properties from an existing
material in the material library. Press the Copy Properties from Material… button in the
General tab, select [Load from Material Library…] in the Copy Properties from Material
dialog box:
Material Properties
The material parameters for a thermal problem can be defined inside the material
parameters dialog box: Modeling: Materials New/Edit New Material . Select the
Thermal tab.
44
If you select a PTC (Perfect Thermal Conductor) type, an infinite thermal conductivity is
assumed. A body with PTC material assigned always has a uniform temperature.
Please note that the conjugate heat transfer solver replaces PTC with copper.
For transient thermal problems (see also the next chapter) and the conjugate heat
transfer solver the heat capacity and the material density must be specified. These
parameters determine how much energy per Kelvin is stored in a certain amount of mass
or volume:
The field Thermal expansion coefficient is available only if the CHT solver type is
selected. It allows to enter the volumetric coefficient of thermal expansion in [1e-6 / k].
The thermal expansion coefficient describes how the volume of a fluid changes with a
change in temperature and must be set when natural convection is simulated i.e. when
the gravity is defined and the fluid flow is simulated. If the value is zero the solver applies
the ideal gas law to calculate the thermal expansion coefficient as the inverse of the
ambient temperature. This law is not valid for liquids therefore it is important to set the
thermal expansion coefficient when a liquid is simulated and gravity is defined.
Because the thermal diffusivity plays an important role for the transient simulation
process, it is shown here as well. The diffusivity can be calculated from the thermal
conductivity, the heat capacity and the material density as follows:
45
k
,
cP
where
: Diffusivity [m² / s]
k: Thermal conductivity [W / K /m]
: Density [kg / m³]
cP: Specific heat capacity [J / K / kg]
Nonlinear heat capacity can be used for simulation of material phase change in transient
computations. This can be achieved by a local increase of heat capacity for a small
interval of temperatures. For more information on simulation of phase changes, please
contact support service.
For simulations which involve biological materials, heating mechanisms of living tissue
can be taken into account (see also: Bioheat Source below). In addition, it is possible to
define a convection coefficient for surface materials of human voxel models (typically:
skin).
The Flow Resistance material parameter is only supported by the conjugate heat
transfer solver. It is used to model the fluid flow behavior across a screen without having
to mesh the screen geometry.
A flow going through a volume resistance experiences a pressure gradient which can
be written as follows:
Pj 0.5 f ij u j u j
j ,
where f ij is the loss coefficient tensor per unit length and u j is a velocity component in
the global coordinate system (X,Y,Z).
The loss coefficient tensor is defined with respect to a local coordinate system (U’,V’,W’)
and transformed into a 3x3 tensor in the global coordinate system.
A Flow Resistance assigned to a surface uses the specifications of the sheet properties
group whereas a Flow Resistance assigned to a solid uses the specifications of the solid
properties.
46
© Dassault Systèmes | CST Studio Suite 2024
Boundary Conditions
The boundary conditions for the thermal and conjugate heat transfer solver can be
defined in the Thermal Boundaries tab of the Boundary Conditions dialog box
(Simulation: Settings Boundaries )
For “isothermal” and “open” boundaries the temperature settings may be assigned by
pressing the corresponding button […]. This button opens the dialog Boundary Settings,
47
in which the temperature value can further be configured, for example, by assigning of
a fixed or floating temperature. By default, the option Unset is selected, which means
the boundary is considered as a PTC surface without sources assigned.
For the "open" boundary condition, it is assumed that the temperature approaches the
predefined value with increasing distance from the structure. Apply this type of boundary
condition if thermal conduction through the surrounding background material plays an
important role for your problem. In order to consider thermal convection effects on the
structure, Thermal surface properties (see p. 53) should be used.
When no heat flow leaves the computational domain through a boundary, use the
"adiabatic" boundary condition. In case the conductive heat flow of an open structure
can be neglected, you can use these boundary conditions instead of “open” boundary
conditions (if radiation or convection effects dominate).
The following table shows an overview, where T is the temperature and Q is the heat
flux density:
The picture below illustrates an example of how thermal fields are influenced by the
different boundary types. It shows a metal sphere at a constant temperature, which is
surrounded by a material with constant thermal conductivity.
48
For Conjugate Heat Transfer Solver:
The conjugate heat transfer solver supports similar types of boundary conditions. It is
however important to note that it interprets these boundary conditions differently, in
particular for the case of open boundaries.
User input for adiabatic walls: friction at the boundary wall (no heat exchange, zero
emissivity).
User input for symmetrical boundaries: none (no friction, no heat exchange, zero
emissivity).
User input for open boundaries: flow temperature, flow velocity or flow gauge
pressure. An open boundary allows flow to enter and leave the domain, which could be
used to model the flow and thermal behavior of an inlet or of an outlet specified by a
pressure gauge or a velocity. If the flow temperature is unknown, which is the case for
outlets or if the flow direction is unknown set the temperature to unset.
The emissivity is set to 1. The reference temperature used for radiation is the radiation
temperature defined in solver parameter dialog.
49
Imports/Sources and Loads
The thermal and conjugate heat transfer solvers can handle several types of sources or
loss mechanisms, which are listed below:
Temperature Source
This source is available via Simulation: Sources and Loads Temperature Source .
This source type can be assigned to a surface of an object with PTC material properties
or any other material with non-zero thermal conductivity. For the transient thermal solver
an initial temperature source can be defined, which is taken into account only for
generation of the initial temperature distribution and ignored during the transient
solution.
When assigned to a solid that is either PTC or PEC it defines the total heat flow
coming from the solid surface. Therefore, a heat source with zero heat flow and a
floating temperature are identical.
50
Thermal Loss Distribution
This source is available via Simulation: Imports Thermal Losses . Thermal losses
can occur inside materials with finite conductivity, on surfaces of good conductors, inside
dispersive materials or at materials where particles hit the surface. These loss
distributions can be imported and used as thermal sources inside thermally conductive
materials. If previously calculated loss distributions are present, you can edit setting by
reopening the dialog box (Simulation: Imports Thermal Losses ).
51
Crashed particles Tracking Solver ( ), PIC Solver ( )
whereas QRadiation stands for the radiated power, T for the surface temperature, TReference
for the reference temperature, which can be equal to ambient or user-defined, for the
Stefan-Boltzmann constant and ASurface for the area for the selected surfaces. An
emissivity value = 0 means that the surface does not lose thermal power by radiation.
A value of 1 means that the thermal power emitted by the surface equals to that of a
black body at the same temperature.
52
where QConvection denotes the power, T the solid surface temperature, TReference the
reference temperature in the fluid and ASurface the area for the selected surfaces.
The thermal surface properties dialog includes additional options for the conjugate heat
transfer solver. The emissivity of the solid defined by the emissivity of its material can
be overwritten by a surface emissivity for the assigned surface. Similarly, the solar
absorptance of the solid defined by its material can be overwritten here for the respective
surface. In addition, the local fluid temperature can be used as the reference
temperature when convection is prescribed by a heat transfer coefficient.
53
The fan characteristics (i.e. fan curve, volume flow rate or stagnation pressure) are given
54
Heat Pipe (not supported by Thermal solver)
Heat Pipes are available via Simulation: Sources and Loads Heat Pipe .
A heat pipe is defined in two steps: In the first step, the shape(s) forming the heat pipe
must be selected. In the second step, the shapes in contact with the heat pipe are
selected. The shapes in contact with the heat pipe generally include a heat source and
a heat sink. After the second step, the properties of the heat pipe can be entered.
55
The junction node of the Delphi model is an internal node which doesn't interact with the
package environment. The top inner, top outer, bottom inner, bottom outer and side
nodes of the Delphi model are called surface or external nodes because they are
associated to certain physical regions or patches located on the package surface,
allowing them to interact with the package environment as shown below:
The top inner node is associated to a rectangular and centered patch on the upper
(package) surface. The top outer node is associated to the patch delimited by the edges
of the top inner node and the edges of the upper surface.
The bottom inner node is associated to a rectangular and centered patch on the bottom
surface. The bottom outer node is associated to the patch delimited by the edges of the
bottom inner node and the edges of the bottom surface.
56
The junction is an internal node. Other internal nodes can be added if necessary by
adding additional nodes to the Node / Resistance table.
Note that the thermal conductivity and specific heat of the shape used to model the
package geometry have no influence on the behavior of the Delphi model.
You can find more detail about Delphi in the online documentation.
When the check box Invert flow direction or the Switch fluid side are visible in the Edit
Lid or Edit Opening dialog, the normal flow direction at the selected face may be
previewed in the 3D main view. Check or uncheck the box to change the arrow direction
57
so that it shows the flow direction (Invert flow direction) or that it points toward the fluid
(Switch fluid side).
A lid specifies the flow properties of a fluid inside a fluid domain. The following boundary
types are available:
58
Open Set the flow temperature. Set the flow velocity,
volume flow rate or gauge pressure.
The opening emissivity is set to one.
The radiation temperature defined in the CHT Solver
Parameters dialog is used as the reference
temperature.
The type opening is used to prescribe flow properties within the same fluid domain.
Therefore, heat and/or mass can be transferred across an opening.
Open Set the flow temperature. Set the flow velocity, volume
flow rate or gauge pressure.
The opening emissivity is set to one.
The radiation temperature defined in the CHT Solver
Parameters dialog is used as the reference
temperature.
Depending if the current temperature value T is higher or lower than the blood
temperature this mechanism cools or heats the surrounding material. The blood
temperature value can be edited inside the Specials dialog box of the thermal solvers
(Simulation: Solver Setup Solver Specials or Simulation: Solver Setup Solver
Specials).
The Basal metabolic rate describes the amount of heat QMetabolic which is produced by
tissue per volume V.
QMetabolic V CMetabolic
59
area [K∙m2/W] as well as thermal capacitance [J/K]), or by its thickness and the thermal
properties of material assigned. Both definitions are equivalent and can be easily
converted into each other:
Here Rθ represents the absolute thermal resistance, rθ the thermal resistance per unit
area, C the thermal capacitance of the contact layer. In the material-based
representation, thermal conductivity k, specific heat capacity cP, material density ρ and
layer thickness l are used. The contact area A is calculated by the solver.
The advantage of contact properties definition through lumped parameters is the ease
and transparency of the parameter values. Besides, the absolute thermal resistance is
independent from the contact area A which may vary in case of solid intersections or
depending on the mesh settings. On the other hand, the material-based definition offers
more flexibility, for example it supports nonlinear material properties.
60
To set a user-defined initial temperature to a shape please create a "New Initial
Temperature On Solid..." and drag and drop the shape on it (or drag the initial condition
on the shape).
To set a user-defined initial temperature and velocity to a fluid domain please create
a "New Initial Condition on Fluid Domain..." and drag and drop the fluid domain on it (or
drag the initial condition on the fluid domain).
This vector defines the velocity with which the material comprising the solid is moving
relatively to the sources and solid geometry. A typical example would be a very long
tube moving through a coil for the purpose of induction heating.
If a velocity vector has been assigned to any solid, the solver saves important
information about the distribution and maximum of Peclet number in order to control the
solution quality.
61
Only tetrahedral-based thermal solvers support this feature. In the transient solution, the
moving media velocity vector may be made time-dependent by assigning Excitation
Signals to its components.
You can find more detail about moving media in the online documentation.
You can find more detail about ECXML import in the online documentation.
Monitors at Points
The monitors of this kind record scalar values that are defined at a point (e.g. the x-
component of the heat current density at a fixed position). You can create these monitors
via Simulation: Monitors Monitor at Point .
Transient thermal solver records the temperature values at the monitor points during the
whole solution time interval.
Two additional types of monitor at point are available for the conjugate heat transfer
solver. The Pressure and Velocity types evaluate, respectively, the pressure and velocity
at the monitor point at each iteration.
The conjugate heat transfer solver saves the values of the monitor points as 1D data
into the Navigation Tree under NT: 1D Results Monitors at Points <monitor name>.
The conjugate heat transfer solver can use the point monitors activated in Simulation:
Setup solver: Accuracy: Custom stop criteria to detect the convergence of the solver.
62
This monitor type is similar, although not identical, to Probes available within CST
MICROWAVE STUDIO.
Monitors on Faces
The monitors of this kind record scalar values defined on a surface. You can create
these monitors via Simulation: Monitors Monitor on Faces .
Two types of monitors on face are available. The type flow flux is used to monitor the
fluid flow, consequently the monitor surfaces must not change the flow and must be
borrowed from a dummy solid. A dummy solid is either a solid whose material is exactly
the same as the background material or a solid not considered for simulation but
considered for the bounding box. If necessary, adjust the local mesh properties of the
dummy solid to match those of the background to avoid unwanted mesh refinements
around the dummy solid.
𝑚̇ = ∯ 𝜌𝐮 ∙ 𝐝𝐀
The type solid flux is used to monitor the heat flux and the heat transfer coefficient at
solid/fluid interfaces, consequently the monitor surfaces must be borrowed from a solid
considered for simulation and considered for the bounding box and whose material is
different from the background material:
𝑃 = ∯ −𝑘 ∙ ∇𝑇 ∙ 𝐝𝐀
The monitors on faces are evaluated at each iteration and the surface quantities are
saved as 1D data into the Navigation Tree under NT: 1D Results Monitors on Faces
<monitor name>.
63
The conjugate heat transfer solver can use the face monitors activated in Simulation:
Setup solver: Accuracy: Custom stop criteria to detect the convergence of the solver.
3D Field Monitors
In contrast to steady state solvers, field distributions delivered by transient solvers need
to be requested by the user in advance by defining Field Monitors by right clicking on
NT: Monitors 3D Field Monitors New Field Monitor . A dialog box opens where
the type of the field, the start time and the sample step width can be defined:
Before starting the solver, it is advisable to look at the Ambient temperature, which is by
default the reference temperature for the radiation and convection models as well as for
the open boundary conditions. Moreover, this temperature may be assigned to PTC
regions without user-defined temperature or heat sources.
If Bioheat properties must be adjusted, one can open the Specials dialog:
64
© Dassault Systèmes | CST Studio Suite 2024
This also applies to the transient thermal solver. For further details, please refer to the
online help.
The conjugate heat transfer solver allows the simulation of stationary or transient heat
transfers between solids immersed in the background or in fluid domains. Radiation,
convection and conduction are simulated according to the settings set in
the Physics section below.
To simulate conduction and convection, please make sure that the density, dynamic
viscosity, thermal specific heat capacity and thermal conductivity are larger than zero.
To simulate only heat transfer (including conduction in the fluids) but without advection
in the fluids, please make sure that the thermal conductivity is larger than zero. This
corresponds to the Fluid flow option turned off.
The simulation takes place in vacuum whenever all thermal parameters of the fluid
material are not defined, i.e. the thermal conductivity, specific heat, material density,
dynamic viscosity and thermal expansion coefficient are set to zero.
To switch between the steady state and the transient conjugate heat transfer solver it is
sufficient to change the simulation regime. By default the simulation regime is steady-
state.
Two different time integration methods are available in the transient regime:
Fixed The time step width is constant during the entire simulation. This type of
time integration will be robust but more computationally expensive than an
adaptive time integration.
Adaptive The time step width is adapted during the simulation process according to
the intermediate results and transient profiles. This type of time integration
should be preferred for most cases.
65
Excitation Signal Settings
For some transient thermal simulations, it is necessary to define time domain excitation
The parameters of the signal depend on the individual signal type and are described in
the online help. The parameter Ttotal must be set for almost all signal types and defines
the size of the definition interval. For time values larger than Ttotal the signal is, in
general, continued by a constant value. It is also possible to import a signal or to create
a user defined signal or to select a pre-defined signal from the signal database.
All defined signals are visible in the Signal folder in the Navigation Tree and can be
displayed by selection in the Navigation Tree:
66
Transient Thermal Solver Settings
You can switch between the steady state and transient thermal solvers by selecting
either
After selecting the transient solver, the solver parameters dialog box can be opened by
clicking on the icon in the Home or the Simulation ribbon (Simulation: Solver Setup
Solver ). Before starting the transient thermal solver, a valid Simulation duration time
must be entered:
Most source types can be weighted with a previously defined excitation function, when
pressing the Excitations button:
67
For each source, a signal can be assigned via a drop down list. The same signal can be
assigned to several sources. Optionally, an individual time delay t can be defined for
each source. The resulting time dependent excitation f is the product of the source
value v (e.g. the temperature) and the (possibly shifted) assigned signal s :
f (t ) s(t t ) v .
The initial temperature distribution can be defined in the Select Start Temperature
dialog, which can be called by pressing the Start temperature: Settings button. The
default setting is to assign the ambient temperature everywhere except regions with
temperature sources. Alternatively, it is possible to assign the solution of the steady-
state problem with initial source values as well as import a temperature distribution from
an external thermal solution.
Result Types
After a steady state thermal simulation run has been completed successfully, new
result entries appear in the navigation tree:
68
The directory 1D Results contains the convergence curve, heat flow values for the heat
sources as well as power scaling values for imported fields.
In the directory 2D/3D Results, beside the scalar temperature field the heat flow density
can be seen, which is a vector field showing the heat flow inside thermally conductive
materials. Moreover, a text file is written where the total heat flow for every source is
listed. In case field losses were imported, further information like interpolated loss
distributions as well as the scaling factor is presented.
The transient thermal solver creates a different output in the navigation tree:
If time domain temperature monitors have been defined for the transient thermal solver,
the associated results will be listed under 2D/3D Results as well. In addition, a couple
of time signals are added to the 1D Results section:
These 1D signals can be updated during the simulation process by selecting the tree
item and pressing 1D Plot: Plot Properties Update Results or the F5 key.
69
The 1D Results contain solution convergence, point and face monitors as well as
performance data, plotted against iterations to give user insights into convergence and
solutions.
70
Chapter 4 – Finding Further Information
After carefully reading this manual, you will already have some idea of how to use CST
MPHYSICS STUDIO efficiently for your own problems. However, when you are creating
your own first models, some questions may arise. In this chapter, we give you a short
overview of the available additional documentation.
The QuickStart Guide is opened automatically on each project start if the checkbox File:
Options Preferences Open QuickStart Guide is checked. Alternatively, you may
start this assistant at any time by selecting QuickStart Guide from the Help button in
the upper right corner.
When the QuickStart Guide is launched, a dialog box opens showing a list of tasks,
where each item represents a step in the model definition and simulation process.
Usually, a project template will already set the problem type and initialize some basic
settings like units and background properties. Otherwise, the QuickStart Guide will first
open a dialog box in which you can specify the type of calculation you wish to analyze
and proceed with the Next button:
In order to access information about the QuickStart Guide itself, click the Help button.
To obtain more information about a particular operation, click on the appropriate item in
the QuickStart Guide.
Online Documentation
The online help system is the primary source of information. You can access the help
system’s overview page at any time by choosing File: Help Help Contents . The
online help system includes a powerful full text search engine.
In each of the dialog boxes, there is a specific Help button, which opens the
corresponding manual page. Additionally, the F1 key gives some context sensitive help
when a particular mode is active. For instance, by pressing the F1 key while a basic
shape generation mode is active, you can get information about the definition of shapes
and possible actions.
71
When no specific information is available, pressing the F1 key will open an overview
page from which you may navigate through the help system.
Please refer to the CST Studio Suite Getting Started manual to find more detailed
explanations about the usage of the CST MPHYSICS STUDIO Online Documentation.
Technical Support
Before contacting Technical Support, you should check the online help system. If this
does not help to solve your problem, you find additional information in the Knowledge
Base and obtain general product support at 3DS.com/Support.
History of Changes
An overview of important changes in the latest version of the software can be obtained
by following the What’s New in this Version link on the help system’s main page or
from the File: Help backstage page. Since there are many new features in each new
version, you should browse through these lists even if you are already familiar with one
of the previous releases.
72