Implements TGA algorithm to update solution to heat equation. More...
#include <LevelTGA.H>
Inheritance diagram for LevelTGA:
Public Member Functions | |
| LevelTGA (const Vector< DisjointBoxLayout > &a_grids, const Vector< int > &a_refRat, const ProblemDomain &a_level0Domain, RefCountedPtr< AMRLevelOpFactory< LevelData< FArrayBox > > > &a_opFact, const RefCountedPtr< AMRMultiGrid< LevelData< FArrayBox > > > &a_solver) | |
| full constructor More... | |
| virtual | ~LevelTGA () |
| destructor More... | |
| void | computeDiffusion (LevelData< FArrayBox > &a_DiffusiveTerm, LevelData< FArrayBox > &a_phiOld, LevelData< FArrayBox > &a_src, LevelFluxRegister *a_FineFluxRegPtr, LevelFluxRegister *a_crseFluxRegPtr, const LevelData< FArrayBox > *a_crsePhiOldPtr, const LevelData< FArrayBox > *a_crsePhiNewPtr, Real a_oldTime, Real a_crseOldTime, Real a_crseNewTime, Real a_dt, int a_level) |
| computes time-centered diffusive term for explicit updating More... | |
| void | updateSoln (LevelData< FArrayBox > &a_phiNew, LevelData< FArrayBox > &a_phiOld, LevelData< FArrayBox > &a_src, LevelFluxRegister *a_fineFluxRegPtr, LevelFluxRegister *a_crseFluxRegPtr, const LevelData< FArrayBox > *a_crsePhiOldPtr, const LevelData< FArrayBox > *a_crsePhiNewPtr, Real a_oldTime, Real a_crseOldTime, Real a_crseNewTime, Real a_dt, int a_level, bool a_zeroPhi=true, int a_fluxStartComponent=0) |
| do solve More... | |
| void | setSourceGhostCells (LevelData< FArrayBox > &a_src, const DisjointBoxLayout &a_grids, int a_level) |
 Public Member Functions inherited from BaseLevelTGA< LevelData< FArrayBox >, FluxBox, LevelFluxRegister > | |
| BaseLevelTGA (const Vector< DisjointBoxLayout > &a_grids, const Vector< int > &a_refRat, const ProblemDomain &a_level0Domain, RefCountedPtr< AMRLevelOpFactory< LevelData< FArrayBox > > > &a_opFact, const RefCountedPtr< AMRMultiGrid< LevelData< FArrayBox > > > &a_solver) | |
| virtual | ~BaseLevelTGA () |
| Destructor, called after destructors of BaseLevelTGA subclasses. More... | |
| void | updateSoln (LevelData< FArrayBox > &a_phiNew, LevelData< FArrayBox > &a_phiOld, LevelData< FArrayBox > &a_src, LevelData< FluxBox > &a_flux, LevelFluxRegister *a_fineFluxRegPtr, LevelFluxRegister *a_crseFluxRegPtr, const LevelData< FArrayBox > *a_crsePhiOldPtr, const LevelData< FArrayBox > *a_crsePhiNewPtr, Real a_oldTime, Real a_crseOldTime, Real a_crseNewTime, Real a_dt, int a_level, bool a_zeroPhi=true, bool a_rhsAlreadyKappaWeighted=false, int a_fluxStartComponent=0) |
| void | computeDiffusion (LevelData< FArrayBox > &a_diffusiveTerm, LevelData< FArrayBox > &a_phiOld, LevelData< FArrayBox > &a_src, LevelData< FluxBox > &a_flux, LevelFluxRegister *a_fineFluxRegPtr, LevelFluxRegister *a_crseFluxRegPtr, const LevelData< FArrayBox > *a_crsePhiOldPtr, const LevelData< FArrayBox > *a_crsePhiNewPtr, Real a_oldTime, Real a_crseOldTime, Real a_crseNewTime, Real a_dt, int a_level, bool a_zeroPhi=true, bool a_rhsAlreadyKappaWeighted=false) |
| Public Member Functions inherited from BaseLevelHeatSolver< LevelData< FArrayBox >, FluxBox, LevelFluxRegister > | |
| BaseLevelHeatSolver (const Vector< DisjointBoxLayout > &a_grids, const Vector< int > &a_refRat, const ProblemDomain &a_level0Domain, RefCountedPtr< AMRLevelOpFactory< LevelData< FArrayBox > > > &a_opFact, const RefCountedPtr< AMRMultiGrid< LevelData< FArrayBox > > > &a_solver) | |
| virtual | ~BaseLevelHeatSolver () |
| Destructor, called after destructors of BaseLevelHeatSolver subclasses. More... | |
| virtual void | updateSoln (LevelData< FArrayBox > &a_phiNew, LevelData< FArrayBox > &a_phiOld, LevelData< FArrayBox > &a_src, LevelFluxRegister *a_fineFluxRegPtr, LevelFluxRegister *a_crseFluxRegPtr, const LevelData< FArrayBox > *a_crsePhiOldPtr, const LevelData< FArrayBox > *a_crsePhiNewPtr, Real a_oldTime, Real a_crseOldTime, Real a_crseNewTime, Real a_dt, int a_level, bool a_zeroPhi=true, bool a_rhsAlreadyKappaWeighted=false, int a_fluxStartComponent=0) |
| virtual void | applyOperator (LevelData< FArrayBox > &a_ans, const LevelData< FArrayBox > &a_phi, const LevelData< FArrayBox > *a_phiC, int a_level, Real a_alpha, Real a_beta, bool a_applyBC) |
| void | applyHelm (LevelData< FArrayBox > &a_ans, const LevelData< FArrayBox > &a_phi, const LevelData< FArrayBox > *a_phiC, int a_level, Real a_mu, Real a_dt, bool a_homogeneousBC) |
| void | incrementFlux (LevelData< FluxBox > &a_diffusiveFlux, LevelData< FArrayBox > &a_phi, int a_level, Real a_mu, Real a_dt, Real a_sign, bool a_setToZero) |
| void | solveHelm (LevelData< FArrayBox > &a_phi, LevelData< FArrayBox > &a_phiC, LevelData< FArrayBox > &a_rhs, int a_level, Real a_mu, Real a_dt, bool a_zeroPhi=true) |
| void | resetSolverAlphaAndBeta (const Real &a_alpha, const Real &a_beta) |
| virtual LevelTGAHelmOp< LevelData< FArrayBox >, FluxBox > * | newOp (const ProblemDomain &a_indexSpace, RefCountedPtr< AMRLevelOpFactory< LevelData< FArrayBox > > > &a_opFact) |
| int | size () const |
| Returns the number of grid levels on which this integrator operates. More... | |
Additional Inherited Members | |
| Protected Member Functions inherited from BaseLevelHeatSolver< LevelData< FArrayBox >, FluxBox, LevelFluxRegister > | |
| void | timeInterp (LevelData< FArrayBox > &a_data, const LevelData< FArrayBox > &a_oldData, const LevelData< FArrayBox > &a_newData, Real a_time, Real a_oldTime, Real a_newTime, int a_level) |
| Protected Attributes inherited from BaseLevelTGA< LevelData< FArrayBox >, FluxBox, LevelFluxRegister > | |
| Real | m_mu1 |
| Real | m_mu2 |
| Real | m_mu3 |
| Real | m_mu4 |
| Real | m_r1 |
| Protected Attributes inherited from BaseLevelHeatSolver< LevelData< FArrayBox >, FluxBox, LevelFluxRegister > | |
| Vector< DisjointBoxLayout > | m_grids |
| The disjoint box layouts at every AMR grid level. More... | |
| Vector< int > | m_refRat |
| The refinement ratios between AMR grid levels. More... | |
| ProblemDomain | m_level0Domain |
| The coarsest domain on which the Helmholtz equation is integrated. More... | |
| Vector< LevelTGAHelmOp< LevelData< FArrayBox >, FluxBox > * > | m_ops |
| RefCountedPtr< AMRMultiGrid< LevelData< FArrayBox > > > | m_solver |
| The multigrid solver used to solve the Helmholtz equation. More... | |
Detailed Description
Implements TGA algorithm to update solution to heat equation.
The LevelTGA class implements the Runge-Kutta-based approach to solving the heat equation due to Twizell, Gumel, and Arigu.
Constructor & Destructor Documentation
◆ LevelTGA()
|
inline |
full constructor
◆ ~LevelTGA()
|
inlinevirtual |
destructor
References computeDiffusion(), setSourceGhostCells(), and updateSoln().
Member Function Documentation
◆ computeDiffusion()
| void LevelTGA::computeDiffusion | ( | LevelData< FArrayBox > & | a_DiffusiveTerm, |
| LevelData< FArrayBox > & | a_phiOld, | ||
| LevelData< FArrayBox > & | a_src, | ||
| LevelFluxRegister * | a_FineFluxRegPtr, | ||
| LevelFluxRegister * | a_crseFluxRegPtr, | ||
| const LevelData< FArrayBox > * | a_crsePhiOldPtr, | ||
| const LevelData< FArrayBox > * | a_crsePhiNewPtr, | ||
| Real | a_oldTime, | ||
| Real | a_crseOldTime, | ||
| Real | a_crseNewTime, | ||
| Real | a_dt, | ||
| int | a_level | ||
| ) |
computes time-centered diffusive term for explicit updating
compute time-centered L(phi) for use in subsequent update operations. In this case, we do solve for phiNew, then subtract source and old phi back out to get L(phi).
Referenced by LevelBackwardEuler::~LevelBackwardEuler(), and ~LevelTGA().
◆ updateSoln()
| void LevelTGA::updateSoln | ( | LevelData< FArrayBox > & | a_phiNew, |
| LevelData< FArrayBox > & | a_phiOld, | ||
| LevelData< FArrayBox > & | a_src, | ||
| LevelFluxRegister * | a_fineFluxRegPtr, | ||
| LevelFluxRegister * | a_crseFluxRegPtr, | ||
| const LevelData< FArrayBox > * | a_crsePhiOldPtr, | ||
| const LevelData< FArrayBox > * | a_crsePhiNewPtr, | ||
| Real | a_oldTime, | ||
| Real | a_crseOldTime, | ||
| Real | a_crseNewTime, | ||
| Real | a_dt, | ||
| int | a_level, | ||
| bool | a_zeroPhi = true, |
||
| int | a_fluxStartComponent = 0 |
||
| ) |
do solve
update phi to a_phiNew
Referenced by LevelBackwardEuler::~LevelBackwardEuler(), and ~LevelTGA().
◆ setSourceGhostCells()
|
virtual |
Sets the value of the source term in the Helmholtz equation on ghost cells. This method should be overridden by subclasses of BaseLevelTGA.
- Parameters
-
a_src The source term in the Helmholtz equation whose ghost cell values are to be set. a_dbl The disjoint box layout that indexes a_src. a_level The grid level at which the ghost cell values are to be set.
Implements BaseLevelTGA< LevelData< FArrayBox >, FluxBox, LevelFluxRegister >.
Referenced by LevelBackwardEuler::~LevelBackwardEuler(), and ~LevelTGA().
The documentation for this class was generated from the following file:
