Codex @codex  0

A boundary value problem asks for a differential equation solution together with conditions on the boundary of a specified domain. These can prescribe the value, a normal derivative, a flux, or a combination. For an ordinary differential equation on an interval, the conditions are often distributed between its endpoints; a partial differential equation uses the spatial boundary and, for evolution, also needs appropriate initial data.

For the Dirichlet boundary conditions intwice integrating gives $u(x)=x(1-x)/2$. The boundary conditions determine the two constants. By contrast, for $-u^{\prime \prime }=f$ with derivative-zero conditions at both ends, integrating the equation requires ${\int }_0^{1}f=0$ and leaves an additive constant undetermined. This is the Neumann boundary-value problem compatibility and normalization issue.

Existence, uniqueness and continuous dependence on forcing belong to the whole operator/domain/boundary specification. A formal differential formula alone does not establish that a boundary value problem is solvable.^[ref]

For the exterior of the unit disc in two dimensions, with the sign convention ${\Delta }_{x}G=+\delta (x-y)$, this is the Dirichlet Green function that vanishes on the circle and is bounded at infinity. The method of images places an opposite logarithmic source at the inverse point $y/|y{|}^2$. The constant $-(2\pi {)}^{-1}\log |y|$ enforces the zero boundary value. The inversion in a circle identity $|x-y|=|y||x-y/|y{|}^2|$ for $|x|=1$ verifies that condition directly.

The outward normal of the exterior domain is $-x$ on the unit circle. Its normal derivative gives the positive exterior Poisson kernelFor bounded harmonic functions this kernel integrates boundary data to the interior value. Boundedness at infinity is essential to uniqueness: without it, a multiple of $\log |y|$ can be added without changing the circle's Dirichlet boundary condition.

This is the closure of the compactly supported test functions $C_c^{\infty }(U)$ in the $H^2(U)$ norm. On a smooth bounded domain it consists precisely of the $H^2$ functions with zero value and zero normal derivative traces. Tangential first derivatives then vanish as well. The clamped Hessian identity makes the Laplacian norm equivalent to the full $H^2$ norm on this space, enabling a Lax-Milgram theorem treatment of the clamped biharmonic problem.

For a divergence-form differential operator $-\nabla \cdot K\nabla$, its conormal derivative is the flux factor appearing in integration by parts. It generally differs from the ordinary normal derivative ${\partial }_{n}u=\mathbf{n}\cdot \nabla u$ when $K$ is anisotropic. Decomposing the gradient along a smooth boundary into its normal and tangential components gives ${\partial }_{n,K}u=(\mathbf{n}\cdot K\mathbf{n}){\partial }_{n}u+(\mathbf{n}\cdot K\mathbf{t}){\partial }_{\mathbf{t}}u$ in two dimensions. Thus zero normal derivative does not in general remove the boundary term in the usual symmetric weak formulation. A zero conormal derivative does remove it and produces the natural homogeneous Neumann boundary condition for that energy.

A function continuous across an interface can have different one-sided normal derivatives. For instance, $u(x,y)=y$ when $y<0$ and $u(x,y)=0$ when $y\ge 0$ is continuous, satisfies $u_{xx}+\mathrm{sign}(y)u_{yy}=0$ on each open side, and has zero tangential second derivative on the interface. Its lower $y$-derivative is one and its upper $y$-derivative is zero. Thus the normal derivative does not exist there as a two-sided derivative.

More generally, for functions smooth on each side with one-sided first derivative traces, global $C^1$ regularity would force both traces to equal the derivative on the interface. Unequal traces therefore preclude a globally classical $C^2$ solution as well. A condition on the function value alone is not a condition on its normal derivative; a specified transmission law or an appropriate weak formulation is additional information.

A diffusion tensor describes direction-dependent diffusive flux. In the local, linear Fickian model,where $C$ is concentration and the entries of $\mathbf{D}$ have the units of a diffusion coefficient. A scalar coefficient is recovered when $\mathbf{D}=D\mathbf{I}$. Material models can instead specify a diagonal matrix or symmetric matrix to represent anisotropic transport.^[ref]

For an original two-dimensional example, chooseA concentration gradient directed along $x$, with components $(g,0)$, produces diffusive flux $(-2g,-g)$ in consistent units. The flux need not be antiparallel to the gradient. Nevertheless,for every nonzero gradient in this positive-definite matrix example.

The eigenvectors $(1,1)/\sqrt{2}$ and $(1,-1)/\sqrt{2}$ have eigenvalues 3 and 1. They identify the principal directions and coefficients. Rotating the coordinate basis changes the matrix entries according to ${\mathbf{D}}^{\prime }=R\mathbf{D}{R}^T$; it does not change these principal coefficients. In an unbounded, homogeneous, drift-free Gaussian model, displacement covariance is $2t\mathbf{D}$, so spreading also differs along those directions.

Symmetry and positive definiteness are assumptions of this ordinary reciprocal diffusion model, not a definition covering every transport law called a tensor. Boundaries also matter: a reflecting boundary condition constrains $\mathbf{n}\cdot \mathbf{D}\nabla C$, which need not be proportional to the normal derivative alone.

For a smooth heat equation solution $w_t=\Delta w$ with $w_t+\alpha {\partial }_{n}w=0$ on the boundary and nonnegative $\alpha$, Green's first identity gives the displayed dissipation identity for $E={\int }_V|\nabla w{|}^2$. The outward normal derivative is used. Writing the boundary loss as $\alpha ({\partial }_{n}w{)}^2$ avoids division by $\alpha$ and remains valid where $\alpha =0$. A stationary solution is harmonic with $\alpha {\partial }_{n}w=0$; it need not be a particular fixed-data Dirichlet problem solution.

With the convention ${\Delta }_z{G}_D=\delta (z-\zeta )$, the displayed Dirichlet Green function vanishes on the unit circle. Indeed $|z-\zeta |=|1-z\bar{\zeta }|$ when $|z|=1$. For $\zeta \ne 0$, the second logarithm corresponds to an opposite image source at $\zeta /|\zeta {|}^2$, outside the unit disc; it is a harmonic function throughout the interior. At $\zeta =0$ the formula becomes $(2\pi {)}^{-1}\log |z|$. Its outward normal derivative in the boundary variable is $(2\pi {)}^{-1}(1-|\zeta {|}^2)/|z-\zeta {|}^2$, the Poisson kernel divided by $2\pi$. Green's second identity therefore gives the Poisson integral on the unit disk.

For equal oriented Dirichlet boundary data on all three sides and a uniquely solvable elliptic problem, rotation invariance makes their outward normal derivatives equal. Pairing the global relation with its Schwarz conjugation of a spectral function at $\ell (k+\lambda /k)=2\pi in$ determines every Fourier coefficient of the common normal trace. The zero mode at $\lambda =0$ follows from zero total flux. At negative parameters, uniqueness and possible removable spectral denominators must be checked rather than assumed.

Given a uniquely solvable Dirichlet problem, this map sends prescribed Dirichlet boundary data to the solution's outward normal derivative. It depends on the differential operator and domain. Spectral global relations offer one way to compute it numerically from boundary traces, while a finite collocation scheme still requires a rank and conditioning check.

Consider the Helmholtz equation $(\Delta +k^2)u=0$ in $z>0$, with $k>0$, boundary value $u(x,0)=g(x)$, and upward outgoing or decaying propagation. Taking the tangential Fourier transform, with frequency $\xi$, givesThus the upward normal derivative is the Fourier multiplier $i\beta (\xi )$. The outward unit normal of the upper half-space points downward, so its Dirichlet-to-Neumann map has multiplier $-i\beta (\xi )$ instead. The normal orientation and the square-root branch determine both signs.

For $|\xi |<k$ the modes propagate; for $|\xi |>k$ they are evanescent waves decaying with height. At high frequency $\beta (\xi )\sim i|\xi |$, so the map grows in magnitude like one tangential derivative. Small-amplitude boundary roughness can therefore have a large derivative effect when it varies on a short spatial scale.

A dynamic boundary condition includes a time derivative of the boundary trace. For a half-line heat equation, $u_t(0,t)+\alpha u_x(0,t)+\beta u(0,t)=f(t)$ gives a separate boundary evolution coupled to the interior flux. Its Laplace transform contains the initial trace $u_0(0)$. The sign of $u_x$ differs from that of the outward normal derivative at a left endpoint.

If a harmonic function on one side of a straight boundary has zero normal derivative there, even reflection extends it harmonically across the boundary. Locally take the boundary to be $x=0$ and assume $u$ and its first derivatives extend continuously to its open segment. Define $\tilde{u}(x,y)=u(|x|,y)$. The zero value of $u_x$ makes the first derivatives match across the join. Applying Green's first identity on the two sides against a smooth compactly supported test function leaves no interface term, so the reflected function is weakly harmonic and hence harmonic by Weyl lemma.

Reflecting across both edges of a right-angle corner gives a harmonic function on a punctured neighbourhood when both normal derivatives vanish. If the original function is bounded at the corner, the removable singularity for a bounded harmonic function makes the extension harmonic there too. Thus a maximum on such an edge or corner can be tested using the interior strong maximum principle for harmonic functions.

A real analytic system solved for the highest normal derivatives, with real analytic Cauchy data, has a unique local real analytic solution. For common order $m$, its normal form is ${\partial }_t^{m}U=F(t,x,\{{\partial }_t^j{\partial }_x^{\alpha }U:j<m,j+|\alpha |\le m\})$. An invertible coefficient/Jacobian matrix for the highest normal derivatives is the system's non-characteristic hypersurface condition. A real analytic coordinate change flattens the initial hypersurface, and adjoining derivatives through order $m-1$ reduces the problem to a first-order real analytic system. Uniqueness here is in the real analytic class; it does not assert Hadamard well-posedness in smooth or Sobolev norms.

A global relation for a linear boundary value problem is an identity connecting transformed initial or boundary traces of a solution. Its form depends on the differential operator, the auxiliary adjoint solutions and the transform convention.

For $L=\Delta -k^2$, suppose $Lu=f$ and $Lv=0$ on a smooth bounded domain $\Omega$. Green's second identity givesThus a homogeneous equation $f=0$ gives a boundary identity with right side zero; an inhomogeneous equation contributes the known volume term. Exponential choices of $v$ turn this into a spectral family connecting the Dirichlet boundary data and the outward normal derivative.

A global relation is a necessary identity. Recovering unknown traces requires enough independent relations and an admissible analytic inversion or a rank-checked numerical scheme; the identity alone is not an existence theorem for arbitrary boundary data.

For a bounded smooth domain $\Omega \subset {\mathbb{R}}^3$ and $u\in C^2(\overline{\Omega })$, let $\Phi (x,y)=1/(4\pi |x-y|)$, normalized by $-{\Delta }_y\Phi (x,y)={\delta }_x(y)$. For an interior point $x$, Green's third identity isIt represents the function through its volume PDE source term, boundary values and outward normal derivative. For a harmonic function the volume term vanishes, so boundary data alone reconstruct the interior.

To derive it, apply Green's second identity on $\Omega \setminus \overline{B_{\epsilon }(x)}$. The inner sphere has the outward normal of the punctured domain, pointing toward $x$. Its singular fundamental solution term tends to $u(x)$ as $\epsilon \downarrow 0$, giving the stated signs. Replacing $\Phi$ by the Dirichlet Green function makes its boundary value zero and leaves the Poisson kernel multiplying the boundary values of $u$. Thus the identity connects Green functions, boundary integral representations and Dirichlet problems.

For the advection-diffusion equation $u_t=u_{xx}+\alpha u_x$ on $x>0$, the spectral divergence form uses $\omega =k^2-i\alpha k$ and factor $e^{-ikx+\omega t}$. With negative-exponential Half-range Fourier transforms, integration gives the displayed global relation for $\mathrm{Im}k\le 0$. The boundary transform $G_1$ refers to $u_x(0,t)$, the negative of the outward normal derivative.

For $L=a^{ij}{D}_{ij}+b^i{D}_i$ with bounded coefficients and uniform ellipticity near a boundary point $y$ satisfying the interior sphere condition, suppose $u\in C^2(\Omega )\cap C^0(\Omega \cup \{y\})$ has $Lu\ge 0$ and $u(x)<u(y)$ in a tangent open ball. With $\mathbf{n}$ the outward unit normal of that ball, ${\mathrm{liminf}}_{t\downarrow 0}[u(y)-u(y-t\mathbf{n})]/t>0$. Consequently the outward normal derivative is positive whenever it exists. The given boundary continuity alone does not assert that the derivative exists. The minimum version follows by applying this to $-u$; operators with a zeroth-order term require additional sign hypotheses.

The Navier slip boundary condition allows a fluid to slide tangentially past a solid wall while resisting that motion through wall friction in the Navier slip boundary condition. For parallel flow past a flat, impermeable wall translating uniformly with tangential velocity $U_w$, write the fluid's tangential velocity as $u$ and its relative tangential velocity as $u_{\mathrm{rel}}=u-U_w$. Thenwhere $b$ is the slip length, also written ${\ell }_s$. The unit normal $\mathbf{n}$ points from the wall into the fluid, and ${\partial }_{n}u=\mathbf{n}\cdot \nabla u$ is the normal derivative. The separate no-penetration boundary condition requires $(\mathbf{u}-{\mathbf{U}}_w)\cdot \mathbf{n}=0$: tangential sliding does not allow flow through the wall.

The sign convention matters. In a channel $0<z<H$, the lower wall has inward unit normal $+{\mathbf{e}}_z$ and the upper wall has inward unit normal $-{\mathbf{e}}_z$. If both walls obey the Navier slip boundary condition, their tangential velocities satisfyHere the prime means the derivative with respect to $z$. Using $+b_H{u}^{\prime }(H)$ at the upper wall would reverse the friction law. The examples Couette flow with a slipping lower wall and planar Poiseuille flow with equal wall slip use different upper-wall boundary conditions.

This flat-wall formula is a local linear model. On a curved wall the Newtonian fluid stress tensor must supply the tangential traction; the normal derivative alone generally misses curvature terms. The author-written review by Lauga, Brenner and Stone states the full strain-based condition with the unit normal directed into the liquid. ^[ref]

The same review credits Navier's 1823 treatise with the linear slip law. Its history is a reminder that the no-slip boundary condition is an experimentally successful model, whose validity depends on the interface and scale. The history here follows that review rather than a new reading of Navier's original memoir. ^[ref]

Neumann boundary data are prescribed values of the outward normal derivative of a function on the boundary of a boundary value problem. They enter a Neumann boundary condition as ${\partial }_{n}u=g_N$. Directional derivatives defined using a fixed coordinate direction must be converted to the outward-normal convention: for example, on the bottom of an upper domain, ${\partial }_{n}u=-u_y$.