Boundary conditions (b.c.) are constraints necessary for the solution of a boundary value problem. A boundary value problem is a differential equation (or system of differential equations) to be solved in a domain on whose boundary a set of conditions is known. It is opposed to the “initial value problem”, in which only the conditions on one extreme of the interval are known. Boundary value problems are extremely important as they model a vast amount of phenomena and applications, from solid mechanics to heat transfer, from fluid mechanics to acoustic diffusion. They rise naturally in every problem based on a differential equation to be solved in space, while initial value problems usually refer to problems to be solved in time. Boundary value problems have been extensively studied by Jacques Charles François Sturm (1803-1855) and Joseph Liouville (1809-1882), who studied the eigenvalues of a linear differential equation of the second order1

${}^{1}$

. They studied the conditions which guarantee the existence and uniqueness of the solution of the differential problem and how it is affected by the boundary conditions2

${}^{2}$

. The Sturm-Liouville theory is extremely important for any computational problem because it enables to understand if a problem is “well-posed” and how it is possible to obtain the solution.

Types of boundary conditions

Both ordinary and partial differential equations need boundary conditions to be solved. Different types of boundary conditions can be imposed on the boundary of the domain (figure 1). The choice of the boundary condition is fundamental for the resolution of the computational problem: a bad imposition of b.c. may lead to the divergence of the solution or to the convergence to a wrong solution.

Dirichlet boundary condition

The Dirichlet boundary condition is a type of boundary condition named after Peter Gustav Lejeune Dirichlet (1805–1859, figure 2)3

${}^{3}$

.

This condition specifies the value that the unknown function needs to take on along the boundary of the domain. Given, for example, the Laplace equation, the boundary value problem with the Dirichlet b.c. is written as:

is the independent variable (e.g. the spatial coordinates), Ω

$\mathrm{\Omega}$

is the function domain, ∂Ω

$\mathrm{\partial}\mathrm{\Omega}$

is the boundary of the domain, and f

$f$

is a given scalar function defined on ∂Ω

$\mathrm{\partial}\mathrm{\Omega}$

. In the framework of numerical simulations, it is usually imposed directly in the algebraic system to be solved. Let’s consider the following algebraic system derived from a numerical algorithm:

The Neumann boundary condition is a type of boundary condition, named after Carl Neumann (1832 – 1925, figure 3)3

${}^{3}$

. When imposed on an ordinary or a partial differential equation, it specifies the values that the derivative of a solution is going to take on the boundary of the domain. Given, for example, the Laplace equation, the boundary value problem with the Dirichlet b.c. is written as:

is the unit normal to the boundary surface, if Ω⊂R3

$\mathrm{\Omega}\subset {R}^{3}$

.

In the case of ordinary differential equations (i.e. Ω⊂R1

$\mathrm{\Omega}\subset {R}^{1}$

), the derivative normal to the boundary coincides with the global derivative φ′

${\phi}^{\prime}$

. In the rare cases in which a temporal dependency is solved through a finite element approach (instead of the more usual finite difference), this type of boundary condition is the most common. Neumann boundary condition is also called “natural” because it naturally appears in the development of the weak formulation in any finite element approach. Let’s consider the following simple equation:

is a given scalar function. This equation rules many phenomena, for instance, the thermal diffusion in 1D and the tension/compression of a beam. The finite element method consists in rewriting the equation from a differential (strong) form to an integral (weak) formulation. This transformation is done through two steps:

Thus, a term including the derivative of the unknown field on the boundary naturally appears; for 1D problems, this term refers to the extremes of the interval, for 2D problems it refers to the contour of the domain, and for 3D problems it refers to the boundary surfaces. The presence of the boundary term on the right-hand side highlights two properties of Neumann boundary conditions:

homogeneous Neumann b.c. are naturally satisfied without any explicit imposition

since Dirichlet b.c. are usually applied by modifying the right-hand term, a homogeneous Neumann condition is applied on all the boundaries where a Dirichlet b.c. is imposed.

Robin boundary condition

The Robin boundary condition is a type of boundary condition named after Victor Gustave Robin (1855–1897)4

${}^{4}$

. It consists of a linear combination of the values of the field and its derivatives on the boundary. Given, for example, the Laplace equation, the boundary value problem with the Dirichlet b.c. is written as:

are real parameters. This condition is also called “impedance condition”.

Mixed boundary condition

It consists of applying different types of boundary conditions in different parts of the domain. It is important to notice that boundary conditions must be applied on the whole boundary: the “free” boundary is anyways subjected to a homogeneous Neumann condition. The mixed boundary condition differs from the Robin condition because the latter consists of different types of boundary conditions applied to the same region of the boundary, while the mixed condition implies different types of b.c. applied to different parts of the boundary.

Cauchy boundary condition

The Cauchy boundary condition is a condition on both the unknown field and its derivatives5

${}^{5}$

. It differs from the Robin condition because the Cauchy condition implies the imposition of two constraints (1 Dirichlet b.c. + 1 Neumann b.c.), while the Robin condition implies only one constraint on the linear combination of the unknown function and its derivatives.

Applications

Structural and Solid mechanics

Dirichlet boundary conditions.

Solid mechanics is usually modelled through a displacement-based model, thus Dirichlet boundary conditions usually consist in imposing the displacement of the structure in given points. Structural mechanics is often based on formulations which include relative rotations, whose numerical resolution requires nonlinear shape functions in the finite element approximation. For instance, frame structures are based on the beam theory and the related finite element has 6 degrees of freedom (3 displacements + 3 rotations in a 3D space). For a 2D problem, each node of the boundary has 3 degrees of freedom on which Dirichlet boundary conditions can be applied: 2 displacements (ux

${u}_{x}$

and uy

${u}_{y}$

) and 1 rotation (ω

$\omega $

). These constraints are usually depicted as follows:

In the table, external constraints to the ground (i.e. set the value to zero) are reported, but the symbols are used also to express a fixed displacement/rotation different than zero.

Neumann boundary conditions

In solid mechanics, the spatial derivatives of displacements are related to the strain tensor. In elasticity, the strain is proportional to the stress, hence the Neumann boundary condition refers to both imposed strains and stresses. Since the stress is also linked to the external forces through the Cauchy’s stress principle, the Neumann condition is also used to apply external loads. As stated in the section dedicated to Neumann boundary conditions, the homogeneous condition is naturally satisfied, so “free” boundaries may not be modelled explicitly.

Robin boundary condition

It is used to model the mechanical impedance of a structure, thus how much it resists to motion when subjected to a harmonic load.

Fluid mechanics

Dirichlet boundary conditions

In computational fluid mechanics, the classical Dirichlet boundary condition consists of the value of velocity and/or pressure to be taken by a certain set of nodes. It is common to refer to some sets of b.c. according to the following terminology:

slip boundary condition: the velocity normal to the boundary is set to zero, while the velocity parallel to the boundary is let free

no-slip boundary condition: both the velocity normal to the boundary and the velocity parallel to the boundary are set equal to zero.

At least one homogeneous b.c. on the pressure (i.e. p=0

$p=0$

) has to be imposed as a reference for open domains, for instance, in the highest boundary of the air domain.

Neumann boundary conditions

Constraints on the derivative of velocity or pressure fields are mainly used in two cases. The first case is the application of a symmetry plane, thus:

Since this condition is always applied in addition to a Dirichlet b.c. on the velocity normal to the boundary, it is naturally satisfied. The second application is the modelling of wall friction, in the case when it is proportional to the strain rate:

It is used to describe semi-reflective walls, which partially absorb waves. It is not a very common application and it can be used only for pressure-based models. It is mostly used for acoustic applications.

Thermodynamics

Dirichlet boundary conditions

In thermodynamics, Dirichlet boundary conditions consist of surfaces (in 3D problems) held at a fixed temperatures.

Neumann boundary conditions

In thermodynamics, the Neumann boundary condition represents the heat flux across the boundaries. The perfect insulator reflects a homogeneous condition (naturally satisfied), while all warmed and cooled boundaries require to explicitly assign the boundary condition. This is normally the case of electronic components (inward heat flux) of external cooling spray/channel (outward heat flux).

Electromagnetism

Maxwell equations are commonly solved through a potential formulation. In this section, the (A,φ

$A,\phi $

) formulation is considered, where A

$A$

is the magnetic vector potential and φ

$\phi $

is the scalar electric potential. Under certain conditions, the (A,φ

$A,\phi $

) can be decoupled and the electric potential can be computed through the Laplace equation:

is usually imposed on the boundary sections of the conductive domain; in the case of wires, the values on one section are normally set to zero, while a fixed value of the electric potential is fixed on the second section. The conditions A

$A$

usually constrain the magnetic field to be tangent to the external boundary, i.e. to have all the magnetic lines inside the computational domain.

Neumann boundary conditions

In electromagnetic modelling, under certain hypotheses, ∇φ

$\mathrm{\nabla}\phi $

is the electric current density. The imposition of a homogeneous Neumann boundary condition (i.e.:math:nablavarphicdot n=0) means forcing the electric current not to cross the boundaries. This condition is also referred as “insulating boundary” and represents the behavior of a perfect insulator.

Robin boundary condition

It is used to model the impedance of an electric circuit, thus the opposition that a circuit presents to a current when a voltage is applied. It is also used to model the impedance of the electromagnetic wave.

: Cheng, A. and D. T. Cheng (2005). Heritage and early history of the boundary element method, Engineering Analysis with Boundary Elements, 29, 268–302.

4

${}^{4}$

: Gustafson, K., (1998). Domain Decomposition, Operator Trigonometry, Robin Condition, Contemporary Mathematics, 218. 432–437.

5

${}^{5}$

: Morse, P. M. and Feshbach, H. Methods of Theoretical Physics, Part I. New York: McGraw-Hill, pp. 678-679, 1953.

Contents

Cookie Policy

SimScale uses cookies to improve your user experience. By using this website you consent to our cookie policy.

This website uses cookies so that we can provide you with the best user experience possible. Cookie information is stored in your browser and performs functions such as recognising you when you return to our website and helping our team to understand which sections of the website you find most interesting and useful.

Strictly Necessary Cookies

Strictly Necessary Cookie should be enabled at all times so that we can save your preferences for cookie settings.

If you disable this cookie, we will not be able to save your preferences. This means that every time you visit this website you will need to enable or disable cookies again.