# /ALE/GRID/STANDARD

Block Format Keyword Describes the standard formulation for ALE grid velocity computation.

## Format

(1) | (2) | (3) | (4) | (5) | (6) | (7) | (8) | (9) | (10) |
---|---|---|---|---|---|---|---|---|---|

/ALE/GRID/STANDARD | |||||||||

$\alpha $ | $\gamma $ | $\eta $ | ${l}_{c}$ | ||||||

Blank Format |

## Definitions

## Comments

- Fictitious springs are introduced on solid elements to control grid
velocities.These springs are nonlinear elastic viscous. To ensure stability, their stiffness is computed from time step. The two types of springs are edge and anti-shear springs.
- Edge springsThe forces for an edge spring are a function of its length variation during time.
(1) $$\text{\Delta}{F}_{edge}=\mathrm{k}\left(h\right)\cdot \left({w}_{2}-{w}_{1}\right)dt$$Where,

${w}_{1},{w}_{2}$ are grid velocities on nodes

`N`_{1}and`N`_{2}, respectively.$h$ is the

`N`_{1}distance from opposite face$dt$ is the time step

and $k\left(h\right)$ is the spring stiffness $\mathrm{k}\left(h\right)={k}_{critical}$

If $h$ is inferior to the characteristic length ${l}_{c}$ and`N`_{1}is moving toward the opposite face then,(2) $$\mathrm{k}\left(h\right)=\frac{1}{{\lambda}^{2}}\left[\gamma +\left(\gamma -1\right){\left(\frac{h-{l}_{c}}{{l}_{c}}\right)}^{3}\right]{k}_{critical}$$$\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{${\lambda}^{2}$}\right.$ is the stability factor taking into account the damping factor $\beta $ , the scale factor $\alpha $ , and time step $dt$ 4otherwise, $\mathrm{k}\left(h\right)=\frac{\gamma}{{\lambda}^{2}}{k}_{critical}$ - Anti-shear springs
The anti-shear forces ${F}_{shear}$ are computed from node penetration. gap is $\raisebox{1ex}{${l}_{c}$}\!\left/ \!\raisebox{-1ex}{$s$}\right.$ from opposite face.

The value of ${F}_{shear}$ is:

(3) $${F}_{shear}=gap\cdot \mathrm{k}\left(h\right)$$and(4) $$\mathrm{k}\left(h\right)=\frac{\gamma}{{\lambda}^{2}}\left[\gamma +\left(\gamma -1\right){\left(\frac{h-{l}_{c}}{{l}_{c}}\right)}^{3}\right]{k}_{critical}$$ - Viscous DampingViscous forces are computed from a critical damping corresponding to the upper bound for stiffness: $\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{${\lambda}^{2}$}\right.$
(5) $${F}_{viscous}=\beta \alpha \left(\sqrt{1+{\beta}^{2}}-\beta \right)\left({w}_{2}-{w}_{1}\right)$$ - Grid VelocityThe grid velocity is then updated according to:
(6) $${w}_{n+1}={w}_{n}+\frac{\left(\text{\Delta}{F}_{edge}+{F}_{shear}+{F}_{viscous}\right)dt}{m}$$Where, $m$ is fictitious mass on node from springs (automatically computed during Starter).

- Edge springs
- Increasing $\alpha =1$ , the maximum stiffness will be increased. The scale factor $\alpha $ determines the maximum stiffness for a given spring at zero length. The scale factor ensures that the critical stability value is not exceeded (to avoid time step decrease).
- This flag is acting on stiffness shape. Stiffness is linear with $\gamma $ = 0. Moreover, increasing $\gamma $ , the lower bound stiffness for edge spring will be increased. Springs have a critical stiffness at zero length (this corresponds to a unitary factor). For a length greater than or equal to the characteristic length, the spring stiffness is the critical stiffness multiplied by $\gamma $ .
- It is recommended to use small
values for
$\beta $
, otherwise damping may become over critical. The stability factor
is:
(7) $$\lambda =\frac{dt}{{\alpha}_{0}\left(\sqrt{1+{\beta}^{2}}-\beta \right)}$$ -
${l}_{c}$
defines the length below which:
- Edge spring stiffness is increased $h<{l}_{c}$
- Anti-shear spring is activated: $h<\raisebox{1ex}{${l}_{c}$}\!\left/ \!\raisebox{-1ex}{$5$}\right.$

- All these parameters can be modified during an Engine restart (/ALE/GRID/STANDARD).
- Mesh auto correction. It is
possible to give more weight to anti-shear forces by either:
- Setting ${l}_{c}$ parameter close to the mesh size
- Setting a negative value for $\alpha $ parameter (elastic forces on edges are set to 0 at the first cycle of current run)

- This method assumes a homogeneous spring repartition around each node. This is not the case when connecting two meshes, where topologies are different.