[zxyf1]

Hi Foamers,

I have been working on a shock tube case for months, using layer addition and removal scheme as my dynamic mesh scheme in order to simulate the motion of diaphragmless shock tube's plate between driven and driver section. The solver I use is rhoCentralDyMFoam which was copied from 3.0.x and modified to implemented in foam-extend-3.0.1.

The problem I have is when the solver try to add an additional layer of cell in the mesh, the Courrant Number will jump from 0.3 to 10 which will cause the solver crashes. Is there any way to solve this kind of problem?

I have been tried increasing the resoultion of the mesh, however it will occur another error, which indicates the temperature is extremely high:

    From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const
    in file /usr/local/OpenFOAM/foam-extend-3.1/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 72
    Maximum number of iterations exceeded.  Rescue by HJ

If you have any idea please tell me, really appreciate that.
Best!