|
An immersed phase field fracture model for microporomechanics with Darcy–Stokes flow Physics of Fluids 33, 016603 (2021); https://doi.org/10.1063/5.0035602 Hyoung Suk Suh (서형석) and WaiChing Sun (孫維正) ABSTRACT This paper presents an immersed phase field model designed to predict the fracture-induced flow due to brittle fracture in vuggy porous media. Due to the multiscale nature of pores in the vuggy porous material, crack growth may connect previously isolated pores, which leads to flow conduits. This mechanism has important implications for many applications such as disposal of carbon dioxide and radioactive materials and hydraulic fracture and mining. To understand the detailed microporomechanics that causes the fracture-induced flow, we introduce a new phase field fracture framework where the phase field is not only used as an indicator function for damage of the solid skeleton but also used as an indicator of the pore space. By coupling the Stokes equation that governs the fluid transport in the voids, cavities, and cracks and Darcy’s flow in the deformable porous media, our proposed model enables us to capture the fluid–solid interaction of the pore fluid and solid constituents during crack growth. Numerical experiments are conducted to analyze how the presence of cavities affects the accuracy of predictions based on the homogenized effective medium during crack growth. 
0 Comments
For some polycrystalline materials such as austenitic stainless steel, magnesium, TATB, and HMX, twinning is a crucial deformation mechanism when the dislocation slip alone is not enough to accommodate the applied strain. To predict this coupling effect between crystal plasticity and deformation twinning, we introduce a mathematical model and the corresponding monolithic and operator splitting solver that couples the crystal plasticity material model with a phase field twining model such that the twinning nucleation and propagation can be captured via an implicit function. While a phase-field order parameter is introduced to quantify the twinning induced shear strain and corresponding crystal reorientation, the evolution of the order parameter is driven by the resolved shear stress on the twinning system. To avoid introducing an additional set of slip systems for dislocation slip within the twinning region, we introduce a Lie algebra averaging technique to determine the Schmid tensor throughout the twinning transformation. Three different numerical schemes are proposed to solve the coupled problem, including a monolithic scheme, an alternating minimization scheme, and an operator splitting scheme. Three numerical examples are utilized to demonstrate the capability of the proposed model, as well as the accuracy and computational cost of the solvers. [preprint] There are many journal articles dedicated to constitutive models that show good matches on selected data. However, how useful are these models if we don't report or know their weakness? Our recently accepted CMAME paper uses deep reinforcement learning and game theory to explore this important question. By conceptualizing the efforts to validating and falsifying the model as a competing activities, we introduce AI agents to compete against each others to explore each others' weakness and through competition improving the resultant model and finding the Nash equilibrium that represents the limits of the state-of-the-art of each model, hand-crafted or AI-generated. The interesting aspect is that this competition can be applied to any material laws and data set. Please consider submit data/model to us and we can validate and falsify the model for you. More information can be found in the preprint [PDF]. The lectures for the entire 3-day workshop can be found at the official website of ALERT. alertgeomaterials.eu/2020/10/videos-of-the-doctoral-school-2020/ Special thanks to the organizers, Professors Wei Wu and Manuel Pastor, and the colleagues who attend the talks. PhD student Hyoung Suk Suh passed the qualification exam. His work is on modeling the thawing and freezing of frozen porous media. In the proposal talk, Hyoung Suk introduces the mathematical framework requires to simulate the freezing and thawing of porous media at the pore scale (see the rehearsal Youtube Video below). Since joining the research group in Fall 2018, Hyoung Suk has published two journal articles and one conference paper.
H.S. Suh, D. O'Conner, W.C. Sun, A phase field model for cohesive fracture in micropolar continua, Computer Methods in Applied Mechanics and Engineering, doi:10.1016/j.cma.2020.113181, 2020. H.S. Suh, W.C. Sun, An open source FEniCS implementation of a phase field fracture model for micropolar continua, International Journal of Multiscale Computational Engineering, doi:10.1615/IntJMultCompEng.2020033422, 2019. [open source code] H.S. Suh, W.C. Sun, An immersed phase field fracture model in fluid-infiltrating porous media with evolving Beavers-Joseph-Saffman condition, 2nd International Conference on Energy Geotechnics, La Jolla, California, USA, 2021. Congratulations, Hyoung Suk!  I am excited to announce that our PhD student Eric Bryant has successfully defended his dissertation "Capturing Evolving Size-Dependent Anisotropy from Brittle Fracture to Plasticity for Geological Materials" and will join Los Alamos National Laboratory as postdoc research associate. His primary appointment will be with the fluid dynamics and solid mechanics group of the Theoretical Division (T-3). During his tenure in the Sun research group, Eric has been awarded the Presidential Fellowship and the Guggenheim Fellowship. His PhD work focuses on modeling the size effect of evolving anisotropy of fracture and plastic deformation. Traditionally, one may introduce scaling law to incorporate how shear and tensile strengths depend on the size of the specimen. In his work, the main focus is also to investigate how does the rotation among the principal directions of the elastic strain, tensorial internal variables and stress tensor evolves across length scales and during the deformation process. He is the authors of 4 journal articles and one manuscript very close to completion. Among them, three of them published in Computer Methods in Applied Mechanics and Engineering, and one in International Journal of Fracture as listed below. Congratulations, Eric! Published work:
Author: Nikolas Vlassis, WaiChing Sun Abstract: We present a machine learning approach that integrates geometric deep learning and Sobolev training to generate a family of finite strain anisotropic hyperelastic models that predict the homogenized responses of polycrystals previously unseen during the training. While hand-crafted hyperelasticity models often incorporate homogenized measures of microstructural attributes, such as the porosity or the averaged orientation of constitutes, these measures may not adequately reflect the topological structures of the attributes. We fill this knowledge gap by introducing the concept of the weighted graph as a new high-dimensional descriptor that represents topological information, such as the connectivity of anisotropic grains in an assemble. By leveraging a graph convolutional deep neural network in a hybrid machine learning architecture previously used in Frankel et al. 2019, the artificial intelligence extracts low-dimensional features from the weighted graphs and subsequently learns the influence of these low-dimensional features on the resultant stored elastic energy functionals. To ensure smoothness and prevent unintentionally generating a non-convex stored energy functional, we adopt the Sobolev training method for neural networks such that a stress measure is obtained implicitly by taking directional derivatives of the trained energy functional. Results from numerical experiments suggest that Sobolev training is capable of generating a hyperelastic energy functional that predicts both the elastic energy and stress measures more accurately than the classical training that minimizes L2 norm. Verification exercises against unseen benchmark FFT simulations and phase-field fracture simulations using the geometric learning generated elastic energy functional are conducted to demonstrate the quality of the predictions. [manuscript]
Presenter: Dr. Steve WaiChing Sun, Columbia University Title: Some applications of graph theory in data-driven multi-scale mechanics Video: Watch a recording of the seminar Abstract: In this talk, we will share our experience in using undirected and directed graphs to solve computational solid mechanics problems with a variety of deep neural networks. In the first half of the talk, we will focus on the usage of undirected weighted graphs that represent the microstructures. We will demonstrate 1) how to effectively represent microstructures such as polycrystals, granular assemblies and composites as node-weighted graphs, a network of nodes with assigned attributes connected by edges, 2) how to create low-dimensional topological descriptors via graph convolution neural network that exhibits spatial and rotational invariance properties and 3) how these topological descriptors can be used to enhance the accuracy and robustness of the forward predictions and generalize the surrogate constitutive models generated via semi-supervised learning. In the second half of the talk, we will examine the application of directed multi-graphs that represent causality/relational knowledge of material laws. By idealizing the process of modeling constitutive laws as a multi-player game, we will examine 4) how the process of writing, validating and falsifying a constitutive law can be formulated as a Markov decision process, and 5) how a model-free deep reinforcement learning paradigm can introduce artificial intelligence (AI) modelers and experimentalists that learn to create hand-crafted-like constitutive models through competitions and repeated trial and errors. Examples will be provide to illustrate how these self-interacting, self-improving AI agents discover new hidden hierarchical structures of mechanics knowledge, spot the weakness of existing models and create new approaches to incorporate non-Euclidean data traditionally excluded in constitutive laws to make predictions more accurate and robust. |
Group NewsNews about Computational Poromechanics lab at Columbia University. Categories
All
Archives
January 2021
|
RSS Feed