May 30, 2022 at 12:18 pmAnanth NarayanSubscriberHello all I am working on post-graduate project on film cooling on flat plate (without rotation). I tried with various pressure-velocity coupling schemes, but even after 1.5 day my result is not attaining steady state, and the contours of temperature is not changing after certain iterations even after 3-4 hours of solving them.
Flow will be turbulent and hence for y+ = 1, first cell spacing from wall has been made, and aspect ratio (max) is around 500, total mesh cell count around 1.3 million. Image of geometry is attached below [Paper]: (here D = 6.35 mm)
I tried almost all possible ways I know to attain it steady state, but nothing worked. The possible combinations which I tried are listed below:
1. When I start to reduce the mesh size, I see that steady state convergence is not reached.
2. For SIMPLE, PISO residuals diverge to 1e5 and I stopped the calculations. But for COUPLED and SIMPLEC, residuals are converging but after continuity residuals 1e-3 they are not converging anymore.
3. For COUPLED with and without Pseudo-Transient, I am not getting steady state convergence. For both my continuity residuals attains 1e-3 and after that convergence doesn't occur. For your reference, I am attaching an image without Pseudo-Transient after 55000 iterations. Also even after so long time, I find the same contour.
4. For SIMPLEC, same thing happens as of COUPLED solver.
5. If I do transient analysis, for time step above 1e-3- solution starts diverging, but when I go below that it takes so so long time.
For your reference, I am attaching similar kind of simulation I did with other geometry:
The paper linked uses temperature piecewise polynomial for Cp, K, Dynamic Viscosity and density is considered as ideal gas. Since both fluid coming at different density and velocity this is one of the way to replicate the phenomenon. Else we can directly use species transport if density and other parameters are not relatable to each other.
I have used piecewise polynomial relation so as to validate Fig 6(a) for Re = 8e4. As every parameters are know we can find mainstream velocity by, (here 6.35 mm is diameter of cooling tube)
How do I achieve the steady state solution, am I missing something or do I need more computational power. I am using Intel i7 with 16 GB, 8 Processor CPU.
For your reference, I am attaching case file. (Please replace .txt to .cas.h5)
Thanks in advance
Regards Ananth Narayan
M.Tech Thermal Engineering
May 31, 2022 at 11:53 amSwathi V. V.Ansys EmployeeHi Ananth Ensure the mesh quality is atleast 0.1. Aong with residual, please check if flux imbalance is low and create monitor points to determine if the solution is converged. Refer to the following documents for more details: 32.21.1. Judging Convergence (ansys.com) and 32.15.3. Monitoring Solution Quantities (ansys.com). You can initally run steady state simulation and then switch to transient to start with a reasonable initial guesses and for fasterconvergence. If you are facing divergence, start your simulation with lower under relaxation factor and then as the solution gets stable increase it to default value.
You can use adaptive time stepping to change time step size as the calculation proceeds and solution attains stability. Please refer to 32.13.1. Inputs for Time-Dependent Problems (ansys.com).
June 1, 2022 at 2:40 amAnanth NarayanSubscriberHello, thanks for initiating answer for this question. Does mesh quality is orthogonal quality?. I am using ICEM CFD for meshing and min skewness is 0.5.
If suppose I reduce the relaxation factor and run steady state and then switch to transient state, do I need to follow adaptive time step or fixed time step. What is the reasonable CFL number I should be using. If not can I use expression for time step as the time increases, time step size should decrease?
- 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.