Species Transport

SpeciesTransport solves the convection-diffusion-reaction transport of each chemical species in a reactive flow. It picks a convective scheme (upwind, central, or van Leer), manages an optional burnt-zone bounding box and a fuel-zone bounding box, and ties into the reaction rate supplied by ThermoChemistry.

Key Classes and Concepts

• SpeciesTransport : public OPObject: per-species transport equation and boundary handling for the reactive-flow sub-tree.

Usage

Input

Defined in the @SpeciesTransport block.

@SpeciesTransport

$BCOrder            Boundary-condition order                 : 0
$TempBurntGas       Temperature of the burnt gas             : 0.0
$Species            Enable multi-species transport           : No
$FuelIndex          Index of the fuel species                : 0

# Advection scheme (pick one) {#advection-scheme-pick-one-2}
$AdvectionUpwind    Upwind convective scheme                 : No
$AdvectionCentral   Central-difference convective scheme     : No
$AdvectionVanLeer   van Leer limited scheme                  : Yes

# Burnt / fuel bounding boxes (physical coordinates; -10 means unused) {#burnt-fuel-bounding-boxes-physical-coordinates--10-means-unused}
$X0BurntZone        Burnt-zone lower X                       : -10
$XNBurntZone        Burnt-zone upper X                       : -10
$Y0BurntZone        Burnt-zone lower Y                       : -10
$YNBurntZone        Burnt-zone upper Y                       : -10
$Z0BurntZone        Burnt-zone lower Z                       : -10
$ZNBurntZone        Burnt-zone upper Z                       : -10

$X0FuelZone         Fuel-zone lower X                        : -10
$XNFuelZone         Fuel-zone upper X                        : -10
$Y0FuelZone         Fuel-zone lower Y                        : -10
$YNFuelZone         Fuel-zone upper Y                        : -10
$Z0FuelZone         Fuel-zone lower Z                        : -10
$ZNFuelZone         Fuel-zone upper Z                        : -10

$IF_REACTION        Enable reaction source term              : No
$nCoeffs            Number of polynomial coefficients        : 15

The defaults for zone boundaries are -10 which the module treats as "unused".

Output

Per-species concentration fields through the standard VTK pipeline.

Example

#include "ReactiveFlows/FlowMixture.h"
#include "ReactiveFlows/SpeciesTransport.h"

FlowMixture      Mix(OPSettings, InputFile);
SpeciesTransport ST(OPSettings, InputFile);

for(RTC.tStep = RTC.tStart; RTC.tStep <= RTC.nSteps; RTC.IncrementTimeStep())
{
    // ... mixture + flow update ...
    ST.Solve(Phi, Mix, BC, RTC.dt);
}

Dependencies

• FlowMixture — global mixture state.
• ThermoChemistry — reaction rates.
• EnergyTransport — temperature field.
• BoundaryConditions.