Microvascular materials containing internal microchannels are able to achieve multi-functionality by flowing different fluids through vasculature. Active cooling is one application to protect structural components and devices from thermal overload, which is critical to modern technology including electric vehicle battery packaging and solar panels on space probes. Creating thermally efficient vascular network designs requires state-of-the-art computational tools. Prior optimization schemes have only considered steady-state cooling, rendering a knowledge gap for time-varying heat transfer behavior. In this study, a transient topology optimization framework is presented to maximize the active-cooling performance and mitigate computational cost. Here, we optimize the channel layout so that coolant flowing within the vascular network can remove heat quickly and also provide a lower steady-state temperature. An objective function for this new transient formulation is proposed that minimizes the area beneath the average temperature versus time curve to simultaneously reduce the temperature and cooling time. The thermal response of the system is obtained through a transient Geometric Reduced Order Finite Element Model (GRO-FEM). The model is verified via a conjugate heat transfer simulation in commercial software and validated by an active-cooling experiment conducted on a 3D-printed microvascular metal. A transient sensitivity analysis is derived to provide the optimizer with analytical gradients of the objective function for further computational efficiency. Example problems are solved demonstrating the method’s ability to enhance cooling performance along with a comparison of transient versus steady-state optimization results. In this comparison, both the steady-state and transient frameworks delivered different designs with similar performance characteristics for the problems considered in this study. This latest computational framework provides a new thermal regulation toolbox for microvascular material designers.
This article introduces a computational design framework for obtaining three‐dimensional (3D) periodic elastoplastic architected materials with enhanced performance, subject to uniaxial or shear strain. A nonlinear finite element model accounting for plastic deformation is developed, where a Lagrange multiplier approach is utilized to impose periodicity constraints. The analysis assumes that the material obeys a von Mises plasticity model with linear isotropic hardening. The finite element model is combined with a corresponding path‐dependent adjoint sensitivity formulation, which is derived analytically. The optimization problem is parametrized using the solid isotropic material penalization method. Designs are optimized for either end compliance or toughness for a given prescribed displacement. Such a framework results in producing materials with enhanced performance through much better utilization of an elastoplastic material. Several 3D examples are used to demonstrate the effectiveness of the mathematical framework.
more » « less- Award ID(s):
- 1926353
- NSF-PAR ID:
- 10452358
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- International Journal for Numerical Methods in Engineering
- Volume:
- 122
- Issue:
- 8
- ISSN:
- 0029-5981
- Page Range / eLocation ID:
- p. 1889-1910
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract A multiphysics phase field framework for coupled electrochemical and elastoplastic behaviors is presented, where the evolution of complex solid-electrolyte is described by the variation of the phase field variable with time. The solid-electrolyte interface kinetics nonlinearly depends on the thermodynamic driving force and can be accelerated by mechanical straining according to the film rupture-dissolution mechanism. A number of examples in two- and three- dimensions are demonstrated based on the finite element-based MOOSE framework. The model successfully captures the pit-to-crack transition under simultaneous electrochemical and mechanical effects. The crack initiation and growth has been demonstrated to depend on a variety of materials properties. The coupled corrosion and crystal plasticity framework also predict the crack initiation away from the perpendicular to the loading direction.
-
In this manuscript, we present a multiscale Adaptive Reduced-Order Modeling (AROM) framework to efficiently simulate the response of heterogeneous composite microstructures under interfacial and volumetric damage. This framework builds on the eigendeformation-based reduced-order homogenization model (EHM), which is based on the transformation field analysis (TFA) and operates in the context of computational homogenization with a focus on model order reduction of the microscale problem. EHM pre-computes certain microstructure information by solving a series of linear elastic problems defined over the fully resolved microstructure (i.e., concentration tensors, interaction tensors) and approximates the microscale problem using a much smaller basis spanned over subdomains (also called parts) of the microstructure. Using this reduced basis, and prescribed spatial variation of inelastic response fields over the parts, the microscale problem leads to a set of algebraic equations with part-wise responses as unknowns, instead of node-wise displacements as in finite element analysis. The volumetric and interfacial influence functions are calculated by using the Interface enriched Generalized Finite Element Method (IGFEM) to compute the coefficient tensors, in which the finite element discretization does not need to conform to the material interfaces. AROM takes advantage of pre-computed coefficient tensors associated with the finest ROM and efficiently computes the coefficient tensors of a series of gradually coarsening ROMs. During the multiscale analysis stage, the simulation starts with a coarse ROM which can capture the initial elastic response well. As the loading continues and response in certain parts of the microstructure starts to localize, the analysis adaptively switches to the next level of refined ROM to better capture those local responses. The performance of AROM is evaluated by comparing the results with regular EHM (no adaptive refinement) and IGFEM under different loading conditions and failure modes for various 2D and 3D microstructures. The proposed AROM provides an efficient way to model history-dependent nonlinear responses for composite materials under localized interface failure and phase damage.more » « less
-
more » « less
ABSTRACT: The chemo-mechanical loading of rocks causes the dissolution and precipitation of multiple phases in the rock. This dissolution and precipitation of load-bearing mineral phases lead to the stress redistribution in neighboring phases, which in turn results in deformational changes of the sample composite. The aim of this study is to investigate the link between microstructural evolution and creep behavior of shale rocks subjected to chemo-mechanical loading through modeling time-dependent deformation induced by the dissolution-precipitation process. The model couples the microstructural evolution of the shale rocks with the stress/strain fields inside the material as a function of time. The modeling effort is supplemented with an experimental study where shale rocks were exposed to CO2-rich brine under high temperature and pressure conditions. 3D snapshots of the sample microstructure were generated using segmented micro-CT images of the shale sample. The time-evolving microstructures were then integrated with the Finite element-based mechanical model to simulate the creep induced by dissolution and precipitation processes independent of the intrinsic viscoelasticity/viscoplasticity of the mineral phases. After computation of the time-dependent viscoelastic properties of the shale composite, the combined microstructure model and finite element model were utilized to predict the time-dependent stress and strain fields in different zones of reacted shale.
1. INTRODUCTION Determination of viscous behavior of shale rocks is key in wide range of applications such as stability of reservoirs, stability of geo-structures subjected to environmental forcing, underground storage of hazardous materials and hydraulic fracturing. Short-term creep strains in hydraulic fracturing can change stress fields and in turn can impact the hydraulic fracturing procedures(H. Sone & Zoback, 2010; Hiroki Sone & Zoback, 2013). While long-term creep strains can hamper the reservoir performance due to the reduction in permeability of the reservoir by closing of fractures and fissures(Du, Hu, Meegoda, & Zhang, 2018; Rybacki, Meier, & Dresen, 2016; Sharma, Prakash, & Abedi, 2019; Hiroki Sone & Zoback, 2014). Owing to these significance of creep strain, it is important to understand the viscoelastic/viscoplastic behavior of shales.
-
Polycrystalline materials consist of grains (crystals) oriented at different angles resulting in a heterogeneous and anisotropic mechanical behavior at that micro-length scale. In this study, a novel method is proposed for the first time to determine the [Formula: see text] crystal orientations of grains in a [Formula: see text] domain, using solely [Formula: see text] deformation fields. The grain boundaries are assumed to be unknown and delineated from the reconstructed changes in the crystallographic orientation. Further, the constitutive equations that describe the mechanical behavior of the domain in [Formula: see text] under plane stress conditions are derived, assuming that the material is transversely isotropic in 3D. Finite element based algorithms are utilized to discretize the inverse problem. The in-house written inverse problem solver is coupled with Matlab-based optimization scripts to solve for the mechanical property distributions. The performance of this method is tested at different noise levels with synthetic displacements that were used as measured data. The reconstructions deteriorate as the noise level is increased. This work presents a first milestone in the verification of this novel technology with synthetic data.more » « less
