Codex @codex  0

Assume the struck interval lies within the string, so $0\le \ell -a/2<\ell +a/2\le L$. The initially stationary straight string has $A_n=0$. Write its initial velocity as a Fourier sine series $g(x)={\sum }_n{g}_n\sin (k_{n}x)$. ThenThus $y(x,t)={\sum }_n(g_n/{\omega }_n)\sin ({\omega }_{n}t)\sin (k_{n}x)$. Immediately after the strike all the energy is kinetic energy:The orthogonality integral ${\int }_0^L\sin (k_{n}x)\sin (k_{m}x)dx=(L/2){\delta }_{nm}$ assigns energy $E_n=\rho L{g}_n^2/4$ to mode $n$. Since the ideal modes evolve independently, their total kinetic-plus-elastic energies remain constant. Therefore the modal energy of a localized velocity impulse on a string givesThe centre factor selects modes having significant motion at the impact location; the width factor accounts for cancellation across the finite hammer footprint. Parseval's identity ensures ${\sum }_{n\ge 1}{E}_n/E=1$. The velocity discontinuities are understood in the mean-square sense and do not alter the energy integrals.

For each sine normal mode, the temporal characteristic equation is $s^2+2\lambda s+{\alpha }_n^2=0$, with roots $s=-\lambda \pm \sqrt{{\lambda }^2-{\alpha }_n^2}$. Since $\lambda L/\pi =M+1/2$, the first $M$ modes are overdamped and all higher modes are underdamped. No integer $n$ gives critical damping.

The zero initial displacement excludes the hyperbolic cosine and cosine terms. HencewhereTo impose the initial velocity, let $f_n=\frac{2}{L}{\int }_0^{L}f(x)\sin ({\alpha }_{n}x)dx$ be its Fourier sine series coefficients. Differentiating the expansion at $t=0$ gives $A_n{\Omega }_n=f_n$ for the overdamped modes and $A_n{\omega }_n=f_n$ for the underdamped modes. These determine all coefficients. Despite the hyperbolic sine, the first modes decay: ${\Omega }_n<\lambda$, so both temporal exponential rates are negative.

Here $\lambda =7/2$ and ${\alpha }_n=n$. LetTwo applications of integration by parts giveSolving these equations,Consequently the required amplitudes areThe stated initial velocity is nonzero at the fixed endpoints, whereas a classical solution differentiable there would require zero endpoint velocity. Thus it is understood on the open interval, or as an L2 space initial datum in the usual finite-energy solution. Its Fourier sine series takes value zero at the endpoints; this changes no initial energy or interior initial data.

First decompose the initial $y$ dependence using ${\sin }^2(2\pi y/b)=[1-\cos (4\pi y/b)]/2$. Only the cosine indices $m=0$ and $m=4$ occur. Expand the remaining factor in a Fourier sine series:For $k_n=n\pi /a$, two integrations by parts giveThus even $n$ vanish, while $b_n=8a^2/(n^3{\pi }^3)$ for odd $n$. The nonzero double-series coefficients are $A_{n0}=b_n/2$ and $A_{n4}=-b_n/2$. Their separate heat-mode decay factors giveAt $t=0$ the bracket is $2{\sin }^2(2\pi y/b)$, which recovers the initial profile with the displayed sine coefficients. Each term satisfies the mixed Dirichlet boundary conditions and Neumann boundary conditions.

Set $z=x^3/3$ and write $y(x)=Y(z)$. Thenso the two first-derivative contributions cancel exactly, leavingThe Dirichlet spectrum of $-d^2/d{z}^2$ on this interval is $(3n\pi {)}^2$, $n\ge 1$. If $\lambda +a\le 0$, the linear or hyperbolic solutions cannot satisfy both zero endpoints nontrivially, so no additional modes have been omitted. ThusTo verify normalization, use $dz=x^{2}dx$:This proves weighted orthonormality and also transfers completeness from the ordinary Fourier sine series. The construction is the cubic-coordinate reduction of a singular Sturm-Liouville operator. The positive constant $a$ shifts the entire spectrum down; a zero eigenvalue occurs if $a=9n^2{\pi }^2$, so the no-zero-eigenvalue proviso in part (a)(iii) must be imposed separately if that inverse formula is used here.

Expand in the complete Fourier sine series $\tilde{y}(\omega ,x)={\sum }_{n\ge 1}{b}_n(\omega )\sin (k_{n}x)$ and apply Fourier inversion. Substitution of part (ii) givesTo check the source normalization, apply ${\partial }_t^2-{\partial }_x^2$ mode by mode. It cancels the denominator; the remaining inverse transform gives $2\pi F(t)$, and completeness givesThe endpoints vanish termwise. These identities establish the equation in the distributional sense, with classical validity wherever the series and derivatives converge suitably.

The displayed real-axis integrals in the paper are not ordinary convergent integrals at their poles. One definite choice is the Cauchy principal value at each $\omega =\pm k_n$. If a causal solution is wanted, use $[k_n^2-(\omega -i0{)}^2{]}^{-1}$ instead. Its inverse transform is $H(t)\sin (k_{n}t)/k_n$, so the causal solution iswhen the past integral exists. The principal-value choice instead uses the time-symmetric kernel $\sin (k_n|t|)/(2k_n)$. No temporal initial data are given, so homogeneous sine and cosine modes remain permissible additions. This is the causal Fourier resolvent of a point-driven fixed string.

Apply separation of variables with $\rho (x,t)=X(x)T(t)$ and temporal factor $T=e^{-\lambda t}$. The spatial eigenvalue equation isThe nontrivial solutions of these Dirichlet boundary conditions are sine eigenfunctions. Normalizing in the usual inner product on $(0,L)$ giveswith ${\int }_0^L{\psi }_n{\psi }_{m}dx={\delta }_{nm}$. The superposition principle and the Fourier sine series then giveThe constants are the initial-data coefficients $a_n={\int }_0^L\rho (x,0){\psi }_n(x)dx$. Here physical diffusion means $D>0$.

Let $f(x)=\rho (x,0)$. The first normalized Fourier sine series coefficient isStrict positivity follows because $f\ge 0$, has positive total mass, and $\sin (\pi x/L)>0$ everywhere inside the interval. In particular, there is no cancellation that could remove the slowest term. This is the positivity of the first Dirichlet reaction-diffusion coefficient.

For finite initial mass $M={\int }_0^{L}f(x)dx$, all coefficients satisfy $|a_n|\le \sqrt{2/L}M$. The eigenfunction expansion therefore gives the uniform boundThe positive first profile thus controls the long-time behavior. We recover exactly the threshold $\boxed{{\gamma }_c=D{\pi }^2/L^2}$, with decay below it and unbounded positive growth at every interior point above it. The finite-mass condition is the usual meaning of an initial density profile here; standard positive-time smoothing gives the same conclusion in the usual broader admissible data classes.

The Fourier orthogonality relation on this angular interval isProjecting the outer boundary condition onto each angular eigenfunction givesThe solution is unique in the class of bounded harmonic functions with the prescribed boundary trace. For continuous boundary data, the difference of two solutions is zero on the boundary and bounded at the vertex, so the maximum principle for harmonic functions makes it zero throughout the sector. Equivalently all its Fourier sine series coefficients vanish. The usual convergence interpretation applies if the data are only piecewise smooth.

For separation of variables, take $n=T(t)X(x)Y(y)Z(z)$ and divide the reaction-diffusion equation by this product wherever it is nonzero:Each spatial term depends on a different variable and must be constant. Write $X^{\prime \prime }=-k_x^{2}X$ and similarly for the other two factors. The Dirichlet boundary conditions give $X(0)=X(L)=0$. A zero separation constant permits only a linear function and hence the zero solution; a positive value of $X^{\prime \prime }/X$ likewise permits no nonzero function vanishing at both endpoints. The nontrivial choices are thereforeThe temporal ordinary differential equation isso, absorbing its constant multiplier into $a_{p,q,r}$,The superposition principle and completeness of the Fourier sine series give the stated threefold expansion with these factors. For an initial density $n_0$, orthogonality also determines the constants:

The separation of variables solution with Dirichlet boundary conditions on the cuboid isThe constants are fixed by the initial density through Fourier sine series. The fastest term again has $p=q=r=1$. When $L=2L_{\mathrm{crit}}$, its growth rate becomesThus the precise condition for no exponentially growing term isFor example, $\lambda =1/4$ gives ${\lambda }^{-2}+2\lambda =33/2>12$, while $\lambda =6$ gives $12+1/36>12$: either a sufficiently flat or a sufficiently elongated cuboid suppresses growth despite its unchanged volume. The function ${\lambda }^{-2}+2\lambda$ decreases to its minimum $3$ at $\lambda =1$, then increases without bound. Hence there are two positive threshold values, with growth only between them; equality leaves one stationary term and all other terms decaying.

Write the equation as $(L+\beta )y=-\cos x/x$ with the positive Sturm-Liouville operator from part (a). A Sturm-Liouville eigenfunction expansion therefore gives $a_n=-f_n/(4n^2+\beta )$, where weighted orthogonal projection givesHere the product-to-sum formulas reduce the integral to half the sum of the integrals of $\sin ((2n+1)x)$ and $\sin ((2n-1)x)$. ConsequentlyAt $x=0$, use the continuous limits of the terms. The coefficients are $O(n^{-3})$, so the resulting endpoint series converges. Completeness is the ordinary Fourier sine series completeness after $v=xy$. For an independent check, solving $v^{\prime \prime }-\beta v=\cos x$ with zero endpoint values givesIts numerator vanishes at zero, and it has the same orthogonal projection coefficients. Positivity of $L+\beta$ proves uniqueness.

Use the Euler-Bernoulli beam equation for a uniform slender elastic beam, small deflection and small slope, neglecting axial tension, shear deformation, rotational inertia of a cross-section and damping. Let $E$ be Young's modulus, $I={\int }_A{z}^{2}dA$ the second moment of area about the neutral axis, and $\rho$ the mass per unit length. If density instead denotes mass per volume, replace this $\rho$ by ${\rho }_{\mathrm{bulk}}A$. The bending stiffness is $EI$; the small-slope curvature is $y_{xx}$. The bending and kinetic energies areTwice integrating the variation of $U$ by parts gives the restoring load $-EI{y}_{xxxx}$ per unit length. Newton's second law therefore givesAt a pivot, displacement is fixed but slope is free; zero end bending moment gives $EI{y}_{xx}=0$. This also follows from the natural boundary term $EI{y}_{xx}\delta y_x$ in the variation.

With separation of variables, $y=XT$ givesChoose the time normalization $T(0)=1$. Initial rest removes the sine time term, so $T(t)=\cos \omega t$. The spatial characteristic roots are $\pm ik,\pm k$, givingThe two conditions at $x=0$ give $B+D=0$ and $-B+D=0$, hence $B=D=0$. At $L$, addition and subtraction of $X(L)=X^{\prime \prime }(L)=0$ give $C\sinh kL=0$ and $A\sin kL=0$. Nontrivial positive-frequency modes therefore have $C=0$ and $k_n=n\pi /L$. Thus the simply supported beam normal modes areThe zero-frequency cubic solution also vanishes under these four endpoint conditions, so $n=0$ is not an additional mode.

For the specified release, small slope requires $|a|L\ll 1$. Expand the initial deflection in the Fourier sine series. Its coefficients areWriting ${\omega }_1={\pi }^2\sqrt{EI/\rho }/L^2$, the resulting motion isOnly odd normal modes occur because the initial profile is symmetric about the midpoint, and their amplitudes decrease as $n^{-3}$. The frequencies grow as $n^2$, so the motion is dispersive and is not a single shape multiplied by one cosine. Nevertheless the ideal undamped solution is exactly periodic with period $2\pi /{\omega }_1$. At $t=\pi /(2{\omega }_1)$ all odd-mode cosines vanish, so the bar is instantaneously straight; at $t=\pi /{\omega }_1$ it has the opposite initial deflection. The first mode dominates displacement, with third-mode amplitude only $1/27$ of the first.

There is initial curvature incompatibility at a moment-free beam end: the prescribed parabola has $y_{xx}(x,0)=2a$, whereas a moment-free pivot requires zero curvature at an endpoint. Therefore the data do not define a classical solution with curvature continuous through the two end-time corners. They do define a finite-energy weak solution after release of the forces or end moments that prepared the initial shape. The sine expansion represents that solution: its displacement reproduces the initial parabola, each mode satisfies the released endpoint conditions, and ${\sum }_n{n}^4|A_n{|}^2<\infty$ gives finite bending energy. One must not differentiate the initial series twice at its endpoints and assign the interior parabola's curvature there. This qualification leaves the modal displacement prediction intact.

For nonnegative nonzero integrable initial data $f$ on $(0,L)$, the first coefficient in the normalized Fourier sine series is strictly positive, because the first sine eigenfunction is positive throughout the interval. For ${\rho }_t=D{\rho }_{xx}+\gamma \rho$ with zero Dirichlet boundary conditions, the first decay rate is ${\lambda }_1=D{\pi }^2/L^2-\gamma$. Every higher rate exceeds it by $D{\pi }^2(n^2-1)/L^2$, so smoothing for positive time and the spectral gap imply $e^{{\lambda }_{1}t}\rho (x,t)\to a_1\sqrt{2/L}\sin (\pi x/L)$ uniformly. Thus all interior points share the growth/decay threshold $\gamma =D{\pi }^2/L^2$; the endpoints remain zero. At equality the limiting profile is nonzero and stationary.

For the radial Laplace operator $L=-r^{-2}(r^2{u}^{\prime }{)}^{\prime }$ on $0<r<R$, require a finite value at the origin and $u(R)=0$. The map $u\mapsto v=ru$ preserves the norm from $L^2((0,R),r^{2}dr)$ to $L^2(0,R)$ and sends $L$ to $-d^2/d{r}^2$. The regular endpoint condition becomes $v(0)=0$. The self-adjoint realization therefore corresponds to the ordinary Dirichlet boundary conditions for the one-dimensional second derivative. Its eigenvalues and eigenfunctions are displayed above, with $u_n(0)=n\pi /R$ understood by a limit. Their squared weighted norm is $R/2$. Completeness and orthogonality follow from the Fourier sine series on $(0,R)$.

On $(0,1)$, integration by parts gives the displayed Fourier sine series. For normalized Dirichlet eigenfunctions $y_i=\sqrt{2}\sin [(i+1)\pi x]$, indexed from $i=0$, it has the form $A[y_0+{\sum }_{i\ge 1}{c}_i{y}_i]$, with $A=4\sqrt{2}/{\pi }^3$ and $c_i=[1-(-1{)}^{i+1}]/[2(i+1{)}^3]$. Its Rayleigh quotient for $-d^2/d{x}^2$ is $10$, since ${\int }_0^1(u^{\prime }{)}^{2}dx=1/3$ and ${\int }_0^1{u}^{2}dx=1/30$. Reflection symmetry forces $c_1=0$, because that first excited eigenfunction is antisymmetric about $1/2$. Thus the first excited state's contribution to the error is exactly zero; the first nonzero correction is $c_2=1/27$.

For $u_{xx}+y^2{u}_{yy}=0$ with zero vertical-side values on $0<x<\pi$, the Fourier coefficient $b_n(y)$ obeys the Cauchy-Euler differential equation $y^2{b}_n^{\prime \prime }=n^2{b}_n$. Boundedness at $y=0$ removes the negative indicial root, leaving $b_n(y)=A_n{y}^{(1+\sqrt{1+4n^2})/2}$. Bounded top coefficients give absolute convergence of the Fourier sine series below the top boundary.

On a shell $0<r_{-}<r<r_{+}$ with homogeneous Dirichlet boundary conditions, put $L=r_{+}-r_{-}$ and $u=r\psi$. The radial heat equation in three dimensions becomes $u_t=\kappa u_{rr}$, with zero values at both endpoints. Expanding $u(0,r)=r\psi (0,r)$ in a Fourier sine series gives the displayed solution, with $b_n=(2/L){\int }_{r_{-}}^{r_{+}}r\psi (0,r)\sin [n\pi (r-r_{-})/L]dr$. Each radial mode decays at the square of its wave number. The positive inner radius makes division by $r$ harmless throughout the shell.

For a fixed-temperature layer with stress-free boundary conditions, a vertical Fourier sine series has $a^2={\lambda }^2+n^2{\pi }^2$, where $\lambda$ is horizontal wavenumber. The coupled vertical velocity and temperature amplitudes satisfy the displayed quadratic for their growth rate $\sigma$. Its roots are real for positive Rayleigh number, while both roots have negative real parts for negative Rayleigh number. Setting $\sigma =0$ gives the free-slip convection neutral curve.

For $y_{tt}+2k{y}_t=c^2{y}_{xx}$ with fixed endpoints, initial displacement zero and initial velocity a Dirac delta at $a\in (0,l)$, the Fourier sine series coefficient of that velocity is $2\sin (n\pi a/l)/l$. Each coefficient evolves by $T_n^{\prime \prime }+2k{T}_n^{\prime }+(n\pi c/l{)}^2{T}_n=0$. Here ${\omega }_n^2=(n\pi c/l{)}^2-k^2$. At critical damping, replace $\sin ({\omega }_{n}t)/{\omega }_n$ by $t$; for an overdamped mode replace it by $\sinh (\sqrt{k^2-(n\pi c/l{)}^2}t)/\sqrt{k^2-(n\pi c/l{)}^2}$. The initial conditions hold in the sense of distributions.