Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “line relaxation”

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 181 records · Page 10

Radiative acceleration in outflows from broad absorption line quasi-stellar objects. 2: Wind models

We investigate the dynamics of radiatively driven broad absorption-line (BAL) outflows in quasi-stellar objects (QSOs) by developing radial and time-independent numerical models. Two limits are explored. The first assumes that the absorbing matter is not forced to comove with the substrate, which provides pressure confinement. This assumption allows us to explore in detail a case in which the acceleration is entirely due to radiation pressure. Using the parameters inferred from observations, we find that under these conditions radiative acceleration (mainly due to resonance line scattering) can readily accelerate the flow to the observed velocities. An important feature of the noncoupled flow is that the line profiles tend to stay relatively flat throughout the velocity interval covered by the line. We discuss how relaxing the assumptions of radial symmetry and time independent may help to explain the structures observed in BALs. In the second class of models, the absorbing flow is assumed to be completely coupled to the substrate in which it is embedded. Aside from being more plausible physically, these models produce line profiles that trail off at higher velocities, a behavior observed in some BALs. We show that, even if the substrate is massless, we have to assume a starting radius very close to the inferred radius of the broad emission-line region (approximately 0.1 pc) in order to obtain a significant contribution from radiative acceleration, given a typical active galactic nucleus (AGN) spectrum. The reason is that the energy input needed to pressurize the substrate, allowing the flow to become supersonic and to retain a reasonable ionization equilibrium, at the same time contributes appreciably to the acceleration. A way to relax the small starting radius constraint is to use a softer ionizing spectrum.

Arav, Nahum↗

Upwind relaxation algorithms for the Navier Stokes equations

The development of upwind relaxation algorithms for obtaining efficient steady-state solutions to the compressible Navier-Stokes equations is described. The method is second-order accurate spatially and naturally disipative, using third-order flux splitting of the pressure and convective terms and second-order central differencing for shear and heat flux terms. A line Gauss-Seidel relaxation approach, shown to be unconditionally stable for model convection and diffusion equations, is used. The algorithm is demonstrated for several flows using the thin-layer form of the equations, including the problem of shock-induced separation over a flat plate.

Thomas, J. L.↗

Upwind relaxation algorithms for Euler/Navier-Stokes equations

A description of and results from a solution algorithm for the compressible Navier-Stokes equations are presented. The main features of the algorithm are second or third order accurate upwind discretization of the convection and pressure derivatives and a relaxation scheme for the unfactored implicit backward Euler time method, implemented in a finite-volume formulation. Upwind methods were successfully used to obtain solutions to the Euler equations for flows with strong shock waves. The particular upwind method being used is based on the flux vector splitting technique developed by Van Leer and both second and third order accurate discretizations were developed. Currently, the most widely used implicit solution technique for the Navier-Stokes equations use approximate factorization (AF) methods to treat multidimensional problems. The time integration scheme being used in the present algorithm corresponds to a line Gauss-Seidel relaxation method. This method produces good convergence rates for steady-state flows, and most of the algorithm was vectorized on the NASA Langley VPS 32 computer. The Navier-Stokes algorithm was tested for several two-dimensional flow problems. Solutions for the problems gave excellent results. The presented effort is directed toward the extension of the scheme to the full three-dimensional Navier-Stokes equations.

Thomas, J. L.↗

Comparative Hydrophobic Core Dynamics Between Wild‐Type Amyloid‐β Fibrils, Glutamate‐3 Truncation, and Serine‐8 Phosphorylation

Abstract Post‐translational modifications (PTMs) of amyloid‐β (Aβ) species are implicated in the modulation of overall toxicities and aggregation propensities. We investigated the internal dynamics in the hydrophobic core of the truncated ΔE3 mutant fibrils of Aβ 1–40 and compared them with prior and new data for wild‐type fibrils as well as with phosphorylated S8 fibrils. Deuteron static solid‐state NMR techniques, spanning line‐shape analysis, longitudinal relaxation, and chemical exchange saturation transfer methods, were employed to assess the rotameric jumps of several methyl‐bearing and aromatic groups in the core of the fibrils. Taken together, the results indicate the rather significant influence of the PTMs on the hydrophobic core dynamics, which propagates far beyond the local site of the chemical modification. The phosphorylated S8 fibrils display an overall rigidifying of the core based on the higher activation barriers of motions than the wild‐type fibrils, whereas the ΔE3 fibrils induce a broader variety of changes, some of which are thermodynamic in nature rather than the kinetic ones.

Vugmeyster, Liliya↗

Efficient Nonthermal Ion and Electron Acceleration Enabled by the Flux-Rope Kink Instability in 3D Nonrelativistic Magnetic Reconnection

The relaxation of field-line tension during magnetic reconnection gives rise to a universal Fermi acceleration process involving the curvature drift of particles. However, the efficiency of this mechanism is limited by the trapping of energetic particles within flux ropes. Using 3D fully kinetic simulations, we demonstrate that the flux-rope kink instability leads to strong field-line chaos in weak-guide-field regimes where the Fermi mechanism is most efficient, thus allowing particles to transport out of flux ropes and undergo further acceleration. As a consequence, both ions and electrons develop clear power-law energy spectra that contain a significant fraction of the released energy. The low-energy bounds are determined by the injection physics, while the high-energy cutoffs are limited only by the system size. Furthermore, these results have strong relevance to observations of nonthermal particle acceleration in space and astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗

Implicit time-marching solution of the Navier-Stokes equations for thrust reversing and thrust vectoring nozzle flows

An implicit finite volume method is investigated for the solution of the compressible Navier-Stokes equations for flows within thrust reversing and thrust vectoring nozzles. Thrust reversing nozzles typically have sharp corners, and the rapid expansion and large turning angles near these corners are shown to cause unacceptable time step restrictions when conventional approximate factorization methods are used. In this investigation these limitations are overcome by using second-order upwind differencing and line Gauss-Siedel relaxation. This method is implemented with a zonal mesh so that flows through complex nozzle geometries may be efficiently calculated. Results are presented for five nozzle configurations including two with time varying geometries. Three cases are compared with available experimental data and the results are generally acceptable.

Imlay, S. T.↗

Magnetic energy dissipation in force-free jets

It is shown that a magnetic pressure-dominated, supersonic jet which expands or contracts in response to variations in the confining external pressure can dissipate magnetic energy through field-line reconnection as it relaxes to a minimum-energy configuration. In order for a continuous dissipation to occur, the effective reconnection time must be a fraction of the expansion time. The dissipation rate for the axisymmetric minimum-energy field configuration is analytically derived. The results indicate that the field relaxation process could be a viable mechanism for powering the synchrotron emission in extragalactic jets if the reconnection time is substantially shorter than the nominal resistive tearing time in the jet.

Choudhuri, Arnab Rai↗

A time accurate finite volume method for propulsion chamber flows

An implicit three-dimensional time-accurate method for propulsion chamber flows is proposed which uses line Gauss-Seidel relaxation and multiple axial sweeps for the convergence of each time step. The general time-integration algorithm employed includes such schemes as the Euler implicit method. The results of spatial and temporal accuracy tests reveal that Roe's (1981) flux difference splitting provides excellent tracking of acoustic wave speeds. In comparison with other methods, no low mean flow Mach number convergence limitation or Courant number stabilization restriction is observed.

Beddini, R. A.↗

Application of multi-grid methods for solving the Navier-Stokes equations

This paper presents the application of a class of multi-grid methods to the solution of the Navier-Stokes equations for two-dimensional laminar flow problems. The methods consists of combining the full approximation scheme-full multi-grid technique (FAS-FMG) with point-, line- or plane-relaxation routines for solving the Navier-Stokes equations in primitive variables. The performance of the multi-grid methods is compared to those of several single-grid methods. The results show that much faster convergence can be procured through the use of the multi-grid approach than through the various suggestions for improving single-grid methods. The importance of the choice of relaxation scheme for the multi-grid method is illustrated.

Demuren, A. O.↗

Application of multi-grid methods for solving the Navier-Stokes equations

The application of a class of multi-grid methods to the solution of the Navier-Stokes equations for two-dimensional laminar flow problems is discussed. The methods consist of combining the full approximation scheme-full multi-grid technique (FAS-FMG) with point-, line-, or plane-relaxation routines for solving the Navier-Stokes equations in primitive variables. The performance of the multi-grid methods is compared to that of several single-grid methods. The results show that much faster convergence can be procured through the use of the multi-grid approach than through the various suggestions for improving single-grid methods. The importance of the choice of relaxation scheme for the multi-grid method is illustrated.

Demuren, A. O.↗

Generalized conjugate-gradient methods for the Navier-Stokes equations

A generalized conjugate-gradient method is used to solve the two-dimensional, compressible Navier-Stokes equations of fluid flow. The equations are discretized with an implicit, upwind finite-volume formulation. Preconditioning techniques are incorporated into the new solver to accelerate convergence of the overall iterative method. The superiority of the new solver is demonstrated by comparisons with a conventional line Gauss-Siedel Relaxation solver. Computational test results for transonic flow (trailing edge flow in a transonic turbine cascade) and hypersonic flow (M = 6.0 shock-on-shock phenoena on a cylindrical leading edge) are presented. When applied to the transonic cascade case, the new solver is 4.4 times faster in terms of number of iterations and 3.1 times faster in terms of CPU time than the Relaxation solver. For the hypersonic shock case, the new solver is 3.0 times faster in terms of number of iterations and 2.2 times faster in terms of CPU time than the Relaxation solver.

Ajmani, Kumud↗

Preconditioned conjugate gradient methods for the Navier-Stokes equations

A preconditioned Krylov subspace method (GMRES) is used to solve the linear systems of equations formed at each time-integration step of the unsteady, two-dimensional, compressible Navier-Stokes equations of fluid flow. The Navier-Stokes equations are cast in an implicit, upwind finite-volume, flux-split formulation. Several preconditioning techniques are investigated to enhance the efficiency and convergence rate of the implicit solver based on the GMRES algorithm. The superiority of the new solver is established by comparisons with a conventional implicit solver, namely line Gauss-Seidel relaxation (LGSR). Computational test results for low-speed (incompressible flow over a backward-facing step at Mach 0.1), transonic flow (trailing edge flow in a transonic turbine cascade), and hypersonic flow (shock-on-shock interactions on a cylindrical leading edge at Mach 6.0) are presented. For the Mach 0.1 case, overall speedup factors of up to 17 (in terms of time-steps) and 15 (in terms of CPU time on a CRAY-YMP/8) are found in favor of the preconditioned GMRES solver, when compared with the LGSR solver. The corresponding speedup factors for the transonic flow case are 17 and 23, respectively. The hypersonic flow case shows slightly lower speedup factors of 9 and 13, respectively. The study of preconditioners conducted in this research reveals that a new LUSGS-type preconditioner is much more efficient than a conventional incomplete LU-type preconditioner.

Ajmani, Kumud↗

Time integration algorithms for the two-dimensional Euler equations on unstructured meshes

Explicit and implicit time integration algorithms for the two-dimensional Euler equations on unstructured grids are presented. Both cell-centered and cell-vertex finite volume upwind schemes utilizing Roe's approximate Riemann solver are developed. For the cell-vertex scheme, a four-stage Runge-Kutta time integration, a fourstage Runge-Kutta time integration with implicit residual averaging, a point Jacobi method, a symmetric point Gauss-Seidel method and two methods utilizing preconditioned sparse matrix solvers are presented. For the cell-centered scheme, a Runge-Kutta scheme, an implicit tridiagonal relaxation scheme modeled after line Gauss-Seidel, a fully implicit lower-upper (LU) decomposition, and a hybrid scheme utilizing both Runge-Kutta and LU methods are presented. A reverse Cuthill-McKee renumbering scheme is employed for the direct solver to decrease CPU time by reducing the fill of the Jacobian matrix. A comparison of the various time integration schemes is made for both first-order and higher order accurate solutions using several mesh sizes, higher order accuracy is achieved by using multidimensional monotone linear reconstruction procedures. The results obtained for a transonic flow over a circular arc suggest that the preconditioned sparse matrix solvers perform better than the other methods as the number of elements in the mesh increases.

Slack, David C.↗

Velocity measurements in the plume of an arcjet engine

A nonintrusive technique has been used to conduct a radial survey in the flow field of an arcjet engine plume. The technique measures the Doppler shift of an optically thin line resulting from recombination and relaxation processes in the high Mach number stream, in order to determine flow velocities. Atom temperature can also be calculated from the same Doppler-broadened line widths, when these shifts are measured with a scanning Fabry-Perot spectrometer whose design is presented in detail.

Pivirotto, T. J.↗

Air- and Self-Broadened Half Widths, Pressure-Induced Shifts, and Line Mixing in the Nu(sub 2) Band of (12)CH4

Lorentz self- and air-broadened half width and pressure-induced shift coefficients and their dependences on temperature have been measured from laboratory absorption spectra for nearly 130 transitions in the nu(sub 2) band of (12)CH4. In addition line mixing coefficients (using the relaxation matrix element formalism) for both self- and airbroadening were experimentally determined for the first time for a small number of transitions in this band. Accurate line positions and absolute line intensities were also determined. These parameters were obtained by analyzing high-resolution (approx. 0.003 to 0.01 per cm) laboratory spectra of high-purity natural CH4 and air-broadened CH4 recorded at temperatures between 226 and 297 K using the McMath-Pierce Fourier transform spectrometer (FTS) located at the National Solar Observatory on Kitt Peak, Arizona. A multispectrum nonlinear least squares technique was used to fit short (5-15 per cm) spectral intervals in 24-29 spectra simultaneously. Parameters were determined for nu(sub 2) transitions up to J" = 16. The variations of the measured broadening and shift parameters with the rotational quantum number index and tetrahedral symmetry species are examined. The present results are also compared with previous measurements available in the literature.

Smith, M. A. H.↗

An efficient iteration strategy for the solution of the Euler equations

A line Gauss-Seidel (LGS) relaxation algorithm in conjunction with a one-parameter family of upwind discretizations of the Euler equations in two-dimensions is described. The basic algorithm has the property that convergence to the steady-state is quadratic for fully supersonic flows and linear otherwise. This is in contrast to the block ADI methods (either central or upwind differenced) and the upwind biased relaxation schemes, all of which converge linearly, independent of the flow regime. Moreover, the algorithm presented here is easily enhanced to detect regions of subsonic flow embedded in supersonic flow. This allows marching by lines in the supersonic regions, converging each line quadratically, and iterating in the subsonic regions, thus yielding a very efficient iteration strategy. Numerical results are presented for two-dimensional supersonic and transonic flows containing both oblique and normal shock waves which confirm the efficiency of the iteration strategy.

Walters, R. W.↗

Efficient solutions to the Euler equations for supersonic flow with embedded subsonic regions

A line Gauss-Seidel (LGS) relaxation algorithm in conjunction with a one-parameter family of upwind discretizations of the Euler equations in two dimensions is described. Convergence of the basic algorithm to the steady state is quadratic for fully supersonic flows and is linear for other flows. This is in contrast to the block alternating direction implicit methods (either central or upwind differenced) and the upwind biased relaxation schemes, all of which converge linearly, independent of the flow regime. Moreover, the algorithm presented herein is easily coupled with methods to detect regions of subsonic flow embedded in supersonic flow. This allows marching by lines in the supersonic regions, converging each line quadratically, and iterating in the subsonic regions, and yields a very efficient iteration strategy. Numerical results are presented for two-dimensional supersonic and transonic flows containing oblique and normal shock waves which confirm the efficiency of the iteration strategy.

Walters, Robert W.↗

Spectral Line Parameters Including Temperature Dependences of Self- and Air-Broadening in the 2 (left arrow) 0 Band of CO at 2.3 micrometers

Temperature dependences of pressure-broadened half-width and pressure-induced shift coefficients along with accurate positions and intensities have been determined for transitions in the 2<--0 band of C-12 O-16 from analyzing high-resolution and high signal-to-noise spectra recorded with two different Fourier transform spectrometers. A total of 28 spectra, 16 self-broadened and 12 air-broadened, recorded using high- purity (greater than or equal to 99.5% C-12-enriched) CO samples and CO diluted with dry air(research grade) at different temperatures and pressures, were analyzed simultaneously to maximize the accuracy of the retrieved parameters. The sample temperatures ranged from 150 to 298K and the total pressures varied between 5 and 700 Torr. A multispectrum nonlinear least squares spectrum fitting technique was used to adjust the rovibrational constants (G, B, D, etc.) and intensity parameters (including Herman-Wallis coefficients), rather than determining individual line positions and intensities. Self-and air-broadened Lorentz half-width coefficients, their temperature dependence exponents, self- and air-pressure-induced shift coefficients, their temperature dependences, self- and air-line mixing coefficients, their temperature dependences and speed dependence have been retrieved from the analysis. Speed-dependent line shapes with line mixing employing off-diagonal relaxation matrix element formalism were needed to minimize the fit residuals. This study presents a precise and complete set of spectral line parameters that consistently reproduce the spectrum of carbon monoxide over terrestrial atmospheric conditions.

Devi, V. Malathy↗