Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonlinear reaction”

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 19 records

EUPDF: Eulerian Monte Carlo Probability Density Function Solver for Applications With Parallel Computing, Unstructured Grids, and Sprays

The success of any solution methodology used in the study of gas-turbine combustor flows depends a great deal on how well it can model the various complex and rate controlling processes associated with the spray's turbulent transport, mixing, chemical kinetics, evaporation, and spreading rates, as well as convective and radiative heat transfer and other phenomena. The phenomena to be modeled, which are controlled by these processes, often strongly interact with each other at different times and locations. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. The influence of turbulence in a diffusion flame manifests itself in several forms, ranging from the so-called wrinkled, or stretched, flamelets regime to the distributed combustion regime, depending upon how turbulence interacts with various flame scales. Conventional turbulence models have difficulty treating highly nonlinear reaction rates. A solution procedure based on the composition joint probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices (such as extinction, blowoff limits, and emissions predictions) because it can account for nonlinear chemical reaction rates without making approximations. In an attempt to advance the state-of-the-art in multidimensional numerical methods, we at the NASA Lewis Research Center extended our previous work on the PDF method to unstructured grids, parallel computing, and sprays. EUPDF, which was developed by M.S. Raju of Nyma, Inc., was designed to be massively parallel and could easily be coupled with any existing gas-phase and/or spray solvers. EUPDF can use an unstructured mesh with mixed triangular, quadrilateral, and/or tetrahedral elements. The application of the PDF method showed favorable results when applied to several supersonic-diffusion flames and spray flames. The EUPDF source code will be available with the National Combustion Code (NCC) as a complete package.

Raju, M. S.↗

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗

Nonlinear structural analysis for fiber-reinforced superalloy turbine blades

A computational capability for predicting the nonlinear thermomechanical structural response of fiber-reinforced superalloy (FRS) turbine blades is described. This capability is embedded in a special purpose computer code (COBSTRAN) developed at the NASA Lewis Research Center. Special features of this computational capability include accounting for: fiber/matrix reaction, nonlinear and anisotropic material behavior, complex stress distribution due to local and global heterogeneity, and residual stresses due to initial fabrication and/or inelastic behavior during subsequent missions. Numerical results are presented from analyses of a hypothetical FRS turbine blade subjected to a fabrication process and subsequent mission cycle. The results demonstrate the capabilities of this computational tool to; predict local stress/strain response and capture trends of local nonlinear and anisotropic material behavior, relate the effects of this local behavior to the global response of a multilayered fiber-composite turbine blade, and trace material history from fabrication through successive missions.

Hopkins, D. A.↗

The role of nonlinear critical layers in boundary layer transition

Asymptotic methods are used to describe the nonlinear self-interaction between pairs of oblique instability modes that eventually develops when initially linear spatially growing instability waves evolve downstream in nominally two-dimensional laminar boundary layers. The first nonlinear reaction takes place locally within a so-called 'critical layer', with the flow outside this layer consisting of a locally parallel mean flow plus a pair of oblique instability waves - which may or may not be accompanied by an associated plane wave. The amplitudes of these waves, which are completely determined by nonlinear effects within the critical layer, satisfy either a single integro-differential equation or a pair of integro-differential equations with quadratic to quartic-type nonlinearities. The physical implications of these equations are discussed.

Goldstein, M.E.↗

Nonlinear a Tollmien-Schlichting/vortex interaction in boundary layers

The nonlinear reaction between two oblique three-dimensional Tollmien-Schlichting (TS) waves and their induced streamwise-vortex flow is considered theoretically for an incompressible boundary layer. The same theory applies to the destabilization of an incident vortex motion by subharmonic TS waves, followed by interaction. The scales and flow structure involved are addressed for high Reynolds numbers. The nonlinear interaction is powerful, starting at quite low amplitudes with a triple-deck structure for the TS waves but a large-scale structure for the induced vortex, after which strong nonlinear amplification occurs. This includes nonparallel-flow effects. The nonlinear interaction is governed by a partial differential system for the vortex flow coupled with an ordinary-differential one for the TS pressure. The solution properties found sometimes produce a breakup within a finite distance and sometimes further downstream, depending on the input amplitudes upstream and on the wave angles, and that then leads to the second stages of interaction associated with higher amplitudes, the main second stages giving either long-scale phenomena significantly affected by nonparallelism or shorter quasi-parallel ones governed by the full nonlinear triple-deck response.

Hall, P.↗

The Effect of Nonlinear Critical Layers on Boundary Layer Transition

Asymptotic methods are used to describe the nonlinear self-interaction between pairs of oblique instability modes that eventually develops when initially linear and spatially growing instability waves evolve downstream in nominally two-dimensional and spanwise periodic laminar boundary layers. The first nonlinear reaction takes place locally within a so-called 'critical layer' with the flow outside this layer consisting of a locally parallel mean flow plus an appropriate superposition of linear instability waves. The amplitudes of these waves are determined by either a single integro-differential equation or by a pair of integro-differential equations with quadratic to quartic-type nonlinearities.

Goldstein, Marvin E.↗

Nonlinear Tollmien-Schlichting/vortex interaction in boundary layers

The nonlinear reaction between two oblique 3-D Tollmein-Schlichting (TS) waves and their induced streamwise-vortex flow is considered theoretically for an imcompressible boundary layer. The same theory applies to the destabilization of an incident vortex motion by subharmonic TS waves, followed by interaction. The scales and flow structure involved are addressed for high Reynolds numbers. The nonlionear interaction is powerful, starting at quite low amplitudes with a triple-deck structure for the TS waves but a large-scale structure for the induced vortex, after which strong nonlinear amplification occurs. This includes nonparallel-flow effects. The nonlinear interaction is governed by a partial differential system for the vortex flow coupled with an ordinary-differential one for the TS pressure. The solution properties found sometimes produce a breakup within a finite distance and sometimes further downstream, depending on the input amplitudes upstream and on the wave angles, and that then leads to the second stages of interaction associated with higher amplitudes, the main second stages giving either long-scale phenomena significantly affected by nonparallelism or shorter quasi-parallel ones governed by the full nonlinear triple-deck response.

Hall, P.↗

Correspondence between discrete and continuous models of excitable media: trigger waves

We present a theoretical framework for relating continuous partial differential equation (PDE) models of excitable media to discrete cellular automata (CA) models on a randomized lattice. These relations establish a quantitative link between the CA model and the specific physical system under study. We derive expressions for the CA model's plane wave speed, critical curvature, and effective diffusion constant in terms of the model's internal parameters (the interaction radius, excitation threshold, and time step). We then equate these expressions to the corresponding quantities obtained from solution of the PDEs (for a fixed excitability). This yields a set of coupled equations with a unique solution for the required CA parameter values. Here we restrict our analysis to "trigger" wave solutions obtained in the limiting case of a two-dimensional excitable medium with no recovery processes. We tested the correspondence between our CA model and two PDE models (the FitzHugh-Nagumo medium and a medium with a "sawtooth" nonlinear reaction source) and found good agreement with the numerical solutions of the PDEs. Our results suggest that the behavior of trigger waves is actually controlled by a small number of parameters.

Non-NASA Center↗

On the numerical treatment of nonlinear source terms in reaction-convection equations

The objectives of this paper are to investigate how various numerical treatments of the nonlinear source term in a model reaction-convection equation can affect the stability of steady-state numerical solutions and to show under what conditions the conventional linearized analysis breaks down. The underlying goal is to provide part of the basic building blocks toward the ultimate goal of constructing suitable numerical schemes for hypersonic reacting flows, combustions and certain turbulence models in compressible Navier-Stokes computations. It can be shown that nonlinear analysis uncovers much of the nonlinear phenomena which linearized analysis is not capable of predicting in a model reaction-convection equation.

Lafon, A.↗

A deterministic particle method for one-dimensional reaction-diffusion equations

We derive a deterministic particle method for the solution of nonlinear reaction-diffusion equations in one spatial dimension. This deterministic method is an analog of a Monte Carlo method for the solution of these problems that has been previously investigated by the author. The deterministic method leads to the consideration of a system of ordinary differential equations for the positions of suitably defined particles. We then consider the time explicit and implicit methods for this system of ordinary differential equations and we study a Picard and Newton iteration for the solution of the implicit system. Next we solve numerically this system and study the discretization error both analytically and numerically. Numerical computation shows that this deterministic method is automatically adaptive to large gradients in the solution.

Mascagni, Michael↗

Model reference adaptive attitude control of spacecraft using reaction wheels

A nonlinear model reference adaptive control law for large angle rotational maneuvers of spacecraft using reaction wheels in the presence of uncertainty is presented. The derivation of control law does not require any information on the values of the system parameters and the disturbance torques acting on the spacecraft. The controller includes a dynamic system in the feedback path. The control law is a nonlinear function of the attitude error, the rate of the attitude error, and the compensator state. Simulation results are prsented to show that large angle rotational maneuvers can be performed in spite of the uncertainty in the system.

Singh, Sahjendra N.↗

The thermal stability of coronal loops by nonlinear diffusion asymptotics

A nonlinear reaction-diffusion equation and some additional constraints are derived which describe the time-dependent behavior of the temperature structure of the plasma in coronal loops. The equation is analyzed using nonlinear diffusion asymptotics, in particular singular perturbation techniques, and the results are interpreted in the context of the physical problem of the thermal stability and temporal behavior of the plasma. The results are consistent with the possibility of cyclic thermal behavior of the plasma, as suggested by Kuin and Martens (1982).

Pakkert, J. W.↗

Meteorological and air pollution modeling for an urban airport

Results are presented of numerical experiments modeling meteorology, multiple pollutant sources, and nonlinear photochemical reactions for the case of an airport in a large urban area with complex terrain. A planetary boundary-layer model which predicts the mixing depth and generates wind, moisture, and temperature fields was used; it utilizes only surface and synoptic boundary conditions as input data. A version of the Hecht-Seinfeld-Dodge chemical kinetics model is integrated with a new, rapid numerical technique; both the San Francisco Bay Area Air Quality Management District source inventory and the San Jose Airport aircraft inventory are utilized. The air quality model results are presented in contour plots; the combined results illustrate that the highly nonlinear interactions which are present require that the chemistry and meteorology be considered simultaneously to make a valid assessment of the effects of individual sources on regional air quality.

Swan, P. R.↗

Post-buckling behavior of a beam-column on a nonlinear elastic foundation with a gap

The structural behavior of an elastic beam-column placed with a gap between two nonlinearity elastic layers each resting on a rigid foundation was examined. The beam-column was laterally supported at both ends and subjected to a uniform transverse load and axial compression. Its slenderness was such that the axial compressive force exceeds the amount that would be necessary to buckle it as a simple supported column. The elastic layers were represented by an elastic foundation with a strongly nonlinear specific reaction taken as a rapidly increasing function of the layer compression. The analytical model developed simulated the entire pattern of the deflection and stress state including layer and end support reactions, under gradually increasing axial force.

Kuznetsov, E. N.↗

Modeling of the aerodynamic response to arbitrary aircraft maneuvres - A numerical validation

The regime of validity of an aerodynamic mathematical model, applicable to describe the nonlinear aerodynamic reactions to a delta wing maneuvering at high angles of attack is investigated. An unsteady vortex-lattice method is used to compute the unsteady flowfields, and thus to evaluate the aerodynamic data required by the model, in terms of specified characteristic motions. Time-histories of the aerodynamic responses to complex motions are generated by means of the model and the evaluated aerodynamic data, and are compared with baseline aerodynamic responses obtained from direct vortex-lattice computations. The validity of the mathematical modeling approach for the maneuvering delta wing is demonstrated by the close agreement of the force and moment responses obtained from the two approaches.

Katz, J.↗