DeltaApproximationHeavisidePhaseWeighted< number > Class Template Reference

Developer Documentation: MeltPoolDG::LevelSet::DeltaApproximationHeavisidePhaseWeighted< number > Class Template Reference
Developer Documentation
MeltPoolDG::LevelSet::DeltaApproximationHeavisidePhaseWeighted< number > Class Template Reference

#include <delta_approximation_phase_weighted.hpp>

Inheritance diagram for MeltPoolDG::LevelSet::DeltaApproximationHeavisidePhaseWeighted< number >:
[legend]
Collaboration diagram for MeltPoolDG::LevelSet::DeltaApproximationHeavisidePhaseWeighted< number >:
[legend]

Public Member Functions

 DeltaApproximationHeavisidePhaseWeighted (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::DeltaApproximationHeavisidePhaseWeighted< number >

Asymmetric, phase weighted Dirac delta approximation function, as proposed in

Yokoi, K. (2014). A density-scaled continuum surface force model within a balanced force formulation. Journal of Computational Physics, 278, 221-228. DOI: 10.1016/j.jcp.2014.08.034

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

    2 ( (1 - φ) w_g + φ w_h )

W(φ) = ------------------------— , 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 they don't sum up to zero.

Constructor & Destructor Documentation

◆ DeltaApproximationHeavisidePhaseWeighted()

template<typename number >
MeltPoolDG::LevelSet::DeltaApproximationHeavisidePhaseWeighted< number >::DeltaApproximationHeavisidePhaseWeighted ( const DeltaApproximationPhaseWeightedData< number > &  data)
inline

Member Function Documentation

◆ compute_weight() [1/2]

template<typename number >
dealii::VectorizedArray< number > MeltPoolDG::LevelSet::DeltaApproximationHeavisidePhaseWeighted< number >::compute_weight ( const dealii::VectorizedArray< number > &  level_set_heaviside) const
inlineoverridevirtual

◆ compute_weight() [2/2]

template<typename number >
number MeltPoolDG::LevelSet::DeltaApproximationHeavisidePhaseWeighted< number >::compute_weight ( const number  level_set_heaviside) const
inlineoverridevirtual

◆ compute_weight_internal()

template<typename number >
template<typename value_type >
value_type MeltPoolDG::LevelSet::DeltaApproximationHeavisidePhaseWeighted< number >::compute_weight_internal ( const value_type &  level_set_heaviside) const
inlineprivate

This function calculates the asymmetric

    2 ( (1 - φ) w_g + φ w_h )

W(φ) = ------------------------— w_g + w_h

Member Data Documentation

◆ correction_factor

template<typename number >
const number MeltPoolDG::LevelSet::DeltaApproximationHeavisidePhaseWeighted< number >::correction_factor
private

◆ w_g

template<typename number >
const number MeltPoolDG::LevelSet::DeltaApproximationHeavisidePhaseWeighted< number >::w_g
private

◆ w_h

template<typename number >
const number MeltPoolDG::LevelSet::DeltaApproximationHeavisidePhaseWeighted< number >::w_h
private

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