adding analysis capabilities for electrochemistry

analysis/ElectroChemistry.cpp

@@ -0,0 +1,29 @@
+#include "analysis/ElectroChemistry.h"
+
+void ElectroChemistryAnalyzer::ElectroChemistryAnalyzer(std::shared_ptr <Domain> Dm):
+	Dm(dm){
+	Nx=dm->Nx; Ny=dm->Ny; Nz=dm->Nz;
+	Volume=(Nx-2)*(Ny-2)*(Nz-2)*Dm->nprocx()*Dm->nprocy()*Dm->nprocz()*1.0;
+	
+	ChemicalPotential.resize(Nx,Ny,Nz);       ChemicalPotential.fill(0);
+	ElectricalPotential.resize(Nx,Ny,Nz);      ElectricalPotential.fill(0);
+	Pressure.resize(Nx,Ny,Nz);     	Pressure.fill(0);
+	Rho.resize(Nx,Ny,Nz);       	Rho.fill(0);
+	Vel_x.resize(Nx,Ny,Nz);         Vel_x.fill(0);	    // Gradient of the phase indicator field
+	Vel_y.resize(Nx,Ny,Nz);         Vel_y.fill(0);
+	Vel_z.resize(Nx,Ny,Nz);         Vel_z.fill(0);
+	SDs.resize(Nx,Ny,Nz);         	SDs.fill(0);
+	
+	DoubleArray Rho;	           // density field 
+	DoubleArray ChemicalPotential;	   // density field 
+	DoubleArray ElectricalPotential;	// density field 
+	DoubleArray Pressure; 	// pressure field
+	DoubleArray Vel_x;		// velocity field
+	DoubleArray Vel_y;
+	DoubleArray Vel_z;
+	DoubleArray SDs;
+}
+
+void ElectroChemistryAnalyzer::~ElectroChemistryAnalyzer(){
+	
+}

analysis/ElectroChemistry.h

@@ -0,0 +1,52 @@
+/*
+ * Sub-phase averaging tools
+ */
+
+#ifndef ElectroChem_INC
+#define ElectroChem_INC
+
+#include <vector>
+#include "common/Domain.h"
+#include "common/Communication.h"
+#include "analysis/analysis.h"
+#include "analysis/distance.h"
+#include "analysis/Minkowski.h"
+#include "common/Utilities.h"
+#include "common/MPI_Helpers.h"
+#include "IO/MeshDatabase.h"
+#include "IO/Reader.h"
+#include "IO/Writer.h"
+
+class ElectroChemistryAnalyzer{
+public:
+	std::shared_ptr <Domain> Dm;
+	double Volume;
+	// input variables
+	double rho_n, rho_w;
+	double nu_n, nu_w;
+	double gamma_wn, beta;
+	double Fx, Fy, Fz;
+
+    //...........................................................................
+    int Nx,Ny,Nz;
+	DoubleArray Rho;	           // density field 
+	DoubleArray ChemicalPotential;	   // density field 
+	DoubleArray ElectricalPotential;	// density field 
+	DoubleArray Pressure; 	// pressure field
+	DoubleArray Vel_x;		// velocity field
+	DoubleArray Vel_y;
+	DoubleArray Vel_z;
+	DoubleArray SDs;
+
+	ElectroChemistryAnalyzer(std::shared_ptr <Domain> Dm);
+	~ElectroChemistryAnalyzer();
+	
+	void SetParams();
+	void Basic();
+	void Write(int time);
+
+private:
+	FILE *TIMELOG;
+}
+#endif
+