Public Types |
Public Member Functions |
Public Attributes |
Protected Member Functions |
List of all members
PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real > Class Template Reference
Navier-Stokes equations. Derived from Euler for the convective terms, which is derived from PhysicsBase. More...
#include <navier_stokes.h>
Inheritance diagram for PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >:
Collaboration diagram for PHiLiP::Physics::NavierStokes< dim, nspecies, 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::Euler< dim, nspecies, nstate, real > | |
| using | two_point_num_flux_enum = Parameters::AllParameters::TwoPointNumericalFlux |
| Public Types inherited from PHiLiP::Physics::PhysicsBase< dim, nspecies, nstate, real > | |
| using | NonPhysicalBehaviorEnum = Parameters::AllParameters::NonPhysicalBehaviorEnum |
Public Member Functions | |
| NavierStokes (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=273.15, 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, nspecies, real > > manufactured_solution_function=nullptr, const two_point_num_flux_enum two_point_num_flux_type=two_point_num_flux_enum::KG, const bool has_nonzero_physical_source=false) | |
| Constructor. | |
| template<typename real2 > | |
| dealii::Tensor< 1, dim, real2 > | compute_temperature_gradient (const std::array< real2, nstate > &primitive_soln, const std::array< dealii::Tensor< 1, dim, real2 >, nstate > &primitive_soln_gradient) const |
| template<typename real2 > | |
| dealii::Tensor< 1, dim, real2 > | compute_velocities_parallel_to_wall (const std::array< real2, nstate > &conservative_soln, const dealii::Tensor< 1, dim, real2 > &normal_vector) const |
| template<typename real2 > | |
| dealii::Tensor< 1, dim, real2 > | compute_wall_tangent_vector (const std::array< real2, nstate > &conservative_soln, const dealii::Tensor< 1, dim, real2 > &normal_vector) const |
| template<typename real2 > | |
| dealii::Tensor< 1, dim, real2 > | compute_wall_tangent_vector_from_velocities_parallel_to_wall (const dealii::Tensor< 1, dim, real2 > &velocities_parallel_to_wall) const |
| template<typename real2 > | |
| real2 | compute_wall_shear_stress (const std::array< real2, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real2 >, nstate > &conservative_soln_gradient, const dealii::Tensor< 1, dim, real2 > &normal_vector) const |
| template<typename real2 > | |
| real2 | compute_viscosity_coefficient (const std::array< real2, nstate > &primitive_soln) const |
| template<typename real2 > | |
| real2 | compute_viscosity_coefficient_from_temperature (const real2 temperature) const |
| template<typename real2 > | |
| real2 | compute_viscosity_coefficient_sutherlands_law (const std::array< real2, nstate > &primitive_soln) const |
| template<typename real2 > | |
| real2 | compute_viscosity_coefficient_sutherlands_law_from_temperature (const real2 temperature) const |
| template<typename real2 > | |
| real2 | scale_viscosity_coefficient (const real2 viscosity_coefficient) const |
| template<typename real2 > | |
| real2 | compute_scaled_viscosity_coefficient (const std::array< real2, nstate > &primitive_soln) const |
| template<typename real2 > | |
| real2 | compute_scaled_heat_conductivity_given_scaled_viscosity_coefficient_and_prandtl_number (const real2 scaled_viscosity_coefficient, const double prandtl_number_input) const |
| template<typename real2 > | |
| real2 | compute_scaled_heat_conductivity (const std::array< real2, nstate > &primitive_soln) const |
| template<typename real2 > | |
| dealii::Tensor< 1, dim, real2 > | compute_heat_flux (const std::array< real2, nstate > &primitive_soln, const std::array< dealii::Tensor< 1, dim, real2 >, nstate > &primitive_soln_gradient) const |
| template<typename real2 > | |
| dealii::Tensor< 1, dim, real2 > | compute_heat_flux_given_scaled_heat_conductivity_and_temperature_gradient (const real2 scaled_heat_conductivity, const dealii::Tensor< 1, dim, real2 > &temperature_gradient) const |
| template<typename real2 > | |
| dealii::Tensor< 1, 3, real2 > | compute_vorticity (const std::array< real2, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real2 >, nstate > &conservative_soln_gradient) const |
| Evaluate vorticity from conservative variables and gradient of conservative variables. | |
| real | compute_vorticity_magnitude_sqr (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &conservative_soln_gradient) const |
| Evaluate vorticity magnitude squared from conservative variables and gradient of conservative variables. | |
| real | compute_vorticity_magnitude (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &conservative_soln_gradient) const |
| Evaluate vorticity magnitude from conservative variables and gradient of conservative variables. | |
| real | compute_second_invariant (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &conservative_soln_gradient) const |
| real | compute_enstrophy (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &conservative_soln_gradient) const |
| Evaluate enstrophy from conservative variables and gradient of conservative variables. | |
| real | compute_incompressible_enstrophy (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &conservative_soln_gradient) const |
| Evaluate incompressible enstrophy from conservative variables and gradient of conservative variables. | |
| real | compute_incompressible_palinstrophy (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, 3 > &vorticity_gradient) const |
| Evaluate incompressible palinstrophy from conservative variables and gradient of vorticity. | |
| real | compute_vorticity_based_dissipation_rate_from_integrated_enstrophy (const real integrated_enstrophy) const |
| real | compute_pressure_dilatation (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &conservative_soln_gradient) const |
| real | compute_dilatation (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &conservative_soln_gradient) const |
| real | compute_density_gradient_magnitude (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &conservative_soln_gradient) const |
| dealii::Tensor< 2, dim, real > | compute_deviatoric_strain_rate_tensor (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &conservative_soln_gradient) const |
| dealii::Tensor< 2, dim, real > | compute_strain_rate_tensor_from_conservative (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &conservative_soln_gradient) const |
| template<typename real2 > | |
| dealii::Tensor< 2, dim, real2 > | compute_strain_rate_tensor_from_conservative_templated (const std::array< real2, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real2 >, nstate > &conservative_soln_gradient) const |
| real | compute_viscosity_times_deviatoric_strain_rate_tensor_magnitude_sqr (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &conservative_soln_gradient) const |
| Evaluate the square of the deviatoric strain-rate tensor magnitude (i.e. double dot product) from conservative variables and gradient of conservative variables. | |
| real | compute_deviatoric_strain_rate_tensor_based_dissipation_rate_from_integrated_viscosity_times_deviatoric_strain_rate_tensor_magnitude_sqr (const real integrated_viscosity_times_deviatoric_strain_rate_tensor_magnitude_sqr) const |
| template<typename real2 > | |
| dealii::Tensor< 2, dim, real2 > | extract_velocities_gradient_from_primitive_solution_gradient (const std::array< dealii::Tensor< 1, dim, real2 >, nstate > &primitive_soln_gradient) const |
| template<typename real2 > | |
| dealii::Tensor< 2, dim, real2 > | compute_strain_rate_tensor (const dealii::Tensor< 2, dim, real2 > &vel_gradient) const |
| real | compute_viscosity_times_strain_rate_tensor_magnitude_sqr (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &conservative_soln_gradient) const |
| Evaluate the square of the strain-rate tensor magnitude (i.e. double dot product) from conservative variables and gradient of conservative variables. | |
| real | compute_strain_rate_tensor_based_dissipation_rate_from_integrated_viscosity_times_strain_rate_tensor_magnitude_sqr (const real integrated_viscosity_times_strain_rate_tensor_magnitude_sqr) const |
| template<typename real2 > | |
| dealii::Tensor< 2, dim, real2 > | compute_viscous_stress_tensor_via_scaled_viscosity_and_strain_rate_tensor (const real2 scaled_viscosity_coefficient, const dealii::Tensor< 2, dim, real2 > &strain_rate_tensor) const |
| dealii::Tensor< 2, dim, real > | compute_viscous_stress_tensor_from_conservative (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &conservative_soln_gradient) const |
| template<typename real2 > | |
| dealii::Tensor< 2, dim, real2 > | compute_viscous_stress_tensor_from_conservative_templated (const std::array< real2, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real2 >, nstate > &conservative_soln_gradient) const |
| template<typename real2 > | |
| dealii::Tensor< 2, dim, real2 > | compute_viscous_stress_tensor (const std::array< real2, nstate > &primitive_soln, const std::array< dealii::Tensor< 1, dim, real2 >, nstate > &primitive_soln_gradient) const |
| real | get_tensor_product_magnitude_sqr (const dealii::Tensor< 2, dim, real > &tensor1, const dealii::Tensor< 2, dim, real > &tensor2) const |
| Tensor product magnitude squared. | |
| dealii::Tensor< 2, dim, real > | compute_germano_idendity_matrix_L_component (const std::array< real, nstate > &conservative_soln) const |
| dealii::Tensor< 2, dim, real > | compute_germano_identity_matrix_M_component (const std::array< real, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real >, nstate > &conservative_soln_gradient) const |
| 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) override |
| virtual std::array< real, nstate > | dissipative_flux_dot_normal_on_adiabatic_boundary (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, const dealii::Tensor< 1, dim, real > &normal) |
| 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 override |
| dealii::Tensor< 1, dim, real > | compute_scaled_viscosity_gradient (const std::array< real, nstate > &primitive_soln, const dealii::Tensor< 1, dim, real > &temperature_gradient) const |
| 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::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 |
| std::array< real, nstate > | dissipative_source_term (const dealii::Point< dim, real > &pos) const |
| Dissipative flux contribution to the source term. | |
| std::array< real, nstate > | source_term (const dealii::Point< dim, real > &pos, const std::array< real, nstate > &conservative_soln, const real current_time) const override |
| Source term is zero or depends on manufactured solution. | |
| dealii::Tensor< 2, nstate, real > | convective_flux_directional_jacobian_via_dfad (std::array< real, nstate > &conservative_soln, const dealii::Tensor< 1, dim, real > &normal) const |
| real | compute_scaled_viscosity_coefficient_derivative_wrt_temperature_via_dfad (std::array< real, nstate > &conservative_soln) const |
| template<typename real2 > | |
| std::array< dealii::Tensor< 1, dim, real2 >, nstate > | dissipative_flux_given_velocities_viscous_stress_tensor_and_heat_flux (const dealii::Tensor< 1, dim, real2 > &vel, const dealii::Tensor< 2, dim, real2 > &viscous_stress_tensor, const dealii::Tensor< 1, dim, real2 > &heat_flux) const |
| dealii::Vector< double > | post_compute_derived_quantities_vector (const dealii::Vector< double > &uh, const std::vector< dealii::Tensor< 1, dim > > &duh, const std::vector< dealii::Tensor< 2, dim > > &dduh, const dealii::Tensor< 1, dim > &normals, const dealii::Point< dim > &evaluation_points) const override |
| For post processing purposes, computes all the quantities we write to the VTK files. | |
| std::vector< std::string > | post_get_names () const override |
| For post processing purposes, sets the base names (with no prefix or suffix) of the computed quantities. | |
| std::vector< dealii::DataComponentInterpretation::DataComponentInterpretation > | post_get_data_component_interpretation () const override |
| For post processing purposes, sets the interpretation of each computed quantity as either scalar or vector. | |
| dealii::UpdateFlags | post_get_needed_update_flags () const override |
| For post processing purposes, updates the required flags for dealii. | |
| real | get_tensor_magnitude_sqr (const dealii::Tensor< 2, dim, real > &tensor) const |
| Returns the square of the magnitude of the tensor (i.e. the double dot product of a tensor with itself) | |
| real | get_tensor_magnitude (const dealii::Tensor< 2, dim, real > &tensor) const |
| Returns the the magnitude of the tensor (i.e. the double dot product of a tensor with itself) | |
| Public Member Functions inherited from PHiLiP::Physics::Euler< dim, nspecies, nstate, real > | |
| Euler (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, std::shared_ptr< ManufacturedSolutionFunction< dim, nspecies, real > > manufactured_solution_function=nullptr, const two_point_num_flux_enum two_point_num_flux_type=two_point_num_flux_enum::KG, const bool has_nonzero_diffusion=false, const bool has_nonzero_physical_source=false) | |
| Constructor. | |
| std::array< dealii::Tensor< 1, dim, real >, nstate > | convective_flux (const std::array< real, nstate > &conservative_soln) const override |
| Convective flux: \( \mathbf{F}_{conv} \). | |
| std::array< real, nstate > | convective_normal_flux (const std::array< real, nstate > &conservative_soln, const dealii::Tensor< 1, dim, real > &normal) const |
| Convective normal flux: \( \mathbf{F}_{conv} \cdot \hat{n} \). | |
| dealii::Tensor< 2, nstate, real > | convective_flux_directional_jacobian (const std::array< real, nstate > &conservative_soln, const dealii::Tensor< 1, dim, real > &normal) const |
| Convective flux Jacobian: \( \frac{\partial \mathbf{F}_{conv}}{\partial w} \cdot \mathbf{n} \). | |
| std::array< real, nstate > | convective_eigenvalues (const std::array< real, nstate > &, const dealii::Tensor< 1, dim, real > &) const override |
| Spectral radius of convective term Jacobian is 'c'. | |
| real | max_convective_eigenvalue (const std::array< real, nstate > &soln) const override |
| Maximum convective eigenvalue. | |
| real | max_convective_normal_eigenvalue (const std::array< real, nstate > &soln, const dealii::Tensor< 1, dim, real > &normal) const override |
| Maximum convective normal eigenvalue (used in Lax-Friedrichs) More... | |
| real | max_viscous_eigenvalue (const std::array< real, nstate > &soln) const override |
| Maximum viscous eigenvalue. | |
| 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 flux: 0. | |
| std::array< real, nstate > | source_term (const dealii::Point< dim, real > &pos, const std::array< real, nstate > &conservative_soln, 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 > | convective_source_term (const dealii::Point< dim, real > &pos) const |
| Convective flux contribution to the source term. | |
| template<typename real2 > | |
| std::array< real2, nstate > | convert_conservative_to_primitive_templated (const std::array< real2, nstate > &conservative_soln) const |
| std::array< real, nstate > | convert_conservative_to_primitive (const std::array< real, nstate > &conservative_soln) const |
| Convert conservative to primitive (real2==real); required by base class. | |
| std::array< real, nstate > | convert_primitive_to_conservative (const std::array< real, nstate > &primitive_soln) const |
| template<typename real2 > | |
| std::array< dealii::Tensor< 1, dim, real2 >, nstate > | convert_conservative_gradient_to_primitive_gradient_templated (const std::array< real2, nstate > &conservative_soln, const std::array< dealii::Tensor< 1, dim, real2 >, nstate > &conservative_soln_gradient) const |
| 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 |
| virtual real | compute_gamma (const std::array< real, nstate > &conservative_soln) const |
| For Euler it returns this->gam. | |
| template<typename real2 > | |
| real2 | compute_pressure_templated (const std::array< real2, nstate > &conservative_soln) const |
| Evaluate pressure from conservative variables. | |
| real | compute_pressure (const std::array< real, nstate > &conservative_soln) const |
| Compute pressure from conservative solution. | |
| template<typename real2 > | |
| real2 | compute_entropy (const real2 density, const real2 pressure) const |
| Evaluate physical entropy = log(p ^{-}) from pressure and density. | |
| real | compute_specific_enthalpy (const std::array< real, nstate > &conservative_soln, const real pressure) const |
| Evaluate pressure from conservative variables. | |
| template<typename real2 > | |
| real2 | compute_entropy_templated (const std::array< real2, nstate > &conservative_soln) const |
| Evaluate physical entropy = log(p ^{-}) from pressure and density. | |
| real | compute_entropy (const std::array< real, nstate > &conservative_soln) const |
| Evaluate physical entropy from conservative variables. | |
| real | compute_numerical_entropy_function (const std::array< real, nstate > &conservative_soln) const |
| Compute numerical entropy function -rho s. | |
| real | compute_sound (const std::array< real, nstate > &conservative_soln) const |
| Evaluate speed of sound from conservative variables. | |
| real | compute_sound (const real density, const real pressure) const |
| Evaluate speed of sound from density and pressure. | |
| template<typename real2 > | |
| dealii::Tensor< 1, dim, real2 > | compute_velocities (const std::array< real2, nstate > &conservative_soln) const |
| Evaluate velocities from conservative variables. | |
| template<typename real2 > | |
| real2 | compute_velocity_squared (const dealii::Tensor< 1, dim, real2 > &velocities) const |
| Given the velocity vector \( \mathbf{u} \), returns the dot-product \( \mathbf{u} \cdot \mathbf{u} \). | |
| template<typename real2 > | |
| dealii::Tensor< 1, dim, real2 > | extract_velocities_from_primitive (const std::array< real2, nstate > &primitive_soln) const |
| Given primitive variables, returns velocities. | |
| real | compute_total_energy (const std::array< real, nstate > &primitive_soln) const |
| Given primitive variables, returns total energy. More... | |
| real | compute_kinetic_energy_from_primitive_solution (const std::array< real, nstate > &primitive_soln) const |
| Given primitive variables, returns kinetic energy. | |
| real | compute_incompressible_kinetic_energy_from_primitive_solution (const std::array< real, nstate > &primitive_soln) const |
| Given primitive variables, returns incompressible kinetic energy. | |
| real | compute_kinetic_energy_from_conservative_solution (const std::array< real, nstate > &conservative_soln) const |
| Given conservative variables, returns kinetic energy. | |
| real | compute_incompressible_kinetic_energy_from_conservative_solution (const std::array< real, nstate > &conservative_soln) const |
| Given conservative variables, returns incompressible kinetic energy. | |
| real | compute_entropy_measure (const std::array< real, nstate > &conservative_soln) const |
| Evaluate entropy from conservative variables. More... | |
| real | compute_entropy_measure (const real density, const real pressure) const |
| Evaluate entropy from density and pressure. | |
| real | compute_mach_number (const std::array< real, nstate > &conservative_soln) const |
| Given conservative variables, returns Mach number. | |
| template<typename real2 > | |
| real2 | compute_temperature (const std::array< real2, nstate > &primitive_soln) const |
| Given primitive variables, returns NON-DIMENSIONALIZED temperature using free-stream non-dimensionalization. More... | |
| real | compute_density_from_pressure_temperature (const real pressure, const real temperature) const |
| Given pressure and temperature, returns NON-DIMENSIONALIZED density using free-stream non-dimensionalization. More... | |
| real | compute_temperature_from_density_pressure (const real density, const real pressure) const |
| Given density and pressure, returns NON-DIMENSIONALIZED temperature using free-stream non-dimensionalization. More... | |
| real | compute_pressure_from_density_temperature (const real density, const real temperature) const |
| Given density and temperature, returns NON-DIMENSIONALIZED pressure using free-stream non-dimensionalization. More... | |
| std::array< dealii::Tensor< 1, dim, real >, nstate > | convective_numerical_split_flux (const std::array< real, nstate > &conservative_soln1, const std::array< real, nstate > &conservative_soln2) const override |
| Evaluates convective flux based on the chosen split form. | |
| std::array< real, nstate > | compute_entropy_variables (const std::array< real, nstate > &conservative_soln) const |
| std::array< real, nstate > | compute_conservative_variables_from_entropy_variables (const std::array< real, nstate > &entropy_var) const |
| std::array< real, nstate > | compute_kinetic_energy_variables (const std::array< real, nstate > &conservative_soln) const |
| Computes the kinetic energy variables. | |
| real | compute_mean_density (const std::array< real, nstate > &conservative_soln1, const std::array< real, nstate > &convervative_soln2) const |
| Mean density given two sets of conservative solutions. More... | |
| real | compute_mean_pressure (const std::array< real, nstate > &conservative_soln1, const std::array< real, nstate > &convervative_soln2) const |
| Mean pressure given two sets of conservative solutions. More... | |
| dealii::Tensor< 1, dim, real > | compute_mean_velocities (const std::array< real, nstate > &conservative_soln1, const std::array< real, nstate > &convervative_soln2) const |
| Mean velocities given two sets of conservative solutions. More... | |
| real | compute_mean_specific_total_energy (const std::array< real, nstate > &conservative_soln1, const std::array< real, nstate > &convervative_soln2) const |
| Mean specific total energy given two sets of conservative solutions. More... | |
| 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 |
| Boundary condition handler. | |
| 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 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 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... | |
| template<typename real2 > | |
| real2 | handle_non_physical_result (const std::string message="") const |
| Function to handle nonphysical results. More... | |
Public Attributes | |
| const double | viscosity_coefficient_inf |
| Nondimensionalized viscosity coefficient at infinity. | |
| const bool | use_constant_viscosity |
| Flag to use constant viscosity instead of Sutherland's law of viscosity. | |
| const double | constant_viscosity |
| Nondimensionalized constant viscosity. | |
| const double | prandtl_number |
| Prandtl number. | |
| const double | reynolds_number_inf |
| Farfield (free stream) Reynolds number. | |
| const double | isothermal_wall_temperature |
| Nondimensionalized isothermal wall temperature. | |
| const thermal_boundary_condition_enum | thermal_boundary_condition_type |
| Thermal boundary condition type (adiabatic or isothermal) | |
| Public Attributes inherited from PHiLiP::Physics::Euler< dim, nspecies, nstate, real > | |
| const double | ref_length |
| Reference length. | |
| const double | gam |
| Constant heat capacity ratio of fluid. | |
| const double | gamm1 |
| Constant heat capacity ratio (Gamma-1.0) used often. | |
| const double | density_inf |
| const double | mach_inf |
| Farfield Mach number. | |
| const double | mach_inf_sqr |
| const double | angle_of_attack |
| Angle of attack. More... | |
| const double | side_slip_angle |
| Sideslip angle. More... | |
| const double | sound_inf |
| Non-dimensionalized sound* at infinity. | |
| const double | pressure_inf |
| Non-dimensionalized pressure* at infinity. | |
| const double | entropy_inf |
| Entropy measure at infinity. | |
| const two_point_num_flux_enum | two_point_num_flux_type |
| Two point numerical flux type (for split form) | |
| double | temperature_inf |
| Non-dimensionalized temperature* at infinity. Should equal 1/density*(inf) | |
| double | dynamic_pressure_inf |
| Non-dimensionalized dynamic pressure* at infinity. | |
| dealii::Tensor< 1, dim, double > | velocities_inf |
| Non-dimensionalized Velocity vector at farfield. 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 | |
| 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 |
| virtual void | boundary_face_values_viscous_flux (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, const std::array< real, nstate > &, const std::array< dealii::Tensor< 1, dim, real >, nstate > &, std::array< real, nstate > &soln_bc, std::array< dealii::Tensor< 1, dim, real >, nstate > &soln_grad_bc) const override |
| Boundary face values for viscous fluxes. | |
| void | boundary_wall_viscous_flux (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 |
| 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 override |
| Evaluate the manufactured solution boundary conditions. | |
| Protected Member Functions inherited from PHiLiP::Physics::Euler< dim, nspecies, nstate, real > | |
| template<typename real2 > | |
| bool | check_positive_quantity (real2 &quantity, const std::string qty_name) const |
| Check positive quantity and modify it according to handle_non_physical_result() More... | |
| 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 |
| void | boundary_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 |
| Wall boundary condition. | |
| void | boundary_pressure_outflow (const real total_inlet_pressure, const real back_pressure, const std::array< real, nstate > &soln_int, std::array< real, nstate > &soln_bc) const |
| void | boundary_inflow (const real total_inlet_pressure, const real total_inlet_temperature, const dealii::Tensor< 1, dim, real > &normal_int, const std::array< real, nstate > &soln_int, std::array< real, nstate > &soln_bc) const |
| 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 |
| void | boundary_farfield (std::array< real, nstate > &soln_bc) const |
| Simple farfield boundary conditions based on freestream values. | |
| void | boundary_p0_extrapolation (const std::array< real, nstate > &soln_int, std::array< real, nstate > &soln_bc, std::array< dealii::Tensor< 1, dim, real >, nstate > &soln_grad_bc) const |
| p0 extrapolation at the boundary | |
| void | boundary_custom (std::array< real, nstate > &soln_bc) const |
| Custom boundary conditions for the left boundary of the astrophysical mach jet case where it is not hypersonic inflow. | |
| void | boundary_astrophysical_inflow (std::array< real, nstate > &soln_bc) const |
| Boundary conditions based on user-defined values. | |
| std::array< real, nstate > | get_manufactured_solution_value (const dealii::Point< dim, real > &pos) const |
| Get manufactured solution value. | |
| std::array< dealii::Tensor< 1, dim, real >, nstate > | get_manufactured_solution_gradient (const dealii::Point< dim, real > &pos) const |
| Get manufactured solution gradient. | |
| std::array< dealii::Tensor< 1, dim, real >, nstate > | convective_numerical_split_flux_kennedy_gruber (const std::array< real, nstate > &conservative_soln1, const std::array< real, nstate > &conservative_soln2) const |
| std::array< real, nstate > | compute_ismail_roe_parameter_vector_from_primitive (const std::array< real, nstate > &primitive_soln) const |
| Compute Ismail-Roe parameter vector from primitive solution. | |
| real | compute_ismail_roe_logarithmic_mean (const real val1, const real val2) const |
| Compute Ismail-Roe logarithmic mean. | |
| std::array< dealii::Tensor< 1, dim, real >, nstate > | convective_numerical_split_flux_ismail_roe (const std::array< real, nstate > &conservative_soln1, const std::array< real, nstate > &conservative_soln2) const |
| std::array< dealii::Tensor< 1, dim, real >, nstate > | convective_numerical_split_flux_chandrashekar (const std::array< real, nstate > &conservative_soln1, const std::array< real, nstate > &conservative_soln2) const |
| Chandrashekar entropy conserving flux. | |
| std::array< dealii::Tensor< 1, dim, real >, nstate > | convective_numerical_split_flux_ranocha (const std::array< real, nstate > &conservative_soln1, const std::array< real, nstate > &conservative_soln2) const |
| Ranocha pressure equilibrium preserving, entropy and energy conserving flux. | |
Protected Attributes | |
| const double | sutherlands_temperature |
| Sutherland's temperature. Units: [K]. More... | |
| const double | freestream_temperature |
| Freestream temperature. Units: [K]. | |
| const double | temperature_ratio |
| Ratio of Sutherland's temperature to freestream temperature. | |
| 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... | |
Detailed Description
template<int dim, int nspecies, int nstate, typename real>
class PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >
Navier-Stokes equations. Derived from Euler for the convective terms, which is derived from PhysicsBase.
Definition at line 12 of file navier_stokes.h.
Member Function Documentation
◆ boundary_wall_viscous_flux()
template<int dim, int nspecies, int nstate, typename real>
|
protected |
No-slip wall boundary conditions
- Given by equations 460-461 of the following paper:
- * Hartmann, Ralf. "Numerical analysis of higher order discontinuous Galerkin finite element methods." (2008): 1-107.
Definition at line 1404 of file navier_stokes.cpp.
◆ compute_density_gradient_magnitude()
template<int dim, int nspecies, int nstate, typename real>
| real PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_density_gradient_magnitude | ( | const std::array< real, nstate > & | conservative_soln, |
| const std::array< dealii::Tensor< 1, dim, real >, nstate > & | conservative_soln_gradient | ||
| ) | const |
Evaluate density gradient magnitude from conservative variables and gradient of conservative variables
Definition at line 567 of file navier_stokes.cpp.
◆ compute_deviatoric_strain_rate_tensor()
template<int dim, int nspecies, int nstate, typename real>
| dealii::Tensor< 2, dim, real > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_deviatoric_strain_rate_tensor | ( | const std::array< real, nstate > & | conservative_soln, |
| const std::array< dealii::Tensor< 1, dim, real >, nstate > & | conservative_soln_gradient | ||
| ) | const |
Evaluate the deviatoric strain-rate tensor from conservative variables and gradient of conservative variables – Reference: de la Llave Plata et al. (2019). "On the performance of a high-order multiscale DG approach to LES at increasing Reynolds number."
Definition at line 612 of file navier_stokes.cpp.
◆ compute_deviatoric_strain_rate_tensor_based_dissipation_rate_from_integrated_viscosity_times_deviatoric_strain_rate_tensor_magnitude_sqr()
template<int dim, int nspecies, int nstate, typename real>
| real PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_deviatoric_strain_rate_tensor_based_dissipation_rate_from_integrated_viscosity_times_deviatoric_strain_rate_tensor_magnitude_sqr | ( | const real | integrated_viscosity_times_deviatoric_strain_rate_tensor_magnitude_sqr | ) | const |
Evaluate non-dimensional theoretical deviatoric strain-rate tensor based dissipation rate from integrated deviatoric strain-rate tensor magnitude squared. – Reference: Cox, Christopher, et al. "Accuracy, stability, and performance comparison between the spectral difference and flux reconstruction schemes." Computers & Fluids 221 (2021): 104922. – Equation (57a) with free-stream nondimensionalization applied
Definition at line 682 of file navier_stokes.cpp.
◆ compute_dilatation()
template<int dim, int nspecies, int nstate, typename real>
| real PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_dilatation | ( | const std::array< real, nstate > & | conservative_soln, |
| const std::array< dealii::Tensor< 1, dim, real >, nstate > & | conservative_soln_gradient | ||
| ) | const |
Evaluate dilatation from conservative variables and gradient of conservative variables
Definition at line 548 of file navier_stokes.cpp.
◆ compute_germano_idendity_matrix_L_component()
template<int dim, int nspecies, int nstate, typename real>
| dealii::Tensor< 2, dim, real > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_germano_idendity_matrix_L_component | ( | const std::array< real, nstate > & | conservative_soln | ) | const |
Nondimensionalized Germano identity tensor, L*, from conservative solution and solution gradient Reference: Flad and Gassner 2017
Definition at line 837 of file navier_stokes.cpp.
◆ compute_germano_identity_matrix_M_component()
template<int dim, int nspecies, int nstate, typename real>
| dealii::Tensor< 2, dim, real > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_germano_identity_matrix_M_component | ( | const std::array< real, nstate > & | conservative_soln, |
| const std::array< dealii::Tensor< 1, dim, real >, nstate > & | conservative_soln_gradient | ||
| ) | const |
Nondimensionalized Germano identity tensor, M*, from conservative solution and solution gradient Reference: Flad and Gassner 2017
Definition at line 852 of file navier_stokes.cpp.
◆ compute_heat_flux()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
| dealii::Tensor< 1, dim, real2 > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_heat_flux | ( | const std::array< real2, nstate > & | primitive_soln, |
| const std::array< dealii::Tensor< 1, dim, real2 >, nstate > & | primitive_soln_gradient | ||
| ) | const |
Nondimensionalized heat flux, q* Reference: Masatsuka 2018 "I do like CFD", p.148, eq.(4.14.13)
Definition at line 365 of file navier_stokes.cpp.
◆ compute_heat_flux_given_scaled_heat_conductivity_and_temperature_gradient()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
| dealii::Tensor< 1, dim, real2 > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_heat_flux_given_scaled_heat_conductivity_and_temperature_gradient | ( | const real2 | scaled_heat_conductivity, |
| const dealii::Tensor< 1, dim, real2 > & | temperature_gradient | ||
| ) | const |
Nondimensionalized heat flux, q*, given the scaled heat conductivity and temperature gradient Reference: Masatsuka 2018 "I do like CFD", p.148, eq.(4.14.13)
Definition at line 382 of file navier_stokes.cpp.
◆ compute_pressure_dilatation()
template<int dim, int nspecies, int nstate, typename real>
| real PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_pressure_dilatation | ( | const std::array< real, nstate > & | conservative_soln, |
| const std::array< dealii::Tensor< 1, dim, real >, nstate > & | conservative_soln_gradient | ||
| ) | const |
Evaluate pressure dilatation from conservative variables and gradient of conservative variables – Reference: Cox, Christopher, et al. "Accuracy, stability, and performance comparison between the spectral difference and flux reconstruction schemes." Computers & Fluids 221 (2021): 104922.
Definition at line 532 of file navier_stokes.cpp.
◆ compute_scaled_heat_conductivity()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
|
inline |
Scaled nondimensionalized heat conductivity, hat{kappa*} Reference: Masatsuka 2018 "I do like CFD", p.148, eq.(4.14.13)
Definition at line 350 of file navier_stokes.cpp.
◆ compute_scaled_heat_conductivity_given_scaled_viscosity_coefficient_and_prandtl_number()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
|
inline |
Scaled nondimensionalized heat conductivity, hat{kappa*}, given scaled nondimensionalized viscosity coefficient and Prandtl number Reference: Masatsuka 2018 "I do like CFD", p.148, eq.(4.14.13)
Definition at line 337 of file navier_stokes.cpp.
◆ compute_scaled_viscosity_coefficient()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
|
inline |
Scaled nondimensionalized viscosity coefficient, hat{mu*} Reference: Masatsuka 2018 "I do like CFD", p.148, eq.(4.14.14)
Definition at line 323 of file navier_stokes.cpp.
◆ compute_scaled_viscosity_coefficient_derivative_wrt_temperature_via_dfad()
template<int dim, int nspecies, int nstate, typename real>
|
inline |
Scaled viscosity coefficient derivative wrt temperature via dfad (automatic differentiation) – Only used for verifying the basic dfad procedure that is extended in convective_flux_directional_jacobian_via_dfad()
Definition at line 1171 of file navier_stokes.cpp.
◆ compute_scaled_viscosity_gradient()
template<int dim, int nspecies, int nstate, typename real>
| dealii::Tensor< 1, dim, real > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_scaled_viscosity_gradient | ( | const std::array< real, nstate > & | primitive_soln, |
| const dealii::Tensor< 1, dim, real > & | temperature_gradient | ||
| ) | const |
Gradient of the scaled nondimensionalized viscosity coefficient Reference: Masatsuka 2018 "I do like CFD", p.148, eq.(4.14.14 and 4.14.17)
Definition at line 911 of file navier_stokes.cpp.
◆ compute_second_invariant()
template<int dim, int nspecies, int nstate, typename real>
| real PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_second_invariant | ( | const std::array< real, nstate > & | conservative_soln, |
| const std::array< dealii::Tensor< 1, dim, real >, nstate > & | conservative_soln_gradient | ||
| ) | const |
Evaluate second invariant (i.e. Q-criterion) from conservative variables and gradient of conservative variables – Reference: Jeong J, Hussain F. On the identification of a vortex. Journal of Fluid Mechanics. 1995;285:69-94. doi:10.1017/S0022112095000462 – Equation (2)
Definition at line 454 of file navier_stokes.cpp.
◆ compute_strain_rate_tensor()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
| dealii::Tensor< 2, dim, real2 > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_strain_rate_tensor | ( | const dealii::Tensor< 2, dim, real2 > & | vel_gradient | ) | const |
Nondimensionalized strain rate tensor, S* Reference: Masatsuka 2018 "I do like CFD", p.148, extracted from eq.(4.14.12)
Definition at line 707 of file navier_stokes.cpp.
◆ compute_strain_rate_tensor_based_dissipation_rate_from_integrated_viscosity_times_strain_rate_tensor_magnitude_sqr()
template<int dim, int nspecies, int nstate, typename real>
| real PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_strain_rate_tensor_based_dissipation_rate_from_integrated_viscosity_times_strain_rate_tensor_magnitude_sqr | ( | const real | integrated_viscosity_times_strain_rate_tensor_magnitude_sqr | ) | const |
Evaluate non-dimensional theoretical strain-rate tensor based dissipation rate from integrated strain-rate tensor magnitude squared. – Reference: Navah, Farshad, et al. "A High-Order Variational Multiscale Approach to Turbulence for Compact Nodal Schemes." – Equation (E.9) with free-stream nondimensionalization applied
Definition at line 745 of file navier_stokes.cpp.
◆ compute_strain_rate_tensor_from_conservative()
template<int dim, int nspecies, int nstate, typename real>
| dealii::Tensor< 2, dim, real > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_strain_rate_tensor_from_conservative | ( | const std::array< real, nstate > & | conservative_soln, |
| const std::array< dealii::Tensor< 1, dim, real >, nstate > & | conservative_soln_gradient | ||
| ) | const |
Nondimensionalized strain rate tensor, S*, from conservative solution and solution gradient Reference: Masatsuka 2018 "I do like CFD", p.148, extracted from eq.(4.14.12)
Definition at line 587 of file navier_stokes.cpp.
◆ compute_strain_rate_tensor_from_conservative_templated()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
| dealii::Tensor< 2, dim, real2 > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_strain_rate_tensor_from_conservative_templated | ( | const std::array< real2, nstate > & | conservative_soln, |
| const std::array< dealii::Tensor< 1, dim, real2 >, nstate > & | conservative_soln_gradient | ||
| ) | const |
Nondimensionalized strain rate tensor, S*, from conservative solution and solution gradient Reference: Masatsuka 2018 "I do like CFD", p.148, extracted from eq.(4.14.12)
Definition at line 597 of file navier_stokes.cpp.
◆ compute_temperature_gradient()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
| dealii::Tensor< 1, dim, real2 > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_temperature_gradient | ( | const std::array< real2, nstate > & | primitive_soln, |
| const std::array< dealii::Tensor< 1, dim, real2 >, nstate > & | primitive_soln_gradient | ||
| ) | const |
Nondimensionalized temperature gradient
Definition at line 62 of file navier_stokes.cpp.
◆ compute_velocities_parallel_to_wall()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
| dealii::Tensor< 1, dim, real2 > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_velocities_parallel_to_wall | ( | const std::array< real2, nstate > & | conservative_soln, |
| const dealii::Tensor< 1, dim, real2 > & | normal_vector | ||
| ) | const |
Nondimensionalized velocities parallel to wall
Definition at line 79 of file navier_stokes.cpp.
◆ compute_viscosity_coefficient()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
|
inline |
Nondimensionalized viscosity coefficient, mu* Based on the use_constant_viscosity flag, it returns a value based on either: (1) Sutherland's viscosity law, or (2) Constant nondimensionalized viscosity value
Definition at line 241 of file navier_stokes.cpp.
◆ compute_viscosity_coefficient_from_temperature()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
|
inline |
Nondimensionalized viscosity coefficient, mu* Based on the use_constant_viscosity flag, it returns a value based on either: (1) Sutherland's viscosity law, or (2) Constant nondimensionalized viscosity value
Definition at line 257 of file navier_stokes.cpp.
◆ compute_viscosity_coefficient_sutherlands_law()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
|
inline |
Nondimensionalized viscosity coefficient, mu* Reference: Masatsuka 2018 "I do like CFD", p.148, eq.(4.14.16)
Based on Sutherland's law for viscosity
- Reference: Sutherland, W. (1893), "The viscosity of gases and molecular force", Philosophical Magazine, S. 5, 36, pp. 507-531 (1893)
- Values: https://www.cfd-online.com/Wiki/Sutherland%27s_law
Definition at line 273 of file navier_stokes.cpp.
◆ compute_viscosity_coefficient_sutherlands_law_from_temperature()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
|
inline |
Nondimensionalized viscosity coefficient, mu* Reference: Masatsuka 2018 "I do like CFD", p.148, eq.(4.14.16)
Based on Sutherland's law for viscosity
- Reference: Sutherland, W. (1893), "The viscosity of gases and molecular force", Philosophical Magazine, S. 5, 36, pp. 507-531 (1893)
- Values: https://www.cfd-online.com/Wiki/Sutherland%27s_law
Definition at line 293 of file navier_stokes.cpp.
◆ compute_viscous_stress_tensor()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
| dealii::Tensor< 2, dim, real2 > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_viscous_stress_tensor | ( | const std::array< real2, nstate > & | primitive_soln, |
| const std::array< dealii::Tensor< 1, dim, real2 >, nstate > & | primitive_soln_gradient | ||
| ) | const |
Nondimensionalized viscous stress tensor, tau* Reference: Masatsuka 2018 "I do like CFD", p.148, eq.(4.14.12)
Definition at line 817 of file navier_stokes.cpp.
◆ compute_viscous_stress_tensor_from_conservative()
template<int dim, int nspecies, int nstate, typename real>
| dealii::Tensor< 2, dim, real > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_viscous_stress_tensor_from_conservative | ( | const std::array< real, nstate > & | conservative_soln, |
| const std::array< dealii::Tensor< 1, dim, real >, nstate > & | conservative_soln_gradient | ||
| ) | const |
Nondimensionalized viscous stress tensor, tau* Reference: Masatsuka 2018 "I do like CFD", p.148, eq.(4.14.12)
Definition at line 788 of file navier_stokes.cpp.
◆ compute_viscous_stress_tensor_via_scaled_viscosity_and_strain_rate_tensor()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
| dealii::Tensor< 2, dim, real2 > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_viscous_stress_tensor_via_scaled_viscosity_and_strain_rate_tensor | ( | const real2 | scaled_viscosity_coefficient, |
| const dealii::Tensor< 2, dim, real2 > & | strain_rate_tensor | ||
| ) | const |
Nondimensionalized viscous stress tensor, tau* Reference: Masatsuka 2018 "I do like CFD", p.148, eq.(4.14.12)
Definition at line 755 of file navier_stokes.cpp.
◆ compute_vorticity_based_dissipation_rate_from_integrated_enstrophy()
template<int dim, int nspecies, int nstate, typename real>
| real PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_vorticity_based_dissipation_rate_from_integrated_enstrophy | ( | const real | integrated_enstrophy | ) | const |
Evaluate non-dimensional theoretical vorticity-based dissipation rate integrated enstrophy. Note: For incompressible flows or when dilatation effects are negligible – Reference: Cox, Christopher, et al. "Accuracy, stability, and performance comparison between the spectral difference and flux reconstruction schemes." Computers & Fluids 221 (2021): 104922. – Equation (56) with free-stream nondimensionalization applied
Definition at line 523 of file navier_stokes.cpp.
◆ compute_wall_shear_stress()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
| real2 PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_wall_shear_stress | ( | const std::array< real2, nstate > & | conservative_soln, |
| const std::array< dealii::Tensor< 1, dim, real2 >, nstate > & | conservative_soln_gradient, | ||
| const dealii::Tensor< 1, dim, real2 > & | normal_vector | ||
| ) | const |
Nondimensionalized wall shear stress
Definition at line 135 of file navier_stokes.cpp.
◆ compute_wall_tangent_vector()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
| dealii::Tensor< 1, dim, real2 > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_wall_tangent_vector | ( | const std::array< real2, nstate > & | conservative_soln, |
| const dealii::Tensor< 1, dim, real2 > & | normal_vector | ||
| ) | const |
Nondimensionalized wall tangent vector
Definition at line 101 of file navier_stokes.cpp.
◆ compute_wall_tangent_vector_from_velocities_parallel_to_wall()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
| dealii::Tensor< 1, dim, real2 > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::compute_wall_tangent_vector_from_velocities_parallel_to_wall | ( | const dealii::Tensor< 1, dim, real2 > & | velocities_parallel_to_wall | ) | const |
Nondimensionalized wall tangent vector from velocities parallel to wall
Definition at line 115 of file navier_stokes.cpp.
◆ convective_flux_directional_jacobian_via_dfad()
template<int dim, int nspecies, int nstate, typename real>
| dealii::Tensor< 2, nstate, real > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::convective_flux_directional_jacobian_via_dfad | ( | std::array< real, nstate > & | conservative_soln, |
| const dealii::Tensor< 1, dim, real > & | normal | ||
| ) | const |
Convective flux Jacobian computed via dfad (automatic differentiation) – Only used for verifying the dfad procedure used in dissipative flux jacobian
Definition at line 1120 of file navier_stokes.cpp.
◆ dissipative_flux()
template<int dim, int nspecies, int nstate, typename real>
|
overridevirtual |
Nondimensionalized viscous flux (i.e. dissipative flux) Reference: Masatsuka 2018 "I do like CFD", p.142, eq.(4.12.1-4.12.4)
Reimplemented from PHiLiP::Physics::Euler< dim, nspecies, nstate, real >.
Definition at line 898 of file navier_stokes.cpp.
◆ dissipative_flux_directional_jacobian()
template<int dim, int nspecies, int nstate, typename real>
| dealii::Tensor< 2, nstate, real > PHiLiP::Physics::NavierStokes< dim, nspecies, 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 |
Dissipative flux Jacobian Note: Only used for computing the manufactured solution source term; computed using automatic differentiation
Definition at line 955 of file navier_stokes.cpp.
◆ dissipative_flux_directional_jacobian_wrt_gradient_component()
template<int dim, int nspecies, int nstate, typename real>
| dealii::Tensor< 2, nstate, real > PHiLiP::Physics::NavierStokes< dim, nspecies, 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 |
Dissipative flux Jacobian wrt gradient component Note: Only used for computing the manufactured solution source term; computed using automatic differentiation
Definition at line 993 of file navier_stokes.cpp.
◆ dissipative_flux_dot_normal()
template<int dim, int nspecies, int nstate, typename real>
|
overridevirtual |
Nondimensionalized viscous flux (i.e. dissipative flux) dot normal vector that accounts for gradient boundary conditions References: (1) Masatsuka 2018 "I do like CFD", p.142, eq.(4.12.1-4.12.4), (2) For the boundary condition case, refer to the equation above equation 458 of the following paper: Hartmann, Ralf. "Numerical analysis of higher order discontinuous Galerkin finite element methods." (2008): 1-107.
If adiabatic on no-slip (1001) wall BCs
Reimplemented from PHiLiP::Physics::PhysicsBase< dim, nspecies, nstate, real >.
Definition at line 1224 of file navier_stokes.cpp.
◆ dissipative_flux_dot_normal_on_adiabatic_boundary()
template<int dim, int nspecies, int nstate, typename real>
|
virtual |
Nondimensionalized viscous flux (i.e. dissipative flux) dot normal vector that accounts for gradient boundary conditions when the on boundary flag is true References: (1) Masatsuka 2018 "I do like CFD", p.142, eq.(4.12.1-4.12.4), (2) For the boundary condition case, refer to the equation above equation 458 of the following paper: Hartmann, Ralf. "Numerical analysis of higher order discontinuous Galerkin finite element methods." (2008): 1-107.
If adiabatic on no-slip (1001) wall BCs
Reimplemented in PHiLiP::Physics::NavierStokes_ChannelFlowConstantSourceTerm_WallModel< dim, nspecies, nstate, real >.
Definition at line 1268 of file navier_stokes.cpp.
◆ dissipative_flux_given_velocities_viscous_stress_tensor_and_heat_flux()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
| std::array< dealii::Tensor< 1, dim, real2 >, nstate > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::dissipative_flux_given_velocities_viscous_stress_tensor_and_heat_flux | ( | const dealii::Tensor< 1, dim, real2 > & | vel, |
| const dealii::Tensor< 2, dim, real2 > & | viscous_stress_tensor, | ||
| const dealii::Tensor< 1, dim, real2 > & | heat_flux | ||
| ) | const |
Nondimensionalized viscous flux (i.e. dissipative flux) computed via given velocities, viscous stress tensor, and heat flux. Reference: Masatsuka 2018 "I do like CFD", p.142, eq.(4.12.1-4.12.4)
Definition at line 1322 of file navier_stokes.cpp.
◆ dissipative_flux_templated()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
|
protected |
Nondimensionalized viscous flux (i.e. dissipative flux) Reference: Masatsuka 2018 "I do like CFD", p.142, eq.(4.12.1-4.12.4)
Definition at line 1198 of file navier_stokes.cpp.
◆ extract_velocities_gradient_from_primitive_solution_gradient()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
| dealii::Tensor< 2, dim, real2 > PHiLiP::Physics::NavierStokes< dim, nspecies, nstate, real >::extract_velocities_gradient_from_primitive_solution_gradient | ( | const std::array< dealii::Tensor< 1, dim, real2 >, nstate > & | primitive_soln_gradient | ) | const |
Extract gradient of velocities
Definition at line 692 of file navier_stokes.cpp.
◆ scale_viscosity_coefficient()
template<int dim, int nspecies, int nstate, typename real >
template<typename real2 >
|
inline |
Scaled nondimensionalized viscosity coefficient, hat{mu*}, given nondimensionalized viscosity coefficient Reference: Masatsuka 2018 "I do like CFD", p.148, eq.(4.14.14)
Definition at line 310 of file navier_stokes.cpp.
Member Data Documentation
◆ sutherlands_temperature
template<int dim, int nspecies, int nstate, typename real>
|
protected |
Sutherland's temperature. Units: [K].
Constants for Sutherland's law for viscosity Reference: Sutherland, W. (1893), "The viscosity of gases and molecular force", Philosophical Magazine, S. 5, 36, pp. 507-531 (1893) Values: https://www.cfd-online.com/Wiki/Sutherland%27s_law
Definition at line 68 of file navier_stokes.h.
The documentation for this class was generated from the following files:
- src/physics/navier_stokes.h
- src/physics/navier_stokes.cpp
