Simulation of long-term transient density-dependent transport in groundwater     

Emil O. Frind 《Advances in water resources》1982,5(2):73-88

A numerical model for the economical simulation of long-term transient response in density-dependent transport problems is introduced. Although a classical Galerkin finite element approach is used, emphasis on optimum efficiency throughout the development results in a scheme that is found to be significantly less costly than comparable existing schemes. This advantage in efficiency increases the scope of simulation problems that can be handled within the constraints of a limited research budget. Some distinctive aspects are the elimination of static quantities in the fluid continuity equation, achieved by the introduction of equivalent freshwater head, and the elimination of numerical integration, achieved by the deliberate choice of linear elements. As a result of this choice, fluid velocities are discontinuous across the element boundaries. It is shown, however, that the solution obtained with discontinuous velocities approaches that obtained with continuous velocities as the grid is refined, and that the two types of solutions give essentially the same results when the elements are in the same size range. The model is applied to simulate the complete transient response for a well-known problem of seawater intrusion in a confined aquifer. The simulation is performed with both the constant dispersion coefficient used by previous researchers, and a more physically realistic velocity-dependent dispersion coefficient. Responses are found to be substantially different for the two types of coefficients, with the velocity-dependent dispersion coefficient producing much slower convergence to a state of dynamic equilibrium, and a much more pointed saltwater toe, which at the bottom of the aquifer tends to a sharp interface at equilibrium. Finally, it is shown by means of large-scale applications that the model is capable of efficiently simulating the long-term transient response in systems of practical significance.  相似文献   

Cenozoic alkali basaltic magmas of western Germany and their products of differentiation   总被引：11，自引：0，他引：11  

K. Hans Wedepohl  Emil Gohn  Gerald Hartmann 《Contributions to Mineralogy and Petrology》1994,115(3):253-278

Analytical data on major elements and 31 trace elements in olivine nephelinites, nepheline basanites, basanitic alkali olivine basalts and their differentiates (tephrites, hawaiites, mugearites, benmoreites, latites, phonolites and trachytes) from Hegau, Kaiserstuhl, Rhön, Hessian Depression, Vogelsberg, Westerwald, Siebengebirge, E Eifel and Hocheifel are evaluated. They were based on 400 samples with new or unpublished data on about one third of the rocks. The Sr–Nd isotopic compositions for 78 rocks are included. The alkali basaltic volcanism is caused by adiabatic decompression of asthenospheric mantle updomed to a minimum depth of 50 km in connection with the Alpine continent collision. The chemical compositions of the primary basaltic melts from the different areas are similar containing about one hundred-fold enrichment of highly incompatible elements relative to the primitive mantle from partial melting of depleted and secondarily enriched peridotite. The elements Cs, K, Pb and Ti are specifically depleted in the basalts partly because of phlogopite being residual at partial melting. The Tertiary alkali basalts range in Nd-isotopic composition from 0.51288 to 0.51273 and in Sr-isotopic ratios from 0.7032 to 0.7042. These ranges indicate mixtures of HIMU, depleted and enriched mantle components in the metasomatically altered peridotite source which resembles that of certain ocean islands. The Nd-Sr-isotopic compositions of the Quaternary E Eifel are close to bulk Earth ratios. East and W Eifel plots differ distinctly from the Tertiary Hocheifel which is geographically intermediate. This isotopic difference, beside specific K/Na ratios, is probably caused by separate metasomatic pulses that immediately preceded the respective periods of volcanism. The metasomatically altered mantle had partly primitive mantle signatures (Nb/Ta, Zr/Sm and Th/U ratios) and partly ocean island (or MORB) source properties (Rb/Cs). A MORB source can be excluded because of the low K/Rb and high Th/U ratios. A correlation of D with ⁸⁷Sr/⁸⁶Sr in amphibole and phlogopite and a slightly larger ¹⁸O than in MORB is conformable with a seawater and crustal impact on the source of alkali basalts. Slightly higher than average water concentrations in the source of certain primary basaltic melts (indicated by amphibole phenocrysts in their basalts) are required for differentiation of these basalts in magma chambers of the upper crust. Model calculations are presented to explain compositions of differentiates which range from about 60% to about 20% residual melt. The latter are represented by phonolites and trachytes. The Nd- and Sr-isotopic signatures of the majority of differentiates indicate contamination by a granitic partial melt from the wall rocks of magma chambers. Olivine nephelinite magma was the common source of contaminated differentiates.  相似文献   

Ein Versuch zur quantitativen Bestimmung der Faltungskraft     

Walther Emil Petrascheck Jr. 《International Journal of Earth Sciences》1950,38(2):133-138

  相似文献   

Machine Learning Can Assign Geologic Basin to Produced Water Samples Using Major Ion Geochemistry     

Shelton  Jenna L.  Jubb  Aaron M.  Saxe  Samuel W.  Attanasi  Emil D.  Milkov  Alexei V.  Engle  Mark  Freeman  Philip A.  Shaffer  Christopher A.  Blondes  Madalyn S. 《Natural Resources Research》2021,30(6):4147-4163

Natural Resources Research - Understanding the geochemistry of waters produced during petroleum extraction is essential to informing the best treatment and reuse options, which can potentially be...  相似文献   

The role of iron in tetrahedrite and tennantite determined by Rietveld refinement of neutron powder diffraction data     

Jens Wenzel Andreasen  Emil Makovicky  Bente Lebech  Sven Karup Møller 《Physics and Chemistry of Minerals》2008,35(8):447-454

Rietveld refinement of neutron powder diffraction data on four samples of synthetic, iron-bearing tetrahedrite (Cu_12?xFe_xSb₄S₁₃) with x = 0.28, 0.69, 0.91, 2.19 and four samples of synthetic tennantite (Cu_12?xFe_xAs₄S₁₃) with x = 0.33, 0.38, 0.86, 1.5 indicate unambiguously that iron is incorporated into tetrahedral M1 (12d) sites and not into triangular M2 (12e) sites in the cubic crystal structure (space group I $ \ifmmode\expandafter\bar\else\expandafter\=\fi{4} Rietveld refinement of neutron powder diffraction data on four samples of synthetic, iron-bearing tetrahedrite (Cu_12−xFe_xSb₄S₁₃) with x = 0.28, 0.69, 0.91, 2.19 and four samples of synthetic tennantite (Cu_12−xFe_xAs₄S₁₃) with x = 0.33, 0.38, 0.86, 1.5 indicate unambiguously that iron is incorporated into tetrahedral M1 (12d) sites and not into triangular M2 (12e) sites in the cubic crystal structure (space group I 3 m). The refinement results also confirm that M2 is a split (24g), flat-pyramidal site situated statistically on both sides of the S1−S1–S2 triangle. In tetrahedrite, this split is about 0.6 ?, in tennantite about 0.7 ?. Trends in bond lengths and magnitude of the M2 split were evaluated by means of linear regression with Fe concentration as the independent variable.  相似文献   

Coastal Ocean Forecasting Science supported by GODAE OceanView Coastal Oceans and Shelf Seas Task Team (COSS-TT)—Part II     

Cirano  Mauro  Charria  Guillaume  De Mey-Frémaux  Pierre  Kourafalou  Vassiliki H.  Stanev  Emil 《Ocean Dynamics》2021,71(6-7):779-783

  相似文献   

Tidal and wind-driven surface currents in the German Bight: HFR observations versus model simulations     

Alexander Port  Klaus-Werner Gurgel  Joanna Staneva  Johannes Schulz-Stellenfleth  Emil V. Stanev 《Ocean Dynamics》2011,61(10):1567-1585

Tidal and wind-driven surface currents in the German Bight between shallow mudflats of the North Frisian islands and the island of Helgoland are studied using coastal high-frequency radar (HFR) observations and hindcasts from a primitive equation numerical model. The setup of the observational system is described, and estimates of expected measurement errors are given. A quantitative comparison of numerical model results and observations is performed. The dominant tidal components are extracted from the two data sources using tidal harmonic analysis and the corresponding tidal ellipses are defined. Results show that the spatial patterns of different tidal ellipse parameters are consistent in the two data sets. Model sensitivity studies with constant and variable salinity and temperature distributions are used to study density-related mechanisms of circulation. Furthermore, the role of the surface wind field in driving the German Bight circulation is investigated using the complex correlation between wind and surface current vectors. The observed change of the respective correlation patterns from the coastal to open ocean is shown to be due to a combination of density effects, the coastline and topography. The overall conclusion is that HFR observations resolve the small-scale and rapidly evolving characteristics of coastal currents well in the studied area and could present an important component for regional operational oceanography when combined with numerical modelling. Some unresolved issues associated with the complex circulation and large instability of circulation in front of the Elbe River Estuary justify further considerations of this area using dedicated surveys and modelling efforts.  相似文献   

Plio-Pleistocene basanitic and melilititic series of the Bohemian Massif: K-Ar ages,major/trace element and Sr–Nd isotopic data     

Jaromír Ulrych  Lukáš Ackerman  Kadosa Balogh  Ernst Hegner  Emil Jelínek  Zoltan Pécskay  Antonín Přichystal  Brian G.J. Upton  Jiří Zimák  Radana Foltýnová 《Chemie der Erde / Geochemistry》2013

The Plio-Pleistocene volcanic rocks of the Bohemian Massif comprise a compositional spectrum involving two series: an older basanitic series (6.0–0.8 Ma) and a younger, melilititic series (1.0–0.26 Ma). The former consists of relatively undifferentiated basaltic rocks, slightly silica-undersaturated, with Mg# ranging from 62 to almost primitive mantle-type values of 74. The major and trace element characteristics correspond to those of primitive intra-plate alkaline volcanic rocks from a common sub-lithospheric mantle source (European Asthenospheric Reservoir – EAR) including positive Nb, and negative K and Pb anomalies. ⁸⁷Sr/⁸⁶Sr ratios of 0.7032–0.7034 and ¹⁴³Nd/¹⁴⁴Nd of 0.51285–0.51288 indicate a moderately depleted mantle source as for other mafic rocks of the central European volcanic province with signs of HIMU-like characteristics commonly attributed to recycling of subducted oceanic crust in the upper mantle during the Variscan orogeny. The melilititic series is characterized by higher degrees of silica-undersaturation, and high Mg# of 68–72 values, compatible with primitive-mantle-derived compositions. The high OIB-like Ce/Pb (19–47) and Nb/U (32–53) ratios indicate that assimilation of crustal material was negligible. In both series, concentrations of incompatible elements are mildly elevated and ⁸⁷Sr/⁸⁶Sr ratios (0.7034–0.7036) and ¹⁴³Nd/¹⁴⁴Nd ratios (0.51285–0.51288) overlap. Variations in incompatible element concentrations and isotopic compositions in the basanitic series and melilititic series can be explained by a lower degree of mantle melting for the latter with preferential melting of enriched mantle domains. The Sr and Nd isotopic compositions of both rock series are similar to those of the EAR. Minor differences in geochemical characteristics between the two series may be attributed to: (i) to different settings with respect to crust and lithospheric mantle conditions in (a) Western Bohemia (WB) and (b) Northeastern Bohemia (NEB) and the Northern Moravia and Silesia (NMS) areas, (ii) a modally metasomatized mantle lithosphere in WB in contrast to cryptically metasomatized domains in the NEB and NMS, (iii) different degrees of partial melting with very low degrees in WB but higher degrees in NEB and NMS. The geochemical and isotopic similarity between the Plio-Pleistocene volcanic rocks and those of the late Cretaceous and Cenozoic (79–6 Ma) suggests that their magmas came from compositionally similar mantle sources, that underwent low degrees of melting over an interval of ∼80 Ma. The Oligocene to Miocene basanitic series that accompanied the Plio-Pleistoicene basanitic series in the NMS region indicate that they shared a common mantle source. There is no geochemical evidence for thermal erosion of the lithospheric mantle or significant changes in mantle compositions within the time of a weak thermal perturbation in the asthenospheric mantle. These perturbations were caused by a dispersed mantle plume or passively upwelling asthenosphere in zones of lithospheric thinning.  相似文献   

Discrete fracture model for simulating waterflooding processes under fracturing conditions          下载免费PDF全文

Emil Gallyamov  Timur Garipov  Denis Voskov  Paul Van den Hoek 《国际地质力学数值与分析法杂志》2018,42(13):1445-1470

In our study, we develop a model for simulating fracturing processes in a poroelastic medium. The proposed approach combines the discrete fracture model enriched with contact plane mechanics. The model captures mechanical interactions of fractures and a deformable medium, fluid, and heat transfer in fractures and in a porous medium. Both effects of poroelasticity and thermoelasticity are accounted in our model. Mass and heat conservation equations are approximated by the finite volume method, and mechanical equilibrium equations are discretized by means of the Galerkin finite element approach. Two‐dimensional grid facets between 3‐dimensional finite elements are considered as possible fracture surfaces. Most of these facets are inactive from the beginning and are activated throughout the simulation. A fracture propagation criterion, based on Irwin's approach, is verified on each nonlinear iteration. When the criterion is satisfied, additional contact elements are added into finite element and discrete fracture model formulations respectively. The proposed approach allows modeling of existing natural and artificially created fractures within one framework. The model is tested on single‐ and multiple‐phase fluid flow examples for both isothermal and thermal conditions and verified against existing semianalytical solutions. The applicability of the approach is demonstrated on an example of practical interests where a sector model of an oil reservoir is simulated with different injection and production regimes.  相似文献   

Tidal wave transformations in the German Bight   总被引：2，自引：2，他引：0  

Emil V. Stanev  Rahma Al-Nadhairi  Joanna Staneva  Johannes Schulz-Stellenfleth  Arnoldo Valle-Levinson 《Ocean Dynamics》2014,64(7):951-968

Mesoscale and submesoscale dynamics associated with tidal wave transformations were addressed in the German Bight using numerical simulations. Tidal gauge and velocity observations in several locations were used to validate the numerical model. A downscaling approach included analysis of simulations with horizontal resolutions of 1, 0.4, and 0.2 km. It was shown that the modified tidal wave lost most of its energy after reflection or refraction over the eastern part of the German Bight. Energy loss resulted in a pronounced change of the wave’s spectral composition and generation of overtides. Tidal oscillations were modified by mesoscale processes associated with bathymetric channels. Semidiurnal and quarterdiurnal tides revealed very different spatial patterns. The former were aligned with the bathymetric channels, while the latter were rather “patchy” and had about half the spatial scales. In numerous areas around the bathymetric channels, the major axis of the M₄ ellipses was normal or at some angle with the major axis of the M₂ ellipses. Thus, higher harmonics developed “orthogonal” patterns that drove secondary circulations. Moreover, the ratio between spring and neap tidal amplitudes was relatively low in the Wadden Sea, showing reduced sensitivity of this very shallow area to fortnightly tidal variations. It was demonstrated that simulated hydrodynamics patterns help explain the physical mechanism shaping the median grain size distribution in the German Bight.  相似文献