Search

This paper describes the experiments carried out on a mediaeval masonry tower in the historic centre of Lucca and some finite element numerical simulations of the tower's experimental response. The Guinigi Tower, one of the most iconic monuments in Lucca, has been continuously monitored by high-sensitivity seismic stations that recorded the structure's response to the dynamic actions of the surrounding environment. The monitoring campaign results have been analysed to show the effectiveness of dynamic monitoring as a valuable source of information on the structural properties of the tower. The dynamic analyses of the tower and the surrounding palace subjected to some seismic events recorded during the experiments have highlighted the capabilities of experiment-based finite element modelling. The calibration of the finite element model and the numerical analysis have been carried out by resorting to procedures developed at ISTI-CNR and able to consider the nonlinear behaviour of masonry materials. [ABSTRACT FROM AUTHOR]

In a recent work, the authors demonstrated the merits of the diametrical core deformation technique (DCDT) for stress determination at an excavation face. An analytical model was developed and verified with laboratory tests. This paper aims to validate the DCDT in the field using a well-documented case study from a Canadian underground mine. Rock cores are extracted from the floor of a footwall drive 1000 m below surface for the purpose of estimating stress magnitude and orientation. The DCDT is applied to the extracted cores and the state of stress in the drift floor is estimated. On the other hand, a mine wide FLAC3D model was developed using the geomechanical data and in situ stress regime obtained from a previous study. The model is a pared-down version of a comprehensive FLAC3D mine wide model that simulates past mine-and-fill sequence and was calibrated with various burst criteria and recorded mining-induced seismic events. A comparison of FLAC3D and DCDT in terms of floor stress estimates and orientation shows good agreement. The DCDT could be used as a ground control tool to help detect stress-related ground control issues. Highlights: Extraction of rock cores from a drift floor 1000 m below surface for the purpose of stress estimation. Application of diametrical core deformation technique (DCDT) to estimate the biaxial state of stress—magnitude and orientation—in the floor. Development of a mine-wide 3D numerical model using geomechanical data and in situ stress regime. Validation of DCDT by comparison of stress orientation and magnitude estimates with calibrated numerical model. The comparison shows good agreement. [ABSTRACT FROM AUTHOR]

The article investigates the impact of applying a dynamic computational model that considers inertia forces on pavement deflections under rapidly changing loads over time. This study is particularly relevant to the modelling of falling weight deflectometer (FWD) testing. Initially, the article examines the deflection values obtained from computational models under loads with varying frequencies. In this context, considering inertia forces was significant for load durations shorter than 0.04 s. In such cases, the results of static and dynamic analyses differed considerably. One application of FWD measurement results is determining the stiffness moduli of pavement layers using backcalculation. The study explored the impact of incorporating inertia forces into the pavement model on the estimated values of stiffness moduli obtained via backcalculation. The results revealed differences of several percent between the stiffness moduli calculated using dynamic and static numerical models. Subsequently, the key pavement deformations and fatigue life were determined using the obtained moduli. Again, significantly different results were observed between dynamic and static cases. Based on these findings, it can be concluded that dynamic effects should not be ignored when using FWD testing for backcalculation. Additionally, the article addresses the sensitivity of backcalculation results, which is crucial for the accurate interpretation of the obtained data. [ABSTRACT FROM AUTHOR]

The mechanical behaviour of Fibre-reinforced sands (FRS) has been extensively studied, presenting improved mechanical properties compared to unreinforced soils. Many models have been developed to predict its general stress-strain behaviour. However, the use of double-phase models in FRS is still incipient. Double-phase models are advantageous because they can simulate the whole FRS and the behaviour of its individual components, soil skeleton and reinforcement. This paper uses a modified model for Municipal Solid Waste to reproduce the FRS mechanical response. Introducing a new hardening parameter and a dilatant zone allowed the model to reproduce FRS dilatancy. The model's variables are easily understood, allowing the reproduction of the mechanical behaviour of FRS formed by sands with void ratios ranging from 0.610 to 0.917 and mean grain size from 0.29 to 0.83 mm. The fibres' lengths varied from 12.5 to 51 mm. The results of triaxial and hollow cylinder torsional tests under different stress paths had their main characteristics (peak strength, post-peak behaviour, dilatancy and reinforcement effectiveness) well captured by the model. Predicted and experimental FRS's deviator stress usually differ by less than 15% and the model performance is equivalent or superior to other available models, even requiring fewer input parameters. [ABSTRACT FROM AUTHOR]

SUMMARY: Using forward mantle convection models starting at 140 Ma, and assimilating plate reconstructions as surface velocity boundary condition, we predict present-day mantle structure and compare them with tomography models, using geoid as an additional constraint. We explore a wide model parameter space, such as different values of Clapeyron slope and density change across 660 km, density and viscosity of the thermochemical piles at the core–mantle boundary (CMB), internal heat generation rate, and model initiation age. We also investigate the effects of different strengths of a weak layer below 660 km and weaker asthenosphere and slabs. Our results suggest that slab structures at different subduction zones are sensitive to the viscosity of the asthenosphere, strength of slabs, values of Clapeyron slope and the density and viscosity of the thermochemical piles, while different internal heat generation rates do not affect the slab structures. We find that with a moderately weak asthenosphere ($10^{20}$ Pa·s) and strong slabs, the predicted slab structures are consistent with the tomography models, and the observed geoid is also matched well. Moreover, our models successfully reproduce the degree-2 structure of the lower mantle beneath Africa and the Pacific, also known as Large Low Shear Velocity provinces (LLSVPs). A moderate Clapeyron slope of −2.5 MPa K−1 at 660 km aids in slab stagnation while higher values result in massive slab accumulation at that depth, ultimately leading to slab avalanches. We also find that the convective patterns in the thermal and thermochemical cases with slightly denser LLSVPs are similar, although the geoid amplitudes are lower for the latter. However, with more dense LLSVPs, the slabs cannot perturb them and no plumes are generated. Plumes arise as thermal instabilities from the edges of the LLSVPs, when cold and viscous slabs perturb them. While our predicted plume locations are consistent with the observed hotspot locations, matching the plume structures in tomography models is difficult. These plumes are essential in fitting the finer features of the observed geoid. In longer-duration models, more voluminous subducted material reaches the CMB, which tends to erode the LLSVPs significantly, and yields a poor fit to the observed geoid. Our results suggest that with the presence of a thin, moderately weak layer below 660 km, a slightly dense LLSVP, and Clapeyron slope of −2.5 MPa K−1, the velocity anomalies in seismic tomography and the long-wavelength geoid can be matched well. One of the limitations of our models is that the assimilated plate motion history may be too short to overcome arbitrary initial conditions effects. Also, assimilated true plate velocities in our models may not represent the true convective vigour of the Earth. [ABSTRACT FROM AUTHOR]

SUMMARY: Large-scale and high-resolution seismic modelling are very significant to simulating seismic waves, evaluating earthquake hazards and advancing exploration seismology. However, achieving high-resolution seismic modelling requires substantial computing and storage resources, resulting in a considerable computational cost. To enhance computational efficiency and performance, recent heterogeneous computing platforms, such as Nvidia Graphics Processing Units (GPUs), natively support half-precision floating-point numbers (FP16). FP16 operations can provide faster calculation speed, lower storage requirements and greater performance enhancement over single-precision floating-point numbers (FP32), thus providing significant benefits for seismic modelling. Nevertheless, the inherent limitation of fewer 16-bit representations in FP16 may lead to severe numerical overflow, underflow and floating-point errors during computation. In this study, to ensure stable wave equation solutions and minimize the floating-point errors, we use a scaling strategy to adjust the computation of FP16 arithmetic operations. For optimal GPU floating-point performance, we implement a 2-way single instruction multiple data (SIMD) within the floating-point units (FPUs) of CUDA cores. Moreover, we implement an earthquake simulation solver for FP16 operations based on curvilinear grid finite-difference method (CGFDM) and perform several earthquake simulations. Comparing the results of wavefield data with the standard CGFDM using FP32, the errors introduced by FP16 are minimal, demonstrating excellent consistency with the FP32 results. Performance analysis indicates that FP16 seismic modelling exhibits a remarkable improvement in computational efficiency, achieving a speedup of approximately 1.75 and reducing memory usage by half compared to the FP32 version. [ABSTRACT FROM AUTHOR]

SUMMARY: We propose a novel method for 3-D magnetotelluric (MT) forward modelling based on hybrid meshless and finite-element (FE) methods. This method divides the earth model into a central computational region and an expansion one. For the central region, we adopt scatter points to discretize the model, which can flexibly and accurately characterize the complex structures without generating unstructured mesh. The meshless method using compact support radial basis function is applied to simulate this area's electromagnetic field. While in the expansion region, to avoid the heavy time consumption and numerical error of the meshless method caused by non-uniform nodes, we adopt a node-based finite-element method with regular hexahedral mesh for stability. Finally, the two discretized systems are coupled at the interface nodes according to the continuity conditions of vector and scalar potentials. Considering that the normal electric field is discontinuous at the interface with resistivity discontinuity, while the shape functions for the meshless method are continuous, we further adopt the visibility criterion in constructing the support region. Numerical experiments on typical models show that using the same degree of freedom, the hybrid meshless-finite element method (FEM) algorithm has higher accuracy than the node-based FEM and meshless method. In addition, the electric field discontinuity at interfaces is well preserved, which proves the effectiveness of the visibility criterion method. In general, compared to the conventional grid-based method, this new approach does not need the complex mesh generation for complex structures and can achieve high accuracy, thus it has the potential to become a powerful 3-D MT forward modelling technique. [ABSTRACT FROM AUTHOR]

In order to predict the ground settlement in a scientific, intuitive, and simple way, based on the theory of Bio-consolidation, a three-dimensional fluid-solid coupled numerical calculation programme FGS-3D for ground settlement was compiled by using the Fortran 95 language, and a front-end operation platform was developed by using Microsoft VisualBasic, so that a three-dimensional variable-parameter fully coupled viscoelastic-plastic model of ground settlement was constructed using the city of Yancheng as an example, and the development of ground settlement and horizontal displacement changes from 2021 to 2030 were predicted. The results show that the three-dimensional fully coupled finite-element numerical model of building load, groundwater seepage, and soil deformation established by the above computer development program can directly create a hydrogeological conceptual model of groundwater mining and predict ground settlement, so as to achieve the visualisation of the three-dimensional seepage of groundwater and the fully coupled simulation of ground subsidence in the whole process of groundwater mining. Under the joint action of construction load and groundwater mining, the water level of the aquifer in Yancheng City rises by 1.26 m on average in the main groundwater mining area of the group III pressurised aquifer, forming two smaller landing funnels, and the lowest water level of the two landing funnels is −15 m. [ABSTRACT FROM AUTHOR]

This paper investigates the effect of innovative cross-sectional configurations on Low-Cost Bamboo Composite (LCBC) structural members. The study employs both experimental and numerical methods with different resin matrices and bamboo species. In this study, LCBC short columns are designed with different innovative cross-sectional configurations in an attempt to overcome the costly production processes of engineered bamboo. This approach uses bundles of bamboo, both in culm and strip forms. A compatible, environmentally responsible, and economically justifiable resin matrix is used to fabricate an LCBC member. The production of LCBC members does not necessitate highly advanced technology. This capability enables the production of LCBC members in custom-designed cross-sectional shapes and lengths. This study introduces the Russian doll (RD), Big Russian doll (BRD), Hawser (HAW), and Scrimber (SCR) cross-sectional configurations. Extra-large, large, medium, and small sizes of bamboo are employed. Synthetic Epoxy (EXP), a Bio-based Experimental soft filler (BE1), Bio-Epoxy (BE2), Furan-based (PF1) matrices are applied. Furthermore, Moso, Guadua, Madake, and Tali bamboo species are incorporated. The results of this study reveal that the most efficient cross-sectional configuration for compressive strength is the HAW configuration, closely followed by the SCR configuration. LCBC members with bio-resins have shown excellent promise in competing in strength with those made with their synthetic counterparts. The maximum compressive strengths (MPa) were achieved by two specimens with synthetic epoxy closely followed by a specimen with bio-epoxy, namely HAW-EPX-M, RD-EPX-M, and RD-BE2-G specimens with 78 MPa, 75 MPa, and 72 MPa, respectively. In terms of the modulus of elasticity of LCBC with different resin matrices, the stiffest specimens were HAW-BE2-M1, HAW-EPX-M, and HAW-BE2-M2 with 3.89 GPa, 3.08 GPa, and 2.54 GPa, respectively. The theoretical and numerical modelling of the LCBC members showed excellent correlation with the experimental results, which provides the capacity to design LCBC for engineering projects. The LCBC design can be further developed with more bamboo and less resin content. [ABSTRACT FROM AUTHOR]

Columns in traditional timber structures are commonly seen resting on stone bases and are very important in resisting lateral loads. This paper numerically modeled the lateral resistance of the columns combined with experimental investigation. Different numerical models were developed, based on which sensitivity analyses were performed. A practical numerical modelling strategy was further proposed and verified. The analysis results indicated that instead of material properties and contact area, the columns' lateral performance was much more sensitive to the variation of surface curvature at bottom surface. Without consideration of the surface curvature, the modeling error in the initial stiffness and peak load of a column was more than 607% and 8%, respectively. By best matching the load–displacement curves of tested column specimens, the optimal surface curvature was identified as 1/20306.5 mm−1. Then, a practical finite numerical model, characterized by the optimal curvature at the bottom surface, was proposed. This modelling method was validated by use of existing test results of traditional timber columns with different diameters and vertical compression loads. The modelling load–displacement curves agreed well with experimental curves both in terms of the initial lateral stiffness and peak load. Detailed simulation results based on the practical modelling strategy were presented and discussed. [ABSTRACT FROM AUTHOR]

This paper pertains to the reliability analysis of transmission tower foundations under variable wind loads. The reliability analysis of the anchor foundation, frequently employed for transmission towers due to their huge uplift capacity at shallow depth, is presented using a systematic probabilistic approach. The paper presents the design of an anchor foundation for a typical transmission tower based on loads obtained by applying lateral wind velocity on the tower. The Multivariate Adaptive Regression Splines (MARS) surrogate model is first developed to construct a deterministic prediction model for the anchor uplift capacity. The developed model is further used in the reliability analysis of the anchor foundation considering variabilities in the soil properties and the uplift load, and the results of the probabilistic analysis are reported in the paper. The results indicate the importance of incorporating variability in design variables for the proper safety of the structure. Soil cohesion is found to be the most sensitive parameter in the uplift capacity of the tower foundation. The comparison of the deterministic factor of safety and failure probability of the foundation brings out the importance of the reliability-based design methodology. [ABSTRACT FROM AUTHOR]

Enriching methane from sources such as natural gas, biogas, and coalbed methane is a valuable gas separation endeavor, given methane as a cleaner and more sustainable fuel resource. Pressure swing adsorption (PSA) is widely regarded as a mature technology for this purpose. The essence of the mathematical model is to capture the fundamental aspects of the problem, enabling an accurate description of the system. This review summarizes recent trends in the mathematical modelling of PSA for methane enrichment. It provides a comprehensive compilation of numerical modelling and simulation studies, with a critical emphasis on methane enrichment applications. A rigorous mathematical model describing the adsorption phenomenology in a packed bed column is presented, incorporating mass, momentum, and energy balances, which serve as the fundamental framework for the cyclic adsorption process. Key analyses and performance metrics of PSA systems are detailed through examples from the literature. A case study based on a single column is presented to exemplify the mathematical modelling pathway and serve as a basis for model validation. The paper concludes by highlighting instructive avenues for future research directions and development in this area. [ABSTRACT FROM AUTHOR]

Most studies on open-pit mine stability have been performed in two-dimensional (2D) mode, which neglects the effect of the circular shape and thus the effect of third-direction boundary conditions. The aim of this work was to study the effect of the third dimension on static and dynamic stability analyses of pit slopes by comparing results from two- and three-dimensional (3D) numerical and limit equilibrium analyses. Geomechanical data from Choghart open-pit mine, located in the seismically active area of central Iran, were considered under a historical earthquake (M = 6.4). Kinematic analyses showed that two tectonic blocks of the mine – BL1 and BL2 – are susceptible to wedge and planar failures, respectively. The 2D and 3D static analyses showed that both blocks are stable, with higher safety factors (SFs) obtained using numerical 3D analyses. In the dynamic analyses, both 2D and 3D analyses showed wedge instability for BL1. However, for BL2, the 2D numerical analyses showed instability whereas the 3D numerical and pseudo-static analyses indicated block stability. The results revealed that the 3D SFs are generally 20–40% greater than those from 2D analyses. This can be explained by the consideration of the third-direction effect, verifying the importance of 3D analyses. [ABSTRACT FROM AUTHOR]

In this study, numerical analyses of connected and disconnected piled-raft foundations (PRFs) comprising a single pile, 2 × 1, 3 × 1, 2 × 2, 3 × 2, 3 × 3 and 5 × 5 pile groups subjected to different lateral loads, with a constant vertical load, have been carried out by three-dimensional finite-element analysis. Results from the analysis have been validated with reported field studies considering vertical load and reported experimental studies considering lateral load. The reported multilayered sub-soil profile has been considered for the analysis. Parametric studies in terms of the effect of variation in pile length (L/d = 20, 30, 39), spacing of piles in the pile group (S/d = 3, 5, 6), number of piles, thickness of raft (0.5, 1.2 and 2 m), variation of axial force, excess pore water pressure, bending moment and shear force along the depth of the piles for connected and disconnected PRFs have been discussed for a time period of 20 years. The horizontal connected piled-raft coefficients have been determined corresponding to the lateral deflections. A generalised equation has been predicted for the horizontal connected PRF coefficient. [ABSTRACT FROM AUTHOR]

This study applies random field theory to simulate the spatial variability of the ground where a tunnel was excavated in a granitic soil mass, in the city of Porto (Portugal). Numerical analyses are performed using the finite-element method combined with random field theory to simulate tunnel excavation. A parametric study to analyse the influence of the correlation coefficient is performed. The results are compared with the measured values and the values obtained with the deterministic computation. The performance of the models is evaluated using two measures of error: root mean squared error (RMSE) and mean absolute deviation (MAD). With the random fields approach (correlation length equal to 0.5) values of MAD equal to 2.18 and RMSE equal to 4.16 were obtained, whereas with the deterministic approach values of MAD equal to 2.29 and RMSE equal to 4.26 were obtained. Furthermore, the performance of the models improves as far as the correlation coefficient increases. Therefore, it is concluded that the results obtained using random fields are better than those obtained using the deterministic analysis. [ABSTRACT FROM AUTHOR]

Underground transportation systems often involve multiple tunnels constructed closely together. Previous studies mainly focus on interaction between circular tunnels; by contrast, interaction mechanisms involving non-circular tunnels are not well understood. In this study, four physical three-dimensional centrifuge tests were performed in dry sand, simulating the response of existing circular and horseshoe-shaped tunnels to a newly excavated tunnel. Two different ratios between pillar depth and tunnel diameter (P/D) of 0.5 and 2.0 were considered. Furthermore, three-dimensional numerical back-analyses considering small-strain stiffness were undertaken. Results reveal that the ground settlement above an existing horseshoe-shaped tunnel is less sensitive to pillar depth than for circular ones. Furthermore, for P/D = 0.5, the existing horseshoe-shaped tunnel experiences both vertical and horizontal compression; more stress reduction occurs vertically than horizontally. A circular tunnel for the same pillar depth becomes compressed vertically but elongated horizontally; stress reduction around the existing circular tunnel is less vertically than horizontally. However, for P/D = 2.0, both types of tunnel become elongated vertically and compressed horizontally because of a larger reduction in vertical stresses than horizontal ones. These results demonstrate that both pillar depth and shape profoundly affect tunnel deformation mechanisms. [ABSTRACT FROM AUTHOR]

The ballast control of a floating dock mainly relies on manual operations, which can be time-consuming and requires skilled workers. This study proposes an automatic ballast control system for floating docks, which improves operational efficiency and safety during the vessel docking process. A numerical model is developed to simulate the dynamic process of the floating dock's operations, which includes a six degrees-of-freedom (6DOF) model, a hydrostatic force model, a hydrodynamic force model, and a hydraulic model. The hydrostatic force model is developed using the Archimedes law and a strip theory along the longitudinal direction. The hydrodynamic model is made based on the effects of added mass and dynamic damping. The hydraulic model is proposed to deal with the hydraulic calculation of the ballast water system. The present automatic ballast control is designed based on a modified proportional controller (P-controller) to control the valve opening angles when the pitch or roll angles are larger than the corresponding threshold values. Without using controllers, the roll angles of the dock can reach 8.9 deg and 13 deg during the ballasting and de-ballasting operations, respectively. The present modified P-controller with optimized control parameters can stabilize the dock during the de-ballasting operation and keep the maximum pitch and roll angles no larger than 0.016 deg and 0.0783 deg, respectively. During the ballasting operation with the same control parameters, the roll and pitch are below 0.0604 deg and 0.0145 deg, respectively. The present automatic control will be further implemented in the vessel docking cases and can significantly improve the stability of the dock. [ABSTRACT FROM AUTHOR]

SUMMARY: Reconciling rock unit boundary geometry is crucial for geological and geophysical studies aiming to achieve a comprehensive 3-D subsurface model. To create a unified 3-D parametrization suitable for both geological modelling and geophysical inversion, an integrated approach utilizing implicit modelling is essential. However, a key challenge lies in encapsulating all pertinent information within the 3-D model, ensuring compatibility with the utilized data sets and existing constraints. In this study, we present a workflow that enables the generation of an integrated 3-D subsurface model primarily using gravity and reflection seismic data sets. Our approach involves a cooperative geophysical inversion workflow, which incorporates the inverted model from the reflection seismic data while leveraging sparse petrophysical information. Despite advances in integrated modelling, the incorporation of implicit modelling approaches in cooperative inversion workflows remains unexplored. In our gravity inversion process, we use a generalized level set method to refine the boundaries of rock units in the prior model. We integrate the inverted model, derived from seismic and other sparse petrophysical data sets, to create a comprehensive 3-D prior model. To enhance the integration of reflection seismic data sets in the level set inversion, we introduce a weighting uncertainty matrix containing constraint terms. This step refines the model's accuracy and ensures greater consistency. Finally, we search for any missing rock units within inverted model through nucleation investigations. The introduced methodology has undergone successful testing in the Boulia region (Southern Mount Isa, Queensland), utilizing two 2-D reflection seismic profiles and regional gravity data sets. This study primarily aims to reconstruct the geometry of major structures within the basement units and the basin at a regional scale. By combining seismic profiles and gravity data sets with constraining information, we are able to create a 3-D model of the area that accurately represents distinct rock units and their boundary geometries. Additionally, relevant legacy data sets and prior modelling results from the region have been incorporated and refined, ensuring that the final model aligns with all available knowledge about the area. [ABSTRACT FROM AUTHOR]

SUMMARY: An accurate magnitude estimation is necessary to properly evaluate seismic hazard, especially in low to moderate seismicity areas such as Metropolitan France. However, magnitudes of small earthquakes are subject to large uncertainties caused by major high-frequency propagation effects which are generally not properly considered. To address this issue, we developed a method to separate source, attenuation and site parameters from the elastic radiative transfer modelling of the full energy envelopes of seismograms. The key feature of our approach is the treatment of attenuation—both scattering and absorption—in a simple but realistic velocity model of the Earth's lithosphere, including a velocity discontinuity at the Moho. To reach this goal, we developed a 2-step inversion procedure, allowing first to extract attenuation parameters for each source-station path from the whole observed energy envelope using the Levenberg–Marquardt and grid-search algorithms, then to determine site amplification and the source displacement spectrum from which the moment magnitude Mw is extracted. In the first step, we use the forward modelling procedure of Heller et al. in order to simulate energy envelopes by taking into account the full treatment of wave polarization, the focal mechanism of the source and the scattering anisotropy. The inversion procedure is then applied to the 2019 ML 5.2 Le Teil and 2014 ML 4.5 Lourdes earthquakes which both occurred in southern France. Data from 6 stations are selected for each event. The inversion results confirm a significant variability in the attenuation parameters (scattering and absorption) at regional scale and a strong frequency dependence. Scattering appears to be stronger towards the French Alps and Western Pyrenees. Absorption is stronger as frequency increases. Although not very resolvable, the mechanism of scattering appears to be forward or very forward. By inverting the source spectrum, we determine moment magnitudes Mw of 5.02 ± 0.17 for the Le Teil earthquake and 4.17 ± 0.15 for the Lourdes earthquake. [ABSTRACT FROM AUTHOR]

SUMMARY: A 3-D finite-difference solver has been developed and implemented for Boussinesq convection in a spherical shell. The solver transforms any complex curvilinear domain into an equivalent Cartesian domain using Jacobi transformation and solves the governing equations in the latter. This feature enables the solver to account for the effects of the non-spherical shape of the convective regions of planets and stars. Apart from parallelization using MPI, implicit treatment of the viscous terms using a pipeline alternating direction implicit scheme and HYPRE multigrid accelerator for pressure correction makes the solver efficient for high-fidelity direct numerical simulations. We have performed simulations of Rayleigh–Bénard convection at two Rayleigh numbers Ra = 105 and 107 while keeping the Prandtl number fixed at unity (Pr = 1). The average radial temperature profile and the Nusselt number match very well, both qualitatively and quantitatively, with the existing literature. Closure of the turbulent kinetic energy budget, apart from the relative magnitude of the grid spacing compared to the local Kolmogorov scales, ensures sufficient spatial resolution. [ABSTRACT FROM AUTHOR]

SUMMARY: Geophysics and Geotechnical Engineering commonly use 1-D wave propagation analysis, simplifying complex scenarios by assuming flat and homogeneous soil layers, vertical seismic wave propagation and negligible pore water pressure effects (total stress analysis). These assumptions are commonly used in practice, providing the basis for applications like analysing site responses to earthquakes and characterizing soil properties through inversion processes. These processes involve various in situ tests to estimate the subsurface soil's material profile, providing insights into its behaviour during seismic events. This study seeks to address the limitations inherent to 1-D analyses by using 3-D physics-based simulations to replicate in situ tests performed in the Argostoli basin, Greece. Active and passive source surveys are simulated, and their results are used to determine material properties at specific locations, using standard geophysical methods. Our findings underscore the potential of 3-D simulations to explore different scenarios, considering different survey configurations, source types and array sets. [ABSTRACT FROM AUTHOR]

Lisbon's historical seismicity, socioeconomic importance and population density contribute to a moderate to high seismic risk. The geological setting of the city includes cases of inclined layers, interbedding sedimentary rock layers in soil deposits, sand and clay layers in the same geological unit, leading to cases of shear wave velocity inversion and a large scatter of geotechnical properties within each geological unit. The morphological setting of the city is characterised by the existence of several hills and relatively shallow, stream-carved valleys filled with alluvial deposits. The seismic site effects in Lisbon were assessed through numerical simulation using the linear equivalent method and adopting the two types of seismic action defined in the Portuguese National Annex of Eurocode 8: (i) one-dimensional subsoil models covering the city, at sites where borehole data and geophysical data were available; (ii) two-dimensional subsoil models along three cross-sections representative of the geological settings and morphology. The distribution of amplification factors in the city revealed a pattern related to ground characteristics that impact seismic soil response, such as the presence of high-thickness cover deposits, significant shear-wave variations, alluvial valleys, a crest or significant slope variations and inclined layers. The 2D/1D spectral ratio highlighted the areas were 2D seismic effects are more important. The soil factor determined in the numerical analyses was consistently greater than the soil factor values indicated in Eurocode 8. [ABSTRACT FROM AUTHOR]

The electrical transport of saturated porous rocks has been studied in 3-dimensions for five different channel geometries. The five selected geometries can form the de-constructed parts of any real transport channel in reservoir rocks. The average cross-sectional area and effective length of each of the channels have been computed to construct the geometric factor. The electrical properties studied showed a quadratic increase with increase in clay fraction in the channel with surface conductivity emerging as the major contributor responsible for conductivity at clay fraction ∼ 0.5 and above. The role of clay conductance in electrical transport was investigated for the different geometries. Variations in pore throat size corresponded to non-linear shifts in effective conductivity with clay conductance, in contrast to the linear increase observed for uniform pore throat size. Channel structures with inverse gradients in flow paths, showed significantly lower electrical conductivity. Structures with more than one parallel path displayed higher electrical conductivity. [ABSTRACT FROM AUTHOR]

Behind a retaining wall, the mean effective stress mainly decreases during an excavation phase following an unloading stress path. The volumetric strains generated by purely elastic soils are systematically dilative which induces aberrant ground uplifts. The introduction of plasticity along with a nonlinear elastic domain turns out to be essential for a realistic prediction of ground movements. In this paper, a numerical analysis is carried out using a finite element code considering an advanced soil constitutive model called "Generalized Hardening Soil" which has been recently developed. This model contains the exact same set of features as the Hardening Soil Small Strain model but with the possibility to activate each of its plastic and nonlinear elastic mechanisms independently. The role of these mechanisms are investigated to assess their impact on the shape and the amplitude of the ground movements. Numerical results demonstrated that plasticity triggers the main contractive volumetric strains leading to settlements. Nevertheless it cannot fully compensate the elastic uplifts due to unloading. The insertion of strain dependent stiffness is essential as well as the stress dependency. A back analysis of the historical excavation of the Taipei National Enterprise Center permitted to validate these findings. [ABSTRACT FROM AUTHOR]

Legacy colliery spoil tip failures pose a significant hazard that can result in harm to persons or damage to property and infrastructure. In this research, the 2020 Wattstown tip landslide caused by heavy rainfall was examined to investigate the likely mechanisms and developmental factors contributing to colliery spoil tip failures in Welsh coalfields. To achieve this, an integrated method was proposed through the combination of remote sensing mapping, stability chart analysis, 2D limit equilibrium (LE) modelling, and 3D finite difference method (FDM) analysis. Various water table geometries were incorporated into these models to ascertain the specific groundwater condition that triggered the occurrence of the 2020 landslide. In addition, sensitivity analyses were carried out to assess the influence of the colliery spoil properties (i.e., cohesion, friction angle, and porosity) on the slope stability analysis. The results indicate that the landslide was characterised by a shallow rotational failure mode and spatially constrained by the critical water table and an underlying geological interface. In addition, the results also imply that the landslide was triggered by the rise of water table associated with heavy rainfall. Through sensitivity analysis, it was found that the properties of the colliery spoil played an important role in confining the extent of the landslide and controlling the process of its development. The findings underscore the detrimental effects of increased pore pressures, induced by heavy rainfall, on the stability of colliery tips, highlighting the urgent needs for local government to enhance water management strategies for this region. [ABSTRACT FROM AUTHOR]

Liquefaction-induced settlement of shallow foundations is the result of bearing capacity failure in undrained conditions and sedimentary settlement during the post-liquefaction process. The bearing capacity of a shallow foundation is highly dependent on the size and dimensions of its footprint. In addition, the reduction in shear strength in liquefiable soil, a key parameter for estimating bearing capacity, depends on the excess pore water pressure generated during an earthquake. This study aims to investigate the impact of earthquake motion on the extent of liquefaction-induced settlement in shallow foundations. A parametric study was conducted by varying the input earthquake motions in a three-dimensional response history analysis to directly consider the interaction between the soil and superstructures. The numerical analysis model constructed for the parametric study was rigorously calibrated using a reference dynamic centrifuge test in a prototype scale. The effects of the horizontal boundary and drainage conditions in the numerical model were closely examined during calibration. The parametric study results indicate that the intensity measures of an earthquake, which quantify the energy associated with the number of reversals, exhibit a close correlation with the resulting liquefaction-induced settlement as opposed to other conventional earthquake motion parameters, such as peak acceleration, magnitude, and frequency. [ABSTRACT FROM AUTHOR]

Integrating ice‐sheet models with empirical data pertaining to palaeo‐ice sheets promotes advances in the models used in sea‐level predictions and can improve our understanding of past ice‐sheet behaviour. The large number of empirical constraints on the last British–Irish Ice Sheet make it ideal for model–data comparison experiments. Here, we present an ensemble of 600 model simulations, which are compared with data on former ice‐flow extent, flow geometry and deglaciation timing. Simulations which poorly recreate data were ruled out, allowing us to examine the remaining physically realistic simulations which capture the ice sheets' behavioural tendencies. Our results led to a novel reconstruction of behaviour in the data‐poor region of the North Sea, insights into the ice stream, potential ice‐shelf and readvance dynamics, and the potential locations of peripheral ice caps. We also propose that the asynchronous behaviour of the British–Irish Ice Sheet is a consequence of the geography of the British Isles and the merging and splitting of different bodies of ice through saddle merger and collapse. Furthermore, persistent model–data mismatches highlight the need for model development, especially regarding the physics of ice–ocean interactions. Thus, this work highlights the power of integrating models and data, a long‐held aim of palaeoglaciology. [ABSTRACT FROM AUTHOR]

Geothermal energy, derived from the Earth's internal heat, can be harnessed due to the geothermal gradient between the Earth's interior and its surface. This heat, sustained by radiogenic decay, varies across regions, and is highest near volcanic areas. In 2020, 108 countries utilised geothermal energy, with an installed capacity of 15,950 MWe for electricity and 107,727 MWt for direct use in 2019. Low-enthalpy sources require binary systems for power production. Open-loop systems face issues like scaling, difficult water treatment, and potential seismicity, while closed-loop systems, using abandoned petroleum or gas wells, reduce costs and environmental impacts greatly. The novel geothermal gravity heat pipe (GGHP) design eliminates parasitic power consumption by using hydrostatic pressure for fluid circulation. Implemented in an abandoned well in north-east (NE) Slovenia, the GGHP uses a numerical finite difference method to model heat flow. The system vaporises the working fluid in the borehole, condenses it at the surface, and uses gravitational flow for circulation, maintaining efficient heat extraction. The model predicts that continuous maximum capacity extraction depletes usable heat rapidly. Future work will explore sustainable heat extraction and potential discontinuous operation for improved efficiency. [ABSTRACT FROM AUTHOR]

Full and partial set-in weld-joints connect the nozzle to pressure vessels to prevent lamellar tearing. This investigation aims to compare full and partial set-in weld-joints through computer simulations, numerical modelling, welding experiments, nondestructive examination and destructive testing. The nozzle-to-shell opening connection's weld-joint design model is DN20, with a 50 mm shell wall thickness. The optimal groove weld depth, at a safety factor of 3.5, ranges from 20 to 50 mm. We investigated weld-joint design variations using design-of-experiment, which included actual welding and mechanical testing. The inter-run temperature and its interaction with the heat input were found to be statistically significant. The mechanical properties of welded joint are better than those of the specified material, A516 Grade 70, which include an average impact toughness of 89.5 J at -50 °C, a microhardness of 195 Hv10, an ultimate tensile strength of 533 MPa, a yield strength of 391.5 MPa and an elongation of 56.5%. The coefficients of determination of impact toughness and microhardness, which have R2 values of 88.73% and 76.65%, respectively, conclude a successful welding experiment. The design and mechanical testing results conclude the reliability of the weld-joints. [ABSTRACT FROM AUTHOR]

Copyright of Hydrogeology Journal is the property of Springer Nature and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)

Geophysical electromagnetic survey methods are particularly effective in locating conductive mineral deposits or mineralization zones in a mineral resource exploration. The forward modelling of the electromagnetic responses over such targets is a fundamental task in quantitatively interpreting the geophysical data into a geological model. Due to the ubiquitous irregular and complex geometries associated with the mineral rock units, it is critical that the numerical modelling approach being used is able to adequately and efficiently incorporate any necessary geometries of the Earth model. To circumvent the difficulties in representing complex but necessary geometry features in an Earth model for the existing mesh-based numerical modelling approaches (e.g., finite element and finite difference methods), I present a meshfree modelling approach that does not require a mesh to solve the Maxwell's equations. The meshfree approach utilizes a set of unconnected points to represent any geometries in the Earth model, allowing for the maximal flexibility to account for irregular surface geometries and topography. In each meshfree subdomain, radial basis functions are used to construct meshfree function approximation in transforming the differential equations in the modelling problem into linear systems of equations. The method solves the potential function equations of the Maxwell's equations in the modelling. The modelling accuracy using the meshfree method is examined and verified using one magnetotelluric model and two frequency-domain controlled-source models. The magnetotelluric model is the well-known Dublin Test Model 2 in which the spherical geometry of the conductor in the shallow subsurface may pose as a challenge for many numerical modelling methods. The first controlled-source model is a simple half-space model with the electric dipole source for which analytical solutions exist for the modelling responses. The second controlled-source model is the volcanic massive sulphide mineral deposit from Voisey's Bay, Labrador, Canada in which the deposit's surface is highly irregular. For all modellings, the calculated electromagnetic responses are found to agree with other independent numerical solutions and the analytical solutions. The advantages of the meshfree method in discretizing the Earth models with complex geometries in the forward modelling of geophysical electromagnetic data is clearly demonstrated. [ABSTRACT FROM AUTHOR]

The purpose of this study is to investigate the effect of the position of the solar collector absorber and stagnant air relative to the airflow path on the thermodynamic properties of traversing air and the overall thermal performance of the collector. A numerical model was developed and validated using five different absorber design positions. Numerical results are validated using experimental data obtained from Epinal (France). The model was subsequently used to parametrically study the influence of these absorber positions on the density, temperature, velocity variations of the airflow through the collector and the thermal efficiency of each case scenario. The study showed the importance of stagnant air in solar collector thermal performance. Also, it showed that when one absorber is used, it is better to locate the absorber on the top of the airflow path than below the airflow path. Sandwiching the airflow path between double absorbers with the upper absorbers having stagnant air between it and the plexiglass cover and the down absorber having stagnant air between it and the bottom of the collector can produce a thermal efficiency of 31% which is 2% higher than the closest case and 13% higher than the worst-case design. [ABSTRACT FROM AUTHOR]

The paper is devoted to analysis of power losses in a resistance welding machine including supplying system and examination of welding conditions of the welding machine current in a case of synchronous (simultaneous) operation of multiple welding machines, i.e., during the conduction of welding current. Analysis of the most important contributors of power losses generated on the current path between a power source and the welded joint is carried out. The analysis is carried out for a DC (direct current) welding machine with power electronics inverter. AC (alternating current) welding machines are also taken into account. The analysis is divided in two parts. The first one is the analysis of single welding machine operation, while in the second part, coexistence (mutual operation) of two synchronous welding machines is considered. The analysis is based on results of numerical and experimental investigations. The first one is focused on calculation of power losses in the energy path (including a Sankey type power loss distribution diagrams), and the second one is based on experimental tests carried out to determine the diameter of the weld nugget and the strength of the joints, for cases of reducing the welding current. An example of simultaneous operation of welding machines was presented and discussed. The percentage voltage drops and power losses in the entire power supply path of the resistance welder are shown. The analysis carried out is extremely important from this point quality of welded joints. [ABSTRACT FROM AUTHOR]

Mine tailings are commonly stored in off-stream reservoirs and are usually composed of water with high concentrations of fine particles (microns). The rupture of a mine-tailings pond promotes, depending on the characteristics of the stored material, the fluidization and release of hyper-concentrated flows that typically behave as non–Newtonian fluids. The simulation of non–Newtonian fluid dynamics using numerical modelling tools is based on the solution of mass and momentum conservation equations, particularizing the shear stress terms by means of a rheological model that accounts for the properties of the fluid. This document presents the extension of Iber, a two-dimensional hydrodynamic numerical tool, for the simulation of non–Newtonian shallow flows, especially those related to mine tailings. The performance of the numerical tool was tested throughout benchmarks and real study cases. The results agreed with the analytical and theoretical solutions in the benchmark tests; additionally, the numerical tool also revealed itself to be adequate for simulating the dynamic and static phases under real conditions. The outputs of this numerical tool provide valuable information, allowing researchers to assess flood hazard and risk in mine-tailings spill propagation scenarios. [ABSTRACT FROM AUTHOR]

Assessment of mining impact on groundwater is one of critical considerations for longwall extension and sustainability, however usually constrained by limited data availability, hydrogeological variation, and the complex coupled hydro-mechanical behaviour. This paper aims to determine the factors and mechanism of groundwater depressurisation and identify knowledge gaps and methodological limitations for improving groundwater impact assessment. Analysis of dewatering cases in Australian, Chinese, and US coalfields demonstrates that piezometric drawdown can further lead to surface hydrology degradation, while the hydraulic responses vary with longwall parameters and geological conditions. Statistical interpretation of 422 height of fracturing datasets indicates that the groundwater impact positively correlates to panel geometry and depth of cover, and more pronounced in panel interaction and top coal caving cases. In situ stress, rock competency, clay mineral infillings, fault, valley topography, and surface–subsurface water interaction are geological and hydrogeological factors influencing groundwater hydraulics and long-term recovery. The dewatering mechanism involves permeability enhancement and extensive flow through fracture networks, where interconnected fractures provide steep hydraulic gradients and smooth flow pathways draining the overlying water to goaf of lower heads. Future research should improve fracture network identification and interconnectivity quantification, accompanied by description of fluid flow dynamics in the high fracture frequency and large fracture aperture context. The paper recommends a research framework to address the knowledge gaps with continuous data collection and field-scale numerical modelling as key technical support. The paper consolidates the understanding of longwall mining impacting mine hydrology and provides viewpoints that facilitate an improved assessment of groundwater depressurisation. [ABSTRACT FROM AUTHOR]

In this study, the effects of soil spatial variability on the behaviour of the embankment supported with a combined retaining structure (CRS) were investigated. The numerical model of the CRS embankment was established and validated with the field data. An application programming interface (API) was developed to deal with the data exchanging issue between the numerical model and the spatial variability characterization model. Based on the verified numerical model and the API, the probabilistic analysis with 500 Monte Carlo simulations was automatically computed. Three influencing factors of the retained soil (the mean of the friction angle, the variation of the friction angle and the vertical correlation length of the random field) are analysed by parametric analysis. The results show that the vertical correlation length of the random field is most important in the earth pressure calculation process, while the mean of the friction angle is the factor with least impact. On the whole, the spatial variability of soil properties has minimal impact on the distribution and magnitude of earth pressure behind the retaining structure. [ABSTRACT FROM AUTHOR]

In the context of a global shift towards clean energy, the integration of hydrogen into existing natural gas pipeline networks is essential for facilitating the energy transition. Specifically, low-pressure natural gas pipeline networks are crucial for transporting gas efficiently and with minimal losses from the source to end-users in nearby cities. The physical properties of hydrogen differ significantly from those of natural gas, which implies that introducing hydrogen into the pipeline network will modify operational parameters and impact the safety and efficiency of gas transmission. Therefore, this study concentrates on the efficiency of transporting and scheduling hydrogen-blended natural gas through existing low-pressure direct pipelines. Utilizing the BWRS (Benedict-Webb-Rubin-Starling) equation of state, mathematical models were developed to calculate gas properties and to simulate both the steady and transient states of hydrogen-blended natural gas networks. These models were then validated with TGNET software. This study examines a low-pressure natural gas transmission and distribution network that transports gas directly from a production site in Sichuan to urban areas. It conducts simulations of hydrogen blending and analyzes the impact of various hydrogen blend ratios on the transport and scheduling responses of the network. Research suggests that maintaining the hydrogen blend below 30% ensures the safe operation of low-pressure gas pipelines and prevents the transmission pressure from exceeding 3 MPa. Additionally, under a fixed energy flow output, blending hydrogen intensifies the challenges associated with pressure drops and temperature increases due to high flows in small diameter pipes, and further complicates the management of interruptions resulting from faults. Therefore, it is recommended that the hydrogen blend ratio be maintained at no more than 30%, and that accident handling procedures be established in advance to control pressure variations following a shutdown. • Development of a model for a low-pressure hydrogen-doped pipeline network for direct delivery from the source to the city. • The minimum inlet pressure requirements for different hydrogen doping ratios are presented. • Elevated hydrogen levels can exacerbate temperature rise and pressure drop phenomena in small diameter, high flow pipes. • With the addition of hydrogen, the incident response time for a transient shutdown of transmission was reduced. [ABSTRACT FROM AUTHOR]

Basalt fibre reactive powder concrete (BFRPC) was prepared by incorporating basalt fibres (BFs) instead of steel fibre into reactive powder concrete. In this study, the impact resistance of BFRPC was experimentally and numerically analysed at different strain rates (101–102 s−1). The 75 mm dia. split-Hopkinson pressure bar was used for impact compression tests and the finite-element software LS- Dyna was used for numerical simulation analysis. The results showed that the high-strength BFRPC has an obvious strain rate effect and the dynamic increase factor (DIF) of compressive strength increases logarithmically with strain rate. Meanwhile, the parameters of the Comité Euro-international du Béton (CEB) model were refitted and the relationship between strain rate and DIF was established. By using the Holmquist–Johnson–Cook material constitutive model (HJC model), the stress–strain curves and failure patterns obtained were consistent with the experimental results. The incorporation of BFs significantly improves the deformation properties of BFRPC. [ABSTRACT FROM AUTHOR]

Overpressure‐driven hydrofracturing pervasively occurs in sedimentary basins worldwide. Hydrofracture zones can vertically penetrate several kilometres of rocks and are dominant pathways for basin‐scale fluid migration and energy circulations. Although hydrofracture zones have been extensively described and analysed in the literature, the mechanisms on how hydrofracture zones form and evolve are still poorly understood. In this study, we explore the formation and evolution of a hydrofracture zone in the northern South China Sea, using numerical models constrained by borehole and seismic data. We show that the radius of hydrofracture zone decreases with the strata permeability. The growth of hydrofracture zone is mainly controlled by rock density (ρ$$ \rho $$), pressure at the origin of hydrofracture zone (pb), Poisson's ratio (v), and the radius of the hydrofracture zone at its origin (r). Moreover, as the hydrofracture zone grows, a transition layer forms between the overpressured hydrofracture zone and the overlying hydrostatic pressure zone. The thickness of this transition layer is controlled by strata permeability, strata thickness, overpressure, and pressure gradient within the hydrofracture zone. This study quantitatively explores the development and evolution of overpressure‐driven hydrofractures for the first time, and has wide applications in geohazard assessment, hydrocarbon exploration, carbon circulation, and climatic change. [ABSTRACT FROM AUTHOR]

Seawalls are constructed by covering and protecting the sloping seashore ground with concrete plates or blocks. Their purpose is to sustain high waves induced by strong winds and prevent ground erosion, but they often collapse, mobilising different modes of failure, including that of the ground. Nevertheless, limited research has been conducted on ground failure caused by high waves. In this study, a series of novel centrifuge model tests was first conducted to investigate the failure mechanisms of seawalls due to wave propagation, focusing on the failure of the ground. Finite-element analyses were subsequently conducted to explore the failures observed in the model tests and to provide further insight as to the state of the ground leading to failure. Two failure modes were demonstrated to prevail: floating of the covering panel and sliding failure of the ground. In addition, of the possible causes of failure, the following three were identified in the current study: increased unit weight and reduced suction from wetting; enhanced seepage force under the panel and around the toe block during backwash; water pressure on the back of the panel and the landward side of the toe block during backwash. [ABSTRACT FROM AUTHOR]

The first and fully validated smoothed particle hydrodynamics (SPH) model is presented to tackle coupled flow–deformation problems in unsaturated porous media that undergo large deformation and post-failure behaviour. Unlike the commonly adopted double-layer SPH framework for saturated soils, this paper presents a three-phase single-layer SPH model capable of predicting anisotropic seepage flows through porous media and their complete time-dependent transition from unsaturated to saturated states, as well as their influence on the mechanical behaviour of the porous media and vice versa. The mathematical framework is developed based on Biot's mixture theory and discretised using the authors' recently developed novel SPH approximation scheme for the second derivatives of a field quantity. The soil is modelled using a suction-dependent elastoplastic constitutive model, expressed in terms of effective stress and suction. In addition, an adaptive two-timescale scheme is proposed for the first time to address existing challenges in solving coupled-flow large-deformation problems that involve a significant difference in the timescale required for the solid and fluid phases. The capability of the proposed SPH model was demonstrated through fundamental consolidation tests and a large-scale rainfall-induced slope failure experiment. Very good agreements with theoretical solutions and experimental results are achieved, suggesting that the proposed SPH model can be readily extended to solve a wide range of large-scale geotechnical applications involving coupled unsaturated seepage–deformation problems. [ABSTRACT FROM AUTHOR]

Fracture structures are important natural gas transport spaces in shale gas reservoirs, and their storage state in shale gas reservoirs seriously affects gas production and extraction efficiency. This work uses numerical modeling techniques to investigate the logging response law of the thermal and epithermal neutrons in the gas-containing fracture environment at various angles, applying neutron logging as a technical method. To increase the precision of the evaluation of the natural gas storage condition in shale gas reservoirs, the angle of the fractures' neutron logging data is analyzed. It is found that even in an environment with the same porosity of the fractures, there are significant differences in the logging results due to the different angles of the fracture alignment: 1. the neutron counts in the high-angle (70–90°) fracture environment are 2.25 times higher than in the low-angle (0–20°), but the diffusion area of the neutrons is only 10.58% of that in the low-angle (0–20°); 2. in the neutron energy spectrum, neutron counts are spreading to the high-energy region (7–13 MeV) along with the increase in the angle of the fracture, and the feature is especially prominent in the approximately vertical (60–90°) fracture environment, which is an increase of 528.12% in comparison with the counts in the approximately horizontal angle (0–30°) environment. The main reason for these differences is the variation in the volume of the fracture within the source radiation. This volumetric difference results from the variation in fracture angles (even though the fracture porosity is the same). In view of the above phenomenon, this paper proposes the concept of "effective fracture volume", which can intuitively reflect the degree of influence of fracture angle on neutron logging results. Further, based on the unique characteristics of shale gas reservoirs and neutrons, this paper provides important theoretical support for the modification of the porosity of the field operation, the evaluation of the physical characteristics of the gas endowment space, and the assessment. [ABSTRACT FROM AUTHOR]

The use of advanced materials, such as ultra-high-performance concrete (UHPC), and slender cross-sectional shapes in bridge girders can result in high susceptibility to shear failure during fire exposure. Shear limit state is usually not considered in fire design, which can be a major risk for I-shaped girders with slender webs. This paper presents the development of an approach to evaluate the shear capacity of UHPC bridge girders during fire through extension of room temperature capacity equations to elevated temperatures. Four types of standard AASHTO concrete bridge girders are analyzed under both standard and hydrocarbon fire scenarios. The output results are utilized to evaluate the progression of temperatures and deflections throughout fire exposure, as well as estimate degradation in flexural and shear capacities. Results show that the fire resistance of UHPC bridge girders are significantly affected by fire severity and concrete mass and that these factors should be taken into consideration in fire design of UHPC bridge girders. [ABSTRACT FROM AUTHOR]

This article presents research results regarding the development of a new manufacturing technology for an element assigned to belt conveyor flights in the extractive industry through hot die forging (of a forging with a double-sided flange) instead of the currently realized process of producing such an element by welding two flanges onto a sleeve or one flange onto a flange forging. The studies were conducted to design an innovative and low-waste technology, mainly with the use of numerical modelling and simulations, partially based on the current technology of producing a flange forging. Additionally, during the development of the forging process, the aspect of robotization was considered, both in respect of the forging tools and the process of transportation and relocation of forging between the impressions and the forging aggregates. A thermo-mechanical model of the process of producing a belt conveyor flight forging with deformable tools was elaborated by means of the Forge 3NxT program. The results of the conducted numerical modelling made it possible, among other things, to develop models of forging tools ensuring the proper manner of material flow and filling of the impressions, as well as temperature and plastic deformation distributions in the forging and also the detection of possible forging defects. For the technology elaborated this way, the tools were built together with a special instrument for flanging in the metal, and technological tests were performed under industrial conditions. The produced forgings were verified through a measurement of the geometry, by way of 3D scanning, as well as the hardness, which definitively confirmed the properness of the developed technology. The obtained technological test results made it possible to confirm that the elaborated construction, as well as the tool impressions, ensure the possibility of implementing the designed technology with the use of robotization and automatization of the forging process. [ABSTRACT FROM AUTHOR]

Numerous studies have looked at rainfall infiltration as a triggering factor for rainfall-induced landslides. In addition to rainfall characteristics, soil hydraulic properties are recognized to play a crucial role in instability mechanisms. This study proposes a methodology for assessing the soil-water characteristic curve and hydraulic conductivity function through finite-element seepage modelling of physical model tests. The approach is applied to an instrumented, small-scale slope model built with uniformly graded sand with a 30 degree inclination, exposed to homogeneous rainfall until the occurrence of failure. In a first stage, a sensitivity analysis is carried out to assess the influence of the Mualem–van Genuchten model parameters on the slope response. For this purpose, six physically based indicators are utilized to compare the numerical modelling results with the experimentally obtained data regarding the hydraulic response of the slope model. In the subsequent stage, a trial-and-error calibration stage is conducted to determine the set of parameters for which the difference between the numerically and experimentally acquired data is minimized. Ultimately, the suggested methodology facilitates the assessment of the optimal set of hydraulic parameters, to which both the sensitivity analysis and the trial-and-error calibration phase are anchored. The approach has demonstrated its effectiveness in calibrating the unsaturated hydraulic properties of the considered soil, as it properly addresses the physical mechanisms associated with rainfall infiltration in a slope. [ABSTRACT FROM AUTHOR]

The equatorial mixed Rossby–gravity wave (MRGW) is an important contributor to tropical variability. Its excitation mechanism capable of explaining the observed MRGW variance peak at synoptic scales in the troposphere remains elusive. This study investigates wave–mean flow interactions as a generation process for the MRGWs using the TIGAR model, which employs Hough harmonics as the basis of spectral expansion on the sphere, thereby representing MRGWs as prognostic variables. Idealized numerical simulations reveal the interactions between waves emanating from a symmetric tropical heat source and an asymmetric subtropical zonal jet as an excitation mechanism for the MRGWs. The excited MRGWs have variance spectra resembling the observed MRGWs in the tropical troposphere. The mixed Rossby–gravity energy spectrum has a maximum at zonal wavenumbers k=4$$ k=4 $$–5 also in the case of an asymmetric forcing that generates MRGWs across large scales. Effects of wave–wave interactions appear of little importance for the MRGW growth compared with wave–mean flow interactions. Application of the zonal‐mean zonal wind profiles from ERA5 reaffirms the importance of the asymmetry of the zonal mean flow. [ABSTRACT FROM AUTHOR]

The dependence of seismic wave speeds on propagation or polarization direction, called seismic anisotropy, is a relatively direct indicator of mantle deformation and flow. Mantle seismic anisotropy is often inferred from measurements of shear-wave splitting. A number of standard techniques to measure shear-wave splitting have been applied globally; for example, *KS splitting is often used to measure upper mantle anisotropy. In order to obtain robust constraints on anisotropic geometry, it is necessary to sample seismic anisotropy from different directions, ideally using different seismic phases with different incidence angles. However, many standard analysis techniques can only be applied for certain epicentral distances and source–receiver geometries. To search for new ways to detect mantle anisotropy, instead of focusing on the sensitivity of individual phases, we investigate the wavefield as a whole: we apply a 'wavefield differencing' approach to (systematically) understand what parts of the seismic wavefield are most affected by splitting due to seismic anisotropy in the mantle. We analyze differences between synthetic global wavefields calculated for isotropic and anisotropic input models, incorporating seismic anisotropy at different depths. Our results confirm that the seismic phases that are commonly used in splitting techniques are indeed strongly influenced by mantle anisotropy. However, we also identify less commonly used phases whose waveforms reflect the effects of anisotropy. For example, PS is strongly affected by splitting due to seismic anisotropy in the upper mantle. We show that PS can be used to fill in gaps in global coverage in shear-wave splitting data sets (for example, beneath ocean basins). We find that PcS is also a promising phase, and present a proof-of-concept example of PcS splitting analysis across the contiguous United States using an array processing approach. Because PcS is recorded at much shorter distances than *KS phases, PcS splitting can therefore fill in gaps in backazimuthal coverage. Our wavefield differencing results further hint at additional potential novel methods to detect and characterize splitting due to mantle seismic anisotropy. [ABSTRACT FROM AUTHOR]

Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS' numerical modelling. [ABSTRACT FROM AUTHOR]

Modelling seismic wavefields in complex 3-D elastic media is the key in many fields of Earth Science: seismology, seismic imaging, seismic hazard assessment and earthquake source mechanism reconstruction. This modelling operation can incur significant computational cost, and its accuracy depends on the ability to take into account the scales of the subsurface heterogeneities varying. The theory of homogenization describes how the small-scale heterogeneities interact with the seismic waves and allows to upscale elastic media consistently with the wave equation. In this study, an efficient and scalable numerical homogenization tool is developed, relying on the similarity between the equations describing the propagation of elastic waves and the homogenization process. By exploiting the optimized implementation of an elastic modelling kernel based on a spectral-element discretization and domain decomposition, a fully scalable homogenization process, working directly on the spectral-element mesh, is presented. Numerical experiments on the entire SEAM II foothill model and a 3-D version of the Marmousi II model illustrate the efficiency and flexibility of this approach. A reduction of two orders of magnitude in terms of absolute computational cost is observed on the elastic wave modelling of the entire SEAM II model at a controlled accuracy. [ABSTRACT FROM AUTHOR]

Earthquakes with moment magnitude (M w) ranging from 6.5 to 7.0 have been observed to generate sufficiently strong acoustic waves (AWs) in the upper atmosphere. These AWs are detectable in Global Navigation Satellite System satellite signals-based total electron content (TEC) observations in the ionosphere at altitudes ∼250–300 km. However, the specific earthquake source parameters that influence the detectability and characteristics of AWs are not comprehensively understood. Here, we extend our approach of coupled earthquake-atmosphere dynamics modelling by combing dynamic rupture and seismic wave propagation simulations with 2-D and 3-D atmospheric numerical models, to investigate how the characteristics of earthquakes impact the generation and propagation of AWs. We developed a set of idealized dynamic rupture models varying faulting types and fault sizes, hypocentral depths and stress drops. We focus on earthquakes of M w 6.0–6.5, which are considered the smallest detectable with TEC, and find that the resulting AWs undergo non-linear evolution and form acoustic shock N waves reaching thermosphere at ∼90–140 km. The results reveal that the magnitude of the earthquakes is not the sole or primary factor determining the amplitudes of AWs in the upper atmosphere. Instead, various earthquake source characteristics, including the direction of rupture propagation, the polarity of seismic wave imprints on the surface, earthquake mechanism, stress drop and radiated energy, significantly influence the amplitudes and periods of AWs. The simulation results are also compared with observed TEC fluctuations from AWs generated by the 2023 M w 6.2 Suzu (Japan) earthquake, finding preliminary agreement in terms of model-predicted signal periods and amplitudes. Understanding these nuanced relationships between earthquake source parameters and AW characteristics is essential for refining our ability to detect and interpret AW signals in the ionosphere. [ABSTRACT FROM AUTHOR]