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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A High-Order Discontinuous Galerkin Spectral Element Method for Compressible Reacting Flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large eddy simulations due to their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reactive Euler equations encountered in high-speed combustion. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of DG approach. Thus, the framework is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. The numerical method is implemented within the spectral element solver Nek5000. Validation cases are conducted for both non-reactive and reactive discontinuous flows to demonstrate the solver capability. In particular, canonical one-dimensional and two-dimensional detonation simulations are performed and the high-order numerical results are validated against available literature data.

computational fluid dynamics (CFD)↗
DOE OSTI2025

A Finite Element Method for Compressible and Turbulent Multiphase Flow Instabilities with Heat Transfer

We present a new finite element framework for modeling compressible, turbulent multiphase flows with heat transfer. For two-fluid systems with a free surface, the Volume of Fluid (VOF) method is implemented without the need for interface reconstruction, while turbulence is resolved using a dynamic Vreman large eddy simulation (LES) model. Unlike most two-phase VOF studies, which neglect heat transfer, the present approach incorporates energy transport equations within the VOF formulation to account for heat exchange, an effect particularly important in turbulent flows. Conjugate heat transfer is often challenging in finite volume methods, which require explicit specification of heat fluxes at the solid–fluid interface, limiting accuracy and predictive capability. By contrast, the finite element formulation does not require heat flux inputs, allowing more accurate and robust simulation of heat transfer between solids and fluids. The method is demonstrated through three representative cases. First, a two-fluid instability with a single-mode perturbation is simulated and validated against analytical growth rates. Second, conjugate heat transfer is examined in a high-temperature flow over a cold metal cylinder, with validation performed both quantitatively—via pressure coefficient comparisons with experimental data—and qualitatively using vector field topology. Finally, compressible spray injection and breakup are modeled, demonstrating the ability of the framework to capture interfacial dynamics and atomization under turbulent, high-speed conditions. In the compressible spray injection and breakup case, the results indicate that the finite element formulation achieved higher predictive accuracy and robustness than the finite-volume method. With the same mesh resolution, the FEM reduced the root mean square error (RMSE) and mean absolute percentage error (MAPE) from 6.96 mm and 26.0% (for the FVM) to 4.85 mm and 12.7%, respectively, demonstrating improved accuracy and robustness in capturing interfacial dynamics and heat transfer. The study also introduced vector field topology to visualize and interpret coherent flow structures and instabilities, offering insights beyond conventional scalar-field analyses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The Advantages of Combining Morphological and Molecular Methods to Characterise Zooplankton Communities: A Case Study of the UNESCO Biosphere Reserve of the Berlengas Archipelago, Portugal

Understanding the ecological dynamics of zooplankton communities is crucial to precisely assessing the health of marine ecosystems and their management. Metabarcoding has contributed to a better understanding of biodiversity in marine environments. However, this methodology still requires protocol optimisation. Here, we used a complementary approach combining molecular and morphological identification methods to identify the zooplankton community inhabiting the Berlengas Archipelago, Portugal. The presence of non-indigenous species was also assessed. The results showed that the metabarcoding approach outperformed the classical morphological identification method, detecting more species with higher resolution. Nevertheless, the classical method was able to identify species that were not detected by the molecular approaches, probably due to a lack of reference data in the databases. The comparison between different molecular approaches showed that COI and bulk DNA gave better results than 18S rRNA and eDNA by detecting higher species diversity. However, complementarities were observed between them. Molecular tools also proved effective in identifying several potential non-indigenous species, identifying, for the first time, several potentially unreported NIS inhabiting the Portuguese marine ecosystems. Overall, our results confirmed the importance of combining both classical and molecular methods to obtain a more refined assessment of the zooplanktonic communities in marine environments.

Simões, Marco (ORCID:0000000179957304)↗
DOE OSTI2024

A comprehensive review of dwell time optimization methods in computer-controlled optical surfacing

Dwell time plays a vital role in determining the accuracy and convergence of the computer-controlled optical surfacing process. However, optimizing dwell time presents a challenge due to its ill-posed nature, resulting in non-unique solutions. To address this issue, several well-known methods have emerged, including the iterative, Bayesian, Fourier transform, and matrix-form methods. Despite their independent development, these methods share common objectives, such as minimizing residual errors, ensuring dwell time's positivity and smoothness, minimizing total processing time, and enabling flexible dwell positions. This paper aims to comprehensively review the existing dwell time optimization methods, explore their interrelationships, provide insights for their effective implementations, evaluate their performances, and ultimately propose a unified dwell time optimization methodology.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Beyond pinball loss: Quantile methods for calibrated uncertainty quantification

Amongthemanywaysofquantifying uncertainty in a regression setting, specifying the full quantile function is attractive, as quantiles are amenable to interpretation and evaluation. A model that predicts the true conditional quantiles for each input, at all quantile levels, presents a correct and efficient representation of the underlying uncertainty. To achieve this, many current quantile-based methods focus on optimizing the pinball loss. However, this loss restricts the scope of applicable regression models, limits the ability to target many desirable properties (e.g. calibration, sharpness, centered intervals), and may produce poor conditional quantiles. In this work, we develop new quantile methods that address these shortcomings. In particular, we propose methods that can apply to any class of regression model, select an explicit balance between calibration and sharpness, optimize for calibration of centered intervals, and produce more accurate conditional quantiles. We provide a thorough experimental evaluation of our methods, which includes a high dimensional uncertainty quantification task in nuclear fusion.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Analysis of Niobium Electropolishing Using a Generalized Distribution of Relaxation Times Method

Using electrochemical impedance spectroscopy, we have devised a method of sensing the microscopic surface conditions on the surface of niobium as it is undergoing an electrochemical polishing (EP) treatment. The method uses electrochemical impedance spectroscopy (EIS) to gather information on the surface state of the electrode without disrupting the polishing reaction. The EIS data is analyzed using a so-called distribution of relaxation times (DRT) method. Using DRT, the EIS data can be deconvolved into discrete relaxation time peaks without any a priori knowledge of the electrode dynamics. By analyzing the relaxation time peaks, we are able to distinguish two distinct modes of the EP reaction. As the polishing voltage is increased, the electrode transitions from the low voltage EP mode, characterized by a single relaxation time peaks, to the high voltage EP mode, characterized by two relaxation time peaks. We theorize that this second peak is caused by the formation of an oxide layer on the electrode. We also find that this oxide induced peak transitions from to a negative relaxation time, which is indicative of a blocking electrode process. By analyzing EPed samples, we show that samples polished in the low voltage mode have significantly higher surface roughness due to grain etching and faceting. We find that the surface roughness of the samples only improves when the oxide film peak is present and in the negative relaxation time region. This shows that EIS combined with DRT analysis can be used to predict etching on EPed Nb. This method can also be performed before or during the EP, which could allow for adjustment of polishing parameters to guarantee a smooth cavity surface finish.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Electromagnetic field-assisted method for chemical conversion

Embodiments relate to methods for enhancing chemical conversions. One or more embodiments relate to a method for enhancing a multi-step chemical conversion reaction. The method includes providing a reactant mixture comprising one or more reacting specie(s); and providing a catalyst or sorbent comprising one or more support materials and one or more deposited catalytically active materials. The method further includes applying an electromagnetic field with a prescribed power, frequency, and pulsing strategy specific to interactions of reactant species and an electromagnetic field with at least one of the support materials, sorbent, and catalytically active materials in a particular chemical reaction.

Shekhawat, Dushyant↗
DOE OSTI2024

Method for electrochemical transformation of amorphous material to crystalline material

A method for converting amorphous boron nitride to crystalline boron nitride, the method comprising immersing the amorphous boron nitride into anhydrous molten magnesium chloride maintained within a temperature range of 720° C.-820° C. while the amorphous boron nitride is cathodically polarized at a voltage within a range of −2.2V to −2.8V for a period of time of at least 2 minutes to result in conversion of the amorphous boron nitride to the crystalline form. Also described herein is a method for converting an amorphous carbon material to a crystalline carbon material, the method comprising immersing said amorphous carbon material into anhydrous molten magnesium chloride maintained within a temperature range of 780° C.-820° C. while the amorphous carbon material is cathodically polarized at a voltage within a range of −2.2V to −2.8V for a period of time of at least 2 minutes to result in conversion of the amorphous carbon material to the crystalline form.

Bagri, Prashant↗
DOE OSTI2025

Method and system embodiments for converting ethanol to para-xylene and ortho-xylene

Disclosed herein are embodiments of a method and system for converting ethanol to para-xylene. The method also provides a pathway to produce terephthalic acid from biomass-based feedstocks. In some embodiments, the disclosed method produces p-xylene with high selectivity over other aromatics typically produced in the conversion of ethanol to xylenes, such as m-xylene, ethyl benzene, benzene, toluene, and the like. And, in some embodiments, the method facilitates the ability to use ortho/para mixtures of methylbenzyaldehyde for preparing ortho/para xylene product mixtures that are amendable to fractionation to separate the para- and ortho-xylene products thereby providing a pure feedstock of para-xylene that can be used to form terephthalic anhydride and a pure feedstock of ortho-xylene that can be used for other purposes, such as phthalic anhydride.

Ramasamy, Karthikeyan Kallupalayam↗
DOE OSTI2025

Methods of improving drought and salt resistance in a plant and genetically engineered plants with improved drought and salt resistance

The present disclosure provides methods for increasing drought resistance, salt resistance, and biomass production of a plant. The methods encompass expression of DiGeorge-Syndrome Critical Region 14 (DGCR14) gene in the plant. In comparison to a plant not manipulated in this manner, the disclosed, genetically-modified, plants display improved drought resistance and salt resistance. Also provided are plants that can be obtained by the method according to the invention, and nucleic acid vectors to be used in the described methods.

Xie, Meng↗
DOE OSTI2025

Opportunities and Challenges of Computational Electromagnetics Methods for Superconducting Circuit Quantum Device Modeling: A Practical Review

High-fidelity numerical methods that model the physical layout of a device are essential for the design of many technologies. For methods that characterize electromagnetic effects, these numerical methods are referred to as computational electromagnetics (CEM) methods. Although the CEM research field is mature, emerging applications can still stress the capabilities of the techniques in use today. The design of superconducting circuit quantum devices falls in this category due to the unconventional material properties and important features of the devices covering nanometer to centimeter scales. Such multiscale devices can stress the fundamental properties of CEM tools which can lead to an increase in simulation times, a loss in accuracy, or even cause no solution to be reliably found. While these challenges are being investigated by CEM researchers, knowledge about them is limited in the broader community of users of these CEM tools. This review is meant to serve as a practical introduction to the fundamental aspects of the major CEM techniques that a researcher may need to choose between to model a device, as well as provide insight into what steps they may take to alleviate some of their challenges. Our focus is on highlighting the main concepts without rigorously deriving all the details, which can be found in many textbooks and articles. After covering the fundamentals, we discuss more advanced topics related to the challenges of modeling multiscale devices with specific examples from superconducting circuit quantum devices. We conclude with a discussion on future research directions that will be valuable for improving the ability to successfully design increasingly more sophisticated superconducting circuit quantum devices. Although our focus and examples are taken from this area, researchers from other fields will still benefit from the details discussed here.

Elkin, Samuel T. [Virginia Tech., Blacksburg]↗
DOE OSTI2024

A high-order Shifted Interface Method for Lagrangian shock hydrodynamics

Here, we present a new method for two-material Lagrangian hydrodynamics, which combines the Shifted Interface Method (SIM) with a high-order Finite Element Method. Our approach relies on an exact (or sharp) material interface representation, that is, it uses the precise location of the material interface. The interface is represented by the zero level-set of a continuous high-order finite element function that moves with the material velocity. This strategy allows to evolve curved material interfaces inside curved elements. By reformulating the original interface problem over a surrogate (approximate) interface, located in proximity of the true interface, the SIM avoids cut cells and the associated problematic issues regarding implementation, numerical stability, and matrix conditioning. Accuracy is maintained by modifying the original interface conditions using Taylor expansions. We demonstrate the performance of the proposed algorithms on established numerical benchmarks in one, two and three dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Method-independent cusps for atomic orbitals in quantum Monte Carlo

Here, we present an approach for augmenting Gaussian atomic orbitals with correct nuclear cusps. Like the atomic orbital basis set itself and unlike previous cusp corrections, this approach is independent of the many-body method used to prepare wave functions for quantum Monte Carlo. Once the basis set and molecular geometry are specified, the cusp-corrected atomic orbitals are uniquely specified, regardless of which density functionals, quantum chemistry methods, or subsequent variational Monte Carlo optimizations are employed. We analyze the statistical improvement offered by these cusps in a number of molecules and find them to offer similar advantages as molecular-orbital-based approaches while remaining independent of the choice of many-body method.

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

Evaluation of Density-Functional Tight-Binding Methods for Simulation of Protic Molecular Ion Pairs

In this work, we benchmark the accuracy of the density-functional tight-binding (DFTB) method, namely the long-range corrected second-order (LC-DFTB2) and third-order (DFTB3) models, for predicting energetics of imidazolium-based ionic liquid (IL) ion pairs. We compare the DFTB models against popular density functionals such as LC-ωPBE and B3LYP, using ab initio domain-based local pair-natural orbital coupled cluster (DLPNO-CC) energies as reference. Calculations were carried out in the gas phase, as well as in aqueous solution using implicit solvent methods. We find that the LC-DFTB2 model shows excellent performance in the gas phase and agrees well with reference energies in implicit solvent, often outperforming DFTB3 predictions for complexation energetics. Our study identifies a range of opportunities for use of the LC-DFTB method and quantifies its sensitivity to protonation states and the types of chemical interactions between ion pairs.

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

Arbitrary Order Virtual Element Methods for High‐Order Phase‐Field Modeling of Dynamic Fracture

ABSTRACT Accurate modeling of fracture nucleation and propagation in brittle and ductile materials subjected to dynamic loading is important in predicting material damage and failure under extreme conditions. Phase‐field fracture models have garnered a lot of attention in recent years due to their success in representing damage and fracture processes in a wide class of materials and under a variety of loading conditions. Second‐order phase‐field fracture models are by far the most popular among researchers (and increasingly, among practitioners), but fourth‐order models have started to gain broader acceptance since their more recent introduction. The exact solution corresponding to these high‐order phase‐field fracture models has higher regularity. Thus, numerical solutions of the model equations can achieve improved accuracy and higher spatial convergence rates. In this work, we develop a virtual element framework for the high‐order phase‐field model of dynamic fracture. The virtual element method (VEM) can be regarded as a generalization of the classical finite element method. In addition to many other desirable characteristics, the VEM allows computing on polytopal meshes. Here, we use ‐conforming virtual elements and the generalized‐ time integration method for the momentum balance equation, and adopt ‐conforming virtual elements for the high‐order phase‐field equation. We verify our virtual element framework using classical quasi‐static benchmark problems and demonstrate its capabilities with the aid of numerical simulations of dynamic fracture in brittle materials.

42 ENGINEERING↗
DOE OSTI2025

Efficient proximal subproblem solvers for a nonsmooth trust-region method

In [R. J. Baraldi and D. P. Kouri, Mathematical Programming, (2022), pp. 1-40], we introduced an inexact trust-region algorithm for minimizing the sum of a smooth nonconvex and nonsmooth convex function. The principle expense of this method is in computing a trial iterate that satisfies the so-called fraction of Cauchy decrease condition—a bound that ensures the trial iterate produces sufficient decrease of the subproblem model. In this paper, we expound on various proximal trust-region subproblem solvers that generalize traditional trust-region methods for smooth unconstrained and convex-constrained problems. We introduce a simplified spectral proximal gradient solver, a truncated nonlinear conjugate gradient solver, and a dogleg method. Finally, we compare algorithm performance on examples from data science and PDE-constrained optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Crack opening calculation in phase-field modeling of fluid-filled fracture: A robust and efficient strain-based method

The phase-field method has become popular for the numerical modeling of fluid-filled fractures, thanks to its ability to represent complex fracture geometry without algorithms. However, the algorithm-free representation of fracture geometry poses a significant challenge in calculating the crack opening (aperture) of phase-field fracture, which governs the fracture permeability and hence the overall hydromechanical behavior. Although several approaches have been devised to compute the crack opening of phase-field fracture, they require a sophisticated algorithm for post-processing the phase-field values or an additional parameter sensitive to the element size and alignment. Here, we develop a novel method for calculating the crack opening of fluid-filled phase-field fracture, which enables one to obtain the crack opening without additional algorithms or parameters. Here we transform the displacement-jump-based kinematics of a fracture into a continuous strain-based version, insert it into a force balance equation on the fracture, and apply the phase-field approximation. Through this procedure, we obtain a simple equation for the crack opening which can be calculated with quantities at individual material points. We verify the proposed method with analytical and numerical solutions obtained based on discrete representations of fractures, demonstrating its capability to calculate the crack opening regardless of the element size or alignment.

58 GEOSCIENCES↗
DOE OSTI2025

A provably stable numerical method for the anisotropic diffusion equation in confined magnetic fields

We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial boundary value problem. A discrete formulation is presented using operator splitting in time with the summation by parts finite difference approximation of spatial derivatives for the perpendicular diffusion operator. Weak penalty procedures are derived for implementing both boundary conditions and parallel diffusion operator obtained by field line tracing. We prove that the fully-discrete approximation is unconditionally stable. Discrete energy estimates are shown to match the continuous energy estimate given the correct choice of penalty parameters. A nonlinear penalty parameter is shown to provide an effective method for tuning the parallel diffusion penalty and significantly minimises rounding errors. Several numerical experiments, using manufactured solutions, the “NIMROD benchmark” problem and a single island problem, are presented to verify numerical accuracy, convergence, and asymptotic preserving properties of the method. Finally, we present a magnetic field with chaotic regions and islands and show the contours of the anisotropic diffusion equation reproduce key features in the field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗