Past exam of the physics B course of the University of Cambridge / 2015

For a smooth isolated conducting sphere, Gauss's law gives the surface electric field $E_s=Q/(4\pi {ϵ}_0{a}^2)$ and electric potential $V=Q/(4\pi {ϵ}_{0}a)=a{E}_s$. At the specified dielectric strength, $V_{\max }=3.0\times 10^3\mathrm{V}$. The capacitance of an isolated conducting sphere is $C=4\pi {ϵ}_{0}a$, so its capacitor energy isThis is the ideal breakdown-limited energy of a conducting sphere, taking air's permittivity approximately equal to vacuum permittivity and neglecting nearby objects and surface irregularities. The vacuum permittivity used numerically is $8.8542\times 10^{-12}\mathrm{F}{\mathrm{m}}^{-1}$.^[ref]

Take the initial thrower and receiver positions to be $(-a,0)$ and $(a,0)$, with counterclockwise angular velocity $\omega$. Interpret $v$ as the launch speed relative to the thrower. The ball inherits the thrower's tangential velocity, so in an inertial frame its initial horizontal velocity is $(v,-a\omega )$. The flight time across the diameter is $t_f\simeq 2a/v$. During that time the receiver moves approximately $+a\omega t_f$ sideways while the ball moves $-a\omega t_f$. Thus the transverse miss on a rotating platform isThe ball passes behind the receiver in the sense of rotation. Equivalently, in the rotating frame its transverse Coriolis acceleration has magnitude $2\omega v$ to leading order; integrating over $t_f$ gives $d=\frac{1}{2}(2\omega v)t_f^2$. Terms neglected here are higher order in $a\omega /v$. Specifying an initially radial launch velocity in the inertial frame instead would define a different throw and would lose the inherited tangential velocity.

The dynamic viscosity measures the transport of tangential momentum between adjacent layers. The kinetic estimates of gas transport coefficients givewhere $n$ is number density, $m$ the molecular mass, $\lambda$ the mean free path, and ${\sigma }_c$ a collision cross section. At fixed temperature, raising the pressure increases the number of momentum carriers but shortens their mean free path by the same factor. Consequently the dynamic viscosity of a dilute gas is approximately independent of pressure.

For a temperature-independent collision cross section, the same estimate gives $\eta \propto \sqrt{m{k}_{B}T}/{\sigma }_c$, so $\eta \propto T^{1/2}$ in the hard-sphere model. Faster molecules carry momentum across layers more rapidly. Real intermolecular interactions make the effective collision cross section temperature-dependent, changing this simple exponent. The pressure-independence also fails when wall collisions, rather than molecule–molecule collisions, set the transport distance; the continuum estimate assumes a small Knudsen number.

The varying magnetic flux through the ring induces an electromotive force by Faraday's law of induction. Its direction follows Lenz's law. With solenoid current $I_s$, ring current $I_r$, ring self-inductance $L_r$, electrical resistance $R$, and mutual inductance $M(z)$,The ring sits in a nonuniform magnetic field: coupling decreases above the solenoid, so $M>0$ and $M^{\prime }<0$ with a common orientation for the two windings. For $I_s=I_0\cos \omega t$, the in-phase part of the induced ring current is negative:Thus alternating-current induction levitation can balance the ring's weight when $⟨F_z⟩=mg$. A falling mean lifting force with height supplies vertical stability. In the strongly inductive limit the two currents are almost opposite: the ring repels the solenoid during both half cycles. Finite self-inductance matters; an ideal purely resistive loop would carry a current exactly in quadrature with $I_s$ and give zero average force in this model. The source supplies the Joule heating in the real ring, and any mechanical work while it rises.

After the photon has thermalized, its energy $E_{\gamma }=hc/\lambda$ belongs to the macroscopic body's internal energy. The body's heat capacity makes the resulting temperature change negligible. Its equilibrium entropy therefore increases byThe Boltzmann entropy $S=k_B\ln \Omega$ then gives the microstate multiplier from photon absorption:For scale, explicitly choosing a representative visible wavelength $\lambda =500\mathrm{nm}$ and room temperature $T=300\mathrm{K}$ gives $\Delta S=1.32\times 10^{-21}\mathrm{J}{\mathrm{K}}^{-1}$, or $\Delta S/k_B=95.9$, and a multiplier approximately $4.5\times 10^{41}$. These are illustrative choices, not additional supplied data: the exact multiplier depends exponentially on the particular visible wavelength and temperature. The exact Planck constant, speed of light, and Boltzmann constant fix the conversion.^[ref]

In electrostatics, the Maxwell equations in matter areHere ${\rho }_f$ is free charge density and $\mathbf{P}$ is polarization density. For a homogeneous isotropic linear dielectric, the constitutive equation is $\mathbf{D}={ϵ}_0{ϵ}_r\mathbf{E}$. Integrating the first equation over a thin pillbox straddling an interface, with unit normal $\mathbf{n}$ from medium 1 to medium 2, gives $\mathbf{n}\cdot ({\mathbf{D}}_2-{\mathbf{D}}_1)={\sigma }_f$. Integrating the second around a vanishingly narrow rectangular loop gives $\mathbf{n}\times ({\mathbf{E}}_2-{\mathbf{E}}_1)=0$. These are the electrostatic interface boundary conditions. When there is no free surface electric charge, the normal electric displacement field and tangential electric field are continuous; the normal electric field generally is not.

Neglect end fringing field of a capacitor. In either axial segment of the coaxial cylindrical capacitor, a cylindrical Gaussian surface gives $D_r={\lambda }_f/(2\pi r)$, where ${\lambda }_f$ is free electric charge per length. Integrating the radial electric field between the conductors givesThe liquid and air segments share the same conductor electric potential difference, so their capacitances add in parallel. Thus the partially filled coaxial capacitor hasHere the question's $ϵ$ is relative permittivity, and air's relative permittivity is approximated by one.

Write $x$ for the rise above the external level. The wide reservoir keeps that level approximately fixed. At constant applied voltage, one must include work done by the source: the appropriate electrical contribution to the mechanical potential is $-C{V}^2/2$. The upward dielectric insertion force at fixed voltage is consequentlyThe excess hydrostatic pressure $\rho gx$ acts over the annular area $A=\pi (b^2-a^2)$. Alternatively its gravitational potential energy is $\rho gA{x}^2/2$, since liquid initially at the reservoir level must be raised through all heights from zero to $x$. Balancing the derivatives givesThis assumes $h+x<L$, a communicating liquid reservoir, negligible capillarity, and a sufficiently long capacitor that end effects are small. Differentiating gives the dielectric-rise voltmeter's sensitivity at a bias $V$:To obtain its numerical value in $\mu \mathrm{m}{\mathrm{V}}^{-1}$, multiply the right-hand side in SI units by $10^6$. A convenient parametric scale, with $\beta =b/a$, isFor a concrete estimate, choose $a=1\mathrm{mm}$, $b=2\mathrm{mm}$, mass density $\rho =10^3\mathrm{kg}{\mathrm{m}}^{-3}$, relative permittivity $ϵ=2$ and operating voltage $V=100\mathrm{V}$. These are illustrative choices, not supplied apparatus data. They giveThus this millimetre-scale example requires fine optical displacement measurement for volt-level changes. The parametric formula states exactly how the estimate changes with the chosen geometry, liquid and bias.

The displacement is quadratic in voltage, cannot distinguish its sign, and has zero first-order sensitivity at zero bias. Smaller gaps and larger relative permittivity increase sensitivity, but capillarity, leakage, dielectric breakdown, temperature-dependent liquid properties, and finite available rise then limit performance. A nonconducting liquid avoids a sustained conduction load, although charging still draws electric charge and leakage is possible. A high voltage can give a readable displacement; a low-voltage instrument would require fine displacement measurement and careful calibration.

The differential form of Faraday's law of induction isBy Stokes theorem, for a stationary oriented circuit it becomes $\mathcal{E}=\oint \mathbf{E}\cdot d\mathbf{\ell }=-d{\Phi }_B/dt$, where ${\Phi }_B=\int \mathbf{B}\cdot d\mathbf{S}$ is magnetic flux and $\mathcal{E}$ is electromotive force. Circuit orientation and surface normal obey the right-hand rule. Lenz's law is the physical meaning of the minus sign: the induced current's magnetic field opposes the change in linked magnetic flux. For a moving circuit, the force per unit electric charge is $\mathbf{E}+\mathbf{v}\times \mathbf{B}$, and its line integral is the corresponding electromotive force. Motion through a static magnetic field can therefore induce a voltage even when ${\partial }_t\mathbf{B}=0$.

Let $\theta$ be the wire's angle from the downward vertical, $r=L\sin \theta$, and choose $\mathbf{\omega }$ and $\mathbf{B}$ upward, as in the figure. The motional term along the wire has radial component $B\omega s$ at cylindrical radius $s$. In steady electrical equilibrium, integrating from the grounded axis to the bob givesThe Lorentz force on this positive electric charge is radially outward. Its ratio to the required centripetal term is $QB/(m\omega )\lesssim C{B}^2{L}^2/(2m)$, so the stated heavy-bob condition makes its effect on the zero-field conical pendulum geometry negligible. With wire tension $T_w$, $T_w\cos \theta =mg$ and $T_w\sin \theta =m{\omega }^{2}r$, giving $\cos \theta =g/({\omega }^{2}L)$. Hence the magnetically charged conical pendulum hasA nonzero-radius orbit requires ${\omega }^{2}L>g$. Reversing the sense of rotation or the magnetic field reverses the sign of $Q_0$.

The response to the added electric field requires an electrical convention. In the usual mechanical approximation, treat the pre-existing electric charge as $Q_0$, include its upward electric force $Q_{0}E$, and neglect all magnetic-force feedback. ThenThis predicts an increased radius for $Q_{0}E>0$. It is the result when the newly applied field is used only as a probe force on the original electric charge.

There is a first-order qualification if the bob remains electrically connected to the grounded pivot and the uniform added field's potential is included consistently. Put $s=\cos \theta$ and set the external electric potential to zero at the pivot. The external potential at the bob is then $ELs$, so the charging relation and full force balance becomeThe second equation follows by eliminating the wire tension between vertical balance and radial balance, including the outward magnetic Lorentz force. Expanding both equations to first order in $E$, and to leading order in $C{B}^2{L}^2/m$, with $s_0=g/({\omega }^{2}L)$, givesThe extra factor arises from the added field changing $Q$ by $-CEL{s}_0$ and therefore changing the radial magnetic force. That term is linear in $E$ too; the heavy-bob assumption alone does not justify retaining $Q_{0}E$ while dropping it. The two expressions refer to the two electrical approximations just specified. The expansion requires a small change relative to the initial radius and is nonuniform as that radius tends to zero.

For an electrical measurement, route the grounded connection through a sensitive electrometer operated near zero input voltage, and integrate the charging electric current when $B$ or $\omega$ is changed. Starting with the field off, the transferred electric charge is $Q_0$ in this ideal capacitance model. Reversing $B$ separates the induced signal from offsets. Alternatively, isolate the charged bob briefly and measure its electric potential with a high-input-impedance electrometer, accounting for the instrument's extra capacitance and for the fact that isolation changes the electrical boundary condition.

The characteristic impedance of a lossless transmission line is the ratio of the voltage to the current amplitude of a single forward-travelling wave. If $\mathcal{L}$ and $\mathcal{C}$ denote inductance and capacitance per unit length, the lossless telegrapher equations giveA backward wave has $I_{-}=-V_{-}/Z$ when current is positive along $+z$. Thus characteristic impedance is not a resistor's electrical resistance: it describes travelling-wave amplitudes and determines impedance matching, power transfer and reflections. A semi-infinite line, or a line terminated by a matching load, presents that impedance at its input.

At a junction with no lumped branch, voltage and electric current are continuous:Eliminating $V_t$ gives the voltage reflection coefficientFor real positive $Z_1$, incident and reflected waves occupy the same line, so their power ratio is $|r{|}^2$.

For the slightly mismatched coaxial cables, the change in characteristic impedance isConsequentlyThe requested power limit $0.01\%=10^{-4}$ implies $|\Delta b|/b<0.02\ln (b/a)$. Since $\ln (b/a)=2\pi nZ/Z_{\mathrm{vac}}=1.918$ for the stated line, the adjacent cables' outer radii must match within approximately $3.84\%$. If a tolerance is quoted separately as a symmetric manufacturing range $\pm \delta b$ and arbitrary worst-case pairs can differ by $2\delta b$, the corresponding individual tolerance is approximately $\pm 1.92\%$. The leading calculation is for small mismatches; the exact logarithmic expression gives slightly asymmetric positive and negative limits. The vacuum wave impedance used here is approximately $376.73\Omega$.^[ref]

A needle is a shunt electrical resistance $R_n$, so the load seen at its location is the continuation of the line in parallel with $R_n$: $Z_L=Z{R}_n/(Z+R_n)$. The local voltage reflection coefficient and transmission amplitude areFor periodic resistive loading of a transmission line, the reflection from the next needle travels an additional round trip $2L$ and crosses each previous needle twice. Taking the first needle as phase reference, the single-return approximation for periodic line reflections therefore sums toHere $v_p$ is the propagation phase velocity in the cable, $c/n$ for the stated nondispersive nonmagnetic filling. Squaring the amplitude gives the requested expressionFor a finite number $N$ of needles the amplitude has the additional numerator $1-[t^2{e}^{i\varphi }{]}^N$. The displayed semi-infinite result assumes sufficiently many needles, no terminal reflection, and neglects intrinsic propagation loss.

For real $r<0$, write $q=t^2$. The denominator of the power coefficient is $1+q^2-2q\cos \varphi$. It has constructive peaks at $\varphi =2\pi j$ and minima at $\varphi =(2j+1)\pi$, repeating with frequency spacing $\Delta \nu =v_p/(2L)$. Within this approximation,Weak individual shunts thus give narrow interference peaks, with phase width of order $|r|$.

The requested sum excludes paths that reflect repeatedly between different needles. It should not be read as an exact multiple-scattering result, particularly at the coherent peaks. For example the exact semi-infinite amplitude satisfiesAt $\varphi =2\pi j$ this equation reduces to $(R+1{)}^2=0$: the exact ideal infinite array has unit power reflection there, rather than the first-return approximation's value near $1/4$. This does not change the requested geometric-series derivation, but identifies its physical limitation.

Choose the guide axis as $z$, with $0<x<a$ and $0<y<b$. The electric field is along $y$, parallel to the walls' $b$ dimension. Take a homogeneous lossless nonmagnetic dielectric with relative permittivity $ϵ$, ideal conducting walls, and time dependence $e^{-i\omega t}$. The lowest transverse electric mode of a rectangular waveguide isAt $x=0,a$, this tangential electric field vanishes. At $y=0,b$ the electric field is normal to the wall; a nonzero value is supported by surface electric charge and violates no conducting boundary condition. The field has zero divergence. The associated normal magnetic field also vanishes at each wall, as verified below. The fundamental sine across the wider dimension minimizes the nonzero transverse wave number; ${\mathrm{TE}}_{00}$ is not a physical propagating mode of a hollow simply connected guide.

Below the waveguide cutoff frequency $f_c=c/(2a\sqrt{ϵ})$, put $\beta =i\kappa$ with $\kappa >0$ to select decay for $z>0$. The evanescent transverse electric mode in a rectangular waveguide then hasFor the fixed frequency in the question,The decay length increases monotonically and diverges as $ϵ\to 4^{-}$; its formal limit at $ϵ\to 0^{+}$ is $a/\pi$. At $ϵ=4$ the mode is at cutoff, with zero axial wave number. For $ϵ>4$, $\beta =(\pi /a)\sqrt{ϵ/4-1}$ is real and the mode propagates without axial attenuation in this ideal lossless model. There is no real finite evanescent decay length to continue beyond cutoff. The paper's graph notation $l$ denotes the same $\ell$ as its preceding formula.

From Faraday's law of induction, $\nabla \times \tilde{\mathbf{E}}=i\omega \tilde{\mathbf{B}}$. Writing $q=\pi /a$ givesThus $B_x$ vanishes on the $x$ walls and $B_y$ vanishes on the $y$ walls, meeting the normal magnetic boundary conditions. The fields also satisfy $\nabla \cdot \mathbf{B}=0$. For real $E_0$ the physical fields are $E_y=E_0\sin (qx)e^{-\kappa z}\cos \omega t$, $B_x=-(\kappa /\omega )E_0\sin (qx)e^{-\kappa z}\sin \omega t$, and $B_z=-(q/\omega )E_0\cos (qx)e^{-\kappa z}\sin \omega t$.

The time-averaged Poynting vector isAll products are imaginary because the electric field and magnetic field are in quadrature. A single decaying evanescent wave stores reactive energy without transmitting average power into a semi-infinite lossless guide. A finite below-cutoff section can transmit power when both independent evanescent solutions are present; that is a different boundary-value problem.

The skin effect also gives exponential field decay, but in a good metal $k$ has both real and imaginary parts and there is a real inward mean Poynting vector. That flux supplies Joule heating $\sigma |\tilde{E}{|}^2/2$. Here attenuation comes from geometric cutoff in a lossless region, with no conductivity and no dissipative power sink; exponential decay alone does not establish absorption.

For the attractive inverse-square force, choose zero potential energy at infinity, so $U=-A/r$. At either radial turning point $r_p,r_a$, the kinetic energy is purely tangential, with conserved angular momentum $J$. ThusIts two roots are the periapsis and apoapsis distances. Their sum is $r_p+r_a=-A/E$, whereas the ellipse has $r_p+r_a=2a$. Therefore $a=-A/(2E)$, with $E<0$ for a nondegenerate bound ellipse. The circular case follows by the double-root limit.

For the two-body problem, define $\mathbf{r}={\mathbf{r}}_B-{\mathbf{r}}_A$, total mass $M=m_A+m_B$, reduced mass $\mu =m_A{m}_B/M$, and center of mass $\mathbf{R}$. The relative equation is $\mu \ddot{\mathbf{r}}=\mathbf{F}$, while the isolated center of mass has $\ddot{\mathbf{R}}=0$. Once the relative orbit is known,The relative energy uses $\mu$ in its kinetic energy and $A=G{m}_A{m}_B$ in the gravitational potential energy.

In the equal-mass case the initial relative speed is $u$, $\mu =m/2$, and the center-of-mass velocity is half B's initial velocity. ThusNegative relative energy requires $u^2<4Gm/a_1$. The escape threshold is $u_E=2\sqrt{Gm/a_1}$; equality is marginal parabolic escape, so a strictly bound pair has $u<u_E$.

For $u=u_E/2=\sqrt{Gm/a_1}$, $E_{\mathrm{rel}}=-3G{m}^2/(4a_1)$ and the relative semimajor axis is $a=2a_1/3$. The initial transverse speed is below the relative circular value $\sqrt{2Gm/a_1}$, so the initial turning point is apoapsis. HenceChoose the initial separation along $+x$ from A to B, with B's initial velocity $+u\hat{\mathbf{y}}$. Conservation of angular momentum gives relative periapsis speed $3u$. At the opposite end of the ellipse the direction is $-\hat{\mathbf{y}}$, soThus B has the same speed as initially but the opposite velocity. A has $2u\hat{\mathbf{y}}$ at that instant. A value $3u/2$ for B would be its center-of-mass-frame speed, not the requested initial laboratory frame result.

The relative orbital period follows from Kepler's third law with total gravitating mass $2m$:In the center-of-mass frame each star follows an ellipse of semimajor axis $a_1/3$ and eccentricity $1/2$, with the center of mass at a focus; they remain opposite on their simultaneous position vectors. In the laboratory frame add the uniform drift $ut/2$ along $y$. The laboratory paths of an equal-mass binary are therefore not closed ellipses, even though their relative motion repeats.

In Lagrangian mechanics, choose independent generalized coordinates $q_i$ compatible with the mechanical constraints, write $L=T-U$, and apply the Euler-Lagrange equationEquivalently, the physical path makes the action stationary under variations with fixed endpoints. For ideal holonomic constraints, their unknown reaction forces need not appear in the reduced equations. Direct use of Newton's laws of motion instead resolves vector forces and reactions and enforces the constraints separately. Both approaches describe the same mechanics; the scalar Lagrangian is often shorter in curvilinear coordinates. Nonconservative generalized forces must be added separately. Here a drive maintains the prescribed azimuthal angular velocity $\Omega$.

The mass has polar speed $a\dot{\theta }$ and azimuthal speed $a\Omega \sin \theta$, so its kinetic energy is $m{a}^2({\dot{\theta }}^2+{\Omega }^2{\sin }^2\theta )/2$. The elastic potential energy is $\mu {\theta }^2/2$, and the height is $a\cos \theta$, giving gravitational potential energy $mga\cos \theta$. Thereforeand the Euler-Lagrange equation givesFor $\mu =2mga$, define $q={\Omega }^{2}a/g$ and $F(\theta )=q\sin \theta \cos \theta -2\theta +\sin \theta$, so $\ddot{\theta }=(g/a)F(\theta )$. Near zero,When $q<1$ the upright equilibrium is linearly stable and the small oscillations have ${\omega }_1^2=g/a-{\Omega }^2$.

For $q>1$, the graph of $F$ initially rises, whereas $F(\pi /2)=1-\pi <0$. There must therefore be a zero inside the interval. More strongly, for $0<\theta <\pi /2$,is strictly decreasing: $\cos \theta$ decreases and $\theta /\sin \theta$ increases because $\sin \theta -\theta \cos \theta >0$. The ratio goes from $q-1>0$ to $1-\pi <0$, giving exactly one zero ${\theta }_0$ in that interval. At the zero, $F^{\prime }({\theta }_0)<0$, so it is a stable equilibrium. The rotating torsion-stabilized inverted pendulum also has a symmetric stable angle $-{\theta }_0$ if both signs of tilt are permitted.

At ${\theta }_0=\pi /4$, setting $F({\theta }_0)=0$ gives $q/2-\pi /2+1/\sqrt{2}=0$, soFor $\theta ={\theta }_0+\xi$, the linearized equation is $\ddot{\xi }=(g/a)F^{\prime }({\theta }_0)\xi$. Since $F^{\prime }=q\cos 2\theta -2+\cos \theta$,At the threshold $q=1$ the linear restoring term vanishes, but the cubic term remains:The upright state remains nonlinearly stable, with no nonzero infinitesimal harmonic frequency. Small finite disturbances produce quartic oscillator motion whose frequency is proportional to its small angular amplitude; its period tends to infinity as that amplitude tends to zero. Zero linear stiffness is not the same as an absent restoring torque.

Use positive normal stress for tension. On a small rectangular material element, ${\tau }_i$ is the normal force per unit area on the face perpendicular to axis $i$, and $e_i=\Delta {\ell }_i/{\ell }_i$ is the corresponding normal strain. Young's modulus $E$ is the ratio of axial normal stress to axial strain in uniaxial stress; Poisson's ratio $\sigma$ is minus the transverse-to-axial strain ratio in the same test. Thus isotropic linear elasticity gives $E{e}_1={\tau }_1-\sigma {\tau }_2-\sigma {\tau }_3$ and cyclic permutations. The stress symbols here follow the question; $\sigma$ denotes Poisson's ratio, not a stress component.

For $e_2=e_3=0$, the two transverse compliance equations imply ${\tau }_2={\tau }_3=\sigma {\tau }_1/(1-\sigma )$. Substitution into the first givesThis uniaxial strain modulus is greater than the uniaxial-stress Young's modulus for the positive Poisson's ratio of interest, because the lateral constraints prevent free Poisson contraction.

For a small twist $\varphi$, the bonded faces' relative tangential displacement at radius $r$ is $r\varphi$. The layer's shear strain is $r\varphi /h$, so the shear stress is $Gr\varphi /h$. Integrating the torque givesThe minus sign denotes a restoring torque; the question quotes its magnitude. This uses the disk's polar second moment of area.

For tilt about H, choose $y$ in the plate perpendicular to H. The vertical displacement of the bonded top face is $y\theta$ to first order, giving normal strain $y\theta /h$. The thin layer is constrained laterally by bonding, so its leading normal traction is $My\theta /h$, not $Ey\theta /h$. Integrating its moment about H givessince the diametric second moment of area is $\pi a^4/4$. This is the specified thin-layer approximation to the torsion and tilt stiffness of a bonded elastic disk. Relaxation near the free cylindrical edge is neglected; for a nearly incompressible layer it can limit quantitative accuracy unless the thinness is sufficient.

The thin square plate's moment of inertia about V is $I_V=2m{b}^2/3$. Its two perpendicular in-plane central moments are equal and sum to $I_V$, so the diagonal H has $I_H=m{b}^2/3$. Dividing the two angular stiffnesses by the appropriate moments of inertia givesThe geometry factors cancel in their ratio. Using the supplied relation between shear modulus and Young's modulus,The result lies just below the incompressible isotropic limit $1/2$, as expected for a rubber with a large compression stiffness compared with its shear modulus. It is an inference within the bonded-layer model, rather than an independently measured material constant.

For a Newtonian fluid with velocity $u(z)\hat{\mathbf{x}}$, $\tau =\eta {\partial }_{z}u$ is the $x$ component of traction on a plane whose normal points along $+z$. The opposite face has the opposite traction. The coefficient $\eta$ is dynamic viscosity, with units $\mathrm{Pa}\mathrm{s}$. Faster layers transfer positive $x$ momentum toward slower neighbours; the internal shear stress resists relative sliding.

Take a thin fluid element with horizontal area $A$ and thickness $dz$. The difference between its upper and lower viscous tractions is $A\eta u_{zz}dz$. The pressure gradient supplies $-A{p}_{x}dz$. Incompressibility is satisfied because $u$ has no $x$ dependence, and the convective acceleration vanishes: $(\mathbf{u}\cdot \nabla )\mathbf{u}=u{\partial }_x\mathbf{u}=0$. Equating net force to mass times acceleration therefore gives the relevant component of the Navier-Stokes equation:For steady flow, integrate $\eta u^{\prime \prime }=p^{\prime }$ twice and impose the no-slip boundary condition at each plate, $u(0)=0$ and $u(a)=u_0$. The Couette-Poiseuille flow profile isThe linear term is Couette flow and the parabolic term is plane Poiseuille flow. A negative $p^{\prime }$ drives flow in $+x$. Integrating the profile across the gap gives the volumetric flow per unit width,It has units ${\mathrm{m}}^2{\mathrm{s}}^{-1}$; multiplying by the plate width gives the actual volume per time.

With zero pressure gradient, substituting $u=\mathrm{Re}[A{e}^{i(kz-\omega t)}]$ gives $-i\rho \omega =-\eta k^2$. Select the root that decays into fluid occupying $z>0$:The other sign applies to decay in the opposite half-space. With a real wall velocity amplitude $U_0$, the established periodic solution of Stokes' second problem isThe oscillation's amplitude decays exponentially and its phase lags with distance from the plate. This is diffusion of tangential momentum, not a freely propagating undamped elastic wave. Initial transients must have decayed before this purely periodic expression applies.

The fluid-displacement amplitude and the wall-displacement amplitude have the same ratio as their velocity amplitudes, since both are divided by the common $\omega$. Thus for the specified amplitudes,Using the standard approximate water density $\rho \simeq 10^3\mathrm{kg}{\mathrm{m}}^{-3}$, explicitly an approximation rather than another value printed on the paper, gives $\delta =0.564\mathrm{mm}$ and $z\simeq 2.21\mathrm{mm}$ from either face of a sufficiently large plate. The viscous penetration depth varies as ${\rho }^{-1/2}$ if a more precise water density is used. The ideal parallel-flow solution ignores finite plate edges and tank walls.

For entropy gain from impact and thermal equilibration, let the large lake be a reservoir at an initially unspecified temperature $T_L$, in kelvin, and take the final meteorite temperature to be $T_L$. Its initial temperature is $T_i=1273.15\mathrm{K}$. The meteorite loses thermal energy $mc(T_i-T_L)$, and its kinetic energy $m{v}^2/2$ is also dissipated into the lake. ThusandBoth heat equilibration and dissipation contribute positively. For the requested estimate, assume a lake temperature $T_L=300\mathrm{K}$. This gives $\Delta S_{\mathrm{universe}}\simeq 2.93\times 10^3\mathrm{J}{\mathrm{K}}^{-1}$. Taking $273$–$300\mathrm{K}$ instead gives about $3.3$–$2.9\mathrm{kJ}{\mathrm{K}}^{-1}$; this shows the estimate's sensitivity without asserting an unprovided lake temperature. The boxed expression retains the full temperature dependence. The model uses the stated constant heat capacity, no phase change of the meteorite, and a lake large enough that its temperature change is negligible.^[ref]