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 613 records · Page 34

Global Energetics of Solar Flares. VII. Aerodynamic Drag in Coronal Mass Ejections

The free energy that is dissipated in a magnetic reconnection process of a solar flare, generally accompanied by a coronal mass ejection (CME), has been considered as the ultimate energy source of the global energy budget of solar flares in previous statistical studies. Here we explore the effects of the aerodynamic drag force on CMEs, which supplies additional energy from the slow solar wind to a CME event, besides the magnetic energy supply. For this purpose, we fit the analytical aerodynamic drag model of Cargill and Vršnak et al. to the height–time profiles r(t) of LASCO/SOHO data in 14,316 CME events observed during the first 8 yr (2010–2017) of the Solar Dynamics Observatory era (ensuring EUV coverage with AIA). Our main findings are (1) a mean solar wind speed of w = 472 ± 414 km s(exp −1), (2) a maximum drag-accelerated CME energy of E(drag) <~2 × 10(exp32) erg, (3) a maximum flare-accelerated CME energy of E(flare <~1.5 × 10(exp33) erg, (4) the ratio of the summed kinetic energies of all flare accelerated CMEs to the drag-accelerated CMEs amounts to a factor of 4, (5) the inclusion of the drag force slightly lowers the overall energy budget of CME kinetic energies in flares from ≈7% to ≈4%, and (6) the arrival times of CMEs at Earth can be predicted with an accuracy of ≈23%.

Aschwanden, Markus J.↗

Three-dimensional MHD modeling of magnetotail dynamics for different polytropic indices

The paper analyzes the 3D resistive tearing instability in a magnetotail configuration by means of a resistive MHD code for various forms of the energy equation, including ohmic and compressional heating. A complete derivation of the set of MHD equations from a kinetic formalism is presented, with all simplifying assumptions noted. Emphasis is placed on the changes in magnetotail evolution caused by a variation of the assumed value of the polytropic exponent gamma, representing the ratio of specific heat in the energy equation which closes the set of MHD equations. Within the gamma range of 2/3-4, a tearing mode is found to develop in all cases. Certain properties, such as the formation and ejection of a plasmoid, associated with the formation of the reversal region of the north-south component of the magnetic field, and strong earthward and tailward flow fields prevail as qualitative features independent of the value of gamma chosen. Ohmic dissipation is found to be the major source of the heating, while adiabatic or nonadiabatic compression may contribute primarily in the closed field region earthward of the reconnection site.

Hesse, Michael↗

Dynamic Models for LES of Turbulent Front Propagation With a Spectral Method

Direct numerical simulation of turbulent reacting flows places extreme demands on computational resources. At the present time, simulations can be performed only for greatly simplified reaction systems and for very low Reynolds numbers. Direct simulation of more realistic cases occurring at higher Reynolds number and including multiple species and numerous chemical reactions will exceed available computational resources far into the future. Because of this, there is a clear need to develop the technique of large eddy simulation for reacting flows. Unfortunately this task is complicated by the fact that combustion arises from chemical reactions that occur at the smallest scales of the flow. Capturing the large-scale behavior without resolving the small-scale details is extremely difficult in combustion problems. Thus LES modeling for turbulent combustion encounters difficulties not present in modeling momentum transport, in which the main effect of the small scales is to provide dissipation. The difficulty is more pronounced in premixed combustion, where detailed chemistry plays an essential role in determining the flame speed (or overall burning rate); in nonpremixed combustion infinite rate chemistry can be assumed, eliminating the small scale features to a first approximation.

Im, H. G.↗

Structural Damage Detection Using Virtual Passive Controllers

This paper presents novel approaches for structural damage detection which uses the virtual passive controllers attached to structures, where passive controllers are energy dissipative devices and thus guarantee the closed-loop stability. The use of the identified parameters of various closed-loop systems can solve the problem that reliable identified parameters, such as natural frequencies of the open-loop system may not provide enough information for damage detection. Only a small number of sensors are required for the proposed approaches. The identified natural frequencies, which are generally much less sensitive to noise and more reliable than the identified natural frequencies, are used for damage detection. Two damage detection techniques are presented. One technique is based on the structures with direct output feedback controllers while the other technique uses the second-order dynamic feedback controllers. A least-squares technique, which is based on the sensitivity of natural frequencies to damage variables, is used for accurately identifying the damage variables.

Lew, Jiann-Shiun↗

Diurnal Forcing of Planetary Atmospheres

This project investigated and developed models which will assist in developing an understanding of the diurnal forcings in the atmospheres of Mars, Venus, Earth and Jupiter. The differences in the temperatures during the day, and night, generate very deep diurnal boundary layers and strong slope-related winds. A general circulation model (GCM) was adapted for computation of Martian atmospheric dynamics, to simulate the global transport of these atmospheric tracers. The model was applied to the study of the Martian seasonal water cycle. The principles have been applied to venus. The region of study is the upper cloud level, and above rather than the area close to the surface. The model was also applied in efforts to study the evolution of volcanic clouds in the Earth's stratosphere. In particular this was used to model the dispersion of sulfuric acid particles in the atmosphere. A brief study of the application of the tidal theory to the planet Jupiter has been undertaken. The goal of the study is to deduce the source of the tidal dissipation which is thought to have led to the orbital resonances amongst the Galilean satellites and considerable heating of the interiors of Io, Europa, and Ganymede.

Houben, Howard C.↗

Dynamics of the Final Stages of Terrestrial Planet Growth and the Formation of the Earth-Moon System

An overview of current theories of star and planet formation, with emphasis on terrestrial planet accretion and the formation of the Earth-Moon system is presented. These models predict that rocky planets should form around most single stars, although it is possible that in some cases such planets are lost to orbital decay within the protoplanetary disk. The frequency of formation of gas giant planets is more difficult to predict theoretically. Terrestrial planets are believed to grow via pairwise accretion until the spacing of planetary orbits becomes large enough that the configuration is stable for the age of the system. Giant impacts during the final stages of growth can produce large planetary satellites, such as Earth's Moon. Giant planets begin their growth like terrestrial planets, but they become massive enough that they are able to accumulate substantial amounts of gas before the protoplanetary disk dissipates.

Lissauer, Jack J.↗

Advancements in engineering turbulence modeling

Some new developments in two-equation models and second order closure models are presented. Two-equation models (k-epsilon models) have been widely used in computational fluid dynamics (CFD) for engineering problems. Most of low-Reynolds number two-equation models contain some wall-distance damping functions to account for the effect of wall on turbulence. However, this often causes the confusion and difficulties in computing flows with complex geometry and also needs an ad hoc treatment near the separation and reattachment points. A set of modified two-equation models is proposed to remove the aforementioned shortcomings. The calculations using various two-equation models are compared with direct numerical simulations of channel flow and flat boundary layers. Development of a second order closure model is also discussed with emphasis on the modeling of pressure related correlation terms and dissipation rates in the second moment equations. All the existing models poorly predict the normal stresses near the wall and fail to predict the 3-D effect of mean flow on the turbulence (e.g. decrease in the shear stress caused by the cross flow in the boundary layer). The newly developed second order near-wall turbulence model is described and is capable of capturing the near-wall behavior of turbulence as well as the effect of 3-D mean flow on the turbulence.

Shih, T.-H.↗

A Dynamic Nonlinear Subgrid-Scale Model for Large-Eddy Simulation of Complex Turbulent Flows

We present a new dynamic nonlinear subgrid-scale (SGS) model for large-eddy simulations (LES) and apply it to compute a flow involving pressure gradients, surface curvature and separation, for which data from a direct numerical simulation are available for comparison. The model, inspired by the triple model idea of Bardina et al. (“Improved Turbulence Models Based on Large Eddy Simulation of Homogeneous, Incompressible, Turbulent Flows,” Report No. TF-19, Thermosciences Division, Department of Mechanical Engineering, Stanford University, 1983), includes a Galilean-invariant term called the modified Leonard stress tensor, and two nonlinear terms comprised of the products of the strain-rate and rotation-rate tensors for an improved representation of the subgrid-scale dissipation, backscatter and anisotropy effects. The model does not employ any ad hoc averaging or clipping procedures, and does not require the specification of a characteristic length scale; hence, it naturally avoids the ambiguities associated with defining a proper length scale for anisotropic grids. Results from the wall-resolved LES of flow past a Gaussian bump using the new model demonstrate improved prediction of skin-friction, flow separation, mean flow profiles and turbulent quantities when compared to implicit LES as well as explicit LES using the Vreman SGS model on the same grid.

Turbulence↗

Evaluation of Limiter Functions for Supersonic Applications

Limiters commonly used in the simulations of flows with discontinuities are compared with the new limiter function proposed by Nishikawa using idealized test cases in two dimensions as well as complex three-dimensional problems. The Nishikawa limiter is observed to be consistently the least dissipative in idealized test cases as well as complex practical problems, for both steady and time-dependent problems. In the case of steady simulations, its iterative convergence characteristics are either similar or better than other limiter functions. The nearfield sonic boom signature of a low-boom demonstrator is computed to demonstrate the utility of the Nishikawa limiter function for realistic supersonic aircraft configurations.

Computational fluid dynamics↗

Kinetic aspects of tail dynamics - Theory and simulation

Kinetic theories relevant to the geomagnetic tail are reviewed. The topics discussed include kinetic instabilities, simulations, and current-sheet particle acceleration. Tearing mode and reconnection theories are emphasized. Kinetic treatment is appropriate for these topics since the tail plasma is collisionless. Fluid calculations are appropriate when stochastic processes dominate and for studies where long wavelengths are important. However, fluid treatments of tearing modes and reconnection require a finite resistivity in the diffusion region. Thus, although 'anomalous resistivity' can be guessed or in some cases calculated, ideally the kinetic treatment is often to be preferred. Particle motion and acceleration in the current sheet can give rise to beam-like distributions in the plasma-sheet boundary layer. Studies of current-sheet particle motion have also been used as the basis for 'kinetic' tail equilibrium models. Furthermore, quite recently current-sheet particle motion is used directly in Coroniti's explosive tail reconnection model. The 'inertial conductivity' from the equilibrium models provides the 'dissipation' necessary for reconnection.

Speiser, T. W.↗

Mechanisms of Interfacial Mixing in the Blown Powder Deposition of Inconel 625 to Copper Alloy for Liquid Rocket Engine Applications

There is an increasing need for functional materials to support high heat flux applications and optimize weight by using bimetallic additive manufacturing (AM). While one promising fabrication technology is the use of blown powder directed energy deposition (DED) AM, challenges remain with bi-metallic combinations. Prominent among these challenges is the need to understand mixing of materials at the bimetallic interface for forming stronger, more durable and reliable joints to extend the versatility of AM into critical applications. One application is the manufacturing of combustion chambers with copper-alloy liners and Superalloy structural jackets. In this application, material combinations are selected for effective heat dissipation, using copper, and retention of strength at high temperatures using Inconel 625. Currently, the dominant problem with this configuration is the highly variable strength of the interface and its correlation with deposition parameters. In this study, samples were fabricated at several vendors with Inconel 625 deposited onto wrought C-18150 (Cu-Cr-Zr) to form a bimetallic interface. Characterization of the interfaces were achieved through optical and electron microscopy, mechanical testing of mini tensile specimens, and electron dispersion spectroscopy to understand diffusion and bulk mixing behaviors at the joint. A theoretical framework was constructed in the context of solidification and mixing mechanisms to explain the resulting differences in the interface. This qualitative approach utilized a fluid dynamics perspective to sum the order of magnitudes of the various acting forces on the weld pool during deposition to predict convection trends based off known parameters. The predicted trends for each vendor are compared to obtain a better understanding of how mixing in the molten pool varied based on deposition parameters.

Additive Manufacturing↗

Quantifying Leak Rates and Vacuum Impact on Granular Soil Column Pneumatic Seal

Introduction: All In-Situ Resource Utilization (ISRU) regolith reactor systems that aim to extract gases or vapors must maintain a seal along the flow path of the regolith from the inlet hopper to where extracted gases or vapors are generated. Using a column of regolith at the inlet has been proposed as an alternative to systems using series of mechanical valves to reduce mass and complexity. Preliminary investigations into the use of regolith itself as a barrier to escaping fluids have been conducted, which this research seeks to expand upon. Different vertical regolith simulant plug diameters and heights were tested under vacuum with dry nitrogen to observe leak rates and approximate steady state pressure gradient equilibria. The effect of absolute ambient pressure was also studied by repeating tests at incremental pressures ranging from 10e-5 Torr up to 1 atm to corroborate ambient bench tests from previous work. The impact of regolith dynamics was addressed to emulate continuous-feed regolith reaction processes by measuring leak rates through continuously draining simulant columns. OB-1 lunar highlands regolith simulant was used. Experimental Setup: Regolith simulant columns were constructed using a motorized butterfly valve separating upper acrylic tubes of various diameters from a larger, lower aluminum pipe plenum (See Figure 1). The butterfly valve was cracked open such that gas could freely pass through whereas regolith would bridge and stay suspended as a cohesive plug above in the acrylic portion. These assemblies were loaded into a vacuum chamber, and a compressed dry nitrogen line was routed into the lower plenum. Pressure transducers measured the absolute pressure of the lower plenum, which was compared to absolute chamber pressure to monitor the pressure gradient across the regolith simulant. A precision mass flow controller enabled fine control of plenum pressure. The plenum pressure was quasi-statically increased to set points informed by previous work [1] before shutting off the gas supply to allow the now pressurized gas volume to dissipate through the simulant into the chamber in what was referred to as a ‘leak test.’ In a separate set of tests, the pressure was allowed to quasi-statically build until visual disturbances in the soil were observed via remote camera feeds. Initial leak rate tests were repeated on the disturbed soil to observe changes in leak rates between initial as-poured and post-bubbling soil conditions. This test schedule was applied to simulant columns of three different diameters (approximately 2, 4, and 6 inches) as well as two different heights (around 8 and 16 inches) and was repeated at least three times for each configuration. Plenums were also pressurized while the butterfly valves were opened slightly further to facilitate slow, controlled, constant flowrate draining of the simulant to observe the impact that continuous-feed systems might have on leak rate through the regolith column. Finally, the regolith columns were reset and refilled and static leak rate tests were repeated approximately every decade in ambient pressure from 10e-5 Torr up to ambient conditions to observe the impact of transitioning out of the molecular flow regime on regolith column seal performance. Analysis: Trends with respect to each factor examined were constructed and reported alongside commentary on the most likely responsible physical phenomena. Confidence intervals were applied to the predictive findings of the data by quantifying instrumentation error and accounting for probabilistic impacts of curve fitting the generic exponential function onto pressure decay data. The leak rate is assumed directly proportional to the exponential term coefficient, while the steady-state static pressure gradient is assumed to be the asymptotic limit of each curve fit. Finally, a generic model developed by Ogino et al. [2] and later tuned for ambient testing was re-tuned using vacuum data to estimate leak rates on the lunar surface. Results: Pressure decay rates decreased up to 27-fold under vacuum when compared to ambient pressure, as shown in Figure 2. Simulant column diameter was found to exponentially correlate to both regolith column leak rate and steady-state equilibrium pressure delta. Height was found to be linearly related to leak rate –doubling the regolith column height consistently yielded around a 40% reduction in leak rate. Dynamic regolith draining was found to slightly increase leak rate, although to a lower degree than expected. Conclusion: Because of the transition into molecular gas flow through regolith, vacuum testing is required to properly estimate performance. Regolith columns should be designed as narrow as possible to avoid bridging and minimize leak rate, and as tall as is feasible to both further reduce leak rate and to help prevent bubbling, which has been shown to degrade the quality of the regolith seal. Disrupted soil was found to exhibit higher leak rates than pristine, as-poured soil, indicating the need to consider the impact of recovery from depressurization events. Regolith columns appear equally feasible for both dosed and continuous-feed reactors.

ISRU↗

Evaluation of Limiter Functions for Supersonic Applications

Limiters commonly used in the simulations of flows with discontinuities are compared with the new limiter function proposed by Nishikawa using idealized test cases in two dimensions as well as complex three-dimensional problems. The Nishikawa limiter is observed to be consistently the least dissipative in idealized test cases as well as complex practical problems, for both steady and time-dependent problems. In the case of steady simulations, its iterative convergence characteristics are either similar or better than other limiter functions. The nearfield sonic boom signature of a low-boom demonstrator is computed to demonstrate the utility of the Nishikawa limiter function for realistic supersonic aircraft configurations.

Computational fluid dynamics↗

The origin of the moon.

Recent studies of the dynamics and thermodynamics of the early solar nebula have provided a basis for a theory of the origin of the moon which is consistent with presently developed views of the origin of the solar system. The hypothesis of inhomogeneous planetary accretion is extended. It is suggested that the anomalous properties of the moon, such as its enrichment in Ca, Al, Ti, and other refractories and its depletion in iron and volatiles can be explained if the bulk of the moon represents a high temperature condensate. It is proposed that the moon is composed chiefly of compounds that condense before iron and that the volatile content of the moon was brought in as a thin veneer after the solar nebula dissipated.

Anderson, D. L.↗

Further numerical calculations of the circulation of the atmosphere of Venus

Results of two-dimensional simulations of the deep circulation of Venus are presented. They prove that the high surface temperature can only be explained by the greenhouse effect, and that Goody and Robinson's (1966) dynamical model is not valid. Analytical radiative equilibrium solutions for a semi-grey atmosphere, both with and without an internal heat source, are presented. It is shown that the greenhouse effect is sufficient to produce the high surface temperature if the thermal optical depth is much greater than 100 and the ratio of solar to thermal optical depths does not exceed 0.005. A two-dimensional version of a three-dimensional model is used to test the validity of the new mechanism proposed by Gierasch (1975) to explain the 4-day circulation. Numerical experiments with large horizontal viscosities failed to show strong zonal velocities even for the case of large Prandtl numbers. It is observed that the dissipation of angular momentum introduced by the strong horizontal diffusion more than compensates for the upward transport of angular momentum due to the Hadley cell.

De Rivas, E. K.↗

Assessment of numerical techniques for unsteady flow calculations

The characteristics of unsteady flow motions have long been a serious concern in the study of various fluid dynamic and combustion problems. With the advancement of computer resources, numerical approaches to these problems appear to be feasible. The objective of this paper is to assess the accuracy of several numerical schemes for unsteady flow calculations. In the present study, Fourier error analysis is performed for various numerical schemes based on a two-dimensional wave equation. Four methods sieved from the error analysis are then adopted for further assessment. Model problems include unsteady quasi-one-dimensional inviscid flows, two-dimensional wave propagations, and unsteady two-dimensional inviscid flows. According to the comparison between numerical and exact solutions, although second-order upwind scheme captures the unsteady flow and wave motions quite well, it is relatively more dissipative than sixth-order central difference scheme. Among various numerical approaches tested in this paper, the best performed one is Runge-Kutta method for time integration and six-order central difference for spatial discretization.

Hsieh, Kwang-Chung↗

Carbon nanotube vacuum gauges with wide-dynamic range and processes thereof

A miniature thermal conductivity gauge employs a carbon single-walled-nanotube. The gauge operates on the principle of thermal exchange between the voltage-biased nanotube and the surrounding gas at low levels of power and low temperatures to measure vacuum across a wide dynamic range. The gauge includes two terminals, a source of constant voltage to the terminals, a single-walled carbon nanotube between the terminals, a calibration of measured conductance of the nanotube to magnitudes of surrounding vacuum and a current meter in electrical communication with the source of constant voltage. Employment of the nanotube for measuring vacuum includes calibrating the electrical conductance of the nanotube to magnitudes of vacuum, exposing the nanotube to a vacuum, applying a constant voltage across the nanotube, measuring the electrical conductance of the nanotube in the vacuum with the constant voltage applied and converting the measured electrical conductance to the corresponding calibrated magnitude of vacuum using the calibration. The nanotube may be suspended to minimize heat dissipation through the substrate, increasing sensitivity at even tower pressures.

Manohara, Harish↗

Dynamical evolution of a solar coronal magnetic field arcade

Calculations of the long-term dynamical evolution of a solar coronal magnetic field arcade which is subjected to shearing photospheric flows are presented. The evolution is obtained by numerical solution of a subset of the resistive magnetohydrodynamic equations. For a simplified model of the bipolar magnetic field observed in the solar corona, it is found that photospheric flow produces a slow evolution of the magnetic field, with a buildup of magnetic energy. For certain photospheric shear profiles, the field configuration produced is linearly unstable to an ideal magnetohydrodynamic mode when the shear exceeds a critical value. The nonlinear evolution of this instability shows the spontaneous formation of current sheets. Reconnection of the magnetic field produces a rapid release of magnetic energy. The major fraction of the energy is dissipated resistively, while a small fraction is converted into kinetic energy of an ejected plasmoid. The relevance of these results to two-ribbon flares is discussed.

Mikic, Z.↗