/ALE/GRID/STANDARD

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

It is an improved /ALE/GRID/SPRING formulation based on edge springs and anti-shear springs. 1

Format

(1)	(3)	(5)	(7)
/ALE/GRID/STANDARD
$α$ $α$	$γ$ $γ$	$η$ $η$	$l_{c}$ $l_{c}$
Blank Format

Definition

Field	Contents	SI Unit Example
$α$ $α$	Scale factor for maximum stiffness. 2 Default = 0.9 (Real)
$γ$ $γ$	Nonlinearity factor for edge spring stiffness. 3 Default = 1e-2 (Real)
$η$ $η$	Damping coefficient. 4 Default = 1e-2 (Real)
$l_{c}$ $l_{c}$	Characteristic length. (Real)	$[m]$ $[m]$

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.
1. Edge springs
  The forces for an edge spring are a function of its length variation during time.(1)
  $Δ F_{e d g e} = k (h) \cdot (w_{2} - w_{1}) d t$ $Δ F_{e d g e} = k (h) \cdot (w_{2} - w_{1}) d t$
  
  Where,
  
  $w_{1}, w_{2}$ $w_{1}, w_{2}$ are grid velocities on nodes N₁ and N₂, respectively.
  
  $h$ $h$ is the N₁ distance from opposite face
  
  $d t$ $d t$ is the time step
  
  and $k (h)$ $k (h)$ is the spring stiffness $k (h) = k_{c r i t i c a l}$ $k (h) = k_{c r i t i c a l}$
  
  If $h$ $h$ is inferior to the characteristic length $l_{c}$ $l_{c}$ and N₁ is moving toward the opposite face then,(2)
  $k (h) = \frac{1}{λ^{2}} [γ + (γ - 1) {(\frac{h - l_{c}}{l_{c}})}^{3}] k_{c r i t i c a l}$ $k (h) = \frac{1}{λ^{2}} [γ + (γ - 1) {(\frac{h - l_{c}}{l_{c}})}^{3}] k_{c r i t i c a l}$
  
  $\frac{1}{λ^{2}}$ $\frac{1}{λ^{2}}$ is the stability factor taking into account the damping factor $β$ $β$ , the scale factor $α$ $α$ , and time step $d t$ $d t$ 4
  
  Figure 2.
  
  otherwise, $k (h) = \frac{γ}{λ^{2}} k_{c r i t i c a l}$ $k (h) = \frac{γ}{λ^{2}} k_{c r i t i c a l}$
  
  Figure 3.
2. Anti-shear springs
  The anti-shear forces $F_{s h e a r}$ $F_{s h e a r}$ are computed from node penetration. gap is $\frac{l_{c}}{s}$ $\frac{l_{c}}{s}$ from opposite face.
  
  The value of $F_{s h e a r}$ $F_{s h e a r}$ is:
  (3)
  $F_{s h e a r} = g a p \cdot k (h)$ $F_{s h e a r} = g a p \cdot k (h)$
  
  and (4)
  $k (h) = \frac{γ}{λ^{2}} [γ + (γ - 1) {(\frac{h - l_{c}}{l_{c}})}^{3}] k_{c r i t i c a l}$ $k (h) = \frac{γ}{λ^{2}} [γ + (γ - 1) {(\frac{h - l_{c}}{l_{c}})}^{3}] k_{c r i t i c a l}$
  
  Figure 4.
3. Viscous Damping
  Viscous forces are computed from a critical damping corresponding to the upper bound for stiffness: $\frac{1}{λ^{2}}$ $\frac{1}{λ^{2}}$ (5)
  $F_{v i s c o u s} = β α (\sqrt{1 + β^{2}} - β) (w_{2} - w_{1})$ $F_{v i s c o u s} = β α (\sqrt{1 + β^{2}} - β) (w_{2} - w_{1})$
4. Grid Velocity
  The grid velocity is then updated according to:(6)
  $w_{n + 1} = w_{n} + \frac{(Δ F_{e d g e} + F_{s h e a r} + F_{v i s c o u s}) d t}{m}$ $w_{n + 1} = w_{n} + \frac{(Δ F_{e d g e} + F_{s h e a r} + F_{v i s c o u s}) d t}{m}$
  
  Where, $m$ $m$ is fictitious mass on node from springs (automatically computed during Starter).
Increasing $α = 1$ $α = 1$ , the maximum stiffness will be increased. The scale factor $α$ $α$ 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 $γ$ $γ$ = 0. Moreover, increasing $γ$ $γ$ , 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 $γ$ $γ$ .
It is recommended to use small values for $β$ $β$ , otherwise damping may become over critical. The stability factor is:(7)
$λ = \frac{d t}{α_{0} (\sqrt{1 + β^{2}} - β)}$ $λ = \frac{d t}{α_{0} (\sqrt{1 + β^{2}} - β)}$
$l_{c}$ defines the length below which:
- Edge spring stiffness is increased $h < l_{c}$
- Anti-shear spring is activated: $h < \frac{l_{c}}{5}$
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 $α$ 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.