共查询到20条相似文献，搜索用时 31 毫秒
1.
《Applied Geochemistry》2001,16(14):1653-1662
Strontium-90 (90Sr) is one of the major radioactive contaminants found in DP Canyon at Los Alamos, New Mexico, USA. Radioactive surveys found that 90Sr is present in surface water and shallow alluvial groundwater environments in Los Alamos National laboratory (LANL). Colloids may influence the transport of this radionuclide in surface and groundwater environments in LANL. In this study, the authors investigated the sorption/desorption behavior of radioactive Sr on Ca-montmorillonite and silica colloids, and the effect of ionic strength of water on the sorption of Sr. Laboratory batch sorption experiments were conducted using 85Sr as a surrogate for 90Sr. Groundwater, collected from Well LAUZ-1 at DP Canyon and from Well J-13 at Yucca Mountain, Nevada, and deionized water, were used. The results show that 92–100% of the 85Sr was rapidly adsorbed onto Ca-montmorillonite colloids in all three waters. Adsorption of 85Sr onto silica colloids varied among the three waters. The ionic strength and Ca2+ concentration in groundwater significantly influence the adsorption of 85Sr onto silica colloids. Desorption of 85Sr from Ca-montmorillonite colloids is slower than from silica colloids. Desorption of 85Sr from silica colloids was faster in LAUZ-1 groundwater than in J-13 groundwater and deionized water. The results suggest that clay and silica colloids may facilitate the transport of Sr along potential flowpaths from DP Canyon to Los Alamos Canyon. 相似文献
2.
A primal C0-conforming virtual element discretization for the approximation of the bidimensional two-phase flow of immiscible fluids in porous media using general polygonal meshes is discussed. This work investigates the potentialities of the Virtual Element Method (VEM) in solving this specific problem of immiscible fluids in porous media involving a time-dependent coupled system of non-linear partial differential equations. The performance of the fully discrete scheme is thoroughly analysed testing it on general meshes considering both a regular problem and more realistic benchmark problems that are of interest for physical and engineering applications.
相似文献3.
Ernest S. GLADNEY David B. CURTIS Daniel R. PERRIN 《Geostandards and Geoanalytical Research》1984,8(1):43-46
Concentrations of boron in 35 silicate reference materials are measured by thermal neutron capture gamma-ray spectrometry. Results are compared for NBS, USGS, and CCRMP reference materials with values from the literature. The use of two prompt gamma facilities at Los Alamos are discussed. 相似文献
4.
Investigation of a hybrid polygonal finite element formulation for confined and unconfined seepage 
下载免费PDF全文
下载免费PDF全文 This study presents a formulation for field problems using hybrid polygonal finite elements, taking steady state seepage through a porous material as the focus. We make comparisons with a conventional finite element formulation based on a single primary variable, focussing on the advantages of the hybrid formulation in terms of flux field accuracy and extension to convex polygonal shaped elements. For the unconfined case, we adopt a head dependent hydraulic conductivity that does not require remeshing. The performance of the hybrid polygonal element formulation is demonstrated through a series of numerical examples. The results show a sensitivity of the location of the free surface in unconfined seepage to mesh configuration for hybrid quadrilateral meshes with various aspect ratios, but not for hybrid polygonal meshes with various orientations and irregularity. Examination of the free surface location results for several conforming shape function options shows an insensitivity to choice of interpolation function, provided that it conforms with the assumptions in the formulation. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
5.
6.
The Pajarito fault forms the western margin of the Rio Grande rift in north-central New Mexico, and lies adjacent to Los Alamos National Laboratory, a major Federal research facility. Vertical displacement on this normal fault over the past 1.2 Ma has created a 50- to 120-m-high fault scarp on Bandelier Tuff (1.2 Ma), yielding a long-term average slip rate of ca. 0.1 mm/yr. In support of a Laboratory-wide seismic hazards assessment, we excavated 14 trenches in the Pajarito fault zone to determine the age of the most recent displacement event, the recurrence interval between events, the displacement per event, and the variability in slip rate and recurrence through time. The large number of trenches was required by the large height of the fault scarp and the complexity of the fault zone. Only about half the trenches contained significant thicknesses of Holocene deposits, but in those trenches there was clear evidence for an early-to-mid-Holocene displacement event. The previous event was at least 20–40 ka, and the average recurrence interval over the past ca. 300 ka was about 20–40 kyr. We infer that much of the structural relief across this fault developed soon after eruption of the Bandelier Tuff between 1.0 and 1.2 Ma, and that slip rate slowed considerably after that time. 相似文献
7.
The Los Alamos National Laboratory conducted a reconnaissance-scale geochemical exploration program on the island of St. Lucia during September 1983. Totals of 307 stream-sediment and 55 beach-sand samples were collected and analyzed for a suite of 44 elements. Data interpretation utilized computer-generated color images. Several anomalies were delineated, the most significant of which indicate previously unrecognized Au mineralization in two areas. 相似文献
8.
9.
Najla Frih Vincent Martin Jean Elizabeth Roberts Ali Saada 《Computational Geosciences》2012,16(4):1043-1060
We consider a model for fluid flow in a porous medium with a fracture. In this model, the fracture is treated as an interface between subdomains, on which specific equations have to be solved. In this article, we analyze the discrete problem, assuming that the fracture mesh and the subdomain meshes are completely independent, but that the geometry of the fracture is respected. We show that despite this nonconformity, first-order convergence is preserved with the lowest-order Raviart–Thomas(-Nedelec) mixed finite elements. Numerical simulations confirm this result. 相似文献
10.
11.
冬季在北白令海陆架区域频繁地出现潜热冰间湖,对当地的生态系统和北极盐跃层贡献很大.将CICE海冰模式应用到该区域,采用高分辨率(6.37 km)网格,模拟2002年11月至2003年4月的海冰变化过程,模拟的海冰总面积和海冰密集度与AMSR-E/Aqua卫星遥感结果吻合很好,其中两者日平均海冰总面积在模拟期间的相关系数达到0.97.模拟结果表明,东北风将海冰向南输运在东西走向的海岸南部形成冰间湖,反映了潜热冰间湖形成和演化的动力过程.对卫星观测数据,将海冰密集度<75%作为冰间湖的判据;而对数值模拟结果,确定海冰密集度<70%为冰间湖的判据.据此讨论白令海4个区域的冰间湖形成过程,与卫星数据进行比较,大部分冰间湖得到很好的模拟.深入讨论了影响冰间湖模拟准确度的主要因素,认为选用恰当的阈值、提高气象强迫场的空间和时间分辨率有助于提高模拟效果.对部分海域的冰间湖模拟效果不佳,需要发展冰海耦合模式才能最终解决. 相似文献
12.
Richards方程常用于非饱和土渗流问题,并且应用广泛。在数值求解中,对Richards方程线性化,进而采用有限差分法进行数值离散以及迭代计算。其中传统的迭代法比如Jacobi迭代、Gauss-Seidel迭代法（GS）和连续超松驰迭代法（successive over-relaxation method,简称SOR）迭代收敛率较慢,尤其在离散空间步长较小以及离散时间步长较大时。因此,采用整体校正法以及多步预处理法对传统迭代法进行改进,提出一种基于整体校正法的多步预处理Gauss-Seidel迭代法（improved Gauss-Seidel iterative method with multistep preconditioner based on the integral correction method,简称ICMP(m)-GS）求解Richards方程导出的线性方程组。通过非饱和渗流算例,并与传统迭代法和解析解对比,对改进算法的收敛率和加速效果进行了验证。结果表明,提出的ICMP(m)-GS可以很大程度地改善线性方程组的病态性,相较于常规方法GS,SOR以及单一改进方法,ICMP(m)-GS具有更快的收敛率,更高的计算效率和计算精度。该方法可以为非饱和土渗流的数值模拟提供一定参考。 相似文献
13.
A general approach for obtaining the consistent tangent operator for constitutive rate equations is presented. The rate equations can be solved numerically by the user's favourite time integrator. In order to obtain reliable results, the substepping in integration should be based on a control of the local error. The main ingredient of the consistent tangent operator, namely the derivative of the stress with respect to the strain increment must be computed simultaneously with the same integrator, applied to a numerical approximation of the variational equations. This information enables finite‐element packages to assemble a consistent tangent operator and thus guarantees quadratic convergence of the equilibrium iterations. Several numerical examples with a hypoplastic constitutive law are given. As numerical integrator we used a second‐order extrapolated Euler method. Quadratic convergence of the equilibrium iteration is shown. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
14.
Numerical modeling of salinity changes in marine environment of Persian Gulf is investigated in this paper. Computer simulation of the problem is performed by the solution of a convection-diffusion equation for salinity concentration coupled with the hydrodynamics equations. The hydrodynamic equations consist of shallow water equations of continuity and motion in horizontal plane. The effects of rain and evaporations are considered in the continuity equation and the effects of bed slope and friction, as well as the Coriolis effects are considered in two equations of motion. The cell vertex finite volume method is applied for solving the governing equations on triangular unstructured meshes. Using unstructured meshes provides great flexibility for modeling the flow problems in arbitrary and complex geometries, such as Persian Gulf domain. The results of evaporation and Coriolis effects, as well as imposing river and tidal boundary conditions to the hydrodynamic model of Persian Gulf (considering variable topology rough bed) are compared with predictions of Admiralty Tide Table, which are obtained from the harmonic analysis. The performance of the developed computer model is demonstrated by simulation of salinity changes due to inflow effects and diffusion effects as well as computed currents. 相似文献
15.
This paper presents a novel mass-conservative mixed multiscale method for solving flow equations in heterogeneous porous media. The media properties (the permeability) contain multiple scales and high contrast. The proposed method solves the flow equation in a mixed formulation on a coarse grid by constructing multiscale basis functions. The resulting velocity field is mass-conservative on the fine grid. Our main goal is to obtain first-order convergence in terms of the mesh size which is independent of local contrast. This is achieved, first, by constructing some auxiliary spaces, which contain global information that cannot be localized, in general. This is built on our previous work on the generalized multiscale finite element method (GMsFEM). In the auxiliary space, multiscale basis functions corresponding to small (contrast-dependent) eigenvalues are selected. These basis functions represent the high-conductivity channels (which connect the boundaries of a coarse block). Next, we solve local problems to construct multiscale basis functions for the velocity field. These local problems are formulated in the oversampled domain, taking into account some constraints with respect to auxiliary spaces. The latter allows fast spatial decay of local solutions and, thus, allows taking smaller oversampled regions. The number of basis functions depends on small eigenvalues of the local spectral problems. Moreover, multiscale pressure basis functions are needed in constructing the velocity space. Our multiscale spaces have a minimal dimension, which is needed to avoid contrast dependence in the convergence. The method’s convergence requires an oversampling of several layers. We present an analysis of our approach. Our numerical results confirm that the convergence rate is first order with respect to the mesh size and independent of the contrast. 相似文献
16.
The Geochemical Evolution of Anorthosite Residual Magmas in the Laramie Anorthosite Complex, Wyoming 总被引：6，自引：4，他引：2
MITCHELL JEREMY N.; SCOATES JAMES S.; FROST CAROL D.; KOLKER ALLAN 《Journal of Petrology》1996,37(3):637-660
Olivine- and pyroxene-bearing Fe-enriched dioritic rocks inthe 1434 Ma Laramie anorthosite complex are interpreted to representvariably fractionated and contaminated magmas residual afterthe crystallization of anorthosite. Geochemical characteristicsof this suite include the following: high contents of TiO2,, and P2O5; high incompatibletrace element contents; rare earth element patterns with a largerange of Eu anomalies; and isotopic compositions that reflectthe geographic location of individual samples, with ISr increasingand Nd decreasing from south to north. After extraction fromanorthosite, fractionation of ferrodioritic residual magmasresulted in secondary residual monzodioritic melts and complementaryoxide-rich ferrodiorite cumulates. Geographic trends in isotopiccomposition reflect an increasing Archean crustal componentfrom south to north. Dioritic dikes and cumulates with isotopiccompositions similar to associated anorthosites were derivedlocally. Large isotopic discrepancies between some dioritesand their hosting anorthosites reflect preferential contaminationof residual magma during ascent and emplacement of mantle-derivedplagioclase-rich diapirs, followed by subsequent extractionand isolation of Fe-enriched interstitial melt. Strong isotopiccontrasts between anorthosite and associated Fe-enriched rocksin anorthosite complexes do not preclude a direct relationshipbetween them and reflect the diversity and complexity of processesduring their petrogenesis. KEY WORDS: anorthosite; ferrodiorite; geochemistry; Laramie anorthosite complex; residual magma
*Corresponding author: Present address: Ceramic Science and Technology, Mail Stop K762, Los Alamos National Laboratory, Los Alamos, NM 87545, USA. Phone: (505) 665-3934. Fax: (505) 665-3935. E-mail: jeremy{at}lanl.gov. 相似文献
17.
Helge Holden Kenneth Hvistendahl Karlsen Knut-Andreas Lie 《Computational Geosciences》2000,4(4):287-322
We present an accurate numerical method for a large class of scalar, strongly degenerate convection–diffusion equations. Important subclasses are hyperbolic conservation laws, porous medium type equations, two-phase reservoir flow equations, and strongly degenerate equations coming from the recent theory of sedimentation–consolidation processes. The method is based on splitting the convective and the diffusive terms. The nonlinear, convective part is solved using front tracking and dimensional splitting, while the nonlinear diffusion part is solved by an implicit–explicit finite difference scheme. In addition, one version of the implemented operator splitting method has a mechanism built in for detecting and correcting unphysical entropy loss, which may occur when the time step is large. This mechanism helps us gain a large time step ability for practical computations. A detailed convergence analysis of the operator splitting method was given in Part I. Here we present numerical experiments with the method for examples modelling secondary oil recovery and sedimentation–consolidation processes. We demonstrate that the splitting method resolves sharp gradients accurately, may use large time steps, has first order convergence, exhibits small grid orientation effects, has small mass balance errors, and is rather efficient. 相似文献
18.
Anne G. Hereford Elizabeth H. Keating George D. GuthrieJr Chen Zhu 《Environmental Geology》2007,52(5):965-977
Reactions and reaction rates within aquifers are fundamental components of critical hydrological processes. However, reactions
simulated in laboratory experiments typically demonstrate rates that are much faster than those observed in the field. Therefore,
it is necessary to conduct more reaction rate analyses in natural settings. This study of geochemical reactions in the regional
aquifer in the Pajarito Plateau near Los Alamos, New Mexico combines modeling with petrographic assessment to further knowledge
and understanding of complex natural hydrologic systems. Groundwater geochemistry shows marked evolution along assumed flow
paths. The flow path chosen for this study was evaluated using inverse mass balance modeling to calculate the mass transfer.
X-ray diffraction and field emission gun scanning electron microscopy were used to identify possible reactants and products.
Considering the mineralogy of the aquifer and saturation indices for the regional water refined initial interpretations. Calculations
yielded dissolution rates for plagioclase on the order of 10−15 mol s−1 m−2 and for K-feldspar on the order of 10−17 mol s−1 m−2, orders of magnitude slower than laboratory rates. While these rates agree with other aquifer studies, they must be considered
in the light of the uncertainty associated with geometric surface area estimates, 14C ages, and aquifer properties. 相似文献
19.
S. R. Sabbagh Yazdi 《International Journal of Environmental Science and Technology》2006,2(4):309-317
In this paper, simulation of oil spill due to tidal currents in Persian Gulf is performed by coupled solution of the hydrodynamics equations and an equation for convection and diffusion of the oil. The hydrodynamic equations utilized in this work consist of depth average equations of continuity and motion in two dimensional horizontal planes. The effect of evaporation is considered in the continuity equation and the effects of bed slope and friction, as well as the Coriolis effects are considered in two equations of motion. The overlapping cell vertex finite volume method is applied for solving the governing equations on triangular unstructured meshes. Using unstructured meshes provides great flexibility for modeling the flow in arbitrary and complex geometries, such as Persian Gulf flow domain. The results of the hydrodynamic model for tidal currents in Persian Gulf domain is examined by imposing tidal fluctuations to the main flow boundary during a limited period of time. Finally, the developed model is used to simulate an accidental oil spill from a point in Persian Gulf. 相似文献
20.
数值模拟是研究撞击坑形成过程的一种主要方法,尤其是认识撞击坑形成机制的重要手段。撞击坑形成过程数值模拟的基本原理是用离散方法描述物质在高速撞击作用下的运动及状态,在模拟中首先将物质与空间划分成离散的网格,在每一次迭代计算中逐步求解各个网格的形变、运动与状态改变的规律。牛顿运动定律、物质的连续体模型与热力学方程是撞击坑形成过程数值模拟的理论基础,牛顿运动定律以偏微分方程的形式贯穿在离散化的网格空间中,物质的连续体模型将物质的屈服强度与破裂、温度、孔隙、振动等联系起来,而热力学方程则通过其他热力学参数计算网格单元内物质的压强与物质所处的状态。 相似文献