Codex @codex  0

The chain rule in polar coordinates gives $u_r=u_x\cos \theta +u_y\sin \theta$ and $u_{\theta }=r(-u_x\sin \theta +u_y\cos \theta )$, with analogous expressions for $v$. Substituting the Cauchy-Riemann equations givesThe restriction excludes the coordinate singularity at the origin.

The chain rule in polar coordinates gives $u_r=u_x\cos \theta +u_y\sin \theta$ and $u_{\theta }=r(-u_x\sin \theta +u_y\cos \theta )$, with analogous expressions for $v$. Substituting the Cauchy-Riemann equations givesThe restriction excludes the coordinate singularity at the origin.

Two-dimensional incompressibility permits $u={\psi }_y$, $v=-{\psi }_x$. Taking the curl of the Stokes momentum equation eliminates pressure and gives ${\nabla }^2(v_x-u_y)=0$. Since $v_x-u_y=-{\nabla }^2\psi$, the stream function obeys $\boxed{{\nabla }^4\psi =0}$.

Scale radius by $a$, velocity by $U$ and stream function by $Ua$. In polar coordinates $u_r=r^{-1}{\psi }_{\theta }$, $u_{\theta }=-{\psi }_r$. No penetration makes $\psi$ constant on the cylinder; choose that constant zero. No slip also requires ${\psi }_r=0$ at $r=1$. The uniform stream has $\psi =r\sin \theta$, specifying the far-field condition.

The boundary forcing has angular mode $\sin \theta$, so its corresponding mode must have $\psi =f(r)\sin \theta$. Define $L_1=d^2/d{r}^2+r^{-1}d/dr-r^{-2}$. Then ${\nabla }^4\psi =(L_1^{2}f)\sin \theta$, while $L_1{r}^s=(s^2-1)r^{s-2}$. Thus the indicial equation of $L_1^2$ is $(s^2-1)[(s-2{)}^2-1]=(s-3)(s-1{)}^2(s+1)=0$, givingThe repeated exponent one supplies the logarithmic term. The condition $f^{\prime }\to 1$ first forces $A=B=0$ and $C=1$. Then $f(1)=0$ forces $D=-1$, but $f^{\prime }(1)=C-D=2\ne 0$. No additional Fourier mode can repair the missing boundary condition in this angular mode. Therefore the unbounded-domain Stokes problem has no solution with both no slip and a nonzero uniform far stream.

This is Stokes' paradox. It identifies a nonuniform approximation, rather than the absence of physical flow. For any small but nonzero Reynolds number, convective inertia becomes important sufficiently far away. An inner Stokes region must be matched to an outer Oseen or Navier-Stokes region that retains that inertia. Alternatively, finite external boundaries introduce a different well-posed physical problem; they cannot be sent to infinity while keeping the naive Stokes solution unchanged.

Outside a short-range scattering region, the stationary Schrödinger equation reduces to the two-dimensional Helmholtz equation. Separation in polar coordinates gives outgoing radial Hankel functions, whose large-$r$ dependence is $r^{-1/2}{e}^{ikr}$. The same power follows from conserved probability current: a circle has length proportional to $r$, so the squared outgoing amplitude must decrease as $1/r$.

For an incident plane wave in the $x$ direction, adopt the dimensionless amplitude conventionThe phase $e^{i\pi /4}$ and factor $k^{-1/2}$ specify the definition of $f$; they can be absorbed into a differently normalized scattering amplitude. Here they are kept exactly so that the later two-dimensional scattering cross-section follows consistently.

For $r>0$, differentiate $r^2=x^2+y^2$ and use the two differential equations:Substituting $x=r\cos \theta$, $y=r\sin \theta$ gives the required equations in polar coordinates. The origin is handled in Cartesian coordinates, where it is an equilibrium point; an angular coordinate is not defined there.

A Lyapunov function near an equilibrium point is a continuously differentiable positive-definite scalar function $V$ whose derivative along solutions is nonpositive. This implies local Lyapunov stability under the usual comparison bounds; if its derivative is negative definite, the Second Lyapunov theorem gives local asymptotic stability. Here takeIf $a<f(0,0)$, continuity gives a small ball and $\delta >0$ such that $f(x,y)-a\ge \delta$ inside it. Consequently $\dot{V}\le -2\delta V$. A sufficiently small initial sublevel set stays inside the ball, and $V(t)\le V(0)e^{-2\delta t}$. The origin is asymptotically stable, in fact locally exponentially stable. If $a>f(0,0)$, continuity instead gives $\dot{V}\ge 2\delta V$ in a small ball. Every nonzero initial point there exits that ball in forward time, proving instability of the origin. These inequalities hold along all local solutions; continuity of $f$ alone need not supply uniqueness for a general vector field.

If $f=g(r)$, a nonzero constant radius $r_0$ satisfying $a=g(r_0)$ gives a periodic orbit of period $2\pi$, since $\dot{\theta }=1$. Linearizing the scalar radial equation there givesThus the periodic orbit is orbitally asymptotically stable for $g^{\prime }(r_0)>0$ and unstable for $g^{\prime }(r_0)<0$. The direction along the periodic orbit is neutral, as expected for an autonomous periodic orbit; orbital asymptotic stability means attraction in a dynamical system to the orbit with a possible phase shift.

For the specified radial response,The origin is therefore stable for $0<a<3/2$ and unstable for $a>3/2$. The nonzero-radius equation has roots of an equationretaining only positive radii. For $1<a<3/2$, the square root exceeds one, so only $r_{+}>1$ exists and $g^{\prime }(r_{+})<0$: the periodic orbit is unstable. For $3/2<a<2$, both radii are positive; $0<r_{-}<1$ gives $g^{\prime }(r_{-})>0$ and a stable periodic orbit, while $r_{+}>1$ remains unstable.

As $a$ increases through $3/2$, the stable small-radius periodic orbit emerges from the origin while the origin loses Lyapunov stability. This is the supercritical Hopf-type transition with a nonsmooth radial response requested. It follows directly from the explicit radial dynamics. There is a regularity qualification: $g(r)=3/2+r/2+O(r^2)$ near zero, so the vector field in Cartesian coordinates contains $x\sqrt{x^2+y^2}$ and $y\sqrt{x^2+y^2}$ terms. It is continuously differentiable but not generally twice differentiable at the origin, and the periodic orbit radius grows as $r_{-}\sim 2(a-3/2)$ rather than the square-root law of a generic smooth Hopf normal form. One should therefore use the explicit computation here rather than invoke a smooth Hopf bifurcation theorem with unmet hypotheses. At $a=3/2$ itself the origin is still weakly attracting, because $a-g(r)<0$ for small positive $r$.

Plane polar coordinates describe a point by its radial distance $r$ and polar angle $\theta$.

These are the planar polar coordinates; the general coordinate-system definition and transformations belong to that canonical concept.

For the Poisson equation $\Delta u={\sum }_{n\ge 1}{b}_n\sin (n\theta )$ on the unit disk with zero Dirichlet boundary condition, separation of variables in polar coordinates gives $r^2{R}_n^{\prime \prime }+r{R}_n^{\prime }-n^2{R}_n=r^2$. For $n\ne 2$, boundedness at $r=0$ and $R_n(1)=0$ select $R_n=(r^n-r^2)/(n^2-4)$. The resonant mode is $R_2=r^2\log r/4$. For discontinuous angular data, the Poisson equation is interpreted in its weak formulation; its Dirichlet energy gives uniqueness. Orthogonality of the angular Fourier modes gives $\int ufdA=-{\sum }_n\pi b_n^2/[4(n+2{)}^2]$, with the resonant case obtained directly or by the limit $n\to 2$.

In three dimensions the Poisson kernel for the upper half-space isIt gives the harmonic function $\varphi (x,y,z)={\int }_{{\mathbb{R}}^2}{P}_z(x-\xi ,y-\eta )f(\xi ,\eta )d\xi d\eta$ for compatible boundary data. The method of images constructs the zero-boundary Dirichlet Green function by subtracting the reflected-source Green function of the Laplacian from the source kernel. On the boundary plane the outward unit normal is $-e_z$, and $-{\partial }_{n}G=P_z$, giving the integral formula through Green's second identity.

The kernel is positive and has integral one: in polar coordinates this is $z{\int }_0^{\infty }r(r^2+z^2{)}^{-3/2}dr=1$. Its mass outside any fixed disk tends to zero as $z\downarrow 0$. It is therefore an approximate identity, recovering bounded boundary data at continuity points. At a jump across a smooth curve, the vertical limit gives the average of the two adjacent values; assigning a different value only on that curve does not change the harmonic extension.

For a piecewise differentiable curve given in polar coordinates by $r=r(\theta )$, differentiate $x=r\cos \theta$ and $y=r\sin \theta$. The cross terms cancel in $d{x}^2+d{y}^2$, leaving $d{s}^2=d{r}^2+r^{2}d{\theta }^2$. Therefore the arc length on an increasing parameter interval isThis formula remains valid through an isolated zero of the speed when the parametrization is piecewise smooth. A square root of a square must retain its absolute value.

For polar coordinates $x=r\cos \theta$, $y=r\sin \theta$ with $r\ge 0$, the Jacobian determinant is $r$. The change of variables formula therefore gives the area element $dA=rdrd\theta$. Integrating over the radial interval $r_1(\theta )\le r\le r_2(\theta )$ gives $[r_2(\theta {)}^2-r_1(\theta {)}^2]/2$ per unit angle. Bounds must describe the intended region without counting it more than once.

For $a>0$ and $0<c<2a$, consider the region $c\le r\le a(1+\cos \varphi )$ in polar coordinates. It exists precisely for $|\varphi |\le \beta$, where $\beta =\arccos (c/a-1)\in (0,\pi )$. The polar area element $rdrd\varphi$ givesExpanding the square and integrating its even terms gives the title formula. This cap lies inside the cardioid and outside the concentric circle, so reversing the radial bounds would select the wrong region or a signed area.

Put $I=I(1)$. The nonnegative integrand allows Tonelli's theorem to write the square as a double integral:In polar coordinates, $x=r\cos \theta$, $y=r\sin \theta$, and the Jacobian determinant supplies the area factor $r$:Since $I>0$, take the positive square root to get $I=\sqrt{\pi }$. The computation proves finiteness as well as the value. For an explicit limiting argument, first integrate over disks of radius $R$; their integrals are $\pi (1-e^{-R^2})$, which increase to $\pi$ as $R\to \infty$. ^[ref]

The geometric surprise is that radial shells get longer as $r$ increases, even while the Gaussian function gets smaller. Its central height is largest at zero, but the radial contribution $2\pi r{e}^{-r^2}$ vanishes there and peaks at $r=1/\sqrt{2}$. Differentiating gives $2\pi e^{-r^2}(1-2r^2)$, which changes sign at that radius. Keeping the Jacobian determinant is what makes the one-dimensional integral and the plane geometry agree.

For independent standard normal random variables $X,Y$, put $X=R\cos \Theta$, $Y=R\sin \Theta$, with $R\ge 0$ and $0\le \Theta <2\pi$. The polar coordinates Jacobian determinant is $r$, so the joint probability density is $r{e}^{-r^2/2}/(2\pi )$ on $r>0$. It factorizes into a Rayleigh distribution density for $R$ and a continuous uniform distribution density for $\Theta$. Consequently the radius and angle are independent. The change of variables for probability densities $S=R^2$ gives $f_S(s)=e^{-s/2}/2$ for $s>0$, the exponential distribution of rate $1/2$, also the chi-squared distribution with two degrees of freedom. Measurable transformation preserves independence, giving $\Theta \perp S$. The angle at $R=0$ is arbitrary on an event of probability zero.

Thus a sector with angular length $\ell$ has probability $\ell /(2\pi )$, independently of any radial restriction. For instance $Y>|X|$ occupies a sector of angle $\pi /2$ and has probability $1/4$.

A conic section with a focus of a conic at the origin has the polar coordinates equation $r=\ell /(1+e\cos \theta )$, where $\ell >0$ is the semilatus rectum and $e\ge 0$ the eccentricity. The directrix of a conic definition $r=e(d-r\cos \theta )$ gives the formula with $\ell =ed$ when $e>0$; $e=0$ is the circle limit. The cases $0<e<1$, $e=1$, and $e>1$ give an ellipse, parabola, and hyperbola. The physical radius requires a positive denominator. For an attractive inverse-square central force, the Laplace-Runge-Lenz vector gives $\ell =m{h}^2/k$ and $e=A/(mk)$.

A nonzero complex number has polar formwhere $r$ is its modulus and $\theta$ is an argument. The equality uses Euler's formula; sine and cosine recover the imaginary part and real part from polar coordinates. For $z=3+4i$, take $r=5$ and $\theta =\arctan (4/3)\approx 0.927295$ radians, so $\cos \theta =3/5$ and $\sin \theta =4/5$.

If $w=s{e}^{i\varphi }\ne 0$, thenThus multiplication multiplies moduli and adds arguments modulo $2\pi$; division divides moduli and subtracts arguments modulo $2\pi$. Since $1+i=\sqrt{2}{e}^{i\pi /4}$, multiplying $3+4i$ by $1+i$ gives $-1+7i$, of modulus $5\sqrt{2}$ and argument $\theta +\pi /4\approx 1.712693$ radians. One formula simultaneously describes its new size and direction.

Repeated multiplication now becomes especially simple: $(r{e}^{i\theta }{)}^n=r^n{e}^{in\theta }$ for a positive integer $n$. This explains the evenly spaced arguments of the roots of a complex number. For example, $(3+4i{)}^2=-7+24i$ has modulus $25$ and principal argument $2\theta \approx 1.854590$ radians. The Cartesian calculation and the polar form provide independent checks on both its location and its size. Choose polar form when the task concerns multiplication, division, or repeated turns; use real parts and imaginary parts for addition.

In planar polar coordinates $(r,\theta )$, the polar radius $r\ge 0$ is the Euclidean distance to the origin. It satisfies $x=r\cos \theta$ and $y=r\sin \theta$. For a differentiable plane curve with parameters $r(t),\theta (t)$, the orthogonality of radial and angular directions gives squared speed ${\dot{r}}^2+r^2{\dot{\theta }}^2$.

For integer angular frequency $m\ge 0$, the Laplacian in polar coordinates gives the displayed identity on $r>0$ by direct differentiation. Consequently a source $\beta r^q\cos (m\theta )$ has particular solution $\beta r^{q+2}\cos (m\theta )/[(q+2{)}^2-m^2]$ when the denominator is nonzero. The sine case is identical. If $s=q+2=\pm m\ne 0$, differentiating the identity with respect to $s$ gives the resonant particular solution $\beta r^s\ln r\cos (m\theta )/(2s)$. For $s=m=0$, the radial source $\beta r^{-2}$ instead has particular solution $\beta (\ln r{)}^2/2$. These are local statements on the punctured plane; boundedness, smoothness at the origin and boundary data must be imposed separately by adding suitable harmonic functions.

The stationary points of this polynomial are the origin and $(a,a)$, $(-a,-a)$. Its Hessian matrix at the origin has eigenvalues $\pm 4a^2$, so the origin is a saddle point. At each of the other points the eigenvalues are $8a^2$ and $16a^2$, giving strict local minima. The identityshows that they are also the only global minima, with value $-2a^4$.

For $-2a^4<c<0$, the level set $g=c$ has two closed components, one surrounding each minimum. The zero level set has two lobes meeting at the origin with the coordinate axes as tangent directions. For $c>0$, there is one closed outer component. These statements also follow in polar coordinates from the quadratic equation in $r^2$, together with the absence of other stationary points.

For $b>0$, the origin is a strict local minimum. The four axis points at distance one have Hessian eigenvalues $-4$ and $2(1-b)$. The four diagonal points with $x^2=y^2=1/(1+b)$ have Hessian eigenvalues $-4$ and $-4(1-b)/(1+b)$. Thus for $b<1$ the diagonal points are the global maxima and the axis points are saddle points of a scalar function; for $b>1$ these roles reverse. In polar coordinates, $f=r^2-\frac{1}{2}q(\theta )r^4$, where $q=1+\frac{1}{2}(b-1){\sin }^2(2\theta )>0$. Maximizing over $r$ gives $1/(2q)$ and proves the global statements. At $b=1$, $f=\frac{1}{2}-\frac{1}{2}(r^2-1{)}^2$, so the entire unit circle is the maximum set.

A Rankine vortex models a finite rotating core joined continuously to an irrotational exterior. In cylindrical polar coordinates, its radial and vertical velocities vanish and its azimuthal velocity has the displayed form, where $a$ is the core radius and $\Omega$ is a signed angular velocity. The axial vorticity is $r^{-1}{\partial }_r(r{u}_{\theta })$, equal to $2\Omega$ inside and zero outside. The velocity is continuous at $r=a$, although its first radial derivative has different one-sided values. The circulation outside the core is $2\pi \Omega a^2$.