January 25, 2021 at 6:00 pmai0013Subscriber
I am simulating an airflow initiated by a pressure discharge from a canister @ 21 bar to a vacuumed vessel at 0.4 bar absolute. At the end of the discharge, I have a uniform atmospheric pressure in all domains, but different temperatures in the reservoirs (as has been already shown in the experiments)
Fluent estimates 145 K in the canister, and from an isoentropic relationship (gamma = 1.4) this should be:January 26, 2021 at 12:16 pmRobAnsys EmployeeLooking at the TUI messages the solver starts to struggle with the flow and then temperature diverges. Chances are the temperature hits a nonphysical value somewhere which then returns a zero Cp causing the failure. If the solver hadn't returned the zero value I suspect you'd see temperature limit warnings followed by divergence. Have a careful look at the flow field about 5-10 iterations before the crash (re-run from a saved data set) and see what's going on. nWhat are the particles for?nJanuary 26, 2021 at 12:41 pmai0013SubscriberDear Thank you for your kind explanation. The particles are because I'm trying to re-create the 20L Siwek Sphere experiment (used for dust explosions). To understand the flow I started with pure gas simulation, with no physics other than K-eps turbulence and temperature equation. I patched the regions accordingly, and I observe that this produces a shockwave entering into the vacuumed sphere. In this case, the simulation successfully represents the pressure-time history for both reservoirs (I'm comparing with OF too)nMy problem comes when I inject the particles alongside the gas. My coal-volatile-mixture was defined using the coal-calculator, with properties of the mixture defined by;ndensity = ideal-gas & Cp = mixing law, I am using the pressure-based solver with default URF's (except for T = 0.6) , 2nd Order Discretization schemes (except MUSCL for Temperature) and fixed dt = 1e-05s. I still don't understand when my flow hits a non-physical temperature if the gas is perfectly OK (Max temp across the shockwave ~ 550 K). Besides, I read somewhere that it could be due to a bad mesh, but I I already refined it and increased the cell quality up to a min of 0.4. Any other suggestion would be much appreciated.January 26, 2021 at 1:48 pmRobAnsys EmployeeHave a look at where the particles are. DPM is intended for a 10-14% max volume fraction, and I suspect your problem is down to having a lot of particle mass somewhere. Are you planning on burning the mixture? If not use inert particles. Note, multiphase and shocks are a very complex subject, and you may find the models don't pick everything up. nJanuary 26, 2021 at 2:07 pmai0013SubscriberThank you Yes, eventually I'll need to burn the mixture. At this stage, I'm looking at the particle dispersion only. As suggested I already tried to use inert instead of combusting particles. I suspected I added all the physics (+reactions, +radiation, +DPM interaction) without any benefit, but when I turned off all the extra things it didn't help too much, as I had the same error.nIn fact, a few authors have used fluent to run the same simulation as I, so I guess I'm still doing something wrong. nJanuary 26, 2021 at 2:51 pmRobAnsys EmployeeCheck mesh resolution, just because the cell quality is good doesn't mean the mesh is suitable. nMarch 3, 2021 at 7:33 amscholarSubscribernI am experiencing the same issue. I have three phases: mixture-air, carbon-solid, and silica-solid. Currently, I am not solving any reactions. I am solving mass, momentum, and energy transport. Due to Cp becoming zero for carbon-solid, I changed Cp to a constant value. The simulation is running fine.nBut when I change the carbon-solid to a mixture-phase (having only carbon-solid as species) and kept Cp as a constant value, the simulation is not running. I see that the temperature is falling to 1 K which is the least set in controls (see the attachment).nPlease let me know your suggestions on this.nnMarch 3, 2021 at 9:42 amRobAnsys EmployeeIf you have a species-mixture with only one material it might be awkward numerically. We solve for n-1 species so if you only have one material I'm not sure how the logic will work out. Add a second species with a small volume fraction, and same properties as the carbon and see what happens. nMarch 3, 2021 at 10:42 amscholarSubscriberHello Rob,nI tried using another species in the mixture-coal. Still, the error persists. Since I kept Cp constant, the Cp value is >0. But the divergence is still there.nThe temperature of phase-2 (mixture-coal) is attached. Also when I change the operating conditions (T = 1023) so that the reactor temperature is 1073 K. But that is not reflected from contours. Should I patch the temperature for the entire domain as an input to the operating temperature of the reactor?nnMarch 3, 2021 at 10:42 amMarch 4, 2021 at 10:43 amRobAnsys EmployeeOperating temperature is for the Bousinesq approximation, click on Help in the panel and it'll explain or link to an explanation of the various fields. nViewing 10 reply threads
Ansys Innovation Space
- You must be logged in to reply to this topic.
Boost Ansys Fluent Simulations with AWS
Computational Fluid Dynamics (CFD) helps engineers design products in which the flow of fluid components is a significant challenge. These different use cases often require large complex models to solve on a traditional workstation. Click here to join this event to learn how to leverage Ansys Fluids on the cloud, thanks to Ansys Gateway powered by AWS.
Earth Rescue – An Ansys Online Series
The climate crisis is here. But so is the human ingenuity to fight it. Earth Rescue reveals what visionary companies are doing today to engineer radical new ideas in the fight against climate change. Click here to watch the first episode.
Subscribe to the Ansys Blog to get great new content about the power of simulation delivered right to your email on a weekly basis. With content from Ansys experts, partners and customers you will learn about product development advances, thought leadership and trends and tips to better use Ansys tools. Sign up here.Trending discussions
- Suppress Fluent to open with GUI while performing in journal file
- Floating point exception in Fluent
- What are the differences between CFX and Fluent?
- Heat transfer coefficient
- Getting graph and tabular data from result in workbench mechanical
- Difference between K-epsilon and K-omega Turbulence Model
- The solver failed with a non-zero exit code of : 2
- Time Step Size and Courant Number
- Mesh Interfaces in ANSYS FLUENT
- error: Received signal SIGSEGV
Top Rated Tags
© 2023 Copyright ANSYS, Inc. All rights reserved.Ansys does not support the usage of unauthorized Ansys software. Please visit www.ansys.com to obtain an official distribution.