Dear Stéphane,
Here are some answers (mine
) to your different questions...
1) Yes. But unfortunately this is completely undocumented stuff... apologize! You'll have to look at the
cs_ventil.c file in the
src/base directory of the kernel package. You'll find there a Fortran API to access some C functions (usually to be called in this order):
- define a fan -> defvtl (to be done only once)
- initialize fans -> inivtl (to be done only once)
- get the number of fans -> tstvtl
- compute the source term associated -> tsvvtl
Other functions exists to compute the mass flux -- both inflow and outflow -- accross the fan (
debvtl), to know which cells are "in" a fan (
numvtl). All of this can be done in the
ustsns subroutine for source term management (or maybe also in
usiniv for fans definition and initialization). And for non-french speakers, fan is translated in French by ventilateur (thus the vtl suffix of all the subroutine API).
ps: if needed, I might dig up an example somewhere...
2) I checked how we compute the face warping coefficient and it seems to be an angle degree (0 being the best value == no warping). How's your mesh ? If there are only tetrahedra, then the warping is (of course) zero ; if there are hexaedra, then it is indeed kind of strange to have so low value (except for brick elements of course).
3) Not really... Let's say that the SOLU scheme is a second-order upwind scheme, but that is generally worse than upwind regarding the min/max principle. Centred scheme is 2nd order as well and may be worse than the SOLU scheme regarding the min/max principle (i.e. more dispersive) but is generally far better regarding energy conservation (i.e. less diffusive). I don't really have advice here... and usually use centred scheme for velocity and scalars (plus clipping to tackle the potential overshoot) and upwind (1st order) for turbulence.
4) Basically, k-epsilon model for historical reasons and habits (some could say bad habits... I don't want to start a flamewar here
) with its linear-production version for jet-impacting flows, but SSG model is being used more often now; k-omega SST is also used but you want to avoid natural convection as it is usually bad for this kind of flows; LES of course, mainly the standard Smagorinsky model, whenever frequential aspects are needed (depending on mesh constraints, wavelength of interest for the physical aspects, ...). v2-f is currently being improved in a new version, not yet integrated within the official version, so you can consider it as in a "stabilization" stage. LRR and mixing length models are rarely used.
Regards,
David