Implementation of advection in alphaEqn.H in settlingFoam

Navip · July 17, 2015, 01:34

Dear Foamers,

I wish to add an advection term (div(U*Alpha)) in the settlingFoam solver, alphaEqn.H.

I try different option but no one works. Does anyone of you have an idea how to do that ? The problem lies on the fact that U and Alpha are matrix I think.

volScalarField divU(fvc::div(phi/fvc::interpolate(rho)));

fvScalarMatrix AlphaEqn
(
fvm::ddt(rho, Alpha)
+ fvm::div(phiAlpha, Alpha)
+ fvm::div(divU,Alpha) //I tried this as divU is showned just up
- fvm::laplacian(mut, Alpha)
);

Help would be appreciate,

Xavier

rapierrz · July 17, 2015, 03:09

Hi Xavier,

I have the same problem with a small difference.when I use

phi/fvc::interpolate(rho), result of alpha become with high gradients in some

locations and I don't know the reason.

But in your case i think if your solver even make succesfully it has dimentional

problem in divU.

July 17, 2015, 01:34	Implementation of advection in alphaEqn.H in settlingFoam	#1
Navip New Member Xavier Pivan Join Date: May 2015 Posts: 10 Rep Power: 12	Dear Foamers, I wish to add an advection term (div(UAlpha)) in the settlingFoam solver, alphaEqn.H. I try different option but no one works. Does anyone of you have an idea how to do that ? The problem lies on the fact that U and Alpha are matrix I think. volScalarField divU(fvc::div(phi/fvc::interpolate(rho))); fvScalarMatrix AlphaEqn ( fvm::ddt(rho, Alpha) + fvm::div(phiAlpha, Alpha) + fvm::div(divU,Alpha) //I tried this as divU is showned just up* - fvm::laplacian(mut, Alpha) ); Help would be appreciate, Xavier

Thread Tools	Search this Thread
Show Printable Version Email this Page	Search this Thread: Advanced Search
Display Modes
Linear Mode Switch to Hybrid Mode Switch to Threaded Mode

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
alphaEqn.H in twoPhaseEulerFoam	cheng1988sjtu	OpenFOAM Bugs	15	May 1, 2016 17:12
ENO advection implementation	mpash	Main CFD Forum	0	January 24, 2014 00:20
Use advection step and Dispersion step to solve advection diffusion equation	cooljd	Main CFD Forum	0	October 15, 2012 16:44
Total variation diminishing scheme for advection or advection dispersion equation	cooljd	Main CFD Forum	7	July 27, 2012 09:32
Implementation issues of fvSchemes / laplacianScheme, in particular gaussLaplacianSch	thomek	OpenFOAM Programming & Development	0	October 18, 2010 06:10

July 17, 2015, 03:09		#2
rapierrz Senior Member Hesam Join Date: Feb 2015 Posts: 139 Rep Power: 12	Hi Xavier, I have the same problem with a small difference.when I use phi/fvc::interpolate(rho), result of alpha become with high gradients in some locations and I don't know the reason. But in your case i think if your solver even make succesfully it has dimentional problem in divU.