SEARCH · Search NASA

Results for “geometric programming”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Evaluation of methods for characterizing the fine particulate matter emissions from aircraft and other diffusion flame combustion aerosol sources

The U. S. Environmental Protection Agency in collaboration with the U. S. Air Force Arnold Engineering Development Complex conducted the VAriable Response In Aircraft nvPM Testing (VARIAnT) 3 and 4 test campaigns to compare nvPM emissions measurements from a variety of diffusion flame combustion aerosol sources (DFCASs), including a Cummins diesel engine, a diesel powered generator, two gas turbine start carts, a J85-GE-5 turbojet engine burning multiple fuels, and a Mini-CAST soot generator. The VARIAnT research program has conducted four test campaigns to date with the VARIAnT 3 and 4 test campaigns focusing on BC mass instrument performance, use of gas turbine engine sources for BC mass instrument calibration, and the characterization of the aerosols produced from the combustion sources. Results from both campaigns revealed agreement of about 20% between the AVL Micro-Soot Sensor, the Cavity Attenuated Phase Shift (CAPS PMSSA) monitor and the thermal-optical reference method for elemental carbon (EC), independent of the calibration source used. For the Artium LII-300, the measured mass concentrations in VARIAnT3 fall within 18% and in VARIAnT4 fall within 32% of the reference EC mass concentration when calibrated on a combustor rig in VARIAnT3 and on a LGT-60 start cart in VARIAnT4, respectively. It was also found that the three mass instrument types (MSS, CAPS PMSSA, and LII-300) can exhibit different BC to reference EC ratios depending on the emission source that appear to correlate to particle geometric mean mobility diameter, morphology, or some other parameter associated with particle geometric mean diameter (GMD) with the LII-300 showing a slightly stronger apparent trend with GMD. Systematic differences in LII-300 measured mass concentrations have been reduced by calibrating with a turbine combustion source (combustor or turbine engine). With respect to the particle size measurements made, the sizing instruments (TSI SMPS, TSI EEPS, and Cambustion DMS 500) were found to be in general agreement in terms of size distributions and concentrations with some exceptions. Teflon filter measurements of the total aerosol mass produced by the various DFCASs differed from the reference EC, BC, and integrated particle sized distribution measured aerosol masses. The Teflon filter measured mass was 30 to 40% greater than the reference EC and 50 to 60% greater than the BC mass for the J85 in both test campaigns. The measurements of particle size distributions and single particle analysis by miniSPLAT indicated the presence of large particles (>100 nm) having more compact morphologies, higher effective density, and a composition dominated by OC and ash. This increased large particle fraction is also associated with higher values of single scattering albedo measured by the CAPS PMSSA instrument and higher OC/EC measurements. These measurements indicate gas turbine engine emissions can be a more heterogeneous mix of particle types beyond the original E31 assumption that it is mainly composed of black carbon.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2024

Property Variations in Modern REBCO Coated Conductors from Multiple Manufacturers

The complex, multilayer structure of REBCO Coated Conductor (CC) poses significant challenges in the fabrication of high magnetic field devices where large stresses may initiate various forms of damage. Our goal is to peer below the cartoon representations of CC so that, amongst other things, we might better understand whether a CC from one manufacturer is interchangeable with that from another. This involves knowledge of a broad range of electromagnetic, geometric, microstructural, and J c (θ,B,T) properties, and their variations that collectively pose challenges for the fault tolerance of REBCO CC devices. Accordingly, comparative measurements of J c , visualization of flux penetration with Magneto-Optical Imaging (MOI), tape geometry from Scanning Electron Microscopy (SEM) of polished cross-sections, and extensive optical microscopy was performed on recently purchased samples from multiple manufacturers. Our analyses reveal many deviations from or characteristics absent from manufacturers’ specifications, while a comparison of different manufacturers’ mechanically and laser slit tapes shows a diverse array of slitting characteristics amongst the manufacturers and variation in properties those made to the same specification. Laser slit tapes from several manufacturers reveal ablated edges with damaged regions extending up to 50 μm, comparable to the damaged region found in the mechanically slit CC of this study. Overall, the aim of this study is to flesh out appropriate ways to understand the real conductor below the manufacturers’ cartoons to avoid surprises in our REBCO CC coil development program. The goal of this work was to perform a broad array of characterizations of the type needed for validation of purpose for making high field magnets: to our surprise we found a wide range of properties which greatly impact the mechanical strength and electromagnetic performance of solenoids composed of these conductors and reinforced the need for a broad characterization program for each conductor prior to its implementation.

Bradford, G.↗

DOE OSTI2025

Non-destructive structural characterization of graphite components using mechanical resonance and deep learning

As compared to conventional nuclear reactors, microreactors have the potential to significantly reduce construction timelines and capital costs, decreasing the barriers for advanced nuclear reactor technologies. However, the lower power output of these microreactors (typically < 20 MWe) creates challenging economics if operation and maintenance costs cannot be sufficiently reduced. The compact size of these designs presents an opportunity for comprehensive in-situ structural health monitoring to provide real-time feedback in order to reduce operational costs associated with maintenance and downtime. Many microreactor concepts use graphite for both in-core neutron moderation and as a structural material, which has typically required some form of periodic and laborious inspection. This report provides a description and assessment of recent work with graphite to couple acoustic-based experimental measurements and characterization with machine learning models to mature structural health monitoring capabilities and generate benefits for the nuclear microreactor industry. With resilient embedded sensors in development in other programs funded by the US Department of Energy’s Office of Nuclear Energy and elsewhere, the work described herein builds upon previously funded efforts to mature non-destructive testing technology that relates measured vibrational signatures to structural changes, using a combination of new experimental measurements and machine learning processing. Building on past successful demonstrations of predictive workflows to identify structural changes in a hexagonal stainless steel test article with excellent acoustic propagation, we first performed baseline characterization on graphite samples with canonical geometries to ensure compatibility and confidence in the applied techniques for a material with distinctly different mechanical properties. In contrast to efforts in previous years, we worked exclusively with unidirectional vibration data that is more comparable to those expected from the existing embedded sensor technologies which are suitable for deployment in a reactor setting. Established acoustic and modern machine-learning-based characterization approaches were applied to the resulting datasets from these simple geometries. Both approaches were found to be highly capable of detecting even small geometric irregularities amongst nominally identical samples. As such, we then moved to testing these approaches for detection of artificial local stress perturbations introduced into a more complex geometry: a hexagonal block with drilled holes. A main outcome of this work is that a generalizable ML workflow can be used to detect and predict the characteristics of small artificial anomalies in a graphite component with a relevant geometry. While this work was performed using surficial vibration data, we expect the approach to be flexible and viable for other monitoring scenarios, such as those with different arrangements or types of sensor arrays. As compared to previously funded efforts, an existing ML workflow based on neural networks was enhanced through the addition of recently developed Fourier neural operators. As applied to previously collected and new vibration datasets, prediction accuracies of anomaly characterizations were greatly improved with minimal added computational cost. As trained on small durations of vibration data (tens of seconds) collected over a realistic number of locations, the model was able to reliably determine the presence of a subtle stress anomaly and begin to provide location estimates. Such an approach is likely to be viable for more relevant reactor damage scenarios for graphite components, such as progressive crack growth or creep.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Post Irradiation Examination Results of Irradiated Yttrium Hydride

Department of Energy’s (DOE’s) Microreactor program (MRP) aims to provide the fundamental data to enable the development of microreactors. As such, material property data of critical materials for microreactor technologies are researched. Because substoichiometric yttrium dihydride (YHx, where x<2) is considered as a potential solid neutron moderator, its material property data has been combined in the Advanced Moderator Material Handbook which includes thermodynamic and thermophysical properties of YHx with the exception of irradiated material’s properties due to limited PIE. To fill the knowledge gap for the irradiated YHx, specimens and irradiation capsules were prepared at Los Alamos National Laboratory (LANL). Specimens were irradiated in the Advanced Test Reactor (ATR) at Idaho National Laboratory’s (INL’s). post-irradiation examination (PIE) was performed at INL’s Materials and Fuels Complex (MFC). This report compiles the PIE results of irradiated YHx specimens through fiscal years 2022 and 2023 (FY22-23) . The PIE data will directly be incorporated into the newer version of the Advanced Moderator Material Handbook. The main takeaways include that (i) the geometrical stability and mechanical integrity of YHx was intact with couple exceptions after high-temperature irradiations (600-800°C), (ii) hydrogen content variation due to manufacturing or irradiation in YHx caused visible surface discoloration, that is also related to the microstructural changes, (iii) qualitative comparisons of PIE methods implied that H retention was significantly higher at 600°C as compared to 800°C, (iv) thermal properties included signatures correlated with the H loss or re-gain, (v) importance of manufacturing readiness and initial as-manufactured specimens history was emphasized, (vi) the needs of targeted irradiations focusing on temperature and time parameters and very targeted PIE were specified.

08 HYDROGEN↗

DOE OSTI2023

Report on the Incorporation of Post-Irradiation Examination Results of Yttrium Hydride for Advanced Moderator Handbook

08 HYDROGEN↗

DOE OSTI2023

Phase-Field Modeling of Materials Interfaces and Nanostructures

Nanostructured materials offer unique properties for a wide range of energy applications. Among various known processing methods, liquid metal dealloying (LMD)—the selective dissolution of a base alloy element into a metallic melt—has emerged as a powerful technique to produce a new class of nano-/meso-scale open porous and bicontinuous composite structures with ultra-high interfacial area. LMD has recently been complemented by the advent of vapor phase dealloying (VPD), a novel technique that exploits the selective evaporation of one element from a parent alloy containing elements with very different vapor pressures, thereby enabling to fabricate open nanoporous structures of various elements from less-noble metals to inorganic elements regardless of their chemical activity without requiring high LMD temperatures and chemical etching. Together LMD and VPD have greatly expanded the scope of dealloying techniques, limited by traditional electrochemical means to noble metals, and boosted the design of new functional and structural materials that combine a wide variety of elements. Topologically-connected open porous structures with ultra-high surface area allow mass transport within the structure while preserving structural integrity, enabling them to serve as catalysts, fuel cells, supercapacitors, or high-capacity battery materials. Bicontinuous composite structures in turn can display high strength and high ductility or superior radiation-damage resistance due to the ultra-high interface area between interpenetrating solid phases. This research program makes use of state-of-the-art computational methods to understand at a basic level the self-organizing dealloying process with main focus on dealloying kinetics and interfacial pattern formation at the dealloying front controlling initial structure size, topology, and phase compositions. Phase-field simulation studies of LMD focus on solid solutions, line compounds, and intermetallic systems that can form ternary composites by nucleation and growth of a new phase. Studies of VPD employ phase-field modeling and a hybrid method combining a kinetic Monte Carlo (KMC) model of evaporation and surface diffusion with molecular dynamics for vapor-phase transport inside nanopores. Simulations explore mechanisms of interface- and diffusion-controlled dealloying kinetics, both observed in VPD but not fundamentally understood. In addition, our newly developed multi-physics phase-field approach of large-volume-change phase transformations is being used to model novel 3D anode geometries based on dealloyed nanoporous structures including novel sandwiched graphene/Si/silica for highrate Li ion battery. Those studies are aimed at elucidating geometric design principles that improve mechanical stability. We expect this research to enhance the capability to tailor nano-/mesoscale structures for a wide range of energy-related materials applications and to yield further advances in computational methodologies that benefit a broad materials research community.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

DOE OSTI2024

Coordinated JWST Imaging of Three Distance Indicators in a Supernova Host Galaxy and an Estimate of the Tip of the Red Giant Branch Color Dependence

Boasting a 6.5 m mirror in space, JWST can increase by several times the number of supernovae (SNe) to which a redshift-independent distance has been measured with a precision distance indicator (e.g., tip of the red giant branch (TRGB) or Cepheids); the limited number of such SN calibrators currently dominates the uncertainty budget in distance ladder Hubble constant (H 0 ) experiments. JWST/NIRCAM imaging of the Virgo Cluster galaxy NGC 4536 is used here to preview JWST program GO-1995, which aims to measure H 0 using three stellar distance indicators (Cepheids, TRGB, and J-branch asymptotic giant branch/carbon stars). Each population of distance indicator was here successfully detected—with sufficiently large number statistics, well-measured fluxes, and characteristic distributions consistent with ingoing expectations—so as to confirm that we can acquire distances from each method precise to about 0.05 mag (statistical uncertainty only). We leverage overlapping Hubble Space Telescope imaging to identify TRGB stars, crossmatch them with the JWST photometry, and present a preliminary constraint on the slope of the TRGB's F115W versus (F115W – F444W) relation equal to -0.99 ± 0.16 mag mag -1 . This slope is consistent with prior slope measurements in the similar Two Micron All-Sky Survey J band, as well as with predictions from the BaSTI isochrone suite. We use the new TRGB slope estimate to flatten the 2D TRGB feature and measure a (blinded) TRGB distance relative to a set of fiducial TRGB colors, intended to represent the absolute fiducial calibrations expected from geometric anchors such as NGC 4258 and the Magellanic Clouds. In doing so, we empirically demonstrate that the TRGB can be used as a standardizable candle at the IR wavelengths accessible with JWST.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Proximal Galerkin: A Structure-Preserving Finite Element Method for Pointwise Bound Constraints

The proximal Galerkin finite element method is a high-order, low iteration complexity, nonlinear numerical method that preserves the geometric and algebraic structure of pointwise bound constraints in infinite-dimensional function spaces. This paper introduces the proximal Galerkin method and applies it to solve free boundary problems, enforce discrete maximum principles, and develop a scalable, mesh-independent algorithm for optimal design with pointwise bound constraints. This paper also introduces the latent variable proximal point (LVPP) algorithm, from which the proximal Galerkin method derives. When analyzing the classical obstacle problem, we discover that the underlying variational inequality can be replaced by a sequence of second-order partial differential equations (PDEs) that are readily discretized and solved with, e.g., the proximal Galerkin method. Throughout this work, we arrive at several contributions that may be of independent interest. These include (1) a semilinear PDE we refer to as the entropic Poisson equation; (2) an algebraic/geometric connection between high-order positivity-preserving discretizations and certain infinite-dimensional Lie groups; and (3) a gradient-based, bound-preserving algorithm for two-field, density-based topology optimization. The complete proximal Galerkin methodology combines ideas from nonlinear programming, functional analysis, tropical algebra, and differential geometry and can potentially lead to new synergies among these areas as well as within variational and numerical analysis. Open-source implementations of our methods accompany this work to facilitate reproduction and broader adoption.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Deciphering and Manipulating Low Dimensional Magnetism

This program investigates the electronic structure and collective quantum phenomena of correlated magnetic materials using advanced angle-resolved photoemission spectroscopy (ARPES) and complementary probes, with a focus on three interrelated material families: (i) the semiconducting van der Waals magnet Cr₂Ge₂Te₆, (ii) metallic Fe-based van der Waals magnets including Fe₃GeTe₂ Fe₃GaTe₂, and Fe5GeTe₂, and (iii) kagome magnets such as FeGe, FeSn, and CsV₃Sb₅/CsCr₃Sb₅. In Cr₂Ge₂Te₆, our work established how spin excitations develop across a dimensional crossover as spins establish correlation to form long range order, providing a clean platform to isolate correlation effects. In metallic Fe-based magnets, we uncovered the dichotomy between flat and dispersive bands, revealed momentum-dependent electronic reconstructions tied to magnetic order, and demonstrated reversible, non-volatile electronic switching near room temperature, highlighting the strong coupling among spin, charge, and lattice degrees of freedom in metallic ferromagnets. In kagome magnets, we identified charge density wave formation, symmetry breaking, flat-band renormalization, and field-induced momentum-dependent electronic anisotropy, elucidating how geometric frustration and electronic correlations conspire to generate emergent quantum states. Collectively, these results establish a unified microscopic framework for understanding correlation-driven electronic reconstruction, symmetry breaking, and collective order across semiconducting and metallic magnetic systems, advancing DOE mission goals in quantum materials discovery and control.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Code Verification and Solution Verification framework in pin-resolved neutron transport code MPACT

Program verification in scientific computing encompasses the application of formal and mathematical techniques to a scientific computing code for its credibility, accuracy, and validity. Code Verification identifies bugs and performance issues in the software development stage. Solution Verification assesses the applicability of the code and the accuracy of the solution to problems of interest. Both activities utilize application cases and quantify the error against prescribed acceptance criteria. However, simply executing more application cases does not guarantee stronger or more comprehensive credibility. Here, we establish a verification framework that involves Code Verification and Solution Verification, both of which work together such that the overarching goal of “converge to the correct answer for the intended application” can be reasonably inferred. The application of such a verification framework is demonstrated using the pin-resolved neutron transport code MPACT, where standard unit tests and regression tests are covered, and where the Method of Exact Solutions and the Method of Manufactured Solutions are successfully used. Additionally, the applicability of Method of Manufactured Solutions is extended to the OECD/NEA C5G7 benchmark problems of practical material and geometric configurations. Solution Verification activities are demonstrated on a practical hierarchy of application models of increasing complexity ranging from 2D pin cell problems to 3D assembly problems. The convergence behavior and rate of convergence with respect to each individual variable are studied and provided. This framework can be adapted broadly to other fields involving scientific computing codes.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

The response field and the saddle points of quantum mechanical path integrals

Highlights: • Moyal quantum mechanics and Marinov’s path integral. • Classical and semiclassical limits of Marinov’s path integral. • Oscillating functional integrals. • Instantons of the Marinov’s path integral. In quantum statistical mechanics, Moyal’s equation governs the time evolution of Wigner functions and of more general Weyl symbols that represent the density matrix of arbitrary mixed states. A formal solution to Moyal’s equation is given by Marinov’s path integral. In this paper we demonstrate that this path integral can be regarded as the natural link between several conceptual, geometric, and dynamical issues in quantum mechanics. A unifying perspective is achieved by highlighting the pivotal role which the response field, one of the integration variables in Marinov’s integral, plays for pure states even. The discussion focuses on how the integral’s semiclassical approximation relates to its strictly classical limit; unlike for Feynman type path integrals, the latter is well defined in the Marinov case. The topics covered include a random force representation of Marinov’s integral based upon the concept of “Airy averaging”, a related discussion of positivity-violating Wigner functions describing tunneling processes, and the role of the response field in maintaining quantum coherence and enabling interference phenomena. The double slit experiment for electrons and the Bohm–Aharonov effect are analyzed as illustrative examples. Furthermore, a surprising relationship between the instantons of the Marinov path integral over an analytically continued (“Wick rotated”) response field, and the complex instantons of Feynman-type integrals is found. The latter play a prominent role in recent work towards a Picard–Lefschetz theory applicable to oscillatory path integrals and the resurgence program.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Evaluating Variable-Impedance Magnetically-Insulated Transmission Lines as a Risk-Mitigation Measure for Next-Generation Pulsed Power

This project has produced the first detailed characterizations of power flow resulting from applying the “variable-impedance MITL” concept to real-life systems in Sandia’s pulsed power program (Z and next-generation pulsed power (NGPP)). We present simulation results and analyses for constant-impedance versions of both Z and NGPP and survey the operational viability of several variable-impedance re-designs in the parameter space of linear tapers. Circuit modeling (SCREAMER/Bertha) was used to pinpoint promising candidate designs, and EM-PIC (Empire) simulations were used to evaluate these candidates more rigorously. This approach was particularly successful in the Z regime which resulted in the identification of several viable variable-impedance MITL designs for each level. The approach was more challenged in the operating space NGPP occupies, producing data points that speak to a more restrictive design space due to anode plasma turn-on. In the end, we were able to converge on one viable variable-impedance design for the highest inductance line (level “F”) and one for the highest current line (level “A”). Altogether, the body of simulation evidence presented in this report suggest there does exist flexibility in operating space for magnetically-insulated transmission lines (MITLs) having variable geometric impedance to be a potential enabling technology for safely increasing current delivery (and potentially lowering stack voltage) in pulsed-power drivers by manipulating electron losses; however, operating points for a particular design must be carefully screened. Circuit and EM-PIC modeling provided consistent verdicts in safe operating regimes for operational viability, but additional physics such as anode plasma turn-on which is included in Empire but not in SCREAMER/Bertha was found to be a critical factor affecting power flow that lead to different assessments between the codes. It is not always the case that the occurrence of anode plasma caused a design to fail (some designs turned on anode plasma yet still delivered load currents meeting design targets); the details matter such as how early in the pulse anode surfaces break down (and how large a region). However, in every case that it did fail it was found that the feedback from anode plasma was the cause (i.e., turning off the anode plasma model in Empire restored agreement with the circuit model prediction). As circuit simulations represent an efficient and practical means of surveying design space compared to more computationally-expensive approaches such as EM-PIC, it could be prudent to invest in the research and development of models to include the effects of anode plasma such as ion emission in circuit codes. The variable-impedance MITL design is a new concept that enables controlled manipulation of the initial electron losses in the outer MITL and can be tested on Z today. We encourage follow-on work to explore further optimization (including alternative variable-impedance profiles, e.g., having constant dZ/dR), and to confirm the major findings presented in this report by fielding test hardware on actual Z shots.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

CO and H 2 adsorption on Au-Ni bimetallic surfaces: a combined experimental and DFT theoretical study

Au-Ni bimetallic thin films were grown on refractory metal substrates. CO and H 2 adsorption on Au-Ni bimetallic surfaces have been studied by a combination of in situ polarization modulation infrared reflection absorption spectroscopy (PM-IRRAS), temperature-programmed desorption (TPD), and density functional theory (DFT) calculations. It is found CO desorption peak shifts from 413 K on pure Ni surfaces to 293 K on the isolated Ni atoms formed by alloying with Au atoms. The sharp decrease of CO desorption temperature on Au-Ni surfaces with increasing Au coverage is caused by the change of the favored CO adsorption sites from bridge/hollow sites on pure Ni surfaces to Ni top sites on Au-Ni bimetallic surfaces. In situ PM-IRRAS shows two CO adsorption bands on Au sites at 2119 cm −1 and 2103 cm −1 on Au-Ni surfaces at 80 K, which are due to CO bound on under-coordinated Au atoms and electron negatively charged Au sites modified with nearby Ni atoms, respectively. Even with the Au-Ni surface temperature at as low as 100 K, CO adsorption induced Ni surface segregation has been observed by in situ PM-IRRAS. Furthermore, DFT calculation results discover the adsorption energy of CO on Ni top sites continues to decrease with increasing Au coverage due to the geometric ensemble effect and the lowered d-band center after Ni alloying with Au. H 2 desorption temperature decreases from 363 K on pure Ni thin films to 302 K with increasing Au coverage to 0.6 ML. A new H 2 peak appears at around 170 K on the Au-Ni surfaces with Au coverages between 0.6 ML and 0.9 ML. This new H 2 TPD peak is assigned to H 2 desorption from the totally isolated Ni sites. With Au coverage above 1.5 ML, there is no any H 2 desorption detected. Finally, the combined surface science studies and DFT calculations provide new insights into the surface structure-activity correlation of Ni-base bimetallic surface alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Shockwave Arrival Times from Operation Redwing and Operation Upshot-Knothole

Back in the 1940s, 50s, and 60s, 210 nuclear tests were conducted in the atmosphere. As part of the testing program, EG&G performed shockwave arrival-time measurements for different shots in different operations using pressure transducers at different distances from ground zero. They then developed correlations unique to each set of data. They did not attempt to correlate all data onto one curve. As a fireball moves, material is picked up and pulled in, which is called entrainment. This extra material increases the mass of the fireball, and as the energy is spread out over a larger and larger mass, the specific internal energy of the fireball gases drops, which causes the temperature to drop as well. When the energy in the solid angle of a spherical shockwave hits the surface, it is reflected. Shortly thereafter, it catches up with the outward-moving shockwave and re-combines with it, causing it to expand at a slightly faster rate. As the shockwave grows, it approaches a hemispherical shape. The purpose of this project is to demonstrate that the shockwave arrival-time data can, in fact, be correlated onto a single function when the shockwave time-dependent solution is modified to include a geometric factor, θ, and is corrected for an entrainment coefficient.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2020

Scalable molecular dynamics on CPU and GPU architectures with NAMD

NAMD is a molecular dynamics program designed for high-performance simulations of very large biological objects on CPU- and GPU-based architectures. NAMD offers scalable performance on petascale parallel supercomputers consisting of hundreds of thousands of cores, as well as on inexpensive commodity clusters commonly found in academic environments. It is written in C++ and leans on Charm++ parallel objects for optimal performance on low-latency architectures. NAMD is a versatile, multipurpose code that gathers state-of-the-art algorithms to carry out simulations in apt thermodynamic ensembles, using the widely popular CHARMM, AMBER, OPLS, and GROMOS biomolecular force fields. Here, we review the main features of NAMD that allow both equilibrium and enhanced-sampling molecular dynamics simulations with numerical efficiency. We describe the underlying concepts utilized by NAMD and their implementation, most notably for handling long-range electrostatics; controlling the temperature, pressure, and pH; applying external potentials on tailored grids; leveraging massively parallel resources in multiple-copy simulations; and hybrid quantum-mechanical/molecular-mechanical descriptions. We detail the variety of options offered by NAMD for enhanced-sampling simulations aimed at determining free-energy differences of either alchemical or geometrical transformations and outline their applicability to specific problems. Last, we discuss the roadmap for the development of NAMD and our current efforts toward achieving optimal performance on GPU-based architectures, for pushing back the limitations that have prevented biologically realistic billion-atom objects to be fruitfully simulated, and for making large-scale simulations less expensive and easier to set up, run, and analyze. NAMD is distributed free of charge with its source code at www.ks.uiuc.edu.

high-performance computing↗

DOE OSTI2024

Geometry of Non-Diffusive Tracer Transport in Gridded Atmospheric Models

A first-order linear and numerically non-diffusive Eulerian transport algorithm, minVAR, was recently developed for preservation of correlations between interrelated tracers during advective transport. The present study extends this work by: (1) providing further investigation of several interesting geometric constructions found in contours of constant minVAR, short for minimum variance, through extension to three coordinate dimensions. These contours capture point-by-point representations of thousands of individual atmospheric aerosol and/or cloud particles as they evolve and are rendered on Eulerian grids at a level of sub-grid resolution limited only by numerical precision; and (2) exploration of geometric similarities between the Arakawa C-grid, used to obtain interpolated values of the wind field at grid scale and minVAR. In particular, we consider interpolation of the u and v horizontal components of wind velocity from grid to sub-grid scales. The last results are motivated by recent applications of the Weather Research and Forecasting (WRF) model applied in the coastal Houston region, where the recent TRacking Aerosol Convection Interactions ExpeRiment (TRACER) field campaign was organized. A unique and fully consistent mapping is obtained between particles moving along meteorological wind trajectories and the non-diffusive, non-dispersive representation of such trajectories on an Eulerian grid.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Extend an innovative HPC-Compatible Multiple Temporal-spatial Resolution Concurrent Finite Element Modeling Approach to Guide Laser Powder Bed Fusion Additive

Laser power bed fusing (PBF) additive manufacturing is a key enabling technology to manufacture highly complex and integrated automotive structures. However, the geometric complexity of PBF-AM technique also leads to highly non-uniform heating and cooling rate in the manufactured part, which may cause flaw formation and produce excessive and nonuniform residual stresses, which increase quality uncertainties and manufacture issues, leading to increases in cost and energy consumption in the form of rejected parts. In this research project, we developed an innovative Multi-Spatial-Temporal-Resolution Finite Element (MUST-FE) method and completed the corresponding high performance computation (HPC) platform-based in-house code, which enables high accuracy prediction of temperature and residual stress fields for component-scale PBF-AM manufacture in efficient computation time. The MUST-FE model is calibrated and validated with a “2D pad” AlSi10Mg experiments by matching the melt pool shape and dimension, and with a “XY-cross” AlSi10Mg experiment by matching the thermal distortion and residual stress. The innovative multi-resolution and concurrent modeling approach adopted in this code ensures accuracy and computational efficiency, which will enable energy-efficient and high-yield, low-cost manufacturing of optimized, qualifiable automotive structures and contribute towards reaching technical targets outlined in AMO’s Program Plan to develop additive manufacturing systems that deliver consistently reliable parts with predictable properties.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Refining perovskite structures to pair distribution function data using collective Glazer modes as a basis

Structural modelling of octahedral tilts in perovskites is typically carried out using the symmetry constraints of the resulting space group. In most cases, this introduces more degrees of freedom than those strictly necessary to describe only the octahedral tilts. It can therefore be a challenge to disentangle the octahedral tilts from other structural distortions such as cation displacements and octahedral distortions. This paper reports the development of constraints for modelling pure octahedral tilts and implementation of the constraints in diffpy-CMI, a powerful package to analyse pair distribution function (PDF) data. The model in the program allows features in the PDF that come from rigid tilts to be separated from non-rigid relaxations, providing an intuitive picture of the tilting. The model has many fewer refinable variables than the unconstrained space group fits and provides robust and stable refinements of the tilt components. It further demonstrates the use of the model on the canonical tilted perovskite CaTiO 3 which has the known Glazer tilt system α + β – β – . The Glazer model fits comparably to the corresponding space-group model Pnma below r = 14 Å and becomes progressively worse than the space-group model at higher r due to non-rigid distortions in the real material.

36 MATERIALS SCIENCE↗