CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Uy not converging in natural convetion peudo 1D melting case

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 11, 2021, 12:29
Default Uy not converging in natural convetion peudo 1D melting case
  #1
New Member
 
Kai Salscheider
Join Date: Aug 2018
Posts: 13
Rep Power: 7
Skaiwalker is on a distinguished road
Hello Forum,

Problem:
I want to simulate the melting of a magnesium ingot in a crucible, so I chose OF8 with its fvOptions "solidificationMeltingSource" to model the 2 phases in 1 fluid. The testcase is a slice (yellow marked area) out of a litte crucible as can be seen in the attached picture.

Convection is important so I want the flow to exit and enter the domain with cell value, as if z would be infinite. As of my understanding this leads to:
- "zeroGradient" on the velocity at top and bottom
- top and bottom: fixedFluxPressure as in chtMultiRegionFoam tutorial heatExchanger.

Edit:
This leads to jumping Uy residuals. The first timestep does not converge, then it runs for several steps just to diverge from 0.027s onwards and converging again from 1.3s onwards.

Question:
Is it crucial that all timesteps converge to a low residual? The pressure gradient in y is marginal and the absolute values of Uz are of magnitude e-09, so i think it's not easy for the numerics?

slice.jpg
testcase_meltingIngot.zip

Last edited by Skaiwalker; February 12, 2021 at 10:20. Reason: Wrong Material Properties in fvOptions and empty testcase, duh!
Skaiwalker is offline   Reply With Quote

Old   February 12, 2021, 10:28
Default
  #2
New Member
 
Kai Salscheider
Join Date: Aug 2018
Posts: 13
Rep Power: 7
Skaiwalker is on a distinguished road
The material properties in fvOptions were wrong, hence the higher velocity. With the right velocity I encountered a new error and edited the original post to suit it. And the testcase isn't empty anymore
Skaiwalker is offline   Reply With Quote

Old   October 1, 2021, 05:15
Default
  #3
Member
 
Join Date: Nov 2020
Posts: 52
Rep Power: 5
mikulo is on a distinguished road
Quote:
Originally Posted by Skaiwalker View Post
The material properties in fvOptions were wrong, hence the higher velocity. With the right velocity I encountered a new error and edited the original post to suit it. And the testcase isn't empty anymore
Hey,

Have you solved your problem? I am running with the same trouble too. It seems like it is not easy to reach convergence using sMS.

Regards,
Mike
mikulo is offline   Reply With Quote

Old   October 1, 2021, 12:11
Default
  #4
New Member
 
Anup Singh
Join Date: Mar 2020
Posts: 22
Rep Power: 5
Anup Singh is on a distinguished road
Just try to solve using the enthalpy porosity method.
Instead of using solidificationmeltingsource just add the source term to the temperature equation and update it according to the enthalpy-porosity or any update scheme.
You can even try to explicitly converge both temperature and liquid fraction at each update using a loop.
Anup Singh is offline   Reply With Quote

Old   November 11, 2021, 04:34
Cool
  #5
Member
 
Join Date: Nov 2020
Posts: 52
Rep Power: 5
mikulo is on a distinguished road
Quote:
Originally Posted by Anup Singh View Post
Just try to solve using the enthalpy porosity method.
Instead of using solidificationmeltingsource just add the source term to the temperature equation and update it according to the enthalpy-porosity or any update scheme.
You can even try to explicitly converge both temperature and liquid fraction at each update using a loop.
Hello,

Thanks for your response. What you said is another option too, though it seems that the fvOption sMS is programmed well. Anyway, do you have some test cases? I would be very glad to see and run it. Thanks!
mikulo is offline   Reply With Quote

Reply

Tags
boundary condition, boussinesq approximation, buoyantpimplefoam

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
Pressure instabilities with interDyMFoam for the floatingObject case nbadano OpenFOAM Running, Solving & CFD 15 October 15, 2021 07:35
Is Playstation 3 cluster suitable for CFD work hsieh OpenFOAM 9 August 16, 2015 15:53
Case (VOF + Solidification&Melting) will not work with Species Transport? ERROR ? andy.eco FLUENT 4 September 14, 2010 18:27
Free surface boudary conditions with SOLA-VOF Fan Main CFD Forum 10 September 9, 2006 13:24
melting case wanghong FLUENT 0 November 21, 2004 22:03


All times are GMT -4. The time now is 07:45.