-
-
November 8, 2023 at 6:04 amae22b001SubscriberHello all, The inlet of my geometry has a supersonic nozzle. I am trying to simulate flow through it in steady state initially and then later inÂunsteady conditions.With the density based solver I was able to get a decent steady state solution. When I started to run the transient solution , my results start diverging near the wall, which I think is caused by high accept ratio cells near the wall/inflation layers and sharp corners. Without the inflation layers, it is working.ÂHowever, when I use the pressure based solver I am able to get really good results with inflation layers as well. I want to recreate the results using density based solver as Mach number in my flow exceeds 3.Please let me know how to fix this issue.I will attach pictures
-
November 10, 2023 at 1:27 pmFedericoAnsys Employee
How are you initializing your solution for the transient case?
Also, what happens if you reduce the Courant number to 1?
-
November 14, 2023 at 9:39 amMatt JonesSubscriber
Hello any ideas on how to fix, thanks
-
-
November 10, 2023 at 2:02 pmMatt JonesSubscriber
Â
Â
Â
Hello, Thanks for replying.
Im using the steady state solution to start the transient solution. I was able to obtain the steady state solution by slowly increasing the courant no from 5 to 15.ÂÂ
Â
-
November 10, 2023 at 2:05 pmMatt JonesSubscriber
Â
after obtaining the steady state solution, without using any udf (so nothing is changing and everything shld instantly converge), I tried out a transient simulation to make sure everything was working, but even with a courant number of 1, the residuals shot up.
Â
-
November 11, 2023 at 4:04 pmMatt JonesSubscriber
Hello any ideas on how to fix?
-
November 14, 2023 at 1:53 pmFedericoAnsys Employee
You mentioned the aspect ratio of your inflation layers. How large is this ratio? What is your y+?
-
November 14, 2023 at 2:07 pmMatt JonesSubscriber
y+ is less than 2 on the walls. Ascept ratio is less than 100.
-
November 15, 2023 at 2:28 pmFedericoAnsys Employee
Looks good to me.Â
What are your solver settings for Transient case?
-
-
-
November 15, 2023 at 2:30 pmae22b001Subscriber
Same as steady state but with timestep of 1e-6, second order timestepping and courant no 5.
Â
-
November 15, 2023 at 3:23 pmFedericoAnsys Employee
I don't see any issues with the information that you are giving me. You can try reducing Turbulent URF and/or the time step to see if this helps with stabilization.
Side note: the pressure-based solver nowadays may be just as accurate even for flows M ~ 3.
-
November 15, 2023 at 4:15 pmMatt JonesSubscriber
Â
I tried using a timestep of 1e-7, and it seemed to be working better. The only problem is that I’ve seen multiple papers use a timestep of 1e-5 and 1e-6. So it seems weird that this problem is happening to me and that it is diverging from the inflation layers.
The bigger problem is that after I get the steady state solution, I simply changed it to transient with a timestep of 1e-5 (nothing else is changing). Ideally it shld converge instantly, but after a few iterations it started diverging.
I saved the residuals for postprocessing and observed where the residuals were high. and it was really high in some parts of the inflation layers.
This is confusing because none of these problems happen when using pressure based solver -
November 15, 2023 at 4:17 pmMatt JonesSubscriber
https://www.researchgate.net/publication/366732294_Effects_of_isolator_length_on_pseudo-shock_wave_in_a_rectangular_duct
for your reference, this is the paper i am talking about -
November 15, 2023 at 4:39 pmFedericoAnsys Employee
time step size will depend on the mesh, so unless you have identical mesh size, the time step size may not be comparable.
-
November 15, 2023 at 5:58 pmae22b001Subscriber
Yes but The bigger problem is that after I get the steady state solution, I simply changed it to transient with a timestep of 1e-5 (nothing else is changing). Ideally it shld converge instantly at every timestep, but after a few iterations it started slowly diverging.
-
November 15, 2023 at 6:22 pmFedericoAnsys Employee
Not sure what you mean by "it should converge instantly at every time step".
As I said, time step size is related to the mesh size. The fact that you see improvement when reducing dt and the fact that the solution is stable when removing inflation layers (which are your smallest cells), really points to the time step being too large for me.
-
November 16, 2023 at 1:42 amae22b001Subscriber
Understood, Thank you for your help
-
-
-
- The topic ‘Diverging Results with density based solver’ is closed to new replies.
- error udf
- Help: About the expression of turbulent viscosity in Realizable k-e model
- Unburnt Hydrocarbons contour in ANSYS FORTE for sector mesh
- Fluent fails with Intel MPI protocol on 2 nodes
- Cyclone (Stairmand) simulation using RSM
- Diesel with Ammonia/Hydrogen blend combustion
- Non-Intersected faces found for matching interface periodic-walls
- Mass Conservation Issue in Methane Pyrolysis Shock Tube Simulation
- Encountering Error in Heterogeneous Surface Reaction
- How to obtain axial and tangential velocity in CFX-post?
-
1131
-
468
-
466
-
225
-
201
© 2024 Copyright ANSYS, Inc. All rights reserved.