bcExtended3d4.h File Reference

Go to the source code of this file.

Functions
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3	ut (also for v and w)!deltaFu
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3 a1	D (i1, i2-2)
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3 a1 a2	D (i1-2, i2)

Variables
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined	!c11
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined	c22
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined	c33
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined	c1
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined	c2
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3	!urr
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3	uss
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3	utt
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3	ur
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3	us
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3	deltaFv
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3	deltaFw
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3	g2f = RHS for extrapolation
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3 a1 a2	!DeltaU
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3 a1 a2	a2
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3 a1 a2 a3 **	a11c
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3 a1 a2 a3 a3 a3 a3 a2 and a3 at the corner	a11
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3 a1 a2 a3 a3 a3 a3 a2 and a3 at the corner a2	Lu
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3 a1 a2 a3 a3 a3 a3 a2 and a3 at the corner a2 a2 D r4	u
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3 a1 a2 a3 a3 a3 a3 a2 and a3 at the corner a2 a2 D r4	x2
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3 a1 a2 a3 a3 a3 a3 a2 and a3 at the corner a2 a2 D r4	x3
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3 a1 a2 a3 a3 a3 a3 a2 and a3 at the corner a2 a2 D r4	x4
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3 a1 a2 a3 a3 a3 a3 a2 and a3 at the corner a2 a2 D r4 i2	is2
!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3 a1 a2 a3 a3 a3 a3 a2 and a3 at the corner a2 a2 D r4 i2 i3	is3

Function Documentation

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 a1 D	(	i1	,
		i2-	2
	)

Referenced by Cgsm::getForcing(), and getRayleighSpeed().

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 a1 a2 D	(	i1-	2,
		i2
	)

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 ut

(

also for v and

w

)

Referenced by Cgad::addForcing(), Cgins::addForcing(), Cgasf::addForcing(), Cgcns::addForcing(), Cgins::determineErrors(), DomainSolver::determineErrors(), Cgmp::getInterfaceResiduals(), Cgcns::getUt(), Cgasf::getUt(), Cgins::getUt(), ogfTaylor(), DomainSolver::takeTimeStepFE(), and Cgins::turbulenceModels().

Variable Documentation

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined !c11

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 a1 a2 !DeltaU

Initial value:

= c11*urr+c22*uss+c33*utt+c1*ur+c2*us+c3*ut - deltaFu
      DeltaV = c11*vrr+c22*vss+c33*vtt+c1*vr+c2*vs+c3*vt - deltaFv
      DeltaW = c11*wrr+c22*wss+c33*wtt+c1*wr+c2*ws+c3*wt - deltaFw

! ** decompose point u(i1-is1,i2-is2,i3-is3) into components along a1

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 !urr

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 a1 a2 a3 a3 a3 a3 a2 and a3 at the corner a11

Referenced by Maxwell::computeTimeStep(), Cgsm::getTimeStep(), Cgmp::initializeInterfaceBoundaryConditions(), MovingGrids::rigidBodyMotion(), RigidBodyMotion::takeStepImplicitRungeKutta(), and PistonMotion::timeStep().

!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3 a1 a2 a3 a3 a3 a3 ** a11c

Initial value:

=A11D3(i1-is1,i2-is2,i3-is3)
      a12c=A12D3(i1-is1,i2-is2,i3-is3)
      a13c=A13D3(i1-is1,i2-is2,i3-is3)
      a21c=A21D3(i1-is1,i2-is2,i3-is3)
      a22c=A22D3(i1-is1,i2-is2,i3-is3)
      a23c=A23D3(i1-is1,i2-is2,i3-is3)
      a31c=A31D3(i1-is1,i2-is2,i3-is3)
      a32c=A32D3(i1-is1,i2-is2,i3-is3)
      a33c=A33D3(i1-is1,i2-is2,i3-is3)

      a1a1=a11c*a11c+a12c*a12c+a13c*a13c
      a1a2=a11c*a21c+a12c*a22c+a13c*a23c
      a1a3=a11c*a31c+a12c*a32c+a13c*a33c
      a2a3=a21c*a31c+a22c*a32c+a23c*a33c
      a3a3=a31c*a31c+a32c*a32c+a33c*a33c
! The tangential component is assumed valid:
      a1Dotu1=a11c*u(i1-is1,i2-is2,i3-is3,ex)+a12c*u(i1-is1,i2-is2,i3-is3,ey)+a13c*u(i1-is1,i2-is2,i3-is3,ez)
      a3Dotu1=a31c*u(i1-is1,i2-is2,i3-is3,ex)+a32c*u(i1-is1,i2-is2,i3-is3,ey)+a33c*u(i1-is1,i2-is2,i3-is3,ez)
 ! u(i1-is1,i2-is2,i3-is3,k) = b1[k]*x1 +g1[k]
      b11 =-a11c*(-a1a2*a3a3+a1a3*a2a3)/(-a1a1*a3a3+a1a3**2)+a21c-a31c*(-a1a1*a2a3+a1a3*a1a2)/(-a1a1*a3a3+a1a3**2)
      b12 =-a12c*(-a1a2*a3a3+a1a3*a2a3)/(-a1a1*a3a3+a1a3**2)+a22c-a32c*(-a1a1*a2a3+a1a3*a1a2)/(-a1a1*a3a3+a1a3**2)
      b13 =-a13c*(-a1a2*a3a3+a1a3*a2a3)/(-a1a1*a3a3+a1a3**2)+a23c-a33c*(-a1a1*a2a3+a1a3*a1a2)/(-a1a1*a3a3+a1a3**2)
      g11 =-(-a11c*a1a3*a3Dotu1+a11c*a1Dotu1*a3a3-a31c*a1a3*a1Dotu1+a31c*a1a1*a3Dotu1)/(-a1a1*a3a3+a1a3**2)
      g12 =-(-a12c*a1a3*a3Dotu1+a12c*a1Dotu1*a3a3-a32c*a1a3*a1Dotu1+a32c*a1a1*a3Dotu1)/(-a1a1*a3a3+a1a3**2)
      g13 =(a13c*a1a3*a3Dotu1-a13c*a1Dotu1*a3a3+a33c*a1a3*a1Dotu1-a33c*a1a1*a3Dotu1)/(-a1a1*a3a3+a1a3**2)

! ** decompose point u(i1-2*is1,i2-2*is2,i3-2*is3) into components along a1

!results from bcExtended3d4 maple!Assign values on the extended boundary next to two PEC boundaries!Here we assume the following are defined c3 a1 a2 a3 a3 a3 a2

Referenced by Cgins::advanceLineSolveOld(), Cgasf::allSpeedImplicitTimeStep(), Cgasf::applyBoundaryConditions(), Maxwell::assignInitialConditions(), Maxwell::computeIntensity(), F90_ID(), Maxwell::getErrors(), Parameters::getUserDefinedKnownSolution(), and DomainSolver::variableTimeStepBoundaryInterpolation().

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c1

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c2

Referenced by Maxwell::advanceC(), Cgins::advanceLineSolveOld(), Maxwell::applyDSIBC(), GridEvolution::getAcceleration(), RigidBodyMotion::getAddedMassMatrices(), RigidBodyMotion::getForce(), RigidBodyMotion::getForceInternal(), Parameters::getUserDefinedKnownSolution(), Maxwell::initializePlaneMaterialInterface(), RigidBodyMotion::integrate(), InterfaceTransfer::internalInterpolate(), DomainSolver::parabolicInflow(), Maxwell::setupDSICoefficients(), RigidBodyMotion::takeStepExtrapolation(), RigidBodyMotion::takeStepImplicitRungeKutta(), PistonMotion::timeStep(), and userDefinedErrorEstimator().

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c22

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c33

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 deltaFv

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 deltaFw

Initial value:

= RHS for Delta(u,v,w)                
 !    g1f

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 g2f = RHS for extrapolation

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 a1 a2 a3 a3 a3 a3 a2 and a3 at the corner a2 a2 D r4 i2 is2

Referenced by Cgins::advanceLineSolveNew(), applyFourthOrderBoundaryConditionsOld(), Maxwell::assignBoundaryConditions(), Cgins::assignLineSolverBoundaryConditions(), Cgins::computeTurbulenceQuantities(), insbe2(), Cgins::lineSolverBoundaryConditions(), and rigidBodyBoundaryProjection().

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 a1 a2 a3 a3 a3 a3 a2 and a3 at the corner a2 a2 D r4 i2 i3 is3

Referenced by Cgins::advanceLineSolveNew(), applyFourthOrderBoundaryConditionsOld(), Cgins::assignLineSolverBoundaryConditions(), Cgins::computeTurbulenceQuantities(), Cgins::lineSolverBoundaryConditions(), and rigidBodyBoundaryProjection().

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 a1 a2 a3 a3 a3 a3 a2 and a3 at the corner a2 Lu

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 a1 a2 a3 a3 a3 a3 a2 and a3 at the corner a2 a2 D r4 u

Initial value:

= g2f  a1.D+s4 u = g1f
 ! e1x := dd11*x1+dd12*x2+dd13*x3+dd14*x4+ f1x
 ! e2x := dd21*x1+dd22*x2+dd23*x3+dd24*x4+ f2x
 ! e3x := dd31*x1+dd32*x2+dd33*x3+dd34*x4+ f3x
 ! e4x := dd41*x1+dd42*x2+dd43*x3+dd44*x4+ f4x
      dd11=cc14a*b11+cc15a*b12+cc16a*b13
      dd12=cc11a*b21+cc12a*b22+cc13a*b23
      dd13=cc14b*b31+cc15b*b32+cc16b*b33
      dd14=cc11b*b41+cc12b*b42+cc13b*b43
      dd21=cc24a*b11+cc25a*b12+cc26a*b13
      dd22=cc21a*b21+cc22a*b22+cc23a*b23
      dd23=cc24b*b31+cc25b*b32+cc26b*b33
      dd24=cc21b*b41+cc22b*b42+cc23b*b43
      dd31=-4*a21*b11-4*a22*b12-4*a23*b13
      dd32=a21*b21+a22*b22+a23*b23
      dd33=0
      dd34=0
      dd41=0
      dd42=0
      dd43=-4*a11*b31-4*a12*b32-4*a13*b33
      dd44=a11*b41+a12*b42+a13*b43

      f1x=cc11a*g21+cc12a*g22+cc13a*g23+cc14a*g11+cc15a*g12+cc16a*g13+cc11b*g41+cc12b*g42+cc13b*g43+cc14b*g31+cc15b*g32+cc16b*g33+f1
      f2x=cc21a*g21+cc22a*g22+cc23a*g23+cc24a*g11+cc25a*g12+cc26a*g13+cc21b*g41+cc22b*g42+cc23b*g43+cc24b*g31+cc25b*g32+cc26b*g33+f2
      f3x=a21*g21+a22*g22+a23*g23-4*a21*g11-4*a22*g12-4*a23*g13+f3
      f4x=a11*g41+a12*g42+a13*g43-4*a11*g31-4*a12*g32-4*a13*g33+f4

!  solution x1

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 ur

Referenced by Cgcns::applyBoundaryConditions(), Cgins::assignLineSolverBoundaryConditions(), Cgins::getAugmentedSolution(), DomainSolver::getAugmentedSolution(), Cgsm::getErrors(), Maxwell::getErrors(), and Cgins::lineSolverBoundaryConditions().

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 us

Referenced by Cgmp::getInterfaceResiduals(), InterfaceTransfer::internalInterpolate(), and CnsParameters::setParameters().

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 uss

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 utt

Referenced by Cgsm::getForcing(), and ogfTaylor().

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 a1 a2 a3 a3 a3 a3 a2 and a3 at the corner a2 a2 D r4 x2

Referenced by DeformingBodyMotion::advanceElasticShell(), applyFourthOrderBoundaryConditionsOld(), Cgsm::assignAnnulusEigenfunction(), Maxwell::assignInterfaceBoundaryConditions(), GridEvolution::getAcceleration(), Maxwell::getChargeDensity(), Cgsm::getInitialConditions(), GridEvolution::getVelocity(), DeformingBodyMotion::initialize(), DeformingBodyMotion::integrate(), main(), DeformingBodyMotion::regenerateComponentGrids(), Cgsm::setupUserDefinedForcing(), and Cgsm::userDefinedForcing().

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 a1 a2 a3 a3 a3 a3 a2 and a3 at the corner a2 a2 D r4 x3

Referenced by DeformingBodyMotion::advanceElasticShell(), Cgsm::assignAnnulusEigenfunction(), GridEvolution::getAcceleration(), and MovingGrids::moveGrids().

! results from bcExtended3d4 maple ! Assign values on the extended boundary next to two PEC boundaries ! Here we assume the following are defined c3 a1 a2 a3 a3 a3 a3 a2 and a3 at the corner a2 a2 D r4 x4

Referenced by Cgsm::assignAnnulusEigenfunction(), and GridEvolution::getAcceleration().