/[escript]/trunk/escript/py_src/flows.py
ViewVC logotype

Diff of /trunk/escript/py_src/flows.py

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

revision 3074 by gross, Tue Jul 27 01:47:45 2010 UTC revision 3360 by jfenwick, Thu Nov 18 00:20:21 2010 UTC
# Line 156  class DarcyFlow(object): Line 156  class DarcyFlow(object):
156     def setValue(self,f=None, g=None, location_of_fixed_pressure=None, location_of_fixed_flux=None, permeability=None):     def setValue(self,f=None, g=None, location_of_fixed_pressure=None, location_of_fixed_flux=None, permeability=None):
157        """        """
158        assigns values to model parameters        assigns values to model parameters
159          
160        :param f: volumetic sources/sinks        :param f: volumetic sources/sinks
161        :type f: scalar value on the domain (e.g. `Data`)        :type f: scalar value on the domain (e.g. `Data`)
162        :param g: flux sources/sinks        :param g: flux sources/sinks
# Line 165  class DarcyFlow(object): Line 165  class DarcyFlow(object):
165        :type location_of_fixed_pressure: scalar value on the domain (e.g. `Data`)        :type location_of_fixed_pressure: scalar value on the domain (e.g. `Data`)
166        :param location_of_fixed_flux:  mask for locations where flux is fixed.        :param location_of_fixed_flux:  mask for locations where flux is fixed.
167        :type location_of_fixed_flux: vector values on the domain (e.g. `Data`)        :type location_of_fixed_flux: vector values on the domain (e.g. `Data`)
168        :param permeability: permeability tensor. If scalar ``s`` is given the tensor with        :param permeability: permeability tensor. If scalar ``s`` is given the tensor with ``s`` on the main diagonal is used.
       ``s`` on the main diagonal is used.  
169        :type permeability: scalar or tensor values on the domain (e.g. `Data`)        :type permeability: scalar or tensor values on the domain (e.g. `Data`)
170    
171        :note: the values of parameters which are not set by calling ``setValue`` are not altered.        :note: the values of parameters which are not set by calling ``setValue`` are not altered.
172        :note: at any point on the boundary of the domain the pressure (``location_of_fixed_pressure`` >0)        :note: at any point on the boundary of the domain the pressure (``location_of_fixed_pressure`` >0)
173        or the normal component of the flux (``location_of_fixed_flux[i]>0`` if direction of the normal        or the normal component of the flux (``location_of_fixed_flux[i]>0``) if direction of the normal is along the *x_i* axis.
174        is along the *x_i* axis.  
175        """        """
176        if f !=None:        if f !=None:
177       f=util.interpolate(f, self.__pde_p.getFunctionSpaceForCoefficient("X"))       f=util.interpolate(f, self.__pde_p.getFunctionSpaceForCoefficient("X"))
# Line 306  class DarcyFlow(object): Line 306  class DarcyFlow(object):
306    
307        :note: The problem is solved as a least squares form        :note: The problem is solved as a least squares form
308        *(K^[-1]+D^* (DKD^*)^[-1] D)u+G p=D^* (DKD^*)^[-1] f + K^[-1]g*        *(K^[-1]+D^* (DKD^*)^[-1] D)u+G p=D^* (DKD^*)^[-1] f + K^[-1]g*
309        *G^*u+G^* K Gp=G^*g*        *G^*u+*G^* K Gp=G^*g*
   
310        where *D* is the *div* operator and *(Gp)_i=p_{,i}* for the permeability *K=k_{ij}*.        where *D* is the *div* operator and *(Gp)_i=p_{,i}* for the permeability *K=k_{ij}*.
311        """        """
312        self.verbose=verbose        self.verbose=verbose

Legend:
Removed from v.3074  
changed lines
  Added in v.3360

  ViewVC Help
Powered by ViewVC 1.1.26