How to solve this coefficient form PDE in OpenFOAM?

lupe123 · January 16, 2020, 03:47

Hi all,

I am new to working with OpenFOAM and am trying to recreate a COMSOL solution for a coefficient form PDE with OpenFOAM.

The equation is:

Code:

div(gamma-grad(u))==0

where u is a scalar that depends on x,y,z and gamma is a constant vector.

I have looked through the existing solvers in OpenFOAM, but could not find one that fits this equation.

So my question is: Can this equation be solved in OpenFOAM with an existing solver or would I have to implement a new solver for this?

Regards,
Lukas

Santiago · January 16, 2020, 06:48

Quote:

Originally Posted by lupe123

Hi all,

I am new to working with OpenFOAM and am trying to recreate a COMSOL solution for a coefficient form PDE with OpenFOAM.

The equation is:

Code:

div(gamma-grad(u))==0

where u is a scalar that depends on x,y,z and gamma is a constant vector.

I have looked through the existing solvers in OpenFOAM, but could not find one that fits this equation.

So my question is: Can this equation be solved in OpenFOAM with an existing solver or would I have to implement a new solver for this?

Regards,
Lukas

If gamma is a constant vector, Wouldn't your PDE reduce just to -div(grad(u)) = 0?
Or did you meant gamma = gamma(x,y,z) constant?

lupe123 · January 16, 2020, 08:16

Hi Santiago,

Thanks for the reply.

I meant gamma = gamma(x,y,z) constant. In COMSOL this is called conservative flux source.

January 16, 2020, 03:47	How to solve this coefficient form PDE in OpenFOAM?	#1
lupe123 New Member Join Date: Jan 2020 Posts: 2 Rep Power: 0	Hi all, I am new to working with OpenFOAM and am trying to recreate a COMSOL solution for a coefficient form PDE with OpenFOAM. The equation is: Code: div(gamma-grad(u))==0 where u is a scalar that depends on x,y,z and gamma is a constant vector. I have looked through the existing solvers in OpenFOAM, but could not find one that fits this equation. So my question is: Can this equation be solved in OpenFOAM with an existing solver or would I have to implement a new solver for this? Regards, Lukas

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
OpenFOAM v3.0.1 Training, London, Houston, Berlin, Jan-Mar 2016	cfd.direct	OpenFOAM Announcements from Other Sources	0	January 5, 2016 04:18
OpenFOAM Training, London, Chicago, Munich, Sep-Oct 2015	cfd.direct	OpenFOAM Announcements from Other Sources	2	August 31, 2015 14:36
Drag Coefficient for Ellipse Form (2D)	Ketut Utama	Main CFD Forum	8	December 11, 2014 12:03
Modified OpenFOAM Forum Structure and New Mailing-List	pete	Site News & Announcements	0	June 29, 2009 06:56
Can I use OpenFOAM to solve unsteady diffusion problem	yongshenglian	OpenFOAM Running, Solving & CFD	1	September 17, 2008 13:03

January 16, 2020, 08:16		#3
lupe123 New Member Join Date: Jan 2020 Posts: 2 Rep Power: 0	Hi Santiago, Thanks for the reply. I meant gamma = gamma(x,y,z) constant. In COMSOL this is called conservative flux source.