半圆拱无力场的稳定性

本文从能量原理出发，导出了半圆拱无力场的能量积分的普遍形式，并证明所有半圆拱无力场都是稳定的。     

Nonlinear Force-Free Modeling of Coronal Magnetic Fields. II. Modeling a Filament Arcade and Simulated Chromospheric and Photospheric Vector Fields

We compare a variety of nonlinear force-free field (NLFFF) extrapolation algorithms, including optimization, magneto-frictional, and Grad – Rubin-like codes, applied to a solar-like reference model. The model used to test the algorithms includes realistic photospheric Lorentz forces and a complex field including a weakly twisted, right helical flux bundle. The codes were applied to both forced “photospheric” and more force-free “chromospheric” vector magnetic field boundary data derived from the model. When applied to the chromospheric boundary data, the codes are able to recover the presence of the flux bundle and the field’s free energy, though some details of the field connectivity are lost. When the codes are applied to the forced photospheric boundary data, the reference model field is not well recovered, indicating that the combination of Lorentz forces and small spatial scale structure at the photosphere severely impact the extrapolation of the field. Preprocessing of the forced photospheric boundary does improve the extrapolations considerably for the layers above the chromosphere, but the extrapolations are sensitive to the details of the numerical codes and neither the field connectivity nor the free magnetic energy in the full volume are well recovered. The magnetic virial theorem gives a rapid measure of the total magnetic energy without extrapolation though, like the NLFFF codes, it is sensitive to the Lorentz forces in the coronal volume. Both the magnetic virial theorem and the Wiegelmann extrapolation, when applied to the preprocessed photospheric boundary, give a magnetic energy which is nearly equivalent to the value derived from the chromospheric boundary, but both underestimate the free energy above the photosphere by at least a factor of two. We discuss the interpretation of the preprocessed field in this context. When applying the NLFFF codes to solar data, the problems associated with Lorentz forces present in the low solar atmosphere must be recognized: the various codes will not necessarily converge to the correct, or even the same, solution. On 07/07/2007, the NLFFF team was saddened by the news that Tom Metcalf had died as the result of an accident. We remain grateful for having had the opportunity to benefit from his unwavering dedication to the problems encountered in attempting to understand the Sun’s magnetic field; Tom had completed this paper several months before his death, leading the team through the many steps described above.     

Energy Buildup,Flux Confinement and Helicity Accumulation in the Solar Corona

Starting from a dipole field and a given distribution of footpoint displace- ment of field lines on the photosphere,we find axisymmetric,force-free field solutions in spherical coordinates that have the same distribution of normal field on the photo- sphere and magnetic topology as the dipole field.A photospheric shear is introduced in the azimuthal direction in a region that strides across the equator and ends at latitude γ_s.The footpoint displacement has a sine distribution in latitude and a peak amplitude of (?)_m.The magnetic energy E,azimuthal flux F_(?),and magnetic helicity H_T in the solar corona are then calculated for each force-free field solution.It is found that for a given shear region range γ_s,all of the three quantities increase monotoni- cally with increasing (?)_m.In particular,both F_(?) and H_T have a linear dependence on (?)m.When (?)_m reaches a certain critical value (?)_(mc),the force-free field loses equilib- rium,leading to a partial opening of the field and the appearance of a current sheet in the equatorial plane.At this point,E,F_(?)and H_T reach their maximum values, E_c,F_((?)c) and H_(Tc).E_c increases,and F_((?)c) and H_(Tc)decrease with decreasing λ_s.It is found that E_c is always smaller than the open field energy,in agreement with the Aly conjecture.Of the three critical parameters,E_c has the weakest dependence on λ_s.Therefore,if one is interested in the transition of a magnetic configuration from a stable state to a dynamic one,the magnetic energy is probably the most appropriate marker of the transition.     

A depth-dependent study of the topographic rectification of tidal currents

Abstract

A depth-dependent model for the topographic rectification of tidal currents in a homogeneous rotating fluid is used to examine the dependence of the rectified mean flow on various tidal, topographic and frictional parameters. Friction is parameterized through a vertically-uniform, time-independent vertical eddy viscosity and a bottom stress law applied near the top of the constant stress layer. The model neglects the interaction of mean and tidal currents, assumes uniformity along isobaths, and is closed with the assumption of zero depth-averaged mean flow across isobaths.

In the limit of depth-independence, the model reduces to that considered by Huthnance (1973) and Loder (1980) which, for weak friction, favours anticyclonic mean circulation around shallow regions and Lagrangian flow which is significantly reduced from the Eulerian. With the inclusion of vertical structure, the magnitude of the anticyclonic flow is amplified suggesting that depth-independent models may underestimate the along-isobath flow. For strong friction the direction of the mean flow depends on the orientation of the tidal ellipse relative to the isobaths. The depthindependent model again underestimates the magnitude of the along-isobath flow, but this can be offset with an appropriate reduction of the bottom friction coefficient.

The cross-isobath mean flows are one to two orders of magnitude weaker than the along-isobath flows and generally have more vertical structure. There is also a significant Stokes drift in the cross-isobath direction. Although there is some tendency for the cross-isobath mean bottom current to be down the cross-isobath mean pressure gradient, it appears that it is not generally possible to infer this current from depth-independent models.     

Laboratory experiments on fronts

Abstract

Supercritically unstable density fronts near a vertical wall in a rotating, two-layer fluid were created on a laboratory turntable by withdrawing the outer wall of an annulus with a narrow gap, and allowing buoyant fluid from within the annulus to collapse toward a state of quasi-geostrophic balance. The resulting “coastal” current has a nearly uniform potential vorticity and is bounded by a front on which ageostrophic, wave-like disturbances grow. If the current width is comparable to the Rossby radius of deformation, the dominant length scale of disturbances is proportional to the width of the current. On the other hand, if the upper layer is much wider than the Rossby radius, then the observed length scale is a constant multiple of the Rossby radius. If the vertical boundary is omitted in the experiments, so that we are left with a circular anticyclonic vortex, the observed length scales and large-amplitude behaviour of disturbances are identical to those for the boundary currents, indicating that the wall has no significant influence on the flow.

At very large amplitude the growing waves lead to the formation of cyclone-anticyclone vortex pairs. For very wide currents, both the mean flow and the disturbances are first confined to a region within a few Rossby radii of the front. However, both the mean flow and the turbulent eddy motions slowly propagate into the previously stationary upper layer until, eventually, the whole of the upper layer is turbulent.     

Tangential discontinuities and the optical analogy for stationary fields. v. formal integration of the force-free field equations

Abstract

This paper demonstrates the appearance of tangential discontinuities in deformed force-free fields by direct integration of the field equation ? x B = αB. To keep the mathematics tractable the initial field is chosen to be a layer of linear force-free field B_x = + B ₀cosqz, B_y = — B ₀sinqz, B_z = 0, anchored at the distant cylindrical surface ? = (x ² + y ²)^1/2 = R and deformed by application of a local pressure maximum of scale l centered on the origin x = y = 0. In the limit of large R/l the deformed field remains linear, with α = q[1 + O(l ²/R ²)]. The field equations can be integrated over ? = R showing a discontinuity extending along the lines of force crossing the pessure maximum. On the other hand, examination of the continuous solutions to the field equations shows that specification of the normal component on the enclosing boundary ? = R completely determines the connectivity throughout the region, in a form unlike the straight across connections of the initial field. The field can escape this restriction only by developing internal discontinuities.

Casting the field equation in a form that the connectivity can be specified explicitly, reduces the field equation to the eikonal equation, describing the optical analogy, treated in papers II and III of this series. This demonstrates the ubiquitous nature of the tangential discontinuity in a force-free field subject to any local deformation.     

Force-free magnetic fields with singular current density surfaces and their stability

Starting from Bernstein's principle of magnetohydrodynamic energy, a general analysis is presented for the stability of a kind of 1-D force-free magnetic fields with singular current density surfaces and a single parameter in cylindrical coordinates. It is found that in the parameter space of this kind of force-free magnetic fields there simultaneously exist stable and unstable regions. Their stability is solely determined by the radial distribution of the magnetic pitch in the neighborhood of the cylinder axis, and is independent of the presence of singular current density surface at the boundary of the field.     

日冕三重无力场电流片的磁场重联

采用二维三分量磁流体力学模型,对日冕三重无力场电流片的磁场重联进行了数值研究,揭示了重联过程的基本物理特征．这类重联过程将加热和加速日冕等离子体,并导致多个高温、高密度、高磁螺度的磁岛的形成和向上喷发．这表明,多重无力场电流片的重联可能在日冕磁能释放、上行等离子体团的形成和太阳磁场螺度向行星际空间的逃逸方面起重要的作用．     

Spontaneous tangential discontinuities and the optical analogy for static magnetic fields. VI. Topology of current sheets

Abstract

The electric surface current in a tangential discontinuity in a force-free magnetic field is conserved. The direction of the current is halfway between the direction of the continuous fields on either side of the surface of discontinuity. Hence the current sheets, i.e. the surface of tangential discontinuity, have a topology that is distinct from the lines of force of the field. The precise nature of the topology of the current sheet depends upon the form of the winding patterns in the field. Hence, invariant winding patterns and random winding patterns are treated separately. Current sheets may have edges, at the junction of two or more topological separatrices. The current lines may, in special cases, be closed on themselves. The lines of force that lie on either side of a current sheet somewhere pass off the sheet across a junction onto another sheet. In most cases the current sheets extending along a field make an irregular honeycomb.

The honeycomb pattern varies along the field if the winding pattern of the field varies. The surface current density in a tangential discontinuity declines inversely, or faster, with distance from its region of origin. The edges of weaker tangential discontinuities (originating in more distant regions) are bounded by the stronger tangential discontinuities (of nearby origin).

An examination of the force-free field equations in a small neighborhood of the line of intersection of two tangential discontinuities shows that the lines of force twist around to cross the line of intersection at right angles. If the angle between the tangential discontinuities exceeds π/2, there is also the possibilitity that the lines twist around so as to come tangent to the line of intersection as they cross it.