Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

A numerical model for Type II supernovae — The collapse stage

We use the following numerical model for the collapse stage of a Type II supernova of 15 M_⊙. Our electron capture rate includes the effects of the inverse reaction and the neutron-proton mass difference. This decreases the electron density at the collapse stage and led to rather large values of the maximum inward velocity and of the corresponding mass (Umax = 3.06 × 10⁹cm/s, Mmax=0.76 M_⊙). These larger values are more favourable for the propagation of shock after the rebounce and the triggering-off of a Type-II supernova explosion. For neutrino transport, we use a leakage model and an equilibrium diffusion model, respectively, for the thin and thick stages and a grey atmosphere model to assess the effect of neutrino precipitation on the collapse. We found this effect to be small, the energy precipitation to be not more than 10^?5 the neutrino energy loss and the momentum precipitation not more than 10^?6 the gravitational acceleration.     

Evolution of Electric Currents Associated with Two M-Class Flares

Using one-minute cadence vector magnetograms from Big Bear Solar Observatory (BBSO), we analyze the temporal behavior of derived longitudinal electric currents associated with two flares on July 26, 2002. One of the events is an M1.0 flare which occurred in active region NOAA 10044, while the other is an M8.7 flare in the adjacent region 10039. Rapid changes of magnetic fields in the form of flux emergence are found to be associated with both of these events. However, the temporal behavior of electric currents are very different. For the M1.0 flare, the longitudinal electric current density drops rapidly near the flaring neutral line; while for the M8.7 flare, the current density rapidly increases, confirming the picture of the current-carrying flux emergence. We offer a possible explanation for such a difference: magnetic reconnection at different heights for the two events, near the photosphere for the M1.0 flare, and higher up for the M8.7 flare.     

Observational evidence for solar coronal decimetric radio pulsations with very short periods

Using the decimetric (700–1500 MHz) radio spectrometer and the synchronous observational system with high temporal resolution at four frequencies (1420, 2130, 2840 and 4260 MHz) of Yunnan Observatory, two rare events were observed on 2001 June 24 and 1990 July 30. The former was a small radio burst exhibiting pulsations with short periods (about 29, 40 and 100 ms) in the impulsive phase. The latter was a large radio burst, which at 2840 MHz produced radio pulsations with period of about 30 ms. This paper focuses on pulsations with very short periods in the range of 29–40 ms. The mechanism of generation of such pulsations may be modulation of radio radiation by the periodic trains of whistler packets originating in unstable regions of the corona. Alternatively, these pulsations can be attributed to wave-wave non-linear interactions of electrostatic upper hybrid waves driven by beams of precipitating electrons in flaring loops.     

Prediction Test for the Two Extremely Strong Solar Storms in October 2003

In late October and early November 2003, a series of space weather hazard events erupted in solar-terrestrial space. Aiming at two intense storm (shock) events on 28 and 29 October, this paper presents a Two-Step method, which combines synoptic analysis of space weather–`observing’ and quantitative prediction – ‘palpating’, and uses it to test predictions. In the first step, ‘observing’, on the basis of observations of the source surface magnetic field, interplanetary scintillation (IPS) and ACE spacecraft, we find that the propagation of the shock waves is asymmetric and northward relative to the normal direction of their solar sources due to the large-scale configuration of the coronal magnetic fields, and the Earth is located near the direction of the fastest speed and greatest energy of the shocks. Being two fast ejection shock events, the fast explosion of extremely high temperature and strong magnetic field, and background solar wind velocity as high as 600 and 1000 km s⁻¹, are also helpful to their rapid propagation. According to the synoptic analysis, the shock travel times can be estimated as 21 and 20 h, which are close to the observational results of 19.97 and 19.63 h, respectively. In the second step, ‘palpating’, we adopt a new membership function of the fast shock events for the ISF method. The predicted results here show that for the onset time of the geomagnetic disturbance, the relative errors between the observational and the predicted results are 1.8 and 6.7%, which are consistent with the estimated results of the first step; and for the magnetic disturbance magnitude, the relative errors between the observational and the predicted results are 4.1 and 3.1%, respectively. Furthermore, the comparison among the predicted results of our Two-Step method with those of five other prevailing methods shows that the Two-Step method is advantageous in predicting such strong shock event. It can predict not only shock arrival time, but also the magnitude of magnetic disturbance. The results of the present paper tell us that understanding the physical features of shock propagation thoroughly is of great importance in improving the prediction efficiency.     

利用自回归模型检测地球极移重力效应

同时采用4个台站的国际超导重力仪长期连续重力观测资料和国际地球自转服务中心提供的同步地球自转参数，研究了极移引起的地球重力场变化特征。利用自回归模型估计了各序列的功率谱密度和积谱密度，结果表明极移导致的重力效应的主要能量集中在Chandler摆动和周年项附近，叠积后实际重力观测与极移重力信号理论值之间的差异分别为0．4％和3．9％，说明超导重力仪可有效监测极移导致的重力变化。     

桐柏─大别山北麓金矿带地质物化探找矿模型

从桐柏─大别山成矿带区域地质、地球物理、地球化学特征的研究入手，分析与金成矿有关的地层、构造、岩浆岩等控矿因素的地球物理场及地球化学特征；建立区域金矿找矿模型。     

西宁盆地泥石流类型及其防治意见

概略地叙述了西宁地区的泥石流形成条件和形成机理,并将西宁地区泥石流划分为稀性泥石流类及泥石质泥石流种类。在此基础上提出防治泥石流灾害的意见。