Laboratory Astrophysics Experiments for Simulation Code Validation: A Case Study

The growing field of Laboratory Astrophysics seeks to study the extreme environments found in many astrophysical events in the controlled setting of a laboratory. In addition to the opportunity to perform basic research into the nature and properties of materials in astrophysical environments, laboratory astrophysics experiments serve beautifully for validating calculations performed by simulation codes designed to model astrophysical phenomena. I present results from our ongoing validation effort for FLASH, a parallel adaptive-mesh hydrodynamics code for the compressible, reactive flows of astrophysical environments. The first test case is a laser-driven shock propagating though a multilayer target introducing Rayleigh–Taylor and Richtmyer–Meshkov fluid instabilities at the material interfaces. The second is an accelerating fluid interface that is subject to the Rayleigh–Taylor instability. We found good agreement between simulations and experiment for the multilayer target case, but disagreement between experiment and simulation in the Rayleigh–Taylor case. I discuss our findings and possible reasons for the disagreement.     

3-D BAROCLINIC NUMERICAL SIMULATION OF THE SOUTH CHINA SEA I. UPPER CIRCULATION

A three-dimensional baroclinic shelf sea model was employed to simulate the seasonal characteristics of the South China Sea (SCS) upper circulation. The results showed that: in summer, an anticyclonic eddy, after its formation between the Bashi Channel and Dongsha Islands in the northeastern SCS, moves southwestward until it disperses slowly. There exists a northward western boundary current along the east shore of the Indo-China Peninsula in the western SCS and an anticyclonic gyre in the southern SCS. But at the end of summer and beginning of autumn, a weak local cyclonic eddy forms in the Nansha Trough, then grows slowly and moves westward till it becomes a cyclonic gyre in the southern SCS in autumn. At the beginning of winter, there exists a cyclonic gyre in the northern and southern SCS, and there is a southward western boundary current along the east shore of the Indo-China Peninsula. But at the end of winter, an anticyclonic eddy grows and moves toward the western boundary after forming in the Nansha Trough. The eddy‘s movement induces a new opposite sign eddy on its eastern side, while the strength of the southward western boundary current gets weakened. This phenomenon continues till spring and causes eddies in the southern SCS.     

An iterative technique for solving equations of statistical equilibrium

Superlevel partitioning is combined with a simple relaxation procedure to construct an iterative technique for solving equations of statistical equilibrium. In treating an N -level model atom, the technique avoids the N ³ scaling in computer time for direct solutions with standard linear equation routines and also does not fail at large N due to the accumulation of round-off errors. As a consequence, the technique allows detailed model atoms with N ≳10³ , such as those required for iron peak elements, to be incorporated into diagnostic codes for analysing astronomical spectra. Tests are reported for a 394-level Fe  ii ion and a 1266-level Ni  i – iv atom.     

The 2dF Galaxy Redshift Survey: luminosity dependence of galaxy clustering

We investigate the dependence of the strength of galaxy clustering on intrinsic luminosity using the Anglo-Australian two degree field galaxy redshift survey (2dFGRS). The 2dFGRS is over an order of magnitude larger than previous redshift surveys used to address this issue. We measure the projected two-point correlation function of galaxies in a series of volume-limited samples. The projected correlation function is free from any distortion of the clustering pattern induced by peculiar motions and is well described by a power law in pair separation over the range     . The clustering of     galaxies in real space is well-fitted by a correlation length     and power-law slope     . The clustering amplitude increases slowly with absolute magnitude for galaxies fainter than M *, but rises more strongly at higher luminosities. At low luminosities, our results agree with measurements from the Southern Sky Redshift Survey 2 by Benoist et al. However, we find a weaker dependence of clustering strength on luminosity at the highest luminosities. The correlation function amplitude increases by a factor of 4.0 between     and −22.5, and the most luminous galaxies are 3.0 times more strongly clustered than L * galaxies. The power-law slope of the correlation function shows remarkably little variation for samples spanning a factor of 20 in luminosity. Our measurements are in very good agreement with the predictions of the hierarchical galaxy formation models of Benson et al.     

Morphology, photometry and kinematics of N-body bars – I. Three models with different halo central concentrations

We discuss the morphology, photometry and kinematics of the bars which have formed in three N -body simulations. These have initially the same disc and the same halo-to-disc mass ratio, but their haloes have very different central concentrations. The third model includes a bulge. The bar in the model with the centrally concentrated halo (model MH) is much stronger, longer and thinner than the bar in the model with the less centrally concentrated halo (model MD). Its shape, when viewed side-on, evolves from boxy to peanut and then to 'X'-shaped, as opposed to that of model MD, which stays boxy. The projected density profiles obtained from cuts along the bar major axis, for both the face-on and the edge-on views, show a flat part, as opposed to those of model MD which are falling rapidly. A Fourier analysis of the face-on density distribution of model MH shows very large  m=2  , 4, 6 and 8 components. Contrary to this, for model MD the components  m=6  and 8 are negligible. The velocity field of model MH shows strong deviations from axial symmetry, and in particular has wavy isovelocities near the end of the bar when viewed along the bar minor axis. When viewed edge-on, it shows cylindrical rotation, which the MD model does not. The properties of the bar of the model with a bulge and a non-centrally concentrated halo (MDB) are intermediate between those of the bars of the other two models. All three models exhibit a lot of inflow of the disc material during their evolution, so that by the end of the simulations the disc dominates over the halo in the inner parts, even for model MH, for which the halo and disc contributions were initially comparable in that region.     

绝热非径向太阳振动的数值计算方法

本文采用改进后的平行Shooting方法，即连续正交归一化方法，通过增加系统的秩和把微分方程从线性变成非线性的方法，使得病态方程变成非病态方程，克服了平行Shooting方法（不连续正交归一化方法）的主要困难：一是选择径向网格点数目的问题，如果被检验振动模的径向节点数大或小，所需要的网格点数则多或少；二是在远离阈值的时候，尚未有合适的方法提高计算精度的问题．数值计算出了标准太阳模型的绝热非径向太阳振动的理论本征频率谱．从数值计算的结果发现改进后的平行Shooting方法比线性化的Henyey方法和平行Shooting方法收敛性更好且精度更高．其误差可以达到1％的精度．