A distributed parallel multiple-relaxation-time lattice Boltzmann method on general-purpose graphics processing units for the rapid and scalable computation of absolute permeability from high-resolution 3D micro-CT images

Digital rock physics (DRP) is a rapidly evolving technology targeting fast turnaround times for repeatable core analysis and multi-physics simulation of rock properties. We develop and validate a rapid and scalable distributed-parallel single-phase pore-scale flow simulator for permeability estimation on real 3D pore-scale micro-CT images using a novel variant of the lattice Boltzmann method (LBM). The LBM code implementation is designed to take maximum advantage of distributed computing on multiple general-purpose graphics processing units (GPGPUs). We describe and extensively test the distributed parallel implementation of an innovative LBM algorithm for simulating flow in pore-scale media based on the multiple-relaxation-time (MRT) model that utilizes a precise treatment of body force. While the individual components of the resulting simulator can be separately found in various references, our novel contributions are (1) the integration of all of the mathematical and high-performance computing components together with a highly optimized code implementation and (2) the delivery of quantitative results with the simulator in terms of robustness, accuracy, and computational efficiency for a variety of flow geometries including various types of real rock images. We report on extensive validations of the simulator in terms of accuracy and provide near-ideal distributed parallel scalability results on large pore-scale image volumes that were largely computationally inaccessible prior to our implementation. We validate the accuracy of the MRT-LBM simulator on model geometries with analytical solutions. Permeability estimation results are then provided on large 3D binary microstructures including a sphere pack and rocks from various sandstone and carbonate formations. We quantify the scalability behavior of the distributed parallel implementation of MRT-LBM as a function of model type/size and the number of utilized GPGPUs for a panoply of permeability estimation problems.     

The First Hop: Use of Beaufort Sea Deltas by Hatch-Year Semipalmated Sandpipers

River deltas along Alaska’s Beaufort Sea coast are used by hatch-year semipalmated sandpipers (Calidris pusilla) after leaving their terrestrial natal sites, but the drivers of their use of these stopover sites on the first “hop” of fall migration are unknown. We quantified sandpiper temporal distribution and abundance as related to food resources at three river deltas during the beginning of their fall migration (post-breeding period) to compare the habitat quality among these deltas. We conducted population counts, sampled invertebrates, and captured birds to collect blood samples from individuals for triglyceride and stable isotope analyses to determine fattening rates and diet. Patterns of sandpiper and invertebrate abundance were complex and varied among deltas and within seasons. River deltas were used by sandpipers from late July to late August, and peak sandpiper counts ranged from 1000 to 4000 individuals, of which 98% were hatch-year semipalmated sandpipers. Isotopic signatures from blood plasma samples indicated that birds switched from a diet of upland tundra to delta invertebrate taxa as the migration season progressed, suggesting a dependence on delta invertebrates. Despite differences in diet among deltas, we found no differences in fattening rates of juvenile sandpipers as indicated by triglyceride levels. The number of sandpipers was positively associated with abundance of Amphipoda and Oligochaeta at the Jago and Okpilak-Hulahula deltas; an isotopic mixing model indicated that sandpipers consumed Amphipoda and Oligochaeta at Jago, mostly Chironomidae at Okpilak-Hulahula and Spionidae at Canning. Regardless of the difference in sandpiper diets at the Beaufort Sea deltas, their similar fattening rates throughout the season indicate that all of these stopover sites provide a critical food resource for hatch-year sandpipers beginning their first migration.     

Great South Bay After Sandy: Changes in Circulation and Flushing due to New Inlet

The coastal ocean model FVCOM is applied to quantify the changes in circulation, flushing, and exposure time in Great South Bay, New York, after Superstorm Sandy breached the barrier island in 2012. Since then, the lagoon system is connected to the Atlantic via five instead of four inlets. The model simulations are run on two high-resolution unstructured grids, one for the pre-breach configuration, one including the new inlet, with tidal-only forcing, and summer and winter forcing conditions. Despite its small cross-sectional size, the breach has a relatively large net inflow that leads to a strengthening of the along-bay through-flow in Great South Bay (GSB); the tidally driven volume transport in central GSB quadrupled. The seasonal forcing scenarios show that the southwesterly sea breeze in summer slows down the tidally driven flow, while the forcing conditions in winter are highly variable, and the circulation is dependent on wind direction and offshore sea level. Changes in flushing and exposure time associated with the modified transport patterns are evaluated using a Eulerian passive tracer technique. Results show that the new inlet produced a significant decrease in flushing time (approximately 35% reduction under summer wind conditions and 20% reduction under winter wind conditions). Maps of exposure time reflect the local changes in circulation and flushing.     

Investigation of metal pollution in Moryayla (Erzurum) and surrounding stream sediments,Turkey

The enrichment factor (EF), geochemical index (Igeo), pollution factor (Cf), and pollutant load index (PLI) were calculated using average shale metal concentration values to determine metal pollution in Moryayla and surrounding stream sediments. It was determined that the average concentration values are above the world average shale values. According to the results obtained, the mean accumulation levels of heavy metals in the study area were Mn > Zn > Cu > Pb > As > Cr > Ni > Sb. Considerable degree of enrichment was observed because the elements As, Cu, and Pb are at EF > 5 at locations 47, 53, 70, and 71. The values obtained from the Igeo, Cf, and PLI are a wide range of pollutants (from unpolluted to moderately polluted to extremely polluted to) for the various elements studied.     

Water stress index and its implication for agricultural land-use policy in Thailand

Land-use change and expansion for agriculture significantly affect freshwater resource availability which in turn results in different levels of water scarcity and competition among users in different areas and time. The study uses the water stress index as an indicator for determining the potential impact of water use considering water deprivation potential. Considering the case of Thailand, the temporal aspects of water withdrawal and availability have been analyzed and characterized as the monthly water stress index for the 25 watersheds in the country. Wide variations of monthly water stress index are obtained; extreme water stress is observed in the areas of the Chao Phraya and Tha Chin watersheds (the central region) during the dry season because of the large areas cultivating second rice (rice planted in dry season). The target of the Thai government on agricultural land-use change and zoning toward the conversion of about 0.37 M·ha of the upland paddy fields with low productivity to sugarcane was evaluated showing serious implications on the monthly water stress index of the relevant watersheds and the water scarcity footprint potentials of rice and sugarcane production. The results reveal that proper policy can help reduce the amount of water requirement for agriculture in June, July, August and September by about 60–220 M·m³, which in turn results in the decrease in monthly water stress index values. Nevertheless, appropriate measures of water resource management for agriculture still need to be designed to avoid water competition as well as protect the ecosystem.     

Determining ammonia nitrogen decay rate of Malaysian river water in a laboratory flume

There is lack of information regarding ammonia nitrogen (AN), (i.e. NH₃-N) decay rate of river water in tropical regions like Malaysia. AN decay rate is a very important kinetic parameter to estimate NH₃-N, nitrate nitrogen (NO₃-N) and dissolved oxygen concentrations of river water by using computer models. This study presents determination of ammonia nitrogen decay rate of river water in the tropical environment of Malaysia. A laboratory flume was used to conduct twelve experiments. The flume was used to represent the turbulent condition of a typical river. Ammonia nitrogen decay rate for the tropical environment of Malaysia was observed to be between 0.194 and 0.554 per day. Median value of AN decay rate was 0.26 per day, which is slightly lower than the global median value of 0.295 per day. To check the accuracy of flume experiments, the AN decay rate of Pusu River obtained from the flume experiment was used to calibrate and validate ammonia nitrogen concentration of the river by using water quality analysis and simulation program (WASP). Very good calibration and validation results were achieved, which substantiated the accuracy of the flume experiments.     

Effect of confining stress and loading frequency on dynamic behavior of municipal solid waste in Kahrizak landfill

In order to analyze the stability of a landfill site, it should analyze some properties of waste introducing the main structural elements. Up to now, it has not been done much research on the properties of municipal solid waste. In addition, due to the differences in waste properties from one country to another and even from one to another landfill site, it is impossible to generalize the results. These conditions caused local research on the evaluation of static and dynamic parameters of municipal solid waste to be done. Because Iran is a seismic country, the short-term behavior of waste controls landfill sites stability during seismic loads; so it is necessary to know the dynamic behavior of these materials. In this paper, about 18 cyclic tests were performed, and in addition to determining the dynamic parameters of municipal solid waste of Tehran Kahrizak Disposal Site using cyclic triaxial test, the effect of confining stress and loading frequency on dynamic properties of these materials was evaluated. The results have shown that with increasing the confining stress and loading frequency, shear modulus and damping ratio was increased and decreased, respectively, and it is related to the composition of the municipal solid waste materials.     

Isolation of <Emphasis Type="Italic">Bacillus cereus</Emphasis> from botanical soil and subsequent biodegradation of waste engine oil

Waste engine oil causes a vital environmental pollution when it spill during change and transportation and products of waste engine oil causes lethal effects to the living systems. Thus, abiotic and biotic approaches are being extensively used for removal of waste engine oil pollution. Therefore in present study, waste engine oil degradation was accomplished by a new bacterial culture, isolated from the soil by an enrichment technique. Morphological, biochemical and gene sequence analysis revealed that isolate was Bacillus cereus. Subsequently, biodegradation potential of B. cereus for waste engine oil was studied. Experimental variables, such as pH, substrate concentration, inoculum size, temperature and time on the biodegradation, were checked in mineral salt medium. The biodegradation efficiency of B. cereus was determined by gravimetry, UV–visible spectrophotometry and gas chromatography. In addition, waste engine oil was also characterized by GC–MS and FTIR for its major constituents, which showed total 38 components in waste engine oil, including hopanes, benzopyrene, long-chain aliphatic hydrocarbons, dibenzothiophenes, biphenyl and their derivatives. Results of successive biodegradation indicated that B. cereus was capable to degrade 1% of waste engine oil with 98.6% degradation potential at pH 7 within 20 days. Hence, B. cereus presents an innovative tool for removing the engine oil from the contaminated area.     

Hydrogeochemical Features of Thermal Waters of South Trungbo (Central Vietnam)

The results of studying the features of the hydrogeological structure and chemical and isotope composition of thermal waters from the central part of Vietnam that are characterized by intense manifestations of intrusive magmatism are presented. It is established that low–and high–thermal waters with temperature varying within 30–85°C are developed in the area under study. The value of total mineralization of the hydrotherms ranges from 0.05 to 10.05 g/dm³. It is assumed that the circulation of thermal waters that are different in temperature and chemical composition occurs at two levels. The regular change of the hydrotherm composition in the direction from mineralized chloride sodium, including with increased Ca content, to fresh sodium bicarbonate is revealed. The ratio of δ¹⁸O–δ²H isotopes indicates that the water component is based on meteoric water. In the coastal areas, there is an isotope shift towards the ocean waters, which is also confirmed by the hydrogeochemical data. The key factors for forming the chemical composition of the thermal waters in South Trungbo are their genetic type, the interaction processes in the “water–rock–gas–organic substance” system, and their equilibrium–nonequilibrium state.     

An orthotropic interface damage model for simulating drying processes in soils

The study of drying process in soils has received an increased attention in the last few years. This is very complex phenomenon that generally leads to the formation and propagation of desiccation cracks in the soil mass. In recent engineering applications, high aspect ratio elements have proved to be well suited to tackle this type of problem using finite elements. However, the modeling of interfaces between materials with orthotropic properties that generally exist in this type of problem using standard (isotropic) constitutive model is very complex and challenging in terms of the mesh generation, leading to very fine meshes that are intensive CPU demanding. A novel orthotropic interface mechanical model based on damage mechanics and capable of dealing with interfaces between materials in which the strength depends on the direction of analysis is proposed in this paper. The complete mathematical formulation is presented together with the algorithm suggested for its numerical implementation. Some simple yet challenging synthetic benchmarks are analyzed to explore the model capabilities. Laboratory tests using different textures at the contact surface between materials were conducted to evaluate the strengths of the interface in different directions. These experiments were then used to validate the proposed model. Finally, the approach is applied to simulate an actual desiccation test involving an orthotropic contact surface. In all the application cases the performance of the model was very satisfactory.