/[escript]/trunk/modellib/py_src/flow.py
ViewVC logotype

Diff of /trunk/modellib/py_src/flow.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 2344 by jfenwick, Mon Mar 30 02:13:58 2009 UTC revision 2474 by gross, Tue Jun 16 06:32:15 2009 UTC
# Line 59  class SteadyIncompressibleFlow(Model): Line 59  class SteadyIncompressibleFlow(Model):
59             self.__p_very_old=None             self.__p_very_old=None
60             self.__dt_old=None             self.__dt_old=None
61             self.__pde=LameEquation(self.domain)             self.__pde=LameEquation(self.domain)
62             self.__pde.setSolverMethod(LameEquation.DIRECT)             self.__pde.getSolverOptions().setSolverMethod(self.__pde.getSolverOptions().DIRECT)
63             if self.location_prescribed_velocity == None: self.location_prescribed_velocit=Data()             if self.location_prescribed_velocity == None: self.location_prescribed_velocit=Data()
64             if self.prescribed_velocity == None: self.prescribed_velocity=Data()             if self.prescribed_velocity == None: self.prescribed_velocity=Data()
65    
# Line 104  class SteadyIncompressibleFlow(Model): Line 104  class SteadyIncompressibleFlow(Model):
104                                sigma=self.pressure*kronecker(self.__pde.getDomain()), \                                sigma=self.pressure*kronecker(self.__pde.getDomain()), \
105                                r=self.prescribed_velocity, \                                r=self.prescribed_velocity, \
106                                q=self.location_prescribed_velocity)                                q=self.location_prescribed_velocity)
107            self.__pde.setTolerance(self.rel_tol/10.)            self.__pde.getSolverOptions().setTolerance(self.rel_tol/10.)
108            self.velocity=self.__pde.getSolution()            self.velocity=self.__pde.getSolution()
109            update=penalty*div(self.velocity)            update=penalty*div(self.velocity)
110            self.pressure=self.pressure-update            self.pressure=self.pressure-update

Legend:
Removed from v.2344  
changed lines
  Added in v.2474

  ViewVC Help
Powered by ViewVC 1.1.26