This section describes how to include a coordinate transformation to the solution of the incompressible Navier-Stokes equations. In some cases, this approach allows a slightly deformed geometry to be mapped into a geometry with a homogeneous direction, which can be treated using a quasi-3D method. It is also useful for problems with a moving body, where otherwise a moving mesh would have to be employed.
To activate the mapping technique, SolverType
needs to be set as:
Also, there are other optional properties in the SolverInfo
section:
the first two options determine if the pressure and viscous terms resulting from the coordinate transformation are treated implicitly using an iterative procedure. If the last option is set to true, the viscous terms in the mapping are not computed. This leads to a faster solution, but the effect on the results need to be determined for the specific case. By default, all mapping terms are computed and treated explicitly.
When treating the mapping terms implicitly, the following parameters can be set:
They determine the tolerance of the iterative solution of the equations, and a relaxation parameter which can improve the numerical stability of the method.
The particular transformation employed is specified by:
where TIMEDEPENDENT
indicates if the transformation varies with time.
The available values for TYPE
, and the transformations they represent, are:
Mapping type | |||
Translation | x + f(t) | y + g(t) | z + h(t) |
XofZ | x + f(z,t) | y | z |
XofXZ | f(x,z,t) | y | z |
XYofZ | x + f(z,t) | y + g(z,t) | z |
XYofXY | f(x,y,t) | g(x,y,t) | z |
General | f(x,y,z,t) | g(x,y,z,t) | h(x,y,z,t) |
where (
, , ) are the Cartesian (physical) coordinates and (x,y,z) are the transformed coordinates. Note that for quasi-3D problems, the z coordinate cannot be transformed.
The function COORDS
(and VEL
for time dependent mappings) indicated in the MAPPING
section
need to be defined, for example as:
the transformation defined by these functions need to be valid (non-zero Jacobian). By
default, any component of COORDS
that is not specified is taken as a trivial transformation, e.g.
= x, and any velocity not specified is considered to be zero.
In case of a time-dependent mapping, a moving body boundary condition is available:
When using the MovingBody
boundary condition, the Dirichlet condition is relative to the
boundary, while the regular Dirichlet boundary condition is taken in an absolute
sense.
All Dirichlet boundary conditions are specified in the Cartesian (physical) space, and are automatically transformed to the computational frame of reference.
NEKTAR/solvers/IncNavierStokesSolver/Tests
. See for example the files
KovaFlow_3DH1D_P8_16modes_Mapping-implicit.xml
and CylFlow_Mov_mapping.xml
.