Codex @codex  0

Write $q_i=a^{ij}{D}_{j}u$. If $u$ is a classical solution, multiply $-{\partial }_i{q}_i=f$ by a smooth function $v$ and apply the divergence theorem. The Neumann boundary condition $q_i{N}_i=0$ removes the boundary term and givesThe density of smooth functions in a Sobolev space and boundedness of the coefficients extend this identity to every $v\in H^1(U)$, so $u$ is a weak solution.

Conversely, take $v$ to be a test function compactly supported in $U$. The weak formulation saysThe fundamental lemma of the calculus of variations gives the equation in $U$. Under the regularity implicit in the stated notion of a classical solution, it holds pointwise. Applying integration by parts again with arbitrary $v\in H^1(U)$ leaveswhere $T$ is the trace operator. Traces of smooth functions can be chosen arbitrarily on the boundary, so the boundary fundamental lemma of the calculus of variations gives $(a^{ij}{\partial }_{j}u)N_i=0$. Thus $u$ is a classical solution.

For $u,g\in L_{\mathrm{loc}}^1(U)$, the function $g$ is the $i$th weak derivative $D_{i}u$ whenfor every test function $\phi \in C_c^{\infty }(U)$. This is the integration by parts identity with no boundary term and agrees with the ordinary derivative whenever $u$ is classically differentiable.

For $1\le p\le \infty$, the first-order Sobolev space iswith norm, for example,Functions equal almost everywhere represent the same Sobolev element.

For $u_0\in L^{\infty }(\mathbb{R})$, call $u\in L^{\infty }({\mathbb{R}}_{+}\times \mathbb{R})$ a weak solution when, for every $\phi \in C_c^1([0,\infty )\times \mathbb{R})$,This follows by multiplying the equation by a test function and applying integration by parts in time and space; $(F\phi {)}_x$ includes the term $F_x\phi$ because the original transport operator is not in divergence form.

If $u$ is $C^1$, test functions supported away from $t=0$ show that $u_t+F{u}_x=g$ as a distributional identity, hence pointwise. Integrating this pointwise equation by parts in the weak identity leavesArbitrary boundary test functions and the fundamental lemma of the calculus of variations give $u(0,x)=u_0(x)$. Thus a $C^1$ weak solution is the unique classical solution from part a.

Write the Inviscid Burgers equation in conservation form asA bounded function $u$ is a weak solution with initial datum $u_0$ whenfor every compactly supported $C^1$ test function $\phi$.

Across a straight discontinuity $x=\sigma t$, integration by parts on its two sides shows that the boundary terms cancel exactly when the Rankine-Hugoniot condition holds:when $u_{-}\ne u_{+}$. For every $c>0$, defineThe three jumps have left and right states $(0,-c)$, $(-c,c)$, and $(c,0)$, so their Rankine-Hugoniot speeds are respectively $-c/2$, $0$, and $c/2$, exactly the speeds of the displayed lines. Hence each $u_c$ satisfies the weak equation away from the origin and across every jump. Moreover, its nonzero support at time $t$ has length $ct$, so $u_c(t,\cdot )\to 0$ in $L_{\mathrm{loc}}^1$ as $t\downarrow 0$; its initial datum is therefore zero in the weak identity.

The zero function and all the distinct functions $u_c$ are bounded weak solutions with the same zero initial datum. Thus weak solutions are not unique. The central jump from $-c$ to $c$ is an expansion shock, which an entropy condition would exclude.

For the p-energythe first variation in the direction $\psi \in C_c^{\infty }(\Omega )$ isAn integration by parts therefore gives the Euler-Lagrange equationwhich is the p-Laplacian equation. In the notation of the question one takes $a^{ij}(x,z,\eta )=|\eta {|}^{p-2}{\delta }^{ij}$.

For $\eta \ne 0$, the principal coefficient matrix isIts eigenvalue in directions orthogonal to $\eta$ is $|\eta {|}^{p-2}$, while its eigenvalue parallel to $\eta$ is $(p-1)|\eta {|}^{p-2}$. The coefficients are $C^{1,\alpha }$ away from $\eta =0$, and the condition number there is at most $p-1$. On every region where $0<m\le |Du|\le M$, this gives uniform ellipticity with constants depending on $m$, $M$, and $p$. At $Du=0$ all principal eigenvalues vanish, so the operator is degenerate there and is not strictly elliptic on a domain containing a critical point.

First choose a test function $v\in C_c^{\infty }(U)$. The weak formulation and integration by parts giveThe fundamental lemma of the calculus of variations implies $-\Delta u=f$ pointwise because $u\in C^2(\overline{U})$ and $f$ is continuous. The Dirichlet boundary condition on ${\Gamma }_1$ already follows from $u\in V$ and continuity of $u$.

For arbitrary $v\in V$, Green's first identity and the interior equation now reduce the weak identity toThe traces of smooth members of $V$ can be chosen freely on compact subsets of ${\Gamma }_2$. Another application of the fundamental lemma of the calculus of variations, now on the boundary, gives ${\partial }_{\nu }u=0$ pointwise on ${\Gamma }_2$. Hence $u$ is a classical solution of the complete mixed boundary value problem.

For every compactly supported $C^1$ test function $\phi$ on $[0,\infty )\times \mathbb{R}$, define a weak solution by the identityThe extra $\phi$ appears because ${\partial }_x(x\phi )=x{\phi }_x+\phi$. This identity is obtained from the linear transport equation by integration by parts in time and space.

Conversely, if $u$ and $u_0$ have the stated $C^1$ regularity, choosing test functions supported away from $t=0$ shows in the distributional sense that $u_t+x{u}_x=0$. Continuity makes the equation pointwise. Integrating that pointwise equation by parts in the displayed identity leavesfor all boundary test functions. The fundamental lemma of the calculus of variations gives $u(0,x)=u_0(x)$, so $u$ is a classical solution.

Because $v$ is a divergence-free vector field, the product rule givesThe integral of a divergence over the periodic domain vanishes. Hence integration by parts yieldsThis is the skew-symmetry of incompressible transport.

Let $q=2p+2$ and multiply the vorticity equation by ${\omega }^{q-1}$, first for an odd power as in the hint. Periodic integration by parts and $\nabla \cdot u=0$ giveThe diffusion term is nonnegative, while the Holder inequality bounds the right-hand side by $\Vert g{\Vert }_q\Vert \omega {\Vert }_q^{q-1}$. ConsequentlyOn the finite-volume torus, $L^q$ norms increase to the essential supremum as $q\to \infty$. Hence the damped-vorticity maximum estimate gives

For $I_1$, the change of variables $t=x^2/\theta$ givesFor $k\ge 1$, use $d(e^{-x^2/\theta })=-(2x/\theta )e^{-x^2/\theta }dx$ and integration by parts:The boundary term vanishes: $x^{k+1}\to 0$ at zero, and the exponential decay dominates the power at infinity. These estimates also establish convergence of the improper integrals used above.

Let $w$ be the difference of two sufficiently regular solutions of the wave equation with the prescribed data. By linearity, $w$ has zero initial displacement, zero initial velocity, and homogeneous Dirichlet boundary conditions. Define its wave energy byFor $\alpha =0$, differentiation under the integral sign and integration by parts yieldThe boundary term vanishes because $w=0$ at both fixed endpoints implies $w_t=0$ there. The zero initial data give $E(0)=0$, so $w_t=w_x=0$ for all times. Thus $w$ is spatially constant, and its zero boundary values give $\boxed{w=0}$. This proves uniqueness for the undamped wave equation.

Choose the positive weightsThese are the volumetric heat capacities. The heat equations and integration by parts giveThe exterior boundary terms vanish because the temperatures there are zero. At the thermal interface, the remaining term isby continuity of temperature and heat flux. Hence the weighted energy estimate isBoth weights are positive, so $I\ge 0$. Zero initial temperatures give $I(0)=0$, forcing $I(t)=0$ for all $t>0$. Each temperature therefore vanishes throughout its slab. Applying this argument to a difference also proves uniqueness for arbitrary common initial data in this thermal interface problem.

For separation of variables, put $U=X(x)H(t)$ and set $X^{\prime \prime }/X=H^{\prime }/(\kappa H)=-\lambda$. The boundary conditions are $X(0)=0$, $X^{\prime }(L)=0$. Multiplying $-X^{\prime \prime }=\lambda X$ by $X$ and using integration by parts givesA nonzero solution therefore has $\lambda >0$: the zero case forces $X$ constant, and then $X(0)=0$ forces $X=0$. Writing $\lambda ={\mu }^2$, the eigenfunctions and their temporal factors areThese sine functions are orthogonal on $[0,L]$, with squared norm $L/2$. Expanding the initial condition and multiplying by $\sin ({\mu }_{m}x)$ gives, using the permitted interchange of the sum and integral,Here $\cos ({\mu }_{n}L)=0$ and $\sin ({\mu }_{n}L)=(-1{)}^n$. Thus the Fourier series solution isThe exponential decay permits differentiation of the series for $t>0$, and each term satisfies both boundary conditions. At $t=0$ it is the mixed-end Fourier expansion of $-T_{0}x/L$.

The gained heat is the increase in internal energy. Since the paper specifies no cross-sectional area, first give the answer per unit cross-sectional area; multiply by $A$ for a rod of cross-sectional area $A$. The initial temperature has integral $T_{0}L/2$, whereas ${\int }_0^L\sin ({\mu }_{n}x)dx=1/{\mu }_n$. At $t_{\ast }=L^2/\kappa$ this givesNumerically $\Delta Q/A\simeq 0.45623855\rho c{T}_{0}L$. This is positive and less than the eventual gain $\rho c{T}_{0}L/2$, as expected from relaxation towards the uniform steady state.

For separation of variables, set $u=X(x)Y(y)$, with $X^{\prime \prime }/X=-Y^{\prime \prime }/Y=-\lambda$. The side Neumann boundary conditions give $X^{\prime }(0)=X^{\prime }(\pi )=0$. As in the standard Sturm-Liouville problem, integration by parts shows $\lambda \ge 0$. The zero eigenvalue gives constant $X$, and $Y^{\prime \prime }=0$ together with $Y^{\prime }(\pi )=0$ makes $Y$ constant too.

For positive eigenvalues, the side conditions give $\lambda =n^2$, $X_n=\cos (nx)$, $n\ge 1$. The top condition chooses the hyperbolic cosine factor $\cosh (n(\pi -y))$. Normalize it to one at $y=0$, obtainingThe orthogonality of integer Fourier modes and the permitted interchange of summation and integration giveFor the remaining Fourier coefficients, integrate by parts twice. The sine boundary term vanishes, and with $f=x(\pi -x)$,The last integral vanishes because $f^{\prime \prime }=-2$ and $n$ is a positive integer. Thus $a_n=-2(1+(-1{)}^n)/n^2$: only even Fourier modes survive, as expected from reflection symmetry about $x=\pi /2$. Writing $n=2m$ givesThe factor involving hyperbolic cosine lies between zero and one on the closed rectangle. Comparison with $\sum m^{-2}$ proves uniform convergence, so the solution is continuous on the entire boundary and takes the prescribed bottom values. On every strip $y\ge \epsilon >0$, the Fourier coefficients have additional exponential decay, allowing differentiation and verifying Laplace's equation and all three Neumann boundary conditions.

To justify the endpoint calculation in b, differentiate at $0<x<\pi$, $y=0$ in the interior-edge sense:Subtracting $1$ from the hyperbolic tangent leaves an absolutely convergent series. For the remaining part, the power series $-\log (1-z)={\sum }_{m\ge 1}{z}^m/m$ gives, by taking real parts at $z=r{e}^{2ix}$ and then $r\uparrow 1$,The same formula follows as the limit of the differentiated interior solution, since its leading multiplier is $e^{-2my}$. This establishes the boundary normal derivative used in b. It also explains why the failure of the corner compatibility condition for mixed boundary data affects first derivatives at the bottom corners without spoiling continuity of $u$ itself.

The homogeneous linear differential equation $y^{\prime \prime }+y=0$ has general solution $C_1\cos x+C_2\sin x$. For a particular solution, use variation of parameters with $y_1=\cos x$, $y_2=\sin x$ and Wronskian one:The double-angle formulas and integration by parts giveSubstitution simplifies to $-x^2\cos x/4+x\sin x/4+\cos x/8$. Absorb the last term into the homogeneous linear differential equation solution. ThereforeIndeed, applying $d^2/d{x}^2+1$ to $x\sin x$ gives $2\cos x$, and applying it to $x^2\cos x$ gives $2\cos x-4x\sin x$; their difference divided by four is the required $x\sin x$.

For the rate-one exponential distribution, the probability density function is $e^{-t}$ for $t\ge 0$. By integration by parts,Thus

On $(-\pi ,\pi )$ the function is odd, so its mean and all its cosine Fourier coefficients vanish:The assigned endpoint value does not affect these integrals. The product of the two odd functions $x|x|$ and $\sin nx$ is even, hence the sine Fourier coefficients areTwo applications of integration by parts giveThereforeThe periodic extension has a jump discontinuity at every odd multiple of $\pi$: its left and right limits are ${\pi }^2$ and $-{\pi }^2$. Its assigned value at each such point is $-{\pi }^2$. Thus it is odd away from these points, but it is not a globally odd function pointwise. The Fourier series will take the midpoint value zero at these jumps.

Subtract the constant steady state profile: set $v=T-T_0$. ThenThe last equation follows from Fourier's law and zero heat flux at the insulated end. For a nonzero separated solution $v=X(x)H(t)$, separation of variables givesTo see that $\mu >0$ for a nonzero eigenfunction, multiply the spatial equation by $X$ and use integration by parts:Both boundary terms vanish. If the right side were zero, $X$ would be constant and $X(0)=0$ would make it zero. Thus write $\mu ={\lambda }^2$ with $\lambda >0$. Solving gives $X=A\sin (\lambda x)$, and $X^{\prime }(L)=0$ requires $\cos (\lambda L)=0$. The mixed Dirichlet-Neumann modes on an interval therefore haveThe eigenvalues of $-d^2/d{x}^2$ are ${\mu }_n={\lambda }_n^2$; ${\lambda }_n$ are their positive square roots. By superposition, with coefficients for which the resulting series and its required derivatives converge at positive times,The eigenfunctions are orthogonal on $[0,L]$ and have squared norm $L/2$. Expanding the initial excess temperature $v(x,0)=T_{0}x/L$ givesThus the requested temperature isAt the insulated endpoint $\sin ({\lambda }_{n}L)=(-1{)}^{n-1}$ cancels the alternating coefficient sign, leavingEvery summand is positive. The endpoint temperature never drops below $T_0$; it decreases towards $T_0$ as $t\to \infty$ and starts at $T(L,0)=2T_0$.

For completeness, the Sturm-Liouville eigenfunction expansion gives the initial profile in the $L^2$ sense. Here the coefficients are absolutely summable, so the Weierstrass M-test also gives uniform convergence on $[0,L]\times [0,\infty )$, and the initial values are recovered everywhere. For positive times the exponential factors justify all required termwise derivatives, giving a classical solution. There is, however, an initial incompatibility with a Neumann boundary condition at $(L,0)$: the initial slope is $T_0/L$, while the insulated Neumann boundary condition requires zero slope for $t>0$. The temperature is continuous to the initial edge, but $T_x$ cannot extend continuously to that corner. One must not differentiate the initial series term by term at $x=L$.

The function $x\sin (px)$ is an even function, so every Fourier sine coefficient is zero. Introduce the auxiliary integralFor a nonzero integer $k$, integration by parts givesConsequentlyFor $n\ge 1$, a product-to-sum formula yieldsIf $n\ne p$, this becomesThe exceptional Fourier coefficient at $n=p$ must be evaluated before dividing by $p-n$:Thus the full Fourier series isThis equality holds on $[-\pi ,\pi ]$ and by periodicity everywhere. Indeed the periodic extension is continuous at the joined endpoints, since both endpoint values are zero, and is piecewise smooth. The theorem on pointwise convergence of a piecewise smooth Fourier series therefore applies, including at the endpoints. The Fourier coefficients also decay as $n^{-2}$, making the resulting series absolutely and uniformly convergent.

Parametrize the surface of revolution byAn oriented vector surface area element isUse the vector field $\mathbf{F}=(x/2,y/2,0)$, whose divergence is one. If this parametrization has the outward orientation on an enclosing surface, the divergence theorem givesThe second equality is integration by parts, using the endpoint values. For example this outward convention applies when the positive-$t$ branch lies above the negative-$t$ branch, $f(t)>f(-t)$ for $0<t<\pi /2$.

The stated $C^1$ regularity and endpoint conditions alone do not guarantee an embedded enclosing surface: branches can cross or coincide. Without an appropriate orientation assumption the flux formula is an algebraic signed volume. If “inside” is defined as the region between the two graph branches for each radius, its nonnegative geometric volume is insteadThis explicitly handles a reversal or change in the order of the two branches.