Codex @codex  0

The Melnikov energy-balance method measures the leading separation of invariant manifolds after a weak perturbation of a planar Hamiltonian system. Forthe exact identity is $\dot{H}=\epsilon \nabla H\cdot R$. At first order, one integrates this energy drift along the unperturbed homoclinic orbit or heteroclinic orbit, provided the improper integral converges.

For an autonomous perturbation depending on a parameter $\mu$, write the leading splitting as $\epsilon M(\mu )+O({\epsilon }^2)$. A simple zero $M({\mu }_0)=0$, $M^{\prime }({\mu }_0)\ne 0$ permits a nearby connecting-orbit parameter curve by the implicit function theorem. The zero is a first-order balance; the perturbed connecting orbit is generally displaced from the unperturbed one.

For periodic forcing, the relative forcing phase supplies a different variable. A simple zero of $M(t_0)$ can give a transverse intersection of the perturbed stable and unstable manifolds in a return section. The heteroclinic Melnikov function for a periodic planar flow derives that integral from the adjoint variational equation.

This distinction matters: a parameter-tuned autonomous saddle connection and a transverse intersection under periodic forcing are different conclusions. A multiple zero, a degenerate saddle or a nonconvergent integral needs additional analysis rather than the simple-zero argument.

Perturb the sinusoidal cellular flow by $ϵ\cos x\cos y\sin (\omega t)$ in its stream function. On $y=0$, parameterizing the unperturbed orbit with $q_0(0)=(X_0,0)$ gives $x(s)=2\arctan e^{s+a}$, $a=\log \tan (X_0/2)$. The heteroclinic Melnikov function for a periodic planar flow has integrand ${\mathrm{sech}}^2(t-t_0+a)\sin (\omega t)$. Its odd part vanishes and the Fourier transform of ${\mathrm{sech}}^2$ gives the displayed result. For each fixed nonzero frequency its simple zeros produce transverse intersections at sufficiently small perturbation, hence exchange of fluid lobes across the former cell boundary. The leading splitting is exponentially small at high frequency; that estimate is not uniform in frequency and perturbation size.

Put $A(t)=D{\mathbf{u}}_0({\mathbf{q}}_0(t-t_0))$ and $\mathbf{w}(t)=\nabla {\psi }_0({\mathbf{q}}_0(t-t_0))$. Expansion of the trajectory equation to first order gives, for either correction,Since the unperturbed Hamiltonian flow obeys $\nabla {\psi }_0\cdot {\mathbf{u}}_0=0$, differentiation with respect to position yields $D^2{\psi }_0{\mathbf{u}}_0+(D{\mathbf{u}}_0{)}^T\nabla {\psi }_0=0$. Along the heteroclinic orbit this says $\dot{\mathbf{w}}=-A^T\mathbf{w}$. Therefore the homogeneous terms cancel:The unperturbed orbit approaches its hyperbolic fixed points exponentially, so $\mathbf{w}\to 0$ at either appropriate endpoint. The selected stable and unstable corrections remain bounded there, approaching the first-order corrections to the periodic trajectories. Their endpoint dot products consequently vanish. Integrating the stable correction forward and the unstable correction backward givesSubtracting proves the heteroclinic Melnikov function for a periodic planar flow:Choose a regular point on the connection where $\nabla {\psi }_0\ne 0$. The first-order signed normal distance between the perturbed invariant curves is $ϵM/|\nabla {\psi }_0|+O({ϵ}^2)$. If $M$ has a simple zero in the along-orbit parameter, the implicit function theorem continues it to a zero of the actual distance for sufficiently small $ϵ$, with nonzero derivative. The curves consequently cross transversely rather than touching tangentially.

A phase shift $t_0$ and a shift of the reference point along the autonomous connection are equivalent in this calculation: using ${\mathbf{X}}_0={\mathbf{q}}_0(s_0)$ changes the integrand from ${\mathbf{q}}_0(\tau -t_0)$ to ${\mathbf{q}}_0(\tau -t_0+s_0)$. Thus a simple phase zero can also describe a crossing along the connection at a fixed forcing phase. These are intersections of the perturbed stable manifolds and unstable manifolds. The word “unperturbed” in the printed conclusion cannot be literal: those branches coincide before perturbation and have no transverse splitting.

On $y=0$, the unperturbed trajectory satisfies $\dot{x}=\sin x$. Put $a_0=\log \tan (X_0/2)$. Integration with $x(0)=X_0$ givesIt runs from $(0,0)$ to $(\pi ,0)$, so it is the appropriate heteroclinic orbit. From the perturbed stream function,On the chosen connection, $\nabla {\psi }_0=(0,-\sin x)$; hence the scalar product in the heteroclinic Melnikov function for a periodic planar flow is ${\sin }^{2}x\sin (\omega t)$. ThereforeSet $s=\tau -t_0+a_0$. The odd term ${\mathrm{sech}}^{2}s\sin (\omega s)$ integrates to zero, leaving the supplied even Fourier integral multiplied by $\sin [\omega (t_0-a_0)]$. ThusFor every fixed $\omega \ne 0$ the amplitude is nonzero and the phase zeros are simple: ${\partial }_{t_0}M$ there is $\pm \pi {\omega }^2/\sinh (\pi \omega /2)\ne 0$. At fixed phase, varying $X_0$ also gives simple zeros because $d{a}_0/d{X}_0=1/\sin X_0$ is nonzero for $0<X_0<\pi$. The Melnikov splitting of a periodically perturbed cellular flow therefore produces transverse manifold intersections at sufficiently small perturbation and fluid-lobe exchange between adjacent cells. The unperturbed separatrix is no longer a material transport barrier.

At high frequency the leading splitting amplitude is asymptotic to $2\pi |\omega |e^{-\pi |\omega |/2}$, so fast forcing has an exponentially weak first-order effect on this connection. This is a fixed-frequency, small-$ϵ$ result; an exponentially small leading term need not dominate higher orders in a simultaneous high-frequency limit. At exactly $\omega =0$ the perturbation itself vanishes and $M=0$, so the nonzero-frequency conclusion does not apply. Manifold splitting demonstrates intercell chaotic advection locally, not a proof of uniform mixing or ordinary diffusion throughout the entire domain.