首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
We develop an approach to deriving the three-dimensional non-force-free coronal magnetic field from vector magnetograms. Based on the principle of minimum dissipation rate, a general non-force-free magnetic field is expressed as the superposition of one potential field and two constant-α (linear) force-free fields. Each is extrapolated from its bottom boundary data, providing the normal component only. The constant-α parameters are distinct and determined by minimizing the deviations between the numerically computed and measured transverse magnetic field at the bottom boundary. The boundary conditions required are at least two layers of vector magnetograms, one at the photospheric level and the other at the chromospheric level, presumably. We apply our approach to a few analytic test cases, especially to two nonlinear force-free cases examined by Schrijver et al. (Solar Phys. 235, 161, 2006). We find that for one case with small α parameters, the quantitative measures of the quality of our result are better than the median values of those from a set of nonlinear force-free methods. The reconstructed magnetic-field configuration is valid up to a vertical height of the transverse scale. For the other cases, the results remain valid to a lower vertical height owing to the limitations of the linear force-free-field solver. Because our method is based on the fast-Fourier-transform algorithm, it is much faster and easy to implement. We discuss the potential usefulness of our method and its limitations.  相似文献   

2.
The announced missions to the Saturn and Jupiter systems renewed the space community interest in simple design methods for gravity assist tours at planetary moons. A key element in such trajectories are the V-Infinity Leveraging Transfers (VILT) which link simple impulsive maneuvers with two consecutive gravity assists at the same moon. VILTs typically include a tangent impulsive maneuver close to an apse location, yielding to a desired change in the excess velocity relative to the moon. In this paper we study the VILT solution space and derive a linear approximation which greatly simplifies the computation of the transfers, and is amenable to broad global searches. Using this approximation, Tisserand graphs, and heuristic optimization procedure we introduce a fast design method for multiple-VILT tours. We use this method to design a trajectory from a highly eccentric orbit around Saturn to a 200-km science orbit at Enceladus. The trajectory is then recomputed removing the linear approximation, showing a Δv change of <4%. The trajectory is 2.7 years long and comprises 52 gravity assists at Titan, Rhea, Dione, Tethys, and Enceladus, and several deterministic maneuvers. Total Δv is only 445 m/s, including the Enceladus orbit insertion, almost 10 times better then the 3.9 km/s of the Enceladus orbit insertion from the Titan–Enceladus Hohmann transfer. The new method and demonstrated results enable a new class of missions that tour and ultimately orbit small mass moons. Such missions were previously considered infeasible due to flight time and Δv constraints.  相似文献   

3.
Wiegelmann  T. 《Solar physics》2004,219(1):87-108
We developed a code for the reconstruction of nonlinear force-free and non-force-free coronal magnetic fields. The 3D magnetic field is computed numerically with the help of an optimization principle. The force-free and non-force-free codes are compiled in one program. The force-free approach needs photospheric vector magnetograms as input. The non-force-free code additionally requires the line-of-sight integrated coronal density distribution in combination with a tomographic inversion code. Previously the optimization approach has been used to compute magnetic fields using all six boundaries of a computational box. Here we extend this method and show how the coronal magnetic field can be reconstructed only from the bottom boundary, where the boundary conditions are measured with vector magnetographs. The program is planed for use within the Stereo mission.  相似文献   

4.
We construct a simple, regularized estimator for the dark energy equation of state by using the recently introduced linear response approximation. We show that even a simple regularization substantially improves the performance of the free-form fitting approach. The use of the linear response approximation allows an analytical construction of the maximum likelihood estimator, in a convenient and easy to use matrix form. We show that, in principle, such regularized free-form fitting can give us an unbiased estimate of the functional form of the equation of state of dark energy. We show the efficacy of this approach on simulated SuperNova Acceleration Probe class data, but it is easy to generalize this method to include other cosmological tests. We provide a possible explanation for the sweet spots seen in other reconstruction methods.  相似文献   

5.
The effect of meridional circulation on spherical shell dynamos is considered in the Parker approximation. We demonstrate that the type of the exited solution crucially depend on the intensity of the meridional circulation. If the circulation is equatorward or if it is polarward however do not exceed some critical value, an oscillating solution in form of an equatorward traveling wave is excited. If the polarward meridional circulation becomes too intensive the solution becomes steady growing. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)  相似文献   

6.
We model the propagation of a coronal shock wave, using nonlinear geometrical acoustics. The method is based on the Wentzel–Kramers–Brillouin (WKB) approach and takes into account the main properties of nonlinear waves: i) dependence of the wave front velocity on the wave amplitude, ii) nonlinear dissipation of the wave energy, and iii) progressive increase in the duration of solitary shock waves. We address the method in detail and present results of the modeling of the propagation of shock-associated extreme-ultraviolet (EUV) waves as well as Moreton waves along the solar surface in the simplest solar corona model. The calculations reveal deceleration and lengthening of the waves. In contrast, waves considered in the linear approximation keep their length unchanged and slightly accelerate.  相似文献   

7.
The minimum dissipative rate (MDR) method for deriving a coronal non-force-free magnetic field solution is partially evaluated. These magnetic field solutions employ a combination of three linear (constant-α) force-free-field solutions with one being a potential field (i.e., α=0). The particular case of the solutions where the other two α’s are of equal magnitude but of opposite sign is examined. This is motivated by studying the SOLIS (Synoptic Optical Long-term Investigation of the Sun (SOLIS), a National Solar Observatory facility) vector magnetograms of AR 10987, which show a global α value consistent with an α=0 value as evaluated by (∇×B) z /B z over the region. Typical of the current state of the observing technology, there is no definitive twist for input into the general MDR method. This suggests that the special α case, of two α’s with equal magnitudes and opposite signs, is appropriate given the data. Only for an extensively twisted active region does a dominant, nonzero α normally emerge from a distribution of local values. For a special set of conditions, is it found that (i) the resulting magnetic field is a vertically inflated magnetic field resulting from the electric currents being parallel to the photosphere, similar to the results of Gary and Alexander (Solar Phys. 186:123, 1999), and (ii) for α≈(α max /2), the Lorentz force per unit volume normalized by the square of the magnetic field is on the order of 1.4×10−10 cm−1. The Lorentz force (F L) is a factor of ten higher than that of the magnetic force d(B 2/8π)/dz, a component of F L. The calculated photospheric electric current densities are an order of magnitude smaller than the maximum observed in all active regions. Hence both the Lorentz force density and the generated electric current density seem to be physically consistent with possible solar dynamics. The results imply that the field could be inflated with an overpressure along the neutral line. However, the implementation of this or any other extrapolation method using the electric current density as a lower boundary condition must be done cautiously, with the current magnetography.  相似文献   

8.
Zhao and Kosovichev (Astrophys. J. 591, 446, 2003) found two opposite sub-photospheric vortical flows in the depth range of 0 – 12 Mm around a fast rotating sunspot. So far there is no theoretical model explaining such flow motions. In this paper, we try to explain this phenomenon from the point of view of magnetic flux tubes interacting with large-scale vortical motions of plasma. In the deeper zone under the photosphere, the magnetic force may be less than the nonmagnetic force of plasma. The vortical flow located there twists the flux tube and magnetic free energy is built up in the tube. In the shallower zone under the photosphere, the magnetic force may be greater than the nonmagnetic force. Thus, part of the stored magnetic free energy is released to drive the plasma to rotate in two opposite directions, e.g., in the depth ranges of 0 – 3(5) and 9 – 12 Mm. In addition, we also define a vector of nonpotential magnetic stress τ, which can be related to flare occurrence. It is calculated for the active region NOAA 10930 on 11 December 2006. We find that: i) the integral of its line-of-sight (LOS) stress successively increases around the magnetic neutral line (MNL) prior to and during the flare and decreases to a minimum after the flare; ii) the integral of its transverse stress exceeds the integral of its LOS component by one order of magnitude over the whole field of view; iii) the transverse stress first points toward the MNL, then along it, and finally it points away from it. We need other data to verify whether or not the magnetic energy is transported in the horizontal direction to the neutral line, and then partly changes into the energy in LOS direction before and during the flare.  相似文献   

9.
The problem of magnetic field generation under screw motion in a toroidal channel is studied numerically. The screw dynamo in the cylinder with periodical boundary conditions was found to be a suitable approximation for generation of the magnetic field by a screw flow in a thin torus. For the thick torus, a principally new solution of the screw dynamo problem was obtained. In this case the growing global magnetic field mode has the scale of a maximal geometrical size of the torus and does not vanish on the axis of the torus (© 2011 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)  相似文献   

10.
Neutral Line associated Sources (NLSs) are quasi-stationary microwave sources projected onto vicinities of the neutral line of the photospheric magnetic field. NLSs are often precursors of powerful flares, but their nature is unclear. We endeavor to reveal the structure of an NLS and to analyze a physical connection between such a source with a site of energy release in the corona above NOAA 10488 (October/November 2003). Evolution of this AR includes emergence and collision of two bipolar magnetic structures, rise of the main magnetic separator, and the appearance of an NLS underneath. The NLS appears at a contact site of colliding sunspots, whose relative motion goes on, resulting in a large shear along a tangent. Then the nascent NLS becomes the main source of microwave fluctuations in the AR. The NLS emission at 17 GHz is dominated by either footpoints or the top of a loop-like structure, an NLS loop, which connects two colliding sunspots. During a considerable amount of time, the emission dominates over that footpoint of the NLS loop, where the magnetic field is stronger. At that time, the NLS resembles a usual sunspot–associated radio source, whose brightness center is displaced towards the periphery of a sunspot. Microwave emission of an X2.7 flare is mainly concentrated in an ascending flare loop, initially coinciding with the NLS loop. The top of this loop is located at the base of a non-uniform bar-like structure visible in soft X-rays and at 34 GHz at the flare onset. We reveal i) upward lengthening of this bar before the flare onset, ii) the motion of the top of an apparently ascending flare loop along the axis of this bar, and iii) a non-thermal microwave source, whose descent along the bar was associated with the launching of a coronal ejection. We connect the bar with a probable position of a nearly vertical diffusion region, a site of maximal energy release inside an extended pre-flare current sheet. The top of the NLS loop is located at the bottom of this region. A combination of the NLS loop and diffusion region constitutes the skeleton of a quasi-stationary microwave NLS.  相似文献   

11.
The objective of this paper is to find periodic solutions of the circular Sitnikov problem by the multiple scales method which is used to remove the secular terms and find the periodic approximated solutions in closed forms. Comparisons among a numerical solution (NS), the first approximated solution (FA) and the second approximated solution (SA) via multiple scales method are investigated graphically under different initial conditions. We observe that the initial conditions play a vital role in the numerical and approximated solutions behaviour. The obtained motion is periodic, but the difference of its amplitude is directly proportional with the initial conditions. We prove that the obtained motion by the numerical or the second approximated solutions is a regular and periodic, when the infinitesimal body starts its motion from a nearer position to the common center of primaries. Otherwise when the start point distance of motion is far from this center, the numerical solution may not be represent a periodic motion for along time, while the second approximated solution may present a chaotic motion, however it is always periodic all time. But the obtained motion by the first approximated solution is periodic and has regularity in its periodicity all time. Finally we remark that the provided solutions by multiple scales methods reflect the true motion of the Sitnikov restricted three–body problem, and the second approximation has more accuracy than the first approximation. Moreover the solutions of multiple scales technique are more realistic than the numerical solution because there is always a warranty that the motion is periodic all time.  相似文献   

12.
We study Parker instability (PI) operating in a non-adiabatic, gravitationally stratified, interstellar medium. We discuss models with two kinds of heating mechanisms. The first one results from photoionization models. The other, relying on supplemental sources, has been postulated by Reynolds, Haffner & Tufte. The cooling rate, corresponding to radiative losses, is an approximation to the one given by Dalgarno & McCray. An unperturbed state of the system represents a magnetohydrostatic and thermal equilibrium. We perform linear stability analysis by solving the boundary value problem. We find that the maximum growth rate of PI rises for increasing magnitudes of non-adiabatic effects. In the pure photoionization model, the maximum growth rate of the general non-adiabatic case coincides with the isothermal limit. Adding other sources of heat leads to a maximum growth rate that is larger than the one corresponding to the isothermal limit. We find that the influence of the supplemental heating on PI also leads to a decrease in temperature in magnetic valleys. Finally, we conclude that the initial gas cooling due to the action of PI may promote a subsequent onset of thermal instability in magnetic valleys and formation of giant molecular clouds.  相似文献   

13.
The radio approach based on the Askaryan effect for detecting the ultra-high energy cosmic neutrinos has become a mature experimental technique. So far the existing calculations of the Cherenkov radiation associated with the Askaryan effect has been mostly based on the far-field approximation, whose validity maybe challenged when the detector is close to the event. In this paper we present an alternative approach to calculate the Cherenkov pulse by a numerical code based on the finite difference time-domain (FDTD) method. This approach has the advantage of providing the solution everywhere in space, contrary to other methods that rely on the far-field approximation. We also present a one-dimensional theoretical model for the shower with analytical solution, which helps to elucidate our nonzero-width simulation results. We show that for a shower with symmetric longitudinal development, the resulting near-field waveform would be asymmetric in time. In addition, we demonstrate that for a shower elongated by the LPM (Landau-Pomeranchuk-Migdal) effect and thus with a multi-peak structure, a bipolar, asymmetric waveform is still preserved in the near-field regime irrespective of the specific variations of the multi-peak structure, which makes it a generic, distinctive feature. This should provide an important characteristic signature for the identification of ultra-high energy cosmogenic neutrinos.  相似文献   

14.
We present an interferometric technique for the reconstruction of ultra-wide band impulsive signals from point sources. This highly sensitive method was developed for the search for ultra-high energy neutrinos with the ANITA experiment but is fully generalizable to any antenna array detecting radio impulsive events. Applications of the interferometric method include event reconstruction, thermal noise and anthropogenic background rejection, and solar imaging for calibrations. We illustrate this technique with applications from the analysis of the ANITA-I and ANITA-II data in the 200–1200 MHz band.  相似文献   

15.
The properties of solar magnetic fields on scales less than the spatial resolution of solar telescopes are studied. A synthetic infrared spectropolarimetric diagnostic based on a 2D MHD simulation of magnetoconvection is used for this. Analyzed are two time sequences of snapshots that likely represent two regions of the network fields with their immediate surroundings on the solar surface with unsigned magnetic flux densities of 300 and 140 G. In the first region from the probability density functions of the magnetic field strength it is found that the most probable field strength at log τ 5=0 is equal to 250 G. Weak fields (B<500 G) occupy about 70% of the surface, whereas stronger fields (B>1000 G) occupy only 9.7% of the surface. The magnetic flux is −28 G and its imbalance is −0.04. In the second region, these parameters are correspondingly equal to 150 G, 93.3%, 0.3%, −40 G, and −0.10. The distribution of line-of-sight velocities on the surface of log τ 5=−1 is estimated. The mean velocity is equal to 0.4 km s−1 in the first simulated region. The average velocity in the granules is −1.2 km s−1 and in the intergranules it is 2.5 km s−1. In the second region, the corresponding values of the mean velocities are equal to 0, −1.8, and 1.5 km s−1. In addition the asymmetry of synthetic Stokes V profiles of the Fe i 1564.8 nm line is analyzed. The mean values of the amplitude and area asymmetry do not exceed 1%. The spatially smoothed amplitude asymmetry is increased to 10% whereas the area asymmetry is only slightly varied.  相似文献   

16.
We present a novel numerical method that allows the calculation of nonlinear force-free magnetostatic solutions above a boundary surface on which only the distribution of the normal magnetic field component is given. The method relies on the theory of force-free electrodynamics and applies directly to the reconstruction of the solar coronal magnetic field for a given distribution of the photospheric radial field component. The method works as follows: we start with any initial magnetostatic global field configuration (e.g. zero, dipole), and along the boundary surface we create an evolving distribution of tangential (horizontal) electric fields that, via Faraday’s equation, give rise to a respective normal-field distribution approaching asymptotically the target distribution. At the same time, these electric fields are used as boundary condition to numerically evolve the resulting electromagnetic field above the boundary surface, modeled as a thin ideal plasma with non-reflecting, perfectly absorbing outer boundaries. The simulation relaxes to a nonlinear force-free configuration that satisfies the given normal-field distribution on the boundary. This is different from existing methods relying on a fixed boundary condition – the boundary evolves toward the a priori given one, at the same time evolving the three-dimensional field solution above it. Moreover, this is the first time that a nonlinear force-free solution is reached by using only the normal field component on the boundary. This solution is not unique, but it depends on the initial magnetic field configuration and on the evolutionary course along the boundary surface. To our knowledge, this is the first time that the formalism of force-free electrodynamics, used very successfully in other astrophysical contexts, is applied to the global solar magnetic field.  相似文献   

17.
Improvements to an existing method for calculating nonlinear force-free magnetic fields (Wheatland, Solar Phys. 238, 29, 2006) are described. In particular a solution of the 3-D Poisson equation using 2-D Fourier transforms is presented. The improved nonlinear force-free method is demonstrated in application to linear force-free test cases with localized nonzero values of the normal component of the field in the boundary. These fields provide suitable test cases for nonlinear force-free calculations because the boundary conditions involve localized nonzero values of the normal components of the field and of the current density, and because (being linear force-free fields) they have more direct numerical solutions. Despite their simplicity, fields of this kind have not been recognized as test cases for nonlinear methods before. The examples illustrate the treatment of the boundary conditions on current in the nonlinear force-free method, and in particular the limitations imposed by field lines that connect outside of the boundary region.  相似文献   

18.
Asteroids and comets 10–100 m in size that collide with Earth disrupt dramatically in the atmosphere with an explosive transfer of energy, caused by extreme air drag. Such airbursts produce a strong blastwave that radiates from the meteoroid's trajectory and can cause damage on the surface. An established technique for predicting airburst blastwave damage is to treat the airburst as a static source of energy and to extrapolate empirical results of nuclear explosion tests using an energy‐based scaling approach. Here we compare this approach to two more complex models using the iSALE shock physics code. We consider a moving‐source airburst model where the meteoroid's energy is partitioned as two‐thirds internal energy and one‐third kinetic energy at the burst altitude, and a model in which energy is deposited into the atmosphere along the meteoroid's trajectory based on the pancake model of meteoroid disruption. To justify use of the pancake model, we show that it provides a good fit to the inferred energy release of the 2013 Chelyabinsk fireball. Predicted overpressures from all three models are broadly consistent at radial distances from ground zero that exceed three times the burst height. At smaller radial distances, the moving‐source model predicts overpressures two times greater than the static‐source model, whereas the cylindrical line‐source model based on the pancake model predicts overpressures two times lower than the static‐source model. Given other uncertainties associated with airblast damage predictions, the static‐source approach provides an adequate approximation of the azimuthally averaged airblast for probabilistic hazard assessment.  相似文献   

19.
The polarization-free (POF) approximation (Trujillo Bueno and Landi Degl'Innocenti, 1996) is capable of accounting for the approximate influence of the magnetic field on the statistical equilibrium, without actually solving the full Stokes vector radiative transfer equation. The method introduces the Zeeman splitting or broadening of the line absorption profile I in the scalar radiative transfer equation, but the coupling between Stokes I and the other Stokes parameters is neglected. The expected influence of the magnetic field is largest for strongly-split strong lines and the effect is greatly enhanced by gradients in the magnetic field strength. Formally the interaction with the other Stokes parameters may not be neglected for strongly-split strong lines, but it turns out that the error in Stokes I obtained through the POF approximation to a large extent cancels the neglect of interaction with the other Stokes parameters, so that the resulting line source functions and line opacities are more accurate than those obtained with the field-free approach. Although its merits have so far only been tested for a two-level atom, we apply the POF approximation to multi-level non-LTE radiative transfer problems on the premise that there is no essential difference between these two cases. Final verification of its validity in multi-level cases still awaits the completion of a non-LTE Stokes vector transfer code.For two realistic multi-level cases (CaII and MgI in the solar atmosphere) it is demonstrated that the POF method leads to small changes, with respect to the field-free method, in the line source functions and emergent Stokes vector profiles (much smaller than for a two-level atom). Real atoms are dominated by strong ultraviolet lines (only weakly split) and continua, and most lines with large magnetic splitting (in the red and the infrared) are at higher excitation energies, i.e. they are relatively weak and unable to produce significant changes in the statistical equilibrium. We find that it is generally unpredictable by how much the POF results will differ from the field-free results, so that it is nearly always necessary to confirm predictions by actual computations.The POF approximation provides more reliable results than the field-free approximation without significantly complicating the radiative transfer problem, i.e. without solving any extra equations and without excessive computational resource requirements, so that it is to be preferred over the field-free approximation.  相似文献   

20.
We develop a radio astronomical approach to 3D‐reconstruction in few projections tomography. It is based on the 2‐CLEAN DSA method which consists of two clean algorithms by using a synthesized beam. In complex cases two extreme solutions are used for the analysis of the image structure. These solutions determine the limits of permissible energy redistribution on the image among the components of small and large scales. Two variants of 3D‐reconstruction proceeding from a set of two‐dimensional projections (3D2D) and from a set of one‐dimensional ones (3D1D) are considered. It is shown that the quality of 3D2D‐reconstruction should be similar to the quality of 2D1D‐reconstruction if the same number of equally spaced scans is used. But a doubled number of projections is required for 3D1D‐reconstruction. We have simulated 3D‐reconstruction of an optically thin emitting object. The present technique is a component of astrotomography and it has good prospects for a wide range of remote sensing. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号