Hi all, I'm totally new to FiPy (and finite volume methods generally) so I apologize if this is a silly question. I need to solve a semilinear wave equation of the form (in 1 dimension):
utt - uxx = f(x, u) where the notation is: u unknown, ux is the first space derivative, utt the second time derivative etc. f(x, u) is a nonlinear function of both x and u which can vary quite rapidly with both variables. I'd also like to solve the analogous problem in 2-D and 3-D and it would be really good to be able to solve it on arbitrary domains. It's easy to get this into first order form, ut = v vt = uxx + f(x, u) But it's still not quite in the standard form for FiPy. I'm guessing I'd have to write the uxx term as a source term? Is FiPy suited to this kind of problem - and if not can you recommend another (preferably open source) tool? I'd appreciate any advice. Cheers, Michael Brown
_______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
