Hello, I need to solve a non-linear PDE of this type:
******************** eq = DiffusionTerm(coeff = numerix.sqrt(PHI)) == Gamma/numerix.sqrt(phi) ******************** where Gamma is a constant. and the mesh is 2D. the problem is de 1/numerix.sqrt(phi), which gives me, depending on baoundary conditions and system size, a lot of divergence problems. is there any way to put it in the other side of the equation, so it is not a division? I do not know how would I express the diffusion term in that case Thanks a lot for your help -- Dr. Francisco Vega Reyes Departamento de Fisica, Universidad de Extremadura _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
