TY - JOUR
T1 - High gradient phenomena in two-dimensional vortex interactions
AU - Yao, H. B.
AU - Zabusky, N. J.
AU - Dritschel, D. G.
PY - 1995/1/1
Y1 - 1995/1/1
N2 - Hyperdiffusion, a simple linear eddy diffusivity scheme, is commonly used in atmospheric and oceanic simulations because it increases the range of inertially behaving spatial scales for a given model resolution. Compared with molecular diffusion (which is utterly negligible in the atmosphere and oceans), hyperdiffusion more sharply confines the dissipation to the smallest scales of the numerical model. But is this all that hyperdiffusion does? In this paper, the inelastic interaction of two distributed vortices of unequal size is examined. Contour surgery (CS) simulations are compared with pseudospectral (PS) simulations employing hyperdiffusion or molecular diffusion. The example illustrates what is believed to be the most fundamental characteristic of two-dimensional (2-D) (and layerwise-2-D) vortex dynamics, namely, the formation of exceedingly high vorticity gradients. There is an excellent agreement between the hyperdiffusive PS and CS calculations at early times (i.e., for a few vortex rotation periods). Thereafter, significant discrepancies develop, beginning abruptly from the time when vorticity-gradient intensification is arrested by diffusion. A rapid inward erosion of the smaller of the two vortices then takes place. This erosion takes place under the joint action of (hyper) diffusion and stripping (the peeling of the vortex periphery by the external flow). With hyperdiffusion, the erosion is accompanied by a serious numerical artifact: a climb in the peak vorticity by 30% in this example. Eventually, the erosion reaches the vortex center and the vortex is sheared into a filament. In the CS calculation, there is no erosion, no climb in peak vorticity, and the vortex appears to last indefinitely. In the PS calculations, the viscosity or hyperdiffusion is adjusted according to the resolution to give the largest possible inertial range while ensuring numerical stability. It is found that vortices that are spanned by fewer than 10-20 grid points are eroded away in only a few vortex rotation periods (a time scale that is very much shorter than one would estimate from pure viscous decay). These findings bring into question the results of many 2-D turbulence simulations using hyperdiffusion, for hyperdiffusion simulates neither inviscid dynamics nor molecular-diffusive dynamics.
AB - Hyperdiffusion, a simple linear eddy diffusivity scheme, is commonly used in atmospheric and oceanic simulations because it increases the range of inertially behaving spatial scales for a given model resolution. Compared with molecular diffusion (which is utterly negligible in the atmosphere and oceans), hyperdiffusion more sharply confines the dissipation to the smallest scales of the numerical model. But is this all that hyperdiffusion does? In this paper, the inelastic interaction of two distributed vortices of unequal size is examined. Contour surgery (CS) simulations are compared with pseudospectral (PS) simulations employing hyperdiffusion or molecular diffusion. The example illustrates what is believed to be the most fundamental characteristic of two-dimensional (2-D) (and layerwise-2-D) vortex dynamics, namely, the formation of exceedingly high vorticity gradients. There is an excellent agreement between the hyperdiffusive PS and CS calculations at early times (i.e., for a few vortex rotation periods). Thereafter, significant discrepancies develop, beginning abruptly from the time when vorticity-gradient intensification is arrested by diffusion. A rapid inward erosion of the smaller of the two vortices then takes place. This erosion takes place under the joint action of (hyper) diffusion and stripping (the peeling of the vortex periphery by the external flow). With hyperdiffusion, the erosion is accompanied by a serious numerical artifact: a climb in the peak vorticity by 30% in this example. Eventually, the erosion reaches the vortex center and the vortex is sheared into a filament. In the CS calculation, there is no erosion, no climb in peak vorticity, and the vortex appears to last indefinitely. In the PS calculations, the viscosity or hyperdiffusion is adjusted according to the resolution to give the largest possible inertial range while ensuring numerical stability. It is found that vortices that are spanned by fewer than 10-20 grid points are eroded away in only a few vortex rotation periods (a time scale that is very much shorter than one would estimate from pure viscous decay). These findings bring into question the results of many 2-D turbulence simulations using hyperdiffusion, for hyperdiffusion simulates neither inviscid dynamics nor molecular-diffusive dynamics.
UR - http://www.scopus.com/inward/record.url?scp=0028976499&partnerID=8YFLogxK
U2 - 10.1063/1.868650
DO - 10.1063/1.868650
M3 - Article
AN - SCOPUS:0028976499
SN - 1070-6631
VL - 7
SP - 539
EP - 548
JO - Physics of Fluids
JF - Physics of Fluids
IS - 3
ER -