|
[Sponsors] |
April 11, 2011, 11:18 |
|
#21 |
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 17 |
Hi V.
where i have to put the reconCentral files? /src/finiteVolume/....? Thanks andrea |
|
April 11, 2011, 13:10 |
|
#22 |
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 17 |
Hi,
another question... both schemes (reconCentral and pointLinear) are interpolation scheme, if i'm not wrong. So if i use one of these for the gradient Should I use the same for interpolation? like gradSchemes ->Gauss pointLinear; interpolationSchemes ->pointLinear; or, with reconCentral: gradSchemes -> Gauss reconCentral; interpolation -> reconCentral; Thanks again andrea |
|
April 11, 2011, 13:10 |
|
#23 | |
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 20 |
Quote:
If my memory does not fail, you can compile it dynamically by typing wmake libso inside the source folder (wherever you choose to place it). After that you have to set properly your fvSchemes dictionary (search for some examples in the forum) and include in your controlDict something like this: libs ( "libOpenFOAM.so" "reconCentral.so" ) Hope this helps V. |
||
April 11, 2011, 13:23 |
|
#24 | |
Senior Member
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 368
Rep Power: 20 |
Quote:
divSchemes { div(phi,U) Gauss reconCentral cellLimited leastSquares 1.0; //... } interpolationSchemes { default linear; interpolate(HbyA) linear; interpolate(U) reconCentral phi cellLimited leastSquares 1.0; } Hope this helps too V. |
||
April 12, 2011, 04:27 |
|
#25 |
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 17 |
Thanks a lot V., you have been very helpfull!!
andrea |
|
April 12, 2011, 12:58 |
|
#26 |
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 17 |
Hi all, just to let you know...
I finally was able to find a good solution with the unstructured mesh. The best choice for my cases was the pointLinear scheme suggested by Markus. The solution now is a (quite) perfect circle and the behavior during the simulation is very similar to the structured case. The drop is not perfectly still in the middle of the domain but is slightly moving in a sort of circle (probably due to the parasite current) and the solution varies slightly around the analytical one, but for the moment i'm quite happy. Hope this will help someone else! andrea |
|
April 14, 2011, 11:01 |
|
#27 |
Senior Member
Markus Rehm
Join Date: Mar 2009
Location: Erlangen (Germany)
Posts: 184
Rep Power: 17 |
Hello Andrea,
could you post your fvSchemes for the best case? Thank you. |
|
April 14, 2011, 12:35 |
|
#28 |
Senior Member
Andrea Ferrari
Join Date: Dec 2010
Posts: 319
Rep Power: 17 |
Hi Markus,
i tried different setting (always using the pointLinear schemes) and i found good results with almost all. here two of these: 1) Code:
ddtSchemes { default Euler; } gradSchemes { default cellMDLimited Gauss pointLinear 0.333; } divSchemes { div(rho*phi,U) Gauss pointLinear 0.333; div(phi,alpha) Gauss vanLeer; div(phirb,alpha) Gauss interfaceCompression; } laplacianSchemes { default Gauss linear limited 0.333; } interpolationSchemes { default pointLinear; } snGradSchemes { default limited 0.333; } fluxRequired { default no; p_rgh; pcorr; alpha1; } Code:
ddtSchemes { default Euler; } gradSchemes { default cellMDLimited Gauss pointLinear 0.333; } divSchemes { div(rho*phi,U) Gauss limitedLinearV 0.333; div(phi,alpha) Gauss vanLeer; div(phirb,alpha) Gauss interfaceCompression; } laplacianSchemes { default Gauss linear limited 0.333; } interpolationSchemes { default pointLinear; } snGradSchemes { default limited 0.333; } fluxRequired { default no; p_rgh; pcorr; alpha1; } hope this help andrea |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|