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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

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.↗
NASA NTRS1995

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.↗
NASA NTRS1979

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.↗
DOE OSTI2020

Simulation Analysis of Geometry and Material Effects for Dropkinson Bar

The reported research is motivated by the need to address a key issue affecting the Dropkinson bar apparatus. This unresolved issue is the interference of the stress wave reflected from the bar-beam boundary with the measurement of the stress-strain response of a material tested in the apparatus. The purpose of the wave beam that is currently connected to the bar is to dissipate the stress wave, but the portion of the wave reflected from the bar-beam boundary is still significant. First, we focused on understanding which parameters affect the reflected wave's arrival time at a strain gauge. Specifically, we used finite-element numerical simulations with the Sierra/SM module to study the effects of various bar-beam connection fixities, alternative wave beam materials, and alternative geometries of the Dropkinson bar system based on a monolithic design. The conclusion of this study is that a partial reflection always occurs at the bar-beam boundary (or, for a monolithic design, at a point where the bar geometry changes). Therefore, given a fixed total length of the bar, it is impossible to increase the reflected wave's arrival time by any significant amount. After reaching this conclusion, we focused instead on trying to minimize the energy of the reflected stress wave circulating up and down through the bar over a relatively long period of time (10 ms). Once again, we used numerical simulations with the Sierra/SM module to investigate the effects of various bar-beam connection fixities, alternative wave beam materials, and parameters of an asymmetric monolithic design of the bar-and-beam system. This study demonstrated that various parameters can significantly affect the energy of the wave reflections, with the difference between best and worst configurations being about one order of magnitude in terms of energy. Based on the obtained results, we conclude with concrete takeaways for Dropkinson bar users and propose potential directions for future research and optimization.

36 MATERIALS SCIENCE↗
NASA NTRS1980

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.↗
NASA NTRS1990

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.↗
NASA NTRS2002

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.↗
DOE OSTI2021

Self-powered Heating: Efficiency Analysis

Conventional fuel-fired heating devices such as furnaces, boilers, and water heaters have fuel efficiency less than 100% on the basis of higher heating value. They also require electricity from the electric grid to power parasitic loads such as blowers, pumps, fans, and ignitors. The primary energy efficiency of the device accounts for both fuel used on-site and primary energy used off-site to produce electric power used by the device. This work compares conventional fuel-fired heating devices to two types of self-powered devices. A self-powered device (SPD) integrates a power cycle onboard to eliminate consumption of grid electricity. We assume that all heat rejected by the onboard power cycle is added to the process fluid, so that, compared with a conventional device, the same amount of heat is provided to the process fluid and the same amount of fuel is consumed, but grid electricity consumption is eliminated. The first SPD type is the basic one: exactly the electricity required is generated. The second type considered is the SPD with heat pump (SPD-HP), in which the power cycle generates more electricity than needed for parasitic loads, and the excess electricity is used to power a heat pump. The heat pump extracts additional heat from the ambient to boost efficiency. Both SPD and SPD-HP self-consume all the generated electricity, in contrast to combined heat and power (CHP) systems that export electricity. In this work, equations are derived to express the efficiency of three classes of heating devices: conventional (consuming grid electricity), self-powered (consuming no grid electricity), and self-powered with heat pump. The efficiency of each is derived as a function of up to six factors: (1) the fraction of combustion heat captured, (2) the rate of parasitic power consumption, (3) the fraction of electric energy dissipated as useful heat, (4) the power cycle conversion efficiency, (5) the grid efficiency, when applicable, and (6) the heat pump COP, when applicable. Scenarios are identified in which it is possible to achieve efficiency greater than 100% on a higher heating value basis. Plausible configurations using existing technology options are outlined.

Gluesenkamp, Kyle↗
NASA NTRS1988

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.↗
NASA NTRS1987

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.↗
DOE OSTI2024

Modeling hyperbranched polymer formation via ATRP using dissipative particle dynamics

Hyperbranched polymers (HBPs) offer distinguishing, advantageous properties that arise from their distinctive complex topology. One of the effective approaches to the synthesis of hyperbranched structures involves the use of a branching initiator (inibramer) that is activated only after incorporation into a polymer chain. There remain, however, challenges in determining and characterizing the structures of the synthesized HBPs. Dissipative particle dynamics (DPD) was used to probe the effects of inibramer concentration, solvent concentration, and inibramer reactivity on the kinetics, molecular weight, and dispersity of HBPs. Additionally, DPD allows for direct observation of branched structures, which was not possible in previously reported Monte Carlo type simulations. It was found that higher inibramer concentrations led to faster monomer consumption while forming more dendritic structures with fewer defects. Additionally, high dispersities characteristic of HBP systems were found to originate from asymmetric propagation rates between inibramer-inibramer and inibramer-monomer reactions.

Biswas, Santidan↗
DOE OSTI2022

Reconciling and Improving Formulations for Thermodynamics and Conservation Principles in Earth System Models (ESMs)

Abstract This paper provides a comprehensive derivation of the total energy equations for the atmospheric components of Earth System Models (ESMs). The assumptions and approximations made in this derivation are motivated and discussed. In particular, it is emphasized that closing the energy budget is conceptually challenging and hard to achieve in practice without resorting to ad hoc fixers. As a concrete example, the energy budget terms are diagnosed in a realistic climate simulation using a global atmosphere model. The largest total energy errors in this example are spurious dynamical core energy dissipation, thermodynamic inconsistencies (e.g., coupling parameterizations with the host model) and missing processes/terms associated with falling precipitation and evaporation (e.g., enthalpy flux between components). The latter two errors are not, in general, reduced by increasing horizontal resolution. They are due to incomplete thermodynamic and dynamic formulations. Future research directions are proposed to reconcile and improve thermodynamics formulations and conservation principles.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1988

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.↗
NASA NTRS1987

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.↗
NASA NTRS1988

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.↗
NASA NTRS1992

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.↗
NASA NTRS2013

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↗
NASA NTRS1989

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.↗