耀斑后环珥系速度场的定量研究

本文研究了在太阳重力、偶极磁压强梯度和气体压力梯度联合作用下环珥系内物质的下落运动,并在等温、定常假设下用数值方法计算了该环系的二维视向速度场。通过将速度场的理论计算跟文[1]中给出的观测速度场拟合,导出了该环系的有关参数及分布。计算表明环系内物质密度和磁场强度对物质下落运动影响较显著。     

Coronal heating mechanisms

As a result of the large body of data available from solar and stellar coronae, our understanding of the mechanisms responsible for the heating of coronal plasmas to temperatures of the order of ～ 10⁸ K has changed. The solar corona is highly structured by magnetic fields and the acoustic shocks which, according to early theories, should have acted as the coronal energy source have not been observed. Einstein Observatory data show moreover that coronae are present in most regions of the H-R diagram. The observed relationship between X-ray luminosity and rotational velocity in dwarf stars from spectral types F to M again suggests an active role for the magnetic fields.The basic picture which is emerging is that coronae in stellar types from F to M are produced because of the interaction of the magnetic field with the convective velocity fields generated in the photosphere resulting in MHD waves or currents which dissipate in the corona. X-ray emission in early type stars cannot be explained with this mechanism and the models which have been proposed for these stars are not yet completely satisfactory.     

日冕冲浪形成的磁流体动力学模拟

应用二维时变可压缩磁流体动力学模拟,数值研究了双极-单极磁场中电阻撕裂模不稳定性引起的磁场重联过程,用于模拟日冕冲浪的形成.结果表明,在包含有三区——双极场、电流片和单极场的磁静力平衡初态下,双极场和单极场中的磁力线将会直接重联,磁场演变成鞭状(whip)结构.由弯曲磁力线支撑的等离子体团向上运动到最高位置后,逐渐下落和弥散.等离子体团上升速度可达到0.10v_A(v_A为双极场中的Alfv'én速度).模拟结果证实日冕冲浪的形成可能与双极-单极场中的磁场重联密切相关.      

Radiation from relativistic shocks in turbulent magnetic fields

Using our new 3-D relativistic particle-in-cell (PIC) code parallelized with MPI, we investigated long-term particle acceleration associated with a relativistic electron–positron jet propagating in an unmagnetized ambient electron–positron plasma. The simulations were performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. Acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value as predicted by hydrodynamic shock compression. In the jet (reverse) shock behind the bow (forward) shock the strongest electromagnetic fields are generated. These fields may lead to time dependent afterglow emission. In order to calculate radiation from first principles that goes beyond the standard synchrotron model used in astrophysical objects we have used PIC simulations. Initially we calculated radiation from electrons propagating in a uniform parallel magnetic field to verify the technique. We then used the technique to calculate emission from electrons in a small simulation system. From these simulations we obtained spectra which are consistent with those generated from electrons propagating in turbulent magnetic fields with red noise. This turbulent magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel instability. A fully developed shock within a larger simulation system may generate a jitter/synchrotron spectrum.     

2011年8月4日螺旋状日冕物质抛射爆发事件

利用多卫星多波段的综合观测数据,通过追踪光球表面等离子体速度分析计算了耀斑爆发前后磁螺度的变化,发现耀斑爆发前活动区中光球表面存在强的水平剪切运动,活动区磁螺度的注入主要由这种剪切运动所产生;使用CESE-MHD-NLFFF重建了耀斑爆发前后活动区的磁场位形,推测出耀斑过程中存在磁绳结构的抛射.基于这些分析,给出了这一螺旋状抛射结构的形成机制:爆发前暗条西侧足点的持续剪切运动驱动磁通量绳增加扭转,高度扭缠的通量绳与东侧足点附近的开放磁力线重联并与东侧足点断开,进而向外抛出并伴随解螺旋运动.另外,利用1AU处WIND卫星的观测数据在对应的行星际日冕物质抛射中找到典型磁云的观测特征.这表明除了传统上双足点均在太阳表面的磁云模型,这种单足点固定于太阳表面的磁通量绳爆发图景同样可能在行星系际空间形成磁云结构.研究结果对进一步认识磁云结构具有重要意义.      

Relation between solar wind velocity and properties of its source region

Different kinds of coronal holes are sources of different kind of solar winds. A successful solar wind acceleration model should be able to explain all those solar winds. For the modeling it is important to find a universal relation between the solar wind physical parameters, such as velocity, and coronal physical parameters such as magnetic field energy. To clarify the physical parameters which control the solar wind velocity, we have studied the relation between solar wind velocity and properties of its source region such as photospheric/coronal magnetic field and the size of each coronal hole during the solar minimum. The solar wind velocity structures were derived by using interplanetary scintillation tomography obtained at Solar-Terrestrial Environment Laboratory, Japan. Potential magnetic fields were calculated to identify the source region of the solar wind. HeI 1083 nm absorption line maps obtained at Kitt Peak National Solar Observatory were used to identify coronal holes. As a result, we found a relation during solar minimum between the solar wind velocity and the coronal magnetic condition which is applicable to different kind of solar winds from different kind of coronal holes.     

电子环束流对航天器表面电位和等离子体鞘层的影响

采用二维三分量静电粒子模拟程序研究了电子环束流对低轨道磁化等离子体中运行的航天器表面电位和等离子体鞘层结构的影响.当有电子束流平行磁场入射时,航天器表面所带负电位的绝对值将增大.电子束流的速度和密度的增加将使航天器表面负电位绝对值增加很快.与此同时等离子体鞘层的尺度也将增大,形状将由向尾部延伸的“泪滴”状变成向两侧展开的“机翼”状.当电子环束流以相对于磁场成一定角度斜入射时,由于磁场的约束作用,航天器表面负电位的绝对值将随入射角增大而减小.      

流向磁场作用下二维磁流体槽道湍流直接模拟

对低磁雷诺数近似下流向磁场作用的二维磁流体槽道湍流进行了直接数值模拟(DNS,Direct Numerical Simulation),给出了 Re=10 000时不同磁相互作用数下二维磁流体槽道湍流的近壁速度分布、湍流脉动速度均方根、雷诺应力等统计量的变化,并与中性流体二维槽道湍流进行了比较.结果表明:流向磁场作用会导致对数区上移;雷诺应力最大值随磁相互作用数呈线性变化;随着磁相互作用数的增大,下壁面平均涡量的时间演化由拟周期性向周期性转变,且脉动周期逐渐增大并当流动层流化后下壁面平均涡量成为常值.      

低电离层elves现象的模拟研究

建立了雷暴放电产生的电磁场在对流层-低电离层之间耦合的时变模式,并用模式研究elves现象的时空特征.模拟结果表明,雷暴放电产生的电场包括准静电场和辐射场,其中由放电电流脉冲产生的辐射场是产生elves现象的直接原因.辐射场的电场强度分布与偶极辐射产生的电场强度的分布相似,在雷暴放电的正上方电场较弱,在地面附近辐射电场较强.模拟结果还表明,elves现象在向外扩张时具有双峰或单峰结构,在摄像仪的底片上相应出现双环或单环结构.elves现象在83 km高度水平向外扩张的过程中,最大光辐射强度开始增大,然后逐渐减小.此外,利用球形电磁波假设,探讨了elves现象超光速水平扩张的原因.      

Analytical-HZETRN model for rapid assessment of active magnetic radiation shielding

The use of active radiation shielding designs has the potential to reduce the radiation exposure received by astronauts on deep-space missions at a significantly lower mass penalty than designs utilizing only passive shielding. Unfortunately, the determination of the radiation exposure inside these shielded environments often involves lengthy and computationally intensive Monte Carlo analysis. In order to evaluate the large trade space of design parameters associated with a magnetic radiation shield design, an analytical model was developed for the determination of flux inside a solenoid magnetic field due to the Galactic Cosmic Radiation (GCR) radiation environment. This analytical model was then coupled with NASA’s radiation transport code, HZETRN, to account for the effects of passive/structural shielding mass. The resulting model can rapidly obtain results for a given configuration and can therefore be used to analyze an entire trade space of potential variables in less time than is required for even a single Monte Carlo run. Analyzing this trade space for a solenoid magnetic shield design indicates that active shield bending powers greater than ∼15 Tm and passive/structural shielding thicknesses greater than 40 g/cm² have a limited impact on reducing dose equivalent values. Also, it is shown that higher magnetic field strengths are more effective than thicker magnetic fields at reducing dose equivalent.     

行星际结构与垂直无碰撞激波的相互作用

应用一维混合模拟方法数值研究了两种行星际结构──反向磁场和高密度等离子团与垂直无碰撞激波的相互作用．结果表明，随着激波上游区磁场的反向，下游区磁场将逐渐改变符号，且等离子体密度和速度分别呈现较强的湍动．激波上游和下游的物理量依然满足Rankine-Hugonoit关系．当高密度的等离子体团通过垂直无碰撞激波时，部分质子被激波反射，部分质子被加速并进入下游区域．由于质子速度分布为非Maxwell分布，在激波下游也激发出较强的湍动．     

Observations of an active region filament

An active region filament was well observed on September 4, 2002 with THEMIS at the Teide observatory and SOHO/MDI. The full Stokes parameters of the filament were obtained in Hα and FeI 6302 Å lines. Using the data, we have studied the fine structure of the filament and obtained the parameters at the barb endpoints, including intensity, velocity and longitudinal magnetic field. Our results indicate: (a) the Doppler velocities are quiet different at barb endpoints; (b) the longitudinal magnetic fields at the barb endpoints are very weak; (c) there is a strong magnetic field structure under the filament spine.     

NUMERICAL MODELLING OF VELOCITY AND TEMPERATURE DISTRIBUTIONS OF THE BUOYANCY CONVECTION EFFECT IN KNbO3 MELT

Numerical modelling of velocity and temperature fields in high-temperature KNbO3 melt of a loop-shaped Pt wire heater is carried out by using the commercial com putational code ANSYS for the mathematical solution of the governing equations.Based on the experimental boundary conditions and the Boussinesq approximation,the numerical modelling of a steady and two-dimensional model is applied to study the process under consideration of the buoyancy-driven convection condition. The result is compared with the previous experimental and theoretical data obtained in our laboratory, and the former is in agreement with the latter. Thus a theoretical guide for reasonable growth conditions is provided by studying in depth the real fluid flow effects in the crystal growth from the melt.     

Magnetogasdynamic shock wave generated by a moving piston in a rotational axisymmetric isothermal flow of perfect gas with variable density

The propagation of a strong cylindrical shock wave in an ideal gas with azimuthal magnetic field, and with or without axisymmetric rotational effects, is investigated. The shock wave is driven out by a piston moving with time according to power law. The ambient medium is assumed to have radial, axial and azimuthal component of fluid velocities. The fluid velocities, the initial density and the initial magnetic field of the ambient medium are assumed to be varying and obey power laws. Solutions are obtained, when the flow between the shock and the piston is isothermal. The gas is assumed to have infinite electrical conductivity and the angular velocity of the ambient medium is assumed to be decreasing as the distance from the axis increases. It is expected that such an angular velocity may occur in the atmospheres of rotating planets and stars. The shock wave moves with variable velocity and the total energy of the wave is non-constant. The effects of variation of the initial density and the Alfven-Mach number on the flow-field are obtained. A comparison is also made between rotating and non-rotating cases.     

Lagrangian analysis of satellite-derived currents: Application to the North Western Mediterranean coastal dynamics

Optimal interpolation methods for improving the reconstruction of coastal dynamics from along-track satellite altimetry measurements have been recently developed over the North Western Mediterranean Sea. Maps of satellite-derived geostrophic current anomalies are generated using these methods, and added to different mean circulation fields in order to obtained absolute geostrophic currents. The resulting fields are then compared to standard AVISO products, and their accuracies are assessed with Lagrangian diagnostics. The trajectories of virtual particle clusters are simulated with a Lagrangian code either with new current fields or with the AVISO ones. The simulated trajectories are then compared to 16 in situ drifter trajectories to evaluate the performance of the different velocity fields. The comparisons show that the new current fields lead to better results than the AVISO one, especially over the shallow continental shelf of the Gulf of Lion. However, despite the use of innovative strategies, some altimetry limitations still persist in the coastal domain, where small scale processes remain sub-sampled by conventional altimetry coverage but will benefit from technological development in the near future. Some of the limitations of the Lagrangian diagnostics presently used are also analyzed, but dedicated studies will be required for future further investigations.     

The atmospheric cosmic- and solar energetic particle radiation environment at aircraft altitudes

Galactic cosmic rays interact with the solar wind, the earth's magnetic field and its atmosphere to produce hadron, lepton and photon fields at aircraft altitudes. In addition to cosmic rays, energetic particles generated by solar activity bombard the earth from time to time. These particles, while less energetic than cosmic rays, also produce radiation fields at aircraft altitudes which have qualitatively the same properties as atmospheric cosmic rays. We have used a code based on transport theory to calculate atmospheric cosmic-ray quantities and compared them with experimental data. Agreement with these data is seen to be good. We have then used this code to calculate equivalent doses to aircraft crews. We have also used the code to calculate radiation doses from several large solar energetic particle events which took place in 1989, including the very large event that occurred on September 29^th and 30^th of that year. The spectra incident on the atmosphere were determined assuming diffusive shock theory.     

歼击机座舱空气流动和传热模拟实验

对某型歼击机座舱的空气流动和传热进行了全尺寸地面模拟实验.模拟了气动热载荷和人体热载荷,并使用ATM2400型36通道多点气流和温度测量系统测量了5个测试面的速度场以及其中两个测试面和邻近服装表面的温度场.观察了供气流量和集中喷嘴开关对舱内空气流动状态和温度分布的影响,评价了该座舱的舒适性.结果表明:舱内速度场 并非相对座舱纵向中心面对称分布,而座舱内壁面温度分布亦极不均匀;集中喷嘴全开和全关时,舱内气流速度场和温度场及邻近服装表面空气温度的分布情况差别显著,但对舱内气流平均速度和平均温度的影响不大;舱内气流平均速度随供气量呈线性变化.实验结果可为数值模拟提供验证依据.     

Unusually strong magnetic fields in Titan’s ionosphere: T42 case study

Observations of unusually large magnetic fields in the ionosphere indicate periods of maximum stress on Titan’s ionosphere and potentially of the strongest loss rates of ionospheric plasma. During Titan flyby T42, the observed magnetic field attained a maximum value of 37 nT between an altitude of 1200 and 1600 km, about 20 nT stronger than on any other Titan pass and close to five times greater in magnetic pressure. The strong fields occurred near the corotation-flow terminator rather than at the sub-flow point, suggesting that the flow which magnetized the ionosphere was from a direction far from corotation and possibly towards Saturn. Extrapolation of solar wind plasma conditions from Earth to Saturn using the University of Michigan MHD code predicts an enhanced solar wind dynamic pressure at Saturn close to this time. Cassini’s earlier exits from Saturn’s magnetosphere support this prediction because the Cassini Plasma Spectrometer instrument saw a magnetopause crossing three hours before the strong field observation. Thus it appears that Titan’s ionosphere was magnetized when the enhanced solar wind dynamic pressure compressed the Saturnian magnetosphere, and perhaps the magnetosheath magnetic field, against Titan. The solar wind pressure then decreased, leaving a strong fossil field in the ionosphere. When observed, this strong magnetic flux tube had begun to twist, further enhancing its strength.     

Global characteristics of ionospheric electric fields and disturbances during the first hours of magnetic storms

The ionospheric plasma density can be significantly disturbed during magnetic storms. In the conventional scenario of ionospheric storms, the negative storm phases with plasma density decreases are caused by neutral composition changes, and the positive storm phases with plasma density increases are often related to atmospheric gravity waves. However, recent studies show that the global redistribution of the ionospheric plasma is dominated primarily by electric fields during the first hours of magnetic storms. In this paper, we present the measurements of ionospheric disturbances by the DMSP satellites and GPS network during the magnetic storm on 6 April 2000. The DMSP measurements include the F region ion velocity and density at the altitude of ∼840 km, and the GPS receiver network provides total electron content (TEC) measurements. The storm-time ionospheric disturbances show the following characteristics. The plasma density is deeply depleted in a latitudinal range of ∼20° over the equatorial region in the evening sector, and the depletions represent plasma bubbles. The ionospheric plasma density at middle latitudes (20°–40° magnetic latitudes) is significantly increased. The dayside TEC is increased simultaneously over a large latitudinal range. An enhanced TEC band forms in the afternoon sector, goes through the cusp region, and enters the polar cap. All the observed ionospheric disturbances occur within 1–5 h from the storm sudden commencement. The observations suggest that penetration electric fields play a major role in the rapid generation of equatorial plasma bubbles and the simultaneous increases of the dayside TEC within the first 2 h during the storm main phase. The ionospheric disturbances at later times may be caused by the combination of penetration electric fields and neutral wind dynamo process.     

太阳风局部加速的磁流体力学过渡特征

本文讨论了轴对称磁力线管中太阳风从亚声速和亚阿尔文速度加速到超声速和超阿尔文速度的过程。考虑到流动和场的二维效应以后,除去临界声速线M2=1以外,还有临界阿尔文速线M_A2=1和M2+M_A2=1。这时,磁力线管中存在三次过渡,使磁流体力学方程组的类型从椭圆型变为双曲型,再变为椭圆型,最后变为双曲型。由于日冕中的典型声速和阿尔文速度具有相同的量级,而且人们对磁力线管的急剧膨胀很有兴趣,这就需要讨论二维磁流体力学方程组的自(冫合)解。太阳风的局部流管中有可能出现多次过渡的流动特征。