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

OF meshless

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

Like Tree6Likes
  • 3 Post By elvis
  • 1 Post By wyldckat
  • 1 Post By wyldckat
  • 1 Post By elvis

Reply
 
LinkBack Thread Tools Display Modes
Old   September 14, 2013, 10:43
Default OF meshless
  #1
Member
 
Join Date: Jun 2011
Posts: 51
Rep Power: 6
cfdivan is on a distinguished road
Dear all,

Considering how time consuming is the mesh process in any of our simulations, is there anybody aware of the implementation of "meshless methodoly" in OF? I am expecting with this approach reduce significantly the hours spent in mesh debug and improve the solution performance by reducing the non-convergence solutions. There are any validation tests together with the code ready to run?

It would be quite interesting try out some performance tests of OF with this mesh methodology.

Regards,
cfdivan is offline   Reply With Quote

Old   September 15, 2013, 06:26
Default
  #2
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,253
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings cfdivan,

It took me a while to remember about this... but I remembered just now. There is a non-official implementation (aka it's not in the official and public OpenFOAM version) of the "Immersed Boundary Method" developed by Wikki: http://wikki.gridcore.se/application...BoundaryMethod

AFAIK, always keep in mind that there is no such thing as a silver bullet in CFD. One meshing strategy might be excellent for one thing, but terrible for another. And the "Immersed Boundary Method" and/or meshless are always implied meshing strategies, namely that even though you do not have to bother to generate a mesh, a mesh does still exist, since the geometry is still analysed in a discrete way, not in a continuous way!

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   September 16, 2013, 07:24
Default
  #3
Senior Member
 
Elvis
Join Date: Mar 2009
Location: Sindelfingen, Germany
Posts: 577
Blog Entries: 5
Rep Power: 13
elvis is on a distinguished road
Hello,

just my 2 cent

you can dig into Immersed Boundary Meth. (IBM) =>Z. Tuković (source code 4 OF-extend) by studing OF-workshop material (OFW7_2012_Darmstadt)
http://sourceforge.net/projects/open...darySlides.tgz
http://sourceforge.net/projects/open...lmMaterial.tgz

I remember two guys that implemented Smoothed Particle Hydrodynamics (SPH)-method with OpenFOAM one presentation at OFW6 @ Penn State
http://sourceforge.net/projects/open...mar_slides.pdf
the second guy Basti Werdelmann was giving a presentation at OF-Stammtisch @ Karlsruhe
http://www.rz.uni-karlsruhe.de/~hj16...Stammtisch.pdf.
You should contact Hendrik Rusche (Wikki) or Ivana Buntic-Ogor (on the OF-Stammtisch List) if you want the Stammtisch-presentation they might have it

But I am not aware of any SPH-code for OF out in the wilderness to the public.

But these examples show it can be done with enough OF-expertise, s.th I do not have
:-(
wyldckat, cfdivan and ali_atrian like this.
elvis is offline   Reply With Quote

Old   September 17, 2013, 17:28
Default
  #4
Member
 
Join Date: Jun 2011
Posts: 51
Rep Power: 6
cfdivan is on a distinguished road
Hi all,
Thank you very much for both posts...they are quite usefull.

Bruno,
Yes, even for "meshless" approach the domain/flow is solved in discrete way but, I meant the immersed BL approach . It is tipically quite easy do a pretty decent surface meshes. Even with the decent surface mesh, problems in "BL zone" (the prism layer extrusion process) are quite often. Those areas are usually the origin for some solution divergence. I knew that wikki already tried that approach. I guess that karalit also implemented that approach. I was just wondering if this approach was already implemented by someone in the forum who can share the code/expereciences.

"AFAIK, always keep in mind that there is no such thing as a silver bullet in CFD. One meshing strategy might be excellent for one thing, but terrible for another". Do you know what kind of applications that immersed BL approach could be terrible?

Elvis,
Thank you very much for all the links. I'll try to dig in a bit more and share the findings/ideas here.

Regards,
cfdivan is offline   Reply With Quote

Old   September 17, 2013, 18:07
Default
  #5
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,253
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Hi cfdivan,

Quote:
Originally Posted by cfdivan View Post
"AFAIK, always keep in mind that there is no such thing as a silver bullet in CFD. One meshing strategy might be excellent for one thing, but terrible for another". Do you know what kind of applications that immersed BL approach could be terrible?
Sorry, I vaguely remember reading somewhere about one or two situations in which it didn't work very well, but I can't remember where I read about it.
Which is why I ended writing about what has been found several times in CFD for many other situations: no matter how many new methods we find for simulating or meshing, they always have some flaw that leads us to go back to the good old standard procedures/methods/models.

Best regards,
Bruno
cfdivan likes this.
wyldckat is offline   Reply With Quote

Old   September 30, 2013, 03:13
Default
  #6
Senior Member
 
Elvis
Join Date: Mar 2009
Location: Sindelfingen, Germany
Posts: 577
Blog Entries: 5
Rep Power: 13
elvis is on a distinguished road
Hi,

I saw http://www.wire.tu-bs.de/workshops/nofun13/d_nofun.html
=>"LBM mesh generation in OpenFoam" a presentation by Andrea Pasquali (TU Braunschweig)
so some kind of Lattice Bolzmann Method might be possible to be handled with as well
elvis is offline   Reply With Quote

Old   October 2, 2013, 17:25
Default
  #7
Member
 
Join Date: Jun 2011
Posts: 51
Rep Power: 6
cfdivan is on a distinguished road
Hi Elvis,

Thk u very much for your link. "LBM mesh generation", means a meshless process? What is the main difference between "LBM mesh generation" and "Immersed Boundary Method"? I checked on the website but I did not find any additional information about the topic. Do you know where I can find more documentation about the topic LBM mesh generation in OpenFoam and/or comparision with classic methods or Immersed Boundary Method?

Many thank
cfdivan is offline   Reply With Quote

Old   October 5, 2013, 01:51
Default
  #8
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,253
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings to all!

@cfdivan:
Quote:
Originally Posted by cfdivan View Post
"LBM mesh generation", means a meshless process?
Nope. LBM is this: http://en.wikipedia.org/wiki/Lattice_Boltzmann_methods
In a nutshell, based on what I studied a few years ago:
  1. It doesn't use FVM (Finite Volume Method), while both OpenFOAM and IBM do use it.
  2. You still have a mesh, usually using uniformly cube-shaped cells.
  3. The data is passed between the cells as if using a particle based mechanism, i.e., one particle can only move from one cell to another. It's not exactly this, but its the mnemonic I use to understand this model.
  4. The advantage of this over FVM CFD is that it's easily parallelizable, namely for using GPGPU hardware.
I think LBM has improved considerably in the last 3-4 years, but it probably still has several limitations.

Best regards,
Bruno
cfdivan likes this.
wyldckat is offline   Reply With Quote

Old   October 5, 2013, 05:43
Default
  #9
Senior Member
 
Elvis
Join Date: Mar 2009
Location: Sindelfingen, Germany
Posts: 577
Blog Entries: 5
Rep Power: 13
elvis is on a distinguished road
Hi,

basicly I wanted to say that OF is a Computational Continuum Mechanics (CCM) toolkit
which can be used for different methods like FVM and FEM which are already implemented.
But other people used it for different methods like SPH
elvis is offline   Reply With Quote

Old   October 5, 2013, 07:12
Default
  #10
Member
 
Join Date: Jun 2011
Posts: 51
Rep Power: 6
cfdivan is on a distinguished road
Hi wyldckat,
Thk u very much for your link. Pretty usefull for a brief review.

Regarding FVM - I already tried to investigate this point some time ago and I found that they likely use FD...is it true? I never was able to confirm that

With FD we have several advantages when we are working with structured meshes....it is very simple and effective and quite easy to obtain higher-order schemes (High-Mach numbers applications for example).

"The data is passed between the cells as if using a particle based mechanism" - This really sounds a restriction for high mach-number flows where you need higher-order schemes .
Considering that you have different sizes of the cells in your domain and the equations are solved in "Lattice world", in one time step you can have one particle go to the neighbor cell while in the other zone of your domain the particle can travel 2 or more cells. Based on this assumption, we will not be able to track the particle distribution. They need to work around with field update, is it right?

The paralelization is in fact an significant advantage of this code. The algorithm is quite simple and only the stream step needs exchange of information with neighbors which does not happen in FVM.

It looks that the efficieny of this code comes from the implementation of simple algorithms where non higher-order schemes are needed. Though, when it comes to high-mach numbers and coupling of density to temperature variations, higher-order schemes are nedded and the simplycity/efficiency of the code decreases.

There is also a point regarding the full unsteadiness of this approach. When compared with DES-SA of NS for example, it might be more fast/easy to obtain a unsteady converged solution. Nevertheless there are an wide range of engineering applications where the RANS/steady state provide us enough accuracy for very low prices in terms of CPh which migh not be possible with LBM. It looks like LBM fits better for applications where the unsteadiness is the king.

In any case, what looks really interesting for me is the "meshless concept" that is very painfull process on my simulations that I am expecting that we do not have neither in IBM nor LBM. From time to time I spend a lot of hours to check why there is no convergence on simulation, ending up play around with disc schemes and mesh.

Regards,
cfdivan is offline   Reply With Quote

Old   October 5, 2013, 07:22
Default
  #11
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,253
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Hi cfdivan,

As far as I can remember, LBM works in a similar way to the Courant number restriction that we're taught in the very first tutorial OpenFOAM User Guide: http://www.openfoam.org/docs/user/cavity.php
In essence, the size of each cube must respect the flow speed limit, otherwise it simply cannot present properly your problem. In addition, the cube sizes also had to take into account the properties of the fluid, mainly the viscosity.
In other words, LBM is almost a "Similitude" modelling strategy: http://en.wikipedia.org/wiki/Similitude_%28model%29

Therefore, with any luck, for supersonic flow, you only need several billion LBM cells, working on a cluster with 8 GPUs per node...


As for "FD": which "FD" are you talking about?

Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   October 9, 2013, 20:10
Default
  #12
Member
 
Join Date: Jun 2011
Posts: 51
Rep Power: 6
cfdivan is on a distinguished road
Hi Bruno,

Sorry for the delay on reply
FD, I meant Finite diffrences.

I also agree that the LBM might works in similiar way of Courant number. To me it sounds that the scalling of the flow properties (like velocity for example) to the "Lattice world" take in consideration the cell size and the speed of the sound, is it right? The "Courant number" in LBM is ajusted at each cell in order to keep the tracking of the particles at each timestep.

Regards,
cfdivan is offline   Reply With Quote

Old   October 9, 2013, 21:47
Default
  #13
Senior Member
 
Julien de Charentenay
Join Date: Jun 2009
Location: Australia
Posts: 229
Rep Power: 9
julien.decharentenay is on a distinguished road
Send a message via Skype™ to julien.decharentenay
Hi cfdivan/Bruno,

Great discussion. The openFoam solvers (simpleFoam, etc) are based on FVM discretisation of the NS equations.

Wouldn't it be great to have access to various methods (FVM, IBM, LBM, SPH) using a unique and compatible set of input data rather than juggling between solvers, pre-processors, etc.... This would make it easier to use the right tool for the right job.
__________________
---
Julien de Charentenay
julien.decharentenay is offline   Reply With Quote

Old   November 9, 2014, 07:49
Default "Basti Werdelmann" presentation at OF-Stammtisch
  #14
New Member
 
Alireza Atrian
Join Date: May 2014
Posts: 26
Rep Power: 3
ali_atrian is on a distinguished road
Quote:
Originally Posted by elvis View Post
the second guy Basti Werdelmann was giving a presentation at OF-Stammtisch @ Karlsruhe
http://www.rz.uni-karlsruhe.de/~hj16...Stammtisch.pdf
Hello elvis
could you please give me a more accurate address of "Basti Werdelmann"?
I'm going to know how to implement SPH in openFoam so I googled the name for hours but noo cunclusion!
is the presentation on 6th openfoam workshop (13-16 June 2011)? or ...?
the link you have inserted does not work any more!
do you know any other presentation, work, coding, tutorials on the issue?
Best Regards
ali_atrian is offline   Reply With Quote

Old   November 10, 2014, 03:33
Default
  #15
Senior Member
 
Elvis
Join Date: Mar 2009
Location: Sindelfingen, Germany
Posts: 577
Blog Entries: 5
Rep Power: 13
elvis is on a distinguished road
Hi Alireza,

the implementation of the SPH presented at Karlsruhe OF-Stammtisch is a different one to that presented at Penn State OF-workshop, but both are closed source to my knowledge.

You should find his address in https://wiki.manchester.ac.uk/spheri..._Abstracts.pdf
Hope you are not a headhunter ;-)
ali_atrian likes this.
elvis 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
Xflow Meshless aproach ?? leo_NM Main CFD Forum 1 June 23, 2013 19:38
meshless Hamou SADAT Main CFD Forum 4 January 3, 2008 17:48
Looking for two Papers on Meshless Gerrit Groot Main CFD Forum 4 September 20, 2007 09:21
Meshless FD method over FEM? Shyam Main CFD Forum 2 May 5, 2004 02:35
Meshless CFD Des Aubery Main CFD Forum 1 August 17, 2003 08:42


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