electrostaticFoam.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2011-2023 OpenFOAM Foundation
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8 License
9  This file is part of OpenFOAM.
10 
11  OpenFOAM is free software: you can redistribute it and/or modify it
12  under the terms of the GNU General Public License as published by
13  the Free Software Foundation, either version 3 of the License, or
14  (at your option) any later version.
15 
16  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
19  for more details.
20 
21  You should have received a copy of the GNU General Public License
22  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
23 
24 Application
25  electrostaticFoam
26 
27 Description
28  Solver for electrostatics.
29 
30 \*---------------------------------------------------------------------------*/
31 
32 #include "argList.H"
33 
34 #include "fvMesh.H"
35 
36 #include "fvcSnGrad.H"
37 
38 #include "fvmDdt.H"
39 #include "fvmDiv.H"
40 #include "fvmLaplacian.H"
41 
42 using namespace Foam;
43 
44 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
45 
46 int main(int argc, char *argv[])
47 {
48  #include "setRootCase.H"
49 
50  #include "createTime.H"
51  #include "createMesh.H"
52  #include "createFields.H"
53 
54  // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
55 
56  Info<< "\nStarting iteration loop\n" << endl;
57 
58  while (runTime.loop())
59  {
60  Info<< "Iteration = " << runTime.name() << nl << endl;
61 
62  solve
63  (
65  );
66 
67  rhoFlux = -k*mesh.magSf()*fvc::snGrad(phi);
68 
69  solve
70  (
71  fvm::ddt(rho) + fvm::div(rhoFlux, rho)
72  );
73 
74  runTime.write();
75 
76  Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
77  << " ClockTime = " << runTime.elapsedClockTime() << " s"
78  << nl << endl;
79  }
80 
81  Info<< "End\n" << endl;
82 
83  return 0;
84 }
85 
86 
87 // ************************************************************************* //
label k
int main(int argc, char *argv[])
Calculate the snGrad of the given volField.
Calculate the matrix for the first temporal derivative.
Calculate the matrix for the divergence of the given field and flux.
Calculate the matrix for the laplacian of the field.
const dimensionedScalar epsilon0
Electric constant: default SI units: [F/m].
tmp< SurfaceField< Type > > snGrad(const VolField< Type > &vf, const word &name)
Definition: fvcSnGrad.C:45
tmp< fvMatrix< Type > > laplacian(const VolField< Type > &vf, const word &name)
Definition: fvmLaplacian.C:47
tmp< fvMatrix< Type > > div(const surfaceScalarField &flux, const VolField< Type > &vf, const word &name)
Definition: fvmDiv.C:48
tmp< fvMatrix< Type > > ddt(const VolField< Type > &vf)
Definition: fvmDdt.C:46
Namespace for OpenFOAM.
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:257
messageStream Info
static const char nl
Definition: Ostream.H:266
SolverPerformance< Type > solve(fvMatrix< Type > &, const word &)
Solve returning the solution statistics given convergence tolerance.