CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development

modify fvMatrix / coupled equations / higher order diffusion

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Display Modes
Old   May 7, 2011, 06:20
Default modify fvMatrix / coupled equations / higher order diffusion
New Member
Jinming Lu
Join Date: May 2011
Posts: 2
Rep Power: 0
Jinming is on a distinguished road
Dear OpenFOAM community,

I've just started using OpenFOAM, trying to solve an equation of the type

df/dt = laplacian(laplacian(f))

for higher order diffusion

To my limited knowledge (please correct if I am wrong), simple diffusion

df/dt = laplacian(f) is already implemented and
fvScalarMatrix fEqn
    fvm::ddt(f) - fvm::laplacian(f)
creates an equation of the type (for the implicit case)

(F - F0)*V/deltaT - A*F

where F is a vector/list containing the cell values of f to be calculated, F0 those of the previous time step, V is the volume, deltaT is the time step, and A denotes the "discretized laplacian" (coefficient matrix).

A naive solution of

df/dt = laplacian(laplacian(f))

would be creating an equation of the type

(F - F0)*V/deltaT - A*A*F,

i.e. multiplying the fvMatrix A obtained from fvm::laplacian(f) by itself.

I would be grateful if someone could give me a hint on the implementation or how to solve the coupled problem

df/dt == laplacian(g)
laplacian(f) == g

Thanks in advance.

Best regards,

Jinming is offline   Reply With Quote


Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Higher order terms in the equations Argyle Main CFD Forum 7 October 31, 2008 08:41
2nd order equations not solving well Paul FLUENT 1 May 16, 2008 17:52
Higher order Hyperbolic systems noName Main CFD Forum 8 August 9, 2005 10:46
Higher order for phase volume fractions? George Bergantz CFX 2 December 11, 2000 12:55
Numerical diffusion error Z.Zeng Main CFD Forum 8 October 22, 1999 09:06

All times are GMT -4. The time now is 21:48.