A hybrid discrete–finite element model for numerical simulation of geomaterials

A hybrid discrete–finite element model is introduced for simulation of mechanical behavior of geomaterials. The soil or rock is modeled as a system of discrete balls that interact through normal and shear springs. The balls can be bonded at the contact points to withstand the applied deviatoric stresses. The important feature of this model is that the confining walls that can be imagined for example as the surrounding membrane or the mold in a physical test are modeled by deformable finite elements. This allows simulation of laboratory test features more realistically compared to the situations where the surrounding walls are rigid. The relationships between micro- and macro-properties are investigated in this paper as well. These relationships and the corresponding curves are helpful tools in calibration of the numerical model for the macroscopic elastic properties.     

Synsedimentary deformations in member 2 of the Mila Formation in the Central Alborz Mountains,Northern Iran

In Alborz Mountains, the thickness of sediments in Member 2 of the Mila Formation (Middle Cambrian) underwent dramatic and abrupt deformations, which were accompanied by normal and reverse faults and asymmetrical folds. These deformed sediments are covered by parallel beds both in the upper and lower sides. The existence of such extension and compression structures adjacent to each other indicates the influence of non-tectonic factors in their evolution. These deformations induced by downslope gliding of sediment packages are fully compatible with the dislocation model of Farrell (J Struct Geol 6:727–736, 1984) and Farrell and Eaton (1987). Moreover, the occurrence of synsedimentary deformations in a vast area in Alborz Mountains reinforces the probability of the impact of seismic shocks in their formation. In fact, incessant seismic events in the Middle Cambrian led to numerous submarine slumpings in the sediments of Member 2 of the Mila Formation.     

Development of a Darcy-Brinkman model to simulate water flow and tracer transport in a heterogeneous karstic aquifer (Val d’Orléans, France)

Darcy’s law is the equation of reference widely used to model aquifer flows. However, its use to model karstic aquifers functioning with large pores is problematic. The physics occurring within the karstic conduits requires the use of a more representative macroscopic equation. A hydrodynamic model is presented which is adapted to the karstic aquifer of the Val d’Orléans (France) using two flow equations: (1) Darcy’s law, used to describe water flow within the massive limestone, and (2) the Brinkman equation, used to model water flow within the conduits. The flow equations coupled with the transport equation allow the prediction of the karst transfer properties. The model was tested by using six dye tracer tests and compared to a model that uses Darcy’s law to describe the flow in karstic conduits. The simulations show that the conduit permeability ranges from 5?×?10^?6 to 5.5?×?10^?5?m² and the limestone permeability ranges from 8?×?10^?11 to 6?×?10^?10?m². The dispersivity coefficient ranges from 23 to 53 m in the conduits and from 1 to 5 m in the limestone. The results of the simulations carried out using Darcy’s law in the conduits show that the dispersion towards the fractures is underestimated.     

The Momentum 4-Vector Imparted by Gravitational Waves in Bianchi-Type Metrics

Considering the Møller, Weinberg and Qadir-Sharif's definitions in general relativity, we find the momentum 4-vector of the closed universe based on the Bianchi-type metrics. The momentum 4-vector (due to matter plus fields) is found to be zero. This result supports the viewpoints of Albrow and Tryon and extends the previous works by Cooperstock–Israelit, Rosen, Johri et al., Banerjee–Sen and Vargas who investigated the problem of the energy in Friedmann–Robertson–Walker universe and Salt?-Havare who studied the problem of the energy-momentum of the viscous Kasner-type space-times.     

Performance assessment of a concept propulsor: the thrust-balanced propeller

This paper presents the results of a numerical performance analysis to demonstrate the worthiness of a recently patented new concept propulsor, the so-called “thrust-balanced propeller (TBP)”. The main advantage of this unconventional propulsor is its inherent ability to reduce the unsteady effect of blade forces and moments when it is operating in a non-uniform wake flow. The propulsor comprises a pair of diametrically opposed blades that are connected to one another and mounted so as to be rotatable together through a limited angle about their spindle axis. A quasi-hydrodynamic approach is described and applied to perform the numerical analysis using a state-of-the-art lifting surface procedure for conventional propellers. Performance comparisons with a conventional fixed-pitch propeller are made for the blade forces and moments, efficiency, cavitation extents and fluctuating hull pressures. Bearing in mind the quasi-static nature of the analyses, the results present favourable performance characteristics for the thrust-balanced propeller and support the worthiness of the concept. However, the concept needs to be proved through physical model tests, which are planned to take in a cavitation tunnel.     

Implementation of baroclinic terms in a three-dimensional hydrodynamic model and its application to Anzali Port, Iran

Baroclinic terms have been implemented in a three-dimensional fully hydrodynamic model developed by Badiei et al. [2008. A three-dimensional non-hydrostatic boundary fitted model for free surface flows. International Journal for Numerical Methods in Fluids, 56(6), 607-627] modifying its momentum equations to account for density gradients and utilizing the scalar (salinity, temperature, etc.) conservation equation (SCE) and a state equation for the calculation of density. In the solution of advection-diffusion terms of the governing Navier-Stokes equations (NSE) and SCE, a symmetric splitting method was applied to ensure the long-term stability of simulations. Correction terms proposed by Ruddic et al. (1995) were applied to SCE to ensure the conservation of the scalar quantity. In the presence of baroclinic terms, the zero gradient pressure in the vertical direction in the vicinity of surface and bottom boundaries assumed by Badiei et al. [2008. A three-dimensional non-hydrostatic boundary fitted model for free surface flows. International Journal for Numerical Methods in Fluids, 56(6), 607-627] created spurious currents. This problem was solved by assuming a hydrostatic pressure variation at those boundaries. The ability of extended model was validated by comparing its results with an experimental test case. The simulation of hydrodynamic and salt intrusion at Anzali Port located at the southern coasts of Caspian Sea in Iran was carried out by the model with both barotropic and baroclinic modes. The simulated results with baroclinic mode show a better agreement with measured data as compared to the results of barotropic mode that clearly demonstrate the significance of baroclinic terms in the simulation of cyclic intrusion of salt wedge into the Port Basin.     

Issues in Eulerian–Lagrangian modeling of sediment transport under saltation regime

The saltation regime is very important for understanding the sediment transport mechanism. However,there is no consensus on a model for the saltation regime. This study answers several questions raised with respect to the Eulerian-Lagrangian modeling of sediment transport. The first question is why the previous saltation models that use different combinations of hydrodynamic forces yielded acceptable results? The second question is which shear lift model(i.e. a shear lift expression and its coefficient) is more appropriate? Another important question is which hydrodynamic forces have greater contributions to the saltation characteristics of a sediment particle? The last question is what are the contributions of the turbulence fluctuations as well as effects of using two-and three-dimensional(2 D and 3 D) models on the simulation results? In order to fairly answer these questions, a systematic study was done by considering different scenarios. The current study is the first attempt to clearly discuss these issues. A comprehensive 3 D saltation model for non-cohesive sediment was developed that includes all the hydrodynamic forces acting on the particle. The random nature of sediment transport was included using turbulent flow and bed-particle collision models. The eddy interaction model was applied to generate a3 D turbulent flow field. Bed-particle collisions were considered using the concept of a contact zone and a corresponding contact point. The validation of the model was done using the available experimental data for a wide range of sediment size(0.03 to 4.8 cm). For the first question, the results indicated that some of the hydrodynamic effects show opposing trends and some have negligible effects. With these opposing effects it is possible to adjust the coefficients of different models to achieve acceptable agreement with the same experimental data while omitting some aspects of the physics of the process. A suitable model for the shear lift force was developed by linking the lift coefficient to the drag coefficient and the contributions of the hydrodynamic forces and turbulence fluctuations as well as the consequences of using of 2 D and 3 D models were studied. The results indicate that the shear lift force and turbulent flow fluctuations are important factors for the saltation of both sand and gravel, and they cannot be ignored.     

Seismic fragility functions for code compliant and non-compliant RC SMRF structures in Pakistan

Fragility functions are derived for low-rise code compliant & non-compliant special moment resisting frames (SMRFs). Non-compliant SMRFs those built in low strength concrete and lacking confining ties in joint panel zones, commonly found in developing countries. Shake table tests were performed on single-storey and two-storey 1:3 reduced scale representative frames to understand the damage mechanism and develop deformation-based damage scale. The non-compliant SMRF experienced column flexure cracking, longitudinal bar-slip in beam and observed with cover concrete spalling from the joint panels. The code compliant SMRF experienced flexure cracks in beam/column, and experienced joint cracking under extreme shaking. Numerical modeling technique is developed for inelastic modeling of reinforced concrete frame with beam bar-slip and joint damageability using SeismoStruct. Natural accelerograms were used to analyze the considered frames through incremental dynamic analyses in SeismoStruct. A probabilistic based approach was used to derive fragility functions for the considered frames. An example case study is presented for damageability evaluation of structures for earthquakes of various return periods (43, 72, 475, 2475 years).     

Speciation and Determination of Inorganic Mercury and Methylmercury by Headspace Single Drop Microextraction and Electrothermal Atomic Absorption Spectrometry in Water and Fish

In this study, headspace single drop microextraction (HS‐SDME) method in combination with electrothermal atomic absorption spectrometry (ETAAS) method was developed and validated for the speciation and determination of inorganic mercury (iHg) and methylmercury (MeHg). MeHg and iHg species were reduced to volatile methylmercury hydride (CH₃HgH) and elemental mercury, respectively, in the presence of NaBH₄ and trapped onto a drop of acceptor phase in the tip of a microsyringe. Thiourea and ammonium pyrrolydinedithiocarbamate (APDC) were tested as the acceptor phase. The experimental parameters of the method such as microextraction time, temperature, NaBH₄ concentration, acceptor phase concentration, and pH of the medium were investigated to obtain distinctive conditions for mercury species. Possible interference effects have also been investigated. In order to validation of the method, analytical figures of merits such as accuracy, precision, limit of detection (LOD), limit of quantitation (LOQ), and linear working range have been evaluated. Accuracy of the method has been verified by analyzing certified reference materials (BCR 453 Tuna fish) and spiked samples. The proposed method was applied for the speciation and determination of mercury species in water and fish samples. Mercury species (MeHg and iHg) have been determined in the real samples with a relative error less than 10%.