Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “compression rate”

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 181 records · Page 10

Implicit Large-Eddy Simulations of Compressible Mixing Layers

Implicit large-eddy simulations of the self-similar regions of two compressible mixing layers at high Reynolds number and convective Mach numbers of 0.381 and 0.690 were carried out. Experimental data was used as the input into a synthetic eddy method turbulent inflow boundary condition to initiate simulations. Computational grids with nearly isotropic cells were used and the grid counts ranged from 24 to 270 million points. Agreement with experimental velocity and Reynolds stress data was good given the simplifications of the simulation. The decrease in mixing layer growth rate with increasing compressibility compares well with the literature. Budgets for the Reynolds stress transport equation were extracted from the simulations. It was demonstrated that the grids used did not resolve the dissipation term at the experimental Reynolds number. However, low Reynolds number simulations indicate that the numerical dissipation replaces the physical dissipation and the dissipation term in the transport equation can be represented by the summation term that quantifies the imbalance in the budget. The budgets for the high Reynolds number mixing layers agree well with previous low Reynolds number studies; indicating that the primary mechanism for the reduced growth rate and increased anistropy of the Reynolds stress tensor due to compressibility is a reduction in the pressure-strain and production terms. A new scaling based on the magnitude of the Reynolds stress tensor is proposed. This scaling provides better relative comparisons between the data and shows that the decrease in pressure-strain in the shear component is the primary driver in reducing the transverse normal stress.

turbulence↗

Shock-induced chemistry and high strain-rate viscoelastic behavior of a phenolic polymer

We use impact experiments and a finite element model (up to 1.2 GPa), and molecular dynamics simulations (up to 60 GPa), to examine the behavior of a phenolic polymer under shock compression, spanning both nonreactive and reactive regimes. In the nonreactive regime, relaxation following compression at strain rates of ∼105 s−1 can be explained by viscoelasticity observed at ordinary laboratory rates (≲1 s−1) by accounting for the temperature dependence of the phenolic β-transition. Reasonable agreement is found between the measured shock Hugoniot up to 1.2 GPa and molecular dynamics simulation for cross-linked structures of comparable density. We also observed a first-order mechanical transition near 0.36 GPa shock stress and estimated a spall strength of 0.102 GPa and Hugoniot elastic limit of 1–2 GPa. The shock stress is found to vary up to 24% among phenolics made with different resin and/or cure processes. Finally, molecular dynamics simulations are used to identify a reactive regime at shock pressures ≳20 GPa that is characterized by chemically driven, rate-dependent relaxation processes, including dehydrogenation and dehydration reactions that promote the formation of a dense, highly cross-linked carbonaceous solid and the release of light volatiles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quasi-static to Dynamic Mechanical Response and Microstructure Development of Tantalum-Tungsten Alloys

Lawrence Livermore National Laboratory (LLNL) is interested the quasi-static to dynamic mechanical response and microstructure evolution of tantalum-tungsten (Ta-W) alloys (Ta-2.5W, Ta-5W, and Ta-10W, wt.%) made by conventional wrought processing and additive manufacturing (AM). This work scope was performed at the Colorado School of Mines (Mines) and included quasi-static (e.g., 10 -3 s -1 ) mechanical testing in tension and compression, along with selected high strain rate (Kolsky) pressure bar testing in compression (e.g., 10 3 s -1 ), with and without temperature variations in some instances. Complementary microstructure characterization was performed on undeformed and deformed samples to understand the role of processing on microstructural evolution and the deformation mechanisms that impact the mechanical response with variations in strain rate, temperature, and strain state (e.g., tension versus compression in selected examples). The wrought material provided by LLNL from Viridis Materials was found to have unrecrystallized regions within the microstructure, which led to unexpected results relative to previously reported properties for Ta-W. The AM Ta-2.5W (wt.%) material provided by LLNL was found to have higher compressive strength than wrought Ta-2.5W (wt.%), which is hypothesized to be due to differences in the crystallographic texture between the two materials. This project partially supported several postdocs and graduate students at Mines.

36 MATERIALS SCIENCE↗

Mixing and chemical reaction in sheared and nonsheared homogeneous turbulence

Direct numerical simulations were made to examine the local structure of the reaction zone for a moderately fast reaction between unmixed species in decaying, homogeneous turbulence and in a homogeneous turbulent shear flow. Pseudospectral techniques were used in domains of 64 exp 3 and higher wavenumbers. A finite-rate, single step reaction between non-premixed reactants was considered, and in one case temperature-dependent Arrhenius kinetics was assumed. Locally intense reaction rates that tend to persist throughout the simulations occur in locations where the reactant concentration gradients are large and are amplified by the local rate of strain. The reaction zones are more organized in the case of a uniform mean shear than in isotropic turbulence, and regions of intense reaction rate appear to be associated with vortex structures such as horseshoe vortices and fingers seen in mixing layers. Concentration gradients tend to align with the direction of the most compressive principal strain rate, more so in the isotropic case.

Leonard, Andy D.↗

Enhancement of Mixing and Reaction in High-speed Combustor Flowfields

Research has been underway for a number of years, both in the United States and abroad, to develop advanced aerospace propulsion systems for use late in this century and beyond. One program is now underway at the NASA Langley Research Research has been underway for a number of years, both in the United States and abroad, to develop advanced aerospace propulsion systems for use late in this century and beyond. One program is now underway at the NASA Langley Research Center to develop a hydrogen-fueled supersonic combustion ramjet (scramjet) that is capable of propelling a vehicle at hypersonic speeds in the atmosphere. A part of that research has been directed toward the optimization of the scramjet combustor and, in particular, the efficiency of fuel-air mixing and reaction taking place in the engine. In the high-speed vehicle configurations currently being considered, achieving a high combustor efficiency becomes particularly difficult. With increasing combustor Mach number, the degree of fuel-air mixing that can be achieved through natural convective and diffusive processes is reduced leading to an overall decrease in combustion efficiency and thrust. Compressible shear layers and jets provide a good model for studying the physical processes occurring in high speed mixing and combustion in a scramjet. Mixing layers are characterized by large scale eddies that form due to the high shear that is present between the fuel and the air streams. These eddies entrain the fuel and air into the mixing region between the fluids leading to increased surface area and locally steep concentration gradients. Molecular diffusion then occurs across the strained interfaces. In an early study of high-speed mixing, Brown and Roshko show that the spreading rate of a supersonic mixing layer decreases with increasing Mach number, exhibiting a factor of three decrease in spread rate as compared with an incompressible mixing layer with the same density ratio. They conclude that the reduced spread rate is primarily due to compressibility. Papamoschou and Roshko also observe that the spreading rate of compressible mixing layers is significantly reduced over that of incompressible layers. To characterize the structure of the flow quantitatively, they define a convective Mach number. The reduction in mixing layer spreading rate is shown in these experiments to correlate well with increasing convective Mach number.

Drummond, J. Philip↗

Analysis of homogeneous turbulent reacting flows

Full turbulence simulations at low Reynolds numbers were made for the single-step, irreversible, bimolecular reaction between non-premixed reactants in isochoric, decaying homogeneous turbulence. Various initial conditions for the scalar field were used in the simulations to control the initial scalar dissipation length scale, and simulations were also made for temperature-dependent reaction rates and for non-stoichiometric and unequal diffusivity conditions. Joint probability density functions (pdf's), conditional pdf's, and various statistical quantities appearing in the moment equations were computed. Preliminary analysis of the results indicates that compressive strain-rate correlates better than other dynamical quantities with local reaction rate, and the locations of peak reaction rates seem to be insensitive to the scalar field initial conditions.

Leonard, A. D.↗

Evolution of microstructure and strength of a high entropy alloy undergoing the strain-induced martensitic transformation

In a recent work, we have reported outstanding strength and work hardening exhibited by a metastable high entropy alloy (HEA), Fe 42 Mn 28 Co 10 Cr 15 Si 5 (in at. %), undergoing the strain-induced martensitic transformation from metastable gamma austenite (γ) to stable epsilon martensite (ε). However, the alloy exhibited poor ductility, which was attributed to the presence of the brittle sigma (σ) phase in its microstructure. The present work reports the evolution of microstructure, strength, and ductility of a similar HEA, Fe 38.5 Mn 20 Co 20 Cr 15 Si 5 Cu 1.5 (in at. %), designed to suppress the formation of σ phase. A cast and then rolled plate of the alloy was processed into four conditions by annealing for 10 and 30 min at 1100 °C and by friction stir processing (FSP) at tool rotation rates of 150 and 400 revolutions per minute (RPM) to facilitate detailed examinations of variable initial grain structures. Neutron diffraction and electron microscopy were employed to characterize the microstructure and texture evolution. The initial materials had variable grain size but nearly 100% γ structure. Diffusionless strain induced γ→ε phase transformation took place under compression with higher rate initially and slower rate at the later stages of deformation, independent on the initial grain size. The transformation facilitated part of plastic strain accommodation and rapid strain hardening owing to a transformation-induced dynamic Hall-Petch-type barrier effect, increase in dislocation density, and texture. Here, the peak strength of nearly 2 GPa was achieved under compression using the structure created by double pass FSP (150 RPM followed by 150 RPM). Remarkably, the tensile elongation exhibited by the alloy was nearly 20% with fracture surfaces featuring a combination of ductile dimples and cleavage.

36 MATERIALS SCIENCE↗

PDF turbulence modeling and DNS

The problem of time discontinuity (or jump condition) in the coalescence/dispersion (C/D) mixing model is addressed in probability density function (pdf). A C/D mixing model continuous in time is introduced. With the continuous mixing model, the process of chemical reaction can be fully coupled with mixing. In the case of homogeneous turbulence decay, the new model predicts a pdf very close to a Gaussian distribution, with finite higher moments also close to that of a Gaussian distribution. Results from the continuous mixing model are compared with both experimental data and numerical results from conventional C/D models. The effect of Coriolis forces on compressible homogeneous turbulence is studied using direct numerical simulation (DNS). The numerical method used in this study is an eight order compact difference scheme. Contrary to the conclusions reached by previous DNS studies on incompressible isotropic turbulence, the present results show that the Coriolis force increases the dissipation rate of turbulent kinetic energy, and that anisotropy develops as the Coriolis force increases. The Taylor-Proudman theory does apply since the derivatives in the direction of the rotation axis vanishes rapidly. A closer analysis reveals that the dissipation rate of the incompressible component of the turbulent kinetic energy indeed decreases with a higher rotation rate, consistent with incompressible flow simulations (Bardina), while the dissipation rate of the compressible part increases; the net gain is positive. Inertial waves are observed in the simulation results.

Hsu, A. T.↗

The role of breakup and coalescence in fine-scale bubble-induced turbulence. II. Kinematics

This second part of our research explores the kinematic aspect of fine-scale bubble-induced turbulence (BIT) to (i) present the effect of bubble breakup and coalescence and (ii) compare it against the universal kinematic fine-scale turbulence characteristics reported in the literature. To this end, we simulate a dilute bubbly system of 0.5% void fraction using two distinct numerical simulations. In the volume-of-fluid (VoF) simulation, bubbles undergo breakup and coalescence. In the immersed boundary method (IBM) simulation, however, they act as rigid spheres. We also perform a simulation of classical homogeneous isotropic turbulence (HIT). The first important outcome of this study is that BIT is radically different from HIT in terms of its kinematic fine-scale characteristics. In the vorticity-dominating regions, BIT exhibits a weak vortex stretching. This weak vortex stretching is due to (a) the intermediate strain-rate eigenvalues skewed weakly to positive and (b) the extensive strain-rate eigenvector aligning perpendicular to the vorticity vector. The BIT has, on average, not only a weak enstrophy production but also a weak strain production in strain-dominating regions. The weak strain production is due to (a) the presence of vortex stretching in highly strained fluid elements and (b) the absolute magnitude of compressive strain-rate eigenvalue being as close to the extensive strain-rate eigenvalue. Thus, none of the kinematic fine-scale HIT characteristics is noted for BIT. The second important conclusion is that bubble breakup and coalescence play little to no influence on the kinematics of fine-scale BIT as VoF and IBM simulations produce similar results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A finite-strain rate- and pressure-dependent constitutive framework for analyzing shock compression behavior of cemented tungsten carbides to 100 GPa

In the present study a thermodynamically-consistent finite-strain rate-and-pressure-dependent constitutive framework is implemented to analyze the shock-compression behavior of cemented tungsten carbides to 100 GPa. Central to this framework is the use of logarithmic strain with a set of invariant basis that allow the Cauchy stress tensor to be expressed as a sum of three response terms that are mutually orthogonal, thus permitting a complete separation of the deviatoric and volumetric (pressure) response. An overstress viscoplasticity model that includes strain and strain rate hardening along with thermal softening is used to represent the deviatoric response, while a complete Mie-Grüneisen equation of state (EoS) is used to obtain the pressure response. Using this formulation, the shock-induced compression behavior of cemented tungsten carbide - obtained from planar plate impact experiments using a 30 mm powder gun to peak stresses of up to ~100 GPa - is analyzed to better understand the structure of the measured shock wave profiles and the associated in-material shock quantities. Of particular interest is the evolution of material inelasticity and strength, and temperature in the tungsten carbide samples during the shock compression process.

Cemented tungsten carbide↗

Unsteady transition measurements on a pitching three-dimensional wing

Boundary layer transition measurements were made during an experimental study of the aerodynamics of a rectangular wing undergoing unsteady pitching motions. The wing was tested at chordwise Mach numbers between 0.2 and 0.6, at sweep angles of 0, 15, and 30 deg, and for steady state, sinusoidal, and constant pitch rate motions. The model was scaled to represent a full size helicopter rotor blade, with chord Reynolds numbers between 2 and 6 x 10(exp 6). Sixteen surface hot-film gages were located along three spanwise stations: 0.08, 0.27, and 0.70 chords from the wing tip. Qualitative heat transfer information was obtained to identify the unsteady motion of the point of transition to turbulence. In combination with simultaneous measurements of the unsteady surface pressure distributions, the results illustrate the effects of compressibility, sweep, pitch rate, and proximity to the wing tip on the transition and relaminarization locations.

Lorber, Peter F.↗

Temperature distortion generator for turboshaft engine testing

The procedures and unique hardware used to conduct an experimental investigation into the response of a small-turboshaft-engine compression system to various hot gas ingestion patterns are presented. The temperature distortion generator described herein uses gaseous hydrogen to create both steady-state and time-variant, or transient, temperature distortion at the engine inlet. The range of transient temperature ramps produced by the distortion generator during the engine tests was from less than 111 deg K/sec (200 deg R/sec) to above 611 deg K/sec (1100 deg R/sec); instantaneous temperatures to 422 deg K (760 deg R) above ambient were generated. The distortion generator was used to document the maximum inlet temperatures and temperature rise rates that the compression system could tolerate before the onset of stall for various circumferential distortions as well as the compressor system response during stall.

Klann, G. A.↗

Temperature distortion generator for turboshaft engine testing

The procedures and unique hardware used to conduct an experimental investigation into the response of a small-turboshaft-engine compression system to various hot gas ingestion patterns are presented. The temperature distortion generator described herein uses gaseous hydrogen to create both steady-state and time-variant, or transient, temperature distortion at the engine inlet. The range of transient temperature ramps produced by the distortion generator during the engine tests was from less than 111 deg K/sec (200 deg R/sec) to above 611 deg K/sec (1100 deg R/sec); instantaneous temperatures to 422 deg K (760 deg R) above ambient were generated. The distortion generator was used to document the maximum inlet temperatures and temperature rise rates that the compression system could tolerate before the onset of stall for various circumferential distortions as well as the compressor system response during stall.

Klann, G. A.↗

Vector-based feedback of continuous wave radiofrequency compression cavity for ultrafast electron diffraction

The temporal resolution of ultrafast electron diffraction at weakly relativistic beam energies (⁠≲100 keV) suffers from space-charge induced electron pulse broadening. We describe the implementation of a radio frequency (RF) cavity operating in the continuous wave regime to compress high repetition rate electron bunches from a 40.4 kV DC photoinjector for ultrafast electron diffraction applications. Active stabilization of the RF amplitude and phase through a feedback loop based on the demodulated in-phase and quadrature components of the RF signal is demonstrated. This scheme yields 144 ± 19 fs RMS temporal resolution in pump–probe studies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quasi-static to Dynamic Response and Microstructure Development of Tantalum-Tungsten Alloys

Lawrence Livermore National Laboratory (LLNL) has a need to understand the quasi-static to dynamic mechanical response and microstructure evolution of tantalum-tungsten (Ta-W) alloys made by conventional wrought processing and additive manufacturing (AM). This two-year work scope will be performed at the Colorado School of Mines (Mines) and will include quasi-static (e.g., 10 -3 s -1 ) mechanical testing in tension and compression, along with high strain rate (Kolsky) pressure bar testing in tension and compression (e.g., 10 3 s -1 ), without and with temperature variations. Gas gun experiments at higher strain rates (e.g., over 10 3 s -1 ) will also be pursued. Complementary microstructure characterization will be performed on undeformed and deformed samples to fundamentally understand the role of processing on the microstructure characteristics and resulting deformation mechanisms that dictate the mechanical response of these alloys with variations in strain rate, temperature, and strain state. Alloys of interest will be Ta-based, containing up to 10 wt.% W. LLNL will supply initial wrought-processed material for sample machining. LLNL will also supply samples built by AM during the project for testing. Arc melting and casting may also be performed at Mines to produce small quantities of different alloys for machining and exploratory testing on an as-needed basis. This project will support multiple postdocs and graduate students at Mines, training them with the knowledge, skills, and abilities needed to study materials under extreme conditions of interest to the National Nuclear Security Agency.

36 MATERIALS SCIENCE↗

Strain-rate/temperature behavior of high density polyethylene in compression

The compressive strain rate/temperature behavior of highly linear, high density polyethylene was analyzed in terms of the predictive relations developed for metals and other crystalline materials. For strains of 5 percent and above, the relationship between applied strain rate, dotted epsilon, and resulting flow stress, sigma, was found to be: dotted epsilon exp times (Q sub f/RT) = k'(sigma/sigma sub c) to the nth power; the left-hand side is the activation-energy-compensated strain rate, where Q sub f is activation energy for flow, R is gas constant, and T is temperature; k is a constant, n is temperature-independent stress exponent, and sigma/sigma sub c is structure-compensated stress. A master curve resulted from a logarithmic plot of activation-energy-compensated strain rate versus structure-compensated stress.

Clements, L. L.↗

Compressibility Effects on the Growth and Structure of Homogeneous Turbulent Shear Flow

Compressibility effects within decaying isotropic turbulence and homogeneous turbulent shear flow have been studied using direct numerical simulation. The objective of this work is to increase our understanding of compressible turbulence and to aid the development of turbulence models for compressible flows. The numerical simulations of compressible isotropic turbulence show that compressibility effects are highly dependent on the initial conditions. The shear flow simulations, on the other hand, show that measures of compressibility evolve to become independent of their initial values and are parameterized by the root mean square Mach number. The growth rate of the turbulence in compressible homogeneous shear flow is reduced compared to that in the incompressible case. The reduced growth rate is the result of an increase in the dissipation rate and energy transfer to internal energy by the pressure-dilatation correlation. Examination of the structure of compressible homogeneous shear flow reveals the presence of eddy shocklets, which are important for the increased dissipation rate of compressible turbulence.

Blaisdell, G. A.↗

Study of $\langle {p}_{\text{T}}\rangle$ and its higher moments, and extraction of the speed of sound in Pb-Pb collisions with ALICE

Ultrarelativistic heavy-ion collisions produce a state of hot and dense strongly interacting QCD matter called quark-gluon plasma (QGP). On an event-by-event basis, the volume of the QGP in ultracentral collisions is mostly constant, while its total entropy can vary significantly due to quantum fluctuations, leading to variations in the temperature of the system. Exploiting this unique feature of ultracentral collisions allows for the interpretation of the correlation of the mean transverse momentum ($\langle$p T $\rangle$) of produced charged hadrons and the number of charged hadrons as a measure for the speed of sound, c s . This speed is related to the rate at which compression waves travel in the QGP and is determined by fitting the relative increase in $\langle$p T $\rangle$ with respect to the relative change in the average charged-particle density ($\langle$dN ch /dη$\rangle$) measured at mid-rapidity. This study reports the event-average $\langle$p T $\rangle$ of charged particles as well as the variance, skewness, and kurtosis of the event-by-event transverse momentum per charged particle ([p T ]) distribution in ultracentral Pb-Pb collisions at a center-of-mass energy of 5.02 TeV per nucleon pair using the ALICE detector. Different centrality estimators based on charged-particle multiplicity or the transverse energy of the event are used to select ultracentral collisions. By ensuring a pseudorapidity gap between the region used to define the centrality and the region used to perform the measurement, the influence of biases and their potential effects on the rise of the mean transverse momentum is tested. The measured c$^{2}_{s}$ is found to strongly depend on the exploited centrality estimator and ranges between 0.1146±0.0028 (stat.)±0.0065 (syst.) and 0.4374±0.0006 (stat.)±0.0184 (syst.) in natural units. The self-normalized variance shows a steep decrease towards ultracentral collisions, while the self-normalized skewness variables show a maximum, followed by a fast decrease. These non-Gaussian features are understood in terms of the vanishing of the impact-parameter fluctuations contributing to the event-to-event [p T ] distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗