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

Issue with MULES alpha

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 16, 2024, 13:37
Default Issue with MULES alpha
  #1
New Member
 
Michael
Join Date: Mar 2024
Posts: 5
Rep Power: 2
GardenLake is on a distinguished road
Hello everyone,

I am having an issue with the interfoam solver. I have a water flowing through the model with an inlet and outlet. I also have air on top of the water. When I have flow=0, I have no issue running. When I start flow (specifying an inlet velocity), I have an issue where my max alpha in MULES goes above 1 (which shouldn't happen, and my pressure goes crazy.

Has anyone had a similar issue and had a fix? I think it may have to do with the interface between the phases (air and water) changing but am unsure.

Thank you for your help!
GardenLake is offline   Reply With Quote

Old   March 17, 2024, 09:36
Default
  #2
Member
 
Michael Sukham
Join Date: Mar 2020
Location: India
Posts: 79
Rep Power: 6
2538sukham is on a distinguished road
Any chance of sharing the case or sharing the U p_rgh (or pd )?
2538sukham is offline   Reply With Quote

Old   March 21, 2024, 19:56
Default
  #3
New Member
 
Michael
Join Date: Mar 2024
Posts: 5
Rep Power: 2
GardenLake is on a distinguished road
Hello,

So I was able to fix the problem, I had forgot to designate a boundary and it was "default face".

Now I am having issues with my alpha value gets very high and low (2*10^18 and -2*10^18) and then crashing. In addition, I have a "wall" with a zero gradient designation getting alpha values below 1 but this is before crash (it is under water by quite a ways, should not be getting anything below 1). I would post a screenshot but can not figure out how to do that, it says I keep getting denied.

Does this problem ring a bell to anyone?

I have one inlet, an outlet, a wall, and an atmospheric top.
Alpha BC's

Inlet:type fixedValue;
value uniform 1;
Outlet: zeroGradient
Wall:zeroGradient
Atmospheric Top: type inletOutlet;
inletValue uniform 0;
value uniform 0;

Do these boundary conditions seem reasonable?
GardenLake is offline   Reply With Quote

Old   March 21, 2024, 21:05
Default
  #4
New Member
 
Michael
Join Date: Mar 2024
Posts: 5
Rep Power: 2
GardenLake is on a distinguished road
It is supposed to be open channel flow.
My U boundary conditions are...

wall: fixed value uniform (0 0 0)
inlet: fixed value uniform (1 0 0)
outlet: zeroGradient
Top: pressureInletOutletVelocity uniform (0 0 0)

My P boundary conditions are...
wall: fixedFluxPressure
inlet: fixedFluxPressure
outlet: fixedFluxPressure
Top: total pressure uniform 0;


Could the issue be with one of these?
GardenLake is offline   Reply With Quote

Old   March 21, 2024, 22:04
Default
  #5
Member
 
Michael Sukham
Join Date: Mar 2020
Location: India
Posts: 79
Rep Power: 6
2538sukham is on a distinguished road
For the pressure in outlet, try the zeroGradient bc. Other bc's look fine.
With regards
2538sukham is offline   Reply With Quote

Old   March 21, 2024, 23:34
Default
  #6
New Member
 
Michael
Join Date: Mar 2024
Posts: 5
Rep Power: 2
GardenLake is on a distinguished road
Hello,

Thank you for the response!

I tried the zeroGradient but now I am getting pressure GAMG value to crash. My alpha values at the crash are still quite large (+- 10^9). My crash occurs after about 0.4 seconds of runtime.

Do you have any other suggestions on what may be the issue? I tried "check mesh" and it came back with everything being ok
GardenLake is offline   Reply With Quote

Old   March 22, 2024, 08:23
Default
  #7
Member
 
Michael Sukham
Join Date: Mar 2020
Location: India
Posts: 79
Rep Power: 6
2538sukham is on a distinguished road
Can you share the case here. I would like to run and check it. I am not an expert but maybe I can look at the case and report it to you.
With regards. (Share zip file). I would put zeroGradient for the inlet as well for pressure.
2538sukham is offline   Reply With Quote

Old   March 25, 2024, 11:35
Default
  #8
New Member
 
Michael
Join Date: Mar 2024
Posts: 5
Rep Power: 2
GardenLake is on a distinguished road
Does anyone have an idea on what my pressure boundary condition should be for a boundary that is immediately above the inlet (which has wetting fluid going through it at a specified velocity), and is just exposed to what should be a normal atmosphere?
GardenLake is offline   Reply With Quote

Reply


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
[OLAFLOW] The OLAFLOW Thread Phicau OpenFOAM Community Contributions 457 March 27, 2024 00:59
Supercavitation / overPimpleDyMFoam canik01@yahoo.com OpenFOAM Running, Solving & CFD 28 September 6, 2021 11:07
multiphaseEulerFoam (OF2.3.0) : Courant number explodes when running in parallel Mehrez OpenFOAM Running, Solving & CFD 10 May 18, 2016 11:44
Stuck in a Rut- interDyMFoam! xoitx OpenFOAM Running, Solving & CFD 14 March 25, 2016 07:09
Alpha to be phased out Heiko Gerhauser Main CFD Forum 17 July 25, 2001 14:40


All times are GMT -4. The time now is 18:59.