[atelcikti1]

Hi CFD fellows,

Currently i am trying to solve mixing (later laminar combustion) on liquid rocket motor. Species are some hydrocarbons, hcl, o2, oh, h2o so on so forth. I am using rhoReactingfoam for this. i am doing unsteady analysis and slowly increasing inlet mass flow rate (flowRateInletVelocity) for U, zeroGradient for P, pressureInletOutletVelocity for U and waveTransmissive for P at outlet then initializing domain with atmospheric conditions. At starting phase of analysis, solution seems like stable, but after some time, time step sizes reduces too much (ie. to 1e-12), so i observe no change in solution, even some variable matricies are not being solved. I am using maxco 0.5, initial stepsize 1e-7,1e-8 whatever doesnt matter at all. I am observing no recirculation at boundaries, no janaf temperature region warnings (i am limiting the numbers and expand the thermo properties until 20000K). Any idea about this strange behavior?

Thanks in advance.