SPH Adaptive Time-Step Providers Reference

Convective CFL, Gravity CFL, Surface Tension CFL, and Viscous CFL time-step providers limit the time-step size to meet the specified target criteria on Courant number conditions. You use these time-step providers to simulate free surface flows using Smoothed Particle Hydrodynamics (SPH).

The SPH method is not fully implicit. Therefore, the time step is restricted by stability conditions. Adaptive time-stepping is used to determine a time-step size $Δ t$ such that the various Courant number (CFL) conditions are met within the computational domain.

Convective CFL: The Convective CFL time-step provider sets the time-step based on the CFL number, the convective time scale. The time-step is adjusted such that the maximum CFL number of all particles is lower than the specified Target Convective CFL Number. The convective CFL number is defined as:
图 1. EQUATION_DISPLAY

${C F L}_{c o n v e c t i v e, i} = Δ t \frac{v_{i}}{Δ x_{i}} < {C F L}_{c o n v e c t i v e, t a r}$

(102)

where $v_{i}$ is the magnitude of the velocity of particle i and $Δ x_{i} = \sqrt[3]{V_{i}}$ is the size of the cube with the same volume $V$ as the particle i.

Gravity CFL

The Gravity CFL time-step provider sets the time-step based on a CFL number that takes into account the resultant of body forces $f_{b, i}$ . The time-step is adjusted such that the maximum gravity CFL number of all particles is lower than the specified Target Gravity CFL Number. The CFL number based on body forces is defined as:

图 2. EQUATION_DISPLAY

{C F L}_{g r a v i t y, i} = \sqrt{\frac{f_{b, i}}{ρ_{0} Δ x_{i}}} < {C F L}_{g r a v i t y, t a r}

(103)

where $ρ_{0}$ is the fluid density. This time-step provider is available only when gravity is selected in the simulation.

Viscous CFL

In viscous flow simulations, the Viscous CFL time-step provider sets the time-step based on a CFL number that takes into account the fluid viscosity. The time-step is adjusted such that the maximum viscous CFL number is lower than the specified Target Viscous CFL Number. The CFL number based on fluid viscosity is defined as:

图 3. EQUATION_DISPLAY

{C F L}_{v i s o c u s, i} = Δ t \frac{μ_{i}}{ρ_{0} {Δ x_{i}}^{2}} < {C F L}_{v i s o c u s, t a r}

(104)

where $μ_{i}$ is the dynamic viscosity of the particle i.

For SPH simulations with the explicit velocity solver, you are advised to use all the SPH compatible time-step providers (Convective CFL, Gravity CFL, and Viscous CFL). The Viscous CFL time-step provider is not required for simulations that use the implicit velocity solver.

表 1. SPH Adaptive Time-Step Providers Reference
Provided By	Right-click Continua > Physics 1 > Models > Adaptive Time-Step > Time-Step Providers and select New, followed by: Convective CFL Gravity CFL Viscous CFL
Example Node Path	Continua > Physics 1 > Models > Adaptive Time-Step > Time-Step Providers > Convective CFL
	Continua > Physics 1 > Models > Adaptive Time-Step > Time-Step Providers > Gravity CFL
	Continua > Physics 1 > Models > Adaptive Time-Step > Time-Step Providers > Viscous CFL
Requires	Time: Implicit Unsteady. Material: Multiphase. Multiphase Models: Smoothed-Particle Hydrodynamics (SPH).
Properties	Enable/Disable See Adaptive Time-Step Provider Common Properties. Target Convective CFL Number and Include Boundaries See Convective CFL Properties. Target Gravity CFL Number See Gravity CFL Properties. Target Viscous CFL Number See Viscous CFL Properties.
Activates	Field Functions	None.

Convective CFL Properties

Target Convective CFL Number: Specifies ${C F L}_{c o n v e c t i v e, t a r}$ in Eqn. (102). The default value is 0.72.

Include Boundaries: When activated, includes the boundary velocities in the calculation of the time-step for the convective CFL condition.

Gravity CFL Properties

Target Gravity CFL Number: Specifies ${C F L}_{g r a v i t y, t a r}$ in Eqn. (103). The default value is 0.073.

Viscous CFL Properties

Target Viscous CFL Number: Specifies ${C F L}_{v i s c o u s, t a r}$ in Eqn. (104). The default value is 0.01.