All Classes Namespaces Files Functions Variables Typedefs Friends Macros Groups
createFaceFields.H
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | www.openfoam.com
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8  Copyright (C) 2011-2016 OpenFOAM Foundation
9  Copyright (C) 2019 OpenCFD Ltd.
10  Copyright (C) 2016-2019 ISP RAS (www.ispras.ru) UniCFD Group (www.unicfd.ru)
11 -------------------------------------------------------------------------------
12 License
13  This file is part of QGDsolver library, based on OpenFOAM+.
14 
15  OpenFOAM is free software: you can redistribute it and/or modify it
16  under the terms of the GNU General Public License as published by
17  the Free Software Foundation, either version 3 of the License, or
18  (at your option) any later version.
19 
20  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
21  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
22  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
23  for more details.
24 
25  You should have received a copy of the GNU General Public License
26  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
27 
28 Global
29  createFaceFields
30 
31 Description
32  Creates the face fields: linear interpolation of fields from volumes
33  to face centers.
34 
35 \*---------------------------------------------------------------------------*/
36 
37 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
38 
39 // Density
40 surfaceScalarField rhof
41 (
42  "rhof",
43  linearInterpolate(rho)
44 );
45 
46 // Velocity
47 surfaceVectorField Uf
48 (
49  "Uf",
50  linearInterpolate(U)
51 );
52 
53 surfaceScalarField pf
54 (
55  "pf",
56 linearInterpolate(p)
57 );
58 
59 //
60 surfaceVectorField Wf
61 (
62  "Wf",
63  linearInterpolate(W)
64 );
65 
66 surfaceVectorField BdFrcf
67 (
68  "BdFrcf",
69  linearInterpolate(BdFrc)
70 );
71 
72 surfaceScalarField Tf
73 (
74  "Tf",
75 linearInterpolate(T)
76 );
77 
78 // Heat capacity at constant pressure
79 surfaceScalarField Cpf
80 (
81  "Cpf",
82  linearInterpolate(thermo.Cp())
83 );
84 
85 surfaceScalarField muf
86 (
87  "muf",
88  linearInterpolate(thermo.mu()*0.0)
89 );
90 
91 /*
92  *
93  * QGD coefficients
94  *
95  */
96 
97 // OpenFOAM delta coeffs can be negative for coupled patches
98 surfaceScalarField taubyrhof
99 (
100  "taubyrhof",
101  tauQGDf/rhof
102 );
103 
104 surfaceScalarField alphaf
105 (
106  "alphaf",
107  linearInterpolate(thermo.alpha())
108 );
109 
110 
111 // ************************************************************************* //
Cpf
Definition: updateFields.H:42
volVectorField W(IOobject("W", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), U)
pf
Definition: updateFields.H:21
BdFrcf
Definition: updateFields.H:36
const surfaceScalarField & tauQGDf
Definition: createFields.H:55
muf
Definition: updateFields.H:38
volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh)
alphaf
Definition: updateFields.H:18
taubyrhof
Definition: updateFluxes.H:11
rhof
Definition: updateFields.H:27
Tf
Definition: updateFields.H:33
BdFrc
Definition: updateFields.H:35
rhoQGDThermo & thermo
Definition: createFields.H:47
Wf
Definition: QHDUEqn.H:31
Uf
Definition: updateFields.H:30
volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho())
volScalarField & T
Definition: createFields.H:53
volScalarField & p
Definition: createFields.H:52