Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

Application of the method of integral relations to boundary layer flows over blunt bodies

Calculations of boundary layer flows past blunt bodies at angles of incidence are presented. Using the method of integral relations together with the method of lines, the full three-dimensional boundary layer equations are reduced to a system of first-order ordinary differential equations. The streamwise shear stress function and the cross-flow velocity component are represented as suitable functions of the streamwise velocity component. The role of the zone of dependence is automatically satisfied by the choice of differencing in the method of lines. Solutions correct to the second order are obtained in the positive shear region for flow over an ellipsoid at 30-deg incidence. The results are compared with corresponding finite difference solutions.

Modarress, D.↗
NASA NTRS1980

The primer vector in linear, relative-motion equations

Primer vector theory is used in analyzing a set of linear, relative-motion equations - the Clohessy-Wiltshire equations - to determine the criteria and necessary conditions for an optimal, N-impulse trajectory. Since the state vector for these equations is defined in terms of a linear system of ordinary differential equations, all fundamental relations defining the solution of the state and costate equations, and the necessary conditions for optimality, can be expressed in terms of elementary functions. The analysis develops the analytical criteria for improving a solution by (1) moving any dependent or independent variable in the initial and/or final orbit, and (2) adding intermediate impulses. If these criteria are violated, the theory establishes a sufficient number of analytical equations. The subsequent satisfaction of these equations will result in the optimal position vectors and times of an N-impulse trajectory. The solution is examined for the specific boundary conditions of (1) fixed-end conditions, two-impulse, and time-open transfer; (2) an orbit-to-orbit transfer; and (3) a generalized rendezvous problem. A sequence of rendezvous problems is solved to illustrate the analysis and the computational procedure.

Source record↗
NASA NTRS1982

A numerical method for integrating the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes

An extension of the method of moments is developed for the numerical integration of the kinetic equations of droplet spectra evolution by condensation/evaporation and by coalescence/breakup processes. The number density function n sub k (x,t) in each separate droplet packet between droplet mass grid points (x sub k, x sub k+1) is represented by an expansion in orthogonal polynomials with a given weighting function. In this way droplet number concentrations, liquid water contents and other moments in each droplet packet are conserved and the problem of solving the kinetic equations is replaced by one of solving a set of coupled differential equations for the number density function moments. The method is tested against analytic solutions of the corresponding kinetic equations. Numerical results are obtained for different coalescence/breakup and condensation/evaporation kernels and for different initial droplet spectra. Also droplet mass grid intervals, weighting functions, and time steps are varied.

Emukashvily, I. M.↗
NASA NTRS1981

Real time estimation of the heaving and pitching motions of a ship, using a Kalman filter

In the present study the estimation of the heave and pitch motion of a ship is considered, using Kalman filtering techniques. A significant part of the study is devoted to constructing appropriate models for the sea and the ship. The governing equations are obtained from hydrodynamic considerations in the form of linear differential equations with frequency dependent coefficients. In addition, nonminimum phase characteristics are obtained due to the spatial integration of the water wave forces. The resulting transfer matrix function is irrational and nonminimum phase. The conditions for a finite-dimensional approximation are considered and the impact of the various parameters is assessed. A numerical application is considered for a DD-963 destroyer.

Triantafyllou, M.↗
NASA NTRS1983

Real time estimation of ship motions using Kalman filtering techniques

The estimation of the heave, pitch, roll, sway, and yaw motions of a DD-963 destroyer is studied, using Kalman filtering techniques, for application in VTOL aircraft landing. The governing equations are obtained from hydrodynamic considerations in the form of linear differential equations with frequency dependent coefficients. In addition, nonminimum phase characteristics are obtained due to the spatial integration of the water wave forces. The resulting transfer matrix function is irrational and nonminimum phase. The conditions for a finite-dimensional approximation are considered and the impact of the various parameters is assessed. A detailed numerical application for a DD-963 destroyer is presented and simulations of the estimations obtained from Kalman filters are discussed.

Triantafyllou, M. S.↗
NASA NTRS1985

The Formation of Preplanetary Disks from the Collapse of Rotating Molecular Cloud Cores

Solutions that describe the collapse of a molecular cloud core that is initially in unstable equilibrium, embedded within an envelope of uniform density, and rotating at the same rate as the envelope are given. Hydrodynamic equations, including self gravity, are deduced to a set of ordinary differential equations, which are solved by the method of matched asymptotic expansions. Results of these calculations are: (1) the range of stellar masses derived seems to correspond to realistic ranges of observed stellar masses and interstellar cloud parameters, (2) the proper measure of dissipation rate is the ratio of accretion time to viscous diffusion time, and (3) the pressure distribution on the surface of an accreting protostar is nonuniform in a way that favors the channeling of a stellar wind into a bipolar flow directed along the rotation axis.

Cassen, P.↗
NASA NTRS1985

Viscous flows in rotating pipes with internal cavities and related vortex breakdowns

The problem of predicting vortex breakdown, which is characterized by a drastic change in core size and trajectory, is still basically unsolved. In the present paper, viscous flow in a long rotating pipe is analyzed, as a model for the type of breakdown observed in wake vortices. The full Navier-Stokes equations are shown to simplify to a set of unlinked partial differential equations. These equations allow for two separate solutions, for both full pipe flows (the classic rotating Hagen-Poiseuille flow) and cavity flows, in which a central area, with no net mass flux or viscous dissipation, is produced. Conditions permitting switching from one type of flow to the other, with no change in the mass flux or velocities at the pipe surface, are found. It is then argued that these conditions indicate when vortex breakdown will occur.

Weihs, D.↗
NASA NTRS1987

Optimal cure cycle design of a resin-fiber composite laminate

A unified computed aided design method was studied for the cure cycle design that incorporates an optimal design technique with the analytical model of a composite cure process. The preliminary results of using this proposed method for optimal cure cycle design are reported and discussed. The cure process of interest is the compression molding of a polyester which is described by a diffusion reaction system. The finite element method is employed to convert the initial boundary value problem into a set of first order differential equations which are solved simultaneously by the DE program. The equations for thermal design sensitivities are derived by using the direct differentiation method and are solved by the DE program. A recursive quadratic programming algorithm with an active set strategy called a linearization method is used to optimally design the cure cycle, subjected to the given design performance requirements. The difficulty of casting the cure cycle design process into a proper mathematical form is recognized. Various optimal design problems are formulated to address theses aspects. The optimal solutions of these formulations are compared and discussed.

Hou, Jean W.↗
NASA NTRS1990

SIAM Conference on Parallel Processing for Scientific Computing, 4th, Chicago, IL, Dec. 11-13, 1989, Proceedings

Attention is given to such topics as an evaluation of block algorithm variants in LAPACK and presents a large-grain parallel sparse system solver, a multiprocessor method for the solution of the generalized Eigenvalue problem on an interval, and a parallel QR algorithm for iterative subspace methods on the CM2. A discussion of numerical methods includes the topics of asynchronous numerical solutions of PDEs on parallel computers, parallel homotopy curve tracking on a hypercube, and solving Navier-Stokes equations on the Cedar Multi-Cluster system. A section on differential equations includes a discussion of a six-color procedure for the parallel solution of elliptic systems using the finite quadtree structure, data parallel algorithms for the finite element method, and domain decomposition methods in aerodynamics. Topics dealing with massively parallel computing include hypercube vs. 2-dimensional meshes and massively parallel computation of conservation laws. Performance and tools are also discussed.

Dongarra, Jack↗
NASA NTRS1951

Analysis of plane-plastic stress problems with axial symmetry in strain-hardening range

A simple method is developed for solving plane-plastic-stress problems with axial symmetry in the strain-hardening range which is based on the deformation theory of plasticity employing the finite-strain concept. The equations defining the problems are first reduced to two simultaneous nonlinear differential equations involving two dependent variables: (a) the octahedral shear strain, and (b) a parameter indicating the ratio of principal stresses. By multiplying the load and dividing the radius by an arbitrary constant, it is possible to solve these problems without iteration for any value of the modified load. The constant is determined by the boundary condition. This method is applied to a circular membrane under pressure, a rotating disk without and with a central hole, and an infinite plate with a circular hole. Two materials, inconel x and 16-25-6, the octahedral shear stress-strain relations of which do not follow the power law, are used. Distributions of octahedral shear strain, as well as of principal stresses and strains, are obtained. These results are compared with the results of the same problems in the elastic range.

Wu, M H Lee↗
NASA NTRS1996

Finite Difference Time Marching in the Frequency Domain: A Parabolic Formulation for Aircraft Acoustic Nacelle Design

An explicit finite difference iteration scheme is developed to study harmonic sound propagation in aircraft engine nacelles. To reduce storage requirements for large 3D problems, the time dependent potential form of the acoustic wave equation is used. To insure that the finite difference scheme is both explicit and stable, time is introduced into the Fourier transformed (steady-state) acoustic potential field as a parameter. Under a suitable transformation, the time dependent governing equation in frequency space is simplified to yield a parabolic partial differential equation, which is then marched through time to attain the steady-state solution. The input to the system is the amplitude of an incident harmonic sound source entering a quiescent duct at the input boundary, with standard impedance boundary conditions on the duct walls and duct exit. The introduction of the time parameter eliminates the large matrix storage requirements normally associated with frequency domain solutions, and time marching attains the steady-state quickly enough to make the method favorable when compared to frequency domain methods. For validation, this transient-frequency domain method is applied to sound propagation in a 2D hard wall duct with plug flow.

Baumeister, Kenneth J.↗
NASA NTRS1996

Finite Difference Time Marching in the Frequency Domain: A Parabolic Formulation for the Convective Wave Equation

An explicit finite difference iteration scheme is developed to study harmonic sound propagation in ducts. To reduce storage requirements for large 3D problems, the time dependent potential form of the acoustic wave equation is used. To insure that the finite difference scheme is both explicit and stable, time is introduced into the Fourier transformed (steady-state) acoustic potential field as a parameter. Under a suitable transformation, the time dependent governing equation in frequency space is simplified to yield a parabolic partial differential equation, which is then marched through time to attain the steady-state solution. The input to the system is the amplitude of an incident harmonic sound source entering a quiescent duct at the input boundary, with standard impedance boundary conditions on the duct walls and duct exit. The introduction of the time parameter eliminates the large matrix storage requirements normally associated with frequency domain solutions, and time marching attains the steady-state quickly enough to make the method favorable when compared to frequency domain methods. For validation, this transient-frequency domain method is applied to sound propagation in a 2D hard wall duct with plug flow.

Baumeister, K. J.↗
NASA NTRS1999

Some Continuum Aspects of Data Assimilation

A long-sought goal in data assimilation is to build more fully the time dimension into the data assimilation process. Truly four-dimensional data assimilation requires evolving second-moment information, namely the estimation error covariance, along with the state estimate itself. The continuum evolution equations for the estimation error covariance constitute a system of partial differential equations in six space dimensions, forced by a model error covariance and by the observations themselves. The high dimensionality of this system poses a difficult computational problem. Numerous methods have been proposed for approximate, discrete solution of this system of equations. Four-dimensional variational schemes solve these equations indirectly, while schemes based on Kalman filtering solve the equations more directly, usually on a small subspace of the full six-dimensional space. In both cases, most proposed solution methods are not derived from the continuum covariance evolution equations themselves. Instead, most methods simply inherit the discretization used for evolving the state estimate, either in the form of a discrete tangent linear model or, in the case of ensemble schemes, in the form of a fully nonlinear discrete model. In this lecture we show that solutions of the continuum covariance evolution equations possess simple properties that are not readily satisfied upon discretization of the equations. we give numerical examples illustrating that "inherited" discretizations sometimes fail to produce meaningfully accurate solutions. Finally, we suggest alternative discretization methods that may yield more faithful approximate solutions of the continuum problem.

Cohn, Stephen E.↗
NASA NTRS2003

Ultrasail

Ultrasail is a complete sail system for the launch, deployment, stabilization and control of very large solar sails enabling reduced mission times for interplanetary and deep space spacecraft. Ultrasail is an innovative, non-traditional approach to propulsion technology achieved by combining propulsion and control systems developed for formation-flying microsatellites with an innovative solar sail architecture to achieve sq km-class controllable sail areas, sail subsystem area densities of 1 gm per sq m, and thrust levels equivalent to 400 kW ion thruster systems used for comparable deep space missions. Ultrasail can conceivably even achieve outer planetary rendezvous, a deep space capability now reserved for high-mass nuclear and chemical systems. Ultrasail is a Delta IV-launched multi-blade spin-stabilized system with blade lengths as long as 50 km, reminiscent of the MacNeal Heliogyro. The primary innovation is the near-elimination of sail supporting structures by attaching the sail tip to a rigid formation-flying microsatellite truss which deploys the sail blade, and which then articulates the blade to provide attitude control, including spin stabilization and precession of the spin axis. These tip microsatellites are controlled by a solar-powered 3-axis microthruster system (electric or cold gas) to maintain proper sail film tension during deployment and spin-up. The satellite mass also provides a stabilizing centrifugal force on the blade while in rotation. Understanding the dynamics of individual blades is key to the overall dynamics of Ultrasail. Forces and torques that must be modeled include those due to solar pressure, those generated by the microsatellite at the blade tip and by torques applied at the blade root. Centrifugal forces also play a significant role in the deployment and maintenance of the sail configuration. To capture the dynamics of the overall system, the equations of motion for the blades have been derived. Using these differential equations, a control law will be derived to maneuver Ultrasail. This law involves the pitching of the individual blades thereby moving the distribution of the radiation pressure on each individual blade and inducing a resultant torque on the system. The direction of the angular momentum vector and its rate of precession can be controlled through the pitch angle of the blades. The Ultrasail trajectory is also being studied. Optimal or near-optimal trajectories are being generated to showcase Ultrasail performance. Various missions, e.g. outer planet and solar polar missions for observation of the Sun, are currently being investigated to demonstrate the performance enhancements generated by Ultrasail technology. Calculus-of-variations-based optimization software is used to produce optimal Ultrasail trajectories. The performance of these trajectories is being compared to optimal results generated with other propulsion models, including chemical propulsion, ion propulsion, and competing solar sail concepts. Results of these studies will quantify the performance of Ultrasail compared to existing solar sail concepts for high energy missions.

Burton, R.↗
NASA NTRS2011

Gravito-Inertial Force Resolution in Perception of Synchronized Tilt and Translation

Natural movements in the sagittal plane involve pitch tilt relative to gravity combined with translation motion. The Gravito-Inertial Force (GIF) resolution hypothesis states that the resultant force on the body is perceptually resolved into tilt and translation consistently with the laws of physics. The purpose of this study was to test this hypothesis for human perception during combined tilt and translation motion. EXPERIMENTAL METHODS: Twelve subjects provided verbal reports during 0.3 Hz motion in the dark with 4 types of tilt and/or translation motion: 1) pitch tilt about an interaural axis at +/-10deg or +/-20deg, 2) fore-aft translation with acceleration equivalent to +/-10deg or +/-20deg, 3) combined "in phase" tilt and translation motion resulting in acceleration equivalent to +/-20deg, and 4) "out of phase" tilt and translation motion that maintained the resultant gravito-inertial force aligned with the longitudinal body axis. The amplitude of perceived pitch tilt and translation at the head were obtained during separate trials. MODELING METHODS: Three-dimensional mathematical modeling was performed to test the GIF-resolution hypothesis using a dynamical model. The model encoded GIF-resolution using the standard vector equation, and used an internal model of motion parameters, including gravity. Differential equations conveyed time-varying predictions. The six motion profiles were tested, resulting in predicted perceived amplitude of tilt and translation for each. RESULTS: The modeling results exhibited the same pattern as the experimental results. Most importantly, both modeling and experimental results showed greater perceived tilt during the "in phase" profile than the "out of phase" profile, and greater perceived tilt during combined "in phase" motion than during pure tilt of the same amplitude. However, the model did not predict as much perceived translation as reported by subjects during pure tilt. CONCLUSION: Human perception is consistent with the GIF-resolution hypothesis even when the gravito-inertial force vector remains aligned with the body during periodic motion. Perception is also consistent with GIF-resolution in the opposite condition, when the gravito-inertial force vector angle is enhanced by synchronized tilt and translation.

Wood, Scott J.↗
DOE OSTI2025

Towards a Quantum Algorithm for the Incompressible Nonlinear Navier-Stokes Equations

In this work, we present novel concepts for quantum algorithms to solve transient, nonlinear partial differential equations (PDEs). The challenge lies in how to effectively represent, encode, process, and evolve the nonlinear system of PDEs on quantum computers. We will discuss the new techniques using the incompressible Navier-Stokes equations as an example, because it represents the fundamental nonlinear feature and yet removes certain complexity in physics, allowing us to focus on the design of quantum algorithms. Previous attempts solving nonlinear PDEs in quantum computation have often involved storing multiple copies of solutions or employing linearizations. Neither is practical due to exponential scaling with evolution time or insufficient solution accuracy. We propose a new framework based on matrix product states (MPSs) and matrix product operators (MPOs), in addition to the Krylov subspace methods. For example, the solution variables of the Navier-Stokes equations are represented by MPSs, and the linear and nonlinear terms are processed by MPOs. The time evolution of the operators is attained by a fast-forwarding algorithm using Krylov subspace methods. Furthermore, we discuss various techniques for efficient encoding of MPSs, measurement reduction for MPOs, and use of tensor operations to treat multi-variate, multi-physics characteristics of Navier-Stokes.

Gopalakrishnan Meena, Murali [ORNL] (ORCID:0000000↗
DOE OSTI2022

Comprehensive uncertainty quantification (UQ) for full engineering models by solving probability density function (PDF) equation

This report details a new method for propagating parameter uncertainty (forward uncertainty quantification) in partial differential equations (PDE) based computational mechanics applications. The method provides full-field quantities of interest by solving for the joint probability density function (PDF) equations which are implied by the PDEs with uncertain parameters. Full-field uncertainty quantification enables the design of complex systems where quantities of interest, such as failure points, are not known apriori. The method, motivated by the well-known probability density function (PDF) propagation method of turbulence modeling, uses an ensemble of solutions to provide the joint PDF of desired quantities at every point in the domain. A small subset of the ensemble is computed exactly, and the remainder of the samples are computed with approximation of the driving (dynamics) term of the PDEs based on those exact solutions. Although the proposed method has commonalities with traditional interpolatory stochastic collocation methods applied directly to quantities of interest, it is distinct and exploits the parameter dependence and smoothness of the dynamics term of the governing PDEs. The efficacy of the method is demonstrated by applying it to two target problems: solid mechanics explicit dynamics with uncertain material model parameters, and reacting hypersonic fluid mechanics with uncertain chemical kinetic rate parameters. A minimally invasive implementation of the method for representative codes SPARC (reacting hypersonics) and NimbleSM (finite- element solid mechanics) and associated software details are described. For solid mechanics demonstration problems the method shows order of magnitudes improvement in accuracy over traditional stochastic collocation. For the reacting hypersonics problem, the method is implemented as a streamline integration and results show very good accuracy for the approximate sample solutions of re-entry flow past the Apollo capsule geometry at Mach 30.

42 ENGINEERING↗