Converting implicit solver to explicit one
Hello,
I'd like to simulate detonation instabilities. The reactingFoam solver is a good basis for this, however I removed the diffusive terms and now need an explicit formulation of the solver. So, my first approach was to use the fvc:: formulation instead of fvm:: wherever possible: fvm::ddt(rho) + fvc::div(phi) fvm::ddt(rho, U) + fvm::div(phi, U) == rho*g - fvc::grad(p) fvm::ddt(rho, h) + mvConvection->fvmDiv(phi, h) == DpDt fvm::ddt(rho, Yi) + mvConvection->fvmDiv(phi, Yi) == chemistry.RR(i) becomes fvm::ddt(rho) + fvc::div(phi) fvm::ddt(rho, U) + fvc::div(phi, U) == rho*g - fvc::grad(p) fvm::ddt(rho, h) + mvConvection->fvcDiv(phi, h) == DpDt fvm::ddt(rho, Yi) + mvConvection->fvcDiv(phi, Yi) == chemistry.RR(i) Is this really all? Do I now have an explicit solver? So far, everything works. Now I have to handle the pressure-velocity coupling in the PISO loop. Initially it looks like this: fvm::ddt(psi, p) + fvm::div(phid, p) - fvm::laplacian(rho*rUA, p) Here, the solver crashes if I replace the fvm::'s in the respective terms. Can anybody explain the problem, please? The error message: Quote:
|
Same error here trying to convert pisoFoam to explicit.
Any clue? |
All times are GMT -4. The time now is 16:51. |