HyperWorks Solvers

/UPWM/SUPG

/UPWM/SUPG

 Previous topic Next topic No expanding text in this topic  

/UPWM/SUPG

 Previous topic Next topic JavaScript is required for expanding text JavaScript is required for the print function  

Engine Keyword

/UPWM/SUPG - Upwind Method for Momentum Advection

Description

This describes the Streamline Upwind Petrov Galerkin formulation.

Format

/UPWM/SUPG

fac

Data

Description

fac

Scale factor applied on Petrov Galerkine weighting function.

Default = 1.00 (recommended value).

Comments

1.The streamline upwind Petrov-Galerkine (SUPG) is an optimal upwinding technique which acts only in the flow direction.
2.If this keyword is used, the usual upwind method for momentum advection is no longer used. As a consequence, the coefficient symbol_n1 from /UPWIND Starter is also not used.
3.The Galerkin (shape function upwnsupg_sym) are modified to integrate the momentum convection term upwnsupg_eq, which leads to a Petrov Galerkin method using weight function W* such as:

upwn_supg_equation

Where, symbol_k_slash is a function of the Reynolds number, material velocity and element length, and fac is a scale factor. fac = 1.00 is the standard value for SUPG formulation (recommended).

4.The /UPWM/SUPG and /UPWM/TG options are available in 3D Euler, 3D ALE and 2D Euler, although they are not recommended in axisymmetric cases.
5.It is available for multi-material law 51.
6.The upwinding technique is introduced to add numerical diffusion to the scheme which, otherwise is generally under-diffusive and thus unstable. The usual method is first: order accurate and strongly depends on the mesh aspect. In some cases (see examples below) the classical method might lead to locally erratic velocities with unsuitable meshes. SUPG method provides more accurate results, if the phenomenon is sensitive to mesh aspect. It is recommended to systematically activate it.

Modeling situation:

blast_wave

Figure 1: Blast Wave modeling
A regular mesh is used although this phenomenon has spherical evolution. Result are obtained with multi-material law 51.

Velocity Result - Classical UPWIND:

classical_upwind

Figure 2: Blast Wave evolution with usual method
Classical UPWIND method for momentum transportation is not suitable with spherical evolution and regular mesh.

Velocity Result - SUPG:

blast_wave_supg

Figure 3: Blast Wave evolution with SUPG method
SUPG method for momentum advection removes mesh sensitivity. It is no longer required to have a spherical mesh to model spherical phenomena.

See Also:

Space Integration in User's Guide