2024-10-30 | | Total: 8
Meso-scale calculations of pore collapse and hotspot formation in energetic crystals provide closure models to macro-scale hydrocodes for predicting the shock sensitivity of energetic materials. To this end, previous works obtained atomistics-consistent material models for two common energetic crystals, HMX and RDX, such that pore collapse calculations adhered closely to molecular dynamics (MD) results on key features of energy localization, particularly the profiles of the collapsing pores, appearance of shear bands, and the transition from viscoplastic to hydrodynamic collapse. However, some important aspects such as the temperature distributions in the hotspot were not as well captured. One potential issue was noted but not resolved adequately in those works, namely the grid resolution that should be employed in the meso-scale calculations for various pore sizes and shock strengths. Conventional computational mechanics guidelines for selecting meshes as fine as possible, balancing computational effort, accuracy and grid independence, were shown not to produce physically consistent features associated with shear localization. Here, we examine the physics of pore collapse, shear band evolution and structure, and hotspot formation, for both HMX and RDX; we then evaluate under what conditions atomistics-consistent models yield physically correct (considering MD as ground truth) hotspots for a range of pore diameters, from nm to microns, and for a wide range of shock strengths. The study provides insights into the effects of pore size and shock strength on pore collapse and hotspots, identifying aspects such as size-independent behaviors, and proportion of energy contained in shear as opposed to jet impact-heated regions of the hotspot. Areas for further improvement of atomistics-consistent material models are also indicated.
Behavioral adoptions are influenced by peers in different ways. While some individuals may change after a single incoming influence, others need multiple cumulated attempts. These two mechanism, known as the simple and the complex contagions, often occur together in social phenomena alongside personal factors determining individual adoptions. Here we aim to identify, which of these contagion mechanism dominate a spreading process propagated by time-varying interactions. We consider three types of spreading scenarios: ones pre-dominated by simple or complex contagion, and mixed dynamics where the dominant mechanism changes during the unfolding of the spreading process. We propose different methods to analytically identify the transitions between these three scenarios and compare them with numerical simulations. This work offers new insights into social contagion dynamics on temporal networks, without assuming prior knowledge about individual's contagion mechanism driving their adoption decisions.
Optimizing the synthesis of zeolites and exploring novel frameworks offer pivotal opportunities and challenges in materials design. While inverse design proves highly effective for simpler crystals, its application to intricate structures like zeolites poses severe challenges. Here, we introduce an innovative inverse design workflow tailored to efficiently reproduce target zeolite frameworks in a binary coarse-grained model using enhanced sampling molecular dynamics simulations. This workflow integrates an evolutionary parameter optimization strategy with a variant of the seeding approach. Using this method, we successfully reproduce Z1 and SGT zeolites, and Type-I clathrates, find new optimal parameters for known phases, such as the SOD and CFI, and even discover novel frameworks, such as Z5. This is done within a simple coarse-grained model for a tetrahedra-forming component and a structure-directing agent. Our methodology not only enables the screening of synthesis protocols but also facilitates the discovery of hypothetical zeolites.
The accurate simulation of complex biochemical phenomena has historically been hampered by the computational requirements of high-fidelity molecular-modeling techniques. Quantum mechanical methods, such as ab initio wave-function (WF) theory, deliver the desired accuracy, but have impractical scaling for modelling biosystems with thousands of atoms. Combining molecular fragmentation with MP2 perturbation theory, this study presents an innovative approach that enables biomolecular-scale ab initio molecular dynamics (AIMD) simulations at WF theory level. Leveraging the resolution-of-the-identity approximation for Hartree-Fock and MP2 gradients, our approach eliminates computationally intensive four-center integrals and their gradients, while achieving near-peak performance on modern GPU architectures. The introduction of asynchronous time steps minimizes time step latency, overlapping computational phases and effectively mitigating load imbalances. Utilizing up to 9,400 nodes of Frontier and achieving 59% (1006.7 PFLOP/s) of its double-precision floating-point peak, our method enables us to break the million-electron and 1 EFLOP/s barriers for AIMD simulations with quantum accuracy.
Data privacy is a significant concern in many environments today. This is particularly true if sensitive information, e.g., engineering, medical, or financial data, is to be processed on potentially insecure systems, as it is often the case in cloud computing. Fully homomorphic encryption (FHE) offers a potential solution to this problem, as it allows for secure computations on encrypted data. In this paper, we investigate the viability of using FHE for privacy-preserving numerical simulations of partial differential equations. We first give an overview of the CKKS scheme, a popular FHE method for computations with real numbers. This is followed by an introduction of our Julia packages OpenFHE.jl and SecureArithmetic.jl, which provide a Julia wrapper for the C++ library OpenFHE and offer a user-friendly interface for secure arithmetic operations. We then present a performance analysis of the CKKS scheme within OpenFHE, focusing on the error and efficiency of different FHE operations. Finally, we demonstrate the application of FHE to secure numerical simulations by implementing two finite difference schemes for the linear advection equation using the SecureArithmetic.jl package. Our results show that FHE can be used to perform cryptographically secure numerical simulations, but that the error and efficiency of FHE operations must be carefully considered when designing applications.
Tropical cyclones (TCs) are powerful, natural phenomena that can severely impact populations and infrastructure. Enhancing our understanding of the mechanisms driving their intensification is crucial for mitigating these impacts. To this end, researchers are pushing the boundaries of TC simulation resolution down to scales of just a few meters. However, higher resolution simulations come with significant computational challenges, increasing both time and energy costs. Adaptive mesh refinement (AMR) is a technique widely used in computational fluid dynamics but has seen limited application in atmospheric simulations. This study explores the use of h-adaptive grids using the spectral element discretization technique to accelerate TC simulations while allowing very high resolutions in certain parts of the domain. By applying AMR to a rapidly intensifying TC test case, we demonstrate that AMR can replicate the results of uniform grid simulations in terms of mean and local wind speed maxima while dramatically reducing computational costs. We show that AMR can speed up dry TC simulations by a factor of 2 - 13 for the set of tested refinement criteria. Additionally, we show that TC intensity changes as resolution is increased and that AMR can deliver high-resolution simulations at the cost of coarser static simulations. Our findings indicate that AMR and spectral element methods are promising tools for enhancing TC simulations.
Computational physics simulation can be a powerful tool to accelerate industry deployment of new scientific technologies. However, it must address the challenge of computationally tractable, moderately accurate prediction at large industry scales, and training a model without data at such large scales. A recently proposed component reduced order modeling (CROM) tackles this challenge by combining reduced order modeling (ROM) with discontinuous Galerkin domain decomposition (DG-DD). While it can build a component ROM at small scales that can be assembled into a large scale system, its application is limited to linear physics equations. In this work, we extend CROM to nonlinear steady Navier-Stokes flow equation. Nonlinear advection term is evaluated via tensorial approach or empirical quadrature procedure. Application to flow past an array of objects at moderate Reynolds number demonstrates $\sim23.7$ times faster solutions with a relative error of $\sim 2.3\%$, even at scales $256$ times larger than the original problem.
Scaling up new scientific technologies from laboratory to industry often involves demonstrating performance on a larger scale. Computer simulations can accelerate design and predictions in the deployment process, though traditional numerical methods are computationally intractable even for intermediate pilot plant scales. Recently, component reduced order modeling method is developed to tackle this challenge by combining projection reduced order modeling and discontinuous Galerkin domain decomposition. However, while many scientific or engineering applications involve nonlinear physics, this method has been only demonstrated for various linear systems. In this work, the component reduced order modeling method is extended to steady Navier-Stokes flow, with application to general nonlinear physics in view. Large-scale, global domain is decomposed into combination of small-scale unit component. Linear subspaces for flow velocity and pressure are identified via proper orthogonal decomposition over sample snapshots collected at small scale unit component. Velocity bases are augmented with pressure supremizer, in order to satisfy inf-sup condition for stable pressure prediction. Two different nonlinear reduced order modeling methods are employed and compared for efficient evaluation of nonlinear advection: 3rd-order tensor projection operator and empirical quadrature procedure. The proposed method is demonstrated on flow over arrays of five different unit objects, achieving $23$ times faster prediction with less than $4\%$ relative error up to $256$ times larger scale domain than unit components. Furthermore, a numerical experiment with pressure supremizer strongly indicates the need of supremizer for stable pressure prediction. A comparison between tensorial approach and empirical quadrature procedure is performed, which suggests a slight advantage for empirical quadrature procedure.