Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear differential equations”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

Program for solution of ordinary differential equations

A program for the solution of linear and nonlinear first order ordinary differential equations is described and user instructions are included. The program contains a new integration algorithm for the solution of initial value problems which is particularly efficient for the solution of differential equations with a wide range of eigenvalues. The program in its present form handles up to ten state variables, but expansion to handle up to fifty state variables is being investigated.

Sloate, H.↗
NASA NTRS1978

A procedure for analyzing transonic flow over harmonically oscillating airfoils

Finite difference procedures were successfully used to solve the steady transonic flow about airfoils and appear to provide a practical means for calculating the corresponding unsteady flow. The purpose of the paper is to describe a finite difference procedure derived from the equations for the potential flow by assuming small perturbations and harmonic motion. The velocity potential is divided into steady and unsteady parts, and the resulting unsteady equation is linearized on the basis of small amplitudes of oscillation. The steady velocity potential, which must be calculated first, is described by the classical nonlinear transonic differential equation.

Weatherill, W. H.↗
NASA NTRS1983

Transient responses of phosphoric acid fuel cell power plant system

An analytical and computerized study of the steady state and transient response of a phosphoric acid fuel cell (PAFC) system was completed. Parametric studies and sensitivity analyses of the PAFC system's operation were accomplished. Four non-linear dynamic models of the fuel cell stack, reformer, shift converters, and heat exchangers were developed based on nonhomogeneous non-linear partial differential equations, which include the material, component, energy balance, and electrochemical kinetic features. Due to a lack of experimental data for the dynamic response of the components only the steady state results were compared with data from other sources, indicating reasonably good agreement. A steady state simulation of the entire system was developed using, nonlinear ordinary differential equations. The finite difference method and trial-and-error procedures were used to obtain a solution. Using the model, a PAFC system, that was developed under NASA Grant, NCC3-17, was improved through the optimization of the heat exchanger network. Three types of cooling configurations for cell plates were evaluated to obtain the best current density and temperature distributions. The steady state solutions were used as the initial conditions in the dynamic model. The transient response of a simplified PAFC system, which included all of the major components, subjected to a load change was obtained. Due to the length of the computation time for the transient response calculations, analysis on a real-time computer was not possible. A simulation of the real-time calculations was developed on a batch type computer. The transient response characteristics are needed for the optimization of the design and control of the whole PAFC system. All of the models, procedures and simulations were programmed in Fortran and run on IBM 370 computers at Cleveland State University and the NASA Lewis Research Center.

Lu, Cheng-Yi↗
NASA NTRS1985

Magnetic and thermodynamic structure of the solar corona during sunspot minimum

A model of the solar corona during sunspot minimum was constructed. It was suggested that pole-equator asymmetry is caused by interaction of global azimuthal electric currents in the corona and potential magnetic fields originating below the corona. The problem was reduced to a nonlinear ordinary differential equation for a structural function. Taking the difference in electron density between the equatorial plane and the polar direction from observations, the basic equation was solved and the global electric current distribution around the Sun along with the magnetic and thermodynamic structures of the solar corona was derived. The asymmetry between the two hemispheres, which can be due to the magnetic quadrupole, is considered in detail. It is shown that an additional quadrupole term in the magnetic field representation affects the position of polar coronal hole boundaries differently in the two hemispheres. As a result, the outflow from the two polar regions is found to be different.

Osherovich, V. A.↗
NASA NTRS1971

Estimation in nonlinear systems with transport delay.

The problem of estimation of state in nonlinear dynamical systems containing time delays is studied. The plant is specified by a set of nonlinear differential-difference equations. Observations are a nonlinear function of current and/or delayed states. Both contain additive disturbances. The criterion used for the optimal estimates is the integral of the weighted squared error. Using the theory of the calculus of variations, equations are developed for the estimation. They are first expressed in the form of a split boundary value problem, which is then converted to an initial value problem for on-line estimation. The result yields a sequential estimation scheme in which filtered and smoothed estimates are computed in a sequential manner. The applicability of the procedure is demonstrated by a practical example.

Stoller, R. L.↗
NASA NTRS1979

Analysis of turbulent free-jet hydrogen-air diffusion flames with finite chemical reaction rates

A numerical analysis is presented of the nonequilibrium flow field resulting from the turbulent mixing and combustion of an axisymmetric hydrogen jet in a supersonic parallel ambient air stream. The effective turbulent transport properties are determined by means of a two-equation model of turbulence. The finite-rate chemistry model considers eight elementary reactions among six chemical species: H, O, H2O, OH, O2 and H2. The governing set of nonlinear partial differential equations was solved by using an implicit finite-difference procedure. Radial distributions were obtained at two downstream locations for some important variables affecting the flow development, such as the turbulent kinetic energy and its dissipation rate. The results show that these variables attain their peak values on the axis of symmetry. The computed distribution of velocity, temperature, and mass fractions of the chemical species gives a complete description of the flow field. The numerical predictions were compared with two sets of experimental data. Good qualitative agreement was obtained.

Sislian, J. P.↗
NASA NTRS1989

Hypervelocity atmospheric flight: Real gas flow fields

Flight in the atmosphere is examined from the viewpoint of including real gas phenomena in the flow field about a vehicle flying at hypervelocity. That is to say, the flow field is subject not only to compressible phenomena, but is dominated by energetic phenomena. There are several significant features of such a flow field. Spatially, its composition can vary by both chemical and elemental species. The equations which describe the flow field include equations of state and mass, species, elemental, and electric charge continuity; momentum; and energy equations. These are nonlinear, coupled, partial differential equations that have been reduced to a relatively compact set of equations in a self-consistent manner (which allows mass addition at the surface at a rate comparable to the free-stream mass flux). The equations and their inputs allow for transport of these quantities relative to the mass-average behavior of the flow field. Thus transport of mass by chemical, thermal, pressure, and forced diffusion; transport of momentum by viscosity; and transport of energy by conduction, chemical considerations, viscosity, and radiative transfer are included. The last of these complicate the set of equations by making the energy equations a partial integrodifferential equation. Each phenomenon is considered and represented mathematically by one or more developments. The coefficients which pertain are both thermodynamically and chemically dependent. Solutions of the equations are presented and discussed in considerable detail, with emphasis on severe energetic flow fields. Hypervelocity flight in low-density environments where gaseous reactions proceed at finite rates chemical nonequilibrium is considered, and some illustrations are presented. Finally, flight where the flow field may be out of equilibrium, both chemically and thermodynamically, is presented briefly.

Howe, John T.↗
NASA NTRS1990

Hypervelocity atmospheric flight: Real gas flow fields

Flight in the atmosphere is examined from the viewpoint of including real gas phenomena in the flow field about a vehicle flying at hypervelocity. That is to say, the flow field is subject not only to compressible phenomena, but is dominated by energetic phenomena. There are several significant features of such a flow field. Spatially, its composition can vary by both chemical and elemental species. The equations which describe the flow field include equations of state and mass, species, elemental, and electric charge continuity; momentum; and energy equations. These are nonlinear, coupled, partial differential equations that were reduced to a relatively compact set of equations of a self-consistent manner (which allows mass addition at the surface at a rate comparable to the free-stream mass flux). The equations and their inputs allow for transport of these quantities relative to the mass-averaged behavior of the flow field. Thus transport of mass by chemical, thermal, pressure, and forced diffusion; transport of momentum by viscosity; and transport of energy by conduction, chemical considerations, viscosity, and radiative transfer are included. The last of these complicate the set of equations by making the energy equation a partial integrodifferential equation. Each phenomenon is considered and represented mathematically by one or more developments. The coefficients which pertain are both thermodynamically and chemically dependent. Solutions of the equations are presented and discussed in considerable detail, with emphasis on severe energetic flow fields. For hypervelocity flight in low-density environments where gaseous reactions proceed at finite rates, chemical nonequilibrium is considered and some illustrations are presented. Finally, flight where the flow field may be out of equilibrium, both chemically and thermodynamically, is presented briefly.

Howe, John T.↗
NASA NTRS2004

Rapid Switching and Modulation by use of Coupled VCSELs

Devices incorporating coupled vertical-cavity surface-emitting lasers (VCSELs) have been proposed as means to effect unprecedentedly rapid modulation and/or switching of light beams in advanced optical communication and information-processing systems. A basic device according to the proposal would comprise two or more VCSELs (1) positioned so close to each other that they are coupled and (2) connected to a current source that biases them at a steady current above the threshold current for lasing. The operation of VCSELs under the conditions stated above has been simulated computationally with an algorithm that solves an approximation of the Maxwell-Bloch equations. (These are nonlinear coupled partial differential equations that model the relevant aspects of the physics of semiconductor lasers.) The results of the simulations show that the intensities of the beams oscillate, out of phase relative to each other, at a frequency that is typically of the order of tens of gigahertz. In particular, one simulation was performed for two coupled VCSELs, each having a circular aperture 5.6 microns in diameter, operating at a wavelength of 980 nm. The results of the simulation showed that the intensities of the two beams oscillated in opposite phase (see Figure 1) at a frequency of 42 GHz. In the far field, the beams were found to be separated by an angle of approx.=8deg. The intensities of the beams were also found to oscillate in the near field, 90deg out of phase with the oscillation in the far field. Figure 2 depicts a simple near-field switching device in which an assembly of microlenses would couple the two oscillating beams to two receivers.

Goorjian, Peter M.↗
DOE OSTI2020

Variational nonlinear WKB in the Eulerian frame

Nonlinear WKB is a multiscale technique for studying locally plane-wave solutions of nonlinear partial differential equations (PDEs). Its application comprises two steps: (1) replacement of the original PDE with an extended system separating the large scales from the small and (2) reduction of the extended system to its slow manifold. In the context of variational fluid theories with particle relabeling symmetry, nonlinear WKB in the mean Eulerian frame is known to possess a variational structure. This much has been demonstrated using, for instance, the theoretical apparatus known as the generalized Lagrangian mean. On the other hand, the variational structure of nonlinear WKB in the conventional Eulerian frame remains mysterious. By exhibiting a variational principle for the extended equations from step (1) above, here we demonstrate that nonlinear WKB in the Eulerian frame is in fact variational. Remarkably, the variational principle for the extended system admits loops of relabeling transformations as a symmetry group. Noether’s theorem therefore implies that the extended Eulerian equations possess a family of circulation invariants parameterized by S 1 . As an illustrative example, we use our results to systematically deduce a variational model of high-frequency acoustic waves interacting with a larger-scale compressible isothermal flow.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1988

Applications of an exponential finite difference technique

An exponential finite difference scheme first presented by Bhattacharya for one dimensional unsteady heat conduction problems in Cartesian coordinates was extended. The finite difference algorithm developed was used to solve the unsteady diffusion equation in one dimensional cylindrical coordinates and was applied to two and three dimensional conduction problems in Cartesian coordinates. Heat conduction involving variable thermal conductivity was also investigated. The method was used to solve nonlinear partial differential equations in one and two dimensional Cartesian coordinates. Predicted results are compared to exact solutions where available or to results obtained by other numerical methods.

Handschuh, Robert F.↗
NASA NTRS1978

A solution procedure for behavior of thick plates on a nonlinear foundation and postbuckling behavior of long plates

Approximate solutions for three nonlinear orthotropic plate problems are presented: (1) a thick plate attached to a pad having nonlinear material properties which, in turn, is attached to a substructure which is then deformed; (2) a long plate loaded in inplane longitudinal compression beyond its buckling load; and (3) a long plate loaded in inplane shear beyond its buckling load. For all three problems, the two dimensional plate equations are reduced to one dimensional equations in the y-direction by using a one dimensional trigonometric approximation in the x-direction. Each problem uses different trigonometric terms. Solutions are obtained using an existing algorithm for simultaneous, first order, nonlinear, ordinary differential equations subject to two point boundary conditions. Ordinary differential equations are derived to determine the variable coefficients of the trigonometric terms.

Stein, M.↗
NASA NTRS1993

Control of micromachined deformable mirrors

A micromachined deformable mirror with pixelated electrostatic actuators is proposed. The paper begins with a physical description of the proposed mirror. Then a mathematical model in the form of a nonlinear partial differential equation describing the mirror surface deformations is derived. This model is used to derive the required voltages for the actuators to achieve a specified static deformation of the mirror surface. This is followed by the derivation of a static nonlinear feedback controller for achieving noninteracting actuation. Then the structure for a complete control system for wavefront correction is proposed. The paper concludes with a discussion of the physical implementation of the proposed control system.

Agronin, M. L.↗
NASA NTRS1992

Control of Micromachined Deformable Mirrors

A micromachined deformable mirror with pixelated electrostatic actuators is proposed. The paper begins with a physical description of the proposed mirror. Then, a mathematical model in the form of a nonlinear partial differential equation describing the mirror surface deformations is derived. This model is used to derive the required voltages for the actuators to achieve a specified static deformation of the mirror surface. This is followed by the derivation of a static nonlinear feedback controller for achieving noninteracting actuation. Then, the structure for a complete control system for wavefront correction is proposed. The paper concludes with a discussion of the physical implementation of the proposed control system.

mirrors↗
NASA NTRS1996

Low-Dimensional Dynamical Models of Thermal Convection

A low-dimensional dynamic model for transitional buoyancy-driven flow in a differentially heated tall enclosure is presented. The full governing partial differential equations with the associated boundary conditions are solved by a spectral element method for a cavity of aspect ratio A=20. Proper orthogonal decomposition is applied to the oscillatory solution at Prandtl number Pr=P tau (omega) = 0.71 and Grashof number G tau (omega) = 3.2 x 10 (exp 4) to construct empirical eigenfunctions. Using the four most energetic empirical eigenfunctions for the velocity and temperature as basis functions and applying Galerkin's method, a reduced model consisting of eight nonlinear ordinary differential equations is obtained. Close to the 'design' conditions (P tau(omega) G tau(omega)), the low-order model (LOM) predictions are in excellent agreement with the predictions of the full model. In particular, the critical Grashof number at the onset of the first temporal flow instability (Hopf bifurcation) was well as the frequency and amplitude of oscillations at supercritical conditions are in excellent agreement with the predictions of the full model. Far from the 'design' conditions, the LOM predicts the existence of multiple stable steady solutions at large values of G tau, and a unique stable steady solution at small values of G tau, and exhibits hysteretic behavior that is qualitatively similar to that observed in direct numerical simulations based on the full model.

Liakopoulos, Anthony↗
DOE OSTI2025

Understanding plasma turbulence through exact coherent structures

Plasma turbulence is a key challenge in understanding transport phenomena in magnetically confined plasmas. This work presents a generalized framework to analyze plasma turbulence that utilizes periodic orbit theory. In periodic orbit theory, doubly periodic solutions (coherent structures) of the governing equation(s) serve as building blocks of the considered turbulent dynamics. To illustrate the concept and method, the particularly simple Kuramoto–Sivashinsky (referred to here as LMRT for the original authors: LaQuey, Mahajan, Rutherford, and Tang) trapped-ion mode toy model is used. By applying numerical optimization techniques to the LMRT equation, we extract coherent spacetime patterns that represent the library of allowable fundamental structures of the equation. These structures provide a framework to systematically describe turbulence as a composition of recurrent solutions, revealing an underlying order within chaotic plasma motion. Although illustrated here using the simplified LMRT model for clarity, this framework provides a general strategy that can be extended to more complex and realistic models of plasma turbulence, including gyrokinetic systems. This offers a new method for predicting and potentially controlling transport processes in fusion plasmas by providing a bridge between nonlinear dynamical systems theory and plasma physics in the form of a generalized framework with which to analyze and understand spatially extended nonlinear partial differential equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗