Codex @codex  0

An additive composite expansion combines an inner expansion and an outer expansion by subtracting their common behavior in the overlap region:The subtraction removes the part counted twice. Near the boundary layer, the outer expansion cancels the overlap contribution; away from the boundary layer, the inner expansion cancels it. An error estimate is needed to establish accuracy across the whole domain.^[ref]

For the singular perturbation problem $ϵy^{\prime \prime }+y^{\prime }=0$, $y(0)=0$, $y(1)=1$, the outer solution is one, the leading inner expansion is $1-e^{-x/ϵ}$, and their overlap is one. Thus $y_{\mathrm{comp}}=1-e^{-x/ϵ}$. With $q=e^{-1/ϵ}$, comparison with the exact normalized solution givesThe supremum norm error is therefore exponentially small on the entire interval, even though the outer solution alone has error one. This is an error bound for function values; errors in derivatives have additional powers of ${ϵ}^{-1}$.

The composite satisfies the differential equation exactly but has right boundary condition $y_{\mathrm{comp}}(1)=1-q$. Dividing it by $1-q$ restores that boundary condition and recovers the exact solution. This example makes both the strength and the limit of matching explicit: it gives a global approximation, while finite-$ϵ$ normalization can still carry information beyond every algebraic order.

An adjoint eigenfunction is a nonzero eigenfunction of the adjoint operator on its adjoint domain. For an inner product conjugate-linear in its first argument, $L^{\ast }g=\overline{\lambda }g$ implies $⟨g,(L-\lambda )f⟩=0$ for every $f$ in the domain of $L$. Consequently a necessary solvability condition for $(L-\lambda )f=h$ is $⟨g,h⟩=0$. If $L-\lambda I$ is a Fredholm operator, orthogonality to every element of its adjoint kernel is also sufficient. Boundary conditions and weights in a coupled system enter the adjoint domain and must be derived by integration by parts; copying the original eigenfunction without checking them can give a wrong projection.

An advanced Green function has support only before its source time. It is selected by a future boundary condition, rather than an initial causal response. For $L={\partial }_t^2+2\lambda {\partial }_t+{\omega }_0^2$ with $\lambda <0$ and ${\Omega }^2={\omega }_0^2-{\lambda }^2>0$, the real-axis inverse Fourier transform of $(-{\omega }^2-2i\lambda \omega +{\omega }_0^2{)}^{-1}$ is $-H(-t)e^{-\lambda t}\sin (\Omega t)/\Omega$. Its first derivative jumps by one at zero. A growing retarded solution exists with a different inverse contour; this illustrates causality versus stability of a Green function.

Suppose a real simple leading eigenfunction is $y_0=C\varphi$ and a first-order correction satisfies $L_{0}v=-\nu w{y}_0-y_0^p$, where $p\ge 2$ is an integer, $L_0$ is self-adjoint with null space spanned by $\varphi$, the boundary conditions are homogeneous and independent of the parameter, and $w>0$. The Fredholm solvability condition givesIndeed, orthogonality of the forcing to $\varphi$ gives $\nu C\int w{\varphi }^2+C^p\int {\varphi }^{p+1}=0$, and $C\ne 0$ permits division. Thus the numerical eigenvalue correction depends on the chosen leading amplitude; fixing an eigenfunction normalization is part of defining the nonlinear branch. This calculation assumes a perturbative branch and does not by itself establish that branch's existence.

For a real one-dimensional potential, a real-axis transmission probability is a modulus square. Its algebraic complex continuation need not have only physical bound-state poles. In a square well, $t=e^{-ika}/D(k)$ while the continued probability is $1/[D(k)D(-k)]$. Zeros of the second factor can encode a different boundary condition. In units ${\hslash }^2/(2m)=1$, a well of width two and depth $1+{\cot }^{2}1$ has a continued-probability pole at $k=i\cot 1$ which satisfies neither decaying even nor decaying odd matching. The physical bound-state criterion uses the outgoing amplitude on its decaying sheet.

If a thin capacitor interior is modeled with axially uniform axisymmetric fields, Maxwell's equations give $r^{-1}(r{E}_z^{\prime }{)}^{\prime }+(\omega /c{)}^2{E}_z=0$. The solution regular at the axis is proportional to $J_0(\omega r/c)$. Normalizing to a central sinusoidal voltage gives $E_z=-(V_0/L)J_0(\omega r/c)\sin \omega t$ in this ideal interior model. It can change sign across the radius when propagation matters. Feed, edge and radiation boundary conditions are needed for a full driven-device solution.

Because one wave travels backward, its incident boundary condition belongs at the opposite end of the medium. With $u(0)=a$ and $v(L)=0$, differentiating once gives $u^{\prime \prime }=-{\kappa }^{2}u$ and the solutionThe invariant $|u{|}^2+|v{|}^2$ is constant along $z$, while the input/output powers obey transmission minus reflection equal to the incident power in this normalization. A proposed solution must satisfy both first-order equations, not just their second derivatives.

For a plate under uniform tension per length $T$, with surface mass density $\sigma$ and flexural rigidity $D_b>0$, the small-displacement equation is $\sigma z_{tt}=T{\nabla }^{2}z-D_b{\nabla }^{4}z$. If the edge conditions admit a transverse sine mode of wavenumber $q_n$, its longitudinal dispersion relation isBending raises mode cutoff frequencies of a membrane waveguide and changes large-$k$ behaviour from $\omega \propto |k|$ to $\omega \propto k^2$. A genuinely clamped plate also fixes the edge slope, so its transverse mode shapes need not remain simple sine functions. The qualitative increase in restoring force survives that change of boundary conditions.

This zero-range boundary condition replaces a short-range interaction by its S-wave scattering length. If $\psi =A/r+B+O(r)$, it requires $B=-A/a$. An outgoing spherical wave also has a constant $ik$ term; retaining that term gives the unitary scattering-length amplitude $f_0=-a/(1+ika)$. The model neglects effective-range corrections and is applied outside the microscopic range, not literally inside a finite potential.

For $L=d^2/d{x}^2-2/x^2$ on $x>0$, the homogeneous solutions $x^2,x^{-1}$ satisfy respectively the zero and infinity boundary conditions. Their Wronskian is $-3$. The displayed Green's function is continuous at $x=\xi$ and its first derivative has jump one, so $L_{x}G=\delta (x-\xi )$ in the distribution sense. Integration against a compactly supported forcing gives the unique solution decaying at both boundaries.

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.

Let $h_1,\dots ,h_n$ be a basis of homogeneous solutions of a linear ordinary differential operator $L$ of order $n$. Let $B_1,\dots ,B_n$ be independent homogeneous linear boundary conditions, which may be coupled. Suppose a fundamental solution $G_0(x,\xi )$ satisfies $L_x{G}_0=\delta (x-\xi )$. Define the boundary matrix $M_{ij}=B_i{h}_j$ and $b_i(\xi )=B_i{G}_0(\cdot ,\xi )$. If $M$ is invertible, thenThe added terms are homogeneous solutions, so the Dirac delta source is preserved. Applying the boundary operators gives $b-M{M}^{-1}b=0$. Any two such Green functions differ, for fixed source position, by a homogeneous solution with coefficient vector in the null space of $M$; invertibility therefore proves uniqueness. A singular boundary matrix signals a nontrivial homogeneous solution satisfying the endpoint conditions, where a compatibility condition and a nonunique inverse may arise. For $L=D^2$, one convenient fundamental solution is $(x-\xi )H(x-\xi )$.

For a sufficiently regular function on a domain $\Omega$, its boundary trace is the restriction to $\partial \Omega$, using the limit from within the domain when necessary. In a half-line evolution problem, the traces $f_j(t)={\partial }_x^{j}u(0,t)$ are boundary values of the spatial derivatives. An initial value does not generally determine every such trace without the partial differential equation and suitable boundary conditions. The Sobolev trace operator extends classical restriction continuously to appropriate Sobolev spaces, where pointwise values need not be defined.

The diagrams solve an explicitly bounded Dirichlet problem in a uniform charge-free medium of constant permittivity. Take a width scale $W=1$ and terminal electric potential difference $V=1$. Two zero-thickness ideal plate segments occupy $|x|\le W/2$ at $y=\pm d/2$, with electric potentials $\pm V/2$. Fix $\varphi =0$ on the entire rectangle $x=\pm 2W$, $y=\pm 1.5W$. This boundary is an additional enclosing equipotential, rather than a condition at infinity. The model is invariant along the omitted third direction. The connection between the electric potential, its gradient, and the charge-free Laplace equation follows from electrostatics.^[ref]

On a square mesh of spacing $h=W/160$, the finite difference method imposesat free nodes, while plate and outer nodes retain their prescribed electric potentials. There are $641\times 481$ nodes. Each plotted gap aligns exactly with the mesh; the smallest gap has 16 mesh intervals. The generators solve the resulting linear system by the conjugate gradient method, stopping when the maximum dimensionless five-point residual is below $2\times 10^{-10}$, and use central finite differences to approximate $\mathbf{E}=-\nabla \varphi$ away from plate nodes.

The discrete maximum principle keeps the computed electric potential between $-V/2$ and $+V/2$. Reflection checks require evenness in $x$ and oddness in $y$; halving $h$ tests discretization effects, while enlarging the outer rectangle tests sensitivity to its artificial boundary conditions. These are different checks: a small residual only establishes that the discrete equations were solved. Differentiating near zero-thickness endpoints is especially sensitive to the mesh, so the plots do not report a converged edge peak. No total electric charge, exact open-domain three-dimensional capacitance, or measured material property is extracted from this model.

For a bounded numerical check, compare $h=W/80$ with the plotted $h=W/160$ at midplane positions $x/W=0,0.25,0.40,0.50,0.60,0.80$. At $d/W=0.10$, the largest change in $-E_{y}d/V$ is approximately $0.0265$, at $x/W=0.50$; the other five sampled changes are below $0.00472$. At $d/W=0.25$, moving the outer rectangle from $\pm 2W\times \pm 1.5W$ to $\pm 3W\times \pm 2W$ while retaining $h=W/80$ changes these samples by at most approximately $0.00622$. These calculated sensitivities support the qualitative comparison and expose its limits; they are not experimental uncertainties or evidence that the artificial boundary has converged to infinity.

In an even, piecewise constant potential, integrate outward from the centre with either even boundary conditions ${\psi }^{\prime }(0)=0$ or odd boundary conditions $\psi (0)=0$. Match the wavefunction and its derivative at every finite step and impose exponential decay outside. The resulting discrete energies alternate in parity, with the fourth bound state odd and containing three nodes of a wavefunction. There is no hard-wall zero at a finite step; the second derivative changes according to ${\psi }^{\prime \prime }=2m(V-E)\psi /{\hslash }^2$.

Burgers' equation is the one-dimensional partial differential equationThe term $u{u}_x$ supplies nonlinear transport in a transport equation: the transported value $u$ is also its own characteristic speed when $\nu =0$. Larger values catch up with smaller values ahead of them, steepening a decreasing profile. The diffusion term $\nu u_{xx}$ opposes this steepening. Positive $\nu$ gives the viscous Burgers equation; setting $\nu =0$ gives the Inviscid Burgers equation.

This small equation makes three different behaviors accessible in one model: Burgers characteristic breaking from continuous decreasing data, an entropy shock for a downward step, and a rarefaction wave for an upward step. The Burgers Riemann problem with positive flux and the Viscous Burgers travelling shock below give exact examples. These are mathematical calculations with chosen data; the figures contain no experimental measurements.

Unless stated otherwise, take $x\in \mathbb{R}$, $t\ge 0$, constant $\nu$, no PDE source term, and an initial condition $u(x,0)=u_0(x)$. There are no finite endpoints requiring boundary conditions. Where a profile has different limiting states at $x=\pm \infty$, those states are specified explicitly. On a finite interval, boundary conditions become part of the problem: positive diffusion normally requires a condition at each endpoint, whereas the Inviscid Burgers equation requires appropriate inflow transport boundary conditions. With dimensional variables, $u$ has units of velocity and $\nu$ of kinematic viscosity; all numerical examples below use dimensionless quantities.

Consider the Riemann problem on $\mathbb{R}$ with initial condition $u_0(x)=u_L$ for $x<0$ and $u_0(x)=u_R$ for $x>0$. No endpoint boundary conditions are imposed. When $u_L>u_R$, its entropy solution is a shock wave with constant speed $s=(u_L+u_R)/2$. In the original example $u_L=2$, $u_R=0$, the Rankine-Hugoniot condition givesThe left characteristic speed is $2$ and the right characteristic speed is $0$, so $2>1>0$: both sides enter the entropy shock. The Viscous Burgers travelling shock gives an exact smooth family tending to this solution.

Reversing the states to $u_L=0$, $u_R=2$ still gives a Rankine-Hugoniot condition speed of $1$, but now characteristic curves move away from the discontinuity. That expansion jump is a weak solution and fails entropy admissibility. For general $u_L<u_R$, the selected entropy solution is the continuous self-similar solutionThis rarefaction wave fills the expanding gap between the two exterior characteristic curves. Inside its fan, $u_t=-x/t^2$ and $u{u}_x=x/t^2$, verifying the Inviscid Burgers equation; its values agree continuously with both exterior states. For $0$ to $2$, the fan occupies $0<x<2t$. If $u_L=u_R$, the solution is constant.^[ref]

The capacitor gap aspect ratio compares the gap $d$ with a lateral plate width $W$. For rectangular plates of widths $W_x,W_z$, require both $d/W_x\ll 1$ and $d/W_z\ll 1$ for the leading parallel-plate capacitor approximation. Large area alone is insufficient: a long narrow strip can have large area and still have $d$ comparable with its smaller width. In addition, the observation point must be many gap widths from an edge. Reducing $d/W$ improves the interior approximation without making the edge itself disappear.

In a two-dimensional straight gap, write the electric potential as the infinite-plate solution plus a perturbation:The perturbation satisfies the Laplace equation and vanishes on the two ideal plates. Separation of variables in a semi-infinite strip, with $s$ the distance inward from its entrance, gives transverse modesThus edge disturbances have a natural penetration scale set by $d$. Modes forbidden by a particular symmetry have zero coefficients; this expression is an explanation of localization in an ideal strip, not a universal numerical correction to finite-plate capacitance. The outside boundary conditions determine the coefficients.

There is no single aspect-ratio cutoff guaranteeing an assigned accuracy: the requested observable, plate shape and thickness, surroundings, dielectric distribution, and guard geometry all matter. The small-gap approximation is local and is not uniform across an edge. It also does not justify extrapolating $C\simeq ϵA/d$ to arbitrarily large $d$, where finite electrical conductors and their external electric field dominate. At fixed $V$, formally shrinking an ideal gap makes the interior electric field grow like $1/d$; real electric currents across an imperfectly insulating gap, surface roughness and contact eventually invalidate that ideal limit.

A change of dependent variable expresses the unknown function of a differential equation in terms of a new unknown, while keeping the independent variable fixed. For example, substituting $y=e^{-ax/2}v$ in $y^{\prime \prime }+a{y}^{\prime }+by=f$ removes the first derivative and gives $v^{\prime \prime }+(b-a^2/4)v=e^{ax/2}f$ when $a,b$ are constant. Transforming the boundary conditions as well is essential, particularly for a Robin boundary condition.

Forwith the boundary conditions $y(0)=y^{\prime }(0)=0$, $y^{\prime \prime }(L)=0$, and $-A{y}^{\prime \prime \prime }(L)=F$, the higher-order Euler-Lagrange equation isIts solution splits as