Robin boundary condition explained

In mathematics, the Robin boundary condition (; properly in French ʁɔbɛ̃/), or third type boundary condition, is a type of boundary condition, named after Victor Gustave Robin (1855 - 1897).[1] When imposed on an ordinary or a partial differential equation, it is a specification of a linear combination of the values of a function and the values of its derivative on the boundary of the domain. Other equivalent names in use are Fourier-type condition and radiation condition.[2]

Definition

Robin boundary conditions are a weighted combination of Dirichlet boundary conditions and Neumann boundary conditions. This contrasts to mixed boundary conditions, which are boundary conditions of different types specified on different subsets of the boundary. Robin boundary conditions are also called impedance boundary conditions, from their application in electromagnetic problems, or convective boundary conditions, from their application in heat transfer problems (Hahn, 2012).

If Ω is the domain on which the given equation is to be solved and ∂Ω denotes its boundary, the Robin boundary condition is:[3]

au+b

\partialu
\partialn

=g    on\partial\Omega

for some non-zero constants a and b and a given function g defined on ∂Ω. Here, u is the unknown solution defined on Ω and denotes the normal derivative at the boundary. More generally, a and b are allowed to be (given) functions, rather than constants.

In one dimension, if, for example, Ω = [0,1], the Robin boundary condition becomes the conditions:

\begin{align}au(0)-bu'(0)&=g(0)\au(1)+bu'(1)&=g(1)\end{align}

Notice the change of sign in front of the term involving a derivative: that is because the normal to [0,1] at 0 points in the negative direction, while at 1 it points in the positive direction.

Application

Robin boundary conditions are commonly used in solving Sturm - Liouville problems which appear in many contexts in science and engineering.

In addition, the Robin boundary condition is a general form of the insulating boundary condition for convection–diffusion equations. Here, the convective and diffusive fluxes at the boundary sum to zero:

ux(0)c(0)-D

\partialc(0)
\partialx

=0

where D is the diffusive constant, u is the convective velocity at the boundary and c is the concentration. The second term is a result of Fick's law of diffusion.

Bibliography

Notes and References

  1. Gustafson, K., (1998). Domain Decomposition, Operator Trigonometry, Robin Condition, Contemporary Mathematics, 218. 432 - 437.
  2. Logan, J. David, (2001). Transport Modeling in Hydrogeochemical Systems. Springer.
  3. Book: Finite Element Analysis with Error Estimators: An Introduction to the FEM and Adaptive Error Analysis for Engineering Students. J. E. Akin. Butterworth-Heinemann. 2005. 9780080472751. 69.