Feature Selection Clustering and Prototype Placement for Turbulence Data Sets.
Abstract not provided.
Barone, Matthew F.↗
SEARCH · Search NASA
Results for “strain tensor”
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.
Uncovering hidden bias in neutron diffraction residual strain measurements
When calculating residual strain via neutron or X-ray diffraction, uncertainties propagated from the peak fit are often inadequate to describe the true scatter of measurements about a singular strain state, such as one that should describe a macroscopic continuum. Because diffraction is inherently a selective process, orientation-dependent scatter arises from the sub-sampling of strong microstructure and strain gradients. This paper investigates the appropriateness of propagated uncertainties with reference to their original intention, i.e. noise about a mean value. Thirty-six unique orientations of strain measurements are taken at multiple locations within an additive friction-stir deposition component with fine-scale gradients (∼200 µm) of plastic strain, texture and residual elastic strain. Multiple strain and stress calculation pathways are compared: direct substitution of three measurements into Hooke's law, direct inversion of any six unique orientations into the strain state tensor and thirty-six measurement least-squares estimation. For the last two cases, the appropriateness of the uncertainty interval is statistically evaluated on the basis of a physical constraint: common agreement under the strain transformation law. For this sample, the direct inversion of six measurements retains a conservative estimate of the uncertainty. However, propagated uncertainties in the least-squares solution greatly underestimate the true experimental scatter. A simple pathway to estimate appropriate uncertainty intervals is suggested. These results demonstrate that the interpretation of uncertainty in residual strain is strongly dependent on intrinsic sample-dependent effects, and that oversampling orientations and statistical analysis can give more accurate results with realistic uncertainties.
36 MATERIALS SCIENCE↗
Using Particle Image Velocimetry to Determine Turbulence Model Parameters
The primary parameter of a standard k-ϵ model, C μ , was calculated from stereoscopic particle image velocimetry (PIV) data for a supersonic jet exhausting into a transonic crossflow. This required the determination of turbulent kinetic energy, turbulent eddy viscosity, and turbulent energy dissipation rate. Image interrogation was optimized, with different procedures used for mean strain rates and Reynolds stresses, to produce useful turbulent eddy viscosity fields. The eddy viscosity was calculated by a least-squares fit to all components of the three-dimensional strain-rate tensor that were available from the PIV data. This eliminated artifacts and noise observed when using a single strain component. Local dissipation rates were determined via Kolmogorov’s similarity hypotheses and the second-order structure function. The eddy viscosity and dissipation rates were then combined to determine C μ . Considerable spatial variation was observed in C μ , with the highest values found in regions where turbulent kinetic energy was relatively low but where turbulent mixing was important, e.g., along the high-strain jet edges and in the wake region. This suggests that use of a constant C μ in modeling may lead to poor Reynolds stress predictions at mixing interfaces. A data-driven modeling approach that can predict this spatial variation of C μ based on known state variables may lead to improved simulation results without the need for calibration.
42 ENGINEERING↗
Active- and transfer-learning applied to microscale-macroscale coupling to simulate viscoelastic flows
Active- and transfer-learning are applied to microscale dynamics of polymer flows for the multiscale discovery of effective constitutive approximations required in viscoelastic flow simulation. The result is macroscopic rheology directly connected to a microstructural model. Micro and macroscale simulations are adaptively coupled by means of Gaussian process regression (GPR) to run the expensive microscale computations only as necessary. This multiscale method is demonstrated with flows of a polymer solution as a model system. At the microscale level dissipative particle dynamics (DPD) is employed to model the fluid as a suspension of bead-spring micro-structures subjected to steady shear flow. The results yield the non-Newtonian viscosity and the first normal stress difference at strain rates as training data used in a GPR model. DPD parameters are calibrated with respect to experimental data for a real polymer solution. Compliance with these data requires adjustment of the DPD model's cutoff radius, which then becomes a function of the second invariant of the strain rate tensor. The FENE-P model is chosen for the macroscale description using the spectral element method (SEM) to simulate channel flow and flow past a circular cylinder. The DPD results at the lowest possible shear strain rate yield an estimate of the zero-shear rate viscosity, which allows the initiation of the macroscale flow by SEM as a Newtonian fluid. The resulting strain-rate field is surveyed to determine additional shear strain rate sampling points for the DPD system. This new information allows an initial fitting of parameters of the constitutive equation followed by new SEM simulations at the macroscale. Additionally, guided by active-learning GPR to select new sampling points, this process continues until convergence is achieved. The effectiveness of this new simulation paradigm for viscoelastic flows is tested with different macroscale operating conditions. The effective closure learned in the channel simulation is then transferred directly to the flow past a circular cylinder at low Reynolds number, where the results show that only two additional DPD simulations are required to achieve a satisfactory constitutive model. With an increase of the Reynolds number, the active-learning scheme automatically detects the inaccuracy of the learned constitutive model, and initiates additional DPD simulations for the extra data needed to once again close the microscale-macroscale coupled system. This new paradigm of active- and transfer-learning for multiscale modeling is readily applicable to other microscale-macroscale coupled simulations of complex fluids and other materials. Furthermore, the coupling between microscale and macroscale solvers can be seamlessly implemented with our open source multiscale universal interface (MUI) library.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Deep learning of structural morphology imaged by scanning X-ray diffraction microscopy
Scanning X-ray nanodiffraction microscopy is a powerful technique for spatially resolving nanoscale structural morphologies by diffraction contrast. One of the critical challenges in experimental nanodiffraction data analysis is posed by the convergence angle of nanoscale focusing optics which creates simultaneous dependency of the far-field scattering data on three independent components of the local strain tensor-corresponding to dilation and two potential rigid body rotations of the unit cell. All three components are in principle resolvable through a spatially mapped sample tilt series; however, traditional data analysis is computationally expensive and prone to artifacts. In this study, we implement NanobeamNN, a convolutional neural network specifically tailored to the analysis of scanning probe X-ray microscopy data. NanobeamNN learns lattice strain and rotation angles from simulated diffraction of a focused X-ray nanobeam by an epitaxial thin film and can directly make reasonable predictions on experimental data without the need for additional fine-tuning. We demonstrate that this approach represents a significant advancement in computational speed over conventional methods, as well as a potential improvement in accuracy over the current standard.
Luo, Aileen [Cornell Univ., Ithaca, NY (United Sta↗
Modal analysis of the behavior of inertial particles in turbulence subjected to Stokes drag
The clustering of small heavy inertial particles subjected to Stokes drag in turbulence, which is the concern of present study, is known to be minimal at small and large Stokes number and substantial at St = $\textit{O}$(1). This nonmonotonic trend, which has been shown computationally and experimentally, is yet to be explained analytically. In this study, we obtain an analytical expression for the Lyapunov exponents that quantitatively predicts this trend. The sum of the exponents, which is the normalized rate of change of the signed volume of a small cloud of particles, is correctly predicted to be negative and positive at small and large Stokes numbers, respectively, asymptoting to $\textit{τQ}$ as $\textit{τ}$ → 0 and $τ^{–1/2}|Q|^{1/4}$ as $\textit{τ} → ∞$, where $\textit{τ}$ is the particle relaxation time and $\textit{Q(τ)}$ is the difference between the norm of the rotation- and strain-rate tensors computed along the particle trajectory, which remains negative at all $\textit{τ}$ in turbulent flows. Additionally, the trajectory crossing is predicted only in hyperbolic flows where Q < 0 for sufficiently inertial particles with a $\textit{τ}$ that scales with $|Q|^{–1/2}$. Following the onset of crossovers, a transition from clustering to dispersion is predicted correctly. We show these behaviors are not unique to three-dimensional isotropic turbulence and can be reproduced closely by a one-dimensional mono-harmonic flow, which appears as a fundamental canonical problem in the study of particle clustering. Analysis of this one-dimensional canonical flow shows that the rate of clustering, quantified as the product of the Lyapunov exponent and particle relaxation time, is bounded by –1/2, behaving with extreme nonlinearity in the hyperbolic flows and always remaining positive in the elliptic flows. Finally, this model problem also confirms that dispersion in a hyperbolic flow will occur only if particles' trajectory cross. These findings, which stem from our analysis, are corroborated by the direct numerical simulations.
42 ENGINEERING↗
A novel transformation of the ice sheet Stokes equations and some of its properties and applications
We introduce a novel transformation of the Stokes equations into a form closely resembling the shallow Blatter–Pattyn equations. The two forms differ by only a few additional terms, while their variational formulations differ only by a single term in each horizontal direction. Specifically, the variational formulation of the Blatter–Pattyn model drops the vertical velocity in the second invariant of the strain rate tensor. Here we make use of the new transformation in two ways. First, we consider incorporating the transformed equations into a code that can be very easily converted from a Stokes to a Blatter–Pattyn model, and vice versa, by switching these terms on or off. This may be generalized so that the Stokes model is switched on adaptively only where the Blatter–Pattyn model loses accuracy. Second, the key role played by the vertical velocity in the Blatter–Pattyn approximation motivates new approximations. Two examples are presented. These require a mesh that enables the discrete continuity equation to be invertible for the vertical velocity in terms of the horizontal velocity components. Examples of such meshes, such as the first-order P1–E0 mesh and the second-order P2–E1 mesh, are given in both 2D and 3D. However, the transformed Stokes model has the same type of gravity forcing as the Blatter–Pattyn model, determined by the ice surface slope, thereby forgoing some of the mesh generality of the traditional formulation of the Stokes model.
58 GEOSCIENCES↗
Quantitative kinetic rules for plastic strain-induced α - ω phase transformation in Zr under high pressure
Plastic strain-induced phase transformations (PTs) and chemical reactions under high pressure are broadly spread in modern technologies, friction and wear, geophysics, and astrogeology. However, because of very heterogeneous fields of plastic strain $E$ p and stress σ tensors and volume fraction c of phases in a sample compressed in a diamond anvil cell (DAC) and impossibility of measurements of σ and $E$ p , there are no strict kinetic equations for them. Here, we develop a kinetic model, finite element method (FEM) approach, and combined FEM-experimental approaches to determine all fields in strongly plastically predeformed Zr compressed in DAC, and specific kinetic equation for α-ω PT consistent with experimental data for the entire sample. Since all fields in the sample are very heterogeneous, data are obtained for numerous complex 7D paths in the space of 3 components of the plastic strain tensor and 4 components of the stress tensor. Kinetic equation depends on accumulated plastic strain (instead of time) and pressure and is independent of plastic strain and deviatoric stress tensors, i.e., it can be applied for various above processes. Our results initiate kinetic studies of strain-induced PTs and provide efforts toward more comprehensive understanding of material behavior in extreme conditions.
36 MATERIALS SCIENCE↗
Shear strain alters the structure and migration mechanism of self-interstitial atoms in copper
Here, we use atomistic modeling to show that externally applied shear strain causes the lowest energy self-interstitial atom (SIA) structure in copper (Cu) to change from a <100>-type dumbbell to a <110>-type dumbbell. Concurrently, SIA migration switches from the 3-D random walk characteristic of <100>-type dumbbells to a 1-D mechanism analogous to that of crowdion SIAs. Furthermore, the relative energies of these two dumbbell structures as a function of strain are well predicted using elastic dipole tensors computed at zero strain, indicating that examination of these tensors may be used to assess the likelihood of strain-induced SIA structure transitions in other materials. Changes in lowest energy SIA structures and associated migration mechanisms stand to impact predictions of SIA behavior in irradiated solids.
36 MATERIALS SCIENCE↗
Low-temperature state in strontium titanate microcrystals using in situ multireflection Bragg coherent x-ray diffraction imaging
Strontium titanate is a classic quantum paraelectric oxide material that has been widely studied in bulk and thin films. It exhibits a well-known cubic-to-tetragonal antiferrodistortive phase transition at 105 K, characterized by the rotation of oxygen octahedra. A possible second phase transition at lower temperature is suppressed by quantum fluctuations, preventing the onset of ferroelectric order. However, recent studies have shown that ferroelectric order can be established at low temperatures by inducing strain and other means. Here, we used in situ multireflection Bragg coherent x-ray diffraction imaging to measure the strain and rotation tensors for two strontium titanate microcrystals at low temperature. We observe strains induced by dislocations and inclusion-like impurities in the microcrystals. Based on radial magnitude plots, these strains increase in magnitude and spread as the temperature decreases. Pearson's correlation heat maps show a structural transition at 50 K, which could possibly be the formation of a low-temperature ferroelectric phase in the presence of strain. We do not observe any change in local strains associated with the tetragonal phase transition at 105 K.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Temporal dynamics of the alignment of the turbulent stress and strain rate
We report that the flux of energy between scales in turbulence can be cast as the result of the interaction between a turbulent stress and a rate of strain. Efficient energy transfer requires that these two tensors be oriented properly relative to each other, but previous work has shown that the instantaneous alignment between them is poor. Here, we consider the temporal dynamics of this alignment in a direct numerical simulation of isotropic turbulence. We show that the orientation of the stress lags behind that of the strain rate, both at a single location and along trajectories. However, the timescale of the reorientation of the stress in these cases does not follow the expected dynamical scaling with length scale. To capture the proper dynamical scaling, we reformulate the energy flux between scales using the right Cauchy-Green strain tensor and the second Piola-Kirchhoff stress tensor. Our results highlight the key role played by the deformation of fluid elements in the physics of the energy cascade, and suggest that their irreversible deformation is a Lagrangian manifestation of the cascade's broken time-reversal symmetry.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Utah FORGE: Phase 1a Tensor Strainmeter Data for the April, 2022 Stimulation of Well 16A(78)-32
Data from two Tensor Optical Fiber Strainmeters that were operational during Stages 1, 2, and 3 of the April, 2022 stimulation of well 16A(78)-32. Each csv file contains data from each stimulation stage (stage1, stage2, stage3) for both Phase 1a strainmeter installations (FS01, formerly FS-C, and FS02, formerly FS1-2) in human-readable comma-separated value text files. There are two header lines in each file describing the data contained in that column along with their units, respectively. Data have been decimated from 2 to 1 Hz to match the Pason data found in the linked GDR dataset below (16A78-32 Stimulation Pason Data). These files contain the time series spanning the same time interval as the Pason data as well as ambient data for 5 hours before the stimulation and 5 hours following shut in. The station locations were chosen based on their proximity to the borehole seismometers owned and operated by the University of Utah. See README.txt for more information.
15 GEOTHERMAL ENERGY↗
Shear-Induced Conformations of Salt-Free Polyelectrolytes in Semidilute Solutions
Here, we present coarse-grained molecular dynamics simulations of salt-free polyelectrolyte chains in semidilute solutions under simple shear flow, with full hydrodynamic interactions and explicit dipolar solvent. At equilibrium, chain orientation statistics follow a pseudo-Voigt distribution, and the structural correlation length and chain end-to-end vector autocorrelation function exhibit scaling behavior consistent with theoretical predictions for polyelectrolytes. Under shear, chains transition from coiled to stretched states and the end-to-end vector autocorrelation function reveals oscillatory dynamics at high Weissenberg numbers. Analysis of the gyration tensor and shear strain distributions identified three distinct chain populations with directional alignment along and against the shear gradient. Compared with their neutral polymer counterparts, polyelectrolytes exhibit stronger shear thinning and enhanced chain alignment under the same Weissenberg number, which is attributed to electrostatic interactions and shear-induced counterion release. These findings provide molecular insight into the distinct flow response of charged polymers and have implications for tailoring the rheological properties of polyelectrolyte-based materials.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Eulerian finite element implementations of a dislocation density-based continuum model
In Eulerian finite element simulations, the mesh moves relative to the material. After every change of position between the mesh and the material, the state variables are interpolated to the new mesh position, which is referred to as advection. Large strain crystal plasticity models are based on the multiplicative decomposition of the total deformation gradient. Further, the stress is evaluated as a function of the thermoelastic strain, temperature, and other state variables. Advection of tensor quantities, such as the strain, is coupled with possibly significant advection errors. In an effort to reduce the advection errors, we develop two rate forms of an established dislocation density-based continuum model. To that end, we replace the multiplicative decomposition of the deformation gradient with the additive decomposition of the velocity gradient, and define the stress rate instead of the total stress. The Eulerian implementation is compared with Lagrangian calculations, and two numerical examples with severe deformation levels are presented.
36 MATERIALS SCIENCE↗
Tensor Basis Gaussian Process Models of Hyperelastic Materials
In this study, we develop Gaussian process regression (GPR) models of isotropic hyperelastic material behavior. First, we consider the direct approach of modeling the components of the Cauchy stress tensor as a function of the components of the Finger stretch tensor in a Gaussian process. We then consider an improvement on this approach that embeds rotational invariance of the stress-stretch constitutive relation in the GPR representation. This approach requires fewer training examples and achieves higher accuracy while maintaining invariance to rotations exactly. Finally, we consider an approach that recovers the strain-energy density function and derives the stress tensor from this potential. Although the error of this model for predicting the stress tensor is higher, the strain-energy density is recovered with high accuracy from limited training data. The approaches presented here are examples of physics-informed machine learning. They go beyond purely data-driven approaches by embedding the physical system constraints directly into the Gaussian process representation of materials models.
36 MATERIALS SCIENCE↗
Strongly Anisotropic Thermomechanical Response to Shock Wave Loading in Oriented Samples of the Triclinic Molecular Crystal 1,3,5-Triamino-2,4,6-trinitrobenzene
All-atom molecular dynamics (MD) simulations were used to study shock wave loading in oriented single crystals of the highly anisotropic triclinic molecular crystal 1,3,5-triamino-2,4,6-trinitrobenzene (TATB). The crystal structure consists of planar hydrogen-bonded sheets of individually planar TATB molecules that stack into graphitic-like layers. Shocks were studied for seven systematically prepared crystal orientations with limiting cases that correspond to shock propagation exactly perpendicular and exactly parallel to the graphitic-like layers. The simulations were performed for initially defect-free crystals using a reverse-ballistic configuration that generates explicit, supported shocks. Final longitudinal stress components are between ≈8.5 and ≈10.5 GPa for the 1.0 km s –1 impact speed studied. Orientation-dependent properties are reported including shock speeds, stresses, temperatures, compression ratios, and local material strain rates. Spatiotemporal maps of the temperature, stress tensor, material flow, and molecular orientations reveal complicated processes that arise for specific shock directions. Furthermore, the results indicate that TATB shock response is highly sensitive to crystal orientation, with significant qualitative differences for the time evolution of the stress tensor and temperature, elastic/inelastic compression response, defect formation and growth, critical von Mises stress, and strain rates during shock rise that span nearly an order of magnitude. A variety of inelastic deformation mechanisms are identified, ranging from crumpling of graphitic-like layers to dislocation-mediated plasticity to intense shear strain localization. To our knowledge, these are the first systematic MD simulations and analysis of explicit shock wave propagation along nontrivial crystal directions in a triclinic molecular crystal.
36 MATERIALS SCIENCE↗
Stress and Strain Heterogeneity and Persistence in Uniaxially‐ and Triaxially‐Loaded Sandstone
Two critical questions in brittle rock mechanics are how rocks developed localized strains and to what extent internal stress heterogeneity controls this localization and subsequent macroscopic failure. Definitive answers have not yet been found, but would provide insight into rock fracture mechanics as relevant to hydrocarbon extraction and sequestration. Here, we use synchrotron X‐ray tomography (XRT) and 3D X‐ray diffraction (3DXRD) during uniaxial and triaxial tests on Nugget and Bentheimer sandstones to examine strain and stress localization prior to mechanical failure. 3DXRD was used to measure intra‐granular lattice strains which were used to compute elastic stress tensors of each grain. Digital volume correlation (DVC) was applied to XRT images to determine the strain field in the sample. Both samples featured marked spatial heterogeneity, localization, and temporal persistence of elevated stresses and strains during their mechanical deformation toward failure. Both samples featured a majority of grains with at least one principal stress component that was tensile, a signature of the influence of heterogeneity on stress transmission. Measurements further revealed that compressive stress orientations and statistics evolved in a similar manner to those of inter‐particle forces in loose granular materials, with triaxially‐compressed rock exhibiting enhanced grain stress heterogeneity compared to uniaxially‐compressed rock. Our results complement recent work by others who employed XRT and scanning 3DXRD to study triaxially‐compressed sandstone, but extend those results to uniaxial compression, sandstones of varied porosity, and grain stress measurements throughout the 3D full extent of the samples rather than in a single layer examined with scanning 3DXRD.
58 GEOSCIENCES↗
A magnetic analog of pressure–strain interaction
We study the evolution equation for magnetic energy density for a non-relativistic magnetized plasma in the (Lagrangian) reference frame comoving with the electron bulk velocity. Analyzing the terms that arise due to the ideal electric field, namely, perpendicular electron compression and magnetic field line bending, we recast them to reveal a quantity with a functional form analogous to the often-studied pressure–strain interaction term that describes one piece of internal energy density evolution of the species in a plasma, except with the species pressure tensor replaced by the magnetic stress tensor. We dub it the “magnetic stress–strain interaction.” We discuss decompositions of the magnetic stress–strain interaction analogous to those used for pressure–strain interaction. These analogies facilitate the interpretation of the evolution of the various forms of energy in magnetized plasmas and should be useful for a wide array of applications, including magnetic reconnection, turbulence, collisionless shocks, and wave–particle interactions. We display and analyze all the terms that can change magnetic energy density in the Lagrangian reference frame of the electrons using a particle-in-cell simulation of magnetic reconnection.
Barbhuiya, M. Hasan (ORCID:0000000163301650)↗