Boundary value problem

A boundary value problem asks for a differential-equation solution satisfying conditions at more than one point or boundary component.

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.

Find functions analytic on the two sides of a contour whose boundary values obey a prescribed jump, with suitable normalization at infinity and other distinguished points. An additive jump $J$ on a counterclockwise closed contour is solved by the Cauchy integral formula $F(k)=(2\pi i{)}^{-1}\int J(\zeta )/(\zeta -k)d\zeta$, subject to the normalization and compatibility conditions. The Sokhotski–Plemelj theorem verifies the jump.

The Wiener-Hopf method solves a boundary value problem with complementary half-line traces by Wiener-Hopf factorization and additive splitting into upper- and lower-half-plane analytic functions. Analytic continuation gives a common entire function, whose growth and edge conditions determine its admissible value.

The plus factor is analytic and nonzero above the transform contour, and the minus factor below it. Their continuations inherit singularities and zeros in the opposite half-planes. Allocation of branch cuts and modal poles must agree with the physical outgoing continuation; no factor may have a zero inside its designated analytic domain.

A scalar Wiener-Hopf equation couples two Half-range Fourier transforms analytic in complementary half-planes. After dividing by a kernel factor, additive splitting separates the known forcing. Source poles, boundary terms and contour prescriptions are part of the problem data.

A forcing pole can be split between analytic half-planes by subtracting a factor's value at the pole. For instance,The second quotient has a removable pole and belongs to the plus expression; the first keeps the prescribed minus-side forcing pole. After Wiener-Hopf factorization, analytic continuation identifies an entire remainder, whose value is fixed by growth and edge conditions. The pole prescription must be stated together with the Fourier transform convention.

The Fokas method uses a global relation for a linear boundary value problem between transformed boundary traces, complex spectral symmetries, and contour integration to eliminate unknown traces or construct a numerical boundary scheme. The distinction between the domains of analyticity of the transforms and the decay sectors of spectral exponentials controls valid contour deformations.

To eliminate an unknown boundary transform by a dispersion relation symmetry, evaluate the global relation at $\nu (k)$ only where the spatial transform remains analytic. Substitute the resulting trace identity into the contour representation. An unwanted transform vanishes only if its transformed argument and exponential admit a valid contour deformation and decay estimate. For constant drift, $\nu (k)=i\alpha -k$ maps the upper domain $\mathrm{Re}(k^2-i\alpha k)<0$ into the lower transform half-plane.

If $P(k)$ is analytic in the upper half-plane with suitable growth, paired integrals along the positive real and imaginary rays cancel when their exponential is $e^{ik(x+iy)}$ with $x,y>0$. For $0<y<\ell$, the analogous pair along the imaginary and negative real rays cancels with $e^{ik(x+i(y-\ell ))}$. These two applications of the Cauchy integral theorem eliminate reflected unknown transforms from a strip global relation; an individual ray integral need not vanish.

A global relation is an identity connecting transforms or integrals of the boundary traces of a solution. For a formally self-adjoint operator $L=\Delta -k^2$, Green second identity givesfor every homogeneous adjoint solution $Lv=0$. Exponential adjoint families turn this into a spectral identity between the Dirichlet boundary condition and unknown normal derivative.

A boundary trace transformed over a finite time interval using the dispersion relation of a linear partial differential equation. For integrable trace data and a polynomial dispersion function, it is an entire function of the complex spectral parameter. Its growth depends on $\mathrm{Re}\omega (k)$; multiplying by $e^{-\omega (k)t}$ changes the integrand to the causal factor $e^{-\omega (k)(t-s)}$.

For $u_{z\bar{z}}={\beta }^{2}u$, pull back its spectral closed differential one-form along each straight side $z=m_j+s{h}_j$. Under counterclockwise traversal and outward normal derivatives $q_j$, the side density is $i{e}^{-i\beta (\lambda z-\bar{z}/\lambda )}[|h_j|q_j+\beta (\lambda h_j+{\bar{h}}_j/\lambda )g_j]$, where $g_j$ is the Dirichlet boundary data. The Generalized Stokes theorem gives the displayed global relation. Clockwise traversal reverses the normal-derivative coefficient; all orientations must be changed consistently.

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.

These are chosen complex values of a spectral parameter for a linear boundary value problem at which a truncated global relation for a linear boundary value problem is enforced. Their directions and magnitudes determine which boundary side and which tangential mode are tested. Suitable row combinations can give much better conditioning than arbitrary raw samples.

The two modified Helmholtz adjoint plane waves $v_{\lambda }$ and ${\tilde{v}}_{\lambda }=\overline{v_{\overline{\lambda }}}$ give boundary identities $\mathcal{G}(\lambda )=0$ and $\tilde{\mathcal{G}}(\lambda )=0$. For real solution traces, $\tilde{\mathcal{G}}(\lambda )=\overline{\mathcal{G}(\overline{\lambda })}$. Reality is needed for this conjugation argument; the parametrization separately yields the reciprocal identity $\tilde{\mathcal{G}}(\lambda )=\mathcal{G}(1/\lambda )$.

On the square $[-1,1{]}^2$, write normal traces in ${\varphi }_n(s)=\sin [n\pi (s+1)/2]$. Set $p_n=n\pi /2$, ${\omega }_n=\sqrt{k^2+p_n^2}$, $r_n=({\omega }_n+p_n)/k$. Paired spectral samples $\{r_n,r_n^{-1}\}$, $\{-r_n,-r_n^{-1}\}$, $\{i{r}_n,i/r_n\}$, $\{-i{r}_n,-i/r_n\}$ give adjoint tests $e^{\mp {\omega }_{n}y}{\varphi }_n(x)$ and $e^{\mp {\omega }_{n}x}{\varphi }_n(y)$. Each test vanishes on the adjacent sides, removing their unknown normal derivatives.

Let $R_{jn}={\int }_{\partial \Omega }f{\partial }_n{v}_{jn}ds$ for a square sine adjoint test, and put ${\rho }_{jn}=e^{-{\omega }_n}{R}_{jn}$, ${\delta }_n=e^{-2{\omega }_n}$. The opposite-side normal-trace coefficients obeywith the same formulas for $L,R$. When the known Dirichlet boundary condition is integrated exactly, these are exact coefficients, since orthogonality eliminates all other unknown modes.

For sine collocation of square modified Helmholtz global relations, each opposite-side mode has a scaled matrixThe assembled system has strict diagonal dominance, and its spectral condition number of a positive-definite matrix is below $(1+e^{-\pi })/(1-e^{-\pi })<1.091$. The conclusion concerns the explicit paired sine rows, not arbitrary uncombined complex rows.

A complex parameter labels auxiliary solutions of a linear differential equation, such as exponential adjoint solutions used in a global relation for a linear boundary value problem. It need not be a discrete eigenvalue or a parameter of a nonlinear Lax pair. Analytic dependence on it allows transformed boundary identities to be sampled or inverted.

Subtract a known function with the required boundary trace so that the remaining unknown has homogeneous boundary conditions. For an interval with Dirichlet boundary conditions $g(t),h(t)$, $\ell =(1-x/L)g+(x/L)h$ is a simple lifting. In $w_t=w_{xx}-cw$, $z=w-\ell$ has forcing $-{\ell }_t-c\ell$ and zero endpoint values. This improves Fourier sine series convergence.

For the free Schrodinger equation with compatible smooth initial and boundary data, subtract $g(t)e^{-x}$ before taking a Fourier sine transform. The lifted initial datum vanishes at the endpoint, making its sine transform $O(k^{-2})$ if its second derivative is integrable. The transformed forcing is $k(ig-g^{\prime })/(1+k^2)$; time integration followed by integration by parts makes its solution contribution uniformly $O(k^{-3})$. An integrable spectral majorant proves uniform convergence even at the compatible initial-boundary corner.