CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Programming & Development

implicit div(k*grad(a*u))

Register Blogs Community New Posts Updated Threads Search

 
 
LinkBack Thread Tools Search this Thread Display Modes
Prev Previous Post   Next Post Next
Old   August 8, 2011, 13:45
Default implicit div(k*grad(a*u))
  #1
New Member
 
Join Date: Aug 2011
Posts: 6
Rep Power: 14
gunnar2 is on a distinguished road
Hello,
I have asked in "Running / Solving / CFD (http://www.cfd-online.com/Forums/ope...-k-grad-u.html)" before, but now this seems to be a programming question.

How can I treat \nabla\cdot k \nabla \alpha u implicitly?

k - scalar field (known diffusivity)
\alpha - scalar field (known mass fraction)
u - scalar field (unknown)

fvm::laplacian(k,alpha*u) generates a compiler error.

I think a possible solution is to right-multiply the fvMatrix fvm::laplacian(k,u) with a diagonal matrix D_\alpha made up of the values of \alpha on it's diagonal.

But how do I write this in C++?

Thank you
Gunnar
gunnar2 is offline   Reply With Quote

 


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
ddt explicit or implicit with courant number!! mrshb4 OpenFOAM 0 December 13, 2010 11:13
Help. implicit discritization dearboy Main CFD Forum 0 November 25, 2010 21:46
Can anyone explain the difference between these implicit concepts? bearcat Main CFD Forum 0 February 20, 2010 18:39
how implicit is implicit? Joe cool Main CFD Forum 7 July 16, 2004 11:10
UDF : Chemical Reaction in Catalytic Converter adhimac FLUENT 1 March 27, 2001 14:04


All times are GMT -4. The time now is 14:45.