CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT > Fluent Multiphase

Results doesn’t converge with fine boundary layer mesh

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By CFDKareem

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 24, 2022, 19:21
Default Results doesn’t converge with fine boundary layer mesh
  #1
New Member
 
Waru
Join Date: Apr 2022
Posts: 7
Rep Power: 3
modorous is on a distinguished road
Hi,
I'm stuck on this for a few weeks now, any help you can give is much appreciated.
As shown below, a liquid coolant is poured onto a structure. Liquid coolant and air are the two phases. Initially, domain is only filled with air. I'm looking at the heat transfer between the liquid and the surface. See attached pictures of the setup "Boundry_conditions" & "geometry".

This is a transient double precision simulation with PISO, Presto!, and Geo-reconstruct as the method. Turbulant model is K omega SST. Multi-phase is solved as a VOF with Implicit formulation, Implicit Body force and bounded second-order implicit transient formulation.
To get an accurate heat transfer coefficient I am using fine boundary layers in my mesh (First height - 0.002mm, layers 7, transition ratio 0.3). Mesh has an orthogonal quality of 0.18 and an aspect ratio of 226 (I'm using Fluent Meshing watertight geometry). See attached "mesh"& "BL_mesh_zoom" to check my mesh.

With these settings, all the parameters converge except for continuity at first. After several iterations, all the parameters started to diverge. I also ran a few simulations after increasing the first height of the boundary layer. With higher "First Height" results converge at first. But after around 100 iterations results start not to converge. I also checked the volume fraction of the domain at this point. This is the point where the secondary phase reaches the boundary layer. Finally, I ran a simulation without any boundary layers. Which ran without any issues, but gave wrong heat transfer results.

Is this due to my boundary layer mesh? as you can see in the above picture, I don't have much room to add a lot of layers. can someone tell me what I am doing wrong here? How can I solve this issue?

Thanks,

Best regards,

Waruna
Attached Images
File Type: jpg Boundry conditions.jpg (40.7 KB, 9 views)
File Type: jpg geometry.jpg (71.0 KB, 12 views)
File Type: jpg Mesh.jpg (88.0 KB, 9 views)
File Type: jpg BL_mesh_zoom.jpg (68.1 KB, 8 views)
modorous is offline   Reply With Quote

Old   December 27, 2022, 13:48
Default
  #2
Senior Member
 
CFDKareem's Avatar
 
Kareem
Join Date: Nov 2022
Location: New York
Posts: 115
Rep Power: 3
CFDKareem is on a distinguished road
Quote:
Originally Posted by modorous View Post
Hi,
I'm stuck on this for a few weeks now, any help you can give is much appreciated.
As shown below, a liquid coolant is poured onto a structure. Liquid coolant and air are the two phases. Initially, domain is only filled with air. I'm looking at the heat transfer between the liquid and the surface. See attached pictures of the setup "Boundry_conditions" & "geometry".

This is a transient double precision simulation with PISO, Presto!, and Geo-reconstruct as the method. Turbulant model is K omega SST. Multi-phase is solved as a VOF with Implicit formulation, Implicit Body force and bounded second-order implicit transient formulation.
To get an accurate heat transfer coefficient I am using fine boundary layers in my mesh (First height - 0.002mm, layers 7, transition ratio 0.3). Mesh has an orthogonal quality of 0.18 and an aspect ratio of 226 (I'm using Fluent Meshing watertight geometry). See attached "mesh"& "BL_mesh_zoom" to check my mesh.

With these settings, all the parameters converge except for continuity at first. After several iterations, all the parameters started to diverge. I also ran a few simulations after increasing the first height of the boundary layer. With higher "First Height" results converge at first. But after around 100 iterations results start not to converge. I also checked the volume fraction of the domain at this point. This is the point where the secondary phase reaches the boundary layer. Finally, I ran a simulation without any boundary layers. Which ran without any issues, but gave wrong heat transfer results.

Is this due to my boundary layer mesh? as you can see in the above picture, I don't have much room to add a lot of layers. can someone tell me what I am doing wrong here? How can I solve this issue?

Thanks,

Best regards,

Waruna
Your settings look good for the boundary conditions/discretization. It is possible that the quality of your boundary layers are causing the divergence in the calculation.

I am not very familiar with Fluent meshing so I am unsure of the algorithm it uses to create the inflation layers. However, looking at your mesh it looks like there are some mesh faces cutting across the inflation that may be creating some small mesh elements. See the picture attached with marks. These small mesh elements can definitely cause some issue with convergence.

If you use the "Check Mesh" in Fluent does it give any warnings about small mesh elements?

If you can choose the inflation algorithm in Fluent meshing try using a "pre inflation" algorithm. This will create the inflation layers on the wall before creating the interior mesh. It will help avoid the cutting of the hex mesh elements in the inflation layer.

Finally, I would try using workbench meshing to create the same mesh. You won't be able to use ploy mesh elements, but would be a good check to confirm that the rest of your setup is working well.

Let me know how it works out.
Attached Images
File Type: jpg VOF_MeshBoundary.jpg (71.8 KB, 7 views)
modorous likes this.
__________________
Please like the answer if it helped!

Video Tutorials and Tips: https://www.youtube.com/@cfdkareem/featured
CFDKareem is offline   Reply With Quote

Reply

Tags
fluent meshing watertight, heat transfer, inflation layers, mesh 3d, vof multiphase

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[snappyHexMesh] No generation of layer onto zerothickness surface crubio.abujas OpenFOAM Meshing & Mesh Conversion 3 October 25, 2022 04:20
Wind turbine simulation Saturn CFX 58 July 3, 2020 02:13
Fine tetrahedral mesh with boundary layer crashes Akanksha90 OpenFOAM Running, Solving & CFD 0 September 17, 2019 05:08
[ICEM] 3D Dynamic Mesh - Boundary layer mesh issues nathanricks ANSYS Meshing & Geometry 0 September 23, 2015 06:14
Question about heat transfer coefficient setting for CFX Anna Tian CFX 1 June 16, 2013 07:28


All times are GMT -4. The time now is 02:31.