Utilization of CAM,CAF, MAX,and FAX for space radiation analyses using HZETRN

To estimate astronaut health risk due to space radiation, one must have the ability to calculate various exposure-related quantities that are averaged over specific organs and tissue types. Such calculations require computational models of the ambient space radiation environment, particle transport, nuclear and atomic physics, and the human body. While significant efforts have been made to verify, validate, and quantify the uncertainties associated with many of these models and tools, relatively little work has focused on the uncertainties associated with the representation and utilization of the human phantoms. In this study, we first examine the anatomical properties of the Computerized Anatomical Man (CAM), Computerized Anatomical Female (CAF), Male Adult voXel (MAX), and Female Adult voXel (FAX) models by comparing the masses of various model tissues used to calculate effective dose to the reference values specified by the International Commission on Radiological Protection (ICRP). The MAX and FAX tissue masses are found to be in good agreement with the reference data, while major discrepancies are found between the CAM and CAF tissue masses and the reference data for almost all of the effective dose tissues. We next examine the distribution of target points used with the deterministic transport code HZETRN (High charge (Z) and Energy TRaNsport) to compute mass averaged exposure quantities. A numerical algorithm is presented and used to generate multiple point distributions of varying fidelity for many of the effective dose tissues identified in CAM, CAF, MAX, and FAX. The point distributions are used to compute mass averaged dose equivalent values under both a galactic cosmic ray (GCR) and solar particle event (SPE) environment impinging isotropically on three spherical aluminum shells with areal densities of 0.4 g/cm², 2.0 g/cm², and 10.0 g/cm². The dose equivalent values are examined to identify a recommended set of target points for each of the tissues and to further assess the differences between CAM, CAF, MAX, and FAX. It is concluded that the previously published CAM and CAF point distributions were significantly under-sampled and that the set of point distributions presented here should be adequate for future studies involving CAM, CAF, MAX, or FAX. It is also found that the errors associated with the mass and location of certain tissues in CAM and CAF have a significant impact on the mass averaged dose equivalent values, and it is concluded that MAX and FAX are more accurate than CAM and CAF for space radiation analyses.     

Differential code bias estimation based on uncombined PPP with LEO onboard GPS observations

Differential Code Bias (DCB) is an essential correction that must be provided to the Global Navigation Satellite System (GNSS) users for precise position determination. With the continuous deployment of Low Earth Orbit (LEO) satellites, DCB estimation using observations from GNSS receivers onboard the LEO satellites is drawing increasing interests in order to meet the growing demands on high-quality DCB products from LEO-based applications, such as LEO-based GNSS signal augmentation and space weather research. Previous studies on LEO-based DCB estimation are usually using the geometry-free combination of GNSS observations, and it may suffer from significant leveling errors due to non-zero mean of multipath errors and short-term variations of receiver code and phase biases. In this study, we utilize the uncombined Precise Point Positioning (PPP) model for LEO DCB estimation. The models for uncombined PPP-based LEO DCB estimation are presented and GPS observations acquired from receivers onboard three identical Swarm satellites from February 1 to 28, 2019 are used for the validation. The results show that the average Root Mean Square errors (RMS) of the GPS satellite DCBs estimated with onboard data from each of the three Swarm satellites using the uncombined PPP model are less than 0.18 ns when compared to the GPS satellite DCBs obtained from IGS final daily Global Ionospheric Map (GIM) products. Meanwhile, the corresponding average RMS of GPS satellite DCBs estimated with the conventional geometry-free model are 0.290, 0.210, 0.281 ns, respectively, which are significantly larger than those obtained with the uncombined PPP model. It is also noted that the estimated GPS satellite DCBs by Swarm A and C satellites are highly correlated, likely attributed to their similar orbit type and space environment. On the other hand, the Swarm receiver DCBs estimated with uncombined PPP model, with Standard Deviation (STD) of 0.065, 0.037 and 0.071 ns, are more stable than those obtained from the official Swarm Level 2 products with corresponding STD values of 0.115, 0.101, and 0.109 ns, respectively. The above indicates that high-quality DCB products can be estimated based on uncombined PPP with LEO onboard observations.     

First results of operational ionospheric dynamics prediction for the Brazilian Space Weather program

It is shown the development and preliminary results of operational ionosphere dynamics prediction system for the Brazilian Space Weather program. The system is based on the Sheffield University Plasmasphere–Ionosphere Model (SUPIM), a physics-based model computer code describing the distribution of ionization within the Earth mid to equatorial latitude ionosphere and plasmasphere, during geomagnetically quiet periods. The model outputs are given in a 2-dimensional plane aligned with Earth magnetic field lines, with fixed magnetic longitude coordinate. The code was adapted to provide the output in geographical coordinates. It was made referring to the Earth’s magnetic field as an eccentric dipole, using the approximation based on International Geomagnetic Reference Field (IGRF-11). During the system operation, several simulation runs are performed at different longitudes. The original code would not be able to run all simulations serially in reasonable time. So, a parallel version for the code was developed for enhancing the performance. After preliminary tests, it was frequently observed code instability, when negative ion temperatures or concentrations prevented the code from continuing its processing. After a detailed analysis, it was verified that most of these problems occurred due to concentration estimation of simulation points located at high altitudes, typically over 4000 km of altitude. In order to force convergence, an artificial exponential decay for ion–neutral collisional frequency was used above mentioned altitudes. This approach shown no significant difference from original code output, but improved substantially the code stability. In order to make operational system even more stable, the initial altitude and initial ion concentration values used on exponential decay equation are changed when convergence is not achieved, within pre-defined values. When all code runs end, the longitude of every point is then compared with its original reference station longitude, and differences are compensated by changing the simulation point time slot, in a temporal adjustment optimization. Then, an approximate neighbor searching technique was developed to obtain the ion concentration values in a regularly spaced grid, using inverse distance weighting (IDW) interpolation. A 3D grid containing ion and electron concentrations is generated for every hour of simulated day. Its spatial resolution is 1° of latitude per 1° of longitude per 10 km of altitude. The vertical total electron content (VTEC) is calculated from the grid, and plotted in a geographic map. An important feature that was implemented in the system is the capacity of combining observational data and simulation outputs to obtain more appropriate initial conditions to the ionosphere prediction. Newtonian relaxation method was used for this data assimilation process, where ionosonde data from four different locations in South America was used to improve the system accuracy. The whole process runs every day and predicts the VTEC values for South America region with almost 24 h ahead.     

中心线变化规律和膨胀比对圆转方过渡段流场形态的影响

为了获得流道中心线变化规律和截面膨胀比ε对圆转方过渡段流场形态的影响规律,在某型发动机的真实工况下,分别对3种中心线变化规律(前缓后急、缓急相当、前急后缓)和5种膨胀比(1.00,1.01,1.03,1.05,1.07)的圆转方过渡段结构进行建模计算,得到了圆转方过渡段内部流场形态、总压恢复系数和流向涡涡量的变化规律。结果表明:不同结构的横截面静压和速度分布差异较大。前急后缓结构在圆转方横截面出现明显的角涡;在三种中心线变化规律结构中总压恢复系数沿流向均呈下降趋势,流向涡涡量沿流向均呈先增大后减小的趋势,三种结构的流向涡涡量最大值从大到小依次为前急后缓结构、前缓后急结构和缓急相当结构,最大最小值相差19%。随着膨胀比的增大,在圆转方横截面会出现角涡;在不同膨胀比结构中总压恢复系数沿流向均呈下降趋势,流向涡涡量沿流向均呈先增大后减小的趋势,膨胀比1.03结构的局部流向涡无量纲涡量比1.07结构小40.9%。缓急相当中心线变化规律和膨胀比ε=1.03为最佳结构参数。     

On the uncertainties in the measurement of absolute (true) TEC over Indian equatorial and low latitude sectors

The Indian sector encompasses the equatorial and low latitude regions where the ionosphere is highly dynamic and is characterized by the equatorial ionization anomaly (EIA) resulting in large latitudinal electron density gradients causing errors and uncertainties in the estimation of range delays in satellite based navigation systems. The diurnal and seasonal variations of standard deviations in the TEC data measured during the low sunspot period 2004–2005 at 10 different Indian stations located from equator to the anomaly crest region and beyond are examined and presented. The day-to-day variability in TEC is found to be lowest at the equatorial station and increases with latitude up to the crest region of EIA and decreases beyond.     

Estimation of Galactic Cosmic Ray exposure inside and outside the Earth’s magnetosphere during the recent solar minimum between solar cycles 23 and 24

The evidently low solar activity observed between solar cycles 23 and 24 during the years 2008–2010 led to a substantial increase in the Galactic Cosmic Ray (GCR) intensity in comparison with preceding solar minima. As the GCRs consist of highly-ionizing charged particles having the potential to cause biological damage, they are a subject of concern for manned missions to space. With the enhanced particle fluxes observed between 2008 and 2010, it is reasonable to assume that the radiation exposure from GCR must have also increased to unusually high levels. In this paper, the GCR exposure outside and inside the Earth’s magnetosphere is numerically calculated for time periods starting from 1970 to the end of 2011 in order to investigate the increase in dose levels during the years 2008–2010 in comparison with the last three solar minima. The dose rates were calculated in a water sphere, used as a surrogate for the human body, either unshielded or surrounded by aluminium shielding of 0.3, 10 or 40 g/cm².     

结合局部次流循环的变几何轴对称进气道研究

针对常规定几何轴对称进气道在低马赫数工作时流量系数低、溢流阻力大的问题,提出了一种结合局部次流循环的变几何轴对称进气道,其通过平移进气道一级锥并引入局部次流循环重构前体激波系相结合的方法,保证了进气道在较宽马赫数范围内的流量捕获能力。通过仿真方法验证了这一设计概念的可行性,并与常规定几何轴对称进气道进行了性能对比。结果表明:该新概念可调轴对称进气道在低马赫数工作时具有良好的流量捕获性能,并且在整个工作范围内保持了较高的总压恢复性能。与按传统方法设计的定几何轴对称进气道相比,其流量系数和总压恢复系数在工作范围内的最大改善幅度分别达到27.45%和14.31%。此外,选择合适的非控制状态贴口马赫数对该设计概念的实现效果具有明显的影响。     

自由流紊流度对串列叶栅性能的影响

通过实验，测量了自由流紊流度（Ｔｕ）在１．４％－６．７％范围内，一双圆弧串列叶栅总压损失系数、气流转角、附面层内速度分布以及吸力面压力系数。结果表明：（１）当Ｔｕ小于大约３％时，随着Ｔｕ增加，总压损失明显降低；当Ｔｕ大于３％时，Ｔｕ变化对总压损失几乎无影响。（２）Ｔｕ在１．４％－６．７％范围内变化对气流转角无影响。（３）Ｔｕ增加，对前排叶片吸力面靠近前缘部分附面层有有利影响。     

Space weather and commercial airlines

Space weather phenomena can effect many areas of commercial airline operations including avionics, communications and GPS navigation systems. Of particular importance at present is the recently introduced EU legislation requiring the monitoring of aircrew radiation exposure, including any variations at aircraft altitudes due to solar activity. With the introduction of new ultra-long-haul “over-the-pole” routes, “more-electric” aircraft in the future, and the increasing use of satellites in the operation, the need for a better understanding of the space weather impacts on future airline operations becomes all the more compelling. This paper will present the various space weather effects, some provisional results of an ongoing 3-year study to monitor cosmic radiation in aircraft, and conclude by summarising some of the identified key operational issues, which must be addressed, with the help of the science community, if the airlines want to benefit from the availability of space weather services.     

Night-time total electron content enhancements at equatorial anomaly region in China

This paper presents small scale (duration ?1 h, ΔTEC ? 1TECU) night-time total electron content (TEC) enhancements observed at the equatorial anomaly region in China, for the first time. The data is from a GPS receiver chain established in 2005 by Institute of Center for Space Science and Applied Research, Chinese Academy of Sciences and a GPS receiver of International GPS Service (IGS), located between Fuzhou (26.1°N, 119.3°E) and Nanning (22.8°N, 108.3°E). Two other GPS observations of IGS taken at higher latitude were also used to investigate the localization of such phenomenon. The characteristics of the night-time TEC enhancement are examined with two case studies. The TEC increases about 1–3TECU, intermittently. While the night-time TEC enhancement mainly occurs at the equatorial anomaly region, it can be observed at middle latitude as well. The spatial size of the enhancement region is less than 5° in longitude. The primary statistical study shows that the TEC enhancement is more often observed in spring and autumn, but rarely in summer. It has no dependence on geomagnetic activity. The enhancement can occur both before and after midnight.