Solves the oxidant diffusion step of the Suvorov et al. (10.1007/s10825-006-0003-z) oxidation model on the Cartesian grid carrying two level sets. More...

#include <lsOxidationDiffusion.hpp>

Inheritance diagram for viennals::OxidationDiffusion< T, D >:

Classes
struct	ReactionBoundarySample
	Sub-grid accurate sample of the reaction boundary crossing closest to a given grid node. Carries the interpolated concentration at the boundary point, plus the axis and offset of the edge on which the crossing lies so that callers can reconstruct the boundary-point position and normal. More...
struct	FrontVelocity
	Si front speed (k C_i / (N alpha), >= 0) and normal (level-set gradient direction) at a point on the front: the area-weighted average over the face crossings within two cells (Gaussian weights, 0.75 cells wide). Not found before the first solve or away from the front. More...

Public Types
using	StressLookup = std::unordered_map<std::size_t, std::array<T, 9>>
using	PressureLookup = std::unordered_map<std::size_t, T>

Public Member Functions
	OxidationDiffusion ()=default
	OxidationDiffusion (SmartPointer< Domain< T, D > > passedReactionInterface, SmartPointer< Domain< T, D > > passedAmbientInterface, OxidationParameters passedParameters={})
void	markGeometryChanged ()
	Call after any in-place modification of the level sets (e.g. after ls::Advect) so that the next apply() rebuilds the Cartesian node grid.
void	setReactionInterface (SmartPointer< Domain< T, D > > passedInterface)
void	setAmbientInterface (SmartPointer< Domain< T, D > > passedInterface)
void	setMaskInterface (SmartPointer< Domain< T, D > > passedInterface, int passedMaskSign=1)
void	clearMaskInterface ()
void	setParameters (OxidationParameters passedParameters)
OxidationParameters	getParameters () const
T	getEffectiveReactionRate (const Vec3D< T > &coordinate) const
void	clearPressureField ()
void	setStress (const IndexType &index, const std::array< T, 9 > &stress)
	Oxide Cauchy stress (Pa, tension positive, row-major 3x3) at a node, for the normal-stress reaction law.
void	getFeedback (const IndexType &index, T &pressure, std::array< T, 9 > &stress) const
	The pressure and stress the last solve used at a node (reference pressure / zero stress where nothing was fed back).
const PressureLookup &	getPressureCache () const
	Pressure and stress feedback carried from one solve to the next, like the concentration cache: lsOxidation builds a new OxidationDiffusion for every coupled solve, and without the carry-over every solve began its stress coupling from zero stress (~7 coupling iterations per solve on a Kao cylinder, residual 1.0 at the first).
const StressLookup &	getStressCache () const
void	setFeedbackCache (PressureLookup pressures, StressLookup stresses)
const std::vector< T > &	getInterfaceNormalStress () const
	The interface normal stress the reaction law used in the last solve (per node; NaN where the law falls back to the unstressed rate).
std::vector< T >	interfaceNormalStressFor (const StressLookup &lookup) const
	The interface normal stress the reaction law would use with lookup: the median sample at every node with a REACTION face, smoothed along the interface (see updateInterfaceNormalStress for why). Non-mutating, so the coupling loop can measure convergence without changing the state.
void	setStress (const Vec3D< T > &coordinate, const std::array< T, 9 > &stress)
void	setPressure (const IndexType &index, T pressure)
void	setPressure (const Vec3D< T > &coordinate, T pressure)
void	setOxideSigns (int passedReactionSign, int passedAmbientSign)
	Set signs defining the oxide band. A node is inside oxide if reactionSign * reactionPhi >= 0 and ambientSign * ambientPhi >= 0.
void	setSolveBounds (const IndexType &passedMinIndex, const IndexType &passedMaxIndex)
	Restrict the dense Cartesian diffusion solve to a finite index box. This is useful for level sets with infinite boundary conditions.
void	clearSolveBounds ()
void	apply ()
T	getScalarVelocity (const Vec3D< T > &coordinate, int material, const Vec3D< T > &normalVector, unsigned long) final
	Should return a scalar value for the velocity at coordinate for a point of material with the given normalVector.
FrontVelocity	frontVelocityAt (const Vec3D< T > &coordinate) const
T	getDissipationAlpha (int, int material, const Vec3D< T > &) final
	If lsLocalLaxFriedrichsAnalytical is used as the spatial discretization scheme, this is called to provide the analytical solution for the alpha values, needed for numerical stability.
T	getConcentration (const Vec3D< T > &coordinate) const
T	getConcentration (const IndexType &index) const
T	getReactionBoundaryConcentration (const Vec3D< T > &coordinate) const
T	getReactionBoundaryConcentration (const IndexType &index) const
unsigned	getIterations () const
T	getResidual () const
T	getNormalizedResidual () const
bool	lastSolveConverged () const
std::size_t	getNumberOfSolutionNodes () const
bool	hasFiniteConcentrationField () const
const std::unordered_map< std::size_t, T > &	getConcentrationCache () const
void	setConcentrationCache (std::unordered_map< std::size_t, T > cache)
void	setGpuMode (GpuMode mode)
	Set the GPU solver selection mode (see GpuMode). On CPU-only builds (VIENNALS_GPU_BICGSTAB not defined) Gpu is an error and Auto the CPU.
void	markSolved ()
	Mark the current solution as valid without re-solving. Call this before any parallel advection (lsAdvect) that uses this field as a velocity source to prevent concurrent apply() calls inside getScalarVelocity().
void	writeConcentrationToLevelSet ()
	Write per-node concentration into ambientInterface->getPointData() so that lsInterior + lsAdvect can carry it across timestep boundaries. Safe to call regardless of the solved flag; uses the most recent nodes.
void	writePressureToLevelSet ()
	Write per-node pressure into ambientInterface->getPointData() so that it survives advection and can warm-start the coupling loop next step.
void	writePersistentFields ()
	Convenience wrapper: persist both concentration and pressure in one call.
ReactionBoundarySample	getReactionBoundarySample (const Vec3D< T > &coordinate) const
	Return the reaction boundary sample for the grid node nearest to coordinate. Used by the deformation solver to get both the sub-grid concentration and the crossing edge so it can compute a sub-grid normal.
T	getScalarVelocityFromSample (const ReactionBoundarySample &sample) const
	Absolute scalar velocity derived from an already-computed boundary sample, avoiding a second call to reactionBoundarySample inside getScalarVelocity.
virtual Vec3D< T >	getVectorVelocity (const Vec3D< T > &, int, const Vec3D< T > &, unsigned long)
	Like getScalarVelocity, but returns a velocity value for each cartesian direction.

Static Public Member Functions
template<class... Args>
static auto	New (Args &&...args)

Protected Member Functions
template<class GridType>
void	setMirrorSides (const GridType &grid)
int	stencilOffset (const IndexType &index, unsigned direction, int offset) const
bool	crosses (T a, T b) const
T	valueAt (ConstSparseIterator &it, const IndexType &index) const
template<class OnGrid>
Vec3D< T >	definedLevelSetGradient (ConstSparseIterator &it, const IndexType &index, OnGrid onGrid) const
bool	inBounds (const IndexType &index) const
void	initNodeLookup ()
std::size_t	lookupNode (const IndexType &index) const
std::size_t	linearIndex (const IndexType &index) const
bool	increment (IndexType &index) const
std::size_t	findNearbyNode (const IndexType &index) const
template<class Accept>
std::vector< std::size_t >	findNearbyNodes (const IndexType &index, Accept accept) const
bool	initializeGridFromInterfaces (SmartPointer< Domain< T, D > > reactionInterface, SmartPointer< Domain< T, D > > ambientInterface, SmartPointer< Domain< T, D > > maskInterface, bool useRequestedBounds, const IndexType &requestedMinIndex, const IndexType &requestedMaxIndex, std::size_t maxGridPoints, const std::string &solverName)
bool	initializeGridFromMask (SmartPointer< Domain< T, D > > maskInterface, bool useRequestedBounds, const IndexType &requestedMinIndex, const IndexType &requestedMaxIndex, std::size_t maxGridPoints, const std::string &solverName)

Protected Attributes
std::vector< std::size_t >	nodeLookupFlat
IndexType	minIndex {}
IndexType	maxIndex {}
std::array< std::size_t, D >	extents {}
std::array< std::size_t, D >	strides {}
T	gridDelta = 1.
std::array< bool, D >	mirrorLow_ {}
std::array< bool, D >	mirrorHigh_ {}

Static Protected Attributes
static constexpr std::size_t	noNode = std::numeric_limits<std::size_t>::max()

Detailed Description

template<class T, int D>
class viennals::OxidationDiffusion< T, D >

Solves the oxidant diffusion step of the Suvorov et al. (10.1007/s10825-006-0003-z) oxidation model on the Cartesian grid carrying two level sets.

The oxide is assumed to lie between the Si/SiO2 reaction interface and the SiO2/O2 ambient interface. Regular Cartesian grid nodes inside this band are solved with finite differences. Whenever an axis-aligned grid edge leaves the oxide, a cross-point is inserted at the zero of the crossed level-set function and a Robin boundary condition is applied there. The diffusion update uses the resulting sub-grid distances in a nonuniform three-point stencil along each Cartesian axis.

The default sign convention matches the multilayer convention described in the paper: oxide is above the reaction interface and below the ambient interface, i.e. reactionPhi >= 0 and ambientPhi <= 0.

Member Typedef Documentation

◆ PressureLookup

template<class T, int D>

using viennals::OxidationDiffusion< T, D >::PressureLookup = std::unordered_map<std::size_t, T>

◆ StressLookup

template<class T, int D>

using viennals::OxidationDiffusion< T, D >::StressLookup = std::unordered_map<std::size_t, std::array<T, 9>>

Constructor & Destructor Documentation

◆ OxidationDiffusion() [1/2]

template<class T, int D>

viennals::OxidationDiffusion< T, D >::OxidationDiffusion ( )

default

◆ OxidationDiffusion() [2/2]

template<class T, int D>

viennals::OxidationDiffusion< T, D >::OxidationDiffusion	(	SmartPointer< Domain< T, D > >	passedReactionInterface,
		SmartPointer< Domain< T, D > >	passedAmbientInterface,
		OxidationParameters	passedParameters = {} )

inline

Member Function Documentation

◆ apply()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::apply ( )

inline

◆ clearMaskInterface()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::clearMaskInterface ( )

inline

◆ clearPressureField()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::clearPressureField ( )

inline

◆ clearSolveBounds()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::clearSolveBounds ( )

inline

◆ crosses()

template<class T, int D>

bool viennals::OxidationSolverBase< T, D >::crosses	(	T	a,
		T	b ) const

inlineprotectedinherited

◆ definedLevelSetGradient()

template<class T, int D>

template<class OnGrid>

Vec3D< T > viennals::OxidationSolverBase< T, D >::definedLevelSetGradient	(	ConstSparseIterator &	it,
		const IndexType &	index,
		OnGrid	onGrid ) const

inlineprotectedinherited

◆ findNearbyNode()

template<class T, int D>

std::size_t viennals::OxidationSolverBase< T, D >::findNearbyNode ( const IndexType & index ) const

inlineprotectedinherited

◆ findNearbyNodes()

template<class T, int D>

template<class Accept>

std::vector< std::size_t > viennals::OxidationSolverBase< T, D >::findNearbyNodes	(	const IndexType &	index,
		Accept	accept ) const

inlineprotectedinherited

◆ frontVelocityAt()

template<class T, int D>

FrontVelocity viennals::OxidationDiffusion< T, D >::frontVelocityAt ( const Vec3D< T > & coordinate ) const

inline

◆ getConcentration() [1/2]

template<class T, int D>

T viennals::OxidationDiffusion< T, D >::getConcentration ( const IndexType & index ) const

inline

◆ getConcentration() [2/2]

template<class T, int D>

T viennals::OxidationDiffusion< T, D >::getConcentration ( const Vec3D< T > & coordinate ) const

inline

◆ getConcentrationCache()

template<class T, int D>

const std::unordered_map< std::size_t, T > & viennals::OxidationDiffusion< T, D >::getConcentrationCache ( ) const

inline

◆ getDissipationAlpha()

template<class T, int D>

T viennals::OxidationDiffusion< T, D >::getDissipationAlpha	(	int	,
		int	material,
		const Vec3D< T > &	)

inlinefinalvirtual

If lsLocalLaxFriedrichsAnalytical is used as the spatial discretization scheme, this is called to provide the analytical solution for the alpha values, needed for numerical stability.

Reimplemented from viennals::VelocityField< T >.

◆ getEffectiveReactionRate()

template<class T, int D>

T viennals::OxidationDiffusion< T, D >::getEffectiveReactionRate ( const Vec3D< T > & coordinate ) const

inline

◆ getFeedback()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::getFeedback	(	const IndexType &	index,
		T &	pressure,
		std::array< T, 9 > &	stress ) const

inline

The pressure and stress the last solve used at a node (reference pressure / zero stress where nothing was fed back).

◆ getInterfaceNormalStress()

template<class T, int D>

const std::vector< T > & viennals::OxidationDiffusion< T, D >::getInterfaceNormalStress ( ) const

inline

The interface normal stress the reaction law used in the last solve (per node; NaN where the law falls back to the unstressed rate).

◆ getIterations()

template<class T, int D>

unsigned viennals::OxidationDiffusion< T, D >::getIterations ( ) const

inline

◆ getNormalizedResidual()

template<class T, int D>

T viennals::OxidationDiffusion< T, D >::getNormalizedResidual ( ) const

inline

◆ getNumberOfSolutionNodes()

template<class T, int D>

std::size_t viennals::OxidationDiffusion< T, D >::getNumberOfSolutionNodes ( ) const

inline

◆ getParameters()

template<class T, int D>

OxidationParameters viennals::OxidationDiffusion< T, D >::getParameters ( ) const

inline

◆ getPressureCache()

template<class T, int D>

const PressureLookup & viennals::OxidationDiffusion< T, D >::getPressureCache ( ) const

inline

Pressure and stress feedback carried from one solve to the next, like the concentration cache: lsOxidation builds a new OxidationDiffusion for every coupled solve, and without the carry-over every solve began its stress coupling from zero stress (~7 coupling iterations per solve on a Kao cylinder, residual 1.0 at the first).

◆ getReactionBoundaryConcentration() [1/2]

template<class T, int D>

T viennals::OxidationDiffusion< T, D >::getReactionBoundaryConcentration ( const IndexType & index ) const

inline

◆ getReactionBoundaryConcentration() [2/2]

template<class T, int D>

T viennals::OxidationDiffusion< T, D >::getReactionBoundaryConcentration ( const Vec3D< T > & coordinate ) const

inline

◆ getReactionBoundarySample()

template<class T, int D>

ReactionBoundarySample viennals::OxidationDiffusion< T, D >::getReactionBoundarySample ( const Vec3D< T > & coordinate ) const

inline

Return the reaction boundary sample for the grid node nearest to coordinate. Used by the deformation solver to get both the sub-grid concentration and the crossing edge so it can compute a sub-grid normal.

◆ getResidual()

template<class T, int D>

T viennals::OxidationDiffusion< T, D >::getResidual ( ) const

inline

◆ getScalarVelocity()

template<class T, int D>

T viennals::OxidationDiffusion< T, D >::getScalarVelocity	(	const Vec3D< T > &	coordinate,
		int	material,
		const Vec3D< T > &	normalVector,
		unsigned long	)

inlinefinalvirtual

Should return a scalar value for the velocity at coordinate for a point of material with the given normalVector.

Reimplemented from viennals::VelocityField< T >.

◆ getScalarVelocityFromSample()

template<class T, int D>

T viennals::OxidationDiffusion< T, D >::getScalarVelocityFromSample ( const ReactionBoundarySample & sample ) const

inline

Absolute scalar velocity derived from an already-computed boundary sample, avoiding a second call to reactionBoundarySample inside getScalarVelocity.

◆ getStressCache()

template<class T, int D>

const StressLookup & viennals::OxidationDiffusion< T, D >::getStressCache ( ) const

inline

◆ getVectorVelocity()

template<class T>

virtual Vec3D< T > viennals::VelocityField< T >::getVectorVelocity	(	const Vec3D< T > &	,
		int	,
		const Vec3D< T > &	,
		unsigned long	)

inlinevirtualinherited

Like getScalarVelocity, but returns a velocity value for each cartesian direction.

Reimplemented in viennals::OxidationDeformation< T, D >, viennals::OxidationMaskCarry< T, D >, and viennals::OxidationMaterialVelocity< T, D >.

◆ hasFiniteConcentrationField()

template<class T, int D>

bool viennals::OxidationDiffusion< T, D >::hasFiniteConcentrationField ( ) const

inline

◆ inBounds()

template<class T, int D>

bool viennals::OxidationSolverBase< T, D >::inBounds ( const IndexType & index ) const

inlineprotectedinherited

◆ increment()

template<class T, int D>

bool viennals::OxidationSolverBase< T, D >::increment ( IndexType & index ) const

inlineprotectedinherited

◆ initializeGridFromInterfaces()

template<class T, int D>

bool viennals::OxidationSolverBase< T, D >::initializeGridFromInterfaces	(	SmartPointer< Domain< T, D > >	reactionInterface,
		SmartPointer< Domain< T, D > >	ambientInterface,
		SmartPointer< Domain< T, D > >	maskInterface,
		bool	useRequestedBounds,
		const IndexType &	requestedMinIndex,
		const IndexType &	requestedMaxIndex,
		std::size_t	maxGridPoints,
		const std::string &	solverName )

inlineprotectedinherited

◆ initializeGridFromMask()

template<class T, int D>

bool viennals::OxidationSolverBase< T, D >::initializeGridFromMask	(	SmartPointer< Domain< T, D > >	maskInterface,
		bool	useRequestedBounds,
		const IndexType &	requestedMinIndex,
		const IndexType &	requestedMaxIndex,
		std::size_t	maxGridPoints,
		const std::string &	solverName )

inlineprotectedinherited

◆ initNodeLookup()

template<class T, int D>

void viennals::OxidationSolverBase< T, D >::initNodeLookup ( )

inlineprotectedinherited

◆ interfaceNormalStressFor()

template<class T, int D>

std::vector< T > viennals::OxidationDiffusion< T, D >::interfaceNormalStressFor ( const StressLookup & lookup ) const

inline

The interface normal stress the reaction law would use with lookup: the median sample at every node with a REACTION face, smoothed along the interface (see updateInterfaceNormalStress for why). Non-mutating, so the coupling loop can measure convergence without changing the state.

◆ lastSolveConverged()

template<class T, int D>

bool viennals::OxidationDiffusion< T, D >::lastSolveConverged ( ) const

inline

◆ linearIndex()

template<class T, int D>

std::size_t viennals::OxidationSolverBase< T, D >::linearIndex ( const IndexType & index ) const

inlineprotectedinherited

◆ lookupNode()

template<class T, int D>

std::size_t viennals::OxidationSolverBase< T, D >::lookupNode ( const IndexType & index ) const

inlineprotectedinherited

◆ markGeometryChanged()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::markGeometryChanged ( )

inline

Call after any in-place modification of the level sets (e.g. after ls::Advect) so that the next apply() rebuilds the Cartesian node grid.

◆ markSolved()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::markSolved ( )

inline

Mark the current solution as valid without re-solving. Call this before any parallel advection (lsAdvect) that uses this field as a velocity source to prevent concurrent apply() calls inside getScalarVelocity().

◆ New()

template<class T, int D>

template<class... Args>

auto viennals::OxidationDiffusion< T, D >::New ( Args &&... args )

inlinestatic

◆ setAmbientInterface()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::setAmbientInterface ( SmartPointer< Domain< T, D > > passedInterface )

inline

◆ setConcentrationCache()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::setConcentrationCache ( std::unordered_map< std::size_t, T > cache )

inline

◆ setFeedbackCache()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::setFeedbackCache	(	PressureLookup	pressures,
		StressLookup	stresses )

inline

◆ setGpuMode()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::setGpuMode ( GpuMode mode )

inline

Set the GPU solver selection mode (see GpuMode). On CPU-only builds (VIENNALS_GPU_BICGSTAB not defined) Gpu is an error and Auto the CPU.

◆ setMaskInterface()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::setMaskInterface	(	SmartPointer< Domain< T, D > >	passedInterface,
		int	passedMaskSign = 1 )

inline

◆ setMirrorSides()

template<class T, int D>

template<class GridType>

void viennals::OxidationSolverBase< T, D >::setMirrorSides ( const GridType & grid )

inlineprotectedinherited

◆ setOxideSigns()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::setOxideSigns	(	int	passedReactionSign,
		int	passedAmbientSign )

inline

Set signs defining the oxide band. A node is inside oxide if reactionSign * reactionPhi >= 0 and ambientSign * ambientPhi >= 0.

◆ setParameters()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::setParameters ( OxidationParameters passedParameters )

inline

◆ setPressure() [1/2]

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::setPressure	(	const IndexType &	index,
		T	pressure )

inline

◆ setPressure() [2/2]

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::setPressure	(	const Vec3D< T > &	coordinate,
		T	pressure )

inline

◆ setReactionInterface()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::setReactionInterface ( SmartPointer< Domain< T, D > > passedInterface )

inline

◆ setSolveBounds()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::setSolveBounds	(	const IndexType &	passedMinIndex,
		const IndexType &	passedMaxIndex )

inline

Restrict the dense Cartesian diffusion solve to a finite index box. This is useful for level sets with infinite boundary conditions.

◆ setStress() [1/2]

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::setStress	(	const IndexType &	index,
		const std::array< T, 9 > &	stress )

inline

Oxide Cauchy stress (Pa, tension positive, row-major 3x3) at a node, for the normal-stress reaction law.

◆ setStress() [2/2]

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::setStress	(	const Vec3D< T > &	coordinate,
		const std::array< T, 9 > &	stress )

inline

◆ stencilOffset()

template<class T, int D>

int viennals::OxidationSolverBase< T, D >::stencilOffset	(	const IndexType &	index,
		unsigned	direction,
		int	offset ) const

inlineprotectedinherited

◆ valueAt()

template<class T, int D>

T viennals::OxidationSolverBase< T, D >::valueAt	(	ConstSparseIterator &	it,
		const IndexType &	index ) const

inlineprotectedinherited

◆ writeConcentrationToLevelSet()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::writeConcentrationToLevelSet ( )

inline

Write per-node concentration into ambientInterface->getPointData() so that lsInterior + lsAdvect can carry it across timestep boundaries. Safe to call regardless of the solved flag; uses the most recent nodes.

◆ writePersistentFields()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::writePersistentFields ( )

inline

Convenience wrapper: persist both concentration and pressure in one call.

◆ writePressureToLevelSet()

template<class T, int D>

void viennals::OxidationDiffusion< T, D >::writePressureToLevelSet ( )

inline

Write per-node pressure into ambientInterface->getPointData() so that it survives advection and can warm-start the coupling loop next step.

Member Data Documentation

◆ extents

template<class T, int D>

std::array<std::size_t, D> viennals::OxidationSolverBase< T, D >::extents {}

protectedinherited

◆ gridDelta

template<class T, int D>

T viennals::OxidationSolverBase< T, D >::gridDelta = 1.

protectedinherited

◆ maxIndex

template<class T, int D>

IndexType viennals::OxidationSolverBase< T, D >::maxIndex {}

protectedinherited

◆ minIndex

template<class T, int D>

IndexType viennals::OxidationSolverBase< T, D >::minIndex {}

protectedinherited

◆ mirrorHigh_

template<class T, int D>

std::array<bool, D> viennals::OxidationSolverBase< T, D >::mirrorHigh_ {}

protectedinherited

◆ mirrorLow_

template<class T, int D>

std::array<bool, D> viennals::OxidationSolverBase< T, D >::mirrorLow_ {}

protectedinherited

◆ nodeLookupFlat

template<class T, int D>

std::vector<std::size_t> viennals::OxidationSolverBase< T, D >::nodeLookupFlat

protectedinherited

◆ noNode

template<class T, int D>

std::size_t viennals::OxidationSolverBase< T, D >::noNode = std::numeric_limits<std::size_t>::max()

staticconstexprprotectedinherited

◆ strides

template<class T, int D>

std::array<std::size_t, D> viennals::OxidationSolverBase< T, D >::strides {}

protectedinherited

The documentation for this class was generated from the following file:

include/viennals/lsOxidationDiffusion.hpp

Classes

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Static Protected Attributes

Detailed Description

Member Typedef Documentation

◆ PressureLookup

◆ StressLookup

Constructor & Destructor Documentation

◆ OxidationDiffusion() [1/2]

◆ OxidationDiffusion() [2/2]

Member Function Documentation

◆ apply()

◆ clearMaskInterface()

◆ clearPressureField()

◆ clearSolveBounds()

◆ crosses()

◆ definedLevelSetGradient()

◆ findNearbyNode()

◆ findNearbyNodes()

◆ frontVelocityAt()

◆ getConcentration() [1/2]

◆ getConcentration() [2/2]

◆ getConcentrationCache()

◆ getDissipationAlpha()

◆ getEffectiveReactionRate()

◆ getFeedback()

◆ getInterfaceNormalStress()

◆ getIterations()

◆ getNormalizedResidual()

◆ getNumberOfSolutionNodes()

◆ getParameters()

◆ getPressureCache()

◆ getReactionBoundaryConcentration() [1/2]

◆ getReactionBoundaryConcentration() [2/2]

◆ getReactionBoundarySample()

◆ getResidual()

◆ getScalarVelocity()

◆ getScalarVelocityFromSample()

◆ getStressCache()

◆ getVectorVelocity()

◆ hasFiniteConcentrationField()

◆ inBounds()

◆ increment()

◆ initializeGridFromInterfaces()

◆ initializeGridFromMask()

◆ initNodeLookup()

◆ interfaceNormalStressFor()

◆ lastSolveConverged()

◆ linearIndex()

◆ lookupNode()

◆ markGeometryChanged()

◆ markSolved()

◆ New()

◆ setAmbientInterface()

◆ setConcentrationCache()

◆ setFeedbackCache()

◆ setGpuMode()

◆ setMaskInterface()

◆ setMirrorSides()

◆ setOxideSigns()

◆ setParameters()

◆ setPressure() [1/2]

◆ setPressure() [2/2]

◆ setReactionInterface()

◆ setSolveBounds()

◆ setStress() [1/2]

◆ setStress() [2/2]

◆ stencilOffset()

◆ valueAt()

◆ writeConcentrationToLevelSet()

◆ writePersistentFields()

◆ writePressureToLevelSet()

Member Data Documentation

◆ extents

◆ gridDelta

◆ maxIndex