How to add this term to openfoam solver

rapierrz · February 25, 2015, 14:59

Hello friends,

how can i add the right hand side equation to icoFoam???

kmooney · February 26, 2015, 20:05

If this is the NS diffusion term its already there.

rapierrz · February 27, 2015, 02:16

yes.This is diffusion term but usually in solvers there is only

div(grad(U)) and i want to know how to write

div(grad(U)+grad(U).tanspose)

I write div(nu,symm(fvc::grad(U))) and Error was :

no matching function for call to ‘div(Foam::volScalarField&...

ssss · February 27, 2015, 04:52

I suppose you are using non constant viscosity so:

Code:

fvc::grad(mu)&fvc::grad(U)+fvc::laplacian(mu,U)

You should change fvc for fvm if you want to use a more implicit calculation in the momentumPredictor.

I recommend yo to take a look at the solver: nonNewtonianIcoFoam

rapierrz · February 27, 2015, 06:20

thank you ssss,

my problem is that how i can define div(mu,grad(U).transpose)

ssss · February 27, 2015, 06:23

Code:

fvc::grad(U).T()

rapierrz · February 27, 2015, 06:35

When i use of fvc::grad(U).T() error is :

has no member named ‘T’

ssss · February 27, 2015, 06:37

Then try:

Code:

T(fvc::grad(U))

rapierrz · February 27, 2015, 06:47

it does not work

my openfoam version is 2.2.1

ssss · February 27, 2015, 09:13

Well then last chance:

Code:

fvc::grad(U)().T()

It compiles for me

rapierrz · February 28, 2015, 13:30

Thank you very much ssss

ut's working correctly.

you was great help for me.

February 25, 2015, 14:59	How can i add this term to icoFoam	#1
rapierrz Senior Member Hesam Join Date: Feb 2015 Posts: 139 Rep Power: 11	Hello friends, how can i add the right hand side equation to icoFoam??? Last edited by rapierrz; February 26, 2015 at 05:11.

February 27, 2015, 04:52		#4
ssss Senior Member anonymous Join Date: Aug 2014 Posts: 205 Rep Power: 13	I suppose you are using non constant viscosity so: Code: fvc::grad(mu)&fvc::grad(U)+fvc::laplacian(mu,U) You should change fvc for fvm if you want to use a more implicit calculation in the momentumPredictor. I recommend yo to take a look at the solver: nonNewtonianIcoFoam

February 27, 2015, 06:23		#6
ssss Senior Member anonymous Join Date: Aug 2014 Posts: 205 Rep Power: 13	Code: fvc::grad(U).T()

February 27, 2015, 06:37		#8
ssss Senior Member anonymous Join Date: Aug 2014 Posts: 205 Rep Power: 13	Then try: Code: T(fvc::grad(U))

February 27, 2015, 09:13		#10
ssss Senior Member anonymous Join Date: Aug 2014 Posts: 205 Rep Power: 13	Well then last chance: Code: fvc::grad(U)().T() It compiles for me rapierrz likes this.

February 26, 2015, 20:05		#2
kmooney Senior Member Kyle Mooney Join Date: Jul 2009 Location: San Francisco, CA USA Posts: 323 Rep Power: 18	If this is the NS diffusion term its already there.

February 27, 2015, 02:16		#3
rapierrz Senior Member Hesam Join Date: Feb 2015 Posts: 139 Rep Power: 11	yes.This is diffusion term but usually in solvers there is only div(grad(U)) and i want to know how to write div(grad(U)+grad(U).tanspose) I write div(nu,symm(fvc::grad(U))) and Error was : no matching function for call to ‘div(Foam::volScalarField&...

February 27, 2015, 06:20		#5
rapierrz Senior Member Hesam Join Date: Feb 2015 Posts: 139 Rep Power: 11	thank you ssss, my problem is that how i can define div(mu,grad(U).transpose)

February 27, 2015, 06:35		#7
rapierrz Senior Member Hesam Join Date: Feb 2015 Posts: 139 Rep Power: 11	When i use of fvc::grad(U).T() error is : has no member named ‘T’

February 27, 2015, 06:47		#9
rapierrz Senior Member Hesam Join Date: Feb 2015 Posts: 139 Rep Power: 11	it does not work my openfoam version is 2.2.1

February 28, 2015, 13:30		#11
rapierrz Senior Member Hesam Join Date: Feb 2015 Posts: 139 Rep Power: 11	Thank you very much ssss ut's working correctly. you was great help for me.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Solve poisson equation just add a source term	nandiganavishal	OpenFOAM Running, Solving & CFD	18	November 14, 2022 10:12
Can't get data from OpenFoam to external solver using externalCoupled	perry	OpenFOAM Running, Solving & CFD	4	May 26, 2014 09:09
Building a custom solver on OpenFOAM 2.0	wschosta	OpenFOAM Programming & Development	1	July 8, 2011 16:07
Add time dependent Source term	libe	OpenFOAM Running, Solving & CFD	13	September 23, 2010 10:12
OpenFOAM Training in Europe and USA	hjasak	OpenFOAM	0	August 8, 2008 06:33