Abstract class for evolving the minimal plasma model in time. More...

#include <CD_CdrPlasmaStepper.H>

Inheritance diagram for Physics::CdrPlasma::CdrPlasmaStepper:

[legend]

Collaboration diagram for Physics::CdrPlasma::CdrPlasmaStepper:

[legend]

Public Member Functions
	CdrPlasmaStepper ()
	Empty constructor.

	CdrPlasmaStepper (RefCountedPtr< CdrPlasmaPhysics > &a_physics)
	Full constructor.

	~CdrPlasmaStepper ()
	Destructor.

virtual void	setupSolvers () override
	Instantiate solvers.

virtual void	allocate () override
	Allocate data for the time stepper and solvers.

virtual void	initialData () override
	Fill all solvers with initial data.

virtual void	postInitialize () override
	Post-initialize operations to be performed after filling solvers with initialData.

virtual void	postCheckpointSetup () override
	Post-initialize operations to be performed after filling solvers with data read from checkpoint files.

virtual void	registerRealms () override
	Register realms to be used for the simulation.

virtual void	registerOperators () override
	Register operators to be used for the simulation.

virtual void	parseRuntimeOptions () override=0
	Parse runtime options. Subclasses must implement this.

virtual void	prePlot () override
	Compute some thing that go into plot files.

virtual int	getNumberOfPlotVariables () const override
	Get the number of plot variables for this time stepper.

virtual Vector< std::string >	getPlotVariableNames () const override
	Get plot variable names.

virtual void	writePlotData (LevelData< EBCellFAB > &a_output, int &a_icomp, const std::string a_outputRealm, const int a_level) const override
	Write plot data to output holder.

virtual Real	computeDt () override=0
	Compute the time step.

virtual Real	advance (const Real a_dt) override=0
	Advance method, advances equations.

virtual void	synchronizeSolverTimes (const int a_step, const Real a_time, const Real a_dt) override
	Synchronize solver times.

virtual void	printStepReport () override
	Print a step report.

virtual void	preRegrid (const int a_lmin, const int a_oldFinestLevel) override
	Pre-regrid function – this is the one called by Driver.

virtual void	preRegridInternals (const int a_lmin, const int a_oldFinestLevel)
	Do a preRegrid operation for internal storage (if needed).

virtual void	regrid (const int a_lmin, const int a_oldFinestLevel, const int a_newFinestLevel) override
	Regrid method.

virtual void	postRegrid () override
	Perform post-regrid operations.

virtual void	computeJ (EBAMRCellData &a_J) const
	Compute the current density.

virtual void	computeElectricField (MFAMRCellData &a_electricFieldCell, const MFAMRCellData &a_potential) const
	Compute the cell-centered electric field on both phases.

virtual void	computeElectricField (EBAMRCellData &a_electricFieldCell, const phase::which_phase a_phase) const
	Compute the cell-centered electric field on a specific phase using whatever is available in m_fieldSolver.

virtual void	computeElectricField (EBAMRCellData &a_E, const phase::which_phase a_phase, const MFAMRCellData &a_potential) const
	Compute the cell-centered electric field on a specific phase.

virtual void	computeElectricField (EBAMRFluxData &a_electricFieldFace, const phase::which_phase a_phase, const EBAMRCellData &a_electricFieldCell) const
	Compute the face-centered electric field. This will only do the normal component.

virtual void	computeElectricField (EBAMRIVData &a_electricFieldEB, const phase::which_phase a_phase, const EBAMRCellData &a_electricFieldCell) const
	Compute the EB-centered electric field.

virtual void	computeMaxElectricField (Real &a_maxElectricField, const phase::which_phase a_phase)
	Compute the maximum of the electric field.

virtual Real	getTime () const
	Return time.

virtual bool	stationaryRTE ()
	Check if RTE solvers are stationary.

virtual void	allocateInternals ()=0
	Allocate internal storage. Must be implemented by subclasses.

virtual void	deallocateInternals ()=0
	Deallocate internals. Must be implemented by subclasses.

virtual void	parseOptions ()=0
	Parse TimeStepper class options. Must be implemented by subclasses.

virtual void	deallocate ()
	Deallocation function. This will deallocate internal storage in the subclasses as well as the solvers.

virtual void	computeSpaceChargeDensity ()
	Compute the centroid-centered space charge density by using the data inside the CDR solvers.

virtual void	computeSpaceChargeDensity (MFAMRCellData &a_rho, const Vector< EBAMRCellData * > &a_cdrDensities)
	Compute the centroid-centered space charge density.

virtual void	computeCellConductivity (EBAMRCellData &a_cellConductivity) const
	Compute the cell-centered conductivity.

virtual void	computeFaceConductivity (EBAMRFluxData &a_conductivityFace, EBAMRIVData &a_conductivityEB, const EBAMRCellData &a_conductivityCell) const
	Compute the face-centered conductivity from a cell-centered conductivity.

virtual void	setupSemiImplicitPoisson (const Real a_dt)
	Set up a semi-implicit Poisson solver.

virtual void	setupSemiImplicitPoisson (const EBAMRFluxData &a_conductivityFace, const EBAMRIVData &a_conductivityEB, const Real a_factor)
	Set up a semi-implicit Poisson solver.

virtual bool	solvePoisson ()
	Solve the electrostatic Poisson equation.

virtual bool	solvePoisson (MFAMRCellData &a_potential, MFAMRCellData &a_rho, const Vector< EBAMRCellData * > a_cdrDensities, const EBAMRIVData &a_sigma)
	General electrostatic Poisson solver routine.

virtual void	advanceReactionNetwork (const Real a_time, const Real a_dt)
	Advance the reaction network. This will compute the electric field (on the appropriate phase) and call the other version.

virtual void	advanceReactionNetwork (Vector< EBAMRCellData * > &a_cdrSources, Vector< EBAMRCellData * > &a_rteSources, const Vector< EBAMRCellData * > &a_cdrDensities, const Vector< EBAMRCellData * > &a_rteDensities, const EBAMRCellData &a_E, const Real &a_time, const Real &a_dt)
	Compute reaction network source terms.

virtual void	advanceReactionNetwork (Vector< EBAMRCellData * > &a_cdrSources, Vector< EBAMRCellData * > &a_rteSources, const Vector< EBAMRCellData * > &a_cdrDensities, const Vector< EBAMRCellData * > &a_cdrGradients, const Vector< EBAMRCellData * > &a_rteDensities, const EBAMRCellData &a_E, const Real &a_time, const Real &a_dt)
	Compute reaction network source terms. This is the AMR version that does all levels.

virtual void	advanceReactionNetwork (Vector< LevelData< EBCellFAB > * > &a_cdrSources, Vector< LevelData< EBCellFAB > * > &a_rteSources, const Vector< LevelData< EBCellFAB > * > &a_cdrDensities, const Vector< LevelData< EBCellFAB > * > &a_cdrGradients, const Vector< LevelData< EBCellFAB > * > &a_rteDensities, const LevelData< EBCellFAB > &a_E, const Real &a_time, const Real &a_dt, const int a_level)
	Compute reaction network source terms. This is the level version.

virtual void	advanceReactionNetworkRegularCells (Vector< FArrayBox * > &a_cdrSources, Vector< FArrayBox * > &a_rteSources, const Vector< FArrayBox * > &a_cdrDensities, const Vector< FArrayBox * > &a_cdrGradients, const Vector< FArrayBox * > &a_rteDensities, const FArrayBox &a_E, const Real &a_time, const Real &a_dt, const Real &a_dx, const Box &a_cellBox)
	Compute reaction network sources. This is the regular version which does not do the cut-cells.

virtual void	advanceReactionNetworkIrreg (Vector< EBCellFAB * > &a_cdrSources, Vector< EBCellFAB * > &a_rteSources, const Vector< EBCellFAB * > &a_cdrDensities, const Vector< EBCellFAB * > &a_cdrGradients, const Vector< EBCellFAB * > &a_cdrVelocities, const Vector< EBCellFAB * > &a_rteDensities, const EBCellFAB &a_E, const Real &a_time, const Real &a_dt, const Real &a_dx, const Box &a_cellBox, const int a_lvl, const DataIndex &a_dit)
	Compute reaction network sources. This is the irregular version which only does the cut-cells.

virtual void	advanceReactionNetworkIrregInterp (Vector< EBCellFAB * > &a_cdrSources, Vector< EBCellFAB * > &a_rteSources, const Vector< EBCellFAB * > &a_cdrDensities, const Vector< EBCellFAB * > &a_cdrGradients, const Vector< EBCellFAB * > &a_cdrVelocities, const Vector< EBCellFAB * > &a_rteDensities, const EBCellFAB &a_E, const Real &a_time, const Real &a_dt, const Real &a_dx, const Box &a_box, const int a_lvl, const DataIndex &a_dit)
	Compute reaction network sources. This version interpolates everything to the centroid and runs with that data.

virtual void	advanceReactionNetworkIrregUpwind (Vector< EBCellFAB * > &a_cdrSources, Vector< EBCellFAB * > &a_rteSources, const Vector< EBCellFAB * > &a_cdrDensities, const Vector< EBCellFAB * > &a_cdrGradients, const Vector< EBCellFAB * > &a_cdrVelocities, const Vector< EBCellFAB * > &a_rteDensities, const EBCellFAB &a_E, const Real &a_time, const Real &a_dt, const Real &a_dx, const Box &a_box, const int a_lvl, const DataIndex &a_dit)
	Compute reaction network sources.

virtual void	computeCdrDiffusion ()
	Compute diffusion things using whatever is available in the solvers.

virtual void	computeCdrDiffusion (const EBAMRCellData &a_electricFieldCell, const EBAMRIVData &a_electricFieldEB)
	Compute diffusion things using whatever is available in the solvers. This version uses the input electric fields.

virtual void	computeCdrDiffusionCell (Vector< EBAMRCellData > &a_cdrDcoCell, const Vector< EBAMRCellData * > &a_cdrDensities, const EBAMRCellData &a_electricFieldCell, const Real &a_time)
	Compute CDR diffusion coefficients on cell centers. This will call the level version.

virtual void	computeCdrDiffusionCell (Vector< LevelData< EBCellFAB > * > &a_cdrDcoCell, const Vector< LevelData< EBCellFAB > * > &a_cdrDensities, const LevelData< EBCellFAB > &a_electricFieldCell, const int a_lvl, const Real &a_time)
	Compute diffusion coefficients. This is the level version and it will iterate through grid patches.

virtual void	computeCdrDiffusionCellRegular (Vector< FArrayBox * > &a_cdrDcoCell, const Vector< FArrayBox * > &a_cdrDensities, const FArrayBox &a_electricFieldCell, const Box a_cellBox, const Real a_dx, const Real a_time)
	Compute diffusion coefficients on cell centers. This does all the regular cells.

virtual void	computeCdrDiffusionCellIrregular (Vector< EBCellFAB * > &a_cdrDcoCell, const Vector< EBCellFAB * > &a_cdrDensities, const EBCellFAB &a_electricFieldCell, const Real a_dx, const Real &a_time, const int a_lvl, const DataIndex &a_dit)
	Compute diffusion coefficients on cell centers. This version does the irregular cells.

virtual void	computeCdrDiffusionFace (Vector< EBAMRFluxData * > &a_cdrDcoFace, const Vector< EBAMRCellData * > &a_cdrDensities, const EBAMRCellData &a_electricFieldCell, const Real &a_time)
	Compute diffusion coefficients on face centers.

virtual void	computeCdrDiffusionEb (Vector< EBAMRIVData * > &a_cdrDcoEB, const Vector< EBAMRIVData * > &a_cdrDensitiesEB, const EBAMRIVData &a_electricFieldEB, const Real &a_time)
	Compute diffusion coefficients on the EB. This is the AMR version – it will call the level version.

virtual void	computeCdrDiffusionEb (Vector< LevelData< BaseIVFAB< Real > > * > &a_cdrDcoEB, const Vector< LevelData< BaseIVFAB< Real > > * > &a_cdrDensitiesEB, const LevelData< BaseIVFAB< Real > > &a_electricFieldEB, const Real &a_time, const int a_lvl)
	Compute diffusion coefficients on the EB. This is a level version but it will go through all patches.

virtual void	computeCdrDriftVelocities ()
	Compute the CDR drift velocities.

virtual void	computeCdrDriftVelocities (Vector< EBAMRCellData * > &a_velocities, const Vector< EBAMRCellData * > &a_cdrDensities, const EBAMRCellData &a_E, const Real &a_time)
	Compute the cell-centered CDR velocities. This will call the level version.

virtual void	computeCdrDriftVelocities (Vector< LevelData< EBCellFAB > * > &a_cdrVelocities, const Vector< LevelData< EBCellFAB > * > &a_cdrDensities, const LevelData< EBCellFAB > &a_electricField, const int a_lvl, const Real &a_time)
	Compute the cell-centered CDR velocities.

virtual void	computeCdrDriftVelocitiesRegular (Vector< FArrayBox * > &a_cdrVelocities, const Vector< FArrayBox * > &a_cdrDensities, const FArrayBox &a_electricField, const Box &a_cellBox, const Real &a_time, const Real &a_dx)
	Compute the cell-centered CDR velocities.

virtual void	computeCdrDriftVelocitiesIrregular (Vector< EBCellFAB * > &a_velocities, const Vector< EBCellFAB * > &a_cdrDensities, const EBCellFAB &a_electricField, const Real &a_time, const Real &a_dx, const int a_lvl, const DataIndex &a_dit)
	Compute the cell-centered CDR velocities.

virtual void	computeCdrFluxes (Vector< EBAMRIVData * > &a_cdrFluxesEB, const Vector< EBAMRIVData * > &a_extrapCdrFluxes, const Vector< EBAMRIVData * > &a_extrapCdrDensities, const Vector< EBAMRIVData * > &a_extrapCdrVelocities, const Vector< EBAMRIVData * > &a_extrapCdrGradients, const Vector< EBAMRIVData * > &a_extrapRteFluxes, const EBAMRIVData &a_electricField, const Real &a_time)
	Compute CDR fluxes on the EB. This is the AMR version – it will call the level version.

virtual void	computeCdrFluxes (Vector< LevelData< BaseIVFAB< Real > > * > &a_cdrFluxesEB, const Vector< LevelData< BaseIVFAB< Real > > * > &a_extrapCdrFluxes, const Vector< LevelData< BaseIVFAB< Real > > * > &a_extrapCdrDensities, const Vector< LevelData< BaseIVFAB< Real > > * > &a_extrapCdrVelocities, const Vector< LevelData< BaseIVFAB< Real > > * > &a_extrapCdrGradients, const Vector< LevelData< BaseIVFAB< Real > > * > &a_extrapRteFluxes, const LevelData< BaseIVFAB< Real > > &a_electricField, const Real &a_time, const int a_lvl)
	Compute CDR fluxes on the EB. This is the level version which will fill EB fluxes on both electrode and dielectric cells.

virtual void	computeCdrDomainFluxes (Vector< EBAMRIFData * > &a_cdrFluxes, const Vector< EBAMRIFData * > &a_extrapCdrFluxes, const Vector< EBAMRIFData * > &a_extrapCdrDensities, const Vector< EBAMRIFData * > &a_extrapCdrVelocities, const Vector< EBAMRIFData * > &a_extrapCdrGradients, const Vector< EBAMRIFData * > &a_extrapRteFluxes, const EBAMRIFData &a_field, const Real &a_time)
	Compute CDR fluxes on domain faces for usage in boundary conditions. This is the AMR version.

virtual void	computeCdrDomainFluxes (Vector< LevelData< DomainFluxIFFAB > * > a_cdrFluxes, const Vector< LevelData< DomainFluxIFFAB > * > &a_extrapCdrFluxes, const Vector< LevelData< DomainFluxIFFAB > * > &a_extrapCdrDensities, const Vector< LevelData< DomainFluxIFFAB > * > &a_extrapCdrVelocities, const Vector< LevelData< DomainFluxIFFAB > * > &a_extrapCdrGradients, const Vector< LevelData< DomainFluxIFFAB > * > &a_extrapRteFluxes, const LevelData< DomainFluxIFFAB > &a_electricField, const Real &a_time, const int a_lvl)
	Compute CDR fluxes on domain faces for usage in boundary conditions. This is the level version.

virtual void	computeExtrapolatedFluxes (Vector< EBAMRIVData * > &a_extrapCdrFluxesEB, const Vector< EBAMRCellData * > a_cdrDensities, const Vector< EBAMRCellData * > a_cdrVelocities, const phase::which_phase a_phase)
	Extrapolate cell-centered fluxes to the EB centroid.

virtual void	computeExtrapolatedDomainFluxes (Vector< EBAMRIFData * > &a_cdrDomainFluxes, const Vector< EBAMRCellData * > a_cdrDensities, const Vector< EBAMRCellData * > a_cdrVelocities, const phase::which_phase a_phase)
	Extrapolate drift fluxes to domain faces.

virtual void	computeExtrapolatedVelocities (Vector< EBAMRIVData * > &a_cdrVelocitiesEB, const Vector< EBAMRCellData * > a_cdrVelocitiesCell, const phase::which_phase a_phase)
	Extrapolated cell-centered velocities to the EB.

virtual void	deallocateSolverInternals ()
	Deallocate internal solver storages.

virtual void	extrapolateToEb (Vector< EBAMRIVData * > &a_ebData, const phase::which_phase a_phase, const Vector< EBAMRCellData * > &a_cellData)
	Extrapolated cell-centered data to the EB. This is the AMR version for all species.

virtual void	extrapolateToEb (EBAMRIVData &a_ebData, const phase::which_phase a_phase, const EBAMRCellData &a_cellData)
	Extrapolated cell-centered data to the EB. This is the AMR version.

virtual void	extrapolateToEb (LevelData< BaseIVFAB< Real > > &a_ebData, const phase::which_phase a_phase, const LevelData< EBCellFAB > &a_cellData, const int a_lvl)
	Extrapolated cell-centered data to the EB. This is the level version.

virtual void	extrapolateToDomainFaces (Vector< EBAMRIFData * > &a_domainData, const phase::which_phase a_phase, const Vector< EBAMRCellData * > &a_cellData)
	Extrapolate cell-centered data to domain faces. This is the AMR version that does all species.

virtual void	extrapolateToDomainFaces (EBAMRIFData &a_domainData, const phase::which_phase a_phase, const EBAMRCellData &a_cellData)
	Extrapolate cell-centered data to domain faces. This is the AMR version.

virtual void	extrapolateToDomainFaces (LevelData< DomainFluxIFFAB > &a_domainData, const phase::which_phase a_phase, const LevelData< EBCellFAB > &a_cellData, const int a_lvl)
	Extrapolate cell-centered data to domain faces. This is the AMR version.

virtual void	extrapolateVectorToDomainFaces (Vector< EBAMRIFData * > &a_domainData, const phase::which_phase a_phase, const Vector< EBAMRCellData * > &a_cellData)
	Extrapolate vector data to domain faces. This is the AMR version that does all species.

virtual void	extrapolateVectorToDomainFaces (EBAMRIFData &a_domainData, const phase::which_phase a_phase, const EBAMRCellData &a_cellData)
	Extrapolate vector data to domain faces. This is the AMR version.

virtual void	extrapolateVelocitiesToDomainFaces (Vector< EBAMRIFData * > &a_domainVelocities, const phase::which_phase a_phase, const Vector< EBAMRCellData * > &a_cellVelocities)
	Extrapolate velocities data to domain faces. This is the AMR version that does all species.

virtual void	getCdrMax (Real &a_cdrMax, std::string &a_solverName) const
	Get maximum density in the CDR species. This will fetch both the value and the solver name.

virtual void	initialSigma ()
	Initialize the surface charge.

virtual void	projectFlux (EBAMRIVData &a_projectedFlux, const EBAMRIVData &a_flux)
	Project a flux onto the EB normal. This is the AMR version.

virtual void	projectFlux (LevelData< BaseIVFAB< Real > > &a_projectedFlux, const LevelData< BaseIVFAB< Real > > &a_flux, const int a_lvl)
	Project a flux onto the EB normal. This is the level version.

virtual void	projectDomain (EBAMRIFData &a_projectedFlux, const EBAMRIFData &a_flux)
	Project flux onto domain edges/faces.

virtual void	regridInternals (const int a_lmin, const int a_oldFinestLevel, const int a_newFinestLevel)=0
	Regrid internal storage. This MUST be overwritten by implementations.

virtual void	regridSolvers (const int a_lmin, const int a_oldFinestLevel, const int a_newFinestLevel)
	Regrid all of our solvers.

virtual void	resetDielectricCells (EBAMRIVData &a_data) const
	Set data defined over dielectric cells to zero.

virtual void	sanityCheck () const
	Do a sanity check. Only works when DEBUG=TRUE.

virtual void	setupCdr ()
	Setup the CDR solvers.

virtual void	setupPoisson ()
	Setup the field solver.

virtual void	setupRadiativeTransfer ()
	Setup the RTE solvers.

virtual void	setupSigma ()
	Set up the surface charge solver.

virtual void	setSolverVerbosity ()
	Set solver verbosity.

virtual void	solverDump ()
	All solvers writes plot files. These are not the same as the plot files written by Driver.

virtual void	setVoltage (std::function< Real(const Real a_time)> a_voltage)
	Set voltage curve.

virtual void	solveRadiativeTransfer (const Real a_dt)
	Solve the radiative transfer problem.

virtual void	setCdrSolvers (RefCountedPtr< CdrLayout< CdrSolver > > &a_cdr)
	Set the CDR layout.

virtual void	setFieldSolver (RefCountedPtr< FieldSolver > &a_fieldSolver)
	Set the field solver.

virtual void	setRadiativeTransferSolvers (RefCountedPtr< RtLayout< RtSolver > > &a_rte)
	Set the RTE solvers.

virtual Real	computeElectrodeCurrent ()
	Compute the ohmic current on electrodes.

virtual Real	computeDielectricCurrent ()
	Compute the ohmic current on dielectrics.

virtual Real	computeDomainCurrent ()
	Compute the ohmic current through the domain faces.

virtual Real	computeOhmicInductionCurrent ()
	Compute induced current in external circuit due to ohmic conduction.

virtual Real	computeRelaxationTime ()
	Compute the relaxation time as dt = eps0/conductivity.

virtual Real	getDt ()
	Get dt.

Public Member Functions inherited from TimeStepper
	TimeStepper ()
	Default constructor (does nothing)

virtual	~TimeStepper ()
	Default destructor (does nothing)

void	setAmr (const RefCountedPtr< AmrMesh > &a_amr)
	Set AmrMesh.

void	setComputationalGeometry (const RefCountedPtr< ComputationalGeometry > &a_computationalGeometry)
	Set the computational geometry.

virtual void	postPlot ()
	An option for calling special functions prior to plotting data. Called by Driver in the IMMEDIATELY after writing the plot file.

virtual Vector< long int >	getCheckpointLoads (const std::string a_realm, const int a_level) const
	Get computational loads to be checkpointed.

virtual bool	needToRegrid ()
	Fuction which can have Driver do regrids at arbitrary points in the simulation.

bool	keepGoing () const
	Query whether the time stepper wants to continue stepping.

virtual bool	loadBalanceThisRealm (const std::string a_realm) const
	Load balancing query for a specified realm. If this returns true for a_realm, load balancing routines will be called during regrids.

virtual void	loadBalanceBoxes (Vector< Vector< int > > &a_procs, Vector< Vector< Box > > &a_boxes, const std::string a_realm, const Vector< DisjointBoxLayout > &a_grids, const int a_lmin, const int a_finestLevel)
	Load balance grid boxes for a specified realm.

Protected Types
enum class	SourceTermComputation { Interpolated , Upwind }
	Enum for switching between source term computations. 'Interpolated' interpolates the cell-centered data to cell centroids. InterpolatedStable is the same as Interpolated but with a stability fix near inflow EBs. CellAverage assumes centroid-centered data, and Upwind uses the upwind formulation in Villa et. al.

enum class	TimeCode { Advection , Diffusion , AdvectionDiffusion , Source , RelaxationTime , Restricted , Hardcap , Error , Physics }
	An enum for encapsulating how time steps were restricted.

Protected Member Functions
virtual void	parseVerbosity ()
	Parse class verbosity.

virtual void	parseSolverVerbosity ()
	Parse solver verbosities.

virtual void	parseCFL ()
	Parse the CFL number.

virtual void	parseRelaxationTime ()
	Parse the relaxation time restriction.

virtual void	parseMinDt ()
	Parse the minimum allowed time step.

virtual void	parseMaxDt ()
	Parse the maximum allowed time step.

virtual void	parseFastPoisson ()
	Parse the "fast poisson" method, i.e. how often we solve the Poisson equation.

virtual void	parseFastRadiativeTransfer ()
	Parse the "fast rte" method, i.e. how often we solve the Poisson equation.

virtual void	parseSourceComputation ()
	Parse how we compute the source terms.

virtual void	writeJ (LevelData< EBCellFAB > &a_output, int &a_comp, const std::string a_outputRealm, const int a_level) const
	Compute and put the current density in output data holder.

virtual void	computePhysicsPlotVars (EBAMRCellData &a_physicsPlotVars) const noexcept
	Calculate the user-supplied plot variables.

virtual void	writeData (LevelData< EBCellFAB > &a_output, int &a_comp, const EBAMRCellData &a_data, const std::string a_outputRealm, const int a_level, const bool a_interpToCentroids, const bool a_interpGhost) const noexcept
	Write data to output. Convenience function.

Protected Attributes
std::string	m_realm
	Realm where the time stepper is registered.

std::string	m_className
	Time stepper class name.

phase::which_phase	m_phase
	Plasma phase.

SourceTermComputation	m_whichSourceTermComputation
	Which source term computation.

RefCountedPtr< MultiFluidIndexSpace >	m_multifluidIndexSpace
	Index space.

RefCountedPtr< CdrPlasmaPhysics >	m_physics
	Plasma kinetics.

RefCountedPtr< CdrLayout< CdrSolver > >	m_cdr
	CDR solvers.

RefCountedPtr< RtLayout< RtSolver > >	m_rte
	Radiative transfer solvers.

RefCountedPtr< FieldSolver >	m_fieldSolver
	Poisson solver.

RefCountedPtr< SurfaceODESolver< 1 > >	m_sigma
	Surface charge solver.

EBAMRCellData	m_currentDensity
	Storage for the current.

EBAMRCellData	m_physicsPlotVars
	Storage for physics-observables.

std::function< Real(const Real a_time)>	m_voltage
	Set the potential.

Real	m_minDt
	Minimum allowed time step.

Real	m_maxDt
	Maximum allowed time step.

Real	m_cfl
	CFL number.

Real	m_relaxTime
	Fudge factor for relaxation time.

Real	m_time
	TIme.

Real	m_dt
	Previous time step size.

Real	m_dtCFL
	Computed CFL time step.

TimeCode	m_timeCode
	Time code for step restriction.

int	m_solverVerbosity
	Verbosity for solvers.

int	m_fastRTE
	Fast RTE solver.

int	m_fastPoisson
	Fast RTE solver.

int	m_upwindFactor
	Upwind factor.

Protected Attributes inherited from TimeStepper
int	m_verbosity
	Class verbosity.

int	m_timeStep
	Time step.

Real	m_time
	TIme.

Real	m_dt
	Previous time step size.

bool	m_keepGoing
	If false, Driver will stop the time loop after the current step.

RefCountedPtr< AmrMesh >	m_amr
	AmrMesh.

RefCountedPtr< ComputationalGeometry >	m_computationalGeometry
	Computational geometry.

Detailed Description

Abstract class for evolving the minimal plasma model in time.

This class implements all the required coupling functions for more easily constructing algorithms. Subclasses must implement the advance method. This class handles all I/O.

Constructor & Destructor Documentation

◆ CdrPlasmaStepper()

CdrPlasmaStepper::CdrPlasmaStepper ( RefCountedPtr< CdrPlasmaPhysics > & a_physics )

Full constructor.

Parameters

[in] a_physics CDR plasma physics.

Member Function Documentation

◆ advance()

virtual Real Physics::CdrPlasma::CdrPlasmaStepper::advance ( const Real a_dt )

overridepure virtual

Advance method, advances equations.

Parameters

[in] a_dt Time step

Returns: Actual time step used (can be different from a_dt for adaptive methods)

Implements TimeStepper.

Implemented in Physics::CdrPlasma::CdrPlasmaGodunovStepper, and Physics::CdrPlasma::CdrPlasmaImExSdcStepper.

◆ advanceReactionNetwork() [1/4]

void CdrPlasmaStepper::advanceReactionNetwork	(	const Real	a_time,
		const Real	a_dt
	)

virtual

Advance the reaction network. This will compute the electric field (on the appropriate phase) and call the other version.

Parameters

[in]	a_time	Time
[in]	a_dt	Time step

◆ advanceReactionNetwork() [2/4]

void CdrPlasmaStepper::advanceReactionNetwork	(	Vector< EBAMRCellData * > &	a_cdrSources,
		Vector< EBAMRCellData * > &	a_rteSources,
		const Vector< EBAMRCellData * > &	a_cdrDensities,
		const Vector< EBAMRCellData * > &	a_cdrGradients,
		const Vector< EBAMRCellData * > &	a_rteDensities,
		const EBAMRCellData &	a_E,
		const Real &	a_time,
		const Real &	a_dt
	)

virtual

Compute reaction network source terms. This is the AMR version that does all levels.

Parameters

[out]	a_cdrSources	Source terms for the CDR solvers.
[out]	a_rteSources	Source terms for the radiative transfer solvers.
[in]	a_cdrDensities	Cell-centered densities for the plasma species (i.e., CDR solvers).
[in]	a_cdrGradients	Cell-centered gradients for the plasma species (i.e., CDR solvers).
[in]	a_rteDensities	Cell-centered densities for the radiative transfer equations (i.e., RTE solvers).
[in]	a_E	Cell-centered electric field.
[in]	a_time	Time
[in]	a_dt	Time step.

◆ advanceReactionNetwork() [3/4]

void CdrPlasmaStepper::advanceReactionNetwork	(	Vector< EBAMRCellData * > &	a_cdrSources,
		Vector< EBAMRCellData * > &	a_rteSources,
		const Vector< EBAMRCellData * > &	a_cdrDensities,
		const Vector< EBAMRCellData * > &	a_rteDensities,
		const EBAMRCellData &	a_E,
		const Real &	a_time,
		const Real &	a_dt
	)

virtual

Compute reaction network source terms.

This will call the other version but will first compute the gradients of the plasma densities.

Parameters

[out]	a_cdrSources	Source terms for the CDR solvers.
[out]	a_rteSources	Source terms for the radiative transfer solvers.
[in]	a_cdrDensities	Cell-centered densities for the plasma species (i.e., CDR solvers).
[in]	a_rteDensities	Cell-centered densities for the radiative transfer equations (i.e., RTE solvers).
[in]	a_E	Cell-centered electric field.
[in]	a_time	Time
[in]	a_dt	Time step.

◆ advanceReactionNetwork() [4/4]

void CdrPlasmaStepper::advanceReactionNetwork	(	Vector< LevelData< EBCellFAB > * > &	a_cdrSources,
		Vector< LevelData< EBCellFAB > * > &	a_rteSources,
		const Vector< LevelData< EBCellFAB > * > &	a_cdrDensities,
		const Vector< LevelData< EBCellFAB > * > &	a_cdrGradients,
		const Vector< LevelData< EBCellFAB > * > &	a_rteDensities,
		const LevelData< EBCellFAB > &	a_E,
		const Real &	a_time,
		const Real &	a_dt,
		const int	a_level
	)

virtual

Compute reaction network source terms. This is the level version.

Parameters

[out]	a_cdrSources	Source terms for the CDR solvers.
[out]	a_rteSources	Source terms for the radiative transfer solvers.
[in]	a_cdrDensities	Cell-centered densities for the plasma species (i.e., CDR solvers).
[in]	a_cdrGradients	Cell-centered gradients for the plasma species (i.e., CDR solvers).
[in]	a_rteDensities	Cell-centered densities for the radiative transfer equations (i.e., RTE solvers).
[in]	a_E	Cell-centered electric field.
[in]	a_time	Time
[in]	a_dt	Time step.
[in]	a_level	Grid level.

◆ advanceReactionNetworkIrreg()

void CdrPlasmaStepper::advanceReactionNetworkIrreg	(	Vector< EBCellFAB * > &	a_cdrSources,
		Vector< EBCellFAB * > &	a_rteSources,
		const Vector< EBCellFAB * > &	a_cdrDensities,
		const Vector< EBCellFAB * > &	a_cdrGradients,
		const Vector< EBCellFAB * > &	a_cdrVelocities,
		const Vector< EBCellFAB * > &	a_rteDensities,
		const EBCellFAB &	a_E,
		const Real &	a_time,
		const Real &	a_dt,
		const Real &	a_dx,
		const Box &	a_cellBox,
		const int	a_lvl,
		const DataIndex &	a_dit
	)

virtual

Compute reaction network sources. This is the irregular version which only does the cut-cells.

Parameters

[out]	a_cdrSources	Source terms for the CDR solvers.
[out]	a_rteSources	Source terms for the radiative transfer solvers.
[in]	a_cdrDensities	Cell-centered densities for the plasma species (i.e., CDR solvers).
[in]	a_cdrGradients	Cell-centered gradients for the plasma species (i.e., CDR solvers).
[in]	a_cdrVelocities	Cell-centered velocities for the CDR solvers.
[in]	a_rteDensities	Cell-centered densities for the radiative transfer equations (i.e., RTE solvers).
[in]	a_E	Cell-centered electric field.
[in]	a_time	Time
[in]	a_dt	Time step.
[in]	a_dx	Grid resolution.
[in]	a_cellBox	Cell-centered computational region.
[in]	a_lvl	Grid level
[in]	a_dit	Grid patch index

Note: This will call one of the other versions.

◆ advanceReactionNetworkIrregInterp()

void CdrPlasmaStepper::advanceReactionNetworkIrregInterp	(	Vector< EBCellFAB * > &	a_cdrSources,
		Vector< EBCellFAB * > &	a_rteSources,
		const Vector< EBCellFAB * > &	a_cdrDensities,
		const Vector< EBCellFAB * > &	a_cdrGradients,
		const Vector< EBCellFAB * > &	a_cdrVelocities,
		const Vector< EBCellFAB * > &	a_rteDensities,
		const EBCellFAB &	a_E,
		const Real &	a_time,
		const Real &	a_dt,
		const Real &	a_dx,
		const Box &	a_box,
		const int	a_lvl,
		const DataIndex &	a_dit
	)

virtual

Compute reaction network sources. This version interpolates everything to the centroid and runs with that data.

Parameters

[out]	a_cdrSources	Source terms for the CDR solvers.
[out]	a_rteSources	Source terms for the radiative transfer solvers.
[in]	a_cdrDensities	Cell-centered densities for the plasma species (i.e., CDR solvers).
[in]	a_cdrGradients	Cell-centered gradients for the plasma species (i.e., CDR solvers).
[in]	a_cdrVelocities	Cell-centered velocities for the CDR solvers.
[in]	a_rteDensities	Cell-centered densities for the radiative transfer equations (i.e., RTE solvers).
[in]	a_E	Cell-centered electric field.
[in]	a_time	Time
[in]	a_dt	Time step.
[in]	a_dx	Grid resolution.
[in]	a_cellBox	Cell-centered computational region.
[in]	a_lvl	Grid level
[in]	a_dit	Grid patch index

◆ advanceReactionNetworkIrregUpwind()

void CdrPlasmaStepper::advanceReactionNetworkIrregUpwind	(	Vector< EBCellFAB * > &	a_cdrSources,
		Vector< EBCellFAB * > &	a_rteSources,
		const Vector< EBCellFAB * > &	a_cdrDensities,
		const Vector< EBCellFAB * > &	a_cdrGradients,
		const Vector< EBCellFAB * > &	a_cdrVelocities,
		const Vector< EBCellFAB * > &	a_rteDensities,
		const EBCellFAB &	a_E,
		const Real &	a_time,
		const Real &	a_dt,
		const Real &	a_dx,
		const Box &	a_box,
		const int	a_lvl,
		const DataIndex &	a_dit
	)

virtual

Compute reaction network sources.

This version is called with the Villa et. al upwind correction. It can also be used with regularly interpolated data, in which case it will turn off reactions in EB cells where the mass flows into the domain. This function assumes that all input is on the cell center.

Parameters

[out]	a_cdrSources	Source terms for the CDR solvers.
[out]	a_rteSources	Source terms for the radiative transfer solvers.
[in]	a_cdrDensities	Cell-centered densities for the plasma species (i.e., CDR solvers).
[in]	a_cdrGradients	Cell-centered gradients for the plasma species (i.e., CDR solvers).
[in]	a_cdrVelocities	Cell-centered velocities for the CDR solvers.
[in]	a_rteDensities	Cell-centered densities for the radiative transfer equations (i.e., RTE solvers).
[in]	a_E	Cell-centered electric field.
[in]	a_time	Time
[in]	a_dt	Time step.
[in]	a_dx	Grid resolution.
[in]	a_cellBox	Cell-centered computational region.
[in]	a_lvl	Grid level
[in]	a_dit	Grid patch index

◆ advanceReactionNetworkRegularCells()

void CdrPlasmaStepper::advanceReactionNetworkRegularCells	(	Vector< FArrayBox * > &	a_cdrSources,
		Vector< FArrayBox * > &	a_rteSources,
		const Vector< FArrayBox * > &	a_cdrDensities,
		const Vector< FArrayBox * > &	a_cdrGradients,
		const Vector< FArrayBox * > &	a_rteDensities,
		const FArrayBox &	a_E,
		const Real &	a_time,
		const Real &	a_dt,
		const Real &	a_dx,
		const Box &	a_cellBox
	)

virtual

Compute reaction network sources. This is the regular version which does not do the cut-cells.

Parameters

[out]	a_cdrSources	Source terms for the CDR solvers.
[out]	a_rteSources	Source terms for the radiative transfer solvers.
[in]	a_cdrDensities	Cell-centered densities for the plasma species (i.e., CDR solvers).
[in]	a_cdrGradients	Cell-centered gradients for the plasma species (i.e., CDR solvers).
[in]	a_rteDensities	Cell-centered densities for the radiative transfer equations (i.e., RTE solvers).
[in]	a_E	Cell-centered electric field.
[in]	a_time	Time
[in]	a_dt	Time step.
[in]	a_dx	Grid resolution.
[in]	a_cellBox	Cell-centered computational region.

◆ allocate()

void CdrPlasmaStepper::allocate ( )

overridevirtual

Allocate data for the time stepper and solvers.

Implements TimeStepper.

◆ allocateInternals()

void CdrPlasmaStepper::allocateInternals ( )

pure virtual

Allocate internal storage. Must be implemented by subclasses.

Implemented in Physics::CdrPlasma::CdrPlasmaGodunovStepper, and Physics::CdrPlasma::CdrPlasmaImExSdcStepper.

◆ computeCdrDiffusion() [1/2]

void CdrPlasmaStepper::computeCdrDiffusion ( )

virtual

Compute diffusion things using whatever is available in the solvers.

This will compute the electric field on the cell center and EB, and then call the other version.

◆ computeCdrDiffusion() [2/2]

void CdrPlasmaStepper::computeCdrDiffusion	(	const EBAMRCellData &	a_electricFieldCell,
		const EBAMRIVData &	a_electricFieldEB
	)

virtual

Compute diffusion things using whatever is available in the solvers. This version uses the input electric fields.

Parameters

[in]	a_electricFieldCell	Electric field on cell centers.
[in]	a_electricFieldEB	Electric field on EB centroids

◆ computeCdrDiffusionCell() [1/2]

void CdrPlasmaStepper::computeCdrDiffusionCell	(	Vector< EBAMRCellData > &	a_cdrDcoCell,
		const Vector< EBAMRCellData * > &	a_cdrDensities,
		const EBAMRCellData &	a_electricFieldCell,
		const Real &	a_time
	)

virtual

Compute CDR diffusion coefficients on cell centers. This will call the level version.

Parameters

[out]	a_cdrDcoCell	CDR diffusion coefficients on cell centers.
[in]	a_cdrDensities	Cell-centered CDR densities
[in]	a_electricFieldCell	Cell-centered electric field
[in]	a_time	Time

◆ computeCdrDiffusionCell() [2/2]

void CdrPlasmaStepper::computeCdrDiffusionCell	(	Vector< LevelData< EBCellFAB > * > &	a_cdrDcoCell,
		const Vector< LevelData< EBCellFAB > * > &	a_cdrDensities,
		const LevelData< EBCellFAB > &	a_electricFieldCell,
		const int	a_lvl,
		const Real &	a_time
	)

virtual

Compute diffusion coefficients. This is the level version and it will iterate through grid patches.

Parameters

[out]	a_cdrDcoCell	CDR diffusion coefficients on cell centers.
[in]	a_cdrDensities	Cell-centered CDR densities
[in]	a_electricFieldCell	Cell-centered electric field
[in]	a_lvl	Grid level
[in]	a_time	Time

◆ computeCdrDiffusionCellIrregular()

void CdrPlasmaStepper::computeCdrDiffusionCellIrregular	(	Vector< EBCellFAB * > &	a_cdrDcoCell,
		const Vector< EBCellFAB * > &	a_cdrDensities,
		const EBCellFAB &	a_electricFieldCell,
		const Real	a_dx,
		const Real &	a_time,
		const int	a_lvl,
		const DataIndex &	a_dit
	)

virtual

Compute diffusion coefficients on cell centers. This version does the irregular cells.

Parameters

[out]	a_cdrDcoCell	CDR diffusion coefficients on cell centers.
[in]	a_cdrDensities	Cell-centered CDR densities
[in]	a_electricFieldCell	Cell-centered electric field
[in]	a_dx	Grid resolution
[in]	a_time	Time
[in]	a_lvl	Grid level
[in]	a_dit	Grid patch index

◆ computeCdrDiffusionCellRegular()

void CdrPlasmaStepper::computeCdrDiffusionCellRegular	(	Vector< FArrayBox * > &	a_cdrDcoCell,
		const Vector< FArrayBox * > &	a_cdrDensities,
		const FArrayBox &	a_electricFieldCell,
		const Box	a_cellBox,
		const Real	a_dx,
		const Real	a_time
	)

virtual

Compute diffusion coefficients on cell centers. This does all the regular cells.

Parameters

[out]	a_cdrDcoCell	CDR diffusion coefficients on cell centers.
[in]	a_cdrDensities	Cell-centered CDR densities
[in]	a_electricFieldCell	Cell-centered electric field
[in]	a_cellBox	Kernel region
[in]	a_dx	Grid resolution
[in]	a_time	Time

◆ computeCdrDiffusionEb() [1/2]

void CdrPlasmaStepper::computeCdrDiffusionEb	(	Vector< EBAMRIVData * > &	a_cdrDcoEB,
		const Vector< EBAMRIVData * > &	a_cdrDensitiesEB,
		const EBAMRIVData &	a_electricFieldEB,
		const Real &	a_time
	)

virtual

Compute diffusion coefficients on the EB. This is the AMR version – it will call the level version.

Parameters

[out]	a_cdrDcoEB	EB-centered_CDR diffusion coefficients
[in]	a_cdrDensitiesEB	EB-centered CDR densities
[in]	a_electricFieldEB	EB-centered electric field
[in]	a_time	Time

◆ computeCdrDiffusionEb() [2/2]

void CdrPlasmaStepper::computeCdrDiffusionEb	(	Vector< LevelData< BaseIVFAB< Real > > * > &	a_cdrDcoEB,
		const Vector< LevelData< BaseIVFAB< Real > > * > &	a_cdrDensitiesEB,
		const LevelData< BaseIVFAB< Real > > &	a_electricFieldEB,
		const Real &	a_time,
		const int	a_lvl
	)

virtual

Compute diffusion coefficients on the EB. This is a level version but it will go through all patches.

Parameters

[out]	a_cdrDcoEB	EB-centered_CDR diffusion coefficients
[in]	a_cdrDensitiesEB	EB-centered CDR densities
[in]	a_electricFieldEB	EB-centered electric field
[in]	a_time	Time
[in]	a_lvl	Grid level.

◆ computeCdrDiffusionFace()

void CdrPlasmaStepper::computeCdrDiffusionFace	(	Vector< EBAMRFluxData * > &	a_cdrDcoFace,
		const Vector< EBAMRCellData * > &	a_cdrDensities,
		const EBAMRCellData &	a_electricFieldCell,
		const Real &	a_time
	)

virtual

Compute diffusion coefficients on face centers.

This will first compute the diffusion coefficients on cell centers, and then interpolate the result to face centers.

Parameters

[out]	a_cdrDcoCell	CDR diffusion coefficients on cell centers.
[in]	a_cdrDensities	Cell-centered CDR densities
[in]	a_electricFieldCell	Cell-centered electric field
[in]	a_time	Time

◆ computeCdrDomainFluxes() [1/2]

void CdrPlasmaStepper::computeCdrDomainFluxes	(	Vector< EBAMRIFData * > &	a_cdrFluxes,
		const Vector< EBAMRIFData * > &	a_extrapCdrFluxes,
		const Vector< EBAMRIFData * > &	a_extrapCdrDensities,
		const Vector< EBAMRIFData * > &	a_extrapCdrVelocities,
		const Vector< EBAMRIFData * > &	a_extrapCdrGradients,
		const Vector< EBAMRIFData * > &	a_extrapRteFluxes,
		const EBAMRIFData &	a_field,
		const Real &	a_time
	)

virtual

Compute CDR fluxes on domain faces for usage in boundary conditions. This is the AMR version.

Parameters

[out]	a_cdrFluxes	CDR domain bounadry fluxes for solvers.
[in]	a_extrapCdrFluxes	CDR fluxes extrapolated to the domain faces.
[in]	a_extrapCdrDensities	CDR densities extrapolated to the domain faces.
[in]	a_extrapCdrVelocities	CDR velocities extrapolated to the domain faces.
[in]	a_extrapCdrGradients	CDR gradients extrapolated to the domain faces.
[in]	a_extrapRteFluxes	RTE fluxes on the domain faces.
[in]	a_electricField	Electric field on the domain faces.
[in]	a_time	Time.

◆ computeCdrDomainFluxes() [2/2]

void CdrPlasmaStepper::computeCdrDomainFluxes	(	Vector< LevelData< DomainFluxIFFAB > * >	a_cdrFluxes,
		const Vector< LevelData< DomainFluxIFFAB > * > &	a_extrapCdrFluxes,
		const Vector< LevelData< DomainFluxIFFAB > * > &	a_extrapCdrDensities,
		const Vector< LevelData< DomainFluxIFFAB > * > &	a_extrapCdrVelocities,
		const Vector< LevelData< DomainFluxIFFAB > * > &	a_extrapCdrGradients,
		const Vector< LevelData< DomainFluxIFFAB > * > &	a_extrapRteFluxes,
		const LevelData< DomainFluxIFFAB > &	a_electricField,
		const Real &	a_time,
		const int	a_lvl
	)

virtual

Compute CDR fluxes on domain faces for usage in boundary conditions. This is the level version.

This will iterate over the patches on each level and compute the domain fluxes.

Parameters

[out]	a_cdrFluxes	CDR domain bounadry fluxes for solvers.
[in]	a_extrapCdrFluxes	CDR fluxes extrapolated to the domain faces.
[in]	a_extrapCdrDensities	CDR densities extrapolated to the domain faces.
[in]	a_extrapCdrVelocities	CDR velocities extrapolated to the domain faces.
[in]	a_extrapCdrGradients	CDR gradients extrapolated to the domain faces.
[in]	a_extrapRteFluxes	RTE fluxes on the domain faces.
[in]	a_electricField	Electric field on the domain faces.
[in]	a_time	Time.
[in]	a_lvl	Grid level

◆ computeCdrDriftVelocities() [1/3]

void CdrPlasmaStepper::computeCdrDriftVelocities ( )

virtual

Compute the CDR drift velocities.

This will compute the cell-centered electric field and then call the other version.

◆ computeCdrDriftVelocities() [2/3]

void CdrPlasmaStepper::computeCdrDriftVelocities	(	Vector< EBAMRCellData * > &	a_velocities,
		const Vector< EBAMRCellData * > &	a_cdrDensities,
		const EBAMRCellData &	a_E,
		const Real &	a_time
	)

virtual

Compute the cell-centered CDR velocities. This will call the level version.

Parameters

[out]	a_cdrVelocities	Cell-centered CDR drift velocities
[out]	a_cdrDensities	Cell-centered CDR densities
[in]	a_electricField	The cell-centered electric field
[in]	a_time	Time

◆ computeCdrDriftVelocities() [3/3]

void CdrPlasmaStepper::computeCdrDriftVelocities	(	Vector< LevelData< EBCellFAB > * > &	a_cdrVelocities,
		const Vector< LevelData< EBCellFAB > * > &	a_cdrDensities,
		const LevelData< EBCellFAB > &	a_electricField,
		const int	a_lvl,
		const Real &	a_time
	)

virtual

Compute the cell-centered CDR velocities.

This is the level version – it will call the patch version for each patch in the grid.

Parameters

[out]	a_cdrVelocities	Cell-centered CDR drift velocities
[out]	a_cdrDensities	Cell-centered CDR densities
[in]	a_electricField	The cell-centered electric field
[in]	a_lvl	Grid level
[in]	a_time	Time

◆ computeCdrDriftVelocitiesIrregular()

void CdrPlasmaStepper::computeCdrDriftVelocitiesIrregular	(	Vector< EBCellFAB * > &	a_velocities,
		const Vector< EBCellFAB * > &	a_cdrDensities,
		const EBCellFAB &	a_electricField,
		const Real &	a_time,
		const Real &	a_dx,
		const int	a_lvl,
		const DataIndex &	a_dit
	)

virtual

Compute the cell-centered CDR velocities.

This is patch version for irregular grid cells.

Parameters

[out]	a_cdrVelocities	Cell-centered CDR drift velocities
[out]	a_cdrDensities	Cell-centered CDR densities
[in]	a_electricField	The cell-centered electric field
[in]	a_time	Time
[in]	a_dx	Grid resolution
[in]	a_lvl	Grid level
[in]	a_dit	Grid index

◆ computeCdrDriftVelocitiesRegular()

void CdrPlasmaStepper::computeCdrDriftVelocitiesRegular	(	Vector< FArrayBox * > &	a_cdrVelocities,
		const Vector< FArrayBox * > &	a_cdrDensities,
		const FArrayBox &	a_electricField,
		const Box &	a_cellBox,
		const Real &	a_time,
		const Real &	a_dx
	)

virtual

Compute the cell-centered CDR velocities.

This is patch version for regular grid cells.

Parameters

[out]	a_cdrVelocities	Cell-centered CDR drift velocities
[out]	a_cdrDensities	Cell-centered CDR densities
[in]	a_electricField	The cell-centered electric field
[in]	a_cellBox	Computational region.
[in]	a_time	Time
[in]	a_dx	Grid resolution

◆ computeCdrFluxes() [1/2]

void CdrPlasmaStepper::computeCdrFluxes	(	Vector< EBAMRIVData * > &	a_cdrFluxesEB,
		const Vector< EBAMRIVData * > &	a_extrapCdrFluxes,
		const Vector< EBAMRIVData * > &	a_extrapCdrDensities,
		const Vector< EBAMRIVData * > &	a_extrapCdrVelocities,
		const Vector< EBAMRIVData * > &	a_extrapCdrGradients,
		const Vector< EBAMRIVData * > &	a_extrapRteFluxes,
		const EBAMRIVData &	a_electricField,
		const Real &	a_time
	)

virtual

Compute CDR fluxes on the EB. This is the AMR version – it will call the level version.

Parameters

[out]	a_cdrFluxesEB	CDR fluxes to be put in CDR solvers as a BC on the EBs.
[in]	a_extrapCdrFluxes	CDR fluxes extrapolated from the inside. Can be used to construct inflow/outflow BCs.
[in]	a_extrapCdrDensities	CDR densities extrapolated from the inside.
[in]	a_extrapCdrVelocities	CDR velocities on the EB, extrapolated from the inside.
[in]	a_extrapCdrGradients	CDR gradients on the EB, extrapolated from the inside.
[in]	a_electricField	Electric field on the EB.
[in]	a_time	Time.

◆ computeCdrFluxes() [2/2]

void CdrPlasmaStepper::computeCdrFluxes	(	Vector< LevelData< BaseIVFAB< Real > > * > &	a_cdrFluxesEB,
		const Vector< LevelData< BaseIVFAB< Real > > * > &	a_extrapCdrFluxes,
		const Vector< LevelData< BaseIVFAB< Real > > * > &	a_extrapCdrDensities,
		const Vector< LevelData< BaseIVFAB< Real > > * > &	a_extrapCdrVelocities,
		const Vector< LevelData< BaseIVFAB< Real > > * > &	a_extrapCdrGradients,
		const Vector< LevelData< BaseIVFAB< Real > > * > &	a_extrapRteFluxes,
		const LevelData< BaseIVFAB< Real > > &	a_electricField,
		const Real &	a_time,
		const int	a_lvl
	)

virtual

Compute CDR fluxes on the EB. This is the level version which will fill EB fluxes on both electrode and dielectric cells.

Parameters

[out]	a_cdrFluxesEB	CDR fluxes to be put in CDR solvers as a BC on the EBs.
[in]	a_extrapCdrFluxes	CDR fluxes extrapolated from the inside. Can be used to construct inflow/outflow BCs.
[in]	a_extrapCdrDensities	CDR densities extrapolated from the inside.
[in]	a_extrapCdrVelocities	CDR velocities on the EB, extrapolated from the inside.
[in]	a_extrapCdrGradients	CDR gradients on the EB, extrapolated from the inside.
[in]	a_electricField	Electric field on the EB.
[in]	a_time	Time.
[in]	a_lvl	Grid level.

◆ computeCellConductivity()

void CdrPlasmaStepper::computeCellConductivity ( EBAMRCellData & a_cellConductivity ) const

virtual

Compute the cell-centered conductivity.

Parameters

[out] a_cellConductivity Cell-centered conductivity.

◆ computeDielectricCurrent()

Real CdrPlasmaStepper::computeDielectricCurrent ( )

virtual

Compute the ohmic current on dielectrics.

This will compute the total charge flux through dielectric surfaces and return it.

◆ computeDomainCurrent()

Real CdrPlasmaStepper::computeDomainCurrent ( )

virtual

Compute the ohmic current through the domain faces.

This will compute the total current through domain faces and return it.

◆ computeDt()

virtual Real Physics::CdrPlasma::CdrPlasmaStepper::computeDt ( )

overridepure virtual

Compute the time step.

Subclasses must implement this one.

Implements TimeStepper.

Implemented in Physics::CdrPlasma::CdrPlasmaGodunovStepper, and Physics::CdrPlasma::CdrPlasmaImExSdcStepper.

◆ computeElectricField() [1/5]

void CdrPlasmaStepper::computeElectricField	(	EBAMRCellData &	a_E,
		const phase::which_phase	a_phase,
		const MFAMRCellData &	a_potential
	)		const

virtual

Compute the cell-centered electric field on a specific phase.

Parameters

[out]	a_electricFieldCell	Cell-centered electric field
[in]	a_phase	Phase
[in]	a_potential	Electrostatic potential (cell-centered)

◆ computeElectricField() [2/5]

void CdrPlasmaStepper::computeElectricField	(	EBAMRCellData &	a_electricFieldCell,
		const phase::which_phase	a_phase
	)		const

virtual

Compute the cell-centered electric field on a specific phase using whatever is available in m_fieldSolver.

This will use the potential from m_fieldSolver.

Parameters

[out]	a_electricFieldCell	Cell-centered electric field
[in]	a_phase	Phase

◆ computeElectricField() [3/5]

void CdrPlasmaStepper::computeElectricField	(	EBAMRFluxData &	a_electricFieldFace,
		const phase::which_phase	a_phase,
		const EBAMRCellData &	a_electricFieldCell
	)		const

virtual

Compute the face-centered electric field. This will only do the normal component.

Parameters

[out]	a_electricFieldFace	Face-centered electric field
[in]	a_phase	Phase
[in]	a_electricFieldCell	Cell-centered electric field (all components).

◆ computeElectricField() [4/5]

void CdrPlasmaStepper::computeElectricField	(	EBAMRIVData &	a_electricFieldEB,
		const phase::which_phase	a_phase,
		const EBAMRCellData &	a_electricFieldCell
	)		const

virtual

Compute the EB-centered electric field.

Parameters

[out]	a_electricFieldEB	EB-centered electric field.
[in]	a_phase	Phase
[in]	a_electricFieldCell	Cell-centered electric field (all components).

◆ computeElectricField() [5/5]

void CdrPlasmaStepper::computeElectricField	(	MFAMRCellData &	a_electricFieldCell,
		const MFAMRCellData &	a_potential
	)		const

virtual

Compute the cell-centered electric field on both phases.

Parameters

[out]	a_electricFieldCell	Cell-centered electric field
[in]	a_potential	Cell-centered electrostatic potential

◆ computeElectrodeCurrent()

Real CdrPlasmaStepper::computeElectrodeCurrent ( )

virtual

Compute the ohmic current on electrodes.

This will compute the total charge flux through electrode surfaces and return it.

Note: In 2D this will yield the current per unit length into the domain.

◆ computeExtrapolatedDomainFluxes()

void CdrPlasmaStepper::computeExtrapolatedDomainFluxes	(	Vector< EBAMRIFData * > &	a_cdrDomainFluxes,
		const Vector< EBAMRCellData * >	a_cdrDensities,
		const Vector< EBAMRCellData * >	a_cdrVelocities,
		const phase::which_phase	a_phase
	)

virtual

Extrapolate drift fluxes to domain faces.

Note that this projects the flux on the domain face vectors.

Parameters

[out]	a_cdrDomainFluxes	CDR fluxes extrapolated and projected onto the domain face.
[in]	a_cdrDensities	Cell-centered CDR solver densities.
[in]	a_cdrVelocities	Cell-centered CDR solver velocities.
[in]	a_phase	Phase

◆ computeExtrapolatedFluxes()

void CdrPlasmaStepper::computeExtrapolatedFluxes	(	Vector< EBAMRIVData * > &	a_extrapCdrFluxesEB,
		const Vector< EBAMRCellData * >	a_cdrDensities,
		const Vector< EBAMRCellData * >	a_cdrVelocities,
		const phase::which_phase	a_phase
	)

virtual

Extrapolate cell-centered fluxes to the EB centroid.

Parameters

[out]	a_extrapCdrFluxesEB	CDR fluxes extrapolated to the EB.
[in]	a_cdrDensities	CDR densities on the cell center
[in]	a_cdrVelocities	CDR velocities on the cell center
[in]	a_phase	Phase

◆ computeExtrapolatedVelocities()

void CdrPlasmaStepper::computeExtrapolatedVelocities	(	Vector< EBAMRIVData * > &	a_cdrVelocitiesEB,
		const Vector< EBAMRCellData * >	a_cdrVelocitiesCell,
		const phase::which_phase	a_phase
	)

virtual

Extrapolated cell-centered velocities to the EB.

This only computes the normal velocity on the EB, i.e. the output is the extrapolated velocity projected along the EB normal.

Parameters

[out]	a_cdrVelocitiesEB	EB-centered velocities
[in]	a_cdrVelocitiesCell	Cell-centered velocities
[in]	a_phase	Phase

◆ computeFaceConductivity()

void CdrPlasmaStepper::computeFaceConductivity	(	EBAMRFluxData &	a_conductivityFace,
		EBAMRIVData &	a_conductivityEB,
		const EBAMRCellData &	a_conductivityCell
	)		const

virtual

Compute the face-centered conductivity from a cell-centered conductivity.

Parameters

[out]	a_conductivityFace	Face-centered conductivity.
[out]	a_conductivityEB	EB-centered conductivity.
[in]	a_conductivityCell	Cell-centered conductivity.

◆ computeJ()

void CdrPlasmaStepper::computeJ ( EBAMRCellData & a_J ) const

virtual

Compute the current density.

Parameters

[out] a_J Current density

◆ computeMaxElectricField()

void CdrPlasmaStepper::computeMaxElectricField	(	Real &	a_maxElectricField,
		const phase::which_phase	a_phase
	)

virtual

Compute the maximum of the electric field.

Parameters

[out]	a_maxElectricField	Maximum electric field magnitude evaluated over the input phase
[in]	a_phase	Input phase

◆ computeOhmicInductionCurrent()

Real CdrPlasmaStepper::computeOhmicInductionCurrent ( )

virtual

Compute induced current in external circuit due to ohmic conduction.

This will return int(E*J dV) over the entire domain.

◆ computeSpaceChargeDensity() [1/2]

void CdrPlasmaStepper::computeSpaceChargeDensity ( )

virtual

Compute the centroid-centered space charge density by using the data inside the CDR solvers.

This will call the other version with the storage holders being the space charge density data holder in m_fieldSolver and the cell-centered densities in the CDR solvers.

◆ computeSpaceChargeDensity() [2/2]

void CdrPlasmaStepper::computeSpaceChargeDensity	(	MFAMRCellData &	a_rho,
		const Vector< EBAMRCellData * > &	a_cdrDensities
	)

virtual

Compute the centroid-centered space charge density.

This will first compute the cell-centered space charge density and then interpolate that to the centroid.

Parameters

[out]	a_rho	Centroid centered space charge density.
[in]	a_cdrDensities	Cell-centered plasma species densities. Must be in the same order as the densities in m_cdr.

◆ deallocateInternals()

void CdrPlasmaStepper::deallocateInternals ( )

pure virtual

Deallocate internals. Must be implemented by subclasses.

Do nothing

Implemented in Physics::CdrPlasma::CdrPlasmaGodunovStepper, and Physics::CdrPlasma::CdrPlasmaImExSdcStepper.

◆ extrapolateToDomainFaces() [1/3]

void CdrPlasmaStepper::extrapolateToDomainFaces	(	EBAMRIFData &	a_domainData,
		const phase::which_phase	a_phase,
		const EBAMRCellData &	a_cellData
	)

virtual

Extrapolate cell-centered data to domain faces. This is the AMR version.

Parameters

[out]	a_domainData	Domain-centered data
[in]	a_phase	Phase
[in]	a_cellData	Cell-centered data.

◆ extrapolateToDomainFaces() [2/3]

void CdrPlasmaStepper::extrapolateToDomainFaces	(	LevelData< DomainFluxIFFAB > &	a_domainData,
		const phase::which_phase	a_phase,
		const LevelData< EBCellFAB > &	a_cellData,
		const int	a_lvl
	)

virtual

Extrapolate cell-centered data to domain faces. This is the AMR version.

Parameters

[out]	a_domainData	Domain-centered data
[in]	a_phase	Phase
[in]	a_cellData	Cell-centered data.
[in]	a_lvl	Grid level

◆ extrapolateToDomainFaces() [3/3]

void CdrPlasmaStepper::extrapolateToDomainFaces	(	Vector< EBAMRIFData * > &	a_domainData,
		const phase::which_phase	a_phase,
		const Vector< EBAMRCellData * > &	a_cellData
	)

virtual

Extrapolate cell-centered data to domain faces. This is the AMR version that does all species.

Parameters

[out]	a_domainData	Domain-centered data
[in]	a_phase	Phase
[in]	a_cellData	Cell-centered data.

◆ extrapolateToEb() [1/3]

void CdrPlasmaStepper::extrapolateToEb	(	EBAMRIVData &	a_ebData,
		const phase::which_phase	a_phase,
		const EBAMRCellData &	a_cellData
	)

virtual

Extrapolated cell-centered data to the EB. This is the AMR version.

Parameters

[out]	a_ebData	EB-centered data
[in]	a_phase	Phase
[in]	a_cellData	Cell-centered data

◆ extrapolateToEb() [2/3]

void CdrPlasmaStepper::extrapolateToEb	(	LevelData< BaseIVFAB< Real > > &	a_ebData,
		const phase::which_phase	a_phase,
		const LevelData< EBCellFAB > &	a_cellData,
		const int	a_lvl
	)

virtual

Extrapolated cell-centered data to the EB. This is the level version.

Parameters

[out]	a_ebData	EB-centered data
[in]	a_phase	Phase
[in]	a_cellData	Cell-centered data
[in]	a_lvl	Grid level

◆ extrapolateToEb() [3/3]

void CdrPlasmaStepper::extrapolateToEb	(	Vector< EBAMRIVData * > &	a_ebData,
		const phase::which_phase	a_phase,
		const Vector< EBAMRCellData * > &	a_cellData
	)

virtual

Extrapolated cell-centered data to the EB. This is the AMR version for all species.

Parameters

[out]	a_ebData	EB-centered data
[in]	a_phase	Phase
[in]	a_cellData	Cell-centered data

◆ extrapolateVectorToDomainFaces() [1/2]

void CdrPlasmaStepper::extrapolateVectorToDomainFaces	(	EBAMRIFData &	a_domainData,
		const phase::which_phase	a_phase,
		const EBAMRCellData &	a_cellData
	)

virtual

Extrapolate vector data to domain faces. This is the AMR version.

This will extrapolate the vector data to domain faces and project it along the domain normal.

Parameters

[out]	a_domainData	Domain-centered data. Must have one component.
[in]	a_phase	Phase
[in]	a_cellData	Cell-centered data. Must have SpaceDim components.

◆ extrapolateVectorToDomainFaces() [2/2]

void CdrPlasmaStepper::extrapolateVectorToDomainFaces	(	Vector< EBAMRIFData * > &	a_domainData,
		const phase::which_phase	a_phase,
		const Vector< EBAMRCellData * > &	a_cellData
	)

virtual

Extrapolate vector data to domain faces. This is the AMR version that does all species.

This will extrapolate the vector data to domain faces and project it along the domain normal.

Parameters

[out]	a_domainData	Domain-centered data. Must have one component.
[in]	a_phase	Phase
[in]	a_cellData	Cell-centered data. Must have SpaceDim components.

◆ extrapolateVelocitiesToDomainFaces()

void CdrPlasmaStepper::extrapolateVelocitiesToDomainFaces	(	Vector< EBAMRIFData * > &	a_domainVelocities,
		const phase::which_phase	a_phase,
		const Vector< EBAMRCellData * > &	a_cellVelocities
	)

virtual

Extrapolate velocities data to domain faces. This is the AMR version that does all species.

This will extrapolate the velocities data to domain faces.

Parameters

[out]	a_domainVelocities	Domain-centered velocities. Must have SpaceDim components.
[in]	a_phase	Phase
[in]	a_cellData	Cell-centered velocities. Must have SpaceDim components.

◆ getCdrMax()

void CdrPlasmaStepper::getCdrMax	(	Real &	a_cdrMax,
		std::string &	a_solverName
	)		const

virtual

Get maximum density in the CDR species. This will fetch both the value and the solver name.

Parameters

[out]	a_cdrMax	Maximum density found in one of the CDR solvers
[out]	a_solverName	The solver name

◆ getNumberOfPlotVariables()

int CdrPlasmaStepper::getNumberOfPlotVariables ( ) const

overridevirtual

Get the number of plot variables for this time stepper.

This is necessary because Driver, not TimeStepper, is responsible for allocating the necessary memory.

Returns: Returns number of plot variables that will be written during writePlotData

Implements TimeStepper.

◆ getPlotVariableNames()

Vector< std::string > CdrPlasmaStepper::getPlotVariableNames ( ) const

overridevirtual

Get plot variable names.

Implements TimeStepper.

◆ getTime()

Real CdrPlasmaStepper::getTime ( ) const

virtual

Return time.

Returns: m_time.

◆ initialData()

void CdrPlasmaStepper::initialData ( )

overridevirtual

Fill all solvers with initial data.

Implements TimeStepper.

◆ parseOptions()

virtual void Physics::CdrPlasma::CdrPlasmaStepper::parseOptions ( )

pure virtual

Parse TimeStepper class options. Must be implemented by subclasses.

Implemented in Physics::CdrPlasma::CdrPlasmaGodunovStepper, and Physics::CdrPlasma::CdrPlasmaImExSdcStepper.

◆ parseRuntimeOptions()

virtual void Physics::CdrPlasma::CdrPlasmaStepper::parseRuntimeOptions ( )

overridepure virtual

Parse runtime options. Subclasses must implement this.

Reimplemented from TimeStepper.

Implemented in Physics::CdrPlasma::CdrPlasmaGodunovStepper, and Physics::CdrPlasma::CdrPlasmaImExSdcStepper.

◆ parseSourceComputation()

void CdrPlasmaStepper::parseSourceComputation ( )

protectedvirtual

Parse how we compute the source terms.

This is for modifications near the EB.

◆ postCheckpointSetup()

void CdrPlasmaStepper::postCheckpointSetup ( )

overridevirtual

Post-initialize operations to be performed after filling solvers with data read from checkpoint files.

Implements TimeStepper.

Reimplemented in Physics::CdrPlasma::CdrPlasmaGodunovStepper.

◆ postInitialize()

void CdrPlasmaStepper::postInitialize ( )

overridevirtual

Post-initialize operations to be performed after filling solvers with initialData.

Implements TimeStepper.

◆ postRegrid()

void CdrPlasmaStepper::postRegrid ( )

overridevirtual

Perform post-regrid operations.

This includes all operations to be done AFTER interpolating data to new grids.

Implements TimeStepper.

Reimplemented in Physics::CdrPlasma::CdrPlasmaGodunovStepper.

◆ prePlot()

void CdrPlasmaStepper::prePlot ( )

overridevirtual

Compute some thing that go into plot files.

Reimplemented from TimeStepper.

◆ preRegrid()

void CdrPlasmaStepper::preRegrid	(	const int	a_lmin,
		const int	a_oldFinestLevel
	)

overridevirtual

Pre-regrid function – this is the one called by Driver.

Parameters

[in]	a_lmin	The coarsest level that changes
[in]	a_oldFinestLevel	The finest level before the regrid.

Implements TimeStepper.

Reimplemented in Physics::CdrPlasma::CdrPlasmaGodunovStepper.

◆ preRegridInternals()

void CdrPlasmaStepper::preRegridInternals	(	const int	a_lmin,
		const int	a_oldFinestLevel
	)

virtual

Do a preRegrid operation for internal storage (if needed).

Parameters

[in]	a_lmin	The coarsest level that changes
[in]	a_oldFinestLevel	The finest level before the regrid.

◆ printStepReport()

void CdrPlasmaStepper::printStepReport ( )

overridevirtual

Print a step report.

This is called by Driver after time step. The routine can be used to display use information about the simulation progress.

Implements TimeStepper.

◆ projectDomain()

void CdrPlasmaStepper::projectDomain	(	EBAMRIFData &	a_projectedFlux,
		const EBAMRIFData &	a_flux
	)

virtual

Project flux onto domain edges/faces.

Parameters

[out]	a_projectedFlux	Projected flux – must have one component.
[in]	a_flux	Vector flux – must have SpaceDim components.

◆ projectFlux() [1/2]

void CdrPlasmaStepper::projectFlux	(	EBAMRIVData &	a_projectedFlux,
		const EBAMRIVData &	a_flux
	)

virtual

Project a flux onto the EB normal. This is the AMR version.

Parameters

[out]	a_projectedFlux	Projected flux – must have one component.
[in]	a_flux	Vector flux – must have SpaceDim components.

◆ projectFlux() [2/2]

void CdrPlasmaStepper::projectFlux	(	LevelData< BaseIVFAB< Real > > &	a_projectedFlux,
		const LevelData< BaseIVFAB< Real > > &	a_flux,
		const int	a_lvl
	)

virtual

Project a flux onto the EB normal. This is the level version.

Parameters

[out]	a_projectedFlux	Projected flux – must have one component.
[in]	a_flux	Vector flux – must have SpaceDim components.

◆ registerOperators()

void CdrPlasmaStepper::registerOperators ( )

overridevirtual

Register operators to be used for the simulation.

Implements TimeStepper.

◆ registerRealms()

void CdrPlasmaStepper::registerRealms ( )

overridevirtual

Register realms to be used for the simulation.

Implements TimeStepper.

◆ regrid()

void CdrPlasmaStepper::regrid	(	const int	a_lmin,
		const int	a_oldFinestLevel,
		const int	a_newFinestLevel
	)

overridevirtual

Regrid method.

Parameters

[in]	a_lmin	The coarsest level that changed.
[in]	a_oldFinestLevel	The finest level before the regrid.
[in]	a_newFinestLevel	The finest level after the regrid.

Implements TimeStepper.

Reimplemented in Physics::CdrPlasma::CdrPlasmaGodunovStepper.

◆ regridInternals()

virtual void Physics::CdrPlasma::CdrPlasmaStepper::regridInternals	(	const int	a_lmin,
		const int	a_oldFinestLevel,
		const int	a_newFinestLevel
	)

pure virtual

Regrid internal storage. This MUST be overwritten by implementations.

Parameters

[in]	a_lmin	The coarsest level that changed.
[in]	a_oldFinestLevel	The finest level before the regrid.
[in]	a_newFinestLevel	The finest level after the regrid.

Implemented in Physics::CdrPlasma::CdrPlasmaGodunovStepper, and Physics::CdrPlasma::CdrPlasmaImExSdcStepper.

◆ regridSolvers()

void CdrPlasmaStepper::regridSolvers	(	const int	a_lmin,
		const int	a_oldFinestLevel,
		const int	a_newFinestLevel
	)

virtual

Regrid all of our solvers.

Parameters

[in]	a_lmin	The coarsest level that changed.
[in]	a_oldFinestLevel	The finest level before the regrid.
[in]	a_newFinestLevel	The finest level after the regrid.

◆ resetDielectricCells()

void CdrPlasmaStepper::resetDielectricCells ( EBAMRIVData & a_data ) const

virtual

Set data defined over dielectric cells to zero.

Parameters

[in] a_data EB-centered data.

◆ setCdrSolvers()

void CdrPlasmaStepper::setCdrSolvers ( RefCountedPtr< CdrLayout< CdrSolver > > & a_cdr )

virtual

Set the CDR layout.

Parameters

[in] a_cdr CDR solvers.

◆ setFieldSolver()

void CdrPlasmaStepper::setFieldSolver ( RefCountedPtr< FieldSolver > & a_fieldSolver )

virtual

Set the field solver.

Parameters

[in] a_fieldSolver Field solver

◆ setRadiativeTransferSolvers()

void CdrPlasmaStepper::setRadiativeTransferSolvers ( RefCountedPtr< RtLayout< RtSolver > > & a_rte )

virtual

Set the RTE solvers.

Parameters

[in] a_rte RTE solvers

◆ setSolverVerbosity()

void CdrPlasmaStepper::setSolverVerbosity ( )

virtual

Set solver verbosity.

This will set the solver verbosity to m_solverVerbosity.

◆ setupSemiImplicitPoisson() [1/2]

void CdrPlasmaStepper::setupSemiImplicitPoisson	(	const EBAMRFluxData &	a_conductivityFace,
		const EBAMRIVData &	a_conductivityEB,
		const Real	a_factor
	)

virtual

Set up a semi-implicit Poisson solver.

This is the version that uses pre-existing conductivities in a_conductivity. The factor dictates the multiplication factor f in the expression div((epsilonr + f*sigma) grad(phi)) where sigma is the conductivity and f is a scaling factor.

Parameters

[in]	a_conductivityFace	Face-centered conductivities (the Poisson solver will interpolate them to face centroids in the multigrid solve).
[in]	a_conductivityEB	EB-centered conductivities.
[in]	a_factor	Factor to multiply conductivities with (they could be pre-scaled by dt/eps0 in which case a_factor = 1).

◆ setupSemiImplicitPoisson() [2/2]

void CdrPlasmaStepper::setupSemiImplicitPoisson ( const Real a_dt )

virtual

Set up a semi-implicit Poisson solver.

This will compute the required conductivities first, using data already present in the CDR solvers. After that, it will set up the Poisson equation with spatially varying conductivities.

Parameters

[in] a_dt Time step.

◆ setupSolvers()

void CdrPlasmaStepper::setupSolvers ( )

overridevirtual

Instantiate solvers.

Implements TimeStepper.

◆ setVoltage()

void CdrPlasmaStepper::setVoltage ( std::function< Real(const Real a_time)> a_voltage )

virtual

Set voltage curve.

Parameters

[in] a_voltage Voltage curve

◆ solvePoisson() [1/2]

bool CdrPlasmaStepper::solvePoisson ( )

virtual

Solve the electrostatic Poisson equation.

This will first compute the centroid-centered space charge density (using the cell-centered data in the CDR solvers) and then ask the field solver to solve.

◆ solvePoisson() [2/2]

bool CdrPlasmaStepper::solvePoisson	(	MFAMRCellData &	a_potential,
		MFAMRCellData &	a_rho,
		const Vector< EBAMRCellData * >	a_cdrDensities,
		const EBAMRIVData &	a_sigma
	)

virtual

General electrostatic Poisson solver routine.

This will compute the space charge density onto a_rho and then solve for the potential.

Parameters

[in,out]	a_potential	Electrostatic potential
[out]	a_rho	Centroid-centered space charge density.
[in]	a_cdrDensities	Plasma species densities.
[in]	a_sigma	Surface charge density on dielectric EBs.

◆ solveRadiativeTransfer()

void CdrPlasmaStepper::solveRadiativeTransfer ( const Real a_dt )

virtual

Solve the radiative transfer problem.

This will call the other version with states taken from the solvers.

Parameters

[in] a_dt Time step.

◆ stationaryRTE()

bool CdrPlasmaStepper::stationaryRTE ( )

virtual

Check if RTE solvers are stationary.

Returns: True if the RTE solvers are stationary.

◆ synchronizeSolverTimes()

void CdrPlasmaStepper::synchronizeSolverTimes	(	const int	a_step,
		const Real	a_time,
		const Real	a_dt
	)

overridevirtual

Synchronize solver times.

Parameters

[in]	a_step	Time step
[in]	a_time	Time (in seconds)
[in]	a_dt	Time step that was used.

Implements TimeStepper.

◆ writeData()

void CdrPlasmaStepper::writeData	(	LevelData< EBCellFAB > &	a_output,
		int &	a_comp,
		const EBAMRCellData &	a_data,
		const std::string	a_outputRealm,
		const int	a_level,
		const bool	a_interpToCentroids,
		const bool	a_interpGhost
	)		const

protectedvirtualnoexcept

Write data to output. Convenience function.

Parameters

[in,out]	a_output	Output data holder.
[in,out]	a_icomp	Starting component where this solver begins writing the output.
[in]	a_data	Data to write.
[in]	a_outputRealm	Realm where a_output lives
[in]	a_level	Grid level
[in]	a_interpToCentroids	If true, a_data will be interpolated to cell centroids before writing to a_output.
[in]	a_interpGhost	If true, interpolate ghost cells

◆ writeJ()

void CdrPlasmaStepper::writeJ	(	LevelData< EBCellFAB > &	a_output,
		int &	a_comp,
		const std::string	a_outputRealm,
		const int	a_level
	)		const

protectedvirtual

Compute and put the current density in output data holder.

Parameters

[in,out]	a_output	Output data holder.
[in,out]	a_icomp	Starting component in a_output to begin at. On output this is incremented by SpaceDim.
[in]	a_outputRealm	Realm where a_output lives
[in]	a_level	Grid level

◆ writePlotData()

void CdrPlasmaStepper::writePlotData	(	LevelData< EBCellFAB > &	a_output,
		int &	a_icomp,
		const std::string	a_outputRealm,
		const int	a_level
	)		const

overridevirtual

Write plot data to output holder.

Parameters

[in,out]	a_output	Output data holder.
[in,out]	a_icomp	Starting component in a_output to begin at.
[in]	a_outputRealm	Realm where a_output belongs
[in]	a_level	Grid level

Implements TimeStepper.

Member Data Documentation

◆ m_realm

std::string Physics::CdrPlasma::CdrPlasmaStepper::m_realm

protected

Realm where the time stepper is registered.

This will always be the primal realm because CdrPlasma has no dual-grid functionality.

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

Physics/CdrPlasma/CD_CdrPlasmaStepper.H
Physics/CdrPlasma/CD_CdrPlasmaStepper.cpp

Public Member Functions

Protected Types

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ CdrPlasmaStepper()

Member Function Documentation

◆ advance()

◆ advanceReactionNetwork() [1/4]

◆ advanceReactionNetwork() [2/4]

◆ advanceReactionNetwork() [3/4]

◆ advanceReactionNetwork() [4/4]

◆ advanceReactionNetworkIrreg()

◆ advanceReactionNetworkIrregInterp()

◆ advanceReactionNetworkIrregUpwind()

◆ advanceReactionNetworkRegularCells()

◆ allocate()

◆ allocateInternals()

◆ computeCdrDiffusion() [1/2]

◆ computeCdrDiffusion() [2/2]

◆ computeCdrDiffusionCell() [1/2]

◆ computeCdrDiffusionCell() [2/2]

◆ computeCdrDiffusionCellIrregular()

◆ computeCdrDiffusionCellRegular()

◆ computeCdrDiffusionEb() [1/2]

◆ computeCdrDiffusionEb() [2/2]

◆ computeCdrDiffusionFace()

◆ computeCdrDomainFluxes() [1/2]

◆ computeCdrDomainFluxes() [2/2]

◆ computeCdrDriftVelocities() [1/3]

◆ computeCdrDriftVelocities() [2/3]

◆ computeCdrDriftVelocities() [3/3]

◆ computeCdrDriftVelocitiesIrregular()

◆ computeCdrDriftVelocitiesRegular()

◆ computeCdrFluxes() [1/2]

◆ computeCdrFluxes() [2/2]

◆ computeCellConductivity()

◆ computeDielectricCurrent()

◆ computeDomainCurrent()

◆ computeDt()

◆ computeElectricField() [1/5]

◆ computeElectricField() [2/5]

◆ computeElectricField() [3/5]

◆ computeElectricField() [4/5]

◆ computeElectricField() [5/5]

◆ computeElectrodeCurrent()

◆ computeExtrapolatedDomainFluxes()

◆ computeExtrapolatedFluxes()

◆ computeExtrapolatedVelocities()

◆ computeFaceConductivity()

◆ computeJ()

◆ computeMaxElectricField()

◆ computeOhmicInductionCurrent()

◆ computeSpaceChargeDensity() [1/2]

◆ computeSpaceChargeDensity() [2/2]

◆ deallocateInternals()

◆ extrapolateToDomainFaces() [1/3]

◆ extrapolateToDomainFaces() [2/3]

◆ extrapolateToDomainFaces() [3/3]

◆ extrapolateToEb() [1/3]

◆ extrapolateToEb() [2/3]

◆ extrapolateToEb() [3/3]

◆ extrapolateVectorToDomainFaces() [1/2]

◆ extrapolateVectorToDomainFaces() [2/2]

◆ extrapolateVelocitiesToDomainFaces()

◆ getCdrMax()

◆ getNumberOfPlotVariables()

◆ getPlotVariableNames()

◆ getTime()

◆ initialData()

◆ parseOptions()

◆ parseRuntimeOptions()

◆ parseSourceComputation()

◆ postCheckpointSetup()

◆ postInitialize()

◆ postRegrid()

◆ prePlot()

◆ preRegrid()

◆ preRegridInternals()