# modify fvMatrix / coupled equations / higher order diffusion

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

 May 7, 2011, 06:20 modify fvMatrix / coupled equations / higher order diffusion #1 New Member   Jinming Lu Join Date: May 2011 Posts: 2 Rep Power: 0 Sponsored Links 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 Code: 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

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Argyle Main CFD Forum 7 October 31, 2008 08:41 Paul FLUENT 1 May 16, 2008 17:52 noName Main CFD Forum 8 August 9, 2005 10:46 George Bergantz CFX 2 December 11, 2000 12:55 Z.Zeng Main CFD Forum 8 October 22, 1999 09:06