Codex @codex  0

Let the fluid approach a uniform, quiescent state $({\rho }_0,p_0,c_0)$ at infinity. Use the outward unit normal vector $N=\nabla f/|\nabla f|$, pointing from the object into the fluid. Write $H=H(f)$ for the Heaviside step function, $q=\rho -{\rho }_0$, and ${\delta }_S=\delta (f)|\nabla f|$. With the viscous stress tensor ${\tau }_{ij}$, defineThus $T_{ij}$ is the Lighthill stress tensor, and $L_i$ is the surface loading exerted on the fluid. The total force on the object has the opposite sign. Extend $q$ by zero inside the object. Since the boundary is fixed and impermeable, $u\cdot N=0$, and the distributional mass conservation and momentum conservation equations areFor example, the surface contribution to the momentum equation is $(\rho u_i{u}_j+P_{ij})N_j{\delta }_S=L_i{\delta }_S$. Differentiating the first equation in time and subtracting the divergence of the second gives the fixed-body Ffowcs Williams-Hawkings equation:There is no acoustic thickness noise: a fixed impermeable boundary injects no mass. Apply the retarded acoustic Green function $\delta (t-R/c_0)/(4\pi c_0^{2}R)$, where $R=|x-y|$. The supplied surface-delta identity converts the loading term to a surface integral. For an exterior observer, the required integral equation isThis acoustic analogy is an exact rearrangement of the fluid equations; source approximations enter later. In the distant linear acoustic region, $p^{\prime }=c_0^{2}q$.

The acoustic far field requires observer distance $r$ much larger than the source extent $\ell$ and $k_{0}r\gg 1$, so derivatives of the retarded phase dominate derivatives of $1/R$. The acoustic compact-source approximation requires $k_0\ell \ll 1$, so the source's differential propagation delays can be neglected. These conditions concern different lengths and neither implies the other. Put $n=x/r$, $t_r=t-r/c_0$, and define the loading and stress momentsThe leading compact acoustic far field isThe surface term is an acoustic dipole; the volume term is an acoustic quadrupole. The plus sign of the loading term follows from the minus sign of its spatial divergence and the derivative of retarded time.

For the power comparison, assume low Mach number $M=U/c_0$, high Reynolds number, eddy size $\ell$, characteristic time $\tau =\ell /U$, and stress magnitude ${\rho }_0{U}^2$. Keep the prescribed turbulent region and its leading volume stress comparable when inserting the object; neglect changes of source statistics and exceptional cancellations of the leading moments. Compactness follows from $\ell /(c_0\tau )=M\ll 1$. Since $Q\sim {\rho }_0{U}^2{\ell }^3$ and the outgoing acoustic intensity is $⟨p^{\prime 2}⟩/({\rho }_0{c}_0)$, the volume power isIf the object's size is comparable to $\ell$ and its fluctuating resultant loading is $F\sim {\rho }_0{U}^2{\ell }^2$, thenAngular factors and dimensionless loading coefficients have been omitted. The acoustic dipole and acoustic quadrupole interference vanishes in the leading power integrated over a full sphere, by odd angular parity; individual directions can show interference.

A size threshold also requires a loading model. For a small nonseparating object of size $a\ll \ell$ in an approximately inviscid, smoothly varying incident eddy, a spatially uniform pressure has zero resultant on its closed surface. The first surviving loading is the pressure-gradient or added mass force, $F\sim {\rho }_0{a}^{3}U/\tau ={\rho }_0{U}^2{a}^3/\ell$. Under these explicitly stated assumptions,The two mechanisms become comparable at $a\sim \ell M^{1/3}$; for smaller objects the cases with and without the object have the same leading order of power. This is the small-body loading-noise threshold. It is not a geometry-independent size law: a separated drag force $F\sim {\rho }_0{U}^2{a}^2$ fluctuating on the original eddy time instead gives ${\mathcal{P}}_D/{\mathcal{P}}_Q\sim (a/\ell {)}^4/M^2$ and $a\sim \ell \sqrt{M}$. If both force and time are set by body-scale eddies, $\tau =a/U$, their loading power is ${\rho }_0{U}^6{a}^2/c_0^3$ and comparison with the original volume source gives $a\sim M\ell$. Each estimate states which source time and force it uses.

The physical distinction is momentum conservation. Unforced turbulence in a fluid supplies internal stresses, whose net-force contribution cancels and leaves acoustic quadrupole radiation. The supported object exchanges momentum with the fluid and the external support, permitting a time-dependent resultant force and acoustic dipole radiation. One fewer retarded derivative makes the loading radiation more efficient at low Mach number.

Take the normal vector outward from the cylinder-and-cone solid. The surface consists of the cylindrical side and conical roof, with no bottom disk; this convention fixes the otherwise unspecified orientation of a surface. The boundary is the radius-two circle at $z=-2$. Reversing every normal vector reverses the final sign.

For $\mathbf{F}=(y{z}^2,0,0)$, its curl is $(0,2yz,-z^2)$. For the cylindrical parametrized surface, use ${\mathbf{r}}_c=(2\cos \theta ,2\sin \theta ,z)$, with $-2\le z\le 2$, and use ${\mathbf{r}}_{c,\theta }\times {\mathbf{r}}_{c,z}=(2\cos \theta ,2\sin \theta ,0)$ as the outward area vector. Its flux isFor the cone, let $q=4-z$, $2\le z\le 4$ and ${\mathbf{r}}_k=(q\cos \theta ,q\sin \theta ,z)$. Its outward area vector is ${\mathbf{r}}_{k,\theta }\times {\mathbf{r}}_{k,z}=(q\cos \theta ,q\sin \theta ,q)$. ThusThe total surface integral is therefore $-16\pi$ with the outward convention, or $+16\pi$ for the opposite convention.

For Stokes theorem, the induced bottom-circle orientation is counterclockwise as viewed from above, opposite to the orientation of an outward bottom cap. Use $\mathbf{r}=(2\cos \theta ,2\sin \theta ,-2)$ with increasing $\theta$. Thenagreeing with the two parametrized fluxes. The seam at $z=2$ is internal and cancels between the two pieces; the cone tip does not supply another boundary curve.

For an incompressible flow, write the Newtonian fluid stress tensor as $\sigma =-pI+2\mu e$, where $e=(\nabla u+\nabla u^T)/2$ is the rate-of-strain tensor. The Stokes equation gives ${\partial }_j{\sigma }_{ij}=0$. Symmetry of the Cauchy stress tensor therefore givesIntegrating and applying the divergence theorem expresses the viscous dissipation as boundary power:At a moving rigid body, the surface velocity is $U+\Omega \times (x-x_c)$. Its boundary-power contribution is $U\cdot F+\Omega \cdot G$, with the force and torque evaluated using the fluid's outward normal vector. Both resultants vanish for a force-free, torque-free inclusion. Thus the new viscous dissipation comes entirely from the unchanged outer boundary velocity. Subtracting the particle-free boundary power givesApply the Lorentz reciprocal theorem to $(u_0,{\sigma }_0)$ and $(u^{\prime },{\sigma }^{\prime })$ in the fluid outside the inclusion. On the outer boundary $u^{\prime }=0$, soWriting $n_p=-n$ for the particle's outward normal vector gives the required extra dissipation due to a rigid inclusion:This subtraction already accounts for the fluid volume displaced by the inclusion; it is not just the integral of the disturbance's local viscous dissipation.

For the sphere, the ambient rate-of-strain tensor $E$ is symmetric and trace free. The sphere in a uniform straining Stokes flow has zero translational and rotational velocity: inversion symmetry eliminates its force, and symmetry of $E$ eliminates its torque. To derive the disturbance, put $r=|x|$, $Q=x\cdot Ex$ and seekThe incompressible flow condition and Stokes equation reduce toA decaying family satisfying these equations is $f=A{r}^{-5}$, $g=B{r}^{-5}-5A{r}^{-7}/2$, $j=2B{r}^{-5}$. The no-slip boundary condition requires $f(a)=-1$ and $g(a)=0$, so $A=-a^5$ and $B=-5a^3/2$. ConsequentlyThe total velocity is zero at $r=a$, and the disturbance decays as $r^{-2}$. The pressure constant has been set to zero. These boundary and far-field conditions, together with Uniqueness of Stokes flow, establish the solution. There is no rotational background in a pure strain flow.

On the sphere $u^{\prime }=-u_0=-aE{n}_p$. Contracting the supplied Newtonian fluid stress tensor with $n_p$ makes its two terms proportional to $(n_p\cdot E{n}_p)n_p$ cancel, leaving $\sigma n_p=5\mu E{n}_p$. Hence the extra dissipation due to a rigid inclusion isHere ${\int }_{S^2}{n}_j{n}_{k}d\Omega =4\pi {\delta }_{jk}/3$. Taking the large outer boundary limit only after using the fixed-boundary power identity avoids replacing that prescribed boundary by an uncontrolled boundary at infinity.

If $N_v$ denotes the number of spheres per unit volume, $N_v(4\pi a^3/3)=\varphi$. ThusAdd this to the particle-free viscous dissipation density $2\mu E:E$. Comparing with $2{\mu }_{\mathrm{eff}}E:E$ gives the Einstein viscosity formula for a dilute suspension:to first order in $\varphi$. The hydrodynamic interactions neglected here contribute beyond that dilute order.

For any vector $\mathbf{a}$, differentiating its components and the rotating basis gives the rotating-frame derivative formulaApply it twice to the position $\mathbf{x}$, with constant angular velocity $\mathbf{\omega }$. The inertial acceleration iswhere dots denote rotating-frame derivatives. Newton's second law therefore becomesThe last terms are the Coriolis force and the force corresponding to centrifugal acceleration.

For the particle on a uniformly rotating inclined plane, choose orthonormal vectors ${\mathbf{e}}_1$ along the horizontal $\xi$ axis and ${\mathbf{e}}_2$ up the slope, and set ${\mathbf{e}}_3={\mathbf{e}}_1\times {\mathbf{e}}_2$, the upward normal vector. Then ${\mathbf{e}}_z=\sin \theta {\mathbf{e}}_2+\cos \theta {\mathbf{e}}_3$ andHere $N$ is the normal reaction. The marble is modeled as the smooth sliding particle specified by the printed equations, without an additional rolling constraint. The rotating-frame Coriolis acceleration and centrifugal acceleration are respectivelyThe normal rotating-frame acceleration is zero, so the normal reaction is not generally $mg\cos \theta$; it satisfiesProjecting the same equation onto the two tangential axes gives exactly

Multiply the equations by $\dot{\xi }$ and $\dot{\eta }$ and add. The Coriolis force terms cancel, as expected because this force does no work relative to the rotating frame. The normal reaction also does no work along the plane. The conserved quantity is rotating-frame kinetic energy plus gravitational and centrifugal potential energy:This need not be the conserved inertial energy, since the externally driven rotating plane can exchange energy with the particle.

Choose the outward normal vector $n=(\cos \phi ,\sin \phi ,0)$ and the positive conventionThe ambient connection vanishes in Cartesian coordinates. Since ${\partial }_{\phi }n=e_{\phi }/\rho$ and ${\partial }_{\hat{z}}n=0$, the extrinsic curvature components are $K_{\phi \phi }=\rho$, $K_{\phi \hat{z}}=K_{\hat{z}\hat{z}}=0$. Taking the trace using the induced metric,Reversing the normal or using the negative extrinsic-curvature convention gives $-1/\rho$. The two principal curvatures are $1/\rho$ and zero in the chosen convention; the averaged mean curvature would be $K/2$, so it must not be confused with the requested trace.

The two displacement vectors $\mathbf{b}-\mathbf{a}$ and $\mathbf{c}-\mathbf{a}$ lie in the required plane. Because the three points are non-collinear, their cross productis nonzero and is a normal vector to that plane. A point with position vector $\mathbf{r}$ lies in the plane exactly whenExpanding the cross product givesThe terms $\mathbf{a}\cdot (\mathbf{a}\times \mathbf{b})$ and $\mathbf{a}\cdot (\mathbf{c}\times \mathbf{a})$ vanish, so the equation of a plane through three points becomesNon-collinearity is important: it ensures this is a genuine plane equation with nonzero normal, rather than the vacuous equality $0=0$.

For the specified coordinates,Their sum is $(3,3,6)$, with Euclidean norm $3\sqrt{6}$. Consequently either orientation of the unit normal is valid:As a check, the scalar triple product is $12$, so the plane equation simplifies to $x+y+2z=4$, which all three supplied points satisfy. The normal direction is $(1,1,2)$.

Write derivatives with respect to $u$. The arc-length condition is $f^{\prime 2}+g^{\prime 2}=1$. The tangent vectors areTheir scalar products give the first fundamental formChoose the unit normal vector $N=(-g^{\prime }\cos v,-g^{\prime }\sin v,f^{\prime })$, since ${\sigma }_u\times {\sigma }_v=fN$ and $f>0$. Scalar products of $N$ with the second partial derivatives giveThe mixed coefficient is zero. Reversing $N$ reverses the second fundamental form but not the Gaussian curvature. Its formula isDifferentiating the arc-length identity gives $f^{\prime }{f}^{\prime \prime }+g^{\prime }{g}^{\prime \prime }=0$. Substitute it directly, without dividing by $g^{\prime }$ (which may vanish), to obtain

Let $\Omega$ be a bounded piecewise smooth domain, $f=u{|}_{\partial \Omega }$ and $q={\partial }_{n}u$, with $n$ the outward normal vector. Integrating the divergence form from part (i) and using the divergence theorem gives two global relations for a linear boundary value problem:These conjugate global relations for the modified Helmholtz equation are identities between the prescribed Dirichlet boundary condition and its unknown normal derivative; neither is a pointwise boundary equation. Both hold for every nonzero boundary spectral parameter under the usual trace regularity assumptions.

If $u$ is real, its two boundary traces $f,q$ are real. Since ${\tilde{v}}_{\lambda }=\overline{v_{\overline{\lambda }}}$ and the outward normal vector is real,Thus apply complex conjugation to the first global relation for a linear boundary value problem, and replace its parameter by $\overline{\lambda }$, to get the second. For complex $u$, conjugation would also replace its traces by their conjugates, so that argument requires the stated reality assumption. There is also an algebraic redundancy in this particular parametrization: ${\tilde{v}}_{\lambda }=v_{1/\lambda }$, hence $\tilde{\mathcal{G}}(\lambda )=\mathcal{G}(1/\lambda )$ even for complex data. This does not change the requested conjugation identity.

Label the sides by $L,R,B,T$, and parametrize the two vertical sides by $s=y\in [-1,1]$, the horizontal sides by $s=x\in [-1,1]$. Their outward normal derivatives are $q_L=-u_x(-1,s)$, $q_R=u_x(1,s)$, $q_B=-u_y(s,-1)$ and $q_T=u_y(s,1)$. Let $f_j$ denote the prescribed side values. Assume compatible, sufficiently regular boundary traces; the square's corners have zero arclength measure and do not require separate normal values.

For sine collocation of square modified Helmholtz global relations, use Legendre polynomials for the known Dirichlet boundary condition and a Fourier sine series for the unknown normal derivative:The known coefficients are $d_{j\ell }=(2\ell +1){\int }_{-1}^1{f}_j{P}_{\ell }ds/2$. The sine functions have ${\int }_{-1}^1{\varphi }_m{\varphi }_{n}ds={\delta }_{mn}$ and form a complete basis in $L^2(-1,1)$. Choosing them for the normal derivative does not impose zero flux at a corner: the expansion is an $L^2$ representation, and endpoint values are not determined by it. If preserving corner values of the approximated Dirichlet trace is necessary, subtract its endpoint-interpolating line before polynomial approximation, and add that line back. All subsequent known-data integrals can alternatively be evaluated with the exact $f_j$.

Define the entire function basis transformsAt $\tau =\pm i{p}_m$ the quotient is evaluated by its removable limit or by the defining integral. For example, ${\mathsf{P}}_0=2\sinh \tau /\tau$, and polynomial expansion of $P_{\ell }$ expresses every ${\mathsf{P}}_{\ell }$ in derivatives of ${\mathsf{P}}_0$. PutFor compactness write $A=A(\lambda )$ and $B=B(\lambda )$. Substituting these expansions into the first global relation for a linear boundary value problem givesThe signs are fixed by the outward normal vectors, rather than by a choice of traversal direction. The second approximate global relation for a linear boundary value problem replaces $A$ by $-A$, leaving $B$ unchanged:Here $\simeq 0$ records omitted boundary-expansion tails. In a finite spectral method the selected equations are set equal to zero to solve for the $4N$ unknown real coefficients $c_{jm}$. For real data the two families obey the same complex conjugation relation as their exact counterparts.

Choose an origin near the bounded loop and write ${\mathbf{r}}^{\prime }$ for the source position. The magnetic vector potential and its far-field multipole expansion arewith the loop held fixed as $r=|\mathbf{r}|\to \infty$. The first term integrates to zero since $C$ is closed. Define its oriented vector area byThe second equality is Stokes theorem, applied componentwise to any oriented spanning surface. In particular the vector area is independent of the choice of spanning surface. From ${\oint }_{C}d(r_i^{\prime }{r}_j^{\prime })=0$, the matrix ${\oint }_C{r}_i^{\prime }d{r}_j^{\prime }$ is antisymmetric, which gives ${\oint }_C(\mathbf{r}\cdot {\mathbf{r}}^{\prime })d{\mathbf{r}}^{\prime }=\mathbf{S}\times \mathbf{r}$. HenceThe orientation follows the current by the right-hand rule. For $r\ne 0$, taking the curl of this leading magnetic vector potential gives the magnetic dipole fieldFor a planar loop, $\mathbf{S}$ is the signed area times the unit normal vector. For a nonplanar loop it is the oriented vector area, rather than the scalar area of an arbitrarily chosen surface.

The positive uniform ellipticity condition from part (i) is needed here. It is not explicitly imposed in the printed preamble. Without it the conclusion is false: on $(0,\pi )$, take $L=u^{\prime \prime }$, $\beta =0$, and $f=0$. The Neumann boundary condition is satisfied by $u(x)=\cos (mx)$, and $(L+m^2)u=0$ for arbitrarily large shifts.

There is also an actual error in the PDF's supplied trace inequality: a nonzero constant has nonzero trace and zero gradient. Use instead the valid multiplicative trace inequality on a bounded smooth domainFor completeness, choose a smooth vector field $X$ agreeing with the outward normal vector on the boundary. The divergence theorem gives, for smooth $u$,Approximation and the Sobolev trace theorem extend it to $H^1(U)$. Thus the intended existence result remains true after correcting the hint.

Set $q=\Vert u{\Vert }_2$, $d=\Vert Du{\Vert }_2$, ${\beta }_{-}=\max \{-\beta ,0\}$, $c_{-}=\max \{-c,0\}$, andUniform ellipticity and the Young inequality giveFor $\lambda \ge C_0+1$, the bilinear form is coercive on the real Hilbert space $H^1(U)$ with coercivity constant $\min \{\theta /2,1\}$.

All coefficients are bounded, and the Sobolev trace theorem gives $\Vert Tu{\Vert }_2\le C_T\Vert u{\Vert }_{H^1}$. The Cauchy-Schwarz inequality applied to each volume and boundary term therefore proves $|B_{\lambda }[u,v]|\le C_{\lambda }\Vert u{\Vert }_{H^1}\Vert v{\Vert }_{H^1}$. Also $\ell (v)={\int }_{U}fv$ is a continuous linear functional, with $|\ell (v)|\le \Vert f{\Vert }_2\Vert v{\Vert }_{H^1}$. Every hypothesis of the Lax-Milgram theorem is now verified, and it gives

For an oriented smooth surface $S$, the Gauss map sends $p$ to the chosen unit normal vector $N(p)\in S^2$. Since $|N{|}^2=1$, differentiation shows that $D{N}_p(X)$ is perpendicular to $N(p)$ and hence lies in $T_{p}S$.

In a local parametrization $\varphi (u,v)$ with $n=N\circ \varphi$, differentiatinggivesThus the bilinear form $(X,Y)\mapsto D{N}_p(X)\cdot Y$ is symmetric, so $D{N}_p$ is self-adjoint. The Gaussian curvature isWriting the coefficients of the first fundamental form asand those of the second fundamental form asone obtains

At an umbilic point, the self-adjoint map $D{N}_p$ has a repeated eigenvalue, so it is a scalar map. If every point is umbilic, there is a function $\lambda$ withEquality of mixed partial derivatives givesThe two tangent vectors are linearly independent, hence ${\lambda }_u={\lambda }_v=0$. Since ${\mathbb{R}}^2$ is connected, $\lambda$ is constant.

If $\lambda =0$, then $n$ is constant andso $S$ lies in a plane. If $\lambda \ne 0$, thenThus $n-\lambda \varphi =c$ is constant, andTherefore $S$ lies in a sphere of radius $1/|\lambda |$. This proves that the surface is part of a plane or part of a sphere.

A unit-speed curve on a smooth embedded surface is a geodesic exactly when its acceleration is normal to the surface, equivalently when its tangential acceleration vanishes.

Every unit-speed geodesic on the cylinder through $(1,0,0)$ has the formIndeed, unrolling the cylinder to its universal cover ${\mathbb{R}}^2$ turns these curves into straight lines. Directly,which is parallel to the cylinder's radial normal vector, verifying the geodesic characterization. Such a geodesic is closed exactly when $b=0$; its image is then the horizontal circle $z=0$.

Yes. In polar coordinates $(r,\theta )$ on ${\mathbb{R}}^2\setminus \{0\}$, use the Riemannian metricThe coordinate $u=\log r$ identifies this surface isometrically with the flat cylinder $\mathbb{R}\times S^1$. Through each point, the circle $u=\text{constant}$ is a closed geodesic, and every other geodesic has nonzero linear motion in $u$ and is not closed. Thus every point lies on a unique closed geodesic.

This is the plane through $a$ perpendicular to the nonzero normal vector $n$.

Take the parallel at height one,It has constant speed $\sqrt{2}$. Its distance from the axis is the constant $\rho =\sqrt{2}$, and its tangent is parallel to the parallel, so $\theta =0$ and it satisfies the Clairaut first integral for a surface of revolution.

However,whereas a normal vector to $x^2+y^2-z^2=1$ at $\gamma (t)$ isThese vectors are not parallel. By the ambient acceleration criterion for a surface geodesic, $\gamma$ is not a geodesic. Thus constant speed together with Clairaut's relation is not sufficient for the geodesic equations.

Write $q^1=u$, $q^2=v$, ${\sigma }_i={\partial }_i\sigma$, and let $g_{ij}={\sigma }_i\cdot {\sigma }_j$ be the first fundamental form. Sincedifferentiating ${\sigma }_u\cdot \dot{\gamma }$ and ${\sigma }_v\cdot \dot{\gamma }$, or equivalently taking the two inner products of the second formula with ${\sigma }_k$, givesThe matrix $(g_{k\ell })$ is invertible because $\sigma$ is a regular embedded surface parametrization. Consequently the two geodesic equations hold exactly when $\ddot{\gamma }$ is orthogonal to both tangent vectors ${\sigma }_u$ and ${\sigma }_v$, which says precisely that $\ddot{\gamma }$ is a normal vector. This is the ambient acceleration criterion for a surface geodesic.