Search NASA⌕ Search

SEARCH · Search NASA

Results for “Methods”

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

Molten Halide Salt Surface Tension: Methods and Correlations

Here, this paper reviews various methods for studying surface tension and their applicability to fluoride and chloride molten salt systems, including a comparison of benefits and drawbacks. Such a comparison aids in experiment design based on desired factors such as scale, accuracy, and repeatability. A detailed review is presented for existing literature data regarding the surface tension of molten fluoride and chloride salts. These reference data were compiled and analyzed to determine cross-validated correlation equations for several alkali and alkaline earth fluoride and chloride salts as functions of temperature. These correlations are necessary for reliable multiphysics modeling approaches as well as accurate design and analysis of multiphase molten salt phenomena such as gas sparging and bubble formation/transport. This analysis supports the development of the thermophysical arm of the Molten Salt Thermal Properties Database (MSTDB-TP) managed by Oak Ridge National Laboratory.

Chloride↗
DOE OSTI2026

A Stochastic Quasi-Newton Method in the Absence of Common Random Numbers

We present Q-SASS, a quasi-Newton method for unconstrained stochastic optimization that does not rely on common random numbers. Most existing quasi-Newton approaches leverage common random numbers to construct second-order updates. However, motivated by challenges in variational quantum algorithms—where such coordination is not possible—we consider the setting in which function values and gradients are accessible only through noisy probabilistic zeroth- and first-order oracles, and no common random numbers can be exploited. We derive high-probability tail bounds on the iteration complexity of our algorithm for nonconvex, convex, and strongly convex (more generally, those satisfying the PL condition) objective functions. Finally, we demonstrate the empirical benefits of our quasi-Newton updating scheme on both synthetic and quantum chemistry problems.

Complexity bound↗
DOE OSTI2026

Adaptive Sampling-Based Bi-Fidelity Stochastic Trust Region Method for Stochastic Derivative-Free Optimization

Bi-fidelity stochastic optimization has gained increasing attention as an efficient approach to reduce computational costs by leveraging a low-fidelity (LF) model to optimize an expensive high-fidelity (HF) objective. In this paper, we propose ASTRO-BFDF, an adaptive sampling trust-region method specifically designed for unconstrained bi-fidelity stochastic derivative-free optimization problems. In ASTRO-BFDF, the LF function serves two purposes: (i) to identify better iterates for the HF function when the optimization process indicates a high correlation between them and (ii) to reduce the variance of the HF function estimates using bi-fidelity Monte Carlo (BFMC). The algorithm dynamically determines sample sizes while adaptively choosing between crude Monte Carlo and BFMC to balance the trade-off between optimization and sampling errors. We prove that the iterates generated by ASTRO-BFDF converge to a first-order stationary point almost surely. Additionally, we demonstrate the effectiveness of the proposed algorithm through numerical experiments on synthetic benchmarks and simulation optimization problems involving discrete event systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Validation of a Tapered Impact Damper for Traffic Signal Structure Vibrations Using the Method of Harmonic Balance

Significant vibration amplitudes and cycles can be produced when traffic signal structures with low inherent damping are excited near one of their natural frequencies. For the mitigation of wind-induced vibrations, dynamic vibration absorbers coupled to the structure are often used. Here, this research investigates the performance of a tapered impact damper, consisting of a hanging spring-mass oscillator inside a housing capable of reducing vibration amplitude over a broader frequency range than the conventional tuned mass damper. A nonlinear, two degree-of-freedom model is developed with coordinates representing the traffic structure and the tapered impact damper. This research focuses on the application of the harmonic balance method to approximate the periodic solutions of the nonlinear equations to compute the nonlinear dynamics of the damped traffic signal structure. After designing and manufacturing a tapered impact damper, the traffic signal structure is tested with and without the damper using free vibration snapback tests. The experimental frequency and damping backbone curves are used to validate the analytical model, and the effectiveness of the damper is discussed.

42 ENGINEERING↗
DOE OSTI2024

Modeling powder spreadability in powder-based processes using the discrete element method

Powder-bed fusion (PBF) processes refer to a subset of Additive Manufacturing (AM) techniques where powder is spread on the build-plate before melting (by a laser or electron beam). While PBF processes are attractive due to their ability for realizing complex structures that are either difficult or impossible to create through conventional means, the parts fabricated with these techniques can exhibit defects such as pores, inclusions, and excessive surface roughness. To minimize these defects, much research has been dedicated towards process maturation by optimizing laser or electron beam parameters. However, these developmental efforts typically do not address the recoating process where achieving dense and uniform layers of powder is a necessity for ensuring process repeatability and part quality. While the recoating process can be studied through experimentation, the dynamics of particle movement are difficult to analyze experimentally. Therefore, here, in this study, powder spreading in PBF was simulated through the Discrete Element Method (DEM) to elucidate the mechanisms that control powder-bed quality. Utilizing the Buckingham Pi theorem, a dimensionless metric referred to as the spreading index is developed that combines powder-bed density, roughness, and particle size to assess the quality of powder layers. The formulated spreading index is then related to several dimensionless quantities that provide insight into the mechanisms dominating powder spreading in PBF. The DEM simulations conducted in this work focused on the scenario where powder is spread onto an existing powder bed and revealed that a reduction in the recoating velocity causes an increase in the spreading index while little to no impact on the spreading index was observed when varying layer thickness from 30 μm to 75 μm.Particle size effects on the powder-bed quality were also investigated.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The rigorous upscaling of advection-dominated transport in heterogeneous porous media via the Method of Finite Averages

Systems involving advection-dominated transport through heterogeneous porous and fractured media are ubiquitous in subsurface engineering applications. However, upscaling such systems continues to challenge rigorous modeling efforts, particularly when advection is stronger than diffusion at fine spatial scales (i.e., when the Péclet number is greater than one at length scales that characterize a system’s unit-cells, representative elementary volumes, or averaging regions). Here, in this work, we propose and validate a strategy for extending the Method of Finite Averages (MoFA), a rigorous upscaling methodology for heterogeneous porous media, to upscale transport systems experiencing stronger advection than diffusion at fine scales (i.e., fine-scale Péclet numbers greater than one). We detail the strategy, the physical conditions under which it can be applied while retaining a priori modeling error guarantees, and implement the strategy to obtain a MoFA model for advective-diffusive transport that accommodates advective physics at fine spatial scales. We then perform two numerical experiments considering systems with system-scale Péclet numbers of 300 and 1000 — which correspond to fine-scale Péclet numbers of 30 and 100, respectively — to verify that the error guarantees are met under the strategy. After, we conduct a numerical study to demonstrate the strategy’s advantages over the original MoFA methodology. The results suggest that rigorously-upscaled transport models for heterogeneous porous media experiencing advective physics at finer spatial scales can be derived through MoFA and resolved orders of magnitude faster than their pore-scale counterparts. The results also suggest that the presented strategy is limited to modeling shallow concentration gradients when there are large differences between the time scales related to advection and a system’s temporally-varying boundary conditions. This limitation hinders the strategy’s practicality in modeling more advective systems, and as such, opportunity exists for developing additional strategies that accommodate rapidly-varying boundary conditions — and consequentially, steeper concentration gradients — while modeling advective systems with MoFA.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Plastic additives in the ocean: Use of a comprehensive dataset for meta-analysis and method development

In excess of 13,000 chemicals are added to plastics (‘additives’) to improve performance, durability, and production of plastic products. They are categorized into numerous chemical classes including flame retardants, light stabilizers, antioxidants, and plasticizers. While research on plastic additives in the marine environment has increased over the past decade, there is a lack of methodological standardization. To direct future measurement of plastic additives, we compiled a first-of-its-kind dataset of literature assessing plastic additives in marine environments, delineated by sample type (plastic debris, seawater, sediment, biota). Using this dataset, we performed a meta-analysis to summarize the state of the science. Currently, our dataset includes 217 publications published between 1978 and May 2023. The majority of publications analyzed plastic additives in biota collected from Europe and Asia. Analyses concentrated on plasticizers, brominated flame retardants, and bisphenols. Common sample preparation techniques included Solvent - Agitation extraction for plastic, sediment, and biota samples, and Solid Phase Extraction for seawater samples with dichloromethane and solvent mixtures including dichloromethane as the organic extraction solvent. Finally, most analyses were performed utilizing gas chromatography/mass spectrometry. There are a variety of data gaps illuminated by this meta-analysis, most notably the small number of compounds that have been targeted for detection compared to the large number of additives used in plastic production. The provided dataset facilitates future investigation of trends in plastic additive concentration data in the marine environment (allowing for comparison to toxicity thresholds) and acts as a starting point for optimizing and harmonizing plastic additive analytical methods.

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

An immersed interface method for microstructure-scale electrochemical battery models: numerical formulation and performance portable implementation

We present the numerical formulation, verification, and performance portable implementation of an immersed interface method for microstructure scale electrochemical modeling of batteries. The innovation in this approach is the resolution of chemical species and electrostatic potential discontinuities at active interfaces without the use of interface conforming unstructured grids. A unified formulation on Cartesian grids for all domains (electrodes and electrolyte) is used with interfacial flux conditions applied using volume fraction or “color” function gradients. We have developed one dimensional and two dimensional test cases with analytic solutions for electrochemical modeling using which we verified the consistency and accuracy of our scheme. Our solver is also validated against solutions from a macroscale model and an unstructured multi-subdomain solver for a full lithium ion cell. We then demonstrated the utility of our solver on an image-based complex battery electrode microstructure at high charging rate. Our technique also exhibits good scalability on distributed memory architectures using central processing units (CPU), with problem sizes up to 1.8 billion degrees of freedom and with 5400 ranks. Initial performance studies of our open-source performance portable solver showed about 70 times speed up using a graphics processing unit (GPU) compared to single compute core for a problem with 4 million cells.

25 ENERGY STORAGE↗
DOE OSTI2025

A new tubular ironless linear induction motor using non-interrupted winding method

The ironless linear induction motor device can be used for any application which requires a large motion span or rotation restrictedly by the ironless feature due to the strong magnetic field environment. One of the possible applications is ITER diagnostics shutters. The tokamak reactor systems have a severe radiation issue during the operation. The traditional Kapton-insulated wires have significantly lower radiation resistance than an expected large radiation dose at ITER Motion Stark Effect (MSE) shutter location. Furthermore, the traditional single-sided linear induction motor winds the wire to make the N-turn pitch coils independently and then connect them in series to make a one-phase winding. Therefore, it is a big challenge to make so many radiation-protection connectors in such a limited space as inside the vacuum vessel of fusion reactors. The approved mineral-insulated kinds of cables need a large bend radius for the end turns, so the traditional single-sided linear motor with two-layer three-phase windings need to be upgraded to the tubular linear motors. The new proposed tubular ironless linear induction motor can not only overcome the large bend radius issue with the mineral-insulated cables but also allow the motor coils wound continuously around the tubular-shaped stator, forward and backward N times to make a N-turn winding for each phase, in order to resolve the radiation protection issue with the winding connection. The 3D transient Maxwell electromagnetic model with a tubular three-phase linear induction motor has been analyzed. Furthermore, this paper will present the detailed design and winding method of ironless linear induction motor actuators, and results of the 3D transient numerical simulation of the actuator performance.

Ironless↗
DOE OSTI2026

Electron-impact ionization from excited states method

We present a distorted-wave with exchange (DWE) formulation for computing electron-impact ionization (EII) cross sections within the recently developed excited states method (ESM) for plasmas. The ESM provides a self-consistent quantum-mechanical description of both bound and continuum electronic states in dense plasmas, incorporating finite-temperature and screening effects absent in isolated-atom approaches. Using this framework, we calculate EII cross sections for lithium across a wide range of temperatures and densities. We show that the ESM reproduces the isolated-atom limit at low densities while capturing strong plasma effects, such as pressure ionization and shape resonances, at solid density. Comparisons with isolated-atom DWE and average-atom (AA) calculations reveal that, although AA cross sections can approximate configuration average behavior, the ESM provides more accurate, state-resolved answers. In conclusion, these results demonstrate that the ESM allows one to obtain rate coefficients for collisional-radiative modeling of plasmas out of local thermodynamical equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Anti-symmetric and positivity preserving formulation of a spectral method for Vlasov-Poisson equations

We analyze the anti-symmetric properties of a spectral discretization for the one-dimensional Vlasov-Poisson equations. The discretization is based on a spectral expansion in velocity with the symmetrically weighted Hermite basis functions, central finite differencing in space, and an implicit Runge Kutta integrator in time. The proposed discretization preserves the anti-symmetric structure of the advection operator in the Vlasov equation, resulting in a stable numerical method. We apply such discretization to two formulations: the canonical Vlasov-Poisson equations and their continuously transformed square-root representation. The latter preserves the positivity of the particle distribution function. We derive analytically the conservation properties of both formulations, including particle number, momentum, and energy, which are verified numerically on the following benchmark problems: manufactured solution, linear and nonlinear Landau damping, two-stream instability, bump-on-tail instability, and ion-acoustic wave.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

An efficient level set method for tracking many materials

Here, we present an efficient level set method to track an arbitrary number of materials. The algorithm is optimal in the sense that it only needs to store a single unsigned distance-like function and a single integer indicator function, independent of the number of materials or distinct regions being tracked. Furthermore, for smooth velocity fields and smooth interface shape, arbitrarily high order solutions can be demonstrated. For interfaces that are or become kinked, the solution is limited to second-order convergence rates in the L 1 norm and first-order in the L ∞ norm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Multi-GPU porting of a phase-change cascaded lattice Boltzmann method for three-dimensional pool boiling simulations

The Lattice Boltzmann method (LBM) has proven effective in simulating phase-change phenomena, such as melting, solidification, evaporation, and boiling. In this work, we develop a highly parallelized multi-GPU implementation of LBM for three-dimensional pool boiling simulations. The code is based on the OpenACC programming model, which enables the code to be deployed efficiently on multi-core CPUs, GPUs, and potentially other accelerators, without the need for architecture-specific rewrites. To support large-scale simulations, the domain is decomposed and distributed across multiple compute nodes using MPI. We demonstrate that the code exhibits excellent scaling properties, with ideal strong-scaling running with up to 256 GPUs on the MareNostrum5 cluster.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Crystal Growth of Quaternary Rare Earth Selenosilicates by Using the Flux-Assisted Boron Chalcogen Mixture Method: Investigation of Their Magnetic and Optical Properties

A series of rare earth magnesium selenosilicates, RE3Mg 0.5 SiSe 7 (RE = Ce, Pr, Nd, Sm, Gd, Tb, Dy) were obtained as single crystal using the flux assisted boron chalcogen mixture (BCM) method. The structures of the crystals were determined by single-crystal X-ray diffraction. The RE3Mg 0.5 SiSe 7 series crystallizes in the hexagonal crystal system in the space group P6 3 . Polycrystalline powders were synthesized to perform physical property measurements. Magnetic measurements over the 2–300 K temperature range reveal that Ce 3 Mg 0.5 SiSe 7 and Gd 3 Mg 0.5 SiSe 7 exhibit paramagnetic behavior with negative Weiss constants (θ W = −14.50, θ W = −6.13 K, respectively). The optical properties of RE 3 Mg 0.5 SiSe 7 (RE = Ce, Pr, Nd, Sm, Gd) were measured by ultraviolet–visible (UV–vis) diffuse reflectance. Density functional theory (DFT) electronic structure calculations were performed. Furthermore, a second harmonic generation measurement was performed on a polycrystalline powder of Ce 3 Mg 0.5 SiSe 7 and was found to be SHG active with an efficiency of 0.11 times the standard potassium dihydrogen phosphate (KDP).

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

Crystal Growth of Quaternary RE 2 EuSi 2 S 8 ( RE = Ce–Nd, Sm, Gd, Tb) Using Flux-Assisted Boron Chalcogen Mixture (BCM) Method: Investigation of Magnetic and Luminescence Properties

A series of quaternary rare-earth containing thiosilicates with the general formula RE 2 EuSi 2 S 8 (RE = Ce–Nd, Sm, Gd, Tb) has been synthesized via the flux-assisted boron chalcogen mixture (BCM) crystal growth method. High-quality single crystals were obtained, and their crystal structures were determined by single-crystal X-ray diffraction. The RE 2 EuSi 2 S 8 series crystallizes in the trigonal system, adopting the space group R-3c. Polycrystalline samples were employed for physical property measurements, including magnetic susceptibility measurements, UV–visible diffuse reflectance, and photoluminescent response. Magnetic data of RE 2 EuSi 2 S 8 ( RE = Ce, Nd, and Gd) were collected over the 2–300 K temperature range. The samples were paramagnetic behavior with negative Weiss constants (θ W = −11.05, −10.55, and −1.35 K respectively). Their thermal stability was investigated using thermogravimetric analysis (TGA). Optical band gaps, estimated from diffuse reflectance spectra, were determined to be 2.2(1) eV for Ce 2 EuSi 2 S 8 , 1.8(1) eV for Nd 2 EuSi 2 S 8 , and 1.7(1) eV for Gd 2 EuSi 2 S 8 respectively. Finally, photoluminescence measurements were collected on Ce 2 EuSi 2 S 8 and Tb 2 EuSi 2 S 8 single crystals.

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

Structural Investigation of Six Quinary Sulfides Synthesized via the Flux-Assisted Boron-Chalcogen Mixture (BCM) Method: Eu 2+ Containing Members of the RE 3 MTQ 7 (M and T = Transition or Main Group Metals, Q = Chalcogens) Family

For this work, a series of six quinary rare-earth sulfides Ce 4+ 1.85 Eu 2+ 1.15 Na 0.30 SiS 7 , Ce 4+ 1.91 Eu 2+ 1.09 K 0.18 SiS 7 , Ce 4+ 1.96 Eu 2+ 1.04 Rb 0.08 SiS 7 , Ce 4+ 1.98 Eu 2+ 1.02 Cs 0.05 SiS 7 , Ce 4+ 1.97 Eu 2+ 1.03 Ag 0.06 SiS 7 , and Ce 4+ 1.50 Eu 2+ 1.50 CuSiS 7 were obtained in an alkali iodide flux using the boron-chalcogen mixture (BCM) method. Single crystal X-ray diffraction was used to determine the structures of the high quality single crystals that were grown; their elemental compositions were confirmed by energy-dispersive spectroscopy (EDS). The compounds crystallize in the hexagonal crystal system in the noncentrosymmetric space group P63. The crystal structure consists of a three-dimensional network composed of mixed cerium and europium bicapped trigonal prisms, isolated SiS4 tetrahedra, and monovalent metals (Na, K, Rb, Cs, Ag, and Cu) located in cavities created by linked Ce/EuS 8 polyhedra. The structures are charge-balanced when Ce and Eu are in their +4 and +2 oxidation states, respectively. The effective magnetic moment of Ce 1.50 4+ Eu 1.50 2+ CuSiS 7 determined from the temperature dependence of the magnetic susceptibility data is consistent with the presence of Ce 4+ and Eu 2+ . Clear correlations between the alkali ion site occupancy, the ionic radius of the alkali cations, and the average bond length of Ce 4+ /Eu 2+ –S, were established. UV–vis diffuse reflectance data were collected for Ce 1.50 4+ Eu 1.50 2+ CuSiS 7 and a band gap of 1.9(1) eV was established.

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

Flux-Assisted Boron Chalcogen Mixture (BCM) Method for Synthesizing Mixed Chalcogenide Semiconductors (AkRE 2 Si 2 Se x S 8– x and CaRE 2 Si 2 Se 8 ) ( Ak = Ca and Sr; RE = La, Ce, Pr, Nd, and Sm): Investigation of Their Magnetic and Optical Properties

We report a detailed structural analysis of a series of ten quaternary rare-earth-containing seleno-thiosilicates AkRE 2 Si 2 Se x S 8-x and selenosilicates, CaRE 2 Si 2 Se 8 (Ak = Ca and Sr; RE = La, Ce, Pr, Nd, and Sm). Single crystals were obtained by using the flux-assisted boron chalcogen mixture (BCM) method and single crystal X-ray diffraction was used to determine their structures. All members of the AkRE 2 Si 2 Se x S 8-x and CaRE 2 Si 2 Se 8 series crystallize in the space group R$\bar{3}$c (space group number 167) of the trigonal crystal system. The single-crystal X-ray diffraction analysis revealed a strong preference for Se/S atoms to occupy one vs. the other of the two available sites. Polycrystalline samples were used for magnetic susceptibility and UV–visible diffuse reflectance measurements. Magnetic measurements show that CaCe₂Si₂Se₁.₇₃S₆.₂₇ and CaNd₂Si₂Se₂.₅S₅.₅ are paramagnetic with negative Weiss constants (θ = –60.1 and –26.2). Diffuse reflectance analysis gives optical band gaps of 2.7(1) eV (CaLa₂Si₂Se₂.₃₈S₅.₆₂), 2.2(1) eV (CaCe₂Si₂Se₁.₇₃S₆.₂₇), 2.5(1) eV (CaNd₂Si₂Se₂.₅S₅.₅), and 2.0(1) eV (CaCe₂Si₂Se₈), consistent with density functional theory calculations. By partially or fully replacing S sites with Se it was possible to achieve band gap tuning. Photoluminescence behavior was also investigated for CaCe 2 Si 2 Se 1.73 S 6.27 via irradiation with 375 nm ultraviolet light.

BCM method↗
DOE OSTI2025

Synthesis of a series of rare-earth-based multi-anion chalcogenide iodides RE 3 Si 2 Se x S 8− x I (RE = La, Ce, Pr, and Nd) using the flux-assisted boron–chalcogen mixture method

Single crystals and polycrystalline powders of rare earth mixed chalcogenide iodides La 3 Si 2 Se 1.21 S 6.79 I, Ce 3 Si 2 Se 1.39 S 6.61 I, Pr 3 Si 2 Se 1.22 S 6.78 I, and Nd 3 Si 2 Se 1.18 S 6.82 I were prepared using the reactive flux-assisted boron–chalcogen mixture (BCM) method at 850 °C. All compounds crystallized in the monoclinic crystal system, space group C2/c (space group number 15). The series adopts the La 3 Si 2 O 8 Cl structure type, containing isolated SiQ 4 tetrahedra connected by REQ 8 (RE = La, Ce, Pr and Nd) polyhedra; this arrangement creates tunnels that are filled by I atoms. The partial substitution of S by Se was carried out to modulate the optical properties. Phase pure samples and uniform solid solutions were obtained for all compositions as determined using powder X-ray diffraction patterns. Polycrystalline powders were used for physical property measurements, including magnetic susceptibility and UV-Vis diffuse reflectance. The solid-state UV-Vis data for the polycrystalline La 3 Si 2 Se 1.21 S 6.79 I, Ce 3 Si 2 Se 1.39 S 6.61 I, and Pr 3 Si 2 Se 1.22 S 6.78 I samples revealed band gaps of E g = 2.5(1), 2.2(1), and 2.3(1) eV, typical of semiconductors. Magnetic measurements indicated that Ce 3 Si 2 Se 1.39 S 6.61 I and Nd 3 Si 2 Se 1.18 S 6.82 I exhibit paramagnetic behavior with slightly negative Weiss constants θ = −25 and −38. The photoluminescence spectrum of Ce 3 Si 2 Se 1.39 S 6.61 I exhibits a broad emission band around ∼493 nm.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗