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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

A note on windowing for the waveform relaxation

The technique of windowing has been often used in the implementation of the waveform relaxations for solving ODE's or time dependent PDE's. Its efficiency depends upon problem stiffness and operator splitting. Using model problems, the estimates for window length and convergence rate are derived. The electiveness of windowing is then investigated for non-stiff and stiff cases respectively. lt concludes that for the former, windowing is highly recommended when a large discrepancy exists between the convergence rate on a time interval and the ones on its subintervals. For the latter, windowing does not provide any computational advantage if machine features are disregarded. The discussion is supported by experimental results.

Zhang, Hong↗
NASA NTRS1993

The piecewise-linear predictor-corrector code - A Lagrangian-remap method for astrophysical flows

We describe a time-explicit finite-difference algorithm for solving the nonlinear fluid equations. The method is similar to existing Eulerian schemes in its use of operator-splitting and artificial viscosity, except that we solve the Lagrangian equations of motion with a predictor-corrector and then remap onto a fixed Eulerian grid. The remap is formulated to eliminate errors associated with coordinate singularities, with a general prescription for remaps of arbitrary order. We perform a comprehensive series of tests on standard problems. Self-convergence tests show that the code has a second-order rate of convergence in smooth, two-dimensional flow, with pressure forces, gravity, and curvilinear geometry included. While not as accurate on idealized problems as high-order Riemann-solving schemes, the predictor-corrector Lagrangian-remap code has great flexibility for application to a variety of astrophysical problems.

Lufkin, Eric A.↗
NASA NTRS1994

The least-squares finite element method for low-mach-number compressible viscous flows

The present paper reports the development of the Least-Squares Finite Element Method (LSFEM) for simulating compressible viscous flows at low Mach numbers in which the incompressible flows pose as an extreme. Conventional approach requires special treatments for low-speed flows calculations: finite difference and finite volume methods are based on the use of the staggered grid or the preconditioning technique; and, finite element methods rely on the mixed method and the operator-splitting method. In this paper, however, we show that such difficulty does not exist for the LSFEM and no special treatment is needed. The LSFEM always leads to a symmetric, positive-definite matrix through which the compressible flow equations can be effectively solved. Two numerical examples are included to demonstrate the method: first, driven cavity flows at various Reynolds numbers; and, buoyancy-driven flows with significant density variation. Both examples are calculated by using full compressible flow equations.

Yu, Sheng-Tao↗
NASA NTRS1994

Multi-partitioning for ADI-schemes on message passing architectures

A kind of discrete-operator splitting called Alternating Direction Implicit (ADI) has been found to be useful in simulating fluid flow problems. In particular, it is being used to study the effects of hot exhaust jets from high performance aircraft on landing surfaces. Decomposition techniques that minimize load imbalance and message-passing frequency are described. Three strategies that are investigated for implementing the NAS Scalar Penta-diagonal Parallel Benchmark (SP) are transposition, pipelined Gaussian elimination, and multipartitioning. The multipartitioning strategy, which was used on Ethernet, was found to be the most efficient, although it was considered only a moderate success because of Ethernet's limited communication properties. The efficiency derived largely from the coarse granularity of the strategy, which reduced latencies and allowed overlap of communication and computation.

Vanderwijngaart, Rob F.↗
NASA NTRS1993

A numerical procedure for analysis of finite rate reacting flows

Combustion processes in rocket propulsion systems are characterized by the existence of multiple, vastly differing time and length scales, as well as flow-speeds at wide variation of Mach numbers. The chemical kinetics processes in the highly active reaction zone are characterized by much smaller scales compared to fluid convective and diffusive time scales. An operator splitting procedure for transient finite rate chemistry problems has been developed using a pressure based method, which can be applied to all speed flows without difficulties. The splitting of chemical kinetics terms formed the fluid-mechanical terms of the species equation ameliorated the difficulties associated with the disparate time scales and stiffness in the set of equations which describes highly exothermic combustion. A combined efficient ordinary differential equations (ODE) solver was used to integrate the effective chemical source terms over the residence time at each grid cell. One and two dimensional reacting flow situations were carried out to demonstrate and verify the current procedure. Different chemical kinetics with different degrees of nonlinearity have also been incorporated to test the robustness and generality of the proposed method.

Shang, H. M.↗
NASA NTRS1995

Convergence properties of finite-difference hydrodynamics schemes in the presence of shocks

We investigate the asymptotic convergence of finite-difference schemes for the Euler equations when the limiting solution contains shocks. The Lax-Wendroff theorem guarantees that certain conservative schemes converge to correct, physically valid solutions. We focus on two one-dimensional operator-split schemes with explicit artificial-viscosity terms. One, an internal-energy scheme, does not satisfy the assumptions of Lax-Wendroff; the other, a conservative total-energy scheme, does. With viscous lengths chosen proportional to the grid size, we find that both schemes converge to their zero-grid-size limits at the theoretically expected rate, but only the conversative scheme converges toward correct solutions of the inviscid fluid equations. We show that the difference in their behaviors results directly from the presence of shocks in the limiting solution. Empirically, we find that when the viscous lenghts tend toward zero more slowly than the grid size, however the nonconservative scheme also converges toward correct solutions. We characterize the asymptotic behavior of the total-energy scheme in a particular problem in which a shock forms. As the grid is refined, a Cauchy error approaches the expected rate of change slowly. We show that the changes in the artificial viscosity alter the diffusion of small-amplitude waves. The differences associated with such waves make the dominant contribution to the Cauchy error. We formulate an analytic model to relate the rate of approach to the effect of varying diffusion in waves and find quantitative agreement with our numerical results.

Kimoto, Paul A.↗
NASA NTRS1995

Nozzle flow with vibrational nonequilibrium

This research concerns the modeling and numerical solutions of the coupled system of compressible Navier-Stokes equations in cylindrical coordinates under conditions of equilibrium and nonequilibrium thermodynamics. The problem considered was the modeling of a high temperature diatomic gas N2 flowing through a converging-diverging high expansion nozzle. The problem was modeled in two ways. The first model uses a single temperature with variable specific heats as functions of this temperature. For the second model we assume that the various degrees of freedom all have a Boltzmann distribution and that there is a continuous redistribution of energy among the various degrees of freedom as the gas passes through the nozzle. Each degree of freedom is assumed to have its own temperature and, consequently, each system state can be characterized by these temperatures. This suggests that formulation of a second model with a vibrational degree of freedom along with a rotational-translation degree of freedom, each degree of freedom having its own temperature. Initially the vibrational degree of freedom is excited by heating the gas to a high temperature. As the high temperature gas passes through the nozzle throat there is a sudden drop in temperature along with a relaxation time for the vibrational degree of freedom to achieve equilibrium with the rotational-translation degree of freedom. That is, we assume that the temperature change upon passing through the throat is so great that the changes in the vibrational degree of freedom occur at a much slower pace and consequently lags behind the rotational-translational energy changes. This lag results in a finite relaxation time. In this context the term nonequilibrium is used to denote the fact that the energy content of the various degrees of freedom are characterized by two temperatures. We neglect any chemical reactions which could also add nonequilibrium effects. We develop the energy equations for the nonequilibrium model from first principles. The resulting equations, which model the nozzle flow, can be expressed in various forms. In most forms the resulting equations are coupled systems of nonlinear partial differential equations subject to certain boundary conditions. To solve the resulting coupled system of nonlinear partial differential equations, several numerical techniques were investigated: (1) the explicit MacCormack method, (2) the explicit-implicit MacCormack method, (3) the method of operator splitting, (4) factorization schemes, and (5) the Steger-Warming scheme.

Heinbockel, J. H.↗
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↗
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.↗
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.↗