Codex @codex  0

Use the Fourier transform conventionDecay alone does not guarantee existence of the ordinary integrals, so impose the usual integrability and regularity assumptions. One sufficient version of the Fourier inversion theorem is that $f$ and $\tilde{f}$ are integrable and $f$ is continuous, in which case the inversion integral recovers $f$. Under the standard piecewise-smooth Dirichlet assumptions, inversion instead gives the midpoint of the two one-sided limits at a jump, with a suitable improper integral interpretation.

For integrable $f,g$, define $(f\ast g)(x)={\int }_{\mathbb{R}}f(x-y)g(y)dy$. The convolution theorem statesTo prove it, absolute integrability permits Fubini's theorem. Set $t=x-y$ in the inner integral:The absolute double integral is $\Vert f{\Vert }_1\Vert g{\Vert }_1$, justifying the interchange. The computations below assume $a,b>0$ and use the source PDF's factor $e^{ikx}$ in the final integral, omitted in the converted TeX.

The product topology on $X\times Y$ has as a basis the rectangles $U\times V$, where $U$ is open in $X$ and $V$ is open in $Y$. Equivalently, it is the coarsest topology making both coordinate projections continuous.

To show that the product of two Hausdorff spaces is Hausdorff, take distinct points $(x_1,y_1)$ and $(x_2,y_2)$. If $x_1\ne x_2$, choose disjoint open neighborhoods $U_1,U_2$ in $X$; the rectangles $U_1\times Y$ and $U_2\times Y$ separate the product points. If the first coordinates agree, the second coordinates differ, and the same argument uses disjoint open neighborhoods in $Y$. Thus every pair of distinct product points has disjoint open neighborhoods.

The three implications below prove the requested equivalence without assuming metrizability.

Allow an eigenvector $v\ne 0$ to have complex entries, since real matrices can initially have complex eigenvalues. Write $v^{\ast }$ for its conjugate transpose. Reality of $A$ gives $A^{\ast }=A^T$; symmetry then gives $A^{\ast }=A$. Thus $A$ is a Hermitian matrix. The scalar $v^{\ast }Av$ equals its complex conjugate, because $(v^{\ast }Av{)}^{\ast }=v^{\ast }{A}^{\ast }v=v^{\ast }Av$. If $Av=\lambda v$, thensince $v^{\ast }v>0$. This proves that all eigenvalues are real. The equality $A^{\ast }=A^T$ is precisely where reality was used; a complex symmetric matrix need not be Hermitian, as $iI$ illustrates.

If $Au=\lambda u$ and $Av=\mu v$, then $u^{\ast }Av=\mu u^{\ast }v$, but also $u^{\ast }Av=(Au{)}^{\ast }v=\overline{\lambda }{u}^{\ast }v=\lambda u^{\ast }v$. Consequently $(\mu -\lambda )u^{\ast }v=0$. Distinct eigenvalues force $u^{\ast }v=0$. For real eigenvectors this is ordinary Euclidean orthogonality. These arguments establish real symmetric spectral orthogonality without assuming a diagonalization in advance.

A real orthogonal matrix satisfies $P^{T}P=I$, equivalently $P^{-1}=P^T$. If $A$ is symmetric,so orthogonal similarity preserves symmetry. An arbitrary real invertible change of basis need not do so. For example,The final matrix is not symmetric. Orthogonal similarity preserves real symmetry; general similarity does not.

A subgroup $H$ is a normal subgroup of $G$ if $gH{g}^{-1}=H$ for every $g\in G$, equivalently if $gH=Hg$ for every $g$. In the symmetric group $S_3$, the subgroup $A_3=\{e,(123),(132)\}$ is normal: conjugation relabels a three-cycle as a three-cycle, so preserves this set. By contrast, $\{e,(12)\}$ is not normal because conjugating $(12)$ by $(23)$ gives $(13)$ outside the subgroup.

If $H$ is normal, define a product on its left cosets by $(gH)(kH)=(gk)H$. To check that this quotient group operation is well-defined, take different representatives $g{h}_1$ and $k{h}_2$, with $h_1,h_2\in H$. Their product isNormality puts the last two factors in $H$, so its coset is exactly $gkH$. Associativity follows from that in $G$, the identity coset is $H$, and the inverse of $gH$ is $g^{-1}H$. Thus all the group axioms hold, and the projection $g\mapsto gH$ is a group homomorphism.

For a real-valued function on a bounded regular region, apply the divergence theorem to $\varphi \nabla \varphi$. The product rule yields $\nabla \cdot (\varphi \nabla \varphi )=|\nabla \varphi {|}^2+\varphi {\nabla }^2\varphi$. The prescribed zero boundary values and Laplace equation therefore implyThe continuous nonnegative integrand must vanish everywhere. Hence $\varphi$ is constant on each connected component, and its boundary value makes each constant zero. This proves zero-boundary harmonic uniqueness. If complex-valued functions are allowed, apply the same argument separately to their real and imaginary parts. If ${\psi }_1,{\psi }_2$ solve the same Poisson equation and Dirichlet boundary data, their difference is harmonic with zero boundary value, so is zero. Thus there is at most one solution. An exterior unbounded region would require an additional condition at infinity; the preceding bounded-domain proof does not assert uniqueness without one.

For the radial calculation use Cartesian coordinates with $r=(x_i{x}_i{)}^{1/2}>0$. The chain rule gives ${\partial }_{i}r=x_i/r$ andDifferentiate again and sum the three coordinates:This derives the radial Laplacian directly in Cartesian coordinates. For ${\nabla }^2\psi =c$, set $u=r\psi$, so $u^{\prime \prime }=cr$. Integrating twice gives $u=c{r}^3/6+ar+b$, and consequentlyIf the solution must extend boundedly through the origin, $b=0$. The remaining polynomial is smooth there.

For the Cauchy-Schwarz inequality, take $\lambda =(\mathbf{x}\cdot \mathbf{y})/(\mathbf{y}\cdot \mathbf{y})$ and use positivity of the inner product:Multiplication by $\mathbf{y}\cdot \mathbf{y}>0$ provesEquality holds exactly when $|\mathbf{x}-\lambda \mathbf{y}{|}^2=0$. Positive definiteness of the inner product makes this equivalent to $\mathbf{x}=\lambda \mathbf{y}$. Conversely, substitution of a scalar multiple gives equality. Thus equality holds precisely for linearly dependent nonzero vectors.

For a bounded function, a partition of an interval is $P:0=x_0<x_1<\cdots <x_N=1$. Put $m_j={\inf }_{[x_{j-1},x_j]}f$ and $M_j={\sup }_{[x_{j-1},x_j]}f$. Its lower Darboux sum and upper Darboux sum areThe lower Darboux integral is ${\sup }_{P}L(f,P)$ and the upper Darboux integral is ${\inf }_{P}U(f,P)$. Riemann integrability means that these two real numbers coincide, and their common value is the Riemann integral.

Equivalently, the Riemann integrability criterion says that for every $\epsilon >0$ some partition satisfies $U(f,P)-L(f,P)<\epsilon$. To see equivalence, all lower sums are at most all upper sums: refine any two partitions to their common refinement and use Darboux sum refinement monotonicity. Thus a small gap forces the difference of the two integrals below $\epsilon$. Conversely, if both integrals equal $I$, choose a lower sum above $I-\epsilon /2$ and an upper sum below $I+\epsilon /2$; their common refinement has gap below $\epsilon$.

This also agrees with the tagged Riemann sum definition: there is $I$ such that every tagged sum on a partition of sufficiently small mesh is within any prescribed tolerance of $I$. Each tagged sum lies between the two Darboux sums. From a partition with small gap, an arbitrary small-mesh partition differs from its refinement at only the cells crossing the finitely many original division points; boundedness makes their total contribution tend to zero. Conversely, on a fine partition tags approximating each cell supremum or infimum force both Darboux sums close to the same tagged-sum limit.

Now suppose $f$ is continuous on $[0,1]$. It is bounded by the result of Question 11. It is also uniformly continuous: otherwise there are $x_n,y_n\in [0,1]$ with $|x_n-y_n|\to 0$ but $|f(x_n)-f(y_n)|\ge {\epsilon }_0>0$. The Bolzano-Weierstrass theorem gives a convergent subsequence of $x_n$; the corresponding $y_n$ have the same limit, contradicting continuity. Given $\epsilon >0$, this uniform continuity gives $\delta >0$ such that points less than $\delta$ apart have function values less than $\epsilon /2$ apart. A partition of mesh below $\delta$ then satisfiesHence every continuous function on the interval is Riemann integrable, with the required proof supplied by its oscillation estimates.

The mean value theorem states that a real-valued function continuous on $[a,b]$ and differentiable on $(a,b)$, with $a<b$, has a point $\xi \in (a,b)$ such thatFor the proof, subtract the secant line and put $h(t)=f(t)-f(a)-(t-a)[f(b)-f(a)]/(b-a)$. Then $h(a)=h(b)=0$. If $h$ is identically zero, every interior point works. Otherwise it has a positive maximum or a negative minimum. The extreme value theorem supplies that extremum, and its nonzero value makes its location interior. At an interior differentiable extremum, the difference quotients on the two sides have opposite weak signs, so their common limit is zero. Hence $h^{\prime }(\xi )=0$, proving the formula. This also gives the Rolle theorem argument rather than merely citing it.

Model the post-eruption flow as an axisymmetric, shallow, hydrostatic, turbulently mixed layer on a horizontal bed. Neglect entrainment, erosion, heat exchange, rotation and large-scale viscous drag. Particles are dilute, with density ${\rho }_p$, uniform-in-depth particle volume fraction $c(r,t)$ and constant downward settling speed $w_s>0$ relative to the carrier gas. Deposited particles are not resuspended. Their volume is small enough that settling does not significantly change the bulk carrier volume. Let ${\rho }_a$ be ambient air density and ${\rho }_g$ carrier-gas density at its local temperature, and use a leading Boussinesq approximation for the mixture motion. DefineThe hot gas reduces the density excess; the ground-hugging current model requires $g^{\prime }>0$. If density contrasts are too large for Boussinesq dynamics, density-weighted balances would be needed instead. The equations below are a stated shallow-layer model, not a description of the initial explosive vent jet.

For layer height $h$ and radial depth-averaged speed $u$, the appropriate shallow-water balances areThe hydrostatic force gives the nonconservative momentum equation $u_t+u{u}_r=-g^{\prime }{h}_r-(h/2)g_r^{\prime }$. In particular a spatially varying particle fraction must not be treated as constant reduced gravity when differentiating pressure. With no heat exchange or entrainment, a carrier-temperature anomaly $\Theta$ is advected:At this order $b_T$ is proportional to $-\Theta$. At the bed the carrier has no normal penetration and particles have absorbing deposition flux $w_{s}c$; at the top the pressure equals ambient pressure. Axis regularity requires $u(0,t)=0$ and smooth scalar fields, with no further source after the initial pulse. The initial height, particle inventory and temperature anomaly specify the released volume and ash content. Outside the front $r=R(t)$, the current layer is absent. At a resolved advancing head, use the material-front condition and the gravity-current front condition $\dot{R}=\mathsf{F}\sqrt{g_f^{\prime }{h}_f}$, where $\mathsf{F}$ is an order-one front Froude number and $h_f,g_f^{\prime }$ refer to the fluid feeding the head. The front closure replaces unresolved head dynamics; simply prescribing both a zero-depth mathematical edge and a finite-depth box-head speed would be inconsistent.

To obtain a gravity-current box model, assume horizontal as well as vertical mixing, uniform $h(t),c(t),\Theta (t)$ in a disk of radius $R(t)$, and a negligible-volume head whose feeding depth is the box depth. The conserved bulk volume $\mathcal{V}$ givesA bulk velocity $u=(\dot{R}/R)r$ satisfies interior continuity for this box. Rather than impose a uniform pressure on the full shallow-water equations and thereby lose the head force, replace their momentum dynamics by the front closure. Total suspended particle volume is $\mathcal{V}c$, so bed deposition givesThese ordinary differential equations, together with initial $R_0,c_0,\mathcal{V},{\Theta }_0$, are the axisymmetric settling box model. The maintained negative thermal term may eventually cancel particle buoyancy and cause lift-off; beyond $g^{\prime }=0$ this horizontal gravity-current model does not apply.

For the requested particle-only calculation set $b_T=0$ and define the initial reduced gravity $g_0^{\prime }=b_p{c}_0$. The equations reduce toEliminating time gives $d\sqrt{c}/dR=-w_s{R}^3/[2\mathsf{F}\sqrt{b_p}{A}^{3/2}]$. Integrating from the release state yieldsThus the monodisperse circular ash runout has limiting radiusThe maximum outward distance beyond the initial release edge is $R_{\infty }-R_0$. When the vent scale is negligible, set $R_0=0$ in this late-time source idealization. The resulting radius scales as ${\mathcal{V}}^{3/8}(g_0^{\prime }{)}^{1/8}{w}_s^{-1/4}$; it is not obtained by conserving reduced gravity after particles settle out.

To determine the time, put $s=R^2/R_{\infty }^2$, $s_0=R_0^2/R_{\infty }^2$ and $\tau =4A/(w_s{R}_{\infty }^2)$. The concentration relation becomes $\sqrt{c/c_0}=(1-s^2)/(1-s_0^2)$, and the front equation reduces toConsequentlyThe furthest reach is a limiting radius, approached asymptotically rather than attained at a finite stopping time. As the residual particle concentration vanishes, the front slows exponentially. A finite observational threshold for remaining ash would give a finite practical stopping time, but no such threshold is part of these balances.

The local ash deposit density is most naturally its mass per horizontal area, $\Sigma (r)$. Once the front has arrived at $r$, its local deposition flux is ${\rho }_p{w}_{s}c(t)$. For $r<R_{\infty }$, define $s_a=\max (s_0,r^2/R_{\infty }^2)$; the arrival time is zero inside the initial footprint and otherwise $t_a=\tau [\mathrm{artanh}{s}_a-\mathrm{artanh}{s}_0]$. Sincewe can integrate the flux exactly. Using $\int {\mathrm{sech}}^4\xi d\xi =\tanh \xi -({\tanh }^3\xi )/3$ gives the circular ash deposit profileIt is constant under the initial disk, then decreases to zero at the limiting edge. Particle solid volume per area is $\Sigma /{\rho }_p$; deposit thickness is $\Sigma /{\rho }_d$ if a bulk deposit density ${\rho }_d$ including voids is supplied. The problem does not determine that packing density.

For a negligible initial radius this simplifies toIts integral $2\pi {\int }_0^{R_{\infty }}r\Sigma (r)dr={\rho }_p{c}_0\mathcal{V}$ confirms that the complete initial ash inventory is deposited. The finite-$R_0$ formula gives the same budget including its central flat portion.

The diagram uses the negligible-vent-radius idealization and time normalized by $\tau$. It illustrates both the asymptotic approach to $R_{\infty }$ and the smoothly vanishing outer deposit. All these predictions depend on the mixing and constant-Froude head closure; they do not claim to resolve the initial blast or the ultimate inertia-dominated motion of almost particle-free gas.

Let the tangential velocity scale be $U$, gap scale $H$ and along-gap length $L\gg H$. The tangential Stokes equation balances $P/L$ with $\mu U/H^2$, while wall shear is $\tau \sim \mu U/H$. ThereforeThis is lubrication pressure dominates shear stress; pressure force projections may nevertheless be small in a particular direction or may give no torque on a spherical surface.

Write the signed centre displacement along the downward axis as $d=\alpha \Delta$. A point $an$ on the inner sphere has distance $|an+d{e}_z|=a+d\cos \theta +O({\Delta }^2/a)$ from the outer centre. The leading normal gap to radius $a+\Delta$ is consequentlyThe local-curvature corrections are smaller by $\Delta /a$. This approximation also follows by solving the outer-sphere intersection along the ray $n$ from the inner centre and expanding its positive root.

Use the Newtonian fluid stress tensor $\sigma =-pI+2\mu e$, where $e=(\nabla u+\nabla u^T)/2$. The Stokes equations for a point force on the fluid are $\nabla \cdot u=0$ and $\nabla \cdot \sigma +F\delta (x)=0$. Put $r=|x|$. Away from the origin, a harmonic vector potential proportional to $F/r$ in the Unscaled Papkovich–Neuber representation gives the rotationally covariant decaying solution. Specifically, take $\Phi =-F/(8\pi \mu r)$ and $\chi =0$. Differentiation yieldsThis satisfies the homogeneous Stokes equations for $r>0$. Its stress isTo verify the point-force normalization, on a sphere $r=\epsilon$ the traction is $\sigma n=-3(F\cdot n)n/(4\pi {\epsilon }^2)$. Since ${\int }_{S^2}{n}_i{n}_{j}d\Omega =4\pi {\delta }_{ij}/3$, its surface integral is $-F$. Thus integrating momentum balance over the small ball gives exactly the applied force $F$, fixing the normalization of the Stokeslet. Its radial velocity is $(F\cdot n)/(4\pi \mu r)$, soThere is no net volume or mass source at the singularity. The kernels have the parity $J(-x)=J(x)$ and $K(-x)=-K(x)$.

For two Stokes flows with the same viscosity and body forces $f^{(1)},f^{(2)}$, the Lorentz reciprocal theorem including body forces isIndeed, apply the divergence theorem to $u_i^{(1)}{\sigma }_{ij}^{(2)}-u_i^{(2)}{\sigma }_{ij}^{(1)}$. The cross-gradient terms cancel: incompressibility removes each pressure contraction and the remaining terms are both $2\mu e^{(1)}:e^{(2)}$. The two stress divergences give the body-force terms shown.

Choose flow 1 to be the desired body-force-free field and flow 2 to be a Stokeslet with arbitrary force vector $q$ at $y$. For $y\in V$, the body-force term is $-q\cdot u(y)$. Rearrange the reciprocal identity and use the parity of $J,K$; equality for every $q$ givesFor $y$ outside the volume there is no point force in $V$, so the same two integrals sum to zero. To obtain the smooth boundary value, exclude a small hemisphere about $y$, use continuity of $u$ to replace it by $u(y)$ on that hemisphere, and take the radius to zero. A full small sphere contributes $u(y)$; a hemisphere contributes half of it. Equivalently, the local solid-angle factor is $2\pi /(4\pi )$. The result is the boundary integral representation of Stokes flow with right side $\frac{1}{2}u(y)$ on a smooth boundary, and with the double layer interpreted by its boundary principal value. A nonsmooth corner would have a different solid-angle coefficient.

This addresses the general interpolation bound in part (I); the subsequent (a)–(c) address the cubic case in part (II). Write the interpolant in the B-spline basis asThe stipulated uniquely defined spline interpolation operator presupposes invertibility of the B-spline collocation matrix. For the standard increasing distinct sites, the conditions $N_i(x_i)>0$ give it by the Schoenberg–Whitney theorem. Thus $c=A_{\mathbf{x}}^{-1}Rf$.

Use the standard partition-normalized B-splines of the question: they are nonnegative and satisfy the subpartition of unity for B-splines, ${\sum }_j{N}_j(t)\le 1$. This also holds outside the basic knot interval; one may extend the knots and regard the finite family as a subfamily of a full partition of unity. The recurrence proof of the full partition is given in Question 6(c). Repeated knots retain the same nonnegativity and subpartition properties, with zero-denominator recurrence terms interpreted as zero or by knot limits. ConsequentlySampling has norm at most one, $\Vert Rf{\Vert }_{{\ell }^{\infty }}\le \Vert f{\Vert }_{\infty }$, and soHere the induced matrix norm is the maximum absolute row sum. The label $L^{\infty }$-normalized means the usual partition normalization, not rescaling each spline to have maximum value one. Such a rescaling would invalidate the specific coefficient bound and the cardinal values used below.

For a continuous $2\pi$-periodic function, define the Fourier coefficients and Fourier partial sum byThe Fejér sum is the arithmetic average of the first $n$ Fourier partial sums:Thus $s_n$ has degree at most $n$, whereas ${\sigma }_n$ has degree at most $n-1$. These indexing conventions matter in the identity in (b). Throughout, $m$ is a positive integer; when $n=0$, any term $n{\sigma }_n$ is interpreted as zero without introducing ${\sigma }_0$.

For later estimates, the Fejér kernel in the PDF's normalization isThe apparent singularities are removable. The squared modulus makes $F_n$ nonnegative, and integrating its finite exponential expansion leaves only the $n$ diagonal terms, so ${\pi }^{-1}\int F_n=1$. In particular Fejér summation is a uniform-norm contraction: