Skip to content
GitLab
Projects Groups Snippets
  • /
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
  • Sign in / Register
  • openfoam openfoam
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
  • Issues 420
    • Issues 420
    • List
    • Boards
    • Service Desk
    • Milestones
  • Merge requests 7
    • Merge requests 7
  • Deployments
    • Deployments
    • Releases
  • Wiki
    • Wiki
  • Activity
  • Graph
  • Create a new issue
  • Commits
  • Issue Boards
Collapse sidebar
  • Development
  • openfoamopenfoam
  • Issues
  • #328
Closed
Open
Issue created Dec 07, 2016 by Admin@OpenFOAM-adminMaintainer

Closed incompressible two-phase system: set pressure reference

I am testing the solver twoPhaseEulerFoam (version 3.0.x) in a closed system to see the visualize the settling process. In my case, I set that fluids are rigorously incompressible. Since the system is closed, I set a no-slip b.c. for velocities (fixedValue (0 0 0)) in all patches and a fixedFluxPressure b.c. for p_rgh. In the dictionary fvSolution, I set pRefCell and pRefValue equal to 0 and 1e5 (as in internalField), respectively. Using GAMG solver for the solution of the pressure-associated linear system, the simulation crashes after 1 time-step. Using PCG solver, the simulation goes on but with an important increase in the number of iterations. The only way that I have to use GAMG for the pressure linear system is to fix a pressure in a patch and avoid to specify pRefCell and pRefValue. But, of course, this is not rigorous since I am considering an incompressible closed system. I attached a case for testing. In the case test, I sett

In the userguide, it is written that "In a closed incompressible system, pressure is relative: it is the pressure range that matters not the absolute values. In these cases, the solver sets a reference level of pRefValue in cell pRefCell where p is the name of the pressure solution variable. Where the pressure is p_rgh, the names are p_rhgRefValue and p_rhgRefCell respectively". Nevertheless, "createSystem.H" requires a pRefCell and pRefValue, not a p_rghRefCell and p_rghRefValue. I have tried to modify the "createSystem.H" but the case asks for a "pRefCell" and "pRefValue".

Thanks for your attention, NicolòcaseTest.tar.gz

## Reattaching the author to the issue ticket: @nic92 ##

Edited Dec 11, 2019 by Kutalmış Berçin
Assignee
Assign to
Time tracking