View Issue Details

IDProjectCategoryView StatusLast Update
0004263OpenFOAMBugpublic2025-08-10 15:27
ReporterWesley Assigned Tohenry  
PrioritynormalSeveritymajorReproducibilityalways
Status closedResolutionno change required 
PlatformRyzen 7 7735HS processorOSKubuntuOS Version24.04
Product Version12 
Summary0004263: incompressibleMultiphaseVoF does not allow use of power law fluid model
DescriptionDescription of the issue: incompressibleMultiphaseVoF solver does not use the powerLaw fluid model values in the momentumTransport.fluid files

Basic simulation:
-T pipe, with 2 fluids entering from the 2 sides of the T, exiting the remaining path
-Both fluids are powerLaw fluids
-Both fluids have the same density
-Both fluids have the same volumetric flow rate
-k values differ by a factor of 10 between the two fluids

Expected results:
-The more viscous fluid will take more of the cross-sectional area of the pipe after the fluids are in contact
-Inlet p and p_rgh are expected to abe about 1.0E+08 with the powerLaw fluid properties (from simulation using incompressibleVoF for a two-phase simulation and twoPhaseTransport in the momentumTransport file)


Simulation results:
-Results are consistent with the rheology in physicalProperties.*, not the rheology in momentumTransport.*
-Each fluid takes 50% of the cross-sectional area
-Inlet p is about 4.9E+02
-Inlet p_rgh is about 2.0E+02

I have not been able to find a tutorial for multi-phase simulation for non-Newtonian fluids, and have been trying get the simulation running. I can get the simulation running in OpenFOAM 9 (before location of fluid properties to changed for non-Newtonian fluids), but have not been able to get it working for OF 10, 11, 12, or 13.

My assumptions that need to be checked:
- incompressibleMultiphaseVoF is supposed to allow for simulation with non-Newtonian fluids
- the files are setup correctly

Desired outcome:
- incompressibleMultiphaseVoF allows non-Newtonian fluids
- if incompressibleMultiphaseVoF does not allow non-Newtonian fluids, suggestion for other solver that does
Steps To ReproduceRun the attached simulation or create your own simulation with a power law fluid. The .stl file is included, and has already been put into the polyMesh directory.

For comparison, I have also included the file with the incompressibleVoF simulation. This simulation solves as I would expect from the planned properties and boundary conditions.
Additional InformationI am marking the severity as major because I can't run the types of simulations I would like, and suspect the issue carries over to all non-Newtonian rheologies.

Including more tutorials for solvers that allow using non-Newtonian fluids would be helpful for some users.

I have reproduced this issue on Debian 12 and on Kubuntu 24.04 on multiple computers.

I have reproduced this issue on OpenFOAM 10, 11, 12, and 13.

I have asked for help and guidance on cfd-online.com and have not received feedback on how to better set up the files or how to get the simulation working.
TagsNo tags attached.

Activities

Wesley

2025-08-10 13:25

reporter  

henry

2025-08-10 15:27

manager   ~0013645

User support request

Issue History

Date Modified Username Field Change
2025-08-10 13:25 Wesley New Issue
2025-08-10 13:25 Wesley File Added: Bug_Submission_incompressibleMultiphaseVoF_Solver.tar.gz
2025-08-10 13:25 Wesley File Added: Transient2PhaseIncompressibleVoF.zip
2025-08-10 15:27 henry Assigned To => henry
2025-08-10 15:27 henry Status new => closed
2025-08-10 15:27 henry Resolution open => no change required
2025-08-10 15:27 henry Note Added: 0013645