-
Notifications
You must be signed in to change notification settings - Fork 3
/
EnergyIntegrationOutput.h
69 lines (49 loc) · 2.21 KB
/
EnergyIntegrationOutput.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
#ifndef _ENERGY_INTEGRATION_OUTPUT_H_
#define _ENERGY_INTEGRATION_OUTPUT_H_
#include <IoData.h>
#include <SpaceVariable.h>
#include <LaserAbsorptionSolver.h>
#include <VarFcnBase.h>
#include<Vector5D.h>
/*********************************************************
* class EnergyIntegrationOperator integrates the energy
* required by user within the specified region
* ********************************************************/
using namespace std;
class EnergyIntegrationOutput
{
MPI_Comm &comm;
OutputData &iod_output;
std::vector<VarFcnBase*> &vf;
MeshData& iod_mesh;
EquationsData& iod_eqs;
LaserAbsorptionSolver* laser;
int frequency;
double frequency_dt;
double last_snapshot_time;
FILE *file[EnergyIntegrationData::SIZE]; //!< one file per solution variable
int numMaterials;
//! Mesh info
SpaceVariable3D& coordinates;
SpaceVariable3D& delta_xyz;
SpaceVariable3D& cell_volume;
public:
EnergyIntegrationOutput(MPI_Comm &comm_, IoData &iod, OutputData &iod_output_,
MeshData &iod_mesh_, EquationsData &iod_eqs_,
LaserAbsorptionSolver* laser_,
std::vector<VarFcnBase*> &vf_, SpaceVariable3D& coordinates_,
SpaceVariable3D& delta_xyz_, SpaceVariable3D& cell_volume_);
~EnergyIntegrationOutput();
void WriteSolutionOfIntegrationEnergy(double time, double dt, int time_step, SpaceVariable3D &V, SpaceVariable3D &ID,
SpaceVariable3D* L, bool force_write);
private:
void IntegrateVolume(SpaceVariable3D &ID, double* vol);
void IntegrateMass(SpaceVariable3D &V, SpaceVariable3D &ID, double* mass);
void IntegrateTotalEnergy(SpaceVariable3D &V, SpaceVariable3D &ID, double* E);
void IntegrateTotalEnthalpy(SpaceVariable3D &V, SpaceVariable3D &ID, double* H);
void IntegrateKineticEnergy(SpaceVariable3D &V, SpaceVariable3D &ID, double* kinetic);
void IntegrateInternalEnergy(SpaceVariable3D &V, SpaceVariable3D &ID, double* internal);
void IntegratePotentialEnergy(SpaceVariable3D &V, SpaceVariable3D &ID, double* potential);
void IntegrateLaserRadiation(SpaceVariable3D &V, SpaceVariable3D &ID, SpaceVariable3D *L, double* radiation);
};
#endif