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

How to Produce a Multiregion Porous Media - simpleFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 3, 2013, 12:19
Default How to Produce a Multiregion Porous Media - simpleFoam
  #1
New Member
 
David Ouellette
Join Date: Feb 2013
Location: Ottawa, Ontario, Canada
Posts: 12
Rep Power: 4
DavidOuellette is on a distinguished road
Hey all,

I recently setup a simple case where I have flow between parallel plates and the media between these plates is a porous material. The flow is also governed by Darcy's Law. Now, I want to expand this case to have flow through multiple porous layers, where each porous layer has different properties. Unfortunately, I'm having trouble with the 'easy case' of creating a second layer.

I played around with chtMultiRegionSimpleFoam as well as the tutorial case of multiRegionHeater. However, when I ran the tutorial as is, without any modification, it did not converge. I also did not get convergence when I play around with mesh size and boundary conditions. I also noticed that a lot of people seem to have trouble with the same solver. My question though, has anyone tried making a multilayer porous media with this solver (chtMultiRegionSimpleFoam)? One of the challenges that I can foresee is that the fluid and solid regions will coincide. Can chtMultiRegionSimpleFoam handle that?

I was also looking into Maria Garcia Camprubi's PhD thesis (http://gfn.cps.unizar.es/?q=node/234) where she modeled a solid oxide fuel cell using OpenFOAM (her version of OF was part of the ext branch, mine is OF 2.1.1). She explained that we can map properties to different regions within a mesh. Does anyone have experience with this technique or know of any good references or start points for this technique? Also, does anyone know of any other possible ways to have multiple regions with different properties?

Thanks,
David
DavidOuellette is offline   Reply With Quote

Old   March 3, 2013, 13:01
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,301
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings David,

If I understand you correctly, what you want are only multiple porous zones, but with not clear separation between fluid zones (i.e. the same fluid goes across all porous zones). If this is the case, check the example I made some time ago: Multiple regions with a porous zone post #5

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   March 3, 2013, 13:14
Default
  #3
New Member
 
David Ouellette
Join Date: Feb 2013
Location: Ottawa, Ontario, Canada
Posts: 12
Rep Power: 4
DavidOuellette is on a distinguished road
Hi Bruno,

Thanks for your quick reply! I just downloaded the suggested files from your post and I will begin playing around with the model soon. I'll let you know how it all turns out.

Thanks again for your help Bruno. Much appreciated!
David
DavidOuellette is offline   Reply With Quote

Old   March 6, 2013, 17:54
Default
  #4
New Member
 
David Ouellette
Join Date: Feb 2013
Location: Ottawa, Ontario, Canada
Posts: 12
Rep Power: 4
DavidOuellette is on a distinguished road
Hey Bruno and fellow Foamers,

I started playing with snappyHexMesh (sHM) where the geometry is a simple cube. The tutorial that I followed is shown in this link (https://sites.google.com/site/snappy...mesh/cube-case), where the files can be found in: https://docs.google.com/file/d/0B2KY...ZDg/edit?pli=1

After running the cases and playing around with the geometry, I now want to create a 2D mesh so it better represents my scenario. Unfortunately, when I try this, my solver doesn't agree and crashes after viewing the results in ParaView. The terminal gives the error below.

HTML Code:
--> FOAM FATAL IO ERROR: 
size 500 is not equal to the given value of 4000

file: /home/fuelcell-laptop/OpenFOAM/fuelcell-laptop-2.1.1/run/FuelCellModels/porousParallelPlatesStl/81/p from line 18 to line 544.

    From function Field<Type>::Field(const word& keyword, const dictionary&, const label)
    in file /home/opencfd/OpenFOAM/OpenFOAM-2.1.1/src/OpenFOAM/lnInclude/Field.C at line 236.
So, my question is, does anyone know of anyway to have sHM create only one cell in the z-direction, and have the required number of cells in the x and y-directions?

Thanks,
David
DavidOuellette is offline   Reply With Quote

Old   March 6, 2013, 18:08
Default
  #5
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,301
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Quick answer: http://openfoamwiki.net/index.php/Sn...rate_2D_meshes
wyldckat is offline   Reply With Quote

Old   March 6, 2013, 18:23
Default
  #6
New Member
 
David Ouellette
Join Date: Feb 2013
Location: Ottawa, Ontario, Canada
Posts: 12
Rep Power: 4
DavidOuellette is on a distinguished road
Awesome! Thanks Bruno! I'm not sure how I missed that page when I originally checked OpenFoam Wiki. But nonetheless, thanks!

David
DavidOuellette 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
Porous media setup issues in Fluent Bernard Van FLUENT 26 July 24, 2015 18:24
How to model granular flow through porous media Axius FLUENT 2 August 7, 2014 10:34
Multiphase Porous Media Flow - Convergence Issues VT_Bromley FLUENT 5 March 23, 2011 11:38
species mass source in porous media ? PK FLUENT 0 February 16, 2007 12:12
porous media: Fluent or Star-CD? Igor Main CFD Forum 0 December 5, 2002 16:16


All times are GMT -4. The time now is 23:20.