Search NASA⌕ Search

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.

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

At least 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ∼1/Np with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly non-equilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ~1/$\sqrt{N_p}$ with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly nonequilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2019

High-Order Finite-Difference Nonlinear Filter Methods for Subsonic Turbulence Simulation with Stochastic Forcing

Numerical stability of high-order filter schemes developed by Yee & Sjogreen is tested on three-dimensional turbulence simulations with stochastic forcing and their performance is compared with that of TVD and WENO schemes. The best­ performing filter method employs an eighth-order central base scheme with the Kennedy & Gruber skew-symmetric splitting of the inviscid flux derivative, a wavelet-based local flow sensor, a nonlinear filter utilizing the dissipative portion of seventh-order \VENO scheme, and an explicit third - or fourth-order Runge­-Kutta time integration. We show that the filter scheme is more computational]y efficient and provides a wider spectral bandwidth compared to the seventh-order WENO scheme. The method also demonstrates robust long-time integration for moderately compressible turbulence. In contrast, the fifth- and seventh­ order WENO schemes show non-trivial evolution of the velocity and density power spectra. over a. few dozen dynamical times, where both TVD and filter schemes recover a so lid statistically stationary turbulent state

Kritsuk, Alexei G.↗
DOE OSTI2023

Novel Chalcopyrites for Advanced Photoelectrochemical Water Splitting

With the support of DoE’s EERE office, our team has established a unique tool-chest of capabilities, including theoretical modeling (Lawrence Livermore National Laboratory: LLNL), state-of-the-art synthesis (Hawaii Natural Energy Institute: HNEI, Stanford, and the National Renewable Energy Laboratory: NREL) and advanced materials and interfaces characterization (University of Nevada, Las Vegas: UNLV, and Lawrence Berkeley National Laboratory: LBNL), to accelerate the development of high efficiency and durable chalcopyrite materials for advanced photoelectrochemical (PEC) water splitting. Using this synergistic approach, we have successfully created new wide bandgap chalcopyrite photocathodes generating over 10 mA/cm 2 , developed innovative strategies to protect them from corrosion, and engineered novel integration methods to circumvent thin film materials mechanical, chemical and thermal incompatibility. In Task 1 “Modeling and synthesis of chalcopyrite photocathodes”, we expanded our library of wide bandgap chalcopyrites for PEC water splitting. With support from LLNL’s “Computational Materials Diagnostics and Optimization of PEC Devices”, LBNL’s “photophysical” and NREL’s “I-III-VI Compound Semiconductors for Water-Splitting” nodes, we investigated two new chalcopyrite candidates for PEC water splitting: Cu(In,Al)Se 2 and Cu(In,B)Se 2 . We also further developed ordered vacancy compounds, such as CuGa 3 Se 5 , with unprecedented durability during PEC waters splitting in acidic solutions. In Task 2 “Interfaces engineering for enhanced efficiency and durability”, we addressed both the non-ideal band-edge positions of chalcopyrites with respect to water redox potentials, as well as their chemical instability under PEC water splitting, with a buried-junctions approach. With help from NREL’s “High-Throughput Experimental Thin Film Combinatorial Capabilities” and “Corrosion Analysis of Materials” nodes, we engineered environmentally friendly n-type buffers, including Mn x Zn 1-x O, to adjust the chalcopyrite band-edge positions and achieved photovoltages as high as 925 mV. Also, we integrated non-precious catalytic-protecting layers, such as WO 3 , to enhance the water splitting long-term stability of chalcopyrite absorbers. Finally, in Task 3 “Hybrid photoelectrode device integration”, we proposed an innovative method to bond wide bandgap photocathodes onto narrow bandgap PV drivers at room temperature using conductive polymers. Our semi-monolithic approach addressed fundamental processing incompatibility issues, as both the photocathode and the PV driver are processed separately. Proof-of-concept whole-chalcopyrite tandems were obtained by consecutive exfoliation and transfer of fully integrated 1.85 eV CuGa 3 Se 5 and 1.13 eV CuInGaSe 2 stacks from their Mo/SLG substrates onto a new single FTO host substrate.

08 HYDROGEN↗
DOE OSTI2024

RESOLVING THE ELECTROCHEMICAL EQUATIONS OF A SOLID OXIDE FUEL CELL FOR USE IN TRANSIENT SIMULATION AND INTEGRATION INTO CYBER-PHYSICAL SYSTEMS

A major challenge with complex cyber-physical systems stems from long model computational time that creates a mismatch between the model system and the physical system. The numerical modeling of solid oxide fuel cells (SOFCs) presents particular challenges due to the highly coupled nature of the underlying equations and the multiphysics needed to fully resolve their behavior during a transient event. To this end current approaches revolve around splitting the computational efforts into resolving temperature effects and resolving electrochemical effects. Current methods employed for the transient simulation of an SOFC for implementation in the Hybrid Performance (HyPer) facility cyberphysical plant at the National Energy Technology Laboratory reveal a distinct need for accelerated results with a high degree of stability. To this aim, an investigation into the computational time for the code reveals that the underlying electrochemical algorithm takes an order of magnitude more time than its thermal counterpart and has a tendency to vary in terms of iteration time and as such a rework of the underlying system is proposed. The primary method for accelerated electrochemical algorithm solutions is to employ higher order root finding recipes for the resolution of the highly coupled electrochemical equations. This is done with the intention to reduce the overall number of subiterations necessary for resolving voltage, current density, and species concentration, properties of the fuel cell that are all directly coupled and require nested iterative approaches. The overall objective of this approach is an order of magnitude reduction in calculation time without sacrificing stability and increasing accuracy. Specific approaches involve using both bounded and unbounded techniques, such as the False Position method and the Secant method (or if applicable Newton-Raphson) respectively, the drawbacks being slower convergence for False Position and instability for the Secant or Newton-Raphson methods. Current preliminary results on simplified versions of the parent functions involved for electrochemical calculations indicate a reduction in computational steps by a factor of two for the secant method and a factor of three for Newton-Raphson. When implemented into new modified electrochemical algorithms, the results indicate a possible order of magnitude reduction in calculation time.

Arias, Jesus↗
DOE OSTI2021

Effect of Ligand Polarity on the Internal Dipoles and Ferroelectric Distortion in BaTiO 3 Nanocubes

Here, surface adsorbates and surrounding matrix species have been demonstrated to affect the properties of nanoscale ferroelectrics and nanoscale ferroelectric composites; potentially counteracting performance losses that can occur in small particle sizes. In this work, the effects of nonpolar oleic acid (OA) and polar tetrafluoroborate (BF 4 - ) ligand capping on the surface of various sizes of BaTiO 3 nanocubes have been investigated with combined neutron diffraction and neutron pair distribution function (PDF), density functional theory (DFT), and ab initio molecular dynamics (AIMD) methods. The low real space PDF region provides an unobstructed view of rhombohedral (split short and long) Ti-O distances in BaTiO 3 nanocubes, mimicking the well-established order-disorder local structure found in bulk BaTiO 3 . Interestingly, the intermediate-range order in nanocubes is found to be orthorhombic, rather than tetragonal. It is concluded that polar ligands adsorbed at BaTiO 3 surfaces stabilize the correlation length scale of local rhombohedral distortions in ferroelectric nanoparticles relative to nonpolar ligands.

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

Visibility of Hypovascularized Liver Tumors during Intra-Arterial Therapy Using Split-Bolus Single-Phase Cone Beam CT

Purpose: To validate a split-bolus contrast injection protocol for single-phase CBCT in terms of detectability of hypovascular liver tumors compared to digital subtraction angiography (DSA). Materials and Methods: In this retrospective, single-center study, 20 consecutive patients with in total 77 hypovascularized tumors referred for intra-arterial therapy received a split-bolus single-phase CBCT. Two readers rated the visibility of the target tumors scheduled for embolization in CBCT and DSA compared to the pre-interventional multiphasic CT or MRI used as reference on a 3-point scoring system (1 = optimal, 3 = not visible) and catheter-associated artifacts (1 = none, 3 = extended). SNR, CNR and contrast values were derived from 37 target tumors in CBCT and MRI. Statistical analysis included the kappa test to determine interrater reliability, the Friedman’s test for the inter-modality comparison evaluating tumor visibility in DSA and CBCT as well as for quantitative assessment. Post hoc analysis included the Wilcoxon signed-rank test. p values < 0.05 were considered significant. Results: Ninety percentage of target tumors were rated as visible in CBCT and 37.5% in DSA (p < 0.001). 70.1% of pre-interventionally detected hypovascularized tumors were depicted with CBCT and 31.2% by DSA (p < 0.001). 7.8% of known tumors were outside the FOV. Quantitative assessment showed higher image contrasts in CBCT (1.91 ± 7.01) compared to hepatobiliary-phase MRI (0.29 ± 0.14, p = 0.003) and to portal-venous (p.v.) MRI (0.31 ± 0.13, p < 0.001), but higher CNR for MRI (1.18 ± 0.80; 13.92 ± 15.82; 13.79 ± 6.65). Conclusion: In conclusion, the split-bolus single-phase CBCT detects significantly more hypovascularized liver tumors compared to DSA performed through the proper hepatic artery with high image contrasts. Level of Evidence: Level III, diagnostic study.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2025

Kβ X-ray Emission Spectra Analysis Using Bayesian Optimization

The Kβ X-ray emission spectrum of 3 d transition metals is rich with electronic and structural information due to strong exchange interactions with the valence shell of the metal, and has become crucial for understanding their spin and oxidation states. The spectrum is commonly treated using crystal-field multiplet theory, a semi-empirical theory that uses tunable parameters to control the strength of the effects present in X-ray emission spectroscopy (XES). However, determining the experimental values of these parameters remains a challenge. We present a methodology that applies Bayesian optimization to crystal-field multiplet theory to determine parameter values. The algorithm is tested on the X-ray emission spectra of a collection of Mn, Co, and Ni oxides. We are able to find optimal values for the four most impactful parameters: Slater−Condon reduction factors F dd , F pd , and G pd , and crystal field splitting 10 Dq . The algorithm produces significantly improved accuracy compared to current analysis methods, and probes interparameter dependencies by modeling the error landscape. This advancement enhances XES analysis by offering an approach of obtaining quantitative electronic structural information on 3 d transition metal valence shells, facilitating applications across various scientific fields.

Bayesian optimization↗
DOE OSTI2021

CATMoS: Collaborative Acute Toxicity Modeling Suite

Background: Humans are exposed to tens of thousands of chemical substances that need to be assessed for their potential toxicity. Acute systemic toxicity testing serves as the basis for regulatory hazard classification, labeling, and risk management. However, it is cost- and time-prohibitive to evaluate all new and existing chemicals using traditional rodent acute toxicity tests. In silico models built using existing data facilitate rapid acute toxicity predictions without using animals. Objectives: The U.S. Interagency Coordinating Committee on the Validation of Alternative Methods Acute Toxicity Workgroup organized an international collaboration to develop in silico models for predicting acute oral toxicity based on five different endpoints: LD50 value, U.S. Environmental Protection Agency hazard categories, Globally Harmonized System for Classification and Labelling hazard categories, very toxic chemicals (LD50 =50 mg/kg), and non-toxic chemicals (LD50 >2000 mg/kg). Methods: An acute oral toxicity data inventory for 11,992 chemicals was compiled, split into training and evaluation sets, and made available to 35 participating international research groups that submitted a total of 139 predictive models. Predictions that fell within the applicability domains of the submitted models were evaluated using external validation sets. These were then combined into consensus models to leverage strengths of individual approaches. Results: The resulting consensus predictions, which leverage the collective strengths of each individual model, form the Collaborative Acute Toxicity Modeling Suite (CATMoS). CATMoS demonstrated high performance in terms of accuracy and robustness when compared to in vivo results. Discussion: CATMoS is being evaluated by regulatory agencies for its utility and applicability as a potential replacement for in vivo rat acute oral toxicity studies. CATMoS predictions for over 800,000 chemicals have been made available via the NTP’s Integrated Chemical Environment. The models are also implemented in a free, standalone open-source tool, OPERA, which allows predictions of new and untested chemicals to be made.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2024

Integration of Concentrating Solar Power with High Temperature Electrolysis for Hydrogen Production: Preprint

Hydrogen (H2) has been identified as a leading sustainable contender to replace fossil fuels in transportation and electricity generation. H2 production can be achieved by concentrating solar thermal power (CSP) systems collecting thermal energy from the sun to various chemical processes for fuel production. Fuel production via solar thermal chemical processes integrated with CSP uses the full spectrum of sunlight compared with photovoltaic power conversion and stores solar energy directly and efficiently [1]. The solar fuel production can be realized by thermochemical processes (e.g., water splitting for H2 production, carbon dioxide reduction, or methane reforming) or thermal electrochemical methods (e.g., integration with solid oxide electrolysis cell). Technology development for CSP-integrated solar fuel production requires broad technological bases from solar energy collection to chemical energy conversion. H2 generated from renewable sources can be an energy carrier for a carbon-free economy. Integrating CSP with high temperature electrolysis (HTE) using solid oxide electrolysis cells (SOEC) provides a renewable path for H2 generation. The CSP-HTE integration approach provides the benefit of thermal energy storage (TES) for continuous operation, improved capacity, and SOEC life. H2 gas has low energy density for transportation, pipeline networks are expensive, and H2 liquefaction is energy intensive. However, an alternative method for H2 distribution is to use carbon dioxide (CO2) capture and liquid hydrocarbon synthesis to convert solar energy into liquid fuels that are compatible with the existing fossil fuel infrastructure.

concentrating solar thermal power↗
NASA NTRS1966

Orbital tube flare for tubing connectors

Marshall engineering personnel have developed and are now refining a new and unique method of mechanically flaring tubing. The system eliminates the use of the conventional split backing die which is used with most tube-flaring equipment. The Marshall-developed tool is adapted for use on standard Leonard 3-CP-HD flaring machines. The stationary split backing die is replaced in function with a die ring, whose i.d.-surface rolls or orbits around the o.d. of the tube flare zone. An i.d.-flaring cone moves axially and forces the metal against the rolling die to form the flare. Major advantages of this system are the elimination of flash, or split die marks, on the back angle of the flare, and improved flare-surface finishes, permitting more uniform sealing characteristics and improvements of the total geometric configuration of the flare as required by the MSFC specification MC-146. The Orbital Adapter system of tooling is currently being used to manufacture flared tube hardware for the Saturn V S-IO space launch vehicle at MSFC, for the Agena space launch vehicle (used for the Gemini project) at Lockheed, Sunnyvale, for LEM checkout equipment for Grumman, and for the Saturn S-IV-B space launch vehicle by Douglas Aircraft Company. Several other aerospace firms are evaluating the system for possible incorporation on their tube flaring requirements. Marshall engineers are now concentrating their efforts on refining and reducing the total cost of the overall system.

Tubing↗
NASA NTRS1973

Prospects for obtaining metallic hydrogen with spherical presses

Description of a split-sphere apparatus modified for use at low temperature and affording a possible method for compressing molecular hydrogen to a pressure in excess of 1 Mbar and for converting it to the metallic state. The construction costs of the apparatus are relatively low and the amount of liquid helium required for low-temperature operation is readily obtainable with modern liquefiers.

Spain, I. L.↗
NASA NTRS1980

Relativistic scattered wave calculations on UF6

Self-consistent Dirac-Slater multiple scattering calculations are presented for UF6. The results are compared critically to other relativistic calculations, showing that the results of all molecular orbital calculations are in qualitative agreement, as measured by energy levels, population analyses, and spin-orbit splittings. A detailed comparison is made to the relativistic X alpha(RX alpha) method of Wood and Boring, which also uses multiple scattering theory, but incorporates relativistic effects in a more approximate fashion. For the most part, the RX alpha results are in agreement with the present results.

Case, D. A.↗
NASA NTRS1986

An LU-SSOR scheme for the Euler and Navier-Stokes equations

A new multigrid relaxation scheme, lower-upper symmetric successive overrelaxation (LU-SSOR) is developed for the steady-state solution of the Euler and Navier-Stokes equations. The scheme, which is based on central differences, does not require flux splitting for approximate Newton iteration. Application to transonic flow shows that the new method is efficient and robust. The vectorizable LU-SSOR scheme needs only scalar diagonal inversions.

Yoon, Seokkwan↗
NASA NTRS1987

An LU-SSOR scheme for the Euler and Navier-Stokes equations

A new multigrid relaxation scheme, lower-upper symmetric successive overrelaxation (LU-SSOR) is developed for the steady-state solution of the Euler and Navier-Stokes equations. The scheme, which is based on central differences, does not require flux splitting for approximate Newton iteration. Application to transonic flow shows that the new method is efficient and robust. The vectorizable LU-SSOR scheme needs only scalar diagonal inversions.

Yoon, Seokkwan↗
NASA NTRS1989

Numerical solution of the incompressible Navier-Stokes equations for steady-state and time-dependent problems

An algorithm for the solution of the incompressible Navier-Stokes equations in three-dimensional generalized curvilinear coordinates is presented. The algorithm can be used to compute both steady-state and time-dependent flow problems. The algorithm is based on the method of artificial compressibility and uses a higher-order flux-difference splitting technique for the convective terms and a second-order central difference for the viscous terms. The steady-state solution of flow through a square duct with a 90 deg bend is computed and the results are compared with experimental data. Good agreement is observed. A comparison with an analytically known exact solution is then performed to verify the time accuracy of the algorithm. Finally, the flow through an artificial heart configuration with moving boundaries is calculated and presented.

Rogers, Stuart E.↗
NASA NTRS1990

A time-dependent dusty gas dynamic model of axisymmetric cometary jets

The present time-dependent, axisymmetric dusty gas dynamical model of inner cometary atmospheres solves the coupled and time-dependent equations of continuity, momentum, and energy for a gas-dust mixture between the surface of the nucleus and 100 km, using an axisymmetric 40 x 40 grid structure. A novel numerical method employing a second-order accurate Godunov-type scheme with dimensional splitting is used to solve the time-dependent pde system. It is established that a subsolar dust spike not predicted by previous calculations is generated by narrow axisymmetric jets, together with a jet cone whose opening angle depends on the jet length.

Korosmezey, A.↗
NASA NTRS1991

Steady and unsteady solutions of the incompressible Navier-Stokes equations

An algorithm for the solution of the incompressible Navier-Stokes equations in three-dimensional generalized curvilinear coordinates is presented. The algorithm can be used to compute both steady-state and time-dependent flow problems. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The equations are solved with a line-relaxation scheme that allows the use of very large pseudotime steps leading to fast convergence for steady-state problems as well as for the subiterations of time-dependent problems. The steady-state solution of flow through a square duct with a 90-deg bend is computed, and the results are compared with experimental data. Good agreement is observed. Computations of unsteady flow over a circular cylinder are presented and compared to other experimental and computational results. Finally, the flow through an artificial heart configuration with moving boundaries is calculated and presented.

Rogers, Stuart E.↗