CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM (http://www.cfd-online.com/Forums/openfoam/)
-   -   Internal Faces to boundary patch conversion (http://www.cfd-online.com/Forums/openfoam/73438-internal-faces-boundary-patch-conversion.html)

robingilbert March 8, 2010 15:46

Internal Faces to boundary patch conversion
 
1 Attachment(s)
Hello Foamers,

I am simulating a room with a raised platform which has an air inlet on the platform. the platform is an internal face. i created the geometry with gmsh and converted it to OpenFoam using gmsh2ToFoam . Since the internal faces are not recognized as boundary patches, I did the following:

1) modified the boundary file,
2) modified all the entries in 0 folder.
3) ran createBaffles.
4) moved the mesh from '0' folder to 'constant' folder

Now the required patch is showing in my boundary file. i applied the BC and on running the simulation, it is showing the following error:

Code:

Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package hPsiThermo<pureMixture<sutherlandTransport<specieThermo<hConstThermo<perfectGas>>>>>
Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RASModel
Selecting RAS turbulence model kEpsilon
Creating field DpDt


Starting time loop

Courant Number mean: 0 max: 0
Time = 1

#0  Foam::error::printStack(Foam::Ostream&) in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#1  Foam::sigFpe::sigFpeHandler(int) in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/lib/linux64GccDPOpt/libOpenFOAM.so"
#2  ?? in "/lib/libc.so.6"
#3  Foam::flowRateInletVelocityFvPatchVectorField::updateCoeffs() in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/lib/linux64GccDPOpt/libfiniteVolume.so"
#4  Foam::fvMatrix<Foam::Vector<double> >::fvMatrix(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&, Foam::dimensionSet const&) in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/lib/linux64GccDPOpt/libfiniteVolume.so"
#5  Foam::fv::gaussLaplacianScheme<Foam::Vector<double>, double>::fvmLaplacianUncorrected(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/lib/linux64GccDPOpt/libfiniteVolume.so"
#6  Foam::fv::gaussLaplacianScheme<Foam::Vector<double>, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/lib/linux64GccDPOpt/libfiniteVolume.so"
#7  Foam::fv::laplacianScheme<Foam::Vector<double>, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/lib/linux64GccDPOpt/libfiniteVolume.so"
#8  Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam::fvm::laplacian<Foam::Vector<double>, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&, Foam::word const&) in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#9  Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam::fvm::laplacian<Foam::Vector<double>, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#10  Foam::compressible::RASModels::kEpsilon::divDevRhoReff(Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>&) const in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#11  main in "/home/robingilbert/OpenFOAM/robingilbert-1.6.x/applications/bin/linux64GccDPOpt/rhoPimpleFoam"
#12  __libc_start_main in "/lib/libc.so.6"
#13  _start at /build/buildd/eglibc-2.10.1/csu/../sysdeps/x86_64/elf/start.S:116
Floating point exception

I followed the steps in:

http://www.cfd-online.com/Forums/ope...-openfoam.html

I cant understand why we need to run 'faceSet'. If at all i need to run it, can you please guide me as to how i should set up the 'faceSetDict'?

I am attaching my .geo file here, just in case

robingilbert March 8, 2010 15:52

Would splitMesh work fine if I use it instead of createBaffles?

steve March 23, 2010 20:40

Hello,

I am facing a similar problem...
I want to calculate the flux through an internal face inside my mesh. In principal I did the same steps as Robin to transform the faceZone into a patch...
I think I will manage to calculate the flux through that patch if only I would know which type I have to assign to this patch. This is where I am stuck. Is there something like a transmissive patch type which does not influence the flow? Are there other ideas how to solve this problem?

Thanks and kind regards,
Steffen

roth March 24, 2010 08:33

calcMassFlow
 
Consider using Bernhard Gschaider's calcMassFlow:

http://openfoamwiki.net/index.php/Contrib_calcMassFlow

This utility operates on faceSets or patches interchangeably. I seem to recall it compiling in 1.5.x without much bother.

Mike

steve March 24, 2010 11:40

Thanks for the hint, Mike - I'll give it a try...

Kind Regards,
Steffen


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