Solve div(c1*phi)?

enoch · March 7, 2012, 12:37

I wanted to solve the following equation for T in OF. T is a scalar quantity.

div(c1*grad(T)) + div(c2*T) = div(c3)

This equation is expressed as the OF language:

fvm::laplacian(c1, T) + fvm::div(c2*T) == fvc::div(c3)

There is an compilation error on div(c2*T). Does OpenFOAM cover only two type of divergence terms: div(phi, T) and div(T)?
How can I deal with div(coeff*T)???

Thanks in advance for your comments.

lichmaster · March 8, 2012, 05:07

Hi!

First of all let's make sure I understand your problem: you want to solve a transport equation for a scalar field T with a diffusive term and a source term, right?

If so, c1 must be a scalar, but c2 and c3 must be vectors, otherwise you are trying to calculate the divergence of a scalar, which doesn't exisit.

Let's assume c2 is the vector transporting T, then this is a convection term that, if I'm not wrong, can be modeled in OpenFOAM with fvm::div(c2,T).

March 7, 2012, 12:37	*Solve div(c1phi)?**	#1
enoch Member Jeong Kim Join Date: Feb 2010 Posts: 42 Rep Power: 16	I wanted to solve the following equation for T in OF. T is a scalar quantity. div(c1grad(T)) + div(c2T) = div(c3) This equation is expressed as the OF language: fvm::laplacian(c1, T) + fvm::div(c2T) == fvc::div(c3) There is an compilation error on div(c2T). Does OpenFOAM cover only two type of divergence terms: div(phi, T) and div(T)? How can I deal with div(coeff*T)??? Thanks in advance for your comments.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Solve Flow or VOF simultaneously ?	Ramsey	FLUENT	1	February 16, 2011 14:16
Linearized NS euqations: how to solve them?(problem with Matrix operations..)	matteoL	OpenFOAM Running, Solving & CFD	0	November 18, 2009 07:58
Solve for two or more "Temperatures"	Rui	CFX	12	September 9, 2008 22:58
How solve UDS explicitely?	Hesam	FLUENT	0	July 1, 2007 10:03
How to solve another continuum and momentum eqn?	west_wing	FLUENT	0	August 25, 2003 11:00

March 8, 2012, 05:07		#2
lichmaster New Member Lorenzo A. Ricciardi Join Date: Mar 2011 Posts: 27 Rep Power: 15	Hi! First of all let's make sure I understand your problem: you want to solve a transport equation for a scalar field T with a diffusive term and a source term, right? If so, c1 must be a scalar, but c2 and c3 must be vectors, otherwise you are trying to calculate the divergence of a scalar, which doesn't exisit. Let's assume c2 is the vector transporting T, then this is a convection term that, if I'm not wrong, can be modeled in OpenFOAM with fvm::div(c2,T).