Comparison between numerical and analytical solution of solute transport models

The fate and movement of dissolved substances in soils and groundwater has generated considerable concern for the quality of the subsurface environment. Many analytical solutions for the partial differential equations that describe solute and pollutant movement exist. Numerical solutions are more general, and often more difficult to verify. In order to determine the model error, the examination of the ability of numerical methods compared to analytical methods is strongly recommended. The objective of the study is to make a comparison between numerical and analytical solution models for solute transport equation. In this study, the numerical solution calculated with the WAVE-model is compared with the analytical solution calculated with CXTFTT-model. The study scenarios considered variables such as compartment depth, applied flux at the top and soil dispersivity under steady-state conditions. The simulations depend on 27 solute infiltration scenarios. The solute concentrations were calculated with the WAVE-model and the CXTFIT-model for each scenario. The WAVE-model error was evaluated with three methods: absolute average maximum error, relative average maximum error and relative average area error. The study implied that the WAVE-model error increased with the increase of the compartment depth, decreasing soil dispersivity, and decrease in flux. The study leads to the recommendation to use compartment depth as thin as possible to minimise the WAVE-model error. Furthermore, it is more useful to use several numerical solution models, such as SWMS-2D model, to evaluate and examine the WAVE-model.     

Li同位素在斜顽辉石、镁铝榴石及镁铝尖晶石晶格内扩散与分馏的理论计算研究

斜顽辉石、镁铝榴石和镁铝尖晶石作为辉石族、石榴石族以及尖晶石族中的重要端元,是地球上地幔主要组成矿物。Li同位素是重要的地幔地球化学示踪剂,其在橄榄石、辉石和石榴石等地幔矿物中的扩散分馏的性质对理解Li同位素作为地幔地球化学示踪剂非常重要。我们通过经典力场经验势方法,从原子尺度上计算研究了不同温压条件下Li同位素在斜顽辉石、镁铝榴石和镁铝尖晶石晶格中分别通过不同的填隙机制和取代空位机制迁移的活化能和其在不同晶格位上的分馏效应。我们发现Li同位素是通过取代空位机制在斜顽辉石、镁铝榴石和镁铝尖晶石中进行迁移扩散。Li同位素在不同晶格位上的分馏作用计算表明,在三种矿物中重同位素7Li会更多地进入晶格填隙位中,而6Li则相对更多进入Mg位。温度是影响这种分馏作用的一个关键因素,相应的结果可用来解释地幔Li同位素组成特征及冷却条件下的同位素分馏等科学问题。     

An algorithm for detecting inter-ellipsoid contacts

An inter-ellipsoid contact detection algorithm was developed and used for simulating the behaviour of assemblies of ellipsoid-shaped particles using the well-known discrete element method (DEM). The contact algorithm was implemented in the modified version of the DEM program TRUBAL originally written to simulate the behaviour of assemblies of spheres. The modified program was used to perform deviatoric and axisymmetric compression tests on a 1000 prolate spheroid assembly in periodic space. The obtained stress–strain–dilation curves conform with the experimental evidence both qualitatively and quantitatively.     

中国区域气候变化的模拟和问题

对中国区域气候变化的模拟和问题作了介绍和讨论。本文还论述了侧边界处理方案对模拟结果的影响。最后对当前区域气候模拟中存在的问题和解决途径进行了必要的讨论。     

CCM3模式中LSM积雪方案的改进研究(Ⅱ):全球模拟试验分析

为了进一步检验新LSM积雪方案的气候模拟性能,经文献[1]单站模拟检验后,在本文中又进一步将其加入CCM3模式中,分别积分近两年。模拟结果表明：原LSM和新LSM积雪方案虽然都能再现全球气候的主要平均特征,但经改进后的新积雪方案对全球积雪气候、降水等要素场及环流等的气候模拟性能比原LSM积雪方案更好些。     

The design of diamond Compton telescope

We have developed radiation detectors using the new synthetic diamonds. The diamond detector has an advantage for observations of “low/medium” energy gamma rays as a Compton telescope. The primary advantage of the diamond detector can reduce the photoelectric effect in the low energy range, which is background noise for tracking of the Compton recoil electron. A concept of the Diamond Compton Telescope (DCT) consists of position sensitive layers of diamond-striped detector and calorimeter layer of CdTe detector. The key part of the DCT is diamond-striped detectors with a higher positional resolution and a wider energy range from 10 keV to 10 MeV. However, the diamond-striped detector is under development. We describe the performance of prototype diamond detector and the design of a possible DCT evaluated by Monte Carlo simulations.      

Kuiper带天体的原始分布模拟

用包括太阳、8颗大行星、冥王星和UB313以及无质量实验粒子在内的N体问题的天体动力学模型,取当前观测的天体轨道根数为初始条件,对具有确定轨道根数的551个Kuiper主带内的小天体进行了10亿年的轨道反演数值模拟.结果显示:当前观测的这些Kuiper 天体中的1/3以上在10亿年前就位于该区域,少部分位于海王星轨道之内,其他在5OAU之外;在4.5亿年前,整个Kuiper主带内的天体呈较好的正态分布,海王星3:2共振带内没有像今天这样的天体聚集现象.     

嵌埋星团质量分层的数值模拟研究

目前,对大质量恒星形成的初始条件不很清楚,特别是大质量恒星是否形成于星团中心仍有争议.有人从时标上考虑,认为嵌埋星团的质量分层现象意味着大质量恒星只能在星团中心诞生.利用Monte Carlo方法对嵌埋星团的动力学演化进行了数值模拟,并与观测进行了比较.假设初始时刻大质量恒星随机分布,一定比例的嵌埋星团因为大质量恒星的随机运动,在演化的某个时刻会呈现暂态质量分层,其中一部分相当明显.这说明,大质量恒星在中心形成并非嵌埋星团质量分层现象的唯一解释.此外,气体的动力学摩擦能有效地减小动力学质量分层的时标,从而增大暂态质量分层的概率.     

Iap general circulation models: a first step towards developing a local area model for weather prediction in nigeria

１．ＩｎｔｒｏｄｕｃｔｉｏｎＮｕｍｅｒｉｃａｌｍｏｄｅｌｓｆｏｒｗｅａｔｈｅｒｐｒｅｄｉｃｔｉｏｎａｒｅｂｅｃｏｍｉｎｇｍｏｒｅｒｅｌｉａｂｌｅａｎｄｐｏｐｕｌａｒｉｎｍｅｔｅｏｒｏｌｏｇｙａｎｄｃｌｉｍａｔｅｒｅｓｅａｒｃｈ．Ｍｕｌｔｉ－ｌａｙｅ...     

Indian summer monsoon rainfall: Dancing with the tunes of the sun

There is strong statistical evidence that solar activity influences the Indian summer monsoon rainfall. To search for a physical link between the two, we consider the coupled cloud hydrodynamic equations, and derive an equation for the rate of precipitation that is similar to the equation of a forced harmonic oscillator, with cloud and rain water mixing ratios as forcing variables. Those internal forcing variables are parameterized in terms of the combined effect of external forcing as measured by sunspot and coronal hole activities with several well known solar periods (9, 13 and 27 days; 1.3, 5, 11 and 22 years). The equation is then numerically solved and the results show that the variability of the simulated rate of precipitation captures very well the actual variability of the Indian monsoon rainfall, yielding vital clues for a physical understanding that has so far eluded analyses based on statistical correlations alone. We also solved the precipitation equation by allowing for the effects of long-term variation of aerosols. We tentatively conclude that the net effects of aerosols variation are small, when compared to the solar factors, in terms of explaining the observed rainfall variability covering the full Indian monsoonal geographical domains.