mathias
Forum Replies Created
-
Posts
-
You need to check yourself. This is too detailed know. If you need further support, please contact me by e-mail.
Best
MathiasIt should work, if you find where it comes from, please, let us know!
Hi steed188,
that is not possible with the actual implementation.
Best
MathiasIt seems to me that you did not compile the code in parallel. You may try 1/ make cleanbuild, 2/ make clean, 3/make. And then try again. Best Mathias
Please have a look at the example cylinder2d, cylinder3d or aorta3d to get an idea of the set-up as well as consider to come to our next spring school in March 2018 in Karlsruhe!
Best
MathiasThese are different Types. Come to the spring school next year and we will have time to explain it in detail. The FluxFunctor should also work for materal numbers or with a greater cylinder including your rextangle plus a material number. Best Mathias
You can check the material numbers in paraview. It looks like the dynamics are not assigned correctly. The best is, if you come to the spring school where we explain it step by step. The forum is not the right place. Best Mathias
You may start with one of the existing OpenLB examples. Then, you set up an indicator of your choice. Then, you use the rename function to rename material numbers which you, finally, use to set dynamics (collision modells) or boundary conditions of your choice. I advice you to go though the examples to see what is possible and register for the next spring school.
Best
MathiasIt all seems correct. Did you check the material numbers for the outflow area?
Changing charL will chance the Re but not the results since charNu and therfore tau stays the same. Best Mathias
The value *uAv should be the actual value of the average velocity at the outflow area. That means you must set it to a desiged value or compute it with the e.g. flux functor every time step.
No, that is not the way. You need to use the define the right value “uAv” by prividing a pointer in the constructor:
virtual void olb::OnLatticeBoundaryCondition2D< T, Lattice >::addConvectionBoundary ( BlockGeometryStructure2D< T > & blockGeometryStructure,
int material,
T omega,
T * uAv = NULL
)Best
MathiasCome to the spring school next year and we will have some time to go through that.
For the convection bc, how did you set the to be imposed velocity?
Dear steed188,
your solotion is not the best in terms of efficiency but I do not see why it should not work. If you find out the problem, please let us know.
Best
Mathias