Temperature equation

IncompressibleNavierStokes.jl supports adding a temperature equation, which is coupled back to the momentum equation through a gravity term [9].

To enable the temperature equation, you need to set the temperature keyword in setup:

setup = Setup(;
    kwargs...,
    temperature = temperature_equation(; kwargs...),
)

where temperature_equation can be configured as follows:

# IncompressibleNavierStokes.temperature_equation — Function.

temperature_equation(
;
    Pr,
    Ra,
    Ge,
    dodissipation,
    boundary_conditions,
    gdir,
    nondim_type
)

Create temperature equation setup (stored in a named tuple).

The equation is parameterized by three dimensionless numbers (Prandtl number, Rayleigh number, and Gebhart number), and requires separate boundary conditions for the temperature field. The gdir keyword specifies the direction gravity, while non_dim_type specifies the type of non-dimensionalization.

source