darcy/examples/twophase_impes/user_fun.hpp

Go to the documentation of this file.

/* -*- mode: c++ -*-

  This file is part of the LifeV Applications.

  Author(s):  A. Fumagalli  <alessio.fumagalli@mail.polimi.it>
       Date: 2010-07-29

  Copyright (C) 2010 EPFL, Politecnico di Milano

  This program is free software; you can redistribute it and/or modify
  it under the terms of the GNU General Public License as published by
  the Free Software Foundation; either version 2.1 of the License, or
  (at your option) any later version.

  This program is distributed in the hope that it will be useful, but
  WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  General Public License for more details.

  You should have received a copy of the GNU General Public License
  along with this program; if not, write to the Free Software
  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
  USA
*/

/**
   @file user_fun.hpp
   @author A. Fumagalli <alessio.fumagalli@mail.polimi.it>
   @date 2010-07-29
*/

#ifndef __user_fun_H
#define __user_fun_H 1

#include <vector>


#include <boost/numeric/ublas/vector.hpp>
#include <boost/numeric/ublas/matrix.hpp>


#include <lifev/core/LifeV.hpp>

using namespace LifeV;

// ===================================================
//! User functions for the pressure equation
// ===================================================
namespace dataProblem
{
typedef boost::numeric::ublas::matrix<Real> Matrix;
typedef boost::numeric::ublas::vector<Real> Vector;

namespace dataPhysical
{
typedef boost::numeric::ublas::matrix<Real> Matrix;
typedef boost::numeric::ublas::vector<Real> Vector;

// Porosity
Real Phi ( const Real& x, const Real& y, const Real& z );

// Dynamic viscosity
const Real mu_w = 1; //1e-3; // [Pa * s]
const Real mu_n = 5; //2*1e-3; // [Pa * s]

// Density
const Real rho_w = 1; //960; // [Kg / m^3]
const Real rho_n = 0.5; //700; // [Kg / m^3]

// Relative permeability
Real k_rw ( const Real& S_w );
Real k_rn ( const Real& S_n );

// Absolute inverse permeability
const Matrix invK ( const Real& t, const Real& x, const Real& y, const Real& z ); // [m^2]

// Gravity acceleration, assumed to be g = [0, 0, -9.7803184]
const Real g = -9.7803184; // [m / s^2]

// Residual saturation
const Real S_wr = 0.1;
const Real S_nr = 0;

// Entry pressure
const Real pd = 1200.; // [Pa]

// Brooks-Corey constant
const Real lambda = 1;

// Capillary pressure
Real pc ( const Real& S_w ); // [Pa]

// First derivative of the relative permeability
Real Dk_rw ( const Real& S_w );
Real Dk_rn ( const Real& S_n );

// First derivative of the capillary pressure
Real Dpc ( const Real& S_w ); // [Pa]

}

// Inverse of permeability matrix
Matrix pressurePermeability ( const Real& /*t*/, const Real& x, const Real& y, const Real& z, const std::vector<Real>& u );

// Source term
Real pressureSource ( const Real& /*t*/, const Real& x, const Real& y, const Real& /*z*/, const ID& /*icomp*/);

// Boundary condition of Dirichlet
Real pressureDirichlet1 ( const Real& /* t */, const Real& x, const Real& y, const Real& z, const ID& /*icomp*/);

// Boundary condition of Dirichlet
Real pressureDirichlet2 ( const Real& /* t */, const Real& x, const Real& y, const Real& z, const ID& /*icomp*/);

// Boundary condition of Dirichlet
Real pressureDirichlet3 ( const Real& /* t */, const Real& x, const Real& y, const Real& z, const ID& /*icomp*/);

// Boundary condition of Neumann
Real pressureNeumann ( const Real& /* t */, const Real& x, const Real& y, const Real& z, const ID& icomp);

// Boundary condition of Robin
Real pressureRobin ( const Real& /* t */, const Real& x, const Real& y, const Real& z, const ID& /*icomp*/);

// ===================================================
//! User functions for the saturation equation
// ===================================================

// Inverse of permeability matrix
Matrix saturationPermeability ( const Real& t, const Real& x, const Real& y, const Real& z, const std::vector<Real>& u );

// Physical flux function
Vector saturationPhysicalFlux ( const Real& /*t*/, const Real& x, const Real& y, const Real& z, const std::vector<Real>& u );

// First derivative in u of the physical flux function
Vector saturationFirstDerivativePhysicalFlux ( const Real& /*t*/, const Real& x, const Real& y, const Real& z, const std::vector<Real>& u );

// Source term
Real saturationSource ( const Real& /*t*/, const Real& x, const Real& y, const Real& /*z*/, const ID& /*icomp*/);

// Initial condition
Real saturationInitialCondition ( const Real& /* t */, const Real& x, const Real& y, const Real& z, const ID& /*icomp*/ );

// Mass function
Real saturationMass ( const Real& /* t */, const Real& x, const Real& y, const Real& z, const ID& /*icomp*/ );

// Boundary condition of Dirichlet
Real saturationDirichlet1 ( const Real& /* t */, const Real& x, const Real& y, const Real& z, const ID& /*icomp*/);

// Boundary condition of Dirichlet
Real saturationDirichlet2 ( const Real& /* t */, const Real& x, const Real& y, const Real& z, const ID& /*icomp*/);

// Boundary condition of Dirichlet
Real saturationDirichlet3 ( const Real& /* t */, const Real& x, const Real& y, const Real& z, const ID& /*icomp*/);

// Boundary condition of Neumann
Real saturationNeumann ( const Real& /* t */, const Real& x, const Real& y, const Real& z, const ID& icomp);

// Boundary condition of Robin
Real saturationRobin ( const Real& /* t */, const Real& x, const Real& y, const Real& z, const ID& /*icomp*/);

}

#endif /* __user_fun_H */