Internal inlet for wave generation and absorption treatment

A new method of implementing, in two-dimensional (2-D) Navier–Stokes equations, a numerical internal wave generation in the finite volume formulation is developed. To our knowledge, the originality of this model is on the specification of an internal inlet velocity defined as a source line for the generation of linear and non-linear waves. The use of a single cell to represent the source line and its transformation to an internal boundary condition proved to be an interesting alternative to the common procedure of adding a mass source term to the continuity equation within a multi-cell rectangular region. Given the reduction of the source domain to a one-dimensional region, this simple new type of source introduced less perturbation than the 2-D source type. This model was successfully implemented in the PHOENICS code (Parabolic Hyperbolic Or Elliptic Numerical Integration Code Series). In addition, the volume of fluid (VOF) fraction was used to describe the free surface displacements. A friction force term was added to the momentum transport equation in the vertical direction, in order to enhance wave damping, within relatively limited number of cells representing the sponge layers at the open boundaries. For monochromatic wave, propagating on constant water depth, numerical and analytical results showed good agreements for free surface profiles and vertical distribution of velocity components. For solitary wave simulation, the wave shape and velocity were preserved; while, small discrepancy in the tailing edge of the free surface profiles was observed. The suitability of this new numerical wave generation model for a two source lines extension was investigated and proven to be innovative. The comparisons between numerical, analytical and experimental results showed that the height of the merging waves was correctly reproduced and that the reflected waves do not interact with the source lines.     

Assessment of reference evapotranspiration across an arid urban environment having poor data monitoring system

Estimation of reference evapotranspiration (ET₀) in urban areas is challenging but essential in arid urban climates. To evaluate ET₀ in an urban environment and non-urban areas, air temperature and relative humidity were measured at five different sites across the arid city of Isfahan, Iran, over 4 years. Wind speed and sunshine hours were obtained from an urban surrounding weather station over the same period and used to estimate ET₀. Calculated ET₀ was compared with satellite-based ET₀ retrieved from the MOD16A2 PET product. Although MODIS PET was strongly correlated with the Valiantzas equation, it overestimated ET₀ and showed average accuracy (r = 0.93–0.94, RMSE = 1.18–1.28 mm/day, MBE = 0.73–0.84 mm/day). The highest ET₀ differences between an urban green space and a non-urban area were 1.1 and 0.87 mm/day, which were estimated by ground measurements and MODIS PET, respectively. The sensitivity of ET₀ to wind speed and sunshine hours indicated a significant effect on cumulative ET₀ at urban sites compared to the non-urban site, which has a considerable impact on the amount of irrigation required in those areas. Although MODIS PET requires improvement to accurately reflect field level microclimate conditions affecting ET₀, it is beneficial to hydrological applications and water resource managers especially in areas where data is limited. In addition, our results indicated that using limited data methods or meteorological data from regional weather stations, leads to incorrect estimation of ET₀ in urban areas. Therefore, decision-makers and urban planners should consider the importance of precisely estimating ET₀ to optimize management of urban green space irrigation, especially in arid and semi-arid climates such as the city of Isfahan.     

An Improved Workflow for Permeability Estimation from Image Logs with Uncertainty Quantification

Natural Resources Research - Production forecasts are extremely sensitive to the spatial distribution of certain reservoir parameters, especially permeability, which controls fluid flow and...     

Wavelet and Gaussian Approaches for Estimation of Groundwater Variations Using GRACE Data

In this study, a scheme is presented to estimate groundwater storage variations in Iran. The variations are estimated using 11 years of Gravity Recovery and Climate Experiments (GRACE) observations from period of 2003 to April 2014 in combination with the outputs of Global Land Data Assimilation Systems (GLDAS) model including soil moisture, snow water equivalent, and total canopy water storage. To do so, the sums of GLDAS outputs are subtracted from terrestrial water storage variations determined by GRACE observations. Because of stripping errors in the GRACE data, two methodologies based on wavelet analysis and Gaussian filtering are applied to refine the GRACE data. It is shown that the wavelet approach could better localize the desired signal and increase the signal‐to‐noise ratio and thus results in more accurate estimation of groundwater storage variations. To validate the results of our procedure in estimation of ground water storage variations, they are compared with the measurements of pisometric wells data near the Urmia Lake which shows favorable agreements with our results.     

Evidence for an underground runoff and soil permeability at the Ouled Fayet (Algiers,Algeria) subsurface landfill pilot project from geophysical investigations

Results from geophysical investigations (electrical resistivity, electromagnetic mapping and seismic refraction) on an excavated cell of the Ouled Fayet (Algiers, Algeria) pilot landfill indicated the presence of an underground runoff and permeable soil underneath the cell. These results contradict those obtained by a feasibility study, based, however, only on the analysis of seventy-six 10-m drilling cores. The 1D boreholes information has been proven to be insufficient and to give biased results. The presence of water at depth is evidenced by lower resistivity, high conductivity anomalies and increase of P-wave velocity. Thus, to the contrary of what is claimed in the feasibility study, a threat of leachate pollution is real. This study shows that landfill construction studies cannot give trustful results without geophysical investigations. More specifically, in Algeria, it is imperative to elaborate a landfill construction code, which should include mandatory geophysical prospecting and deeper drilling cores.     

Dynamic earth pressure on rigid retaining walls induced by a neighboring machine foundation,by the meshless local Petrov-Galerkin method

Dynamic earth pressure induced by machine foundations on a neighboring retaining wall is analyzed with emphasis on factors which control the intensity and location of the design forces. The meshless local Petrov-Galerkin(MLPG) method is used to analyze the problem for a variety of retaining wall and machine foundation geometries. The soil medium is assumed to be homogeneous and visco-elastic. The machine foundation is idealized as a harmonic sinusoidal dynamic force often encountered in practice. A number of analyses have been made to reveal the effect of the loading frequency, the location and size of the foundation and the soil shear wave velocity on the distribution and magnitude of the dynamic earth pressure. Results indicate that there is a critical frequency and a critical location for which the passive pressure takes the maxima in the entire duration of the dynamic load.     

New structural implications for the central Sahara (Algeria), from the revisited Upper Carboniferous “Hassi Bachir” Formation: Paleomagnetic constraints

Paleomagnetic investigations of the folded Upper Namurian–Lower Moscovian “Hassi Bachir” Formation cropping out in the “Ahnet” basin (Central Sahara, Algeria) yield two magnetic components. A pre-folding primary magnetization (D = 136.1°, I = 22.0°, k = 217, α₉₅ = 2.6°) enables us to define a paleomagnetic pole (32.8°S, 55.7°E, K = 328 and A₉₅ = 2.0°) which better constrains a paleopole that was determined by Daly and Irving [Daly, L., Irving, E., 1983. Paléomagnétisme des roches carbonifères du Sahara central; analyse des aimantations juxtaposées; configurations de la Pangée. Ann. Geophys. 1, 207–216] for the same formation. A secondary component consists in a synfolding remagnetization and shows that post-Permian tectonics account for at least about half of the total folding in the studied area. This indicates that Mesozoic folding noted 150 km to the West in the Reggane basin [Smith, B., Derder, M.E.M., Henry, B., Bayou, B., Amenna, M., Djellit, H., Yelles, A.K., Garces, M., Beamud, E., Callot, J.P., Eschard, R., Chambers, A., Aifa, T., Ait Ouali, R., Gandriche, H., 2006. Relative importance of the Hercynian and post-Jurassic tectonic phases in the Saharan platform: a palaeomagnetic study of Jurassic sills in the Reggane basin (Algeria). Geophys. J. Int. 167, 380–396] is not local and affected at least the entire north-western part of the Hoggar area. This reconfirms that the folding of the Paleozoic cover in the Sahara platform should not be restricted to the Hercynian orogeny.     

Bearing capacity of non-associative coaxial granular materials by upper bound limit analysis and finite elements

In this paper, an upper bound estimate of the limit load on non-associative coaxial granular materials is presented. The kinematic approach of the upper bound limit analysis has been utilised. The failure mechanism is assumed to coincide with the direction of the shear bands at every point throughout the body. The shear band orientation in non-associative coaxial materials, i.e. those with the same major principal stress and major principal strain increment directions, can be found based on the angle of dilation and the major principal stress direction. Therefore, having known the stress field at limiting equilibrium, the orientation of the shear bands and hence, the failure mechanism can be obtained. In this study, the stress field is first determined by the method of stress characteristics. Then, the finite element interpolation technique is used to interpolate the stress field and to find the orientation of the shear bands at every point within the field. Once the failure mechanism and the stress state at every point along velocity discontinuities have been found, the upper bound limit analysis has been performed to estimate the limit load.