
[Sponsors] 
September 7, 2011, 08:27 

#21 
Senior Member

Thank you Alberto, I also found it and now trying to understand this expression.
Please see attached file Volumefield_reconstruction.pdf. I am looking for the original reference where this formulation is established. Could you please give some tips?
__________________
Best regards, Dr. Alexander VAKHRUSHEV Christian Doppler Laboratory for "Advanced Process Simulation of Solidification and Melting" Simulation and Modelling of Metallurgical Processes Department of Metallurgy University of Leoben FranzJosefStr. 18 A  8700 Leoben Österreich / Austria Tel.: +43 3842  402  3125 http://smmp.unileoben.ac.at 

September 7, 2011, 21:16 

#22 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,910
Rep Power: 28 
I do not have a reference for this, but I think you got it.
__________________
Alberto Passalacqua GeekoCFD  A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats. OpenQBMM  An opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

September 8, 2011, 03:37 

#23 
Senior Member
Anton Kidess
Join Date: May 2009
Location: Germany
Posts: 1,215
Rep Power: 23 
Alexander, I think it would be nice to have this on the wiki. If you agree, do you have the time to copy it to an article there?


September 8, 2011, 03:39 

#24 
Senior Member

Hi Anton!
Ok, I will put it there. And if Alberto doesn't mind, it would look better with his description of the body forces and the pressure gradient treatment technique.
__________________
Best regards, Dr. Alexander VAKHRUSHEV Christian Doppler Laboratory for "Advanced Process Simulation of Solidification and Melting" Simulation and Modelling of Metallurgical Processes Department of Metallurgy University of Leoben FranzJosefStr. 18 A  8700 Leoben Österreich / Austria Tel.: +43 3842  402  3125 http://smmp.unileoben.ac.at Last edited by makaveli_lcf; September 8, 2011 at 03:43. Reason: Additions 

September 13, 2011, 01:20 

#25  
Senior Member
Santiago Marquez Damian
Join Date: Aug 2009
Location: Santa Fe, Santa Fe, Argentina
Posts: 430
Rep Power: 16 
Alexander,
Quote:
Regards.
__________________
Santiago MÁRQUEZ DAMIÁN, Ph.D. Research Scientist Research Center for Computational Methods (CIMEC)  CONICET/UNL Tel: 543424511594 Int. 7032 Colectora Ruta Nac. 168 / Paraje El Pozo (3000) Santa Fe  Argentina. http://www.cimec.org.ar 

July 26, 2012, 04:28 
Could you be so kind to write down the corresponding code for the equation for p?

#26  
New Member
Join Date: Nov 2010
Posts: 2
Rep Power: 0 
Dear Alberto,
thank you very much for this very detailed explanation about the issue how to add correctly a body force to the (incompressible) momentum equation! I would be grateful if you would still specify the corresponding code for the p equation. Best regards, Vladimir Quote:


July 29, 2012, 15:20 

#27  
Member
Join Date: May 2012
Location: Dresden, Germany
Posts: 32
Rep Power: 7 
Hallo Vladimir,
I don't understand everything yet, but based on interFoam it should be sth like that: Quote:
Quote:
Working on pisoFoam, I think rho has to be removed. I do not know if phi has to contain only the velocity flux, or the 'force flux', too. Any idea? For my understanding: Can someone give me a hint, what 'reconstruct' does? It 'adds' terms to the equation and changes A and H? But there must be a difference between Quote:
Quote:


July 30, 2012, 01:33 
plz help me

#28 
Member
vahid
Join Date: Feb 2012
Location: MashhadIran
Posts: 80
Rep Power: 6 
Hi Dear alberto again:
I want to add surface tension(sigma) in one solver,for this reason I added : please help me,and tell me ,How to correct this problem???#include ''fvCFD.H'' fvc::interpolate(interface.sigma()) in this code: Foam::tmp<Foam::volScalarField> Foam:haseChangeTwoPhaseMixtures::SchnerrSauer: Coeff ( const volScalarField& p ) const { volScalarField limitedAlpha1(min(max(alpha1_, scalar(0)), scalar(1))); volScalarField rho ( limitedAlpha1*rho1() + (scalar(1)  limitedAlpha1)*rho2() ); return //......I want to change it( <<sigma>> surface tension multiple in it): (3*rho1()*rho2())*sqrt(2/(3*rho1()))*(fvc::interpolate(interface.sigma())) *rRb(limitedAlpha1)/(rho*sqrt(mag(p  pSat()) + 0.01*pSat())); //.................................................. ...... } dont successful wmake, and seen(was not declared ): phaseChangeTwoPhaseMixtures/SchnerrSauer/SchnerrSauer.C:113: error: 'interface' was not declared in this scope make: *** [Make/linux64GccDPOpt/SchnerrSauer.o] Error 1 

June 4, 2014, 10:51 

#29 
Senior Member
M. Montero
Join Date: Mar 2009
Location: Madrid
Posts: 123
Rep Power: 10 
HI all,
How could it be done in case it is not a body force but a vectorField? All is in relation to the rotorDiskSource and how introducing this source in the momentum equation can create wiggles in the pressure field and velocity field. Thank you Last edited by be_inspired; June 6, 2014 at 07:14. 

February 23, 2015, 10:16 

#30  
New Member
Peter
Join Date: May 2012
Location: New York
Posts: 18
Rep Power: 7 
Quote:
Hello Harry, Did you finally solve your problem? I have met a similar problem regarding the bodyforceweighted scheme for pressure interpolation. Do you have any ideas on this? Thanks a lot for your help and time! Best, Peter 

May 12, 2017, 09:26 

#31 
New Member
Roman G.
Join Date: Apr 2017
Posts: 9
Rep Power: 2 
Hi,
I know this treat is old but I'm struggling with the implementation of the force term by the use of the RhiChow interpolation. I'm using pimpleFoam. So far I added in the momentum predictor my force term like this: Code:
surfaceScalarField B_fs = fvc::interpolate(lorentz/rho) & mesh.Sf() if (pimple.momentumPredictor()) { solve ( UEqn == fvc::reconstruct (  fvc::snGrad(p)* mesh.magSf() + B_fs ) ); fvOptions.correct(U); } I added a flux term in the phiHbyA equation: Code:
surfaceScalarField rAUf("rAUf", fvc::interpolate(rAU)); surfaceScalarField phil(rAUf*B_fs); surfaceScalarField phiHbyA ( "phiHbyA", fvc::flux(HbyA) + fvc::interpolate(rAU)*fvc::ddtCorr(U, phi) + phil ); Code:
U = HbyA + rAU*fvc::reconstruct((phil  pEqn.flux())/rAUf); Code:
Continuity error cannot be removed by adjusting the outflow. Please check the velocity boundary conditions and/or run potentialFoam to initialise the outflow. Total flux : 8.226416123979263e05 Specified mass inflow : 0.0005613679754126387 Specified mass outflow : 6.154738583391114e17 Adjustable mass outflow : 0 Thanks in advance 

May 18, 2017, 05:15 

#32 
New Member
Roman G.
Join Date: Apr 2017
Posts: 9
Rep Power: 2 
The problem is solved. I had to add the flux from the force after adjusting phi.
Code:
surfaceScalarField phil(rAUf*B_fs); surfaceScalarField phiHbyA ( "phiHbyA", fvc::flux(HbyA) + fvc::interpolate(rAU)*fvc::ddtCorr(U, phi) ); MRF.makeRelative(phiHbyA); adjustPhi(phiHbyA, U, p); phiHbyA += phil; 

Tags 
body force, force 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
urgent help needed (rhiechow interpolation problem)  Ardalan  Main CFD Forum  2  March 18, 2011 16:22 
Help on 2D interpolation in StarCCM+  madhuri  CDadapco  0  November 3, 2010 15:21 
Surface interpolation schemes and parallelization  jutta  OpenFOAM Running, Solving & CFD  0  February 25, 2010 15:32 
momentum interpolation for collocated grid  Hadian  Main CFD Forum  4  December 25, 2009 08:25 
spline interpolation  bajjal  Main CFD Forum  0  May 29, 2006 08:27 