Search NASA⌕ Search

SEARCH · Search NASA

Results for “Operator splitting”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Establishing the Role of Operating Potential and Mass Transfer in Multicarbon Product Generation for Photoelectrochemical CO 2 Reduction Cells Using a Cu Catalyst

There is increasing interest in the possibility of photoelectrochemical (PEC) reduction of CO 2 to C 2+ products; however, the criteria for maximizing PEC solar-to-C 2+ (STC 2+ ) rates are not well understood. We report here a continuum-scale model of PEC CO 2 reduction (CO 2 R) on Cu in 0.1 M CsHCO 3 and use it to optimize the design and operating conditions for generating C 2+ products. Furthermore, we demonstrate that the potential-dependent product distribution of CO 2 R on Cu requires operating near the potential that maximizes C 2+ generation rates ($V$ id ), unlike PEC water splitting, which desires operation at the maximum photocurrent density. Because of this requirement, the criterion for a high STC 2+ rate includes high-photocurrent semiconductors with photovoltages near $V$ id and low series resistance. The STC 2+ rate in these systems is enhanced by optimal CO 2 transport and exhibits low sensitivity to dirunal solar irradiance variations.

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

Energy Use in Quantum Data Centers: Scaling the Impact of Computer Architecture, Qubit Performance, Size, and Thermal Parameters

As quantum computers increase in size, the total energy used by a quantum data center, including the cooling, will become a greater concern. The cooling requirements of quantum computers, which operate at temperatures near absolute zero, are determined by computing system parameters, including the number and type of physical qubits, the packaging efficiency of the system, and the split between circuits operating at cryogenic temperatures and those operating at room temperature. When combined with thermal system parameters such as cooling efficiency and cryostat heat transfer, the total energy use can be determined using a first-principles energy model. These models show that cooling of quantum computers differs in two fundamental ways from conventional data centers: (1) the energy required for cooling is much greater than the energy required for computation, and (2) the cooling loads are sensitive to the computational architecture. The temperature requirements for different qubit types can change energy requirements by orders of magnitude. Power use and computational power, as quantified by quantum volume, are analytically correlated. Approaches are identified for minimizing energy use in integrated quantum systems relative to computational power. Furthermore, designing a sustainable quantum computer will require both efficient cooling and system design that minimizes cooling requirements.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1996

Exponential Formulae and Effective Operations

One of standard methods to predict the phenomena of squeezing consists in splitting the unitary evolution operator into the product of simpler operations. The technique, while mathematically general, is not so simple in applications and leaves some pragmatic problems open. We report an extended class of exponential formulae, which yield a quicker insight into the laboratory details for a class of squeezing operations, and moreover, can be alternatively used to programme different type of operations, as: (1) the free evolution inversion; and (2) the soft simulations of the sharp kicks (so that all abstract results involving the kicks of the oscillator potential, become realistic laboratory prescriptions).

Mielnik, Bogdan↗
DOE OSTI2023

Oxynitrides enabled photoelectrochemical water splitting with over 3,000 hrs stable operation in practical two-electrode configuration

Solar photoelectrochemical reactions have been considered one of the most promising paths for sustainable energy production. To date, however, there has been no demonstration of semiconductor photoelectrodes with long-term stable operation in a two-electrode configuration, which is required for any practical application. Herein, we demonstrate the stable operation of a photocathode comprising Si and GaN, the two most produced semiconductors in the world, for 3,000 hrs without any performance degradation in two-electrode configurations. Measurements in both three- and two-electrode configurations suggest that surfaces of the GaN nanowires on Si photocathode transform in situ into Ga-O-N that drastically enhances hydrogen evolution and remains stable for 3,000 hrs. First principles calculations further revealed that the in-situ Ga-O-N species exhibit atomic-scale surface metallization. This study overcomes the conventional dilemma between efficiency and stability imposed by extrinsic cocatalysts, offering a path for practical application of photoelectrochemical devices and systems for clean energy.

14 SOLAR ENERGY↗
DOE OSTI2021

Operational Limits of Redox Metal Oxides Performing Thermochemical Water Splitting

Solar thermochemical hydrogen production is an attractive technology that stores intermittent solar energy in the form of chemical bonds. Efficient operation requires the identification of a redox-active metal oxide (MO x ) material that can achieve high conversion of water to hydrogen at minimal energy input. Water splitting occurs by consecutive reduction and reoxidation reactions of MO x . MO x is reduced to MO x-δ and, in the second step, is reoxidized by water recovering the initial MO x and generate H 2 . The material must reduce at temperatures achievable in concentrated solar receiver/reactors, while maintaining a thermodynamic driving force to split water. At equilibrium, extent of reduction depends on temperature and oxygen partial pressure, and in this analysis, a set of thermodynamic properties, namely, enthalpy and entropy of oxygen vacancy formation, is sufficient to represent MO x . Herein, a method to easily classify materials based on these thermodynamic properties under any condition of oxygen partial pressure and temperature is presented. This method is based on fundamental thermodynamic principles and is applicable for any redox material with known thermodynamic properties. Despite the simplicity of the method, it is believed that this analysis will support future research in targeting thermodynamic properties of redox-active metal oxides.

water splitting↗
NASA NTRS2019

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Olsson & Oliger, Gerritsen & Olsson, and Yee et al. (15, 7, 37) is revisited. The objective of this paper is to prove for the first time that the entropy split scheme is an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for nonlinear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. The proof is to replace the spatial derivatives by summation-by-parts (SBP) difference operators in the entropy split form of the equations using the physical entropy of the Euler equations. The numerical boundary closure follows directly from the SBP operator. No additional numerical boundary procedure is required. In contrast, Tadmor-type entropy conserving schemes (31) using mathematical entropies and more recently in (35], do not naturally come with a numerical boundary closure and a generalized SBP operator has to be developed (18). Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grubber splitting.

Sjogreen, Bjorn↗
NASA NTRSDate not supplied

Demand Capacity Balancing at Vertiports for Initial Strategic Conflict Management of Urban Air Mobility Operations

Urban Air Mobility (UAM) is a new transportation concept that enables highly automated, cooperative, passenger or cargo-carrying air transportation services in and around urban areas. To achieve the high level of operational density and complexity desired by the UAM community, an airspace system that allows UAM operators to readily access and operate safely and efficiently in the airspace is needed. This airspace system will require air traffic management designed to reduce the risk of conflicts and loss of separation between UAM flights. In general, strategic conflict management is considered as the first layer of conflict management for safe flight operations to condition the traffic to reduce the need for airborne separation provision, the second layer of conflict management. Demand Capacity Balancing (DCB) is one of the concept components to achieve strategic conflict management. DCB strategically evaluates traffic demand and resource capacities to allow UAM operators to determine when, where and how they operate, while mitigating conflicting needs for airspace and vertiport capacity. DCB can be applied whenever UAM demand exceeds the capacity in airspace or at vertiports. As the UAM ecosystem evolves with advanced technologies and matured operational procedures, more complicated conflict management will likely be needed. In the current UAM ‘Concept of Operation (ConOps) 1.0’ operational stage defined by FAA, however, it will be meaningful to explore the demand capacity balancing at vertiports only, as an initial strategic conflict management approach for UAM operations because vertiport capacity seems to be a bottleneck of UAM traffic. For this research, we developed a demand-capacity imbalance detection and resolution service for UAM. This DCB service identifies the demand from operators and compares the demand to a given capacity at the shared resources (i.e., vertiports) over the upcoming time horizon which is divided into time bins having a constant interval. When a new flight plan is submitted, the algorithm embedded in the DCB service checks the available time bins based on the desired departure time and estimated arrival time at origin and destination vertiports, respectively. If the time bins for the originally desired times are already occupied by other flights (i.e., demand is at or above capacity), the algorithm finds the next available time bins for takeoff and landing and shifts the conflicting departure time to the earliest time that satisfies the capacity constraints at both origin and destination vertiports. The details of the algorithm will be described in the final manuscript. Figure 1 shows that the proposed DCB algorithm works well for a sample traffic scenario. In this example, a total of 144 flights, split between two operators, are planned over 2 hours, traveling 10 routes between five vertiports. In the heatmaps, the horizontal axis shows 12 time bins where each bin represents a 12-minute interval, and the vertical axis shows five vertiports. The number in each cell shows the number of operations, counting both departures and arrivals, at a specific vertiport in each time bin. For the given capacity of 2 operations/vertiport/bin, Figure 1 shows that the original demand sometimes exceeds the capacity, but the modified demand is reduced to the given capacity after resolving demand-capacity imbalances. When UAM flights are operated, it is expected that many practical issues would arise in the federated system architecture with multiple operators. UAM operators may experience a time synchronization issue due to communication delay between operator and vehicle. UAM vehicles would fly at different flight speeds, depending on vehicle models. Actual departure and arrival times can have large variations, compared to the schedule. The lead time from flight plan submission to desired departure time can vary by service type (e.g., regular shuttle service vs. on-demand service). Using the proposed DCB algorithm, we also investigated how the actual flight schedule and DCB performance are affected by these uncertainties such as unsynchronized times between operators, flight speed differences, lead time differences, and departure time errors. The final manuscript will include the background of this research work, the description of the DCB algorithm and its use cases with traffic scenarios. It will also provide the analytical results about the impact of various uncertainties that can occur in actual UAM operations on the DCB at vertiports, in terms of demand distribution changes, number of simultaneous operations, and delay propagation.

Urban Air Mobility↗
DOE OSTI2024

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

ZPR/ZPPR Critical Experiment Facilities and As-Built ZPR/ZPPR Models

Argonne National Laboratory (ANL) operated four split-table critical experiment facilities, ZPR 3, ZPR-6, ZPR-9 and ZPPR, between 1955 and 1990. ZPR is the abbreviation for Zero Power Reactor, and ZPPR was the abbreviation for Zero Power Plutonium Reactor which was changed to Zero Power Physics Reactor in the 1980s due to negative views on Plutonium breeding and proliferation concerns. ZPR-3, ZPR-6, ZPR-9 and ZPPR were four reactors in a series of zero power (or lower power reactors) designed at Argonne in the 1950s and 1960s.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

L'Arlesienne de ROOT

Over many years, ROOT users have repeatedly stumbled over—and loudly rediscovered—the infamous 1 GB limit on individual I/O operations, a constraint that somehow survived long past the era when anyone thought a gigabyte was “a lot.” As experiments embraced ever-larger objects and collections, this limit became an increasingly unavoidable rite of passage. This contribution recounts the sustained, multi-year quest by ROOT I/O developers to finally retire this relic, navigating a maze of legacy APIs, memory-management assumptions, and integer boundaries that seemed determined to preserve the status quo. We describe how internal interfaces were carefully modernized to introduce fully 64-bit–capable code paths without breaking the mountains of existing user code that would definitely have noticed. With the limit now lifted, ROOT can finally handle multi-gigabyte objects in a single read or write operation, even when splitting them into an RNTuple is not an option (we’re looking at you, large RooWorkspaces and giant histograms), liberating users from yet another “fun” debugging adventure and clearing the way for the massive analyses of the HL-LHC and beyond.

Canal, Philippe G. [Fermilab] (ORCID:0000000277487↗
DOE OSTI2026

L'Arlesienne de ROOT

Over many years, ROOT users have repeatedly stumbled over—and loudly rediscovered—the infamous 1 GB limit on individual I/O operations, a constraint that somehow survived long past the era when anyone thought a gigabyte was “a lot.” As experiments embraced ever-larger objects and collections, this limit became an increasingly unavoidable rite of passage. This contribution recounts the sustained, multi-year quest by ROOT I/O developers to finally retire this relic, navigating a maze of legacy APIs, memory-management assumptions, and integer boundaries that seemed determined to preserve the status quo. We describe how internal interfaces were carefully modernized to introduce fully 64-bit–capable code paths without breaking the mountains of existing user code that would definitely have noticed. With the limit now lifted, ROOT can finally handle multi-gigabyte objects in a single read or write operation, even when splitting them into an RNTuple is not an option (we’re looking at you, large RooWorkspaces and giant histograms), liberating users from yet another “fun” debugging adventure and clearing the way for the massive analyses of the HL-LHC and beyond.

Canal, Philippe G. [Fermilab] (ORCID:0000000277487↗
NASA NTRS1979

Hydrogen as the solar energy translator

Many concepts are being investigated to convert sunlight to workable energy forms with emphasis on electricity and thermal energy. The electrical alternatives include direct conversion of photons to electricity via photovoltaic solar cells and solar/thermal production of electricity via heat-energy cycles. Solar cells, when commercialized, are expected to have efficiencies of about 12 to 14 percent. The cells would be active about eight hours per day. However, solar-operated water-splitting process research, initiated through JPL, shows promise for direct production of hydrogen from sunlight with efficiencies of up to 35 to 40 percent. The hydrogen, a valuable commodity in itself, can also serve as a storable energy form, easily and efficiently converted to electricity by fuel cells and other advanced-technology devices on a 24-hour basis or on demand with an overall efficiency of 25 to 30 percent. Thus, hydrogen serves as the fundamental translator of energy from its solar form to electrical form more effectively, and possibly more efficiently, than direct conversion. Hydrogen also can produce other chemical energy forms using solar energy.

Kelley, J. H.↗
DOE OSTI2024

Speckle contrast from the split-and-delay unit with seeded X-ray pulses of the MID instrument at European XFEL

We report on the coherence properties and characteristics of the split- and-delay unit at the Materials Imaging and Dynamics instrument of the European XFEL under seeded-beam conditions. Our investigation focuses on the speckle contrast extracted from the scattering patterns from static scatterers and pulse splitting characteristics. Seeded-beam operation enabled a high throughput of the split-and-delay unit. We highlight the invaluable potential of the split-and-delay unit for experimental investigations for enhancing our understanding of ultrafast phenomena in molecular liquids, such as water and aqueous solutions.

Fuoss, Paul↗
DOE OSTI2023

Effect of Split-Injection Strategies on Engine Performance and Emissions under Cold-Start Operation

The recently concluded partnership for advancing combustion engines (PACE) was a US Department of Energy consortium involving multiple national laboratories focused on addressing key efficiency and emission barriers in light-duty engines. Generation of detailed experimental data and modeling capabilities to understand and predict cold-start behavior was a major pillar in this program. Cold-start, as defined by the time between first engine crank and three-way catalyst light-off, is responsible for a large percentage of NOx, unburned hydrocarbon, and particulate matter emissions in light-duty engines. Minimizing emissions during cold-start is a trade-off between achieving faster three-way catalyst light-off, and engine out emissions during that period. In this study, engine performance, emissions, and catalyst warmup potential were monitored while the engine was operated using a single direct injection (baseline case) as well as a two-way-equal-split direct injection strategy. These injection strategies were analyzed at a range of cold-start-operation relevant retarded spark timings of up to 25 degrees after top dead center of firing (dATDCf). A stoichiometric 2-bar NIMEP steady-state condition was used for all cases to simulate cold-start operation. Significant improvement in engine stability was observed with the two-way-split injection strategy at the retarded spark timings allowing for up to 2.5x increase in exhaust heat rate when engine operation is stability constrained. Similar fuel-loss-to-oil trends with exhaust heat rate were observed for both single and two-way-split injection strategies. However, the two-way split injection was observed to produce higher NOx emissions per unit exhaust heat rate. A single data point run with three-way-split direct injection at a very retarded spark-timing of 30 dATDCf pointed to further improvements in engine stability and reduction in fuel-loss-to-oil as compared to single injection strategy. Engine stability decreased as spark timing was initially retarded with a single injection but was observed to plateau and stabilize beyond spark timing of 10 dATDCf. Finally, for the two-way-split-injection strategy, retarding the start of injection (SOI) timing of the second injection led to a decrease in engine stability as well as an increase in soot emissions.

42 ENGINEERING↗
NASA NTRS1972

Lunar rock splitter/can sealer

A manually powered and operated force generating device fractures lunar material as received from the moon and seals cans containing lunar samples. The rock splitting and can sealing operations are carried out in a vacuum environment. Can sealing is a cold welding process that requires a pressure of 125,000 psi (generating 32,000 lb for the largest can used). The design and modification of the device and the operating experience are described.

Johnson, K. G.↗
NASA NTRS1993

Point form relativistic quantum mechanics and relativistic SU(6)

The point form is used as a framework for formulating a relativistic quantum mechanics, with the mass operator carrying the interactions of underlying constituents. A symplectic Lie algebra of mass operators is introduced from which a relativistic harmonic oscillator mass operator is formed. Mass splittings within the degenerate harmonic oscillator levels arise from relativistically invariant spin-spin, spin-orbit, and tensor mass operators. Internal flavor (and color) symmetries are introduced which make it possible to formulate a relativistic SU(6) model of baryons (and mesons). Careful attention is paid to the permutation symmetry properties of the hadronic wave functions, which are written as polynomials in Bargmann spaces.

Klink, W. H.↗
NASA NTRSDate not supplied

A New Approach for a Wider Class of Entropy Split Methods for Compressible Gas Dynamics and MHD

The high order entropy split methods of Sjögreen & Yee [1, 2] by entropy splitting of the compressible Euler (inviscid) flux derivatives for a thermally-perfect gas are based on Harten’s entropy function [3, 4, 5]. Their derivation takes advantage of the homogeneity property of Euler flux, symmetrizable Euler flux derivatives and energy-norm stability in conjunction with high order classical spatial central, DRP (dispersion relation-preserving) [6, 7, 8] or Padé (compact) spatial discretizations [9] with summation-by-parts (SBP) operators [10]. Our entropy split methods have been proven entropy conserving and stable [1, 11, 12]. Our proofs do not rely on a two-point numerical flux, but rather only a linear difference operator is required to derive these methods. To extend the entropy split method for the MHD, we used the Godunov symmetrizable non-conservative MHD form [12, 13, 14]. These high order entropy split methods not only preserve certain physical properties of the chosen governing equations but are also known to either improve numerical stability, and/or minimize aliasing errors in long time integration of turbulent flow computations without the aid of added numerical dissipation. In our previous published work, extensive error norm comparison with grid refinement was performed to show the high accuracy performance of these methods. These studies also showed how well the entropy split methods conserve the entropy, momentum and mass, and preserve the kinetic energy for long time integration of the various flows [1, 2, 12, 13, 14]. The objective of the present work is to use a new approach to obtain a wider class of entropy split methods consisting of a two-point numerical flux portion and a non-conservative portion in such a way that the homogeneity property of the compressible Euler flux is not required. For high order classical spatial central, DRP (dispersion relation-preserving) or Padé (compact) spatial discretizations, this new approach can be proven to be entropy conservative with conservative spatial dsicretizations while at the same time allowing a wider class of symmetrizable inviscid flux derivatives. We also use this generalization to derive an entropy split scheme that is entropy conserving for the equations of MHD without the homogeneity property using the Godunov symmetrizable ideal MHD formulation [15].

High Order Physical Preserving Methods↗