Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.     

Developing constitutive models from EPR‐based self‐learning finite element analysis

A constitutive model that captures the material behavior under a wide range of loading conditions is essential for simulating complex boundary value problems. In recent years, some attempts have been made to develop constitutive models for finite element analysis using self‐learning simulation (SelfSim). Self‐learning simulation is an inverse analysis technique that extracts material behavior from some boundary measurements (eg, load and displacement). In the heart of the self‐learning framework is a neural network which is used to train and develop a constitutive model that represents the material behavior. It is generally known that neural networks suffer from a number of drawbacks. This paper utilizes evolutionary polynomial regression (EPR) in the framework of SelfSim within an automation process which is coded in Matlab environment. EPR is a hybrid data mining technique that uses a combination of a genetic algorithm and the least square method to search for mathematical equations to represent the behavior of a system. Two strategies of material modeling have been considered in the SelfSim‐based finite element analysis. These include a total stress‐strain strategy applied to analysis of a truss structure using synthetic measurement data and an incremental stress‐strain strategy applied to simulation of triaxial tests using experimental data. The results show that effective and accurate constitutive models can be developed from the proposed EPR‐based self‐learning finite element method. The EPR‐based self‐learning FEM can provide accurate predictions to engineering problems. The main advantages of using EPR over neural network are highlighted.     

Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Is the age problem resolved?

The cosmological, astrophysical, and nucleocosmochronological methods for estimating the age of the universe and the corresponding uncertainties are comparatively studied in the present paper. We are led to the conclusion that the new measurements of cosmological parameters, and the recent estimates of the age of globular clusters have led to the gradual disappearance of the age problem from the arena of modern cosmology.     

SHRIMP zircon dating and Sm/Nd isotopic investigations of Neoproterozoic granitoids, Eastern Desert, Egypt

There is an increasing evidence for the involvement of pre-Neoproterozoic zircons in the Arabian–Nubian Shield, a Neoproterozoic crustal tract that is generally regarded to be juvenile. The source and significance of these xenocrystic zircons are not clear. In an effort to better understand this problem, older and younger granitoids from the Egyptian basement complex were analyzed for chemical composition, SHRIMP U–Pb zircon ages, and Sm–Nd isotopic compositions. Geochemically, the older granitoids are metaluminous and exhibit characteristics of I-type granites and most likely formed in a convergent margin (arc) tectonic environment. On the other hand, the younger granites are peraluminous and exhibit the characteristics of A-type granites; these are post-collisional granites. The U–Pb SHRIMP dating of zircons revealed the ages of magmatic crystallization as well as the presence of slightly older, presumably inherited zircon grains. The age determined for the older granodiorite is 652.5 ± 2.6 Ma, whereas the younger granitoids are 595–605 Ma. Xenocrystic zircons are found in most of the younger granitoid samples; the xenocrystic grains are all Neoproterozoic, but fall into three age ranges that correspond to the ages of other Eastern Desert igneous rocks, viz. 710–690, 675–650 and 635–610 Ma. The analyzed granitoids have (+3.8 to +6.5) and crystallization ages, which confirm previous indications that the Arabian–Nubian Shield is juvenile Neoproterozoic crust. These results nevertheless indicate that older Neoproterozoic crust contributed to the formation of especially the younger granite magmas.     

Travel time based approach for the assessment of vulnerability of karst groundwater: the Transit Time Method

Karst aquifers represent important water resources in many parts of the world. Unfortunately, karst aquifers are characterised by high contamination risks. This paper presents a travel time based method for the estimation of karst groundwater vulnerability. It considers (1) physics-based lateral flow within the uppermost weathered zone (epikarst) in a limestone-dominated region and (2) high velocities of vertical infiltration at discrete infiltration points (e.g. sinkholes) or lines (e.g. dry valleys, faults). Consequently, the Transit Time Method honours the actual flow path within the unsaturated zone of a karst aquifer system. A test site in Northern Jordan was chosen for the demonstration of the assessment technique, i.e. the catchment area of the Qunayyah Spring north of the capital Amman. The results demonstrate that zones of highest vulnerability lie within valleys and nearby main fault zones. It also reveals that regions, categorised as protected areas by other methods due to thick unsaturated zones, contribute to a major degree to the total risk.     

Geomorphic changes in Ras Al-Subiyah area,Kuwait

The Ras Al-Subiyah area is considered one of the most promising areas in Kuwait for future development. This development will include a new town called Subiyah and its associated infrastructure. This area is also being considered as the location for connection between Boubyan Island, which is now undergoing major development and the Kuwait mainland. The present study investigates the geomorphology of the Ras Al-Sabiyah area in the northern sector of Kuwait. The study area is generally flat, and it is located west of the Jal Az-Zor escarpment. It is bordered on the east by the Khor Al-Sabiyah tidal channel and on the south by Kuwait Bay. The area receives sediments from several sources; currently the most important are aeolian sediments and the deposition of mud delivered through the Khor Al-Sabiyah from the Iraqi marshes. The study area has been subjected to severe environmental changes due to the Gulf wars and the drainage of Iraqi marshes and the associated artificial changes in fluvial system. Twenty-two surface sediments were collected from the Ras Al-Subiyah area. Samples were collected to include the main geomorphologic characteristic features of the study area. Field observations and remote sensing images from 1990 and 2001 were used to produce an updated geomorphologic map for the Ras Al-Subiyah and a map showing geomorphic changes between 1990 and 2001. Grain size of the surface sediment ranges from gravel to medium sand. In general, grain size statistical analysis indicates that most of the areas are composed of two or more classes of sands transported and deposited from different sources including aeolian, sabkhas, river and the bays. The variability in the grain size statistical parameters may be attributed to the complexity of surface morphology as well as the diversity in the type of depositional environment in the Ras Al-Subiyah area. The total area subjected to change during the 12-year period (1990–2001) is about 32 km² as calculated using GIS techniques.