SEARCH · Engineering Papers

Results for “linear equation systems”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,675 records · Page 93

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Unraveling Electronic and Vibrational Coherences Following a Charge Transfer Process in a Photosystem II Reaction Center

A reaction center is a unique biological system that performs the initial charge separation within a Photosystem II (PSII) multiunit enzyme, which eventually drives the catalytic water-splitting in plants and algae. The possible role of quantum coherences coinciding with the energy and charge transfer processes in PSII reaction center is one of the active areas of research. Here, we study these quantum coherences by using a numerically exact method on an excitonic dimer model, including linear vibronic coupling and employing optimal parameters from experimental two-dimensional coherent spectroscopic measurements. This enables us to precisely capture the excitonic interaction between pigments and the dissipation of the energy from electronic and charge-transfer (CT) states to the protein environment. We employ the time nonlocal (TNL) quantum master equation to calculate the population dynamics, which yields numerically reliable results. The calculated results show that, due to the strong dissipation, the lifetime of electronic coherence is too short to have direct participation in the charge transfer processes. However, there are long-lived vibrational coherences present in the system at frequencies close to the excitionic energy gap. These are strongly coupled with the electronic coherences, which makes the detection of the electronic coherences with conventional techniques very challenging. Additionally, we unravel the strong excitonic interaction of radical pair (PD1 and PD2) in the reaction center, which results in a long-lived electronic coherence of >100 fs, even at room temperature. Our work provide important physical insight to the charge separation process in PSII reaction center, which may be helpful for better understanding of photophysical processes in other natural and artificial light-harvesting systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

NASA NTRS2001

Thermal Analysis of a Carbon Fiber Rope Barrier for Use in the Reusable Solid Rocket Motor Nozzle Joint-2

This study provides for development and verification of analysis methods used to assess performance of a carbon fiber rope (CFR) thermal barrier system that is currently being qualified for use in Reusable Solid Rocket Motor (RSRM) nozzle joint-2. Modeled geometry for flow calculations considers the joint to be vented with the porous CFR barriers placed in the "open' assembly gap. Model development is based on a 1-D volume filling approach where flow resistances (assembly gap and CFRs) are defined by serially connected internal flow and the porous media "Darcy" relationships. Combustion gas flow rates are computed using the volume filling code by assuming a lumped distribution total joint fill volume on a per linear circumferential inch basis. Gas compressibility, friction and heat transfer are included in the modeling. Gas-to-wall heat transfer is simulated by concurrent solution of the compressible flow equations and a large thermal 2-D finite element (FE) conduction grid. The derived numerical technique loosely couples the FE conduction matrix with the compressible gas flow equations, Free constants that appear in the governing equations are calibrated by parametric model comparison to hot fire subscale test results. The calibrated model is then used to make full-scale motor predictions using RSRM aft dome environments. Model results indicate that CFR thermal barrier systems will provide a thermally benign and controlled pressurization environment for the RSRM nozzle joint-2 primary seal activation.

Clayton, J. Louie↗

NASA NTRS2002

Thermal Analysis of a Carbon Fiber Rope Barrier for Use in the Reusable Solid Rocket Motor Nozzle Joint-2

This study provides development and verification of analysis methods used to assess performance of a carbon fiber rope (CFR) thermal barrier system that is currently being qualified for use in Reusable Solid Rocket Motor (RSRM) nozzle joint-2. Modeled geometry for flow calculations considers the joint to be vented with the porous CFR barriers placed in the 'open' assembly gap. Model development is based on a 1-D volume filling approach where flow resistances (assembly gap and CFRs) are defined by serially connected internal flow and the porous media 'Darcy' relationships. Combustion gas flow rates are computed using the volume filling code by assuming a lumped distribution total joint fill volume on a per linear circumferential inch basis. Gas compressibility, friction and heat transfer are included in the modeling. Gas-to-wall heat transfer is simulated by concurrent solution of the compressible flow equations and a large thermal 2-D finite element (FE) conduction grid. The derived numerical technique loosely couples the FE conduction matrix with the compressible gas flow equations. Free constants that appear in the governing equations are calibrated by parametric model comparison to hot fire subscale test results. The calibrated model is then used to make full-scale motor predictions using RSRM aft dome environments. Model results indicate that CFR thermal barrier systems will provide a thermally benign and controlled pressurization environment for the RSRM nozzle joint-2 primary seal activation.

Clayton, J. Louie↗

DOE OSTI2022

Combining DPG in space with DPG time-marching scheme for the transient advection–reaction equation

In this article, we present a general methodology to combine the Discontinuous Petrov–Galerkin (DPG) method in space and time in the context of methods of lines for transient advection–reaction problems. We first introduce a semidiscretization in space with a DPG method redefining the ideas of optimal testing and practicality of the method in this context. Then, we apply the recently developed DPG-based time-marching scheme, which is of exponential-type, to the resulting system of Ordinary Differential Equations (ODEs). Further, we also discuss how to efficiently compute the action of the exponential of the matrix coming from the space semidiscretization without assembling the full matrix. Finally, we verify the proposed method for 1D+time advection–reaction problems showing optimal convergence rates for smooth solutions and more stable results for linear conservation laws comparing to the classical exponential integrators.

42 ENGINEERING↗

DOE OSTI2022

Combining DPG in space with DPG time-marching scheme for the transient advection-reaction equation

In this article, we present a general methodology to combine the Discontinuous PetrovGalerkin (DPG) method in space and time in the context of methods of lines for transient advection-reaction problems. We first introduce a semidiscretization in space with a DPG method redefining the ideas of optimal testing and practicality of the method in this context. Then, we apply the recently developed DPG-based time-marching scheme, which is of exponential-type, to the resulting system of Ordinary Differential Equations (ODEs). We also discuss how to efficiently compute the action of the exponential of the matrix coming from the space semidiscretization without assembling the full matrix. Finally, we verify the proposed method for 1D+time advection-reaction problems showing optimal convergence rates for smooth solutions and more stable results for linear conservation laws comparing to the classical exponential integrators.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1971

A multifilter phase-lock loop.

The phase model for the generalized multifilter phase-lock loop (M PLL) is considered and state equations for this model are derived. A linear analysis is presented to aid in the preliminary design of an M PLL and to indicate the noise improvement over a conventional phase-lock loop (PLL). Performance characteristics are examined for an M PLL with low-pass and bandpass characteristics used in a specific FM communication system. Both single and double sinusoidal FM are used and a region of proper operation of the M PLL is determined in terms of modulation index and modulati ng frequency. These results are obtained from both analog and digital computer simulation of the nonlinear system.

Carden, F.↗

NASA NTRS1971

Verification of theory on weak turbulence relating to the sequence of diffuse plasma resonances in space.

An interpretation of the sequence of diffuse plasma resonances observed by space probes (Alouette 2 and ISIS-1 satellites) is developed in terms of wave-particle nonlinear interaction in a weakly turbulent plasma including the electrostatic electron cyclotron harmonic wave instability. The longest time duration of the center frequency of the diffuse plasma resonance is found to coincide with the most favorable condition for the electrostatic electron cyclotron harmonic wave instability which is obtained by solving the dispersion equation obtained for a linear approximation of the kinetic wave equation for the warm magnetoactive plasma. The electrostatic field due to the transmission of the intense rf pulse produces plasma turbulence involving nonlinear wave-wave interaction and temperature anisotropy which leads to instability. This instability supplies energy to the turbulence. The process can be thought of as a feedback system.

Oya, H.↗

NASA NTRS1982

An approximation theory for nonlinear partial differential equations with applications to identification and control

Approximation results from linear semigroup theory are used to develop a general framework for convergence of approximation schemes in parameter estimation and optimal control problems for nonlinear partial differential equations. These ideas are used to establish theoretical convergence results for parameter identification using modal (eigenfunction) approximation techniques. Results from numerical investigations of these schemes for both hyperbolic and parabolic systems are given.

Banks, H. T.↗

NASA NTRS1986

Numerical study of hydrogen-air supersonic combustion by using elliptic and parabolized equations

The two-dimensional Navier-Stokes and species continuity equations are used to investigate supersonic chemically reacting flow problems which are related to scramjet-engine configurations. A global two-step finite-rate chemistry model is employed to represent the hydrogen-air combustion in the flow. An algebraic turbulent model is adopted for turbulent flow calculations. The explicit unsplit MacCormack finite-difference algorithm is used to develop a computer program suitable for a vector processing computer. The computer program developed is then used to integrate the system of the governing equations in time until convergence is attained. The chemistry source terms in the species continuity equations are evaluated implicitly to alleviate stiffness associated with fast chemical reactions. The problems solved by the elliptic code are re-investigated by using a set of two-dimensional parabolized Navier-Stokes and species equations. A linearized fully-coupled fully-implicit finite difference algorithm is used to develop a second computer code which solves the governing equations by marching in spce rather than time, resulting in a considerable saving in computer resources. Results obtained by using the parabolized formulation are compared with the results obtained by using the fully-elliptic equations. The comparisons indicate fairly good agreement of the results of the two formulations.

Chitsomboon, T.↗

NASA NTRS1986

Parallel scheduling of recursively defined arrays

A new method of automatic generation of concurrent programs which constructs arrays defined by sets of recursive equations is described. It is assumed that the time of computation of an array element is a linear combination of its indices, and integer programming is used to seek a succession of hyperplanes along which array elements can be computed concurrently. The method can be used to schedule equations involving variable length dependency vectors and mutually recursive arrays. Portions of the work reported here have been implemented in the PS automatic program generation system.

Myers, T. J.↗

NASA NTRS1988

An approximation theory for the identification of nonlinear distributed parameter systems

An abstract approximation framework for the identification of nonlinear distributed parameter systems is developed. Inverse problems for nonlinear systems governed by strongly maximal monotone operators (satisfying a mild continuous dependence condition with respect to the unknown parameters to be identified) are treated. Convergence of Galerkin approximations and the corresponding solutions of finite dimensional approximating identification problems to a solution of the original finite dimensional identification problem is demonstrated using the theory of nonlinear evolution systems and a nonlinear analog of the Trotter-Kato approximation result for semigroups of bounded linear operators. The nonlinear theory developed here is shown to subsume an existing linear theory as a special case. It is also shown to be applicable to a broad class of nonlinear elliptic operators and the corresponding nonlinear parabolic partial differential equations to which they lead. An application of the theory to a quasilinear model for heat conduction or mass transfer is discussed.

Banks, H. T.↗

NASA NTRS1990

An approximation theory for the identification of nonlinear distributed parameter systems

An abstract approximation framework for the identification of nonlinear distributed parameter systems is developed. Inverse problems for nonlinear systems governed by strongly maximal monotone operators (satisfying a mild continuous dependence condition with respect to the unknown parameters to be identified) are treated. Convergence of Galerkin approximations and the corresponding solutions of finite dimensional approximating identification problems to a solution of the original finite dimensional identification problem is demonstrated using the theory of nonlinear evolution systems and a nonlinear analog of the Trotter-Kato appproximation result for semigroups of bounded linear operators. The nonlinear theory developed here is shown to subsume an existing linear theory as a special case. It is also shown to be applicable to a broad class of nonlinear elliptic operators and the corresponding nonlinear parabolic partial differential equations to which they lead. An application of the theory to a quasilinear model for heat conduction or mass transfer is discussed.

Banks, H. T.↗

NASA NTRS1991

Alleviation of whirl-flutter on a joined-wing tilt-rotor aircraft configuration using active controls

The feasibility of using active controls to delay the onset of whirl-flutter on a joined-wing tilt rotor aircraft was investigated. The CAMRAD/JA code was used to obtain a set of linear differential equations which describe the motion of the joined-wing tilt-rotor aircraft. The hub motions due to wing/body motion is a standard input to CAMRAD/JA and were obtained from a structural dynamics model of a representative joined-wing tilt-rotor aircraft. The CAMRAD/JA output, consisting of the open-loop system matrices, and the airframe free vibration motion were input to a separate program which performed the closed-loop, active control calculations. An eigenvalue analysis was performed to determine the flutter stability of both open- and closed-loop systems. Sensor models, based upon the feedback of pure state variables and based upon hub-mounted sensors, providing physically measurable accelerations, were evaluated. It was shown that the onset of tilt-rotor whirl-flutter could be delayed from 240 to above 270 knots by feeding back vertical and span-wise accelerations, measured at the rotor hub, to the longitudinal cyclic pitch. Time response calculations at a 270-knot cruise condition showed an active cyclic pitch control level of 0.009 deg, which equates to a very acceptable 9 pound active-control force applied at the rotor hub.

Vanaken, Johannes M.↗

NASA NTRS1999

A Parallel Compact Multi-Dimensional Numerical Algorithm with Aeroacoustics Applications

In this study we propose a novel method to parallelize high-order compact numerical algorithms for the solution of three-dimensional PDEs (Partial Differential Equations) in a space-time domain. For this numerical integration most of the computer time is spent in computation of spatial derivatives at each stage of the Runge-Kutta temporal update. The most efficient direct method to compute spatial derivatives on a serial computer is a version of Gaussian elimination for narrow linear banded systems known as the Thomas algorithm. In a straightforward pipelined implementation of the Thomas algorithm processors are idle due to the forward and backward recurrences of the Thomas algorithm. To utilize processors during this time, we propose to use them for either non-local data independent computations, solving lines in the next spatial direction, or local data-dependent computations by the Runge-Kutta method. To achieve this goal, control of processor communication and computations by a static schedule is adopted. Thus, our parallel code is driven by a communication and computation schedule instead of the usual "creative, programming" approach. The obtained parallelization speed-up of the novel algorithm is about twice as much as that for the standard pipelined algorithm and close to that for the explicit DRP algorithm.

Povitsky, Alex↗

NASA NTRS1974

A method for determining optimum phasing of a multiphase propulsion system for a single-stage vehicle with linearized inert weight

A general analytical treatment is presented of a single-stage vehicle with multiple propulsion phases. A closed-form solution for the cost and for the performance and a derivation of the optimal phasing of the propulsion are included. Linearized variations in the inert weight elements are included, and the function to be minimized can be selected. The derivation of optimal phasing results in a set of nonlinear algebraic equations for optimal fuel volumes, for which a solution method is outlined. Three specific example cases are analyzed: minimum gross lift-off weight, minimum inert weight, and a minimized general function for a two-phase vehicle. The results for the two-phase vehicle are applied to the dual-fuel rocket. Comparisons with single-fuel vehicles indicate that dual-fuel vehicles can have lower inert weight either by development of a dual-fuel engine or by parallel burning of separate engines from lift-off.

Martin, J. A.↗

DOE OSTI2021

Driven-dissipative phase transition in a Kerr oscillator: From semiclassical PT symmetry to quantum fluctuations

We study a minimal model that has a driven-dissipative quantum phase transition, namely a Kerr non-linear oscillator subject to driving and dissipation. Using mean-field theory, exact diagonalization, and the Keldysh formalism, we analyze the critical phenomena in this system, showing which aspects can be captured by each approach and how the approaches complement each other. Then critical scaling and finite-size scaling are calculated analytically using the quantum Langevin equation. Furthermore, the physics contained in this simple model is surprisingly rich: it includes a continuous phase transition, Z 2 symmetry breaking, PT symmetry, state squeezing, and critical fluctuations. Due to its simplicity and solvability, this model can serve as a paradigm for exploration of open quantum many-body physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1973

Perform - A performance optimizing computer program for dynamic systems subject to transient loadings

A description and applications of a computer capability for determining the ultimate optimal behavior of a dynamically loaded structural-mechanical system are presented. This capability provides characteristics of the theoretically best, or limiting, design concept according to response criteria dictated by design requirements. Equations of motion of the system in first or second order form include incompletely specified elements whose characteristics are determined in the optimization of one or more performance indices subject to the response criteria in the form of constraints. The system is subject to deterministic transient inputs, and the computer capability is designed to operate with a large linear programming on-the-shelf software package which performs the desired optimization. The report contains user-oriented program documentation in engineering, problem-oriented form. Applications cover a wide variety of dynamics problems including those associated with such diverse configurations as a missile-silo system, impacting freight cars, and an aircraft ride control system.

Pilkey, W. D.↗

NASA NTRS1979

A fast semi-implicit algorithm for problems of mixed type

Certain physical processes are modeled by partial differential equations which are parabolic over part of the domain and elliptic over the remainder. A family of semi-implicit algorithms which are well suited to initial-boundary value problems of this mixed type is discussed. One important feature of these algorithms is the use of an approximate inverse for the solution of the implicit linear system. A strong error analysis results in an estimate of the total error as a function of approximate inverse error e and time step h.

Frederickson, P. O.↗