CSourceVorticityConfinement Class Reference

Class for a source term due to vorticity confinement. More...

#include <flow_sources.hpp>

Inheritance diagram for CSourceVorticityConfinement:

Public Member Functions
	CSourceVorticityConfinement (unsigned short val_nDim, unsigned short val_nVar, const CConfig *config)
	Constructor of the class.
ResidualType	ComputeResidual (const CConfig *config) override
	Residual of the rotational frame source term.
Public Member Functions inherited from CSourceBase_Flow
	~CSourceBase_Flow () override
	Destructor of the class.
void	SetStreamwisePeriodicValues (const StreamwisePeriodicValues SolverSPvals) final
	Set massflow, heatflow & inlet temperature for streamwise periodic flow.
Public Member Functions inherited from CNumerics
	CNumerics (void)
	Constructor of the class.
	CNumerics (unsigned short val_nDim, unsigned short val_nVar, const CConfig *config)
virtual	~CNumerics (void)
	Destructor of the class.
void	SetTimeStep (su2double val_timestep)
	Set the time step.
void	SetVelocity2_Inf (su2double val_velocity2)
	Set the freestream velocity square.
void	SetVorticity (const su2double *val_vorticity_i, const su2double *val_vorticity_j)
	Set the value of the vorticity.
void	SetStrainMag (su2double val_strainmag_i, su2double val_strainmag_j)
	Set the value of the rate of strain magnitude.
void	SetConservative (const su2double *val_u_i, const su2double *val_u_j)
	Set the value of the conservative variables.
void	SetPrimitive (const su2double *val_v_i, const su2double *val_v_j)
	Set the value of the primitive variables.
void	SetSecondary (const su2double *val_s_i, const su2double *val_s_j)
	Set the value of the primitive variables.
void	SetConsVarGradient (CMatrixView< const su2double > val_consvar_grad_i, CMatrixView< const su2double > val_consvar_grad_j)
	Set the gradient of the conservative variables.
void	SetConsVarGradient (CMatrixView< const su2double > val_consvar_grad)
	Set the gradient of the conservative variables.
void	SetPrimVarGradient (CMatrixView< const su2double > val_primvar_grad_i, CMatrixView< const su2double > val_primvar_grad_j)
	Set the gradient of the primitive variables.
void	SetAdjointVar (const su2double *val_psi_i, const su2double *val_psi_j)
	Set the value of the adjoint variable.
void	SetAdjointVarGradient (CMatrixView< const su2double > val_psivar_grad_i, CMatrixView< const su2double > val_psivar_grad_j)
	Set the gradient of the adjoint variables.
void	SetScalarVar (const su2double *val_scalarvar_i, const su2double *val_scalarvar_j)
	Set the value of the scalar variable.
void	SetTransVar (const su2double *val_transvar_i, const su2double *val_transvar_j)
	Set the value of the turbulent variable.
void	SetScalarVarGradient (CMatrixView< const su2double > val_scalarvar_grad_i, CMatrixView< const su2double > val_scalarvar_grad_j)
	Set the gradient of the scalar variables.
void	SetTransVarGradient (CMatrixView< const su2double > val_transvar_grad_i, CMatrixView< const su2double > val_transvar_grad_j)
	Set the gradient of the turbulent variables.
void	SetLocalGridLength (const su2double val_localGridLength_i)
	Set the value of the turbulent variable.
void	SetTurbAdjointVar (const su2double *val_turbpsivar_i, const su2double *val_turbpsivar_j)
	Set the value of the adjoint turbulent variable.
void	SetTurbAdjointGradient (CMatrixView< const su2double > val_turbpsivar_grad_i, CMatrixView< const su2double > val_turbpsivar_grad_j)
	Set the gradient of the adjoint turbulent variables.
template<class Vec , class Mat >
void	ComputeStochReynStress (su2double density, su2double tke, const Vec &rndVec, su2double Cmag, Mat &stochReynStress)
	Compute a random contribution to the Reynolds stress tensor (Stochastic Backscatter Model).
su2double	ComputeStochRelaxFactor (const CConfig *config)
	Compute relaxation factor for stochastic source term in momentum equations (Stochastic Backscatter Model).
virtual void	SetF1blending (su2double val_F1_i, su2double val_F1_j)
	Set the value of the first blending function.
virtual void	SetF2blending (su2double val_F2_i)
	Set the value of the second blending function.
virtual void	SetCrossDiff (su2double val_CDkw_i)
	Set the value of the cross diffusion for the SST model.
void	SetIntermittencyEff (su2double val_intermittency_eff_i)
	Set the value of the effective intermittency for the LM model.
void	SetIntermittency (su2double val_intermittency_i)
	Set the value of the intermittency for the LM model.
su2double	GetIntermittencyEff () const
	Get the value of the effective intermittency for the transition model.
void	SetAuxVarGrad (CMatrixView< const su2double > val_auxvar_grad_i, CMatrixView< const su2double > val_auxvar_grad_j)
	Set the gradient of the auxiliary variables.
void	SetDiffusionCoeff (const su2double *val_diffusioncoeff_i, const su2double *val_diffusioncoeff_j)
	Set the diffusion coefficient.
void	SetHeatFluxDiffusion (su2double val_heatfluxdiffusion)
	Set the heat flux due to enthalpy diffusion.
void	SetJacHeatFluxDiffusion (su2double val_jacheatfluxdiffusion)
	Set Jacobian of the heat flux due to enthalpy diffusion.
void	SetLaminarViscosity (su2double val_laminar_viscosity_i, su2double val_laminar_viscosity_j)
	Set the laminar viscosity.
void	SetThermalConductivity (su2double val_thermal_conductivity_i, su2double val_thermal_conductivity_j)
	Set the thermal conductivity (translational/rotational)
void	SetThermalConductivity_ve (su2double val_thermal_conductivity_ve_i, su2double val_thermal_conductivity_ve_j)
	Set the thermal conductivity (translational/rotational)
void	SetSpecificHeat (su2double val_specific_heat_i, su2double val_specific_heat_j)
	Set the specifc heat c_p.
void	SetEddyViscosity (su2double val_eddy_viscosity_i, su2double val_eddy_viscosity_j)
	Set the eddy viscosity.
void	SetTurbKineticEnergy (su2double val_turb_ke_i, su2double val_turb_ke_j)
	Set the turbulent kinetic energy.
void	SetStochVar (unsigned short iDim, su2double val_stochvar_i, su2double val_stochvar_j)
	Set the stochastic variables from Langevin equations (Stochastic Backscatter Model).
void	SetSbsInBoxSensor (int8_t val_sbsInBox_i)
	Set the sensor to locate the box where the Stochastic Backscatter Model is active.
void	SetLES_Mode (su2double val_lesMode_i, su2double val_lesMode_j)
	Set the LES sensor for hybrid RANS-LES methods.
void	SetDistance (su2double val_dist_i, su2double val_dist_j)
	Set the value of the distance from the nearest wall.
void	SetStochSource (su2double val_stoch_source, unsigned short iDim)
	Set the stochastic source term for the Langevin equations (Backscatter Model).
void	SetRoughness (su2double val_roughness_i, su2double val_roughness_j)
	Set the value of the roughness from the nearest wall.
void	SetCoord (const su2double *val_coord_i, const su2double *val_coord_j)
	Set coordinates of the points.
void	SetGridVel (const su2double *val_gridvel_i, const su2double *val_gridvel_j)
	Set the velocity of the computational grid.
void	SetPressure (su2double val_pressure_i, su2double val_pressure_j)
	Set the value of the pressure.
void	SetDensity (su2double val_densityinc_i, su2double val_densityinc_j)
	Set the value of the density for the incompressible solver.
void	SetBetaInc2 (su2double val_betainc2_i, su2double val_betainc2_j)
	Set the value of the beta for incompressible flows.
void	SetSoundSpeed (su2double val_soundspeed_i, su2double val_soundspeed_j)
	Set the value of the sound speed.
void	SetEnthalpy (su2double val_enthalpy_i, su2double val_enthalpy_j)
	Set the value of the enthalpy.
void	SetLambda (su2double val_lambda_i, su2double val_lambda_j)
	Set the value of the spectral radius.
void	SetUndivided_Laplacian (const su2double *val_und_lapl_i, const su2double *val_und_lapl_j)
	Set the value of undivided laplacian.
void	SetSensor (su2double val_sensor_i, su2double val_sensor_j)
	Set the value of the pressure sensor.
void	SetNeighbor (unsigned short val_neighbor_i, unsigned short val_neighbor_j)
	Set the number of neighbor to a point.
void	SetNormal (const su2double *val_normal)
	Set the value of the normal vector to the face between two points.
void	SetVolume (su2double val_volume)
	Set the value of the volume of the control volume.
void	SetAvgVolume (su2double val_avgvolume)
	Set the value of AvgVolume Variable.
void	SetDissipation (su2double diss_i, su2double diss_j)
	Sets the values of the roe dissipation.
su2double	GetDissipation () const
	Get the final Roe dissipation factor.
void	GetInviscidProjFlux (const su2double *val_density, const su2double *val_velocity, const su2double *val_pressure, const su2double *val_enthalpy, const su2double *val_normal, su2double *val_Proj_Flux) const
	Compute the projected inviscid flux vector.
void	GetInviscidIncProjFlux (const su2double *val_density, const su2double *val_velocity, const su2double *val_pressure, const su2double *val_betainc2, const su2double *val_enthalpy, const su2double *val_normal, su2double *val_Proj_Flux) const
	Compute the projected inviscid flux vector for incompresible simulations.
void	GetInviscidProjJac (const su2double *val_velocity, const su2double *val_energy, const su2double *val_normal, su2double val_scale, su2double **val_Proj_Jac_tensor) const
	Compute the projection of the inviscid Jacobian matrices.
void	GetInviscidIncProjJac (const su2double *val_density, const su2double *val_velocity, const su2double *val_betainc2, const su2double *val_normal, su2double val_scale, su2double **val_Proj_Jac_tensor)
	Compute the projection of the inviscid Jacobian matrices (incompressible).
void	GetInviscidIncProjJac (const su2double *val_density, const su2double *val_velocity, const su2double *val_betainc2, const su2double *val_enthalpy, const su2double *val_dRhodh, const su2double *val_normal, su2double val_scale, su2double **val_Proj_Jac_Tensor) const
	Compute the projection of the inviscid Jacobian matrices (overload for low speed preconditioner version).
void	GetPreconditioner (const su2double *val_density, const su2double *val_velocity, const su2double *val_betainc2, const su2double *val_enthalpy, const su2double *val_drhodh, su2double **val_Precon) const
	Compute the low speed preconditioning matrix.
void	GetPreconditionedProjJac (const su2double *val_density, const su2double *val_velocity, const su2double *val_betainc2, const su2double *val_normal, su2double **val_Proj_Jac_Tensor) const
	Compute the projection of the preconditioned inviscid Jacobian matrices.
void	GetInviscidProjJac (const su2double *val_velocity, const su2double *val_enthalphy, const su2double *val_chi, const su2double *val_kappa, const su2double *val_normal, su2double val_scale, su2double **val_Proj_Jac_tensor) const
	Compute the projection of the inviscid Jacobian matrices for general fluid model.
void	GetPrimitive2Conservative (const su2double *val_Mean_PrimVar, const su2double *val_Mean_SecVar, su2double **val_Jac_PC) const
	Mapping between primitives variables P and conservatives variables C.
void	GetPMatrix (const su2double *val_density, const su2double *val_velocity, const su2double *val_soundspeed, const su2double *val_enthalpy, const su2double *val_chi, const su2double *val_kappa, const su2double *val_normal, su2double **val_p_tensor) const
	Computation of the matrix P for a generic fluid model.
template<typename Matrix >
void	GetPMatrix (const su2double &density, const su2double *velocity, const su2double &soundspeed, const su2double *normal, Matrix &p_tensor) const
	Computation of the matrix P, this matrix diagonalize the conservative Jacobians in the form $P^{-1}(A.Normal)P=Lambda$.
void	GetinvRinvPe (su2double Beta2, su2double val_enthalpy, su2double val_soundspeed, su2double val_density, const su2double *val_velocity, su2double **val_invR_invPe) const
	Computation of the matrix Rinv*Pe.
void	GetRMatrix (su2double val_pressure, su2double val_soundspeed, su2double val_density, const su2double *val_velocity, su2double **val_invR_invPe) const
	Computation of the matrix R.
void	GetRMatrix (su2double val_soundspeed, su2double val_density, su2double **R_Matrix) const
	Computation of the matrix R.
void	GetLMatrix (su2double val_soundspeed, su2double val_density, su2double **L_Matrix) const
	Computation of the matrix R.
void	ComputeResJacobianGiles (CFluidModel *FluidModel, su2double pressure, su2double density, const su2double *turboVel, su2double alphaInBC, su2double gammaInBC, su2double **R_c, su2double **R_c_inv)
	Computation of the flow Residual Jacobian Matrix for Non Reflecting BC.
void	InvMatrix3D (su2double **matrix, su2double **invMatrix)
	Computate the inverse of a 3x3 matrix.
void	InvMatrix4D (su2double **matrix, su2double **invMatrix)
	Computate the inverse of a 4x4 matrix.
void	GetCharJump (su2double val_soundspeed, su2double val_density, const su2double *prim_jump, su2double *char_jump) const
	Computation of the matrix R.
void	GetPrecondJacobian (su2double Beta2, su2double r_hat, su2double s_hat, su2double t_hat, su2double rB2a2, const su2double *val_Lambda, const su2double *val_normal, su2double **val_absPeJac) const
	Computation of the matrix Td, this matrix diagonalize the preconditioned conservative Jacobians in the form $Tg |Lambda| Td = Pc{-1}|Pc (A.Normal)|$.
void	GetPMatrix_inv (su2double **val_invp_tensor, const su2double *val_density, const su2double *val_velocity, const su2double *val_soundspeed, const su2double *val_chi, const su2double *val_kappa, const su2double *val_normal) const
	Computation of the matrix P^{-1}, this matrix diagonalize the conservative Jacobians in the form $P^{-1}(A.Normal)P=Lambda$.
template<typename Matrix >
void	GetPMatrix_inv (const su2double &density, const su2double *velocity, const su2double &soundspeed, const su2double *normal, Matrix &inv_p_tensor) const
	Computation of the matrix P^{-1}, this matrix diagonalize the conservative Jacobians in the form $P^{-1}(A.Normal)P=Lambda$.
void	GetAdjViscousFlux_Jac (su2double Pressure_i, su2double Pressure_j, su2double Density_i, su2double Density_j, su2double ViscDens_i, su2double ViscDens_j, const su2double *Velocity_i, const su2double *Velocity_j, su2double sq_vel_i, su2double sq_vel_j, su2double XiDens_i, su2double XiDens_j, su2double **Mean_GradPhi, const su2double *Mean_GradPsiE, su2double dPhiE_dn, const su2double *Normal, const su2double *Edge_Vector, su2double dist_ij_2, su2double *val_residual_i, su2double *val_residual_j, su2double **val_Jacobian_ii, su2double **val_Jacobian_ij, su2double **val_Jacobian_ji, su2double **val_Jacobian_jj, bool implicit) const
	Compute viscous residual and jacobian.
void	GetJacInviscidLambda_fabs (const su2double *val_velocity, su2double val_soundspeed, const su2double *val_normal, su2double *val_Lambda_Vector) const
	Computation of the projected inviscid lambda (eingenvalues).
virtual void	ComputeResidual (su2double *val_residual, CConfig *config)
	Compute the numerical residual.
virtual void	ComputeResidual (su2double *val_residual, su2double **val_Jacobian_i, su2double **val_Jacobian_j, CConfig *config)
virtual void	ComputeResidual (su2double *val_residual_i, su2double *val_residual_j, su2double **val_Jacobian_ii, su2double **val_Jacobian_ij, su2double **val_Jacobian_ji, su2double **val_Jacobian_jj, CConfig *config)
virtual void	ComputeResidual (su2double *val_resconv_i, su2double *val_resvisc_i, su2double *val_resconv_j, su2double *val_resvisc_j, su2double **val_Jacobian_ii, su2double **val_Jacobian_ij, su2double **val_Jacobian_ji, su2double **val_Jacobian_jj, CConfig *config)
virtual void	ComputeResidual (su2double *val_resconv, su2double *val_resvisc, su2double **val_Jacobian_i, su2double **val_Jacobian_j, CConfig *config)
virtual void	ComputeResidual (su2double *val_residual, su2double **val_Jacobian_i, CConfig *config)
virtual void	ComputeResidual_TransLM (su2double *val_residual, su2double **val_Jacobian_i, su2double **val_Jacobian_j, CConfig *config, su2double &gamma_sep)
	Residual for transition problems.
virtual ResidualType	ComputeAxisymmetric (const CConfig *config)
	Residual for source term integration.
virtual ResidualType	ComputeVibRelaxation (const CConfig *config)
virtual ResidualType	ComputeChemistry (const CConfig *config)
	Calculation of the chemistry source term.
virtual su2double	GetGammaBC (void) const
	A virtual member.
virtual void	Compute_Tangent_Matrix (CElement *element_container, const CConfig *config)
	A virtual member to compute the tangent matrix in structural problems.
virtual void	Compute_NodalStress_Term (CElement *element_container, const CConfig *config)
	A virtual member to compute the nodal stress term in non-linear structural problems.
virtual void	SetMeshElasticProperties (unsigned long iElem, su2double val_E)
	Set the element-based local Young's modulus in mesh problems.
virtual void	Set_DV_Val (unsigned short i_DV, su2double val_DV)
	A virtual member to set the value of the design variables.
virtual su2double	Get_DV_Val (unsigned short i_DV) const
	A virtual member to retrieve the value of the design variables.
virtual void	Set_ElectricField (unsigned short i_DV, su2double val_EField)
	A virtual member to set the electric field.
virtual void	SetMaterial_Properties (unsigned short iVal, su2double val_E, su2double val_Nu)
	A virtual member to set the material properties.
virtual void	SetMaterial_Density (unsigned short iVal, su2double val_Rho, su2double val_Rho_DL)
	A virtual member to set the material properties.
virtual void	Compute_Mass_Matrix (CElement *element_container, const CConfig *config)
	A virtual member to compute the mass matrix.
virtual void	Compute_Body_Forces (CElement *element_container, const CConfig *config)
	A virtual member to compute the residual component due to body forces.
virtual su2double	Compute_Averaged_NodalStress (CElement *element_container, const CConfig *config)
	A virtual member to compute the averaged nodal stresses.
void	CreateBasis (const su2double *val_Normal, su2double *l, su2double *m)
	Computes a basis of orthogonal vectors from a supplied vector.
virtual void	SetTau_Wall (su2double val_tauwall_i, su2double val_tauwall_j)
	Set the value of the Tauwall.
su2double	GetRoe_Dissipation (const su2double Dissipation_i, const su2double Dissipation_j, const su2double Sensor_i, const su2double Sensor_j, const CConfig *config) const
	Calculate the central/upwind blending function for a face
virtual void	SetRadVar (const su2double *val_radvar_i, const su2double *val_radvar_j)
	Set the value of the radiation variable.
virtual void	SetRadVarGradient (CMatrixView< const su2double > val_radvar_grad_i, CMatrixView< const su2double > val_radvar_grad_j)
	Set the gradient of the radiation variables.
void	SetRadVarSource (const su2double *val_radvar_source)
	Set the gradient of the radiation variables.
virtual void	SetdPdU (su2double *val_dPdU_i, su2double *val_dPdU_j)
	Set the pressure derivatives.
virtual void	SetdTdU (su2double *val_dTdU_i, su2double *val_dTdU_j)
	Set the temperature derivatives.
virtual void	SetdTvedU (su2double *val_dTvedU_i, su2double *val_dTvedU_j)
	Set the vib-el temperture derivatives.
virtual void	SetEve (su2double *val_Eve_i, su2double *val_Eve_j)
	Set the vib-elec energy.
virtual void	SetCvve (su2double *val_Cvve_i, su2double *val_Cvve_j)
	Set the vib-elec specific heat.
virtual void	SetGamma (su2double val_Gamma_i, su2double val_Gamma_j)
	Set the ratio of specific heats.
void	SetMassFlux (const su2double val_MassFlux)
	SetMassFlux.
bool	GetBoundedScalar () const
	Obtain information on bounded scalar problem.

Additional Inherited Members
Static Public Member Functions inherited from CNumerics
template<class Mat1 , class Mat2 >
static void	ComputeMeanRateOfStrainMatrix (size_t nDim, Mat1 &rateofstrain, const Mat2 &velgrad)
	Compute the mean rate of strain matrix.
template<class Mat1 , class Mat2 , class Scalar >
static void	ComputeStressTensor (size_t nDim, Mat1 &stress, const Mat2 &velgrad, Scalar viscosity, Scalar density=0.0, Scalar turb_ke=0.0, bool reynolds3x3=false)
	Compute the stress tensor from the velocity gradients.
template<class Mat1 , class Mat2 , class Scalar2 >
static void	AddQCR (size_t nDim, const Mat1 &gradvel, Mat2 &tau, Scalar2 turb_fraction)
	Add a correction using a Quadratic Constitutive Relation to the stress tensor.
template<class Mat1 , class Mat2 , class Scalar >
static void	ComputePerturbedRSM (size_t nDim, size_t uq_eigval_comp, bool uq_permute, su2double uq_delta_b, su2double uq_urlx, const Mat1 &velgrad, Scalar density, Scalar viscosity, Scalar turb_ke, Mat2 &MeanPerturbedRSM)
	Perturb the Reynolds stress tensor based on parameters.
template<class Vec1 , class Vec2 , class Mat >
static su2double	ComputeProjectedGradient (int nDim, int nVar, const Vec1 &normal, const Vec1 &coord_i, const Vec1 &coord_j, const Mat &grad_i, const Mat &grad_j, bool correct, const Vec2 &var_i, const Vec2 &var_j, su2double *projNormal, su2double *projCorrected)
	Project average gradient onto normal (with or w/o correction) for viscous fluxes of scalar quantities.
static bool	CheckResidualNaNs (bool implicit, int nVar, const ResidualType<> residual)
	Check if residual constains a NaN value.
Protected Member Functions inherited from CSourceBase_Flow
	CSourceBase_Flow (unsigned short val_nDim, unsigned short val_nVar, const CConfig *config)
	Constructor of the class.
Protected Attributes inherited from CSourceBase_Flow
su2double *	residual = nullptr
su2double **	jacobian = nullptr
struct StreamwisePeriodicValues	SPvals
Protected Attributes inherited from CNumerics
unsigned short	nDim
unsigned short	nVar
	Number of dimensions and variables.
su2double	Gamma
	Fluid's Gamma constant (ratio of specific heats).
su2double	Gamma_Minus_One
	Fluids's Gamma - 1.0 .
su2double	Minf
	Free stream Mach number .
su2double	Gas_Constant
	Gas constant.
su2double	Prandtl_Turb
	Turbulent Prandtl's number.
su2double	MassFlux
	Mass flux across edge.
su2double *	Proj_Flux_Tensor
	Flux tensor projected in a direction.
su2double **	tau
	Viscous stress tensor.
const su2double	delta [3][3] = {{1.0, 0.0, 0.0},{0.0, 1.0, 0.0}, {0.0, 0.0, 1.0}}
	Identity matrix.
const su2double *	Diffusion_Coeff_i
	Species diffusion coefficients at point i.
const su2double *	Diffusion_Coeff_j
	Species diffusion coefficients at point j.
su2double	Laminar_Viscosity_i
	Laminar viscosity at point i.
su2double	Laminar_Viscosity_j
	Laminar viscosity at point j.
su2double	Thermal_Conductivity_i
	Thermal conductivity at point i.
su2double	Thermal_Conductivity_j
	Thermal conductivity at point j.
su2double	Thermal_Conductivity_ve_i
	vibrational-electronic Thermal conductivity at point i.
su2double	Thermal_Conductivity_ve_j
	vibrational-electronic Thermal conductivity at point j.
su2double	Cp_i
	Cp at point i.
su2double	Cp_j
	Cp at point j.
su2double	Eddy_Viscosity_i
	Eddy viscosity at point i.
su2double	Eddy_Viscosity_j
	Eddy viscosity at point j.
su2double	turb_ke_i
	Turbulent kinetic energy at point i.
su2double	turb_ke_j
	Turbulent kinetic energy at point j.
su2double	intermittency_eff_i
	effective intermittency at point i.
su2double	intermittency_i
	intermittency at point i.
su2double	Pressure_i
	Pressure at point i.
su2double	Pressure_j
	Pressure at point j.
su2double	Density_i
	Density at point i.
su2double	Density_j
	Density at point j.
su2double	DensityInc_i
	Incompressible density at point i.
su2double	DensityInc_j
	Incompressible density at point j.
su2double	BetaInc2_i
	Beta incompressible at point i.
su2double	BetaInc2_j
	Beta incompressible at point j.
su2double	Lambda_i
	Spectral radius at point i.
su2double	Lambda_j
	Spectral radius at point j.
su2double	SoundSpeed_i
	Sound speed at point i.
su2double	SoundSpeed_j
	Sound speed at point j.
su2double	Enthalpy_i
	Enthalpy at point i.
su2double	Enthalpy_j
	Enthalpy at point j.
su2double	dist_i
	Distance of point i to the nearest wall.
su2double	dist_j
	Distance of point j to the nearest wall.
const su2double *	Und_Lapl_i
	Undivided laplacians at point i.
const su2double *	Und_Lapl_j
	Undivided laplacians at point j.
su2double	Sensor_i
	Pressure sensor at point i.
su2double	Sensor_j
	Pressure sensor at point j.
const su2double *	GridVel_i
	Grid velocity at point i.
const su2double *	GridVel_j
	Grid velocity at point j.
const su2double *	U_i
	Vector of conservative variables at point i.
const su2double *	U_j
	Vector of conservative variables at point j.
const su2double *	V_i
	Vector of primitive variables at point i.
const su2double *	V_j
	Vector of primitive variables at point j.
const su2double *	S_i
	Vector of secondary variables at point i.
const su2double *	S_j
	Vector of secondary variables at point j.
const su2double *	Psi_i
	Vector of adjoint variables at point i.
const su2double *	Psi_j
	Vector of adjoint variables at point j.
const su2double *	ScalarVar_i
	Vector of scalar variables at point i.
const su2double *	ScalarVar_j
	Vector of scalar variables at point j.
su2double	HeatFluxDiffusion
	Heat flux due to enthalpy diffusion for multicomponent.
su2double	JacHeatFluxDiffusion
	Heat flux jacobian due to enthalpy diffusion for multicomponent.
const su2double *	TransVar_i
	Vector of turbulent variables at point i.
const su2double *	TransVar_j
	Vector of turbulent variables at point j.
const su2double *	TurbPsi_i
	Vector of adjoint turbulent variables at point i.
const su2double *	TurbPsi_j
	Vector of adjoint turbulent variables at point j.
CMatrixView< const su2double >	ConsVar_Grad_i
	Gradient of conservative variables at point i.
CMatrixView< const su2double >	ConsVar_Grad_j
	Gradient of conservative variables at point j.
CMatrixView< const su2double >	ConsVar_Grad
	Gradient of conservative variables which is a scalar.
CMatrixView< const su2double >	PrimVar_Grad_i
	Gradient of primitive variables at point i.
CMatrixView< const su2double >	PrimVar_Grad_j
	Gradient of primitive variables at point j.
CMatrixView< const su2double >	PsiVar_Grad_i
	Gradient of adjoint variables at point i.
CMatrixView< const su2double >	PsiVar_Grad_j
	Gradient of adjoint variables at point j.
CMatrixView< const su2double >	ScalarVar_Grad_i
	Gradient of scalar variables at point i.
CMatrixView< const su2double >	ScalarVar_Grad_j
	Gradient of scalar variables at point j.
CMatrixView< const su2double >	TransVar_Grad_i
	Gradient of turbulent variables at point i.
CMatrixView< const su2double >	TransVar_Grad_j
	Gradient of turbulent variables at point j.
CMatrixView< const su2double >	TurbPsi_Grad_i
	Gradient of adjoint turbulent variables at point i.
CMatrixView< const su2double >	TurbPsi_Grad_j
	Gradient of adjoint turbulent variables at point j.
CMatrixView< const su2double >	AuxVar_Grad_i
	Gradient of an auxiliary variable at point i.
CMatrixView< const su2double >	AuxVar_Grad_j
	Gradient of an auxiliary variable at point i.
su2double	LocalGridLength_i
	Local grid length at point i.
const su2double *	RadVar_Source
	Source term from the radiative heat transfer equation.
const su2double *	Coord_i
	Cartesians coordinates of point i.
const su2double *	Coord_j
	Cartesians coordinates of point j.
unsigned short	Neighbor_i
	Number of neighbors of the point i.
unsigned short	Neighbor_j
	Number of neighbors of the point j.
const su2double *	Normal = nullptr
	Normal vector, its norm is the area of the face.
su2double	UnitNormal [MAXNDIM] = {0.0}
	Unitary normal vector.
su2double	TimeStep
	Time step useful in dual time method.
su2double	Area
	Area of the face i-j.
su2double	Volume
	Volume of the control volume around point i.
su2double	AvgVolume
	Average of the control Volume around point i for vorticity confinement parameter correction.
su2double	vel2_inf
	value of the square of freestream speed.
const su2double *	WindGust_i
	Wind gust at point i.
const su2double *	WindGust_j
	Wind gust at point j.
const su2double *	Vorticity_i
const su2double *	Vorticity_j
	Vorticity.
su2double	StrainMag_i
su2double	StrainMag_j
	Strain rate magnitude.
su2double	Dissipation_i
su2double	Dissipation_j
	Dissipation.
su2double	Dissipation_ij
su2double	roughness_i = 0.0
	Roughness of the wall nearest to point i.
su2double	roughness_j = 0.0
	Roughness of the wall nearest to point j.
su2double	MeanPerturbedRSM [3][3]
	Perturbed Reynolds stress tensor
su2double	stochReynStress [3][3] = {{0.0}}
	Stochastic contribution to Reynolds stress tensor for Backscatter Model.
su2double	stochSource [3] = {0.0}
	Source term for Langevin equations in Stochastic Backscatter Model.
su2double	stochVar_i [3] = {0.0}
	Stochastic variables at point i for Stochastic Backscatter Model.
su2double	stochVar_j [3] = {0.0}
	Stochastic variables at point j for Stochastic Backscatter Model.
su2double	lesMode_i = 0.0
	LES sensor at point i for hybrid RANS-LES methods.
su2double	lesMode_j = 0.0
	LES sensor at point j for hybrid RANS-LES methods.
int8_t	sbsInBox_i = 0
	Sensor to assess if point i lies inside the box where the Stochastic Backscatter Model is active.
SST_ParsedOptions	sstParsedOptions
	additional options for the SST turbulence model
unsigned short	Eig_Val_Comp
	Component towards which perturbation is perfromed.
su2double	uq_delta_b
	Magnitude of perturbation.
su2double	uq_urlx
	Under-relaxation factor for numerical stability.
bool	uq_permute
	Flag for eigenvector permutation.
bool	nemo
	Flag for NEMO problems
bool	energy_multicomponent = false
	Flag for multicomponent and reacting flow
bool	bounded_scalar = false
	Flag for bounded scalar problem.
Static Protected Attributes inherited from CNumerics
static constexpr size_t	MAXNDIM = 3
	Max number of space dimensions, used in some static arrays.

Detailed Description

Class for a source term due to vorticity confinement.

Vorticity Confinement (VC) technique to counter the numerical diffusion offered by the numerical scheme

Author

: Y Chandukrishna, Josy P Pullockara, T N Venkatesh. Computational and Theoretical Fluid Dynamics division, CSIR-National Aerospace Laboratories (NAL), Bangalore. Academy of Scientific and Innovative Research (AcSIR), Ghaziabad. First release date : 23 December 2022 modified on:

VC technique introduces an additional term to the N-S equations that counters the numerical difusion offerd by the numerical schemes. The additional term is introduced as a source term. VC technique requires an input confinement parameter that controls the magnitude of the source term. A suitable Confinement parameter is problem, scheme and grid dependant.

Required changes in config file - VORTICITY_CONFINEMENT = YES CONFINEMENT_PARAMETER = <0.0(default)> Currently have support for Viscous, Inviscid cases for both 2D and 3D problems. Can be used along with other source terms also.

The current implementation follows, R. Loehner and C. Yang, Vorticity confinement on unstructured grids, 2002. N. Butsuntorn and A. Jameson, Time Spectral Method for Rotorcraft Flow, 2008.

Constructor & Destructor Documentation

CSourceVorticityConfinement()

CSourceVorticityConfinement::CSourceVorticityConfinement	(	unsigned short	val_nDim,
		unsigned short	val_nVar,
		const CConfig *	config
	)

Constructor of the class.

Parameters

[in]	val_nDim	- Number of dimensions of the problem.
[in]	val_nVar	- Number of variables of the problem.
[in]	config	- Definition of the particular problem.

Member Function Documentation

ComputeResidual()

CNumerics::ResidualType CSourceVorticityConfinement::ComputeResidual ( const CConfig *  config )

overridevirtual

Residual of the rotational frame source term.

Parameters

[in]

config

- Definition of the particular problem.

Returns

Lightweight const-view of residual and Jacobian.

Reimplemented from CNumerics.

---

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

• SU2_CFD/include/numerics/flow/flow_sources.hpp
• SU2_CFD/src/numerics/flow/flow_sources.cpp