Article topics

We develop a PyTorch-based Graph Network Simulator (GNS) that learns physics and predicts the flow behavior of particulate and fluid systems. GNS discretizes the domain with nodes representing a collection of material points and the links connecting the nodes representing the local interaction between particles or clusters of particles. The GNS learns the interaction laws through message passing on the graph. GNS has three components: (a) Encoder, which embeds particle information to a latent graph, the edges are learned functions; (b) Processor, which allows data propagation and computes the nodal interactions across steps; and (c) Decoder, which extracts the relevant dynamics (e.g., particle acceleration) from the graph. We introduce physics-inspired simple inductive biases, such as an inertial frame that allows learning algorithms to prioritize one solution (constant gravitational acceleration) over another, reducing learning time. The GNS implementation uses semi-implicit Euler integration to update the next state based on the predicted accelerations. GNS trained on trajectory data is generalizable to predict particle kinematics in complex boundary conditions not seen during training. The trained model accurately predicts within a 5% error of its associated material point method (MPM) simulation. The predictions are 5,000x faster than traditional MPM simulations (2.5 hours for MPM simulations versus 20 s for GNS simulation of granular flow). GNS surrogates are popular for solving optimization, control, critical-region prediction for in situ viz, and inverse-type problems. The GNS code is available under the open-source MIT license at https://github.com/geoelements/gns.

We present a novel method of using explainability techniques to design physics-aware neural networks. We demonstrate our approach by developing a convolutional neural network (CNN) for solving an inverse problem for shallow subsurface imaging. Although CNNs have gained popularity in recent years across many fields, the development of CNNs remains an art, as there are no clear guidelines regarding the selection of hyperparameters that will yield the best network. While optimization algorithms may be used to select hyperparameters automatically, these methods focus on developing networks with high predictive accuracy while disregarding model explainability (descriptive accuracy). However, the field of Explainable Artificial Intelligence (XAI) addresses the absence of model explainability by providing tools that allow developers to evaluate the internal logic of neural networks. In this study, we use the explainability methods Score-CAM and Deep SHAP to select hyperparameters, such as kernel sizes and network depth, to develop a physics-aware CNN for shallow subsurface imaging. We begin with a relatively deep Encoder-Decoder network, which uses surface wave dispersion images as inputs and generates 2D shear wave velocity subsurface images as outputs. Through model explanations, we ultimately find that a shallow CNN using two convolutional layers with an atypical kernel size of 3x1 yields comparable predictive accuracy but with increased descriptive accuracy. We also show that explainability methods can be used to evaluate the network’s complexity and decision-making. We believe this method can be used to develop neural networks with high predictive accuracy while also providing inherent explainability.

Individual particle rotation and displacement were measured in triaxial tests on transparent sand stabilized with geogrid simulants. The Cellpose U-Net model, originally developed to segment biological cells, was trained to segment images of fused quartz particles. The Score-CAM metric from the field of Explainable AI was used to validate the application of Cellpose to segment particles of fused quartz. These segmented particles were characterized in terms of Fourier shape descriptors and tracked across images. The measured particle displacements in the monotonic triaxial tests correlated with displacement fields from Digital Image Correlation (DIC). In contrast to DIC, the new technique also allows for the measurement of individual particle rotation. The particle rotation measurements were found to be repeatable across different specimens. A state boundary line between probable and improbable particle motions could be identified for a given test based on the measured particle displacements and rotations. The size of the zone of probable motions was used to quantify the effectiveness of the stabilizing inclusions. The results of repeated load tests revealed that the honeycomb inclusions used stabilized the specimens by reducing both particle displacements and rotations.

Inductive heating with high-voltage cables reduces the risk of hydrate formation by raising the temperature of the production fluid in pipelines. Heating the pipeline results in losing a certain fraction of the heat to the surrounding soil through conduction or convection-dominated flow through the soil. However, the amount of heat lost in conduction versus convection and the transition from conduction to convection-dominated heat loss remains unknown. Soil permeability, temperature gradient between cable and mudline, and burial depth influence the mode of heat transfer and the amount of heat lost. We study the dominant mode of heat transfer in pipelines with inductive heating using 2D Finite Difference analysis under different soil and environmental conditions. Low permeability soils primarily exhibit conductive heat transfer, thus losing minimum heat to the surrounding soil. In contrast, convective flow drives a significant fraction of the heat away from the pipeline and towards the ground surface for highly permeable soils, barely heating the fluid in the pipe. We identify a critical Rayleigh-Darcy number of 1 as the controlling value separating conduction and convection-dominated heat transfer. An increase in burial depth deteriorates the heating efficiency in convection-dominated high permeability soils, while it remains unaffected in conduction-dominated low permeability soils.

Visualizing regional-scale landslides is the key to conveying the threat of natural hazards to stakeholders and policymakers. Traditional visualization techniques are restricted to post-processing a limited subset of simulation data and are not scalable to rendering exascale models with billions of particles. In-situ visualization is a technique of rendering simulation data in real-time, i.e., rendering visuals in tandem while the simulation is running. In this study, we develop a scalable N:M interface architecture to visualize regional-scale landslides. We demonstrate the scalability of the architecture by simulating the long runout of the 2014 Oso landslide using the Material Point Method coupled with the Galaxy ray tracing engine rendering 4.2 million material points as spheres. In-situ visualization has an amortized runtime increase of 2\% compared to non-visualized simulations. The developed approach can achieve in-situ visualization of regional-scale landslides with billions of particles with minimal impact on the simulation process.

Compared to most other primates, humans are characterized by a tight fit between the maternal birth canal and the fetal head, leading to a relatively high risk of neonatal and maternal mortality and morbidities. Obstetric selection is thought to favor a spacious birth canal, whereas the source for opposing selection is frequently assumed to relate to bipedal locomotion. An alternative, yet under-investigated, hypothesis is that a more expansive birth canal suspends the soft tissue of the pelvic floor across a larger area, which is disadvantageous for continence and support of the weight of the inner organs and fetus. To test this “pelvic floor hypothesis” we generated a finite element model of the human female pelvic floor and varied its radial size and thickness while keeping all else constant. This allowed us to study the effect of pelvic geometry on pelvic floor deflection (i.e., the amount of bending from the original position) and tissue stresses and stretches. Deflection grew disproportionately fast with increasing radial size, and stresses and stretches also increased. By contrast, an increase in thickness increased pelvic floor stiffness - i.e. the resistance to deformation - which reduced deflection but was unable to fully compensate for the effect of increasing radial size. Moreover, larger thicknesses increase the intra-abdominal pressure necessary for childbirth. Our results support the pelvic floor hypothesis and evince functional trade-offs affecting not only the size of the birth canal but also the thickness and stiffness of the pelvic floor.

Non-invasive subsurface imaging using full waveform inversion (FWI) has the potential to fundamen-tally change engineering site characterization by enabling the recovery of high resolution 2D/3D maps ofsubsurface stiffness. Yet, the accuracy of FWI remains quite sensitive to the choice of the initial start-ing model due to the complexity and non-uniqueness of the inverse problem. In response, we presentthe novel application of convolutional neural networks (CNNs) to transform an experimental seismicwavefield acquired using a linear array of surface sensors directly into a robust starting model for 2DFWI. We begin by describing three key steps used for developing the CNN, which include: selectionof a network architecture, development of a suitable training set, and performance of network training.The ability of the trained CNN to predict a suitable starting model for 2D FWI was compared againstother commonly used starting models for a classic near-surface imaging problem; the identification of anundulating, two-layer, soil-bedrock interface. The CNN developed during this study was able to predictcomplex 2D subsurface images of the testing set directly from their seismic wavefields with an averagemean absolute percent error of 6%. When compared to other common approaches, the CNN approachwas able to produce starting models with smaller seismic image and waveform misfits, both before andafter FWI. The ability of the CNN to generalize to subsurface models which were dissimilar to the onesupon which it was trained was assessed using a more complex, three-layered model. While the predictiveability of the CNN was slightly reduced, it was still able to achieve seismic image and waveform misfitscomparable to the other commonly used starting models. This study demonstrates that CNNs have greatpotential as a tool for developing good starting models for FWI, which are critical for producing accurateFWI results.

A new thixotropic model is developed integrating the Papanastasiou-Bingham model with thixotropy equations to simulate the flow behaviour of Tremie Concrete in the Material Point Method framework. The effect of thixotropy on the rheological behaviour of fresh concrete is investigated by comparing field measurements with numerical simulations. The comparison yields new insights into a critical and often overlooked behaviour of concrete. A parametric study is performed to understand the effect of model parameters and rest-time on the shear stress response of fresh concrete. The Material Point Method with the Papanastasiou-Bingham model reproduces slump-flow measurements observed in the field. The novel model revealed a decline in concrete workability during the Slump-flow test after a period of rest due to thixotropy, which the physical version of the test fails to capture. This reduction in workability significantly affects the flow behaviour and the effective use of fresh concrete in construction operation.

Abstract:

In situ visualization techniques are hampered by a lack of foresight: crucial simulation phenomena can be missed due to a poor sampling rate or insufficient detail at critical timesteps. Keeping a human in the loop is impractical, and defining statistical triggers can be difficult. This paper demonstrates the potential for using a machine-learning-based simulation surrogate as an oracle to identify expected critical regions of a large-scale simulation. These critical regions are used to drive the in situ analysis, providing greater data fidelity and analysis resolution with an equivalent I/O budget to a traditional in situ framework. We develop a distributed asynchronous in situ visualization by integrating TACC Galaxy with CB-Geo MPM for material point simulation of granular flows. We employ a PyTorch-based 3D Graph Network Simulator (GNS) trained on granular flow problems as an oracle to predict the dynamics of granular flows. Critical regions of interests are manually tagged in GNS for in situ rendering in MPM.

Shear stress history controls the pore pressure response in liquefiable soils. The excess pore pressure does not increase under cyclic loading when shear stress amplitude is lower than the peak prior amplitude – the shielding effect. Many sophisticated constitutive models fail to capture the shielding effect observed in the cyclic liquefaction experiments. We develop a data-driven machine learning model based on the LSTM neural network to capture the liquefaction response of soils under cyclic loading. The LSTM model is trained on 12 laboratory cyclic simple shear tests on Nevada sand in loose and dense conditions subjected to different cyclic simple shear loading conditions. The LSTM model features include the relative density of soil and the previous stress history to predict the pore water pressure response. The LSTM model successfully replicates the pore pressure response for three cyclic simple test results considering the shielding and density effects.

The performance evaluation of a potentially unstable slope involves two key components: the initiation of the slope failure and the post-failure runout. The Finite Element Method (FEM) excels at modeling the initiation of instability but quickly loses accuracy in modeling large-deformation problems due to mesh distortion. Hence, the FEM is unable to accurately model post-failure slope runout. Hybrid Eulerian-Lagrangian methods, such as the Material Point Method (MPM), offer a promising alternative for solving large-deformation problems, because particles can move freely across a background mesh, allowing for large deformation without computational issues. However, the use of moving material points in MPM for integration rather than the fixed Gauss points of the FEM reduces the accuracy of MPM in predicting stress distribution and thus failure initiation. We have created a hybrid method by initiating a failure simulation in FEM and subsequently transferring the coordinates, velocities, and stresses to MPM particles to model the runout behavior, combining the strength of both methods. We demonstrate the capability of the hybrid approach by simulating the collapse of a frictional granular column, comparing it to an empirical solution, and evaluating a suitable time to transfer from FEM to MPM by trialing multiple iterations with transfers at different stages of the collapse.

In nature, submarine slope failures usually carry thousands of cubic-meters of sediments across extremely long distances and cause tsunamis and damages to offshore structures. This paper uses the granular column collapse experiment to investigate the effect of slope angle on the runout behavior of submarine granular landslides for different initial volumes. A two-dimensional coupled lattice Boltzman and discrete element method (LBM-DEM) approach is adopted for numerically modeling the granular column collapse. Columns with four different slope angles and six different volumes are modelled under both dry and submerged conditions. The effects of hydrodynamic interactions, including the generation of excess pore pressures, hydroplaning, and drag forces and formation of turbulent vortices, are used to explain the difference in the runout behavior of the submerged columns compared to the dry columns. The results show that at any given slope angle, there is a threshold volume above which the submerged columns have a larger final runout compared to their dry counterpart, and this threshold volume decreases with slope angle.

When power electronics are deployed under the road surface as part of a wireless system it is important to know that their packaging provides adequate heat extraction as well as the required environmental protection – often conflicting requirements. Presently very little can be found in wireless charging standards and literature on the topic of thermal modelling for in-ground components. Yet, this is a topic of great practical significance especially for in-road systems. Traditional cooling methods are not readily applicable underground. This paper uses finite element thermal modelling to investigate the cooling of a representative medium-power in-road wireless system, housed in a sealed ground assembly (GA) chamber and installed to UK requirements (HAUC). The paper quantitatively compares design options and provides practical recommendations for in-road installation thermal management.

Submarine landslides transport thousands of cubic meters of sediment across continental shelves even at slopes as low as 1° and can cause significant casualty and damage to infrastructure. The run-out mechanism in a submarine landslide is affected by factors such as the initial packing density, permeability, slope angle, and initial volume. While past studies have focused on the influence of density, permeability, and slope angle on the granular column collapse, the impact of volume on the run-out characteristics has not been investigated. This study aims to understand how the initial volume affects the run-out using a two-dimensional coupled lattice Boltzman and discrete element (LBM-DEM) method. The coupled LBM-DEM approach allows simulating fluid flow at the pore-scale resolution to understand the grain-scale mechanisms driving the complex continuum-scale response in the granular column collapse. For submerged granular column collapse, the run-out mechanism is heavily influenced by the interaction between the grains and the surrounding fluid. The development of negative pore pressures during shearing and hydrodynamic drag forces inhibit the flow. On the other hand, entrainment of water resulting in hydroplaning enhances the flow. With an increase in volume, the interaction between the grains and the surrounding fluid varies, causing changes in the run-out behavior. For smaller volumes, the forces inhibiting the underwater flow predominates, resulting in shorter run-outs than their dry counterparts. At large volumes, hydroplaning results in larger run-out than the dry cases, despite the inhibiting effects of drag forces and negative pore pressures.

The soil water retention curve (SWRC) is the most commonly used relationship in the study of unsaturated soil. In this paper, the effect of porosity on the SWRC is investigated by numerically modeling unsaturated soil using the Shan-Chen multiphase Lattice Boltzmann Method. The shape of simulated SWRCs are compared against that predicted by the van Genuchten model, demonstrating a good fit except at low degrees of saturation. The simulated SWRCs show an increase in the air-entry value as porosity decreases.

We develop a PyTorch-based Graph Network Simulator (GNS) that learns physics and predicts the flow behavior of particulate and fluid systems. GNS discretizes the domain with nodes representing a collection of material points and the links connecting the nodes representing the local interaction between particles or clusters of particles. The GNS learns the interaction laws through message passing on the graph. GNS has three components: (a) Encoder, which embeds particle information to a latent graph, the edges are learned functions; (b) Processor, which allows data propagation and computes the nodal interactions across steps; and (c) Decoder, which extracts the relevant dynamics (e.g., particle acceleration) from the graph. We introduce physics-inspired simple inductive biases, such as an inertial frame that allows learning algorithms to prioritize one solution (constant gravitational acceleration) over another, reducing learning time. The GNS implementation uses semi-implicit Euler integration to update the next state based on the predicted accelerations. GNS trained on trajectory data is generalizable to predict particle kinematics in complex boundary conditions not seen during training. The trained model accurately predicts within a 5% error of its associated material point method (MPM) simulation. The predictions are 5,000x faster than traditional MPM simulations (2.5 hours for MPM simulations versus 20 s for GNS simulation of granular flow). GNS surrogates are popular for solving optimization, control, critical-region prediction for in situ viz, and inverse-type problems. The GNS code is available under the open-source MIT license at https://github.com/geoelements/gns.

The performance evaluation of a potentially unstable slope involves two key components: the initiation of the slope failure and the post-failure runout. The Finite Element Method (FEM) excels at modeling the initiation of instability but quickly loses accuracy in modeling large-deformation problems due to mesh distortion. Hence, the FEM is unable to accurately model post-failure slope runout. Hybrid Eulerian-Lagrangian methods, such as the Material Point Method (MPM), offer a promising alternative for solving large-deformation problems, because particles can move freely across a background mesh, allowing for large deformation without computational issues. However, the use of moving material points in MPM for integration rather than the fixed Gauss points of the FEM reduces the accuracy of MPM in predicting stress distribution and thus failure initiation. We have created a hybrid method by initiating a failure simulation in FEM and subsequently transferring the coordinates, velocities, and stresses to MPM particles to model the runout behavior, combining the strength of both methods. We demonstrate the capability of the hybrid approach by simulating the collapse of a frictional granular column, comparing it to an empirical solution, and evaluating a suitable time to transfer from FEM to MPM by trialing multiple iterations with transfers at different stages of the collapse.

Visualizing regional-scale landslides is the key to conveying the threat of natural hazards to stakeholders and policymakers. Traditional visualization techniques are restricted to post-processing a limited subset of simulation data and are not scalable to rendering exascale models with billions of particles. In-situ visualization is a technique of rendering simulation data in real-time, i.e., rendering visuals in tandem while the simulation is running. In this study, we develop a scalable N:M interface architecture to visualize regional-scale landslides. We demonstrate the scalability of the architecture by simulating the long runout of the 2014 Oso landslide using the Material Point Method coupled with the Galaxy ray tracing engine rendering 4.2 million material points as spheres. In-situ visualization has an amortized runtime increase of 2\% compared to non-visualized simulations. The developed approach can achieve in-situ visualization of regional-scale landslides with billions of particles with minimal impact on the simulation process.

A new thixotropic model is developed integrating the Papanastasiou-Bingham model with thixotropy equations to simulate the flow behaviour of Tremie Concrete in the Material Point Method framework. The effect of thixotropy on the rheological behaviour of fresh concrete is investigated by comparing field measurements with numerical simulations. The comparison yields new insights into a critical and often overlooked behaviour of concrete. A parametric study is performed to understand the effect of model parameters and rest-time on the shear stress response of fresh concrete. The Material Point Method with the Papanastasiou-Bingham model reproduces slump-flow measurements observed in the field. The novel model revealed a decline in concrete workability during the Slump-flow test after a period of rest due to thixotropy, which the physical version of the test fails to capture. This reduction in workability significantly affects the flow behaviour and the effective use of fresh concrete in construction operation.

Geoelements group wins TACC pathways proposal to simulate the Oso landslide.

The water retention curve (WRC) defines the relationship between matric suction and saturation and is a key function for determining the hydro-mechanical behavior of unsaturated soils. We investigate possible microscopic origins of the water retention behavior of granular soils using both Computed Tomography (CT) experiment and multiphase lattice Boltzmann Method (LBM). We conduct a CT experiment on Hamburg sand to obtain its WRC and then run LBM simulations based on the CT grain skeleton. The multiphase LBM simulations capture the hysteresis and pore-scale behaviors of WRC observed in the CT experiment. Using LBM, we observe that the spatial distribution and morphology of gas clusters varies between drainage and imbibition paths and is the underlying source of the hysteresis. During drainage, gas clusters congregate at the grain surface; the local suction increases when gas clusters enter through small pore openings and decreases when gas clusters enter through large pore openings. Whereas, during imbibition, gas clusters disperse in the liquid; the local suction decreases uniformly. Large pores empty first during drainage and small pores fill first during imbibition. The pore-based WRC shows that an increase in pore size causes a decrease in suction during drainage and imbibition, and an increase in hysteresis.

The soil-water characteristic curve (SWCC) is the most fundamental relationship in unsaturated soil mechanics, relating the amount of water in the soil to the corresponding matric suction. From experimental evidence, it is known that SWCC exhibits hysteresis (i.e. wetting/drying path dependence). Various factors have been proposed as contributors to SWCC hysteresis, including air entrapment, contact angle hysteresis, ink-bottle effect, and change of soil fabric due to swelling and shrinkage, however, the significance of their contribution is debated. From our pore-scale numerical simulations, using the multiphase lattice Boltzmann method, we see that even when controlling for all these factors SWCC hysteresis still occurs, indicating that there is some underlying source that is not accounted for in these factors. We find this underlying source by comparing the liquid/gas phase distributions for simulated wetting and drying experiments of 2D and 3D granular packings. We see that during wetting (i.e. pore filling) many liquid bridges expand simultaneously and join together to fill the pores from the smallest to the largest, allowing menisci with larger radii of curvature (lower matric suction). Whereas, during drying (i.e. pore emptying), only the limited existing gas clusters can expand, which become constrained by the size of the pore openings surrounding them and result in menisci with smaller radii of curvature (higher matric suction).

In nature, submarine slope failures usually carry thousands of cubic-meters of sediments across extremely long distances and cause tsunamis and damages to offshore structures. This paper uses the granular column collapse experiment to investigate the effect of slope angle on the runout behavior of submarine granular landslides for different initial volumes. A two-dimensional coupled lattice Boltzman and discrete element method (LBM-DEM) approach is adopted for numerically modeling the granular column collapse. Columns with four different slope angles and six different volumes are modelled under both dry and submerged conditions. The effects of hydrodynamic interactions, including the generation of excess pore pressures, hydroplaning, and drag forces and formation of turbulent vortices, are used to explain the difference in the runout behavior of the submerged columns compared to the dry columns. The results show that at any given slope angle, there is a threshold volume above which the submerged columns have a larger final runout compared to their dry counterpart, and this threshold volume decreases with slope angle.

Submarine landslides transport thousands of cubic meters of sediment across continental shelves even at slopes as low as 1° and can cause significant casualty and damage to infrastructure. The run-out mechanism in a submarine landslide is affected by factors such as the initial packing density, permeability, slope angle, and initial volume. While past studies have focused on the influence of density, permeability, and slope angle on the granular column collapse, the impact of volume on the run-out characteristics has not been investigated. This study aims to understand how the initial volume affects the run-out using a two-dimensional coupled lattice Boltzman and discrete element (LBM-DEM) method. The coupled LBM-DEM approach allows simulating fluid flow at the pore-scale resolution to understand the grain-scale mechanisms driving the complex continuum-scale response in the granular column collapse. For submerged granular column collapse, the run-out mechanism is heavily influenced by the interaction between the grains and the surrounding fluid. The development of negative pore pressures during shearing and hydrodynamic drag forces inhibit the flow. On the other hand, entrainment of water resulting in hydroplaning enhances the flow. With an increase in volume, the interaction between the grains and the surrounding fluid varies, causing changes in the run-out behavior. For smaller volumes, the forces inhibiting the underwater flow predominates, resulting in shorter run-outs than their dry counterparts. At large volumes, hydroplaning results in larger run-out than the dry cases, despite the inhibiting effects of drag forces and negative pore pressures.

The soil water retention curve (SWRC) is the most commonly used relationship in the study of unsaturated soil. In this paper, the effect of porosity on the SWRC is investigated by numerically modeling unsaturated soil using the Shan-Chen multiphase Lattice Boltzmann Method. The shape of simulated SWRCs are compared against that predicted by the van Genuchten model, demonstrating a good fit except at low degrees of saturation. The simulated SWRCs show an increase in the air-entry value as porosity decreases.

We develop a PyTorch-based Graph Network Simulator (GNS) that learns physics and predicts the flow behavior of particulate and fluid systems. GNS discretizes the domain with nodes representing a collection of material points and the links connecting the nodes representing the local interaction between particles or clusters of particles. The GNS learns the interaction laws through message passing on the graph. GNS has three components: (a) Encoder, which embeds particle information to a latent graph, the edges are learned functions; (b) Processor, which allows data propagation and computes the nodal interactions across steps; and (c) Decoder, which extracts the relevant dynamics (e.g., particle acceleration) from the graph. We introduce physics-inspired simple inductive biases, such as an inertial frame that allows learning algorithms to prioritize one solution (constant gravitational acceleration) over another, reducing learning time. The GNS implementation uses semi-implicit Euler integration to update the next state based on the predicted accelerations. GNS trained on trajectory data is generalizable to predict particle kinematics in complex boundary conditions not seen during training. The trained model accurately predicts within a 5% error of its associated material point method (MPM) simulation. The predictions are 5,000x faster than traditional MPM simulations (2.5 hours for MPM simulations versus 20 s for GNS simulation of granular flow). GNS surrogates are popular for solving optimization, control, critical-region prediction for in situ viz, and inverse-type problems. The GNS code is available under the open-source MIT license at https://github.com/geoelements/gns.

We present a novel method of using explainability techniques to design physics-aware neural networks. We demonstrate our approach by developing a convolutional neural network (CNN) for solving an inverse problem for shallow subsurface imaging. Although CNNs have gained popularity in recent years across many fields, the development of CNNs remains an art, as there are no clear guidelines regarding the selection of hyperparameters that will yield the best network. While optimization algorithms may be used to select hyperparameters automatically, these methods focus on developing networks with high predictive accuracy while disregarding model explainability (descriptive accuracy). However, the field of Explainable Artificial Intelligence (XAI) addresses the absence of model explainability by providing tools that allow developers to evaluate the internal logic of neural networks. In this study, we use the explainability methods Score-CAM and Deep SHAP to select hyperparameters, such as kernel sizes and network depth, to develop a physics-aware CNN for shallow subsurface imaging. We begin with a relatively deep Encoder-Decoder network, which uses surface wave dispersion images as inputs and generates 2D shear wave velocity subsurface images as outputs. Through model explanations, we ultimately find that a shallow CNN using two convolutional layers with an atypical kernel size of 3x1 yields comparable predictive accuracy but with increased descriptive accuracy. We also show that explainability methods can be used to evaluate the network’s complexity and decision-making. We believe this method can be used to develop neural networks with high predictive accuracy while also providing inherent explainability.

Individual particle rotation and displacement were measured in triaxial tests on transparent sand stabilized with geogrid simulants. The Cellpose U-Net model, originally developed to segment biological cells, was trained to segment images of fused quartz particles. The Score-CAM metric from the field of Explainable AI was used to validate the application of Cellpose to segment particles of fused quartz. These segmented particles were characterized in terms of Fourier shape descriptors and tracked across images. The measured particle displacements in the monotonic triaxial tests correlated with displacement fields from Digital Image Correlation (DIC). In contrast to DIC, the new technique also allows for the measurement of individual particle rotation. The particle rotation measurements were found to be repeatable across different specimens. A state boundary line between probable and improbable particle motions could be identified for a given test based on the measured particle displacements and rotations. The size of the zone of probable motions was used to quantify the effectiveness of the stabilizing inclusions. The results of repeated load tests revealed that the honeycomb inclusions used stabilized the specimens by reducing both particle displacements and rotations.

This paper presents a frequency-velocity convolutional neural network (CNN) for rapid, non-invasive 2D shear wave velocity (Vs) imaging of near-surface geo-materials. Operating in the frequency-velocity domain allows for significant flexibility in the linear-array, active-source experimental testing configurations used for generating the CNN input, which are normalized dispersion images. Unlike wavefield images, normalized dispersion images are relatively insensitive to the experimental testing configuration, accommodating various source types, source offsets, numbers of receivers, and receiver spacings. We demonstrate the effectiveness of the frequency-velocity CNN by applying it to a classic near-surface geophysics problem, namely, imaging a two-layer, undulating, soil-over-bedrock interface. This problem was recently investigated in our group by developing a time-distance CNN, which showed great promise but lacked flexibility in utilizing different field-testing configurations. Herein, the new frequency-velocity CNN is shown to have comparable accuracy to the time-distance CNN while providing greater flexibility to handle varied field applications. The frequency-velocity CNN was trained, validated, and tested using 100,000 synthetic near-surface models. The ability of the proposed frequency-velocity CNN to generalize across various acquisition configurations is first tested using synthetic near-surface models with different acquisition configurations from that of the training set, and then applied to experimental field data collected at the Hornsby Bend site in Austin, Texas, USA. When fully developed for a wider range of geological conditions, the proposed CNN may ultimately be used as a rapid, end-to-end alternative for current pseudo-2D surface wave imaging techniques or to develop starting models for full waveform inversion.

In situ visualization techniques are hampered by a lack of foresight: crucial simulation phenomena can be missed due to a poor sampling rate or insufficient detail at critical timesteps. Keeping a human in the loop is impractical, and defining statistical triggers can be difficult. This paper demonstrates the potential for using a machine-learning-based simulation surrogate as an oracle to identify expected critical regions of a large-scale simulation. These critical regions are used to drive the in situ analysis, providing greater data fidelity and analysis resolution with an equivalent I/O budget to a traditional in situ framework. We develop a distributed asynchronous in situ visualization by integrating TACC Galaxy with CB-Geo MPM for material point simulation of granular flows. We employ a PyTorch-based 3D Graph Network Simulator (GNS) trained on granular flow problems as an oracle to predict the dynamics of granular flows. Critical regions of interests are manually tagged in GNS for in situ rendering in MPM.

Shear stress history controls the pore pressure response in liquefiable soils. The excess pore pressure does not increase under cyclic loading when shear stress amplitude is lower than the peak prior amplitude – the shielding effect. Many sophisticated constitutive models fail to capture the shielding effect observed in the cyclic liquefaction experiments. We develop a data-driven machine learning model based on the LSTM neural network to capture the liquefaction response of soils under cyclic loading. The LSTM model is trained on 12 laboratory cyclic simple shear tests on Nevada sand in loose and dense conditions subjected to different cyclic simple shear loading conditions. The LSTM model features include the relative density of soil and the previous stress history to predict the pore water pressure response. The LSTM model successfully replicates the pore pressure response for three cyclic simple test results considering the shielding and density effects.

Non-invasive subsurface imaging using full waveform inversion (FWI) has the potential to fundamen-tally change engineering site characterization by enabling the recovery of high resolution 2D/3D maps ofsubsurface stiffness. Yet, the accuracy of FWI remains quite sensitive to the choice of the initial start-ing model due to the complexity and non-uniqueness of the inverse problem. In response, we presentthe novel application of convolutional neural networks (CNNs) to transform an experimental seismicwavefield acquired using a linear array of surface sensors directly into a robust starting model for 2DFWI. We begin by describing three key steps used for developing the CNN, which include: selectionof a network architecture, development of a suitable training set, and performance of network training.The ability of the trained CNN to predict a suitable starting model for 2D FWI was compared againstother commonly used starting models for a classic near-surface imaging problem; the identification of anundulating, two-layer, soil-bedrock interface. The CNN developed during this study was able to predictcomplex 2D subsurface images of the testing set directly from their seismic wavefields with an averagemean absolute percent error of 6%. When compared to other common approaches, the CNN approachwas able to produce starting models with smaller seismic image and waveform misfits, both before andafter FWI. The ability of the CNN to generalize to subsurface models which were dissimilar to the onesupon which it was trained was assessed using a more complex, three-layered model. While the predictiveability of the CNN was slightly reduced, it was still able to achieve seismic image and waveform misfitscomparable to the other commonly used starting models. This study demonstrates that CNNs have greatpotential as a tool for developing good starting models for FWI, which are critical for producing accurateFWI results.

In nature, submarine slope failures usually carry thousands of cubic-meters of sediments across extremely long distances and cause tsunamis and damages to offshore structures. This paper uses the granular column collapse experiment to investigate the effect of slope angle on the runout behavior of submarine granular landslides for different initial volumes. A two-dimensional coupled lattice Boltzman and discrete element method (LBM-DEM) approach is adopted for numerically modeling the granular column collapse. Columns with four different slope angles and six different volumes are modelled under both dry and submerged conditions. The effects of hydrodynamic interactions, including the generation of excess pore pressures, hydroplaning, and drag forces and formation of turbulent vortices, are used to explain the difference in the runout behavior of the submerged columns compared to the dry columns. The results show that at any given slope angle, there is a threshold volume above which the submerged columns have a larger final runout compared to their dry counterpart, and this threshold volume decreases with slope angle.

Submarine landslides transport thousands of cubic meters of sediment across continental shelves even at slopes as low as 1° and can cause significant casualty and damage to infrastructure. The run-out mechanism in a submarine landslide is affected by factors such as the initial packing density, permeability, slope angle, and initial volume. While past studies have focused on the influence of density, permeability, and slope angle on the granular column collapse, the impact of volume on the run-out characteristics has not been investigated. This study aims to understand how the initial volume affects the run-out using a two-dimensional coupled lattice Boltzman and discrete element (LBM-DEM) method. The coupled LBM-DEM approach allows simulating fluid flow at the pore-scale resolution to understand the grain-scale mechanisms driving the complex continuum-scale response in the granular column collapse. For submerged granular column collapse, the run-out mechanism is heavily influenced by the interaction between the grains and the surrounding fluid. The development of negative pore pressures during shearing and hydrodynamic drag forces inhibit the flow. On the other hand, entrainment of water resulting in hydroplaning enhances the flow. With an increase in volume, the interaction between the grains and the surrounding fluid varies, causing changes in the run-out behavior. For smaller volumes, the forces inhibiting the underwater flow predominates, resulting in shorter run-outs than their dry counterparts. At large volumes, hydroplaning results in larger run-out than the dry cases, despite the inhibiting effects of drag forces and negative pore pressures.

We present a novel method of using explainability techniques to design physics-aware neural networks. We demonstrate our approach by developing a convolutional neural network (CNN) for solving an inverse problem for shallow subsurface imaging. Although CNNs have gained popularity in recent years across many fields, the development of CNNs remains an art, as there are no clear guidelines regarding the selection of hyperparameters that will yield the best network. While optimization algorithms may be used to select hyperparameters automatically, these methods focus on developing networks with high predictive accuracy while disregarding model explainability (descriptive accuracy). However, the field of Explainable Artificial Intelligence (XAI) addresses the absence of model explainability by providing tools that allow developers to evaluate the internal logic of neural networks. In this study, we use the explainability methods Score-CAM and Deep SHAP to select hyperparameters, such as kernel sizes and network depth, to develop a physics-aware CNN for shallow subsurface imaging. We begin with a relatively deep Encoder-Decoder network, which uses surface wave dispersion images as inputs and generates 2D shear wave velocity subsurface images as outputs. Through model explanations, we ultimately find that a shallow CNN using two convolutional layers with an atypical kernel size of 3x1 yields comparable predictive accuracy but with increased descriptive accuracy. We also show that explainability methods can be used to evaluate the network’s complexity and decision-making. We believe this method can be used to develop neural networks with high predictive accuracy while also providing inherent explainability.

Individual particle rotation and displacement were measured in triaxial tests on transparent sand stabilized with geogrid simulants. The Cellpose U-Net model, originally developed to segment biological cells, was trained to segment images of fused quartz particles. The Score-CAM metric from the field of Explainable AI was used to validate the application of Cellpose to segment particles of fused quartz. These segmented particles were characterized in terms of Fourier shape descriptors and tracked across images. The measured particle displacements in the monotonic triaxial tests correlated with displacement fields from Digital Image Correlation (DIC). In contrast to DIC, the new technique also allows for the measurement of individual particle rotation. The particle rotation measurements were found to be repeatable across different specimens. A state boundary line between probable and improbable particle motions could be identified for a given test based on the measured particle displacements and rotations. The size of the zone of probable motions was used to quantify the effectiveness of the stabilizing inclusions. The results of repeated load tests revealed that the honeycomb inclusions used stabilized the specimens by reducing both particle displacements and rotations.

This paper presents a frequency-velocity convolutional neural network (CNN) for rapid, non-invasive 2D shear wave velocity (Vs) imaging of near-surface geo-materials. Operating in the frequency-velocity domain allows for significant flexibility in the linear-array, active-source experimental testing configurations used for generating the CNN input, which are normalized dispersion images. Unlike wavefield images, normalized dispersion images are relatively insensitive to the experimental testing configuration, accommodating various source types, source offsets, numbers of receivers, and receiver spacings. We demonstrate the effectiveness of the frequency-velocity CNN by applying it to a classic near-surface geophysics problem, namely, imaging a two-layer, undulating, soil-over-bedrock interface. This problem was recently investigated in our group by developing a time-distance CNN, which showed great promise but lacked flexibility in utilizing different field-testing configurations. Herein, the new frequency-velocity CNN is shown to have comparable accuracy to the time-distance CNN while providing greater flexibility to handle varied field applications. The frequency-velocity CNN was trained, validated, and tested using 100,000 synthetic near-surface models. The ability of the proposed frequency-velocity CNN to generalize across various acquisition configurations is first tested using synthetic near-surface models with different acquisition configurations from that of the training set, and then applied to experimental field data collected at the Hornsby Bend site in Austin, Texas, USA. When fully developed for a wider range of geological conditions, the proposed CNN may ultimately be used as a rapid, end-to-end alternative for current pseudo-2D surface wave imaging techniques or to develop starting models for full waveform inversion.

Non-invasive subsurface imaging using full waveform inversion (FWI) has the potential to fundamen-tally change engineering site characterization by enabling the recovery of high resolution 2D/3D maps ofsubsurface stiffness. Yet, the accuracy of FWI remains quite sensitive to the choice of the initial start-ing model due to the complexity and non-uniqueness of the inverse problem. In response, we presentthe novel application of convolutional neural networks (CNNs) to transform an experimental seismicwavefield acquired using a linear array of surface sensors directly into a robust starting model for 2DFWI. We begin by describing three key steps used for developing the CNN, which include: selectionof a network architecture, development of a suitable training set, and performance of network training.The ability of the trained CNN to predict a suitable starting model for 2D FWI was compared againstother commonly used starting models for a classic near-surface imaging problem; the identification of anundulating, two-layer, soil-bedrock interface. The CNN developed during this study was able to predictcomplex 2D subsurface images of the testing set directly from their seismic wavefields with an averagemean absolute percent error of 6%. When compared to other common approaches, the CNN approachwas able to produce starting models with smaller seismic image and waveform misfits, both before andafter FWI. The ability of the CNN to generalize to subsurface models which were dissimilar to the onesupon which it was trained was assessed using a more complex, three-layered model. While the predictiveability of the CNN was slightly reduced, it was still able to achieve seismic image and waveform misfitscomparable to the other commonly used starting models. This study demonstrates that CNNs have greatpotential as a tool for developing good starting models for FWI, which are critical for producing accurateFWI results.

The performance evaluation of a potentially unstable slope involves two key components: the initiation of the slope failure and the post-failure runout. The Finite Element Method (FEM) excels at modeling the initiation of instability but quickly loses accuracy in modeling large-deformation problems due to mesh distortion. Hence, the FEM is unable to accurately model post-failure slope runout. Hybrid Eulerian-Lagrangian methods, such as the Material Point Method (MPM), offer a promising alternative for solving large-deformation problems, because particles can move freely across a background mesh, allowing for large deformation without computational issues. However, the use of moving material points in MPM for integration rather than the fixed Gauss points of the FEM reduces the accuracy of MPM in predicting stress distribution and thus failure initiation. We have created a hybrid method by initiating a failure simulation in FEM and subsequently transferring the coordinates, velocities, and stresses to MPM particles to model the runout behavior, combining the strength of both methods. We demonstrate the capability of the hybrid approach by simulating the collapse of a frictional granular column, comparing it to an empirical solution, and evaluating a suitable time to transfer from FEM to MPM by trialing multiple iterations with transfers at different stages of the collapse.

NSF OAC Awards Cognitasium project to develop new data-driven discovery workflows in natural hazards

Civil Engineering analysis technique of finite elements is used for the first time to answer an evolutionary question

Compared to most other primates, humans are characterized by a tight fit between the maternal birth canal and the fetal head, leading to a relatively high risk of neonatal and maternal mortality and morbidities. Obstetric selection is thought to favor a spacious birth canal, whereas the source for opposing selection is frequently assumed to relate to bipedal locomotion. An alternative, yet under-investigated, hypothesis is that a more expansive birth canal suspends the soft tissue of the pelvic floor across a larger area, which is disadvantageous for continence and support of the weight of the inner organs and fetus. To test this “pelvic floor hypothesis” we generated a finite element model of the human female pelvic floor and varied its radial size and thickness while keeping all else constant. This allowed us to study the effect of pelvic geometry on pelvic floor deflection (i.e., the amount of bending from the original position) and tissue stresses and stretches. Deflection grew disproportionately fast with increasing radial size, and stresses and stretches also increased. By contrast, an increase in thickness increased pelvic floor stiffness - i.e. the resistance to deformation - which reduced deflection but was unable to fully compensate for the effect of increasing radial size. Moreover, larger thicknesses increase the intra-abdominal pressure necessary for childbirth. Our results support the pelvic floor hypothesis and evince functional trade-offs affecting not only the size of the birth canal but also the thickness and stiffness of the pelvic floor.

The soil water retention curve (SWRC) is the most commonly used relationship in the study of unsaturated soil. In this paper, the effect of porosity on the SWRC is investigated by numerically modeling unsaturated soil using the Shan-Chen multiphase Lattice Boltzmann Method. The shape of simulated SWRCs are compared against that predicted by the van Genuchten model, demonstrating a good fit except at low degrees of saturation. The simulated SWRCs show an increase in the air-entry value as porosity decreases.

Geoelements group wins TACC pathways proposal to simulate the Oso landslide.

Individual particle rotation and displacement were measured in triaxial tests on transparent sand stabilized with geogrid simulants. The Cellpose U-Net model, originally developed to segment biological cells, was trained to segment images of fused quartz particles. The Score-CAM metric from the field of Explainable AI was used to validate the application of Cellpose to segment particles of fused quartz. These segmented particles were characterized in terms of Fourier shape descriptors and tracked across images. The measured particle displacements in the monotonic triaxial tests correlated with displacement fields from Digital Image Correlation (DIC). In contrast to DIC, the new technique also allows for the measurement of individual particle rotation. The particle rotation measurements were found to be repeatable across different specimens. A state boundary line between probable and improbable particle motions could be identified for a given test based on the measured particle displacements and rotations. The size of the zone of probable motions was used to quantify the effectiveness of the stabilizing inclusions. The results of repeated load tests revealed that the honeycomb inclusions used stabilized the specimens by reducing both particle displacements and rotations.

The water retention curve (WRC) defines the relationship between matric suction and saturation and is a key function for determining the hydro-mechanical behavior of unsaturated soils. We investigate possible microscopic origins of the water retention behavior of granular soils using both Computed Tomography (CT) experiment and multiphase lattice Boltzmann Method (LBM). We conduct a CT experiment on Hamburg sand to obtain its WRC and then run LBM simulations based on the CT grain skeleton. The multiphase LBM simulations capture the hysteresis and pore-scale behaviors of WRC observed in the CT experiment. Using LBM, we observe that the spatial distribution and morphology of gas clusters varies between drainage and imbibition paths and is the underlying source of the hysteresis. During drainage, gas clusters congregate at the grain surface; the local suction increases when gas clusters enter through small pore openings and decreases when gas clusters enter through large pore openings. Whereas, during imbibition, gas clusters disperse in the liquid; the local suction decreases uniformly. Large pores empty first during drainage and small pores fill first during imbibition. The pore-based WRC shows that an increase in pore size causes a decrease in suction during drainage and imbibition, and an increase in hysteresis.

The soil-water characteristic curve (SWCC) is the most fundamental relationship in unsaturated soil mechanics, relating the amount of water in the soil to the corresponding matric suction. From experimental evidence, it is known that SWCC exhibits hysteresis (i.e. wetting/drying path dependence). Various factors have been proposed as contributors to SWCC hysteresis, including air entrapment, contact angle hysteresis, ink-bottle effect, and change of soil fabric due to swelling and shrinkage, however, the significance of their contribution is debated. From our pore-scale numerical simulations, using the multiphase lattice Boltzmann method, we see that even when controlling for all these factors SWCC hysteresis still occurs, indicating that there is some underlying source that is not accounted for in these factors. We find this underlying source by comparing the liquid/gas phase distributions for simulated wetting and drying experiments of 2D and 3D granular packings. We see that during wetting (i.e. pore filling) many liquid bridges expand simultaneously and join together to fill the pores from the smallest to the largest, allowing menisci with larger radii of curvature (lower matric suction). Whereas, during drying (i.e. pore emptying), only the limited existing gas clusters can expand, which become constrained by the size of the pore openings surrounding them and result in menisci with smaller radii of curvature (higher matric suction).

The soil water retention curve (SWRC) is the most commonly used relationship in the study of unsaturated soil. In this paper, the effect of porosity on the SWRC is investigated by numerically modeling unsaturated soil using the Shan-Chen multiphase Lattice Boltzmann Method. The shape of simulated SWRCs are compared against that predicted by the van Genuchten model, demonstrating a good fit except at low degrees of saturation. The simulated SWRCs show an increase in the air-entry value as porosity decreases.

NSF OAC Awards Cognitasium project to develop new data-driven discovery workflows in natural hazards

Civil Engineering analysis technique of finite elements is used for the first time to answer an evolutionary question

Compared to most other primates, humans are characterized by a tight fit between the maternal birth canal and the fetal head, leading to a relatively high risk of neonatal and maternal mortality and morbidities. Obstetric selection is thought to favor a spacious birth canal, whereas the source for opposing selection is frequently assumed to relate to bipedal locomotion. An alternative, yet under-investigated, hypothesis is that a more expansive birth canal suspends the soft tissue of the pelvic floor across a larger area, which is disadvantageous for continence and support of the weight of the inner organs and fetus. To test this “pelvic floor hypothesis” we generated a finite element model of the human female pelvic floor and varied its radial size and thickness while keeping all else constant. This allowed us to study the effect of pelvic geometry on pelvic floor deflection (i.e., the amount of bending from the original position) and tissue stresses and stretches. Deflection grew disproportionately fast with increasing radial size, and stresses and stretches also increased. By contrast, an increase in thickness increased pelvic floor stiffness - i.e. the resistance to deformation - which reduced deflection but was unable to fully compensate for the effect of increasing radial size. Moreover, larger thicknesses increase the intra-abdominal pressure necessary for childbirth. Our results support the pelvic floor hypothesis and evince functional trade-offs affecting not only the size of the birth canal but also the thickness and stiffness of the pelvic floor.

We present a novel method of using explainability techniques to design physics-aware neural networks. We demonstrate our approach by developing a convolutional neural network (CNN) for solving an inverse problem for shallow subsurface imaging. Although CNNs have gained popularity in recent years across many fields, the development of CNNs remains an art, as there are no clear guidelines regarding the selection of hyperparameters that will yield the best network. While optimization algorithms may be used to select hyperparameters automatically, these methods focus on developing networks with high predictive accuracy while disregarding model explainability (descriptive accuracy). However, the field of Explainable Artificial Intelligence (XAI) addresses the absence of model explainability by providing tools that allow developers to evaluate the internal logic of neural networks. In this study, we use the explainability methods Score-CAM and Deep SHAP to select hyperparameters, such as kernel sizes and network depth, to develop a physics-aware CNN for shallow subsurface imaging. We begin with a relatively deep Encoder-Decoder network, which uses surface wave dispersion images as inputs and generates 2D shear wave velocity subsurface images as outputs. Through model explanations, we ultimately find that a shallow CNN using two convolutional layers with an atypical kernel size of 3x1 yields comparable predictive accuracy but with increased descriptive accuracy. We also show that explainability methods can be used to evaluate the network’s complexity and decision-making. We believe this method can be used to develop neural networks with high predictive accuracy while also providing inherent explainability.

This paper presents a frequency-velocity convolutional neural network (CNN) for rapid, non-invasive 2D shear wave velocity (Vs) imaging of near-surface geo-materials. Operating in the frequency-velocity domain allows for significant flexibility in the linear-array, active-source experimental testing configurations used for generating the CNN input, which are normalized dispersion images. Unlike wavefield images, normalized dispersion images are relatively insensitive to the experimental testing configuration, accommodating various source types, source offsets, numbers of receivers, and receiver spacings. We demonstrate the effectiveness of the frequency-velocity CNN by applying it to a classic near-surface geophysics problem, namely, imaging a two-layer, undulating, soil-over-bedrock interface. This problem was recently investigated in our group by developing a time-distance CNN, which showed great promise but lacked flexibility in utilizing different field-testing configurations. Herein, the new frequency-velocity CNN is shown to have comparable accuracy to the time-distance CNN while providing greater flexibility to handle varied field applications. The frequency-velocity CNN was trained, validated, and tested using 100,000 synthetic near-surface models. The ability of the proposed frequency-velocity CNN to generalize across various acquisition configurations is first tested using synthetic near-surface models with different acquisition configurations from that of the training set, and then applied to experimental field data collected at the Hornsby Bend site in Austin, Texas, USA. When fully developed for a wider range of geological conditions, the proposed CNN may ultimately be used as a rapid, end-to-end alternative for current pseudo-2D surface wave imaging techniques or to develop starting models for full waveform inversion.

Non-invasive subsurface imaging using full waveform inversion (FWI) has the potential to fundamen-tally change engineering site characterization by enabling the recovery of high resolution 2D/3D maps ofsubsurface stiffness. Yet, the accuracy of FWI remains quite sensitive to the choice of the initial start-ing model due to the complexity and non-uniqueness of the inverse problem. In response, we presentthe novel application of convolutional neural networks (CNNs) to transform an experimental seismicwavefield acquired using a linear array of surface sensors directly into a robust starting model for 2DFWI. We begin by describing three key steps used for developing the CNN, which include: selectionof a network architecture, development of a suitable training set, and performance of network training.The ability of the trained CNN to predict a suitable starting model for 2D FWI was compared againstother commonly used starting models for a classic near-surface imaging problem; the identification of anundulating, two-layer, soil-bedrock interface. The CNN developed during this study was able to predictcomplex 2D subsurface images of the testing set directly from their seismic wavefields with an averagemean absolute percent error of 6%. When compared to other common approaches, the CNN approachwas able to produce starting models with smaller seismic image and waveform misfits, both before andafter FWI. The ability of the CNN to generalize to subsurface models which were dissimilar to the onesupon which it was trained was assessed using a more complex, three-layered model. While the predictiveability of the CNN was slightly reduced, it was still able to achieve seismic image and waveform misfitscomparable to the other commonly used starting models. This study demonstrates that CNNs have greatpotential as a tool for developing good starting models for FWI, which are critical for producing accurateFWI results.

Inductive heating with high-voltage cables reduces the risk of hydrate formation by raising the temperature of the production fluid in pipelines. Heating the pipeline results in losing a certain fraction of the heat to the surrounding soil through conduction or convection-dominated flow through the soil. However, the amount of heat lost in conduction versus convection and the transition from conduction to convection-dominated heat loss remains unknown. Soil permeability, temperature gradient between cable and mudline, and burial depth influence the mode of heat transfer and the amount of heat lost. We study the dominant mode of heat transfer in pipelines with inductive heating using 2D Finite Difference analysis under different soil and environmental conditions. Low permeability soils primarily exhibit conductive heat transfer, thus losing minimum heat to the surrounding soil. In contrast, convective flow drives a significant fraction of the heat away from the pipeline and towards the ground surface for highly permeable soils, barely heating the fluid in the pipe. We identify a critical Rayleigh-Darcy number of 1 as the controlling value separating conduction and convection-dominated heat transfer. An increase in burial depth deteriorates the heating efficiency in convection-dominated high permeability soils, while it remains unaffected in conduction-dominated low permeability soils.

Facebook awards Dr Chadi El Mohtar and Dr Krishna Kumar an industry researh project to improve the heat dissipation in duct banks.

When power electronics are deployed under the road surface as part of a wireless system it is important to know that their packaging provides adequate heat extraction as well as the required environmental protection – often conflicting requirements. Presently very little can be found in wireless charging standards and literature on the topic of thermal modelling for in-ground components. Yet, this is a topic of great practical significance especially for in-road systems. Traditional cooling methods are not readily applicable underground. This paper uses finite element thermal modelling to investigate the cooling of a representative medium-power in-road wireless system, housed in a sealed ground assembly (GA) chamber and installed to UK requirements (HAUC). The paper quantitatively compares design options and provides practical recommendations for in-road installation thermal management.

The water retention curve (WRC) defines the relationship between matric suction and saturation and is a key function for determining the hydro-mechanical behavior of unsaturated soils. We investigate possible microscopic origins of the water retention behavior of granular soils using both Computed Tomography (CT) experiment and multiphase lattice Boltzmann Method (LBM). We conduct a CT experiment on Hamburg sand to obtain its WRC and then run LBM simulations based on the CT grain skeleton. The multiphase LBM simulations capture the hysteresis and pore-scale behaviors of WRC observed in the CT experiment. Using LBM, we observe that the spatial distribution and morphology of gas clusters varies between drainage and imbibition paths and is the underlying source of the hysteresis. During drainage, gas clusters congregate at the grain surface; the local suction increases when gas clusters enter through small pore openings and decreases when gas clusters enter through large pore openings. Whereas, during imbibition, gas clusters disperse in the liquid; the local suction decreases uniformly. Large pores empty first during drainage and small pores fill first during imbibition. The pore-based WRC shows that an increase in pore size causes a decrease in suction during drainage and imbibition, and an increase in hysteresis.

The soil-water characteristic curve (SWCC) is the most fundamental relationship in unsaturated soil mechanics, relating the amount of water in the soil to the corresponding matric suction. From experimental evidence, it is known that SWCC exhibits hysteresis (i.e. wetting/drying path dependence). Various factors have been proposed as contributors to SWCC hysteresis, including air entrapment, contact angle hysteresis, ink-bottle effect, and change of soil fabric due to swelling and shrinkage, however, the significance of their contribution is debated. From our pore-scale numerical simulations, using the multiphase lattice Boltzmann method, we see that even when controlling for all these factors SWCC hysteresis still occurs, indicating that there is some underlying source that is not accounted for in these factors. We find this underlying source by comparing the liquid/gas phase distributions for simulated wetting and drying experiments of 2D and 3D granular packings. We see that during wetting (i.e. pore filling) many liquid bridges expand simultaneously and join together to fill the pores from the smallest to the largest, allowing menisci with larger radii of curvature (lower matric suction). Whereas, during drying (i.e. pore emptying), only the limited existing gas clusters can expand, which become constrained by the size of the pore openings surrounding them and result in menisci with smaller radii of curvature (higher matric suction).

The soil water retention curve (SWRC) is the most commonly used relationship in the study of unsaturated soil. In this paper, the effect of porosity on the SWRC is investigated by numerically modeling unsaturated soil using the Shan-Chen multiphase Lattice Boltzmann Method. The shape of simulated SWRCs are compared against that predicted by the van Genuchten model, demonstrating a good fit except at low degrees of saturation. The simulated SWRCs show an increase in the air-entry value as porosity decreases.

We develop a PyTorch-based Graph Network Simulator (GNS) that learns physics and predicts the flow behavior of particulate and fluid systems. GNS discretizes the domain with nodes representing a collection of material points and the links connecting the nodes representing the local interaction between particles or clusters of particles. The GNS learns the interaction laws through message passing on the graph. GNS has three components: (a) Encoder, which embeds particle information to a latent graph, the edges are learned functions; (b) Processor, which allows data propagation and computes the nodal interactions across steps; and (c) Decoder, which extracts the relevant dynamics (e.g., particle acceleration) from the graph. We introduce physics-inspired simple inductive biases, such as an inertial frame that allows learning algorithms to prioritize one solution (constant gravitational acceleration) over another, reducing learning time. The GNS implementation uses semi-implicit Euler integration to update the next state based on the predicted accelerations. GNS trained on trajectory data is generalizable to predict particle kinematics in complex boundary conditions not seen during training. The trained model accurately predicts within a 5% error of its associated material point method (MPM) simulation. The predictions are 5,000x faster than traditional MPM simulations (2.5 hours for MPM simulations versus 20 s for GNS simulation of granular flow). GNS surrogates are popular for solving optimization, control, critical-region prediction for in situ viz, and inverse-type problems. The GNS code is available under the open-source MIT license at https://github.com/geoelements/gns.

In situ visualization techniques are hampered by a lack of foresight: crucial simulation phenomena can be missed due to a poor sampling rate or insufficient detail at critical timesteps. Keeping a human in the loop is impractical, and defining statistical triggers can be difficult. This paper demonstrates the potential for using a machine-learning-based simulation surrogate as an oracle to identify expected critical regions of a large-scale simulation. These critical regions are used to drive the in situ analysis, providing greater data fidelity and analysis resolution with an equivalent I/O budget to a traditional in situ framework. We develop a distributed asynchronous in situ visualization by integrating TACC Galaxy with CB-Geo MPM for material point simulation of granular flows. We employ a PyTorch-based 3D Graph Network Simulator (GNS) trained on granular flow problems as an oracle to predict the dynamics of granular flows. Critical regions of interests are manually tagged in GNS for in situ rendering in MPM.

In nature, submarine slope failures usually carry thousands of cubic-meters of sediments across extremely long distances and cause tsunamis and damages to offshore structures. This paper uses the granular column collapse experiment to investigate the effect of slope angle on the runout behavior of submarine granular landslides for different initial volumes. A two-dimensional coupled lattice Boltzman and discrete element method (LBM-DEM) approach is adopted for numerically modeling the granular column collapse. Columns with four different slope angles and six different volumes are modelled under both dry and submerged conditions. The effects of hydrodynamic interactions, including the generation of excess pore pressures, hydroplaning, and drag forces and formation of turbulent vortices, are used to explain the difference in the runout behavior of the submerged columns compared to the dry columns. The results show that at any given slope angle, there is a threshold volume above which the submerged columns have a larger final runout compared to their dry counterpart, and this threshold volume decreases with slope angle.

Submarine landslides transport thousands of cubic meters of sediment across continental shelves even at slopes as low as 1° and can cause significant casualty and damage to infrastructure. The run-out mechanism in a submarine landslide is affected by factors such as the initial packing density, permeability, slope angle, and initial volume. While past studies have focused on the influence of density, permeability, and slope angle on the granular column collapse, the impact of volume on the run-out characteristics has not been investigated. This study aims to understand how the initial volume affects the run-out using a two-dimensional coupled lattice Boltzman and discrete element (LBM-DEM) method. The coupled LBM-DEM approach allows simulating fluid flow at the pore-scale resolution to understand the grain-scale mechanisms driving the complex continuum-scale response in the granular column collapse. For submerged granular column collapse, the run-out mechanism is heavily influenced by the interaction between the grains and the surrounding fluid. The development of negative pore pressures during shearing and hydrodynamic drag forces inhibit the flow. On the other hand, entrainment of water resulting in hydroplaning enhances the flow. With an increase in volume, the interaction between the grains and the surrounding fluid varies, causing changes in the run-out behavior. For smaller volumes, the forces inhibiting the underwater flow predominates, resulting in shorter run-outs than their dry counterparts. At large volumes, hydroplaning results in larger run-out than the dry cases, despite the inhibiting effects of drag forces and negative pore pressures.

In situ visualization techniques are hampered by a lack of foresight: crucial simulation phenomena can be missed due to a poor sampling rate or insufficient detail at critical timesteps. Keeping a human in the loop is impractical, and defining statistical triggers can be difficult. This paper demonstrates the potential for using a machine-learning-based simulation surrogate as an oracle to identify expected critical regions of a large-scale simulation. These critical regions are used to drive the in situ analysis, providing greater data fidelity and analysis resolution with an equivalent I/O budget to a traditional in situ framework. We develop a distributed asynchronous in situ visualization by integrating TACC Galaxy with CB-Geo MPM for material point simulation of granular flows. We employ a PyTorch-based 3D Graph Network Simulator (GNS) trained on granular flow problems as an oracle to predict the dynamics of granular flows. Critical regions of interests are manually tagged in GNS for in situ rendering in MPM.

Visualizing regional-scale landslides is the key to conveying the threat of natural hazards to stakeholders and policymakers. Traditional visualization techniques are restricted to post-processing a limited subset of simulation data and are not scalable to rendering exascale models with billions of particles. In-situ visualization is a technique of rendering simulation data in real-time, i.e., rendering visuals in tandem while the simulation is running. In this study, we develop a scalable N:M interface architecture to visualize regional-scale landslides. We demonstrate the scalability of the architecture by simulating the long runout of the 2014 Oso landslide using the Material Point Method coupled with the Galaxy ray tracing engine rendering 4.2 million material points as spheres. In-situ visualization has an amortized runtime increase of 2\% compared to non-visualized simulations. The developed approach can achieve in-situ visualization of regional-scale landslides with billions of particles with minimal impact on the simulation process.

Abstract:

In situ visualization techniques are hampered by a lack of foresight: crucial simulation phenomena can be missed due to a poor sampling rate or insufficient detail at critical timesteps. Keeping a human in the loop is impractical, and defining statistical triggers can be difficult. This paper demonstrates the potential for using a machine-learning-based simulation surrogate as an oracle to identify expected critical regions of a large-scale simulation. These critical regions are used to drive the in situ analysis, providing greater data fidelity and analysis resolution with an equivalent I/O budget to a traditional in situ framework. We develop a distributed asynchronous in situ visualization by integrating TACC Galaxy with CB-Geo MPM for material point simulation of granular flows. We employ a PyTorch-based 3D Graph Network Simulator (GNS) trained on granular flow problems as an oracle to predict the dynamics of granular flows. Critical regions of interests are manually tagged in GNS for in situ rendering in MPM.

Visualizing regional-scale landslides is the key to conveying the threat of natural hazards to stakeholders and policymakers. Traditional visualization techniques are restricted to post-processing a limited subset of simulation data and are not scalable to rendering exascale models with billions of particles. In-situ visualization is a technique of rendering simulation data in real-time, i.e., rendering visuals in tandem while the simulation is running. In this study, we develop a scalable N:M interface architecture to visualize regional-scale landslides. We demonstrate the scalability of the architecture by simulating the long runout of the 2014 Oso landslide using the Material Point Method coupled with the Galaxy ray tracing engine rendering 4.2 million material points as spheres. In-situ visualization has an amortized runtime increase of 2\% compared to non-visualized simulations. The developed approach can achieve in-situ visualization of regional-scale landslides with billions of particles with minimal impact on the simulation process.

A new thixotropic model is developed integrating the Papanastasiou-Bingham model with thixotropy equations to simulate the flow behaviour of Tremie Concrete in the Material Point Method framework. The effect of thixotropy on the rheological behaviour of fresh concrete is investigated by comparing field measurements with numerical simulations. The comparison yields new insights into a critical and often overlooked behaviour of concrete. A parametric study is performed to understand the effect of model parameters and rest-time on the shear stress response of fresh concrete. The Material Point Method with the Papanastasiou-Bingham model reproduces slump-flow measurements observed in the field. The novel model revealed a decline in concrete workability during the Slump-flow test after a period of rest due to thixotropy, which the physical version of the test fails to capture. This reduction in workability significantly affects the flow behaviour and the effective use of fresh concrete in construction operation.

Individual particle rotation and displacement were measured in triaxial tests on transparent sand stabilized with geogrid simulants. The Cellpose U-Net model, originally developed to segment biological cells, was trained to segment images of fused quartz particles. The Score-CAM metric from the field of Explainable AI was used to validate the application of Cellpose to segment particles of fused quartz. These segmented particles were characterized in terms of Fourier shape descriptors and tracked across images. The measured particle displacements in the monotonic triaxial tests correlated with displacement fields from Digital Image Correlation (DIC). In contrast to DIC, the new technique also allows for the measurement of individual particle rotation. The particle rotation measurements were found to be repeatable across different specimens. A state boundary line between probable and improbable particle motions could be identified for a given test based on the measured particle displacements and rotations. The size of the zone of probable motions was used to quantify the effectiveness of the stabilizing inclusions. The results of repeated load tests revealed that the honeycomb inclusions used stabilized the specimens by reducing both particle displacements and rotations.

Abstract:

Shear stress history controls the pore pressure response in liquefiable soils. The excess pore pressure does not increase under cyclic loading when shear stress amplitude is lower than the peak prior amplitude – the shielding effect. Many sophisticated constitutive models fail to capture the shielding effect observed in the cyclic liquefaction experiments. We develop a data-driven machine learning model based on the LSTM neural network to capture the liquefaction response of soils under cyclic loading. The LSTM model is trained on 12 laboratory cyclic simple shear tests on Nevada sand in loose and dense conditions subjected to different cyclic simple shear loading conditions. The LSTM model features include the relative density of soil and the previous stress history to predict the pore water pressure response. The LSTM model successfully replicates the pore pressure response for three cyclic simple test results considering the shielding and density effects.

Abstract:

The soil-water characteristic curve (SWCC) is the most fundamental relationship in unsaturated soil mechanics, relating the amount of water in the soil to the corresponding matric suction. From experimental evidence, it is known that SWCC exhibits hysteresis (i.e. wetting/drying path dependence). Various factors have been proposed as contributors to SWCC hysteresis, including air entrapment, contact angle hysteresis, ink-bottle effect, and change of soil fabric due to swelling and shrinkage, however, the significance of their contribution is debated. From our pore-scale numerical simulations, using the multiphase lattice Boltzmann method, we see that even when controlling for all these factors SWCC hysteresis still occurs, indicating that there is some underlying source that is not accounted for in these factors. We find this underlying source by comparing the liquid/gas phase distributions for simulated wetting and drying experiments of 2D and 3D granular packings. We see that during wetting (i.e. pore filling) many liquid bridges expand simultaneously and join together to fill the pores from the smallest to the largest, allowing menisci with larger radii of curvature (lower matric suction). Whereas, during drying (i.e. pore emptying), only the limited existing gas clusters can expand, which become constrained by the size of the pore openings surrounding them and result in menisci with smaller radii of curvature (higher matric suction).

When power electronics are deployed under the road surface as part of a wireless system it is important to know that their packaging provides adequate heat extraction as well as the required environmental protection – often conflicting requirements. Presently very little can be found in wireless charging standards and literature on the topic of thermal modelling for in-ground components. Yet, this is a topic of great practical significance especially for in-road systems. Traditional cooling methods are not readily applicable underground. This paper uses finite element thermal modelling to investigate the cooling of a representative medium-power in-road wireless system, housed in a sealed ground assembly (GA) chamber and installed to UK requirements (HAUC). The paper quantitatively compares design options and provides practical recommendations for in-road installation thermal management.