Search NASA⌕ Search

SEARCH · Search NASA

Results for “force matrix”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Degradation science: Integrating modeling and experiments to predict localized corrosion processes (Annual Progress Report)

Additively manufactured (AM) eutectic high-entropy alloys (EHEAs), such as nano-lamellar AlCoCrFeNi 2.1 , have excellent strength, ductility, and wear resistance even at elevated temperatures, but their corrosion behavior in aggressive acids at different length scales remain poorly understood. This work investigates the corrosion behavior of laser powder bed–fused (L-PBF) AlCoCrFeNi 2.1 as a function of annealing temperatures, probing degradation mechanisms from nanoscopic to macroscopic length scales. The alloy is dual phase consisting of a ductile FCC L1 2 phase and a high-strength BCC B2 phase. Rapid solidification during L-PBF produces a far-from-equilibrium nano-lamellar structure with nearly homogeneous elemental distribution, which tends to evolve upon annealing toward Cr/Co/Fe-enriched FCC and Al/Ni-enriched B2. Three conditions were studied: as-printed, 600 °C/5 h, and 1000 °C/1 h, over which B2 lamellae coarsen, lamellar spacing increases, and elemental segregation becomes more prominent. Microstructure and chemistry were characterized by scanning electron microscopy (SEM) and energy-dispersive spectroscopy (EDS), while in-situ electrochemical atomic force microscopy (EC-AFM) was used to link early (<5 h) local dissolution to microstructure after exposure in sulfuric acid. EC-AFM highlights preferential dissolution of the BCC/B2 phase where the surrounding matrix is Cr-depleted and directly quantifies the dissolution rates within individual phases, tracks the transition from early nano-scale attack to partial repassivation, to correlate height differences with current and impedance responses. To monitor longer-term behavior (up to 96 h), ex-situ AFM, SEM, and confocal imaging were combined with conventional bulk electrochemical tests, bridging nanoscale observations to micro/meso-scale damage morphologies. At the meso-scale, the deepest dissolution channels align with the build-direction lamellae and melt-pool boundaries, indicating that printing directionality guides the propagation of these localized corrosion sites. Annealing modifies corrosion by restructuring BCC/FCC phase fractions, lamellar spacing, and Cr/Al segregation, thereby changing the cathode/anode ratio and passive film stability. The results clarify how as-printed and annealed nano-lamellar architectures differ in their susceptibility to selective dissolution; how elemental segregation competes with residual stresses along the build direction. With these insights, future work will use CALPHAD-guided alloy modification to stabilize higher Cr contents in the B2 phase while retaining a dominant FCC+B2/BCC microstructure, with the goal of designing mechanically robust, corrosion-resistant EHEAs for safety-critical applications to leverage the LLNL’s broader national and global security mission.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Forced Oscillation Source Location of Bulk Power Systems Using Synchrosqueezing Wavelet Transform

Forced oscillation source location (FOSL) plays a significant role in mitigating forced oscillations (FOs), which threaten the power system stability. Here, this paper proposes a data-driven approach for FOSL in power systems using synchrosqueezing wavelet transform (SWT). The proposed approach conducts SWT on the measured system responses to obtain the SWT matrix. Then, the SWT-based dissipating energy flow (DEF) model in time-frequency domain and dissipating energy spectrum (DES) model in frequency domain are derived from the traditional DEF model. Further, the characteristics of SWT-based DEF and DES are revealed by referring to the traditional DEF, and the FOSL criteria of the SWT-based DEF and DES can be hereby obtained. Using the obtained FOSL criteria, the FO source can be located from the measured responses. The performance of the proposed FOSL method is evaluated by simulation data of the WECC 179-bus test system and field-measurement PMU data of the ISO New England. The results confirm the accuracy and efficiency of the proposed method in the FOSL.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Structure optimization with stochastic density functional theory

Linear-scaling techniques for Kohn–Sham density functional theory are essential to describe the ground state properties of extended systems. Still, these techniques often rely on the localization of the density matrix or accurate embedding approaches, limiting their applicability. In contrast, stochastic density functional theory (sDFT) achieves linear- and sub-linear scaling by statistically sampling the ground state density without relying on embedding or imposing localization. In return, ground state observables, such as the forces on the nuclei, fluctuate in sDFT, making optimizing the nuclear structure a highly non-trivial problem. In this work, we combine the most recent noise-reduction schemes for sDFT with stochastic optimization algorithms to perform structure optimization within sDFT. We compare the performance of the stochastic gradient descent approach and its variations (stochastic gradient descent with momentum) with stochastic optimization techniques that rely on the Hessian, such as the stochastic Broyden–Fletcher–Goldfarb–Shanno algorithm. In conclusion, we further provide a detailed assessment of the computational efficiency and its dependence on the optimization parameters of each method for determining the ground state structure of bulk silicon with varying supercell dimensions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

A family of tubular pili from harmful algal bloom forming cyanobacterium Microcystis aeruginosa

Cyanobacteria are vital photosynthetic prokaryotes, but some form harmful algal blooms (cyanoHABs) that disrupt ecosystems and produce toxins. The mechanisms by which these blooms form have yet to be fully understood, particularly the role of extracellular components. Here, we present a 2.4 Å cryo-EM structure of a pilus, termed the cyanobacterial tubular (CT) pilus, found in the cyanoHAB-forming Microcystis aeruginosa. The pilin exhibits a unique protein fold, forming a tubular pilus structure with tight, double-layer anti-parallel β-sheet interactions. We show that CT pili are essential for buoyancy by facilitating the formation of micro-colonies, which increases drag force and prevents sinking. The CT pilus surface is heavily glycosylated with ten monosaccharide modifications per pilin. Furthermore, CT pili can enrich microcystin, potentially enhancing cellular resilience, and co-localize with iron-enriched extracellular matrix components. Thus, we propose that this pilus plays an important role in the proliferation of cyanoHABs. This just discovered pilus family appears to be widely distributed across several cyanobacterial orders. Our structural and functional characterization of CT pili provide insights into cyanobacterial cell morphology, physiology, and toxin interactions, and identify potential targets for disrupting bloom formation.

Cryoelectron microscopy↗
DOE OSTI2025

A Poromechanical Framework for Internal Interactions Induced by Solid Inclusions

The framework of poromechanics is generalized to simulate the multiscale behavior of porous media subjected to internal loadings stemming from the growth of solid inclusions. This generalization is designed to enable the study of anisotropic internal stress generation from solid growth within the pores, while recovering isotropic fluid-induced loading as a particular case. For this purpose, a mathematical strategy to define constitutive tensors in a thermodynamically consistent form is proposed, thus offering new opportunities for determining the poromechanical properties of a porous solid through advanced experimentation or micromechanical models. The framework is specialized by means of established elastic solutions for single pore–matrix interaction, as well as through homogenization schemes considering the interaction among congruent pores. In particular, the second Eshelby solution and the Tanaka–Mori–Benveniste homogenization scheme are used to derive a microporoelastic model. At an elemental scale, the model is tested under mixed control conditions by replicating different scenarios of geomaterial testing. In addition, the model characteristics are outlined with reference to inelastic microscopic loadings replicating chemo-mechanical forcing, such as expansive crystal formation. Through a series of parametric analyses, it is shown that the microstructure of the pores significantly influences the properties of porous media. Most notably, it is shown that the effects of a solid forming within the pores depend in a highly nonlinear fashion on the constitutive characteristics of the inhomogeneities and can therefore not be readily quantified or predicted without models capturing the diverse multiscale interactions among pores, inhomogeneities, and matrix.

chemomechanics↗
DOE OSTI2020

How equation of state selection impacts accuracy near the critical point: Forced convection supercritical CO 2 flow over a cylinder

Increasing application of supercritical fluids leads to increased use of numerical simulations to guide designs. However, the impact of Equations of State (EoS) selection on computational cost and accuracy is not well understood. Laminar flow around a heated cylinder is used as a test case. A robust matrix of 99-simulations is used to analyze the flow of carbon dioxide just above the critical point at the isobaric condition of 8 MPa and temperatures between 305 and 390 K. Four EoS were studied: Ideal Gas Law, cubic EoS’s Peng-Robinson and Soave-Redlich-Kwong, and the Span-Wagner EoS. Results were identical for isothermal cases, but for heated or cooled cases, the high-fidelity Span-Wagner EoS simulations presented drastically different wake structures. Compared to the high-fidelity SW EoS, cEoS’s underestimated Nu values by up to 38% and cEoS’s underestimated CD values by 18–42%. In conclusion, the Span-Wagner EoS required only 17% more computation time then the cubic EoS simulations, making the Span-Wagner EoS recommended for engineering design.

42 ENGINEERING↗
DOE OSTI2026

Single-particle detection of enhanced polarizability in Au-decorated semiconducting nanorods via scanning dielectric microscopy

Hybrid nanostructures that combine semiconducting and metallic components offer great potential for photothermal therapy, optoelectronics, and sensing, by integrating tunable optical properties with enhanced light absorption and charge transport. Boosting the integrated performance of these hybrid systems demands techniques capable of probing local variations of the physical properties inaccessible to bulk analysis. Here, we report the single-particle dielectric characterization of hybrid, semiconducting bismuth sulfide (Bi 2 S 3 ) nanorods (NR) decorated with metallic Au nanoparticles (NP), employing scanning dielectric microscopy, which uses electrostatic force microscopy in combination with finite-element numerical simulations. We reveal a pronounced enhancement in the local dielectric response of Bi2S3 upon Au decoration, attributed to interfacial polarization and electron transfer from Au to the Bi 2 S 3 matrix, thus suggesting a enhanced metallic-like polarizability at the single-particle level. Numerical simulations show that the response is dominated by the vertical component of the permittivity and that the decorating metallic Au NP produce only moderate shielding of the semiconductor Bi 2 S 3 NR core, indicating that the large increase in the dielectric response originates primarily from intrinsic modifications within the NR. Overall, these findings provide direct insight into structure–property relationships at the single-particle level, supporting the rational design of advanced hybrid nanostructures with tailored electronic functionalities.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Long short-term memory embedded nudging schemes for nonlinear data assimilation of geophysical flows

Reduced rank nonlinear filters are increasingly utilized in data assimilation of geophysical flows, but often require a set of ensemble forward simulations to estimate forecast covariance. On the other hand, predictor-corrector type nudging approaches are still attractive due to their simplicity of implementation when more complex methods need to be avoided. However, optimal estimate of nudging gain matrix might be cumbersome. In this paper, we put forth a fully nonintrusive recurrent neural network approach based on a long short-term memory (LSTM) embedding architecture to estimate the nudging term, which plays a role not only to force the state trajectories to the observations but also acts as a stabilizer. Furthermore, our approach relies on the power of archival data and the trained model can be retrained effectively due to power of transfer learning in any neural network applications. In order to verify the feasibility of the proposed approach, we perform twin experiments using Lorenz 96 system. Our results demonstrate that the proposed LSTM nudging approach yields more accurate estimates than both extended Kalman filter (EKF) and ensemble Kalman filter (EnKF) when only sparse observations are available. With the availability of emerging AI-friendly and modular hardware technologies and heterogeneous computing platforms, we articulate that our simplistic nudging framework turns out to be computationally more efficient than either the EKF or EnKF approaches.

42 ENGINEERING↗
DOE OSTI2026

Rotation of Crystal Seed During Early Stages of Growth Reveals the Anisotropy of Glass Matrix

Glass-ceramics are a highly valuable class of materials that are needed for realizing many high-tech as well as consumer products. They are made by controlling crystallization of glass. Theoretical atomistic insights obtained using molecular dynamics simulations in this work reveal that crystal seeds can spontaneously rotate within glass at high temperatures due to anisotropic forces at the glass-crystal interface. This rotation challenges the long‐held view of glass as isotropic solid, and highlights how atomic‐scale inhomogeneities influence crystal growth and microstructure, hence the physical properties glass‐ceramics.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Power capture and power take-off load of a self-balanced dual-flap oscillating surge wave energy converter

Wave energy converters are an important part of future renewable energy infrastructure. Predicting their power matrix, capture width ratio, and power take-off loads at a targeted site is required for performance assessment before deployment. Because their testing is very expensive, numerical modeling and simulations play a significant role in those assessments. Linear potential flow theory has limited accuracy under large amplitude wave forcing. More accurate predictions can be obtained by using higher-fidelity models, which are computationally expensive. We present a framework for multi-fidelity numerical simulations to determine the hydrodynamic response, wave capture capability, and power take-off load of a full-scale dual-flap oscillating surge wave energy converter. This design exploits out-of-phase motion by setting the distance between the flaps to half the wavelength of the most occurring wave. The simulations are validated using a 1:10 model experiments in a wave tank. Based on these validations, it was determined that Euler simulations provide an acceptable prediction with 90% reduction in computational time with only 11% error. Utilizing Euler simulations at full-scale, the results demonstrate that the annual electrical energy output is 1.79 GWh under regular wave conditions. Here, one significant improvement over single-flap designs is the capture width ratio which exceeds unity.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Elastic Antiproton-Nucleus Scattering from Chiral Forces

Elastic scattering of antiprotons off He 4 , C 12 , and O 16 , 18 is described for the first time with a consistent microscopic approach based on the calculation of an optical potential (OP) describing the antiproton-target interaction. The OP is derived using the recent antiproton-nucleon ( p ¯ N ) chiral interaction to calculate the p ¯ N t matrix, while the target densities are computed with the ab initio no-core shell model using chiral interactions as well. Here, our results are in good agreement with the existing experimental data and the results computed at different chiral orders of the p ¯ N interaction display a well-defined convergence pattern.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

On the melt pool dynamic of voxel-controlled metal matrix composites via hybrid additive manufacturing: Laser powder bed fusion and ink-jetting

In this study, the effect of the addition of reinforcement nanoparticles to the 316L matrix by adopting ex-situ and in-situ method (drop on demand jetting) to produce 316L/Al 2 O 3 nanocomposite was investigated. In the ex-situ method, the Al 2 O 3 nanoparticles (NPs) were lightly mixed with 316L powder and processed by laser powder bed fusion. In the in-situ method, an ethanol-based ink containing Al 13 nanoclusters (NCs) was added to 316L powder and then processed by laser. Both ex-situ and in-situ method produced nanocomposites with Al-Si-Mn-O-enriched precipitations within the 316L matrix. The addition of NPs/NCs to the 316L matrix, altered the geometrical characteristic of the single-track melt pools. At the same laser power, with increasing the amount of Al 2 O 3 NPs and Al 13 NCs the melt pool deepened due to reduced thermal conductivity and prolonged liquid presence. Further, as a result, 316L/1 wt% Al 13 NCs deposited single track showed larger grains in comparison to 316L single track. At a high laser power of 150W, the Marangoni flow and the buoyancy force caused the nanoparticles to agglomerate and float to the top surface of tracks; therefore, the wt% fraction of precipitation was drastically reduced due to the loss of Al. The 316L/Al 2 O 3 NPs and 316L/Al 13 NCs exhibited the microhardness of 285 ± 13 HV and 293 ± 7 HV, respectively, higher than the deposited 316L single track, 265 ± 15 HV. Lastly, a hybrid LPBF+ink-jet printer was adopted to selectively change the composition of different zones by adding Al 13 NCs ink to 316L and producing a voxel-controlled metal matrix composite.

316L↗
DOE OSTI2021

Thermodynamic properties of Gd-Bi alloys determined by emf measurements in LiCl-KCl-GdCl 3 electrolyte

Thermodynamic properties of binary Gd-Bi alloys (mole fraction, x Gd = 0.02–0.35) were determined using electromotive force (emf) measurements in molten LiCl-KCl-GdCl 3 electrolyte in complement with structural, microstructural, and thermal characterization. X-ray diffractometry (XRD) and scanning electron microscopy (SEM) identified the phase constituents of Gd-Bi alloys (x Gd = 0.02–0.35) as Bi matrix and GdBi compound. According to differential scanning calorimetry (DSC) measurements up to 1250 K, a eutectic transition [L = Bi + GdBi] was observed at 543 K but the reported peritectic transition [L + GdBi = GdBi 2 ] was not detected, suggesting the absence of the GdBi 2 compound. Based on structural, microstructural, and thermal characterization of Gd-Bi alloys (x Gd = 0.02–0.35), the GdBi 2 compound was not detected, and its thermodynamic stability is questionable. The emf of Gd-Bi alloy (x Gd = 0.16) at 700–1048 K relative to Gd(s) was measured by electrodepositing pure Gd metal at 25 K increments. Further, the emf values of Gd-Bi alloys were determined via coulometric titration of Gd into liquid Bi using a two-phase [L + GdBi] alloy as a reference electrode at 773–973 K. The solubility of Gd in liquid Bi was also estimated to be 0.50 mol% (773 K), 0.99 mol% (873 K), and 2.21 mol% (973 K) and the excess partial molar Gibbs energy (Δ$\mathrm{\bar{G}^{ex}_{Gd}}$) of liquid Gd-Bi alloys was as low as 65 kJ mol –1 , indicating strong chemical interactions between Gd and Bi.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Causality and dimensionality in geometric scattering

We report the scattering matrix which describes low-energy, non-relativistic scattering of spin-1/2 fermions interacting via finite-range potentials can be obtained from a geometric action principle in which space and time do not appear explicitly. In the case of zero-range forces, causality leads to constraints on scattering trajectories in the geometric picture. The effect of spatial dimensionality is also investigated by considering scattering in two and three dimensions. In the geometric formulation it is found that dimensionality is encoded in the phase of the harmonic potential that appears in the geometric action.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Two-dimensional patterning of mesoscale fibers using acoustophoresis

The performance of functional composites can rely critically on the arrangement of secondary phases; for example, patterned networks of conductive particles can impart anisotropic thermal, electric or ionic conductivity while preserving flexibility in the matrix. We demonstrate the use of standing acoustic waves to generate periodic patterns of short fibers. We extend the range of possible patterns with the first demonstration of both rectangular grids and arrays of octagons interspersed with rectangles. These newly demonstrated patterns are rationalized using theoretical models of acoustic forces and torques on fibers that account for two-dimensional spatial variations arising from applied acoustic fields. The models enable simulations of fiber motion, which are used to (i) map out final fiber positions as a function of initial position and orientation, and (ii) corroborate experiments visualizing fiber motion and final patterns. This approach provides a fast and accurate way to predict emergent fiber patterns as a function of excitation modality and fiber length. The theory and experiments clearly indicate strong coupling between the length of the fibers and the spacing of the acoustic nodes. Furthermore, this coupling is used to estimate reductions in percolation thresholds associated with the ratio of fiber length and acoustic wavelength.

25 ENERGY STORAGE↗
DOE OSTI2022

Colossal piezoresistance in narrow-gap Eu 5 In 2 Sb 6

Piezoresistance, the change of the electrical resistance (R) of a material in response to an applied mechanical stress (σ), is the driving principle of electromechanical devices such as strain gauges, accelerometers, and cantilever force sensors. Enhanced piezoresistance has been traditionally observed in two classes of uncorrelated materials: nonmagnetic semiconductors and composite structures. In this study, we report the discovery of a remarkably large piezoresistance in Eu 5 In 2 Sb 6 single crystals, wherein anisotropic metallic clusters naturally form within a semiconducting matrix due to electronic interactions. Eu 5 In 2 Sb 6 shows a highly anisotropic piezoresistance, and uniaxial pressure along [001] of only 0.4 GPa leads to a resistivity drop of >99.95%, which results in a colossal piezoresistance factor of 5000×10 –11 Pa –1 . Our result not only reveals the role of interactions and phase separation in the realization of colossal piezoresistance, but it also highlights a route to multifunctional devices with large responses to both pressure and magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Seismic Bracing Development for the Moore 3 Diamond Turning Machine at Lawrence Livermore National Laboratory

To allow for on-site Potassium Dihydrogen Phosphate (KDP) crystal finishing capabilities, Lawrence Livermore National Laboratory (LLNL) acquired the Moore 3 diamond turning machine. Equipment such as the Moore 3 is required by LLNL to have adequate seismic bracing before permanent installation to reduce the risk of personnel injury or machine damage during an earthquake. The purpose of this project was to design a bracket that can withstand seismic forces induced by the machine without hindering function. To develop a solution, background research—including stakeholder identification, stakeholder interviews, and a literature survey—was conducted. Based on LLNL site and stakeholder requirements, several approaches were developed for bracing the Moore 3. Ideations were evaluated based on a weighted decision matrix and a final design strategy was selected. This design was an L-shaped weldment that featured provisions for four HILTI concrete anchors, three triangular gussets, and a steel construction.

42 ENGINEERING↗
DOE OSTI2021

Quantum-Based Molecular Dynamics Simulations Using Tensor Cores

Tensor cores, along with tensor processing units, represent a new form of hardware acceleration specifically designed for deep neural network calculations in artificial intelligence applications. Tensor cores provide extraordinary computational speed and energy efficiency but with the caveat that they were designed for tensor contractions (matrix–matrix multiplications) using only low-precision floating-point operations. Despite this perceived limitation, we demonstrate how tensor cores can be applied with high efficiency to the challenging and numerically sensitive problem of quantum-based Born–Oppenheimer molecular dynamics, which requires highly accurate electronic structure optimizations and conservative force evaluations. The interatomic forces are calculated on-the-fly from an electronic structure that is obtained from a generalized deep neural network, where the computational structure naturally takes advantage of the exceptional processing power of the tensor cores and allows for high performance in excess of 100 Tflops on a single Nvidia A100 GPU. Stable molecular dynamics trajectories are generated using the framework of extended Lagrangian Born–Oppenheimer molecular dynamics, which combines computational efficiency with long-term stability, even when using approximate charge relaxations and force evaluations that are limited in accuracy by the numerically noisy conditions caused by the low-precision tensor core floating-point operations. A canonical ensemble simulation scheme is also presented, where the additional numerical noise in the calculated forces is absorbed into a Langevin-like dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗