public ProjectFile()
{
fields = new List <FieldFile>();
boundaryNames = new List <string>();
otherFiles = new List <FoamFile>();
controlDict = new ControlDict();
transportProperties = new TransportProperties();
turbulenceProperties = new TurbulenceProperties(new KEpsilon());
material = new Material("Steel");
material.Density = 6940.0;
material.DynamicViscosity = 0.006293;
fvSchemes = new FvSchemes();
fvSolution = new FvSolution(new OpenCFD.Solutions.Simple());
}
static void test1()
{
FieldFile ff = new FieldFile("U");
ff.Dimensions = Dimension.U;
ff.InternalField = new InternalField(new VectorVar(1, 0, 0, true));
ff.BoundaryField = new BoundaryField();
ff.BoundaryField.AddBoundary(new FixedValue("in", new VectorVar(1, 0, 0, true)));
ff.BoundaryField.AddBoundary(new Symmetry("sym"));
//ff.Write("e:\\hope\\test");
ControlDict cd = new ControlDict();
//cd.Write("e:\\hope\\test");
FvSchemes fs = new FvSchemes();
Interpolation i = new Interpolation(Interpolation.InterpolationTypes.linear);
Grad g = new Grad(Grad.GradTypes.Gauss, i);
SnGrad s = new SnGrad(SnGrad.SnGradTypes.corrected);
Laplacian l = new Laplacian(Laplacian.LaplacianTypes.Gauss, i, s);
Div d = new Div(Div.DivTypes.Gauss, i);
fs.AddDdtSchemes(new Ddt(Ddt.DdtTypes.steadyState));
fs.AddGradSchemes(g);
fs.AddDivSchemes(d);
fs.AddDivSchemes("div(phi,U)", d);
fs.AddLaplacianSchemes(l);
fs.AddInterpolationSchemes(i);
fs.AddSnGradSchemes(s);
//fs.Write("e:\\hope\\test");
FvSolution fsn = new FvSolution(new Simple());
PBiCG psu = new PBiCG("U", 0.001, 0.001, PBiCG.Preconditioners.DIC);
PBiCG psp = new PBiCG("p", 0.001, 0.001, PBiCG.Preconditioners.DIC);
fsn.AddSolver(psu);
fsn.AddSolver(psp);
fsn.SetFieldRelaxtionFactors("U", 1);
fsn.SetFieldRelaxtionFactors("U", 0.8);
//fsn.Write("e:\\hope\\test");
}
