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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Modeling Wave Driven Non-linear Flow Oscillations: The Terrestrial QBO and a Solar Analog

The Quasi Biennial Oscillation (QBO) of the zonal circulation observed in the terrestrial atmosphere at low latitudes is driven by wave mean flow interaction as was demonstrated first by Lindzen and Holton (1968), shown in a laboratory experiment by Plumb and McEwan (1978), and modeled by others (e.g., Plumb, Dunkerton). Although influenced by the seasonal cycle of solar forcing, the QBO, in principle, represents a nonlinear flow oscillation that can be maintained by a steady source of upward propagating waves. The wave driven non-linearity is of third or odd order in the flow velocity, which regenerates the fundamental harmonic itself to keep the oscillation going - the fluid dynamical analog of the displacement mechanism in the mechanical clock. Applying Hines' Doppler Spread Parameterization (DSP) for gravity waves (GW), we discuss with a global-scale spectral model numerical experiments that elucidate some properties of the QBO and its possible effects on the climatology of the atmosphere. Depending on the period of the QBO, wave filtering can cause interaction with the seasonal variations to produce pronounced oscillations with beat periods around 10 years. Since the seasonal cycle and its variability influence the period of the QBO, it may also be a potent conduit of solar activity variations to lower altitudes. Analogous to the terrestrial QBO, we propose that a flow oscillation may account for the 22-year periodicity of the solar magnetic cycle, potentially answering Dicke (1978) who asked, "Is there a chronometer hidden deep inside the Sun?" The oscillation would occur below the convection region, where gravity waves can propagate. Employing a simplified, analytic model, Hines' DSP is applied to estimate the flow oscillation. Depending on the adopted horizontal wavelengths of GW's, wave amplitudes less than 10 m/s can be made to produce oscillating zonal flows of about 20 m/s that should be large enough to generate a significant oscillation in the magnetic field. For the large length scales of the Sun, the flow cycle period tends to be very long. The period, however, can be made to be 22 years, provided the buoyancy frequency (stability) is sufficiently small, thus placing the proposed flow near the base of the convection zone where a dynamo is now believed to operate.

Mayr, Hans G.↗
NASA NTRS1990

Surface tension and viscosity from damped free oscillations of viscous droplets

Damped oscillations of a viscous droplet in vacuum or in an inert gas of negligible density are considered. The dependence of the complex decay factor on the properties of the liquid is investigated for the first time, and numerical results are compared with earlier studies for special cases. A new method is developed to determine both surface tension and viscosity from a single experiment in which the damping rate and frequency of oscillations are measured. The procedure to determine surface tension and viscosity from oscillating levitated liquids is outlined, and results presented for various modes of shape oscillations.

Suryanarayana, P. V. R.↗
NASA NTRS2002

Analysis of High Order Difference Methods for Multiscale Complex Compressible Flows

Accurate numerical simulations of complex multiscale compressible viscous flows, especially high speed turbulence combustion and acoustics, demand high order schemes with adaptive numerical dissipation controls. Standard high resolution shock-capturing methods are too dissipative to capture the small scales and/or long-time wave propagations without extreme grid refinements and small time steps. An integrated approach for the control of numerical dissipation in high order schemes with incremental studies was initiated. Here we further refine the analysis on, and improve the understanding of the adaptive numerical dissipation control strategy. Basically, the development of these schemes focuses on high order nondissipative schemes and takes advantage of the progress that has been made for the last 30 years in numerical methods for conservation laws, such as techniques for imposing boundary conditions, techniques for stability at shock waves, and techniques for stable and accurate long-time integration. We concentrate on high order centered spatial discretizations and a fourth-order Runge-Kutta temporal discretizations as the base scheme. Near the bound-aries, the base scheme has stable boundary difference operators. To further enhance stability, the split form of the inviscid flux derivatives is frequently used for smooth flow problems. To enhance nonlinear stability, linear high order numerical dissipations are employed away from discontinuities, and nonlinear filters are employed after each time step in order to suppress spurious oscillations near discontinuities to minimize the smearing of turbulent fluctuations. Although these schemes are built from many components, each of which is well-known, it is not entirely obvious how the different components be best connected. For example, the nonlinear filter could instead have been built into the spatial discretization, so that it would have been activated at each stage in the Runge-Kutta time stepping. We could think of a mechanism that activates the split form of the equations only at some parts of the domain. Another issue is how to define good sensors for determining in which parts of the computational domain a certain feature should be filtered by the appropriate numerical dissipation. For the present study we employ a wavelet technique introduced in as sensors. Here, the method is briefly described with selected numerical experiments.

Sjoegreen, Bjoern↗
NASA NTRSDate not supplied

Optical Time Transfer for Bistatic SAR Spacecraft

A spacecraft-to-spacecraft optical time-transfer simulation has been developed as a tool for informing NASA’s Surface Deformation and Change (SDC) mission architecture. The SDC mission will combine radar images from multiple spacecraft to improve understanding of the Earth’s sea-level and landscape changes. Spacecraft must be precisely synchronized in order to create sharp radar images. Simulation of multiple spacecraft time-synchronizing via laser communication can inform technology choices of a mission by providing a picosecond-precision level estimate of clock error. This timing and ranging simulation has been combined with a radar system performance analysis pipeline. The simulated timing errors are used to predict performance of bistatic SAR systems in the presence of oscillator noise and time synchronization in accuracy. This analysis includes both analytic approximation equations from existing literature, and a numerical radar simulation to extract key system performance parameters like phase error and signal-to-noise ratio (SNR)degradation. Precision time-transfer techniques facilitate the accurate synchronization of clocks between any combination of terminals. Most time-transfer technology for comparing two clocks at different terminals use radio frequencies (RF) to measure the time delay between the sending and receiving of signals. Laser technology offers the capability to transmit high data rates with systems that are of smaller size and lower power than comparable RF systems. The clocks on independent spacecraft will have some phase and frequency errors between them that result in clock drift. The two clock models that are included in this bi-directional MATLAB simulation are a cesium-based Chip-Scale Atomic Clock (CSAC) and a rubidium-based Miniature Atomic Clock (MAC). The CSAC has flown as hardware for small satellite missions such as the University of Florida’s CHOMPTT mission. A study of example orbits, including that of NASA NASA-ISRO Synthetic Aperture Radar Mission (NISAR) mission, and lasing rates demonstrate the impact of flight configuration parameters on the synchronization error between two spacecraft. The MATLAB timing simulation uses a Runge-Kutta 4th-order method to propagate spacecraft orbits and computes the light-travel time estimate between them. The simulation outputs the estimated range and estimated clock error based on a user-defined spacecraft cluster configuration. The radar simulation and analytic approximations are applied to evaluate a potential future NASA bistatic SAR constellation architecture. In the proposed architecture, satellites follow each other in the same orbit at 800 km altitude, with a 210 km baseline. We also baseline the CSAC as an ultra stable oscillator, and use NASA’s NISAR for baseline radar system parameters to compute a clock-system introduced phase error of 5.6 degrees without synchronization by frequent time transfer. We build on this base case with a sensitivity analysis of radar performance over a proposed range of constellation and radar system parameters. With this analysis pipeline, we comment on which radar parameters should or should not be changed to minimize synchronization requirements. This analysis technique could be extended or modified to evaluate the timing requirements of other geometries for other future multistatic SAR missions, or other interferometric satellite missions.

Surface Deformation and Change↗
NASA NTRS1982

Band-pass filtering of one year of daily mean pressures on Mars

The sol-mean atmospheric pressures measured at the two Viking landers exhibit fluctuations ranging in period from seasonal to a few days. The time series are highly nonstationary, which severely limits standard methods such as power-spectrum analysis. Advantage is taken of the rather precise annual periodicity of these data to design sharp numerical filters and use them to split the series at each lander into six time series, each containing oscillations associated with a discrete frequency band. This analysis reveals 2- to 4-sol waves during the cold seasons at both landers, and 8-sol waves during the spring at Lander 2, confirming previous results from spectral analysis. As previously predicted, high-frequency waves appear suppressed during the second dust storm. Certain low-frequency signals, probably associated with the global dust storms, are detected and demonstrated to be significant for the first time.

Niver, D. S.↗
NASA NTRS1985

Three-dimensional unsteady lifting surface theory in the subsonic range

The methods of the unsteady lifting surface theory are surveyed. Linearized Euler's equations are simplified by means of a Galileo-Lorentz transformation and a Laplace transformation so that the time and the compressibility of the fluid are limited to two constants. The solutions to this simplified problem are represented as integrals with a differential nucleus; these results in tolerance conditions, for which any exact solution must suffice. It is shown that none of the existing three-dimensional lifting surface theories in subsonic range satisfy these conditions. An oscillating elliptic lifting surface which satisfies the tolerance conditions is calculated through the use of Lame's functions. Numerical examples are calculated for the borderline cases of infinitely stretched elliptic lifting surfaces and of circular lifting surfaces. Out of the harmonic solutions any such temporal changes of the down current are calculated through the use of an inverse Laplace transformation.

Kuessner, H. G.↗
NASA NTRS2003

A Cartesian Adaptive Level Set Method for Two-Phase Flows

In the present contribution we develop a level set method based on local anisotropic Cartesian adaptation as described in Ham et al. (2002). Such an approach should allow for the smallest possible Cartesian grid capable of resolving a given flow. The remainder of the paper is organized as follows. In section 2 the level set formulation for free surface calculations is presented and its strengths and weaknesses relative to the other free surface methods reviewed. In section 3 the collocated numerical method is described. In section 4 the method is validated by solving the 2D and 3D drop oscilation problem. In section 5 we present some results from more complex cases including the 3D drop breakup in an impulsively accelerated free stream, and the 3D immiscible Rayleigh-Taylor instability. Conclusions are given in section 6.

Ham, F.↗
NASA NTRSDate not supplied

A Multivibrational Mode Relaxation Scheme for Two-Temperature Thermochemical State Models

An approach to estimate the effective vibrational relaxation time for polyatomic molecules under the context of two-temperature thermochecmial state model is proposed. Landau-Teller vibrational relaxation equation is used as a basis, and vibrational energies are estimated with partition energy functions assuming harmonic oscillator model. A special weighting procedure is presented, and effective vibrational relaxation time is derived. A series of numerical experiments with CO2 and NH3 mixtures is performed, and the aeroheating results are presented. The results of CO2 environment simulations are compared with the vibrational relaxation model of Park et al. [JTHT, 8(1994), pp. 9–23]. The simulations are also conducted with a modified Landau-Teller expression.

Vibrational relaxation time↗
NASA NTRS2022

Seasonal and Interannual Variability of Tidal Mixing Signatures in Indonesian Seas from High-Resolution Sea Surface Temperature

With their complex narrow passages and vigorous mixing, the Indonesian seas provide the only low-latitude pathway between the Pacific and Indian Oceans and thus play an essential role in regulating Pacific-Indian Ocean exchange, regional air-sea interaction, and ultimately, global climate phenomena. While previous investigations using remote sensing and numerical simulations strongly suggest that this mixing is tidally driven, the impacts of monsoon and El Niño Southern Oscillation (ENSO) on tidal mixing in the Indonesian seas must play an important role. Here we use high-resolution sea surface temperature from June 2002 to June 2021 to reveal monsoon and ENSO modulations of mixing. The largest spring-neap (fortnightly) signals are found to be localized in the narrow passages/straits and sills, with more vigorous tidal mixing during the southeast (boreal summer) monsoon and El Niño than that during the northwest (boreal winter monsoon) and La Niña. Therefore, tidal mixing, which necessarily responds to seasonal and interannual changes in stratification, must also play a feedback role in regulating seasonal and interannual variability of water mass transformations and Indonesian throughflow. The findings have implications for longer-term variations and changes of Pacific–Indian ocean water mass transformation, circulation, and climate.

Indonesian seas↗
NASA NTRS1997

Zonal and Meridional Mass Redistribution in the Tropical Pacific During El Nino Southern Oscillation Events

This note aims at documenting the ENSO-related sea level changes for the whole tropical Pacific, and for the 1979-1996 period which covers numerous El Nino and La Nina events. For this purpose, we rely on (a) 0/450 dbar surface dynamic height anomaly (an alias for sea level) derived from XBT and TOGA-TAO measurements during 1979-1996, (b) TOPEX/Poseidon derived sea level anomaly during 1993-1996, and (c) Murtugudde et al. [1996] model sea level output during 1979-1995.

Delcroix, Thierry↗
NASA NTRS1975

Vaporization response of evaporating drops with finite thermal conductivity

A numerical computing procedure was developed for calculating vaporization histories of evaporating drops in a combustor in which travelling transverse oscillations occurred. The liquid drop was assumed to have a finite thermal conductivity. The system of equations was solved by using a finite difference method programmed for solution on a high speed digital computer. Oscillations in the ratio of vaporization of an array of repetitivity injected drops in the combustor were obtained from summation of individual drop histories. A nonlinear in-phase frequency response factor for the entire vaporization process to oscillations in pressure was evaluated. A nonlinear out-of-phase response factor, in-phase and out-of-phase harmonic response factors, and a Princeton type 'n' and 'tau' were determined. The resulting data was correlated and is presented in graphical format. Qualitative agreement with the open literature is obtained in the behavior of the in-phase response factor. Quantitatively the results of the present finite conductivity spray analysis do not correlate with the results of a single drop model.

Agosta, V. D.↗
NASA NTRS1986

Transitions to chaos in two-dimensional double-diffusive convection

The partial differential equations governing two-dimensional thermosolutal convection in a Boussinesq fluid with free boundary conditions have been solved numerically in a regime where oscillatory solutions can be found. A systematic study of the transition from nonlinear periodic oscillations to temporal chaos has revealed sequences of period-doubling bifurcations. Overstability occurs if the ratio of the solutal to the thermal diffusivity tau is less than 1 and the solutal Rayleigh number Rs is sufficiently large. Solutions have been obtained for two representative values of tau. For tau = 0.316, R(s) = 10,000, symmetrical oscillations undergo a bifurcation to asymmetry, followed by a cascade of period-doubling bifurcations leading to aperiodicity, as the thermal Rayleigh number R(T) is increased. At higher values of R(T), the bifurcation sequence is repeated in reverse, restoring simple periodic solutions. As R(T) is further increased more period-doubling cascades, followed by chaos, can be identified. Within the chaotic regions there are narrow periodic windows, and multiple branches of oscillatory solutions coexist. Eventually the oscillatory branch ends and only steady solutions can be found. The development of chaos has been investigated for tau = 0.1 by varying R(T) for several different values of R(s). When R(s) is sufficiently small there are periodic solutions whose period becomes infinite at the end of the oscillatory branch. As R(s) is increased, chaos appears in the neighborhood of these heteroclinic orbits. At higher values of R(s), chaos is found for a broader range in R(T). A truncated fifth-order model suggest that the appearance of chaos is associated with heteroclinic bifurcations.

Knobloch, E.↗
NASA NTRS1986

Experimental forecasts of El Nino

A deterministic numerical model of the coupled evolution of the tropical ocean and atmosphere was used to forecast all El Nino/Southern Oscillation (ENSO) events from 1970 to 1986. More particularly, the model, originally developed for studying large-scale ocean-atmosphere interactions in the tropics, successfully predicted the characteristics of the spatial and temporal structure of ENSO observed in the study interval. The model indicated that rainfall moving eastward over the Pacific slackens the surface winds that would otherwise cool the eastern Pacific by drawing up cooler subsurface waters. The oceanic thermocline increases, a poleward flow of westerly flowing warn waters deplets the equatorial warm water reservoir, and sea surface temperatures decline. These ENSO conditions are statistically tractable with the model several months in advance, provided upper ocean layer thermal data are available.

Cane, M. A.↗
NASA NTRS1988

Stability analysis of large space structure control systems with delayed input

Large space structural systems, due to their inherent flexibility and low mass to area ratio, are represented by large dimensional mathematical models. For implementation of the control laws for such systems a finite amount of time is required to evaluate the control signals; and this time delay may cause instability in the closed loop control system that was previously designed without taking the input delay into consideration. The stability analysis of a simple harmonic oscillator representing the equation of a single mode as a function of delay time is treated analytically and verified numerically. The effect of inherent damping on the delay is also analyzed. The control problem with delayed input is also formulated in the discrete time domain.

Reddy, A. S. S. R.↗
NASA NTRS1988

Stability analysis of large space structure control systems with delayed input

Large space structural systems, due to their inherent flexibility and low mass to area ratio, are represented by large dimensional mathematical models. For implementation of the control laws for such systems a finite amount of time is required to evaluate the control signals; and this time delay may cause instability in the closed loop control system that was previously designed without taking the input delay into consideration. The stability analysis of a simple harmonic oscillator representing the equation of a single mode as a function of delay time is analyzed analytically and verified numerically. The effect of inherent damping on the delay is also analyzed. The control problem with delayed input is also formulated in the discrete time domain.

Reddy, A. S. S. R.↗
NASA NTRS1989

Nonlinear filters for efficient shock computation

A new type of methods for the numerical approximation of hyperbolic conservation laws with discontinuous solution is introduced. The methods are based on standard finite difference schemes. The difference solution is processed with a nonlinear conservation form filter at every time level to eliminate spurious oscillations near shocks. It is proved that the filter can control the total variation of the solution and also produce sharp discrete shocks. The method is simpler and faster than many other high resolution schemes for shock calculations. Numerical examples in one and two space dimensions are presented.

Engquist, Bjorn↗
NASA NTRS1993

Modeling of Non-Spherical Droplet Dynamics

A two-dimensional time-dependent computer code based on the modified Arbitrary Lagrangian Eulerian (ALE) technique, has been developed to simulate non-spherical droplet dynamics and evaporation under convective flows at real rocket combustion chamber conditions. The equations of mass, momentum, energy and species are simultaneously solved for both liquid and gas phases with an accurate dynamic interface tracking. The jump boundary conditions across the deforming droplet surface are obtained by applying the integral forms of conservation of mass, momentum, and energy. At each time step, the interface geometry and flow properties at the droplet surface are implicitly solved by satisfying the interface boundary conditions. A Lagrangian technique was developed to track the arbitrarily moving interface between the liquid droplet and the external gas. An elliptic grid generator is adopted to dynamically reconstruct grids both inside and outside the droplet surface. This code has been used to study droplet oscillation, droplet deformation/breakup, nonspherical droplet evaporation in both low and high pressure convective flows. This presentation briefly describes the numerical algorithm for modeling of the nonspherical droplet dynamics and demonstrates the representative simulation results of nonspherical droplet evaporation at low and high pressure convective flows. Potential applications of this code to rocket combustor design and performance predictions are discussed.

Deng, Zheng-Tao↗
NASA NTRS2003

Inertio Gravity Waves in the Upper Mesosphere

In the polar region of the upper mesosphere, horizontal wind oscillations have been observed with periods around 10 hours (Hernandez et al., 1992). Such waves are generated in our Numerical Spectral Model (NSM) and appear to be inertio gravity waves (IGW). Like the planetary waves (PW) in the model, the IGWs are generated by instabilities that arise in the mean zonal circulation. In addition to stationary waves for m = 0, eastward and westward propagating waves for m = 1 to 4 appear above 70 km that grow in magnitude up to about 110 km, having periods between 9 and 11 hours. The m = 1 westward propagating IGWs have the largest amplitudes, which can reach at the poles 30 m/s. Like PWs, the IGWs are intermittent but reveal systematic seasonal variations, with the largest amplitudes occurring generally in winter and spring. The IGWs propagate upward with a vertical wavelength of about 20 km.

Mayr, H. G.↗