Form a couple equation eqn = eqn_ux & eon_uy
and then call eqn.sweep(...) on that. With coupled equations, you don't pass a variable to sweep. > On Oct 20, 2016, at 6:56 AM, Francisco Vega Reyes <[email protected]> wrote: > > Hello, > > Is it possible to express in a simple way the sweep method for more than > one variable? > > More concretely, in a similar problem to the Stokes Cavity example I got > an equation with DiffusionTerm both for the X velocity and the Y > velocity (in the same equation, due to tensorial viscosity). > > However, following the SIMPLE algorithm procedure in the example, the > sweep method is used, and this only seems to accept one variable per > equation, > > How could I solve the Stokes Cavity an equation (actually, two) of this > form?: > > mux = numerix.array( ((2.0, 0.), (0.0, 1.0)) ) ; > muy = numerix.array( ((1.0, 0.), (0.0, 2.0)) ) ; > muyx = numerix.array( ((0., 0.), (1.0, 0.)) ) ; > muxy = numerix.array( ((0., 1.), (0., 0.)) ) ; > > eqn_ux = DiffusionTerm(coeff=mux,var=ux) + DiffusionTerm(coeff=muyx , > var =uy) - p.grad.dot([1.,0.]) > > eqn_uy = DiffusionTerm(coeff=muy, var=uy) + DiffusionTerm(coeff=muxy , > var =ux) - p.grad.dot([0.,1.]) > > thanks a lot > > Francisco Vega > > _______________________________________________ > fipy mailing list > [email protected] > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] _______________________________________________ fipy mailing list [email protected] http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
