[shota]

Simply, call the coupling operator every n-th time step using if(iT%n == 0) or something alike.

[shota]

Hi,

I think when you are new to OpenLB, the best would be to check out some example cases which are already implemented. To match your topic, I would start with those in the examples/advectionDiffusionReaction folder. There you can experiment with few apps, choose one which fits your needs best. I would then recommend to transform the example to solve your problem.

Best wishes,
Shota

[shota]

Dear Hervé Mutelle,

in the squareCavity2d example the wall is a simple straight wall. In this case, you can define your own AnalyticalFunctor to prescribe a function for the temperature depending on the position.
Please refer to the example/advectionDiffusionReaction/convectedPlate3d. There a custom AnalyticalFunctor with the name ConvPlate3D which computes the analytical solution as a function of the spatial coordinates. This is done by the operator() function taking the spatial coordinates via input[] and passes the results via output[].
Using the same technique you can define a function for the wall and pass that functor to initEquilibrium and defineRho for the material numbers of the wall.

Hope this helps.

Best wishes,
Shota

[shota]

Dear Mutelle,

sorry for the late response.
I looked into the implementation of the case and it seems that that example uses a deprecated way of assigning material numbers. I would refer to the laminar/cavity2d example regarding this issue. We will work on updating the deprecated case in the next release.
If you need more specific assist, please let us know.

Best wishes,
Shota Ito

[shota]

Dear participant,

We investigated your concern about the cylinder3d example. However, we were not able to reconstruct your problem. Please provide more input in order to further understand your concerns.

For this purpose, we created an issue regarding your concerns on our release repository on GitLab.
There, it is possible to discuss more in detail about your issue:
https://gitlab.com/openlb/release/-/issues/2

Best wishes,
Shota

[shota]

Dear Danial,

thank you for your notice. We will have a look on your issue and reply on gitlab.

[shota]

Please send us the pictures via mail or upload them and send us a link, then we can have more insight into your issue.

[shota]

Dear Danial,

indeed the “setAdvectionDiffusionTemperatureBoundary” can not be used on curved walls.
You can try the “setBouzidiAdeDirichlet” BC which is a similar scheme to the Bouzidi scheme but for prescribing a fixed concentration on curved walls for ADE. However, this BC is not fully validated yet. The implemented scheme can be found here https://doi.org/10.1016/j.jcp.2012.11.027.

[shota]

Dear Massload,

the boundary conditions are applied depending on being a local or non-local operation (e.g., the interpolatedVelocityboundary is non-local as it requires interpolation from neighboring cells) at different stages during collideAndStream() managed by STAGE (e.g., PreCollide(), PostCollide(), etc.).

When you look into src/boundary/setInterpolatedVelocityBoundary3D.hh, the used dynamics and postprocessors (for the local and non-local operations respectively) are listed depending on the alignment of the boundary. In the definitions of the applied dynamics and post-processors, you can find the actually applied boundary treatments. At this point, I kindly recommend you to refer to our User-Guide Chap. 2 about core concepts of OpenLB.

[Author]

Posts