compressibleInterFoam diverging Temperature/Velocity problem
Hi to all,
I am trying to run laval Nozzle simulations using compressibleInterFoam. Nozzle has 3 inlets, 2 of them are liquid nitrogen whereas other one gas nitrogen. I tried several boundary condition, however finally always getting number of iterations exceeded based on final Temp became negative whereas the velocity is too high. My boundary conditions are as follows: alphaLiquid: Code:
boundaryField Code:
inletliq1 Code:
boundaryField Code:
inletliq1 Code:
divSchemes Code:
PIMPLE: iteration 1 Is there anyone that can give some suggestion. Thank you in advance. |
First and most important step.
You should run checkMesh. If your nonOrthogonality is above 50 add nonOrthogonalCorrectors. If it is above 60 consider remeshing. The aspect ratio is another big problem for the interfoam solvers from my experience. nOuterCorrectors can be increased and relaxation factors added. However those are all just measures to counteract bad meshes. Please post the checkMesh log. Are you using a turbulence model? Some things that might help, but probably won't. You could try adding this to your grad schemes: Code:
grad(U) cellLimited Gauss linear 1; |
Quote:
Thank you for your suggestions. I am running the case using 2D geometry and turbulence model is not so big deal. first, I just want to fix problem which I mentioned above. (Velocity is extremely increasing almost 15~20 times larger than inlet velocity, whereas Temp is going to negative...) This is checkMesh: Code:
Create time I set also the fvSolution as you suggested: Code:
PIMPLE Any other request? Thank you. |
As an adiditon, I also commented out TEqn.H to see whether it affects or not velocity increasing. However, no solution. still velocity is increasing extremely high and giving same error which i post at #1.
Any suggestion will be appreaciated. thank you |
Your mesh seems alright. Therefore it is most likely a boundary condition.
You set T at your outlet to fixedValue. This needs to be zeroGradient. You can not define it as fixedValue on both inlet and outlet. Did not even see this before. For p_rgh you might try setting hotPlate to fixedFluxPressure. Everything within the p file should be set to calculated. |
Quote:
your reply is most appreciated. |
You should add the changes I talked about. Change the Temperature outlet boundary condition to zeroGradient. Change p_rgh at walls to fixedFluxPressure etc.
|
Hello everyone,
I am using CompressibleInterFoam for a test case of a rectangular computational domain. In my case instead of what Shipman had, I have two walls on the left boundary and an Inlet in the middle. Top and Bottom Patches are free surfaces. And the Right patch is outlet. And as a test case I initialize all the regions to be air. I was getting negative temperature before but using zeroGradiaent at the inlet fixed that problem. (Thanks for all the useful information here) The temperature decreases rapidly but does not become negative anymore. However, I am getting the following error: Code:
Courant Number mean: 0.017525563 max: 0.36055112 alpha.water Code:
object alpha.water; Code:
dimensions [0 1 -1 0 0 0 0]; Code:
dimensions [1 -1 -2 0 0 0 0]; Code:
dimensions [1 -1 -2 0 0 0 0]; Code:
dimensions [0 0 0 1 0 0 0]; I would appreciate if someone could give me some guidance on how to resolve the issue and help me understand why this error is happening. Thanks to all, |
Have you solved your problem. I meet the same problem and cannot figure where the problem is? Your reply is highly appreciated.
|
Hi everyone,
I am having the same problem with diverging temperature and velocity using compressibleInterFoam. I am using OF v5. I am using uniformFixedValue for p_rgh, calculated for p, zeroGradient for T, U and alpha but the simulation still crashes after few time-steps. I have tried several schemes but the simulation fails at more or less the same point in time. Has anyone managed to solve this problem? |
Hey guys,
I know it is a old thread but be careful when you run simulation with Liquid N2. Check in the following path according to your version : /src/thermophysicalModels/thermophysicalProperties/liquidProperties/N2/N2.C If you have: mu_(32.165, 496.9, 3.9069, -1.08e-21, 10.0), you need to add a minus to the coeff in bold because without this correction it compute a dynamic viscosity over 1e+24 Pa.s which is wrong. See this link for more information: https://bugs.openfoam.org/view.php?id=3136 Maybe it can solve your problems with the velocity and temperature. Cheers, |
I am now running a different simulation for internal flow in compressibleInterFoam using air and water.
I still end up having the same problem, i.e. velocity and pressure keep on increasing and increasing until the simulation diverges :mad: I've tried all sorts of boundary conditions and changed solver settings but it still keeps having the same issues although at different point in the simulation. I've also suppressed solution of the temperature equation to try and single out the issue but it seems that temperature is NOT causing the problem... Can someone please tell me how they solved this problem?? :confused: :( |
Hi, I've been modifying the compressibleInterFoam solver recently as well, and ran into the problem you raised. I found that the calculations are correct when the grid is perfectly orthogonal, but quickly diverge when there are non-orthogonal parts of the grid. Did you solve this problem?
|
All times are GMT -4. The time now is 21:39. |