Write , with . Applying polar coordinates and Fubini's theorem gives, for bounded nonnegative measurable ,
Also by absolute continuity of measures. Extending by zero on the negative half-line proves
Both radial pushforward measures have compact support of a measure, so Brenier's theorem on the line supplies a unique optimal transport plan and a transport map , unique -almost everywhere, with
In one dimension, the gradient of a convex function is a nondecreasing function. Nonnegativity holds -almost everywhere because is carried by . To choose a nonnegative, nondecreasing representative everywhere, let be the radial cumulative distribution functions and use the monotone rearrangement
At set , and at use the largest point of the support of a measure . These choices preserve monotonicity and boundedness. Since is continuous, has the uniform distribution on when has law , and the quantile function sends this uniform distribution to . Thus this representative has the required pushforward measure, even if the target has atoms.
Polar coordinates give speed squared , so unit mass has kinetic energy . The radial force is , giving after choosing the additive constant zero. Therefore is the stated expression.
The two Euler-Lagrange equations are
The second gives conservation of angular momentum . The mechanical energy is
Direct differentiation expresses its derivative as
Thus energy is conserved for the motion wherever .
Write the complex potential of a two-dimensional potential flow as . Its complex velocity is . Here
so the far-field motion has velocity and speed . In polar coordinates the stream function, shifted by a harmless constant, is
It is zero on , proving the cylinder is a streamline.
Stagnation points are zeros of , hence roots of . Put . The two roots are
Both have modulus and are distinct. Their angles satisfy , namely and . The tangential velocity on the cylinder is , confirming the same zeros. Streamlines are level curves of ; the positive circulation places both stagnation points on the upper side and increases the flow speed below the cylinder.
Figure 1.
Streamlines and the two surface stagnation points for positive subcritical circulation
. Details.
The annulus is open and hence Lebesgue measurable; multiplying its indicator by the continuous Gaussian gives a measurable function. It is bounded by the indicator of , whose Lebesgue integral is . Thus it is Lebesgue integrable by comparison.
Use the change of variables formula with polar coordinates , . Removing one radial cut makes this a smooth one-to-one coordinate map; the cut is a null set and does not affect the integral. The absolute Jacobian determinant is . By Fubini theorem,
The antiderivative is , which also makes the sign of the answer transparent.
For the isolated line vortex, use polar coordinates to obtain
Thus for nonzero the streamlines are circles and their counterclockwise circulation is . If the fluid is stationary and the streamline description degenerates.
Place an opposite image vortex at . The method of images gives
On , the distances to the two vortex positions agree, so . Hence the wall is a streamline. Directly the wall velocity has and
For the stipulated steady flow, Bernoulli's equation in the irrotational fluid outside the vortex core gives
The force exerted by the fluid on this wall points in the direction for positive pressure. Relative to uniform pressure , however, its change is in the direction. Using the supplied integral, the force change per unit span is
Thus the wall is attracted towards the vortex, with this magnitude if the other side remains at . Subtracting the uniform pressure is essential for an infinite wall. A freely moving vortex would translate parallel to the wall under its image-induced velocity; the calculation here is for the stated steady vortex configuration, with its position held fixed.
For variations vanishing at the endpoints, the Euler-Lagrange equations are
Differentiate the stated energy expression using the chain rule. The acceleration terms cancel, leaving
Along a stationary trajectory the two brackets vanish. Thus , and in the absence of explicit time dependence .
In plane polar coordinates, the squared speed is , so the kinetic energy is . An inverse-square central gravitational force has potential energy proportional to . With the sign convention printed in this paper, write , giving
For attraction to a star of mass , this convention has . Keeping real and signed is essential; no sign is changed in differentiating the printed functional.
Since , , and , the equations of motion are
The second equation gives the conserved angular momentum ; equivalently . The autonomous-functional identity gives a second conserved quantity,
This is the total mechanical energy; for attractive gravity the potential contribution is negative because .
The function is nonnegative and continuous on the positive quadrant. In polar coordinates its product integral is, by the nonnegative change of variables formula and Tonelli's theorem,
The angular factor is at most one, while the supplied radial integral is , so the product integral is at most and is finite. This establishes integrability before any signed interchange is used.
Let denote the angular integral. The polar evaluation is . In Cartesian coordinates Fubini's theorem instead factors the integral as
using , . All factors converge: their powers exceed at zero and the exponential controls infinity. Equating the evaluations gives
Take a smooth variation with on the boundary. Differentiating under the integral gives the first variation
Integration by parts removes the derivatives of , and the boundary contribution vanishes. The fundamental lemma of the calculus of variations therefore gives
For a regular constraint , the Lagrange multiplier rule says that an extremal satisfies for a constant . Equivalently, apply the preceding Euler–Lagrange equation to . Regularity means that is not identically zero as a linear functional on the admissible variations. For the quadratic constraint below this holds at every nonzero admissible ; a zero constraint leaves only .
In polar coordinates the constrained integrand is
Its Euler–Lagrange equation is , or
The field is single-valued and periodic in . Substitution with gives
The centre is part of the disk, so one must also select the regular radial solution: for , smoothness gives and the regular derivative behaviour; for , is bounded and . A singular solution at zero does not define an admissible smooth field on the disk.
For fixed , let , so and for . The angular derivative term has integral , including zero when . Consequently
The one-dimensional Euler–Lagrange equation for is
exactly the same radial equation. The common factor cancels.
For two regular solutions with the same and different multipliers, multiply the equation for by , multiply that for by , subtract and integrate. The result is
At , both functions vanish; at zero, regularity makes the boundary term tend to zero. Since the multipliers are distinct,
This is weighted Sturm–Liouville orthogonality with weight .
Finally, multiply one radial equation by and integrate by parts. The same endpoint conditions yield
Multiplication by proves
For any nonzero admissible solution the energy on the left before subtraction is positive, so its multiplier is the positive Rayleigh quotient .
Use the Fourier transform convention , initially in a class where the integral transforms and subsequent limits are defined. Transforming Laplace's equation gives on real . Decay for selects
The boundary trace is the sum of its negative-half-line unknown and the given positive-half-line datum. The normal derivative is zero on the negative side. Therefore
For the convention, negative support gives a integral transform analytic for , and positive support gives one analytic for , subject to the appropriate endpoint interpretation. The normal trace in this particular problem requires a finite-part interpretation, detailed below.
For the Wiener-Hopf method, regularize the real-axis kernel by
Take the branch cut of vertically downward from and that of vertically upward from . Their product on the real axis is the positive square root . The branch of has its cut downward from zero. Dividing the regularized equation by gives
The forcing is analytic in , between its two cuts. Decompose by Cauchy integrals on two horizontal lines within that strip. Separation gives
when the endpoint growth class makes the common entire remainder bounded. For one can take the minus Cauchy integral on , with , obtaining
The integrand is at infinity, so the Cauchy integral converges there. In the limit choose , for example ; a fixed positive cannot remain between zero and . With denoting the limiting lower-analytic branch and with the prescribed limiting interpretation,
The numerator square root has a cut going upward from ; the power has a cut going downward from zero. The contour passes between them. The pole is below this contour when is in the lower integral transform domain. These branch choices are essential, rather than independent principal powers of the factors.
Figure 1.
The Wiener-Hopf contour lies between the upward square-root cut and the downward quarter-power cut
. Details.
There is an important mathematical limitation in the source's ordinary integral definition of . An explicit harmonic function particular solution in polar coordinates is
Each separated polar power is harmonic function; at it has the required Dirichlet boundary condition value, and at its angular derivative vanishes, giving the required zero normal derivative. It also decays at infinity. Its unknown boundary traces are
The last function is not locally integrable at zero. Therefore the printed ordinary integral diverges, even for . The integral transform calculation is valid for this solution only after specifying a Hadamard finite-part integral, or equivalently a distribution boundary trace.
For the canonical homogeneous extension of this power, the two transformed traces are
The first is analytic below and the second above. The second formula follows from and the finite-part continuation of the Laplace transform. Substituting their real-axis boundary values directly verifies the boxed integral transform relation for both signs of , so the qualification can be checked without relying only on formal factorization.
There is also homogeneous corner freedom:
satisfies zero Dirichlet boundary condition data on the positive axis, zero Neumann boundary condition data on the negative axis and decay at infinity. Its negative trace has integral transform , corresponding to the term . Higher singular modes would also satisfy the homogeneous boundary conditions for integers , but are excluded if the negative trace is required to be locally integrable. Thus the one-constant formula belongs to a specified corner-growth class. The original boundary conditions alone, without that class, do not assert uniqueness. These facts form the finite-part normal trace for mixed quarter-power boundary data.
Use a reference frame moving with Jupiter, neglect the spacecraft's effect on the planet and neglect the Sun during the short encounter. The Newtonian central force makes the orbit planar and conserves the angular momentum . In plane polar coordinates, the kinetic energy is . Choosing zero gravitational potential energy at infinity gives
The terms are respectively radial kinetic energy, transverse kinetic energy, and attractive gravitational potential energy. The second term is often combined with the third into the radial effective potential. Far from Jupiter,
Here is the perpendicular impact parameter of the incoming asymptote, not the eventual periapsis distance.
For a differentiable central potential, the force is . The angular momentum is , so
This proves conservation of angular momentum. For , the trajectory lies in the fixed plane perpendicular to . In plane polar coordinates, . Substituting this into the kinetic energy gives the radial energy equation
The second term is the centrifugal contribution to the effective potential. Its radial derivative reproduces . Purely radial motion is the case.
For the satellite, write , with Earth's mass and the satellite's mass. At apogee and perigee, the radial velocity is zero, so angular momentum conservation and conservation of energy give
Eliminating yields
For an elliptic orbit, . Therefore the vis-viva equation at apogee is
Substitution into the orbital energy gives . The circular case follows by continuity, or directly from the radial force balance.
An instantaneous tangential retrograde burn at apogee leaves the position unchanged. For the speed change , the orbital energy change is
Differentiating gives to first order. Since , the vis-viva equation implies . Hence
The factor of two distinguishes the major axis from the semimajor axis. As a finite-change check, the new semimajor axis is exactly
within the instantaneous-burn, unchanged-mass model. Expanding this expression gives the stated result.
The massless inextensible string has hanging length , measured downward from the table. Thus . In polar coordinates, the tabletop mass has squared speed . Adding the two kinetic energies and setting the hanging gravitational potential energy to gives
and the Lagrangian is
The canonical momenta are and . Performing the Legendre transform gives the Hamiltonian
It equals the total mechanical energy and is conserved because has no explicit time dependence. Also is a cyclic coordinate, so is conserved angular momentum about the vertical through the hole. The table and the string guide supply external forces; conservation of the system's total linear momentum does not follow. These are the two independent generic integrals needed for the reduced tabletop mass coupled to a hanging mass.
For motion in a central potential, conservation of angular momentum gives . In polar coordinates, the conserved mechanical energy is
Consequently, for ,
The sign in distinguishes inward and outward motion. Choosing the orientation so , it is also the sign of . The expression under the square root must be nonnegative. If , the motion is radial and is constant, so cannot be parametrized by this way.
Choose the orientation . To leading order in the radial disturbance, the azimuthal angular frequency is , so
For either orientation the positive frequency ratio is . At , the ratio is : there is one radial oscillation per revolution. The bound Kepler orbit is an ellipse with the force centre at a focus, with one closest and one furthest approach per revolution. Its exact polar coordinates equation is
At , the ratio is : there are two radial oscillations per revolution. The isotropic harmonic oscillator trajectory is an ellipse centred on the force centre. An exact parametrization is
Here repeats twice during one full circuit. The two different placements of the force centre distinguish these ellipses. These are the two all-bound-orbit cases of Bertrand's theorem.
Figure 1.
Nearly circular ellipses: one radial cycle about a focus for the inverse-distance potential, and two about the centre for the quadratic potential
. Details.
For image charges in a sixty-degree conducting wedge, use polar coordinates in a cross-section normal to the planes' intersection, putting the physical wedge at and the real point charge at . Repeated reflection in the two grounded faces produces five image charges, all at radius and the same coordinate along the intersection:
Equivalently, the positive copies are at and the negative copies at , for . Only the original point charge lies inside the physical wedge.
Figure 1.
Five image charges outside the physical sixty-degree conducting wedge
. Details.
For any observation point on either face, reflection pairs every positive point charge with an equidistant negative one. Their Coulomb potentials cancel, so both grounded faces have zero electric potential. The electric potential also tends to zero at infinity. The uniqueness theorem for the electrostatic potential then identifies this method of images construction with the physical solution.
Point-charge field and superposition. Coulomb's law gives the Coulomb force on a test electric charge a displacement from a source point charge :
The electric field is force per infinitesimal test charge, so
It points outward for positive and inward for negative . For several point charges, the superposition principle adds their fields as vectors:
This formula excludes the source points themselves, where the ideal point-charge field is singular.
Exact field components of the pair. Put
The displacement vectors from the positive and negative charges to the observation point are and . Including the signs of the charges gives
Reading off both components,
For example, on the perpendicular bisector , the two components cancel and the components add toward the negative charge, checking both relative signs.
Far-field polar components. The pair's electric dipole moment points from negative to positive charge:
Use polar coordinates , , with measured from that positive dipole axis. The supplied leading electric dipole field becomes
Taking the dot products with the radial and angular unit vectors gives
Therefore
The unit vectors for polar coordinates are orthonormal, so the magnitude is
The units are , since the electric dipole moment has units .
Angular sketch and comparison. At fixed large , the magnitude has period . It is largest on the dipole axis, , with value , and smallest on the perpendicular directions, , with value . The field magnitude never vanishes at finite in this dipole approximation:
Figure 1.
The angular dependence of the far-field electric dipole magnitude at a fixed radius
. Details.
The neutral charge pair has no net electric charge, so its leading far field falls as and depends on direction. An isolated nonzero point charge has an field whose magnitude is isotropic, with purely radial direction. The dipole also generally has a nonzero component.
A molecular example is the interaction between permanent electric dipoles of hydrogen chloride molecules, whose unequal charge distributions give each molecule a nonzero electric dipole moment. The electric dipole-dipole interaction depends on their relative orientations:
For two parallel dipoles pointing along their line of separation, , illustrating an attractive molecular configuration.
Set . In polar coordinates, the event is , the angular sector
Because has the continuous uniform distribution on an interval of length , the probability is the sector's angular length divided by :
The symmetry of the standard normal distribution in makes the lower sector equally likely. Their disjoint union is , apart from , an event of probability zero. Therefore the requested cumulative distribution function is
This is the unit-scale half-Cauchy distribution. At there is no atom, and as the cumulative distribution function tends to one.
Figure 1.
Gaussian angular probability and the support of a transformed exponential pair
. Details.
The left panel uses , giving a sector of angle and probability . The right panel displays the domain for the next change of variables for probability densities; a vertical section through it determines the conditional range.
The joint probability density of the two independent standard normal random variables is . The polar coordinates map has absolute Jacobian determinant , so the change of variables for probability densities gives
At the origin the angle may be assigned arbitrarily: that point has probability zero. The joint probability density factorizes as a Rayleigh distribution density and a continuous uniform distribution density . Hence and are independent.
For , substitute to find
Therefore , equivalently , and applying a measurable function to preserves its independence from . Explicitly, with ,
Thus and independently. This polar decomposition of a planar standard normal vector separates the uniformly distributed direction from its random squared distance.
Write . The first Hermite polynomials give and proportional, with the same constant, to and . Thus in Cartesian coordinates the new wavefunctions are
Using polar coordinates , gives
The orbital angular momentum is , which immediately verifies its eigenvalues . Multiplying these spatial wavefunctions by for and for gives the corresponding full stationary states.
Use the Christoffel symbol formula with its first index raised:
The diagonal metric tensor has and , with inverse components and . Since and all angular derivatives vanish,
For example, . The other components vanish. The formula supplied in the paper gives the first-kind coefficients ; multiplication by the inverse metric tensor is needed for the requested second-kind coefficients. These polar coordinates are understood away from .