List of forthcoming papers

List of forthcoming papers     

The meteoritic origin of Tutankhamun's iron dagger blade

Scholars have long discussed the introduction and spread of iron metallurgy in different civilizations. The sporadic use of iron has been reported in the Eastern Mediterranean area from the late Neolithic period to the Bronze Age. Despite the rare existence of smelted iron, it is generally assumed that early iron objects were produced from meteoritic iron. Nevertheless, the methods of working the metal, its use, and diffusion are contentious issues compromised by lack of detailed analysis. Since its discovery in 1925, the meteoritic origin of the iron dagger blade from the sarcophagus of the ancient Egyptian King Tutankhamun (14th C. BCE) has been the subject of debate and previous analyses yielded controversial results. We show that the composition of the blade (Fe plus 10.8 wt% Ni and 0.58 wt% Co), accurately determined through portable x‐ray fluorescence spectrometry, strongly supports its meteoritic origin. In agreement with recent results of metallographic analysis of ancient iron artifacts from Gerzeh, our study confirms that ancient Egyptians attributed great value to meteoritic iron for the production of precious objects. Moreover, the high manufacturing quality of Tutankhamun's dagger blade, in comparison with other simple‐shaped meteoritic iron artifacts, suggests a significant mastery of ironworking in Tutankhamun's time.     

Interactive visualization of the largest radioastronomy cubes

3D visualization is an important data analysis and knowledge discovery tool, however, interactive visualization of large 3D astronomical datasets poses a challenge for many existing data visualization packages. We present a solution to interactively visualize larger-than-memory 3D astronomical data cubes by utilizing a heterogeneous cluster of CPUs and GPUs. The system partitions the data volume into smaller sub-volumes that are distributed over the rendering workstations. A GPU-based ray casting volume rendering is performed to generate images for each sub-volume, which are composited to generate the whole volume output, and returned to the user. Datasets including the HI Parkes All Sky Survey (HIPASS – 12 GB) southern sky and the Galactic All Sky Survey (GASS – 26 GB) data cubes were used to demonstrate our framework’s performance. The framework can render the GASS data cube with a maximum render time <0.3 s with 1024 × 1024 pixels output resolution using three rendering workstations and eight GPUs. Our framework will scale to visualize larger datasets, even of Terabyte order, if proper hardware infrastructure is available.     

Prediction of Land Subsidence Under Cyclic Pumping Based on Laboratory and Numerical Simulations

In this study subsidence due to groundwater withdrawal was investigated. Kerman Province in Iran is struggling with land subsidence problem due to extensive groundwater withdrawal mainly for farming. The rate and type of groundwater withdrawal has very important impact on settlement rate. In this research, effective parameters on land subsidence caused by groundwater withdrawal were determined by laboratory tests. Sampling had done up to depth of 300 m mainly with remolded specimens from Shams-abad, Nouq plain in Kerman province. Similar to the field preconsolidation pressure was applied on specimens in the laboratory. Rate of applied stress on prepared specimens was similar to effect of oscillation of groundwater level. In order to model the actual soil behavior in the laboratory, one-dimensional consolidation device (odometer) was adopted for testing. In these tests, the effect of loading caused by seasonal oscillation of groundwater table is considered by means of cyclic loading in the testing which has great effect on rate of settlements. The results of tests show that when the water table level periodically increases and decreases the amount of settlement decrease, comparing with the case when the groundwater table drop to a constant level. In order to predict the further effects of groundwater level oscillation and actual field condition on land subsidence, a finite element model based on Biots’ three-dimensional consolidation theory was developed. After calibration of finite element model with laboratory tests, this model was used for prediction the effect of groundwater level oscillation on actual field conditions.     

Geochemical Characteristics,Origin, and Evolution of Ore‐Forming Fluids of the Khut Copper Skarn Deposit,West of Yazd in Central Iran

The Khut copper skarn deposit is located at about 50 km northwest of Taft City in Yazd province in the middle part of the Urumieh‐Dokhtar magmatic arc. Intrusion of granitoid of Oligocene–Miocene age into carbonate rocks of the Triassic Nayband Formation led to the formation of marble and a calcic skarn. The marble contains high grade Cu mineralization that occurs mainly as open space filling and replacement. Cu‐rich sulfide samples from the mineralized marble are also anomalous in Au, Zn, and Pb. In contrast, the calcic skarn is only weakly anomalous in Cu and W. The calcic skarn is divided into garnet skarn and garnet–pyroxene skarn zones. Paragenetic relationships and microthermometric data from fluid inclusions in garnet and calcite indicate that the compositional evolution of skarn minerals occurred in three main stages as follows. (i) The early prograde stage, which is characterized by Mg‐rich hedenbergite (Hd_53.7Di_42.3–Hd_86.1Di_9.5) with Al‐bearing andradite (69.8–99.5 mol% andradite). The temperature in the early prograde skarn varies from 400 to 500°C at 500 bar. (ii) The late prograde stage is manifested by almost pure andradite (96.2–98.4 mol% andradite). Based on the fluid inclusion data from garnet, fluid temperature and salinity in this stage is estimated to vary from 267 to 361°C and from 10.1 to 21.1 wt% NaCl equivalent, respectively. Pyrrhotite precipitation started during this stage. (iii) The retrograde stage occurs in an exoskarn, which consists of an assemblage of ferro‐actinolite, quartz, calcite, epidote, chlorite, sphalerite, pyrite, and chalcopyrite that partially replaces earlier mineral assemblages under hydrostatic conditions during fracturing of the early skarn. Fluids in calcite yielded lower temperatures (T < 260°C) and fluid salinity declined to ～8 wt% NaCl equivalent. The last stage mineralization in the deposit is supergene weathering/alteration represented by the formation of iron hydroxide, Cu‐carbonate, clay minerals, and calcite. Sulfur isotope data of chalcopyrite (δ³⁴S of +1.4 to +5.2‰) show an igneous sulfur source. Mineralogy and mineral compositions of the prograde assemblage of the Khut skarn are consistent with deposition under intermediately oxidized and slightly lower fS₂ conditions at shallow crustal levels compared with those of other typical Fe‐bearing Cu–Au skarn systems.     

Mixed precision parallel preconditioner and linear solver for compositional reservoir simulation

Reservoir simulation role in value creation and strategic management decisions cannot be over emphasized. Simulation of complex challenging reservoirs with millions of grid blocks especially in compositional mode is very time-consuming even with fast modern computers. On the other hand, high price of cluster supercomputers prevents them for being commonly used for fast simulation of such reservoirs. In recent years, the development of many-core processors like cell processors, DSPs, and graphical processing units (GPUs) has provided a very cost-effective hardware platform for fast computational operations. However, programming for such processors is much more difficult than conventional CPUs, and new parallel algorithm design and special parallel implementation methods are needed. Using the computational power of CPUs, GPUs, and/or any other processing unit, Open Computing Language (OpenCL) provides a framework for programming for heterogeneous platforms. In this paper, OpenCL is used to employ the computational power of a GPU to build a preconditioner and solve the linear system arising from compositional formulation of multiphase flow in porous media. The proposed parallel preconditioner is proved to be quite effective, even in heterogeneous porous media. Using data-parallel modules on GPU, the preconditioner/solver runtime reduced at least 1 order of magnitude compared to their serial implementation on CPU.     

Removal of copper ions from aqueous solution using silica derived from rice straw: comparison with activated charcoal

The batch removal of copper(II) ions from aqueous solution under different experimental conditions using alkali-leached silica and activated charcoal was investigated in this study. The copper(II) uptake was dependent on varying time, pH, copper concentration and temperature. Copper sorption was found fast reaching equilibrium within 1 h with better performance for alkali-leached silica than charcoal. Copper sorption was low at low pH values and increased with rise in initial pH-value until 6.7. Sorption fits well the Langmuir and Freundlich equations with higher uptake by increasing temperature. According to Langmuir equation, the maximum uptake of Cu(II) ions by alkali-leached SiO₂ and charcoal was found to be 242.5 and 94.4 mmol/g at temperature 60 °C and pH 6. Thermodynamic studies confirm that the process was spontaneous and endothermic nature. Kinetic data for Cu(II) sorption was found to follow pseudo-second-order model.