

November 19, 2023 at 2:54 pmjasonsum25Subscriber
Hi, I am simulating the wind flowing across a Yshape building to evaulate the wake effect behind.
Mesh was done in Fluent Meshing, with minimum orthogonality of 0.2, and 5 inflation laters added to ground and building surface.
with BOI refinement in volume near the building
To be more accurate towards reality, I want to take the atmospheroc boundry layer for the wind flow inlet condition.
However, divergence and float points popped up on console when initiaizing. (while just using constant input to inlet velocity did not cause such issues), it is earnestly hoped for some guidance on this and some additional tips on using UDF, thanks!!!!
Error codes:
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0
Divergence detected in AMG solver: ads0 10 0.000000e+00
Hybrid initialization is done.
Error at host: floating point exception
Error at Node 4: floating point exception
Error at Node 0: floating point exception
Error at Node 1: floating point exception
Error at Node 2: floating point exception
Error at Node 3: floating point exception
Error at Node 5: floating point exception
Error at Node 6: floating point exception
Error at Node 7: floating point exception
Error at Node 8: floating point exception
Error at Node 9: floating point exception
Error at Node 10: floating point exception
Error at Node 11: floating point exception
Error at Node 12: floating point exception
Error at Node 13: floating point exception
Error at Node 14: floating point exception
Error at Node 15: floating point exception
Error at Node 16: floating point exception
Error at Node 17: floating point exception
Error at Node 18: floating point exception
Error at Node 19: floating point exception
===============Message from the Cortex Process================================
Compute processes interrupted. Processing can be resumed.
==============================================================================
Error: floating point exception
Error Object: #f

Below is the UDF I taken reference from
&t=2395s:#include "udf.h" #define UREF 4.491 // reference wind speed at ZREF #define ZREF 0.16 // height of the building, 2b #define LREF 0.08 // characteristic length, b #define CMU 0.09 // DEFINE_PROFILE(velocity_profile, thread, position) { float x[ND_ND]; float y; float u; face_t f; begin_f_loop(f, thread) { F_CENTROID(x,f,thread); y=x[1]; u = UREF*pow(y/ZREF,0.27); F_PROFILE(f,thread,position) = u; } end_f_loop(f, thread) } /* profile for kinetic energy, k */ DEFINE_PROFILE(k_profile, thread, position) { float x[ND_ND]; face_t f; float A1 = 0.0042 ; float A2 = 0.0719 ; float A3 = 0.2692 ; float A4 = 0.3671 ; float k, y, yb; begin_f_loop(f, thread) { F_CENTROID(x,f,thread); y = x[1]; yb = y/LREF; k = A1*pow(yb,3)+A2*pow(yb,2)+A3*yb+A4 ; F_PROFILE(f,thread,position)=k; } end_f_loop(f, thread) } /* profile for dissipation rate, epsilon */ DEFINE_PROFILE(dissip_profile, thread, position) { float x[ND_ND]; face_t f; float A1 = 0.0042 ; float A2 = 0.0719 ; float A3 = 0.2692 ; float A4 = 0.3671 ; float k, y, yb; begin_f_loop(f, thread) { F_CENTROID(x,f,thread); y=x[1]; yb = y/LREF; k = A1*pow(yb,3)+A2*pow(yb,2)+A3*yb+A4 ; F_PROFILE(f,thread,position)=pow(CMU,0.75)*pow(k,1.5)/LREF; } end_f_loop(f,thread) } /* References: [1] Yan MENG and Kazuki HIBI. Turbulent measurements of the flow field around a highrise building. Journal of wind engineering. No 76, July 1998. */

November 20, 2023 at 12:31 pmRobForum Moderator
You may want to check your direction vector for "up". If z is up it's not [1] in the vector definition!

 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.
Ansys Blog
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.
 legend min and max
 Ensight hot iron palette from an image
 Streamlines in EnSight using MRI data
 Import MRI data into Ensight
 FLUENT APPLICATIION ERROR
 Total Surface Heat Flux Calculation in Fluent
 Drop Test of a WaterFilled Tube
 Difference between “total pressure” and “absolute pressure”?
 Minimum Orthogonal Quality Less than 0.01 For Transonic Airfoil Flow Analysis
 obtaining pressure distribution by making points in ansys

8808

4658

3153

1688

1474
© 2023 Copyright ANSYS, Inc. All rights reserved.