sfepy.solvers.oseen module

class sfepy.solvers.oseen.Oseen(conf, context=None, **kwargs)[source]

The Oseen solver for Navier-Stokes equations.

Kind: ‘nls.oseen’

For common configuration parameters, see Solver.

Specific configuration parameters:


stabil_mat : str

The name of stabilization material.

adimensionalize : bool (default: False)

If True, adimensionalize the problem (not implemented!).

check_navier_stokes_residual : bool (default: False)

If True, check the Navier-Stokes residual after the nonlinear loop.

i_max : int (default: 1)

The maximum number of iterations.

eps_a : float (default: 1e-10)

The absolute tolerance for the residual, i.e. ||f(x^i)||.

eps_r : float (default: 1.0)

The relative tolerance for the residual, i.e. ||f(x^i)|| /

macheps : float (default: 2.2204460492503131e-16)

The float considered to be machine “zero”.

lin_red : float (default: 1.0)

The linear system solution error should be smaller than (eps_a * lin_red), otherwise a warning is printed.

lin_precision : float or None

If not None, the linear system solution tolerances are set in each nonlinear iteration relative to the current residual norm by the lin_precision factor. Ignored for direct linear solvers.

name = ‘nls.oseen’
class sfepy.solvers.oseen.StabilizationFunction(name_map, gamma=None, delta=None, tau=None, tau_red=1.0, tau_mul=1.0, delta_mul=1.0, gamma_mul=1.0, diameter_mode=’max’)[source]

Definition of stabilization material function for the Oseen solver.


  • tau_red <= 1.0; if tau is None: tau = tau_red * delta
  • diameter mode: ‘edge’: longest edge ‘volume’: volume-based, ‘max’: max. of previous

Get the maps of names and indices of variables in state vector.


Setup common problem-dependent data.

sfepy.solvers.oseen.are_close(a, b, rtol=0.2, atol=1e-08)[source]
sfepy.solvers.oseen.scale_matrix(mtx, indx, factor)[source]