SD-WACCM模式对平流层化学组分的模拟研究

利用美国大气研究中心开发的全球气候模式,对2008年平流层的化学组分(HNO₃、HCl和O₃)进行了模拟研究,并结合了MLS卫星资料进行了对比分析。结果表明,模式可以较好地再现平流层的各化学组分的时空分布状况。并选取了8个区域,将模拟的2008—2009年O₃柱浓度与臭氧监测仪资料对比,结果表明,模式可以较好地再现全球O₃总量的季节变化情况。     

Finite deformation analysis of geomaterials

The mathematical structure and numerical analysis of classical small deformation elasto–plasticity is generally well established. However, development of large deformation elastic–plastic numerical formulation for dilatant, pressure sensitive material models is still a research area. In this paper we present development of the finite element formulation and implementation for large deformation, elastic–plastic analysis of geomaterials. Our developments are based on the multiplicative decomposition of the deformation gradient into elastic and plastic parts. A consistent linearization of the right deformation tensor together with the Newton method at the constitutive and global levels leads toward an efficient and robust numerical algorithm. The presented numerical formulation is capable of accurately modelling dilatant, pressure sensitive isotropic and anisotropic geomaterials subjected to large deformations. In particular, the formulation is capable of simulating the behaviour of geomaterials in which eigentriads of stress and strain do not coincide during the loading process. The algorithm is tested in conjunction with the novel hyperelasto–plastic model termed the B material model, which is a single surface (single yield surface, affine single ultimate surface and affine single potential surface) model for dilatant, pressure sensitive, hardening and softening geomaterials. It is specifically developed to model large deformation hyperelasto–plastic problems in geomechanics. We present an application of this formulation to numerical analysis of low confinement tests on cohesionless granular soil specimens recently performed in a SPACEHAB module aboard the Space Shuttle during the STS‐89 mission. We compare numerical modelling with test results and show the significance of added confinement by the thin hyperelastic latex membrane undergoing large stretching. Copyright © 2001 John Wiley & Sons, Ltd.     

‘Surface Drag in the Arctic Marginal Sea-ice Zone: A Comparison of Different Parameterisation Concepts’

Two parameterisation schemes for the turbulent surface fluxes and drag coefficients over the Arctic marginal sea-ice zone (MIZ) are (further) developed, and their results are compared with each other. Although the schemes are based on different principles (flux averaging and parameter averaging), the resulting drag coefficients differ only slightly in the case of neutral and stable stratification. For unstable stratification and sea-ice conditions being typical for the north-eastern Fram Strait, the drag coefficients resulting from the parameter-averaging concept are 5–10% larger than those of the flux-averaging concept. At a sea-ice concentration of 45%, the parameter-averaging method overestimates the heat fluxes by a factor of 1.2. An inclusion in the schemes of form drag caused by floe edges and ridges has a much larger effect on the drag coefficient, and on the momentum fluxes, than the choice between the parameter-averaging or flux-averaging methods. Based on sensitivity studies with the flux-averaging scheme, a simple formula for the effective drag coefficient above the Arctic MIZ is derived. It reduces the computational costs of the more complex parameterisations and could also be used in larger scale models. With this simple formula, the effective drag coefficient can be calculated as a function of the sea-ice concentration and skin drag coefficients for water and ice floes. The results obtained with this parameterisation differ only slightly from those using the more complex schemes. Finally, it is shown that in the MIZ, drag coefficients for sea-ice models may differ significantly from the effective drag coefficients used in atmospheric models.     

The effect of random recharge on uniform steady free-surface flow in heterogeneous porous formations

The effect of parametric uncertainty in recharge rate and spatial variability of hydraulic conductivity upon free-surface flow is investigated in a stochastic framework. We examine the three-dimensional free-surface gravitational flow problem for sloped mean uniform flow in a randomly heterogeneous porous medium under the influence of random recharge. We develop analytic solutions for the variance of free-surface position, head, and specific discharge on the free surface. Additionally, we obtain semi-analytic solutions for the statistical moments of head and specific discharge beneath the free-surface. Statistical moments are derived using a first-order approximation and then compared with their parallel in an unbounded medium. The effect of recharge mean and variability on the statistical moments is analyzed. Results can be applied to more complex flows, slowly varying in the mean.     

Stochastic analysis of spatial variability in unconfined groundwater flow

In this paper, spatial variability in steady one-dimensional unconfined groundwater flow in heterogeneous formations is investigated. An approach to deriving the variance of the hydraulic head is developed using the nonlinear filter theory. The nonlinear governing equation describing the one-dimensional unconfined groundwater flow is decomposed into three linear partial differential equations using the perturbation method. The linear and quadratic frequency response functions are obtained from the first- and second-order perturbation equations using the spectral method. Furthermore, under the assumption of the exponential covariance function of log hydraulic conductivity, the analytical solutions of both the spectrum and the variance of the hydraulic head produced from the linear system are derived. The results show that the variance derived herein is less than that of Gelhar (1977). The reason is that the log transmissivity is linearized in Gelhars work. In addition, the analytical solutions of both the spectrum and the variance of the hydraulic head produced from the quadratic system are derived as well. It is found that the correlation scale and the trend in mean of log hydraulic conductivity are important to the dimensionless variance ratio.     

Bias and uncertainty in regression-calibrated models of groundwater flow in heterogeneous media

Groundwater models need to account for detailed but generally unknown spatial variability (heterogeneity) of the hydrogeologic model inputs. To address this problem we replace the large, m-dimensional stochastic vector β that reflects both small and large scales of heterogeneity in the inputs by a lumped or smoothed m-dimensional approximation γθ_∗, where γ is an interpolation matrix and θ_∗ is a stochastic vector of parameters. Vector θ_∗ has small enough dimension to allow its estimation with the available data. The consequence of the replacement is that model function f(γθ_∗) written in terms of the approximate inputs is in error with respect to the same model function written in terms of β, f(β), which is assumed to be nearly exact. The difference f(β) − f(γθ_∗), termed model error, is spatially correlated, generates prediction biases, and causes standard confidence and prediction intervals to be too small. Model error is accounted for in the weighted nonlinear regression methodology developed to estimate θ_∗ and assess model uncertainties by incorporating the second-moment matrix of the model errors into the weight matrix. Techniques developed by statisticians to analyze classical nonlinear regression methods are extended to analyze the revised method. The analysis develops analytical expressions for bias terms reflecting the interaction of model nonlinearity and model error, for correction factors needed to adjust the sizes of confidence and prediction intervals for this interaction, and for correction factors needed to adjust the sizes of confidence and prediction intervals for possible use of a diagonal weight matrix in place of the correct one. If terms expressing the degree of intrinsic nonlinearity for f(β) and f(γθ_∗) are small, then most of the biases are small and the correction factors are reduced in magnitude. Biases, correction factors, and confidence and prediction intervals were obtained for a test problem for which model error is large to test robustness of the methodology. Numerical results conform with the theoretical analysis.     

Application of Dynamic Subgrid-scale Models for Large-eddy Simulation of the Daytime Convective Boundary Layer over Heterogeneous Surfaces

The sensitivity of large-eddy simulation (LES) to the representation of subgrid-scale (SGS) processes is explored for the case of the convective boundary layer (CBL) developing over surfaces with varying degrees of spatial heterogeneity. Three representations of SGS processes are explored: the traditional constant Smagorinsky–Lilly model and two other dynamic models with Lagrangian averaging approaches to calculate the Smagorinsky coefficient (C _S ) and SGS Prandtl number (Pr). With initial data based roughly on the observed meteorology, simulations of daytime CBL growth are performed over surfaces with characteristics (i.e. fluxes and roughness) ranging from homogeneous, to striped heterogeneity, to a realistic representation of heterogeneity as derived from a recent field study. In both idealized tests and the realistic case, SGS sensitivities are mostly manifest near the surface and entrainment zone. However, unlike simulations over complex domains or under neutral or stable conditions, these differences for the CBL simulation, where large eddies dominate, are not significant enough to distinguish the performance of the different SGS models, irrespective of surface heterogeneity.     

The fast Ice Nucleus chamber FINCH

We present first results of our new developed Ice Nucleus (IN) counter FINCH from the sixth Cloud and Aerosol Characterization Experiment (CLACE 6) campaign at Jungfraujoch station, 3571 m asl. Measurements were made at the total and the ICE CVI inlet. Laboratory measurements of ice onset temperatures by FINCH are compared to those of the static diffusion chamber FRIDGE (FRankfurt Ice Deposition Freezing Experiment). Within the errors of both new instruments the results compare well to published data.     

Integration schemes with substepping algorithms for creep analysis in geomaterials

This paper describes a robust and efficient methodology for predicting displacements, deformations, and stresses in geomaterials that are susceptible to creep. The methodology is based on two integration schemes, which consider substepping algorithms. The first scheme is used for integrating space-time relations in a global sense, whereas the second scheme is used for integrating stress-stain relations in a local sense. Different from previous studies, both integration schemes are easy to implement and general in the sense that they can be applied to any type of creep law. Through an in-house finite element simulator, several numerical tests are performed. They include triaxial and wellbore closure analyses considering soft soils and salt rocks. The results show that the combination of both schemes leads to stable and accurate solutions with reduced computational time.     

卫星遥感结合地面观测估算非均匀地表区域能量通量

卫星遥感在研究非均匀陆面上地－气间能量和水循环过程时有其独到的作用。文中介绍了一种利用ＬＡＮＤＳＡＴＴＭ资料进行非均匀陆面上区域能量平衡研究的参数化方案。并以夏季和深秋两个景的ＴＭ资料为个例,结合“黑河实验”（ＨＥＩＦＥ）期间的地面观测资料分析了实验区非均匀陆面上地表特征参数及能量平衡各分量的区域分布及季节差异,同时将所得的结果与地面观测的“真值”作了比较,所得的这些结果是基本可信的。这种方法仍然处于发展的初始阶段,对此也进行了讨论。