SEARCH · Engineering Papers

Results for “partial differential equation”

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 829 records · Page 46

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Three-dimensional adaptive grid generation for body-fitted coordinate system

This report describes a numerical method for generating 3-D grids for general configurations. The basic method involves the solution of a set of quasi-linear elliptic partial differential equations via pointwise relaxation with a local relaxation factor. It allows specification of the grid spacing off the boundary surfaces and the grid orthogonality at the boundary surfaces. It includes adaptive mechanisms to improve smoothness, orthogonality, and flow resolution in the grid interior.

Chen, S. C.↗

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 NTRS1990

Run-time scheduling and execution of loops on message passing machines

Sparse system solvers and general purpose codes for solving partial differential equations are examples of the many types of problems whose irregularity can result in poor performance on distributed memory machines. Often, the data structures used in these problems are very flexible. Crucial details concerning loop dependences are encoded in these structures rather than being explicitly represented in the program. Good methods for parallelizing and partitioning these types of problems require assignment of computations in rather arbitrary ways. Naive implementations of programs on distributed memory machines requiring general loop partitions can be extremely inefficient. Instead, the scheduling mechanism needs to capture the data reference patterns of the loops in order to partition the problem. First, the indices assigned to each processor must be locally numbered. Next, it is necessary to precompute what information is needed by each processor at various points in the computation. The precomputed information is then used to generate an execution template designed to carry out the computation, communication, and partitioning of data, in an optimized manner. The design is presented for a general preprocessor and schedule executer, the structures of which do not vary, even though the details of the computation and of the type of information are problem dependent.

Saltz, Joel↗

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 NTRS1990

An abbreviated Reynolds stress turbulence model for airfoil flows

An abbreviated Reynolds stress turbulence model is presented for solving turbulent flow over airfoils. The model consists of two partial differential equations, one for the Reynolds shear stress and the other for the turbulent kinetic energy. The normal stresses and the dissipation rate of turbulent kinetic energy are computed from algebraic relationships having the correct asymptotic near wall behavior. This allows the model to be integrated all the way to the wall without the use of wall functions. Results for a flat plate at zero angle of attack, a NACA 0012 airfoil and a RAE 2822 airfoil are presented.

Gaffney, R. L., Jr.↗

NASA NTRS1990

Partitioning of regular computation on multiprocessor systems

Problem partitioning of regular computation over two dimensional meshes on multiprocessor systems is examined. The regular computation model considered involves repetitive evaluation of values at each mesh point with local communication. The computational workload and the communication pattern are the same at each mesh point. The regular computation model arises in numerical solutions of partial differential equations and simulations of cellular automata. Given a communication pattern, a systematic way to generate a family of partitions is presented. The influence of various partitioning schemes on performance is compared on the basis of computation to communication ratio.

Lee, Fung F.↗

NASA NTRS1990

Distributed parameter estimation for large beam-like structures for general inputs

In this paper, maximum likelihood estimation for distributed parameter models of large flexible structures has been formulated. Distributed parameter models involve far fewer unknown parameters than independent modal characteristics of finite element models. The closed-form solutions for the partial differential equations with general forced inputs have been derived. The closed-form expressions of the sensitivity functions lead to highly efficient algorithms for analyzing ground or on-orbit test results. Numerical simulations with impulse and step inputs to the NASA Mini-MAST trust have been demonstrated. The estimations of modal properties involve its longitudinal elongation modes, lateral bending modes and torsional modes. The results show that distributed parameter models are promising in the parameter estimation of large flexible structures.

Shen, Ji-Yao↗

NASA NTRS1990

A manual for PARTI runtime primitives

Primitives are presented that are designed to help users efficiently program irregular problems (e.g., unstructured mesh sweeps, sparse matrix codes, adaptive mesh partial differential equations solvers) on distributed memory machines. These primitives are also designed for use in compilers for distributed memory multiprocessors. Communications patterns are captured at runtime, and the appropriate send and receive messages are automatically generated.

Berryman, Harry↗

NASA NTRS1990

Spurious frequencies as a result of numerical boundary treatments

The stability theory for finite difference Initial Boundary-Value approximations to systems of hyperbolic partial differential equations states that the exclusion of eigenvalues and generalized eigenvalues is a sufficient condition for stability. The theory, however, does not discuss the nature of numerical approximations in the presence of such eigenvalues. In fact, as was shown previously, for the problem of vortex shedding by a 2-D cylinder in subsonic flow, stating boundary conditions in terms of the primitive (non-characteristic) variables may lead to such eigenvalues, causing perturbations that decay slowly in space and remain periodic time. Characteristic formulation of the boundary conditions avoided this problem. A more systematic study of the behavior of the (linearized) one-dimensional gas dynamic equations under various sets of oscillation-inducing legal boundary conditions is reported.

Abarbanel, Saul↗

NASA NTRS1990

Reduced-order filtering for flexible space structures

There is a need for feedback control of the large flexible space structures which are going to be increasingly important in the future of the space program. These structures are very lightly damped, and vibrations may persist for a long time when the system is disturbed unless an active feedback control strategy is used to damp out the vibrations. The system is best described by a partial differential equation description, but the more common approach is to use a large set of second order differential equations, where a large number of modes must be retained if the mathematics is to provide an adequate description of the dynamical process. Sensors, such as accelerometers and rate gyros, may provide data to the feedback controller so that it may respond appropriately to control the system. The data from the sensors is not perfect, but is subject to noise, called measurement noise, and the dynamical process itself is subject to disturbances referred to as process noise. Filtering the sensor signals to remove the measurement noise, and using the resulting state estimates to control the system are investigated.

Sims, Craig S.↗

NASA NTRS1989

Stochastic histories of dust grains in the interstellar medium

The purpose is to study an evolving system of refractory dust grains within the Interstellar Medium (ISM). This is done via a combination of Monte Carlo processes and a system of partial differential equations, where refractory dust grains formed within supernova remnants and ejecta from high mass loss stars are subjected to the processes of sputtering and collisional fragmentation in the diffuse media and accretion within the cold molecular clouds. In order to record chemical detail, the authors take each new particle to consist of a superrefractory core plus a more massive refractory mantle. The particles are allowed to transfer to and fro between the different phases of the interstellar medium (ISM) - on a time scale of 10(exp 8) years - until either the particles are destroyed or the program finishes at a Galaxy time of 6x10(exp 9) years. The resulting chemical and size spectrum(s) are then applied to various astrophysical problems with the following results. For an ISM which has no collisional fragmentation of the dust grains, roughly 10 percent by mass of the most refractory material survives the rigors of the ISM intact, which leaves open the possibility that fossilized isotopically anomalous material may have been present within the primordial solar nebula. Stuctured or layered refractory dust grains within the model cannot explain the observed interstellar depletions of refractory material. Fragmentation due to grain-grain collisions in the diffuse phase plus the accretion of material in the molecular cloud phase can under certain circumstances cause a bimodal distribution in grain size.

Liffman, Kurt↗

NASA NTRS1990

On the Goertler instability in hypersonic flows: Sutherland law fluids and real gas effects

The Goertler vortex instability mechanism in a hypersonic boundary layer on a curved wall is investigated. The precise roles of the effects of boundary layer growth, wall cooling, and gas dissociation is clarified in the determination of stability properties. It is first assumed that the fluid is an ideal gas with viscosity given by Sutherland's law. It is shown that when the free stream Mach number M is large, the boundary layer divides into two sublayers: a wall layer of O(M sup 3/2) thickness over which the basic state temperature is O(M squared) and a temperature adjustment layer of O(1) thickness over which the basic state temperature decreases monotonically to its free stream value. Goertler vortices which have wavelengths comparable with the boundary layer thickness are referred to as wall modes. It is shown that their downstream evolution is governed by a set of parabolic partial differential equations and that they have the usual features of Goertler vortices in incompressible boundary layers. As the local wavenumber increases, the neutral Goertler number decreases and the center of vortex activity moves towards the temperature adjustment layer. Goertler vortices with wavenumbers of order one or larger must necessarily be trapped in the temperature adjustment layer and it is this mode which is most dangerous. For this mode, it was found that the leading order term in the Goertler number expansion is independent of the wavenumber and is due to the curvature of the basic state. This term is also the asymptotic limit of the neutral Goertler numbers of the wall mode. To determine the higher order corrections terms in the Goertler number expansion, two wall curvature cases are distinguished. Real gas effects were investigated by assuming that the fluid is an ideal dissociating gas. It was found that both gas dissociation and wall cooling are destabilizing for the mode trapped in the temperature adjustment layer, but for the wall mode trapped near the wall the effect of gas dissociation can be either destabilizing or stabilizing.

Fu, Yibin B.↗

NASA NTRS1990

On the receptivity and non-parallel stability of travelling disturbances in rotating disk flow

The generation and evolution of small amplitude wavelength traveling disturbances in rotating disk flow is discussed. The steady rotational speed of the disk is perturbed so as to introduce high frequency oscillations in the flow field. Secondly, surface imperfections are introduced on the disk such as roughness elements. The interaction of these two disturbances will generate the instability waves whose evolution is governed by parabolic partial differential equations that are solved numerically. For the class of disturbances considered (wavelength on the order of Reynolds number), it is found that eigensolutions exist which decay or grow algebraically in the radial direction. However, these solutions grow only for frequencies larger than 4.58 times the steady rotational speed of the disk. The computed receptivity coefficient shows that there is an optimum size of roughness for which these modes are excited the most. The width of these roughness elements in the radial direction is about .1 r(sub 0) where r(sub 0) is the radial location of the roughness. It is also found that the receptivity coefficient is larger for a negative spanwise wavenumber than for a positive one. Typical wave angles found for these disturbances are about -26 degrees.

Balakumar, P.↗

NASA NTRS1991

One shot methods for optimal control of distributed parameter systems 1: Finite dimensional control

The efficient numerical treatment of optimal control problems governed by elliptic partial differential equations (PDEs) and systems of elliptic PDEs, where the control is finite dimensional is discussed. Distributed control as well as boundary control cases are discussed. The main characteristic of the new methods is that they are designed to solve the full optimization problem directly, rather than accelerating a descent method by an efficient multigrid solver for the equations involved. The methods use the adjoint state in order to achieve efficient smoother and a robust coarsening strategy. The main idea is the treatment of the control variables on appropriate scales, i.e., control variables that correspond to smooth functions are solved for on coarse grids depending on the smoothness of these functions. Solution of the control problems is achieved with the cost of solving the constraint equations about two to three times (by a multigrid solver). Numerical examples demonstrate the effectiveness of the method proposed in distributed control case, pointwise control and boundary control problems.

Taasan, Shlomo↗

NASA NTRS1990

The method of lines in analyzing solids containing cracks

A semi-numerical method is reviewed for solving a set of coupled partial differential equations subject to mixed and possibly coupled boundary conditions. The line method of analysis is applied to the Navier-Cauchy equations of elastic and elastoplastic equilibrium to calculate the displacement distributions in various, simple geometry bodies containing cracks. The application of this method to the appropriate field equations leads to coupled sets of simultaneous ordinary differential equations whose solutions are obtained along sets of lines in a discretized region. When decoupling of the equations and their boundary conditions is not possible, the use of a successive approximation procedure permits the analytical solution of the resulting ordinary differential equations. The use of this method is illustrated by reviewing and presenting selected solutions of mixed boundary value problems in three dimensional fracture mechanics. These solutions are of great importance in fracture toughness testing, where accurate stress and displacement distributions are required for the calculation of certain fracture parameters. Computations obtained for typical flawed specimens include that for elastic as well as elastoplastic response. Problems in both Cartesian and cylindrical coordinate systems are included. Results are summarized for a finite geometry rectangular bar with a central through-the-thickness or rectangular surface crack under remote uniaxial tension. In addition, stress and displacement distributions are reviewed for finite circular bars with embedded penny-shaped cracks, and rods with external annular or ring cracks under opening mode tension. The results obtained show that the method of lines presents a systematic approach to the solution of some three-dimensional mechanics problems with arbitrary boundary conditions. The advantage of this method over other numerical solutions is that good results are obtained even from the use of a relatively coarse grid.

Gyekenyesi, John P.↗

NASA NTRS1991

On the anomaly of velocity-pressure decoupling in collocated mesh solutions

The use of various pressure correction algorithms originally developed for fully staggered meshes can yield a velocity-pressure decoupled solution for collocated meshes. The mechanism that causes velocity-pressure decoupling is identified. It is shown that the use of a partial differential equation for the incremental pressure eliminates such a mechanism and yields a velocity-pressure coupled solution. Example flows considered are a three dimensional lid-driven cavity flow and a laminar flow through a 90 deg bend square duct. Numerical results obtained using the collocated mesh are in good agreement with the measured data and other numerical results.

Kim, Sang-Wook↗

NASA NTRS1991

Spillover, nonlinearity, and flexible structures

Many systems whose evolution in time is governed by Partial Differential Equations (PDEs) are linearized around a known equilibrium before Computer Aided Control Engineering (CACE) is considered. In this case, there are infinitely many independent vibrational modes, and it is intuitively evident on physical grounds that infinitely many actuators would be needed in order to control all modes. A more precise, general formulation of this grave difficulty (spillover problem) is due to A.V. Balakrishnan. A possible route to circumvention of this difficulty lies in leaving the PDE in its original nonlinear form, and adding the essentially finite dimensional control action prior to linearization. One possibly applicable technique is the Liapunov Schmidt rigorous reduction of singular infinite dimensional implicit function problems to finite dimensional implicit function problems. Omitting details of Banach space rigor, the formalities of this approach are given.

Bass, Robert W.↗

NASA NTRS1991

Likelihood estimation for distributed parameter models for NASA Mini-MAST truss

A maximum likelihood estimation for distributed parameter models of large flexible structures was formulated. Distributed parameter models involve far fewer unknown parameters than independent modal characteristics or finite element models. The closed form solutions for the partial differential equations with corresponding boundary conditions were derived. The closed-form expressions of sensitivity functions led to highly efficient algorithms for analyzing ground or on-orbit test results. For an illustration of this approach, experimental data of the NASA Mini-MAST truss was used. The estimations of modal properties involve lateral bending modes and torsional modes. The results show that distributed parameter models are promising in the parameter estimation of large flexible structures.

Shen, Ji-Yao↗