CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

Solving stationary two-phase Laplace-equation

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

LinkBack Thread Tools Display Modes
Old   June 3, 2009, 18:39
Default Solving stationary two-phase Laplace-equation
New Member
Martin Tek
Join Date: May 2009
Posts: 7
Rep Power: 8
Martin123 is on a distinguished road
Sorry first for my beginner question, I am quite new to CFD. I would like to solve the Laplace equation for inhomogeneous media:

k div grad U + grad k grad U = 0

For starting, I tried to solve the first term only and managed to setup the k-field (volscalarfield) on run-time. However, I don't see any influence of varying k when using fvm:: discretization - I solved

fvm::laplacian(k, U)

Am I on the right way or should I rethink completely my approach? I am very unsure how OF deals with the extra volScalarField from K: I want to solve certainly for U but I wonder what OF would see with my equation above? Should I better use somehow dictionary values for k and grad k?
Martin123 is offline   Reply With Quote

Old   June 5, 2009, 11:51
New Member
Martin Tek
Join Date: May 2009
Posts: 7
Rep Power: 8
Martin123 is on a distinguished road
OMG I am really lost. What I tried:

volVectorField gradepsilon = fvc::grad(epsilon);
volVectorField gradU = fvc::grad(U);

fvm::laplacian(U) + ((gradepsilon & gradU)/epsilon)

But still, I don't get correct results. Is my assumption right, that I may not use a second geometric field here inside since it is considered as an additional degree of freedom? Did I understand correctly that: if I use fvc:: this is added as a source term whereas all fvm:: operations cause the DOF (e.g. U) to be considered?

Is it alternatively possible to construct a dictionary from a volScalarField? Or taking dictionary values from a run-time function?

Many thanks in advance for your help.
Martin123 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
Moving mesh Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 06:20
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
Differences between serial and parallel runs carsten OpenFOAM Bugs 11 September 12, 2008 11:16
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07

All times are GMT -4. The time now is 00:59.