resetBoundaries.H
Go to the documentation of this file.
1{
2 // Keep standard formulation on domain boundaries to ensure compatibility
3 // with existing boundary conditions
5 (
6 fvc::interpolate(rho.boundaryField()*rAU.boundaryField())
7 );
8
10 (
11 fvc::interpolate(rho.boundaryField()*HbyA.boundaryField())
12 );
13
14 surfaceScalarField::Boundary& rhorAUfbf = rhorAUf.boundaryFieldRef();
15 surfaceVectorField::Boundary& rhoHbyAfbf = rhoHbyAf.boundaryFieldRef();
16
17 forAll(U.boundaryField(), patchi)
18 {
19 if (!U.boundaryField()[patchi].coupled())
20 {
21 rhorAUfbf[patchi] = rhorAUf_orig[patchi];
22 rhoHbyAfbf[patchi] = rhoHbyA_orig[patchi];
23 }
24 }
25}
A field of fields is a PtrList of fields with reference counting.
Definition: FieldField.H:80
U
Definition: pEqn.H:72
surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho *rAU))
HbyA
Definition: pcEqn.H:74
surfaceVectorField rhoHbyAf("rhoHbyAf", fvc::interpolate(rho) *fvc::interpolate(U)+rhorAUf *fvc::interpolate(fvc::grad(p)))
tmp< volScalarField > rAU
Definition: initCorrectPhi.H:1
surfaceVectorField::Boundary & rhoHbyAfbf
const Foam::FieldField< Foam::fvsPatchField, vector > rhoHbyA_orig(fvc::interpolate(rho.boundaryField() *HbyA.boundaryField()))
surfaceScalarField::Boundary & rhorAUfbf
#define forAll(list, i)
Loop across all elements in list.
Definition: stdFoam.H:333