Hi Fabien, I looked over your problem and I haven't been able to debug it. I tried various combinations of constraints and I'm not sure why it's becoming nonphysical (dT >1 or dT < -1). It would require getting into the code and looking at the numbers in the offending cells.
One thing that worries me is that the equations don't have pressure or the continuity equation, which seems strange. I would be rather suspicious of this set of equations because there is no strong coupling between the x and z momentum and also because there they are pure convection without diffusion, which requires special schemes to solve correctly and which FiPy doesn't handle very well. Cheers, Daniel On Thu, Sep 20, 2018 at 8:09 AM fgendr01 <fabien.gend...@univ-lr.fr> wrote: > > Hi Jon, Thibaut and the community ! > Thanks for your help. > I change my parameters back and my problem seems to be stable now. > But I think, there’s still a problem. > On my simulation, there is a continuity problem at the top left and bottom > right. I do not understand why a bubble of average temperature appears! It’s > a bit strange ! > I tested with several boundary conditions but it does not change anything. > Yet I impose a temperature of 1 to the left and -1 to the right and adiabatic > walls at the top and bottom. Do you have an idea ? > Here is some screenshots : -- Daniel Wheeler _______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
