Comparison of theory and experiment for solute transport in weakly heterogeneous bimodal porous media

In this work, the influence of non-equilibrium effects on solute transport in a weakly heterogeneous medium is discussed. Three macro-scale models (upscaled via the volume averaging technique) are investigated: (i) the two-equation non-equilibrium model, (ii) the one-equation asymptotic model and (iii) the one-equation local equilibrium model. The relevance of each of these models to the experimental system conditions (duration of the pulse injection, dispersivity values…) is analyzed. The numerical results predicted by these macroscale models are compared directly with the experimental data (breakthrough curves). Our results suggest that the preasymptotic zone (for which a non-Fickian model is required) increases as the solute input pulse time decreases. Beyond this limit, the asymptotic regime is recovered. A comparison with the results issued from the stochastic theory for this regime is performed. Results predicted by both approaches (volume averaging method and stochastic analysis) are found to be consistent.     

Modelling the tidal mixing fronts and seasonal stratification of the Northwest European Continental shelf

We investigate mixing processes under stratified conditions on the Northwest European Continental shelf using a numerical model (POLCOMS). Our results indicate that convection induced by vertical shearing of horizontal density gradients (‘shear-induced convection’) is a regularly occurring feature in the bottom and surface boundary layers in this open shelf-sea situation. Two types of turbulence models are investigated to study their capability for reproducing the observed location of tidal mixing fronts, and the physical processes occurring in seasonally stratified waters. The first model is a one-equation variant of the Mellor–Yamada model, whereas the second model combines a more recent second-momentum closure with a two-equation model. It is found that generally mean frontal positions (as estimated from ICES data) are predicted more accurately by the two-equation model. The one-equation model reproduces the mean frontal locations to 18.1 km (<3 grid spacings) and the two-equation model to 17.1 km; although in the Celtic Sea the accuracy is ∼33 and ∼12 km, respectively. Comparison with historical tide gauges, current metres, CTD stations, and thermistor chain data from the North Sea Project all show an improvement in accuracy when the two-equation model is used. This is particularly apparent in the model's ability to reproduce the spring–neap variability during stratification. We find that in the presence of shear-induced convection the routinely applied clipping of the turbulent length-scale, previously thought to be a minor ingredient in a turbulence model, has a dramatic effect on the results: if the length-scale clipping is not applied, substantial over-mixing is observed to occur. The causes and possible remedies of this effect are investigated. Overall our results demonstrate a sensitivity to the details of the turbulence model that is significantly greater than previously thought.     

基于高阶双渐近透射边界的重力坝-层状地基动力相互作用分析

时域高阶双渐近透射边界能够同时模拟层状介质中行波和快衰波的传播,具有很高的计算精度和计算效率.本文将高阶双渐近透射边界推广应用到多层层状地基系统弹性波传播问题的模拟,采用广义特征值分解分析该透射边界的数值稳定性,通过移谱法消除导致数值不稳定的虚假模态.将高阶双渐近透射边界以超单元的形式直接嵌入到近场有限元方程,建立了有限元-高阶双渐近透射边界时域耦合分析模型,并将其应用到重力坝-层状地基动力相互作用分析.数值算例分析结果表明,该时域耦合分析模型具有很高的精度和计算效率,适用于实际重力坝工程的地震响应分析.     

Modelling of solute and water transport in semi-permeable clay membranes: comparison with experiments

Theories of osmosis in groundwater flow are increasingly used to explain anomalies of salinity in clayey environments. However, predictive modelling through mathematical analysis can hardly be found in literature. In this paper, a model of chemical osmosis based on non-equilibrium thermodynamics, is used to predict the evolution of pressure and salinity in a clay membrane. Analysis of this model reveals simplifications that hold for specific situations. Two experiments from literature serve to show that the analytical modelling solution agrees with numerical and experimental results. Moreover, it is shown that the commonly applied Boussinesq approximation necessarily does not hold when osmosis is involved. Indeed, the clay system must be able to store the excess flow of water induced by osmosis.     

Modeling non-equilibrium mass transport in biologically reactive porous media

We develop a one-equation non-equilibrium model to describe the Darcy-scale transport of a solute undergoing biodegradation in porous media. Most of the mathematical models that describe the macroscale transport in such systems have been developed intuitively on the basis of simple conceptual schemes. There are two problems with such a heuristic analysis. First, it is unclear how much information these models are able to capture; that is, it is not clear what the model's domain of validity is. Second, there is no obvious connection between the macroscale effective parameters and the microscopic processes and parameters. As an alternative, a number of upscaling techniques have been developed to derive the appropriate macroscale equations that are used to describe mass transport and reactions in multiphase media. These approaches have been adapted to the problem of biodegradation in porous media with biofilms, but most of the work has focused on systems that are restricted to small concentration gradients at the microscale. This assumption, referred to as the local mass equilibrium approximation, generally has constraints that are overly restrictive. In this article, we devise a model that does not require the assumption of local mass equilibrium to be valid. In this approach, one instead requires only that, at sufficiently long times, anomalous behaviors of the third and higher spatial moments can be neglected; this, in turn, implies that the macroscopic model is well represented by a convection–dispersion–reaction type equation. This strategy is very much in the spirit of the developments for Taylor dispersion presented by Aris (1956). On the basis of our numerical results, we carefully describe the domain of validity of the model and show that the time-asymptotic constraint may be adhered to even for systems that are not at local mass equilibrium.     

Study of a smart platform based on backstepping control method

A structural model is significant for the verification of structural control algorithms. However, for nonlinear behavior, experiments are mostly destructive tests that are costly, and conducting repetitive structural experiments is difficult. Therefore, a repetitive structural vibration model is important for structural vibration control. In this study, a smart platform to realize different structural behaviors is developed based on the backstepping control algorithm. Lyapunov functions are used to derive the control law. Simulations show that the designed model can track the structural responses of different arbitrary linear structures very well. In addition, the proposed platform can track responses of different piecewise linear structures and desired models with various hysteresis very well. Numerical results verify the effectiveness of the proposed tracking controller through the backstepping method for the established platform.     

三维热导方程求解的波数域方法(MWN3D)及模型实验研究

给出了三维热导方程波数域的基本解和算法（MWN3D）,并通过数值模型的方式验证、讨论了这一方法的可靠性和适应能力. 模型实验表明这一方法具有计算速度快,处理问题灵活的特点. 模拟实验还显示,只有在横向较为均匀的三维构造条件下才能用一维模型替代三维模型. 事实上,由于在地壳和上地幔中物质的横向变化不可忽略,通过本方法构建一个真三维的热结构模型是必要的.     

EFFECT OF COARSENING OF SURFACE BED MATERIAL ON NON-EQUILIBRIUM SEDIMENT TRANSPORT PROCESS DURING RIVER DEGRADATION

1INTRODUCTIONNon-equilibriumsedimenttransportina"at'Uralstreamwithnon-uniformbedmaterialisasubjectofilltensiveresearch.ThesedimentconcentrationmaybenotequaltothetransportcapacityofflowinanalluvialrivedItgraduallyapproachesequilibriumbydegradationoraggradationinalongdistance.usuallyover100kilometers.Thisproblemisofgreatimportancetothepredictionofthedistributionoferosionordepositionalongachannel,especiallyforlargeriversinChina.Scientistsstudiednonequilibriumsedimenttransportandthemainresul…     

一种新的针对反射波的逆时偏移真振幅成像条件

真振幅成像是一种代表性的定量估计模型参数扰动高波数部分的地震波成像方法.经典的真振幅成像方法在高频近似和理想照明假设条件下求取显式对角Hessian逆矩阵作为偏移振幅加权算子,用以校正波传播过程中的几何扩散效应,得到模型参数扰动的带限估计.真振幅保真成像方法在利用逆时偏移（RTM）框架实现时会产生低波数噪声,影响对高波数参数估计的精度.本文给出了一种新的基于RTM框架的真振幅保真成像条件,该成像条件针对反射波数据,在高频近似下散射模式对应正问题及Bayes反问题框架下导出.与传统基于高频渐进反演的波动方程成像方法类似,利用本文提出RTM成像条件能够保证计算结果与高频近似下反演结果的一致性.同时,利用本文提出RTM真振幅成像条件能够在成像过程中自动保真的消除传统真振幅RTM算法中存在低波数噪声,模型数值实验结果验证了本文方法的正确性和有效性.     

Transport in chemically and mechanically heterogeneous porous media. I: Theoretical development of region-averaged equations for slightly compressible single-phase flow

Transport processes in heterogeneous porous media are often treated in terms of one-equation models. Such treatment assumes that the velocity, pressure, temperature, and concentration can be represented in terms of a single large-scale averaged quantity in regions having significantly different mechanical, thermal, and chemical properties. In this paper we explore the process of single-phase flow in a two-region model of heterogeneous porous media. The region-averaged equations are developed for the case of a slightly compressible flow which is an accurate representation for a certain class of liquid-phase flows. The analysis leads to a pair of transport equations for the region averaged pressures that are coupled through a classic exchange term, in addition to being coupled by a diffusive cross effect. The domain of validity of the theory has been identified in terms of a series of length and timescale constraints.In Part II the theory is tested, in the absence of adjustable parameters, by comparison with numerical experiments for transient, slightly compressible flow in both stratified and nodular models of heterogeneous porous media. Good agreement between theory and experiment is obtained for nodular and stratified systems, and effective transport coefficients for a wide range of conditions are presented on the basis of solutions of the three closure problems that appear in the theory. Part III of this paper deals with the principle of large-scale mechanical equilibrium and the region-averaged form of Darcy's law. This form is necessary for the development and solution of the region-averaged solute transport equations that are presented in Part IV. Finally, in Part V we present results for the dispersion tensors and the exchange coefficient associated with the two-region model of solute transport with adsorption.     

Determination of recovery factor for simulation of non-equilibrium sedimentation in reservoir

It is generally acceptable to assume that bed material load is equal to sediment transport capacity, if the exchange between bed and flow occurs instantaneously. However, for non-equilibrium sediment t...     

Groundwater flow to a system of drainage canals

A classical model of groundwater flow to a system of drainage canals is considered within a framework of a two-dimensional steady-state problem. The solution to this problem is derived based on the Riemann-Schwarts principle of symmetry. This solution yields simple analytical relationships expressed in terms of special or elementary functions. Numerical calculations are used to analyze in detail the effect of all physical characteristics of the model on the flow pattern. In particular, it was established that the presence of water in channels has a significant effect on the flow regime. The limiting cases of the scheme are considered, and simple approximated formulas are derived in these cases for the flow and discharge components.     

A probabilistic approach to structural model updating

The problem of updating a structural model and its associated uncertainties by utilizing measured dynamic response data is addressed. A Bayesian probabilistic formulation is followed to obtain the posterior probability density function (PDF) of the uncertain model parameters for given measured data. The present paper discusses the issue of identifiability of the model parameters and reviews existing asymptotic approximations for identifiable cases. The focus of the paper is on the treatment of the general unidentifiable case where the earlier approximations are not applicable. In this case the posterior PDF of the parameters is found to be concentrated in the neighborhood of an extended and extremely complex manifold in the parameter space. The computational difficulties associated with calculating the posterior PDF in such cases are discussed and an algorithm for an efficient approximate representation of the above manifold and the posterior PDF is presented. Numerical examples involving noisy data are presented to demonstrate the concepts and the proposed method.     

The use of a spatially heterogeneous simulation model for studying biotransformation processes of nitrogen and phosphorus compounds and the dynamics of oxygen dissolved in water in the ecosystem of Neva Bay, the Gulf of Finland: 2. Input data for calculations, modeling results, and their analysis

The data collected for model calculations is systematized for natural and design conditions. The natural conditions did not take into account the effect of the complex of water-protection structures on the hydrological regime of Neva Bay and on the biogenic load onto this water area, while the design conditions, conversely, reflected the possible impact of those structures on the hydrology and ecology of the area. Numerical experiments were used to study the processes of transformation of N and P compounds and the dynamics of dissolved O₂ in Neva Bay water area. In the comparison of the calculated and observed concentrations of biogenic substances, Theil criterion was evaluated to assess the adequacy of the model in reproducing the concentration fields of the distribution of biogenic substances over Neva Bay water area. The major qualitative and quantitative features of the formation of the spatial heterogeneity and the time variations in the concentrations of biogenic element compounds over Neva Bay water area are identified. Possible improvements of the model in the reproduction of the complex of processes that are of particular importance for the development of substance transformations in shallow ecosystems are considered.     

Modelling dam-break flows over mobile beds using a 2D coupled approach

Dam-break flows usually propagate along rivers and floodplains, where the processes of fluid flow, sediment transport and bed evolution are closely linked. However, the majority of existing two-dimensional (2D) models used to simulate dam-break flows are only applicable to fixed beds. Details are given in this paper of the development of a 2D morphodynamic model for predicting dam-break flows over mobile beds. In this model, the common 2D shallow water equations are modified, so that the effects of sediment concentrations and bed evolution on the flood wave propagation can be considered. These equations are used together with the non-equilibrium transport equations for graded sediments and the equation of bed evolution. The governing equations are solved using a matrix method, thus the hydrodynamic, sediment transport and morphological processes can be jointly solved. The model employs an unstructured finite volume algorithm, with an approximate Riemann solver, based on the Roe-MUSCL scheme. A predictor–corrector scheme is used in time stepping, leading to a second-order accurate solution in both time and space. In addition, the model considers the adjustment process of bed material composition during the morphological evolution process. The model was first verified against results from existing numerical models and laboratory experiments. It was then used to simulate dam-break flows over a fixed bed and a mobile bed to examine the differences in the predicted flood wave speed and depth. The effects of bed material size distributions on the flood flow and bed evolution were also investigated. The results indicate that there is a great difference between the dam-break flow predictions made over a fixed bed and a mobile bed. At the initial stage of a dam-break flow, the rate of bed evolution could be comparable to that of water depth change. Therefore, it is often necessary to employ the turbid water governing equations using a coupled approach for simulating dam-break flows.     

An asymptotic model for resistive instability in the Earth's outer core

Abstract

Numerical work indicates that resistive instability may be the dominant mode of instability in the Earth's outer core for realistic core parameter regimes. In this paper, we assume that the Elsasser number is large in order to obtain an asymptotic analysis of resistive instability in an electrically conducting fluid confined to a rotating cylindrical shell of infinite extent in the axial direction. The dimensionless equations of motion are linearized about an ambient magnetic field which is purely azimuthal and depends only on the cylindrical radial variable. Applying the theory of ordinary differential equations with a large parameter, we obtain an asymptotic approximation to the solution. Relatively simple analytic expressions for the complex frequencies are obtained by applying the boundary conditions for insulating boundaries at the cylindrical sidewalls and then assuming that the ambient magnetic field vanishes at one or both of those sidewalls. The results appear to be consistent with previous numerical work.