数值逼近法在丹东市大气降尘量计算中应用

把被定义区域的污染密度看成一张超曲面,再把该超曲面划分为形状相同相互连结的小单元,如三角形,在每个三角形顶点上设一个监测点,测该点的污染密度,以此来推导任一三角形的污染密度和整个区域的平均污染密度公式。把它应用到大气降尘上,就可求出大气降尘量。为大气降尘量和其他污染密度的计算,提供了一个新的方法。      

水质管理规划中BOD5响应矩阵的一种改型模式

本文推导对水质管理规划中经典BOD5响应矩阵的一种改型模式,并对数值例加以阐述.该模式以BOD5总量为控制对象,任一河流断面的BOD5输入量在此响应矩阵中仅出现一次,从而使响应关系由下三角矩阵改变为每一行均含三个变量的稀疏矩阵。这一改型不仅提高了计算速度,更重要的是,由于改型模式具有运输问题类型良好的特殊结构,可以方便地判断一个水质管理规划问题是否可行和给出可行的充分必要条件,这一特性具有较大的实际经济意义.改型模式与经典模式所需数据完全相同。      

周边式二次沉淀池流态的计算模拟研究的数值方法

利用一种半解析半离散的数值方法来求解偏微分方程,从而确定二次沉淀池的流态。该方法具有精度高、速度快、占用计算机内存小的特点。所编制软件可在486微机上运行。计算结果与文献相比较,一致性很好,从而验证了本文方法的正确性。还给出了池型尺寸变化对池中心上部水流死区的影响,为具体设计者提供了参考依据。为进一步的优化设计创造了条件。      

一维平流输送问题的欧拉-泰勒-伽辽金算法及数值试验

介绍了一维常系数平流输送问题的Euler-Taylor-Galerkin(ETG)算法?为获得精确的时间导数,该算法采用了包括二阶和三阶时间导数的前向泰勒级数展开,这些导数值可从空间控制微分方程求出?由此产生一般的时间分离方程?该方程采用标准的Bubnov-Galerkin有限元方法在空间离散化?在使用线性元(帽子函数)和Euler时间差分格式时,可获得精确的Taylor-Galerkin算法?当采用2步显示步骤时,算法为计算效率很高的显式?数值实验结果显示,ETG算法的误差很小,而相误差更可忽略不计,结果令人满意      

Modelling of CO2 arrival time at Ketzin – Part I

The injection of CO₂ at the Ketzin storage site and the chemical detection of its arrival in the observation well allowed testing of different numerical simulation codes. ECLIPSE 100 (E100, black-oil simulator), ECLIPSE 300 (E300, compositional CO₂STORE) and MUFTE-UG were used for predictive modelling applying a constant injection rate of 1 kg s^?1 CO₂ and for a history match applying the actual variable injection rate which ranged from 0 to 0.7 kg s^?1 and averaged 0.23 kg s^?1. The geological model applied, is based on all available geophysical and geological information and has been the same for all programs.The results of the constant injection regime show a good agreement among the programs with a discrepancy of 21–33% for the CO₂ arrival times. However, it is determined from the comparison of the cumulative mass of CO₂ at the time of CO₂ arrival that the injection regime is an important factor for the accurate prediction of CO₂ migration within a saline aquifer. Comparing the actual variable injection regime with the simulations applying a constant injection rate the results are relatively inaccurate.Regarding the actual variable injection regime, which was evaluated using all three simulators, the computational results show a good agreement with the data actually measured at the first observation well. Here, the calculated arrival times exceeded the actual ones by 8.1% (E100), 9.2% (E300) and 17.7% (MUFTE-UG).It can be concluded that irrespective of the deviations of the simulations, due to combinations of different codes and slight differences in input parameters, all three programs are well equipped to give a reliable estimate of the arrival of CO₂. Deviations in the results mainly occur due to different input data and grid size choices done by the different modelling teams working independently of each other. Deviations of the simulations results compared to the actual CO₂ arrival time result from uncertainties in the implementation of the geological model, which was set up based on well log data and analogue studies.     

Indoor Aerosol Modeling: Basic Principles and Practical Applications

The type and amount of indoor air pollutants affects the comfort and quality of indoor environments. Therefore, indoor air quality is an important issue with different social, economic, and health aspects because people in developing countries spend most of their time indoors being exposed to different kinds of indoor pollutants. The indoor air quality can be assessed empirically by measuring the pollutant concentrations or can be predicted by means of mathematical models. An indoor aerosol model describes the dynamic behavior of indoor air pollutants. The basic concept of indoor air models is the mass-balance-conservation where several factors that govern the indoor particle concentrations can be described. These factors may include direct emissions from indoor sources, outdoor aerosol particles penetrating indoors as a result of the ventilation and filtration processes, deposition onto indoor surfaces, and removal from indoor air by means of ventilation. Here we present principles of indoor aerosol models and we also give examples of different kind of models.     

A screening model for fate and transport of biodegradable polyesters in soil

A numerical model for predicting the fate and transport of biodegradable polyester residues in soil, following successive applications of mulch film, was developed and applied. The polymer, applied on surface soil, was assumed to be converted into by-products (monomers), according to a first order kinetics with constant K(1deg). The monomers released were assumed to sorb on soil organic matter (according to a first-order kinetics with constant K(s)), to be leached with the seepage water, through vertical advection and hydrodynamic dispersion, and biodegraded (according to a first-order kinetics with constant K(b)). Results suggested that, to assess a possible build-up of mulch film (as a polymer) on the surface soil, the degradation constant K(1deg) relating the polymer conversion to by-products should be known, whereas the biodegradation constant K(b) indicates there is no danger of groundwater pollution. Likewise, on the basis of by-product concentration in deep soil, soil pollution should not occur.     

Monte Carlo based calibration of scintillation detectors for laboratory and in situ gamma ray measurements

The calibration of scintillation detectors for gamma radiation in a well characterized setup can be transferred to other geometries using Monte Carlo simulations to account for the differences between the calibration and the other geometry. In this study a calibration facility was used that is constructed from bricks of well-known activity concentrations of ⁴⁰K and of radionuclides from the ²³⁸U- and ²³²Th-series. Transfer of the calibration was attempted to a Marinelli beaker geometry with the detector inside a lead shield and to an in situ application with the detector positioned on a sand bed. In general this resulted in good correspondence (within 5-10%) between the activity concentrations derived using the transferred calibration and activities that were derived by independent measurements. Some discrepancies were identified that were attributed to coincident summing in the natural decay series and interference of radon.     

H2O/CO2污染对煤油燃料超声速燃烧影响数值研究

在纯净空气与H₂O/ CO₂污染空气来流对比试验结果基础上,采用数值计算方法和化学动力学方法,研究了H₂O和CO₂污染组分对煤油燃料超声速燃烧的影响,获得了试验手段难以得到的燃烧室流场参数和性能数据。完成了相应的煤油燃料超声速燃烧室二维数值计算,其中匹配了进口总温、总压、马赫数、氧气摩尔分数和工作当量油气比。将数值计算结果与相应试验测量值进行了对比分析,并结合燃烧室流场数据、性能参数分析了H₂O和CO₂污染的动力学影响、以及对燃烧室性能的影响。研究表明:(1)数值计算结果与实验测量值总体上吻合,两种手段均体现了纯净空气来流时不同煤油当量油气比的燃烧室性能,并反映了一致的“污染效应”影响趋势;(2) H₂O污染、H₂O+ CO₂污染的存在降低了煤油燃料超声速燃烧室性能,体现在燃烧诱导压升、燃烧效率、流向冲量增量的下降,而且随着污染组分含量的增加,燃烧室性能下降越加显著。      

Modeling colloid-facilitated transport of multi-species contaminants in unsaturated porous media

Colloid-facilitated transport has been recognized as a potentially important and overlooked contaminant transport process. In particular, it has been observed that conventional two phase sorption models are often unable to explain transport of highly sorbing compounds in the subsurface appropriately in the presence of colloids. In this study a one-dimensional model for colloid-facilitated transport of chemicals in unsaturated porous media is developed. The model has parts for simulating coupled flow, and colloid transport and dissolved and colloidal contaminant transport. Richards' equation is solved to model unsaturated flow, and the effect of colloid entrapment and release on porosity and hydraulic conductivity of the porous media is incorporated into the model. Both random sequential adsorption and Langmuir approaches have been implemented in the model in order to incorporate the effect of surface jamming. The concept of entrapment of colloids into the air-water interface is used for taking into account the effect of retardation caused due to existence of the air phase. A non-equilibrium sorption approach with options of linear and Langmuir sorption assumptions are implemented that can represent the competition and site saturation effects on sorption of multiple compounds both to the solid matrix and to the colloidal particles. Several demonstration calculations are performed and the conditions in which the non-equilibrium model can be approximated by an equilibrium model are also studied.