# Trouble understanding the terms of UEqns.H in twoPhaseEulerFoam

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 February 12, 2019, 18:49 Trouble understanding the terms of UEqns.H in twoPhaseEulerFoam #1 Member   Join Date: Apr 2016 Posts: 30 Rep Power: 10 Hi all, I am trying to add variable source terms in the twoPhaseEulerFoam. However I am unable to understand the various terms of the momentum equation. I am trying to add two force terms fsid and fl in the momentum equation. However don't know where to add the terms. Could somebody explain what the bold terms physically mean? I also cannot find the laplacian term for viscous stresses. Thanks and regards, Shantanu Code: ```//#include "Random.H" //#include "StochasticDispersionRAS.H" #include "contErrs.H" Info<< "Constructing momentum equations" << endl; MRF.correctBoundaryVelocity(U1); MRF.correctBoundaryVelocity(U2); MRF.correctBoundaryVelocity(U); fvVectorMatrix U1Eqn(U1, rho1.dimensions()*U1.dimensions()*dimVol/dimTime); fvVectorMatrix U2Eqn(U2, rho2.dimensions()*U2.dimensions()*dimVol/dimTime); dimensionedScalar dimen ( "dimen", dimensionSet(0,0,1,0,0,0,0), scalar(1.0) ); dimensionedScalar dimen1 ( "dimen1", dimensionSet(0,0,-1,0,0,0,0), scalar(1.0) ); dimensionedScalar dimen2 ( "dimen1", dimensionSet(0,0,-1,0,0,0,0), scalar(1.0) ); volVectorField gr_alpha = fvc::grad(alpha2); volTensorField grU_2 = (fvc::grad(U2))*dimen; volScalarField sc = grU_2 && (grU_2 + grU_2.T()); const scalar csid = 25; const scalar a = 4e-6; //const scalar k = .1; const scalar eps = 0.206; //Lift force //dimensionedScalar up = ("up", [0 0 -1 0 0 0 0], 1.0); //const scalar up =1.0; volScalarField up = .1629+ (.2088/(.1476+(Foam::sqrt(sc)))); //volScalarField up = /*1.0*sc/(sc+ROOTSMALL)*/ sc*dimen1; const scalar cl = 1.2; volVectorField fsid = U2*dimen1; // copy of the field; volVectorField fl = U2*dimen1; volScalarField dis = wallDist(mesh).y(); volVectorField ncap = wallDist(mesh).n(); //Info<<"data="<owner().rndGen(); scalar rn = rnm.GaussNormal(); //Info< wpd; //scalar dis = wpd.updateCell(); //wallDist wpd; // volScalarField dis = wallDist(mesh).y(); // Info<<"data="<

 November 19, 2021, 10:27 #2 New Member   Mehrdad Join Date: May 2020 Posts: 1 Rep Power: 0 Did you manage to understand the equations? where's the viscous stress term?

 Tags fvm, momentum, source terms, twophaseeulerfoam, ueqn.h

 Thread Tools Search this Thread Search this Thread: Advanced Search Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are Off Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post AMK53 OpenFOAM Running, Solving & CFD 5 September 18, 2019 01:50 ARUN K RAJ OpenFOAM Programming & Development 1 February 19, 2015 05:43 JMDag2004 OpenFOAM Running, Solving & CFD 2 January 16, 2015 08:10 titio OpenFOAM Running, Solving & CFD 0 March 19, 2009 16:02 Brindaban Ghosh Main CFD Forum 2 June 24, 2000 04:22

All times are GMT -4. The time now is 23:08.

 Contact Us - CFD Online - Privacy Statement - Top