How to stop the oversetMesh from leaving the main domain
Hello!
I am trying to simulate something like a stone falling down in a vessel. I use OF v2312, oversetMesh and sixDoFRigidBodyMotion without any constraints and restraints (because I want the solid object to move according to the flow). I include the libs (overset fvMotionSolvers) in the controlDict file. But the problem is that the object falls down and leaves the main domain. How can I make it "feel" there is a bottom? Maybe I should include some restraints, but I don't know what exactly. I have another issue, which is related to the complex geometry of the main vessel. I use the adjustTimeStep and maxCo options, but the result is a very huge Courant Number and very small deltaT - smaller than 1e-20 or 1e-30. Is there any possibility to fix this besides making finer mesh, which would probably lead to very small deltaT again? I would highly appreciate any advice! Best regards! |
Hey!
Without going as far as modelling proper impact, there is a constraint that acts as a soft wall: https://www.openfoam.com/news/main-n...sics#soft-wall Maybe it can suit your needs! |
Thank you very much for your reply, Alczem! It really suits my case and helped me a lot! The only problem now is that whenever the object reaches the bottom, this leads to error and the simulation stops. Maybe I did not use proper values for psi and C (first I used 0 for both). Actually, I tried different values, but then it behaved very strange - the object goes far away outside the domain in the side and upper direction for example. Do you have any recommendations for these parameters?
Thanks a lot again and best regards! |
I never used it myself, but you can check out the source code of the constraint. There we can see this:
Code:
{ It also seems to me that the anchor point should be away from the wall you want to avoid, in order to let d become negative. For instance, if your "ground" is at z=0m, maybe have the anchor at z=0.5m and wallNormal should be (0,0,1). I don't know how you first set up the constraint, but that is my interpretation of the code. |
Yes, I will try to calculate some reasonable values. Тhank you very much for your help and advices! :)
|
All times are GMT -4. The time now is 20:08. |