|
[Sponsors] |
![]() |
![]() |
#1 |
Member
Elwardi Fadeli
Join Date: Dec 2016
Location: Boumerdes, Algeria
Posts: 41
Rep Power: 9 ![]() |
Hello,
I'm stuck with a good-for-nothing-but-working-fine solver with no gravity capabilities for around a week now. I can't find where the problem is, so, please, please help me out ![]() I'm solving two-phase compressible flow in porous media, so, my first attempt to calculate the (water) flux was: Code:
// Kwf is a surfaceScalarField as it's important to do this // because of very specific "interpolation" schemes that // have to be used phiwG = (Kwf * rhowf * g) & mesh.Sf(); // Then, in the equation, I do fvScalarMatrix wEqn ( Coeff1*fvm::ddt(water.alpha()) + fvc::div(phiwG) + sources ); So, when I set g = (0 0 0), I get fine results. But when g = (0 0 -9.81), everything gets messed up, ie. bounded variables get out of their bounds. In my many attempts to debug the problem, I guessed either: - My calculation of phiwG is wrong (You can't mess up g, so Kwf or rhowf, but they are fine). - My calculation of other fluxes is wrong, and after inspection Code:
// This is not valid OpenfOAM code // units missmatch fvc::laplacian(-Kw, p)/(fvc::div(phiwG)+VSMALL) - rhow strongly depends on the pressure (which is a main variable), and I don't have enough loops ![]() Most of the other approaches, that are used by standard solver, ex. taking out rho (with snGrad) are not suitable in this case because even if the flow is reduced to be incompressible, difference of rho of both fluids still matters. If anyone has even a hint of a suggestion, please don't hold back. Best wishes, and thanks in advance. |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,978
Blog Entries: 45
Rep Power: 128 ![]() ![]() ![]() ![]() ![]() ![]() |
Quick answer: I came here after receiving the PM you sent me, but I'm not understanding what's going on either.
Because:
The other thing that comes to mind is that the use o "mesh.Sf()" might be wrong, given that it's the product of face area with face normal, so perhaps you were not meant to have the face area? Perhaps the "div(phiwG)" term is utterly incorrect, in the sense that maybe it should be something like "div(phi,g)". Or perhaps the wrong discretization scheme is being used... from the description you've given, nothing other than "upwind" should perhaps be used? Last but not least, perhaps you're not familiar with this page: https://openfoamwiki.net/index.php/HowTo_debugging ?
__________________
|
|
![]() |
![]() |
![]() |
![]() |
#3 |
Member
Elwardi Fadeli
Join Date: Dec 2016
Location: Boumerdes, Algeria
Posts: 41
Rep Power: 9 ![]() |
Thanks @wyldckat for your informative input; I will try to clarify things a bit; I know I've been over protective about my solver
![]() 1. The problematic solver is a coupled solver, similar, to some degree, to pUcoupledFoam in foam-extend. 2. "Kwf" is the face interpolation of what you, CFD guys, may call diffusion. It's closer to a "transmissiblity" really. The important thing that it has some factors that must be interpolated using a specific scheme (upwind & harmonic). 3. Well, I got ambitious and tried to solve the most complex form of equations there is. Then I went back and simplified them to two-phase, incompressible BlackOil equations: ![]() ![]() where ![]() ![]() ![]() ![]() ![]() 4. I tried neutralizing "rho" effect before by setting both densities to 1 but the solver fails anyway. 5. 'rhowf' is the face interpolation of water density (which I keep constant through the domain for the purpose of testing). ... And of course I know about GDB, and gdbFoam; I don't think I would get this far if I didn't use them ![]() Your suggestion that "div(phiwG)" term should be expressed in an other form got me thinking, so I went back and tested it on a "segregated" version of the solver (But the timeStep is so small It's painful to watch) and it worked fine. That's why I assumed It would work in a coupled solver in the first place. I'm very careful with discretization schemes and all my face fluxes need to be multiplied by mesh.magSf() somehow if I want to "div" them as the end result should be in dimless/dimTime dimensions ![]() Thanks for reading this long thread, but I'm starting to think my coupling of these equations is not entirely correct so I'm re-building on top of the working segregated solver; I'll let you know how that goes. |
|
![]() |
![]() |
![]() |
![]() |
#4 |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,978
Blog Entries: 45
Rep Power: 128 ![]() ![]() ![]() ![]() ![]() ![]() |
Quick answer: Unfortunately this really goes outside of my current level of experience.
One of the things that popped into my mind was that perhaps "phiwG" might be accounting as many times as faces per cell, for the actual divergence value... although it doesn't make much sense, given that there are other discretization terms handled this way... Take a look into the "nextRelease" branch in foam-extend 4.0, because they have additional solvers there for coupling... And I took a second look to the equations you provided and I'm confused as to which variables you're actually solving for here... and I got the feeling that you need to restructure how the coupling is being done. From what I can briefly see in the coupled solver you mentioned, there are terms that have to be coupled manually because they are present in more than one equation, so they have to be handled in a special way; others had to be constructed based on the other equation. This may be too late for you to attent, but just in case it's still on time, you might want to attend this, since the people that worked on this kind of problem-solving will likely be there or nearby: NUMAP FOAM Summer School 2019 The other group of people I believe have worked with this kind of thing are the ones who develop Caelus-CML. |
|
![]() |
![]() |
![]() |
![]() |
#5 | |
Member
Elwardi Fadeli
Join Date: Dec 2016
Location: Boumerdes, Algeria
Posts: 41
Rep Power: 9 ![]() |
Thanks a lot Mr. @wyldckat for your insightful assistance on this subject.
As we suspected, the coupling of the Code:
fvc::div(phiwG) Without your assistance, I know I would never caught the bug ![]() Now, things work as I expected, well, I still have some issues, but nothing that GDB can't handle. And by the way, Quote:
![]() |
||
![]() |
![]() |
![]() |
![]() |
#6 | |
Retired Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 10,978
Blog Entries: 45
Rep Power: 128 ![]() ![]() ![]() ![]() ![]() ![]() |
Quote:
I'm not entirely certain how at 1AM I was able to pinpoint that being one of the likely issues... we both got lucky ![]() |
||
![]() |
![]() |
![]() |
Tags |
gravity flow, gravity force |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
NS Equations for multiphase flow in porous media | quir | OpenFOAM Running, Solving & CFD | 3 | December 26, 2020 11:25 |
Multiphase Homogeneous Model equations | marcoymarc | CFX | 3 | July 3, 2013 19:11 |
Question about heat transfer coefficient setting for CFX | Anna Tian | CFX | 1 | June 16, 2013 06:28 |
Flux splitting | Dr B.M. Smith (Smith) | OpenFOAM Running, Solving & CFD | 19 | January 9, 2013 04:07 |
Replace periodic by inlet-outlet pair | lego | CFX | 3 | November 5, 2002 20:09 |