SEARCH · Search NASA

Results for “splitting 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.

SourceResearch connectionFrom yearThrough year

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

At least 1,225 records · Page 68

All sources NASA NTRS DOE OSTI

NASA NTRS1986

An isoparametric spectral element method for solution of the Navier-Stokes equations in complex geometry

High-order (p-type) finite element methods combine the geometric flexibility of standard low-order finite element (or finite volume) techniques with the rapid convergence properties of spectral methods. Various p-type schemes have been proposed for elliptic problems. However, an application of p-type methods to more complex equations, in particular the passive scalar or Navier-Stokes equations, is complicated by the nature of the equations and by efficiency considerations as regards the solution procedure. Patera (1984) has proposed a p-type method for the Navier-Stokes equations, taking into account the arising additional complications. The present study is concerned with an isoparametric spectral element discretization for the solution of flow problems in fully general, curvy geometries. Attention is given to the isoparametric formulation for elliptic equations and a time-splitting scheme for the Navier-Stokes equations.

Korczak, K. Z.↗

NASA NTRS1987

A new algorithm for the Navier-Stokes equations applied to transonic flows over wings

A new algorithm has been developed for the Navier-Stokes equations. For the convective terms, this algorithm employs flux vector splitting together with a locally rotated coordinate system that is aligned with the streamwise direction. The algorithm has been implemented into a Navier-Stokes code, the NASA Ames Research Center's TNS code, and several cases of steady flow have been calculated. The results show improvements over the original method in the code. First, in a case of separated flow, there is better agreement with the experimental results and also there is an absence of numerical oscillations that occur with the original method. Second, in a case with a strong supersonic to subsonic shock wave, comparisons show an improvement in the resolution of the shock wave and the reexpansion singularity.

Goorjian, Peter M.↗

NASA NTRS1989

A dual potential formulation of the Navier-Stokes equations

A dual potential formulation for numerically solving the Navier-Stokes equations is developed and presented. The velocity field is decomposed using a scalar and vector potential. Vorticity and dilatation are used as the dependent variables in the momentum equations. Test cases in two dimensions verify the capability to solve flows using approximations from potential flow to full Navier-Stokes simulations. A three-dimensional incompressible flow formulation is also described. An interesting feature of this approach to solving the Navier-Stokes equations is the decomposition of the velocity field into a rotational part (vector potential) and an irrotational part (scalar potential). The Helmholtz decomposition theorem allows this splitting of the velocity field. This approach has had only limited use since it increases the number of dependent variables in the solution. However, it has often been used for incompressible flows where the solution scheme is known to be fast and accurate. This research extends the usage of this method to fully compressible Navier-Stokes simulations by using the dilatation variable along with vorticity. A time-accurate, iterative algorithm is used for the uncoupled solution of the governing equations. Several levels of flow approximation are available within the framework of this method. Potential flow, Euler and full Navier-Stokes solutions are possible using the dual potential formulation. Solution efficiency can be enhanced in a straightforward way. For some flows, the vorticity and/or dilatation may be negligible in certain regions (e.g., far from a viscous boundary in an external flow). It is possible to drop the calculation of these variables then and optimize the solution speed. Also, efficient Poisson solvers are available for the potentials. The relative merits of non-primitive variables versus primitive variables for solution of the Navier-Stokes equations are also discussed.

Gegg, S. G.↗

NASA NTRSDate not supplied

Pilot Workload Rating Predictions Using Image Data and Recurrent Neural Networks

In this work, we augmented existing methods for estimating pilot workload ratings with deep neural networks trained using data from simulated flight tests in the Vertical Motion Simulator (VMS). We used an existing method, Spare Capacity Operations Estimator (SCOPE), along with a recurrent neural network and conducted comparison studies between the two methods individually, and when used together. We found that using both methods together can improve the result over using either approach alone. In our first test case, we achieved an improved linear correlation coefficient of 0.409 over that of SCOPE alone at 0.352 on the training dataset. Through cross validation, we also found that the results may be dependent on the split of training vs. validation data, and that further investigation should be conducted to understand what additional inputs to the neural network model should be made.

Image Data↗

DOE OSTI2020

Performance of position-sensitive resistive silicon detectors in the Forward Array Using Silicon Technology (FAUST)

The silicon telescopes of the Forward Array Using Silicon Technology (FAUST) have been recently upgraded with silicon detectors which use resistive charge-splitting to measure the position of charged particles. This is in addition to the measurement of the total energy that these charged particles deposit in the detector. The upgrade results in increased angular resolution with a much smaller number of signals than silicon strip detectors with similar resolution. Here, a novel method of calibration has been used to demonstrate the efficacy of these position-sensitive detectors, and the detectors have been used in a commissioning experiment.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2023

Identifying native point defect configurations in α-alumina

Intimately intertwined atomic and electronic structures of point defects govern diffusion-limited corrosion and underpin the operation of optoelectronic devices. For some materials, complex energy landscapes containing metastable defect configurations challenge first-principles modeling efforts. Here, we thoroughly reevaluate native point defect geometries for the illustrative case of α-Al 2 O 3 by comparing three methods for sampling candidate geometries in density functional theory calculations: displacing atoms near a naively placed defect, initializing interstitials at high-symmetry points of a Voronoi decomposition, and Bayesian optimization. We find symmetry-breaking distortions for oxygen vacancies in some charge states, and we identify several distinct oxygen split-interstitial geometries that help explain literature discrepancies involving this defect. We also report a surprising and, to our knowledge, previously unknown trigonal geometry favored by aluminum interstitials in some charge states. Importantly, these new configurations may have transformative impacts on our understanding of defect migration pathways in aluminum-oxide scales protecting metal alloys from corrosion. Overall, the Voronoi scheme appears most effective for sampling candidate interstitial sites because it always succeeded in finding the lowest-energy geometry identified in this study, although no approach found every metastable configuration. Finally, we show that the position of defect levels within the band gap can depend strongly on the defect geometry, underscoring the need to conduct careful searches for ground-state geometries in defect calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1984

Chemical application of diffusion quantum Monte Carlo

The diffusion quantum Monte Carlo (QMC) method gives a stochastic solution to the Schroedinger equation. This approach is receiving increasing attention in chemical applications as a result of its high accuracy. However, reducing statistical uncertainty remains a priority because chemical effects are often obtained as small differences of large numbers. As an example, the single-triplet splitting of the energy of the methylene molecule CH sub 2 is given. The QMC algorithm was implemented on the CYBER 205, first as a direct transcription of the algorithm running on the VAX 11/780, and second by explicitly writing vector code for all loops longer than a crossover length C. The speed of the codes relative to one another as a function of C, and relative to the VAX, are discussed. The computational time dependence obtained versus the number of basis functions is discussed and this is compared with that obtained from traditional quantum chemistry codes and that obtained from traditional computer architectures.

Reynolds, P. J.↗

NASA NTRS1987

Three-dimensional unsteady Euler solutions for propfans and counter-rotating propfans in transonic flow

An Euler code designed for computing the unsteady, three-dimensional, transonic flow about single-rotating and counter-rotating propfans using dynamic blocked-grids is presented. The algorithm is a finite volume, flux-split, upwind, implicit scheme and solves the equations which have been written in a time-dependent curvilinear coordinate system. Relative motion of the blades for counter-rotating configurations is handled by requiring that grid lines be aligned after each discrete rotation of fore and aft rotor grid blocks. The method by which information is passed across block interfaces, as well as how downstream characteristic outflow boundary conditions which enforce simple radial equilibrium are implemented, is discussed. Comparisons of computed flow-field parameters and propfan performance with experimental data indicate good overall agreement between predictions and measurements.

Whitfield, D. L.↗

NASA NTRS1973

A Method of Monitoring Non-Resonant Raman Lidar Returns During Daylight Hours

A method for using non-resonant Raman lidar during daytime hours is presented. A dual (or quadruple) arrangement of electro-optical channels is set up in the lidar receiver. Each channel contains its own optical spectral response as determined by an appropriate set of interference filters. The incoming signal is split and fed into two channels, filtered (simultaneously) and then inserted into the input terminals of an operational amplifier. The difference signal is selected and amplified whereas the common portion of the signal is highly attenuated (Common Mode Re.jection Ratio = 75-80db). For the present hardware, daylight intensities 10(exp 2)-10(exp 3) times Raman signals are received simultaneously. Mean values of intensities of daylight signals which give rise to a constant differences merely shift the value of the measured variable at the ground. These mean value signals have their common part highly attenuated. The differences of the fluctuation portions of the daylight signal represents a noise signal. Identical (space-time) optical paths for each component of the daylight signal produce the same fluctuational spectra in each channel. Hence noise levels, much less than signal levels, are anticipated. First a N2-H2 ) signal pair is differenced. Then an N2-O2 pair is differenced and used for instrumental normalization purposes. Signal levels to be encountered are quite well known as a result of prior field work. Improved optical design of laser system can reduce mean value of daylight intensities by 10(exp 2). Thus the mean daylight intensity can be set equal to the expected signal levels by more advanced optical design. Performance characteristics of photomultiplier pairs are discussed in this connection.

Cooney, John↗

DOE OSTI2024

Solar Thermochemical Carbon Dioxide Splitting Using Ceria and Iron Aluminate Foam Devices and Simulation of a Plant System for Demonstration

An international research project has been undertaken to integrate a unique solar thermal processing reactor system with ceria and iron aluminate as active redox materials for CO2 splitting. Experimental investigations for CO2 splitting were conducted using a solar simulator and tube furnace at Niigata University, followed by demonstrations using a high-flux solar furnace (HFSF) at the National Renewable Energy Laboratory (NREL) in Golden, CO. Each experimental setup consisted of foam devices composed of reticulated porous ceramic (RPC). The RPC has a full ceria or iron aluminate body. It fabricated using the replica method and subjected to a two-step redox reaction, which iteratively separated a stream of CO2 into O2 and CO. Reactivity was evaluated using CO production per mass of the reactive material. The tubular furnace yielded a CO production of 6.41 mL/g at a reduction temperature of 1600degrees C, showcasing a higher CO production rate and total amount than those obtained from experiments conducted with solar simulators and solar furnace setups. For iron aluminate RPC, the productivity was measured as 3.57 mL/g using HFSF at a reduction temperature of 1450degrees C. These results are somewhat higher than those of the previous experiment at lower reduction temperatures of 1400degrees C-1500degrees C. Additionally, the production of CO in the case of ceria RPC was compared with the steady flow model simulation, which assumed chemical equilibrium at various levels of oxygen partial pressure during the reduction process. On the basis of these results, this study proposes a solar fuel system with an open receiver that uses a high-temperature heat transfer fluid.

carbon dioxide thermochemical splitting↗

NASA NTRS1985

Implicit finite-difference methods for the Euler equations

The present paper is concerned with two-dimensional Euler equations and with schemes which are in use of the time of this writing. Most of the development presented carries over directly to three dimensions. The characteristics of the two-dimensional Euler equations in Cartesian coordinates are considered along with generalized curvilinear coordinate transformations, metric relations, invariants of the transformation, flux Jacobian matrices and eigensystems, numerical algorithms, flux split algorithms, implicit and explicit nonlinear control (smoothing), upwind differencing in supersonic regions, unsteady and steady-state computation, the diagonal form of implicit algorithm, metric differencing and invariants, boundary conditions, geometry and mesh generation, and sample solutions.

Pulliam, T. H.↗

DOE OSTI2024

Highly confined epsilon-near-zero and surface phonon polaritons in SrTiO 3 membranes

Recent theoretical studies have suggested that transition metal perovskite oxide membranes can enable surface phonon polaritons in the infrared range with low loss and much stronger subwavelength confinement than bulk crystals. Such modes, however, have not been experimentally observed so far. Here, using a combination of far-field Fourier-transform infrared (FTIR) spectroscopy and near-field synchrotron infrared nanospectroscopy (SINS) imaging, we study the phonon polaritons in a 100 nm thick freestanding crystalline membrane of SrTiO 3 transferred on metallic and dielectric substrates. We observe a symmetric-antisymmetric mode splitting giving rise to epsilon-near-zero and Berreman modes as well as highly confined (by a factor of 10) propagating phonon polaritons, both of which result from the deep-subwavelength thickness of the membranes. Theoretical modeling based on the analytical finite-dipole model and numerical finite-difference methods fully corroborate the experimental results. Our work reveals the potential of oxide membranes as a promising platform for infrared photonics and polaritonics.

36 MATERIALS SCIENCE↗

NASA NTRS1974

Lunar surface dynamics: Some general conclusions and new results from Apollo 16 and 17

Exposure ages of Apollo 17 rocks as measured by tracks and the Kr-Kr rare gas method are reported. Concordant ages of 22 - or + 1 million year (my) are obtained for the station 6 boulder sample 76315. This value is interpreted as the time when the station 6 boulder was emplaced in its present position. Reasonable agreement is also obtained by the two methods for another station 6 boulder, sample 76015. Discordant ages (respectively 5 and 28 my by the track and rare gas methods) are obtained for the station 7 boulder sample, 77135, indicating that the boulder was emplaced at least 5 my ago. The 72 my exposure age of 75035, in general agreement with previous measurements of approximately 85 my for another Camelot boulder, may well date the formation of Camelot. Rock 76015 was split and one surface exposed to the sky through a very small solid angle.

Crozaz, G.↗

NASA NTRS1976

Thermodynamic analysis of alternate energy carriers, hydrogen and chemical heat pipes

The paper discusses the production concept and efficiency of two new energy transmission and storage media intended to overcome the disadvantages of electricity as an overall energy carrier. These media are hydrogen produced by water-splitting and the chemical heat pipe. Hydrogen can be transported or stored, and burned as energy is needed, forming only water and thus obviating pollution problems. The chemical heat pipe envisions a system in which heat is stored as the heat of reaction in chemical species. The thermodynamic analysis of these two methods is discussed in terms of first-law and second-law efficiency. It is concluded that chemical heat pipes offer large advantages over thermochemical hydrogen generation schemes on a first-law efficiency basis except for the degradation of thermal energy in temperature thus providing a source of low-temperature (800 K) heat for process heat applications. On a second-law efficiency basis, hydrogen schemes are superior in that the amount of available work is greater as compared to chemical heat pipes.

Cox, K. E.↗

DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗

DOE OSTI2021

Practical Effects of Integrating Temperature with Strang Split Reactions

For astrophysical reacting flows, operator splitting is commonly used to couple hydrodynamics and reactions. Each process operates independent of one another, but by staggering the updates in a symmetric fashion (via Strang splitting) second order accuracy in time can be achieved. However, approximations are often made to the reacting system, including the choice of whether or not to integrate temperature with the species. Here we demonstrate through a simple convergence test that integrating an energy equation together with reactions achieves the best convergence when modeling reactive flows with Strang splitting. Additionally, second order convergence cannot be achieved without integrating an energy or temperature equation.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Expansion dynamics and chemistry evolution in ultrafast laser filament produced plasmas

Laser ablation in conjunction with optical emission spectroscopy is a potential non-contact, standoff detection method for all elements in the periodic table and certain isotopes such as radionuclides. Currently, significant development efforts are on-going to use ultrafast laser filaments for remote sensing. The application of filaments is of particular interest in extending the range of stand-off capability associated with elemental and isotopic detection via laser-induced breakdown spectroscopy. In this study, we characterize the expansion dynamics and chemical evolution of filament-produced uranium (U) plasmas. Laser filaments are generated in the laboratory by loosely focusing 35 femtosecond (fs), 6 milli Joule (mJ) pulses in air. Two-dimensional spectrally-integrated and time-resolved imaging was performed to study hydrodynamics and evolution of U atomic and UO molecular emission in filament-produced U plasmas. Our results highlight that filament ablation of U plasmas gives a cylindrical plasma plume morphology with an appearance of plume splitting into slow and fast moving components at later times of its evolution. Emission from the slow-moving region shows no distinct spectral features (i.e. broadband-like) and is contributed in part by nanoparticles generated during ultrafast laser ablation. Additionally, we find U atoms and U oxide molecules (i.e. UO, UxOy) co-exist in the filament produced plasma, which can be attributed to the generation of low-temperature plasma conditions during filament ablation

(140.3440) Laser-induced breakdown, (280.5395) Pla↗

DOE OSTI2025

Material Fracturing and Failure Simulation Datasets

Fracturing is a fundamental physics phenomena with broad relevance across multiple domains, ranging from infrastructure integrity, aerospace durability, reservoir production, and seismic events. We present a diverse dataset of simulated fracture evolution and material failure generated from two numerical solvers: the phase-field method and the combined finite-discrete element method (FDEM). These solvers differ in formulation, physical fidelity, and computational efficiency. The dataset includes five materials: PBX, anisotropic shale, tungsten, aluminum, and steel. For each, phase-field simulations span 400,000 cases: 200,000 under uniaxial tension and 200,000 under biaxial tension. The computationally expensive FDEM simulations include 90,000 split evenly among PBX, shale, and tungsten under uniaxial loading. All simulations begin with randomized initial fracture patterns. Each entry includes temporal data capturing fracture propagation dynamics. This comprehensive dataset is designed to support the development of foundational or surrogate machine learning approaches for predicting material failure. While no such models are introduced here, the dataset lays a robust foundation for advancing future research and innovation in these areas.

36 MATERIALS SCIENCE↗