|
[Sponsors] |
June 11, 2014, 10:33 |
InterDyMFoam modification - alpha1
|
#1 |
New Member
Join Date: Jun 2014
Posts: 7
Rep Power: 11 |
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 |
|
June 11, 2014, 11:45 |
internal face in blockMesh
|
#2 |
New Member
arough
Join Date: Mar 2013
Location: tunisia
Posts: 24
Rep Power: 13 |
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? |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
alpha1 grows unbounded in interDyMFoam | tayo | OpenFOAM Running, Solving & CFD | 15 | May 8, 2015 04:35 |
convection term alpha1 interDyMFoam | Elisabeth_ofoam | OpenFOAM Running, Solving & CFD | 0 | June 10, 2014 16: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 03:51 |