On Mon, Feb 9, 2009 at 5:33 PM, Daniel Wheeler <[email protected]> wrote: > > Hi Ionut, > > I'm not sure why your script is not converging. However, I do have > some queries below. >>> fSecondDerivative = -exp(-faceVar) + (3 * b / faceVar**4) >>> DiffusionTerm(coeff=fSecondDerivative) > > Have you tried the above? This seems like the most sensible approach > for the second order part of the equation.
yes, I tried also this possibility, but I have the same result as before. >> fDerivative = exp(-var)-(b / var**3) >> #fSecondDerivative = -exp(-faceVar) + (3 * b / faceVar**4) >> >> diffTerm4 = - DiffusionTerm((h**3, 1)) >> eq = TransientTerm() == -DiffusionTerm((h**3, 1)) + >> (fDerivative.getFaceGrad() * h**3).getDivergence() >> #eq = TransientTerm() == -DiffusionTerm((h**3, 1)) + >> DiffusionTerm(coeff=fSecondDerivative) > > What happens when you uncomment this? Actually, what is h? Shouldn't > "h" be "var"? If I uncomment and I am using the second derivative, I have the same error/result as before. Thanks, Ionut
