非等间距显式算法在重力波非线性传播数值模拟中的初步应用

采用3阶TVD型显式时间积分方法,建立了2维可压缩大气中重力波非线性传播的非等间距显式数值模式.对小振幅重力波传播过程的模拟结果表明,该模式能够很好地再现小振幅重力波的传播过程,并能保持能量守恒关系和各扰动量之间的相位关系,与线性重力波理论预测的结果吻合很好.对有限振幅重力波的模拟结果表明,与格点数相同的等间距网格中的模拟结果相比,采用非等间距网格能够以更高的分辨率模拟重力波的不稳定传播直到破碎的整个过程;在非等间距网格中的模拟结果与在加密一倍的等间距网格中的模拟结果相同,说明了用非等间距网格模拟重力波的饱和与破碎过程是正确的,并且采用非等间距网格能够大大减少计算量.      

一种全球临近空间大气密度建模方法及应用

基于TIMED/SABER卫星2002—2018年观测的20~100 km大气密度数据,统计获得多年月平均值和标准偏差的全球网格数据。利用网格数据,分析了大气密度的变化特征。以网格数据为基准,计算了USSA76的相对偏差,分析了USSA76相对偏差的分布特征。以网格数据为驱动,将大气密度表征为平均值与大尺度扰动量和小尺度扰动量的加和,大尺度扰动和小尺度扰动分别采用余弦函数和一阶自回归模型表征,初步建立了全球临近空间大气密度模型。通过对比模型仿真值与激光雷达观测值,表明模型仿真值与观测值具有较好的吻合度,验证了建模方法的可行性。利用蒙特卡罗方法可再现给定轨迹上所有可能的大气状态。      

空间电子辐射环境探测载荷测试定标试验平台研制

介绍了中国科学院国家空间科学中心新建成的空间电子辐射环境探测载荷测试定标试验平台.该平台由中、高能极弱流电子加速器以及内置多维真空转台的真空靶室试验终端组成，用于对星载空间电子辐射探测器进行地面加速器测试定标.重点描述了为得到中能极弱流均匀平行束，采用电子轨迹程序Egun对中能极弱流电子加速器进行的物理设计和模拟计算，给出球形结构电子枪在栅网孔不加栅网、加理想栅网和直径1mm孔栅网以及在不同加速管出口能量情况下，初聚系统和加速管以及经过二次扩束时输运段中电子轨迹的模拟结果.最终得出能够实现电子枪初始束流减弱8个数量级，获得满足测试定标试验需求的极弱流均匀平行电子束（在试验终端直径50mm靶上束流面密度为105~109cm-2·s-1）的结论.      

多种观测数据驱动的三维行星际太阳风MHD模拟

三维磁流体力学（MHD）数值模拟是行星际太阳风研究的重要手段.本文发展了一种由多种观测数据驱动的三维行星际太阳风MHD数值模型.模型的计算区域为0.1AU到1AU附近，使用Lax-Friedrich差分格式在六片网格系统中进行数值求解.边界条件中磁场使用GONG台站观测的光球磁图外推获得，密度通过LASCO观测的白光偏振亮度反演得到，速度根据以上两种观测数据并利用一种基于人工神经网络技术（ANN）的方法得到，温度通过自洽方法根据磁场和密度导出.利用该模型模拟了第2062卡灵顿周（CR2062）时期的行星际太阳风，模拟结果显示出丰富的观测特征，并与OMNI以及Ulysses的实际观测值符合得较好.该模型可用于提供接近真实的行星际太阳风，有助于提高空间天气预报的精度.      

嵌套闭磁场日冕背景结构对触发CME特征的影响

在子午面内,研究具有嵌套闭磁场结构冕流背景对触发日冕物质抛射(CME)特征的影响.在冕流背景磁场结构内有三个小尺度的闭磁场结构,其中中间小尺度闭磁场结构的磁场方向和冕流整体偶极磁场方向相反.CME触发模型在这个小尺度闭磁场结构下方浮出,它具有同心圆形磁场结构,半径为a=0.1Rs(Rs为太阳半径);CME触发模型前半部分磁场方向和触发位置处小尺度闭磁场结构的磁场方向相反,与冕流整体偶极磁场方向相同.数值模拟结果表明,当CME触发模型中心等离子压力与边界压力之比m≥2时,上浮CME触发模型可以触发CME;当m＜2时,上浮CME触发模型不能触发CME,计算结论的误差小于1%.      

利用栅网控制空间等离子体环境效应模拟中的等离子体参数

通常采用ECR等离子体源产生的等离子体的温度和密度都比较大，通过附加适当目数的栅网，并在栅网上加一定的偏置电位来对等离子体参数，尤其是温度，密度进行调整，满足空间等离子体环境模拟要求，本文利用14目，25目的栅网，对已有的地面实验室空间等离子体环境进行了改进，得到了更加接近空间等离子体参数的一个地面模拟环境。     

Field-aligned current signature in the dusk magnetosphere by INTERBALL-AU and POLAR

The POLAR and INTERBALL-AU spacecrafts magnetic field experiments allow investigation of the time-spatial variation of field-aligned current structures at mid and high altitudes. The study includes two intervals over north high latitudes in the dusk magnetosphere: (1) 10 January 09 – 11 UT (consecutive transits) and (2) 11 January, 03 – UT (nearly conjugate situation at 04:30 UT). By using both spacecraft, the identification of current regions is more confident. On January 11 cycles of the polarization reversal of small-scale disturbances in the poleward portion of region R1 have been detected at high and mid altitudes.     

等离子体密度对栅极系统束流影响分析

针对现有离子推力器栅极系统束流引出分析中未考虑放电室等离子体密度均匀性影响的不足,基于数值方法对该问题进行了研究。采用单元内粒子 （Particle In Cell,PIC） 方法对栅极系统束流引出进行了数值模拟。PIC方法中电场求解采用有限差分方法,离子加速方法采用跳蛙格式。以国产20cm Xe离子推力器作为算例,计算得到了栅极中心和边缘孔引出束流离子和电场的分布。结果显示栅极中心孔相对边缘孔束流聚焦性好,束流发散角小,但中心孔轴线的鞍点电势较边缘孔高24.8V。因此,放电室等离子体密度分布对束流引出具有重要影响,电子返流现象最先出现在栅极中心位置。     

球坐标系六片网格下三维定态行星际太阳风模拟

采用二阶MacCormack差分格式, 利用稳态的磁流体(MHD)方程组在球坐标系六片网格下模拟研究了行星际太阳风. 六片网格系统能有效避免极区奇性和网格收敛性. 迭代按径向方向推进求解, 很大程度上减少了计算量, 节约了计算时间. 内边界条件根据太阳与行星际观测确定, 比较测试了5种内边界条件, 模拟给出了1922卡林顿周的背景太阳风结构. 几种内边界条件所得模拟结果与行星际观测基本吻合. 太阳风速度采用McGregor 等的经验公式给出, 磁场由水平电流片(HCCS)模型得到, 密度和温度分别根据动量守恒和气压守恒得到, 研究表明采用这样的边界条件模拟结果最佳.      

Transition radiation mechanism for broad-band component of stimulated electromagnetic emission

As is well known in the F-region of the ionosphere modified by high power HF radio waves broad-band electromagnetic stimulated emission (SEE) is observed. It was discovered both the beams of superthermal electrons and intensive small-scale irregularities in modified region. These magnetic field aligned irregularities have caviton shape — deep electron density holes. The presence of such irregularities and superthermal electrons create conditions for generation of transition emission. We consider this radiation mechanism for interpretation the SEE broad-band component.     

On high-altitude electric fields in the auroral downward current region and their coupling to ionospheric electric fields

The downward field-aligned current region plays an active role in magnetosphere–ionosphere coupling processes associated with aurora. A quasi-static electric field structure with a downward parallel electric field forms at altitudes between 800 km and 5000 km, accelerating ionospheric electrons upward, away from the auroral ionosphere. Other phenomena including energetic ion conics, electron solitary waves, low-frequency wave activity, and plasma density cavities occur in this region, which also acts as a source region for VLF saucers. Results are presented from high-altitude Cluster observations with particular emphasis on the characteristics and dynamics of quasi-static electric field structures. These, extending up to altitudes of at least 4–5 Earth radii, appear commonly as monopolar or bipolar electric fields. The former occur at sharp boundaries, such as the polar cap boundary whereas the bipolar fields occur at softer boundaries within the plasma sheet. The temporal evolution of quasi-static electric field structures, as captured by the pearls-on-a-string configuration of the Cluster spacecraft, indicates that the formation of electric field structures and of ionospheric plasma density cavities are closely coupled processes. A related feature of the downward current is a broadening of the current sheet with time, possibly related to the depletion process. Preliminary studies of the coupling of electric fields in the downward current region, show that small-scale structures are typically decoupled from the ionosphere, similar to what has been found for the upward current region. However, exceptions are also found where small-scale electric fields couple perfectly between the ionosphere and Cluster altitudes. Recent FAST results indicate that the degree of coupling differs between sheet-like and curved structures, and that it is typically partial. The electric field coupling further depends on the current–voltage relationship, which is highly non-linear in the downward current region, and still unrevealed, as to its specific form.     

模拟三维重力波传播过程的并行数值模式

从三维大气运动的基本控制方程出发, 基于MPI消息传递接口和区域分解的思想, 建立了模拟三维可压缩大气中重力波传播过程的并行数值模式. 在对垂直方向进行区域分解的基础上, 针对跳点网格的特点, 以垂直速度主格点线为子区域的下边界, 状态变量主格点线为子区域的上边界进行区域分解. 利用MPI消息传递接口来传递计算各子区域中变量时所要用到的相邻子区域上相应变量, 从而顺利地完成对整个计算域的并行计算. 根据线性重力波理论, 通过模拟小振幅重力波的传播过程对所建模式进行了验证. 结果表明, 模式可以很好地模拟小振幅重力波在三维空间中的传播过程, 模拟的重力波振幅随着传播高度的增加以指数形式增长, 对能量传播路径以及有效扰动位能和扰动动能的模拟结果均与理论预测的结果吻合很好. 另外, 随着进程个数的增加, 完成相同的计算量所用的计算时间也显著减少. 这表明本文建立的并行数值模式不但能够很好的模拟重力波的传播过程, 而且能够有效地节约计算时间.      

Regional 4-D modeling of the ionospheric electron density from satellite data and IRI

Accurate knowledge of the electron density is the key point in correcting ionospheric delays of electromagnetic measurements and in studying ionosphere physics. During the last decade Global Navigation Satellite Systems (GNSS) have become a promising tool for monitoring ionospheric parameters such as the total electron content (TEC). In this contribution we present a four-dimensional (4-D) model of the electron density consisting of a given reference part, i.e., the International Reference Ionosphere (IRI), and an unknown correction term expanded in terms of multi-dimensional base functions. The corresponding series coefficients are calculable from the satellite measurements by applying parameter estimation procedures. Since satellite data are usually sampled between GPS satellites and ground stations, finer structures of the electron density are modelable just in regions with a sufficient number of ground stations. The proposed method is applied to simulated geometry-free GPS phase measurements. The procedure can be used, for example, to study the equatorial anomaly.     

Spectra and time variability of black-hole binaries in the low/hard state

We propose a jet model for the low/hard state of galactic black-hole X-ray sources which explains the energy spectra from radio to X-rays and a number of timing properties in the X-ray domain such as the time lag spectra, the hardening of the power density spectra and the narrowing of the autocorrelation function with increasing photon energy. The model assumes that (i) there is a magnetic field along the axis of the jet, (ii) the electron density in the jet drops inversely proportional to distance, (iii) the jet is “hotter” near its center than at its periphery, and (iv) the electrons in the jet follow a power-law distribution function. We have performed Monte Carlo simulations of Compton upscattering of soft photons from the accretion disk and have found power-law high-energy spectra with photon-number index in the range 1.5–2 and cutoff at a few hundred keV, power-law time lags versus Fourier frequency with index 0.8, and an increase of the rms amplitude of variability and a narrowing of the autocorrelation function with increasing photon energy as they have been observed in Cygnus X-1. The spectrum at long wavelengths (radio, infrared, optical) is modeled to come from synchrotron radiation of the energetic electrons in the jet. We find flat to inverted radio spectra that extend from the radio up to about the optical band. For magnetic field strengths of the order 10⁵–10⁶ G at the base of the jet, the calculated spectra agree well in slope and flux with the observations.     

SIP-AMR-CESE模式对CR2055背景太阳风的模拟

使用三维太阳行星际自适应网格守恒元解元太阳风模型（SIP-AMR-CESE MHD）,模拟从太阳表面到地球轨道附近的太阳风.该模型使用六片网格技术,同时利用PARAMESH软件包实现网格自适应.在该模型的基础上,通过增加广义拉格朗日乘子（GLM）磁场散度误差消去方法,完善网格加密放粗判据,微调加速加热形式等方法,使模拟结果与观测更好地符合.另外,通过控制不同时刻的计算区域,显著提高了模型的计算效率.在此基础上,给出了模型改进后模拟得到的CR2055太阳风稳态解与观测的对比分析.      

Analysis of coastal altimetry in the Mexican Caribbean

A reprocessing of sea-level anomalies (SLA) resulting from X-TRACK coastal altimetry was carried out for the ENVISAT (2002–2010) and TOPEX/POSEIDON-Jason (1992–2019) satellite missions in the coastal area of the Mexican Caribbean. This consisted of applying a tidal correction to coastal altimetry sea level observations. Harmonic analysis of five coastal tide gauge records was performed to estimate the most important tidal components of the area, resulting on M2, N2, O1, S2, K1, MF, and MM. The tidal signal was reconstructed with the seven tidal components using the TPXO9 model. The SLA signals corrected with the seven tidal components were validated with in situ data from coastal tide gauges. The validation showed that the TPXO9 tidal barotropic model (1/30° grid) used to reconstruct the tidal signal with the seven representative tidal components performed better than the FES2012 global model (1/16° grid) that uses 33 tidal components. The reprocessed SLAs showed clear seasonality with significant signals at 4, 6, and 12 months, with the annual signal being the dominant one. In the Mexican Caribbean coastal zone, oceanographic processes with different scales (from coastal to mesoscale) converge, showing their complexity in the different SLA signals observed. The aim of this work is to contribute to the analysis of coastal altimetry data and understanding the sea level variations in the Mexican Caribbean. This work is the first step in the implementation of methodologies that take advantage of coastal satellite altimetry in the Caribbean Sea.     

小型平板CPL实现高热流密度散热的研究

针对小型平板毛细抽吸两相流体回路(CPL)在高热流密度下的特点,分析了不同工质时系统的压力损失与毛细芯的毛细抽吸力,得出采用氨工质有着较好的传热性能和更高的毛细限,同时得出影响系统毛细限的主要因素是蒸汽联管管径和工质的蒸气密度,提出了工质传输系数作为选取工质的重要指标。建立了蒸发器多孔芯,金属壁面以及工质气、液空间区域的耦合数学模型,并运用SIMPLE算法进行求解,得出蒸发器内的温度分布及由于侧壁效应对多孔芯传热传质与传热极限的影响,同时提出小型平板CPL系统存在侧壁效应传热极限,它是影响系统最大传热量的一个重要极限,在设计小型平板CPL必须予以考虑。     

Computerized ionospheric tomography with the IRI model

Computerized ionospheric tomography (CIT) is a method to estimate ionospheric electron density distribution by using the global positioning system (GPS) signals recorded by the GPS receivers. Ionospheric electron density is a function of latitude, longitude, height and time. A general approach in CIT is to represent the ionosphere as a linear combination of basis functions. In this study, the model of the ionosphere is obtained from the IRI in latitude and height only. The goal is to determine the best representing basis function from the set of Squeezed Legendre polynomials, truncated Legendre polynomials, Haar Wavelets and singular value decomposition (SVD). The reconstruction algorithms used in this study can be listed as total least squares (TLS), regularized least squares, algebraic reconstruction technique (ART) and a hybrid algorithm where the reconstruction from the TLS algorithm is used as the initial estimate for the ART. The error performance of the reconstruction algorithms are compared with respect to the electron density generated by the IRI-2001 model. In the investigated scenario, the measurements are obtained from the IRI-2001 as the line integral of the electron density profiles, imitating the total electron content estimated from GPS measurements. It has been observed that the minimum error between the reconstructed and model ionospheres depends on both the reconstruction algorithm and the basis functions where the best results have been obtained for the basis functions from the model itself through SVD.     

Mathematical modeling of plasma drifts over equatorial low latitude regions

This paper presents a mathematical model to simulate ionospheric plasma drifts at equatorial low latitude regions by coupling of E- and F-regions. The governing non-linear differential equations (of elliptic and parabolic nature) are solved numerically through finite-difference schemes and obtained neutral winds and electric fields. The temperature and electron density profiles are generated utilizing MSIS-86 atmospheric model. The continuity equation is employed to obtain night-time E-region density profile using measured ionograms at Trivandrum (India). The computed vertical and zonal plasma drifts are comparable with measured Jacamarca plasma drifts with little variations during noon and evening times. The plasma drifts at Trivandrum (8.5° N, 76.5° E, dip 0.5° N) are compared with those of Jicamarca (12° S, 76.9° W, dip 2° N). Neutral wind simulations of present model agree well with those of horizontal wind model (HWM-93). The post-sunset enhancement and its reversal are also discussed.