Solving stationary two-phase Laplace-equation
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
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?
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.
|All times are GMT -4. The time now is 08:00.|