# Models

The following models are currently implemented in OpenCMP. In general, all models are expected to include a standard Galerkin finite element method formulation, a discontinuous Galerkin formulation, and a diffuse interface formulation. Nonlinear models are further expected to include formulations for Oseen-style linearization and IMEX linearization. The following sections provide the weak forms for all implemented formulations of the available models.

## Poisson Equation

The Poisson equation takes the following form on a simulation domain . There are three possible types of boundary condition - Dirichlet, Neumann, and Robin - applied on boundary regions , , and respectively.

where is a constant typically taken to be the diffusivity and is some source function.

The different formulations of the Poisson equation finite element weak form are given below. Note that in all cases is the trial function. In the case of the discontinuous Galerkin method, and refer to the average and jump operators respectively and is the penalty coefficient. In the case of the diffuse interface method, is the phase field, are masks for different boundary condition regions, and is the Nitsche method penalty parameter. Also note that for the diffuse interface method all integrals become volume integrals over the enclosing simple domain .

### Standard Galerkin Finite Element Formulation

The standard Galerkin finite element method formulation is as follows:

Dirichlet boundary conditions are imposed strongly by setting values at applicable boundary degrees of freedom.

### Discontinuous Galerkin Formulation

The discontinuous Galerkin formulation is as follows:

where volume integrals are over individual mesh elements then summed over the entire triangulation . Surface integrals are over mesh element facets and summed over either interior facets or boundary facets , , and .

### Diffuse Interface Formulation

The diffuse interface formulation is as follows for the standard Galerkin finite element method:

and for the discontinuous Galerkin method:

## Stokes Equations

The Stokes equations take the following form on a simulation domain . There are two possible types of boundary condition - velocity Dirichlet and normal stress - applied on boundary regions and respectively.

where is the velocity and is the pressure. Furthermore, is the constant kinematic viscosity and is some body force.

The different formulations of the Stokes equations finite element weak form are given below. Note that in all cases and are the trial functions for velocity and pressure respectively. In the case of the discontinuous Galerkin method, and refer to the average and jump operators respectively and is the penalty coefficient. In the case of the diffuse interface method, is the phase field, are masks for different boundary condition regions, and is the Nitsche method penalty parameter. Also note that for the diffuse interface method all integrals become volume integrals over the enclosing simple domain .

### Standard Galerkin Finite Element Formulation

The standard Galerkin finite element method formulation is as follows:

Dirichlet boundary conditions are imposed strongly by setting values at applicable boundary degrees of freedom.

### Discontinuous Galerkin Formulation

The discontinuous Galerkin formulation is as follows:

where volume integrals are over individual mesh elements then summed over the entire triangulation . Surface integrals are over mesh element facets and summed over either interior facets or boundary facets and .

### Diffuse Interface Formulation

The diffuse interface formulation is as follows for the standard Galerkin finite element method:

and for the discontinuous Galerkin method:

## Multi-Component Flow

Note

To be written once discontinuous Galerkin and diffuse interface formulations are available for the multi-component flow model.