Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shear flow”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Deformation behavior of nanoscale Al–Al 2 Cu eutectics studied by in situ micropillar compression

We report deformation behavior of nanoscale laser processed Al–Al 2 Cu eutectics at room temperature is characterized through in situ micro-pillar compression testing in a scanning transmission microscope. Interlamellar spacing of Al–Al 2 Cu eutectics varies from hundreds of nanometers to 20 nm. Three different sizes of micro-pillars are fabricated in order to study the deformation behaviors of single colony and multiple colonies, corresponding to the single crystal and polycrystal respectively. For single colonies, lamellar orientations parallel, normal or inclined to the loading direction were tested. The main findings are: 1) the plasticity mechanisms strongly depend on loading orientation: buckling and kinking in the parallel-loaded eutectics, planar sliding along Al–Al 2 Cu lamellar interfaces in the incline-loaded eutectics and localized shearing in the normal-loaded eutectics. 2) the incline-loaded eutectics exhibits the lowest compression flow strength, and the normal-loaded eutectics has the highest compression flow strength. 3) with decreasing inter-lamellar spacing, the strength increases and plasticity is uniformly distributed, as opposed to shear localization. Highest compressive plasticity is observed in polycrystalline eutectics with an ensemble of lamellar orientations with ~20 nm average spacing: 17.9% at flow stress of 1.63 GPa, and degenerate, bimodal morphology: 11.1% at flow stress of 1.36 GPa.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Classical and quantum computing of shear viscosity for ( 2 + 1 ) D SU(2) gauge theory

We perform a nonperturbative calculation of the shear viscosity for ( 2 + 1 )-dimensional SU(2) gauge theory by using the lattice Hamiltonian formulation. The retarded Green’s function of the stress-energy tensor is calculated from real time evolution via exact diagonalization of the lattice Hamiltonian with a local Hilbert space truncation, and the shear viscosity is obtained via the Kubo formula. When taking the continuum limit, we account for the renormalization group flow of the coupling but no additional operator renormalization. We find the ratio of the shear viscosity and the entropy density η s is consistent with a well-known holographic result 1 4 π at several temperatures on a 4 × 4 honeycomb lattice with the local electric representation truncated at j max = 1 2 . We also find the ratio of the spectral function and frequency ρ x y ( ω ) ω exhibits a peak structure when the frequency is small. Both the exact diagonalization method and simple matrix product state classical simulation method beyond j max = 1 2 on bigger lattices require exponentially growing resources. So we develop a quantum computing method to calculate the retarded Green’s function and analyze various systematics of the calculation including j max truncation and finite size effects, Trotter errors and the thermal state preparation efficiency. Our thermal state preparation method still requires resources that grow exponentially with the lattice size, but with a very small prefactor at high temperature. We test our quantum circuit on both the Quantinuum emulator and the IBM simulator for a small lattice and obtain results consistent with the classical computing ones. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Scaling of Shear Rheology of Concentrated Charged Colloidal Suspensions across Glass Transition

Electrostatic interparticle interactions are a key component in controlling and designing rheological characteristics of concentrated charged colloidal suspensions. Herein, we investigate electroviscous effects on shear rheology using highly charged silica particles. By fixing the volume fraction but varying the salinity, the system undergoes a glass transition as evidenced by the evolution of the yield stress and zero-shear viscosity. We show that the steady shear viscosities obey a critical scaling relation that scales the flow curves into a super- and a sub-critical branch with glass transition salinity serving as the bifurcation point; we also demonstrate an isoviscosity scaling that collapses all isoviscosity lines into a single master curve that exhibits no singularity. Based on each scaling relation, in conjunction with common modeling equations, the quantitative relationships between the shear viscosity, stress, and salinity are established. This study demonstrates a new framework to model the steady shear rheology of concentrated charged colloids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Structural uniformity and compositional homogeneity of solid-phase alloyed rod

Solid-phase processes have emerged as an alternative to fusion-based alloying to avoid coarse microstructures, undesirable phase formation, and high energy consumption. However, achieving uniform distribution of alloying elements during friction-based processing remains challenging due to highly heterogeneous thermomechanical conditions. This work evaluates the structural uniformity and compositional homogeneity of Al–Cu–Zn alloyed rods produced by friction extrusion (FE) and establishes the role of the rotational speed to feed rate ratio (N/V) on alloying effectiveness. A systematic matrix of FE experiments was conducted at constant extrusion ratio with N/V values ranging from 3.7 to 300. Compositional uniformity was assessed along the rod length (ICP-OES), in three dimensions (X-ray computed tomography), and at the microscale (SEM–EDS), supported by a gray-level co-occurrence matrix (GLCM)–based homogeneity metric. Smoothed particle hydrodynamics (SPH) simulations were used to reveal material flow and thermomechanical fields. Results show that N/V = 100 produces a high-shear mixing zone that eliminates the unmixed core and enables near-full dissolution and dispersion of Cu and Zn. At lower N/V, a laminar flow region persists at the rod center, causing segregation and large composition gradients. The combined experimental–computational analysis provides mechanistic insight into the transition from fragmented particle dispersion to thermomechanically assisted metallurgical mixing. This study establishes processing–structure relationships for solid-phase alloying and provides guidance for achieving homogenized compositions comparable to wrought alloys via rapid, scalable FE processing.

Aluminum↗
DOE OSTI2022

The radial phase variation of reversed-shear and toroidicity-induced Alfvén eigenmodes in DIII-D

The eigenfunction of an instability contains information about energy flow in the wave. Here, the amplitude and phase of electron cyclotron emission radiometer data from hundreds of DIII-D reversed shear Alfvén eigenmodes (RSAE) and toroidicity-induced Alfvén eigenmodes (TAE) are analyzed along the outboard horizontal midplane. The radial phase profile can be flat, linearly rising or falling, convex or concave; in other words, a wide variety of shapes is observed. For a particular mode, often the radial phase profile remains approximately constant as the mode evolves in time but sometimes it changes rapidly. Many TAEs and some RSAEs have phase profiles that are rather flat where the mode amplitude is largest but rise steadily by ~2π at large major radius. Rapid phase changes are observed when the frequencies of an RSAE and TAE overlap and the modes couple. The phase profile depends weakly on the fast-ion gradient that would appear in the absence of wave-induced transport. Linear and quadratic fits to the phase profiles, together with many plasma parameters, are assembled into RSAE and TAE databases. In both cases, large variability is observed. For RSAEs, the strongest phase dependencies are on electron temperature T e , RSAE mode frequency, and the density of carbon impurities. For TAEs, the strongest dependencies are on beam power and major radius of the mode. In general, the average RSAE radial phase profile is essentially flat but the TAE profile has non-zero slope and curvature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Dynamics of Three-Dimensional Shock-Wave/Boundary-Layer Interactions

Advances in measuring and understanding separated, nominally two-dimensional (2D) shock-wave/turbulent-boundary-layer interactions (STBLI) have triggered recent campaigns focused on three-dimensional (3D) STBLI, which display far greater configuration diversity. Nonetheless, unifying properties emerge for semi-infinite interactions, taking the form of conical asymptotic behavior where shock-generator specifics become insignificant. The contrast between 2D and 3D separation is substantial; the skewed vortical structure of 3D STBLI reflects the essentially 2D influence of the boundary layer on the 3D character of the swept shock. As with 2D STBLI, conical interactions engender prominent spectral content below that of the turbulent boundary layer. However, the uniform separation length scale, which is crucial to normalizing the lowest-frequency dynamics in 2D STBLI, is absent. Comparatively, the spectra of 3D STBLI are more representative of the mid-frequency, convective, shear-layer dynamics in 2D, while phenomena associated with 2D separation-shock breathing are muted. Asymptotic behavior breaks down in many regions important to 3D-STBLI dynamics, occurring in a configuration-dependent manner. Aspects of inceptive regions near shock generators and symmetry planes are reviewed. Focused efforts toward 3D modal and nonmodalanalyses, moving-shock/boundary-layer interactions, fluid/structure interactions, and flow control are suggested as directions for future work.

42 ENGINEERING↗
DOE OSTI2022

Calibrating hypersonic turbulence flow models with the HIFiRE-1 experiment using data-driven machine-learned models.

In this paper we study the efficacy of combining machine-learning methods with projection-based model reduction techniques for creating data-driven surrogate models of computationally expensive, high-fidelity physics models. Such surrogate models are essential for many-query applications e.g., engineering design optimization and parameter estimation, where it is necessary to invoke the high-fidelity model sequentially, many times. Surrogate models are usually constructed for individual scalar quantities. However there are scenarios where a spatially varying field needs to be modeled as a function of the model’s input parameters. Here we develop a method to do so, using projections to represent spatial variability while a machine-learned model captures the dependence of the model’s response on the inputs. The method is demonstrated on modeling the heat flux and pressure on the surface of the HIFiRE-1 geometry in a Mach 7.16 turbulent flow. The surrogate model is then used to perform Bayesian estimation of freestream conditions and parameters of the SST (Shear Stress Transport) turbulence model embedded in the high-fidelity (Reynolds-Averaged Navier–Stokes) flow simulator, using shock-tunnel data. The paper provides the first-ever Bayesian calibration of a turbulence model for complex hypersonic turbulent flows. We find that the primary issues in estimating the SST model parameters are the limited information content of the heat flux and pressure measurements and the large model-form error encountered in a certain part of the flow.

42 ENGINEERING↗
DOE OSTI2026

Investigating shock-induced chemical reactions in Ni+Al multilayers: A continuum-based mesoscale approach with Arrhenius kinetics and artificial thermal conduction

A new continuum-based mesoscale modeling approach for shock-induced chemical reactions (SICRs) in Ni+Al multilayers is demonstrated in Sandia’s shock physics hydrocode, CTH. The approach utilizes Arrhenius-type kinetics and artificial thermal conduction. Our work builds upon previous efforts to parameterize equations of state for Ni x Al y [J. Appl. Phys. 137, 075102 (2025)], as well as simulations of inert shocks in realistic 2D microstructures [J. Appl. Phys. 137, 225301 (2025)]. To calibrate the reaction kinetics, pairs of the reaction coordinate, R′, vs time are extracted from the molecular dynamics (MD) literature. Here, the MD-informed kinetics are used to simulate the dynamic evolution of pressure and temperature in 2D mesoscale simulations. Overall, the MD-informed kinetics obtained for planar interfaces are too slow, as initial reaction is not observed on a nanosecond time scale. Even with quasi-periodic shock focusing leading to the formation of so-called “hot-spots,” the hot spots are unable to grow and coalesce using the fitted Arrhenius rate constants. However, by increasing the rate constants by two orders of magnitude, SICRs are observed at a shock pressure near 30 GPa, which is supported by experiments. Consequently, these mesoscale simulations suggest that unresolved shear-based mechanical mixing might possibly account for the discrepancies in kinetic rates, with shock-generated intense perturbations, interfacial vortical flows, and elevated temperatures serving as favorable reaction conditions. Future work will calibrate a shear-dependent reaction rate from the MD simulations with realistic interfaces that are wavy, diffuse, and disordered.

Kittell, David E. [Sandia National Laboratories (S↗
DOE OSTI2020

Hydrodynamic electron transport near charge neutrality

We develop the theory of hydrodynamic electron transport in a long-range disorder potential for conductors in which the underlying electron liquid lacks Galilean invariance. For weak disorder, we express the transport coefficients of the system in terms of the intrinsic kinetic coefficients of the electron liquid and the correlation function of the disorder potential. We apply these results to analyze the doping and temperature dependence of transport coefficients of graphene devices. In this work, we show that at charge neutrality, long-range disorder increases the conductivity of the system above the intrinsic value. The enhancement arises from the predominantly vortical hydrodynamic flow caused by local deviations from charge neutrality. Its magnitude is inversely proportional to the shear viscosity of the electron liquid and scales as the square of the disorder correlation radius. This is qualitatively different from the situation away from charge neutrality. In that case, the flow is predominantly potential and produces negative viscous contributions to the conductivity, which are proportional to the sum of shear and bulk viscosities and inversely proportional to the square of disorder correlation radius.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Evaluation of intracoronary hemodynamics identifies perturbations in vorticity

Coronary artery disease (CAD) is highly prevalent and associated with adverse events. Challenges have emerged in the treatment of intermediate coronary artery stenoses. These lesions are often interrogated with fractional flow reserve (FFR) testing to determine if a stenosis is likely to be causative for ischemia in a cardiac territory. This invasive test requires insertion of a pressure wire into a coronary vessel. Recently computational fluid dynamics (CFD) has been used to noninvasively assess fractional flow reserve in vessels reconstructed from medical imaging data. However, many of these simulations are unable to provide additional information about intravascular hemodynamics, including velocity, endothelial shear stress (ESS), and vorticity. We hypothesized that vorticity, which has demonstrated utility in the assessment of ventricular and aortic diseases, would also be an important hemodynamic factor in CAD. Three-dimensional (3D), patient-specific coronary artery geometries that included all vessels >1 mm in diameter were created from angiography data obtained from 10 patients who underwent diagnostic angiography and FFR testing (n = 9). A massively parallel CFD solver (HARVEY) was used to calculate coronary hemodynamic parameters including pressure, velocity, ESS, and vorticity. These simulations were validated by comparing velocity flow fields from simulation to both velocities derived from in vitro particle image velocimetry and to invasively acquired pressure wire-based data from clinical testing. There was strong agreement between findings from CFD simulations and particle image velocimetry experimental testing (p < 0.01). CFD-FFR was also highly correlated with invasively measured FFR (ρ = 0.77, p = 0.01) with an average error of 5.9 ± 0.1%. CFD-FFR also had a strong inverse correlation with the vorticity (ρ = -0.86, p = 0.001). Simulations to determine the effect of the coronary stenosis on intravascular hemodynamics demonstrated significant differences in velocity and vorticity (both p < 0.05). Further evaluation of an angiographically normal appearing non-FFR coronary vessel in patients with CAD also demonstrated differences in vorticity when compared with FFR vessels (p < 0.05). The use of highly accurate 3D CFD-derived intravascular hemodynamics provides additional information beyond pressure measurements that can be used to calculate FFR. Vorticity is one parameter that is modified by a coronary stenosis and appears to be abnormal in angiographically normal vessels in patients with CAD, highlighting a possible use-case in preventative screening for early coronary disease.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The Straightening of a River Meander Leads to Extensive Losses in Flow Complexity and Ecosystem Services

To assist river restoration efforts we need to slow the rate of river degradation. This study provides a detailed explanation of the hydraulic complexity loss when a meandering river is straightened in order to motivate the protection of river channel curvature. We used computational fluid dynamics (CFD) modeling to document the difference in flow dynamics in nine simulations with channel curvature (C) degrading from a well-established tight meander bend (C = 0.77) to a straight channel without curvature (C = 0). To control for covariates and slow the rate of loss to hydraulic complexity, each of the nine-channel realizations had equivalent bedform topography. The analyzed hydraulic variables included the flow surface elevation, streamwise and transverse unit discharge, flow velocity at streamwise, transverse, and vertical directions, bed shear stress, stream function, and the vertical hyporheic flux rates at the channel bed. The loss of hydraulic complexity occurred gradually when initially straightening the channel from C = 0.77 to C = 0.33 (i.e., the radius of the channel is three-times the channel width), and additional straightening incurred rapid losses to hydraulic complexity. Other studies have shown hydraulic complexity provides important riverine habitat and is positively correlated with biodiversity. This study demonstrates how hydraulic complexity can be gradually and then rapidly lost when unwinding a river, and hopefully will serve as a cautionary tale.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Shear strength and permeability of the cement-casing interface

Here, the shear strength and hydraulic permeability of the interface between well cement and casing was investigated using a triaxial direct shear apparatus. For the first time, these experiments provide measurements under controlled stress conditions with fluid flow measurements along the interface. The low cohesion (1.1 ± 1.1 MPa) and the high friction angle (43.4 ± 2.0°) indicates that the shear strength of the interface is provided by friction. This implies that the state of stress of the cement is critical to well integrity. The hydraulic aperture of the undamaged cement-steel samples was 6.8 ± 1.0 microns. Shear damage to the interface caused a decrease (-20 %) in hydraulic aperture for samples aged up to 1 month, and an increase (+300 %) for samples cured for two years. We performed numerical simulations to estimate the leakage potential from a carbon storage operation. This model predicts negligible leakage amounts (47 tonnes) in a shear-damaged well for the modeled injection of ~1.26 million tonnes of CO 2 . Thus, our measurements indicate that the cement-casing interface is not a significant leakage pathway in its intact or damaged state, and that shear-driven failure scenarios for this interface are not a significant risk to CO 2 storage security.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Polyelectrolyte complex scaffoldings for photocrosslinked hydrogels

Photocrosslinkable precursors (small molecules or polymers) undergo rapid crosslinking upon photoirradiation, forming covalently crosslinked hydrogels. The spatiotemporally controlled crosslinking, which can be achieved in situ, encourages the utility of photocrosslinked hydrogels in biomedicine as bioadhesives, bioprinting inks, and extracellular matrix mimics. However, the low viscosity of the precursor solutions results in unwanted flows and dilution, leading to handling difficulties and compromised strength of the photocrosslinked hydrogels. Here, we introduce oppositely charged triblock polyelectrolytes as additives for precursor solutions that transform them into self-assembled polyelectrolyte complex (PEC) hydrogels with enhanced shear strength and viscosity, providing interim protection against precursor dilution and mitigating secondary flows. The PEC network also augments the properties of the photocrosslinked hydrogels. Crosslinking of the precursors upon photoirradiation results in the formation of interpenetrating polymer network hydrogels with PEC and covalently-linked networks that exhibit shear moduli exceeding the linear combination of the moduli of the constituent networks and overcome the tensile strength–extensibility tradeoff that restricts the performance of covalently-linked hydrogels. In conclusion, the reinforcement approach is shown to be compatible with four types of photocrosslinkable precursors, does not require any modification of the precursors, and introduces minimal processing steps, paving the way for a broader translation of photocrosslinkable materials for biomedical applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Effect of outer-cylinder rotation on the radially heated Taylor–Couette flow

A Taylor–Couette setup with radial heating is considered where a Boussinesq fluid is sheared in the annular region between two concentric, independently rotating cylinders maintained at different temperatures. Linear stability analysis is performed to determine the Taylor number for the onset of instability. Two radius ratios corresponding to wide and thin gaps with several rotation rate ratios are considered. The rotation of the outer cylinder is found to have a general stabilizing effect on the stability threshold as compared to pure inner-cylinder rotation, with a few exceptions. The radial heating sets up an axial flow which breaks the reflection symmetry of isothermal Taylor–Couette flow in the axial coordinate. This symmetry breaking separates linear stability thresholds, and we find the fastest growing modes with both positive and negative azimuthal numbers for different parameters. Another important finding of the current study is the discovery of unstable modes in the Rayleigh-stable regime. Furthermore, closed disconnected neutral curves (CDNCs) are observed for both wide and thin gaps which can separate from or merge into open neutral–stability curves. Alternatively, CDNCs can also morph into open neutral stability curves as the rotation rate ratio is changed. CDNCs are observed to be sensitive to changes in control parameters, and their appearance/disappearance is shown to induce discontinuous jumps in the critical Taylor number. Finally, for both wide and thin gaps, the fastest-growing modes found in the pure corotation case are shown to have their origins in the instability islands at smaller values of rotation rate ratios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Reynolds stress scaling in the near-wall region of wall-bounded flows

A new scaling is derived that yields a Reynolds-number-independent profile for all components of the Reynolds stress in the near-wall region of wall-bounded flows, including channel, pipe and boundary layer flows. The scaling demonstrates the important role played by the wall shear stress fluctuations and how the large eddies determine the Reynolds number dependence of the near-wall turbulence behaviour.

42 ENGINEERING↗
DOE OSTI2023

Microscopic activated dynamics theory of the shear rheology and stress overshoot in ultradense glass-forming fluids and colloidal suspensions

We formulate a particle and force level, activated dynamics-based statistical mechanical theory for the continuous startup nonlinear shear rheology of ultradense glass-forming hard sphere fluids and colloidal suspensions in the context of the elastically collective nonlinear Langevin equation approach and a generalized Maxwell model constitutive equation. Activated structural relaxation is described as a coupled local-nonlocal event involving caging and longer range collective elasticity which controls the characteristic stress relaxation time. Theoretical predictions for the deformation-induced enhancement of mobility, the onset of relaxation acceleration at remarkably low values of stress, strain, or shear rate, apparent power law thinning of the steady-state structural relaxation time and viscosity, a nonvanishing activation barrier in the shear thinning regime, an apparent Herschel–Buckley form of the shear rate dependence of the steady-state shear stress, exponential growth of different measures of a yield or flow stress with packing fraction, and reduced fragility and dynamic heterogeneity under deformation were previously shown to be in good agreement with experiments. The central new question we address here is the defining feature of the transient response—the stress overshoot. In contrast to the steady-state flow regime, understanding the transient response requires an explicit treatment of the coupled nonequilibrium evolution of structure, elastic modulus, and stress relaxation time. We formulate a new quantitative model for this aspect in a physically motivated and computationally tractable manner. Theoretical predictions for the stress overshoot are shown to be in good agreement with experimental observations in the metastable ultradense regime of hard sphere colloidal suspensions as a function of shear rate and packing fraction, and accounting for deformation-assisted activated motion appears to be crucial for both the transient and steady-state responses.

Mechanics↗
DOE OSTI2023

Large Eddy Simulation of Rotationally Induced Ingress and Egress around an Axial Seal between Rotor and Stator Disks

In gas turbines, the hot gas exiting the combustor can have temperatures as high as 2000 °C, and some of this hot gas enter into the space between the stator and rotor disks (wheelspace). Since the entering hot gas could damage the disks, its ingestion must be minimized. This is carried out by rim seals and by introducing a cooler flow from the compressor (sealing flow) into the wheelspace. Ingress and egress into rim seals are driven by the stator vanes, the rotor and its rotation, and the rotor blades. This study focuses on the ingress and egress driven by the rotor and its rotation. This is carried out by performing wall-resolved large eddy simulation (LES) around an axial seal in a rotor–stator configuration without vanes and blades. Results obtained show the mechanisms by which the rotor and its rotation induce ingress, egress, and flow trajectories. Kelvin–Helmholtz instability was found to create a wavy shear layer and displacement thickness that produces alternating regions of high and low pressures around the rotor side of the seal. Vortex shedding on the backward-facing side of the seal and its impingement on the rotor side of the seal also produces alternating regions of high and low pressures. The locations of the alternating regions of high and low pressures were found to be statistically stationary and to cause ingress to start on the rotor side of the seal. Vortex shedding and recirculating flow in the seal clearance also cause ingress by entrainment. With the effects of the rotor and its rotation on ingress and egress isolated, this study enables the effects of stator vanes and rotor blades to be assessed.

42 ENGINEERING↗
DOE OSTI2023

Shear band velocity and activation volume during shear deformation by acoustic emission in a Zr-based bulk metallic glass

Here, an analysis is presented for the shear band velocity of plastic deformation and activation volume of the shear transformation zone (STZ) in a Zr 60 Cu 30 Al 10 bulk metallic glass (BMG). Compared to the values calculated by the stick-slip model, the values of shear band velocity calculated after optimization by the acoustic emission technique are much larger than those reported previously. Plastic flow in BMGs is found to be described by a power-law correlation between the maximum velocity of the shear band and the strain rates. Combining the free volume theory and the STZ theory, the relationship between the volume change in the STZ and the shear band velocity is derived as $Ω=A$(27–ln$v/l$). The value of A, which is a material dependent constant, is 5.43 × 10 –3 for the current BMG at room temperature.

36 MATERIALS SCIENCE↗