|
[Sponsors] |
September 23, 2016, 10:36 |
what is the meaning of p_rgh?
|
#1 |
New Member
Simone Colucci
Join Date: Mar 2016
Location: Pisa (Italy)
Posts: 23
Rep Power: 10 |
Hi,
I have a simple question. In many solvers (for example in reactingTwoPhaseEulerFoam, OpenFoam 4.0) the PEqn solves an equation for p_rgh, instead of p. p_rgh is defined as: Code:
p_rgh = p - rho*gh; Code:
volScalarField gh("gh", (g & mesh.C()) - ghRef); Thank you in advance Simone |
|
September 26, 2016, 12:45 |
|
#2 |
Member
Hasan Celik
Join Date: Sep 2016
Posts: 64
Rep Power: 9 |
Can you clarify your point little bit more? I mean, I guess h is the difference in height in both cases you mentioned, so isn't it 10?
|
|
September 26, 2016, 22:15 |
|
#3 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 |
Hi Simone,
the truth is that it has no physical meaning, it is just a convenient mathematical artifact. It is explained in my thesis (page 111): https://repositorio.unican.es/xmlui/...=1&isAllowed=y Best, Pablo |
|
June 14, 2020, 04:07 |
|
#4 | |
Member
Himanshu
Join Date: Jan 2017
Posts: 34
Rep Power: 9 |
Quote:
So, it mean no matter what the co-ordinate system our model is in. the result will be same as long as our BC remain same. To elaborate my question, if I did simulation for one geometry and now I translate that geometry by 10 cm. Will the result be same if I perform simulation on this new translated geometry? |
||
June 14, 2020, 16:37 |
|
#5 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 |
Hi,
not necessarily. The main advantage is to avoid large round-off errors if numbers are large, as P can be. Best, Pablo |
|
June 14, 2020, 16:55 |
|
#6 | |
Senior Member
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 932
Rep Power: 12 |
Quote:
- The governing equations of fluid dynamics are Galilean invariant, which means that making any transformation or rotation on the computational domain within an inertial frame of reference will not change the predictions if the boundary conditions are modified accordingly. - Therefore, translating the domain by 10 cm in a given direction will not change the numerical predictions if the control volume keeps modelling the same volume.
__________________
The OpenFOAM community is the biggest contributor to OpenFOAM: User guide/Wiki-1/Wiki-2/Code guide/Code Wiki/Journal Nilsson/Guerrero/Holzinger/Holzmann/Nagy/Santos/Nozaki/Jasak/Primer Governance Bugs/Features: OpenFOAM (ESI-OpenCFD-Trademark) Bugs/Features: FOAM-Extend (Wikki-FSB) Bugs: OpenFOAM.org How to create a MWE New: Forkable OpenFOAM mirror |
||
June 15, 2020, 08:49 |
|
#7 |
Member
Himanshu
Join Date: Jan 2017
Posts: 34
Rep Power: 9 |
Thank You for the replies.
|
|
June 15, 2020, 09:07 |
|
#8 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 |
Hi HPE,
despite the equations being Galilean invariant, once you discretize them it is another story. The round off errors and the solution will be different. To what extent is a question of the processes and how long you run the simulation for, but "chaos" plays a role here and results may be noticeably different. For example, you can observe this in my simulations when waves break. Best, Pablo |
|
June 15, 2020, 11:18 |
|
#9 |
Senior Member
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 932
Rep Power: 12 |
Hi,
There are circumstances where Galilean invariance breaks down for the finite volume method in Eulerian frame, true; but from an engineering perspective, the reason in the current context is the acute sensitivity of the governing equations of particular physics on initial/boundary conditions. What you have observed can also happen for a restarted simulation wherein the restart is not bitwise, which is very common in any simulation of any software. Personally, I would not be concerned with such ultra academic break downs of Galilean invariance in discretised domains since even the Navier-Stokes equations have only weak solutions, which means there might be some flow behaviour not obeying N-S in some part of universe at some point. But, who cares. Thank you for your kind remarks.
__________________
The OpenFOAM community is the biggest contributor to OpenFOAM: User guide/Wiki-1/Wiki-2/Code guide/Code Wiki/Journal Nilsson/Guerrero/Holzinger/Holzmann/Nagy/Santos/Nozaki/Jasak/Primer Governance Bugs/Features: OpenFOAM (ESI-OpenCFD-Trademark) Bugs/Features: FOAM-Extend (Wikki-FSB) Bugs: OpenFOAM.org How to create a MWE New: Forkable OpenFOAM mirror |
|
November 6, 2020, 21:21 |
|
#10 | |
Member
Hasan Celik
Join Date: Sep 2016
Posts: 64
Rep Power: 9 |
Quote:
|
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Who can tell me what is the meaning of the macro in fluent UDF:C_PREMIXC_RATE(c,t) | ChenZhan | Fluent UDF and Scheme Programming | 1 | May 12, 2017 03:40 |
What is the meaning of pd in foam-extended 3.1 | lxwd | OpenFOAM Running, Solving & CFD | 2 | February 3, 2017 11:54 |
Meaning of Ychar and Ypmma | aylalisa | OpenFOAM Pre-Processing | 2 | October 20, 2013 05:49 |
What's meaning of UDF FUNCTION | zhaoxinyu | Fluent UDF and Scheme Programming | 0 | March 31, 2010 08:04 |
want to know meaning | Sangamesh | Siemens | 0 | May 15, 2007 05:15 |