
[Sponsors] 
March 6, 2015, 10:24 
DPMFoam  Interpretation of UcEqn.H

#1 
New Member
Jeroen Hofman
Join Date: Feb 2015
Location: The Netherlands
Posts: 15
Rep Power: 15 
Dear Foamers,
I am currently looking into the DPMFoam code and I have trouble understanding the UcEqn.H file: 1. Given the UcEqn: fvVectorMatrix UcEqn ( fvm::ddt(alphac, Uc) + fvm::div(alphaPhic, Uc)  fvm::Sp(fvc::ddt(alphac) + fvc::div(alphaPhic), Uc) + continuousPhaseTurbulence>divDevRhoReff(Uc) == (1.0/rhoc)*cloudSU ); Why is the interface momentum transfer term (cloudSU) already included? Also, I would guess the '==' is already solving the equation, but this doesn't make sense to me. 2. Directly under the definition of UcEqn: UcEqn.relax(); What does it mean to 'relax' an equation? I understand the general principle of relaxation but fail to understand how an equation statement can be relaxed. 3. Regarding solving the momentum predictor: if (pimple.momentumPredictor()) { solve ( UcEqn == fvc::reconstruct ( phicForces/rAUcf  fvc::snGrad(p)*mesh.magSf() ) ); } If the 'cloudSU' term from question 1 would already be included, it is now included twice in the solver, as phicForces also represents interphase momentum transfer. The answer to question 1 will probably also explain the answer to this question. Any help or reference to proper documentation regarding DPMFoam/MPPIC would be highly appreciated. Kind regards, Jeroen 

March 10, 2015, 09:49 
Answer question 2

#2 
New Member
Jeroen Hofman
Join Date: Feb 2015
Location: The Netherlands
Posts: 15
Rep Power: 15 
Dear Foamers,
To answer my own question regarding relaxation of the momentum equation (question 2): It is described in Jasak's thesis on Finite Volume Method, equation 3.96: In stead of relaxing the solution, which is done at a later stage in DPMFoam for the pressure, the coefficients of the momentum equation are relaxed to increase the diagonal dominance and hence the convergence rate. As far as I understood this is particularly important for steady state solutions, where the temporal derivative does not contribute to the diagonal. Hope this helps anyone struggling as well 

March 20, 2015, 07:12 
answer question 1/3

#3 
New Member
Jeroen Hofman
Join Date: Feb 2015
Location: The Netherlands
Posts: 15
Rep Power: 15 
Again to answer my own question:
DPMFoam adds a semiimplicit treatment for the Uc term, meaning that cloudSU contains:  drag (explicitly calculated, so handled as source term): drag_coefficient * (U_P  Uc_P)  pressure gradient force (if added, also a source term)  extra implicit part of the drag (drag_coefficient * Uc)  extra MINUS explicit part of the drag (drag_coefficient * Uc) The latter two sum to zero but is included to increase stability, as the drag itself is fully explicit. Notice the difference between the drag, which is calculated with the interpolated Uc value to the particle location, whereas the extra terms added contain simply the cellmid Uc value. Any further explanation on why this is so would be helpful . cloudSU has its source term part set to zero in the DPMFoam code, so eventually it only contains the third bullet from above (drag_coefficient * Uc, implicitly), this explains why it's included in the Uc equation, as it is an implicit term with Uc. cloudSUSu then contains the source part, which is drag + pressure gradient force  explicit drag part. cloudSUSu is included in the phicforces and ends up in solving the pEqn. Hope this clarifies things a bit. Jeroen 

March 24, 2015, 04:18 
Thanks！

#4 
New Member
xue
Join Date: Mar 2015
Posts: 5
Rep Power: 10 
your clarification is very helpful to me


January 19, 2016, 05:37 

#5 
Senior Member
Albrecht vBoetticher
Join Date: Aug 2010
Location: Zürich, Swizerland
Posts: 237
Rep Power: 16 
Hi Jeroen,
I am not so familiar with DPMFoam and my modified version is not stable yet, but I think the extra terms account for the numerical oscillations you may get in the pressure equation when simply adding source terms in the cell center. The PISO algorithm goes for solving U at the faces to avoid the pressurevelocity coupling problem on collocated meshes (Ferziger & Peric 2002). But as I saw you already took a look at the checkerboarding at http://www.cfdonline.com/Forums/ope...tml#post547057 

January 14, 2021, 14:28 

#6 
New Member
Mingze Gao
Join Date: Sep 2020
Posts: 7
Rep Power: 5 
Hi Jeroen,
I am also interested in how to calculate the source term. I know it is related to UTrans, but I could not find where the OpenFOAM defines the UTrans, could you please give me any clue. Thanks. Ming 

Tags 
dpmfoam 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Collision spring stiffness value in DPMFoam  maysmech  OpenFOAM Programming & Development  11  March 25, 2022 23:53 
DPMFoam  Serious Error particleladen flow in simple geometric config  benz25  OpenFOAM Running, Solving & CFD  27  December 19, 2017 21:47 
Unrealistic results using DPMFoam  Mojtaba.a  OpenFOAM Running, Solving & CFD  4  January 24, 2015 11:28 
Crash when using DPMfoam for LPT of objects contacting a vibrating wall  ansubru  OpenFOAM Running, Solving & CFD  0  May 1, 2014 04:24 
Riemann invariants....Any physical interpretation?  Farouk  Main CFD Forum  3  July 11, 2013 08:58 