Triply periodic minimal surface (TPMS) metamaterials characterized by mathematically-controlled topologies exhibit better mechanical properties compared to uniform structures. The unit cell topology of such metamaterials can be further optimized to improve a desired mechanical property for a specific application. However, such inverse design involves multiple costly 3D finite element analyses in topology optimization and hence has not been attempted. Data-driven models have recently gained popularity as surrogate models in the geometrical design of metamaterials. Gyroid-like unit cells are designed using a novel voxel algorithm, a homogenization-based topology optimization, and a Heaviside filter to attain optimized densities of 0-1 configuration. Few optimization data are used as input-output for supervised learning of the topology optimization process from a 3D CNN model. These models could then be used to instantaneously predict the optimized unit cell geometry for any topology parameters, thus alleviating the need to run any topology optimization for future design. The high accuracy of the model was demonstrated by a low mean square error metric and a high dice coefficient metric. This accelerated design of 3D metamaterials opens the possibility of designing any computationally costly problems involving complex geometry of metamaterials with multi-objective properties or multi-scale applications.

Deep Operator Network (DeepONet), a recently introduced deep learning operator network, approximates linear and nonlinear solution operators by taking parametric functions (infinite-dimensional objects) as inputs and mapping them to solution functions in contrast to classical neural networks (NNs) that need re-training for every new set of parametric inputs. In this work, we have extended the classical formulation of DeepONets by introducing recurrent neural networks (RNNs) in its branch in so-called sequential DeepONets (S-DeepONets) thus allowing accurate solution predictions in the entire domain for parametric and time-dependent loading histories. We have demonstrated this novel formulation’s generality and exceptional accuracy with thermal and mechanical random loading histories applied to highly nonlinear thermal solidification and plastic deformation use cases. We show that once S-DeepONet is properly trained, it can accurately predict the final solutions in the entire domain and is several orders of magnitude more computationally efficient than the finite element method for arbitrary loading histories without additional training.

Abstract The paper explores the possibility of using the novel Deep Operator Networks (DeepONet) for forward analysis of numerically intensive and challenging multiphysics designs and optimizations of advanced materials and processes. As an important step towards that goal, DeepONet networks were devised and trained on GPUs to solve the Poisson equation (heat-conduction equation) with the spatially variable heat source and highly nonlinear stress distributions under plastic deformation with variable loads and material properties. Since DeepONet can learn the parametric solution of various phenomena and processes in science and engineering, it was found that a properly trained DeepONet can instantly and accurately inference thermal and mechanical solutions for new parametric inputs without re-training and transfer learning and several orders of magnitude faster than classical numerical methods.

In hydrology, projected climate change impact assessment studies typically rely on ensembles of downscaled climate model outputs. Due to large modeling uncertainties, the ensembles are often averaged to provide a basis for studying the effects of climate change. A key issue when analyzing averages of a climate model ensemble is whether to weight all models in the ensemble equally, often referred to as the equal‐weights or unweighted approach, or to use a weighted approach, where, in general, each model would have a different weight. Many studies have advocated for the latter, based on the assumption that models that are better at simulating the past, that is, the models with higher hindcast accuracy, will give more accurate forecasts for the future and thus should receive higher weights. To examine this issue, observed and modeled daily precipitation frequency (PF) estimates for three urban areas in the United States, namely Boston, Massachusetts; Houston, Texas; and Chicago, Illinois, were analyzed. The comparison used the raw output of 24 Coupled Model Intercomparison Project Phase 5 (CMIP5) models. The PFs from these models were compared with the observed PFs for a specific historical training period to determine model weights for each area. The unweighted and weighted averaged model PFs from a more recent testing period were then compared with their corresponding observed PFs to determine if weights improved the estimates. These comparisons indeed showed that the weighted averages were closer to the observed values than the unweighted averages in nearly all cases. The study also demonstrated how weights can help reduce model spread in future climate projections by comparing the unweighted and weighted ensemble standard deviations in these projections. In all studied scenarios, the weights actually reduced the standard deviations compared to the equal‐weights approach. Finally, an analysis of the results' sensitivity to the areal reduction factor used to allow comparisons between point station measurements and grid‐box averages is provided.

A graph convolutional network (GCN) is employed in the deep energy method (DEM) model to solve the momentum balance equation in three‐dimensional space for the deformation of linear elastic and hyperelastic materials due to its ability to handle irregular domains over the traditional DEM method based on a multilayer perceptron (MLP) network. The method's accuracy and solution time are compared to the DEM model based on a MLP network. We demonstrate that the GCN‐based model delivers similar accuracy while having a shorter run time through numerical examples. Two different spatial gradient computation techniques, one based on automatic differentiation (AD) and the other based on shape function (SF) gradients, are also accessed. We provide a simple example to demonstrate the strain localization instability associated with the AD‐based gradient computation and show that the instability exists in more general cases by four numerical examples. The SF‐based gradient computation is shown to be more robust and delivers an accurate solution even at severe deformations. Therefore, the combination of the GCN‐based DEM model and SF‐based gradient computation is potentially a promising candidate for solving problems involving severe material and geometric nonlinearities.

Using recent advancements in high-performance computing data assimilation to combine satellite InSAR data with numerical models, the prolonged unrest of the Sierra Negra volcano in the Galápagos was tracked to provide a fortuitous, but successful, forecast 5 months in advance of the 26 June 2018 eruption. Subsequent numerical simulations reveal that the evolution of the stress state in the host rock surrounding the Sierra Negra magma system likely controlled eruption timing. While changes in magma reservoir pressure remained modest (<15 MPa), modeled widespread Mohr-Coulomb failure is coincident with the timing of the 26 June 2018 moment magnitude 5.4 earthquake and subsequent eruption. Coulomb stress transfer models suggest that the faulting event triggered the 2018 eruption by encouraging tensile failure along the northern portion of the caldera. These findings provide a critical framework for understanding Sierra Negra’s eruption cycles and evaluating the potential and timing of future eruptions.

Physics‐informed neural networks have gained growing interest. Specifically, they are used to solve partial differential equations governing several physical phenomena. However, physics‐informed neural network models suffer from several issues and can fail to provide accurate solutions in many scenarios. We discuss a few of these challenges and the techniques, such as the use of Fourier transform, that can be used to resolve these issues. This paper proposes and develops a physics‐informed neural network model that combines the residuals of the strong form and the potential energy, yielding many loss terms contributing to the definition of the loss function to be minimized. Hence, we propose using the coefficient of variation weighting scheme to dynamically and adaptively assign the weight for each loss term in the loss function. The developed PINN model is standalone and meshfree. In other words, it can accurately capture the mechanical response without requiring any labeled data. Although the framework can be used for many solid mechanics problems, we focus on three‐dimensional (3D) hyperelasticity, where we consider two hyperelastic models. Once the model is trained, the response can be obtained almost instantly at any point in the physical domain, given its spatial coordinates. We demonstrate the framework's performance by solving different problems with various boundary conditions.