Density estimation of two-phase flow with multiscale and randomly perturbed data

In this paper, we describe an efficient approach for quantifying uncertainty in two-phase flow applications due to perturbations of the permeability in a multiscale heterogeneous porous medium. The method is based on the application of the multiscale finite element method within the framework of Monte Carlo simulation and an efficient preprocessing construction of the multiscale basis functions. The quantities of interest for our applications are the Darcy velocity and breakthrough time and we quantify their uncertainty by constructing the respective cumulative distribution functions. For the Darcy velocity we use the multiscale finite element method, but due to lack of conservation, we apply the multiscale finite volume element method as an alternative for use with the two-phase flow problem. We provide a number of numerical examples to illustrate the performance of the method.     

Neumann stochastic finite element method for calculating temperature field of frozen soil based on random field theory

To study the effect of uncertain factors on the temperature field of frozen soil, we propose a method to calculate the spatial average variance from just the point variance based on the local average theory of random fields. We model the heat transfer coefficient and specific heat capacity as spatially random fields instead of traditional random variables. An analysis for calculating the random temperature field of seasonal frozen soil is suggested by the Neumann stochastic finite element method, and here we provide the computational formulae of mathematical expectation, variance and variable coefficient. As shown in the calculation flow chart, the stochastic finite element calculation program for solving the random temperature field, as compiled by Matrix Laboratory (MATLAB) software, can directly output the statistical results of the temperature field of frozen soil. An example is presented to demonstrate the random effects from random field parameters, and the feasibility of the proposed approach is proven by comparing these results with the results derived when the random parameters are only modeled as random variables. The results show that the Neumann stochastic finite element method can efficiently solve the problem of random temperature fields of frozen soil based on random field theory, and it can reduce the variability of calculation results when the random parameters are modeled as spatially random fields.     

冻融作用对系统与环境间能量交换的影响

通过求解古典Neumann问题,研究土壤在冻结和融化过程的热交换特征,从理论上分析了冻融作用对系统与环境间能量交换的影响.以草炭亚粘土和亚粘土为例,计算了广泛环境条件下的冻融作用对系统温度变化和与环境间热交换的影响.结果表明,冻融作用使相应过程表层温度梯度加大,热交换强度增强.     

地下水模型的Neumann展开Monte-Carlo随机有限元法

讨论了求解地下水模型(水流模型和水质模型)的Neumann展开Monte-Carlo随机有限元法。从基本的随机变量入手,避免了过程中随机变量的增多问题,给出结点水头(浓度)的均值、方差和水头在某区间的概率计算方法;改进了矩阵求逆的效率,对输入随机变量较多、随机变量变异较大的非稳定地下水问题特别有效。同时选取二维承压地下水水流问题(有解析解)作为例子,进行了随机数值模拟实验。      

Data assimilation methods for estimating a heterogeneous conductivity field by assimilating transient solute transport data via ensemble Kalman filter

An ensemble Kalman filter (EnKF) is developed to identify a hydraulic conductivity distribution in a heterogeneous medium by assimilating solute concentration measurements of solute transport in the field with a steady‐state flow. A synthetic case with the mixed Neumann/Dirichlet boundary conditions is designed to investigate the capacity of the data assimilation methods to identify a conductivity distribution. The developed method is demonstrated in 2‐D transient solute transport with two different initial instant solute injection areas. The influences of the observation error and model error on the updated results are considered in this study. The study results indicate that the EnKF method will significantly improve the estimation of the hydraulic conductivity field by assimilating solute concentration measurements. The larger area of the initial distribution and the more observed data obtained, the better the calculation results. When the standard deviation of the observation error varies from 1% to 30% of the solute concentration measurements, the simulated results by the data assimilation method do not change much, which indicates that assimilation results are not very sensitive to the standard deviation of the observation error in this study. When the inflation factor is more than 1.0 to enlarge the model error by increasing the forecast error covariance matrix, the updated results of the hydraulic conductivity by the data assimilation method are not good at all. Copyright © 2012 John Wiley & Sons, Ltd.     

The GPS-gravimetry boundary value problem

Boundary value problem (BVP) plays a funda-mental role in physical geodesy that aims at determin-ing the earth’s shape and its external gravity field. TheMolodensky BVP and the Stokes BVP are typical inphysical geodesy, and the gravity anomaly is a kind ofbasic data. With the wide use of GPS, measurementaccuracy of the earth’s surface can reach one centime-ter, while that of the gravity measurement can reachμgals. Hence, it is necessary to establish a new kind ofBVP which can satisfy…     

重力矢量及张量信息在地球重力场中的应用

随着重力测量技术日新月异的发展，尤其是空间重力测量技术的发展，使反映地球中、高频信息的重力矢量及张量的获取成为可能。文中建立了基于重力矢量及张量信息求解重力场中扰动场元的数学模型，对于将来重力场的精化及其时变性的研究具有一定的理论意义。     

BED CONDITIONS OF NON-EQUILIBRIUM TRANSPORT OF SUSPENDED SEDIMENT

I.INTRODUCTIONSuspendedsedimenttransportissaidbeinanycaseaphenomenaofa3Dprocess.Becauseofthelimitationsincomputercapacity,mostoftheexistingmodelstreateditasatlvoorevenonedimensionalproblem.WiththerapidadvancesofcomputersinbothCPUandmemorycapacity,andbecauseoftheavailabilityof3Dmodelsforturbulentflowsimulationsinpracticalapplications,thefully3Dsimulationofsuspendedsedimentbecomesafeasibleandurgenttaskforhydraulicengineers.TheprocessofsuspendedsedimenttransportisdescribedbytheadVection-d…     

对流弥散方程不同有限元处理方法的比较分析

Galerkin有限元在处理含第二类边界条件的对流弥散方程时,针对对流项和弥散项有两种不同的格林积分变换,所得数值结果的精度也不同。一种方法是把对流和弥散项整体考虑实施格林积分转换(降低微分阶数,由二阶降成一阶),应用边界条件,得出变分方程;另一种处理方法是只对弥散项实施积分变换,应用边界条件,得出变分方程。以一维问题为参考,对两种方法的数值结果与解析解进行比较分析。