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

Boundary condition for fourth order differential equation

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

Like Tree3Likes
  • 1 Post By K.C.
  • 2 Post By K.C.

LinkBack Thread Tools Search this Thread Display Modes
Old   March 23, 2015, 11:26
Default Boundary condition for fourth order differential equation
Join Date: Mar 2015
Posts: 36
Rep Power: 11
K.C. is on a distinguished road
Dear Foamers,

actually i am a bit confused of OpenFoam's implementation of boundary conditions.
Actually I have to solve the followong pDGL for my concentration profile:

fvm::ddt(c) == fvm::laplacian(D, "laplacian(c)") (+ other terms)
with D as diffusion coefficient in form of a volScalarField

In this pDGL fourth order derivative with respect to space (and some mixed derivatives) are present, but OpenFoan uses only one boundary condition for the c-Field.
How can this work and is there an easy way to implement the second boundary condition?
Or more detailed: What happens in the fvm::laplacian( ... "laplacian")?
The best way would be to output the CEqn in matrix-form to an external file like the c-Field. Has anybody of you ever done this?

(Btw.: An explicit treatment is already working like A = fvc::laplacian(C) and fvm::ddt(c) == fvc::laplacian(D,A) with seperated BC to c- and A-field.
The question is about the implicit treatment to use Crank-Nicolson in the end and become second order acurate in time.)

Thanks for any comments.
tom_flint2012 likes this.
K.C. is offline   Reply With Quote

Old   May 12, 2015, 09:25
Join Date: Mar 2015
Posts: 36
Rep Power: 11
K.C. is on a distinguished road
I think, I "solved" the problem on my own.

The matrix structure of OpenFoam (LDU-matrix) does not allow to use the neighbour of a neighbour for discretisation. Only the direct neighbour can be used, so the normal "laplace"-operator is the highest possible implicit derivative.

Therefore a higher implicit derivative is only possible by rewriting all matrix classes and the used solver.
tom_flint2012 and jobin2600 like this.
K.C. is offline   Reply With Quote


boundary-conditions, crank-nicolson, fourth order pdgl

Thread Tools Search this Thread
Search this Thread:

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

Similar Threads
Thread Thread Starter Forum Replies Last Post
mixed inflow/outflow downstream boundary condition question peob OpenFOAM Running, Solving & CFD 3 February 3, 2017 10:54
Question about heat transfer coefficient setting for CFX Anna Tian CFX 1 June 16, 2013 06:28
CFX Solver : Sudden crash Hervé CFX 2 June 16, 2008 06:40
boundary condition for species equation zhou1 FLUENT 1 June 3, 2004 16:15
Boundary Conditions Jan Ramboer Main CFD Forum 11 August 16, 1999 08:59

All times are GMT -4. The time now is 20:38.