Have not looked at your code. But a quick note -- Burger's equation is prone to forming shocks, whose width is set by the diffusive terms. In numerical solutions, the solution then is often determined by numerical diffusivity, which is non-optimal. Play around with analytical solutions with shocks and make sure your code is doing something sensible. It is a hard one.
Jamie On Mon, Aug 29, 2016 at 9:09 AM, Sebastian Glane <gl...@tu-berlin.de> wrote: > Dear fipy-users and developers, > > I am new to fipy and would like to solve Burgers' equation as a first > step. This problem involves a non-linear convection term. According to > the FAQs, for the diffusion term it is to program the nonlinearity using > coeff. I did the same for the convection term: > fipy.VanLeerConvectionTerm( coeff = (0.5 * phi, ), var = phi ). But > somehow it does not work out. I also tried with phi.faceValue. > > Best regards > > Sebastian >
_______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]