Forum Replies Created
July 7, 2021 at 8:01 pm in reply to: Poiseuille 2D Center Velocity Error (Forced,bounceBack) #5811
LBM is a second order scheme, i.e. the numerical solution will approach the analytical one up to an numerical error. Best MathiasJuly 7, 2021 at 6:02 pm in reply to: Poiseuille 2D Center Velocity Error (Forced,bounceBack) #5809
Standard bounce back in general work like this with the wall assumed to be located 1/2 off the nodes. In the mentioned example we have show how to realize it. Best MathiasJune 28, 2021 at 10:48 am in reply to: [BlockGeometryStatistics3D] No discreteNormal is found #5756
You may consider taking part in our next spring school, where we explain how to set up a geometry step by step. It is hard to see what your problem here is.. Basically, in OpenLB you set up a geometry by indicator functions which can be given by stl files or geoetry primitives like cylinders, cubes, spheres, etc.
It scales the desity to e.g. physical units while the computation always uses lattice units where who is scaled to rho.
my team helped you to get a basic version running. Finxing stability issues by means of our forum, however, is more that we can do.
I have no idea about that topic? Which equation system do you want to solve?
You can see that in our multi-component examples. Best Mathias
it is basically there ot distigish betwwen bounce back boundaries (0.5) and other (1.0).
Using a pressure boundary means fixing the density which is not a good idea if you are interested in the density evolutions. You should set a convetion boundary condition instead, which is a Neuman condition for density and velocity!
You can do, yes.
Yes, we do have a fre engery model implemented. Have a look at our examples. Best Mathias
If you refer to the OpenLB examples, that is a Free Energy Model and not SC. BEst Mathias
You need to write a new method which starts the reagion growing for any point within an idicator, lets say asphere, which intersections with a particular material number, e.g. 2, to start the algirthm. But keep in mind that you need to program the parallel region growing algorithm first. Within a colaboration we could help you further but that needs some efford.
you can work with intersections of geometry primitives (indicators) like cylinders, cuboids, etc. with the actualy generated voxel mesh by and stl-file or you may extend the region growing apporach (https://www.openlb.net/DoxyGen/html/d3/d2a/classolb_1_1BlockGeometryStructure3D.html#aa1a30cd27881b163d0486556769fcc70) to the superGeometry implementation, which is only available for blockGeometries at the momenent.
MathiasApril 19, 2021 at 7:57 pm in reply to: Visualization of output for domains created with imported STL file #5614
In allof the OpenLB examples a file “geometry” is written as vtk file, which you can be used for it.