Hi, I am trying to solve the following nonlinear conservation law in one spatial dimension in the variable u(t,x):
u_x + f(u)_t = 0 with initial condition u(.,0) = u0. Here f is a nonlinear function of u (f is a local operator). I want to solve this on the domain [-tmin, tmax] x [0,L] (this is just a conservation law with variables x,t interchanged...), so the mesh is on the interval [-tmin, tmax]. So far I am using periodic boundary conditions for simplicity. If f were of the form f(u) = f_0(u)*u I could simply use a Convection term and update the coefficient f_0(0) at each step, I have done that successfully. But for the above problem I am not sure how to employ the .divergence or getGradient() methods and I am struggling to find examples. Can someone please help? I am new to FiPy, so this might be an easy thing to do, but I haven't found a solution in the manual / FAQ. Thanks, Max _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
