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

Internal Faces to boundary patch conversion

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 8, 2010, 15:46
Default Internal Faces to boundary patch conversion
  #1
Member
 
Robin Gilbert
Join Date: Jan 2010
Posts: 66
Rep Power: 7
robingilbert is on a distinguished road
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
Attached Files
File Type: gz trialgeo2.geo.tar.gz (1.5 KB, 20 views)
robingilbert is offline   Reply With Quote

Old   March 8, 2010, 15:52
Default
  #2
Member
 
Robin Gilbert
Join Date: Jan 2010
Posts: 66
Rep Power: 7
robingilbert is on a distinguished road
Would splitMesh work fine if I use it instead of createBaffles?
robingilbert is offline   Reply With Quote

Old   March 23, 2010, 20:40
Default
  #3
New Member
 
Steffen
Join Date: Oct 2009
Posts: 7
Rep Power: 7
steve is on a distinguished road
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
steve is offline   Reply With Quote

Old   March 24, 2010, 08:33
Default calcMassFlow
  #4
Member
 
Michael Roth
Join Date: Mar 2009
Location: Guelph, Ontario, Canada
Posts: 46
Rep Power: 8
roth is on a distinguished road
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
roth is offline   Reply With Quote

Old   March 24, 2010, 11:40
Default
  #5
New Member
 
Steffen
Join Date: Oct 2009
Posts: 7
Rep Power: 7
steve is on a distinguished road
Thanks for the hint, Mike - I'll give it a try...

Kind Regards,
Steffen
steve is offline   Reply With Quote

Reply

Tags
boundary, gmsh, internal faces

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
Wind turbine simulation Saturn CFX 34 October 16, 2014 05:27
BlockMeshmergePatchPairs hjasak OpenFOAM Native Meshers: blockMesh 11 August 15, 2008 07:36
Only one boundary patch after gambitToFoam conversion bewuethr OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 3 August 22, 2007 10:30
boundary conditions for internal faces Tom Pendrey FLUENT 1 March 5, 2007 19:25
Trimmed cell and embedded refinement mesh conversion issues michele OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 2 July 15, 2005 04:15


All times are GMT -4. The time now is 15:02.