chombo-discharge
Loading...
Searching...
No Matches
CD_CdrPlasmaStepper.H
Go to the documentation of this file.
1/*
2 * SPDX-FileCopyrightText: 2021-2026 SINTEF Energy Research
3 *
4 * SPDX-License-Identifier: GPL-3.0-or-later
5 */
6
13#ifndef CD_CDRPLASMASTEPPER_H
14#define CD_CDRPLASMASTEPPER_H
15
16// Std includes
17#include <functional>
18
19// Our includes
20#include <CD_ComputationalGeometry.H>
21#include <CD_CdrPlasmaPhysics.H>
22#include <CD_AmrMesh.H>
23#include <CD_MultiFluidIndexSpace.H>
24#include <CD_CdrLayout.H>
25#include <CD_RtLayout.H>
26#include <CD_FieldSolver.H>
27#include <CD_SurfaceODESolver.H>
28#include <CD_TimeStepper.H>
29#include <CD_NamespaceHeader.H>
30
31namespace Physics {
32namespace CdrPlasma {
33
45class CdrPlasmaStepper : public TimeStepper
46{
47public:
51 CdrPlasmaStepper();
52
57 CdrPlasmaStepper(RefCountedPtr<CdrPlasmaPhysics>& a_physics);
58
62 ~CdrPlasmaStepper();
63
67 virtual void
68 setupSolvers() override;
69
73 virtual void
74 allocate() override;
75
79 virtual void
80 initialData() override;
81
85 virtual void
86 postInitialize() override;
87
91 virtual void
92 postCheckpointSetup() override;
93
97 virtual void
98 registerRealms() override;
99
103 virtual void
104 registerOperators() override;
105
109 virtual void
110 parseRuntimeOptions() override = 0;
111
115 virtual void
116 prePlot() override;
117
118#ifdef CH_USE_HDF5
124 virtual void
125 writeCheckpointData(HDF5Handle& a_handle, const int a_lvl) const override;
126#endif
127
128#ifdef CH_USE_HDF5
134 virtual void
135 readCheckpointData(HDF5Handle& a_handle, const int a_lvl) override;
136#endif
137
143 virtual int
144 getNumberOfPlotVariables() const override;
145
150 virtual Vector<std::string>
151 getPlotVariableNames() const override;
152
160 virtual void
161 writePlotData(LevelData<EBCellFAB>& a_output,
162 int& a_icomp,
163 const std::string& a_outputRealm,
164 const int a_level) const override;
165
171 virtual Real
172 computeDt() override = 0;
173
179 virtual Real
180 advance(const Real a_dt) override = 0;
181
188 virtual void
189 synchronizeSolverTimes(const int a_step, const Real a_time, const Real a_dt) override;
190
196 virtual void
197 printStepReport() override;
198
204 virtual void
205 preRegrid(const int a_lmin, const int a_oldFinestLevel) override;
206
212 virtual void
213 preRegridInternals(const int a_lmin, const int a_oldFinestLevel);
214
221 virtual void
222 regrid(const int a_lmin, const int a_oldFinestLevel, const int a_newFinestLevel) override;
223
228 virtual void
229 postRegrid() override;
230
235 virtual void
236 computeJ(EBAMRCellData& a_J) const;
237
243 virtual void
244 computeElectricField(MFAMRCellData& a_electricFieldCell, const MFAMRCellData& a_potential) const;
245
252 virtual void
253 computeElectricField(EBAMRCellData& a_electricFieldCell, const phase::which_phase a_phase) const;
254
261 virtual void
262 computeElectricField(EBAMRCellData& a_E, const phase::which_phase a_phase, const MFAMRCellData& a_potential) const;
263
270 virtual void
271 computeElectricField(EBAMRFluxData& a_electricFieldFace,
272 const phase::which_phase a_phase,
273 const EBAMRCellData& a_electricFieldCell) const;
274
281 virtual void
282 computeElectricField(EBAMRIVData& a_electricFieldEB,
283 const phase::which_phase a_phase,
284 const EBAMRCellData& a_electricFieldCell) const;
285
291 virtual void
292 computeMaxElectricField(Real& a_maximumElectricField, const phase::which_phase a_phase);
293
298 virtual Real
299 getTime() const;
300
305 virtual bool
306 stationaryRTE();
307
311 virtual void
312 allocateInternals() = 0;
313
317 virtual void
318 deallocateInternals() = 0;
319
323 virtual void
324 parseOptions() = 0;
325
329 virtual void
330 deallocate();
331
337 virtual void
338 computeSpaceChargeDensity();
339
347 virtual void
348 computeSpaceChargeDensity(MFAMRCellData& a_rho, const Vector<EBAMRCellData*>& a_cdrDensities);
349
354 virtual void
355 computeCellConductivity(EBAMRCellData& a_cellConductivity) const;
356
363 virtual void
364 computeFaceConductivity(EBAMRFluxData& a_conductivityFace,
365 EBAMRIVData& a_conductivityEB,
366 const EBAMRCellData& a_conductivityCell) const;
367
374 virtual void
375 setupSemiImplicitPoisson(const Real a_dt);
376
388 virtual void
389 setupSemiImplicitPoisson(const EBAMRFluxData& a_conductivityFace,
390 const EBAMRIVData& a_conductivityEB,
391 const Real a_factor);
392
399 virtual bool
400 solvePoisson();
401
411 virtual bool
412 solvePoisson(MFAMRCellData& a_potential,
413 MFAMRCellData& a_rho,
414 const Vector<EBAMRCellData*>& a_cdrDensities,
415 const EBAMRIVData& a_sigma);
416
423 virtual void
424 advanceReactionNetwork(const Real a_time, const Real a_dt);
425
437 virtual void
438 advanceReactionNetwork(Vector<EBAMRCellData*>& a_cdrSources,
439 Vector<EBAMRCellData*>& a_rteSources,
440 const Vector<EBAMRCellData*>& a_cdrDensities,
441 const Vector<EBAMRCellData*>& a_rteDensities,
442 const EBAMRCellData& a_E,
443 const Real& a_time,
444 const Real& a_dt);
445
457 virtual void
458 advanceReactionNetwork(Vector<EBAMRCellData*>& a_cdrSources,
459 Vector<EBAMRCellData*>& a_rteSources,
460 const Vector<EBAMRCellData*>& a_cdrDensities,
461 const Vector<EBAMRCellData*>& a_cdrGradients,
462 const Vector<EBAMRCellData*>& a_rteDensities,
463 const EBAMRCellData& a_E,
464 const Real& a_time,
465 const Real& a_dt);
466
479 virtual void
480 advanceReactionNetwork(Vector<LevelData<EBCellFAB>*>& a_cdrSources,
481 Vector<LevelData<EBCellFAB>*>& a_rteSources,
482 const Vector<LevelData<EBCellFAB>*>& a_cdrDensities,
483 const Vector<LevelData<EBCellFAB>*>& a_cdrGradients,
484 const Vector<LevelData<EBCellFAB>*>& a_rteDensities,
485 const LevelData<EBCellFAB>& a_E,
486 const Real& a_time,
487 const Real& a_dt,
488 const int a_lvl);
502 virtual void
503 advanceReactionNetworkRegularCells(Vector<FArrayBox*>& a_cdrSources,
504 Vector<FArrayBox*>& a_rteSources,
505 const Vector<FArrayBox*>& a_cdrDensities,
506 const Vector<FArrayBox*>& a_cdrGradients,
507 const Vector<FArrayBox*>& a_rteDensities,
508 const FArrayBox& a_E,
509 const Real& a_time,
510 const Real& a_dt,
511 const Real& a_dx,
512 const Box& a_cellBox);
513
531 virtual void
532 advanceReactionNetworkIrreg(Vector<EBCellFAB*>& a_cdrSources,
533 Vector<EBCellFAB*>& a_rteSources,
534 const Vector<EBCellFAB*>& a_cdrDensities,
535 const Vector<EBCellFAB*>& a_cdrGradients,
536 const Vector<EBCellFAB*>& a_cdrVelocities,
537 const Vector<EBCellFAB*>& a_rteDensities,
538 const EBCellFAB& a_E,
539 const Real& a_time,
540 const Real& a_dt,
541 const Real& a_dx,
542 const Box& a_cellBox,
543 const int a_lvl,
544 const DataIndex& a_dit);
545
563 virtual void
564 advanceReactionNetworkIrregInterp(Vector<EBCellFAB*>& a_cdrSources,
565 Vector<EBCellFAB*>& a_rteSources,
566 const Vector<EBCellFAB*>& a_cdrDensities,
567 const Vector<EBCellFAB*>& a_cdrGradients,
568 const Vector<EBCellFAB*>& a_cdrVelocities,
569 const Vector<EBCellFAB*>& a_rteDensities,
570 const EBCellFAB& a_E,
571 const Real& a_time,
572 const Real& a_dt,
573 const Real& a_dx,
574 const Box& a_cellBox,
575 const int a_lvl,
576 const DataIndex& a_dit);
577
597 virtual void
598 advanceReactionNetworkIrregUpwind(Vector<EBCellFAB*>& a_cdrSources,
599 Vector<EBCellFAB*>& a_rteSources,
600 const Vector<EBCellFAB*>& a_cdrDensities,
601 const Vector<EBCellFAB*>& a_cdrGradients,
602 const Vector<EBCellFAB*>& a_cdrVelocities,
603 const Vector<EBCellFAB*>& a_rteDensities,
604 const EBCellFAB& a_E,
605 const Real& a_time,
606 const Real& a_dt,
607 const Real& a_dx,
608 const Box& a_box,
609 const int a_lvl,
610 const DataIndex& a_dit);
611
616 virtual void
617 computeCdrDiffusion();
618
625 virtual void
626 computeCdrDiffusion(const EBAMRCellData& a_electricFieldCell, const EBAMRIVData& a_electricFieldEB);
627
635 virtual void
636 computeCdrDiffusionCell(Vector<EBAMRCellData>& a_cdrDcoCell,
637 const Vector<EBAMRCellData*>& a_cdrDensities,
638 const EBAMRCellData& a_electricFieldCell,
639 const Real& a_time);
640
649 virtual void
650 computeCdrDiffusionCell(Vector<LevelData<EBCellFAB>*>& a_cdrDcoCell,
651 const Vector<LevelData<EBCellFAB>*>& a_cdrDensities,
652 const LevelData<EBCellFAB>& a_electricFieldCell,
653 const int a_lvl,
654 const Real& a_time);
655
665 virtual void
666 computeCdrDiffusionCellRegular(Vector<FArrayBox*>& a_cdrDcoCell,
667 const Vector<FArrayBox*>& a_cdrDensities,
668 const FArrayBox& a_electricFieldCell,
669 const Box a_cellBox,
670 const Real a_dx,
671 const Real a_time);
672
683 virtual void
684 computeCdrDiffusionCellIrregular(Vector<EBCellFAB*>& a_cdrDcoCell,
685 const Vector<EBCellFAB*>& a_cdrDensities,
686 const EBCellFAB& a_electricFieldCell,
687 const Real a_dx,
688 const Real& a_time,
689 const int a_lvl,
690 const DataIndex& a_dit);
691
701 virtual void
702 computeCdrDiffusionFace(Vector<EBAMRFluxData*>& a_cdrDcoFace,
703 const Vector<EBAMRCellData*>& a_cdrDensities,
704 const EBAMRCellData& a_electricFieldCell,
705 const Real& a_time);
706
714 virtual void
715 computeCdrDiffusionEb(Vector<EBAMRIVData*>& a_cdrDcoEB,
716 const Vector<EBAMRIVData*>& a_cdrDensitiesEB,
717 const EBAMRIVData& a_electricFieldEB,
718 const Real& a_time);
719
728 virtual void
729 computeCdrDiffusionEb(Vector<LevelData<BaseIVFAB<Real>>*>& a_cdrDcoEB,
730 const Vector<LevelData<BaseIVFAB<Real>>*>& a_cdrDensitiesEB,
731 const LevelData<BaseIVFAB<Real>>& a_electricFieldEB,
732 const Real& a_time,
733 const int a_lvl);
734
739 virtual void
740 computeCdrDriftVelocities();
741
749 virtual void
750 computeCdrDriftVelocities(Vector<EBAMRCellData*>& a_cdrVelocities,
751 const Vector<EBAMRCellData*>& a_cdrDensities,
752 const EBAMRCellData& a_E,
753 const Real& a_time);
754
764 virtual void
765 computeCdrDriftVelocities(Vector<LevelData<EBCellFAB>*>& a_cdrVelocities,
766 const Vector<LevelData<EBCellFAB>*>& a_cdrDensities,
767 const LevelData<EBCellFAB>& a_electricField,
768 const int a_lvl,
769 const Real& a_time);
770
781 virtual void
782 computeCdrDriftVelocitiesRegular(Vector<FArrayBox*>& a_cdrVelocities,
783 const Vector<FArrayBox*>& a_cdrDensities,
784 const FArrayBox& a_electricField,
785 const Box& a_cellBox,
786 const Real& a_time,
787 const Real& a_dx);
788
800 virtual void
801 computeCdrDriftVelocitiesIrregular(Vector<EBCellFAB*>& a_cdrVelocities,
802 const Vector<EBCellFAB*>& a_cdrDensities,
803 const EBCellFAB& a_electricField,
804 const Real& a_time,
805 const Real& a_dx,
806 const int a_lvl,
807 const DataIndex& a_dit);
808
821 virtual void
822 computeCdrFluxes(Vector<EBAMRIVData*>& a_cdrFluxesEB,
823 const Vector<EBAMRIVData*>& a_extrapCdrFluxes,
824 const Vector<EBAMRIVData*>& a_extrapCdrDensities,
825 const Vector<EBAMRIVData*>& a_extrapCdrVelocities,
826 const Vector<EBAMRIVData*>& a_extrapCdrGradients,
827 const Vector<EBAMRIVData*>& a_extrapRteFluxes,
828 const EBAMRIVData& a_electricField,
829 const Real& a_time);
830
845 virtual void
846 computeCdrFluxes(Vector<LevelData<BaseIVFAB<Real>>*>& a_cdrFluxesEB,
847 const Vector<LevelData<BaseIVFAB<Real>>*>& a_extrapCdrFluxes,
848 const Vector<LevelData<BaseIVFAB<Real>>*>& a_extrapCdrDensities,
849 const Vector<LevelData<BaseIVFAB<Real>>*>& a_extrapCdrVelocities,
850 const Vector<LevelData<BaseIVFAB<Real>>*>& a_extrapCdrGradients,
851 const Vector<LevelData<BaseIVFAB<Real>>*>& a_extrapRteFluxes,
852 const LevelData<BaseIVFAB<Real>>& a_electricField,
853 const Real& a_time,
854 const int a_lvl);
855
867 virtual void
868 computeCdrDomainFluxes(Vector<EBAMRIFData*>& a_cdrFluxes,
869 const Vector<EBAMRIFData*>& a_extrapCdrFluxes,
870 const Vector<EBAMRIFData*>& a_extrapCdrDensities,
871 const Vector<EBAMRIFData*>& a_extrapCdrVelocities,
872 const Vector<EBAMRIFData*>& a_extrapCdrGradients,
873 const Vector<EBAMRIFData*>& a_extrapRteFluxes,
874 const EBAMRIFData& a_electricField,
875 const Real& a_time);
876
890 virtual void
891 computeCdrDomainFluxes(Vector<LevelData<DomainFluxIFFAB>*> a_cdrFluxes,
892 const Vector<LevelData<DomainFluxIFFAB>*>& a_extrapCdrFluxes,
893 const Vector<LevelData<DomainFluxIFFAB>*>& a_extrapCdrDensities,
894 const Vector<LevelData<DomainFluxIFFAB>*>& a_extrapCdrVelocities,
895 const Vector<LevelData<DomainFluxIFFAB>*>& a_extrapCdrGradients,
896 const Vector<LevelData<DomainFluxIFFAB>*>& a_extrapRteFluxes,
897 const LevelData<DomainFluxIFFAB>& a_electricField,
898 const Real& a_time,
899 const int a_lvl);
900
908 virtual void
909 computeExtrapolatedFluxes(Vector<EBAMRIVData*>& a_extrapCdrFluxesEB,
910 const Vector<EBAMRCellData*>& a_cdrDensities,
911 const Vector<EBAMRCellData*>& a_cdrVelocities,
912 const phase::which_phase a_phase);
913
922 virtual void
923 computeExtrapolatedDomainFluxes(Vector<EBAMRIFData*>& a_cdrDomainFluxes,
924 const Vector<EBAMRCellData*>& a_cdrDensities,
925 const Vector<EBAMRCellData*>& a_cdrVelocities,
926 const phase::which_phase a_phase);
927
936 virtual void
937 computeExtrapolatedVelocities(Vector<EBAMRIVData*>& a_cdrVelocitiesEB,
938 const Vector<EBAMRCellData*>& a_cdrVelocitiesCell,
939 const phase::which_phase a_phase);
940
944 virtual void
945 deallocateSolverInternals();
946
953 virtual void
954 extrapolateToEb(Vector<EBAMRIVData*>& a_ebData,
955 const phase::which_phase a_phase,
956 const Vector<EBAMRCellData*>& a_cellData);
957
964 virtual void
965 extrapolateToEb(EBAMRIVData& a_ebData, const phase::which_phase a_phase, const EBAMRCellData& a_cellData);
966
974 virtual void
975 extrapolateToEb(LevelData<BaseIVFAB<Real>>& a_ebData,
976 const phase::which_phase a_phase,
977 const LevelData<EBCellFAB>& a_cellData,
978 const int a_lvl);
979
986 virtual void
987 extrapolateToDomainFaces(Vector<EBAMRIFData*>& a_domainData,
988 const phase::which_phase a_phase,
989 const Vector<EBAMRCellData*>& a_cellData);
990
997 virtual void
998 extrapolateToDomainFaces(EBAMRIFData& a_domainData,
999 const phase::which_phase a_phase,
1000 const EBAMRCellData& a_cellData);
1001
1009 virtual void
1010 extrapolateToDomainFaces(LevelData<DomainFluxIFFAB>& a_domainData,
1011 const phase::which_phase a_phase,
1012 const LevelData<EBCellFAB>& a_cellData,
1013 const int a_lvl);
1014
1022 virtual void
1023 extrapolateVectorToDomainFaces(Vector<EBAMRIFData*>& a_domainData,
1024 const phase::which_phase a_phase,
1025 const Vector<EBAMRCellData*>& a_cellData);
1026
1034 virtual void
1035 extrapolateVectorToDomainFaces(EBAMRIFData& a_domainData,
1036 const phase::which_phase a_phase,
1037 const EBAMRCellData& a_cellData);
1038
1046 virtual void
1047 extrapolateVelocitiesToDomainFaces(Vector<EBAMRIFData*>& a_domainVelocities,
1048 const phase::which_phase a_phase,
1049 const Vector<EBAMRCellData*>& a_cellVelocities);
1050
1056 virtual void
1057 getCdrMax(Real& a_cdrMax, std::string& a_solverName) const;
1058
1062 virtual void
1063 initialSigma();
1064
1070 virtual void
1071 projectFlux(EBAMRIVData& a_projectedFlux, const EBAMRIVData& a_flux);
1072
1079 virtual void
1080 projectFlux(LevelData<BaseIVFAB<Real>>& a_projectedFlux, const LevelData<BaseIVFAB<Real>>& a_flux, const int a_lvl);
1081
1087 virtual void
1088 projectDomain(EBAMRIFData& a_projectedFlux, const EBAMRIFData& a_flux);
1089
1096 virtual void
1097 regridInternals(const int a_lmin, const int a_oldFinestLevel, const int a_newFinestLevel) = 0;
1098
1105 virtual void
1106 regridSolvers(const int a_lmin, const int a_oldFinestLevel, const int a_newFinestLevel);
1107
1112 virtual void
1113 resetDielectricCells(EBAMRIVData& a_data) const;
1114
1118 virtual void
1119 sanityCheck() const;
1120
1124 virtual void
1125 setupCdr();
1126
1130 virtual void
1131 setupPoisson();
1132
1136 virtual void
1137 setupRadiativeTransfer();
1138
1142 virtual void
1143 setupSigma();
1144
1149 virtual void
1150 setSolverVerbosity();
1151
1155 virtual void
1156 solverDump();
1157
1162 virtual void
1163 setVoltage(std::function<Real(const Real a_time)> a_voltage);
1164
1170 virtual void
1171 solveRadiativeTransfer(const Real a_dt);
1172
1177 virtual void
1178 setCdrSolvers(RefCountedPtr<CdrLayout<CdrSolver>>& a_cdr);
1179
1184 virtual void
1185 setFieldSolver(RefCountedPtr<FieldSolver>& a_fieldSolver);
1186
1191 virtual void
1192 setRadiativeTransferSolvers(RefCountedPtr<RtLayout<RtSolver>>& a_rte);
1193
1200 virtual Real
1201 computeElectrodeCurrent();
1202
1208 virtual Real
1209 computeDielectricCurrent();
1210
1216 virtual Real
1217 computeDomainCurrent();
1218
1224 virtual Real
1225 computeOhmicInductionCurrent();
1226
1231 virtual Real
1232 computeRelaxationTime();
1233
1238 virtual Real
1239 getDt();
1240
1241protected:
1247 enum class SourceTermComputation
1248 {
1249 Interpolated,
1250 Upwind
1251 };
1252
1256 enum class TimeCode
1257 {
1258 Advection,
1259 Diffusion,
1260 AdvectionDiffusion,
1261 Source,
1262 RelaxationTime,
1263 Restricted,
1264 Hardcap,
1265 Error,
1266 Physics
1267 };
1268
1273 std::string m_realm;
1274
1278 std::string m_className;
1279
1283 phase::which_phase m_phase;
1284
1288 SourceTermComputation m_whichSourceTermComputation;
1289
1293 RefCountedPtr<MultiFluidIndexSpace> m_multifluidIndexSpace;
1294
1298 RefCountedPtr<CdrPlasmaPhysics> m_physics;
1299
1303 RefCountedPtr<CdrLayout<CdrSolver>> m_cdr;
1304
1308 RefCountedPtr<RtLayout<RtSolver>> m_rte;
1309
1313 RefCountedPtr<FieldSolver> m_fieldSolver;
1314
1318 RefCountedPtr<SurfaceODESolver<1>> m_sigma;
1319
1323 EBAMRCellData m_currentDensity;
1324
1328 EBAMRCellData m_physicsPlotVars;
1329
1333 std::function<Real(const Real a_time)> m_voltage;
1334
1338 Real m_minDt;
1339
1343 Real m_maxDt;
1344
1348 Real m_cfl;
1349
1353 Real m_relaxTime;
1354
1358 Real m_time;
1359
1363 Real m_dt;
1364
1368 Real m_dtCFL;
1369
1373 TimeCode m_timeCode;
1374
1378 int m_solverVerbosity;
1379
1383 int m_fastRTE;
1384
1388 int m_fastPoisson;
1389
1393 int m_upwindFactor;
1394
1398 virtual void
1399 parseVerbosity();
1400
1404 virtual void
1405 parseSolverVerbosity();
1406
1410 virtual void
1411 parseCFL();
1412
1416 virtual void
1417 parseRelaxationTime();
1418
1422 virtual void
1423 parseMinDt();
1424
1428 virtual void
1429 parseMaxDt();
1430
1434 virtual void
1435 parseFastPoisson();
1436
1440 virtual void
1441 parseFastRadiativeTransfer();
1442
1447 virtual void
1448 parseSourceComputation();
1449
1457 virtual void
1458 writeJ(LevelData<EBCellFAB>& a_output, int& a_icomp, const std::string& a_outputRealm, const int a_level) const;
1459
1464 virtual void
1465 computePhysicsPlotVars(EBAMRCellData& a_plotVars) const noexcept;
1466
1478 virtual void
1479 writeData(LevelData<EBCellFAB>& a_output,
1480 int& a_comp,
1481 const EBAMRCellData& a_data,
1482 const std::string& a_outputRealm,
1483 const int a_level,
1484 const bool a_interpToCentroids,
1485 const bool a_interpGhost) const noexcept;
1486};
1487} // namespace CdrPlasma
1488} // namespace Physics
1489
1490#include <CD_NamespaceFooter.H>
1491
1492#endif
CD_AmrMesh.H
Declaration of core class for handling AMR-related operations (with embedded boundaries)
CD_CdrLayout.H
Declaration of a class that holds a set of CdrSolvers (to cut down on typing).
CD_CdrPlasmaPhysics.H
Declaration of the Physics::CdrPlasma::CdrPlasmaPhysics interface class.
CD_ComputationalGeometry.H
Declaration of base class for defining geometries.
CD_FieldSolver.H
Contains declaration of a base electrostatics solver class.
CD_MultiFluidIndexSpace.H
Multi-fluid index space.
CD_RtLayout.H
Declaration of a class that holds a set of RtSolvers.
CD_SurfaceODESolver.H
Declaration of a cut-cell ODE solver.
CD_TimeStepper.H
Declaration of main (abstract) time stepper class.
CdrLayout
Class for holding a set of CdrSolvers.
Definition CD_CdrLayout.H:40
Physics::CdrPlasma::CdrPlasmaStepper
Abstract base class for advancing the coupled CDR plasma equations in time.
Definition CD_CdrPlasmaStepper.H:46
Physics::CdrPlasma::CdrPlasmaStepper::getTime
virtual Real getTime() const
Return the current simulation time.
Definition CD_CdrPlasmaStepper.cpp:4322
Physics::CdrPlasma::CdrPlasmaStepper::computeExtrapolatedFluxes
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.
Definition CD_CdrPlasmaStepper.cpp:2577
Physics::CdrPlasma::CdrPlasmaStepper::setupSigma
virtual void setupSigma()
Set up the surface charge solver.
Definition CD_CdrPlasmaStepper.cpp:3901
Physics::CdrPlasma::CdrPlasmaStepper::computeFaceConductivity
virtual void computeFaceConductivity(EBAMRFluxData &a_conductivityFace, EBAMRIVData &a_conductivityEB, const EBAMRCellData &a_conductivityCell) const
Compute the face-centered conductivity from a cell-centered conductivity.
Definition CD_CdrPlasmaStepper.cpp:234
Physics::CdrPlasma::CdrPlasmaStepper::setFieldSolver
virtual void setFieldSolver(RefCountedPtr< FieldSolver > &a_fieldSolver)
Set the field solver.
Definition CD_CdrPlasmaStepper.cpp:3286
Physics::CdrPlasma::CdrPlasmaStepper::computeExtrapolatedDomainFluxes
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.
Definition CD_CdrPlasmaStepper.cpp:2643
Physics::CdrPlasma::CdrPlasmaStepper::printStepReport
virtual void printStepReport() override
Print a step report.
Definition CD_CdrPlasmaStepper.cpp:4811
Physics::CdrPlasma::CdrPlasmaStepper::extrapolateToEb
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.
Definition CD_CdrPlasmaStepper.cpp:3072
Physics::CdrPlasma::CdrPlasmaStepper::extrapolateToDomainFaces
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.
Definition CD_CdrPlasmaStepper.cpp:3123
Physics::CdrPlasma::CdrPlasmaStepper::m_whichSourceTermComputation
SourceTermComputation m_whichSourceTermComputation
Which source term computation.
Definition CD_CdrPlasmaStepper.H:1288
Physics::CdrPlasma::CdrPlasmaStepper::projectFlux
virtual void projectFlux(EBAMRIVData &a_projectedFlux, const EBAMRIVData &a_flux)
Project a flux onto the EB normal. This is the AMR version.
Definition CD_CdrPlasmaStepper.cpp:3412
Physics::CdrPlasma::CdrPlasmaStepper::m_dtCFL
Real m_dtCFL
Computed CFL time step.
Definition CD_CdrPlasmaStepper.H:1368
Physics::CdrPlasma::CdrPlasmaStepper::setSolverVerbosity
virtual void setSolverVerbosity()
Set solver verbosity.
Definition CD_CdrPlasmaStepper.cpp:3838
Physics::CdrPlasma::CdrPlasmaStepper::computeCdrFluxes
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.
Definition CD_CdrPlasmaStepper.cpp:2152
Physics::CdrPlasma::CdrPlasmaStepper::m_solverVerbosity
int m_solverVerbosity
Verbosity for solvers.
Definition CD_CdrPlasmaStepper.H:1378
Physics::CdrPlasma::CdrPlasmaStepper::computeJ
virtual void computeJ(EBAMRCellData &a_J) const
Compute the current density.
Definition CD_CdrPlasmaStepper.cpp:2834
Physics::CdrPlasma::CdrPlasmaStepper::m_className
std::string m_className
Time stepper class name.
Definition CD_CdrPlasmaStepper.H:1278
Physics::CdrPlasma::CdrPlasmaStepper::m_relaxTime
Real m_relaxTime
Scaling factor for the dielectric relaxation time step restriction.
Definition CD_CdrPlasmaStepper.H:1353
Physics::CdrPlasma::CdrPlasmaStepper::parseFastPoisson
virtual void parseFastPoisson()
Parse the "fast poisson" method, i.e. how often we solve the Poisson equation.
Definition CD_CdrPlasmaStepper.cpp:3764
Physics::CdrPlasma::CdrPlasmaStepper::computeOhmicInductionCurrent
virtual Real computeOhmicInductionCurrent()
Compute induced current in external circuit due to ohmic conduction.
Definition CD_CdrPlasmaStepper.cpp:4238
Physics::CdrPlasma::CdrPlasmaStepper::m_sigma
RefCountedPtr< SurfaceODESolver< 1 > > m_sigma
Surface charge solver.
Definition CD_CdrPlasmaStepper.H:1318
Physics::CdrPlasma::CdrPlasmaStepper::advanceReactionNetworkRegularCells
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.
Definition CD_CdrPlasmaStepper.cpp:800
Physics::CdrPlasma::CdrPlasmaStepper::allocate
virtual void allocate() override
Allocate data for the time stepper and solvers.
Definition CD_CdrPlasmaStepper.cpp:3824
Physics::CdrPlasma::CdrPlasmaStepper::m_cdr
RefCountedPtr< CdrLayout< CdrSolver > > m_cdr
CDR solvers.
Definition CD_CdrPlasmaStepper.H:1303
Physics::CdrPlasma::CdrPlasmaStepper::parseMinDt
virtual void parseMinDt()
Parse the minimum allowed time step.
Definition CD_CdrPlasmaStepper.cpp:3716
Physics::CdrPlasma::CdrPlasmaStepper::regridInternals
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.
Physics::CdrPlasma::CdrPlasmaStepper::computeCdrDomainFluxes
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_electricField, const Real &a_time)
Compute CDR fluxes on domain faces for usage in boundary conditions. This is the AMR version.
Definition CD_CdrPlasmaStepper.cpp:2385
Physics::CdrPlasma::CdrPlasmaStepper::parseRuntimeOptions
virtual void parseRuntimeOptions() override=0
Parse runtime options. Subclasses must implement this.
Physics::CdrPlasma::CdrPlasmaStepper::m_multifluidIndexSpace
RefCountedPtr< MultiFluidIndexSpace > m_multifluidIndexSpace
Index space.
Definition CD_CdrPlasmaStepper.H:1293
Physics::CdrPlasma::CdrPlasmaStepper::getCdrMax
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.
Definition CD_CdrPlasmaStepper.cpp:3246
Physics::CdrPlasma::CdrPlasmaStepper::deallocate
virtual void deallocate()
Deallocation function. This will deallocate internal storage in the subclasses as well as the solvers...
Definition CD_CdrPlasmaStepper.cpp:4334
Physics::CdrPlasma::CdrPlasmaStepper::extrapolateVelocitiesToDomainFaces
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.
Definition CD_CdrPlasmaStepper.cpp:2777
Physics::CdrPlasma::CdrPlasmaStepper::sanityCheck
virtual void sanityCheck() const
Do a sanity check. Only works when DEBUG=TRUE.
Definition CD_CdrPlasmaStepper.cpp:3636
Physics::CdrPlasma::CdrPlasmaStepper::advance
virtual Real advance(const Real a_dt) override=0
Advance method, advances equations.
Physics::CdrPlasma::CdrPlasmaStepper::parseVerbosity
virtual void parseVerbosity()
Parse class verbosity.
Definition CD_CdrPlasmaStepper.cpp:3660
Physics::CdrPlasma::CdrPlasmaStepper::m_timeCode
TimeCode m_timeCode
Time code for step restriction.
Definition CD_CdrPlasmaStepper.H:1373
Physics::CdrPlasma::CdrPlasmaStepper::computeRelaxationTime
virtual Real computeRelaxationTime()
Compute the relaxation time as dt = eps0/conductivity.
Definition CD_CdrPlasmaStepper.cpp:4283
Physics::CdrPlasma::CdrPlasmaStepper::deallocateInternals
virtual void deallocateInternals()=0
Deallocate internals. Must be implemented by subclasses.
Definition CD_CdrPlasmaStepper.cpp:418
Physics::CdrPlasma::CdrPlasmaStepper::initialSigma
virtual void initialSigma()
Initialize the surface charge.
Definition CD_CdrPlasmaStepper.cpp:3353
Physics::CdrPlasma::CdrPlasmaStepper::computeCdrDriftVelocitiesRegular
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.
Definition CD_CdrPlasmaStepper.cpp:1987
Physics::CdrPlasma::CdrPlasmaStepper::setupCdr
virtual void setupCdr()
Setup the CDR solvers.
Definition CD_CdrPlasmaStepper.cpp:3852
Physics::CdrPlasma::CdrPlasmaStepper::computeDomainCurrent
virtual Real computeDomainCurrent()
Compute the ohmic current through the domain faces.
Definition CD_CdrPlasmaStepper.cpp:4153
Physics::CdrPlasma::CdrPlasmaStepper::solvePoisson
virtual bool solvePoisson()
Solve the electrostatic Poisson equation.
Definition CD_CdrPlasmaStepper.cpp:355
Physics::CdrPlasma::CdrPlasmaStepper::parseRelaxationTime
virtual void parseRelaxationTime()
Parse the relaxation time restriction.
Definition CD_CdrPlasmaStepper.cpp:3700
Physics::CdrPlasma::CdrPlasmaStepper::m_dt
Real m_dt
Previous time step size.
Definition CD_CdrPlasmaStepper.H:1363
Physics::CdrPlasma::CdrPlasmaStepper::advanceReactionNetworkIrregInterp
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_cellBox, const int a_lvl, const DataIndex &a_dit)
Compute reaction network sources. This version interpolates everything to the centroid and runs with ...
Definition CD_CdrPlasmaStepper.cpp:997
Physics::CdrPlasma::CdrPlasmaStepper::computeElectrodeCurrent
virtual Real computeElectrodeCurrent()
Compute the ohmic current on electrodes.
Definition CD_CdrPlasmaStepper.cpp:3968
Physics::CdrPlasma::CdrPlasmaStepper::~CdrPlasmaStepper
~CdrPlasmaStepper()
Destructor.
Definition CD_CdrPlasmaStepper.cpp:93
Physics::CdrPlasma::CdrPlasmaStepper::advanceReactionNetworkIrregUpwind
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.
Definition CD_CdrPlasmaStepper.cpp:1109
Physics::CdrPlasma::CdrPlasmaStepper::getDt
virtual Real getDt()
Get dt.
Definition CD_CdrPlasmaStepper.cpp:4328
Physics::CdrPlasma::CdrPlasmaStepper::CdrPlasmaStepper
CdrPlasmaStepper()
Empty constructor.
Definition CD_CdrPlasmaStepper.cpp:30
Physics::CdrPlasma::CdrPlasmaStepper::advanceReactionNetworkIrreg
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.
Definition CD_CdrPlasmaStepper.cpp:931
Physics::CdrPlasma::CdrPlasmaStepper::solverDump
virtual void solverDump()
All solvers writes plot files. These are not the same as the plot files written by Driver.
Definition CD_CdrPlasmaStepper.cpp:3918
Physics::CdrPlasma::CdrPlasmaStepper::parseFastRadiativeTransfer
virtual void parseFastRadiativeTransfer()
Parse the "fast rte" method, i.e. how often we solve the Poisson equation.
Definition CD_CdrPlasmaStepper.cpp:3748
Physics::CdrPlasma::CdrPlasmaStepper::computeCdrDiffusionCell
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.
Definition CD_CdrPlasmaStepper.cpp:1326
Physics::CdrPlasma::CdrPlasmaStepper::m_maxDt
Real m_maxDt
Maximum allowed time step.
Definition CD_CdrPlasmaStepper.H:1343
Physics::CdrPlasma::CdrPlasmaStepper::computeDt
virtual Real computeDt() override=0
Compute the time step.
Physics::CdrPlasma::CdrPlasmaStepper::advanceReactionNetwork
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 cal...
Definition CD_CdrPlasmaStepper.cpp:429
Physics::CdrPlasma::CdrPlasmaStepper::writeData
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.
Definition CD_CdrPlasmaStepper.cpp:4507
Physics::CdrPlasma::CdrPlasmaStepper::m_physicsPlotVars
EBAMRCellData m_physicsPlotVars
Storage for physics-observables.
Definition CD_CdrPlasmaStepper.H:1328
Physics::CdrPlasma::CdrPlasmaStepper::m_upwindFactor
int m_upwindFactor
Upwind factor.
Definition CD_CdrPlasmaStepper.H:1393
Physics::CdrPlasma::CdrPlasmaStepper::writePlotData
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.
Definition CD_CdrPlasmaStepper.cpp:4468
Physics::CdrPlasma::CdrPlasmaStepper::setupPoisson
virtual void setupPoisson()
Setup the field solver.
Definition CD_CdrPlasmaStepper.cpp:3868
Physics::CdrPlasma::CdrPlasmaStepper::computeCdrDiffusionCellRegular
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.
Definition CD_CdrPlasmaStepper.cpp:1450
Physics::CdrPlasma::CdrPlasmaStepper::m_fastRTE
int m_fastRTE
Solve RTE every m_fastRTE steps instead of every step. 0 or 1 means every step.
Definition CD_CdrPlasmaStepper.H:1383
Physics::CdrPlasma::CdrPlasmaStepper::m_currentDensity
EBAMRCellData m_currentDensity
Storage for the current.
Definition CD_CdrPlasmaStepper.H:1323
Physics::CdrPlasma::CdrPlasmaStepper::parseSourceComputation
virtual void parseSourceComputation()
Parse how we compute the source terms.
Definition CD_CdrPlasmaStepper.cpp:3780
Physics::CdrPlasma::CdrPlasmaStepper::setRadiativeTransferSolvers
virtual void setRadiativeTransferSolvers(RefCountedPtr< RtLayout< RtSolver > > &a_rte)
Set the RTE solvers.
Definition CD_CdrPlasmaStepper.cpp:3297
Physics::CdrPlasma::CdrPlasmaStepper::computeCdrDriftVelocitiesIrregular
virtual void computeCdrDriftVelocitiesIrregular(Vector< EBCellFAB * > &a_cdrVelocities, 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.
Definition CD_CdrPlasmaStepper.cpp:2056
Physics::CdrPlasma::CdrPlasmaStepper::parseMaxDt
virtual void parseMaxDt()
Parse the maximum allowed time step.
Definition CD_CdrPlasmaStepper.cpp:3732
Physics::CdrPlasma::CdrPlasmaStepper::initialData
virtual void initialData() override
Fill all solvers with initial data.
Definition CD_CdrPlasmaStepper.cpp:3330
Physics::CdrPlasma::CdrPlasmaStepper::m_realm
std::string m_realm
Realm where the time stepper is registered.
Definition CD_CdrPlasmaStepper.H:1273
Physics::CdrPlasma::CdrPlasmaStepper::computeCdrDiffusionFace
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.
Definition CD_CdrPlasmaStepper.cpp:1624
Physics::CdrPlasma::CdrPlasmaStepper::m_cfl
Real m_cfl
CFL number.
Definition CD_CdrPlasmaStepper.H:1348
Physics::CdrPlasma::CdrPlasmaStepper::computeElectricField
virtual void computeElectricField(MFAMRCellData &a_electricFieldCell, const MFAMRCellData &a_potential) const
Compute the cell-centered electric field on both phases.
Definition CD_CdrPlasmaStepper.cpp:2912
Physics::CdrPlasma::CdrPlasmaStepper::computeCdrDiffusion
virtual void computeCdrDiffusion()
Compute diffusion things using whatever is available in the solvers.
Definition CD_CdrPlasmaStepper.cpp:1241
Physics::CdrPlasma::CdrPlasmaStepper::preRegrid
virtual void preRegrid(const int a_lmin, const int a_oldFinestLevel) override
Pre-regrid function – this is the one called by Driver.
Definition CD_CdrPlasmaStepper.cpp:2807
Physics::CdrPlasma::CdrPlasmaStepper::postInitialize
virtual void postInitialize() override
Post-initialize operations to be performed after filling solvers with initialData.
Definition CD_CdrPlasmaStepper.cpp:48
Physics::CdrPlasma::CdrPlasmaStepper::m_fieldSolver
RefCountedPtr< FieldSolver > m_fieldSolver
Poisson solver.
Definition CD_CdrPlasmaStepper.H:1313
Physics::CdrPlasma::CdrPlasmaStepper::computeExtrapolatedVelocities
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.
Definition CD_CdrPlasmaStepper.cpp:2695
Physics::CdrPlasma::CdrPlasmaStepper::parseSolverVerbosity
virtual void parseSolverVerbosity()
Parse solver verbosities.
Definition CD_CdrPlasmaStepper.cpp:3672
Physics::CdrPlasma::CdrPlasmaStepper::computeDielectricCurrent
virtual Real computeDielectricCurrent()
Compute the ohmic current on dielectrics.
Definition CD_CdrPlasmaStepper.cpp:4061
Physics::CdrPlasma::CdrPlasmaStepper::computeMaxElectricField
virtual void computeMaxElectricField(Real &a_maximumElectricField, const phase::which_phase a_phase)
Compute the maximum of the electric field.
Definition CD_CdrPlasmaStepper.cpp:3034
Physics::CdrPlasma::CdrPlasmaStepper::m_physics
RefCountedPtr< CdrPlasmaPhysics > m_physics
Plasma kinetics.
Definition CD_CdrPlasmaStepper.H:1298
Physics::CdrPlasma::CdrPlasmaStepper::setupSolvers
virtual void setupSolvers() override
Instantiate solvers.
Definition CD_CdrPlasmaStepper.cpp:3308
Physics::CdrPlasma::CdrPlasmaStepper::stationaryRTE
virtual bool stationaryRTE()
Check if RTE solvers are stationary.
Definition CD_CdrPlasmaStepper.cpp:102
Physics::CdrPlasma::CdrPlasmaStepper::resetDielectricCells
virtual void resetDielectricCells(EBAMRIVData &a_data) const
Set data defined over dielectric cells to zero.
Definition CD_CdrPlasmaStepper.cpp:3589
Physics::CdrPlasma::CdrPlasmaStepper::computeCdrDiffusionCellIrregular
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.
Definition CD_CdrPlasmaStepper.cpp:1530
Physics::CdrPlasma::CdrPlasmaStepper::m_minDt
Real m_minDt
Minimum allowed time step.
Definition CD_CdrPlasmaStepper.H:1338
Physics::CdrPlasma::CdrPlasmaStepper::projectDomain
virtual void projectDomain(EBAMRIFData &a_projectedFlux, const EBAMRIFData &a_flux)
Project flux onto domain edges/faces.
Definition CD_CdrPlasmaStepper.cpp:3475
Physics::CdrPlasma::CdrPlasmaStepper::m_rte
RefCountedPtr< RtLayout< RtSolver > > m_rte
Radiative transfer solvers.
Definition CD_CdrPlasmaStepper.H:1308
Physics::CdrPlasma::CdrPlasmaStepper::setVoltage
virtual void setVoltage(std::function< Real(const Real a_time)> a_voltage)
Set voltage curve.
Definition CD_CdrPlasmaStepper.cpp:3649
Physics::CdrPlasma::CdrPlasmaStepper::SourceTermComputation
SourceTermComputation
Enum for switching between source term computations. 'Interpolated' interpolates the cell-centered da...
Definition CD_CdrPlasmaStepper.H:1248
Physics::CdrPlasma::CdrPlasmaStepper::postRegrid
virtual void postRegrid() override
Perform post-regrid operations.
Definition CD_CdrPlasmaStepper.cpp:70
Physics::CdrPlasma::CdrPlasmaStepper::setupRadiativeTransfer
virtual void setupRadiativeTransfer()
Setup the RTE solvers.
Definition CD_CdrPlasmaStepper.cpp:3884
Physics::CdrPlasma::CdrPlasmaStepper::m_phase
phase::which_phase m_phase
Plasma phase.
Definition CD_CdrPlasmaStepper.H:1283
Physics::CdrPlasma::CdrPlasmaStepper::computeSpaceChargeDensity
virtual void computeSpaceChargeDensity()
Compute the centroid-centered space charge density by using the data inside the CDR solvers.
Definition CD_CdrPlasmaStepper.cpp:113
Physics::CdrPlasma::CdrPlasmaStepper::computePhysicsPlotVars
virtual void computePhysicsPlotVars(EBAMRCellData &a_plotVars) const noexcept
Calculate the user-supplied plot variables.
Definition CD_CdrPlasmaStepper.cpp:4610
Physics::CdrPlasma::CdrPlasmaStepper::synchronizeSolverTimes
virtual void synchronizeSolverTimes(const int a_step, const Real a_time, const Real a_dt) override
Synchronize solver times.
Definition CD_CdrPlasmaStepper.cpp:3950
Physics::CdrPlasma::CdrPlasmaStepper::registerOperators
virtual void registerOperators() override
Register operators to be used for the simulation.
Definition CD_CdrPlasmaStepper.cpp:79
Physics::CdrPlasma::CdrPlasmaStepper::writeJ
virtual void writeJ(LevelData< EBCellFAB > &a_output, int &a_icomp, const std::string &a_outputRealm, const int a_level) const
Compute and put the current density in output data holder.
Definition CD_CdrPlasmaStepper.cpp:4569
Physics::CdrPlasma::CdrPlasmaStepper::m_voltage
std::function< Real(const Real a_time)> m_voltage
Applied voltage as a function of time. Used to set the Poisson solver boundary condition.
Definition CD_CdrPlasmaStepper.H:1333
Physics::CdrPlasma::CdrPlasmaStepper::postCheckpointSetup
virtual void postCheckpointSetup() override
Post-initialize operations to be performed after filling solvers with data read from checkpoint files...
Definition CD_CdrPlasmaStepper.cpp:4783
Physics::CdrPlasma::CdrPlasmaStepper::registerRealms
virtual void registerRealms() override
Register realms to be used for the simulation.
Definition CD_CdrPlasmaStepper.cpp:57
Physics::CdrPlasma::CdrPlasmaStepper::extrapolateVectorToDomainFaces
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.
Definition CD_CdrPlasmaStepper.cpp:2734
Physics::CdrPlasma::CdrPlasmaStepper::m_fastPoisson
int m_fastPoisson
Solve Poisson every m_fastPoisson steps instead of every step. 0 or 1 means every step.
Definition CD_CdrPlasmaStepper.H:1388
Physics::CdrPlasma::CdrPlasmaStepper::computeCdrDriftVelocities
virtual void computeCdrDriftVelocities()
Compute the CDR drift velocities.
Definition CD_CdrPlasmaStepper.cpp:1837
Physics::CdrPlasma::CdrPlasmaStepper::parseOptions
virtual void parseOptions()=0
Parse TimeStepper class options. Must be implemented by subclasses.
Physics::CdrPlasma::CdrPlasmaStepper::prePlot
virtual void prePlot() override
Pre-plot hook. Computes current density and physics plot variables for output.
Definition CD_CdrPlasmaStepper.cpp:4600
Physics::CdrPlasma::CdrPlasmaStepper::computeCdrDiffusionEb
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.
Definition CD_CdrPlasmaStepper.cpp:1702
Physics::CdrPlasma::CdrPlasmaStepper::allocateInternals
virtual void allocateInternals()=0
Allocate internal storage. Must be implemented by subclasses.
Definition CD_CdrPlasmaStepper.cpp:395
Physics::CdrPlasma::CdrPlasmaStepper::setCdrSolvers
virtual void setCdrSolvers(RefCountedPtr< CdrLayout< CdrSolver > > &a_cdr)
Set the CDR layout.
Definition CD_CdrPlasmaStepper.cpp:3276
Physics::CdrPlasma::CdrPlasmaStepper::regridSolvers
virtual void regridSolvers(const int a_lmin, const int a_oldFinestLevel, const int a_newFinestLevel)
Regrid all of our solvers.
Definition CD_CdrPlasmaStepper.cpp:3575
Physics::CdrPlasma::CdrPlasmaStepper::deallocateSolverInternals
virtual void deallocateSolverInternals()
Deallocate internal solver storages.
Definition CD_CdrPlasmaStepper.cpp:3058
Physics::CdrPlasma::CdrPlasmaStepper::computeCellConductivity
virtual void computeCellConductivity(EBAMRCellData &a_cellConductivity) const
Compute the cell-centered conductivity.
Definition CD_CdrPlasmaStepper.cpp:165
Physics::CdrPlasma::CdrPlasmaStepper::m_time
Real m_time
Current simulation time in seconds.
Definition CD_CdrPlasmaStepper.H:1358
Physics::CdrPlasma::CdrPlasmaStepper::solveRadiativeTransfer
virtual void solveRadiativeTransfer(const Real a_dt)
Solve the radiative transfer problem.
Definition CD_CdrPlasmaStepper.cpp:3931
Physics::CdrPlasma::CdrPlasmaStepper::getNumberOfPlotVariables
virtual int getNumberOfPlotVariables() const override
Get the number of plot variables for this time stepper.
Definition CD_CdrPlasmaStepper.cpp:4398
Physics::CdrPlasma::CdrPlasmaStepper::preRegridInternals
virtual void preRegridInternals(const int a_lmin, const int a_oldFinestLevel)
Do a preRegrid operation for internal storage (if needed).
Definition CD_CdrPlasmaStepper.cpp:2825
Physics::CdrPlasma::CdrPlasmaStepper::TimeCode
TimeCode
An enum for encapsulating how time steps were restricted.
Definition CD_CdrPlasmaStepper.H:1257
Physics::CdrPlasma::CdrPlasmaStepper::setupSemiImplicitPoisson
virtual void setupSemiImplicitPoisson(const Real a_dt)
Set up a semi-implicit Poisson solver.
Definition CD_CdrPlasmaStepper.cpp:281
Physics::CdrPlasma::CdrPlasmaStepper::getPlotVariableNames
virtual Vector< std::string > getPlotVariableNames() const override
Get plot variable names.
Definition CD_CdrPlasmaStepper.cpp:4433
Physics::CdrPlasma::CdrPlasmaStepper::regrid
virtual void regrid(const int a_lmin, const int a_oldFinestLevel, const int a_newFinestLevel) override
Regrid method.
Definition CD_CdrPlasmaStepper.cpp:3534
Physics::CdrPlasma::CdrPlasmaStepper::parseCFL
virtual void parseCFL()
Parse the CFL number.
Definition CD_CdrPlasmaStepper.cpp:3684
RtLayout
Class for holding a set of RtSolvers. T must derive from RtSolver.
Definition CD_RtLayout.H:28
TimeStepper
Base class for advancing equations.
Definition CD_TimeStepper.H:32
Physics
Namespace containing physics models for use with chombo-discharge.
Definition CD_AdvectionDiffusion.H:16
phase::which_phase
which_phase
Enumeration of supported phases.
Definition CD_MultiFluidIndexSpace.H:38
Generated by doxygen 1.9.8