June 26, 2020 at 2:50 amhaitengSubscriber
I am using ANSYS Fluent to simulate the experiment for natural convection of a single heated cylinder in a large water tank made of Plexiglas (approximated as thermal-insulating material), which is 1200 mm wide with a water depth of 1600 mm. Computational domain is a rectangular with the same width and depth as the experiment. The thrid dimension is neglected because it is assumed that spanwise flow field is uniform, which is reasonable. I use steady, pressure-based solver and Boussinesq approximation for buoyancy. Turbulence is modelled by transition SST. Operating temperature is set to 293 K. Density is Boussinesq at 998.92 kg/m3, which is evaluated at the operating temperature. Gravity is specified verticall downward.
The major problem I have is how to set the boundary condition correctly on the four sides of the rectangular domain. I have tried two options.
Option 1: No-slip, adiabatic wall on left, right & bottom boundary. Zero-shear, isothermal (Tw=293 K) wall on top boundary, to simulate the free water surface. This boundary condition setting is most representative of real scenario in experiment. However, the result is non-sensible. Contour of velocity magnitude suggests buoyant plume is completely disrupted, as shown below.
Option 2: No-slip, adiabatic wall on left & right boundary. Pressure inlet at bottom boudanry. Pressure outlet at top boundary. This setting, although not representative of experimental scenario, produces sensible result. Buoyant plume show reasonable profile and velocity magnitude, as shown below.
Could anyone tell me, why Option 1 generates the non-sensible result, although it is exactly the same as experiment? Why Option 2 seems so good, although the boundary condition does not reflect the real scenario?...
June 26, 2020 at 5:20 amDrAmineAnsys EmployeeWhat about the temperature plots especially in the corners at top? Anyhow you need to run this case transient.
June 26, 2020 at 3:18 pmRKAnsys Employee
Please refer to the Fluent User's guide and make sure you have evaluated the density based on the model mentioned.
June 27, 2020 at 3:34 amhaitengSubscriber
Temperature contour is as follows and really weird. Temeperature in the ambient is around 300 K, although it is set to 293 K in the "Operating Temperature". Near the corner of the top, it decreases to 293 K (specified temperature at top boundary) within 1 grid point.
I don't think the transient solver helps, given that steady solver, as an approximation of the system, performs so badly...
June 27, 2020 at 3:36 amhaitengSubscriber
Thanks. I have read the guideline previously and followed it. I used Boussinesq approximation (commonly adopted in literature) so only operating temperature needs to be specified.
June 27, 2020 at 8:11 amDrAmineAnsys EmployeeThe natural boundary condition at top should be an opening. I expect a singularity at a adiabatic/fixed temp wall. Transient solver is the way to go for further analysis.
June 27, 2020 at 8:53 amhaitengSubscriber
Did you mean trying transient solver (URANS or LES) by using the same boundary condition (zero-shear isothermal wall at the top & no-slip adiabatic wall on other three sides)?
Singularity at the top corner might exist. But I don't understand why changing the solver from steady to transient help to solve this problem. Shouldn't we first debug boundary condition setting, or mesh grading near top boundary, using the less-expensive RANS steady solver?
June 28, 2020 at 2:42 am
July 15, 2020 at 12:29 pmhaitengSubscriber
The problem is solved by switching to transient solver. thanks!
- 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.
- 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.