Engineering Papers⌕ Search

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.

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

At least 847 records · Page 47

A model for the three-dimensional spacecraft control laboratory experiment

A model for the three-dimensional Spacecraft Control Laboratory Experiment (SCOLE) is developed. The objective behind this method of modelling is to utilize the basic partial differential equations of motion for this distributed parameter system and not to use the modal expansion in developing the model. The final model obtained is in terms of a transfer function matrix which relates the flexible mast parameters like displacement, slope, shear stress, etc. to external forces and moments.

Kakad, Yogendra P.↗

Comparison of SMAC, PISO, and iterative time-advancing schemes for unsteady flows

Calculations of unsteady flows using a simplified marker and cell (SMAC), a pressure implicit splitting of operators (PSIO), and an iterative time advancing scheme (ITA) are presented. A partial differential equation for incremental pressure is used in each time advancing scheme. Example flows considered are a polar cavity flow starting from rest and self-sustained oscillating flows over a circular and a square cylinder. For a large time step size, the SMAC and ITA are more strongly convergent and yield more accurate results than PSIO. The SMAC is the most efficient computationally. For a small time step size, the three time advancing schemes yield equally accurate Strouhal numbers. The capability of each time advancing scheme to accurately resolve unsteady flows is attributed to the use of new pressure correction algorithm that can strongly enforce the conservation of mass. The numerical results show that the low frequency of the vortex shedding is caused by the growth time of each vortex shed into the wake region.

Kim, Sang-Wook↗

Sinc-Galerkin estimation of diffusivity in parabolic problems

A fully Sinc-Galerkin method for the numerical recovery of spatially varying diffusion coefficients in linear partial differential equations is presented. Because the parameter recovery problems are inherently ill-posed, an output error criterion in conjunction with Tikhonov regularization is used to formulate them as infinite-dimensional minimization problems. The forward problems are discretized with a sinc basis in both the spatial and temporal domains thus yielding an approximate solution which displays an exponential convergence rate and is valid on the infinite time interval. The minimization problems are then solved via a quasi-Newton/trust region algorithm. The L-curve technique for determining an approximate value of the regularization parameter is briefly discussed, and numerical examples are given which show the applicability of the method both for problems with noise-free data as well as for those whose data contains white noise.

Smith, Ralph C.↗

A manual for PARTI runtime primitives, revision 1

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.

Das, Raja↗

Mapping implicit spectral methods to distributed memory architectures

Spectral methods were proven invaluable in numerical simulation of PDEs (Partial Differential Equations), but the frequent global communication required raises a fundamental barrier to their use on highly parallel architectures. To explore this issue, a 3-D implicit spectral method was implemented on an Intel hypercube. Utilization of about 50 percent was achieved on a 32 node iPSC/860 hypercube, for a 64 x 64 x 64 Fourier-spectral grid; finer grids yield higher utilizations. Chebyshev-spectral grids are more problematic, since plane-relaxation based multigrid is required. However, by using a semicoarsening multigrid algorithm, and by relaxing all multigrid levels concurrently, relatively high utilizations were also achieved in this harder case.

Overman, Andrea L.↗

Mapping unstructured grid problems to the connection machine

We present a highly parallel graph mapping technique that enables one to solve unstructured grid problems on massively parallel computers. Many implicit and explicit methods for solving discretizated partial differential equations require each point in the discretization to exchange data with its neighboring points every time step or iteration. The time spent communicating can limit the high performance promised by massively parallel computing. To eliminate this bottleneck, we map the graph of the irregular problem to the graph representing the interconnection topology of the computer such that the sum of the distances that the messages travel is minimized. We show that, in comparison to a naive assignment of processors, our heuristic mapping algorithm significantly reduces the communication time on the Connection Machine, CM-2.

Hammond, Steven W.↗

Long-time behaviour of absorbing boundary conditions

A new class of computational far-field boundary conditions for hyperbolic partial differential equations was recently introduced by the authors. These boundary conditions combine properties of absorbing conditions for transient solutions and properties of far-field conditions for steady states. This paper analyses the properties of the wave equation coupled with these new boundary conditions: well-posedness, dissipativity and convergence in time.

Engquist, B.↗

Simulation of the collapse and dissipation of Langmuir wave packets

Particle-in-cell (PIC) simulations and Zakharov's partial differential equations (PDEs) are used to investigate the collapse of isolated Langmuir wave packets in two dimensions. Collapse thresholds are determined numerically, and the roles of enhanced Langmuir wave damping and nonlinearities not included in the standard Zakharov equations are discussed. The Langmuir wave and ion dynamics in PIC simulations are compared with the predictions of PDE simulations incorporating enhanced Langmuir damping. Electron heating and coherent acceleration in the PIC simulations are discussed and compared with predictions of the transit-time theory.

Newman, D. L.↗

Statistical mechanics models for motion and force planning

The models of statistical mechanics provide an alternative to the methods of classical mechanics more traditionally used in robotics. They have a potential to: improve analysis of object collisions; handle kinematic and dynamic contact interactions within the same frmework; and reduce the need for perfect deterministic world model information. The statistical mechanics models characterize the state of the system as a probability density function (p.d.f.) whose time evolution is governed by a partial differential equation subject to boundary and initial conditions. The boundary conditions when rigid objects collide reflect the conservation of momentum. The models are being developed to embedd in remote semi-autonomous systems with a need to reason and interact with a multiobject environment.

Rodriguez, G.↗

A three-dimensional approach for analysis of sidewall injector mixing and combustion

A three-dimensional (3-D) analytical approach to calculate the mixing and combustion efficiency along engine combustors that use sidewall fuel injectors is described. The analysis consists of three parts: (1) an empirical correlation for fuel penetration, (2) application of integral conservation laws following jet turning and pressure equilibration, and (3) an analytical solution for the simplified three-dimensional partial differential equations to describe the downstream mixing and combustion process. Numerical results are compared to some empirical and experimental work, and favorable agreements are shown. The importance of including anisotropic turbulence and three-dimensional effects are illustrated by comparison with data and other less comprehensive analyses. The analytical solution is suitable for rapid estimates of engine performance and parametric studies of combustor designs containing transverse or angled injector ports arranged in rows and/or columns.

Quan, Victor↗

Statistical-mechanical selection of the shapes of disk galaxies

A new method is proposed for selecting steady state shapes of disk galaxies as 'most probable states' of a large number of stars, given only the total energy and total angular momentum. A partial differential equation is derived for the mean gravitational potential; it is closely related to the 'sinh-Poisson' equation for the mean-field description of a line vortex system or electrostatic guiding-center plasma. A 'water bag' approximation to the distribution function for bound stars renders the equation analytically tractable, but accurate solution of it may require numerical integration, in view of its general nonlinearity.

Montgomery, David↗

Time domain convergence properties of Lyapunov stable penalty methods

Linear hyperbolic partial differential equations are analyzed using standard techniques to show that a sequence of solutions generated by the Liapunov stable penalty equations approaches the solution of the differential-algebraic equations governing the dynamics of multibody problems arising in linear vibrations. The analysis does not require that the system be conservative and does not impose any specific integration scheme. Variational statements are derived which bound the error in approximation by the norm of the constraint violation obtained in the approximate solutions.

Kurdila, A. J.↗

Modeling flow at the nozzle of a solid rocket motor

The mechanical behavior of a rocket motor internal flow field results in a system of nonlinear partial differential equations which can be solved numerically. The accuracy and the convergence of the solution of the system of equations depends largely on how precisely the sharp gradients can be resolved. An adaptive grid generation scheme is incorporated into the computer algorithm to enhance the capability of numerical modeling. With this scheme, the grid is refined as the solution evolves. This scheme significantly improves the methodology of solving flow problems in rocket nozzle by putting the refinement part of grid generation into the computer algorithm.

Chow, Alan S.↗

A theoretical model for intense microwave pulse propagation in an air breakdown environment

A theoretical model based on two coupled partial differential equations is established to describe the propagation of an intense microwave pulse in air breakdown environment. One is derived from the Poynting theorem, and the other one is the rate equation of electron density. A semiempirical formula of the ionization frequency is adopted for this model. A transformation of these two equations to local time frame of reference is introduced so that they can be solved numerically with considerably reduced computation time.

Kuo, S. P.↗

On damping mechanisms in beams

A partial differential equation model of a cantilevered beam with a tip mass at its free end is used to study damping in a composite. Four separate damping mechanisms consisting of air damping, strain rate damping, spatial hysteresis and time hysteresis are considered experimentally. Dynamic tests were performed to produce time histories. The time history data is then used along with an approximate model to form a sequence of least squares problems. The solution of the least squares problem yields the estimated damping coefficients. The resulting experimentally determined analytical model is compared with the time histories via numerical simulation of the dynamic response. The procedure suggested here is compared with a standard modal damping ratio model commonly used in experimental modal analysis.

Banks, H. T.↗

Scheme For Finite-Difference Computations Of Waves

Compact algorithms generating and solving finite-difference approximations of partial differential equations for propagation of waves obtained by new method. Based on concept of discrete dispersion relation. Used in wave propagation to relate frequency to wavelength and is key measure of wave fidelity.

Davis, Sanford↗

Stability of the flow around a cylinder: The spin-up problem

A concern is the flow around an infinite cylinder, which at a certain instant impulsively starts to spin. The growth of vortices in the resulting boundary layer occurring outside the cylinder is investigated. This layer is essentially a Rayleigh layer which grows with time, so the mechanism involved is similar to that studied in Hall (1983). Vortices with wavenumber comparable to the layer thickness are shown to be described by partial differential equations that govern the system numerically. It is assumed that the Rayleigh layer is thin, so particles are confined to move in a path with radius of curvature the same as the cylinder. The Goertler number is a function of time, so the time scale which produces an order, is considered one Goertler number. The right hand branch calculation is considered by letting the time tend to infinity, also inviscid Goertler modes are considered.

Otto, S. R.↗

Modelling crystal growth: Convection in an asymmetrically heated ampoule

The objective was to develop and implement a numerical method capable of solving the nonlinear partial differential equations governing heat, mass, and momentum transfer in a 3-D cylindrical geometry in order to examine the character of convection in an asymmetrically heated cylindrical ampoule. The details of the numerical method, including verification tests involving comparison with results obtained from other methods, are presented. The results of the study of 3-D convection in an asymmetrically heated cylinder are described.

Alexander, J. Iwan D.↗