November 15, 2019 at 3:11 amgabrielcamethruSubscriber
I am trying to simulate category 5 hurricane winds on a roof using ANSYS 2019 R2 Fluent. My project is to design roof-edge devices that can reduce suction pressure on the roof.
I am concerned about the accuracy of the results I have been obtaining, particularly because I'm uncertain about the wall functions and the correctness of my inflation layers to represent the actual boundary layers.
As you can see in the pressure contour below, there is suction pressure in a flow separation region above the roof. Note that the right border of the air domain is the inlet, with air going left towards the house at 70 m/s.
I have chosen to use non-equilibrium wall functions (NWF) for the near-wall treatment since I read that these wall functions consider flow separation. Use of NWF requires y+ values between 30 and 300 so that the centroids of the cells adjacent to walls may lie in the log layer of the boundary layer. I have assigned each edge a y+ value (ranging from 30 for Re = 500,000 and 300 for Re = 10^. When calculating the Re of each edge, I used corresponding edge length as the characteristic length.
To refine the mesh, I have been applying edge sizing and inflation layers on each edge other than the top border of the domain (the sky). I have been using equations for turbulent flow over a flat plate to calculate boundary layer thickness, δ, for each edge, as well as to calculate the wall shear stress needed to find y, the distance from the wall to the centroid of the cells in the adjacent layer. Obviously, these calculations can only give me very rough estimates for δ and y.
Then, I have used the first layer thickness inflation option and set the first layer height equal to 2y (so that the centroid of the cells in the first layer equals y), and each edge has 20 inflation layers whose combined thickness meet or exceed the calculated δ of that edge.
I've noticed that, except for the ground, my approach has not been accurate in keeping y+ values between 30 and 300.
So, to sum up my questions:
1) Is it okay to be using non-equilibrium wall functions for this project? Or should another near-wall treatment be used?
2) How much should I concern myself with the y+ values? What tips would you recommend?
3) I was wondering how the reference values in setup change my results... Since my simulation is 2D, does changing Area or Depth do anything? Or do they only affect drag/lift coefficients (which I have no interest in) ?
Any other information or suggestions would also be greatly appreciated. Thank you beforehand
November 15, 2019 at 11:47 amRobAnsys Employee
Have a look at Enhanced Wall Treatment.
Also, think about the flow in a hurricane: is the velocity constant or does it fluctuate?
- 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
- 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
© 2023 Copyright ANSYS, Inc. All rights reserved.