Search NASA⌕ Search

SEARCH · Search NASA

Results for “fractional mechanical dissipation”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

A Themodynamically Consistent Fractional Visco-Elasto-Plastic Model with Memory-Dependent Damage for Anomalous Materials

We develop a thermodynamically consistent, fractional visco-elastoplastic model coupled with damage for anomalous materials. The model utilizes Scott-Blair rheological elements for both visco- elastic/plastic parts. The constitutive equations are obtained through Helmholtz free-energy potentials for Scott-Blair elements, together with a memory-dependent fractional yield function and dissipation inequalities. A memory-dependent Lemaitre-type damage is introduced through fractional damage energy release rates. For time-fractional integration of the resulting nonlinear system of equations, we develop a first-order semi-implicit fractional return-mapping algorithm. We also develop a finite-difference discretization for the fractional damage energy release rate, which results into Hankel-type matrix-vector operations for each time-step, allowing us to reduce the computational complexity from O(N 3 ) to O(N 2 ) through the use of Fast Fourier Transforms. Our numerical results demonstrate that the fractional orders for visco-elasto-plasticity play a crucial role in damage evolution, due to the competition between the anomalous plastic slip and bulk damage energy release rates.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Adiabatic Heating in Modular Plasticity Models

Plastic deformations in metals are dissipative. Some fraction of the dissipated mechanical energy (plastic work) is converted into thermal energy and serves as a heat source. In cases where the heat cannot be readily transferred to the environment, the local temperature will increase thereby producing variations in mechanical behaviors associated with temperature-dependent properties (e.g. thermal softening due to decreasing yield strengths). This issue is often referred to as "adiabatic heating as an adiabatic temperature condition corresponds to the limiting case where no heat transfer takes place. The impact of converting plastic work into heat on the mechanical response of metals has been long studied. Nonetheless, it still remains an issue. For instance, with respect to ductile failure, the second Sandia Fracture Challenge noted that accounting for plastic heat generation was necessary for predictions under dynamic loading conditions. Furthermore, both experimental and modeling efforts continue to be pursued to better describe and understand the effect of plastic work conversion into heat on structural responses. Noting the need for capturing plastic work conversion into heat in structural analyses, a simple and fairly traditional representation of these responses has been added into existing modular plasticity models in the Library of Advanced Materials for Engineering (LAME). Here, these capabilities are briefly described with the underlying theory and numerical implementation discussed in Sections 2 and 3, respectively. Examples of syntax are given in Section 4 and some verification exercises are found in Section 5. Simple structural analyses are presented in Section 6 to briefly highlight the impact of these features and concluding thoughts are given.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Fast Dissipation of Colliding Alfvén Waves in a Magnetically Dominated Plasma

Magnetic energy around compact objects often dominates over plasma rest mass, and its dissipation can power the object's luminosity. We describe a dissipation mechanism that works faster than magnetic reconnection. The mechanism involves two strong Alfvén waves with anti-aligned magnetic fields B 1 and B 2 that propagate in opposite directions along the background magnetic field B 0 and collide. The collision forms a thin current sheet perpendicular to B 0 , which absorbs the incoming waves. The current sheet is sustained by an electric field E breaking the magnetohydrodynamic condition E < B and accelerating particles to high energies. We demonstrate this mechanism with kinetic plasma simulations using a simple setup of two symmetric plane waves with amplitude A = B 1 / B 0 = B 2 /B 0 propagating in a uniform B 0 . The mechanism is activated when A > 1/2. It dissipates a large fraction of the wave energy, † = (2A - 1)/A 2 , reaching 100% when A = 1. The plane geometry allows one to see the dissipation process in a one-dimensional simulation. We also perform two-dimensional simulations, enabling spontaneous breaking of the plane symmetry by the tearing instability of the current sheet. At moderate A of main interest, the tearing instability is suppressed. Dissipation transitions to normal, slower, magnetic reconnection at A $\gg$ 1. The fast dissipation described in this paper may occur in various objects with perturbed magnetic fields, including magnetars, jets from accreting black holes, and pulsar wind nebulae.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1985

Ion reflection, gyration, and dissipation at supercritical shocks

Topics on collisionless plasmas are discussed and include: (1) Ion reflection is the dominant ion dissipation mechanism at nearly perpendicular, supercritical shocks. (2) An increasing fraction of the ions incident on a supercritical shock is reflected as the Mach number increases. The actual fraction reflected can be predicted using the Rankine Hugoniot conservation relations. (3) The effective temperature associated with the dispersion in velocity space associated with ion reflection accounts for a large fraction of the temperature rise observed across supercritical, quasi perpendicular shocks.

Gosling, J. T.↗
DOE OSTI2019

New Method for Describing Damped Vibrations of a Beam with a Built-in End

We propose a method for describing damped vibrations of a beam with a built-in end considering the dynamic hysteresis that determines mechanical energy dissipation due to viscoelasticity. As the mathematical basis, we have used the fractional integro-differentiation apparatus. Rapidly damped vibrations of a foamed polypropylene beam have been studied experimentally. It is shown that the theoretical model successfully describes experimental data.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in planar geometry are derived in the small Atwood number (Boussinesq) limit. The models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar (heavy-fluid mass fraction) variance S and its dissipation rate Χ modeled either differentially or algebraically. The models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing, respectively. Mixing layer growth parameters and other physical observables relevant to each instability are obtained explicitly as functions of the model coefficients. The turbulent fields are also expressed in terms of the model coefficients, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. The model calibration methodology is described and discussed. Expressions for a subset of the various physical observables are used to calibrate each of the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to these values of the observables and to specific canonical Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz turbulent flows. A calibrated four-equation model is then used to reconstruct the mean and turbulent fields, and late-time turbulent equation budgets for each instability-induced flow across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗
DOE OSTI2024

Intermittent cluster dynamics and temporal fractional diffusion in a bulk metallic glass

Glassy solids evolve towards lower-energy structural states by physical aging. This can be characterized by structural relaxation times, the assessment of which is essential for understanding the glass’ time-dependent property changes. Conducted over short times, a continuous increase of relaxation times with time is seen, suggesting a time-dependent dissipative transport mechanism. By focusing on micro-structural rearrangements at the atomic-scale, we demonstrate the emergence of sub-diffusive anomalous transport and therefore temporal fractional diffusion in a metallic glass, which we track via coherent x-ray scattering conducted over more than 300,000 s. At the longest probed decorrelation times, a transition from classical stretched exponential to a power-law behavior occurs, which in concert with atomistic simulations reveals collective and intermittent atomic motion. Our observations give a physical basis for classical stretched exponential relaxation behavior, uncover a new power-law governed collective transport regime for metallic glasses at long and practically relevant time-scales, and demonstrate a rich and highly non-monotonous aging response in a glassy solid, thereby challenging the common framework of homogeneous aging and atomic scale diffusion.

42 ENGINEERING↗
DOE OSTI2026

A physically based mechanical model for Mullins effect in thermoplastic polyurethanes

Despite decades of research, connecting the chemical and physical structure of thermoplastic polyurethanes to their mechanical properties remains highly challenging. Of particular note are their large-deformation and rate-dependent behaviors, which vary greatly with molecular chemistry, including the type and relative content of soft and hard segments. In this work, we develop a physically motivated mechanical theory for predicting the behavior of thermoplastic polyurethanes. The theory incorporates a representation of microstructural evolution during mechanical deformation, which captures the signatures of stress softening over cyclic loading (commonly referred to as the Mullins effect). There are only eight physically motivated fitting parameters, including a direct dependence on the hard segment fraction. The model predicts that increasing the hard segment fraction leads to higher stiffness and greater energy dissipation, in quantitative agreement with published experimental data. Furthermore, we provide a comprehensive analysis of the model and validate its predictions across several independent datasets focused on mechanical characterization. Direct comparisons to experimental data demonstrate its predictive capability on the effect of loading rate, cyclic deformations, and applied tension or compression. Altogether, this work establishes a predictive framework that connects polymer chemistry and microstructure to emergent mechanical behaviors.

36 MATERIALS SCIENCE↗
NASA NTRS1978

Reverse Current in Solar Flares

An idealized steady state model of a stream of energetic electrons neutralized by a reverse current in the pre-flare solar plasma was developed. These calculations indicate that, in some cases, a significant fraction of the beam energy may be dissipated by the reverse current. Joule heating by the reverse current is a more effective mechanism for heating the plasma than collisional losses from the energetic electrons because the Ohmic losses are caused by thermal electrons in the reverse current which have much shorter mean free paths than the energetic electrons. The heating due to reverse currents is calculated for two injected energetic electron fluxes. For the smaller injected flux, the temperature of the coronal plasma is raised by about a factor of two. The larger flux causes the reverse current drift velocity to exceed the critical velocity for the onset of ion cyclotron turbulence, producing anomalous resistivity and an order of magnitude increase in the temperature. The heating is so rapid that the lack of ionization equilibrium may produce a soft X-ray and EUV pulse from the corona.

Knight, J. W., III↗
NASA NTRS1994

Unsteady outflow models for cosmological gamma-ray bursts

The 'event' that triggers a gamma-ray burst cannot last for more than a few seconds. This is, however, long compared with the dynamical timescale of a compact stellar-mass object (approximately 10 (exp-3) s). Energy is assumed to be released as an outflow with high mean Lorentz factor Gamma. But a compact stellar-mass collapse or merger is, realistically, likely to generate a mass (or energy) flux that is unsteady on some timescales in the range 10(exp -3) - 10 s. If Gamma fluctuates by a factor of approximately 2 around its mean value, relative motions within the outflowing material will themselves (in the comoving frame) be relativistic, and can give rise to internal shocks. For Gamma approximately 10(exp 2), the resultant dissipation occurs outside the 'photosphere' and can convert a substantial fraction of the overall outflow energy into nonthermal radiation. This suggests a mechanism for cosmological bursts that demands less extreme assumptions (in respect of Gamma-values, freedom from baryonic contamination, etc.) than earlier proposals.

Rees, M. J.↗
DOE OSTI2023

Thermomechanical Modeling of Copper

The plastic deformation of metals is a dissipative process. Some fraction of the plastic work is converted to heat which, given the temperature dependent response of metals, produces a thermal-mechanical coupling. In various cases, for instance when the loading is dynamic, this interaction can impact the resulting response of a material and/or system. Thus, appropriately capturing the heat generation from plastic work is necessary for various solid mechanics analysis. Determination of the fraction of work converted to heat has been long studied. Recent developments have demonstrated that the fraction is not constant but depends on various state variables. Resolving these features requires combined modeling and experimental studies. To this end, 304L stainless steel – a poor thermal conductor – was recently subjected to such an investigation. Advanced modeling capabilities were deployed to assess novel thermomechanically coupled experiments. As a complement to that study, in the current work a similar investigation is performed on copper – a good thermal conductor – to assess performance on the opposite end of the spectrum. The current document discusses these modeling efforts.

36 MATERIALS SCIENCE↗
DOE OSTI2023

The Seismic Signature of a High-Energy Density Physics Laboratory and Its Potential for Measuring Time-Dependent Velocity Structure

The Z Machine at Sandia National Laboratories is a pulsed power facility for high-energy density physics experiments that can shock materials to extreme temperatures and pressures through a focused energy release of up to ~ 25 MJ in < 100 nanoseconds. It has been in operation for more than two decades and conducts up to ~ 100 experiments, or “shots,” per year. Based on a set of 74 known shot times from 2018, we determined that Z Machine shots produce detectable ~ 3–17 Hz ground motion 12 km away at the Albuquerque Seismological Laboratory, New Mexico (ANMO), borehole seismograph, with peak signal at ~ 7 Hz. The known shot waveforms were used to create a three-component template, leading to the detection of 2339 Z Machine shots since 1998 through single-station cross-correlation. Local seismic magnitude estimates range from local magnitude (M L ) —2 to —1.3 and indicate that only a small fraction of the shot energy is transmitted by seismic phases observable at 12 km distance. The most recent major facility renovation, which was intended to decrease mechanical dissipation, is associated with an abrupt decrease in observed seismic amplitudes at ANMO despite stable maximum shot energy. The highly repetitive impulsive sources are well suited to coda-wave interferometry to investigate time-dependent velocity structures. Relative velocity variations (dv/v) show an annual cycle with amplitude of ~ 0.2%. Local minima are observed in the late spring, and dv/v increases through the summer monsoon rainfall, possibly reflecting patchy saturation as rainfall infiltrates near the eastern edge of the Albuquerque basin. Importantly, the cumulative results demonstrate that forensic seismology can provide insight into long-term operation of facilities such as pulsed-power laboratories, and that their recurring signals may be valuable for studies of time-dependent structure

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1986

Magnetic energy dissipation in force-free jets

It is shown that a magnetic pressure-dominated, supersonic jet which expands or contracts in response to variations in the confining external pressure can dissipate magnetic energy through field-line reconnection as it relaxes to a minimum-energy configuration. In order for a continuous dissipation to occur, the effective reconnection time must be a fraction of the expansion time. The dissipation rate for the axisymmetric minimum-energy field configuration is analytically derived. The results indicate that the field relaxation process could be a viable mechanism for powering the synchrotron emission in extragalactic jets if the reconnection time is substantially shorter than the nominal resistive tearing time in the jet.

Choudhuri, Arnab Rai↗
DOE OSTI2021

Toroidal modeling of runaway avalanche in DIII-D discharges

Abstract A toroidal modeling tool is developed to study the runaway electron (RE) avalanche production process in tokamak plasmas, by coupling the Rosenbluth–Putvinski avalanche model (Rosenbluth and Putvinski 1997 Nucl. Fusion 37 1355) with an n = 0 magneto-hydrodynamic (MHD) solver. Initial value numerical simulations are carried out for two DIII-D discharges with different plasma shapes (one near circular, and the other with high elongation). It is found that, assuming the same level of about 1% seed current level, the Rosenbluth–Putvinski model somewhat underestimates the RE plateau current for the circular-shaped plasma, as compared with that measured in DIII-D experiments. For an elongated, higher current plasma, simulations find strong runaway current avalanche production despite the lack of measured plateau RE current in experiments. A possible reason for this discrepancy is a lack of additional RE dissipation physics in the present two-dimensional model. Systematic scans of the plasma boundary shape, at fixed pre-disruption plasma current, find that the plasma elongation helps to reduce the RE avalanche production, confirming recent results obtained with an analytic model (Fülöp et al 2020 J. Plasma Phys. 86 474860101). The effect of the plasma triangularity (either positive or negative), on the other hand, has a minor effect. On the physics side, the avalanche process involves two competing mechanisms associated with the electric field. On the one hand, a stronger electric field produces a higher instantaneous avalanche growth rate. On the other hand, a fast growing RE current quickly reduces the fraction of the conduction current together with the electric field, and hence a faster dissipation of the poloidal flux. As a final result of these two dynamic processes, the runaway plateau current is not always the largest with the strongest initial electric field. These results lay the foundation for future self-consistent inclusion of the MHD dynamics and the RE amplification process.

Physics↗
NASA NTRS1972

Free oscillations of the sun and their possible stimulation by solar flares.

The main characteristics of free oscillations of the sun are described, and the sources and sinks of energy are estimated. Adiabatic oscillations and damping are considered. It is shown how a large solar flare can exert a significant mechanical impulse on the sun by causing a wave of compression to move subsonically into the interior. The solar interior does not easily dissipate low-amplitude acoustical energy, and therefore a relatively large fraction of the energy should remain available to go into the normal modes. An estimate is made of how this energy might be distributed among the more interesting modes.

Wolff, C. L.↗
DOE OSTI2022

Spatio-Temporal Clustering of Seismicity Enabled by Off-Fault Plasticity

While significant progress has been made in understanding earthquake source processes in linear elastic domains, the effect of more realistic rheologies including plasticity is poorly understood. Here, in this work, we simulate the sequence of earthquake and aseismic slip of a 2D antiplane rate-and-state fault embedded in a full-space elastic-plastic bulk. We show that off-fault plasticity may lead to partial ruptures as well as temporal clustering of seismic events. Furthermore, the interaction of fault slip and off-fault plasticity results in pockets of slip deficit. While the energy dissipated through plastic deformation remains a small fraction of the total energy budget, its impact on the source characteristics is disproportionally large through the redistribution of stresses and viscous relaxation. Our results suggest a new mechanism of dynamic heterogeneity in earthquake physics that may have important implications on earthquake size distribution and energy budget.

58 GEOSCIENCES↗
NASA NTRS1995

Electron energy transport in the solar wind: Ulysses observations

The electron heat flux in the solar wind has been measured by the Ulysses solar wind plasma experiment in the ecliptic from 1 to 5 AU and out of the ecliptic during the recently completed pass over the solar south pole and the ongoing pass over the solar north pole. Although the electron heat flux contains only a fraction of the kinetic energy of the solar wind. the available energy is sufficient to account for the non-adiabatic expansion of the solar wind electrons. The Ulysses measurements indicate that the electron heat flux is actively dissipated in the solar wind. The exact mechanism or mechanisms is unknown. but a model based on the whistler heat flux instability predicts radial gradients for the electron heat flux in good agreement with the data. We will present measurements of the correlation between wave activity measured by the unified radio and plasma experiment (URAP) and the electron heat flux throughout the Ulysses mission. The goal is to determine if whistler waves are a good candidate for the observed electron heat flux dissipation. The latitudinal gradients of the electron heat flux. wave activity. and electron pressure will be discussed in light of the changes in the magnetic field geometry from equator to poles.

Scime, Earl↗
DOE OSTI2023

The Dynamics of Megafire Smoke Plumes in Climate Models: Why a Converged Solution Matters for Physical Interpretations

As the climate system warms, megafires have become more frequent with devastating effects. A byproduct of these events is the creation of smoke plumes that can rise into the stratosphere and spread across the globe where they reside for many months. To gain a deeper understanding of the plume dynamics, global climate simulations of a megafire were performed at a wide range of grid spacings from 2.0° down to 7 km, including a 7 km nonhydrostatic experiment. The analysis focuses on how the resolved dynamics affects the specification of the plume characteristics such as injection height and black carbon (BC) mass. Prior studies initialize the smoke plume at one or a few grid points and this is shown here to produce severely dissipative dynamics. In order to validate such simulations with observations, enhancements of the plume characteristics to offset the dissipation is necessary. Using a numerically converged simulation, sensitivity tests show that to approximate the observed stratospheric lifetime, a reduction in BC fraction by 50% is necessary for external mixtures. The vorticity dynamics of the plume is also analyzed with a Lagrangian budget to understand the mechanisms responsible for the evolution of a collocated anticyclonic vortex. The results can be distilled down into a simple conceptual model. As the plume rises, the air diverges at the top of the updraft where the largest concentrations of smoke are found. This divergence induces a dilution of the background cyclonic absolute vorticity producing an anticyclonic vortex. Vortex decay occurs from opposite arguments.

54 ENVIRONMENTAL SCIENCES↗