Difference between revisions of "Thermal Calculations"

From OSUPDOCS
Jump to navigation Jump to search
 
(30 intermediate revisions by 2 users not shown)
Line 21: Line 21:
== Conduction ==
== Conduction ==


The <tt>Conduction</tt> command activates thermal conductivity calculations and also activates various optionals involving heat flow in the simulations or other mechanisms that cause temperature changes. These settings can change and select the [[#Thermodynamics Modes in MPM|thermodynamics mode]] of the simulation.
The <tt>Conduction</tt> command activates thermal conductivity calculations and also activates various optionals involving heat flow in the simulations or other mechanisms that cause temperature changes. These settings can change and select the [[#Thermodynamics Modes in MPM|thermodynamics mode]] of the simulation. By default, conduction uses update methods analogous to FLIP methods used in mechanics. This update, however, sometimes results in temperature oscillations on particles within one cell. Conduction simulations with oscillations can be improved by using [[PeriodicXPIC Custom Task#Using FMPM(k) For Transport Properties|periodic FMPM(k)]] for conduction updates.


=== Scripted Input Commands ===
=== Scripted Input Commands ===
Line 33: Line 33:
<ul>
<ul>


<li><tt>(isOn)</tt> is "<tt>Yes</tt>" or "<tt>Yes</tt>" to include or to not include thermal conductivity calculations in the MPM analysis.
<li><tt>(isOn)</tt> is "<tt>Yes</tt>" or "<tt>No</tt>" to include or to not include thermal conductivity calculations in the MPM analysis.
</li>
</li>


</ul>
</ul>


<p>The optional parameters <tt>(option1)</tt>, <tt>(option2)</tt>, <i>etc</i>. can activate various physical mechanisms as heat sources that may cause temperature rises and may change the [[#Thermodynamics Modes in MPM|thermodynamics]] of the analysis. The options are
<p>The optional parameters <tt>(option1)</tt>, <tt>(option2)</tt>, <i>etc</i>. can activate various physical mechanisms as heat sources that may cause temperature rises and may change the [[#Thermodynamics Modes in MPM|thermodynamics]] of the analysis. The first two are complementary options (''i.e''., <tt>Adiabatic</tt> and <tt>Isothermal</tt> are complementary options, which means only one should be used. If both are omitted, the default setting is <tt>Isothermal</tt>.) are:


<ul>
<ul>
<li>"<tt>Adiabatic</tt>" (used to be <tt>Mechanical energy</tt>, which is still allowed) to be [[#Thermodynamics Modes in MPM|locally adiabatic]], which means any dissipated energy in a material results in a heat source that causes a temperature rise. This option can be used with conduction on or off. If conduction is off, the particle temperature will rise and will remain on the particle. When conduction is on, the heat may transfer to neighboring particles meaning with time, the particle will not be locally adiabatic. The entire problem, however, will be [[#Thermodynamics Modes in MPM|globally adiabatic]] unless there are [[#Thermal Boundary Conditions|thermal boundary conditions]].</li>
<li>"<tt>Adiabatic</tt>" (used to be <tt>Mechanical energy</tt>, which is still allowed) to be [[#Thermodynamics Modes in MPM|locally adiabatic]], which means any elastic thermal energy or dissipated energy in a material results in a corresponding change in the particle temperature, the  magnitude of which is determined by the material's heat capacity. This option can be used with conduction on or off. If conduction is off, the particle temperature rise will remain on the particle or the particle will be locally adiabatic. When conduction is on, the temperature rise can conduct as heat to neighboring particles meaning with time, the particle will not be locally adiabatic. The entire problem, however, will be [[#Thermodynamics Modes in MPM|globally adiabatic]] unless there are [[#Thermal Boundary Conditions|thermal boundary conditions]]. Finally, the available mechanics to generate thermal energy depend on [[#Material Support for Thermal Energy|material model being used]].</li>


<li>"<tt>Isothermal</tt>" to be [[#Thermodynamics Modes in MPM|locally isothermal]], which means the material response will never cause a temperature rise. Implicitly all dissipated energy is expelled to the exterior. The problem will be [[#Thermodynamics Modes in MPM|globally isothermal]] unless there are [[#Thermal Boundary Conditions|thermal boundary conditions]].</li>
<li>"<tt>Isothermal</tt>" to be [[#Thermodynamics Modes in MPM|locally isothermal]], which means the material response will never cause a temperature rise. Implicitly all dissipated energy is expelled to the exterior. The problem will be [[#Thermodynamics Modes in MPM|globally isothermal]] unless there are [[#Thermal Boundary Conditions|thermal boundary conditions]].</li>


<li>"<tt>Crack Tips</tt>" to have cracks convert total energy released into heat at the crack tips. For crack-tip heating to work correctly, you have to set the crack thickness to define its thickness. This option has no affect unless <tt>(isOn)</tt> is &quot;Yes.&quot; (<i>i.e.</i>, the temperature rises results from a heat source term in the conduction equations).</li>
</ul>
 
The following options activate conversion of various physical mechanisms into heat. Note that these mechanisms inject heat into the problem though heat flux term in the conduction equation. In other words, they have no affect unless conduction is activated by setting <tt>(isOn)</tt> is &quot;Yes&quot;:
 
<ul>
 
<li>"<tt>Crack Tips</tt>" to have cracks convert total energy released into heat at the crack tips. For crack-tip heating to work correctly, you have to set the crack thickness to define its thickness. </li>


<li>"<tt>Friction</tt>" to convert work of [[Friction|frictional sliding]] between materials in [[Multimaterial MPM|multimaterial mode]] into heat.</li>
<li>"<tt>Friction</tt>" to convert work of [[Friction|frictional sliding]] between materials in [[Multimaterial MPM|multimaterial mode]] into heat.</li>
Line 52: Line 58:


</ul>
</ul>
Note that <tt>Adiabatic</tt> and <tt>Isothermal</tt> are complementary options, which means only one should be used. If both are omitted, the default setting is <tt>Isothermal</tt>.


=== XML Input Commands ===
=== XML Input Commands ===
Line 85: Line 89:
: This command sets up a [[#Thermal Ramp|thermal ramp]] (see next section).
: This command sets up a [[#Thermal Ramp|thermal ramp]] (see next section).


== Thermal Ramp ==
=== Material Support for Thermal Energy ===
 
Mechanics calculations of thermal conduction are coupled phenomena, but the degree of coupling depends on the analysis method. The two ways they are coupled are:


One way to run simulations with initial thermal stresses is set the temperature to all particles at the start of the calculations to a temperature that differs from the [[#Stress Free Temperature|stress free temperature]]. Such an instantaneous temperature change is analogous to impact loading and may cause stress and strain oscillations. To avoid these dynamic effects, it is better to ramp up particle temperature difference by using a thermal ramp. In scripted input files, the commands are:
# Thermal expansion - a change in temperature causes most material the develop thermal strains that will change the particle stress.
# Thermal energy - various mechanical process can generate thermal energy that changes the particle temperature


ThermalRamp (diff),&lt;(time)&gt;
Most analytical methods (especially for elastic solids) include thermal expansion, but ignore thermal energy. This approximation is done because it decouples the thermal conduction equation for the mechanics equation and thereby simplifies the analysis. It is a good approximation for elastic solids because small-stain, elastic deformation does not generate much thermal energy, But when modeling large strains or modeling [[Material Models|materials]] that generate or dissipate more thermal energy (plastic or viscoelastic materials, liquids, and gases), thermal analyses should always be coupled, which involves adding heat source terms to the thermal conduction equation. Furthermore, when solving the two equations simultaneously by numerical methods, there is virtually no performance penalty to including all coupling terms (just the slight work to calculate the thermal energy).
RampStart (start time)


In <tt>XML</tt> input files, the command, which must be in the <tt><Thermal></tt> block is:
Nevertheless, proper modeling of coupled mechanical and thermal problems requires that the [[Material Models|material model]] implemented in <tt>NairnMPM</tt> correctly calculates all thermal energy terms. That calculation is not part of generic MPM because it depends on the [[Material Models|material model]] being used. The following lists explain which materials support various types of thermal coupling. Materials not in these lists may or may support these features as could be checked by looking at source code.


<Isothermal time="(time)" start="(start time)">(diff)</Isothermal>
* Thermal Expansion - All materials
* Dilational Energy - By thermodynamics, all materials heat when compressed and cool when expanded (gases a lot; solids very little). The follow materials support dilational energy:
** [[Isotropic Material|Isotropic]]
** [[Isotropic, Elastic-Plastic Material|IsoPlasticity]]
** [[Isotropic Softening Material|IsoSoftening]]
** [[Isotropic Plastic Softening Material|IsoPlasticSoftening]]
** [[Viscoelastic Material|Viscoelastic]]
** [[Mooney Material|Mooney]]
** [[Neo-Hookean Material|Neohookean]]
** [[Isotropic, Hyperelastic-Plastic Material|HEIsotropic]]
** [[Isotropic, Hyperelastic-Plastic Mie-Gr&#252;neisen Material|HEMGEOSMaterial]]
** [[Tait Liquid Material|TaitLiquid]]
** [[Ideal Gas Material|IdealGas]] (heat of expansion and compression, which is analog of elastic thermal energy in solids)
* Dissipated Energy - Some materials dissipate energy and calculations can cause that energy to be converted into heat. The following materials calculate their dissipated energy:
** [[Isotropic, Elastic-Plastic Material|IsoPlasticity]] (from plastic energy)
** [[Anisotropic, Elastic-Plastic Material|HillPlastic]] (from plastic energy)
** [[Viscoelastic Material|Viscoelastic]] (from viscoelastic dissipation and [[Common Material Properties#Artificial Viscosity|artificial viscosity]])
** [[Isotropic Softening Material|IsoSoftening]]
** [[Isotropic Plastic Softening Material|IsoPlasticSoftening]]
** [[Transversely Isotropic Softening Material|TransIsoSoftening]]
** [[Mooney Material|Mooney]] (from [[Common Material Properties#Artificial Viscosity|artificial viscosity]])
** [[Neo-Hookean Material|Neohookean]] (from [[Common Material Properties#Artificial Viscosity|artificial viscosity]])
** [[Isotropic, Hyperelastic-Plastic Material|HEIsotropic]] (from plastic energy and [[Common Material Properties#Artificial Viscosity|artificial viscosity]])
** [[Isotropic, Hyperelastic-Plastic Mie-Gr&#252;neisen Material|HEMGEOSMaterial]] (from plastic energy and [[Common Material Properties#Artificial Viscosity|artificial viscosity]])
** [[Tait Liquid Material|TaitLiquid]] (from shear work)


where
== Thermal Ramp ==
* <tt>(diff)</tt> is the final temperature difference to apply to all particles in degrees C (or K) after the ramp is completed.
 
* The optional <tt>(time)</tt> is the total time (in ms) to reach the final temperature difference. Enter a time less than zero to apply the entire temperature difference on the first time step. The default is to ramp in one step. It is usually better to ramp the temperature difference over a finite time determined by wave speed of the materials in the model.
To avoid thermal shocks, it is usually better to apply initial thermal fields using a [[PropertyRamp Custom Task|PropertyRamp custom task]].
* The optional <tt>(start time)</tt> is the time (in ms) to start the thermal ramp; the default is to start at time zero.


Thermal ramps can be used for residual stress calculations. Because they create uniform temperature changes they do not lead to any conduction. Thus a ramp can be used without doing conduction calculations or it can be used in combination with conduction calculations and other [[#Thermal Boundary Conditions|thermal boundary conditions]].
The scripting commands <tt>ThermalRamp</tt> and <tt>RampStart</tt> and the <tt>XML</tt> command <tt>&lt;Isothermal&gt;</tt> have been deleted and replaced by the new [[PropertyRamp Custom Task|PropertyRamp custom task]].


== Thermodynamics Modes in MPM ==
== Thermodynamics Modes in MPM ==


By using settings for [[#Conduction|conduction calculations]] along with optional [[#Thermal Boundary Conditions|thermal boundary conditions]], [[NairnMPM]] can simulate various thermodynamics conditions, and it important to be aware of the thermodynamics that is being models. The two particle modes are "Adiabatic" and "Isothermal" (set using a [[#Conduction|conduction option]]). These modes control how each particle handles localized energy flow. The full thermodynamics, however, also depends on interactions between particles and with the exterior. These options are:
By using settings for [[#Conduction|conduction calculations]] along with optional [[#Thermal Boundary Conditions|thermal boundary conditions]], [[NairnMPM]] can simulate various thermodynamics conditions, and it is important to be aware of the thermodynamics that is being modeled. The two basic thermodynamics modes refer to instantaneous temperature change on an isolated particle and are the  "Adiabatic" and "Isothermal" modes (set using a [[#Conduction|conduction option]]). Formally, the "Adiabatic" model is modeling the full dynamic problem while "Isothermal" mode is used to emulate quasi-static simulations. These modes control how each particle handles localized energy flow. The full thermodynamics, however, also depends on interactions between particles through conduction and on exchange of heat with the exterior through boundary conditions. These options are:


# [[#Conduction|Conduction]]: when conduction is one, heat can flow between particles; when it is off no heat flows between particles.
# [[#Conduction|Conduction]]: when conduction is on, heat can flow between particles; when it is off no heat flows between particles.
#  [[#Thermal Boundary Conditions|Thermal Boundary Conditions]]: an "Isolated" system means there are no [[#Thermal Boundary Conditions|thermal boundary conditions]] that can heat or change temperature of particles; a "Nonisolated" system has [[#Thermal Boundary Conditions|thermal boundary conditions]]
#  [[#Thermal Boundary Conditions|Thermal Boundary Conditions]]: an "Isolated" system means there are no [[#Thermal Boundary Conditions|thermal boundary conditions]] that can heat or change temperature of particles; a "Nonisolated" system has [[#Thermal Boundary Conditions|thermal boundary conditions]]


The following grid explains all combinations of these settings and describes the behavior of some thermodynamics quantities. In the table <b>dq</b> is cumulative heat flow on a single particle, <b>dQ</b> is heat flow for all particles, <b>dS</b> is total entropy change, and <b>dT</b> is change in average temperature.
The following grid explains all combinations of these settings and describes the behavior of some thermodynamics quantities. In the table <b>dq</b> is cumulative heat flow on a single particle, <b>dQ</b> is sum of heat flow for all particles, <b>dS</b> is total entropy change, and <b>dT</b> is change in average temperature.


<table border="1" cellpadding="2" cellspacing="0" width="90%" align="center">
<table border="1" cellpadding="2" cellspacing="0" width="90%" align="center">
<tr><th></th><th></th><th colspan="4">Adiabatic</th><th colspan="4">Isothermal</th>
<tr><th></th><th></th><th colspan="4">Adiabatic Mode</th><th colspan="4">Isothermal Mode</th>
</tr>
</tr>
<tr><th>System</th><th>Conduction</th><th>dq</th><th>dS</th><th>dQ</th><th>dT</th>
<tr><th>System (Thermal BCs)</th><th>Conduction (Particles)</th><th>dq</th><th>dQ</th><th>dS</th><th>dT</th>
<th>dq</th><th>dS</th><th>dQ</th><th>dT</th>
<th>dq</th><th>dQ</th><th>dS</th><th>dT</th>
</tr>
</tr>
<tr><th>Isolated</th><th>No</th>
<tr><th>Isolated</th><th>No (Isolated)</th>
<td align="center">0<sup>1</sup></td><td align="center">0</td><td align="center">0</td><td align="center">≠0</td>
<td align="center">0<sup>1</sup></td><td align="center">0</td><td align="center">0</td><td align="center">≠0</td>
<td align="center">≠0</td><td align="center">≠0</td><td align="center">≠0</td><td align="center">0</td>
<td align="center">≠0</td><td align="center">≠0</td><td align="center">≠0</td><td align="center">0</td>
</tr>
</tr>
<tr><th>Isolated</th><th>Yes</th>
<tr><th>Isolated</th><th>Yes (Nonisolated)</th>
<td align="center">≠0<sup>1</sup></td><td align="center">0</td><td align="center">≥0<sup>2</sup></td><td align="center">≠0</td>
<td align="center">≠0<sup>1</sup></td><td align="center">0</td><td align="center">≥0<sup>2</sup></td><td align="center">≠0</td>
<td align="center">≠0<sup>3</sup></td><td align="center">≠0<sup>3</sup></td>
<td align="center">≠0<sup>3</sup></td><td align="center">≠0<sup>3</sup></td>
<td align="center">≠0<sup>3</sup></td><td align="center">0<sup>3</sup></td>
<td align="center">≠0<sup>3</sup></td><td align="center">0<sup>3</sup></td>
</tr>
</tr>
<tr><th>Nonisolated</th><th>No</th>
<tr><th>Nonisolated</th><th>No (Isolated)</th>
<td rowspan="2" align="center">≠0</td><td  rowspan="2" align="center">≠0</td>
<td rowspan="2" align="center">≠0</td><td  rowspan="2" align="center">≠0</td>
<td  rowspan="2" align="center">≠0</td><td  rowspan="2" align="center">≠0<sup>4</sup></td>
<td  rowspan="2" align="center">≠0</td><td  rowspan="2" align="center">≠0<sup>4</sup></td>
Line 133: Line 162:
<td  rowspan="2" align="center">≠0</td><td  rowspan="2" align="center">≠0<sup>4</sup></td>
<td  rowspan="2" align="center">≠0</td><td  rowspan="2" align="center">≠0<sup>4</sup></td>
</tr>
</tr>
<tr><th>Nonisolated</th><th>Yes</th></tr>
<tr><th>Nonisolated</th><th>Yes (Nonisolated)</th></tr>
</table>
</table>


<ol>
<ol>
<li>If any particles start with a temperature that is different then the stress free temperature, the first time step will add dq = CV (dT<sub>i</sub> −dT<sub>SF</sub>) to the particle heat energy. The above conditions will hold thereafter, but a constant will be added to dS and dQ.</li>
<li>If any particles start with a temperature that is different then the stress free temperature, the first time step will add dq = C<sub>V</sub> (dT<sub>i</sub> −dT<sub>SF</sub>) to the particle heat energy. The above conditions will hold thereafter, but a constant will be added to dS and dQ.</li>


<li>The dq will be due to conduction only and entropy will increase because conduction is irreversible.</li>
<li>The dq will be due to conduction only and entropy will increase because conduction is irreversible.</li>
Line 164: Line 193:
The possible thermal boundary conditions are:
The possible thermal boundary conditions are:


* You can set [[Setting Velocity, Temperature, and Concentration#Temperature Conditions|temperature on the grid]].
* You can set [[Setting Velocity and Transport Values#Temperature Conditions|temperature on the grid]].
* [[Rigid Material|Rigid particles]] can provide moving temperature boundary conditions.
* [[Rigid Material|Rigid particles]] can provide moving temperature boundary conditions.
* You can set a [[Setting Forces and Fluxes#Heat Flux Conditions|heat flux on particle surfaces]].
* You can set a [[Setting Forces and Fluxes#Heat Flux Conditions|heat flux on particle surfaces]].
* A [[#Thermal Ramp|thermal ramp]] and ramp of the temperature of all particles
* A [[#Thermal Ramp|thermal ramp]] can ramp of the temperature of all particles


== Thermal Material Properties ==
== Thermal Material Properties ==


The results of thermal calculations will depend on the thermal properties of the materials in the model, which should all be set when doing such calculations. The relevant thermal properties are thermal expansion coefficients, thermal conductivity, and heat capacity. See the [[Material Models|material definitions]] for how to set these properties. A few [[Material Models|material models]] and [[Hardening Laws|hardening laws]] also depend on thermal transitions, such as a melting point, and on the [[#Stress Free Temperature|stress free temperature]].
The results of thermal calculations will depend on the thermal properties of the materials in the model, which should all be set when doing such calculations. The relevant thermal properties are thermal expansion coefficients, thermal conductivity, and heat capacity. See the [[Material Models|material definitions]] for how to set these properties. A few [[Material Models|material models]] and [[Hardening Laws|hardening laws]] also depend on thermal transitions, such as a melting point, and on the [[#Stress Free Temperature|stress free temperature]].

Latest revision as of 21:17, 30 January 2023

Thermal calculations in MPM include stresses and strains induced by temperature differentials and conduction induced by thermal gradients and boundary conditions.

Introduction

NairnMPM can do coupled elasticity-thermal conductivity calculations, provide thermal boundary conditions, convert various material processes into heat, and/or or can apply a constant temperature difference to all particles. In the presence of these temperature changes and heat input or output, the calculations will find thermal stresses and tracks heat flow for valid thermodynamics or, when needed, for tracking thermodynamic quantities. This section is for MPM only; see separate help on thermal calculations when doing FEA calculations with NairnFEA.

Stress Free Temperature

Thermal strains and stresses are always calculated relative to the stress free temperature, which is set in the StressFreeTemp command. In scripted input files, the command is:

StressFreeTemp (temp)

In XML input files, the command, which must be in the <MPMHeader> block is

<StressFreeTemp>(temp)</StressFreeTemp>

where (temp) is the stress free temperature in degrees. The default value is 0. Note that some Material Models and Traction Laws depend on the stress free temperature on an absolute temperature scale and therefore calculations that use those models should always set this temperature in degrees Kelvin.

Note that setting initial particle temperature different than the stress free temperature will cause strains to immediately evolve toward the changed state. The net effect will be an instantaneous "impact" of temperature change that might cause undesirable dynamic effects.

Conduction

The Conduction command activates thermal conductivity calculations and also activates various optionals involving heat flow in the simulations or other mechanisms that cause temperature changes. These settings can change and select the thermodynamics mode of the simulation. By default, conduction uses update methods analogous to FLIP methods used in mechanics. This update, however, sometimes results in temperature oscillations on particles within one cell. Conduction simulations with oscillations can be improved by using periodic FMPM(k) for conduction updates.

Scripted Input Commands

In scripted input files, the command is

Conduction (isOn),<(option1)>,<(option2)>...

where

  • (isOn) is "Yes" or "No" to include or to not include thermal conductivity calculations in the MPM analysis.

The optional parameters (option1), (option2), etc. can activate various physical mechanisms as heat sources that may cause temperature rises and may change the thermodynamics of the analysis. The first two are complementary options (i.e., Adiabatic and Isothermal are complementary options, which means only one should be used. If both are omitted, the default setting is Isothermal.) are:

  • "Adiabatic" (used to be Mechanical energy, which is still allowed) to be locally adiabatic, which means any elastic thermal energy or dissipated energy in a material results in a corresponding change in the particle temperature, the magnitude of which is determined by the material's heat capacity. This option can be used with conduction on or off. If conduction is off, the particle temperature rise will remain on the particle or the particle will be locally adiabatic. When conduction is on, the temperature rise can conduct as heat to neighboring particles meaning with time, the particle will not be locally adiabatic. The entire problem, however, will be globally adiabatic unless there are thermal boundary conditions. Finally, the available mechanics to generate thermal energy depend on material model being used.
  • "Isothermal" to be locally isothermal, which means the material response will never cause a temperature rise. Implicitly all dissipated energy is expelled to the exterior. The problem will be globally isothermal unless there are thermal boundary conditions.

The following options activate conversion of various physical mechanisms into heat. Note that these mechanisms inject heat into the problem though heat flux term in the conduction equation. In other words, they have no affect unless conduction is activated by setting (isOn) is "Yes":

  • "Crack Tips" to have cracks convert total energy released into heat at the crack tips. For crack-tip heating to work correctly, you have to set the crack thickness to define its thickness.
  • "Friction" to convert work of frictional sliding between materials in multimaterial mode into heat.
  • "Crack Friction" to convert work of frictional sliding between crack surfaces into heat.

XML Input Commands

In XML input files, conduction is activated with a <Thermal> block:

<Thermal>
  <Conduction/>
  <EnergyCoupling/>
  <CrackTipHeating/>
  <CrackContactHeating/>
  <ContactHeating/>
  <Isothermal time="5" start="2">-200</Isothermal>
</Thermal>

where

<Conduction>
This tag activates condition calculations. Its omission turns conduction off. When conducton is activated, you should set all needed thermal material properties.
<EnergyCoupling/>
The presence of this tag is equivalent to the "Adiabatic" option while its omission is equivalent to the "Isothermal" option (see above).
<CrackTipHeating/>
The presence of this tag is equivalent to the "Crack Tips" option (see above).
<CrackContactHeating/>
The presence of this tag is equivalent to the "Crack Friction" option (see above).
<ContactHeating/>
The presence of this tag is equivalent to the "Friction" option (see above).
<Isothermal>
This command sets up a thermal ramp (see next section).

Material Support for Thermal Energy

Mechanics calculations of thermal conduction are coupled phenomena, but the degree of coupling depends on the analysis method. The two ways they are coupled are:

  1. Thermal expansion - a change in temperature causes most material the develop thermal strains that will change the particle stress.
  2. Thermal energy - various mechanical process can generate thermal energy that changes the particle temperature

Most analytical methods (especially for elastic solids) include thermal expansion, but ignore thermal energy. This approximation is done because it decouples the thermal conduction equation for the mechanics equation and thereby simplifies the analysis. It is a good approximation for elastic solids because small-stain, elastic deformation does not generate much thermal energy, But when modeling large strains or modeling materials that generate or dissipate more thermal energy (plastic or viscoelastic materials, liquids, and gases), thermal analyses should always be coupled, which involves adding heat source terms to the thermal conduction equation. Furthermore, when solving the two equations simultaneously by numerical methods, there is virtually no performance penalty to including all coupling terms (just the slight work to calculate the thermal energy).

Nevertheless, proper modeling of coupled mechanical and thermal problems requires that the material model implemented in NairnMPM correctly calculates all thermal energy terms. That calculation is not part of generic MPM because it depends on the material model being used. The following lists explain which materials support various types of thermal coupling. Materials not in these lists may or may support these features as could be checked by looking at source code.

Thermal Ramp

To avoid thermal shocks, it is usually better to apply initial thermal fields using a PropertyRamp custom task.

The scripting commands ThermalRamp and RampStart and the XML command <Isothermal> have been deleted and replaced by the new PropertyRamp custom task.

Thermodynamics Modes in MPM

By using settings for conduction calculations along with optional thermal boundary conditions, NairnMPM can simulate various thermodynamics conditions, and it is important to be aware of the thermodynamics that is being modeled. The two basic thermodynamics modes refer to instantaneous temperature change on an isolated particle and are the "Adiabatic" and "Isothermal" modes (set using a conduction option). Formally, the "Adiabatic" model is modeling the full dynamic problem while "Isothermal" mode is used to emulate quasi-static simulations. These modes control how each particle handles localized energy flow. The full thermodynamics, however, also depends on interactions between particles through conduction and on exchange of heat with the exterior through boundary conditions. These options are:

  1. Conduction: when conduction is on, heat can flow between particles; when it is off no heat flows between particles.
  2. Thermal Boundary Conditions: an "Isolated" system means there are no thermal boundary conditions that can heat or change temperature of particles; a "Nonisolated" system has thermal boundary conditions

The following grid explains all combinations of these settings and describes the behavior of some thermodynamics quantities. In the table dq is cumulative heat flow on a single particle, dQ is sum of heat flow for all particles, dS is total entropy change, and dT is change in average temperature.

Adiabatic ModeIsothermal Mode
System (Thermal BCs)Conduction (Particles)dqdQdSdT dqdQdSdT
IsolatedNo (Isolated) 0100≠0 ≠0≠0≠00
IsolatedYes (Nonisolated) ≠010≥02≠0 ≠03≠03 ≠0303
NonisolatedNo (Isolated) ≠0≠0 ≠0≠04 ≠0≠0 ≠0≠04
NonisolatedYes (Nonisolated)
  1. If any particles start with a temperature that is different then the stress free temperature, the first time step will add dq = CV (dTi −dTSF) to the particle heat energy. The above conditions will hold thereafter, but a constant will be added to dS and dQ.
  2. The dq will be due to conduction only and entropy will increase because conduction is irreversible.
  3. This mode is identical to system isolated with no conduction if the initial temperature of all particles is equal to the stress-free temperature. Thus, such a simulation should turn off conduction for efficiency because it is not needed.
  4. When conduction is "No", the only possible thermal boundary condition is a thermal ramp that applies uniform temperature change to all particles. If there is no thermal ramp, then the system is isolated and refer to first mode instead.

Tracking Thermodynamic Quantities

In any thermodynamics mode, you can easily track various thermodynamics quantities:

Thermal Boundary Conditions

The possible thermal boundary conditions are:

Thermal Material Properties

The results of thermal calculations will depend on the thermal properties of the materials in the model, which should all be set when doing such calculations. The relevant thermal properties are thermal expansion coefficients, thermal conductivity, and heat capacity. See the material definitions for how to set these properties. A few material models and hardening laws also depend on thermal transitions, such as a melting point, and on the stress free temperature.