Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dissipative dynamics”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Developments in the computation of turbulent boundary layers

Computational techniques applicable to turbulent boundary layers are classified into solutions of Reynolds-averaged equations, in which all the effects of the turbulence are modelled, and solutions of three-dimensional, time dependent Navier-Stokes equations, in which the large eddies are calculated and only the turbulence at scales smaller than the computational mesh spacings has to be modelled. Current computation costs place engineering computations in the first of these categories; large eddy simulations are appropriate currently for special studies of the dynamical processes of turbulence in idealized flow fields. It is shown that the two methods are interrelated and that each can gain from advances in the other. The degree of success of a pair of increasingly complex Reynolds stress models to broaden their range of applicability is examined through comparisons with experimental data for a variety of flow conditions. An example of a large-eddy simulation is presented, compared with experimental results, and used to evaluate the models for pressure rate-of-strain correlations and dissipation in the Reynolds-averaged equations.

Rubesin, M. W.↗
NASA NTRS2024

Inertia of Ionospheric Conductance During Electron Precipitation Events

Using the SuperThermal Electron Transport (STET) model coupled with the Super‐thermal Proton Electron Atomic Hydrogen—tRansport in the Ionosphere and Thermosphere (SPEAH‐RIT) codes, we demonstrate that temporal variability of ionospheric conductance is defined by several time scales: the temporal variations of the magnetospheric source, the start time of electron precipitation, and the termination of the corresponding source. In these cases, the time scales are defined by dissipation of energetic electrons and effective recombination processes. The results presented in this paper were applied in the regions of pulsating aurora and polararcs, demonstrating the fact that ionospheric conductance requires some time to form and decay. These time delays constitute an effective “inertia” in the conductance calculation which is- not accounted for in many global models which assume an instantaneous connection between precipitation and conductance. Ionospheric conductance inertia influences the temporal variation in ionospheric and magnetospheric electric fields, and as a result, impacts magnetosphere‐ionosphere (MI) dynamics on a global scale.

George V Khazanov↗
NASA NTRS1994

The tropical water and energy cycles in a cumulus ensemble model. Part 1: Equilibrium climate

A cumulus ensemble model is used to study the tropical water and energy cycles and their role in the climate system. The model includes cloud dynamics, radiative processes, and microphysics that incorporate all important production and conversion processes among water vapor and five species of hydrometeors. Radiative transfer in clouds is parameterized based on cloud contents and size distributions of each bulk hydrometeor. Several model integrations have been carried out under a variety of imposed boundary and large-scale conditions. In Part 1 of this paper, the primary focus is on the water and heat budgets of the control experiment, which is designed to simulate the convective - radiative equilibrium response of the model to an imposed vertical velocity and a fixed sea surface temperature at 28 C. The simulated atmosphere is conditionally unstable below the freezing level and close to neutral above the freezing level. The equilibrium water budget shows that the total moisture source, M(sub s), which is contributed by surface evaporation (0.24 M(sub s)) and the large-scale advection (0.76 M(sub s)), all converts to mean surface precipitation bar-P(sub s). Most of M(sub s) is transported verticaly in convective regions where much of the condensate is generated and falls to surface (0.68 bar-P(sub s)). The remaining condensate detrains at a rate of 0.48 bar-P(sub s) and constitutes 65% of the source for stratiform clouds above the melting level. The upper-level stratiform cloud dissipates into clear environment at a rate of 0.14 bar-P(sub s), which is a significant moisture source comparable to the detrained water vapor (0.15 bar-P(sub s)) to the upper troposphere from convective clouds. In the lower troposphere, stratiform clouds evaporate at a rate of 0.41 bar-P(sub s), which is a more dominant moisture source than surface evaporation (0.22 bar-P(sub s)). The precipitation falling to the surface in the stratiform region is about 0.32 bar-P(sub s). The associated latent heating in the water cycle is the dominant source in the heat budget that generates a net upward motion in convective regions, upper stratiform regions (above the freezing level), and a downward motion in the lower stratiform regions. The budgets reveal a cycle of water and energy resulted from radiation-dynamic-convection interactions that maintain equilibrium of the atmosphere.

Sui, C. H.↗
NASA NTRS1993

A nonlinear control design for energy sink simulation in the Euler-Poinsot problem

A nonlinear control design is presented for the purpose of quantitatively simulating the effects of internal damping mechanisms modeled as energy sinks on the attitude dynamics of rigid body spacecraft. Damping is important because it is often the driving mechanism behind passive attitude acquisition maneuvers. Introduction of the controller into the Euler attitude equations of motion allows for the explicit representation of damping without the introduction of additional degrees of freedom required for a physical damping mechanism. This result is significant because perturbation techniques which rely on the closed form solution of the unperturbed problem can then be used to analyze the effects of perturbations upon a damped system. The controller is designed to dissipate kinetic energy while maintaining the magnitude of the angular momentum vector. Control torques are nonlinear functions of the angular momentum components expressed in a body-fixed frame. A numerical simulation of an actual damping mechanism during a decay from minor axis spin into a flat spin is presented showing that the nonlinear controller gives a good qualitative representation and, in many instances, a good quantitative approximation of the attitude motion of a representative spacecraft containing a damping mechanism.

Kammer, Daniel C.↗
NASA NTRS2010

Time Frequency Analysis of Spacecraft Propellant Tank Spinning Slosh

Many spacecraft are designed to spin about an axis along the flight path as a means of stabilizing the attitude of the spacecraft via gyroscopic stiffness. Because of the assembly requirements of the spacecraft and the launch vehicle, these spacecraft often spin about an axis corresponding to a minor moment of inertia. In such a case, any perturbation of the spin axis will cause sloshing motions in the liquid propellant tanks that will eventually dissipate enough kinetic energy to cause the spin axis nutation (wobble) to grow further. This spinning slosh and resultant nutation growth is a primary design problem of spinning spacecraft and one that is not easily solved by analysis or simulation only. Testing remains the surest way to address spacecraft nutation growth. This paper describes a test method and data analysis technique that reveal the resonant frequency and damping behavior of liquid motions in a spinning tank. Slosh resonant frequency and damping characteristics are necessary inputs to any accurate numerical dynamic simulation of the spacecraft.

Green, Steven T.↗
NASA NTRS1991

Supercritical baroclinic disturbances under the influence of topography

An Eady model modified to include Ekman dissipation and sloping horizontal boundaries is used to investigate the nonlinear evolution of baroclinic waves in the presence of surface topography. The effectiveness of the topography to phase lock and equilibrate a given mode vs that of baroclinity to propagate and vacillate that mode is found to depend on the topographic height, its zonal structure, and the level of supercriticality. When topography is sinusoidal and of the same wavelength as the baroclinically most unstable mode, it induces a short-period amplitude modulation whose envelope represents the baroclinic evolution. The dynamics of this modulation is explored analytically. When the sinusoidal topography is of a longer zonal wavelength, a 'mixed' wave state persists with amplitude dominance shifting between the topographic mode and the baroclinic mode, depending on the relative strength of form drag and baroclinicity. Model results are also compared, in appropriate parameter limits, with near-resonance asymptotic results and annulus experiments.

Chou, Shih-Hung↗
NASA NTRS1998

Numerical Modeling Studies of Wake Vortex Transport and Evolution Within the Planetary Boundary Layer

In support of the wake vortex effect of the Terminal Area Productivity program, we have put forward four tasks to be accomplished in our proposal. The first task is validation of two-dimensional wake vortex-turbulence interaction. The second task is investigation of three-dimensional interaction between wake vortices and atmospheric boundary layer (ABL) turbulence. The third task is ABL studies. The, fourth task is addition of a Klemp-Durran condition at the top boundary for TASS model. The accomplishment of these tasks will increase our understanding of the dynamics of wake vortex and improve forecasting systems responsible for air safety and efficiency. The first two tasks include following three parts: (a) Determine significant length scale for vortex decay and transport, especially the length scales associated with the onset of Crow instability (Crow, 1970); (b) Study the effects of atmospheric turbulence on the decay of the wake vortices; and (c) Determine the relationships between decay rate, transport properties and atmospheric parameters based on large eddy simulation (LES) results and the observational data. These parameters may include turbulence kinetic energy, dissipation rate, wind shear and atmospheric stratification. The ABL studies cover LES modeling of turbulence structure within planetary boundary layer under transition and stable stratification conditions. Evidences have shown that the turbulence in the stable boundary layer can be highly intermittent and the length scales of eddies are very small compared to those in convective case. We proposed to develop a nesting grid mesh scheme and a modified Klemp-Durran conditions (Klemp and Wilhelmson, 1978) at the top boundary for TASS model to simulate planetary boundary layer under stable stratification conditions. During the past year, our group has made great efforts to carry out the above mentioned four tasks simultaneously. The work accomplished in the last year will be described in the next section.

Lin, Yuh-Lang↗
NASA NTRS1992

Chromospheric dynamics based on infrared solar brightness variations

The NASA Kuiper Airborne Observatory was used to observe far-infrared continuum brightness fluctuations in the lower chromosphere due to solar 5 minute oscillations on the quiet sun. Brightness measurements made at 50, 100, 200, and 400 microns show a strong correlation with visible-line Doppler measurements from photospheric and chromospheric altitudes. The motion of the chromosphere is nearly in phase over a large range of heights, while the infrared brightness lags the Doppler velocity by phases varying from significantly less than 90 deg at low altitudes to nearly 90 deg at higher altitudes. It is proposed that this is the result of a nonadiabatic response of the chromospheric gas to compression and may indicate an important mechanism for wave dissipation. Thermal relaxation times ranging from about 40 s at 340 km above the tau(5000) = 1 photosphere to about 300 s at 600 km are proposed.

Kopp, G.↗
NASA NTRS2003

Advanced Computational Aeroacoustics Methods for Fan Noise Prediction

Direct computation of fan noise is presently not possible. One of the major difficulties is the geometrical complexity of the problem. In the case of fan noise, the blade geometry is critical to the loading on the blade and hence the intensity of the radiated noise. The precise geometry must be incorporated into the computation. In computational fluid dynamics (CFD), there are two general ways to handle problems with complex geometry. One way is to use unstructured grids. The other is to use body fitted overset grids. In the overset grid method, accurate data transfer is of utmost importance. For acoustic computation, it is not clear that the currently used data transfer methods are sufficiently accurate as not to contaminate the very small amplitude acoustic disturbances. In CFD, low order schemes are, invariably, used in conjunction with unstructured grids. However, low order schemes are known to be numerically dispersive and dissipative. dissipative errors are extremely undesirable for acoustic wave problems. The objective of this project is to develop a high order unstructured grid Dispersion-Relation-Preserving (DRP) scheme. would minimize numerical dispersion and dissipation errors. contains the results of the funded portion of the project. scheme on an unstructured grid has been developed. constructed in the wave number space. The characteristics of the scheme can be improved by the inclusion of additional constraints. Stability of the scheme has been investigated. Stability can be improved by adopting the upwinding strategy.

Envia, Edmane↗
NASA NTRSDate not supplied

A New Approach for a Wider Class of Entropy Split Methods for Compressible Gas Dynamics and MHD

The high order entropy split methods of Sjögreen & Yee [1, 2] by entropy splitting of the compressible Euler (inviscid) flux derivatives for a thermally-perfect gas are based on Harten’s entropy function [3, 4, 5]. Their derivation takes advantage of the homogeneity property of Euler flux, symmetrizable Euler flux derivatives and energy-norm stability in conjunction with high order classical spatial central, DRP (dispersion relation-preserving) [6, 7, 8] or Padé (compact) spatial discretizations [9] with summation-by-parts (SBP) operators [10]. Our entropy split methods have been proven entropy conserving and stable [1, 11, 12]. Our proofs do not rely on a two-point numerical flux, but rather only a linear difference operator is required to derive these methods. To extend the entropy split method for the MHD, we used the Godunov symmetrizable non-conservative MHD form [12, 13, 14]. These high order entropy split methods not only preserve certain physical properties of the chosen governing equations but are also known to either improve numerical stability, and/or minimize aliasing errors in long time integration of turbulent flow computations without the aid of added numerical dissipation. In our previous published work, extensive error norm comparison with grid refinement was performed to show the high accuracy performance of these methods. These studies also showed how well the entropy split methods conserve the entropy, momentum and mass, and preserve the kinetic energy for long time integration of the various flows [1, 2, 12, 13, 14]. The objective of the present work is to use a new approach to obtain a wider class of entropy split methods consisting of a two-point numerical flux portion and a non-conservative portion in such a way that the homogeneity property of the compressible Euler flux is not required. For high order classical spatial central, DRP (dispersion relation-preserving) or Padé (compact) spatial discretizations, this new approach can be proven to be entropy conservative with conservative spatial dsicretizations while at the same time allowing a wider class of symmetrizable inviscid flux derivatives. We also use this generalization to derive an entropy split scheme that is entropy conserving for the equations of MHD without the homogeneity property using the Godunov symmetrizable ideal MHD formulation [15].

High Order Physical Preserving Methods↗
NASA NTRS1992

Kinetic structure of slow shocks - Effects of the electromagnetic ion/ion cyclotron instability

The structure of slow magnetosonic shocks in the low beta regime is analyzed with attention given to ion heating and the effects of waves upstream of the electromagnetic ion/ion cyclotron (EMIIC) instability. Shock formation is assessed by means of three methods - a relaxation method and two based on dynamic flow interactions - to determine the effects of initialization and boundary conditions on the formation. Good solutions are found with the piston method and the similar flow-flow method in which the plasma is injected from two boundaries to form two slow shocks. Plasma parameters and shock normal angle are found to be the key variables dictating the structure of the magnetosonic shocks. Four unique classes of resultant shock structures are described in which classical, steady, or nonsteady behavior is found. The analysis also yields insight into the relationship between EMIIC instability and ion dissipation.

Omidi, N.↗
NASA NTRS1990

Magnetohydrodynamic modeling of the solar corona

The ideal and resistive magnetohydrodynamic (MHD) model is used to examine the dynamics and structure of the solar corona. When the coronal magnetic field is deformed by photospheric flow it can evolve to states that become unstable to ideal MHD modes. The nonlinear evolution of these instabilities can lead to the generation of current sheets, field line reconnection, and energy release. The disruption of an arcade field and the kinking of coronal loops is described. The braiding of the large-scale coronal field by convective photospheric motions develops fine-scale structure in the magnetic field and leads to the development of intense current filaments. The resistive dissipation of these currents can provide an efficient coronal heating mechanism.

Mikic, Zoran↗
NASA NTRS1993

Upwind MacCormack Euler solver with non-equilibrium chemistry

A computer code, designated UMPIRE, is currently under development to solve the Euler equations in two dimensions with non-equilibrium chemistry. UMPIRE employs an explicit MacCormack algorithm with dissipation introduced via Roe's flux-difference split upwind method. The code also has the capability to employ a point-implicit methodology for flows where stiffness is introduced through the chemical source term. A technique consisting of diagonal sweeps across the computational domain from each corner is presented, which is used to reduce storage and execution requirements. Results depicting one dimensional shock tube flow for both calorically perfect gas and thermally perfect, dissociating nitrogen are presented to verify current capabilities of the program. Also, computational results from a chemical reactor vessel with no fluid dynamic effects are presented to check the chemistry capability and to verify the point implicit strategy.

Sherer, Scott E.↗
NASA NTRS1987

Dynamical modeling of a planetary wave mechanism for a Martian polar warming

The mechanisms involved in the global dust storm and polar warming seen in the Martian atmosphere by the Viking IRTM during the winter solstice of 1977 are investigated theoretically by means of numerical simulations. A two-component dynamical model (based on the combined action of a zonally symmetric 'Hadley' circulation at low and middle latitudes and a planetary-wave circulation at middle and high latitudes) is constructed by analogy to the model of Holton and Mass (1976) for terrestrial sudden stratospheric warmings. The Viking data and simulation results are presented in extensive graphs and characterized in detail. It is demonstrated that a planetary-wave mechanism, based primarily on wavenumber 1 and including a high degree of topographical or thermal wave forcing, can reproduce the observed polar warming. The roles of radiative damping, dissipation, and the transport of dust and water are explored.

Barnes, Jeffrey R.↗
NASA NTRS2004

Simulating Virtual Terminal Area Weather Data Bases for Use in the Wake Vortex Avoidance System (Wake VAS) Prediction Algorithm

During the research project, sounding datasets were generated for the region surrounding 9 major airports, including Dallas, TX, Boston, MA, New York, NY, Chicago, IL, St. Louis, MO, Atlanta, GA, Miami, FL, San Francico, CA, and Los Angeles, CA. The numerical simulation of winter and summer environments during which no instrument flight rule impact was occurring at these 9 terminals was performed using the most contemporary version of the Terminal Area PBL Prediction System (TAPPS) model nested from 36 km to 6 km to 1 km horizontal resolution and very detailed vertical resolution in the planetary boundary layer. The soundings from the 1 km model were archived at 30 minute time intervals for a 24 hour period and the vertical dependent variables as well as derived quantities, i.e., 3-dimensional wind components, temperatures, pressures, mixing ratios, turbulence kinetic energy and eddy dissipation rates were then interpolated to 5 m vertical resolution up to 1000 m elevation above ground level. After partial validation against field experiment datasets for Dallas as well as larger scale and much coarser resolution observations at the other 8 airports, these sounding datasets were sent to NASA for use in the Virtual Air Space and Modeling program. The application of these datasets being to determine representative airport weather environments to diagnose the response of simulated wake vortices to realistic atmospheric environments. These virtual datasets are based on large scale observed atmospheric initial conditions that are dynamically interpolated in space and time. The 1 km nested-grid simulated datasets providing a very coarse and highly smoothed representation of airport environment meteorological conditions. Details concerning the airport surface forcing are virtually absent from these simulated datasets although the observed background atmospheric processes have been compared to the simulated fields and the fields were found to accurately replicate the flows surrounding the airport where coarse verification data were available as well as where airport scale datasets were available.

Kaplan, Michael L.↗
NASA NTRS1984

Formation and dynamics of large-scale magnetic structures in the ionosphere of Venus

The formation and dynamics of large-scale magnetic structures in the ionosphere of Venus are examined. It is shown that such structures must be the result of steady state convection of interplanetary field lines into the ionosphere by the small amount of solar wind plasma (less than or approximately equal to 1-5 percent) absorbed by the planetary atmosphere below the ionopause, rather than isolated remnants of large fields persisting for long periods without connection to the solar wind induced current and convection pattern. In particular, it is demonstrated that the magnetic diffusion of such structures would result in their dissipation with time scales of 1-10 min, if they were not steady state structures in convective and diffusive equilibriuim. It is shown that the equations governing the diffusion of these magnetic structures are similar to those governing diffusion of a gas out of an enclosed chamber with a porous wall, and a simple analog is illustrated. The application of these results to magnetic fields of astrophysical plasmas is discussed.

Cloutier, P. A.↗
NASA NTRS1999

The Role of Gravity Waves in Generating Equatorial Oscillations in Modulating Atmospheric Tides

We discuss a Numerical Spectral Mode (NSM) that extends from the ground up into the thermosphere and incorporates Hines' Doppler spread parameterization (DSP) for small scale gravity waves (GW). This model is applied to describe the seasonal variations in the mean zonal circulation, the semi-annual and quasi-biennial oscillations (SAO and QBO), as well as the tides and planetary waves in the middle atmosphere. Initial results showed that this model can reproduce the salient features observed, including the QBO extending into the upper mesosphere inferred from UARS measurements. The model has now been extended to simulate also: (a) the zonal circulation of the lower stratosphere and upper troposphere, and (b) the upwelling at equatorial latitudes associated with the Brewer Dobsen circulation that affects the dynamics significantly as pointed out by Dunkerton. Upward vertical winds increase the period of the QBO observed from the ground. To compensate for that, one needs to increase in the model the eddy diffusivity and the GW momentum flux, bringing the latter closer to values recommended in the DSP. This development is conducive to extending the QBO and SAO to higher latitudes through global scale momentum redistribution. Multi-year interannual oscillations are generated through wave filtering by the solar driven annual oscillation in the zonal circulation. In a 3D version of the model, wave momentum is absorbed and dissipated by tides and planetary waves. A somewhat larger GW source (well within the DSP range) is then required to generate realistic QBO and SAO amplitudes. Since GW momentum is deposited in the altitude regime of increasing winds, the amplitude of the diurnal tide is amplified and its vertical wavelength is reduced at altitudes between 70 and 120 km. Wave filtering by the mean zonal circulation causes the GW flux to peak during equinox, and this produces a large semi-annual variation in the tide that has been observed on UARS. Without the diurnal tide, the semidiurnal tide would also be modulated in this way. But the diurnal tide filters out the GW preferentially during equinox, so that the semidiurnal tide tends to peak during solstice. Under the influence of GW, the tides are modulated significantly by planetary waves that are generated preferentially during solstice in part due to baroclinic instability.

Mayr, H. G.↗
NASA NTRS1985

Dynamics of the earth's ring current - Theory and observation

The development of currents within an arbitrary distribution of particles trapped in the geomagnetic field is described. These currents combine to form the earth's ring current and thus are responsible for the worldwide depressions of surface magnetic field strength during periods of magnetic activity known as magnetic storms. Following a brief review of trapped particle motion in magnetic fields, ring current development is described and presented in terms of basic field and particle distribution parameters. Experimental observations then are presented and discussed within the theoretical framework developed earlier. New results are presented which, in the area of composition and charge state observations, hold high promise in solving many long standing ring current problems. Finally, available experimental results will be used to assess the present understanding as to ring current sources, generation, and dissipation.

Williams, D. J.↗