GASDYN_CorrectEnergy

correct total energy in GASDYN application

GASDYN_CorrectEnergy = 0.01 # allow 1% of the defect energy to be corrected during a time step
Default: GASDYN_CorrectEnergy = 0.0 The correction preocedure is \begin{align} \rho E_i^{corrected} = \rho E_i + \gamma W^{\rho E}_i\end{align}

Classical case:

Choose GASDYN_CorrectEnergy > 0 In this case, we define \( W^{\rho E}_i := {\rho E}_i\) and determine \( \gamma\) such that \begin{align} \int \limits_{\Omega} \gamma W^{\rho E} dV \approx \gamma \sum \limits_{\Omega} W^{\rho E}_i V_i = \text{GasdynCorrectTotalEnergy} \cdot \text{DefectTotalEnergy}\end{align} Remark: This option scales the total energy by a factor such that the correct analytical energy is assumed.

Expert option:

Choose GASDYN_CorrectEnergy < 0 In this case, we define \( W^{\rho E}_i := {\rho E}_i \vert \nabla^T \mathbf{v} \vert \Delta t\) and determine \( \gamma\) such that \begin{align} \int \limits_{\Omega} \gamma W^{\rho E} dV \approx \gamma \sum \limits_{\Omega} W^{\rho E}_i V_i = \text{DefectTotalEnergy}\end{align} with the requirement \( \vert \gamma \vert < \vert \text{GasdynCorrectTotalEnergy} \vert\) Remark: This option scales the current dynamic changes of the total energy by a factor such that the correct analytical energy is assumed.
This item is referenced in:
GASDYN_CorrectEnergy correct total energy in GASDYN application