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 343 records · Page 19

Effects of material non-linearity on the residual stresses in a dendritic silicon crystal ribbon

Thermal stresses developed in a dendritic silicon crystal ribbon have been shown to cause plastic deformation and residual stresses in the ribbon. This paper presents an implementation of a numerical model proposed for thermoelastoplastic behavior of a material. The model has been used to study the effects of plasticity of silicon on the residual stresses. The material properties required to implement this model are all assumed, and the response of the material to the variations in these assumed parameters of the constitutive law and in the finite element mesh is investigated. The steady state growth process is observed to be periodic with nonzero residual stresses. Numerical difficulties are also encountered in the computer solution process, resulting in sharp jumps and large oscillations in the stress responses.

Ray, Sujit K.↗

Numerical simulation of shock-induced combustion/detonation

A numerical study is conducted to simulate the shock-induced combustion in premixed H2-air mixtures. Two types of bodies, blunt (spherical projectile) and sharp (wedge), are considered in the study. A nine-species, 18-step finite-rate H2-air chemical reaction mechanism coupled with the Navier-Stokes equations is solved. The flow field over the blunt body is found to be unsteady, when the projectile velocity is same as the Chapman-Jouget velocity of the mixture. The unsteadiness is caused by the periodic instabilities originating in the stagnation zone. Numerical results show good qualitative agreement with the ballistic range shadowgraph. In addition, the frequency of oscillations, determined by using the Fourier power spectrum, is found to be in good agreement with the experiment. The flow field over the wedge is found to be stable for the conditions considered in this study. The oblique detonation wave structure is investigated and the important flow features are discussed.

Singh, D. J.↗

Pressure dither in venting analyses

Venting analyses simulate the time-dependent flow of air or another gas from one compartment to another, or to an ambient profile. The flow through the vent connecting the compartments is assumed to be one-dimensional, quasi-steady, and isentropic with an empirical discharge coefficient included to account for irreversible flow losses. When a venting analysis is performed numerically, an upper limit on the accuracy exists due to the effects of the numerical integration method. Because this manifests itself in a first-order method by a rapid oscillation about the true answer, this phenomenon is referred to as pressure dither. This paper explores the origin of the phenomenon and gives simple formulae for predicting when it will occur and the greatest achievable accuracy for a given time step.

Fay, John F.↗

A flux splitting scheme with high-resolution and robustness for discontinuities

A flux splitting scheme is proposed for the general nonequilibrium flow equations with an aim at removing numerical dissipation of Van-Leer-type flux-vector splittings on a contact discontinuity. The scheme obtained is also recognized as an improved Advection Upwind Splitting Method (AUSM) where a slight numerical overshoot immediately behind the shock is eliminated. The proposed scheme has favorable properties: high-resolution for contact discontinuities; conservation of enthalpy for steady flows; numerical efficiency; applicability to chemically reacting flows. In fact, for a single contact discontinuity, even if it is moving, this scheme gives the numerical flux of the exact solution of the Riemann problem. Various numerical experiments including that of a thermo-chemical nonequilibrium flow were performed, which indicate no oscillation and robustness of the scheme for shock/expansion waves. A cure for carbuncle phenomenon is discussed as well.

Wada, Yasuhiro↗

Nonperiodic variations in astrophysical systems: Investigating frequency evolution

We present a method related to the wavelet transform, the Gabor transform, for investigating astronomical time series containing nonconsistent frequencies. Instances in which such data sets may arise include variable star light curves, numerical studies of the gravitational three-body problem, X-ray binaries, and signals from more exotic objects such as planets around pulsars and mass infall from accretion structures onto compact objects. As an illustration of its power, we apply the technique to a numerical data set of a gravitational three-body interaction and to photometry of the rapidly oscillating peculiar A star HD 60435. In the three-body example, the method provides an insightful shorthand that allows for the determination of episodes where the system behaves as two nearly Keplerian orbits. For HD 60435, the power in the main frequency exhibits unusual evolution over the duration of the observation.

Boyd, Patricia T.↗

Model of Wave Driven Flow Oscillation for Solar Cycle

At low latitudes in the Earth's atmosphere, the observed zonal flow velocities are dominated by the semi-annual and quasi-biennial oscillations with periods of 6 months and 20 to 32 months respectively. These terrestrial oscillations, the SAO and QBO respectively, are driven by wave-mean flow interactions due to upward propagating planetary-scale waves (periods of days) and small-scale gravity waves (periods of hours). We are proposing (see also Mayr et al., GRL, 2001) that such a mechanism may drive long period oscillations (reversing flows) in stellar and planetary interiors, and we apply it to the Sun. The reversing flows would occur below the convective envelope where waves can propagate. We apply a simplified, one dimensional, analytical flow model that incorporates a gravity wave parameterization due to Hines (1997). Based on this analysis, our estimates show that relatively small wave amplitudes less than 10 m/s can produce zonal flow amplitudes of 20 m/s, which should be sufficient to generate the observed variations in the magnetic field. To produce the 22-year period of oscillation, a low buoyancy frequency must be chosen, and this places the proposed flow in a region that is close to (and below) the base of the convective envelope. Enhanced turbulence associated with this low stability should help to generate the dynamo currents. With larger stability at deeper levels in the solar interior, the model can readily produce also oscillations with much longer periods. To provide an understanding of the fluid dynamics involved, we present numerical results from a 2D model for the terrestrial atmosphere that exemplify the non-linear nature of the wave interaction for which a mechanical analog is the escapement mechanism of the clock.

Mayr, Hans G.↗

Parametric analysis of a passive cyclic control device for helicopters

A parametric study of a passive device which provides a cyclic longitudinal control moment for a helicopter rotor was performed. It utilizes a rotor blade tip which is structurally decoupled from the blade inboard section. This rotor configuration is generally called the Free-Tip Rotor. A two dimensional numerical model was used to review the Constant Lift Tip Rotor, a predecessor of the current configuration, and then the same model was applied to the Passive Cyclic Control Device. The Constant Lift Tip was proven to have the ability to suppress the vibratory lift loading on the tip around the azimuth and to eliminate a significant negative lift peak on the advancing tip. The Passive Cyclic Control Device showed a once-per-revolution lift oscillation with a large amplitude, while minimizing the higher harmonic terms of the lift oscillation. This once-per-revolution oscillation results in the cyclic moment to trim the rotor longitudinally. A rotor performance analysis was performed by a three dimensional numerical model. It indicated that the vortices shed from the junction between the tip and the inboard section has a strong influence on the tip, and it may severely limit the tip performance. It was also shown that the Free-Tip allows the inboard section to have a larger twist, which results in a better performance.

Kumagai, H.↗

Theory and Simulation of A Novel Viscosity Measurement Method for High Temperature Semiconductor

The properties of molten semiconductors are good indicators for material structure transformation and hysteresis under temperature variations. Viscosity, as one of the most important properties, is difficult to measure because of high temperature, high pressure, and vapor toxicity of melts. Recently, a novel method was developed by applying a rotating magnetic field to the melt sealed in a suspended quartz ampoule, and measuring the transient torque exerted by rotating melt flow on the ampoule wall. The method was designed to measure viscosity in short time period, which is essential for evaluating temperature hysteresis. This paper compares the theoretical prediction of melt flow and ampoule oscillation with the experimental data. A theoretical model was established and the coupled fluid flow and ampoule torsional vibration equations were solved numerically. The simulation results showed a good agreement with experimental data. The results also showed that both electrical conductivity and viscosity could be calculated by fitting the theoretical results to the experimental data. The transient velocity of the melt caused by the rotating magnetic field was found reach equilibrium in about half a minute, and the viscosity of melt could be calculated from the altitude of oscillation. This would allow the measurement of viscosity in a minute or so, in contrast to the existing oscillation cup method, which requires about an hour for one measurement.

Lin, Bochuan↗

Nonradial oscillations of neutron stars

Linear, adiabatic, Newtonian, nonradial pulsation analyses have been performed for finite-temperature neutron star models with a fluid core, solid crust, and thin surface fluid 'ocean', including the effects of the neutron star crust. The pulsation equations are considered, including the spheroidal and toroidal modes. A local analysis is performed to provide information about the pulsation modes in the short-wavelength limit. Numerical calculations are made on the mode spectrum and systematic properties. Damping mechanisms are investigated, including gravitational radiation damping, neutrino emission damping, electromagnetic radiation from an oscillating stellar magnetic field, nonadiabatic effects, and internal friction and viscosity.

Mcdermott, P. N.↗

Analysis of cavitation bubble dynamics in a liquid

General differential equations governing the dynamics of the cavitation bubbles in a liquid were derived. With the assumption of spherical symmetry the governing equations were simplified. Closed form solutions were obtained for simple cases, and numerical solutions were calculated for complicated ones. The growth and the collapse of the bubble were analyzed, oscillations of the bubbles were studied, and the stability of the cavitation bubbles were investigated. The results show that the cavitation bubbles are unstable, and the oscillation is not sinusoidal.

Fontenot, L. L.↗

Non-linear analysis of vibrations of irregular plates

A numerical perturbation method is used to investigate the forced vibrations of irregular plates. Nonlinear terms associated with the midplane stretching are retained in the analysis. The numerical part of the method involves the use of linear, finite element techniques to determine the free oscillation mode shapes and frequencies and to obtain the linear midplane stress resultants caused by the midplane stretching. Representing the solution as an expansion in terms of these linear mode shapes, these modes and the resultants are used to determine the equations governing the time-dependent coefficients of this expansion. These equations are solved by using the method of multiple scales. Specific solutions are given for the main-resonant vibrations of an elliptical plate in the presence of internal resonances. The results indicate that modes other than the driven mode can be drawn into the steady state response. Though the excitation is composed of a single harmonic, the response may not be periodic. Moreover, the particular types of responses that can occur are highly dependent on the mode being excited and are sensitive to small geometrical changes.

Lobitz, D. W.↗

Dynamics of charged particles in a Paul radio-frequency quadrupole trap

A molecular-dynamics simulation of hundreds of ions confined in a Paul trap has been performed. The simulation includes the trapped particles' micromotion and interparticle Coulomb interactions. A random walk in velocity was implemented to bring the secular motion to a given temperature which was numerically measured. When the coupling Gamma is large the ions from concentric shells which undergo a quadrupole oscillation at the RF frequency, while the ions within a shell form a 2D hexagonal lattice. Ion clouds at 5 mK show no RF heating for q(z) less than about 0.6, whereas rapid heating is seen for qz = 0.8.

Prestage, J. D.↗

Satellite observations of variations in tropical cyclone convection caused by upper-tropospheric troughs

Satellite observations and numerical model results have been used to study the relationship between upper-tropospheric forcing and the oscillation of convection of tropical cyclones Florence (1988) and Irene (1981) during their mature stage over open warm oceans (SST greater than or equal to 26 C). It is suggested that the initiation and maintenance of intense convective outbreaks in tropical cyclones are related to the channeling and strengthening of their outflow by upper-tropospheric troughs. It is possible to enhance the convection in response to the outflow jet-induced import of eddy relative angular momentum and ascending motion associated with the thermally direct circulation. Both Florence and Irene are found to intensify after the onset of these convective episodes. It is also suggested that the cessation in the convection of the two tropical cyclones occurs when the upper-tropospheric troughs move near or over the tropical cyclones, resulting in the weakening of their outflow and the entrainment of dry upper-tropospheric air into their inner core.

Rodgers, Edward B.↗

Experimental investigation of the effect of an oscillating airstream (Katzmayr effect) on the characteristics of airfoils

A series of experiments were conducted related to the action of an airstream oscillating vertically on supporting surfaces. The object of the experiments was to verify the very interesting results of Mr. Katzmayr, Director of the Vienna Aerodynamics Laboratory, and, if possible, to obtain more complete data on the effect of the amplitude and velocity of the oscillations of the airstream. The results obtained by Mr. Katzmayr are briefly summarized. The conduct of the numerous experiments to verify his results are described in detail. Experimental results are given in tabular and graphical form.

TOUSSAINT↗

Jupiter's Great Red Spot and other vortices

A theoretical explanation of Jupiter's Great Red Spot (GRS) as the self-organization of vorticity in turbulence is presented. A number of properties of the GRS and other Jovian vortices that are unambiguous from the data are listed. The simplest possible model that explains these properties one at a time rather than in a difficult all-encompassing planetary global circulation model is presented. It is shown that Jovian vortices reflect the behavior of quasi-geostrophic (QG) vortices embedded in an east-west wind with bands of uniform potential vorticity. It is argued that most of the properties of the Jovian vortices can be easily explained and understood with QG theory. Many of the signatures of QG vortices are apparent on Voyager images. In numerical and laboratory experiments, QG vortices relax to approximately steady states like the Jovian vortices, rather than oscillating or rotating Kida ellipses.

Marcus, Philip S.↗

Direct Numerical Simulation of Automobile Cavity Tones

The Navier Stokes equation is solved computationally by the Dispersion-Relation-Preserving (DRP) scheme for the flow and acoustic fields associated with a laminar boundary layer flow over an automobile door cavity. In this work, the flow Reynolds number is restricted to R(sub delta*) < 3400; the range of Reynolds number for which laminar flow may be maintained. This investigation focuses on two aspects of the problem, namely, the effect of boundary layer thickness on the cavity tone frequency and intensity and the effect of the size of the computation domain on the accuracy of the numerical simulation. It is found that the tone frequency decreases with an increase in boundary layer thickness. When the boundary layer is thicker than a certain critical value, depending on the flow speed, no tone is emitted by the cavity. Computationally, solutions of aeroacoustics problems are known to be sensitive to the size of the computation domain. Numerical experiments indicate that the use of a small domain could result in normal mode type acoustic oscillations in the entire computation domain leading to an increase in tone frequency and intensity. When the computation domain is expanded so that the boundaries are at least one wavelength away from the noise source, the computed tone frequency and intensity are found to be computation domain size independent.

Kurbatskii, Konstantin↗

Non-Migrating Diurnal Tides Generated with Planetary Waves in the Mesosphere

We report here the results from a modeling study with our Numerical Spectral Model (NSM) that extends from the ground into thermosphere. The NSM incorporates Hines Doppler Spread Parameterization for small-scale gravity waves (GWs) and describes the major dynamical features of the atmosphere, including the wave driven equatorial oscillations (QBO and SAO), and the seasonal variations of tides and planetary waves. Accounting solely for the solar migrating tidal excitation sources, the NSM generates through dynamical interactions also non-migrating tides in the mesosphere that have amplitudes comparable to those observed. The model produces the diurnal (and semidiurnal) oscillations of the zonal mean (m = 0), and eastward and westward propagating tides for zonal wave numbers m = 1 to 4. To identify the mechanism of excitation for these tides, a numerical experiment is performed. The NSM is run without the heat source for the zonal-mean circulation and temperature variation, and the amplitudes of the resulting nonmigrating tides are then negligibly small. This leads to the conclusion that the planetary waves, which normally are excited in the NSM by instabilities but are suppressed in this case, generate the nonmigrating tides through nonlinear interactions with the migrating tides.

Mayr, H. G.↗

Entropy Generation/Availability Energy Loss Analysis Inside MIT Gas Spring and "Two Space" Test Rigs

The results of the entropy generation and availability energy loss analysis under conditions of oscillating pressure and oscillating helium gas flow in two Massachusetts Institute of Technology (MIT) test rigs piston-cylinder and piston-cylinder-heat exchanger are presented. Two solution domains, the gas spring (single-space) in the piston-cylinder test rig and the gas spring + heat exchanger (two-space) in the piston-cylinder-heat exchanger test rig are of interest. Sage and CFD-ACE+ commercial numerical codes are used to obtain 1-D and 2-D computer models, respectively, of each of the two solution domains and to simulate the oscillating gas flow and heat transfer effects in these domains. Second law analysis is used to characterize the entropy generation and availability energy losses inside the two solution domains. Internal and external entropy generation and availability energy loss results predicted by Sage and CFD-ACE+ are compared. Thermodynamic loss analysis of simple systems such as the MIT test rigs are often useful to understand some important features of complex pattern forming processes in more complex systems like the Stirling engine. This study is aimed at improving numerical codes for the prediction of thermodynamic losses via the development of a loss post-processor. The incorporation of loss post-processors in Stirling engine numerical codes will facilitate Stirling engine performance optimization. Loss analysis using entropy-generation rates due to heat and fluid flow is a relatively new technique for assessing component performance. It offers a deep insight into the flow phenomena, allows a more exact calculation of losses than is possible with traditional means involving the application of loss correlations and provides an effective tool for improving component and overall system performance.

Ebiana, Asuquo B.↗