createFaFields.H
Go to the documentation of this file.
1Info<< "Reading field Cs" << endl;
3(
4 IOobject
5 (
6 "Cs",
7 runTime.timeName(),
8 mesh,
9 IOobject::MUST_READ,
10 IOobject::AUTO_WRITE
11 ),
12 aMesh
13);
14
15dimensionedScalar Cs0("Cs0", dimMass/dimArea, 1.0);
16
17const areaVectorField& R = aMesh.areaCentres();
18
19Cs = Cs0*(1.0 + R.component(vector::X)/mag(R));
20
21dimensionedScalar Ds("Ds", dimViscosity, 1.0);
22
24(
25 IOobject
26 (
27 "Us",
28 runTime.timeName(),
29 mesh,
30 IOobject::NO_READ,
31 IOobject::NO_WRITE
32 ),
33 aMesh,
34 dimensionedVector(dimVelocity, Zero)
35);
36
37dimensioned<scalar> Uinf("Uinf", dimVelocity, 1.0);
38
39forAll(Us, faceI)
40{
41 Us[faceI].x() =
42 Uinf.value()*(0.25*(3.0 + sqr(R[faceI].x()/mag(R[faceI]))) - 1.0);
43
44 Us[faceI].y() =
45 Uinf.value()*0.25*R[faceI].x()*R[faceI].y()/sqr(mag(R[faceI]));
46
47 Us[faceI].z() =
48 Uinf.value()*0.25*R[faceI].x()*R[faceI].z()/sqr(mag(R[faceI]));
49}
50
51Us -= aMesh.faceAreaNormals()*(aMesh.faceAreaNormals() & Us);
52
53
54edgeScalarField phis
55(
56 IOobject
57 (
58 "phis",
59 runTime.timeName(),
60 mesh,
61 IOobject::NO_READ,
62 IOobject::NO_WRITE
63 ),
64 linearEdgeInterpolate(Us) & aMesh.Le()
65);
#define R(A, B, C, D, E, F, K, M)
dynamicFvMesh & mesh
engineTime & runTime
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
messageStream Info
Information stream (stdout output on master, null elsewhere)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:372
GeometricField< vector, faPatchField, areaMesh > areaVectorField
Definition: areaFieldsFwd.H:79
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
GeometricField< scalar, faPatchField, areaMesh > areaScalarField
Definition: areaFieldsFwd.H:78
dimensioned< vector > dimensionedVector
Dimensioned vector obtained from generic dimensioned type.
const scalarField & Cs
edgeScalarField phis(IOobject("phis", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), linearEdgeInterpolate(Us) &aMesh.Le())
Us
Info<< "Reading field Cs"<< endl;areaScalarField Cs(IOobject("Cs", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), aMesh);dimensionedScalar Cs0("Cs0", dimMass/dimArea, 1.0);const areaVectorField &R=aMesh.areaCentres();Cs=Cs0 *(1.0+R.component(vector::X)/mag(R));dimensionedScalar Ds("Ds", dimViscosity, 1.0);areaVectorField Us(IOobject("Us", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), aMesh, dimensionedVector(dimVelocity, Zero));dimensioned< scalar > Uinf("Uinf", dimVelocity, 1.0)
faMesh aMesh(mesh)
#define forAll(list, i)
Loop across all elements in list.
Definition: stdFoam.H:333