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 1,099 records · Page 61

Calculation of collisionless pitch-angle scattering of runaway electrons with synchrotron radiation via high-order guiding-centre equation

Recently, the collisionless pitch-angle scattering for relativistic runaway electrons (REs) in toroidal geometries such as tokamaks was discovered through a full orbit simulation approach (Liu et al., Nucl. Fusion, vol. 56, 2016, p. 064002), and it was then theoretically investigated that a new expression for the magnetic moment, including the second-order corrections, could essentially reproduce the so-called collisionless pitch-angle scattering process (Liu et al., Nucl. Fusion, vol. 58, 2018, p. 106018). In this paper, with synchrotron radiation, extensive numerical verification of the validity of the high-order guiding-centre theory is given for simulations involving REs by incorporating such an expression for the magnetic moment into our particle tracing code. A high-order guiding-centre simulation approach with synchrotron radiation (HGSA) is applied. Synchrotron radiation plays an essential role in the life cycle of REs. The energy of REs first increases and then becomes saturated until the electric field acceleration is balanced by the radiation dissipation. Unfortunately, the process cannot be simulated accurately with the standard guiding-centre model, i.e. the first-order guiding-centre model. Remarkably, it is found that the HGSA can effectively produce the fundamental process of REs. Since the time scale of the energy saturation of REs is close to seconds, the computational cost becomes significant. In order to save costs, it is necessary to estimate the time of energy saturation. An analytical estimate is derived for the time it takes for synchrotron drag to balance an accelerating electric field and the provided formula has been numerically verified. Finally, test calculations reveal that HGSA is favourable for exploiting the dynamics of REs in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Asynchronous Ballistic Reversible Computing using Superconducting elements

Computing uses energy. At the bare minimum, erasing information in a computer increases the entropy. Landauer has calculated %7E k B T log(2) Joules is dissipated per bit of energy erased. While the success of Moores law has allowed increasing computing power and efficiency for many years, these improvements are coming to an end. This project asks if there is a way to continue those gains by circumventing Landauer through reversible computing. We explore a new reversible computing paradigm, asynchronous ballistic reversible computing or ABRC. The ballistic nature of data in ABRC matches well with superconductivity which provides a low-loss environment and a quantized bit encoding the fluxon. We discuss both these and our development of a superconducting fabrication process at Sandia. We describe a fully reversible 1-bit memory cell based on fluxon dynamics. Building on this model, we propose several other gates which may also offer reversible operation.

97 MATHEMATICS AND COMPUTING↗

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↗

A new effective theory for stochastic relativistic hydrodynamics

Thermal fluctuations are a fundamental feature of dissipative systems that are essential for understanding physics near the expected critical point of QCD and in small systems. When such fluctuations are modeled naively in relativistic systems, strange features can appear such as negative self-correlation functions. We construct an effective theory for nonlinear stochastic relativistic hydrodynamics that ensure a well-posed mathematical formulation. Using Crooks fluctuation theorem, we derive a symmetry of the effective action that incorporates fluctuations through a suitable free energy functional. For divergence type theories, the action can then be fully specified using a single vector generating current. The equations of motion obtained using this procedure are guaranteed to be flux conservative and symmetric hyperbolic when the dynamics is causal. This ensures that these equations are well-posed (for suitable initial data) and are in a form that can easily be simulated, including with Metropolis techniques.

Mullins, Nicki [University of Illinois at Urbana-C↗

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

Constraints on dark matter from dynamical heating of stars in ultrafaint dwarfs. II. Substructure and the primordial power spectrum

There is a large and growing interest in observations of small-scale structure in dark matter. We propose a new way to probe dark matter structures in the ∼ 10 – 10 8 M ⊙ range. This allows us to constrain the primordial power spectrum over shorter distances scales than possible with direct observations from the CMB. For k in the range ∼ 10 – 1000 Mpc − 1 our constraints on the power spectrum are orders of magnitude stronger than previous bounds. We also set some of the strongest constraints on dark matter isocurvature perturbations. Our method relies on the heating effect such dark matter substructures would have on the distribution of stars in an ultrafaint dwarf galaxy. Many models of inflation produce enhanced power at these short distance scales and can thus be constrained by our observation. Further, many dark matter models such as axion dark matter, self-interacting dark matter and dissipative dark matter, produce dense structures which could be constrained this way. Published by the American Physical Society 2024

Graham, Peter W. (ORCID:0000000216001601)↗

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

Dynamic Molecular Investigation of the Solid-Electrolyte Interphase of an Anode-Free Lithium Metal Battery Using In Situ Liquid SIMS and Cryo-TEM

Solid electrolyte interphase (SEI) has been widely perceived to play a critical role in the stable cycling of rechargeable batteries. However, associated with the fragile and air-sensitive nature of the SEI layer, delineation of the formation process and the nature of SEI remains a big challenge. Here, we use in situ liquid time-of-flight secondary ion mass spectroscopy (TOF-SIMS), cryo- transmission electron microscope (TEM) and density functional theory (DFT) calculation to delineate molecular process on the formation of SEI layer under the dynamic operating condition. We discover that the onset potential for SEI layer formation and the thickness of the SEI show dependence on the solvation shell structure. Using LiCoO 2 as a cathode and Cu film as an anode, the SEI is noticed to start to form at around 2.0 V and reach its final thickness (irreversible part, ~ 40-50 nm) at about 3.0 V in the 1 M LiPF 6 –EC/DMC electrolyte, while for the case of 1 M LiFSI–DME, the SEI starts to form at around 1.5 V and reaches its final thickness (~ 20 nm) at about 2.0 V. The in situ TOF-SIMS clearly indicates the outer SEI layer formation and dissipation upon charging and discharging, implying a continued evolution of electrolyte structure with extended cycling. In conclusion, the present work establishes a direct correlation between the molecular signature of SEI layer with solvation feature of electrolytes in lithium batteries, providing insights for tailoring SEI layer toward improved electrochemical properties of lithium batteries.

25 ENERGY STORAGE↗

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

Measurements of the imploding plasma sheath in triple-nozzle gas-puff z pinches

Gas-puff z-pinch implosions are characterized by the formation of a dense annular plasma shell, the sheath, that is driven to the axis by magnetic forces and therefore subject to the magneto-Rayleigh–Taylor instability. Here, the conditions within these sheaths are measured on the 1-MA COBRA generator at Cornell University for various gas species and initial fill densities. The gas-puff loads are initialized by a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet. Thomson scattering and laser interferometry provide spatially resolved flow, temperature, and electron density profiles midway through the implosion, while extreme ultraviolet pinhole cameras record the evolution of the plasma column and photoconducting diodes measure x-ray emission. Analysis of the scattering spectra includes a means of discriminating between thermal and non-thermal broadening to test for the presence of hydrodynamic turbulence. Two types of sheath profiles are observed, those with sharp discontinuities at the leading edge and those with smooth gradients. In both cases, non-thermal broadening is generally peaked at the front of the sheath and exhibits a characteristic decay length that roughly scales with the sheath ion mean free path. Here, we demonstrate that this non-thermal broadening term is inconsistent with laminar velocity gradients and is more consistent with dissipative turbulence driven by unstable plasma waves in a collisionless shock. The resulting differences in sheath profile are then set by the sheath ion collisionality in a manner consistent with recent 1D kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗