Solving stationary two-phase Laplace-equation

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

 June 3, 2009, 18:39 Solving stationary two-phase Laplace-equation #1 New Member   Martin Tek Join Date: May 2009 Posts: 7 Rep Power: 8 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?

 June 5, 2009, 11:51 #2 New Member   Martin Tek Join Date: May 2009 Posts: 7 Rep Power: 8 OMG I am really lost. What I tried: ... volVectorField gradepsilon = fvc::grad(epsilon); volVectorField gradU = fvc::grad(U); solve ( 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.

 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 Niklas Wikstrom (Wikstrom) OpenFOAM Running, Solving & CFD 122 June 15, 2014 06:20 xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33 carsten OpenFOAM Bugs 11 September 12, 2008 11:16 msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58 liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07

All times are GMT -4. The time now is 10:58.