Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical differentiation”

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 793 records · Page 44

An improved numerical process for solution of solid mechanics problems

This paper gives an overview of the development and status of an improved numerical process for the solution of solid mechanics problems. The proposed process uses a mixed formulation with the fundamental unknowns consisting of both stress and displacement parameters. The problem is formulated either by means of first-order partial differential equations or in a variational form by using a Hellinger-Reissner-type mixed variational principle. For presentation purposes, the components of a numerical process are characterized and the criteria for an ideal process are outlined. Commonly used finite-difference and finite-element procedures are examined in the light of these criteria and it is shown that they fall short in a number of ways. The proposed numerical process, on the other hand, satisfies most of the optimality criteria and appears to be particularly suited for use with the forthcoming generation computers.

Noor, A. K.↗

Ephemeris of a highly eccentric orbit - Explorer 28.

An ephemeris has been obtained for Explorer 28 (IMP 3) which agrees well with 2 years of radio observations and with SAO observations a year later. This ephemeris is generated over the 3 year lifetime by a numerical integration method utilizing a set of initial conditions at launch and without requiring further differential correction. Because highly eccentric orbits are difficult to compute with acceptable accuracy and because a long continuous arc has been obtained which compares with actual data to a known precision, this ephemeris may be used as a standard for computing highly eccentric orbits in the earth-moon system. Orbit improvement was used to obtain the initial conditions which generated the ephemeris. This improvement was based on correcting the energy by adjusting the semimajor axis to match computed times of perigee passage with the observed. This procedure may generate errors in semimajor axis to compensate for model errors in the energy; however, this compensation error is also implicit in orbit determination itself.

Lowrey, B. E.↗

L-orthogonal signal transmission and detection.

Investigation of the detailed capabilities and performance characterization of systems that employ L-orthogonal signaling techniques. L-orthogonal signals represent a unified set of signals wherein the polyphase and orthogonal (bi-orthogonal) signal sets are included as special cases. This fact is important, since orthogonal (bi-orthogonal) and polyphase signaling sets represent opposing forces as far as tradeoffs between error probability, energy-to-noise ratio, and bandwidth requirements are concerned. Bounds on the performance of the optimum receiver and the performance of various suboptimum (practical) receiver structures are given. Coherent and differentially coherent detection of differentially encoded signals are also pursued. Various comparisons and tradeoffs are made by means of numerical evaluation of the error probability expressions.

Lindsey, W. C.↗

The stability of motion of satellites with cavities partially filled with liquid

The stability and time dependent motion of a spinning satellite, simulated by a rigid body with a cavity partially filled with liquid is examined. The problem formulation, consisting of the boundary-value problem for the liquid and moment equations for the entire system is presented. Because of large Reynold's numbers involved, viscosity effects are negligible everywhere except for a thin boundary layer near the wetted surface. Using a boundary-layer analysis, the effect of the boundary layer is replaced by modified boundary conditions for the liquid. The solution of the differential equations for the inviscid problem is solved in closed form. A semi-analytical numerical solution of the inviscid equations subject to the viscous boundary condition has proved unsucessful.

Nayfeh, A. H.↗

Symbolic-numeric interface: A review

A survey of the use of a combination of symbolic and numerical calculations is presented. Symbolic calculations primarily refer to the computer processing of procedures from classical algebra, analysis, and calculus. Numerical calculations refer to both numerical mathematics research and scientific computation. This survey is intended to point out a large number of problem areas where a cooperation of symbolic and numerical methods is likely to bear many fruits. These areas include such classical operations as differentiation and integration, such diverse activities as function approximations and qualitative analysis, and such contemporary topics as finite element calculations and computation complexity. It is contended that other less obvious topics such as the fast Fourier transform, linear algebra, nonlinear analysis and error analysis would also benefit from a synergistic approach.

Ng, E. W.↗

Computation of the modes and polar moment of inertial of the blades of an HAWT

The coupled differential equations of motion of the blades of a horizontal axis wind turbine are solved numerically, permitting the optimization of the design at relatively low cost. The equation of motion is transformed into a set of first order equations and solved with fourth order Runge-Kutta integrators. This technique is applied to a twisted, tapered blade of variable cross section and stiffness including discontinuities. The first six natural frequencies and mode shapes are obtained. The polar moment of inertia of the blades is obtained as a function of frequency and rotational speed.

Beaulieu, G.↗

Methods of sequential estimation for determining initial data in numerical weather prediction

Numerical weather prediction (NWP) is an initial-value problem for a system of nonlinear differential equations, in which initial values are known incompletely and inaccurately. Observational data available at the initial time must therefore be supplemented by data available prior to the initial time, a problem known as meteorological data assimilation. A further complication in NWP is that solutions of the governing equations evolve on two different time scales, a fast one and a slow one, whereas fast scale motions in the atmosphere are not reliably observed. This leads to the so called initialization problem: initial values must be constrained to result in a slowly evolving forecast. The theory of estimation of stochastic dynamic systems provides a natural approach to such problems. For linear stochastic dynamic models, the Kalman-Bucy (KB) sequential filter is the optimal data assimilation method, for linear models, the optimal combined data assimilation-initialization method is a modified version of the KB filter.

Cohn, S. E.↗

Radiative Amplification of Acoustic Waves in Hot Stars

The discovery of broad P Cygni profiles in early type stars and the detection of X-rays emitted from the envelopes of these stars made it clear, that a considerable amount of mechanical energy has to be present in massive stars. An attack on the problem, which has proven successful when applied to late type stars is proposed. It is possible that acoustic waves form out of random fluctuations, amplify by absorbing momentum from stellar radiation field, steepen into shock waves and dissipate. A stellar atmosphere was constructed, and sinusoidal small amplitude perturbations of specified Mach number and period at the inner boundary was introduced. The partial differential equations of hydrodynamics and the equations of radiation transfer for grey matter were solved numerically. The equation of motion was augmented by a term which describes the absorption of momentum from the radiation field in the continuum and in lines, including the Doppler effect and allows for the treatment of a large number of lines in the radiative acceleration term.

Wolf, B. E.↗

Quasi-geostrophic free mode models of long-lived Jovian eddies: Forcing mechanisms and crucial observational tests

Observations of Jupiter and Saturn long-lived eddies, such as Jupiter's Great Red Spot and White Ovals, are presently compared with laboratory experiments and corresponding numerical simulations for free thermal convection in a rotating fluid that is subject to horizontal differential heating and cooling. Difficulties in determining the essential processes maintaining and dissipating stable eddies, on the basis of global energy budget studies, are discussed; such difficulties do not arise in considerations of the flow's potential vorticity budget. On Jupiter, diabatically forced and transient eddy-driven flows primarily differ in the implied role of transient eddies in transporting potential vorticity across closed geostrophic streamlines in the time mean.

Read, P. L.↗

A second-order modeling study of confined swirling flow

A numerical study of a confined strong swirling flow is presented. Computations are performed using a differential second-moment closure. The effect of inlet dissipation rate on calculated mean and turbulence fields is investigated. Two pressure-strain models with various model constants are examined to demonstrate their influences on the predicted results. Finally, comparison of the differential second-moment calculation is made with algebraic second-moment results to determine a better suited model for complex swirling flows.

Nikjooy, M.↗

Comparison Of Two Viscous-Flow Computer Codes

Two viscous-flow computer codes compared by applying them to five test cases of steady-state transonic viscous flows about transonic airfoils. Two codes were: FLOMG, which solves Navier-Stokes differential equations of flow by implementing explicit, Runge-Kutta, finite-volume, multigrid numerical-integration algorithm; and ARC2D, which implements implicit, finite-difference, approximate-factorization, eigenvector-diagonalization numerical-integration algorithm.

Maksymiuk, C. M.↗

On the nonlinear stability of the unsteady, viscous flow of an incompressible fluid in a curved pipe

The stability of the flow of an incompressible, viscous fluid through a pipe of circular cross-section curved about a central axis is investigated in a weakly nonlinear regime. A sinusoidal pressure gradient with zero mean is imposed, acting along the pipe. A WKBJ perturbation solution is constructed, taking into account the need for an inner solution in the vicinity of the outer bend, which is obtained by identifying the saddle point of the Taylor number in the complex plane of the cross-sectional angle co-ordinate. The equation governing the nonlinear evolution of the leading order vortex amplitude is thus determined. The stability analysis of this flow to periodic disturbances leads to a partial differential system dependent on three variables, and since the differential operators in this system are periodic in time, Floquet theory may be applied to reduce this system to a coupled infinite system of ordinary differential equations, together with homogeneous uncoupled boundary conditions. The eigenvalues of this system are calculated numerically to predict a critical Taylor number consistent with the analysis of Papageorgiou. A discussion of how nonlinear effects alter the linear stability analysis is also given, and the nature of the instability determined.

Shortis, Trudi A.↗

Dynamics of emerging active region flux loops

The buoyant rise of a magnetic flux loop arising from a single perturbed segment of a toroidal flux ring lying slightly beneath the base of the convection zone is studied by way of numerical simulations. We have considered flux loop evolution assuming both solid-body rotation, and differential rotation consistent with recent results from helioseismology. Our major results are presented, and we offer some speculations on the decay of active regions, based on the results of our studies. We speculate that as plasma in the tube attempts to establish hydrostatic equilibrium along the field lines after the flux emergence has taken place, the tube field strength at some intermediate depths below the surface becomes sufficiently small at the surface portions of the tube (which have cooled and undergone convective collapse) become dynamically disconnected from those portions near the base of the convection zone. The surface proportions of the emerged flux tubes are then transported by motions near the photosphere, such as supergranular convection and meridional flow.

Fan, Y.↗

Temperature structure of active regions deduced from helium-like sulphur lines

Solar active-region temperatures have been determined from the full-Sun spectra of helium-like sulfur (S XV) observed by the Bragg Crystal Spectrometer on Board the Yohkoh satellite. The average temperature deduced from S XV is demonstrated to vary with the solar activity level: A temperature of 2.5 x 10(exp 6) K is derived from the spectra taken during low solar activity, similar to the general corona, while 4 x 10(exp 6) K is obtained during a higher activity phase. For the latter, the high- temperature tail of the differential emission measure of active regions is found most likely due to the superposition of numerous flare-like events (micro/nano-flares).

Watanabe, Tetsuya↗

Experimental Study of Large-Amplitude Faraday Waves in Rectangular Cylinders

Experiment on single-mode Faraday waves having two, thee, and four wavelengths across a rectangular cylinder of high aspect ratio is the subject of discussion. Previous experiments recently done by Henderson & Miles (1989) and by Lei Jiang et. a1 (1996) focused on Faraday waves with one and two wavelengths across rectangular cylinders. In this experimental study the waves steepness ranges from small at threshold levels to a large amplitude which according to Penny & Price theory (1952) approaches the maximum sustainable amplitude for a standing wave. The waves characteristics for small amplitudes are evaluated against an existing well known linear theory by Benjamin & Ursell (l954) and against a weakly nonlinear theory by J. Miles (1984) which includes the effect of viscous damping. The evaluation includes the wave neutral stability and damping rate. In addition, a wave amplitude differential equation of a linear theory including viscous effect by Cerda & Tirapegui (1998) is solved numerically to yield prediction of temporal profiles of both wave damping and wave formation at the threshold. An interesting finding from this exercise is that the fluid kinematic viscosity needs to increase ten times in order to obtain good agreement between the theoretical prediction and the experimental data for both wave damping and wave starting. For large amplitude waves, the experimental data are evaluated against the theory of Penny & Price which predicts wave characteristics of any amplitude up to the point at which the wave reaches its maximum amplitude attainable for a standing wave. The theory yields two criteria to show the maximum wave steepness, the vertical acceleration at the wave crest of half the earth gravity field acceleration and the including angle at the crest of 90 degrees. Comparison with experimental data shows close agreement for the wave crest acceleration but a large discrepancy for the including angle. Additional information is included in the original extended abstract.

Iek, Chanthy↗

Towards Formalization of Advanced Linear Algebra with Applications to Dynamical Systems using PVS

Linear Algebra is essential for numerous aerospace problems of interest. Formal reasoning about hybrid systems that contain variables modeled by differential equations rely on concepts from Linear Algebra such as eigenvalues, matrix decompositions, and matrix valued functions. For example, the long-term dynamics of a system of differential equations depend on the stability/instability of its equilibrium points, which often reduces to an eigenvalue problem. This talk will embark on a quest to formalize theorems and results about eigenvalues and eigenvectors using PVS. We shall start our journey with 2 x 2 complex matrices, where we will apply our PVS code to a simple example of a dynamical system. Since it can be difficult or impossible to give simple expressions of eigenvalues for larger matrices (i.e. 5 x 5 or higher), we then move towards specifying the power method for verified computation of eigenvalue approximations in PVS. This effort requires development of multivariate complex arithmetic. At the end of the day, having such additions to the PVS NASA libraries will help move towards the use of formal methods to verify concepts of control theory and system level verification.

Linear Algebra↗

Communication-Aware Orbit Design for Small Spacecraft Swarms around Small Bodies

Exploration of small Solar System bodies has traditionally been performed by single monolithic spacecraft carrying a number of science instruments. However, science instruments typically cannot be operated simultaneously due to the instrument requirements including optimal viewing angle, surface illumination, altitude and ground resolution, power, and data constraints. This observation has motivated interest in multi-spacecraft architectures where a swarm of small spacecraft, each carrying a single science instrument, studies a small body after being deployed by a carrier spacecraft, which then collects data from the vehicles and relays it to Earth. Such architectures hold promise to yield significant improvements in mission efficiency, increases in data quality, and shorter mission duration. A key difficulty in the design of such missions is the selection of orbits for the small spacecraft, which must satisfy not only instrument requirements, but also strict inter-spacecraft communication and on-board storage constraints. To address this, in this paper, we present a novel computationally-efficient optimization algorithm for \emph{communication-aware design} of the orbits of a small spacecraft swarm orbiting a small body. The proposed approach captures constraints including instrument requirements, inter-spacecraft communication bandwidths, and on-board memory usage, and it can accommodate highly irregular gravity field models and surface geometries. We propose an efficient algorithm for optimization of instrument observations and inter-spacecraft communications; we then leverage the differentiable nature of the proposed algorithm to accelerate a gradient-based global search algorithm. Numerical simulations of a six-spacecraft swarm studying 433 Eros show that the proposed approach successfully identifies high-quality orbits, and significantly outperform communication-agnostic optimization techniques, resulting in a 10% increase in scientific returns and a 30% increase in the quality of the collected data.

Rahmani, Amir↗

Flow Applications of the Least Squares Finite Element Method

The main thrust of the effort has been towards the development, analysis and implementation of the least-squares finite element method (LSFEM) for fluid dynamics and electromagnetics applications. In the past year, there were four major accomplishments: 1) special treatments in computational fluid dynamics and computational electromagnetics, such as upwinding, numerical dissipation, staggered grid, non-equal order elements, operator splitting and preconditioning, edge elements, and vector potential are unnecessary; 2) the analysis of the LSFEM for most partial differential equations can be based on the bounded inverse theorem; 3) the finite difference and finite volume algorithms solve only two Maxwell equations and ignore the divergence equations; and 4) the first numerical simulation of three-dimensional Marangoni-Benard convection was performed using the LSFEM.

Jiang, Bo-Nan↗