|
[Sponsors] |
January 19, 2013, 11:43 |
polyMesh files, Elements?
|
#1 |
New Member
---
Join Date: Jan 2013
Posts: 17
Rep Power: 13 |
Hi you all!
Pleasure to post here (my first post ). I'm trying to do some freeSurface simulations with OpenFoam, I do the set-up as well as I can but the thing comes when I do have to set where there's VOF (alpha = 1). I'm following the damBreak example which is more or less the same, but the problem comes up when I want to generate the alpha list. what I want to do is to set some elements on a certain zone with fluid. to do so what I need to do is to know which are the elements on this zone and set the alpha value equal to 1. am I right? So the question to come is, where the jeck is the list of elements inside polyMesh folder?, the list of alphas in damBreak is as long as nFaces, but I'm working on a 3D problem, there's anyway to know how many alphas do I have to set? are they related with the elements? Thank you, |
|
January 19, 2013, 17:44 |
|
#2 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,981
Blog Entries: 45
Rep Power: 128 |
Greetings gironidegirona and welcome to the forum!
Look into the file "system/setFieldsDict" and you'll learn how you can select a zone in the mesh domain. When you run the command setFields it will use that dictionary file for applying your requests. For more example files, run these command: Code:
find $FOAM_TUTORIALS $FOAM_UTILITIES -name setFieldsDict Bruno
__________________
|
|
January 29, 2013, 06:10 |
|
#3 |
New Member
---
Join Date: Jan 2013
Posts: 17
Rep Power: 13 |
Thanks wyldckat ,
Thanks for the fast reply, a couple of hours after posting here I realized how to use the setFieldsDict, but the thing I'm facing now is different. I can start a simulation with some predefined alphas, but during the simulation I would like to reduce the number of alphas or to increase it in determined places. The easiest thing for me is to know the relation with the alpha list and the points of that determine the cell, that way I'll do my check "is this element in the redion I'm checking" and do the operations I need to do? I hope I've explained myself. Thanks, Sergi |
|
January 29, 2013, 19:11 |
|
#4 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,981
Blog Entries: 45
Rep Power: 128 |
Greetings Sergi,
That's an interesting kind of simulation you're trying to perform. But AFAIK, not very simple to put into practice. From what I know, you have a few possible solutions:
Good luck! Best regards, Bruno
__________________
|
|
January 30, 2013, 03:56 |
|
#5 |
New Member
---
Join Date: Jan 2013
Posts: 17
Rep Power: 13 |
Pretty interesting Bruno thank you very much!
Actually yesterday I was writing a python program to stop the simulation (seting and endTime after writing data and the continue with and endTime = endTime + writeInterval) I'll try the FunkysetFields, I think It'll meet my purpose. Moreover is it simple to change boundary conditions (considering that I stop the simulation and I'm able to write the time folders). Imagine I want to set a PID for the inlet and change the flow depending on the results I get. Thanks, Sergi |
|
January 30, 2013, 17:28 |
|
#6 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,981
Blog Entries: 45
Rep Power: 128 |
Hi Sergi,
implementing a PID... nice! Although I'm not familiar with any such implementation on OpenFOAM... I think swak4Foam has some functionality for using Python controls directly in OpenFOAM... then again, it might need pythonFlu for the ability to control boundary conditions directly with Python code. You'll have to research more on this topic. You can find references to pythonFlu in swak4Foam's wiki page. Either way, it should be possible to do what you want, but coding is definitely involved! Either in Python or C++. Best regards, Bruno
__________________
|
|
Tags |
cell, element, freesurface, mesh, polymesh |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
[ICEM] Overlapping elements | Anorky | ANSYS Meshing & Geometry | 6 | November 6, 2017 09:22 |
[ICEM] Mesh along axis for a wedge gives elements with extreme aspect ratio | sangrampp | ANSYS Meshing & Geometry | 0 | November 17, 2010 07:57 |
[ICEM] Negative Volumes & 7 node HEXA elements | mr_stoked | ANSYS Meshing & Geometry | 1 | September 21, 2010 13:45 |
[Gmsh] Import problem | ARC | OpenFOAM Meshing & Mesh Conversion | 0 | February 27, 2010 11:56 |
Results saving in CFD | hawk | Main CFD Forum | 16 | July 21, 2005 21:51 |