code_saturne User's Forum

Hello, The variation in density is not necessarily linked to the compresibility of the fluids. As you see in your case, the density is only or almost only a function of temperature, sometimes the density is a function of mixture of two gases too, though the Mach number is actually low(even smaller t...

So you just need to use the idilat option to represent the variable density.
Feel free if you have any questions about the initialization

Li

I made a mistake of structure name.
It should be
cs_velocity_pressure_model_t *vp_model
= cs_get_glob_velocity_pressure_model();
then you can change it by
vp_model->idilat = 2;

It seems there were errors of dimension when accessing the velocity array. under fortran if you get the pointer of velocity vel, then the initialisation for the cell iel is vel(1, iel) = u0 vel(2, iel) = v0 vel(3, iel) = w0. Notice now the compressible module is not compatible with combustion module...

you can change te idilat to 2 in cs_user_parameters-base.c by setting void cs_user_parameters(cs_domain_t *domain) { cs_glob_velocity_pressure_param->idilat = 2 ; } In addition, it is possible that density is not postprocessed in visualization output. You can enforce the option by setting void cs_us...

Hello, irovar = 1 doesn't necessarily mean that the density variation is taken into account in the fluids equations, it indicates the density may be a function of a state law. You may find the idilat option in setup.log of your run folder that indicates which dilatation model you used. Could you ple...

Hello, Sorry for the response a little bit late. I tested this case rapidly. It is true that the ignite process is not observed whatever the inlet temperature is set. I got some information that may help you from one researcher who used Libby Williams model (maybe the only one at EDF R&D for now...