Public Member Functions |
Public Attributes |
Protected Member Functions |
Protected Attributes |
List of all members
PHiLiP::Physics::ConvectionDiffusion< dim, nspecies, nstate, real > Class Template Reference
Convection-diffusion with linear advective and diffusive term. Derived from PhysicsBase. More...
#include <convection_diffusion.h>
Inheritance diagram for PHiLiP::Physics::ConvectionDiffusion< dim, nspecies, nstate, real >:
Collaboration diagram for PHiLiP::Physics::ConvectionDiffusion< dim, nspecies, nstate, real >:
Public Member Functions | |
| ConvectionDiffusion (const Parameters::AllParameters *const parameters_input, const bool convection=true, const bool diffusion=true, const dealii::Tensor< 2, 3, double > input_diffusion_tensor=Parameters::ManufacturedSolutionParam::get_default_diffusion_tensor(), const dealii::Tensor< 1, 3, double > input_advection_vector=Parameters::ManufacturedSolutionParam::get_default_advection_vector(), const double input_diffusion_coefficient=Parameters::ManufacturedSolutionParam::get_default_diffusion_coefficient(), std::shared_ptr< ManufacturedSolutionFunction< dim, nspecies, real > > manufactured_solution_function=nullptr, const Parameters::AllParameters::TestType parameters_test=Parameters::AllParameters::TestType::run_control, const bool has_nonzero_physical_source=false) | |
| Constructor. | |
| std::array< dealii::Tensor< 1, dim, real >, nstate > | convective_flux (const std::array< real, nstate > &solution) const |
| Convective flux: \( \mathbf{F}_{conv} = u \). | |
| std::array< dealii::Tensor< 1, dim, real >, nstate > | convective_numerical_split_flux (const std::array< real, nstate > &soln1, const std::array< real, nstate > &soln2) const override |
| Convective numerical split flux for split form. | |
| std::array< real, nstate > | convert_conservative_to_primitive (const std::array< real, nstate > &conservative_soln) const |
| Convert conservative variables to primitive variables. | |
| std::array< real, nstate > | convert_primitive_to_conservative (const std::array< real, nstate > &primitive_soln) const |
| Convert primitive solution to conservative solution. | |
| std::array< dealii::Tensor< 1, dim, real >, nstate > | convert_conservative_gradient_to_primitive_gradient (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &conservative_soln_gradient) const |
| std::array< dealii::Tensor< 1, dim, real >, nstate > | convert_primitive_gradient_to_conservative_gradient (const std::array< real, nstate > &primitive_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &primitive_soln_gradient) const |
| std::array< real, nstate > | compute_entropy_variables (const std::array< real, nstate > &conservative_soln) const |
| Computes the entropy variables. | |
| std::array< real, nstate > | compute_conservative_variables_from_entropy_variables (const std::array< real, nstate > &entropy_var) const |
| Computes the conservative variables from the entropy variables. | |
| std::array< real, nstate > | convective_eigenvalues (const std::array< real, nstate > &, const dealii::Tensor< 1, dim, real > &) const |
| Spectral radius of convective term Jacobian is 'c'. | |
| real | max_convective_eigenvalue (const std::array< real, nstate > &soln) const |
| Maximum convective eigenvalue. | |
| real | max_viscous_eigenvalue (const std::array< real, nstate > &soln) const |
| Maximum viscous eigenvalue. | |
| std::array< dealii::Tensor< 1, dim, real >, nstate > | dissipative_flux (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 |
| Dissipative flux: u. | |
| std::array< dealii::Tensor< 1, dim, real >, nstate > | dissipative_flux (const std::array< real, nstate > &solution, const std::array< dealii::Tensor< 1, dim, real >, nstate > &solution_gradient) const |
| (function overload) Dissipative flux: u | |
| 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 is zero or depends on manufactured solution. | |
| std::array< real, nstate > | source_term (const dealii::Point< dim, real > &pos, const std::array< real, nstate > &solution, const real current_time) const |
| (function overload) Source term is zero or depends on manufactured solution | |
| void | boundary_face_values (const int, const dealii::Point< dim, real > &, const dealii::Tensor< 1, dim, real > &, const std::array< real, nstate > &, const std::array< dealii::Tensor< 1, dim, real >, nstate > &, std::array< real, nstate > &, std::array< dealii::Tensor< 1, dim, real >, nstate > &) const |
| If diffusion is present, assign Dirichlet boundary condition. More... | |
 Public Member Functions inherited from PHiLiP::Physics::PhysicsBase< dim, nspecies, nstate, real > | |
| PhysicsBase (const Parameters::AllParameters *const parameters_input, const bool has_nonzero_diffusion_input, const bool has_nonzero_physical_source_input, const dealii::Tensor< 2, 3, double > input_diffusion_tensor=Parameters::ManufacturedSolutionParam::get_default_diffusion_tensor(), std::shared_ptr< ManufacturedSolutionFunction< dim, nspecies, real > > manufactured_solution_function_input=nullptr) | |
| Default constructor that will set the constants. | |
| PhysicsBase (const Parameters::AllParameters *const parameters_input, const bool has_nonzero_diffusion_input, const bool has_nonzero_physical_source_input, std::shared_ptr< ManufacturedSolutionFunction< dim, nspecies, real > > manufactured_solution_function_input=nullptr) | |
| Constructor that will call default constructor. | |
| virtual | ~PhysicsBase ()=default |
| Virtual destructor required for abstract classes. | |
| virtual real | compute_pressure (const std::array< real, nstate > &conservative_soln) const |
| Compute pressure from conservative solution. | |
| virtual real | compute_entropy (const std::array< real, nstate > &conservative_soln) const |
| Compute entropy from conservative solution. | |
| virtual real | compute_gamma (const std::array< real, nstate > &conservative_soln) const |
| Compute gamma from conservative solution. | |
| virtual std::array< real, nstate > | compute_kinetic_energy_variables (const std::array< real, nstate > &conservative_soln) const |
| Compute kinetic energy variables from conservative solution. | |
| virtual 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) | |
| virtual std::array< real, nstate > | dissipative_flux_dot_normal (const std::array< real, nstate > &solution, const std::array< dealii::Tensor< 1, dim, real >, nstate > &solution_gradient, const std::array< real, nstate > &filtered_solution, const std::array< dealii::Tensor< 1, dim, real >, nstate > &filtered_solution_gradient, const bool on_boundary, const dealii::types::global_dof_index cell_index, const dealii::Tensor< 1, dim, real > &normal, const int boundary_type) |
| Dissipative fluxes dot normal vector. | |
| virtual std::array< dealii::Tensor< 1, dim, real >, nstate > | dissipative_flux (const std::array< real, nstate > &solution, const std::array< dealii::Tensor< 1, dim, real >, nstate > &solution_gradient, const std::array< real, nstate > &filtered_solution, const std::array< dealii::Tensor< 1, dim, real >, nstate > &filtered_solution_gradient, const dealii::types::global_dof_index cell_index) |
| Dissipative fluxes that will be differentiated ONCE in space. | |
| 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 term that does require differentiation. | |
| virtual std::array< real, nstate > | artificial_source_term (const real viscosity_coefficient, const dealii::Point< dim, real > &pos, const std::array< real, nstate > &solution) const |
| Artificial source term that does not require differentiation stemming from artificial dissipation. | |
| virtual void | boundary_face_values (const int, const dealii::Point< dim, real > &, const dealii::Tensor< 1, dim, real > &, const std::array< real, nstate > &, const std::array< dealii::Tensor< 1, dim, real >, nstate > &, const std::array< real, nstate > &, const std::array< dealii::Tensor< 1, dim, real >, nstate > &, std::array< real, nstate > &, std::array< dealii::Tensor< 1, dim, real >, nstate > &) const |
| Evaluates boundary values and gradients on the other side of the face for the convective flux. | |
| virtual void | boundary_face_values_viscous_flux (const int, const dealii::Point< dim, real > &, const dealii::Tensor< 1, dim, real > &, const std::array< real, nstate > &, const std::array< dealii::Tensor< 1, dim, real >, nstate > &, const std::array< real, nstate > &, const std::array< dealii::Tensor< 1, dim, real >, nstate > &, std::array< real, nstate > &, std::array< dealii::Tensor< 1, dim, real >, nstate > &) const |
| Evaluates boundary values and gradients on the other side of the face for the viscous flux. | |
| 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 dealii::Vector< double > | post_compute_derived_quantities_scalar (const double &uh, const dealii::Tensor< 1, dim > &, const dealii::Tensor< 2, dim > &, const dealii::Tensor< 1, dim > &, const dealii::Point< dim > &) const |
| Returns current scalar 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... | |
| 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 dealii::UpdateFlags | post_get_needed_update_flags () const |
| Returns required update flags of the solution to be used by PhysicsPostprocessor to output current solution. More... | |
| template<typename real2 > | |
| real2 | handle_non_physical_result (const std::string message="") const |
| Function to handle nonphysical results. More... | |
Public Attributes | |
| const bool | hasConvection |
| Turns ON/OFF convection term. | |
| const bool | hasDiffusion |
| const Parameters::AllParameters::TestType | test_type |
| Allows convection diffusion to distinguish between different unsteady test types. More... | |
| Public Attributes inherited from PHiLiP::Physics::PhysicsBase< dim, nspecies, nstate, real > | |
| const bool | has_nonzero_diffusion |
| Flag to signal that diffusion term is non-zero. | |
| const bool | has_nonzero_physical_source |
| Flag to signal that physical source term is non-zero. | |
| const Parameters::AllParameters *const | all_parameters |
| Pointer to parameters object. | |
| const NonPhysicalBehaviorEnum | non_physical_behavior_type |
| Determines type of nonphysical behavior. | |
| std::shared_ptr< ManufacturedSolutionFunction< dim, nspecies, real > > | manufactured_solution_function |
| Manufactured solution function. | |
| const double | BIG_NUMBER = 1e100 |
| BIG_NUMBER which is returned in place of NaN according to handle_non_physical_result() More... | |
Protected Member Functions | |
| dealii::Tensor< 1, dim, real > | advection_speed () const |
| Linear advection speed: c. | |
| real | diffusion_coefficient () const |
| Diffusion coefficient. | |
Protected Attributes | |
| double | linear_advection_velocity [3] = { 1.1, -atan(1)*4.0 / exp(1), exp(1)/(atan(1)*4.0) } |
| Linear advection velocity in x, y, and z directions. | |
| double | diffusion_scaling_coeff = 0.1*atan(1)*4.0/exp(1) |
| Diffusion scaling coefficient in front of the diffusion tensor. | |
| Protected Attributes inherited from PHiLiP::Physics::PhysicsBase< dim, nspecies, nstate, real > | |
| dealii::ConditionalOStream | pcout |
| ConditionalOStream. More... | |
| dealii::Tensor< 2, dim, double > | diffusion_tensor |
| Anisotropic diffusion matrix. More... | |
Additional Inherited Members | |
| Public Types inherited from PHiLiP::Physics::PhysicsBase< dim, nspecies, nstate, real > | |
| using | NonPhysicalBehaviorEnum = Parameters::AllParameters::NonPhysicalBehaviorEnum |
Detailed Description
template<int dim, int nspecies, int nstate, typename real>
class PHiLiP::Physics::ConvectionDiffusion< dim, nspecies, nstate, real >
Convection-diffusion with linear advective and diffusive term. Derived from PhysicsBase.
State variable: \( u \)
Convective flux \( \mathbf{F}_{conv} = u \)
Dissipative flux \( \mathbf{F}_{diss} = -\boldsymbol\nabla u \)
Source term \( s(\mathbf{x}) \)
Equation:
\[ \boldsymbol{\nabla} \cdot ( \mathbf{F}_{conv}( u ) + \mathbf{F}_{diss}( u, \boldsymbol{\nabla}(u) ) = s(\mathbf{x}) \]
Definition at line 29 of file convection_diffusion.h.
Member Function Documentation
◆ boundary_face_values()
template<int dim, int nspecies, int nstate, typename real>
|
virtual |
If diffusion is present, assign Dirichlet boundary condition.
Using Neumann boundary conditions might need to modify the functional in order to obtain the optimal 2p convergence of the functional error
Implements PHiLiP::Physics::PhysicsBase< dim, nspecies, nstate, real >.
Definition at line 19 of file convection_diffusion.cpp.
◆ convert_conservative_gradient_to_primitive_gradient()
template<int dim, int nspecies, int nstate, typename real>
|
virtual |
Obtain gradient of primitive variables from gradient of conservative variables
Implements PHiLiP::Physics::PhysicsBase< dim, nspecies, nstate, real >.
Definition at line 118 of file convection_diffusion.cpp.
◆ convert_primitive_gradient_to_conservative_gradient()
template<int dim, int nspecies, int nstate, typename real>
|
virtual |
Obtain gradient of conservative variables from gradient of primitive variables
Implements PHiLiP::Physics::PhysicsBase< dim, nspecies, nstate, real >.
Definition at line 109 of file convection_diffusion.cpp.
Member Data Documentation
◆ hasDiffusion
template<int dim, int nspecies, int nstate, typename real>
| const bool PHiLiP::Physics::ConvectionDiffusion< dim, nspecies, nstate, real >::hasDiffusion |
Turns ON/OFF diffusion term.
Definition at line 50 of file convection_diffusion.h.
◆ test_type
template<int dim, int nspecies, int nstate, typename real>
| const Parameters::AllParameters::TestType PHiLiP::Physics::ConvectionDiffusion< dim, nspecies, nstate, real >::test_type |
Allows convection diffusion to distinguish between different unsteady test types.
Pointer to all parameters
Definition at line 52 of file convection_diffusion.h.
The documentation for this class was generated from the following files:
- src/physics/convection_diffusion.h
- src/physics/convection_diffusion.cpp
