|
[Sponsors] |
May 4, 2015, 17:50 |
cyclic BC - modify boundary face values
|
#1 |
Member
Luca
Join Date: Mar 2011
Location: Italy
Posts: 62
Rep Power: 15 |
Dear All,
consider a straight duct with inlet and outlet defined as being cyclic patches. Inside of my own solver I am solving for a scalar transport equation for, say phi. I also define another volScalarField, say phi1, which is not solved but it is derived from the values of phi: Code:
forAll(mesh.cells(),cellI) phi1[cellI] = ( phi[cellI] + (B * mesh.C().internalField()[cellI]) + C ); Moreover I modify the boundary values of phi1 at the inlet and outlet patch faces as follows: Code:
const fvPatch& currPatchInlet = mesh.boundaryField()[patchInlet] forAll(currPatchInlet,faceI) phi1.boundaryField()[patchInlet][faceI] = phi.boundaryField()[patchID][faceI] + C const fvPatch& currPatchOutlet = mesh.boundaryField()[patchOutlet] forAll(currPatchOutlet,faceI) phi1.boundaryField()[patchOutlet][faceI] = phi.boundaryField()[patchOutlet][faceI] + C Can someone explain me why? Is this related to the cyclic BC or am I missing something? I hope my question is clear enough to get an answer. Thank you in advance, Luca |
|
Tags |
cyclic bc, face values |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[Other] mesh airfoil NACA0012 | anand_30 | OpenFOAM Meshing & Mesh Conversion | 13 | March 7, 2022 18:22 |
snappyhexmesh remove blockmesh geometry | philipp1 | OpenFOAM Running, Solving & CFD | 2 | December 12, 2014 11:58 |
Radiation interface | hinca | CFX | 15 | January 26, 2014 18:11 |
Pressure instability with rhoSimpleFoam | daniel_mills | OpenFOAM Running, Solving & CFD | 44 | February 17, 2011 18:08 |
[blockMesh] BlockMeshmergePatchPairs | hjasak | OpenFOAM Meshing & Mesh Conversion | 11 | August 15, 2008 08:36 |