Codex @codex  0

At $a=\pi$, the nonzero homogeneous solution $e^{-x}\sin x$ satisfies both zero boundary conditions. Thus the boundary value problem has a nontrivial kernel and cannot have an inverse for arbitrary forcing.

More directly, the adjoint operator has the zero-endpoint homogeneous solution $w(x)=e^x\sin x$. If an ordinary Green's function existed, the integration by parts identity would implyBut $w(\xi )\ne 0$ for $0<\xi <\pi$. This contradiction rules out the Green's function for the unmodified operator. It is the Fredholm alternative obstruction.

For a sufficiently regular solution, put $w=-\Delta u$, with $\Delta$ the Laplacian. Then $-\Delta u-w=0$ and $u-\Delta w=u+{\Delta }^{2}u=f$. The zero Dirichlet boundary conditions on $u$ and $\Delta u$ become zero Dirichlet boundary conditions on $u$ and $w$. Conversely, eliminating $w$ from these two equations recovers the fourth-order equation with the biharmonic operator. For the weak solution we will construct, these interior identities hold as distributional identities; zero Sobolev traces are imposed on $u$ and $w$, without assuming classical fourth derivatives of $u$.

Let $V=H_0^1(\Omega )\times H_0^1(\Omega )$, the Cartesian product of two zero-boundary Sobolev spaces, with normThis is a Hilbert space. Applying integration by parts separately to the two equations gives the weak formulationEquivalently, for $U=(u,w)$ and $T=(v,z)$, define the bilinear form and linear functionalWe seek $U\in V$ with $b(U,T)=F(T)$ for every $T\in V$.

The skew-coupled elliptic system has a useful bounded bilinear form estimate with constant one. Indeed, write the integrand as the dot product of $(\nabla p,\nabla q,-q,p)$ and $(\nabla v,\nabla z,v,z)$. The Cauchy-Schwarz inequality givesFor equal arguments, the two zeroth-order coupling terms cancel. If $K$ is the constant in the supplied Poincaré inequality, thenThus $b$ is a coercive bilinear form, although it is not a symmetric bilinear form. The Lax-Milgram theorem applies and yieldsIn particular, $\Vert (u,w){\Vert }_V\le (1+K)\Vert f{\Vert }_0$. Testing with compactly supported smooth test functions identifies $-\Delta u=w$ and $-\Delta w=f-u$ as distributional identities. This gives the asserted relation to the original boundary value problem at the available regularity.

Choose the positive $z$ direction along the sphere's laboratory velocity $v\hat{\mathbf{z}}$. In its rest frame the far-field fluid velocity is $-v\hat{\mathbf{z}}$. Therefore the far-field potential is $-vr\cos \theta$. The impermeable sphere imposes ${\partial }_r\varphi =0$ at $r=R$. The monopole coefficient vanishes because the sphere is not a source or sink, and the dipole coefficient follows from the wall condition:ThusThe regular exterior boundary value problem has no forced higher multipoles; uniqueness fixes their coefficients to zero. Differentiating givesHere $u_{\theta }$ is the polar velocity component. An equivalent Cartesian form isOnly no-penetration, not no-slip, is appropriate to this ideal inviscid solution. The sphere density ${\rho }_0$ does not enter the velocity field for an externally prescribed constant speed.^[ref]

Here a prime denotes $d/d\sigma$. It is useful to retain terms through $O({ϵ}^2)$ in the normal mechanical equilibrium equation, even though only first-order shape corrections are sought. The undeformed elastic membrane is tensionless, so a first-order change of curvature contributes to normal force only at second order.

The geometric equations written in reference arc length areTheir first-order Taylor expansions giveand henceThe principal stretches and the first-order constitutive equation areBecause $t_s=t_{\phi }=t$ exactly, the tangential mechanical equilibrium equation gives ${\partial }_{s}t=0$. Thus both $t_1$ and $t_2$ are spatial constants: $t_1^{\prime }=t_2^{\prime }=0$.

Expand the geometric factors in the normal balance:ConsequentlyMatching the first-order coefficient gives the requested relationSince the prescribed pressure has no second-order term, and ${\theta }_1^{\prime }=-r_1^{\prime \prime }$, its second-order coefficient givesA closed boundary value problem for the three requested unknowns is thereforeThe first equation has order two, the second order one, and the last order one, so the four boundary conditions match the total order. The unknown constant $t_2$ is selected by the axial endpoint constraint, rather than by an additional specified traction.

For $p\ne 0$, an equivalent formulation eliminates $z_1$ and displays that constraint as an integral:After solving this system, recover $z_1$ by integration from $z_1(0)=0$. No second derivatives of the unspecified function $f$ are required: they enter the next-order shape corrections through the second-order constitutive equation, whereas the displayed normal balance already determines the leading shape. The case $p=0$ has the undeformed solution $r_1=z_1=t_2=0$.

For a static base state under dead loading, put $v_i=u_i(\xi )e^{i\omega t}$ and take zero incremental body force and traction. The admissible squared frequencies are the eigenvalues of the weighted boundary value problemMultiply by $\overline{u_i}$ and integrate over the reference body. The boundary term in integration by parts vanishes, giving the elastic normal-mode energy identityFor a nonzero mode and positive mass density, the denominator is real and positive. The numerator is real by major symmetry of the incremental elastic moduli. Explicitly, if $u=p+iq$ with real $p,q$, the two imaginary cross terms cancel, leaving the sum of the quadratic energies of $\nabla p$ and $\nabla q$. Consequently every admissible ${\omega }^2$ is real. A negative value, when allowed by indefinite moduli, corresponds to exponential growth or decay rather than an oscillation with real frequency.

If the moduli quadratic form is strictly positive on every nonzero real matrix, that same decomposition shows that the numerator is positive whenever $\nabla u$ is nonzero. Thus ${\omega }^2>0$ for every nontranslation mode. Under the usual assumptions of a bounded connected body, regular coefficients and uniform positivity, the Rayleigh quotient defines a self-adjoint operator with discrete positive squared frequencies on the mass-weighted mean-zero displacement space. Its eigenfunctions give the nontrivial small harmonic oscillations.

There is an unavoidable zero-frequency qualification for the pure traction problem as printed: $u_i=b_i$, for any nonzero constant vector $b$, has zero gradient, zero incremental nominal stress and ${\omega }^2=0$. These three constant translations obey all the homogeneous incremental boundary conditions. Under the strict quadratic-form hypothesis they are the entire zero eigenspace on a connected body. Therefore the unrestricted conclusion isStrict positivity for every nonzero displacement requires excluding these modes, for example by imposing ${\int }_{{\Omega }_0}{\rho }_{0}udV=0$, or by adding a displacement constraint. If harmonic perturbations are defined to have nonzero frequency from the outset, the requested positivity holds for all those perturbations.

The fixed-displacement boundary value problem now has $u=0$ on the whole boundary. The elastic normal-mode energy identity still holds, because its boundary term contains $\overline{u}$. DefineFor real constant $K$, $B$ has major symmetry: both Levi-Civita symbols change sign upon exchanging the two index pairs. Its quadratic density is a quadratic minor null Lagrangian. Indeed,The boundary integral vanishes by the Dirichlet boundary condition, and the volume integral vanishes because mixed derivatives are symmetric in $\alpha ,\beta$ whereas ${ϵ}_{\alpha \beta \gamma }$ is antisymmetric. Constancy of $K$ is essential: otherwise its derivatives would produce extra terms. This calculation initially applies to smooth fields and extends to fields with zero Sobolev boundary trace by density.

It follows that the Rayleigh quotient numerator can be replaced by the quadratic form with $H$:The assumed positivity of $H$ on all nonzero real matrices also gives positivity on nonzero complex matrices by splitting real and imaginary parts. A nonzero $u$ satisfying $u=0$ on the boundary cannot have identically zero gradient. Hence the right side is strictly positive, andThis works even when the original moduli quadratic form is not pointwise positive on all matrices: its integrated energy on fixed-boundary gradients agrees with a positive density after adding the quadratic minor null Lagrangian.

Work at one radius and write $\mathcal{R}=k/({\mu }_m{m}_p)$. Use the surface density to fix the mass density normalization and introduce the height scaleChoose the dimensionless variablesThe pressure and temperature scales make the hydrostatic and ideal gas equations parameter-free. The flux scale fixes the heating coefficient at $9/4$, and the choice of $h$ makes the diffusion coefficient unity. Direct substitution gives the universal ordinary differential equationsFor example, the dimensionless coefficient in the temperature equation is $3\alpha \kappa {\Sigma }^2{\mathcal{R}}^4/(16\sigma {\Omega }^5{h}^6)=1$.

Let the upper free surface be at $\zeta ={\zeta }_s$, so its physical semithickness is $H=h{\zeta }_s$. Midplane symmetry and the radiative-zero disk surface conditions giveThe last condition is precisely $\Sigma ={\int }_{-H}^H\rho dz$. The unknown surface position and the two positive midplane values $P(0),\Theta (0)$ are determined as part of the normalized boundary value problem. Symmetry makes $P,D,\Theta$ even and $Q$ odd; $P^{\prime }(0)={\Theta }^{\prime }(0)=0$ already follow from the equations and are not additional independent boundary conditions.

The zero conditions do not mean zero surface flux. Integrating the heating equation gives $Q({\zeta }_s)=(9/4){\int }_0^{{\zeta }_s}Pd\zeta >0$. Near a regular surface with finite outgoing flux, dividing the hydrostatic equation by the temperature equation yields $dP/d\Theta \simeq {\zeta }_s{\Theta }^3/Q_s$. Thus $P\propto {\Theta }^4$, $D\propto {\Theta }^3$ and $\Theta \propto {\zeta }_s-\zeta$, so mass density also vanishes. The apparent singular quotient in the diffusion equation has a finite regular limit. All boundary conditions and equations are independent of the physical parameters: this is homologous vertical structure of a constant-opacity alpha disk.

Use the pressure scale $\rho {\hat{U}}^2$, and denote the scaled transverse velocity by $v$. The steady incompressible flow equations under the prescribed coordinate scaling areThe no-slip boundary condition is $u=v=0$ at $y=1$ and at the lower wall $y=ϵF(x)$, taking $F=0$ upstream. Far upstream, $u\to U_0=6y(1-y)$, $v\to 0$, and $p_x\to -12\lambda /R$. These are the dimensional Navier-Stokes equations divided by the streamwise and transverse inertial scales, respectively; the transverse equation's ${\lambda }^{-2}$ is essential.

A consistent nonlinear wall-layer regime can be chosen asNear a wall the oncoming plane Poiseuille flow is linear in distance, so a layer of thickness $ϵ$ has $u=O(ϵ)$ and $v=O({ϵ}^2)$. Its inertia is $O({ϵ}^2)$; its transverse viscous term is $(\lambda /R)O({ϵ}^{-1})=O({ϵ}^2)$, explaining the distinguished balance. The resulting pressure perturbation is $O({ϵ}^2)$. In the core the induced transverse pressure variation is $O(ϵ/{\lambda }^2)$, negligible compared with that pressure perturbation when ${\lambda }^2ϵ\gg 1$, equivalently $\lambda \gg R^{1/7}$. Also $\lambda \ll R$ makes both the indentation and the viscous boundary layer thin. Streamwise viscous derivatives are smaller still. Thus the core perturbation is inviscid at first order and pressure is common to the two wall layers at their leading order.

Write $u=U_0+ϵu_1+\cdots$ and $v=ϵv_1+\cdots$. The core pressure gradient only enters at order ${ϵ}^2$, so first-order continuity equation and streamwise momentum giveEliminating $u_{1x}$ shows $(v_1/U_0{)}_y=0$. Let that ratio be $-A^{\prime }(x)$; integration in $x$, with the unperturbed upstream condition, then givesThis displacement form automatically has zero first-order change in total core flux, because ${\int }_0^1{U}_0^{\prime }dy=0$.

For the lower boundary layer, usewhere $p_0^{\prime }=-12{ϵ}^3$ is the smaller upstream pressure gradient. The extra term in $v$ removes the velocity induced by the moving coordinate surface. The chain rule gives ${\partial }_x{|}_y={\partial }_x{|}_z-F^{\prime }{\partial }_z$ and ${\partial }_y={ϵ}^{-1}{\partial }_z$; the wall-slope terms cancel from both continuity equation and convection. At leading order the Prandtl boundary-layer equation isMatching to the core near the lower wall gives $U\sim 6[z+F+A]$, so $H_{\mathrm{lower}}=F+A$.

For the upper boundary layer, set $y=1-ϵz$, $u=ϵU$, $v=-{ϵ}^{2}V$. The same equations, wall data and common $P(x)$ result, but the core matching gives $U\sim 6[z-A]$, so $H_{\mathrm{upper}}=-A$. The two pressure-driven wall-layer problems have the same upstream shear, wall conditions and pressure forcing. Under the question's uniqueness assumption their solutions coincide, and their displacement limits must coincide. Thus the symmetric viscous interaction in an indented Poiseuille channel gives

For $F=b{x}^{1/3}$, one has $A=-(b/2)x^{1/3}$ and $H=(b/2)x^{1/3}$. Matching $U\sim 6[z+H(x)]$ in a similarity solution $U=x^{\beta }{G}^{\prime }(z{x}^{-\gamma })$ requires $\beta =\gamma =1/3$. Equating the powers in convection $x^{2\beta -1}$, viscosity $x^{\beta -2\gamma }$ and pressure gradient $x^{\sigma -1}$ gives $\sigma =2/3$. HenceUse the stream function $\psi =x^{2/3}G(\eta )$, with $\eta =z{x}^{-1/3}$. ThenDirect differentiation gives $U{U}_x+V{U}_z=x^{-1/3}[(G^{\prime }{)}^2/3-2G{G}^{\prime \prime }/3]$, whereas $-P^{\prime }+U_{zz}=x^{-1/3}[-2\bar{p}/3+G^{\prime \prime \prime }]$. The one-third-power similarity of an indented channel boundary layer is therefore governed byThe original PDF has $(G^{\prime }{)}^2$ in this equation; the converted TeX loses the derivative and gives $G^2$. Existence remains conditional on a value of $\bar{p}$ for which this boundary value problem has a solution, as stipulated in the question. The similarity describes the stated thin-indentation regime on $x>0$.

Take $L=-d^2/d{x}^2+x$. The Airy equation is $Lu=0$. To discuss positivity of this linear operator, use its homogeneous variation domain $D(L)=\{v\in H^2(0,1):v(0)=0,v^{\prime }(1)=0\}$. Integration by parts givesThe boundary term is zero at both ends. The same calculation with two different functions proves that the associated form is a symmetric bilinear form. Positivity here concerns the homogeneous domain; the condition $u(0)=1$ instead defines an affine set of admissible solutions.

For the variational problem defineand minimizeThe Sobolev trace at zero is meaningful in $H^1$. No derivative boundary value is imposed on this trial space: the Neumann boundary condition will emerge naturally. For $v\in V$, Cauchy-Schwarz inequality and integration giveso $a(v,v)\ge \frac{2}{3}\Vert v{\Vert }_{H^1}^2$. Thus $a$ is a bounded, symmetric, coercive bilinear form on the Hilbert space $V$. The Lax-Milgram theorem says that a bounded coercive bilinear form and a bounded linear functional determine a unique weak solution. Apply it toand set $u=1+v$. This is precisely the first variation condition $a(u,w)=0$ for $J$.

Compactly supported test functions first yield $u^{\prime \prime }=xu$ in the sense of distributional derivatives. Since $xu\in L^2$, the solution is in $H^2$. Integration by parts then leaves $u^{\prime }(1)w(1)=0$ for every $w\in V$, hence $u^{\prime }(1)=0$. The other boundary condition is built into $1+V$. Conversely a solution of this boundary value problem satisfies the first variation condition. Finally,This proves existence, uniqueness and the global minimum characterization, including the natural boundary condition. It is the mixed-boundary Airy energy principle.

For the Dirichlet Green's function, solve $-D{G}_{xx}=\delta (x-\xi )$ with zero endpoint values. Away from $\xi$, $G$ is linear. Continuity at $\xi$ and the derivative jump $G_x(\xi +,\xi )-G_x(\xi -,\xi )=-1/D$ giveFor nonzero Dirichlet boundary conditions, subtract their linear interpolant. The resulting boundary value problem has solutionIn the final heat equation, a steady state has zero time derivative. Here $l=1$, $f(x)=x$, $\alpha =1/D$, $\beta =2/D$. Integrating $u^{\prime \prime }=-x/D$ twice and imposing the endpoint values givesIts second derivative and both boundary conditions verify the answer directly.

Substitute part (b) into the spectral representation. The terms containing $\psi (ik)$ areThe transform $\psi (ik)$ is analytic for $\mathrm{Im}k>0$. Write $k=a+ib$. In the first quadrant,which decays since $x,y>0$. The Cauchy integral theorem in that quadrant equates the two integrals in the first parentheses. In the second quadrant,which decays since $a<0$, $b>0$, and $y<\ell$. The same contour integration equates the two integrals in the second parentheses. Therefore $\psi (ik)$ makes no contribution. This is upper-quadrant cancellation of reflected boundary transforms, and it uses the whole paired expression rather than attempting to discard an individual integral.

The terms involving $c$ cancel by the same two quadrant arguments with the analytic factors $k\pm \gamma$ and $\psi$ omitted. There is a sign defect in the printed reconstruction prefactor. The sides specified in part (a) form a clockwise boundary, so the Cauchy integral theorem gives a negative reconstruction sign. Indeed, integrating each spectral ray first yields $i/(z-\zeta )$, and hence the sum of the three ray integrals is $i{\oint }_{\mathrm{clockwise}}{q}_z(\zeta )/(z-\zeta )d\zeta =-2\pi q_z(z)$. Thus the printed positive prefactor reconstructs $-q_z$. The requested cancellation above holds with either overall sign, but the actual data-dependent derivative isDecay at infinity fixes the integration constant when recovering $q$ from its Wirtinger derivative.

For completeness, the missing condition at infinity and the unrestricted printed $\gamma$ deserve an explicit solvability of a decaying Robin strip check. Let $Y_n$ be orthonormal eigenfunctions from a Sturm-Liouville problem of $-d^2/d{y}^2$ with $Y_n^{\prime }(0)=\gamma Y_n(0)$ and $Y_n^{\prime }(\ell )=-\gamma Y_n(\ell )$, and eigenvalues ${\lambda }_n$. With $g_n={\int }_0^{\ell }g{Y}_{n}dy$, separation of variables givesIndeed, each coefficient obeys $q_n^{\prime \prime }-{\lambda }_n{q}_n=0$ and $q_n^{\prime }(0)=g_n$. A positive eigenvalue has one decaying exponential; a zero eigenvalue has only affine solutions and a negative eigenvalue only oscillatory solutions, neither of which decays unless its coefficient vanishes. This also proves uniqueness in the decay class and justifies reflection symmetry there. For $\gamma >0$, all eigenvalues are positive. For $\gamma =0$, the constant eigenfunction requires ${\int }_0^{\ell }gdy=0$; the apparent zero of $D$ at the origin is then removable. For $\gamma <0$, compatibility with every nonpositive eigenmode is necessary. These restrictions cannot be inferred from smoothness and reflection symmetry alone. The printed boundary value problem without any far-field condition permits additional growing or nondecaying homogeneous solutions, whereas the spectral transforms select the decaying branch.

Fix the convention $\Delta q-m^{2}q=0$, $m>0$, and let $c_m=m^2/4$. With $z=x+iy$, two families of modified Helmholtz adjoint plane waves areBoth satisfy $\Delta E_{\pm }=m^2{E}_{\pm }$, since the products of their $z$ and $\bar{z}$ exponents are $m^2/4$. They obey $E_{-}(z,k)=\overline{E_{+}(z,\bar{k})}$. Green second identity gives the two conjugate global relations for the modified Helmholtz equationFor real boundary traces, ${\mathcal{G}}_{-}(k)=\overline{{\mathcal{G}}_{+}(\bar{k})}$. Thus the second identity is the conjugate spectral companion, rather than an unrelated extra boundary condition. In complex differential form the first relation is equivalentlyIts integrand is a closed one-form: differentiating its coefficients gives $2E_{+}(q_{z\bar{z}}-c_{m}q)=0$. This supplies a direct derivation of the global relation as well as its outward-normal sign convention.

Write the known Dirichlet boundary condition on side $j$ as $f_j(s)$ and the unknown outward normal derivative as $h_j(s)$, with $-1<s<1$. For any linear combination $v$ of the adjoint waves, the global relation isExpanding the $h_j$ in a finite basis and enforcing these identities at selected collocation points for a global relation gives a linear system with a known right side. The square permits especially simple paired tests.

Let $p_n=n\pi /2$, ${\eta }_n=\sqrt{m^2+p_n^2}$, and ${\phi }_n(s)=\sin [p_n(s+1)]$ for $n\ge 1$. Define four real adjoint solutionsThey satisfy the modified Helmholtz equation because ${\eta }_n^2-p_n^2=m^2$. Each is a linear combination of two adjoint exponentials, hence is obtained from the two spectral global relations. More explicitly, put ${\rho }_n^{\pm }=({\eta }_n\pm p_n)/2$, so ${\rho }_n^{+}{\rho }_n^{-}=c_m$. For the $E_{+}$ convention above, the required spectral pairs areAppropriate phase-weighted differences produce $\sin [p_n(s+1)]$. For example the top test is $e^{-{\eta }_n}[e^{i{p}_n}{E}_{+}({\rho }_n^{-})-e^{-i{p}_n}{E}_{+}({\rho }_n^{+})]/(2i)$. The conjugate relation ensures a real system for real data. This is sine collocation of square modified Helmholtz global relations.

The functions ${\phi }_n$ are an orthonormal Fourier sine basis on $[-1,1]$. Set $h_{j,n}={\int }_{-1}^1{h}_j{\phi }_{n}ds$ and compute the known quantity $R_{j,n}={\int }_{\partial \Omega }f{\partial }_n{v}_{j,n}ds$ by numerical integration. Since ${\phi }_n(\pm 1)=0$, each adjoint test vanishes on both adjacent sides. On its own side it equals ${\phi }_n$, and on the opposite side it equals ${\delta }_n{\phi }_n$, where ${\delta }_n=e^{-2{\eta }_n}$. Thus the square modified Helmholtz Dirichlet-to-Neumann coefficients obeyEach block is inverted explicitly: for opposite sides $j,j^{\prime }$,Compute these coefficients for $1\le n\le N$ and reconstruct $h_j^{(N)}={\sum }_{n=1}^N{h}_{j,n}{\phi }_n$. This is a semi-analytical scheme: the spectral basis integrals and two-by-two inverses are explicit, while the known boundary integrals are evaluated numerically. Increase $N$ and the quadrature resolution until the desired convergence is observed; compatible smooth side data have convergent normal-trace expansions, while corner singularities require the usual weaker trace interpretation and more careful quadrature.

The requested weak interaction between sides is particularly clear: the adjacent-side unknown traces contribute exactly zero, and the opposite-side coefficient is $e^{-2\sqrt{m^2+(n\pi /2{)}^2}}\to 0$. The own-side coefficient stays equal to one. Thus diagonal dominance of paired square global-relation collocation is genuine after pairing; individual uncombined exponential samples need not have the same conditioning. The eigenvalues of each block are $1\pm {\delta }_n$, so its condition number is $(1+{\delta }_n)/(1-{\delta }_n)$.

There is also a direct localization check for an individual adjoint wave. On side $j$ with outward unit normal $n_j$ and tangent $t_j$, the normalized test $e^{\eta (n_j\cdot r-1)+ip{t}_j\cdot r}$ has $\eta =\sqrt{m^2+p^2}$. Its opposite-side magnitude is $e^{-2\eta }$, and its magnitude integrates to at most $1/\eta$ along either adjacent side. This side localization of modified Helmholtz plane waves explains why suitable large spectral parameters suppress remote-side effects even before the exact sine cancellation.

Finally the interior solution can be evaluated from the computed traces using the two-dimensional modified Helmholtz fundamental solutionHere $K_0$ is the Modified Bessel function of the second kind, and $(-\Delta +m^2){\Gamma }_m=\delta$. Replace $h$ by its computed finite Fourier sine series and use numerical integration; at interior points the boundary kernels are smooth. The sign follows from Green second identity with this fundamental-solution convention. This completes the numerical integration of the boundary value problem with a Dirichlet boundary condition, rather than stopping at an equation for its unknown boundary derivatives.