CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

InterDyMFoam modification - alpha1

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   June 11, 2014, 09:33
Default InterDyMFoam modification - alpha1
  #1
New Member
 
Join Date: Jun 2014
Posts: 7
Rep Power: 4
Elisabeth_ofoam is on a distinguished road
Dear OpenFoam users,

I am a new OpenFoam user and I need your help.
I try to remove the convective term on the air side for the floatingObject case, like in the example here:
http://www.imac.unavarra.es/Cursos/A...AMWorkshop.pdf
I cannot simply add alpha1 like in the example, because I use openfoam 2.3.0 and the files for interFoam and interDyMFoam look differently.

old version 2.1.1
surfaceScalarField muEff
(
"muEff",
twoPhaseProperties.muf()
+ fvc::interpolate(rho*turbulence->nut())
);

fvVectorMatrix UEqn
(
fvm::ddt(rho, U)
+ fvm::div(rhoPhi, U)
- fvm::laplacian(muEff, U)
- (fvc::grad(U) & fvc::grad(muEff))
//- fvc::div(muEff*(fvc::interpolate(dev(fvc::grad(U)) ) & mesh.Sf()))
);

UEqn.relax();

if (pimple.momentumPredictor())
{
solve
(
UEqn
==
fvc::reconstruct
(
(
fvc::interpolate(interface.sigmaK())*fvc::snGrad(a lpha1)
- ghf*fvc::snGrad(rho)
- fvc::snGrad(p_rgh)
) * mesh.magSf()
)
);
}

new version 2.3.0

fvVectorMatrix UEqn
(
fvm::ddt(rho, U)
+ fvm::div(rhoPhi, U)
+ turbulence->divDevRhoReff(rho, U)
);

UEqn.relax();

if (pimple.momentumPredictor())
{
solve
(
UEqn
==
fvc::reconstruct
(
(
fvc::interpolate(interface.sigmaK())*fvc::snGrad(a lpha1)
- ghf*fvc::snGrad(rho)
- fvc::snGrad(p_rgh)
) * mesh.magSf()
)
);
}

How can I solve this problem?
Thank you in advance for your help!
Elisabeth
Elisabeth_ofoam is offline   Reply With Quote

Old   June 11, 2014, 10:45
Default internal face in blockMesh
  #2
rou
New Member
 
arough
Join Date: Mar 2013
Location: tunisia
Posts: 13
Rep Power: 5
rou is on a distinguished road
Hi Foamers,

I have tried to run two concentric cylinders but when I define patches I have an error message.
Error messages:
Trying to specify a boundary face 4(4 5 17 16) on the face on cell 1 which is either an internal face or already belongs to some other patch. This is face 0 of patch 5 named evap.
From function polyMesholyMesh
(
const IOobject& io,
const pointField& points,
const cellShapeList& cellsAsShapes,
const faceListList& boundaryFaces,
const wordList& boundaryPatchTypes,
const wordList& boundaryPatchNames,
const word& defaultBoundaryPatchType
)
in file C:\tmp\OpenFOAM-1.5\src\OpenFOAM\meshes\polyMesh\polyMeshFromShape Mesh.C at line 483.
FOAM aborting


could anybody help me to know how we can define an internal face in blockMesh?
rou is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
alpha1 grows unbounded in interDyMFoam tayo OpenFOAM Running, Solving & CFD 15 May 8, 2015 03:35
convection term alpha1 interDyMFoam Elisabeth_ofoam OpenFOAM Running, Solving & CFD 0 June 10, 2014 15:57
Negative alpha1 using interDyMFoam jrrygg OpenFOAM Running, Solving & CFD 14 March 2, 2013 11:26
same geometry,structured and unstructured mesh,different behaviour. sharonyue OpenFOAM Running, Solving & CFD 13 January 2, 2013 23:40
alpha or alpha1 in interFoam & interDyMFoam wavytracy OpenFOAM Bugs 3 September 10, 2009 02:51


All times are GMT -4. The time now is 17:58.