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

moveDynamicMesh

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 10, 2010, 10:40
Default moveDynamicMesh
  #1
Member
 
Sylvain Aguinaga
Join Date: Feb 2010
Posts: 32
Rep Power: 6
Sylvain is on a distinguished road
Hi,

I just took a look on the SnakeRiverCanyon tutorial for moveDynamicMesh, and it looks to me a very interesting way to mesh large topological areas for atmospheric simulations.

In the triSurface folder one can see the .stl file which is actually used but also an .XYZ file. I suppose it means that we can also use a XYZ file instead of a stl file. But I didn't succeed in setting the .0/pointDisplacement file in this way.

Anyone has an idea???

Thank you in advance!!!

Sylvain
Sylvain is offline   Reply With Quote

Old   September 15, 2010, 18:56
Default Test MoveDynamicMesh
  #2
Member
 
ubald's Avatar
 
Nicolas Lussier Clément
Join Date: Apr 2009
Location: Montréal, Qc, Canada
Posts: 43
Rep Power: 7
ubald is on a distinguished road
I don't know for the xyz file but in the snakeRiverCanon it look for a stl file: AcrossRiver.stl
in the 0/ bundary

for my part I'm trying to understand how it work and what are the limit of this option for a mesh on complex terrain !

I tried a 2d case of a bump and realize that the transportProperties nu play a great role in this solver at least with displacementSBRStress that I try to understand better !

There is also one thing I like to do with moveDynamicMesh. Its to set a boundary condition that would set the distance equal on the stl boundary. Ultimately I would like to ave the distance of the point on the stl depend on the curvature of this stl.

But No such boundary condition seam to exist on this solver and for me right now this seam a huge task for me a beginner !

Please Look at the case and the png and if you know how to make it better you are more then welcome

Regards
Attached Images
File Type: jpg Capture.jpg (54.8 KB, 167 views)
Attached Files
File Type: gz MoveDynamicMeshTestBump.tar.gz (7.8 KB, 36 views)
ubald is offline   Reply With Quote

Old   September 23, 2010, 10:07
Default
  #3
Member
 
ubald's Avatar
 
Nicolas Lussier Clément
Join Date: Apr 2009
Location: Montréal, Qc, Canada
Posts: 43
Rep Power: 7
ubald is on a distinguished road
I change parameters, use diffusivity exponential 2 inverseDistance 1(minZ), 2 blocks and I got better results I think !

But it still could be improve at the boundary and on top. On the ground the BC dos not allow the nod to move to relax the mesh close to the ground. Once the point is on the stl it stay right there. Witch it is a problem I think! If there is a way to do it I don’t know it for now !
Attached Images
File Type: jpg moveDynamicMeshTestExp.jpg (81.0 KB, 100 views)
File Type: jpg moveDynamicMeshTestExpZoom.jpg (74.6 KB, 98 views)
ubald is offline   Reply With Quote

Old   September 21, 2011, 11:08
Default
  #4
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 447
Rep Power: 8
samiam1000 is on a distinguished road
Hi guys,

I am a step behind you and I decided to write to you.

I can not find a guide to use this tutorial. But I really wish to be able to run it, since I need to perform a wind simulation over a terrain in a certain area of Europe.

What I want to do is to download a DTED map and then perform a simulation.

The point is that I am pretty new, here and I need help.

Could you explain how to solve and post-process this SnakeRiverCanyon tutorial?

Thanks a lot for help,

Samuele
samiam1000 is offline   Reply With Quote

Old   September 21, 2011, 15:45
Default
  #5
Member
 
ubald's Avatar
 
Nicolas Lussier Clément
Join Date: Apr 2009
Location: Montréal, Qc, Canada
Posts: 43
Rep Power: 7
ubald is on a distinguished road
I can but you need to tel me where you are stok!!

blockMesh
moveDynamicMesh

set the case:
move the mesh you just did to the new case ...

I do wind simulation to. Just to know.

Did you change the wall function ?
Do you use fixed gradient or fixed value on top ?

Regards
ubald is offline   Reply With Quote

Old   September 21, 2011, 17:47
Default
  #6
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 447
Rep Power: 8
samiam1000 is on a distinguished road
Dear Ubald,

thanks a lot for answering.

Well, the 1st thing that I do is to give the 2 commands that you mentioned:
- blockMesh
- moveDynamicMesh (or snappyHexMesh)

Now I should get a mesh (in which folder?) and my question is: which is a good case to perform a wind simulation? I mean: in which case I have to copy my new mesh?

Furthermore, could you share with my your case?

Actually, so far, I have not changed any wall function.

If you agree, we can share a folder on dropbox (samuele.zampini@gmail.com)

Thanks for helping,

Samuele
samiam1000 is offline   Reply With Quote

Old   September 21, 2011, 23:59
Default
  #7
Member
 
ubald's Avatar
 
Nicolas Lussier Clément
Join Date: Apr 2009
Location: Montréal, Qc, Canada
Posts: 43
Rep Power: 7
ubald is on a distinguished road
the final mesh will be in the final folder exemple "1123/" ther you'll find polyMesh folder

You put the content of the final foalder in constant foalder after backing up constant/polyMesh/

mv constant/polyMesh/ constant/polyMeshOrig ...

I will look how to add you to dropbox

regards
ubald is offline   Reply With Quote

Old   September 22, 2011, 04:18
Default
  #8
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 447
Rep Power: 8
samiam1000 is on a distinguished road
Dear Nicolas,

first of all, thanks for answering.

Well, to add me in dropbox is very simple: just upload the case folder in dropbox and then share it with me.
I will receive a notification and everything will be ok!

Also, if I check the folfer 25/polyMesh (being 25 my last folder!), I have only a file inside. It is "point".
On the other hand, if I check the folder 0/polyMesh, I find more files in it.

Is that correct?

Furthermore, do you know why if I use snappyHexMesh instead of moveDynamicMesh I get the following error:

Code:
sammy@nash:~/OpenFOAM/sammy-2.0.1/run/tutorials/mesh/moveDynamicMesh/SnakeRiverCanyon$ snappyHexMesh 
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.0.1-51f1de99a4bc
Exec   : snappyHexMesh
Date   : Sep 22 2011
Time   : 10:17:54
Host   : nash
PID    : 3525
Case   : /home/sammy/OpenFOAM/sammy-2.0.1/run/tutorials/mesh/moveDynamicMesh/SnakeRiverCanyon
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Read mesh in = 0.05 s


--> FOAM FATAL IO ERROR: 
cannot open file

file: /home/sammy/OpenFOAM/sammy-2.0.1/run/tutorials/mesh/moveDynamicMesh/SnakeRiverCanyon/system/snappyHexMeshDict at line 0.

    From function regIOobject::readStream()
    in file db/regIOobject/regIOobjectRead.C at line 87.

FOAM exiting
samiam1000 is offline   Reply With Quote

Old   September 22, 2011, 10:14
Default
  #9
Member
 
ubald's Avatar
 
Nicolas Lussier Clément
Join Date: Apr 2009
Location: Montréal, Qc, Canada
Posts: 43
Rep Power: 7
ubald is on a distinguished road
to run snappyHexmesh you need a dictionary as you do for moveDynamicMesh that is why you get this error there is no dictionary for snappy in SnakeRiverCanyon !

You can use both on the same case but then you do it one after the other !!

I put an exemple on the dropbox this one use snappy end moveDynamicMesh maybe I should right a script to help you but at lest its an other exemple !!


From what I understand, OpenFOAM write in the solution folder only the thing that change from original so it would be normal to ave only points file in your case.

Each tool in OpenFOAM is different you might want to learn one at the time !

I also find an other tool to split a subset of cells in the direction you want not necessary in the three direction at the same time but that is an other story.

ps: please dont unzip the foalder in the drupbox !

Regards
Attached Images
File Type: jpg meshExemple.jpg (96.5 KB, 105 views)
ubald is offline   Reply With Quote

Old   September 22, 2011, 11:22
Default
  #10
Senior Member
 
Samuele Z
Join Date: Oct 2009
Location: Mozzate - Co - Italy
Posts: 447
Rep Power: 8
samiam1000 is on a distinguished road
Thanks a lot.

I'll work on it, soon and I'll let you know.

Keep in touch,

Samuele
samiam1000 is offline   Reply With Quote

Old   September 22, 2011, 14:39
Default
  #11
Member
 
ubald's Avatar
 
Nicolas Lussier Clément
Join Date: Apr 2009
Location: Montréal, Qc, Canada
Posts: 43
Rep Power: 7
ubald is on a distinguished road
welcome
ubald is offline   Reply With Quote

Old   October 4, 2011, 23:12
Default convert .xyz to .stl
  #12
New Member
 
Olivier Benichou
Join Date: Sep 2010
Location: australia but want to move back to europe
Posts: 8
Rep Power: 5
olivier78 is on a distinguished road
Hello openFoamer,

I just had a look at the SnakeRiverCanyon tutorial, it is an interesting tutorial and I agree that the SBR solver is using the .stl surface. However a .xyz file is also present in the tutorial directory.
my question is: how can the .xyz file be converted into and .stl surface ?

Olivier
olivier78 is offline   Reply With Quote

Old   October 5, 2011, 00:58
Default
  #13
Member
 
Sylvain Aguinaga
Join Date: Feb 2010
Posts: 32
Rep Power: 6
Sylvain is on a distinguished road
Hi Olivier,

You can convert un xyz to a STL using paraview:

1) change the xyz extension to a csv
2) open the csv in paraview using the csv reader, choose 1 column for each coordinate
3) use the "TableToPoints' filter to obtain an array of points.
4)use the delaunay tool to "map" a suface from the point
5) save the data, you 'll be able to save it as an stl

TADAAA!!! paraview is your friend. Of course if you have a CAD software as Rhinoceros, it's even easier


Sincerely

Sylvain
Sylvain is offline   Reply With Quote

Old   November 10, 2011, 23:22
Default
  #14
Member
 
Join Date: Sep 2011
Posts: 45
Rep Power: 4
ic3wall is on a distinguished road
Quote:
Originally Posted by Sylvain View Post
Hi Olivier,

You can convert un xyz to a STL using paraview:

1) change the xyz extension to a csv
2) open the csv in paraview using the csv reader, choose 1 column for each coordinate
3) use the "TableToPoints' filter to obtain an array of points.
4)use the delaunay tool to "map" a suface from the point
5) save the data, you 'll be able to save it as an stl

TADAAA!!! paraview is your friend. Of course if you have a CAD software as Rhinoceros, it's even easier


Sincerely

Sylvain
I just tried your procedure and got the following warning:

Delaunay3D Unable to factor linear system


Then I create the surface (Extract Surface) but it's impossible to see the results. I save it as .stl and can't see the content of the file.

EDIT: Use Global Mapper instead http://www.globalmapper.com/ !

Last edited by ic3wall; November 11, 2011 at 13:03.
ic3wall 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
moveDynamicMesh broken ? podallaire OpenFOAM Bugs 2 November 11, 2009 13:44


All times are GMT -4. The time now is 10:54.