sphereEigenmode.h File Reference

Go to the source code of this file.

Macros
#define	U1(x, y, z)   (amp*cost*psi*( a0xy*(y) -a0zx*(z) ))
#define	V1(x, y, z)   (amp*cost*psi*(-a0xy*(x)+a0yz*(z) ))
#define	W1(x, y, z)   (amp*cost*psi*( -a0yz*(y)+a0zx*(x) ))
#define	U2(x, y, z)   (amp*cost*psi*( a0xy*(y)*(z) -a0zx*(z)*(y) ))
#define	V2(x, y, z)   (amp*cost*psi*(-a0xy*(x)*(z)+a0yz*(z)*(x) ))
#define	W2(x, y, z)   (amp*cost*psi*( -a0yz*(y)*(x) +a0zx*(x)*(y) ))
#define	VB2A(xa, wxa, pxa)
#define	VB2(xa, wxa, pxa)
#define	VB2X(xj, wj, pj,xa, wa, pa, wja, pja, deltaja)

Functions
	if (pdeVariation==SmParameters::hemp)
	assert (mg.numberOfDimensions()==3)
	if (vibrationClass==1)

Variables
const int	v2c = parameters.dbase.get<int >("v2c")
const int	v3c = parameters.dbase.get<int >("v3c")
bool	assignVelocities = v1c>=0
const int	s11c = parameters.dbase.get<int >("s11c")
const int	s12c = parameters.dbase.get<int >("s12c")
const int	s13c = parameters.dbase.get<int >("s13c")
const int	s21c = parameters.dbase.get<int >("s21c")
const int	s22c = parameters.dbase.get<int >("s22c")
const int	s23c = parameters.dbase.get<int >("s23c")
const int	s31c = parameters.dbase.get<int >("s31c")
const int	s32c = parameters.dbase.get<int >("s32c")
const int	s33c = parameters.dbase.get<int >("s33c")
const int	pc = parameters.dbase.get<int >("pc")
const real &	rho =parameters.dbase.get<real>("rho")
const real &	mu = parameters.dbase.get<real>("mu")
const real &	lambda = parameters.dbase.get<real>("lambda")
const RealArray &	muGrid = parameters.dbase.get<RealArray>("muGrid")
const RealArray &	lambdaGrid = parameters.dbase.get<RealArray>("lambdaGrid")
bool	assignStress = s11c >=0
std::vector< real > &	data = parameters.dbase.get<std::vector<real> >("sphereEigenmodeData")
const int	vibrationClass = (int)data[0]
const int	n = (int)data[1]
const int	m = (int)data[2]
const real	rad = data[3]
const real	cp =sqrt( (lambda+2.*mu)/rho )
const real	cs = sqrt( mu/rho )
real	omega
real	cPhi =1.
real	cOmega =1.
	else
real	kappa = omega/cs
real	h = omega/cp
const real	h2 = h*h
const real	h3 =h*h*h
const real	h4 =h*h*h*h
const real	kappa2 =kappa*kappa
const real	kappa3 =kappa*kappa*kappa
const real	kappa4 =kappa*kappa*kappa*kappa

Macro Definition Documentation

#define U1	(		x,
			y,
			z
	)		(amp*cost*psi*( a0xy*(y) -a0zx*(z) ))

#define U2	(		x,
			y,
			z
	)		(amp*cost*psi*( a0xy*(y)*(z) -a0zx*(z)*(y) ))

#define V1	(		x,
			y,
			z
	)		(amp*cost*psi*(-a0xy*(x)+a0yz*(z) ))

#define V2	(		x,
			y,
			z
	)		(amp*cost*psi*(-a0xy*(x)*(z)+a0yz*(z)*(x) ))

#define VB2	(		xa,
			wxa,
			pxa
	)

Value:

( (-1./(h*h))*( psi2hr*(wxa) + h*h*(xa)*psi3hr*w )   \
            + psi1kr*(pxa) - (n/(n+1.))*psi3kr*kappa*kappa*( r*r*(pxa) - n2p1*(xa)*p ) )

#define VB2A	(		xa,
			wxa,
			pxa
	)

Value:

( (-1./(h*h))*( psi2hr + (hr2/n2p1)*psi3hr )*(wxa)   \
            + (1./n2p1)*psi3hr*( r*r*(wxa) - n2p1*(xa)*w )  \
            + psi1kr*(pxa) - (n/(n+1.))*psi3kr*kappa*kappa*( r*r*(pxa) - n2p1*(xa)*p ) )

#define VB2X	(		xj,
			wj,
			pj,
			xa,
			wa,
			pa,
			wja,
			pja,
			deltaja
	)

Value:

( (-1./(h2))*( h2*xa*psi3hr*wj + psi2hr*wja + h2*deltaja*psi3hr*w + h4*xj*xa*psi4hr*w + h2*xj*psi3hr*wa )\
            + kappa2*xa*psi2kr*pj + psi1kr*pja \
            - (n/(n+1.))*kappa4*xa*psi4kr*( r*r*pj - n2p1*xj*p )\
            - (n/(n+1.))*kappa2*psi3kr*( 2.*xa*pj - n2p1*(deltaja)*p + r*r*pja - n2p1*xj*pa ) )

#define W1	(		x,
			y,
			z
	)		(amp*cost*psi*( -a0yz*(y)+a0zx*(x) ))

#define W2	(		x,
			y,
			z
	)		(amp*cost*psi*( -a0yz*(y)*(x) +a0zx*(x)*(y) ))

Function Documentation

assert

(

mg.

numberOfDimensions() = =3

)

if

(

pdeVariation

= = SmParameters::hemp

)

References printF().

else if

(

vibrationClass

= =1

)

References assert(), OV_ABORT(), and rad.

Variable Documentation

bool assignStress = s11c >=0

bool assignVelocities = v1c>=0

real cOmega =1.

Referenced by Cgsm::assignSpecialInitialConditions().

const real cp =sqrt( (lambda+2.*mu)/rho )

real cPhi =1.

Referenced by Cgsm::assignSpecialInitialConditions().

const real cs = sqrt( mu/rho )

std::vector<real>& data = parameters.dbase.get<std::vector<real> >("sphereEigenmodeData")

else

Initial value:

{
    OV_ABORT("ERROR: vibrationClass")

h = omega/cp

Referenced by Cgcns::addForcing(), Cgsm::assignSpecialInitialConditions(), MaterialProperties::hF(), MaterialProperties::hMinusRef(), MaterialProperties::lnKp(), main(), Maxwell::setupGrids(), TimeFunction::update(), and DomainSolver::userDefinedInitialConditions().

const real h2 = h*h

Referenced by Cgsm::assignSpecialInitialConditions().

const real h3 =h*h*h

Referenced by Cgsm::assignSpecialInitialConditions().

const real h4 =h*h*h*h

Referenced by Cgsm::assignSpecialInitialConditions().

kappa = omega/cs

Referenced by Cgad::addForcing(), Cgins::addForcing(), Cgins::addForcingToPressureEquation(), Cgad::advanceADI(), Cgmp::assignInterfaceBoundaryConditions(), Cgsm::assignSpecialInitialConditions(), AdParameters::displayPdeParameters(), Cgad::formMatrixForImplicitSolve(), PenaltyWallFunctionBC::getShearStresses(), getSpalartAllmarasParameters(), Cgad::getTimeSteppingEigenvalue(), Cgad::getUt(), Maxwell::initializePlaneMaterialInterface(), AdParameters::saveParametersToShowFile(), Cgad::saveShowFileComments(), AdParameters::setParameters(), AdParameters::setPdeParameters(), Cgad::setPlotTitle(), Cgins::turbulenceModels(), and Cgad::writeParameterSummary().

const real kappa2 =kappa*kappa

Referenced by Cgsm::assignSpecialInitialConditions().

const real kappa3 =kappa*kappa*kappa

Referenced by Cgsm::assignSpecialInitialConditions().

const real kappa4 =kappa*kappa*kappa*kappa

Referenced by Cgsm::assignSpecialInitialConditions().

const real& lambda = parameters.dbase.get<real>("lambda")

Referenced by Cgsm::advance(), Cgsm::advanceFOS(), Cgsm::advanceSOS(), Cgsm::assignAnnulusEigenfunction(), Cgsm::assignBoundaryConditionsFOS(), Cgsm::assignBoundaryConditionsSOS(), Cgsm::assignGaussianPulseInitialConditions(), Cgsm::assignHempInitialConditions(), Maxwell::assignInitialConditions(), Cgsm::assignParabolicInitialConditions(), Cgsm::assignSpecialInitialConditions(), Cgsm::assignTwilightZoneInitialConditions(), Parameters::buildMaterialParametersDialog(), Cgsm::Cgsm(), Cgsm::checkDisplacementAndStress(), Maxwell::computeTimeStep(), FOR_3D(), Cgsm::getEnergy(), Cgsm::getErrors(), Maxwell::getErrors(), Cgsm::getForcing(), Cgsm::getInitialConditions(), AsfParameters::getNormalForce(), CnsParameters::getNormalForce(), FlowSolutions::getShockElasticPistonFlow(), Cgsm::getTimeStep(), RigidBodyMotion::getTimeStepEstimate(), Parameters::getUserDefinedKnownSolution(), Cgasf::getUt(), Cgsm::getVelocityAndStress(), if(), Cgmp::interfaceProjection(), Cgsm::plot(), Parameters::readFromAShowFile(), SmParameters::saveParametersToShowFile(), SmParameters::setPdeParameters(), Cgsm::setupGrids(), SmParameters::updateToMatchGrid(), Parameters::updateUserDefinedKnownSolution(), Cgsm::userDefinedForcing(), and Cgsm::userDefinedInitialConditions().

const RealArray& lambdaGrid = parameters.dbase.get<RealArray>("lambdaGrid")

Referenced by Cgsm::advance(), Cgsm::advanceFOS(), Cgsm::advanceSOS(), Cgsm::assignAnnulusEigenfunction(), Cgsm::assignBoundaryConditionsFOS(), Cgsm::assignBoundaryConditionsSOS(), Cgsm::assignGaussianPulseInitialConditions(), Cgsm::assignHempInitialConditions(), Cgsm::assignParabolicInitialConditions(), Cgsm::assignSpecialInitialConditions(), Cgsm::assignTwilightZoneInitialConditions(), Cgsm::Cgsm(), Cgsm::checkDisplacementAndStress(), Cgsm::computeDissipation(), Cgsm::getEnergy(), Cgsm::getErrors(), Cgsm::getForcing(), Cgsm::getTimeStep(), Cgsm::getVelocityAndStress(), SmParameters::saveParametersToShowFile(), SmParameters::setPdeParameters(), Cgsm::setupGrids(), SmParameters::updateToMatchGrid(), and Cgsm::userDefinedForcing().

const int m = (int)data[2]

const real& mu = parameters.dbase.get<real>("mu")

Referenced by Cgasf::addForcing(), Cgcns::addForcing(), Cgasf::addForcingToPressureEquation(), Cgsm::advance(), Cgsm::advanceFOS(), Cgsm::advanceSOS(), Cgcns::applyBoundaryConditions(), Cgcns::applyBoundaryConditionsForImplicitTimeStepping(), Cgsm::assignAnnulusEigenfunction(), RadiationBoundaryCondition::assignBoundaryConditions(), Cgsm::assignBoundaryConditionsFOS(), Cgsm::assignBoundaryConditionsSOS(), Cgsm::assignGaussianPulseInitialConditions(), Cgsm::assignHempInitialConditions(), Cgsm::assignParabolicInitialConditions(), Cgsm::assignSpecialInitialConditions(), Cgsm::assignTwilightZoneInitialConditions(), Parameters::buildMaterialParametersDialog(), Cgsm::Cgsm(), Cgsm::checkDisplacementAndStress(), Cgcns::formImplicitTimeSteppingMatrix(), Cgsm::getEnergy(), Cgsm::getErrors(), Cgsm::getForcing(), Cgsm::getInitialConditions(), AsfParameters::getNormalForce(), CnsParameters::getNormalForce(), getRayleighSpeed(), FlowSolutions::getShockElasticPistonFlow(), FlowSolutions::getSupersonicExpandingFlow(), Cgsm::getTimeStep(), Cgcns::getTimeSteppingEigenvalue(), Cgasf::getTimeSteppingEigenvalue(), Parameters::getUserDefinedKnownSolution(), Cgcns::getUt(), Cgasf::getUt(), Cgsm::getVelocityAndStress(), if(), Cgmp::interfaceProjection(), muFromPmu(), Cgsm::plot(), Parameters::readFromAShowFile(), SmParameters::saveParametersToShowFile(), Cgcns::saveShowFileComments(), SmParameters::setPdeParameters(), Cgcns::setPlotTitle(), Cgsm::setupGrids(), Cgasf::solveForAllSpeedPressure(), Cgasf::solveForTimeIndependentVariables(), SmParameters::updateToMatchGrid(), Parameters::updateUserDefinedKnownSolution(), Cgsm::userDefinedForcing(), and Cgsm::userDefinedInitialConditions().

const RealArray& muGrid = parameters.dbase.get<RealArray>("muGrid")

Referenced by Cgsm::advance(), Cgsm::advanceFOS(), Cgsm::advanceSOS(), Cgsm::assignAnnulusEigenfunction(), Cgsm::assignBoundaryConditionsFOS(), Cgsm::assignBoundaryConditionsSOS(), Cgsm::assignGaussianPulseInitialConditions(), Cgsm::assignHempInitialConditions(), Cgsm::assignParabolicInitialConditions(), Cgsm::assignSpecialInitialConditions(), Cgsm::assignTwilightZoneInitialConditions(), Cgsm::Cgsm(), Cgsm::checkDisplacementAndStress(), Cgsm::computeDissipation(), Cgsm::getEnergy(), Cgsm::getErrors(), Cgsm::getForcing(), Cgsm::getTimeStep(), Cgsm::getVelocityAndStress(), SmParameters::saveParametersToShowFile(), SmParameters::setPdeParameters(), Cgsm::setupGrids(), SmParameters::updateToMatchGrid(), and Cgsm::userDefinedForcing().

const int n = (int)data[1]

real omega

Referenced by Cgins::advanceLineSolveNew(), Cgins::advanceLineSolveOld(), Cgasf::allSpeedImplicitTimeStep(), Cgsm::assignAnnulusEigenfunction(), Cgsm::assignBoundaryConditionsFOS(), Cgsm::assignBoundaryConditionsSOS(), Cgsm::assignHempInitialConditions(), Cgmp::assignInterfaceRightHandSide(), Cgmp::assignInterfaceRightHandSideOld(), Cgsm::assignSpecialInitialConditions(), Cgsm::buildForcingOptionsDialog(), Cgsm::buildVariableDissipation(), RigidBodyMotion::getAngularVelocities(), MovingGrids::getBoundaryAcceleration(), RigidBodyMotion::getCoordinates(), Cgsm::getForcing(), Cgsm::getForcingOption(), MovingGrids::getGridVelocity(), Cgsm::getInitialConditions(), RigidBodyMotion::getRotationMatrix(), DomainSolver::getTimeDependentBoundaryConditions(), MovingGrids::getUserDefinedBoundaryAcceleration(), MovingGrids::getUserDefinedGridVelocity(), Parameters::getUserDefinedKnownSolution(), MovingGrids::gridAccelerationBC(), if(), main(), MovingGrids::moveGrids(), DeformingBodyMotion::regenerateComponentGrids(), SmParameters::setTwilightZoneFunction(), CnsParameters::setTwilightZoneFunction(), InsParameters::setTwilightZoneFunction(), Parameters::setTwilightZoneParameters(), Cgsm::setupUserDefinedForcing(), Cgasf::solveForAllSpeedPressure(), MovingGrids::update(), Parameters::updateUserDefinedKnownSolution(), MovingGrids::updateUserDefinedMotion(), Cgsm::userDefinedBoundaryValues(), Cgsm::userDefinedForcing(), MovingGrids::userDefinedGridAccelerationBC(), and MovingGrids::userDefinedTransformMotion().

const int pc = parameters.dbase.get<int >("pc")

const real rad = data[3]

Referenced by Cgcns::addForcing(), Cgsm::assignSpecialInitialConditions(), DomainSolver::computeBodyForcing(), Maxwell::defineRegionsAndBodies(), Parameters::defineVariableBoundaryValues(), Maxwell::getChargeDensity(), Cgsm::getInitialConditions(), if(), DeformingBodyMotion::initialize(), main(), DeformingBodyMotion::regenerateComponentGrids(), DomainSolver::setVariableBoundaryValues(), DomainSolver::setVariableMaterialProperties(), Cgsm::userDefinedBoundaryValues(), DomainSolver::userDefinedBoundaryValues(), userDefinedErrorEstimator(), DomainSolver::userDefinedForcing(), and DomainSolver::userDefinedInitialConditions().

const real& rho =parameters.dbase.get<real>("rho")

Referenced by CnsParameters::assignParameterValues(), Parameters::assignParameterValues(), Parameters::buildMaterialParametersDialog(), Reactions::checkChemicalSourceJacobian(), Reactions::computeEigenvaluesOfTheChemicalSourceJacobian(), CnsParameters::conservativeToPrimitive(), Maxwell::getChargeDensity(), FlowSolutions::getSupersonicExpandingFlow(), CnsParameters::primitiveToConservative(), Parameters::readFromAShowFile(), SmParameters::setPdeParameters(), DomainSolver::sizeOf(), and Reactions::solveImplicitForYGivenRTP().

const int s11c = parameters.dbase.get<int >("s11c")

const int s12c = parameters.dbase.get<int >("s12c")

const int s13c = parameters.dbase.get<int >("s13c")

const int s21c = parameters.dbase.get<int >("s21c")

const int s22c = parameters.dbase.get<int >("s22c")

const int s23c = parameters.dbase.get<int >("s23c")

const int s31c = parameters.dbase.get<int >("s31c")

const int s32c = parameters.dbase.get<int >("s32c")

const int s33c = parameters.dbase.get<int >("s33c")

const int v2c = parameters.dbase.get<int >("v2c")

const int v3c = parameters.dbase.get<int >("v3c")

const int vibrationClass = (int)data[0]

Referenced by Cgsm::assignSpecialInitialConditions(), and Cgsm::getInitialConditions().