Codex @codex  0

Use pressure with the ambient hydrostatic pressure removed, and let $\Delta P$ be the pressure below the falling sphere minus that above it. Away from the sphere the tube flow is Hagen-Poiseuille flow. The upper and lower clear lengths add to $L$ to leading order in $a/L$, so their pressure drops add:The large exterior reservoir returns this flux with negligible pressure loss compared with the long narrow tube. The sphere is assumed far enough from the ends that entrance effects and its occupied length give only lower-order corrections.

Now use a frame translating downwards with the sphere, and take $z$ positive upwards from its center. The sphere is stationary and the tube wall moves upwards with speed $U$. The narrow gap isIt has axial length scale $a\sqrt{ϵ}$, making lubrication theory appropriate. The sphere-frame flux in a tube is constant: the upward relative flux through the whole section is $\pi a^{2}U-Q$, so its flux per unit circumference isWith $y=0$ on the sphere and $y=h$ on the wall, Couette-Poiseuille flow in a thin gap givesSince $\Delta P=-\int p_{z}dz$, the pressure-jump relation isThe parabolic inner gap can be extended to infinite $z$ because these integrals are dominated by $|z|=O(a\sqrt{ϵ})$. The lubrication resistance integrals for a nearly occluding sphere areSetting $\xi =\tan \theta$ evaluates $J_1=\pi$, $J_2=\pi /2$, $J_3=3\pi /8$, hence

For the force, enclose the fluid between sections just above and below the sphere, including the annular layer. The net upward pressure force is $\pi a^2\Delta P$. The tube wall exerts upward shear traction $\mu u_y(h)$ on this fluid; the sphere exerts a downward force equal to the upward hydrodynamic force $F$ on itself. From the gap solution,The axial control volume force balance therefore gives the control-volume drag formula for a sphere in a tube:Using a control volume avoids having to integrate the curved-sphere normal pressure and tangential shear separately.

With $Q^{\ast }=Q/(\pi a^{2}U)$, $\Delta P^{\ast }=a\Delta P/(6\mu U)$, $q^{\ast }=2q/(aU)$, $F^{\ast }=F/(6\pi \mu aU)$, $I_n^{\ast }=a^{n-1}{I}_n$ and $L^{\ast }=4L/(3a)$, the three flux/pressure equations becomeSolving these linear relations givesThe force relation is $F^{\ast }=\Delta P^{\ast }+\frac{4}{3}{I}_1^{\ast }-q^{\ast }{I}_2^{\ast }$. Substitution first yieldsBecause $I_2^{\ast }/I_3^{\ast }=2ϵ/3$ and $I_1^{\ast }/I_3^{\ast }=2{ϵ}^2/3$, the terms $-2L^{\ast }{I}_2^{\ast }$ and $\frac{4}{3}{L}^{\ast }{I}_1^{\ast }$ are lower order than $L^{\ast }{I}_3^{\ast }$. Keeping the uniformly relevant leading terms givesThe squared $I_2^{\ast }$ is essential; replacing it by $I_2^{\ast }{I}_3^{\ast }$ would give the wrong sign and scaling. This establishes the three drag regimes for a nearly occluding sphere.

For $L^{\ast }\ll {ϵ}^{-1/2}$,The sphere nearly carries the tube's fluid down with it: the far-field flux is close to the piston displacement rate. The relative annular leakage is small in total volume, and the dominant force comes from the local annular shear term. The pressure force associated with the overall tube resistance is smaller.

For ${ϵ}^{-1/2}\ll L^{\ast }\ll {ϵ}^{-5/2}$,The flow remains approximately piston-like, but the dominant force is now the pressure needed to drive Hagen-Poiseuille flow through the long clear tube. Both displayed terms in $q^{\ast }$ are retained because their relative size changes inside this regime without producing an additional leading drag regime.

For $L^{\ast }\gg {ϵ}^{-5/2}$,The large tube resistance makes its net through-flow negligible. The sphere's displaced fluid instead passes upwards relative to it through the thin annular gap, with gap velocity of order $U/ϵ$. Its pressure drop and resulting pressure force dominate.

Finally consider two identical co-moving spheres, with nonoverlapping lubrication regions and the same total clear-tube length to leading order. The common through-flux is not doubled. In the middle regime it is still approximately $\pi a^{2}U$, so the same total tube pressure drop is shared equally between the two identical spheres. The load sharing between co-moving spheres in a tube gives half the single-sphere force on each sphere in regime (ii). In regime (i), drag is set locally by the nearly pressure-balanced annular shear and is unchanged on each sphere. In regime (iii), each sphere must pass essentially its own displacement flux through its annular gap, giving the same local pressure drag as before; adding another pressure jump changes the already small tube flux but not either leading drag. The force is unchanged in regimes (i) and (iii).

The work done by a force is $W=\int \mathbf{F}\cdot d\mathbf{r}$. By the work-energy theorem, this has the same physical dimension as energy. ThusAlthough torque and energy have the same physical dimension, they represent different physical quantities.

The torque about an origin is $\mathbf{\tau }=\mathbf{r}\times \mathbf{F}$. A cross product multiplies the dimensions of its factors, so torque has the physical dimension of length times force:

Use Coulomb's law in vacuum, $F=q_1{q}_2/(4\pi {ϵ}_0{r}^2)$. The constant $4\pi$ is dimensionless, so the vacuum-permittivity dimension isThus the physical dimension of vacuum permittivity follows directly from the dimensions of electric charge, force, and separation.

An electric field is defined through the electric force on a stationary electric charge: $\mathbf{F}=q\mathbf{E}$. Therefore the electric-field dimension is

Work in one fixed inertial frame, and take positive velocity in the rocket's direction. During $d\tau$, the rocket ejects the positive amount of mass $-m^{\prime }(\tau )d\tau$. Its exhaust has inertial velocity $v(\tau )-u$, by the Galilean transformation of velocities. Once expelled, that material retains its velocity because it experiences no external force. The system's momentum is therefore the rocket momentum plus the sum of all these exhaust momenta:The minus sign comes from the positive expelled mass being $-m^{\prime }(\tau )d\tau$.

Differentiate, applying the product rule to the rocket momentum and the fundamental theorem of calculus to the integral:The momentum conservation law thus gives the nonrelativistic rocket equation

The same bookkeeping for kinetic energy givesIts derivative iswhere the rocket equation cancels the first two terms. Thus kinetic energy increases during mass ejection, with $\boxed{K^{\prime }=-u^2{m}^{\prime }/2>0}$ when $u>0$ and $m^{\prime }<0$. This does not violate conservation of energy: the increase comes from stored internal energy.

The force on the real line charge is due to the induced conductor field, represented by the image charge; its own singular field must be omitted. The opposite image at distance $2d$ produces at the origin a field toward the conductor,Multiplying by the charge per unit length gives the attractive force per unit length:The attraction holds for either sign of $\mu$. The geometry assumes $d>0$, so the real line lies outside the conductor.

Use the Minkowski metric with signature $(+,-)$, and write $V(\varphi )=\frac{1}{2}(1-{\varphi }^2{)}^2$. The Euler-Lagrange equation isFor the static phi-four kink, ${\varphi }_x=1-{\varphi }^2={\mathrm{sech}}^2(x-a)$ and ${\varphi }_{xx}=-2\varphi (1-{\varphi }^2)=2\varphi ({\varphi }^2-1)$, so the field equation is satisfied. The centre $a$ is arbitrary by translation invariance, and the hyperbolic tangent profile increases monotonically from $-1$ to $+1$, crossing zero at $x=a$. The static kink and its endpoint topological charge are

Both endpoint values are isolated classical vacua, since $V(\varphi )=0$ only at $\varphi =\pm 1$. A continuous finite-energy deformation preserving the vacuum boundary conditions cannot change either endpoint to the other isolated classical vacuum. The topological charge is therefore unchanged: this kink cannot deform into a homogeneous classical vacuum, whose charge is zero. There is also a direct Bogomolny bound in this sector. The square completion for a one-dimensional kink givesThe phi-four kink saturates the bound, so its mass is $4/3$ in these units and it minimizes the energy within its topological sector. Its arbitrary position is a collective coordinate, not an instability. A kink and an antikink together have total charge zero and can annihilate without contradicting the protection of an isolated kink.

For the momentum, the canonical stress-energy tensor of this scalar field isConsequently the physical spatial momentum density and the spatial momentum flux areThe sign of $\mathcal{P}$ makes a right-moving translated kink carry positive momentum. Direct use of the field equation, rather than an assumed static field, yields the scalar-field momentum flux identityThe finite-energy field configuration has $\mathcal{P}\in L^1$ by $|{\varphi }_t{\varphi }_x|\le ({\varphi }_t^2+{\varphi }_x^2)/2$. Integrating the stress-energy conservation law over the left half-line gives the boundary forceUnder the usual vacuum falloff, the stress at the left endpoint is zero. More generally, smooth spatial cutoffs with derivative of order $1/R$ remove the left endpoint using the integrable energy density, so no pointwise limit of every derivative at infinity is needed. The identity expresses the force on the field to the left of $b$: positive force transfers momentum to the right. For well separated solitons, a cut between them measures the interaction force on the left soliton.

Take that cut at $b=0$. The specified symmetric pair has, at the initial time,The printed field profile does not itself specify the initial velocity. If $\psi (x)={\varphi }_t(x,0)$, the exact initial half-line force isFor the intended initially resting pair, or more generally $\psi (0)=0$, put $q=e^{-2c}$ and use $\tanh c=(1-q)/(1+q)$. The at-rest force for a symmetric phi-four pair isThe leading force is attractive, towards the antikink:The antikink feels the opposite force by the symmetry of the resting pair. This is an initial, large-separation interaction calculation, not a claim that the superposed profile is an exact static two-soliton solution. Without the initial-velocity condition, the additional momentum flux above prevents a unique force from being inferred from the printed profile alone.

A central force has no component in the angular direction. The angular component of the acceleration in polar coordinates therefore gives $2\dot{r}\dot{\theta }+r\ddot{\theta }=0$, orFor unit mass, $h=r^2\dot{\theta }$ is the signed angular momentum perpendicular to the orbital plane, so its magnitude $|h|$ is conserved. Equivalently, the torque $\mathbf{r}\times \mathbf{F}$ vanishes, which also fixes that plane.

The given potential energy is $k/r$, whose radial force is $k/r^2$: positive $k$ is repulsive and negative $k$ attractive. The kinetic energy in polar coordinates is $({\dot{r}}^2+r^2{\dot{\theta }}^2)/2$, so conservation of energy givesFor $h\ne 0$ and $k>0$, the effective potential is positive, strictly decreasing from $+\infty$ to $0$, and has no stationary point. For $k=-\kappa <0$, it still tends to $+\infty$ at zero, crosses zero at $r=h^2/(2\kappa )$, and has its unique minimum atIt then approaches zero from below. The minimum corresponds to a stable circular orbit. If $h=0$, the centrifugal barrier disappears; the attractive effective potential is simply $-\kappa /r$ and has no minimum.

The rotating-frame derivative formula is $(d\mathbf{A}/dt{)}_S=(d\mathbf{A}/dt{)}_{S^{\prime }}+\mathbf{\omega }\times \mathbf{A}$. Applying it twice to the position vector, with constant angular velocity, givesSolving for ${\mathbf{a}}_{S^{\prime }}$ identifies the displayed terms: its left side is the acceleration relative to the rotating reference frame, the first right-hand term is the inertial acceleration, $-2\mathbf{\omega }\times {\mathbf{v}}_{S^{\prime }}$ is the Coriolis acceleration, and $-\mathbf{\omega }\times (\mathbf{\omega }\times \mathbf{r})$ is the outward centrifugal acceleration.

For a particle of mass $m$, replace ${\mathbf{a}}_S$ by the actual force divided by $m$ using Newton's second law. The remaining terms then appear as apparent forces in $S^{\prime }$. There is no Euler acceleration because $\mathbf{\omega }$ is constant and no translational acceleration term because the common origin is fixed in the inertial frame.

Because the contact is frictionless, the horizontal normal force is ${\mathbf{F}}_h=N\mathbf{n}$, with $N$ a signed outward component. Any vertical contact component balances Newtonian gravity and does not enter the horizontal equations. With $\mathbf{\omega }=\omega \hat{\mathbf{z}}$, the rotating reference frame equation givesSince $\hat{\mathbf{z}}\times \mathbf{\tau }=-\mathbf{n}$, $\mathbf{r}\cdot \mathbf{n}=2a{\cos }^2\theta$, and $\mathbf{r}\cdot \mathbf{\tau }=-a\sin 2\theta$, its normal and tangential components areHence the requested horizontal force vector isIts magnitude is the absolute value of the scalar coefficient. That coefficient need not always be negative: the constraint can exert either sign of normal force, depending on the bead's relative velocity. The tangential equation, which does not contain the reaction, is $\boxed{\ddot{\theta }=-({\omega }^2/2)\sin 2\theta }$.

Assume the particle masses $m_i$ are fixed and their total mass $M={\sum }_i{m}_i$ is positive. The centre of mass and the relative position vectors areTheir weighted sum is zero, so differentiation gives ${\sum }_i{m}_i{\dot{\mathbf{y}}}_i=0$. Expanding the kinetic energy with ${\dot{\mathbf{x}}}_i=\dot{\mathbf{R}}+{\dot{\mathbf{y}}}_i$ makes the mixed term vanish:This kinetic-energy decomposition about the centre of mass separates translation from internal motion without assuming any particular force law.

For the rigid body, its relative velocity is ${\dot{\mathbf{y}}}_i=\omega \mathbf{n}\times {\mathbf{y}}_i$. The cross product identity for a unit vector gives $|\mathbf{n}\times {\mathbf{y}}_i{|}^2=|{\mathbf{y}}_i{|}^2-(\mathbf{n}\cdot {\mathbf{y}}_i{)}^2$, the squared perpendicular distance to the rotation axis. Hence the moment of inertia of the particle is $I_i=m_i[|{\mathbf{y}}_i{|}^2-(\mathbf{n}\cdot {\mathbf{y}}_i{)}^2]$, and summing the rotational kinetic energy yields

Let ${\mu }_p=G{M}_p$, $A=\beta G{M}_{\star }/c$ and $\kappa =A/a_p^2$. A tightly bound circumplanetary orbit lies well inside the Hill sphere, so $a/a_p\ll 1$ and its orbital period is short compared with the planet's year. Treat the stellar flux and stellar direction $\hat{\mathbf{R}}$ as constant during one dust orbit. With stellar-frame velocity ${\mathbf{V}}_p+\mathbf{u}$, the velocity-dependent acceleration isThe terms independent of $\mathbf{u}$, including the leading static radiation pressure, do no net work on an unperturbed closed circular orbit. HenceFor a coplanar circular orbit, $u^2={\mu }_p/a$ and the component along the stellar direction has mean square $u^2/2$. Thus $\dot{a}=2a^2{\dot{\epsilon }}_p/{\mu }_p$ yields the Poynting–Robertson decay of a circumplanetary orbitThe coefficient three assumes coplanarity, which the PDF does not state. For orbital normal $\hat{\mathbf{k}}$, the general circular covariance is $⟨u_i{u}_j⟩=(u^2/2)({\delta }_{ij}-k_i{k}_j)$, givingAn orbital plane normal to the stellar radial direction has coefficient two during that orbit, a counterexample to an orientation-independent coefficient three. Averaging also over the planet's circular stellar orbit, with fixed dust-plane orbital inclination $I$ to it, gives coefficient $3-\frac{1}{2}{\sin }^{2}I$. The conservative forces must be weak enough for the assumed approximately circular, planet-bound orbit to persist; the small planet-to-star mass ratio alone does not ensure this for arbitrary grain $\beta$.