MAIA bb96820c
Multiphysics at AIA
Loading...
Searching...
No Matches
Level Set (LS)

levelSetCtrlPntMethod

MFloat LsCartesianSolver::m_GCtrlPntMethod
default = 2

Choose control point methods.
Possible values are:

  • 1 for single body.
  • 2 for multiple bodies.

Keywords: LEVELSET

nogRKSteps

MInt LsCartesianSolver::initializeIntegrationScheme::m_nogRKSteps
default = 0

Sets the number of generalized Runge Kutta steps.

  • non-negative integers

Keywords: LEVELSET, INTEGRATION, RUNGE KUTTA

grkalpha-step

MFloat LsCartesianSolver::m_gRKalpha
Sets the alpha coefficient for the Runge Kutta steps used in the integration scheme. Should contain m_nogRKSteps values. Keywords: LEVELSET, RUNGEKUTTA

radius

MFloat R0
default: 0.5

Radius of moving cylinder in the the level-set test-case 450.
Keywords: LEVEL_SET, BODY

GFieldFromSTLInitMode

default = 1

Choose method how to generate Levelset from STL.
In the Initialization-part, possible values are:

  • 0: only G0-cells are initialized by the stl
  • 1: G-cells specified by the ContstructFlag are initialized by the stl (for ContstructFlag=1: all G-cells, for ContstructFlag=2/3: all Band-cells,)
  • 2: G-cells specified by the ContstructFlag are initialized by a sphere (with radius 0.5 in the origin). Very special testcase, barely useful!
  • 3: only G0-cells are initialized by the stl

gReinitIterationsForGFieldFromSTL

MInt FvCartesianSolver::m_gReinitIterations
default = no default value

This property triggers the maximum number of reinitialization steps to reach the convergence criterion same functionality as gReinitIterations in lssolver.cpp: m_gReinitIterations reads in the property before calling the function levelSetConstrainedReinitialization after the "Reinit", return the old value to m_gReinitIterations possible values are:

  • Any non-negative integer value

Keywords: LEVELSET, REINITIALIZATION, ITERATION, STEPS

levelSetMaxCorrectionThresholdLevels

MInt LsCartesianSolver::spatiallyAdaptiveCorrectionFromSTL::maxThresholdLevels
default = 3
spatiallyAdaptiveCorrectionFromSTL does a correction based on the curvature of the geometry. The high curvature regions are corrected, while the low curvature regions are left untouched. The number of the corrected is controlled by the multiple level thresholds. remarks: this works only with control 2 (moving STL).
Possible values are:

  • 1 to ?

Keywords: LEVEL SET, STL, CURVATURE, GEOMETRY

amplitudes

MFloat* LsCartesianSolver::m_static_computeBodyProperties_amplitude
default = none

Amplitude for body motion of embedded bodies.
NOTE: also used in FV-MB solver for some special cases.

Possible values are:

  • list of floating point numbers

Keywords: LEVELSET, MOVING, BODY, BODY_MOTION

freqFactors

MFloat* LsCartesianSolver::m_static_computeBodyProperties_freqFactor
default = none

Set the frequency factors for prescribing body motion for all embedded bodies.
NOTE: also used in FV-MB solver for some special cases.

Possible values are:

  • list of positive floating point numbers

Keywords: LEVELSET, MOVING, BODY, BODY_MOTION

bodyMovementFunctions

MInt* LsCartesianSolver::m_static_computeBodyProperties_bodyToFunction
default = 1

Prescribes the functions for the body movement. Check the switch case in computeBodyProperties() for what each case actually does.

Possible values are:

  • integer from 1 to 7

Keywords: LEVELSET, BODY, BODY_MOTION, MOVING

initialBodyCenters

MFloat LsCartesianSolver::initialBodyCenter
For each body, nDim Float values. With this property, one can move an STL around to its initial position. In initialInsidePoints is the "real" center of the stl files. This property, is the movement from the "real" center to the initial center used during calculation. Keywords: LEVELSET, MULTILEVELSET, MB

liftStartAngles1

MFloat LsCartesianSolver::liftStartAngle1
default = 0.0

For each body, sets the start angle of the translation. The translation is described by a:

  • bodyToFunction case 1: cosine function
  • bodyToFunction case 2: valve lift shifted quadratic sine function (first angle)
    • Between 0 and 4.0 * PI
    Keywords: LEVELSET, MULTILEVELSET, MB

liftStartAngles2

MFloat LsCartesianSolver::liftStartAngle2
default = 3.0 * PI

For each body, sets the start angle of the translation. The translation is described by a:

  • bodyToFunction case 2: valve lift shifted quadratic sine function (second angle)
    • Between 0 and 4.0 * PI
    Keywords: LEVELSET, MULTILEVELSET, MB

liftEndAngles1

MFloat LsCartesianSolver::liftEndAngle1
default = PI

For each body, sets the end angle of the translation.

  • Between 0 and 4.0 * PI

Keywords: LEVELSET, MULTILEVELSET, MB

liftEndAngles2

MFloat LsCartesianSolver::liftEndAngle2
default = 4.0 * PI

For each body, sets the end angle of the translation.

  • Between 0 and 4.0 * PI

Keywords: LEVELSET, MULTILEVELSET, MB

circleStartAngles

MFloat LsCartesianSolver::circleStartAngle
default = 0.0

For each body, sets the start angle for the circular motion case. (i.e. bodyToFunction case 5).

  • Between 0 and 4.0 * PI

Keywords: LEVELSET, MULTILEVELSET, MB

rotAngle

MFloat LsPar::rotAngle
default = "0.0"

Used to rotate the primary direction of the body movement. information in the property-file is expected in degree!

  • Any positive Float

Keywords: LEVELSET EMBEDED BOUNDARY, MOVEMENT FUNCTIONS

initialInsidePoints

MFloat LsCartesianSolver::initialInsidePoints[s_maxNoEmbeddedBodies*3]
default = 0.0

Defines a reference point that is located inside the body for each embedded body (in x-,y-,z-coordinates) at the beginning of the calculation.
Keywords: LEVEL_SET

gapRegionNormals

Default: 0, 0, 0
Set the normal vectors of the gap region.
Possible values are:

  • Coordinates [noGapRegions * nDim]

Keywords: LEVEL-SET, GAP REGION

gapRegionCenters

Default: 0, 0, 0
Set the center of the gap region.
Possible values are:

  • Coordinates [m_noGapRegions * nDim]

Keywords: DISCONTINUOUS_GALERKIN, APE

gapRegionRadii

MFloat LsCartesianSolver::setUpPotentialGapCells::radius
default: 1.0

Set the radius for each gap region, where any cell within the radius around the gap center will be treated as a potencial gap cell.
Keywords: FINITE_VOLUME, GAP CELL

noGapRegionsClose

MInt LsCartesianSolver::setUpPotentialGapCells::noGapRegionsClose
Default: 1
Sets the number of gap regions, which will be closed by the solver.
Possible values are:

  • Any positive integer.

Keywords: DISCONTINUOUS_GALERKIN, APE, GAP CELL

gapRegionRadiiClose

MFloat LsCartesianSolver::setUpPotentialGapCells::radiusClose
default: 1.0

Set the radius for each gap closure region, where any cell within the radius around the gap closure center will be treated as a potencial gap closure cell.
Keywords: FINITE_VOLUME, GAP CELL

gapRegionHeightsClose

MFloat LsCartesianSolver::setUpPotentialGapCells::gapRegionHeightsClose
Default: none
Sets the heights of gap regions, which will be closed by the solver.
Possible values are:

  • floating point values

Keywords: DISCONTINUOUS_GALERKIN, APE, GAP CELL

gapRegionBodyClose

MInt LsCartesianSolver::setUpPotentialGapCells::gapRegionBodyClose
Default: none
Sets the gap regions on the body which will be closed by the solver.
Possible values are:

  • integers

Keywords: DISCONTINUOUS_GALERKIN, APE, GAP CELL

referenceLength

MFloat LsCartesianSolver::m_referenceLength
default = 1.0

WARNING: Do NOT use any value different than 1.0 - The correct implementation of this is not checked, so it probably will not do what you think it does/should do. Don't use it unless you REALLY know what you are doing. Reference Length L - The length = 1.0 of the grid is scaled with L. Possible values are:

  • 1.0 +- eps

Keywords: LEVELSET, VARIABLES

referenceLength

MFloat FvCartesianSolver::m_Pr
default = sqrt(1/(1+0.5*(1.4-1)/2*Ma^2

Reference velocity used for non-dimensionalisation of some of the movement functions! NOTE: this is not consistently used, some versions use the velocity, other the speed of sound (both in the infinity state) as reference Velocity! Keywords: LEVELSET, VARIABLES

virtualSurgery

MBool LsCartesianSolver::m_virtualSurgery
default = false

Activates the virtual surgery environment of the level-set solver. When activated, the solver modifies the 0th level-set field by interpolating between the 1st and 2nd level-set field. The 1st and 2nd level-set field are static. The 0th level-set field must not be a collected level-set. Keywords: LEVELSET, VARIABLES

sphereRadiusLimit

MFloat LsCartesianSolver::m_sphereRadiusLimit
default = 5.0

The sphereRadiusLimit defines the maximum area in which the level-set is defined. This factor can be set, if the band should be small to save cells, but the level-set value outside of the band should still be calculated. Is only used in mode 4 of constructGFieldFromSTL. Keywords: LEVELSET, VARIABLES

approxNoInterpRegions

MInt LsCartesianSolver::m_approxNoInterpRegions
default =

Number of Interpolations regions in virtual surgery environment. The number can only be approximated a priori. Keywords: LEVELSET, VARIABLES

interpStartTimes

MInt LsCartesianSolver::m_interpStartTime
default =

This property defines the start time step of a virtual surgery for each interpolation region. Keywords: LEVELSET, VARIABLES

noInterpTimeSteps

MInt LsCartesianSolver::m_noInterpTimeSteps
default =

This property defines the number of time steps used for the interpolation in a virtual surgery for each interpolation region. Keywords: LEVELSET, VARIABLES

forcing

MInt LsCartesianSolver::m_forcing
default = no default

Activates the excitation / forcing via sponge in 2D.

Possible values are:

  • 0: inactive
  • 1: active

Keywords: LEVELSET, FLAME, SPONGE

jetHalfLength

MFloat FvCartesianSolver::m_jetHalfLength
default = 4.165

parameter used in generation of various jet inflows possible values are:

  • any positive float

Keywords: jet

steadyFlameLength

MFloat LsCartesianSolver::m_steadyFlameLength
no default

This variable defines the steady flame front amplitude
possible values are:

  • Non-negative floating point values

Keywords: COMBUSTION, STEADY, FLAME, FRONT, AMPLITUDE

maxNoLevelSets

MInt LsCartesianSolver::m_maxNoSets
default = 1

This property fixes the maximum number of separate level set functions to be computed cells
possible values are:

  • Positive integer values, should be O(1) -> more than 10 separate sets will become very slow and should be unnencessary

Keywords: LEVELSET, MULTIPLE LEVEL SET FUNCTIONS

GFieldFromSTLReInitMode

default = 2

Choose method how to Reinit Levelset from STL.
In the Reinitialization part: 0: low order reinitialization 1: low order reinitialization on all levels is included, to avoid initialization errors! 2: no reinitialisation at all (recommended version!) 3: uses high order reinitialization Keywords: LEVELSET

levelSetSign

MInt* LsCartesianSolver::m_levelSetSign
default = 1

This property specifies if the sign of the level-set function during initialization or reconstruction should be assumed to be as computed e.g. from the stl geometry (1) or inverted (-1).
Possible values are:

  • 1
  • -1

Keywords: LEVELSET, MULTIPLE LEVEL SET FUNCTIONS, MOVING_BOUNDARY

computeSet

MBool* LsCartesianSolver::m_computeSet
default = 1

This property specifies if the level set function should be transported in time or if initialization at startup is sufficient.
Possible values are:

  • 1
  • 0

Keywords: LEVELSET

GFieldInitFromSTL

MInt LsCartesianSolver::m_GFieldInitFromSTL
default = 0

This property triggers if the level-set function(s) are initialized from an .stl-geometry.

  • off - level-set functions are not initialized from .stl-geometries
  • on - level-set functions are initialized from .stl-geometries

Keywords: LEVELSET, LEVEL_SET_FROM_STL

GFieldFromSTLInitCheck

MBool LsCartesianSolver::m_GFieldFromSTLInitCheck
default = 0

This property triggers if an output is generated after the initialization. Careful: produces valgrind errors as only the level-set function is initialized, the other variables which are written out are not.

  • off - no check output is generated
  • on - check output is generated after the initialization

Keywords: LEVELSET, LEVEL_SET_FROM_STL

levelSetWithSTLCorrection

MBool LsCartesianSolver::m_GWithReConstruction
default = 0

This property triggers if the level-set function(s) are corrected using the .stl-geometry

  • off - level-set stl correction is off
  • on - level-set stl correction is on

Keywords: LEVELSET, LEVEL_SET_FROM_STL

buildCollectedLevelSetFunction

MInt LsCartesianSolver::m_buildCollectedLevelSetFunction
default = 0

This property triggers if the level-set functions from m_startSet to m_noSets are collected in set 0
possible values are:

  • off - 0th level-set function is computed normally and not collected
  • on - 0th level-set function is collected based on the other level-set functions

Keywords: LEVELSET, MULTILEVELSET

movingBndryCndId

MInt* LsCartesianSolver::movingBndryCndId
default = 3006

This property specifies the standard moving boundary condition Id -> required for the identification of moving objects. If no other bndryCndIds are specified in bodyBndryCndIds, this Id is assumed to be the Id of all moving bodies (in this case, just one body).
Possible values are:

  • Valid boundary condition ids (see boundary condition ids in the code) -> not necessarily a moving bndry cnd Id, this is only required if you work with the MB part of the code!

Keywords: LEVELSET, MULTIPLE LEVEL SET FUNCTIONS

noEmbeddedBodies

MInt* LsCartesianSolver::noEmbeddedBodies
default = 1

This property specifies the number of bodies released into the flow field.
Possible values are:

  • 1

Keywords: LEVELSET, MULTIPLE LEVEL SET FUNCTIONS

bodyBndryCndIds

MInt* LsCartesianSolver::m_bodyBndryCndIds
default = movingBndryCndId

This property specifies the moving boundary condition Ids -> required for the identification of moving objects. If no bndryCndIds are specified here, the Id specified in movingBndryCndId is assumed to be the Id of all moving bodies (in this case, just one body).
Possible values are:

  • all valid Boundary Condition Ids -> not necessarily a moving bndry cnd Id, this is only required if you work with the MB part of the code!

Keywords: LEVELSET, MULTIPLE LEVEL SET FUNCTIONS, MOVING_BOUNDARY

highOrderDeltaFunction

MBool LsCartesianSolver::m_highOrderDeltaFunction
default = 0

test trigger for high order arclength calculcation which equals the flame surface area
possible values are:

  • 0 - off
  • 1 - on

Keywords: LEVELSET, COMBUSTION, DELTA, FUNCTION, FLAME, SURFACE, AREA

fourthOrderNormalCurvatureComputation

MBool LsCartesianSolver::m_fourthOrderNormalCurvatureComputation
default = 0

test trigger for fourth order normal and curvature computation, should not be used for combustion
possible values are:

  • 0 - off
  • 1 - on

Keywords: LEVELSET, COMBUSTION, FOURTH, ORDER, NORMAL, CURVATURE

curvatureDamp

MBool LsCartesianSolver::m_curvatureDamp
default = 0

temporarly trigger for damping the curvature and heat release
possible values are:

  • 0 - off
  • 1 - on

Keywords: COMBUSTION, LEVELSET, DAMPING, CURVATURE

curvatureDampFactor

MFloat LsCartesianSolver::curvatureDampFactor
default = 0

temporarly trigger for damping the curvature by a factor
possible values are:

  • Non-negative floating point values of the order of 3.0

Keywords: COMBUSTION, LEVELSET, DAMPING, CURVATURE

sharpDamp

MBool LsCartesianSolver::m_sharpDamp
default = 0

temporarly trigger for sharp damping the curvature and heat release
possible values are:

  • 0 - off
  • 1 - on

Keywords: COMBUSTION, LEVELSET, SHARP, DAMPING, CURVATURE

useLocalMarksteinLength

MBool LsCartesianSolver::m_useLocalMarksteinLength
default = 0

temporarly trigger for using the local markstein length when computing the neutral markstein length, (see testcases 1990, 19901)
possible values are:

  • 0 - off
  • 1 - on

Keywords: COMBUSTION, LEVELSET, LOCAL, MARKSTEIN, LENGTH

hyperbolicCurvature

MBool LsCartesianSolver::m_hyperbolicCurvature
default = 0

test trigger for the use of the hyperbolic extension of the curvature itselfs.
possible values are:

  • 0 - off
  • 1 - on

Keywords: LEVELSET, COMBUSTION, CURVATURE, HYPERBOLIC, EXTENSION

gRKMethod

MInt LsCartesianSolver::m_gRKMethod
default = no default value

This property sets the level set runge kutta method
possible values are:

  • 0 -
  • 1 - three step third order TVD Runge-kutta scheme
  • 2 -
  • 3 -
  • 4 -
  • 5 - no runge kutta is solved (fixes the level set function)

Keywords: LEVELSET, RUNGE, KUTTA, SCHEME

levelSetDiscretizationScheme

MString LsCartesianSolver::m_levelSetDiscretizationScheme
default = no default value

This property sets the level set discretization scheme of the level set equation.
possible values are:

  • US1 - first order upwind scheme
  • UC3 - third order upwind scheme
  • UC3_SB - third order upwind scheme
  • UC5 - fifth order upwind scheme
  • UC5_SB - fifth order upwind scheme
  • WENO5 - fifth order upwind scheme
  • WENO5_SB - fifth order upwind scheme
  • UC11 - higher order upwind scheme

Keywords: LEVELSET, DISCRETIZATION, METHOD, EQUATION, RHS

LsRotate

MString LsCartesianSolver::m_LsRotate
default = flase

This property needs to be true for rotating levelset

gBandWidth

MInt LsCartesianSolver::m_gBandWidth
default = no default value

This property sets the level set band width in which the level set equation is solved. The outer band cells are reset to a constant level set function value G and are not taken into account for the level set computation. E.g. gBandWidth = 15 means 15 cells in each direction of the G=0 cells are put into the level set band
possible values are:

  • Any non-negative floating point smaller than the number of time steps timeSteps

Keywords: LEVELSET, BAND, WIDTH

gShadowWidth

MInt LsCartesianSolver::m_gShadowWidth
default = no default value is set

This property sets the level set shadow width in which the level set cells are refined to the maximum G-Cell refinement level maxGCellLevel.
possible values are:

  • Any non-negative floating point number which greater than the level set band width m_gBandWith

Keywords: LEVELSET, SHADOW, WIDTH, MAXIMUM, REFINEMENT, LEVEL

gShadowWidthRans

MInt LsCartesianSolver::m_gShadowWidthRans
default = no default value is set

This property sets the level set shadow width for levelSetRans in which the level set cells are refined to the maximum G-Cell refinement level maxGCellLevel.
possible values are:

  • Any non-negative floating point number which greater than the level set band width m_gBandWith

Keywords: LEVELSET, SHADOW, WIDTH, MAXIMUM, REFINEMENT, LEVEL

gInnerBound

MInt LsCartesianSolver::m_gInnerBound
default = no default value

.
possible values are:

Keywords: LEVELSET, METHOD

maxGCellLevel

MInt LsCartesianSolver::m_maxGCellLevel
default = no default value

This property sets the maximum level set refinement level of the level set grid. The grid width is calculated via \( m_gCellDistance = c_cellLengthAtLevel(m_maxGCellLevel) \).
The maximum level set refinement lebel is limited by the maximum refinement level maxRfnmtLevel
possible values are:

  • Any non-negative integer value smaller than the maximum refinement level

Keywords: LEVELSET, MAXIMUM, GCELL, REFINEMENT, LEVEL

computeExtVel

MInt LsCartesianSolver::m_computeExtVel
default = 1

This property sets the level set method to compute the corrected burning velocity and the extension velocity of the level set cells.
possible values are:

  • 1 - computes extension velocity (curvature effect controlled by the Markstein length), marksteinLength
  • 10 - computes extension velocity (curvature effect controlled by the Markstein length), marksteinLength
  • 2 - computes extension velocity for the combustion solver FVCombstnSolver (curvature effect controlled by the Markstein length), marksteinLength
  • 5 - computes extension velocity for the G-Equation Progress Variable Model FVGequPvSolver (cuvrature effect controlled by the molecular diffusivity)
  • 6 - computes extension velocity for the G-Equation Progress Variable Model FVGequPvSolver (no curvature effect)
  • 60 - computes extension velocity for the G-Equation Progress Variable Model FVGequPvSolver (no curvature effect)
  • 7 - computes extension velocity for the G-Equation Progress Variable Model FVGequPvSolver (curvature effect controlled by the Markstein length), marksteinLength
  • 70 - computes extension velocity for the G-Equation Progress Variable Model FVGequPvSolver (curvature effect controlled by the Markstein length + stretch effects), marksteinLength
  • 71 - computes extension velocity for the G-Equation Progress Variable Model FVGequPvSolver (curvature effect controlled by the Markstein length for symmetric flames), marksteinLength
  • 8 - computes extension velocity for the G-Equation Progress Variable Model FVGequPvSolver (computation of neutral Markstein length and growth rate of perturbed flame surfaces), marksteinLength
  • 81 - computes extension velocity for the G-Equation Progress Variable Model FVGequPvSolver (stability test with use of neutral Markstein length (enforced version of case 8), neutral Markstein length is computed for every time step), marksteinLength

Keywords: LEVELSET, COMBUSTION, EXTENSION, VELOCITY, NEUTRAL, MARKSTEIN, LENGTH

smoothExtVel

MInt LsCartesianSolver::m_smoothExtVel
default = 1

This property triggers the smoothed computation of the corrected burning velocity in the band layer which is done by different computations, see computeExtVel
possible values are:

  • 0 - computation of extension velocity only at the G zero front cells
  • 1 - computation of smoothed extension velocity for the band layer cells, see computeExtVel

Keywords: LEVELSET, SMOOTHED, EXTENSION, VELOCITY, BAND, LAYER, G0, ZERO, FRONT, CELLS

extVelIterations

MInt LsCartesianSolver::m_extVelIterations
default = 10000

controls the maximum iterations to reach the convergence criteria of the hyperbolic extension, extVelCFL, extVelConvergence
possible values are:

  • Non-negative integer values of the order of 1000

Keywords: LEVELSET, COMBUSTION, CURVATURE, HYPERBOLIC, EXTENSION

extVelConvergence

MInt LsCartesianSolver::m_extVelConvergence
no default valaue used

controls the convergence of the hyperbolic extension, extVelCFL, extVelConvergence
possible values are:

  • Non-negative floating point values of the order of 0.1

Keywords: LEVELSET, COMBUSTION, COMBUSTION, HYPERBOLIC, EXTENSION

extVelCFL

MFloat LsCartesianSolver::m_extVelCFL
default = 0.5

controls the CFL number of the pseudo time discretization of the hyperbolic extension, extVelConvergence, extVelIterations
possible values are:

  • Non-negative floating point values of the order of 0.1

Keywords: LEVELSET, COMBUSTION, COMBUSTION, HYPERBOLIC, EXTENSION

reinitMethod

MString LsCartesianSolver::m_reinitMethod
default = no default value

This property sets the reinitialization method, for more information see Dissertation D. Hartmann
possible values are:

  • CR1
  • CR2
  • HCR1
  • HCR2
  • HCR2_LIMITED
  • HCR2_FULLREINIT
  • SUS5CR
  • CR2PLUS
  • RSU
  • DL1
  • DL2
  • SUS_1
  • SUS_1PLUS
  • SUS_2
  • SUS_WENO5
  • SUS_WENO5PLUS

Keywords: LEVELSET, REINITIALIZATION, METHOD

gapReinitMethod

MString LsCartesianSolver::m_gapReinitMethod
default = m_reinitMethod

This property sets the reinitialization method for the gap reinitialization for multiple level-set functions, for more information see Dissertation D. Hartmann
possible values are:

  • CR1
  • CR2
  • HCR1
  • HCR2
  • HCR2_LIMITED
  • HCR2_FULLREINIT
  • SUS5CR
  • CR2PLUS
  • RSU
  • DL1
  • DL2
  • SUS_1
  • SUS_1PLUS
  • SUS_2
  • SUS_WENO5
  • SUS_WENO5PLUS

Keywords: LEVELSET, REINITIALIZATION, METHOD

gapReinitMethod

MBool LsCartesianSolver::m_closeGaps
default = false

This property triggers if the reinitialization based closure procedure for narrow gaps between different bodies is called for the collected level-set function

  • 0
  • 1

Keywords: LEVELSET, MULTILEVELSET

determineG0Mode

MBool LsCartesianSolver::m_determineG0CellsMode
default = false

Triggers if G0 cells of combined level-set function are determined from phi^0 (mode=0) or if they are inherited from the individual level-set functions (mode=1)

  • 0
  • 1

Keywords: LEVELSET, MULTILEVELSET

noGapRegions

MFloat noGapRegions
Default: -1
Set how many regions can be gap regions.
Possible values are:

  • Any positive integer number.

Keywords: DISCONTINUOUS_GALERKIN, APE

G0regionId

MInt LsCartesianSolver::m_G0regionId
default = -1
defines a region for which the combined levelSet is used to smooth the overlapping of different bodies. However no gapClosing where the combined levelSet is additionally modified is applied! NOTE: possibly a sponge needs to a applied in this region if any of the bodies is moving!

  • -1: not used (1)
  • positiv initger: regionId for which this is used!

Keywords: MOVING BOUNDARY, GAP, FINITE_VOLUME, LEVELSET

gapInitMethod

MInt LsCartesianSolver::m_gapInitMethod
default = 2
Possible values are:

  • 0: all gap Cells are reseted and initialized (older version from Claudia) (1)
  • 1: only new arising gap Cells are reseted! (2)
  • 2: only new arising gap Cells are reseted and mew epsilons are used

Keywords: MOVING BOUNDARY, GAP, FINITE_VOLUME

forceNoGaps

MInt LsCartesianSolver::m_forceNoGaps
default = 0
triggers the surpression of gap cells

  • 0: do not surpress gap cells (1)
  • 1: surpress gap cells for multi-valve engine (2)
  • 2: surpress gap cells for single-valve engine

Keywords: MOVING BOUNDARY, GAP, LS

deltaMin

MBool LsCartesianSolver::m_gapDeltaMin
default = 1

  • 0
  • 1

Keywords: LEVELSET, MULTILEVELSET

gReinitIterations

MInt LsCartesianSolver::m_gReinitIterations
default = no default value

This property triggers the maximum number of reinitialization steps to reach the convergence criterion reinitConvergence
possible values are:

  • Any non-negative integer value

Keywords: LEVELSET, REINITIALIZATION, ITERATION, STEPS

minReinitializationSteps

MInt LsCartesianSolver::m_minReinitializationSteps
default = m_gReinitIterations

This property triggers the minimum number of reinitialization steps to reach the convergence criterion reinitConvergence
possible values are:

  • Any non-negative integer value

Keywords: LEVELSET, MINIMUM REINITIALIZATION, ITERATION, STEPS

maintenanceIterations

MInt LsCartesianSolver::m_maintenanceIterations
default = gReinitIterations

This property triggers the maximum number of iterations for maintaining the outer band layers
possible values are:

  • Any non-negative integer value

Keywords: LEVELSET, MAINTENANCE, ITERATION, STEPS

guaranteeReinit

MBool LsCartesianSolver::m_guaranteeReinit
default = 0

this property guarantees the reinitialization at every time step
possible values are:

  • 0 - off
  • 1 - on

Keywords: LEVELSET, GUARANTEE, REINITIALIZATION

reinitCFL

MInt LsCartesianSolver::m_reinitCFL
default = 1

This property sets the reinitialization cfl number which is used to calculate the pseudo time step for the reinitialization \( \Delta t_{pseudo} = CFL_{reinit} * \Delta x \), with \( \Delta x \) being the smallest G cell width, see maxGCellLevel
possible values are:

  • Any non-negative floating point number smaller than 1

Keywords: LEVELSET, REINITIALIZATION, CFL, PSEUDO, TIME, STEP

intermediateReinitIterations

MInt LsCartesianSolver::m_intermediateReinitIterations
default = 1

This property triggers the maximum number of iterations for reinitializing the level set function at every runge kutta step. Only used for the runge kutta method 3, gRKMethod.
possible values are:

  • Any non-negative integer value

Keywords: LEVELSET, REINITIALIZATION, ITERATION, STEPS, RUNGE, KUTTA

reinitConvergence

MFloat LsCartesianSolver::m_reinitConvergence
default = no default value

This property sets the convergence criterion of the reinitialization. The criterion should be of the order of 10^{-3} - 10^{-10}, see also gReinitIterations
possible values are:

  • Any non-negative floating point number smaller than 1

Keywords: LEVELSET, REINITIALIZATION, CONVERGENCE, CRITERION, ITERATION, STEPS

reinitThreshold

MFloat LsCartesianSolver::m_reinitThreshold
default = 0

This property sets the reinitialization threshold. The important parameter for the reinitialization is the maintenance of the property \( \left| \nabla G \right| = 1\).
The deviation to this property is the so called reinitialization threshold calculated for all level set front cells (G0 cells) via
\( \Delta_{thresh} = \left( \left| \nabla G \right| - 1 \right)^{2} \). If the maximum threshold and the averaged threshold, see reinitThresholdAvg,
is reached the level set will be reinitialized otherwise maintained.
possible values are:

  • Any non-negative floating point number

Keywords: LEVELSET, REINITIALIZATION, THRESHOLD, G0, CELLS

reinitThresholdAvg

MFloat LsCartesianSolver::m_reinitThresholdAvg
default = 0

This property sets the averaged reinitialization threshold. The important parameter for the reinitialization is the maintenance of the property \( \left| \nabla G \right| = 1\).
The deviation to this property is the so called reinitialization threshold calculated for all level set front cells (G0 cells), see reinitThreshold, and the averaged threshold
\( \overline{\Delta_{threshAvg}} = \overline{\left( \left| \nabla G \right| - 1 \right)^{2}} \) of the G0 cells.
If both criterions are reached the level set will be reinitialized otherwise maintained.
possible values are:

  • Any non-negative floating point number

Keywords: LEVELSET, REINITIALIZATION, THRESHOLD, G0, CELLS

omegaReinit

MInt LsCartesianSolver::m_omegaReinit
default = 1

This property is used for the constrained reinitialization methods Cr1, CR2, CR2PLUS and RSU, see reinitMethod.
possible values are:

  • Any non-negative floating point number

Keywords: LEVELSET, REINITIALIZATION, OMEGA

relaxationFactor

MInt LsCartesianSolver::m_relaxationFactor
default = 0.5

This property forces the reinitialization.
possible values are:

  • Any non-negative floating point number smaller than one

Keywords: LEVELSET, REINITIALIZATION, RELAXATION, FACTOR, CELLS, TO, CORRECT

levelSetTestCase

MInt LsCartesianSolver::m_levelSetTestCase
default = 0

This property triggers the computation of the propagation speed for a few test cases and triggers the calculation of a few level set ouput informations, see file levelSetData2 in LsCartesianSolver::writeLevelSet().
possible values are:

  • 0 - no test case chosen
  • 1 - ?
  • 200 - diagonal convection
  • 201 - merge problem
  • 202 - oscillating circle problem, sine wave around a cirlce
  • 203 - Zalesak's problem
  • 204 - Enright's test
  • 205 - Harmonic forcing problem (Preetham07)
  • 206 - simple expansion/compression?
  • 207 - periodic array of 5 vortices in a 3x1 domain
  • 208 - rotation
  • 1001 - ?
  • 1002 - ?

3D test cases:

  • 300 - convection
  • 301,30100 - merge/expand
  • 302 - merge/shrink
  • 303 - Oscillating circle problem
  • 305 - Enright test
  • 306 - Cold vortex-flame interaction
  • 366 - Oscillating sphere
  • 302 - Pancake sphere
  • 3005 - Enright's test by STL
  • 3006 - Zalesak Disk by STL, refered to Dupont's and Back compensation
  • 3009 - General rigid body motion template

Keywords: LEVELSET, CASES, PROPAGATION, SPEED, OUTPUT

levelSetBoundaryCondition

MInt LsCartesianSolver::m_levelSetBoundaryCondition
default = m_levelSetTestCase

This property triggers the use of special boundary conditions for each specific level set case, see also levelSetTestCase
possible values are:

  • 1 - 2D bunsen flame
  • 1949 - ?
  • 1957 - ?
  • 1966,1967,1968 - ?
  • 17511,175110 - forced response
  • 17512 - forced response bunsen flame
  • 17513,17514,17515,17516,17518 - forced response bunsen flame -> fixes the level set function on the flame tube edges, radiusFlameTube, xOffsetFlameTube(2) for one or two flames, twoFlames

  • 17517 - forced response symmetric bunsen flame -> fixes the level set function on the flame tube edges, radiusFlameTube, xOffsetFlameTube

Keywords: LEVELSET, SMOOTHED, EXTENSION, VELOCITY, BAND, LAYER, G0, ZERO, FRONT, CELLS

levelSetBC

MString FvCartesianSolverPar::m_levelSetBC
default = "SYMMETRIC"

This property triggers the use of simplifications in the GField-calculations for special testcases and setups, possible values are:

  • SYMMETRIC
  • PERIODIC
  • NONE

If the Levelset-Field is neither symmetric nor periodic, the NONE specification has to be selected! Keywords: LEVELSET BAND, LAYER, G0, ZERO, FRONT, CELLS

noHaloLayers

MInt LsCartesianSolver::m_noHaloLayers
default = m_noHaloLayers

This property sets the number of halo layers for the level set grid.
possible values are:

  • Any non-negative integer number greater than one for parallel computation

Keywords: LEVELSET, NUMBER, HALO, LAYERS, PARALLEL, COMPUTATION

reinitInterval

MInt LsCartesianSolver::m_reinitInterval
default = 1

This property triggers the reinitialization interval, but it doesn't ensures the reinitialization, this is depend on wether the thresholds, reinitThreshold, reinitThresholdAvg
or could be ensured by the property guaranteeReinit.
possible values are:

  • Any non-negative integer number smaller than the maximum number of time steps timeSteps

Keywords: LEVELSET, REINITIALIZATION, INTERVAL

maintainOuterBandLayers

MBool LsCartesianSolver::m_maintainOuterBandLayers
default = 0

This property guarantees that the outer band layers are maintained, see also reinitInterval and guaranteeReinit.
possible values are:

  • 0 - off
  • 1 - on

Keywords: LEVELSET, MAINTENANCE, GUARANTEE, OUTER, BAND, LAYERS

writeReinitializationStatistics

MInt LsCartesianSolver::m_writeReinitializationStatistics
default = 0

This property triggers the reinitialization statistics. If the statistics are turned on you get the ASCII files deviation and avgGradientAfter containing the statistics for the G0 cells (level set front cells):
deviation:

avgGradientAfter:

  • globalTimeStep
  • averaged deviation (here called: averaged gradient), see reinitThresholdAvg
  • variance \( \sigma = \overline{\Delta_{threshAvg}}^{2} \)
  • mean level set function gradient \( \overline{\left| \nabla G \right| }\)
  • maximum level set function gradient \( \left| \nabla G \right|_{MAX}\)
  • minimum level set function gradient \( \left| \nabla G \right|_{MIN}\)

The property is only meaningful for the renitialization methods, reinitMethod, HCR1,HCR2,HCR2_LIMITED,HCR2_FULLREINIT. The statistics should be included also for the other reinitialization methods.
possible values are:

  • 0 - off
  • 1 - on

Keywords: LEVELSET, REINITIALIZATION, STATISTICS, DEVIATION, AVGGRADIENTAFTER, G0, CELLS, GRADIENT

interpolateFlowFieldToFlameFront

MBool LsCartesianSolver::m_interpolateFlowFieldToFlameFront
default = 0

This property triggers the transfer of the velocity v from the flow to the G-grid and interpolates the flow field to the flame front, see computeExtVel (cases 5..81)
possible values are:

  • off - transfer is done but without interpolation (possible for all test cases)
  • on - transfer is done with interpolations (probably not possible for all test cases)

Keywords: LEVELSET, TRANSFER, FLOW, FIELD, INTERPOLATION

writeOutAllLevelSetFunctions

MBool LsCartesianSolver::m_writeOutAllLevelSetFunctions
default = 0

This property triggers if all level-set functions are written out (.vtk-file), or just the zeroth level-set function
possible values are:

  • off - output of the zeroth level-set function only
  • on - output of all level-set functions (file may become large!)

Keywords: LEVELSET, MULTILEVELSET

writeOutAllExtensionVelocities

MBool LsCartesianSolver::m_writeOutAllExtensionVelocities
default = 0

This property triggers if all level-set extension velocities are written out (.vtk-file), or just respect to the zeroth level-set function
possible values are:

  • off - output of the extension velocities with respecot to the zeroth level-set function only
  • on - output of all level-set extension velocities (file may become large!)

Keywords: LEVELSET, MULTILEVELSET

writeOutAllCurvatures

MBool LsCartesianSolver::m_writeOutAllCurvatures
default = 0

This property triggers if all level-set function curvatures are written out (.vtk-file), or just the zeroth level-set function curvature
possible values are:

  • off - output of the zeroth level-set function curvature only
  • on - output of all level-set function curvatures (file may become large!)

Keywords: LEVELSET, MULTILEVELSET

writeOutAllCorrectedBurningVelocity

MBool LsCartesianSolver::m_writeOutAllCorrectedBurningVelocity
default = 0

This property triggers if all level-set function corrected burning velocities are written out (.vtk-file), or just the zeroth level-set function corrected burning velocity
possible values are:

  • off - output of the zeroth level-set function corrected burning velocity only
  • on - output of all level-set function corrected burning velocities (file may become large!)

Keywords: LEVELSET, MULTILEVELSET

writeOutAllFlameSpeeds

MBool LsCartesianSolver::m_writeOutAllFlameSpeeds
default = 0

This property triggers if all level-set function flame speeds are written out (.vtk-file), or just the zeroth level-set function flame speed
possible values are:

  • off - output of the zeroth level-set function flame speed only
  • on - output of all level-set function flame speeds (file may become large!)

Keywords: LEVELSET, MULTILEVELSET

writeOutAllNormalVectors

MBool LsCartesianSolver::m_writeOutAllNormalVectors
default = 0

This property triggers if all level-set function normal vectors are written out (.vtk-file)
possible values are:

  • off - output of the normal vectors
  • on - output of the normal vectors (file may become large!)

Keywords: LEVELSET, MULTILEVELSET

engineSetup

MBool LsCartesianSolver::m_engineSetup
Triggers specific stuff for TINA or other engine applications default = false

Keywords: LEVELSET, ADAPTATION

refineDiagonals

MBool LsCartesianSolver::m_refineDiagonals
default = true

Determines whether the diagonal cells for the levelSet band should be refined as well!

  • 1.0 +- eps

Keywords: LEVELSET, ADAPTATION

bodyIdOutput

MBool LsCartesianSolver::m_bodyIdOutput
default = true

Determines whether the additionale bodyIdOutput output should be written! Default is true for semiLagrange, but the output is not necessary for a restart and should be set to false for large applications, to reduce the amount of data! Keywords: LEVELSET, SEMI-LAGRANGE

LsGeometryChange

MBool LsCartesianSolver::m_geometryChange
default = false

Used to trigger an adaptation right after the restart, which can be used to change the geometry within the ls-BandWidth. Useful for setup changes. Keywords: LEVELSET, SEMI-LAGRANGE, GEOMETRY

weightBaseCell

MFloat LS::m_weightBaseCell
default = 1.0

Weight applied for any level-set cell during static weight computation for domain decomposition during balance, good value could be 0.01 Keywords: LS, WEIGHTING, BALANCE

weightLeafCell

MFloat LS::m_weightLeafCell
default = 1.0

Weight applied for any level-set leaf-cell during static weight computation for domain decomposition during balance, good value could be 0.05 Keywords: LS, WEIGHTING, BALANCE

weightBandCell

MFloat LS::m_weightBandCell
default = 1.0

Weight applied for any level-set leaf-cell band cell during static weight computation for domain decomposition during balance, good value could be 0.1 Keywords: LS, WEIGHTING, BALANCE

weightMulitSolverFactor

MFloat LS::m_weightMulitSolverFactor
default = 1.0

Mutli-solver weight factor applied to all level-set cell weights for static weight computation for domain decomposition during balance. 1.0 for single solver application , otherwise setup dependent. Keywords: LS, WEIGHTING, BALANCE

limitWeights

MBool LS::m_limitWeights
default = false

Limit weight of level-set cells by a factor of the largest weight, to ensure a more even distribution of solver memory across ranks. Keywords: LS, WEIGHTING, BALANCE

weightLevelSet

MBool LS::m_weightLevelSet
default = true

Triggers if level-set cells should be considered as DLB weights. Can be deactivated if the level-set is not updated each time-step can no overhead exists i.e., for static boundaries. Keywords: LS, WEIGHTING, BALANCE