SEARCH · Engineering Papers
Results for “coarse-grained simulation”
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.
A lumped particle direct simulation Monte-Carlo method combined with the collisional-radiative model for simulations of non-equilibrium laser-induced plasma plumes
Collisional plasma plumes induced by laser irradiation of material targets exhibit large variations in local density as well as ionization and excitation states, making purely hydrodynamic or kinetic simulations inaccurate or infeasible. To address this challenge and capture non-equilibrium effects in laser-induced plasma plumes at arbitrary degrees of ionization, we develop a hybrid computational approach that combines the kinetic direct simulation Monte Carlo (DSMC) method with a collisional-radiative model (CRM). This ℓDSMC-CRM approach utilizes a lumped particle method to represent minor fractions of excited ions in particle-based simulations and a special coarse-graining technique for atomic spectra and photoionization rates, ensuring numerical convergence at reduced computational cost. The hybrid approach is applied to simulate spatially homogeneous relaxation as well as one- and two-dimensional expansions of plasma plumes induced by irradiation of a copper target by a nanosecond laser pulse in a vacuum or background gas. The comparison with an equilibrium model, where local Saha-Boltzmann equilibrium is enforced, shows that the non-equilibrium effects play a dominant role. The equilibrium model can fail to predict the flow structure and strongly underestimate the degree of absorption of laser radiation by the plume. The ℓDSMC-CRM approach is validated against experimental data demonstrating reasonable agreement with the experimental electron density and temperature, while the equilibrium model is found to dramatically underestimate electron density and temperature. The flexibility of the ℓDSMC-CRM approach allows for its seamless integration into existing DSMC frameworks, making it a valuable tool for high-fidelity plasma modeling in laser-material interactions, laser-based manufacturing, and beyond.
Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse-Grained Model with Explicit Counterions
Clay-rich geomaterials play a critical role in many subsurface systems. The macroscale properties of these materials (low permeability, high ionic conductivity, high swelling pressure, etc.) are sensitive to molecular-level adsorption and hydration interactions at clay−water interfaces. Efforts to develop multiscale simulation approaches to predict these properties reveal a scale gap between atomistic simulations (typically limited to systems smaller than 10 nm) and continuum-scale models (which use computational grid elements with dimensions ≳ 10 μm). In this study, we present a coarse-grained (CG) framework that partly bridges this gap by simulating compacted smectite clay assemblages with dimensions of 0.1 μm containing 1,400 clay particles across a range of dry densities (1,050 to 1,850 kg·m −3 ) and Na/Ca counterion compositions (Na fraction ranging from 0.2 to 1). The simulated systems, along with their reconstructed binary three-dimensional pore networks, are used to evaluate the microstructure, pore size distribution, tortuosity, ion diffusivity, and swelling pressure of compacted smectite clay. Results show that our approach captures important features of the mesoscale heterogeneity of compacted clays, including tactoid formation, hierarchical porosity, and anisotropic pore networks. Results also reveal how compaction and counterion composition govern emergent behaviors, including dominant pore sizes, directional transport, and electrochemical response. This work highlights the potential of CG simulations to bridge molecular and continuum scales and to advance geotechnical and environmental applications involving clay-rich materials as well as related nanoporous media such as geopolymers and calcium-silicate-hydrate. However, the results also suggest that accurate prediction of certain microstructural and mechanical properties (e.g., swelling pressure) may require even larger-scale systems on the order of 1 μm.
The Complexation Properties of Self-Defensive Microgel-Modified Antimicrobial Surfaces
The complexation of cationic antimicrobials with polyanionic microgels on a biomaterial surface can render that surface self-defensive against bacteria by killing those bacteria which physically contact the antimicrobial-loaded microgels. This killing has been attributed to the contact-driven transfer of antimicrobial from a microgel to a challenging bacterium, though much remains unknown about this process. Here, in this study, we use a combination of experiments and computational modeling to identify key aspects of the complexation phenomena which influence the self-defensive properties. We synthesize poly(acrylic acid) (PAA) microgels (∼2–5 μm diameter) via membrane emulsification and electrostatically deposit them onto polycaprolactone (PCL) coupons or onto glass to form a discontinuous submonolayer. Subsequent microgel loading with colistin or with Sub5 antimicrobial peptide (AMP) causes microgel deswelling. Under physiological conditions Sub5 remains stably sequestered whereas colistin is quickly released. Coarse-grained molecular dynamics (CGMD) simulations confirm stronger Sub5/PAA complexation. CGMD calculations also indicate that Sub5 forms dimers and higher-order structures, a prediction confirmed experimentally by Small-Angle X-ray Scattering (SAXS). Supramolecular structure entropically enhances the complexation strength because of enhanced counterion release per complexation event, and this finding can help identify other antimicrobials well suited for such a nonelutive yet self-defensive strategy. CGMD simulations also show that Sub5 has a higher complexation strength with the Staphylococcus aureus membrane than it does with PAA, confirming that there is a thermodynamic driving force for antimicrobial transfer. Such self-defensive surfaces significantly reduce S. aureus colonization (over 90% reduction relative to unmodified controls) in an in vitro hematogenous contamination model and remain cyto-compatible as evidenced by mesenchymal stem cell spreading and proliferation.
Closed-Loop Control of Active Nematic Flows
Stabilizing and shaping autonomous flows of active fluids is a fundamental challenge and a prerequisite for applications. We embed a light-responsive microtubule-based nematic in a proportional-integral control loop that adjusts the applied light intensity in response to real-time measurements of the spatially averaged flow speed. The self-regulating hardware-software-wetware system maintains a target flow speed against external or internal perturbations, including protein aging and aggregation, sample-to-sample variability, and temperature variation. Varying the controller’s gains reveals antagonistic roles between feedback and intrinsic processes, leading to nontrivial dynamics observed in fluctuation spectra. In particular, oscillations emerge from the interplay between the controller, motor binding kinetics, and active hydrodynamic relaxation. Accounting for the underlying binding timescale, our coarse-grained model and nematohydrodynamics simulations corroborate these observations. This work provides insight into the coupled dynamics of controlled active matter, laying the foundation for spatiotemporal patterning of active stress to generate and stabilize new dynamical configurations.
Challenges of conventional iterative all-atom and coarse-grained multiscale molecular dynamics
In this work, we evaluate the biomolecular dynamics behaviors when conventionally iterating between all-atom (AA) and coarse-grained (CG) molecular dynamics (MD) simulations over multiple cycles. We implemented the workflow to iterate between AA and CG in OpenMM, namely the iterative multiscale MD (iMMD) simulation workflow. In particular, we aim to identify practical applications for iterating between AA and CG simulations in a conventional manner without any constraints or model modifications. We evaluate the iMMD workflow on four representative systems, spanning folding of two soluble proteins and protein-protein as well as protein-lipid interactions of two membrane proteins. We observe that iteration between AA and CG representations could help the soluble proteins exit undesirable metastable states to fold, resulting from random protein structural distortions due to cycling. Consequently, the most reliable use of iterative AA and CG simulations appears to be to accelerating complex lipid mixing for membrane-bound protein systems rather than sampling protein conformational space. Our work explores the practical usages and limitations for iterative AA and CG simulations using readily available AA and CG force fields. The evaluated iMMD workflow in OpenMM is made available at https://github.com/lanl/iMMD.
Molecular Mechanisms Underlying Surfactant-Based Plastics De-Inking
Surfactant-mediated binder removal is critical for de-inking pretreatment in plastics recycling, yet the molecular mechanisms governing surfactant performance remain poorly understood. We used all-atom (AA) and coarse-grained (CG) molecular dynamics (MD) simulations alongside alkaline surfactant washing experiments to investigate interactions between a series of surfactants and a polyether urethane (PEU) binder in solution and on a polyethylene (PE) surface. Experiments reveal a range of de-inking efficiencies, ranging from <25 to >95% depending upon surfactant headgroup charge and tail length. AA simulations reveal that charged surfactants reach stable levels of surfactant coverage, while nonionic surfactants aggregate on the binder. CG umbrella sampling calculations quantify the thermodynamics of binder desorption in water. In ∼0.25 M surfactant solutions, up to a 52% reduction in the free energy barrier is computed, with trends in good agreement (R 2 = 0.92, Pearson’s r = –0.96, Spearman’s ρ = –0.83) with experimental de-inking efficiencies. We find that charged surfactants are more effective than nonionic surfactants for de-inking and propose three regimes of surfactant de-inking processes: good de-inking occurs in surfactants that promote PEU desorption with a low radius of gyration (R g ); moderate de-inking occurs when surfactants stabilize PEU but increase R g ; poor de-inking occurs in surfactants that aggregate on the binder and promote extensive anchoring to the surface. Overall, these molecular-level insights have the potential to guide the design of surfactant formulations for plastics recycling applications.
Correction to “COCOMO2: A Coarse-Grained Model for Interacting Folded and Disordered Proteins”
Biomolecular interactions are essential in many biological processes, including complex formation and phase separation processes. Coarse-grained computational models are especially valuable for studying such processes via simulation. Here, we present COCOMO2, an updated residue-based coarse-grained model that extends its applicability from intrinsically disordered peptides to folded proteins. This is accomplished with the introduction of a surface exposure scaling factor, which adjusts interaction strengths based on solvent accessibility, to enable the more realistic modeling of interactions involving folded domains without additional computational costs. COCOMO2 was parametrized directly with solubility and phase separation data to improve its performance on predicting concentration-dependent phase separation for a broader range of biomolecular systems compared to the original version. COCOMO2 enables new applications including the study of condensates that involve IDPs together with folded domains and the study of complex assembly processes. COCOMO2 also provides an expanded foundation for the development of multiscale approaches for modeling biomolecular interactions that span from residue-level to atomistic resolution.
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.
Structure and Flow-Viscosity of Filled-Polymer-Based 3D Printing Ink: Exploration through Coarse-Grained Molecular Dynamics
The addition of nanofiller particles to a polymer matrix has long been known to enhance or modify the composite’s mechanical and rheological properties. However, quantitatively capturing such changes with molecular level simulations remains computationally challenging. Toward that goal, we performed coarse-grained molecular dynamics of a nanocomposite system at a fixed (25 vol %) filler loading under nonspecific, weak polymer–filler interactions representative of a broad class of technologically important materials. We report several interesting results, including: (1) the equilibrium chain-configuration remains Gaussian-like as in an unfilled melt; (2) smaller filler particles display a stronger tendency to cluster; (3) larger fillers act as plasticizers by reducing the entanglement density and accelerating the chain mobility; and (4) fillers enhance the tensile response modulus, with the effect being stronger for larger particles. We also simulate cluster breakup, yielding, and elongational flow under an applied time-linear tensile strain and study the flow viscosity as a function of filler-size and chain-length.
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.
Temporal Coarse Graining for Classical Stochastic Noise in Quantum Systems
Simulations of quantum systems with Hamiltonian classical stochastic noise can be challenging when the noise exhibits temporal correlations over a multitude of time scales, such as for 1/f noise in solid-state quantum information processors. Here we present an approach for simulating Hamiltonian classical stochastic noise that performs temporal coarse-graining by effectively integrating out the high-frequency components of the noise. We focus on the case where the stochastic noise can be expressed as a sum of Ornstein-Uhlenbeck processes. Temporal coarse-graining is then achieved by conditioning the stochastic process on a coarse realization of the noise, expressing the conditioned stochastic process in terms of a sum of smooth, deterministic functions and bridge processes with boundaries fixed at zero, and performing the ensemble average over the bridge processes. For Ornstein-Uhlenbeck processes, the deterministic components capture all dependence on the coarse realization, and the stochastic bridge processes are not only independent but taken from the same distribution with correlators that can be expressed analytically, allowing the associated noise propagators to be precomputed once for all simulations. This combination of noise trajectories on a coarse time grid and ensemble averaging over bridge processes has practical advantages, such as a simple concatenation rule, that we highlight with numerical examples.
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.
Seamlessly joining length scales: From atomistic thermal graphs to anisotropic continuum conductivity
Thermal transport in complex solids is governed by local structure, defects, and anisotropy, yet most continuum models still rely on oversimplified and homogenized conductivities. Here, we bridge atomistic and continuum descriptions by building finite element (FE) models directly from the site-projected thermal conductivity (SPTC), an atomic-level decomposition of the Green–Kubo thermal conductivity. We introduce a toolkit, the “Simulator Collection for Atomic-to-Continuum Scales (SCACS)”, which uses a graph neural network to predict SPTC on large atomic structures, coarse-grains these fields into anisotropic conductivity tensors, and embeds them into the heat-flow FE equation with a customized, anisotropy-aware adaptive mesh refinement scheme. Applied to silicon nanostructures, the resulting FE models act as representative volume elements, reproduce bulk conductivities, and capture interfacial and defect-driven anisotropy while maintaining thermodynamic consistency. Additionally, SCACS predicts experimental conductance trends and fields. This work demonstrates a general route for transferring atomistic transport information into device-scale thermal simulations with physics-based approximations.
High-Order Methods in NASA’s Next Generation of Computational Fluid Dynamics Tools
The missions of the National Aeronautics and Space Administration (NASA) routinely produce unique requirements and challenges for development and application of Computational Fluid Dynamics (CFD) methods. NASA presently embodies four distinct Mission Directorates: Aeronautics Research, Exploration Systems, Science, and Space Operations. These missions generate requirements for systems that operate in a wide variety of environments. They range from the high-speed flight of aerodynamically optimized vehicles operating in the earth’s atmosphere to spacecraft designed for missions that don’t favor aerodynamic optimization, some operating in the atmosphere of planets and planetary moons such as Mars and Venus or Saturn’s moon Titan. Systems supporting these vehicles, such as rocket and jet propulsion, reaction control systems, fluid and thermal transfer systems, etc. can also generate their own unique set of flow phenomena that challenge today’s CFD methodology. Through the NASA Engineering and Safety Center (NESC), NASA annually conducts state-of-the-discipline assessments in fifteen distinct engineering disciplines. These assessments are performed by the NASA Technical Fellows that lead Technical Discipline Teams (TDT) of recognized experts in these fifteen areas. In the Aerosciences discipline, three topics have been identified as the top challenges for the discipline: aero-plume interaction prediction, unsteady separated flows, and aerothermodynamic prediction. These challenge areas are defined by the Agency’s high-risk projects and problems on which the NESC is requested to perform independent testing, analysis, and assessments. When viewed as a whole, these tests, analyses, and assessments provide a clear view of the recurring technical challenges facing Agency engineers and researchers and can be used to guide future research and technology development. The present state-of-the-art in the application of CFD at NASA is the use of Reynolds-Averaged Navier- Stokes (RANS) solvers, primarily executed in a steady-state mode of operation. In isolated cases, Unsteady RANS (URANS) solvers have been employed when steady RANS solutions produce poorly converging or oscillating results or in cases, such as aeroelastic analysis, which require unsteady aerodynamic simulation. For most traditional external and internal aerodynamic flows, structured overset grids or unstructured grids are employed to minimize geometric modeling and grid generation times. Grid adaptation, primarily as a series of coarse-grain intermediate processing steps is also seeing use on particularly complex flow problems and configurations. In the case of aerothermodynamic flows, engineers have been forced to continue to employ structured grid techniques as the present unstructured grid methodology has proven inadequate in the prediction of surface heating. In the area of aero-plume interaction modeling, two-gas, frozen chemistry simulation is generally the state-of-the- art, with some production solvers capable of predicting flows with only a single gas component. Prediction of flows falling into the afore-mentioned top Aerosciences technical challenges have severely stressed the present state-of-the-art in CFD prediction and for some problems, such as unsteady separated flows and aero-plume interaction cases, engineers have begun employing Large Eddy Simulation (LES) and Hybrid RANS/LES techniques. In some isolated aero-propulsion interaction cases, chemically reacting flow simulations have been applied. These methods are highly evolutionary and engineers have little experience in their application, so they cannot be heavily relied upon in today’s application environment. Therefore, this leads one to muse over which numerical technologies will be included in the CFD tools that will be employed 30 years in the future. This presentation will describe specific technical problems that have stressed NASA’s traditional CFD methods to their breaking point and will link these issues to the Agency’s top Aerosciences technical challenges. The discussion will then shift to the characteristics of future CFD solvers that will be required to attack these challenges and how these characteristics differ from the present state-of-the- art. High-order methods certainly appear to have a place in the development of future CFD tools and some of the physical characteristics of our most challenging problems suggest that high-order methods are the only way to effectively solve them. But there are some relatively severe implementation issues that face these methods, particularly in the area of general applicability and robust operation as an engineering tool. Desired characteristics of next-generation CFD solvers will be discussed and the author’s view of which emerging numerical technologies might be employed to address these attributes will also be presented
Development of an Amine Oxide Polyzwitterion Brush Martini Model with Polarizable Water and Ions
Abstract In order to accurately simulate the fouling process of proteins onto polyzwitterion brushes, models that accurately capture the hydration properties and chain conformations of such brushes must first be established. We developed a Martini coarse-grained (CG) model for amine oxide polyzwitterion (PNOMA) brushes, a promising class of antifouling materials, in polarizable water and ions by fitting to all-atom bond and angle distributions, monomer hydration free energy, monomer–monomer distance potential of mean force (PMF), and monomer–salt radial distribution functions (RDFs). Martini 2.2P was selected for compatibility with the established polarizable water and ion models. For comparison with PNOMA, we also constructed models for conventional sulfobetaine (PSBMA) and phosphorylcholine (PMPC) polyzwitterions and the polycation PMETAC using established nonbonded bead types from the literature and refitting bond and angle potentials. We simulated each polymer brush chemistry for varying grafting density and chain length, validating brush height scaling relations against experimental data. The CG models captured the relative hydration strengths among different polyzwitterion chemistries, and brush heights extrapolated to higher molecular weights are in agreement with experimental ellipsometry data. We find that chain swelling of the superhydrophilic PNOMA brushes lies between that of the traditional polyzwitterions PSBMA/PMPC and the polycation PMETAC. For PNOMA brushes in NaCl solution, simulated brush height decreases with salt concentration due to the selectively strong interactions between amine oxide and sodium ions.
Long-Range Dispersion Governs the Behavior of Near-Critical Fluids: Universal Scaling and Implications for Accurate Molecular Simulation
Computer simulations of near-critical and supercritical fluids often deviate from experimental results, a discrepancy commonly attributed to force-field resolutions and inaccuracies. We demonstrate that these errors primarily stem from using finite cutoff schemes for dispersion interactions, which become effectively long-range as the correlation length grows near the critical point. By employing the smooth particle-mesh Ewald method to account for full long-range dispersion, we show that coarse-grained models can also achieve high quantitative accuracy for n-undecane and propylbenzene. This approach enables the precise determination of density fluctuations and correlation lengths, which exhibit universal critical scaling and depend on the molecular size and shape. Our findings resolve a longstanding debate in molecular simulations of near-critical fluids, establishing that capturing long-wavelength fluctuations is essential for bridging the gap between the microscopic behavior and macroscopic critical phenomena.
Complex orders and chirality in the classical Kitaev-Γ model
It is well recognized that the low-energy physics of many Kitaev materials is governed by two dominant energy scales, the Ising-type Kitaev coupling 𝐾 and the symmetric off-diagonal Γ coupling. An understanding of the interplay between these two scales is therefore the natural starting point toward a quantitative description that includes subdominant perturbations that are inevitably present in real materials. This study focuses on the classical 𝐾−Γ model on the honeycomb lattice, with a specific emphasis on the region 𝐾< 0 and Γ > 0 , which is the most relevant for the available materials and which remains enigmatic in both quantum and classical limits, despite much effort. We employ large-scale Monte Carlo simulations on specially designed finite-size clusters and unravel the presence of a complex multisublattice magnetic order in a wide region of the phase diagram, whose structure is characterized in detail. We show that this order can be quantified in terms of a coarse-grained scalar-chirality order, featuring a counterrotating modulation on the two spin sublattices. Here, we also provide a comparison to previous studies and discuss the impact of quantum fluctuations on the phase diagram.