[nimasam]

hi dear friends
density is a variable of p,T, X
rho = f(X,T,p)
p = pressure
T = temperature
X = conversion

i changed the pEqn like below to consider compressiblity
fvScalarMatrix pEqn
( fvm::laplacian(rhof*rUAf, p) == fvc::dd(rho) + fvc::div(rho,phi)

);

but i enocunter very small time step 10 ^-11 OR less
do u have any clue to increase time step with keeping courant number bounded? is using simple algorithm help to keep courant number bounded?
can i use relaxation for transient solution?