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
A diffuse modeling approach for embedded interfaces in linear elasticity
In this contribution, we present a diffuse modeling approach to embed material interfaces into nonconforming meshes with a focus on linear elasticity. For this purpose, a regularized indicator function is employed that describes the distribution of the different materials by a scalar value. The material in the resulting diffuse interface region is redefined in terms of this indicator function and recomputed by a homogenization of the adjacent material parameters. The applied homogenization method fulfills the kinematic compatibility across the interface and the static equilibrium at the interface. In addition, an
- NSF-PAR ID:
- 10116199
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- GAMM-Mitteilungen
- Volume:
- 43
- Issue:
- 1
- ISSN:
- 0936-7195
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Over the two decades, amorphous oxide semiconductors (AOSs) and their thin film transistor (TFT) channel application have been intensely explored to realize high performance, transparent and flexible displays due to their high field effect mobility (μFE=5-20 cm2/Vs), visible range optical transparency, and low temperature processability (25-300 °C).[1-2] The metastable amorphous phase is to be maintained during operation by the addition of Zn and additional third cation species (e.g., Ga, Hf, or Al) as an amorphous phase stabilizer.[3-5] To limit TFT off-state currents, a thin channel layer (10-20 nm) was employed for InZnO (IZO)-based TFTs, or third cations were added to suppress carrier generations in the TFT channel. To resolve bias stress-induced instabilities in TFT performance, approaches to employ defect passivation layers or enhance channel/dielectric interfacial compatibility were demonstrated.[6-7] Metallization contact is also a dominating factor that determines the performance of TFTs. Particularly, it has been reported that high electrical contact resistance significantly sacrifices drain bias applied to the channel, which leads to undesirable power loss during TFT operation and issues for the measurement of TFT field effect mobilities. [2, 8] However, only a few reports that suggest strategies to enhance contact behaviors are available in the literature. Furthermore, the previous approaches (1) require an additional fabrication complexity due to the use of additional treatments at relatively harsh conditions such as UV, plasma, or high temperatures, and (2) may lead to adverse effects on the channel material attributed to the chemical incompatibility between dissimilar materials, and exposures to harsh environments. Therefore, a simple and easy but effective buffer strategy, which does not require any additional process complexities and not sacrifice chemical compatibility, needs to be established to mitigate the contact issues and therefore achieve high performance and low power consumption AOS TFTs. The present study aims to demonstrate an approach utilizing an interfacial buffer layer, which is compositionally homogeneous to the channel to better align work functions between channel and metallization without a significant fabrication complexity and harsh treatment conditions. Photoelectron spectroscopic measurements reveal that the conducting IZO buffer, of which the work function (Φ) is 4.37 eV, relaxes a relatively large Φ difference between channel IZO (Φ=4.81 eV) and Ti (Φ=4.2-4.3 eV) metallization. The buffer is found to lower the energy barrier for charge carriers at the source to reach the effective channel region near the dielectric. In addition, the higher carrier density of the buffer and favorable chemical compatibility with the channel (compositionally the same) further contribute to a significant reduction in specific contact resistance as much as more than 2.5 orders of magnitude. The improved contact and carrier supply performance from the source to the channel lead to an enhanced field effect mobility of up to 56.49 cm2/Vs and a threshold voltage of 1.18 V, compared to 13.41 cm2/Vs and 7.44 V of IZO TFTs without a buffer. The present work is unique in that an approach to lower the potential barrier between the source and the effective channel region (located near the channel/dielectric interface, behaving similar to a buried-channel MOSFET [9]) by introducing a contact buffer layer that enhances the field effect mobility and facilitates carrier supply from the source to the effective channel region.more » « less
-
more » « less
Abstract Magnetic fields are dynamically important in the diffuse interstellar medium. Understanding how gravitationally bound, star-forming clouds form requires modeling of the fields in a self-consistent, supernova-driven, turbulent, magnetized, stratified disk. We employ the FLASH magnetohydrodynamics code to follow the formation and early evolution of clouds with final masses of 3–8 × 103
M ⊙within such a simulation. We use the code’s adaptive mesh refinement capabilities to concentrate numerical resolution in zoom-in regions covering single clouds, allowing us to investigate the detailed dynamics and field structure of individual self-gravitating clouds in a consistent background medium. Our goal is to test the hypothesis that dense clouds are dynamically evolving objects far from magnetohydrostatic equilibrium. We find that the cloud envelopes are magnetically supported with field lines parallel to density gradients and flow velocity, as indicated by the histogram of relative orientations and other statistical measures. In contrast, the dense cores of the clouds are gravitationally dominated, with gravitational energy exceeding internal, kinetic, or magnetic energy and accelerations due to gravity exceeding those due to magnetic or thermal pressure gradients. In these regions, field directions vary strongly, with a slight preference toward being perpendicular to density gradients, as shown by three-dimensional histograms of relative orientation. -
more » « less
Interactions between an evolving solid and inviscid flow can result in substantial computational complexity, particularly in circumstances involving varied boundary conditions between the solid and fluid phases. Examples of such interactions include melting, sublimation, and deflagration, all of which exhibit bidirectional coupling, mass/heat transfer, and topological change of the solid–fluid interface. The diffuse interface method is a powerful technique that has been used to describe a wide range of solid-phase interface-driven phenomena. The implicit treatment of the interface eliminates the need for cumbersome interface tracking, and advances in adaptive mesh refinement have provided a way to sufficiently resolve diffuse interfaces without excessive computational cost. However, the general scale-invariant coupling of these techniques to flow solvers has been relatively unexplored. In this work, a robust method is presented for treating diffuse solid–fluid interfaces with arbitrary boundary conditions. Source terms defined over the diffuse region mimic boundary conditions at the solid–fluid interface, and it is demonstrated that the diffuse length scale has no adverse effects. To show the efficacy of the method, a one-dimensional implementation is introduced and tested for three types of boundaries: mass flux through the boundary, a moving boundary, and passive interaction of the boundary with an incident acoustic wave. Two-dimensional results are presented as well these demonstrate expected behavior in all cases. Convergence analysis is also performed and compared against the sharp-interface solution, and linear convergence is observed. This method lays the groundwork for the extension to viscous flow and the solution of problems involving time-varying mass-flux boundaries.
-
Abstract Adaptive mesh refinement (AMR) is the art of solving PDEs on a mesh hierarchy with increasing mesh refinement at each level of the hierarchy. Accurate treatment on AMR hierarchies requires accurate prolongation of the solution from a coarse mesh to a newly defined finer mesh. For scalar variables, suitably high-order finite volume WENO methods can carry out such a prolongation. However, classes of PDEs, such as computational electrodynamics (CED) and magnetohydrodynamics (MHD), require that vector fields preserve a divergence constraint. The primal variables in such schemes consist of normal components of the vector field that are collocated at the faces of the mesh. As a result, the reconstruction and prolongation strategies for divergence constraint-preserving vector fields are necessarily more intricate. In this paper we present a fourth-order divergence constraint-preserving prolongation strategy that is analytically exact. Extension to higher orders using analytically exact methods is very challenging. To overcome that challenge, a novel WENO-like reconstruction strategy is invented that matches the moments of the vector field in the faces, where the vector field components are collocated. This approach is almost divergence constraint-preserving, therefore, we call it WENO-ADP. To make it exactly divergence constraint-preserving, a touch-up procedure is developed that is based on a constrained least squares (CLSQ) method for restoring the divergence constraint up to machine accuracy. With the touch-up, it is called WENO-ADPT. It is shown that refinement ratios of two and higher can be accommodated. An item of broader interest in this work is that we have also been able to invent very efficient finite volume WENO methods, where the coefficients are very easily obtained and the multidimensional smoothness indicators can be expressed as perfect squares. We demonstrate that the divergence constraint-preserving strategy works at several high orders for divergence-free vector fields as well as vector fields, where the divergence of the vector field has to match a charge density and its higher moments. We also show that our methods overcome the late time instability that has been known to plague adaptive computations in CED.more » « less
