Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PERIODIC FUNCTION”

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 37 records · Page 2

Aeroelastic stability of periodic systems with application to rotor blade flutter

The dynamics of a helicopter blade in forward flight are described by a system of linear differential equations with periodic coefficients. The stability of this periodic aeroelastic system is determined, using multivariable Floquet-Liapunov theory. The transition matrix at the end of the period is evaluated by: (1) direct numerical integration, and (2) a new, approximate method, which consists in approximating a periodic function by a series of step functions. The numerical accuracy and efficiency of the methods is compared, and the second method is shown to be superior by far. Results illustrating the effect of the periodic coefficients and various blade parameters are presented.

Friedmann, P.↗

Quadrature methods for periodic singular and weakly singular Fredholm integral equations

High-accuracy numerical quadrature methods for integrals of singular periodic functions are proposed. These methods are based on the appropriate Euler-Maclaurin expansions of trapezoidal rule approximations and their extrapolations. They are subsequently used to obtain accurate quadrature methods for the solution of singular and weakly singular Fredholm integral equations. Throughout the development the periodic nature of the problem plays a crucial role. Such periodic equations are used in the solution of planar elliptic boundary value problems such as those that arise in elasticity, potential theory, conformal mapping, and free surface flows. The use of the quadrature methods is demonstrated with numerical examples.

Sidi, Avram↗

"Dust Devils": Gardening Agents on the Surface of Mars, and Hidden Hazards to Human Exploration?

Dust devils are familiar sites in the and regions of the world: they can produce quite spectacular displays of dust lofting when the vortices scavenge very loose dust from a dry lake bed or from recently disturbed agricultural fields. If one were to arrive at the center of an arid region, take one photograph, or even a series of photographs over a period of several days, then return the images for laboratory analysis, it would be most likely concluded that the region was inactive from an aeolian perspective. No images of general dust movement were obtained, nor were any dust devils "caught on camera" owing to their ephemeral and unpredictable appearance, and the fact that there was deceptively little residue of their actions. If, however, a camera were to take a 360 degree continuous recording over a period of a year, and the film were then to be shown at high speed over a period a several minutes, the impression might be that of a region ravaged by air vorticity and dust movement. Extrapolate this over geological time, and it is possible to visualize dust devils as prime aeolian agents, rather than insignificant vagaries of nature, On Mars, the thin atmosphere permits the surface of the planet to be heated but it does not itself retain heat with the capacity of the earth's atmosphere. This gives rise to greater thermal instability near the surface of Mars as "warm" air pockets diapiritically inject themselves into higher atmospheric layers. Resulting boundary-layer vorticity on Mars might therefore be expected to produce dust devils in abundance, if only seasonally. The spectacular images of dust devils obtained by Pathfinder within its brief functional period on the planet testify to the probability of highly frequent surface vorticity in light of the above reasoning about observational probability. Notably, the Pathfinder devils appeared to be at least a kilometer in height. There are several consequences for the geology of Mars, and for human exploration, if dust devils are to be expected in reasonable abundance. First, from a geological perspective, the vortices will act as "gardening" agents for the top few centimeters of entrainable material. Over time (hundreds of millions, or billions of years being available), they will cover the surface with scouring paths, and the grain sizes that can be lofted by a vortex probably extends over the whole sand to dust range. The depositional paths are, of course, much larger, so that vortex-induced deposition is more widespread than vortex-induced erosion, and will without doubt, affect the whole region in which the dust devils occur (this might explain why rocks at the Viking site seemed oddly capped with dust in a region apparently subject to general aeolian scouring). On Mars, the lift forces in dust devils might be less than on earth owing to the much thinner atmosphere, but this may be counterbalanced by lower gravity and greater vortex velocities. Certainly, when active, other aeolian phenomena on Mars --sand motion and dust storms, seem no less energetic and no less capable of lofting sediments than equivalent terrestrial aeolian phenomena. Every several years, within the current climatic regime, the surface of Mars is subject to light dust fall from global dust storms. Over time, this should develop a very uniform surface layer, with commensurate uniformity in grain size, mineralogy, albedo, color, and general spectroscopic properties. Dust devils will disturb this situation by continually mixing the surface dust with underlying layers, perhaps composed of silt and sand. This size mixing will also involve compositional mixing. After some years, the thin layer of dust that may be difficult to entrain alone, becomes progressively mixed with coarser materials that could reduce the general aeolian threshold of the soil. Certainly the continual disturbance by vorticity will prevent surface stabilization that may bind or indurate grains (caused by slow cementation or ice welding at grain boundaries). If dust devils continually loft dust to kilometer heights, and the dust is sprayed into many cubic kilometers of atmosphere each time, could the devils produce a continual background of atmospheric dust that might be mistaken for the fallout of a distant large-scale dust storm? From a human exploration perspective, dust devils are unlikely to pose any, life- threatening situation for an astronaut unfortunate enough to encounter a momentary swirling cloud of loose soil. However, it is noted that pervasive dust is probably one of the greatest long-term hazards for a human encampment. The fineness and penetration capabilities of the dust, its electrostatic adhesive properties, and its complete ubiquity, render the material a persistent nuisance at best, but at worst, over a period of many months it is possible that space suits, machinery, habitat interiors, air filters, and so forth, could become jeopardized. Owing to dust penetration, the space suits used in the Apollo landings were rendered unusable after a few EVA activities. There will be a definite attempt to situate a human colony on Mars in an area that is far removed from the regions of the planet known for being the centers of major dust storms. At the heart of these storm systems, the dust lofting mechanics are unknown, but they are energetic and perhaps potentially life-threatening for an astronaut. Locating a colony in a region that appears from space to be meteorologically benign may lead to colony placement in a region prone to dust devils, but dust devils are not (or have not been) detectable from orbital observations: the region surveyed for placement will appear like the apparently inactive and area referred to earlier. The region may be spared from highly energetic weather systems, but it may not be necessarily immune from continual dust disturbance. Additional information is contained in the original.

Marshall, J.↗

Quadrature solution for the general relativistic motion of a satellite or a planet.

The general relativistic motion of a particle in a bounded orbit around a central mass M is derived with the aid of the Schwarzschild field of M. A quadrature gives the central angle as a quasi-periodic function of an effective true anomaly. The mean rate of advance of the pericenter is investigated. Classical and general relativistic orbiters with the same initial conditions are compared, taking into account the treatment of the general relativistic effect by means of a differential correction. General relativistic variations are discussed, giving attention to the variation of the eccentric anomaly and the variation of the central angle.

Vinti, J. P.↗

Real-time computer data system for the 40- by 80-foot wind tunnel facility at Ames Research Center

The background material and operational concepts of a computer-based system for an operating wind tunnel are described. An on-line real-time computer system was installed in a wind tunnel facility to gather static and dynamic data. The computer system monitored aerodynamic forces and moments of periodic and quasi-periodic functions, and displayed and plotted computed results in real time. The total system is comprised of several off-the-shelf, interconnected subsystems that are linked to a large data processing center. The system includes a central processor unit with 32,000 24-bit words of core memory, a number of standard peripherals, and several special processors; namely, a dynamic analysis subsystem, a 256-channel PCM-data subsystem and ground station, a 60-channel high-speed data acquisition subsystem, a communication link, and static force and pressure subsystems. The role of the test engineer as a vital link in the system is also described.

Cambra, J. M.↗

Theory of orthodontic motions

A general theory of orthodontic motion is developed that can be applied to determine the forces necessary to induce a given tooth to move to the predetermined desirable position. It is assumed that the natural (nonorthodontic) forces may be represented by a periodic function and the orthodontic forces may be superimposed upon the natural forces. A simple expression is derived for the applied stress.

Pepe, S.↗

The response of a constant-volume balloon to periodic three-dimensional flow

A theoretical treatment is presented of the motion of a constant-volume balloon immersed in a three-dimensional spatially nonhomogeneous periodic flow field with a mean translational flow velocity. The governing equations are the equations of conservation of momentum of the balloon, coupled with the equations of conservation of momentum of the fluid. In the mathematical model employed, the flow field is represented by a periodic function. The numerical results are in general agreement with observations. Noteworthy is that the balloon does not move with the same mean velocity as the wind and, thus, when Taylor's hypothesis for the fluid is satisfied, the balloon does not 'lock in' with the wind velocity. When Taylor's hypothesis for the parcel is satisfied, however, the balloon does match the mean wind direction (but not its speed).

Tatom, F. B.↗

Method for predicting impulsive noise generated by wind turbine rotors

Large wind turbines can generate both broad band and impulsive noises. These noises can be controlled by proper choice of rotor design parameters such as rotor location with respect to the supporting tower, tower geometry and tip speed. A method was developed to calculate the impulsive noise generated when the wind turbine blade experiences air forces that are periodic functions of the rotational frequency. This phenomenon can occur when the blades operate in the wake of the support tower and the nonuniform velocity field near the ground due to wind shear. Results from this method were compared with measured sound spectra taken at locations of one to two rotor diameters from the DOE/NASA Mod-1 wind turbine. The calculated spectra generally agreed with the measured data in both the amplitude of the predominant harmonics and the roll of rate with frequency. Measured sound pressure levels far from the Mod-1 (15 rotor diameters), however, were higher than predicted. Simultaneous measurements in the near and far field indicated that the propagation effects could enhance the sound levels by more than 10 dB above that expected by spherical dispersion. These propagation effects are believed to be due to terrain and atmospheric characteristics of the Mod-1 site.

Viterna, L. A.↗

Design of a rotary stepped auger for a lunar environment

A lunar outpost will have need for deep drilling operations for both explorative and practical purposes. As in any drilling operation, the cuttings must be cleared from the hole. The hard vacuum of the lunar environment renders conventional flushing methods of cutting removal unfeasible, and requires a new system of removal. A rotary stepped auger (RSA) is a simple mechanical method of removing dry cuttings from a deep hole, and is ideally suited to the lunar environment. The RSA consists of a helical auger with stepped ramps which allow cuttings to slide up the helix, but will prevent them from sliding back down. The auger is driven in a pulsed manner by applying a periodic function of acceleration to the auger shaft. These pulses will compel the cuttings to slide up the auger's helix while the stepped ramps prevent the cuttings from backsliding while the auger accelerates. A mathematical model of the RSA was developed and experimentally evaluated. The math model produced a good baseline design, but the experimental model required some tuning to account for the approximations made in the math model. This design is suited for lunar drilling because it is mechanically simple, integral to the drill string, requires no fluids, is suited to the dry soil, and has relatively low weight and power requirements.

Dardet, Eduardo↗

A perturbation method and some applications

For differential equations with one fast variable, a perturbation method is introduced that transforms a solution valid over only a short time interval to a new solution composed of averaged variables plus a periodic function of the averaged variables. The averaged variables are governed by a set of differential equations where the fast variable has been removed and thus can be numerically integrated quickly or solved directly. This method is applied to a perturbed harmonic oscillator with a cubic perturbation, van der Pol's equation, coorbital motion in the restricted three-body problem, and to nearly circular motion of a particle near one of the primaries in the restricted three-body problem.

Konopliv, Alex↗

Rational trigonometric approximations using Fourier series partial sums

A class of approximations (S(sub N,M)) to a periodic function f which uses the ideas of Pade, or rational function, approximations based on the Fourier series representation of f, rather than on the Taylor series representation of f, is introduced and studied. Each approximation S(sub N,M) is the quotient of a trigonometric polynomial of degree N and a trigonometric polynomial of degree M. The coefficients in these polynomials are determined by requiring that an appropriate number of the Fourier coefficients of S(sub N,M) agree with those of f. Explicit expressions are derived for these coefficients in terms of the Fourier coefficients of f. It is proven that these 'Fourier-Pade' approximations converge point-wise to (f(x(exp +))+f(x(exp -)))/2 more rapidly (in some cases by a factor of 1/k(exp 2M)) than the Fourier series partial sums on which they are based. The approximations are illustrated by several examples and an application to the solution of an initial, boundary value problem for the simple heat equation is presented.

Geer, James F.↗

A spectral boundary integral equation method for the 2-D Helmholtz equation

In this paper, we present a new numerical formulation of solving the boundary integral equations reformulated from the Helmholtz equation. The boundaries of the problems are assumed to be smooth closed contours. The solution on the boundary is treated as a periodic function, which is in turn approximated by a truncated Fourier series. A Fourier collocation method is followed in which the boundary integral equation is transformed into a system of algebraic equations. It is shown that in order to achieve spectral accuracy for the numerical formulation, the nonsmoothness of the integral kernels, associated with the Helmholtz equation, must be carefully removed. The emphasis of the paper is on investigating the essential elements of removing the nonsmoothness of the integral kernels in the spectral implementation. The present method is robust for a general boundary contour. Aspects of efficient implementation of the method using FFT are also discussed. A numerical example of wave scattering is given in which the exponential accuracy of the present numerical method is demonstrated.

Hu, Fang Q.↗

Annual and Semi-Annual Temperature Oscillations in the Upper Mesosphere

Fourier transform spectrometer observations of the mesosphere have been performed at the University of Michigan (latitude: 42.5 N) on a long term basis. A database of near infrared Meinel hydroxyl spectra has been accumulated from which rotational temperatures have been determined. Harmonic analysis of one-day averaged temperatures for the period 1992.0 to 1994.5 has shown a distinct annual and semi-annual variation. Subsequent fitting of a five term periodic function characterizing the annual and semi-annual temperature oscillations to the daily averaged temperatures was performed. The resultant mean temperature and the amplitudes and phases of the annual and semi-annual variations are shown to coincide with an emission height slightly above 85 km which is consistent with the mean rocket derived altitude for peak nocturnal hydroxyl emission.

Niciejewski, R. J.↗

Exponential Approximations Using Fourier Series Partial Sums

The problem of accurately reconstructing a piece-wise smooth, 2(pi)-periodic function f and its first few derivatives, given only a truncated Fourier series representation of f, is studied and solved. The reconstruction process is divided into two steps. In the first step, the first 2N + 1 Fourier coefficients of f are used to approximate the locations and magnitudes of the discontinuities in f and its first M derivatives. This is accomplished by first finding initial estimates of these quantities based on certain properties of Gibbs phenomenon, and then refining these estimates by fitting the asymptotic form of the Fourier coefficients to the given coefficients using a least-squares approach. It is conjectured that the locations of the singularities are approximated to within O(N(sup -M-2), and the associated jump of the k(sup th) derivative of f is approximated to within O(N(sup -M-l+k), as N approaches infinity, and the method is robust. These estimates are then used with a class of singular basis functions, which have certain 'built-in' singularities, to construct a new sequence of approximations to f. Each of these new approximations is the sum of a piecewise smooth function and a new Fourier series partial sum. When N is proportional to M, it is shown that these new approximations, and their derivatives, converge exponentially in the maximum norm to f, and its corresponding derivatives, except in the union of a finite number of small open intervals containing the points of singularity of f. The total measure of these intervals decreases exponentially to zero as M approaches infinity. The technique is illustrated with several examples.

Banerjee, Nana S.↗

X-38 V201 Avionics Architecture

The X-38 is an experimental NASA project developing a core human capable spacecraft at a fraction of the cost of any previous human rated vehicle. The first operational derivative developed from the X-38 program will be the International Space Station (ISS) Crew Return Vehicle (CRV). Although the current X-38 vehicles are designed as re-entry vehicles only, the option exists to modify the vehicle for uses as an upward vehicle launched from an expendable launch vehicle or from the X-33 operational derivative. The Operational CRV, that will be derived from the X-38 spaceflight vehicle, will provide an emergency return capability from the International Space Station (ISS). The spacecraft can hold a crew of up to seven inside a pressurized cabin. The CRV is passively delivered to ISS, stays up to three year on-orbit attached to ISS in a passive mode with periodic functional checkout, before separation from ISS, de-orbit, entry and landing. The X-38 Vehicle 201 (V201) is being developed at NASA/JSC to demonstrate key technologies associated with the development of the CRV design. The X-38 flight test will validate the low cost development concept by demonstrating the entire station departure, re-entry, guidance and landing portions of the CRV mission. All new technologies and subsystems proposed for CRV will be validated during either the on orbit checkout or flight phases of the X-38 space flight test. The X-38 subsystems are required to be similar to those subsystems required for the CRV to the greatest extent possible. In many cases, the subsystems are identical to those that will be utilized on the Operational CRV.

Bedos, Thierry↗

QRCODE: Massively parallelized real-time time-dependent density functional theory for periodic systems

We present a new software module, QRCODE (Quantum Research for Calculating Optically Driven Excitations), for massively parallelized real-time time-dependent density functional theory (RT-TDDFT) calculations of periodic systems in the open-source Qbox software package. Our approach utilizes a custom implementation of a fast Fourier transformation scheme that significantly reduces inter-node message passing interface (MPI) communication of the major computational kernel and shows impressive scaling up to 16,344 CPU cores. In addition to improving computational performance, QRCODE contains a suite of various time propagators for accurate RT-TDDFT calculations. As benchmark applications of QRCODE, we calculate the current density and optical absorption spectra of hexagonal boron nitride (h-BN) and photo-driven reaction dynamics of the ozone-oxygen reaction. We also calculate the second and higher harmonic generation of monolayer and multi-layer boron nitride structures as examples of large material systems. Our optimized implementation of RT-TDDFT in QRCODE enables large-scale calculations of real-time electron dynamics of chemical and material systems with enhanced computational performance and impressive scaling across several thousand CPU cores.

97 MATHEMATICS AND COMPUTING↗