dynamicRefineMent in compressible thermo cases
4 Attachment(s)
Hello everyone,
I've been doing some experimenting with dynamic refinement in compressible solvers and have found a problem I am not sure how to deal with. I am using the dynamicRefineFvMesh on the aachenBomb which has been modified with a coarser starting geometry and a sprayDyMFoam solver which is a synthesis of rhoPimpleDyMFoam and sprayFoam. For some reason, after solving rhoEqn, the density in the refined cells increases dramatically to the extent that the the solver eventually crashes. I have attached an image of the rho, p, U and C7H16 fields immediately after rhoEqn.H. I confirmed that rhoEqn.H is the cause as I also wrote out the fields just after mapping onto the refined mesh and before rhoEqn and the fields don't suffer the discontinuity. Anyone have any ideas as to what could be causing this? It seems like the ddt(rho) term in the equation is what is causing the problem. I've included the dynamicMeshDict, fvSchemes and fvSolutions files for reference. Thanks, |
Hello mturcios777:
Recently, I am trying to implement dynamic mesh refinement in sprayFoam solver,but I faid. As you have succeded sprayDyMFoam,I was wondering could you tell me in volScalarField YtfSource of YEqn.H file, what should I add, dimensionedScalar("zero", parcels.*****().dimensions(),0.0)??????? great thanks in advance! |
Hi Tjucus,
This post should help you out: http://www.cfd-online.com/Forums/ope...tml#post340938 Good luck! Marco |
Thanks a lot! I will try it!
Regards, |
Hello mturcios 777
Sorry for bothering you again! I want to make some modification in the scalar Urel in the SprayParcel.C,to do this, the injection direction of injector a SprayParcel belongs to is necessary? Is there any function avaliable to get the inj direction in the SprayParcel.C, I also tried to use dictionary to look up the directions in the case files,but it seems that dictionary does not work in the SprayParcel? Could you help me solve this problem? Great thanks in advance! |
All times are GMT -4. The time now is 19:38. |