Value of gravx

[Madhurjya814]

Incompact3d/src/transeq.f90

Line 776 in 519d615

dux1(i, j, k, 1) = dux1(i, j, k, 1) + peculiar_density1(i, j, k) * richardson * gravx

In this line the term richardson already contains gravity value g=9.81 m/s2. So the gravx value should be put as -1 if my flow is upward along x direction and buoyancy is acting in the negative x direction? We should not use -9.81 as gravx value, as richardson already contains g=9.81? Is it correct?

[mathrack]

Hello,

You can check the 2D cavity for an example of a buoyant flow : https://github.com/xcompact3d/Incompact3d/blob/master/examples/Cavity/input_2D_LR_Re14084.i3d

The exact definition of dimensionless numbers (such as richardson) is case-specific. In the example above, the gravity vector has a unit length yes.

Cheers,
Mathrack