IMEXSfloW  1.0
templategithubproject
 All Classes Files Functions Variables Pages
Public Member Functions | Public Attributes | List of all members
constitutive Module Reference

Constitutive equations. More...

Public Member Functions

subroutine init_problem_param
 Initialization of relaxation flags. More...
 
subroutine phys_var (Bj, r_qj, c_qj)
 Physical variables. More...
 
subroutine eval_local_speeds (qj, Bj, vel_min, vel_max)
 Local Characteristic speeds. More...
 
subroutine eval_local_speeds2 (qj, Bj, vel_min, vel_max)
 Local Characteristic speeds. More...
 
subroutine qc_to_qp (qc, B, qp)
 Conservative to physical variables. More...
 
subroutine qp_to_qc (qp, B, qc)
 Physical to conservative variables. More...
 
subroutine qc_to_qp2 (qc, B, qp)
 Conservative to 2nd set of physical variables. More...
 
subroutine qp2_to_qc (qp, B, qc)
 From 2nd set of physical to conservative variables. More...
 
subroutine eval_fluxes (Bj, c_qj, r_qj, c_flux, r_flux)
 Hyperbolic Fluxes. More...
 
subroutine eval_nonhyperbolic_terms (Bj, Bprimej, c_qj, c_nh_term_impl, r_qj, r_nh_term_impl)
 Non-Hyperbolic terms. More...
 
subroutine eval_explicit_forces (Bj, Bprimej, qj, expl_forces_term)
 Explicit Forces term. More...
 

Public Attributes

logical, dimension(:), allocatable implicit_flag
 
character(len=20) phase1_name
 
character(len=20) phase2_name
 
complex *16 h
 height More...
 
complex *16 u
 velocity (x direction) More...
 
real *8 grav
 gravitational acceleration More...
 

Detailed Description

Constitutive equations.

Definition at line 4 of file constitutive.f90.

Member Function/Subroutine Documentation

subroutine constitutive::eval_explicit_forces ( real*8, intent(in)  Bj,
real*8, intent(in)  Bprimej,
real*8, dimension(n_eqns), intent(in)  qj,
real*8, dimension(n_eqns), intent(out)  expl_forces_term 
)

Explicit Forces term.

This subroutine evaluates the forces to be treated explicitely in the DIRK numerical scheme (e.g. gravity)

Date
01/06/2012
Parameters
[in]qjconservative variables
[out]expl_forces_termforces term
[in]qjconservative variables
[out]expl_forces_termexplicit forces

Definition at line 453 of file constitutive.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

subroutine constitutive::eval_fluxes ( real*8, intent(in)  Bj,
complex*16, dimension(n_vars), intent(in), optional  c_qj,
real*8, dimension(n_vars), intent(in), optional  r_qj,
complex*16, dimension(n_eqns), intent(out), optional  c_flux,
real*8, dimension(n_eqns), intent(out), optional  r_flux 
)

Hyperbolic Fluxes.

This subroutine evaluates the numerical fluxes given the conservative variables qj, accordingly to the equations for the single temperature model introduced in Romenki et al. 2010.

Date
22/11/2016
Parameters
[in]c_qjcomplex conservative variables
[in]r_qjreal conservative variables
[out]c_fluxcomplex analytical fluxes
[out]r_fluxreal analytical fluxes

Definition at line 304 of file constitutive.f90.

Here is the caller graph for this function:

subroutine constitutive::eval_local_speeds ( real*8, dimension(n_vars), intent(in)  qj,
real*8, intent(in)  Bj,
real*8, dimension(n_vars), intent(out)  vel_min,
real*8, dimension(n_vars), intent(out)  vel_max 
)

Local Characteristic speeds.

This subroutine evaluates an the largest pos and neg characteristic speeds from the conservative variables qj, without any change on u and h, changing u when h is small (Eq. 2.17 K&P 2007).

Date
10/04/2012

Definition at line 105 of file constitutive.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

subroutine constitutive::eval_local_speeds2 ( real*8, dimension(n_vars), intent(in)  qj,
real*8, intent(in)  Bj,
real*8, dimension(n_vars), intent(out)  vel_min,
real*8, dimension(n_vars), intent(out)  vel_max 
)

Local Characteristic speeds.

This subroutine evaluates an the largest pos and neg characteristic speeds from the conservative variables qj, without any change on u and h.

Date
10/04/2012

Definition at line 129 of file constitutive.f90.

Here is the caller graph for this function:

subroutine constitutive::eval_nonhyperbolic_terms ( real*8, intent(in)  Bj,
real*8, intent(in)  Bprimej,
complex*16, dimension(n_vars), intent(in), optional  c_qj,
complex*16, dimension(n_eqns), intent(out), optional  c_nh_term_impl,
real*8, dimension(n_vars), intent(in), optional  r_qj,
real*8, dimension(n_eqns), intent(out), optional  r_nh_term_impl 
)

Non-Hyperbolic terms.

This subroutine evaluates the non-hyperbolic terms (relaxation terms and forces) of the system of equations, both for real or complex inputs. These terms are treated implicitely in the DIRK numerical scheme.

Date
01/06/2012
Parameters
[in]c_qjcomplex conservative variables
[in]r_qjreal conservative variables
[out]c_nh_term_implcomplex non-hyperbolic terms
[out]r_nh_term_implreal non-hyperbolic terms

Definition at line 378 of file constitutive.f90.

Here is the caller graph for this function:

subroutine constitutive::init_problem_param ( )

Initialization of relaxation flags.

This subroutine set the number and the flags of the non-hyperbolic terms.

Date
07/09/2012

Definition at line 34 of file constitutive.f90.

Here is the caller graph for this function:

subroutine constitutive::phys_var ( real*8, intent(in)  Bj,
real*8, dimension(n_vars), intent(in), optional  r_qj,
complex*16, dimension(n_vars), intent(in), optional  c_qj 
)

Physical variables.

This subroutine evaluates from the conservative local variables qj the local physical variables ( $\alpha_i, u_i, \rho_i, x_i, $) for the two phases and the mixture denisity $ \rho $.

Parameters
[in]r_qjreal conservative variables
[in]c_qjcomplex conservative variables
Date
15/08/2011

Definition at line 58 of file constitutive.f90.

Here is the caller graph for this function:

subroutine constitutive::qc_to_qp ( real*8, dimension(n_vars), intent(in)  qc,
real*8, intent(in)  B,
real*8, dimension(n_vars), intent(out)  qp 
)

Conservative to physical variables.

This subroutine evaluates from the conservative variables qc the array of physical variables qp:

  • qp(1) = $ h+B $
  • qp(2) = $ u $

The physical variables are those used for the linear reconstruction at the cell interfaces (limiters are applied to the reconstructed slopes) and for the boundary conditions. Please note that the physical variable u is modified according to Eq. 2.17 of K&P 2007 by the call to the subroutine phys_var.

Parameters
[in]qcconservative variables
[out]qpphysical variables
Date
15/08/2011

Definition at line 177 of file constitutive.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

subroutine constitutive::qc_to_qp2 ( real*8, dimension(n_vars), intent(in)  qc,
real*8, intent(in)  B,
real*8, dimension(n_vars), intent(out)  qp 
)

Conservative to 2nd set of physical variables.

This subroutine evaluates from the conservative variables qc the array of physical variables qp:

  • qp(1) = $ h $
  • qp(2) = $ u $

The physical variables are those used for the linear reconstruction at the cell interfaces. Limiters are applied to the reconstructed slopes.

Parameters
[in]qcconservative variables
[out]qpphysical variables
Date
15/08/2011

Definition at line 241 of file constitutive.f90.

Here is the call graph for this function:

Here is the caller graph for this function:

subroutine constitutive::qp2_to_qc ( real*8, dimension(n_vars), intent(in)  qp,
real*8, intent(in)  B,
real*8, dimension(n_vars), intent(out)  qc 
)

From 2nd set of physical to conservative variables.

This subroutine evaluates the conservative variables qc from the array of physical variables qp:

  • qp(1) = $ h $
  • qp(2) = $ u $
Parameters
[in]qpphysical variables
[out]qcconservative variables
Date
15/08/2011

Definition at line 271 of file constitutive.f90.

Here is the caller graph for this function:

subroutine constitutive::qp_to_qc ( real*8, dimension(n_vars), intent(in)  qp,
real*8, intent(in)  B,
real*8, dimension(n_vars), intent(out)  qc 
)

Physical to conservative variables.

This subroutine evaluates the conservative variables qc from the array of physical variables qp:

  • qp(1) = $ h + B $
  • qp(2) = $ hu $
Parameters
[in]qpphysical variables
[out]qcconservative variables
Date
15/08/2011

Definition at line 206 of file constitutive.f90.

Here is the caller graph for this function:

Member Data Documentation

real*8 constitutive::grav

gravitational acceleration

Definition at line 22 of file constitutive.f90.

complex*16 constitutive::h

height

Definition at line 18 of file constitutive.f90.

logical, dimension(:), allocatable constitutive::implicit_flag

Definition at line 10 of file constitutive.f90.

character(len=20) constitutive::phase1_name

Definition at line 12 of file constitutive.f90.

character(len=20) constitutive::phase2_name

Definition at line 13 of file constitutive.f90.

complex*16 constitutive::u

velocity (x direction)

Definition at line 19 of file constitutive.f90.

The documentation for this module was generated from the following file:
Generated on Tue Nov 22 2016 16:23:06 for IMEXSfloW by   doxygen 1.8.7