January 11, 2023 at 4:10 pmandrewcaratenutoSubscriber
I am trying to simulate a distillation process between two columns of water at different temperatures which are connected by a partition/valve. The system takes advantage of hydrostatic boiling, where tubes of water are raised up quite high (~10 m) to lower the pressures and produce phase change at close-to-ambient temperatures. My geometry represents the top ~0.1 m of the system to minimize complexity.
Gravity is enabled, and I have open pressure boundaries at the bottom of each tube. I've attached a schematic for clarity.
At the start, I initialize the system with only liquid water. Based on the low operating pressure, the effects of gravity, and my multiphase model (Lee model for boiling/condensation and species model for evaporation/concentration-driven mass transfer, with temperature-dependent/pressure-dependent properties), the water boils on both sides. After about 15 seconds, I get a layer of vapor at the top of each column. The pressure on each side corresponds to the saturation pressure of water at the given temperatures (which I believe means they are in the saturated vapor state).
Up to this point I believe the model is correct, but after I open the mesh interface to start the diffusion process, I get lost. From theory and experiments on these systems, vapor should diffuse from the hot side to the cold side because the vapor pressure of the hot side is greater, creating a concentration gradient. This is only possible due to the saturated vapor conditions on each side.
However, when I open the interface, the pressures equalize instantaneously. Their value is about the average of the two temperatures, but the temperature of the vapor region does not support this.
I'm hoping to get some ideas about how to treat saturated vapor in Fluent multiphase. I know I start with saturated vapor, but after I open the interface, the pressures and temperatures no longer correspond to saturation conditions. Any advice on why the pressures are equalizing, or how to preserve the saturation conditions, would be greatly appreciated.
Images of density and pressure state before partition opens are attached.
Model details on model:
Eulerian multiphase; Lee model for boiling/condensation (pressure-dependent saturation temperature via UDF); species model for evaporation/concentration-driven mass transfer (temperature-dependent saturation pressure via UDF); Pressure-based solver; dt = 0.001 s; All URFs default besides energy (0.1);
January 23, 2023 at 11:15 amPrashanthAnsys Employee
I'm writing as nobody else responsed. Have you tried switching off flow equations after remoning the partition? I mean just solve the species transport equations to model pure diffusion.
January 23, 2023 at 12:36 pmandrewcaratenutoSubscriber
Hi Prashanth, thanks for your reply, this is a good suggestion.
How exactly can I switch off the flow equations though? Is it a setting in the viscous model window, or somewhere else?
January 24, 2023 at 4:53 pmandrewcaratenutoSubscriber
Nevermind, I found it in controls > equations. I'll report back based on what I find with this strategy, thanks again.
- You must be logged in to reply to this topic.
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.
- 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
- The solver failed with a non-zero exit code of : 2
- Difference between K-epsilon and K-omega Turbulence Model
- Time Step Size and Courant Number
- Mesh Interfaces in ANSYS FLUENT
- error in cfd post
© 2023 Copyright ANSYS, Inc. All rights reserved.