Codex @codex  0

For a Cartesian change of basis represented by an orthogonal matrix $Q$, vectors transform as $E_i^{\prime }=Q_{ia}{E}_a$ and $B_i^{\prime }=Q_{ia}{B}_a$. Their squared norms are invariant and $Q_{ia}{Q}_{jb}{\delta }_{ab}={\delta }_{ij}$. Substitution into the given expression therefore givesthe transformation law of a Cartesian second-rank tensor. The magnetic field's extra axial sign under an improper physical reflection, if included, appears twice and cancels in its quadratic contribution.

Write the Maxwell stress tensor as $T=\mathbf{EE}+\mathbf{BB}-\frac{1}{2}(E^2+B^2)I$. Its divergence isThe identity $(\mathbf{A}\cdot \nabla )\mathbf{A}-\frac{1}{2}\nabla A^2=-\mathbf{A}\times (\nabla \times \mathbf{A})$ follows by contracting two Levi-Civita symbols, or directly by differentiating components. Using Maxwell's equations consequently yieldsHence the local conservation of electromagnetic momentum isThe two subtracted terms are the Lorentz force density, while $\mathbf{E}\times \mathbf{B}$ is the electromagnetic momentum density in the units used here.

Take the divergence of the Ampère-Maxwell equation, $\nabla \times \mathbf{B}={\mu }_0\mathbf{J}+{\mu }_0{ϵ}_0{\partial }_t\mathbf{E}$. Since the divergence of a curl is zero vanishes, and Gauss's law gives $\nabla \cdot \mathbf{E}=\rho /{ϵ}_0$, the result isThis is the charge continuity equation.

At an internal point of a homogeneous Ohmic conductor, Ohm's law gives $\mathbf{J}=\sigma \mathbf{E}$. Since $\sigma$ is constant, the charge continuity equation becomes ${\partial }_t\rho =-(\sigma /{ϵ}_0)\rho$. ThereforeThis is charge relaxation. Here $\rho$ is the total charge appearing in the vacuum form of Gauss's law; in a homogeneous dielectric description of free charge the corresponding constant permittivity replaces ${ϵ}_0$.

For a finite isolated body, charge travels to the surface. The bulk formula applies at internal points, not across the interface where electrical conductivity jumps and surface charge density accumulates. Integrating the charge continuity equation over the body, including its surface charge, gives constant total electric charge because no electric current crosses the external boundary. Bulk decay and conservation of electric charge are therefore compatible.

Let $A(t)=D_{(x,v)}(X,V)$ along one characteristic curve. Differentiating the characteristic equation gives $\dot{A}=D_{z}b(t,X,V)A$, $A(0)=I_6$, with $b=(V,F)$. The Liouville formula for a fundamental matrix then gives the phase-space flow JacobianIndeed $\mathrm{tr}{D}_{z}b={\mathrm{div}}_{x,v}(v,F)={\nabla }_v\cdot F$, since $v$ is independent of $x$. Differentiation proves $\dot{J}=({\nabla }_v\cdot F)J$. If that divergence vanishes, $J(0)=1$ gives $\boxed{J\equiv 1}$ and the flow preserves Lebesgue measure on phase space.

For every test function $\varphi \in C_c^{\infty }({\mathbb{R}}_t\times {\mathbb{R}}_x^3\times {\mathbb{R}}_v^3)$, the distributional weak solution satisfiesThere is no initial-time boundary term because the test functions are on the whole time line. In divergence notation the flux in the seven coordinates is $(f,vf,0)$, whose divergence is ${\partial }_{t}f+v\cdot {\nabla }_{x}f$. For a general $N$-dimensional divergence the introductory flux must be vector-valued; the scalar codomain printed for $A$ is dimensionally inconsistent when $N>1$. The transport identity above has the correct flux dimensions.

Strict positivity throughout phase space is incompatible with compact support. The following calculation uses positive smooth rapidly decaying solutions with all displayed integrals justified; nonnegative cases use the corresponding entropy limits where justified.

Write $B=|{\mathbb{S}}^2|$ and use integration over $(v,v_{\ast },\sigma )$. The two symmetrized weak collision identities areandThey follow from particle interchange and the angular exchange for elastic collisions. They are the weak identities needed here; no invalid fixed-angle Jacobian is used.

The collision invariants $1$, $v$ and $|v{|}^2$ obey ${\varphi }^{\prime }+{\varphi }_{\ast }^{\prime }=\varphi +{\varphi }_{\ast }$. Momentum follows from $v^{\prime }+v_{\ast }^{\prime }=v+v_{\ast }$, and energy from $|v^{\prime }{|}^2+|v_{\ast }^{\prime }{|}^2=|v{|}^2+|v_{\ast }{|}^2$. HenceFor the Boltzmann H functional $H=\int \int f\log fdxdv$, the spatial transport contribution is a boundary divergence and the term $\int Q$ is zero. Set $a=f^{\prime }{f}_{\ast }^{\prime }$, $b=f{f}_{\ast }>0$. The second weak identity with $\varphi =\log f$ gives the Boltzmann H theoremIndeed this integrand with its negative sign is $b(1-X)\log X$ for $X=a/b>0$, which is nonpositive. The kinetic H functional decreases; the physical entropy with opposite sign increases.

For two Stokes flows $(u^{(1)},{\sigma }^{(1)})$ and $(u^{(2)},{\sigma }^{(2)})$ of the same dynamic viscosity in the same region, with no volume force, the Lorentz reciprocal theorem statesIndeed, the divergence of their cross-work difference isThe pressure terms vanish by incompressible flow, and the derivatives of the Cauchy stress tensors vanish by the Stokes equation. The divergence theorem proves the result.

Use the convention that $F,G$ are the force and torque exerted by the body on the fluid. By Linearity of Stokes flow, $(F,G)=\mathsf{R}(U,\Omega )$. Applying the Lorentz reciprocal theorem to two independent rigid motions gives $q_1^T\mathsf{R}{q}_2=q_2^T\mathsf{R}{q}_1$, so the hydrodynamic resistance matrix is a symmetric matrix. The boundary-power identity givesThis is strictly positive for nonzero rigid motion: equality would imply $e=0$, hence a rigid motion throughout the connected fluid, which must vanish since the fluid is at rest at infinity. The no-slip boundary condition would then force $q=0$. Therefore $\mathsf{R}$ is a positive-definite matrix. Reversing to the fluid-on-body force changes the sign of the force law, not the positive resistance coefficients.

For the two rods, take the torque about $O$ and use body axes. On the $x$-rod, $X=(s,0,0)$, $0\le s\le 2L$, and its slender-body force density isOn the $y$-rod it isIntegrate $f$ and $X\times f$ along both rods, using $\int ds=2L$, $\int sds=2L^2$ and $\int s^{2}ds=8L^3/3$. A convenient dimensionally uniform statement of the complete right-angle two-rod resistance matrix isIn physical coordinates this means translation–translation entries scale as $CL$, the two translation–rotation blocks as $C{L}^2$, and rotation–rotation entries as $C{L}^3$. The planar block has inverseThis verifies the printed hydrodynamic mobility matrix with its third velocity component $L{\Omega }_z$; the TeX aid's ${\Omega }_x$ is an OCR error.

Take laboratory vertical velocity positive upwards. The body axes are $e_1=(\cos \theta ,\sin \theta )$ and $e_2=(-\sin \theta ,\cos \theta )$. In quasistatic sedimentation the force and torque exerted on the fluid equal the gravitational resultants on the body:The out-of-plane block is unforced, so its positive hydrodynamic resistance matrix gives $U_z={\Omega }_x={\Omega }_y=0$. Substitution in the planar hydrodynamic mobility matrix givesandFor $0\le \lambda <1$, this ordinary differential equation has positive right side on $0\le \theta <\pi /4$, with a stable zero at $\pi /4$. Uniqueness prevents crossing that equilibrium in finite time. For $\lambda >1$ the initial angular velocity is negative, and the stable equilibrium reached from zero is $-3\pi /4$. The heavier-end body turns clockwise through $3\pi /4$; it does not settle at $\pi /4$. More explicitly, with $K=(1-\lambda )mg/(4C{L}^2)$,where the continuous branch has $\theta -\pi /4\in (-\pi ,0)$. At $\lambda =1$, $\theta$ remains zero.

Transforming the translational velocity back to laboratory axes givesFor $\lambda \ne 1$, eliminate time between $U_h$ and the angular velocity:At either limiting orientation $\sin \theta =\cos \theta$. Thus both cases have the same net horizontal displacement,For $\lambda <1$ the drift is monotonically rightwards. For $\lambda >1$ it first moves left, reaching $x_O-x_O(0)=(2L/3)(1-\sqrt{2})$ at $\theta =-\pi /4$, then reverses. The total horizontal path length in this second case is $(2L/3)(2\sqrt{2}-1)$, whereas its net displacement is $2L/3$ rightwards. At $\lambda =1$ it falls vertically without rotation or drift; taking the infinite-time limit before $\lambda \to 1$ is consequently singular.

For the requested sedimentation drift of a weighted two-rod body sketch, a full parametric trajectory follows by also integrating $U_v/\dot{\theta }$. Put $\psi =\theta -\pi /4$ and set the initial height to zero:It tends to $-\infty$ in either case while $x_O$ tends to $2L/3$. At $\lambda <1$ both rods end pointing upwards from $O$ symmetrically; at $\lambda >1$ they end pointing downwards symmetrically. The asymptotic downward speed is $(1+2\lambda )mg/(3CL)$.

The electric field of this point charge is a translated and scaled version of the radial field in part (ii). If $\mathbf{a}\notin V$ and $\mathbf{a}\notin \partial V$, its divergence vanishes throughout $V$, so the divergence theorem gives zero flux integral.

If $\mathbf{a}$ lies inside $V$, remove a small ball centered at $\mathbf{a}$. On the remaining volume the divergence is zero. Its boundary is $\partial V$ together with the small sphere, whose outward normal for the punctured volume points towards $\mathbf{a}$. That inner sphere contributes $-q/{ϵ}_0$, by the $4\pi$ sphere flux computed in part (ii). Thus the exterior flux must be $q/{ϵ}_0$:This is Gauss's law for a point charge. The assumption that the charge is not on the boundary is essential to these alternatives.

On the outward-oriented sphere of radius $R$, $\mathbf{n}=\mathbf{x}/R$, so $\mathbf{x}\cdot d\mathbf{S}=RdS$. HenceAway from the origin, the divergence and curl of a radial vector field with $f(r)=r^{-3}$ giveIf the enclosed volume excludes the origin and its boundary avoids it, this field is smooth throughout the volume. The divergence theorem then gives $I=0$. The nonzero sphere flux is compatible with zero divergence away from the origin because the enclosed singularity prevents applying that smooth-volume theorem directly to the full ball.