Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Flow Applications of the Least Squares Finite Element Method

The main thrust of the effort has been towards the development, analysis and implementation of the least-squares finite element method (LSFEM) for fluid dynamics and electromagnetics applications. In the past year, there were four major accomplishments: 1) special treatments in computational fluid dynamics and computational electromagnetics, such as upwinding, numerical dissipation, staggered grid, non-equal order elements, operator splitting and preconditioning, edge elements, and vector potential are unnecessary; 2) the analysis of the LSFEM for most partial differential equations can be based on the bounded inverse theorem; 3) the finite difference and finite volume algorithms solve only two Maxwell equations and ignore the divergence equations; and 4) the first numerical simulation of three-dimensional Marangoni-Benard convection was performed using the LSFEM.

Jiang, Bo-Nan↗
NASA NTRS2013

Modeling of Particle Acceleration at Multiple Shocks Via Diffusive Shock Acceleration: Preliminary Results

We present preliminary results from a model that diffusively accelerates particles at multiple shocks. Our basic approach is related to box models (Protheroe and Stanev, 1998; Moraal and Axford, 1983; Ball and Kirk, 1992; Drury et al., 1999) in which a distribution of particles is diffusively accelerated inside the box while simultaneously experiencing decompression through adiabatic expansion and losses from the convection and diffusion of particles outside the box (Melrose and Pope, 1993; Zank et al., 2000). We adiabatically decompress the accelerated particle distribution between each shock by either the method explored in Melrose and Pope (1993) and Pope and Melrose (1994) or by the approach set forth in Zank et al. (2000) where we solve the transport equation by a method analogous to operator splitting. The second method incorporates the additional loss terms of convection and diffusion and allows for the use of a variable time between shocks. We use a maximum injection energy (Emax) appropriate for quasi-parallel and quasi-perpendicular shocks (Zank et al., 2000, 2006; Dosch and Shalchi, 2010) and provide a preliminary application of the diffusive acceleration of particles by multiple shocks with frequencies appropriate for solar maximum (i.e., a non-Markovian process).

Parker, Linda Neergaard↗
NASA NTRS2013

Modeling of Particle Acceleration at Multiple Shocks via Diffusive Shock Acceleration: Preliminary Results

Successful forecasting of energetic particle events in space weather models require algorithms for correctly predicting the spectrum of ions accelerated from a background population of charged particles. We present preliminary results from a model that diffusively accelerates particles at multiple shocks. Our basic approach is related to box models in which a distribution of particles is diffusively accelerated inside the box while simultaneously experiencing decompression through adiabatic expansion and losses from the convection and diffusion of particles outside the box. We adiabatically decompress the accelerated particle distribution between each shock by either the method explored in Melrose and Pope (1993) and Pope and Melrose (1994) or by the approach set forth in Zank et al. (2000) where we solve the transport equation by a method analogous to operator splitting. The second method incorporates the additional loss terms of convection and diffusion and allows for the use of a variable time between shocks. We use a maximum injection energy (E(sub max)) appropriate for quasi-parallel and quasi-perpendicular shocks and provide a preliminary application of the diffusive acceleration of particles by multiple shocks with frequencies appropriate for solar maximum (i.e., a non-Markovian process).

Parker, L. Neergaard↗
NASA NTRSDate not supplied

Wear Test of the 12.5-kW Advanced Electric Propulsion System Engineering Test Unit Hall Thruster

This work presents a summary of the first wear test of the 12.5 kW Advanced Electric Propulsion System Engineering Test Unit 2 (AEPS ETU-2) thruster produced by Aerojet Rocketdyne. The ETU-2 Wear Test accumulated approximately 730 hours of operation split between the nominal 600 V/12.5 kW condition and the 300 V/6.25 kW condition previously identified as the worst-case erosion condition.

Jason D. Frieman↗
NASA NTRSDate not supplied

Discretization Error Estimation and Control for Farfield Acoustic Signatures

We investigate the utility of adjoint-based error estimates for sonic boom farfield simulations governed by solutions of the augmented Burgers’ equation. Solution of this nonlinear system uses operator splitting with a second-order finite volume discretization in space and second-order Runge-Kutta time marching, while the absorption and molecular relaxation are solved using second-order central differencing. The discretization error in selected ground sonic boom cost functionals is estimated using the method of adjoint-weighted residuals. Key elements of the implementation process are emphasized with details provided on the practical aspects as appliedto the sonic boom farfield propagation. We establish the accuracy of the adjoint solutions usingcomplex step and finite difference approaches, and examine the accuracy of the error estimates using analytical N-wave solutions. We then apply it to a pressure waveform corresponding to the X-59 research aircraft. The investigations demonstrate that the method of adjoint-weighted residuals accurately predicts the level of discretization error present in sonic boom farfield simulations while offering insight into which features of the near field signal are the primary drivers of ground noise metrics. The numerical results indicate that at sampling frequencies as low as50kHz, discretization error in the propagation is under 0.01 dB[A] for realistically complex examples.

CST↗
DOE OSTI2020

Multiphysics coupling plan

In this chapter, multiphysics coupling simulations of nuclear reactor systems are reviewed. On the one hand, the mainstream methods of multiphysics coupling are demonstrated. The fundamental theory and coupling scheme of operator splitting methods, Jacobian-free Newton–Krylov methods, and approximate block Newton methods are summarized and presented. On the other hand, some significant works on neutronic and thermal-hydraulic codes coupling, and large research project, including NURESAFE European project, Multiphysics Object-Oriented Simulation Environment plan, and Consortium for the Advanced Simulation of Light Water Reactors project, are reviewed. This chapter helps the readers have a better understanding on the worldwide current status of multiphysics coupling research.

Gui, Miao↗
NASA NTRSDate not supplied

Extended Wear Testing of the 12-kW Advanced Electric Propulsion System Engineering Test Unit Hall Thruster

This work presents a summary of the first wear test of the 12 kW Advanced Electric Propulsion System (AEPS) intended to isolate the impact of discharge current on component wear rates. Testing was led by Aerojet Rocketdyne and performed with the AEPS Engineering Test Unit 2 (ETU-2) thruster. In total, approximately 860 hours of operation were accumulated split between operating conditions of 600 V/9 kW, 600 V/11 kW, and 600 V/12 kW. Thruster performance and stability were invariant throughout the wear test for all thruster throttle conditions and shown to be equal to the values previously measured with ETU-2. Inner front pole cover erosion rates were found to be invariant to discharge current as the measurements at all three operating conditions were equal to within the empirical uncertainty. Outer front pole cover erosion rates were found to be equal for the 600 V/9 kW and 600 V/12 kW conditions with operation at 600 V/11 kW yielding elevated erosion rates. Overall, the data shows that the AEPS thruster design has a high probability of meeting lifetime and performance requirements as the design proceeds to flight production and qualification.

Jason D Frieman↗
NASA NTRSDate not supplied

Extended Wear Testing of the 12-kW Advanced Electric Propulsion System Engineering Test Unit Hall Thruster

This work presents a summary of the first wear test of the 12 kW Advanced Electric Propulsion System (AEPS) intended to isolate the impact of discharge current on component wear rates. Testing was led by Aerojet Rocketdyne and performed with the AEPS Engineering Test Unit 2 (ETU-2) thruster. In total, approximately 860 hours of operation were accumulated split between operating conditions of 600 V/9 kW, 600 V/11 kW, and 600 V/12 kW. Thruster performance and stability were invariant throughout the wear test for all thruster throttle conditions and shown to be equal to the values previously measured with ETU-2. Inner front pole cover erosion rates were found to be invariant to discharge current as the measurements at all three operating conditions were equal to within the empirical uncertainty. Outer front pole cover erosion rates were found to be equal for the 600 V/9 kW and 600 V/12 kW conditions with operation at 600 V/11 kW yielding rates that were comparable to those on the inner front pole cover. Overall, the data shows that the AEPS thruster design has a high probability of meeting lifetime and performance requirements as the design proceeds to flight production and qualification.

HERMeS↗
NASA NTRS1986

The effects of various implicit operators on a flux vector splitting method

Three different implicit operators in a numerical method for solving the two-dimensional steady Euler equations using flux vector splitting are investigated. These include the implementation of the true Jacobian matrices in the implicit part, the use of their splitting approximate form, and a modification of the implicit part of the scheme that uses the scalar form of the implicit part based on the spectral radii of the split Jacobian matrices. The three versions of the basic algorithm are compared with the results given by two common numerical methods using several two-dimensional test cases. Special attention is paid to the quality of results as well as computational efficiency and convergence properties.

Von Lavante, E.↗
DOE OSTI2021

Substation-Level Grid Topology Optimization Using Bus Splitting: Preprint

Operations of substation circuit breakers are of high significance for performing system maintenance and topology optimization. Bus splitting is one type of topology changes where the two bus-bars at a substation become electrically disconnected after certain actions of circuit breakers. As these events involve detailed substation modeling, they are not typically considered in power system routine operation and control. In this paper, an improved substation-level topology optimization is developed by expanding traditional line switching with breaker-level bus splitting, which can further reduce grid congestion and generation costs. A tight McCormick relaxation is proposed to reformulate the bi-linear terms in the resultant topology optimization model. Thus, a tractable mixed-integer linear program formulation is presented which can be efficiently solved for real-time control. Numerical studies on the IEEE 14-bus and 118-bus systems demonstrate the performance and economic benefits of the proposed topology optimization approach.

bus split↗
DOE OSTI2023

Protecting Customer Privacy Through Distributed Energy Resource Anonymization

Due to their stochastic nature, the increase of Renewable Energy Resources (RERs) as a primary source of energy for power grids creates challenges regarding the reliability and resilience of the system. In order to combat these obstacles, expansion of Distributed Energy Resources (DERs) and their participation in Demand Response (DR) programs is necessary. Widespread participation requires prioritizing customer privacy and addressing concerns that may arise regarding communication between DERs and the Grid Service Provider (GSP). This paper discusses the use of flow reservation resources to split the operating cycles of DER load profiles into unique phases. The splitting of phases increases anonymization of the DERs by making it more difficult to determine the individual characteristics of the device. We discuss an example of this using simulated DER load profile data and examine the resulting effectiveness by using a machine learning algorithm for classification, called Support Vector Machine (SVM).

Distributed Energy Resource, Anonymization, Renewa↗
NASA NTRS1934

Preliminary Wind-Tunnel and Flight Tests of a Balanced Split Flap, Special Report

One disadvantage that has been apparent in the operation of split flaps as used to date is the time and effort required to operate them. In this communication an investigation is being made of possible means for balancing them aerodynamically to make their operation easier. Several arrangements have been tested in the 7 by 210 foot wind tunnel, and the results of the wind-tunnel tests as well as preliminary flight tests on one of the more promising forms are given in this paper.

Weick, Fred E.↗
NASA NTRS1972

Computational efficiency achieved by time splitting of finite difference operators.

A technique is presented for computing multidimensional time-dependent flow fields that avoids much of the inefficiency typically found in finite difference calculations. The technique initially divides the flow field into regions, each containing a mesh of general quadrilateral cells chosen to provide spatial resolution of the local features of the flow. A finite difference operator of second order accuracy, consisting of a sequence of one-dimensional operators (each operating at near maximum Courant-Friedrich-Lewy number) is then constructed for each region. Numerical results illustrating the technique for inviscid flows about simple bodies that generate shock waves, embedded shock waves, and expansion fans are presented and compared with exact theory.

Maccormack, R. W.↗
DOE OSTI2021

Split representation of adaptively compressed polarizability operator

The polarizability operator plays a central role in density functional perturbation theory and other perturbative treatment of first principle electronic structure theories. The cost of computing the polarizability operator generally scales as O(Ne4) where Ne is the number of electrons in the system. The recently developed adaptively compressed polarizability operator (ACP) formulation [L. Lin, Z. Xu and L. Ying, Multiscale Model. Simul. 2017] reduces such complexity to O(Ne 3 ) in the context of phonon calculations with a large basis set for the first time, and demonstrates its effectiveness for model problems. In this paper, we improve the performance of the ACP formulation by splitting the polarizability into a near singular component that is statically compressed, and a smooth component that is adaptively compressed. The new split representation maintains the O(Ne 3 ) complexity, and accelerates nearly all components of the ACP formulation, including Chebyshev interpolation of energy levels, iterative solution of Sternheimer equations, and convergence of the Dyson equations. For simulation of real materials, we discuss how to incorporate nonlocal pseudopotentials and finite temperature effects. In this work, we demonstrate the effectiveness of our method using one-dimensional model problem in insulating and metallic regimes, as well as its accuracy for real molecules and solids.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Stochastic projective splitting

Here, we present a new, stochastic variant of the projective splitting (PS) family of algorithms for inclusion problems involving the sum of any finite number of maximal monotone operators. This new variant uses a stochastic oracle to evaluate one of the operators, which is assumed to be Lipschitz continuous, and (deterministic) resolvents to process the remaining operators. Our proposal is the first version of PS with such stochastic capabilities. We envision the primary application being machine learning (ML) problems, with the method’s stochastic features facilitating “mini-batch” sampling of datasets. Since it uses a monotone operator formulation, the method can handle not only Lipschitz-smooth loss minimization, but also min–max and noncooperative game formulations, with better convergence properties than the gradient descent-ascent methods commonly applied in such settings. The proposed method can handle any number of constraints and nonsmooth regularizers via projection and proximal operators. We prove almost-sure convergence of the iterates to a solution and a convergence rate result for the expected residual, and close with numerical experiments on a distributionally robust sparse logistic regression problem.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

The suppression of pulsar and gamma-ray burst annihilation lines by magnetic photon splitting

Neutron stars, relativistic and compact by nature, show great potential for the copious creation of electron-positron pairs in the magnetospheres; these rapidly cool, thermalize, and then annihilate. It is therefore expected that many neutron sources might display evidence of pair annihilation lines in the 400-500 keV range. It is shown that magnetic photon splitting, which operates effectively at these energies and in the enormous neutron star magnetic fields, can destroy an annihilation feature by absorbing line photons and reprocessing them to lower energies. In so doing, photon splitting creates a soft gamma-ray bump and a broad quasi-power-law contribution to the X-ray continuum, which is too flat to conflict with the observed X-ray paucity in gamma-ray bursts. The destruction of the line occurs in neutron stars with surface fields of 5 x 10 exp 12 G or maybe even less, depending on the size of the emission region.

Baring, Matthew G.↗
DOE OSTI2021

Design Principles for Efficient and Stable Water Splitting Photoelectrocatalysts

Photoelectrochemical water splitting is a promising avenue for sustainable production of hydrogen used in the chemical industry and hydrogen fuel cells. The basic components of most photoelectrochemical water splitting systems are semiconductor light absorbers coupled to electrocatalysts, which perform the desired chemical reactions. A critical challenge for the design of these systems is the lack of stability for the majority of desired semiconductors under operating water splitting conditions. One strategy to address this issue is to protect the semiconductor by covering it with a stabilizing insulator layer, creating a metal–insulator–semiconductor (MIS) architecture, which has demonstrated improved stability. In addition to enhanced stability, the insulator layer may significantly affect the electron and hole transfer, which governs the recombination rates. Furthermore, the insertion of an insulator layer leads to the introduction of additional insulator/electrocatalyst and insulator/semiconductor interfaces. These interfaces can impact the system’s performance significantly, and they need to be carefully engineered to optimize the efficiencies of MIS systems. In this Account, we describe our recent progress in shedding light on the critical role of the insulator and the interfaces on the performance of MIS systems. We discuss our findings by focusing on the concrete example of planar n-type Si protected by a HfO 2 insulator layer and coupled to a Ni or Ir electrocatalyst that performs the oxygen evolution reaction, one of the water splitting half-reactions. To improve our fundamental understanding of the insulator layer, we precisely control the HfO 2 insulator thickness using atomic layer deposition (ALD), and we perform a series of rigorous electrochemical experiments coupled with theory and modeling. We demonstrate that by tuning the insulator thickness, we can control the flux and recombination of photogenerated electrons and holes to optimize the generated photovoltage. Despite optimizing the thickness, we find that the maximum generated photovoltage in MIS systems is often significantly lower than the upper performance limit, i.e., there are additional losses in the system that could not be addressed by optimizing the insulator thickness. We identify the sources of these losses and describe strategies to minimize them by a combination of improving the semiconductor light absorption, removing nonidealities associated with interfacial defects, and finding alternative insulators with improved charge carrier selectivity. Finally, we quantify the improvements that can be obtained by implementing these specific strategies. Our collective work outlines strategies to analyze MIS systems, identify the sources of efficiency losses, and optimize the design to approach the fundamental performance limits. Furthermore, these general approaches are broadly applicable to photoelectrochemical materials that utilize sunlight to produce value-added chemicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗