OpenFOAM Ship Hydrodynamics
I have found OpenFOAM model files of a wigley hull example with an experimental validation report. The model was programmed to be running in version 1.7.x. To be able to make it run in the 2.2.0 version, I had to make the following changes.
In the "fvSchemes" file, "divSchemes" section: Old code: divSchemes { div(rho*phi,U) Gauss linearUpwind cellLimited Gauss linear 1.0; div(phi,alpha) Gauss vanLeer; div(phirb,alpha) Gauss interfaceCompression; div(phi,k) Gauss upwind; div(phi,omega) Gauss upwind; div(phi,R) Gauss upwind; div(R) Gauss linear; div(phi,nuTilda) Gauss upwind; div((nuEff*dev(grad(U).T()))) Gauss linear; } New code: divSchemes { div(rho*phi,U) Gauss linearUpwindV grad(U); div(phi,alpha) Gauss vanLeer; div(phirb,alpha) Gauss interfaceCompression; div(phi,k) Gauss upwind; div(phi,omega) Gauss upwind; div(phi,R) Gauss upwind; div(R) Gauss linear; div(phi,nuTilda) Gauss upwind; div((muEff*dev(T(grad(U))))) Gauss linear; div((nuEff*dev(T(grad(U))))) Gauss linear; } In the "fvSolution" file, the following section is added: "(U|k|omega)Final" { $U; tolerance 1e-06; relTol 0; } In the "nut" file, "boundary" section: Old code hull { type nutWallFunction; value uniform 0; } New code: hull { type nutkWallFunction; value uniform 0; } Now the file is running perfectly on version 2.2.0. But the drag forces I get from it does not match with the values presented in the report. The values I read is 20% over. Do you think the changes I made can effect the the drag output of the model? Thanks, |
Hi, kilroy, I hope to get an OpenFOAM model files of a wigley hull example. Do you know where to get it? Thanks
|
All times are GMT -4. The time now is 22:29. |