[Saleh Abuhanieh]

Hi,


I am trying to solve a simple problem (it has analytical solution) using OpenFOAM:


d/dx[ 2x*y(x) + (epsilon + x^2)*dy(x)/dx] = 0

The boundary values of y are given

I've converted the equation to OpenFOAM as follow:




fvm::div(phi, y) + fvm::laplacian(D, y)


where

phi -> linearInterpolate(U) & mesh.Sf()
U -> the evaluation of 2x multiplied by x unit vector
D -> (epsilon + x^2)


The result is zero everywhere and y doesn't change at all (zero iteration)




Any help in highlighting my mistake(s) will be appreciated


Saleh