1 Attachment(s)
Thank you Alberto, I also found it and now trying to understand this expression.
Please see attached file Attachment 9124. I am looking for the original reference where this formulation is established. Could you please give some tips? |
I do not have a reference for this, but I think you got it.
|
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?
|
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. |
Alexander,
Quote:
Regards. |
Could you be so kind to write down the corresponding code for the equation for p?
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:
|
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:
|
plz help me
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 |
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 |
Quote:
Hello Harry, Did you finally solve your problem? I have met a similar problem regarding the body-force-weighted scheme for pressure interpolation. Do you have any ideas on this? Thanks a lot for your help and time! Best, Peter |
Hi,
I know this treat is old but I'm struggling with the implementation of the force term by the use of the Rhi-Chow 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() I added a flux term in the phiHbyA equation: Code:
surfaceScalarField rAUf("rAUf", fvc::interpolate(rAU)); Code:
U = HbyA + rAU*fvc::reconstruct((phil - pEqn.flux())/rAUf); Code:
Continuity error cannot be removed by adjusting the outflow. Thanks in advance |
The problem is solved. I had to add the flux from the force after adjusting phi.
Code:
surfaceScalarField phil(rAUf*B_fs); |
Quote:
Have you also modified the line Code:
// Calculate the conservative fluxes |
Hi gaza,
I had problems when adding a source term without adjusting the flux to ensure mass conservation. Try to add the source terms after adjusting phi to your flux phiHybA. Quote:
|
Hi Novel,
Thank you for your tip however in my case I do not see the difference after adjustPhi. |
All times are GMT -4. The time now is 04:05. |