SEARCH · Search NASA

Results for “Implementation”

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 595 records · Page 33

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Biogenic carbon capture at pulp mills via sodium spiking and oxy-fuel calcination

Over 13 million metric tons of biogenic carbon dioxide (CO 2 ) are mineralized yearly in United States (US) pulp mill recovery boilers as molten sodium carbonate. The mineralized CO 2 is released downstream in the rotary lime kiln where it can be captured at a relatively low cost for permanent sequestration. The use of biocarbon and bioenergy at pulp mills enables the possibility for atmospheric carbon removal when biogenic CO2 is captured and sequestered. We demonstrate the feasibility of capturing biogenic CO 2 via sodium spiking coupled with oxy-fuel calcination in the rotary lime kiln at existing kraft pulp mills. Sodium (Na) spiking elevates the Na ion loading in the kraft chemical looping process by replacing chlorine-based bleaching with a highly alkaline bleaching sequence. Chemical pulping processes are modeled and simulated to understand the technical limits of implementing sodium spiking in existing pulp mills. Each 1 % of sodium added to the bleaching operations increases the rate of CO 2 mineralization by 4 %, and the maximum increase in sodium content in the kraft process is 9 %. Without sodium spiking, estimated CO 2 capture costs are $131 and $107/mt CO 2 for air and oxy-fuel combustion, respectively. With the implementation of sodium spiking, the cost of CO 2 captured decreases by 27 % and 31 %, showing costs of $96 and $74/mt CO 2 for air and oxy-fuel combustion, respectively. Oxy-fuel combustion reduces the costs of CO 2 capture compared to air combustion, but merged deployment with sodium spiking seems to be the most cost-effective pathway to integrate carbon capture.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2026

A decade of progress in understanding and managing legacy well integrity for geologic carbon storage

This study reviews a decade of research progress in legacy well integrity and risk management for geologic carbon storage (GCS) to commemorate the 20 th anniversary of the Intergovernmental Panel on Climate Change’s 2005 Special Report on Carbon Capture and Storage. In the past ten years, legacy well research has benefited from global efforts to constrain emissions from abandoned oil and gas wells, a continued focus on well materials performance in the presence of CO 2 -rich fluids, and practical experience gained through GCS implementation. Field measurements of abandoned well emissions show that leakage is not universal or catastrophic but forms a continuum of low-to-moderate fluxes that depend on isolation integrity and environmental attenuation. Materials research has constrained the conditions under which Portland cements exhibit self-sealing and non-sealing behaviors, and has identified the impact of geomechanical properties, non-uniform pathway apertures, multi-phase flow, and impurities in the CO 2 stream, on leakage pathways as important new areas for investigation. GCS projects at brownfield sites have inspired the creation of new workflows that integrate various tools and technologies to manage legacy well leakage risks. GCS implementation has also motivated a push towards scenario-based well modeling that directly informs permit applications. These advances inspire new research questions for the coming decade, particularly around the level of legacy well leakage risk that is environmentally acceptable and tolerable to stakeholders when sequestering millions of tonnes of CO 2 annually.

Carbon capture and storage↗

DOE OSTI2024

Modeling inter- and intra-granular dislocation transport using crystal plasticity

Here, this work presents the development of a crystal plasticity material model that incorporates both dislocation transport within grains and dislocation transfer across grain boundaries. This model has been implemented in the open-source finite element code MOOSE. In addition, a novel geometry-based criterion is developed to determine the direction of dislocation transfer across grain boundaries. The transfer criterion incorporates the geometric features of the grain boundary, such as the grain boundary plane normal, and its misorientation, which is accounted for through the orientation of the incoming and outgoing slip systems. The model is tested with several cases, including a copper single crystal, bi-crystal, and polycrystal. The development of the transfer criterion, implementation of the model, and its application to these test cases are discussed in detail.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Manufactured solutions for an electromagnetic slot model

The accurate modeling of electromagnetic penetration is an important topic in computational electromagnetics. Electromagnetic penetration occurs through intentional or inadvertent openings in an otherwise closed electromagnetic scatterer, which prevent the contents from being fully shielded from external fields. To efficiently model electromagnetic penetration, aperture or slot models can be used with surface integral equations to solve Maxwell's equations. A necessary step towards establishing the credibility of these models is to assess the correctness of the implementation of the underlying numerical methods through code verification. Surface integral equations and slot models yield multiple interacting sources of numerical error and other challenges, which render traditional code-verification approaches ineffective. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources for the method-of-moments implementation of the electric-field integral equation with a slot model. Finally, we demonstrate the effectiveness of these approaches for a variety of cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Asymptotic-preserving semi-implicit finite volume scheme for extended magnetohydrodynamics

A Finite Volume (FV) scheme is developed for solving the extended magnetohydrodynamic (XMHD) equations, yielding accurate results in the ideal, resistive, and Hall MHD limits. This is accomplished by first re-writing the XMHD equations such that it allows the algorithm to retain the use of ideal MHD Riemann solvers and the constrained transport method to preserve divergence-free magnetic fields. Incorporation of electron inertia and displacement current introduces additional numerical stiffness which motivates a semi-implicit FV scheme that re-formulates the XMHD model as a relaxation system. The equations are then advanced in time using an explicit 2nd-order Runge–Kutta scheme with operator splitting applied to the implicit source term updates at each sub-stage. For additional numerical stability, a density-dependent slope limiter is implemented to increase flux diffusivity at low density regions where non-ideal effects become significant. The algorithm is subsequently implemented in a scalable adaptive mesh refinement (AMR) framework. As the new algorithm retains many aspects of the ideal MHD formulations, it asymptotes naturally to the ideal MHD limit. Moreover, it shows promising results at the resistive and Hall MHD limits. This is verified against reference test problems for ideal, resistive and Hall MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Assessing the viability of a new subsize tensile specimen geometry for evaluation of structural nuclear and additively manufactured materials

The use of subsize specimens in nuclear materials testing has been a subject of ongoing interest due to radiation safety concerns and resource conservation needs. It is also of interest in additive manufacturing (AM), also known as 3D-printing, as subsize specimens can more accurately represent the behavior of the small, intricate geometries often produced using AM. A novel, extremely small geometry, called the Subsize Teeny (SST) was recently developed and is of interest for implementation. Given its extraordinarily small size and the complexities associated with subsize specimen testing, adequate vetting of this geometry is necessary to ensure data quality. Here, a variety of unirradiated nuclear structural materials were tested in the SST geometry and compared against the well-established SSJ3 geometry. In addition, two case studies implementing the SST as a screening geometry for AM materials were also conducted. The question of SST viability was found to be highly nuanced and will often be dependent on the context or application in question. It was determined, however, that the SST is a largely invalid geometry for exceptionally coarse-grained materials or in cases where the physical defect volume equals or exceeds 0.1 % or where the specimen machining parameters result in significant surface alterations. On the other hand, it was determined that the SST may be employed with confidence if the test material is nearly or totally free of physical defects, isotropic, demonstrates homogeneous plastic deformation, and possesses a fine-grained, nearly or totally homogeneous microstructure with at least twelve slip systems.

Additive manufacturing↗

DOE OSTI2026

Documentation, assessment, and improvements of fission gas modeling capabilities in BISON

Fission gas behavior critically affects nuclear fuel performance, influencing thermal conductivity, microstructure, swelling, and fuel-cladding interactions. The Simple Integrated Fission Gas Release and Swelling (Sifgrs) model, implemented within the BISON fuel performance code, provides advanced capabilities for modeling these phenomena. This manuscript reviews and assesses Sifgrs’ modular structure and submodels describing intragranular, intergranular, and release mechanisms. The modular implementation enables rapid improvements as our understanding advances and new fuel forms emerge. After an overview of capabilities, we present a comprehensive validation suite, comparing predictions to experimental data from tests under normal and transient conditions, including reproduction of the Vitanza curve. Results demonstrate Sifgrs’ accuracy, robustness, and applicability. Current limitations and ongoing efforts are discussed, including mechanistic, multiscale modeling, uncertainty quantification, and validation within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, supporting the safe and efficient use of nuclear fuels.

BISON↗

DOE OSTI2025

A product data network to enable faster, easier, and better planning of building envelopes

The building envelopes contributes significantly to the energy-efficiency of the building. Building performance simulation has made it possible to compare façade technologies regarding energy demand, daylighting, thermal and visual comfort in detail. Planners, such as architects and engineers, need experience to find product data with the right quality and level of detail, and to process the data to fit the calculation and the application. In the available time, planners can compare only a limited number of products, which means that better solutions could go unnoticed. This paper presents a new concept for making product data easily accessible for building façade planning. The concept consists of a network of databases for the efficient exchange and use of optical and calorimetric data of glazing units, shading devices, and combinations of both. The paper presents the research questions, an analysis of the current challenges, six design goals for the product data network and its implementation together with a discussion. Many product data sources can be connected to many planning software applications via the specified application programming interface. When planning software connects to the product data network, the planning of building envelopes can be much faster because planners do not need to spend so much time to search and process product data manually. The planning of building envelopes can also become much easier, especially for planners with limited experience. They do not need to understand all the details about which data fits which calculation if the software company implements this. The planning of building envelopes can become much more reliable when software companies validate their use of the product data network, because the current manual process is prone to errors. The planning of building envelopes can also improve because more products can be compared in the available time, allowing better solutions to be found.

Maurer, Christoph↗

DOE OSTI2024

Accelerating Multivariate Functional Approximation Computation with Domain Decomposition Techniques⋆

Modeling large datasets through Multivariate Functional Approximations (MFA) provide an elegant way to handle many visualization and scientific analysis workflows. The process necessitates scalable data partitioning methods to compute MFA representations efficiently without compromising the accuracy or continuity of the reconstructed solution. We propose a domain -decomposed method for computing the MFA with B -spline bases, which reduces the total work per task and uses a restricted Additive Schwarz (RAS) method to converge the control point data degrees -of -freedom along subdomain boundaries. We provide an in-depth analysis of the parallel approach with domain decomposition solvers, aiming to minimize local subdomain error residuals and recover high -order continuity at subdomain interfaces with appropriate choices of knot overlaps. The communication cost, determined by the overlap regions in the RAS implementation, is optimized to recover the numerical error profile of the single subdomain case. Our proposed method stands in contrast to previous methods, which typically only recover either C 0 or at best C 1 continuity for arbitrary B -spline degree expansions, or those that require post -processing to blend discontinuities in the reconstructed data. We demonstrate the effectiveness of our approach using analytical and real -world datasets in 1D, 2D, and 3D through both strong and weak scaling studies. The performance results indicate that the overall cost of computing the approximation is directly proportional to the underlying nearest -neighbor communication implementation, and is only weakly dependent on the overlap region size that determines the size of the messages. This finding underscores the efficiency and scalability of our proposed method, making it a promising solution for handling large datasets in scientific workflows.

additive Schwarz solvers↗

DOE OSTI2024

Development of a multi-layer canopy model for E3SM Land Model with support for heterogeneous computing

The vertical structure of vegetation canopies creates micro-climates. However, the land components of most Earth System Models, including the Energy Exascale Earth System Model (E3SM), typically neglect vertical canopy structure by using a single layer big-leaf representation to simulate water, CO 2 , and energy exchanges between the land and the atmosphere. In this study, we developed a Multi-Layer Canopy Model for the E3SM Land Model to resolve the micro-climate created by vegetation canopies. The model developed in this study re-implements the CLM-ml_v1 to support heterogeneous computing architectures consisting of CPUs and GPUs and includes three additional optimization-based stomatal conductance models. The use of Portable, Extensible Toolkit for Scientific Computation provides a speedup of 25–50 times on a GPU relative to a CPU. The numerical implementation of the model was verified against CLM-ml_v1 for a month-long simulation using data from the Ameriflux US-University of Michigan Biological Station site. Model structural uncertainty was explored by performing control simulations for five stomatal conductance models that exclude and include the control of plant hydrodynamics (PHD) on photosynthesis. The bias in simulated sensible and latent heat fluxes was lower when PHD was accounted for in the model. Additionally, six idealized simulations were performed to study the impact of three environmental variables (i.e. air temperature, atmospheric CO 2 , and soil moisture) on canopy processes (i.e. net CO 2 assimilation, leaf temperature, and leaf water potential). Increasing air temperature reduced net CO 2 assimilation and increased air temperature. Net CO 2 assimilation increased at higher atmospheric CO 2 , while decreasing soil moisture resulted in lower leaf water potential.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Design of an alumina forming coating for Nb-base refractory alloys

Refractory multi-principal element alloys (RMPEAs) promise to significantly enhance gas turbine engine efficiency, but their poor oxidation performance inhibits their implementation. Alumina-forming bond coat alloys can provide oxidation protection, but discovering suitable chemistries remains a challenge. We employed a design methodology that screens for alumina-formation capability using Al activity and phase constitution predictions from CalPhaD (Thermo-Calc). Alloy down-selection from approximately 7,800 alloys in the Nb-Si-Ti-Al-Hf system was conducted via analysis of calculated thermodynamic properties with number-density topology style maps. This approach is validated by creating and testing the composition Nb 12 Si 23 Ti 24 Al 36 Hf 5 , which forms protective alumina scales up to 1400 °C and resists pesting at 800 °C. Further, the alloy has an average coefficient of thermal expansion of ~10.1 ppm/K, making it well matched to Nb-based refractory alloys. The methodology will be useful for the design of coatings for RMPEAs, enabling their implementation and significant efficiency benefits sooner.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Generating Buoyancy-Driven Convection in Membrane Distillation

Membrane distillation (MD) is a thermally-driven desalination process that can treat hypersaline brines. Considerable MD literature has focused on mitigating temperature and concentration polarization. This literature largely neglects that temperature and concentration polarization increase the feed density near the membrane. With gravity properly oriented, this increase in density could trigger buoyancy-driven convection and increase permeate production. Convection could also be strengthened by heating the feed channel wall opposite the membrane. To investigate that possibility, we perform a series of experiments using a plate-and-frame direct contact MD system with an active membrane area of 300 cm2 and a feed channel wall heated using a resistive heater. The experiments measure the average transmembrane permeate flux for two gravitational orientations, feed Reynolds numbers between 128 and 1128, and wall heat fluxes up to 12 kW/m2. The results confirm that with gravity properly oriented, wall-heating can trigger buoyancy-driven convection for a wide range of feed Reynolds numbers, and increase permeate production between roughly 20 and 130 %. We estimate, however, that at high Reynolds numbers (Re > 800), more than 70 % of the wall heat is carried out of the MD system by the feed flow, without contributing to permeate production. This suggests the need for longer membranes and heat recovery steps in any future practical implementation.

buoyancy-driven convection↗

DOE OSTI2025

Monochromatization interaction region optics design for direct s-channel Higgs production at FCC-ee

The FCC-ee offers the potential to measure the electron Yukawa coupling via direct s-channel Higgs production, e + e − → H , at a centre-of-mass (CM) energy of ∼125 GeV. This measurement is significantly facilitated if the CM energy spread of e + e − collisions can be reduced to a level comparable to the natural width of the Higgs boson, Γ H = 4.1 MeV, without substantial loss in luminosity. Achieving this reduction in collision-energy spread is possible through the “monochromatization” concept. The basic idea is to create opposite correlations between spatial position and energy deviation within the colliding beams, which can be accomplished in beam optics by introducing a nonzero dispersion function with opposite signs for the two beams at the interaction point. Since the first proposal in 2016, the implementation of monochromatization at the FCC-ee has been continuously improved, starting from preliminary parametric studies. In this paper, we present a detailed study of the interaction region optics design for this newly proposed collision mode, exploring different potential configurations and their implementation in the FCC-ee global lattice, along with beam dynamics simulations and performance evaluations including the impact of “beamstrahlung.”

43 PARTICLE ACCELERATORS↗

DOE OSTI2025

Development and validation of a software for simulating γ-γ coincidence emission and detection probabilities

Gamma-gamma coincidence spectrometers have the potential to significantly enhance detection sensitivity for ultra-trace radionuclide measurements. The implementation of these spectrometers, however, is limited by the complexity of acquisition hardware, data processing and quantification. This work reports development of a novel radionuclide quantification software for γ-γ coincidence measurements. For any radionuclide, the software parses the Evaluated Nuclear Structure Data File (ENSDF) database, recursively simulating all possible γ-γ coincidence signatures and their respective emission and detection probabilities. Implemented using Python programming language, the software employs several strategies to boost overall computational performance. Since coincidence-based spectrometers are of notable interest in monitoring compliance for the Comprehensive Nuclear-Test-Ban Treaty (CTBT), the software’s execution was tested for 84 CTBT-relevant radionuclides. To date, the software has been experimentally validated for 15 radionuclides using the Advanced Radionuclide Gamma spectrOmeter (ARGO) at Pacific Northwest National Laboratory, USA (PNNL). Notably, the software can be operated in convergence mode, whereby coincidence detection efficiency’s convergence behavior can help avoid unreliable radionuclide activity estimates. With growing number of coincidence spectrometers worldwide, this paper aims to assist the radiation metrology community in developing similar software for their system.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2026

Design and fabrication of a four-rod radio frequency quadrupole with a variable aperture

Radiofrequency quadrupole (RFQ) linear accelerators with variable aperture cell parameters offer an effective means of accelerating high-intensity ion beams exceeding 100 mA and have been implemented in various RFQ designs. However, in a four-rod RFQ, the use of dielectric perturbation methods to measure the voltage distribution behind the rod electrodes is not applicable when implementing variable aperture geometries, because changes in the rod shape influence the electric field at the location of the perturbation element. Therefore, it was necessary to design the rod profile for each cell such that the electric field distribution behind the rods remains unaffected by changes in the aperture. In this study, we present a method to adjust the voltage distribution by modifying the geometry of the rod electrodes according to each cell's aperture parameters. This design approach enables the integration of variable aperture cells into a four-rod RFQ linac, significantly enhancing the achievable beam current. To ensure uniform capacitance and maintain a constant rod-electrode voltage despite variations in the average aperture radius, we developed a systematic method for shaping the rod electrodes. This method was refined using two-dimensional (2D) and three-dimensional (3D) electromagnetic field simulations to optimize the electrode cross-sectional profiles. The fabricated rod electrodes were then installed in a four-rod RFQ, and the electric field distribution was measured using the perturbation method and tuned by adjusting the base plates. As a result, the measured variation of ±1.1 % confirms that this design methodology can produce a sufficiently uniform voltage distribution for high-current beam acceleration. These results demonstrate the practical feasibility of a variable-aperture four-rod RFQ linac.

43 PARTICLE ACCELERATORS↗

DOE OSTI2026

A 24-channel ultra-low-noise preamplifier for dN/dx measurements with drift tube detectors

Cluster counting $(dN/dx)$ is a promising method to enhance particle identification for gaseous detectors, especially in next-generation collider experiments like the FCC-ee, where good $π/K$ separation over a broad momentum range is essential. However, its implementation in large-scale systems has been limited by the challenging requirements for high-resolution signal amplification and readout. This paper presents a 24-channel ultra-low-noise preamplifier board designed for drift tube detectors to enable $dN/dx$ measurements. The three-stage amplification topology employs SiGe transistors and integrates dedicated noise-minimization techniques, achieving a charge gain of 21.11 mV/fC from 0.3 fC to 50 fC, a bandwidth of 542 MHz, and a voltage gain of 47.8 dB. The measured voltage noise density is 0.35 nV/$\sqrt{\textrm{Hz}}$ , surpassing most of the state-of-the-art preamplifiers for gaseous and silicon detectors. Validation tests conducted on the sMDT chambers at the CERN Proton Synchrotron test beam facility demonstrate that the proposed design meets the stringent preamplifier requirements for implementing the $dN/dx$ method in drift-tube detector systems, achieving an equivalent noise charge of 0.14 fC and a signal-to-noise ratio of 73 when operated with a He:iC 4 H 10 (90:10) gas mixture. The design also shows promise for broader application in other gaseous or semiconductor detectors.

Drift tube detector↗

DOE OSTI2026

Design of a high voltage delivery system for noble liquid time projection chambers

Noble liquid time projection chambers (TPCs) are a leading technology in the detection of ionizing radiation, particularly in applications such as accelerator neutrino physics, dark matter detection, and neutrinoless double beta decay. This work addresses the design considerations for implementing stable high voltage (HV) systems within large noble liquid TPCs, with a focus on the nEXO experiment. Utilizing insights from prior HV research and experimental investigations, we outline factors influencing HV stability and discuss design choices to improve stability and prevent electrical discharges. A novel HV delivery system concept is presented, tailored for the nEXO TPC, which incorporates these design considerations while also meeting the stringent radiopurity requirements of the nEXO neutrinoless double beta decay search. These design considerations and their specific implementation towards a HV delivery system offer guidance to future experiments applying high voltage in noble liquid environments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗