MeltPoolDG::LevelSet::DeltaApproximationReciprocalPhaseWeighted< number > Class Template Reference

#include <delta_approximation_phase_weighted.hpp>

Inheritance diagram for MeltPoolDG::LevelSet::DeltaApproximationReciprocalPhaseWeighted< number >:

Collaboration diagram for MeltPoolDG::LevelSet::DeltaApproximationReciprocalPhaseWeighted< number >:

Public Member Functions
	DeltaApproximationReciprocalPhaseWeighted (const DeltaApproximationPhaseWeightedData< number > &data)
number	compute_weight (const number level_set_heaviside) const override
dealii::VectorizedArray< number >	compute_weight (const dealii::VectorizedArray< number > &level_set_heaviside) const override
Public Member Functions inherited from MeltPoolDG::LevelSet::DeltaApproximationBase< number >
virtual	~DeltaApproximationBase ()=default

Private Member Functions
template<typename value_type >
value_type	compute_weight_internal (const value_type &level_set_heaviside) const

Private Attributes
const number	w_g
const number	w_h
const number	correction_factor

Detailed Description

template<typename number>
class MeltPoolDG::LevelSet::DeltaApproximationReciprocalPhaseWeighted< number >

Asymmetric, Dirac delta approximation function, phase weighted with the density distribution consistent with the mass flux due to evaporation

This function can be used to approximate the Dirac delta function for diffuse interfaces. The approximation is asymmetric and you can weigh the two phases individually. The function is defined as

δ_w(φ) = δ(φ) * W(φ).

The function is dependent only on the heaviside representation of the level set φ (=indicator) The symmetric delta function δ(φ) is defined as the norm of the indicator gradient ∇φ:

δ = ||∇φ||

The weight function W(φ) is defined as

               1                       w_g - w_h

W(φ) = ----------------------— * ----------------------— φ / w_h + (1 - φ) / w_g w_g w_h ln( w_g / w_h )

where w_g is the weight of the gaseous phase (at level set = -1) and w_h is the weight of the heavy phase (at level set = 1). The weights can be chosen arbitrarily, as long as

w_g > 0 and w_h > 0 and w_g != w_h

Note

If the density is determined consistent with mass flux due to evaporation, this Dirac delta approximation can be used to scale interface quantities with the density distribution by choosing the phases' weights equal to their densities.

Constructor & Destructor Documentation

DeltaApproximationReciprocalPhaseWeighted()

template<typename number >

MeltPoolDG::LevelSet::DeltaApproximationReciprocalPhaseWeighted< number >::DeltaApproximationReciprocalPhaseWeighted ( const DeltaApproximationPhaseWeightedData< number > &  data )

inline

Member Function Documentation

compute_weight() [1/2]

template<typename number >

dealii::VectorizedArray< number > MeltPoolDG::LevelSet::DeltaApproximationReciprocalPhaseWeighted< number >::compute_weight ( const dealii::VectorizedArray< number > &  level_set_heaviside ) const

inlineoverridevirtual

Implements MeltPoolDG::LevelSet::DeltaApproximationBase< number >.

compute_weight() [2/2]

template<typename number >

number MeltPoolDG::LevelSet::DeltaApproximationReciprocalPhaseWeighted< number >::compute_weight ( const number  level_set_heaviside ) const

inlineoverridevirtual

Implements MeltPoolDG::LevelSet::DeltaApproximationBase< number >.

compute_weight_internal()

template<typename number >

template<typename value_type >

value_type MeltPoolDG::LevelSet::DeltaApproximationReciprocalPhaseWeighted< number >::compute_weight_internal ( const value_type &  level_set_heaviside ) const

inlineprivate

This function calculates the asymmetric

               1                       w_g - w_h

W(φ) = ----------------------— * ----------------------— φ / w_h + (1 - φ) / w_g w_g w_h ln( w_g / w_h )

Member Data Documentation

correction_factor

template<typename number >

const number MeltPoolDG::LevelSet::DeltaApproximationReciprocalPhaseWeighted< number >::correction_factor

private

w_g

template<typename number >

const number MeltPoolDG::LevelSet::DeltaApproximationReciprocalPhaseWeighted< number >::w_g

private

w_h

template<typename number >

const number MeltPoolDG::LevelSet::DeltaApproximationReciprocalPhaseWeighted< number >::w_h

private

---

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

• include/meltpooldg/level_set/delta_approximation_phase_weighted.hpp