Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coarse grained simulations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

Dynamics and lipid membrane coupling of the RAS-RAF complex revealed via multiscale simulations

To gain molecular and mechanistic insights into initiation of the RAS-RAF signaling cascade, we developed and used a combination of multiscale simulation and experimental approaches. The influence and impact of the membrane on RAS and RAF proteins is a factor we are just beginning to understand and appreciate in more detail. Molecular simulation is an ideal methodology to further study this complicated relationship between the membrane and associated proteins. Our previous work using Multiscale Machine-learned Modeling Infrastructure investigated different lipid compositions solely around the KRAS4b protein and the interplay between protein behavior and these membrane environments. Multiscale Machine-learned Modeling Infrastructure uses machine learning to couple adjacent simulation scales and has been efficiently scaled across some of the world’s largest high-performance computers. Recently, we have expanded this multiresolution framework to include the all-atom simulation scale and to incorporate the RAF RBDCRD domains. Here, we present the overall analysis results from this new simulation campaign comprising a mixture of RAS and RAF RBDCRD proteins. Approximately 35,000 coarse-grained and 10,000 all-atom molecular dynamics simulations were completed, sampled from a variety of protein/lipid composition configurations that were generated from a micron-scale continuum simulation containing hundreds of copies of the proteins. Our studies suggest that orientations of the RAS-RBDCRD complex on the membrane occupy distinct configurational states, and the spatial patterns of lipid arrangements around these different protein states are unique to each state. The extent and size of lipid “fingerprints” imposed on the membrane by the RAS-RBDCRD protein complex are significantly larger than observed for just the RAS protein on its own. These protein complexes strongly associate, but we do not observe statistically significant preferred protein-protein orientations. These observations indicate that spatial colocalization of RAS-RBDCRD proteins in the same vicinity may be assisted by specific membrane environments, acting to increase the probability of signaling complex formation.

Carpenter, Timothy S. [Lawrence Livermore National↗

Uncovering Backbone Conformation for Rigid DPP-Based Donor–Acceptor Conjugated Polymer Using Deuterium Labeling and Neutron Scattering

The conjugated polymer’s backbone conformation dictates the delocalization of electrons, ultimately affecting its optoelectronic properties. Most conjugated polymers can be viewed as semirigid rods with their backbone embedded among long alkyl side chains. Thus, it is challenging to experimentally quantify the conformation of a conjugated backbone. Here, we performed contrast variation neutron scattering on rigid conjugated donor–acceptor (D–A) diketopyrrolopyrrole (DPP) polymers with selectively deuterated side chains to measure the conjugated backbone conformation. We first synthesized DPPbased polymers with deuterated side chains, confirmed by NMR and FTIR. Using contrast variation neutron scattering, we found that the DPP-based conjugated polymers are much more rigid than poly(3-alkylthiophenes), with persistence length (L p ) at 16–18 nm versus 2–3 nm. More importantly, in contrast to the relatively flexible poly(3-alkylthiophenes) whose backbone is more flexible than the whole polymer, we found that the backbone of DPP-based polymers has the same L p value compared to the whole polymer chain. This indicates that side chain interference on backbone conformation is not present for the semirigid polymer, which is further confirmed by coarse-grained molecular dynamics (CG-MD) simulations. Our work provides a novel protocol to probe polymer’s backbone conformation and paradigm-shifting understanding of the backbone conformation of semirigid conjugated polymers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Implementation and Characterization of Three-Dimensional Particle-in-Cell Codes on Multiple-Instruction-Multiple-Data Massively Parallel Supercomputers

A three-dimensional electrostatic particle-in-cell (PIC) plasma simulation code has been developed on coarse-grain distributed-memory massively parallel computers with message passing communications. Our implementation is the generalization to three-dimensions of the general concurrent particle-in-cell (GCPIC) algorithm. In the GCPIC algorithm, the particle computation is divided among the processors using a domain decomposition of the simulation domain. In a three-dimensional simulation, the domain can be partitioned into one-, two-, or three-dimensional subdomains ("slabs," "rods," or "cubes") and we investigate the efficiency of the parallel implementation of the push for all three choices. The present implementation runs on the Intel Touchstone Delta machine at Caltech; a multiple-instruction-multiple-data (MIMD) parallel computer with 512 nodes. We find that the parallel efficiency of the push is very high, with the ratio of communication to computation time in the range 0.3%-10.0%. The highest efficiency (> 99%) occurs for a large, scaled problem with 64(sup 3) particles per processing node (approximately 134 million particles of 512 nodes) which has a push time of about 250 ns per particle per time step. We have also developed expressions for the timing of the code which are a function of both code parameters (number of grid points, particles, etc.) and machine-dependent parameters (effective FLOP rate, and the effective interprocessor bandwidths for the communication of particles and grid points). These expressions can be used to estimate the performance of scaled problems--including those with inhomogeneous plasmas--to other parallel machines once the machine-dependent parameters are known.

Lyster, P. M.↗

Multiscale Modeling of Vinyl-Addition Polynorbornenes: The Effect of Stereochemistry

Vinyl-addition polynorbornenes are candidates for designing high-performance polymers due to unique characteristics, which include a high glass transition temperature associated with a rigid backbone. Recent studies have established that the processability and properties of these polymers can be fine-tuned by using targeted substitutions. However, synthesis with different catalysts results in materials with distinct properties, potentially due to the presence of various stereoisomers that are difficult to quantify experimentally. Herein, we develop all-atom models of polynorbornene oligomers based on classical force fields and density functional theory. To establish the relationship between chemical architecture, chain conformations, and melt structure, we perform detailed molecular dynamics simulations with the fine-tuned atomistic force field and propose simpler coarse-grained descriptions to address the high molecular weight limit. All-atom simulations of oligomers suggest high glass transition temperatures in the range of 550–600 K. In the melt state (800 K), meso chains form highly rigid extended coils (C∞≈11) with amorphous structural characteristics similar to the X-ray diffraction data observed in the literature. In contrast, simulations with racemo chains predict highly helical tubular chain conformations that could promote assembly into crystalline structures.

Polymer Science↗

Data-driven particle dynamics: Structure-preserving coarse-graining for emergent behavior in non-equilibrium systems

Multiscale systems are ubiquitous in science and technology, but are notoriously challenging to simulate as short spatiotemporal scales must be appropriately linked to emergent bulk physics. When expensive high-dimensional dynamical systems are coarse-grained into low-dimensional models, the entropic loss of information leads to emergent physics which are dissipative, history-dependent, and stochastic. To machine learn coarse-grained dynamics from time-series observations of particle trajectories, we propose a framework using the metriplectic bracket formalism that preserves these properties by construction; most notably, the framework guarantees discrete notions of the first and second laws of thermodynamics, conservation of momentum, and a discrete fluctuation-dissipation balance crucial for capturing non-equilibrium statistics. We introduce the mathematical framework abstractly before specializing to a particle discretization. As labels are generally unavailable for entropic state variables, we introduce a novel self-supervised learning strategy to identify emergent structural variables. We validate the method on benchmark systems and demonstrate its utility on two challenging examples: (1) coarse-graining star polymers at challenging levels of coarse-graining while preserving non-equilibrium statistics, and (2) learning models from high-speed video of colloidal suspensions that capture coupling between local rearrangement events and emergent stochastic dynamics. We provide open-source implementations in both PyTorch and LAMMPS, enabling large-scale inference and extensibility to diverse particle-based systems.

Computational Engineering, Finance, and Science (c↗

The Martini 3 Lipidome: Expanded and Refined Parameters Improve Lipid Phase Behavior

Lipid membranes are central to cellular life. Complementing experiments, computational modeling has been essential in unraveling complex lipid-biomolecule interactions, crucial in both academia and industry. The Martini model, a coarse-grained force field for efficient molecular dynamics simulations, is widely used to study membrane phenomena but has faced limitations, particularly in capturing realistic lipid phase behavior. Here, we present refined Martini 3 lipid models with a mapping scheme that distinguishes lipid tails that differ by just two carbon atoms, enhancing the structural resolution and thermodynamic accuracy of model membrane systems including ternary mixtures. The expanded Martini lipid library includes thousands of models, enabling simulations of complex and biologically relevant systems. These advancements establish Martini as a robust platform for lipid-based simulations across diverse fields.

Lipids↗

Thermomechanical Modeling of Woven Materials With Particle-Based, Explicit-Fiber Simulations

Fiber-based materials are extensively used to protect spacecraft during entry. Insulative fibers, often in a fiber network or woven, provide rigidity, strength, and control of material anisotropy and density. Woven thermal protection materials, such as ADEPT (Adaptable, Deployable Entry and Placement Technology), 3D-MAT (3-Dimensional Multifunctional Ablative Thermal Protection), and 3MDCP (3D Woven Mid-Density Carbon Phenolic), enable missions with stronger and denser materials for entry profiles with high shear and heat flux. Vulnerabilities to woven thermal protection materials include manufacturing-induced material property variation, and impact from micrometeoroids. Simulating woven materials under these conditions require models that can resolve hierarchal structures, thermomechanical behavior, and failure. To address this, we simulate weave thermal conduction and mechanical deformation. We simulate the full weave with a coarse-grained yarn model is presented. The model combines a validated, high-resolution single 3MDCP yarn model and phenolic resin model. Instead of modeling every fiber, each yarn ply with order 10, instead of order 1000, fibers. The discrete element bonded particle model (DEM-BPM) of fibers captures the thermal and mechanical behavior within and between fibers. We study the proportion of heat transfer and stress via the contact network, fiber bonds, and overall weave geometry.

bonded particle↗

Perpendicular crossing chains enable high mobility in a noncrystalline conjugated polymer

The nature of interchain π-system contacts, and their relationship to hole transport, are elucidated for the high-mobility, noncrystalline conjugated polymer C16-IDTBT by the application of scanning tunneling microscopy, molecular dynamics, and quantum chemical calculations. The microstructure is shown to favor an unusual packing motif in which paired chains cross-over one another at near-perpendicular angles. By linking to mesoscale microstructural features, revealed by coarse-grained molecular dynamics and previous studies, and performing simulations of charge transport, it is demonstrated that the high mobility of C16-IDTBT can be explained by the promotion of a highly interconnected transport network, stemming from the adoption of perpendicular contacts at the nanoscale, in combination with fast intrachain transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

3D pattern formation of a protein–membrane suspension

Many essential cellular processes, including cell division and the establishment of cell polarity during embryogenesis, are regulated by pattern-forming proteins. These proteins often need to bind to a substrate, such as the cell membrane, onto which they interact and form two-dimensional (2D) patterns. It is unclear how the membrane’s continuity and dimensionality impact pattern formation. Here, we address this gap using the MinDE system, a prototypical example of pattern-forming membrane proteins. We show that when the lipid substrate is fragmented into submicrometer-sized diffusive liposomes, adenosine triphosphate-driven protein–protein interactions generate three-dimensional (3D) spatially extended patterns, despite the complete loss of membrane continuity. Remarkably, these 3D patterns emerge at scales four orders of magnitude larger than the individual liposomes. By systematically varying protein concentration, liposome size, and density, we observed and characterized a variety of 3D dynamical patterns not seen on continuous 2D membranes, including traveling waves, dynamical spirals, and a coexistence phase. Simulations and linear stability analysis of a coarse-grained model revealed that the physical properties of the dispersed membrane effectively rescale both the protein–membrane binding rates and diffusion, two key parameters governing pattern formation and wavelength selection. These findings highlight the robustness of Min’s pattern-forming ability, suggesting that protein–membrane suspensions could serve as an adaptable template for studying out-of-equilibrium self-organization in 3D, beyond in vivo contexts.

36 MATERIALS SCIENCE↗

Direct particle simulation on the Connection Machine CM-2

Particle simulation is a useful technique for analyzing low density flows. The Connection Machine CM-2 is a useful test bed for studying the fine-grain data objects decomposition and the coarse-grain domain decomposition single instruction multiple datastream (SIMD) approaches to particle simulation. Both approaches are investigated for the model problem of uniform flow through a channel and the algorithms required for the SIMD domain decomposition approach are presented. An unresolved issue with the domain decomposition approach is the effect of a poor partitioning on flows with real geometries. Initial results with the channel flow problem indicate that a poor partitioning has only a small detrimental effect on the overall performance.

Dagum, L.↗

A polymer, random walk model for the size-distribution of large DNA fragments after high linear energy transfer radiation

DNA double-strand breaks (DSBs) produced by densely ionizing radiation are not located randomly in the genome: recent data indicate DSB clustering along chromosomes. Stochastic DSB clustering at large scales, from > 100 Mbp down to < 0.01 Mbp, is modeled using computer simulations and analytic equations. A random-walk, coarse-grained polymer model for chromatin is combined with a simple track structure model in Monte Carlo software called DNAbreak and is applied to data on alpha-particle irradiation of V-79 cells. The chromatin model neglects molecular details but systematically incorporates an increase in average spatial separation between two DNA loci as the number of base-pairs between the loci increases. Fragment-size distributions obtained using DNAbreak match data on large fragments about as well as distributions previously obtained with a less mechanistic approach. Dose-response relations, linear at small doses of high linear energy transfer (LET) radiation, are obtained. They are found to be non-linear when the dose becomes so large that there is a significant probability of overlapping or close juxtaposition, along one chromosome, for different DSB clusters from different tracks. The non-linearity is more evident for large fragments than for small. The DNAbreak results furnish an example of the RLC (randomly located clusters) analytic formalism, which generalizes the broken-stick fragment-size distribution of the random-breakage model that is often applied to low-LET data.

DNA/radiation effects↗

Probabilistic Design of a Wind Tunnel Model to Match the Response of a Full-Scale Aircraft

approach is presented for carrying out the reliability-based design of a plate-like wing that is part of a wind tunnel model. The goal is to design the wind tunnel model to match the stiffness characteristics of the wing box of a flight vehicle while satisfying strength-based risk/reliability requirements that prevents damage to the wind tunnel model and fixtures. The flight vehicle is a modified F/A-18 aircraft. The design problem is solved using reliability-based optimization techniques. The objective function to be minimized is the difference between the displacements of the wind tunnel model and the corresponding displacements of the flight vehicle. The design variables control the thickness distribution of the wind tunnel model. Displacements of the wind tunnel model change with the thickness distribution, while displacements of the flight vehicle are a set of fixed data. The only constraint imposed is that the probability of failure is less than a specified value. Failure is assumed to occur if the stress caused by aerodynamic pressure loading is greater than the specified strength allowable. Two uncertain quantities are considered: the allowable stress and the thickness distribution of the wind tunnel model. Reliability is calculated using Monte Carlo simulation with response surfaces that provide approximate values of stresses. The response surface equations are, in turn, computed from finite element analyses of the wind tunnel model at specified design points. Because the response surface approximations were fit over a small region centered about the current design, the response surfaces were refit periodically as the design variables changed. Coarse-grained parallelism was used to simultaneously perform multiple finite element analyses. Studies carried out in this paper demonstrate that this scheme of using moving response surfaces and coarse-grained computational parallelism reduce the execution time of the Monte Carlo simulation enough to make the design problem tractable. The results of the reliability-based designs performed in this paper show that large decreases in the probability of stress-based failure can be realized with only small sacrifices in the ability of the wind tunnel model to represent the displacements of the full-scale vehicle.

Mason, Brian H.↗

Universal Spreading of Conditional Mutual Information in Noisy Random Circuits

For this work, we study the evolution of conditional mutual information (CMI) in generic open quantum systems, focusing on one-dimensional random circuits with interspersed local noise. Unlike in noiseless circuits, where CMI spreads linearly while being bounded by the light cone, we find that noisy random circuits with an error rate 𝑝 exhibit superlinear propagation of CMI, which diverges far beyond the light cone at a critical circuit depth 𝑡 𝑐 ∝ 𝑝 −1 . We demonstrate that the underlying mechanism for such rapid spreading is the combined effect of local noise and a scrambling unitary, which selectively removes short-range correlations while preserving long-range correlations. To analytically capture the dynamics of CMI in noisy random circuits, we introduce a coarse-graining method, and we validate our theoretical results through numerical simulations. Furthermore, we identify a universal scaling law governing the spreading of CMI.

decoherence↗

KRAS4a and KRAS4b show distinct lipid-dependent regulation of RAS-RAF membrane dynamics

KRAS4a and KRAS4b are important regulators of signaling, and their interactions with the plasma membrane are dynamic and influenced by lipid composition. KRAS 4a and 4b have nearly identical globular domains but differ in their membrane-associated hyper variable region (HVR). The functional distinctions between these isoforms remain unclear, particularly with regards to their dependence on specific lipids and the membrane environment. Previous work showed that the membrane orientation of KRAS4b affects its ability to bind to RAF kinase RBDCRD and that the KRAS–RBDCRD complex adopts different poses on the membrane as well as influences the size and composition of the lipid environment. To model differences between KRAS 4a and 4b protein–lipid interactions, we extended the Multiscale Machine-Learned Modeling Infrastructure (MuMMI) to incorporate continuum simulations in the grand canonical ensemble, enabling sampling across macroscopic, coarse-grained, and all-atom resolutions. Using this framework, we systematically altered PIP2 concentrations, KRAS 4a versus 4b, and RAF RBDCRD complexation to assess impacts on membrane–protein interactions and dynamics. Our results reveal that reducing PIP2 shifts and broadens the membrane orientational preference of both KRAS 4b and 4a, with stronger effects on 4b HVR localization versus 4a. We demonstrate that with depletion of the strong negatively charged PIP2 lipid, the less charged phosphatidylserine replaces PIP2. Our findings highlight similarities and distinctions in the dynamics and lipid dependency of KRAS isoforms and suggest that ordering of the local lipid composition by HVRs is a shared property and key modulator of RAS-mediated signaling at the plasma membrane.

Biological and medical sciences↗

Group implicit concurrent algorithms in nonlinear structural dynamics

During the 70's and 80's, considerable effort was devoted to developing efficient and reliable time stepping procedures for transient structural analysis. Mathematically, the equations governing this type of problems are generally stiff, i.e., they exhibit a wide spectrum in the linear range. The algorithms best suited to this type of applications are those which accurately integrate the low frequency content of the response without necessitating the resolution of the high frequency modes. This means that the algorithms must be unconditionally stable, which in turn rules out explicit integration. The most exciting possibility in the algorithms development area in recent years has been the advent of parallel computers with multiprocessing capabilities. So, this work is mainly concerned with the development of parallel algorithms in the area of structural dynamics. A primary objective is to devise unconditionally stable and accurate time stepping procedures which lend themselves to an efficient implementation in concurrent machines. Some features of the new computer architecture are summarized. A brief survey of current efforts in the area is presented. A new class of concurrent procedures, or Group Implicit algorithms is introduced and analyzed. The numerical simulation shows that GI algorithms hold considerable promise for application in coarse grain as well as medium grain parallel computers.

Ortiz, M.↗

ezAlign: A Tool for Converting Coarse-Grained Molecular Dynamics Structures to Atomistic Resolution for Multiscale Modeling

Soft condensed matter is challenging to study due to the vast time and length scales that are necessary to accurately represent complex systems and capture their underlying physics. Multiscale simulations are necessary to study processes that have disparate time and/or length scales, which abound throughout biology and other complex systems. Herein we present ezAlign, an open-source software for converting coarse-grained molecular dynamics structures to atomistic representation, allowing multiscale modeling of biomolecular systems. The ezAlign v1.1 software package is publicly available for download at github.com/LLNL/ezAlign. Its underlying methodology is based on a simple alignment of an atomistic template molecule, followed by position-restraint energy minimization, which forces the atomistic molecule to adopt a conformation consistent with the coarse-grained molecule. The molecules are then combined, solvated, minimized, and equilibrated with position restraints. Validation of the process was conducted on a pure POPC membrane and compared with other popular methods to construct atomistic membranes. Additional examples, including surfactant self-assembly, membrane proteins, and more complex bacterial and human plasma membrane models, are also presented. By providing these examples, parameter files, code, and an easy-to-follow recipe to add new molecules, this work will aid future multiscale modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A coarse-grained model of clay colloidal aggregation and consolidation with explicit representation of the electrical double layer

The aggregation of clay minerals in liquid water exemplifies colloidal self-assembly in nature. These negatively charged aluminosilicate platelets interact through multiple mechanisms with different sensitivities to particle shape, surface charge, aqueous chemistry, and interparticle distance and exhibit complex aggregation structures. Experiments have difficulty resolving the associated colloidal assemblages at the scale of individual particles. Conversely, all-atom molecular dynamics (MD) simulations provide detailed insight on clay colloidal interaction mechanisms, but they are limited to systems containing a few particles. We develop a new coarse-grained (CG) model capable of representing assemblages of hundreds of clay particles with accuracy approaching that of MD simulations, at a fraction of the computational cost. Our CG model is parameterized based on MD simulations of a pair of smectite clay particles in liquid water. A distinctive feature of our model is that it explicitly represents the electrical double layer (EDL), i.e., the cloud of charge-compensating cations that surrounds the clay particles. Our model captures the simultaneous importance of long-range colloidal interactions (i.e., interactions consistent with simplified analytical models, already included in extant clay CG models) and short-range interactions such as ion correlation and surface and ion hydration effects. The resulting simulations correctly predict, at low solid-water ratios, the existence of ordered arrangements of parallel particles separated by water films with a thickness up to ~10 nm and, at high solid-water ratios, the coexistence of crystalline and osmotic swelling states, in agreement with experimental observations.

54 ENVIRONMENTAL SCIENCES↗

Versatile Turbine Disk Alloy Designed and Processed for Higher Temperature Applications

Advanced turbine engine configurations using higher temperature combustor and airfoil concepts require compressor and turbine disks that can withstand temperatures higher than the 650 C typical of current engines. This requires disk alloy and processing improvements. A versatile disk alloy was needed that could be given either fine grain heat treatments for high strength and fatigue and creep resistance up to 704 C or given coarse grain heat treatments for high strength and creep and dwell crack-growth resistance at higher temperatures. This alloy could produce disks with uniform microstructures and properties, as well as with varied, optimized microstructures at disk bore and rim sections. A series of experimental disk alloys were designed and processed to give subscale disks about 13 cm in diameter and 4 cm thick. These alloys had varying levels of key elements to affect mechanical properties and change the "solvus" solution heat-treatment temperature necessary to produce coarse grain microstructures. Disks were given coarse grain heat treatments followed by rapid oil quenching and slower fan air quenching. These heat treatments were intended to simulate the cooling paths of rapidly cooled full-scale disks at the outermost rim and interior bore locations, respectively. Preliminary quench tests of tensile specimens and coin-size minidisks had indicated that alloys having high solvus temperatures were more prone to cracking during rapid quenching from coarse grain heat treatments. These findings were confirmed in the subscale disks, where such alloys did form undesirable quench cracks. Mechanical tests were performed on specimens from subscale disks given these coarse grain heat treatments, as well as on specimens given a fine grain heat treatment. Tensile, creep, and crack growth tests were performed at 704 C and higher temperatures. A versatile alloy was identified that had a low solvus temperature for resistance to quench cracking as well as an optimal combination of high levels of strengthening refractory elements that produced balanced high mechanical properties for both fine grain and coarse grain microstructures. This low-solvus, high-refractory (LSHR) alloy has been scaled-up to produce prototype full-scale turbine disks typical of regional jet turbofan engines. Disks were successfully heat treated to give uniform coarse grain and uniform fine grain microstructures. Additional disks were given a NASA dual microstructure heat treatment (DMHT) that intentionally varied the solution heat-treatment temperatures between the disk rim and bore (re f. 1). The disk rim was heated to a high enough temperature to produce a coarse grain microstructure, while the bore was maintained at a lower temperature to produce a fine grain microstructure (see the figure). This DMHT can produce optimal high strength, fatigue, and creep resistance up to 704 C in the cooler running disk bore, and high strength, creep resistance, and dwell crack growth resistance at higher temperatures for the hotter disk rim. Extensive mechanical testing is being initiated to compare the mechanical properties of the uniform and DMHT disks.

Gabb, TImothy P.↗