Can OpenFoam solve this problem?
Hello everybody
I am trying to solve the next problem that I have been struggling with trying to implement it by myself. I was wondering if I could use OpenFoam instead and how long it could take me in order to figure out how to solve it (I'm pretty good with C++ and my knowledge on turbulence models is decent at this point) I am asking this question because I have a time constraint (1 month) and I don't know if I will have time to switch and solve the problem with OpenFoam or try to debug my implementation. So basically I'm asking for two simple questions: Can OpenFoam solve this problem? how long would it take me given my current knowledge on turbulence models and C++? The problem is this described in this paper (http://onlinelibrary.wiley.com/doi/1...20126/abstract) I will describe it here as well, but you can find the equations there. I'm solving a turbidity current vertical profile. From the 2D RANS equations, I neglect the convective terms because I assume the flow has developed and does not change in the longitudinal direction. With that assumption and the conservation equation, I am left with just the streamwise velocity component's equation. This mean flow equation has a source term which is the sediment concentration. The concentration follows a conservation equation. As a turbulence model, I'm using k-epsilon. One thing different from the paper mentioned above is that in my case there is no roof, the turbidity current is in the lower part of a larger domain that it is mostly at rest. The tricky part of the implementation is the wall function. Special care needs to be taken because of the buoyant forces produced by the concentration. Can OpenFoam take this into account? The common log-law is not accurate since it ignores buoyant forces. So, can OpenFoam solve this? If so, how long will it take me? Please let me know if you might need more details. Thanks in advance. Edit: High Reynolds number is assumed. Miguel |
Hi Miguel, did you solve your problem ?
If there's settling within your turbidity current you can use settlingFoam. I use it and reproduce well experiment found in literature except that the advective term in the AlphaEqn.H are missing. It could explain the early the stop of my turbidity current. I try to implement the advectiveFvPatchField.H within this solver and see how it does affect the motion. Xavier |
All times are GMT -4. The time now is 13:50. |