Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonlinear partial 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

A unified mathematical theory of electrophoretic processes

A mathematical theory is presented which shows that each of the four classical electrophoretic modes (zone electrophoresis, moving boundary electrophoresis, isotachophoresis, and isoelectric focusing) is based on the same general principles and can collectively be described in terms of a single set of equations. This model can predict the evolution of the four electrophoretic modes as a function of time. The model system is one-dimensional, neglecting the effects of electroosmosis, temperature gradients, and any bulk flows of liquid. The model is based on equations which express the components' dissociation equilibria, the mass transport due to electromigration and diffusion, electroneutrality, and the conservation of mass and charge. The model consists of a system of coupled partial differential and nonlinear algebraic equations which can be solved numerically by use of a computer. The versatility of this model was verified using an example of a three-component system containing cacodylate, tris hydroxylmethylaminomethane, and histidine. Results show that this model not only correctly predicts the characteristic features of each electrophoretic mode, but also gives details of the concentration, pH, and conductivity profiles not easily amenable to direct experimental measurement.

Bier, M.↗
NASA NTRS1978

Numerical method for solution of systems of non-stationary spatially one-dimensional nonlinear differential equations

A computational scheme and a standard program is proposed for solving systems of nonstationary spatially one-dimensional nonlinear differential equations using Newton's method. The proposed scheme is universal in its applicability and its reduces to a minimum the work of programming. The program is written in the FORTRAN language and can be used without change on electronic computers of type YeS and BESM-6. The standard program described permits the identification of nonstationary (or stationary) solutions to systems of spatially one-dimensional nonlinear (or linear) partial differential equations. The proposed method may be used to solve a series of geophysical problems which take chemical reactions, diffusion, and heat conductivity into account, to evaluate nonstationary thermal fields in two-dimensional structures when in one of the geometrical directions it can take a small number of discrete levels, and to solve problems in nonstationary gas dynamics.

Morozov, S. K.↗
NASA NTRS1981

Galactic civilizations - Population dynamics and interstellar diffusion

A model is developed of the interstellar diffusion of galactic civilizations which takes into account the population dynamics of such civilizations. The problem is formulated in terms of potential theory, with a family of nonlinear partial differential and difference equations specifying population growth and diffusion for an organism with advantageous genes that undergoes random dispersal while increasing in population locally, and a population at zero population growth. In the case of nonlinear diffusion with growth and saturation, it is found that the colonization wavefront from the nearest independently arisen galactic civilization can have reached the earth only if its lifetime exceeds 2.6 million years, or 20 million years if discretization can be neglected. For zero population growth, the corresponding lifetime is 13 billion years. It is concluded that the earth is uncolonized not because interstellar spacefaring civilizations are rare, but because there are too many worlds to be colonized in the plausible colonization lifetime of nearby civilizations, and that there exist no very old galactic civilizations with a consistent policy of the conquest of inhabited worlds.

Newman, W. I.↗
NASA NTRS1995

Numerical solution of differential equations by artificial neural networks

Conventionally programmed digital computers can process numbers with great speed and precision, but do not easily recognize patterns or imprecise or contradictory data. Instead of being programmed in the conventional sense, artificial neural networks (ANN's) are capable of self-learning through exposure to repeated examples. However, the training of an ANN can be a time consuming and unpredictable process. A general method is being developed by the author to mate the adaptability of the ANN with the speed and precision of the digital computer. This method has been successful in building feedforward networks that can approximate functions and their partial derivatives from examples in a single iteration. The general method also allows the formation of feedforward networks that can approximate the solution to nonlinear ordinary and partial differential equations to desired accuracy without the need of examples. It is believed that continued research will produce artificial neural networks that can be used with confidence in practical scientific computing and engineering applications.

Meade, Andrew J., Jr.↗
NASA NTRS1974

Quantum mechanical streamlines. I - Square potential barrier

Exact numerical calculations are made for scattering of quantum mechanical particles hitting a square two-dimensional potential barrier (an exact analog of the Goos-Haenchen optical experiments). Quantum mechanical streamlines are plotted and found to be smooth and continuous, to have continuous first derivatives even through the classical forbidden region, and to form quantized vortices around each of the nodal points. A comparison is made between the present numerical calculations and the stationary wave approximation, and good agreement is found between both the Goos-Haenchen shifts and the reflection coefficients. The time-independent Schroedinger equation for real wavefunctions is reduced to solving a nonlinear first-order partial differential equation, leading to a generalization of the Prager-Hirschfelder perturbation scheme. Implications of the hydrodynamical formulation of quantum mechanics are discussed, and cases are cited where quantum and classical mechanical motions are identical.

Hirschfelder, J. O.↗
NASA NTRS1990

Two-dimensional convection and radiation with scattering from a Poiseuille flow

Two-dimensional combined convection and radiation heat transfer from a gray scattering fluid in a reflecting channel is considered. The model, represented by a set of simultaneous nonlinear integro-partial differential equations, is solved numerically. The effects of aspect ratio, conduction-radiation parameter, scattering albedo, and wall emissivity, are systematically investigated. It is found that these parameters have a significant influence on the temperature field and alter the radiative and convective fluxes at the hot and cold walls. In particular, when radiation effects are considerable, the heat-transfer characteristics of the fluid at the hot and cold walls are very different.

Kassemi, M.↗
NASA NTRS1989

TDFAC - A composite grid method for time dependent problems

This paper presents the Fast Adaptive Composite Grid Method for Time Dependent Problems (TDFAC), a local grid refinement scheme for nonlinear, time dependent partial differential equations. The method is presented in a one dimensional setting to illustrate the basic features and components. Computational results which show the effectiveness of TDFAC as a computational tool are also included.

Heroux, Michael↗
NASA NTRS1993

Development of programmable artificial neural networks

Conventionally programmed digital computers can process numbers with great speed and precision, but do not easily recognize patterns or imprecise or contradictory data. Instead of being programmed in the conventional sense, artificial neural networks are capable of self-learning through exposure to repeated examples. However, the training of an ANN can be a time consuming and unpredictable process. A general method is being developed to mate the adaptability of the ANN with the speed and precision of the digital computer. This method was successful in building feedforward networks that can approximate functions and their partial derivatives from examples in a single iteration. The general method also allows the formation of feedforward networks that can approximate the solution to nonlinear ordinary and partial differential equations to desired accuracy without the need of examples. It is believed that continued research will produce artificial neural networks that can be used with confidence in practical scientific computing and engineering applications.

Meade, Andrew J.↗
NASA NTRS1992

Simulation of particle mixing by turbulent convective flows on the connection machine

Mixing of particles by chaotic flow fields was simulated on the Connection Machine. Each cell was assigned to the processor, and the coordinates of particles residing on the cell were kept in the local memory of the processor. This approach implies the exchange between the local memories, when a particle moves from one cell to another. Approximately 10 exp 5 particles were injected into a time-dependent flow field obtained by solving the nonlinear system of partial differential equations describing turbulent thermal convection. The flow field was calculated on a CRAY, and data were transferred to a CM-200 through a high-speed HIPPI channel.

Malevsky, A. V.↗
NASA NTRS1993

Analysis of Radiation-natural Convection Interactions in 1-g and low-g Environments using the Discrete Exchange Factor Method

A new numerical method is presented for the analysis of combined natural convection and radiation heat transfer with applications in many engineering situations such as materials processing, combustion and fire research. Because of the recent interest in the low gravity environment of space, attention is devoted to both 1-g and low-g applications. The two-dimensional mathematical model is represented by a set of coupled nonlinear integro-partial differential equations. Radiative exchange is formulated using the Discrete Exchange Factor method (DEF). This method considers point to point exchange and provides accurate results over a wide range of radiation parameters. Numerical results show that radiation significantly influences the flow and heat transfer in both low-g and 1-g applications. In the low-g environment, convection is weak, and radiation can easily become the dominant heat transfer mode. It is also shown that volumetric heating by radiation gives rise to an intricate cell pattern in the top heated enclosure.

Kassemi, M.↗
NASA NTRS1972

A computer program for the geometrically nonlinear static and dynamic analysis of arbitrarily loaded shells of revolution, theory and users manual

A digital computer program known as SATANS (static and transient analysis, nonlinear, shells) for the geometrically nonlinear static and dynamic response of arbitrarily loaded shells of revolution is presented. Instructions for the preparation of the input data cards and other information necessary for the operation of the program are described in detail and two sample problems are included. The governing partial differential equations are based upon Sanders' nonlinear thin shell theory for the conditions of small strains and moderately small rotations. The governing equations are reduced to uncoupled sets of four linear, second order, partial differential equations in the meridional and time coordinates by expanding the dependent variables in a Fourier sine or cosine series in the circumferential coordinate and treating the nonlinear modal coupling terms as pseudo loads. The derivatives with respect to the meridional coordinate are approximated by central finite differences, and the displacement accelerations are approximated by the implicit Houbolt backward difference scheme with a constant time interval. The boundaries of the shell may be closed, free, fixed, or elastically restrained. The program is coded in the FORTRAN 4 language and is dimensioned to allow a maximum of 10 arbitrary Fourier harmonics and a maximum product of the total number of meridional stations and the total number of Fourier harmonics of 200. The program requires 155,000 bytes of core storage.

Ball, R. E.↗
NASA NTRS1990

Krylov subspace methods - Theory, algorithms, and applications

Projection methods based on Krylov subspaces for solving various types of scientific problems are reviewed. The main idea of this class of methods when applied to a linear system Ax = b, is to generate in some manner an approximate solution to the original problem from the so-called Krylov subspace span. Thus, the original problem of size N is approximated by one of dimension m, typically much smaller than N. Krylov subspace methods have been very successful in solving linear systems and eigenvalue problems and are now becoming popular for solving nonlinear equations. The main ideas in Krylov subspace methods are shown and their use in solving linear systems, eigenvalue problems, parabolic partial differential equations, Liapunov matrix equations, and nonlinear system of equations are discussed.

Sad, Youcef↗
NASA NTRS1983

Postbuckling of long orthotropic plates in combined shear and compression

The nonlinear large-deflection partial differential equations of von karman for orthotropic plates loaded in combined shear and compression are converted into a set of first-order nonlinear ordinary differential equations by assuming trigonometric functions in one direction. These equations are solved numerically using a two point boundary problem solver which makes use of Newton's method. Results are obtained which determine the postbuckling behavior of rectangular plates with loading up to about three times the buckling load. Both isotropic and orthotropic composite plates are considered. Results show that orthotropic plates may behave quite differently than isotropic plates and that in-plane boundary conditions are important for plates loaded in shear.

Stein, M.↗
NASA NTRS1985

A theory of post-stall transients in axial compression systems. I - Development of equations

An approximate theory is presented for post-stall transients in multistage axial compression systems. The theory leads to a set of three simultaneous nonlinear third-order partial differential equations for pressure rise, and average and disturbed values of flow coefficient, as functions of time and angle around the compressor. By a Galerkin procedure, angular dependence is averaged, and the equations become first order in time. These final equations are capable of describing the growth and possible decay of a rotating-stall cell during a compressor mass-flow transient. It is shown how rotating-stall-like and surgelike motions are coupled through these equations, and also how the instantaneous compressor pumping characteristic changes during the transient stall process.

Moore, F. K.↗
NASA NTRS1977

Effect of flow on quasi-one-dimensional acoustic wave propagation in a variable area duct of finite length

The general equation for the velocity potential of quasi-one-dimensional acoustic wave motion in a variable area, finite duct with one-dimensional flow is derived by using a perturbation technique. The nonlinear second-order partial differential equation is linearized and then solved, by either a power series expansion method or the Runge-Kutta fourth-order method, for harmonic time dependence. The boundary condition taken at the duct mouth is that of matching the impedance of the duct sound field to that of the radiation field at the duct opening. Three axial Mach number variations along the duct axis are considered and the results obtained are compared with those for the case of constant Mach number, to determine the influence of the axial velocity gradient on sound propagation. The effect of flow on the radiation impedance is also considered.

Lumsdaine, E.↗
NASA NTRS1972

The responses of the thermosphere due to a geomagnetic storm: A MHD model

A magnetohydrodynamics theory was used to study the dynamic response of the neutral atmosphere to a geomagnetic storm. A full set of magnetohydrodynamic equations appropriate for the present problem is derived and their various orders of approximation are discussed in some detail. In order to demonstrate the usefulness of this theoretical model, the May 1967 geomagnetic storm data were used in the resulting set of nonlinear, time dependent, partial differential magnetohydrodynamic equations to calculate variations of the thermosphere due to the storm. The numerical results are presented for wind speeds, electric field strength, and amount of joule heating at a constant altitude for the data recorded. Data show that the strongest thermospheric responses are at the polar region becoming weaker in the equatorial region. This may lead to the speculation that a thermospheric wave is generated in the polar region due to the geomagnetic storm which propagates towards the equator.

Wu, S. T.↗
NASA NTRS1987

Selection principles and pattern formation in fluid mechanics and nonlinear shell theory

Research accomplishments are summarized and publications generated under the contract are listed. The general purpose of the research was to investigate various symmetry breaking problems in fluid mechanics by the use of structure parameters and selection principles. Although all of the nonlinear problems studied involved systems of partial differential equations, many of these problems led to the study of a single nonlinear operator equation of the form F(w, lambda, gamma) = 0, (w is an element of H), (lambda is an element of R1), (gamma is an element of R1). Instead of varying only the load parameter lambda, as is often done in the study of such equations, one of the main ideas used was to vary the structure parameter gamma in such a way that stable solutions were obtained. In this way one determines detailed stability results by making use of the structure of the model equations and the known physical parameters of the problem. The approach was carried out successfully for Benard-type convection problems, Taylor-like problems for short cylinders, rotating Couette-Poiseuille channel flows, and plane Couette flows. The main focus of the research was on wave theory of vortex breakdown in a tube. A number of preliminary results for inviscid axisymmetric flows were obtained.

Sather, Duane P.↗