Physical interpretation of Robin boundary conditions

Question

In a (bounded) domain $\Omega \subset \mathbb{R}^n$, if we're studying the Laplace equation or heat equation or such PDE's we can impose the Dirichlet $u|_{\partial\Omega} \equiv 0$, Neumann $D_{\nu} u|_{\partial\Omega}\equiv 0$ or Robin (for $\alpha \in \mathbb{R}$) $(D_{\nu} u + \alpha u)|_{\partial \Omega} \equiv 0$.

I know that, for example for the heat equation, Dirichlet eigenvalues correspond physically to the boundary being in contact with a (large) heat bath at $T=0$. Or, in the Laplace equation, if we're interested in the modes supported by $\Omega$ (as a drum), Dirichlet boundary conditions can be thought of keeping the boundary from moving.

Neumann boundary conditions, for the heat flow, correspond to a perfectly insulated boundary. For the Laplace equation and drum modes, I think this corresponds to allowing the boundary to flap up and down, but not move otherwise.

---

My question is: what sort of physical interpretations are there for the Robin boundary conditions? Wikipedia says that they are related to electromagnetic problems, but gives no details. I'd be happy with answers that are not necessarily physics-related, for example, if there was somewhere that Robin boundary conditions naturally arise in a mathematical context, I'd be interested to know about that as well.

Answer 1

Here is an example where $\Omega = \mathbb{R}^3$. One way to establish dispersion for the wave equation involves taking a temporal Fourier transform. In order to do this one has to multiply by a cutoff function supported in $t \in [0,\infty)$. You then get the equation

$(\Delta+\omega^2)\psi = F$

where $\psi$ is the temporal Fourier transform of the product of the original solution with the cutoff, $\omega$ is the Fourier variable, and $F$ is a function controllable by initial data via a finite time energy inequality. If this plan of attack is going to work, we need to make sure that $\psi$ is uniquely determined by $F$. This of course requires appropriate boundary conditions at $\infty$. These turn out to be

1) $\psi = O\left(|x|^{-1}\right)$

2) $\frac{\partial\psi}{\partial r} - i\omega\psi = O\left(|x|^{-2}\right)$

This is a sort of Robin condition at infinity. See http://terrytao.wordpress.com/2011/04/21/the-limiting-absorption-principle/ for more details.

Answer 2

In the case of the heat equation there is a very typical example.

If you have heat transfer by convection in one of the boundaries of your domain.

Imagine you have a solid where you are solving the heat equation, and in some of the boundaries, you have a liquid in contact at temperature $T_l$. If this liquid is moving (forced convection) or if you let it move by buoyancy of the hotter parts vs the colder ones, then the boundary condition you must impose on the solid on that boundary is:

$$\vec q\cdot \vec{n} = h\left(T-T_l\right)$$ where $h$ is the convection coefficient, and $\vec q$ is the heat loss by that boundary:

$$\vec q = -\kappa \boldsymbol{\nabla}T$$

so finally you have

$$\kappa\frac{\partial T}{\partial\vec{n}}+h\left(T-T_l\right)=0$$

and if you do a variable change like $\theta = T - T_l$

you get your homogeneous robin condition

$$\kappa\frac{\partial \theta}{\partial\vec{n}}+h\theta=0$$

Answer 3

Check out Section II.1.7 of Tikhonov & Samarskii's text Equations of Mathematical Physics for a nice discussion of the physical interpretation of Dirichlet, Neumann and Robin boundary conditions for the simple example of the 1+1 wave equation ($u_{tt} = u_{xx}$) describing the transverse vibrations of a spring on the interval $x\in[0,l]$.

Here's a brief summary. The value of $u(t,0)$ is the transverse position of the spring at $x=0$. The value of $u_x(t,0)$ is the vertical component of the tension, so anything connected to the spring at this end will experience this vertical force, and by Newton's third law apply the same force to the end of the spring. Here are the interpretations. Dirichlet, $u(t,0)=0$: the end of the spring is transversally clamped or fixed. Neumann, $u_x(t,0)=0$: the end of the spring is undergoes free transverse motion, conversely no external transverse force acts on this end. Robin, $u_x(t,0) = k u(t,0)$: a linearly restorative transverse force is applied to the end of the spring, that is the end is transversally restrained, but elastically rather than rigidly. Actually, for the force to be restorative, one must pick a particular sign of $k$, which I don't feel like figuring out at the moment.

Answer 4

There is a book by Daniel J. Hoppe, "Impedance Boundary Conditions In Electromagnetics", CRC, 1995. I give here its summary from Amazon.com which seems to answer the question:

Electromagnetic scattering from complex objects has been an area of in-depth research for many years. A variety of solution methodologies have been developed and utilised for the solution of ever increasingly complex problems. Among these methodologies, the subject of impedance boundary conditions has interested the authors for some time. In short, impedance boundary conditions allow one to replace a complex structure with an appropriate impedance relationship between the electric and magnetic fields on the surface of the object. This simplifies the solution of the problem considerably, allowing one to ignore the complexity of the internal structure beneath the surface. This book examines impedance boundary conditions in electromagnetics. The introductory chapter provides a presentation of the role of the impedance boundary conditions in solving practical electromagnetic problems and some historical background. One of the main objectives of this book is to present a unified and thorough discussion of this important subject. A method based on a spectral domain approach is presented to derive the Higher Order Impedance Boundary Conditions (HOIBC). The method includes all of the existing approximate boundary conditions, such as the Standard Impedance Boundary Condition, the Tensor Impedance Boundary Condition and the Generalised Impedance Boundary Conditions, as special cases. The special domain approach is applicable to complex coatings and surface treatments as well as simple dielectric coatings. The spectral domain approach is employed to determine the appropriate boundary conditions for planar dielectric coatings, chiral coatings and corregated conductors. The accuracy of the proposal boundary conditions is discussed. The approach is then extended to include the effects of curvature and is applied to curved dielectric and chiral coatings. Numerical data is presented to critically assess the accuracy of the results obtained using various forms of the impedance boundary conditions. A number of appendices that provide more detail on some of the topics addressed in the main body of the book and a selective list of references directly related to the topics addressed in this book are also included.

Answer 5

Robbins conditions correspond to immersing one end of the rod in a constant temperature bath. The temperature bath does not have to match the initial temperature at the end of the rod. The immersed end of rod may come into thermal equilibrium with the bath, or a constant temperature difference may come about (to ensure a constant heat flux).

In thermodynamics, reference is always made to work reservoirs (ideal mechanical energy sources that can do work reversibly without any entropy increase) and reference is made to heat reservoirs that can transfer heat without any change in temperature (of the reservoir). Such a heat reservoir can be attained by making it a water/ice mixture - so that the temperature of the bath is maintained due to the phase change occurring.

The Robbins condition itself can be interpreted as a statement of Newton's law of heating (cooling): the rate of heat transfer is directly proportional to the difference in temperature (between the hot and cold object).

dQ/dt ̇ = h(To - T) Newton’s law of heat transfer

Here dQ/dt represents the rate of heat transfer between the end of rod and the temperature bath.

The constant h is a thermal boundary layer coefficient whose reciprocal represents IMPEDANCE to the transfer of heat.

To is the temperature of the bath.

T is the temperature of the end of the rod.

We can get a Robbins condition of this by invoking Fourier's first law of heat conduction:

q = -k dT/dx

Here q is the heat flux and equals (1/A)(dQ/dt) where dQ/dt is the rate of heat transfer and A is the cross sectional area of the solid rod.

Combining Newton's law of heat transfer with Fourier's first law of conduction yields

-kA dT/dx = h(To - T) right end immersed

kA dT/dx = h(To - T) left end immersed

An equivalent form is dT/dx - aT = b or dT/dx + aT = b where a and b are constant dependent on k, A, h and To.

The concept of impedance is hidden in a Robbins condition because Newton's law of heating (cooling) is hidden. The Robbins condition comes about from applying Newton's law.

Newton's law of cooling is an impedance equation and the reciprocal of h (the thermal boundary layer coefficient) acts like a resistance. A good analogy is to compare the rate of heat transfer to current, the temperature difference to voltage, and the thermal boundary layer to a resistor.