Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dissipative systems”

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 865 records · Page 48

On-orbit performance testing of the Pointing Calibration & Reference Sensor for the Spitzer Space Telescope

We present the on-orbit performance results of the Pointing Calibration and Reference Sensor (PCRS) for the Spitzer Space Telescope. A cryogenic optical (center wavelength 0.55 mu) imager, the PCRS serves as the Observatory's fine guidance sensor by providing an alignment reference between the telescope boresight and the external spacecraft attitude determination system. The PCRS makes precision measurements of the positions of known guide stars; these are used to calibrate measurements from Spitzer's star tracker and gyroscopes to obtain the actual pointing of the Spitzer telescope. The PCRS calibrates out thermomoechanical drifts between the 300 K spacecraft bus and the 5.5 K telescope. By using only 16 pixels, the PCRS provides high precision centroiding with extremely low (`64 mu W) power dissipation, resulting in minimal impact to Spritzer's helium lifetime. We have demonstrated that the PCRS meets its centroiding accuracy requirement of 0.14 arcsec 1-sigma radial, which represents about1/100 pixel centroiding. The Spitzer Space Telescope was launched in 25 August, 2003 and completed its In-Orbit Checkout phase two months later; the PCRS has been operating failure free ever since.

Spitzer↗

Sea State Based Estimation of White Cap Fraction: Implications for Primary Marine Aerosol Fluxes

Oceanic whitecaps (hereafter, W) or the characteristic whiteness of the sea foam is an important feature for predicting exchange of gases, sea spray aerosols (SSAs), heat and momentum transfer between the ocean and the atmosphere at the air-sea interface. Due to its increased surface emission and brightness temperature, whitecaps are critical for satellite retrievals of ocean albedo, ocean color, ocean surface wind vectors from satellite borne radiometer and microwave instruments. Most of the existing models predict W using wind speed and sea surface temperature (SST). However, numerous publications have pointed out that there are large uncertainties in the predicted W and using parameterizations based on wind-wave state can improve the precision of the predicted W. Here, we integrate the University of Miami Wave Model - 2.0 (UMWM) in Goddard Earth Observing System (GEOS) and use wave diagnostics to predict W. We choose the year 2006 for our global UMWM/GEOS runs because of the availability of W dataset from satellite observations. We run UMWM/GEOS at 0.5o x 0.5o by replaying to MERRA2 meteorology and evaluate the wave diagnostics using measurements from fixed buoys and satellite altimeters. We use three different parameterizations for W based on: 1) Reynolds number, 2) wave dissipation energy, and 3) volume of air entrained by breaking waves. We compare our results of W with previous studies and also with the satellite based observational dataset. Predicting W is important for understanding the processes at the air-sea interface. Therefore, this work is a step further in improving the uncertainties in the aerosol and atmospheric chemistry modules of the global models.

Oceanic Whitecaps↗

PROGRESS REPORT ON THE DEVELOPMENT OF PROTECTED CONSTRUCTION FOR HYPERSONIC VEHICLES

The structural problems of re-entry occasioned by aerodynamic heating are now generally well known. Typical of these heating problems is the environment experienced by the manned lifting vehicle re-entering from a low altitude orbit. Figure 1 shows typical curves of lower surface temperature as a function of time for a wing loading of 25 Ibs/ft(exp 2), and a lift-drag ratio of 2.5. The two curves cover a practical range of re-entry angles, the lower curve representing an ideal re-entry path with zero dive angle and the upper curve assuming a 2° error in the re-entry angle. Maneuvers for course change or correction are also added to the upper curve. Significant factors from this curve are the relationship between the maximum temperatures and the capabilities of available metallic material, and also the long re-entry time and its effect on total heat load. If a very shallow re-entry is made, maximum lower surface temperatures reach 2000°F which is just within the range of conventional superalloys. To accommodate practical re-entry angles and maneuvers, however, the temperature capability must reach about 2500°F which requires refractory metals. The re-entry time may be as high as 100 minutes, which gives total heat loads of approximately 40,000 BTU/ft(exp 2). A heat load of this magnitude, with equilibrium surface temperatures of the values shown, suggests that a lighter airframe can be constructed by dissipating the heat by radiation from the surface, rather than by absorbing it with a heat sink, or surface cooling, or ablation. Air Force programs to provide airframes for this type of environment have involved the parallel development of a number of different structural concepts. The development to be discussed here was carried out by Bell Aerosystems Company for the Fabrication Branch, Manufacturing Technology Laboratory, Directorate of Materials and Processes, Aeronautical Systems' Division, Wright-Patterson Air Force Base, Ohio, under Contract AF33(600)-40100 (Double-Wall Cooled Structure).

Reentry vehicle↗

Geometric radiating surfaces for spacecraft thermal control

Various geometric radiating surface (GRS) configurations are evaluated analytically and experimentally for potential use in spacecraft thermal control systems. Egg-crate, honeycomb, flat-parallel, curved-parallel, chevron-shape, and composite fin radiator panel designs were evaluated in terms of mass, apparent emissivity, meteoroid protection, and ease of fabrication; the parallel-finned GRS in composite form was chosen as the most efficient design. The performance of the composite-finned GRS was compared with that of a silver-backed Teflon radiator. The emissivity and dissipative power capabilities of the GRS and radiator are measured. It is observed that the composite-finned GRS has 34 percent less surface area, 34 percent higher heat flux capability, greater meteoroid protection, and 2.8 times higher mass than the Teflon radiator. Future developments and applications for GRCs are discussed.

Keddy, M. D.↗

An active homopolar magnetic bearing with high temperature superconductor (HTS) coils and ferromagnetic cores

A proof-of-feasibility demonstration showed that high temperature superconductor (HTS) coils can be used in a high-load, active magnetic bearing in liquid nitrogen. A homopolar radial bearing with commercially wound HTS (Bi 2223) bias and control coils produced over 200 lb (890 N) radial load capacity (measured non-rotating) and supported a shaft to 14000 rpm. The goal was to show that HTS coils can operate stably with ferromagnetic cores in a feedback controlled system at a current density similar to that in Cu in liquid nitrogen. Design compromises permitted use of circular coils with rectangular cross section. Conductor improvements will eventually permit coil shape optimization, higher current density and higher bearing load capacity. The bias coil, wound with non-twisted, multifilament HTS conductor, required negligible power to carry its direct current. The control coils were wound with monofilament HTS sheathed in Ag. These dissipated negligible power for direct current (i.e. for steady radial load components). When an alternating current (AC) was added, the AC component dissipated power which increased rapidly with frequency and quadratically with AC amplitude. In fact at frequencies above about 2 hz, the effective resistance of the control coil conductor actually exceeds that of the silver which is in electrical parallel with the oxide superconductor. This is at least qualitatively understandable in the context of a Bean-type model of flux and current penetration into a Type II superconductor. Fortunately the dynamic currents required for bearing stability are of small amplitude. These results show that while twisted multifilament conductor is not needed for stable levitation, twisted multifilaments will be required to reduce control power for sizable dynamic loads, such as those due to unbalance.

Brown, G. V.↗

Diagrammatic theory of transition of pendulum like systems

Orbit-orbit and spin-orbit gravitational resonances are analyzed using the model of a rigid pendulum subject to both a time-dependent periodic torque and a constant applied torque. First, a descriptive model of passage through resonance is developed from an examination of the polynomial equation that determines the extremes of the momentum variable. From this study, a probability estimate for capture into libration is derived. Second, a lowest order solution is constructed and compared with the solution obtained from numerical integration. The steps necessary to systematically improve this solution are also discussed. Finally, the effect of a dissipative term in the pendulum equation is analyzed.

Yoder, C. F.↗

Active nutation controller

An apparatus is described for controlling nutation motion in a spinning body, comprised of an angular accelerometer with its input axis perpendicular to the spin axis of the body, a flywheel with an axis of rotation perpendicular to the axis of the accelerometer and to the spin axis of the body, and a motor for driving the flywheel to attenuate or build nutation. The motor is controlled by circuitry that monitors the output of the angular accelerometer and drives the motor clockwise or counterclockwise during predetermined nutation angles synchronized to the zero crossover points of the accelerometer signal centered about the nutation peaks. The motor drive is phased to damp nutation motion to zero for stabilization. To increase the noise immunity of the system, when the output of the accelerometer falls below a threshold level, the circuitry operates in an open loop, beat mode where data representing the last accelerometer signal that exceeded that threshold level is stored, and the motor drive is controlled by the stored data. In a second version, the motor is controlled to supply a predetermined amount of nutation motion to a body undergoing testing on a spin table for energy dissipation evaluation. In each version, the use of an angular accelerometer rather than a linear accelerometer or gyro to monitor nutation enables placement of the nutation control apparatus at any location relative to the spin axis of the body requiring only crude orientation and no calibration.

Hoffman, H. C.↗

Statistical energy analysis of a geared rotor system

The vibroacoustic response of a generic geared rotor system is analyzed on an order of magnitude basis utilizing an approximate statistical energy analysis method. This model includes a theoretical coupling loss factor for a generic bearing component, which properly accounts for the vibration transmission through rolling element bearing. A simplified model of a NASA test stand that assumes vibratory energy flow from the gear mesh source to the casing through shafts and bearings is given as an example. Effects of dissipation loss factor and gearbox radiation efficiency models are studied by comparing predictions with NASA test results.

Lim, Teik C.↗

Wake Vortex Algorithm Scoring Results

This report compares the performance of two models of trailing vortex evolution for which interaction with the ground is not a significant factor. One model uses eddy dissipation rate (EDR) and the other uses the kinetic energy of turbulence fluctuations (TKE) to represent the effect of turbulence. In other respects, the models are nearly identical. The models are evaluated by comparing their predictions of circulation decay, vertical descent, and lateral transport to observations for over four hundred cases from Memphis and Dallas/Fort Worth International Airports. These observations were obtained during deployments in support of NASA's Aircraft Vortex Spacing System (AVOSS). The results of the comparisons show that the EDR model usually performs slightly better than the TKE model.

Robins, R. E.↗

Observed tidal braking in the earth/moon/sun system

The low degree and order terms in the spherical harmonic model of the tidal potential were observed through the perturbations which are induced on near-earth satellite orbital motions. Evaluations of tracking observations from 17 satellites and a GEM-T1 geopotential model were used in the tidal recovery which was made in the presence of over 600 long-wavelength coefficients from 32 major and minor tides. Wahr's earth tidal model was used as a basis for the recovery of the ocean tidal terms. Using this tidal model, the secular change in the moon's mean motion due to tidal dissipation was found to be -25.27 + or - 0.61 arcsec/century-squared. The estimation of lunar acceleration agreed with that observed from lunar laser ranging techniques (-24.9 + or - 1.0 arcsec/century-squared), with the corresponding tidal braking of earth's rotation being -5.98 + or - 0.22 X 10 to the -22 rad/second-squared. If the nontidal braking of the earth due to the observed secular change in the earth's second zonal harmonic is considered, satellite techniques yield a total value of the secular change in the earth's rotation rate of -4.69 + or - 0.36 X 10 to the -22 rad/second-squared.

Christodoulidis, D. C.↗

Observed tidal braking in the earth/moon/sun system

The low degree and order terms in the spherical harmonic model of the tidal potential were observed through the perturbations which are induced on near-earth satellite orbital motions. Evaluations of tracking observations from 17 satellites and a GEM-T1 geopotential model were used in the tidal recovery which was made in the presence of over 600 long-wavelength coefficients from 32 major and minor tides. Wahr's earth tidal model was used as a basis for the recovery of the ocean tidal terms. Using this tidal model, the secular change in the moon's mean motion due to tidal dissipation was found to be -25.27 + or - 0.61 arcsec/century squared. The estimation of lunar acceleration agreed with that observed from lunar laser ranging techniques (-24.9 + or - 1.0 arcsec/century squared), with the corresponding tidal braking of earth's rotation being -5.98 + or - 0.22 x 10 to the minus 22 rad/second squared. If the nontidal braking of the earth due to the observed secular change in the earth's second zonal harmonic is considered, satellite techniques yield a total value of the secular change of the earth's rotation rate of -4.69 + or - 0.36 x 10 to the minus 22 rad/second squared.

Christodoulidis, D. C.↗

A streamwise upwind algorithm for the Euler and Navier-Stokes equations applied to transonic flows

A new algorithm was developed for the Euler and Navier-Stokes equations that uses upwind differencing based on the streawise direction. This algorithm is time accurate and can be used in codes for calculating unsteady transonic flows over wings. Such codes can be used for the flutter analysis of wings. In this algorithm, the coordinate system is locally rotated to align with the streamwise direction. For differencing the convective terms in the streamwise direction, a new form of flux splitting is employed, in which the biasing depends on the local Mach number. In the plane perpendicular to the stream direction, the new flux splitting uses the condition of no flow in that local plane. By using a locally rotated coordinate system, the convective flux vector biasing depends on the total Mach number. Hence, the switching of the flux vector biasing occurs across shock waves and the proper domain of dependence is used in supersonic regions. For comparison, many other upwind methods switch differencing based on Mach number of shock waves in multidimensional flows. The formulas for the convective flux vector differencing do not contain any user specified parameters. So, the amount of numerical dissipation is automatically determined.

Goorjian, Peter M.↗

Europa, tidally heated oceans, and habitable zones around giant planets

Tidal dissipation in the satellites of a giant planet may provide sufficient heating to maintain an environment favorable to life on the satellite surface or just below a thin ice layer. Europa could have a liquid ocean which may occasionally receive sunlight through cracks in the overlying ice shell. In such a case, sufficient solar energy could reach liquid water that organisms similar to those found under Antarctic ice could grow. In other solar systems, larger satellites with more significant heat flow could represent environments that are stable over an order of eons and in which life could perhaps evolve. A zone around a giant planet is defined in which such satellites could exist as a tidally-heated habitable zone. This zone can be compared to the habitable zone which results from heating due to the radiation of a central star. In this solar system, this radiatively-heated habitable zone contains the earth.

Reynolds, Ray T.↗

A streamwise upwind algorithm for the Euler and Navier-Stokes equations applied to transonic flows

A new algorithm was developed for the Euler and Navier-Stokes equations that uses upwind differencing based on the streamwise direction. This algorithm is time accurate and can be used in codes for calculating unsteady transonic flows over wings. Such codes can be used for the flutter analysis of wings. In this algorithm, the coordinate system is locally rotated to align with the streamwise direction. For differencing the convective terms in the streamwise direction, a new form of flux splitting is employed, in which the biasing depends on the local Mach number. In the plane perpendicular to the stream direction, the new flux splitting uses the condition of no flow in that local plane. By using a locally rotated coordinate system, the convective flux vector biasing depends on the total Mach number. Hence, the switching of the flux vector biasing occurs across shock waves and the proper domain of dependence is used in supersonic regions. For comparison, many other upwind methods switch differencing based on Mach number of shock wave in multidimensional flows. The formulas for the convective flux vector differencing do not contain any user specified parameters. So, the amount of numerical dissipation is automatically determined.

Goorjian, P. M.↗

Thermal and dynamical regimes of single- and two-phase magmatic flow in dikes

The coupling between thermal and dynamical regimes of single- and two-phase magmatic flow in dikes, due to temperature-dependent viscosity and dissipation, was investigated using finite element calculations of magma flow in dikelike channels with length-to-width ratios of 1000:1 or more. Solutions of the steady state equations governing magma flow are obtained for a variety of conditions ranging from idealized plane-parallel models to cases involving nonparallel geometry and two-phase flows. The implications of the numerical simulations for the dynamics of flow in a dike-reservoir system and the consequences of dike entrance conditions on magmatic storage are discussed. Consideration is also given to an unmixing/self-lubrication mechanism which may be important for the lubrication of silicic magmas rising to the earth's surface in mixed magma ascent scenarios, which naturally segregates magma mixtures of two components with differing viscosities to minimize the driving pressure gradient.

Carrigan, Charles R.↗

Comparative Study on High-Order Positivity-preserving WENO Schemes

In gas dynamics and magnetohydrodynamics flows, physically, the density and the pressure p should both be positive. In a standard conservative numerical scheme, however, the computed internal energy is obtained by subtracting the kinetic energy from the total energy, resulting in a computed p that may be negative. Examples are problems in which the dominant energy is kinetic. Negative may often emerge in computing blast waves. In such situations the computed eigenvalues of the Jacobian will become imaginary. Consequently, the initial value problem for the linearized system will be ill posed. This explains why failure of preserving positivity of density or pressure may cause blow-ups of the numerical algorithm. The adhoc methods in numerical strategy which modify the computed negative density and/or the computed negative pressure to be positive are neither a conservative cure nor a stable solution. Conservative positivity-preserving schemes are more appropriate for such flow problems. The ideas of Zhang & Shu (2012) and Hu et al. (2012) precisely address the aforementioned issue. Zhang & Shu constructed a new conservative positivity-preserving procedure to preserve positive density and pressure for high-order WENO schemes by the Lax-Friedrichs flux (WENO/LLF). In general, WENO/LLF is too dissipative for flows such as turbulence with strong shocks computed in direct numerical simulations (DNS) and large eddy simulations (LES). The new conservative positivity-preserving procedure proposed in Hu et al. (2012) can be used with any high-order shock-capturing scheme, including high-order WENO schemes using the Roe's flux (WENO/Roe). The goal of this study is to compare the results obtained by non-positivity-preserving methods with the recently developed positivity-preserving schemes for representative test cases. In particular the more difficult 3D Noh and Sedov problems are considered. These test cases are chosen because of the negative pressure/density most often exhibited by standard high-order shock-capturing schemes. The simulation of a hypersonic nonequilibrium viscous shock tube that is related to the NASA Electric Arc Shock Tube (EAST) is also included. EAST is a high-temperature and high Mach number viscous nonequilibrium flow consisting of 13 species. In addition, as most common shock-capturing schemes have been developed for problems without source terms, when applied to problems with nonlinear and/or sti source terms these methods can result in spurious solutions, even when solving a conservative system of equations with a conservative scheme. This kind of behavior can be observed even for a scalar case (LeVeque & Yee 1990) as well as for the case consisting of two species and one reaction (Wang et al. 2012). For further information concerning this issue see (LeVeque & Yee 1990; Griffiths et al. 1992; Lafon & Yee 1996; Yee et al. 2012). This EAST example indicated that standard high-order shock-capturing methods exhibit instability of density/pressure in addition to grid-dependent discontinuity locations with insufficient grid points. The evaluation of these test cases is based on the stability of the numerical schemes together with the accuracy of the obtained solutions.

WENO↗

An approximate analytic solution of a set of nonlinear model alpha-omega-dynamo equations for marginally unstable systems

An approximate analytic solution of a set of nonlinear model alpha-omega-dynamo equations is obtained. The reaction of the Lorentz force on the velocity shear which stretches and, hence, amplifies the magnetic field is incorporated into the model. To single out the effect of the Lorentz force on the omega-effect, the effect of the Lorentz force on the alpha-effect is neglected in this study. The solution represents a nonlinear oscillation with the amplitude and period determined by the dynamo number N. The amplitude is proportional to N - 1, while the period is almost exactly the same as the dissipation time of the unstable mode (proportional to N).

Hinata, S.↗

The interaction between giant gaseous protoplanets and the primitive solar nebula

The manner in which a giant gaseous protoplanet becomes embedded in the primitive solar nebula determines surface boundary conditions which must be used in studying the evolution of such objects. On the one hand, if the system resembles a contact binary system, then the envelope of the protoplanet should approach the entropy of the surrounding nebula. On the other hand angular momentum transfer by resonance and tidal effects between the nebula and the protoplanet may cause the nebula to exhibit a zone of avoidance near the protoplanet, thus inhibiting exchange of material. This problem has been studied with a computer program developed by D. N. C. Lin which simulates disk hydrodynamics by particle motions with dissipation. These studies suggest that for expected values of the protoplanet/protosun mass ratios, significant inhibition of mass exchange is likely, so that it is a reasonable next step to undertake protoplanet evolution studies with the imposition of minimum protoplanet surface temperatures.

Cameron, A. G. W.↗