WALE (Wall-Adapting Local Eddy-viscosity) eddy viscosity model. Derived from LargeEddySimulation_Smagorinsky for only modifying compute_eddy_viscosity. More...

#include <large_eddy_simulation.h>

Inheritance diagram for PHiLiP::Physics::LargeEddySimulation_WALE< dim, nstate, real >:

Collaboration diagram for PHiLiP::Physics::LargeEddySimulation_WALE< dim, nstate, real >:

Public Types
using	thermal_boundary_condition_enum = Parameters::NavierStokesParam::ThermalBoundaryCondition

using	two_point_num_flux_enum = Parameters::AllParameters::TwoPointNumericalFlux

Public Types inherited from PHiLiP::Physics::LargeEddySimulation_Smagorinsky< dim, nstate, real >
using	thermal_boundary_condition_enum = Parameters::NavierStokesParam::ThermalBoundaryCondition

using	two_point_num_flux_enum = Parameters::AllParameters::TwoPointNumericalFlux

Public Types inherited from PHiLiP::Physics::LargeEddySimulationBase< dim, nstate, real >
using	thermal_boundary_condition_enum = Parameters::NavierStokesParam::ThermalBoundaryCondition

using	two_point_num_flux_enum = Parameters::AllParameters::TwoPointNumericalFlux

Public Member Functions
	LargeEddySimulation_WALE (const Parameters::AllParameters *const parameters_input, const double ref_length, const double gamma_gas, const double mach_inf, const double angle_of_attack, const double side_slip_angle, const double prandtl_number, const double reynolds_number_inf, const bool use_constant_viscosity, const double constant_viscosity, const double temperature_inf, const double turbulent_prandtl_number, const double ratio_of_filter_width_to_cell_size, const double model_constant, const double isothermal_wall_temperature=1.0, const thermal_boundary_condition_enum thermal_boundary_condition_type=thermal_boundary_condition_enum::adiabatic, std::shared_ptr< ManufacturedSolutionFunction< dim, real > > manufactured_solution_function=nullptr, const two_point_num_flux_enum two_point_num_flux_type=two_point_num_flux_enum::KG)

real	compute_eddy_viscosity (const std::array< real, nstate > &primitive_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &primitive_soln_gradient, const dealii::types::global_dof_index cell_index) const override

FadType	compute_eddy_viscosity_fad (const std::array< FadType, nstate > &primitive_soln, const std::array< dealii::Tensor< 1, dim, FadType >, nstate > &primitive_soln_gradient, const dealii::types::global_dof_index cell_index) const override

Public Member Functions inherited from PHiLiP::Physics::LargeEddySimulation_Smagorinsky< dim, nstate, real >
	LargeEddySimulation_Smagorinsky (const Parameters::AllParameters *const parameters_input, const double ref_length, const double gamma_gas, const double mach_inf, const double angle_of_attack, const double side_slip_angle, const double prandtl_number, const double reynolds_number_inf, const bool use_constant_viscosity, const double constant_viscosity, const double temperature_inf, const double turbulent_prandtl_number, const double ratio_of_filter_width_to_cell_size, const double model_constant, const double isothermal_wall_temperature=1.0, const thermal_boundary_condition_enum thermal_boundary_condition_type=thermal_boundary_condition_enum::adiabatic, std::shared_ptr< ManufacturedSolutionFunction< dim, real > > manufactured_solution_function=nullptr, const two_point_num_flux_enum two_point_num_flux_type=two_point_num_flux_enum::KG)

double	get_model_constant_times_filter_width (const dealii::types::global_dof_index cell_index) const
	Returns the product of the eddy viscosity model constant and the filter width.

dealii::Tensor< 2, dim, real >	compute_SGS_stress_tensor (const std::array< real, nstate > &primitive_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &primitive_soln_gradient, const dealii::types::global_dof_index cell_index) const
	Nondimensionalized sub-grid scale (SGS) stress tensor, (tau^sgs)*.

dealii::Tensor< 1, dim, real >	compute_SGS_heat_flux (const std::array< real, nstate > &primitive_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &primitive_soln_gradient, const dealii::types::global_dof_index cell_index) const
	Nondimensionalized sub-grid scale (SGS) heat flux, (q^sgs)*.

dealii::Tensor< 2, dim, FadType >	compute_SGS_stress_tensor_fad (const std::array< FadType, nstate > &primitive_soln, const std::array< dealii::Tensor< 1, dim, FadType >, nstate > &primitive_soln_gradient, const dealii::types::global_dof_index cell_index) const
	Nondimensionalized sub-grid scale (SGS) stress tensor, (tau^sgs)* (Automatic Differentiation Type: FadType)

dealii::Tensor< 1, dim, FadType >	compute_SGS_heat_flux_fad (const std::array< FadType, nstate > &primitive_soln, const std::array< dealii::Tensor< 1, dim, FadType >, nstate > &primitive_soln_gradient, const dealii::types::global_dof_index cell_index) const
	Nondimensionalized sub-grid scale (SGS) heat flux, (q^sgs)* (Automatic Differentiation Type: FadType)

Public Member Functions inherited from PHiLiP::Physics::LargeEddySimulationBase< dim, nstate, real >
	LargeEddySimulationBase (const Parameters::AllParameters *const parameters_input, const double ref_length, const double gamma_gas, const double mach_inf, const double angle_of_attack, const double side_slip_angle, const double prandtl_number, const double reynolds_number_inf, const bool use_constant_viscosity, const double constant_viscosity, const double temperature_inf, const double turbulent_prandtl_number, const double ratio_of_filter_width_to_cell_size, const double isothermal_wall_temperature=1.0, const thermal_boundary_condition_enum thermal_boundary_condition_type=thermal_boundary_condition_enum::adiabatic, std::shared_ptr< ManufacturedSolutionFunction< dim, real > > manufactured_solution_function=nullptr, const two_point_num_flux_enum two_point_num_flux_type=two_point_num_flux_enum::KG)
	Constructor.

std::array< dealii::Tensor< 1, dim, real >, nstate >	convective_flux (const std::array< real, nstate > &conservative_soln) const
	Convective flux: $ \mathbf{F}_{conv} $.

std::array< dealii::Tensor< 1, dim, real >, nstate >	dissipative_flux (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &solution_gradient, const dealii::types::global_dof_index cell_index) const
	Dissipative (i.e. viscous) flux: $ \mathbf{F}_{diss} $.

std::array< real, nstate >	convective_eigenvalues (const std::array< real, nstate > &, const dealii::Tensor< 1, dim, real > &) const override
	Convective eigenvalues of the additional models' PDEs. More...

real	max_convective_eigenvalue (const std::array< real, nstate > &soln) const
	Maximum convective eigenvalue of the additional models' PDEs. More...

real	max_convective_normal_eigenvalue (const std::array< real, nstate > &soln, const dealii::Tensor< 1, dim, real > &normal) const
	Maximum convective normal eigenvalue (used in Lax-Friedrichs) of the additional models' PDEs. More...

std::array< real, nstate >	source_term (const dealii::Point< dim, real > &pos, const std::array< real, nstate > &solution, const real current_time, const dealii::types::global_dof_index cell_index) const
	Source term for manufactured solution functions.

double	get_filter_width (const dealii::types::global_dof_index cell_index) const
	Compute the nondimensionalized filter width used by the SGS model given a cell index. More...

Public Member Functions inherited from PHiLiP::Physics::ModelBase< dim, nstate, real >
	ModelBase (std::shared_ptr< ManufacturedSolutionFunction< dim, real > > manufactured_solution_function_input=nullptr)
	Constructor.

virtual	~ModelBase ()=default
	Virtual destructor required for abstract classes.

virtual std::array< real, nstate >	physical_source_term (const dealii::Point< dim, real > &pos, const std::array< real, nstate > &solution, const std::array< dealii::Tensor< 1, dim, real >, nstate > &solution_gradient, const dealii::types::global_dof_index cell_index) const
	Physical source terms additional to the baseline physics (including physical source terms in additional PDEs of model)

void	boundary_face_values (const int boundary_type, const dealii::Point< dim, real > &pos, const dealii::Tensor< 1, dim, real > &normal, const std::array< real, nstate > &soln_int, const std::array< dealii::Tensor< 1, dim, real >, nstate > &soln_grad_int, std::array< real, nstate > &soln_bc, std::array< dealii::Tensor< 1, dim, real >, nstate > &soln_grad_bc) const
	Boundary condition handler.

virtual dealii::Vector< double >	post_compute_derived_quantities_vector (const dealii::Vector< double > &uh, const std::vector< dealii::Tensor< 1, dim > > &, const std::vector< dealii::Tensor< 2, dim > > &, const dealii::Tensor< 1, dim > &, const dealii::Point< dim > &) const
	Returns current vector solution to be used by PhysicsPostprocessor to output current solution. More...

virtual std::vector< dealii::DataComponentInterpretation::DataComponentInterpretation >	post_get_data_component_interpretation () const
	Returns DataComponentInterpretation of the solution to be used by PhysicsPostprocessor to output current solution. More...

virtual std::vector< std::string >	post_get_names () const
	Returns names of the solution to be used by PhysicsPostprocessor to output current solution. More...

Private Member Functions
template<typename real2 >
real2	compute_eddy_viscosity_templated (const std::array< real2, nstate > &primitive_soln, const std::array< dealii::Tensor< 1, dim, real2 >, nstate > &primitive_soln_gradient, const dealii::types::global_dof_index cell_index) const

Additional Inherited Members
Public Attributes inherited from PHiLiP::Physics::LargeEddySimulation_Smagorinsky< dim, nstate, real >
const double	model_constant
	SGS model constant.

Public Attributes inherited from PHiLiP::Physics::LargeEddySimulationBase< dim, nstate, real >
const double	turbulent_prandtl_number
	Turbulent Prandtl number.

const double	ratio_of_filter_width_to_cell_size
	Ratio of filter width to cell size.

std::unique_ptr< NavierStokes< dim, nstate, real > >	navier_stokes_physics
	Pointer to Navier-Stokes physics object.

Public Attributes inherited from PHiLiP::Physics::ModelBase< dim, nstate, real >
std::shared_ptr< ManufacturedSolutionFunction< dim, real > >	manufactured_solution_function
	Manufactured solution function.

dealii::LinearAlgebra::distributed::Vector< int >	cellwise_poly_degree
	Cellwise polynomial degree.

dealii::LinearAlgebra::distributed::Vector< double >	cellwise_volume

Protected Member Functions inherited from PHiLiP::Physics::LargeEddySimulation_Smagorinsky< dim, nstate, real >
template<typename real2 >
dealii::Tensor< 2, dim, real2 >	compute_SGS_stress_tensor_templated (const std::array< real2, nstate > &primitive_soln, const std::array< dealii::Tensor< 1, dim, real2 >, nstate > &primitive_soln_gradient, const dealii::types::global_dof_index cell_index) const
	Templated nondimensionalized sub-grid scale (SGS) stress tensor, (tau^sgs)*.

template<typename real2 >
dealii::Tensor< 1, dim, real2 >	compute_SGS_heat_flux_templated (const std::array< real2, nstate > &primitive_soln, const std::array< dealii::Tensor< 1, dim, real2 >, nstate > &primitive_soln_gradient, const dealii::types::global_dof_index cell_index) const
	Templated nondimensionalized sub-grid scale (SGS) heat flux, (q^sgs)*.

template<typename real2 >
real2	scale_eddy_viscosity_templated (const std::array< real2, nstate > &primitive_soln, const real2 eddy_viscosity) const
	Templated scale nondimensionalized eddy viscosity for Smagorinsky model.

Protected Member Functions inherited from PHiLiP::Physics::LargeEddySimulationBase< dim, nstate, real >
template<typename real2 >
real2	get_tensor_magnitude_sqr (const dealii::Tensor< 2, dim, real2 > &tensor) const
	Returns the square of the magnitude of the tensor (i.e. the double dot product of a tensor with itself)

template<typename real2 >
std::array< dealii::Tensor< 1, dim, real2 >, nstate >	dissipative_flux_templated (const std::array< real2, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real2 >, nstate > &solution_gradient, const dealii::types::global_dof_index cell_index) const
	Templated dissipative (i.e. viscous) flux: $ \mathbf{F}_{diss} $.

dealii::Tensor< 2, nstate, real >	dissipative_flux_directional_jacobian (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &solution_gradient, const dealii::Tensor< 1, dim, real > &normal, const dealii::types::global_dof_index cell_index) const

dealii::Tensor< 2, nstate, real >	dissipative_flux_directional_jacobian_wrt_gradient_component (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &solution_gradient, const dealii::Tensor< 1, dim, real > &normal, const int d_gradient, const dealii::types::global_dof_index cell_index) const

std::array< real, nstate >	get_manufactured_solution_value (const dealii::Point< dim, real > &pos) const
	Get manufactured solution value (repeated from Euler)

std::array< dealii::Tensor< 1, dim, real >, nstate >	get_manufactured_solution_gradient (const dealii::Point< dim, real > &pos) const
	Get manufactured solution value (repeated from Euler)

std::array< real, nstate >	dissipative_source_term (const dealii::Point< dim, real > &pos, const dealii::types::global_dof_index cell_index) const
	Dissipative flux contribution to the source term (repeated from NavierStokes) More...

Protected Member Functions inherited from PHiLiP::Physics::ModelBase< dim, nstate, real >
virtual void	boundary_manufactured_solution (const dealii::Point< dim, real > &pos, const dealii::Tensor< 1, dim, real > &normal_int, const std::array< real, nstate > &soln_int, const std::array< dealii::Tensor< 1, dim, real >, nstate > &soln_grad_int, std::array< real, nstate > &soln_bc, std::array< dealii::Tensor< 1, dim, real >, nstate > &soln_grad_bc) const
	Evaluate the manufactured solution boundary conditions.

virtual void	boundary_wall (std::array< real, nstate > &soln_bc, std::array< dealii::Tensor< 1, dim, real >, nstate > &soln_grad_bc) const
	Wall boundary condition.

virtual void	boundary_outflow (const std::array< real, nstate > &soln_int, const std::array< dealii::Tensor< 1, dim, real >, nstate > &soln_grad_int, std::array< real, nstate > &soln_bc, std::array< dealii::Tensor< 1, dim, real >, nstate > &soln_grad_bc) const
	Outflow Boundary Condition.

virtual void	boundary_inflow (const std::array< real, nstate > &soln_int, const std::array< dealii::Tensor< 1, dim, real >, nstate > &soln_grad_int, std::array< real, nstate > &soln_bc, std::array< dealii::Tensor< 1, dim, real >, nstate > &soln_grad_bc) const
	Inflow boundary conditions.

virtual void	boundary_farfield (std::array< real, nstate > &soln_bc) const
	Farfield boundary conditions based on freestream values.

virtual void	boundary_riemann (const dealii::Tensor< 1, dim, real > &normal_int, const std::array< real, nstate > &soln_int, std::array< real, nstate > &soln_bc) const

virtual void	boundary_slip_wall (const dealii::Tensor< 1, dim, real > &normal_int, const std::array< real, nstate > &soln_int, const std::array< dealii::Tensor< 1, dim, real >, nstate > &soln_grad_int, std::array< real, nstate > &soln_bc, std::array< dealii::Tensor< 1, dim, real >, nstate > &soln_grad_bc) const

Protected Attributes inherited from PHiLiP::Physics::ModelBase< dim, nstate, real >
const MPI_Comm	mpi_communicator
	MPI communicator.

dealii::ConditionalOStream	pcout
	Parallel std::cout that only outputs on mpi_rank==0.

Detailed Description

template<int dim, int nstate, typename real>
class PHiLiP::Physics::LargeEddySimulation_WALE< dim, nstate, real >

WALE (Wall-Adapting Local Eddy-viscosity) eddy viscosity model. Derived from LargeEddySimulation_Smagorinsky for only modifying compute_eddy_viscosity.

Definition at line 255 of file large_eddy_simulation.h.

Constructor & Destructor Documentation

◆ LargeEddySimulation_WALE()

template<int dim, int nstate, typename real >

PHiLiP::Physics::LargeEddySimulation_WALE< dim, nstate, real >::LargeEddySimulation_WALE	(	const Parameters::AllParameters *const	parameters_input,
		const double	ref_length,
		const double	gamma_gas,
		const double	mach_inf,
		const double	angle_of_attack,
		const double	side_slip_angle,
		const double	prandtl_number,
		const double	reynolds_number_inf,
		const bool	use_constant_viscosity,
		const double	constant_viscosity,
		const double	temperature_inf,
		const double	turbulent_prandtl_number,
		const double	ratio_of_filter_width_to_cell_size,
		const double	model_constant,
		const double	isothermal_wall_temperature = `1.0`,
		const thermal_boundary_condition_enum	thermal_boundary_condition_type = `thermal_boundary_condition_enum::adiabatic`,
		std::shared_ptr< ManufacturedSolutionFunction< dim, real > >	manufactured_solution_function = `nullptr`,
		const two_point_num_flux_enum	two_point_num_flux_type = `two_point_num_flux_enum::KG`
	)

Constructor for the sub-grid scale (SGS) model: WALE Reference 1: de la Llave Plata et al. (2019). "On the performance of a high-order multiscale DG approach to LES at increasing Reynolds number." Reference 2: Nicoud & Ducros (1999) "Subgrid-scale stress modelling based on the square of the velocity gradient tensor"

Definition at line 638 of file large_eddy_simulation.cpp.

Member Function Documentation

◆ compute_eddy_viscosity()

template<int dim, int nstate, typename real >

real PHiLiP::Physics::LargeEddySimulation_WALE< dim, nstate, real >::compute_eddy_viscosity	(	const std::array< real, nstate > &	primitive_soln,
		const std::array< dealii::Tensor< 1, dim, real >, nstate > &	primitive_soln_gradient,
		const dealii::types::global_dof_index	cell_index
	)		const

overridevirtual

Nondimensionalized eddy viscosity for the WALE model. Reference: Nicoud & Ducros (1999) "Subgrid-scale stress modelling based on the square of the velocity gradient tensor"

Reimplemented from PHiLiP::Physics::LargeEddySimulation_Smagorinsky< dim, nstate, real >.

Definition at line 679 of file large_eddy_simulation.cpp.

◆ compute_eddy_viscosity_fad()

template<int dim, int nstate, typename real >

FadType PHiLiP::Physics::LargeEddySimulation_WALE< dim, nstate, real >::compute_eddy_viscosity_fad	(	const std::array< FadType, nstate > &	primitive_soln,
		const std::array< dealii::Tensor< 1, dim, FadType >, nstate > &	primitive_soln_gradient,
		const dealii::types::global_dof_index	cell_index
	)		const

overridevirtual

Nondimensionalized eddy viscosity for the WALE model. (Automatic Differentiation Type: FadType) Reference: Nicoud & Ducros (1999) "Subgrid-scale stress modelling based on the square of the velocity gradient tensor"

Reimplemented from PHiLiP::Physics::LargeEddySimulation_Smagorinsky< dim, nstate, real >.

Definition at line 689 of file large_eddy_simulation.cpp.

◆ compute_eddy_viscosity_templated()

template<int dim, int nstate, typename real >

template<typename real2 >

real2 PHiLiP::Physics::LargeEddySimulation_WALE< dim, nstate, real >::compute_eddy_viscosity_templated	(	const std::array< real2, nstate > &	primitive_soln,
		const std::array< dealii::Tensor< 1, dim, real2 >, nstate > &	primitive_soln_gradient,
		const dealii::types::global_dof_index	cell_index
	)		const

private

Templated nondimensionalized eddy viscosity for the WALE model. Reference: Nicoud & Ducros (1999) "Subgrid-scale stress modelling based on the square of the velocity gradient tensor"

Get traceless symmetric square of velocity gradient tensor, i.e. ${S}^{d}$ Reference: Nicoud and Ducros (1999) - Equation (10)

Eddy viscosity is zero in the absence of turbulent fluctuations, i.e. zero strain rate and zero rotation rate. See Nicoud and Ducros (1999). Since the denominator in this eddy viscosity model will go to zero, we must explicitly set the eddy viscosity to zero to avoid a division by zero. Or equivalently, update it from its zero initialization only if there is turbulence.

Definition at line 700 of file large_eddy_simulation.cpp.

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

src/physics/large_eddy_simulation.h
src/physics/large_eddy_simulation.cpp

Public Types

Public Member Functions

Private Member Functions

Additional Inherited Members

Detailed Description

template<int dim, int nstate, typename real> class PHiLiP::Physics::LargeEddySimulation_WALE< dim, nstate, real >

Constructor & Destructor Documentation

◆ LargeEddySimulation_WALE()

Member Function Documentation

◆ compute_eddy_viscosity()

◆ compute_eddy_viscosity_fad()

◆ compute_eddy_viscosity_templated()

template<int dim, int nstate, typename real>
class PHiLiP::Physics::LargeEddySimulation_WALE< dim, nstate, real >