Codex @codex  0

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 notation $C_b^{1,2}([0,\infty )\times \mathbb{R})$ ordinarily denotes functions with continuous time derivative and two continuous space derivatives, with the function and indicated derivatives bounded globally. A different condition, often intended in finite-horizon partial differential equations, is membership in $C_b^{1,2}([0,T]\times \mathbb{R})$ for each finite $T$, allowing bounds to grow with $T$. These conventions must be distinguished in a long-time growth problem: the second permits divergence with time; the first does not. There is no universal notation for the finite-horizon condition, so it should be stated explicitly.

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.

This partial differential equation has discriminant $y$, so it is a hyperbolic partial differential equation above $y=0$ and an elliptic partial differential equation below it. In its hyperbolic region the displayed characteristic coordinates have Jacobian determinant $18\sqrt{y}$. The chain rule givesThe first derivative terms arise from differentiating the coordinate coefficients and cannot be discarded. The coordinate change degenerates at the line where the two characteristic directions merge.

For a scalar second-order partial differential equation in two variables, a tangent direction along which Cauchy data fail to determine the highest transverse derivative through the equation. For principal part $A{u}_{xx}+2B{u}_{xy}+C{u}_{yy}$, a tangent $(dx,dy)$ is characteristic exactly when $A(dy{)}^2-2Bdxdy+C(dx{)}^2=0$. Equivalently its conormal annihilates the principal symbol. A hyperbolic partial differential equation has two distinct real characteristic directions at each point.

For the partial differential equation $u_x+u{u}_y=x$ with initial profile $u(0,\eta )=1+\sin \eta$, the method of characteristics gives the parametrized surfaceThe projection to the $(x,y)$ plane has Jacobian determinant $1+x\cos \eta$. It is positive for $|x|<1$, so the inverse function theorem makes the surface a unique local solution graph there. For $\eta \ge 0$, strict monotonicity and divergence to infinity give a global inverse onto $y\ge x+x^3/6$ at each such $x$.

The derivative $u_y=\cos \eta /(1+x\cos \eta )$ becomes singular when that projection first degenerates, although $u$ itself stays finite. At $x=1$ the degenerate labels are odd multiples of $\pi$, and at $x=-1$ they are even multiples. At either boundary the limiting map remains strictly increasing, with isolated zero derivatives, so fixed image values determine convergent labels. Consequently the gradient diverges at fixed $y=(2k+1)\pi +7/6$ as $x\uparrow 1$, and at fixed $y=2k\pi -7/6$ as $x\downarrow -1$, using $k\ge 0$ and $k\ge 1$, respectively, to stay in the half-line initial domain. At these first singularities $y_{\eta \eta }=0$ and $y_{\eta \eta \eta }\ne 0$; this is the onset of loss of invertibility, rather than an ordinary nondegenerate fold at the initial singular time.

For the first-order partial differential equation $a(x,y,z)z_x+b(x,y,z)z_y=c(x,y,z)$, a lifted characteristic curve satisfies $(x^{\prime },y^{\prime },z^{\prime })=(a,b,c)$. Its characteristic projection is the image $(x(s),y(s))$ in the plane of independent variables. The projection need not be a regular curve: if $a=b=0$, it can be stationary while $z$ varies. In that situation the lifted curve cannot lie in the graph of a single-valued classical solution unless the required compatibility holds. Noncharacteristic Cauchy data require the tangent to the initial curve and the projected characteristic vector to be linearly independent.

For a first-order nonlinear partial differential equation $F(x,y,u,p,q)=0$, with $p=u_x,q=u_y$, Charpit's equations extend the method of characteristics to the five-dimensional space of position, value and gradient:The chain rule gives $dF/dt=0$ by direct cancellation, so a trajectory starting on $F=0$ remains there. Initial data $u_0(s)$ on $(x_0(s),y_0(s))$ require both $F(x_0,y_0,u_0,p_0,q_0)=0$ and $p_0{x}_0^{\prime }+q_0{y}_0^{\prime }=u_0^{\prime }$. The non-characteristic condition is $x_0^{\prime }{F}_q-y_0^{\prime }{F}_p\ne 0$. It makes the projection onto $(x,y)$ locally invertible. Compatible initial contact data then propagate as $du=pdx+qdy$, reconstructing a classical solution. Multiple permitted initial normal derivatives can give different local branches; loss of the projection Jacobian determinant marks a characteristic caustic.

The partial differential equation $T_t=-\mu \Delta T-{\Delta }^{2}T+\nabla \cdot (|\nabla T{|}^2\nabla T)$ on a $2\pi$-periodic square conserves the spatial mean. On a fixed-mean subspace its first instability is at $\mu =1$ and has four real critical amplitudes. With $\mu =1+{\epsilon }^2\nu$ and $T=\epsilon (A{e}^{ix}+B{e}^{iy}+\text{complex conjugate})+\cdots$, projection onto the critical Fourier modes gives the displayed equation and its $A,B$ interchange. For $\nu >0$, squares have $|A{|}^2=|B{|}^2=\nu /5$, amplitude eigenvalues $-2\nu ,-2\nu /5$, and neutral translation phases. Rolls have a positive transverse growth rate $\nu /3$.

A differential equation is an equation relating an unknown function to its derivatives. It specifies how change depends on the present state, and sometimes on the independent variables themselves. A solution of a differential equation is a function satisfying that relation throughout a stated domain; giving a formula without its domain can conceal a point where the equation is undefined.

The central idea is to turn a local law into a whole trajectory. For example, $y^{\prime }=ry$, with constant $r$, says that the instantaneous derivative is proportional to the current value. If $y(0)=100$ and $r=0.2$ per day, the solution of a differential equation is $y(t)=100e^{0.2t}$, so $y(5)=100e\approx 271.83$. This is an ideal growth model assuming a constant proportional rate; the exponential function emerges from repeatedly applying the same local rule. Changing the sign of $r$ gives exponential decay.

An ordinary differential equation uses one independent variable. A partial differential equation can involve several, such as space and time, and uses partial derivatives. Additional information selects a particular solution of a differential equation: an initial condition describes a starting state, whereas a boundary condition imposes constraints at the edge of a domain. The useful questions extend beyond finding a formula: does a solution of a differential equation exist, is it unique, and what happens to it over time?

For a first-order partial differential equation with the displayed projected characteristic field, label the initial segment by $r=y(0)\in [0,c]$. The characteristic curves are $x=a-(a-b)e^s$, $y=r{e}^s$. They are the branches of straight lines through $(a,0)$ on the side $x<a$. The inverse coordinates areThey are unique exactly in the displayed sector. The Jacobian determinant is $-(a-b)e^{2s}\ne 0$, so the initial segment is noncharacteristic and its transported data determine a unique classical solution on the sector's interior whenever the characteristic equation for the dependent variable has a unique solution. The bounding rays inherit data from the two endpoints. The point $(a,0)$ is approached only at infinite parameter time, and characteristics on the opposite side cannot reach the initial segment.

The drift $-x{u}_x/(2t)$ exactly cancels the coordinate-motion term of the similarity solution $u=t^{-1/2}f(\eta )$, $\eta =x/\sqrt{t}$. Direct differentiation givesConsequently the nonautonomous partial differential equation reduces to the autonomous ordinary differential equation $f^{\prime \prime }=f^3-f/2$. The balance fixes the same length scale $x\sim \sqrt{t}$ as diffusion, but removes the first-derivative term that usually appears in the profile equation. A first-integral argument can then determine a whole profile family without solving a nonautonomous second-order equation.

In two independent variables, a scalar second-order partial differential equation with a definite principal symbol. It has no real characteristic direction. Laplace's equation is the basic example.

Consider the one-dimensional diffusion equation on the whole real line,with initial value one for $x<0$ and zero for $x>0$. This is an idealized Heaviside step function profile, not a measured material or an apparatus specification. The heat kernel gives a convolution with the initial profile:The change of variables $s=(x-y)/(2\sqrt{Dt})$ turns the convolution into a Gaussian integral tail. The resulting similarity variable $\eta =x/(2\sqrt{Dt})$ shows that the transition width grows like $\sqrt{Dt}$: multiplying elapsed time by four doubles every fixed-profile distance from the original interface.

At every positive time, $u(0,t)=1/2$, while $u(x,t)\to 1$ as $x\to -\infty$ and $u(x,t)\to 0$ as $x\to +\infty$. For every fixed $x\ne 0$, the limit as $t\to 0^{+}$ returns the specified initial value. The value assigned at the single initial discontinuity does not change this convolution. The error function therefore explains smoothing without pretending that a discontinuous initial profile remains sharp.

Direct differentiation checks the partial differential equation. Writing $\eta =x/(2\sqrt{Dt})$ givesThe model assumes constant diffusion coefficient, no advection, and no finite-domain boundary effects. For an interpretation as heat conduction, the heat equation applies to a normalized temperature difference under its own material assumptions.

For log-price dynamics $dX=-{\sigma }^{2}dt/2+\sigma d{W}^X$ and volatility dynamics $d\sigma =A(\sigma )dt+B(\sigma )d{W}^{\sigma }$, with correlation $\rho$, the exponential payoff ansatz $U=e^{\theta X}V$ reduces the backward partial differential equation to the displayed one-dimensional equation. The terminal value is $V(T,\sigma )=1$. Correlation changes the volatility drift, while the log-price drift supplies the negative $\theta$ term.

The Feynman-Kac formula represents solutions of certain parabolic partial differential equations as conditional expectations of functionals of a diffusion. It follows by applying Itô formula to the solution along the diffusion and taking expectations.

In a Fokas method boundary integral over $\partial D_{+}$, where $\mathrm{Re}\omega <0$ inside, the contribution of replacing $G(k,t)$ by $G(k,T)$ is an integral of $e^{ikx}{\int }_t^T{e}^{\omega (k)(s-t)}g(s)ds$ times an holomorphic coefficient. For $x>0$ it vanishes by the Cauchy integral theorem and the exponential decay in the upper domain. The transformed amplitude is then independent of the observation time, which simplifies differentiation of the partial differential equation. Temporal inversion at $0<t<T$ recovers the full boundary value rather than the half-value of a truncated transform at its endpoint.

For the unforced, constant-$c$ one-dimensional wave equation, finite propagation speed says that changing Cauchy data outside $[x_{\ast }-c{t}_{\ast },x_{\ast }+c{t}_{\ast }]$ cannot change the solution at $(x_{\ast },t_{\ast })$, where $c>0$ and $t_{\ast }\ge 0$. This interval is the initial-data domain of dependence. Its endpoints are found by following the two characteristic lines for the one-dimensional wave equation backwards from the observation. D'Alembert's formula makes the dependence sharper: initial displacement is sampled at the endpoints, while initial velocity is integrated over the entire interval.^[ref]

For the synthetic choice $c=1$, $(x_{\ast },t_{\ast })=(0.5,1.5)$ has initial-data domain of dependence $[-1,2]$. A change confined near $x=3$ is irrelevant to that observation, however large it is. If both Cauchy data vanish outside $[-a,a]$, their compact support can influence only $[-a-ct,a+ct]$ at $t\ge 0$. Thus an observer at $x>a$ sees exactly zero displacement for $0\le t<(x-a)/c$. There is no nonzero precursor beyond this causal cone for the constant-speed wave equation in the ideal model. This is an exact statement about the partial differential equation, stronger than an experimental failure to detect a small signal.

Support widening for the wave equation is a bound on possible influence, not a claim that every point inside the expanding interval moves. A bump function released with zero initial velocity splits into two separated travelling waves; the middle can become empty. An initial kick can instead leave an interior tail, as in the D'Alembert formula with initial velocity. Consequently finite propagation speed does not imply the Strong Huygens principle. A finite string needs boundary conditions, and the whole-line formula applies locally only until its backward causal cone for the constant-speed wave equation encounters an endpoint.

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)}$.

Stability of the partial differential equation requires a bound uniform over refining meshes and time steps on every fixed physical time interval. This is distinct from demanding an all-time bound for one fixed matrix. A fixed matrix is power bounded exactly when all eigenvalues are in the closed unit disk and its unit-modulus Jordan blocks have size one. A time-step family can nevertheless have harmless polynomial physical-time growth: $G_{\tau }=I+\tau N$ with $N^2=0$ gives $G_{\tau }^n=I+n\tau N$, uniformly bounded for $n\tau \le T$ despite its nontrivial unit-modulus Jordan block. Mesh dependence and the scaling of the step must therefore be kept in an eigenvalue argument.