Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Oscillation”

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 181 records · Page 10

Transonic airfoil analysis and design in nonuniform flow

A nonuniform transonic airfoil code is developed for applications in analysis, inverse design and direct optimization involving an airfoil immersed in propfan slipstream. Problems concerning the numerical stability, convergence, divergence and solution oscillations are discussed. The code is validated by comparing with some known results in incompressible flow. A parametric investigation indicates that the airfoil lift-drag ratio can be increased by decreasing the thickness ratio. A better performance can be achieved if the airfoil is located below the slipstream center. Airfoil characteristics designed by the inverse method and a direct optimization are compared. The airfoil designed with the method of direct optimization exhibits better characteristics and achieves a gain of 22 percent in lift-drag ratio with a reduction of 4 percent in thickness.

Chang, J. F.↗

Essentially nonoscillatory spectral Fourier methods for shock wave calculations

An essentially nonoscillatory spectral Fourier method for the solution of hyperbolic partial differential equations is presented. The method is based on adding a nonsmooth function to the trigonometric polynomials which are the usual basis functions for the Fourier method. The high accuracy away from the shock is enhanced by using filters. Numerical results confirm that essentially no oscillations develop in the solution.

Cai, Wei↗

Satellite Altimetry, Ocean Circulation, and Data Assimilation

Ocean circulation is a critical factor in determining the Earth's climate. Satellite altimetry has been proven a powerful technique for measuring the height of the sea surface for the study of global ocean circulation dynamics. A major objective of my research is to investigate the utility of altimeter data for ocean circulation studies. The 6 years' data record of TOPEX/POSEIDON have been analyzed to study the spatial and temporal characteristics of large-scale ocean variability. A major result obtained in 1998 is the discovery of large-scale oscillations in sea level with a period of 25 days in the Argentine Basin of the South Atlantic Ocean (see diagram). They exhibit a dipole pattern with counterclockwise rotational propagation around the Zapiola Rise (centered at 45S and 317E), a small seamount in the abyssal plain of the basin. The peak-to-trough amplitude is about 10 cm over a distance of 500-1000 km. The amplitude of these oscillations has large seasonal-to-interannual variations. The period and rotational characteristics of these oscillations are remarkably similar to the observations made by two current meters deployed near the ocean bottom in the region. What TOPEX/POSEIDON has detected apparently are manifestations of the movement of the entire water column (barotropic motion). The resultant transport variation is estimated to be about 50 x 10(exp 6) cubic M/S, which is about 50% of the total water transport in the region. Preliminary calculations suggest that these oscillations are topographically trapped waves. A numerical model of the South Atlantic is used to investigate the nature of and causes for these waves. A very important property of sea surface height is that it is directly related to the surface geostrophic velocity, which is related to deep ocean circulation through the density field. Therefore altimetry observations are not only useful for determining the surface circulation but also for revealing information about the deep ocean. Another objective of my research is to investigate the extent of information one can infer about the deep ocean from altimetric observations. The approach is data assimilation by ocean models. The objective is to develop a modeling/data assimilation system that will produce estimates of the 3-dimensional state of the ocean for the entire duration of the TOPEX/POSEIDON Mission.

Fu, Lee-Lueng↗

Planetary-Scale Inertio Gravity Waves in the Numerical Spectral Model

In the polar region of the upper mesosphere, horizontal wind oscillations have been observed with periods around 10 hours. Waves with such a period are generated in our Numerical Spectral Model (NSM), and they are identified as planetary-scale inertio gravity waves (IGW). These IGWs have periods between 9 and 11 hours and appear above 60 km in the zonal mean (m = 0), as well as in zonal wavenumbers m = 1 to 4. The waves can propagate eastward and westward and have vertical wavelengths around 25 km. The amplitudes in the wind field are typically between 10 and 20 m/s and can reach 30 m/s in the westward propagating component for m = 1 at the poles. In the temperature perturbations, the wave amplitudes above 100 km are typically 5 K and as large as 10 K for m = 0 at the poles. The IGWs are intermittent but reveal systematic seasonal variations, with the largest amplitudes occurring generally in late winter and spring. In the NSM, the IGW are generated like the planetary waves (PW). They are produced apparently by the instabilities that arise in the zonal mean circulation. Relative to the PWs, however, the IGWs propagate zonally with much larger velocities, such that they are not affected much by interactions with the background zonal winds. Since the IGWs can propagate through the mesosphere without much interaction, except for viscous dissipation, one should then expect that they reach the thermosphere with significant and measurable amplitudes.

Mayr, H. G.↗

Numerical simulation of unsteady viscous flows

Most unsteady viscous flows may be grouped into two categories, i.e., forced and self-sustained oscillations. Examples of forced oscillations occur in turbomachinery and in internal combustion engines while self-sustained oscillations prevail in vortex shedding, inlet buzz, and wing flutter. Numerical simulation of these phenomena was achieved due to the advancement of vector processor computers. Recent progress in the simulation of unsteady viscous flows is addressed.

Hankey, Wilbur L.↗

The obliquity of Mars and 'climate friction'

A mathematical theory is presented which explains the increase of the mean obliquity of Mars over geologic time due to the so called 'climate friction' (i.e., the climatic changes associated with obliquity oscillations of Mars). The theory is compared with a 10 m.y. numerical integration of the equations performed for a hypothetically large amount of climate friction for two cases of the obliquity oscillations: (1) a single sinusoid and (2) a sum of three sinusoids. The theory and numerics agree for both cases within about 12 percent on the size of the secular increase in obliquity. One possible mechanism of climate friction investigated is 'postglacial rebound' on Mars. According to this theory, giant polar caps form when the obliquity is low, and slowly squeeze out an equatorial bulge. When the obliquity is high, the caps disappear, but the bulge takes some time to collapse, due to mantle viscosity, causing it to oscillate out of phase with the obliquity oscillations. This causes a secular increase in the average obliquity.

Rubincam, David P.↗

A numerical method for unsteady aerodynamics via acoustics

Formal solutions to the wave equation may be conveniently described within the framework of generalized function theory. A generalized function theory is used to yield a formulation and formal solution of a wave equation describing oscillation of a flat plate from which a numerical method may be derived.

Hodge, Steve↗

Unsteady two dimensional airloads acting on oscillating thin airfoils in subsonic ventilated wind tunnels

The numerical calculation of unsteady two dimensional airloads which act upon thin airfoils in subsonic ventilated wind tunnels was studied. Neglecting certain quadrature errors, Bland's collocation method is rigorously proved to converge to the mathematically exact solution of Bland's integral equation, and a three way equivalence was established between collocation, Galerkin's method and least squares whenever the collocation points are chosen to be the nodes of the quadrature rule used for Galerkin's method. A computer program displayed convergence with respect to the number of pressure basis functions employed, and agreement with known special cases was demonstrated. Results are obtained for the combined effects of wind tunnel wall ventilation and wind tunnel depth to airfoil chord ratio, and for acoustic resonance between the airfoil and wind tunnel walls. A boundary condition is proposed for permeable walls through which mass flow rate is proportional to pressure jump.

Fromme, J.↗

Noise in Josepson effect mixers and the RSJ model

Josephson effect mixers have previously been observed to display 'excess' noise both in experiments with point contacts and in numerical simulations using the resistively shunted junction (RSJ) model. This excess noise causes the mixer noise temperature to be a factor of typically 20-100 times the physical temperature of the device. Previously, this excess was ascribed to conversion from unwanted sidebands of the local oscillator and Josephson frequencies and their harmonics. Our numerical modeling of the RSJ equations has led to a new understanding of the excess noise, which is simply due to the intrinsic Josephson oscillations of the device. In addition, we have extended the modeling to include the previously ignored case of finite device capacitance (i.e. RSJ capacitance parameter beta(sub c) does not equal 0, which is more realistic for lithographically defined Josephson such as shunted tunnel junctions or SNS bridges. For some cases, this yields an improvement of a factor of two in noise temperature from the zero capacitance models. We will discuss the device parameters which optimize the mixer performance for frequencies approaching the characteristic frequency of the device, which is given by the Josephson frequency at the I(sub c)R(sub n) voltage (nu = 2eI(sub c)R(sub n)/h). These modeling results predict good conversion efficiency and a noise temperature within a factor of a few of the physical temperature. Experiments are in progress to determine the accuracy of this modeling using a waveguide mixer at 100 GHz with optimized, resistively shunted Nb tunnel junctions. If the modeling results are valid, they are particularly encouraging for mixers in the submillimeter regime, given the possibility of obtaining non-hysteretic Josephson devices with I(sub c)R(sub n) products in excess of a millivolt, using for instance, high-T(sub c) SNS bridges. We discuss the modifications to the classical RSJ model which are necessary in the quantum regime (h nu greater than kT), and conclude the Josephson mixers may attain noise temperatures less than ten times the quantum limit at high frequencies.

Schoelkopf, R.↗

Seasonal Tidewater Glacier Terminus Oscillations Bias Multi-Decadal Projections of Ice Mass Change

Numerical, process-based simulations of tidewater glacier evolution are necessary to project future sea-level change under various climate scenarios. Previous work has shown that nonlinearities in tidewater glacier and ice stream dynamics can lead to biases in simulated ice mass change in the presence of noisy forcings. Ice sheet modeling projections that are used in the IPCC Assessment Report 6 (AR6) utilize atmospheric and oceanic forcings at annual temporal resolution, omitting any higher frequency forcings. Here, we quantify the effect of seasonal (<1 year) tidewater glacier terminus oscillations on decadal-scale (30 years) mass change. We use an idealized geometry to mimic realistic tidewater glacier geometries and investigate the impact of the magnitude of seasonal oscillations, bed slope at the glacier terminus, and basal friction law. We find that omitting seasonal terminus motion results in biased mass change projections, with up to an 39% underestimate or 25% overestimate of mass loss when seasonality is neglected, depending on bed slope and magnitude of oscillations. The bias is most sensitive to the magnitude of the seasonal terminus oscillations and exhibits very little sensitivity to choice of friction law. Although omitting seasonal terminus motion may bias mass change projections of individual glaciers, further work is needed to investigate the impact on ice sheet mass change projections. In order to achieve this, seasonality in atmospheric and oceanic forcings must be adequately represented and observations of seasonal terminus positions and tidewater glacier thickness changes must be acquired to evaluate numerical models. Note: Zip folder contains all videos listed in the supplemental document and can be downloaded for viewing.

D. Felikson↗

Magnetic Damping of g-Jitter Driven Flows: 3-D Calculations

A 3-D numerical model is developed to represent the oscillating natural convection induced in a cylindrical cavity filled with Ga-doped germanium with and without the presence of an external magnetic field. The model is developed based on the penalty-finite element solution of the equations describing the transport of momentum, heat and solutal element as well as the electromagnetic field distribution in the melt pool. Automatic time step control is applied to help speed up the calculations. Numerical simulations are conducted to study the convection and magnetic damping effects as a function of frequency, directions and amplitudes of g-jitter and also the direction and magnitudes of the applied magnetic fields. The results show that the g-jitter driven flow is time dependent and exhibits a complex recirculating convection pattern in three dimensions and that an applied magnetic field can be employed to suppress this deleterious convective flow and both magnitude and orientation of the applied field are important in magnetic damping of the g-jitter induced convective flows.

Shang, D. Y.↗

An investigation of the information propagation and entropy transport aspects of Stirling machine numerical simulation

Aspects of the information propagation modeling behavior of integral machine computer simulation programs are investigated in terms of a transmission line. In particular, the effects of pressure-linking and temporal integration algorithms on the amplitude ratio and phase angle predictions are compared against experimental and closed-form analytic data. It is concluded that the discretized, first order conservation balances may not be adequate for modeling information propagation effects at characteristic numbers less than about 24. An entropy transport equation suitable for generalized use in Stirling machine simulation is developed. The equation is evaluated by including it in a simulation of an incompressible oscillating flow apparatus designed to demonstrate the effect of flow oscillations on the enhancement of thermal diffusion. Numerical false diffusion is found to be a major factor inhibiting validation of the simulation predictions with experimental and closed-form analytic data. A generalized false diffusion correction algorithm is developed which allows the numerical results to match their analytic counterparts. Under these conditions, the simulation yields entropy predictions which satisfy Clausius' inequality.

Goldberg, Louis F.↗

On the 'flip-flop' instability of Bondi-Hoyle accretion flows

A simple physical interpretation is advanced by means of an analysis of the shock cone in the accretion flows past a compact object and with an examination of the accretion-line stability analyses. The stability of the conical shock is examined against small angular deflections with attention given to several simplifying assumptions. A line instability is identified in the Bondi-Hoyle accretion flows that leads to the formation of a large opening-angle shock. When the opening angle becomes large the instability becomes irregular oscillation. The analytical methodology is compared to previous numerical configurations that demonstrate different shock morphologies. The Bondi-Hoyle accretion onto a compact object is concluded to generate a range of nonlinear instabilities in both homogeneous and inhomogeneous cases with a quasiperiodic oscillation in the linear regime.

Livio, Mario↗

Two dimensional aerodynamic interference effects on oscillating airfoils with flaps in ventilated subsonic wind tunnels

The numerical computation of unsteady airloads acting upon thin airfoils with multiple leading and trailing-edge controls in two-dimensional ventilated subsonic wind tunnels is studied. The foundation of the computational method is strengthened with a new and more powerful mathematical existence and convergence theory for solving Cauchy singular integral equations of the first kind, and the method of convergence acceleration by extrapolation to the limit is introduced to analyze airfoils with flaps. New results are presented for steady and unsteady flow, including the effect of acoustic resonance between ventilated wind-tunnel walls and airfoils with oscillating flaps. The computer program TWODI is available for general use and a complete set of instructions is provided.

Fromme, J.↗

Velocity-coupled flow oscillations in a simulated solid-propellant rocket environment

A comprehensive numerical analysis has been carried out to study the unsteady flowfields in a simulated rocket-motor environment. The model is based on the time-dependent compressible Navier-Stokes equations with a two-equation turbulence closure scheme. Various important aspects of the coupling between acoustic oscillations and mean flowfields, including flow reversal, modification of transport properties, etc., are addressed. Results indicate that multi-dimensional effects play important roles in determining local flow structures and wave characteristics. In much of the domain, acoustic velocity nodal points are observed in the near-wall region. The classical one-dimensional theory fails to describe several important mechanisms associated with velocity-induced flow instabilities.

Yang, Vigor↗

Artificial Boundary Conditions for Computation of Oscillating External Flows

In this paper, we propose a new technique for the numerical treatment of external flow problems with oscillatory behavior of the solution in time. Specifically, we consider the case of unbounded compressible viscous plane flow past a finite body (airfoil). Oscillations of the flow in time may be caused by the time-periodic injection of fluid into the boundary layer, which in accordance with experimental data, may essentially increase the performance of the airfoil. To conduct the actual computations, we have to somehow restrict the original unbounded domain, that is, to introduce an artificial (external) boundary and to further consider only a finite computational domain. Consequently, we will need to formulate some artificial boundary conditions (ABC's) at the introduced external boundary. The ABC's we are aiming to obtain must meet a fundamental requirement. One should be able to uniquely complement the solution calculated inside the finite computational domain to its infinite exterior so that the original problem is solved within the desired accuracy. Our construction of such ABC's for oscillating flows is based on an essential assumption: the Navier-Stokes equations can be linearized in the far field against the free-stream back- ground. To actually compute the ABC's, we represent the far-field solution as a Fourier series in time and then apply the Difference Potentials Method (DPM) of V. S. Ryaben'kii. This paper contains a general theoretical description of the algorithm for setting the DPM-based ABC's for time-periodic external flows. Based on our experience in implementing analogous ABC's for steady-state problems (a simpler case), we expect that these boundary conditions will become an effective tool for constructing robust numerical methods to calculate oscillatory flows.

Tsynkov, S. V.↗

The condensational instability in the solar transition region and corona

The stability of plasmas at temperatures and densities typical of the solar transition region and corona is investigated using both a linear analysis and nonlinear time-dependent numerical simulations. Growth rates, decay rates, and oscillation frequencies of the perturbations determined from the linear analysis are in excellent agreement with the simulations. The nonlinear regime is characterized by a bifurcation of the plasma into a cool dense condensation surrounded by a hot tenuous corona. The condensation may then be accelerated by forces in the plasma such as those arising from gravity or differential heating. Finally, the results of the detailed simulation show that the transition region is a dynamically stable structure which is the result of the nonlinear evolution of the condensational instability.

Oran, E. S.↗

A least squares finite element scheme for transonic flow around harmonically oscillating airfoils

The present investigation shows that a finite element scheme with a weighted least squares variational principle is applicable to the problem of transonic flow around a harmonically oscillating airfoil. For the flat plate case, numerical results compare favorably with the exact solution. The obtained numerical results for the transonic problem, for which an exact solution is not known, have the characteristics of known experimental results. It is demonstrated that the performance of the employed numerical method is independent of equation type (elliptic or hyperbolic) and frequency. The weighted least squares principle allows the appropriate modeling of singularities, which such a modeling of singularities is not possible with normal least squares.

Cox, C. L.↗