Search NASA⌕ Search

SEARCH · Search NASA

Results for “STRESS 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.

At least 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Non-Boussinesq subgrid-scale model with dynamic tensorial coefficients

A major drawback of Boussinesq-type subgrid-scale stress models used in large-eddy simulations is the inherent assumption of alignment between large-scale strain rates and filtered subgrid-stresses. A priori analyses using direct numerical simulation (DNS) data have shown that this assumption is invalid locally as subgrid-scale stresses are poorly correlated with the large-scale strain rates [J. Bardina, J. Ferziger, and W. Reynolds, Improved subgrid-scale models for large-eddy simulation, in Proceedings of the 13th Fluid and Plasmadynamics Conference, AIAA (1980); C. Meneveau and K. Katz, Scale-invariance and turbulence models for large-eddy simulation, Ann. Rev. Fluid Mech. 32, 1 (2000)]. In the present work, a new, non-Boussinesq subgrid-scale model is presented where the model coefficients are computed dynamically. Some previous non-Boussinesq models have observed issues in providing adequate dissipation of turbulent kinetic energy [e.g., Bardina et al., Proceedings of the 13th Fluid and Plasmadynamics Conference (1980); R. A. Clark, J. Ferziger, and W.C. Reynolds. Evaluation of subgrid-scale models using an accurately simulated turbulent flow, J. Fluid Mech. 91, 1 (1979); S. Stolz and N. A. Adams, An approximate deconvolution procedure for large-eddy simulation, Phys. Fluids 11, 1699 (1999)]; however, the present model is shown to provide sufficient dissipation using dynamic coefficients. Modeled subgrid-scale Reynolds stresses satisfy the consistency requirements of the governing equations for large-eddy simulation (LES), vanish in laminar flow and at solid boundaries, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The new model, referred to as the dynamic tensor-coefficient Smagorinsky model (DTCSM), has been tested in simulations of canonical flows: decaying and forced homogeneous isotropic turbulence, and wall-modeled turbulent channel flow at high Reynolds numbers. The results show favorable agreement with DNS data. It has been shown that DTCSM offers similar predictive capabilities as the dynamic Smagorinsky model for canonical flows. In order to assess the performance of DTCSM in more complex flows, wall-modeled simulations of high Reynolds number flow over a Gaussian bump (Boeing speed bump) exhibiting smooth-body flow separation are performed. Predictions of surface pressure and skin friction, compared against DNS and experimental data, show improved accuracy from DTCSM in comparison to existing static coefficient (Vreman) and dynamic Smagorinsky model. The computational cost of performing LES with this model is up to 15% higher than the dynamic Smagorinsky model.

42 ENGINEERING↗
DOE OSTI2022

Nonlinear deformation and elasticity of BCC refractory metals and alloys

Application of isotropic pressure or uniaxial strain alters the elastic properties of materials; sufficiently large strains can drive structural transformations. Linear elasticity describes stability against infinitesimal strains, while nonlinear elasticity describes the response to finite deformations. Here, it was previously shown that uniaxial strain along [100] drives refractory metals and alloys towards mechanical instabilities. These include an extensional instability, and a symmetry-breaking orthorhombic distortion caused by a Jahn-Teller-Peierls instability that splays the cubic lattice vectors. Here we analyze these transitions in depth. Eigenvalues and eigenvectors of the Wallace tensor identify and classify linear instabilities in the presence of strain. We show that both instabilities are discontinuous, leading to discrete jumps in the lattice parameters. We provide physical intuition for the instabilities by analyzing the changes in first-principles energy, stress, bond lengths, and angles upon application of strain. Electronic band structure calculations show differential occupation of bonding and antibonding orbitals, driven by the changing bond lengths and leading to the structural transformations. Strain thresholds for these instabilities depend on the valence electron count.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Research Introduction [Slides]

In the topic of Reverse Time Imaging, we proposed a new IC to reduce computation cost but reserve image resolution for distributed sensor networks. For Induced Seismicity in Oklahoma, we analyzed fault stress state analysis at state scale, and we applied machine learning techniques to polarity picking and seismicity rate forecasting. The results provide better understanding of fault properties, stress field, and the relationships among fault, stress state, injections, and potential seismic hazards. Lastly, for Microseismic Monitoring, we detected and located 770 low-frequency events (5-50 Hz): (1) Shallow events are highly clustered, consistent with the pathway from injection well 13-10A to monitoring well; moment tensor inversion shows dominant tensile cracking; (2) Deep events are scattered and show migration pattern to the basement; moment tensors show that most events are shear cracks.

58 GEOSCIENCES↗
DOE OSTI2023

Multi-scale investigation of short-range order and dislocation glide in MoNbTi and TaNbTi multi-principal element alloys

Refractory multi-principal element alloys (RMPEAs) are promising materials for high-temperature structural applications. Here, we investigate the role of short-range ordering (SRO) on dislocation glide in the MoNbTi and TaNbTi RMPEAs using a multi-scale modeling approach. Monte carlo/molecular dynamics simulations with a moment tensor potential show that MoNbTi exhibits a much greater degree of SRO than TaNbTi and the local composition has a direct effect on the unstable stacking fault energies (USFEs). From mesoscale phase-field dislocation dynamics simulations, we find that increasing SRO leads to higher mean USFEs and stress required for dislocation glide. The gliding dislocations experience significant hardening due to pinning and depinning caused by random compositional fluctuations, with higher SRO decreasing the degree of USFE dispersion and hence, amount of hardening. Finally, we show how the morphology of an expanding dislocation loop is affected by the applied stress.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Sprain energy consequences for damage localization and fracture mechanics

The 2023 smooth Lagrangian Crack-Band Model (slCBM), inspired by the 2020 invention of the gap test, prevented spurious damage localization during fracture growth by introducing the second gradient of the displacement field vector, named the “sprain,” as the localization limiter. The key idea was that, in the finite element implementation, the displacement vector and its gradient should be treated as independent fields with the lowest ( C 0 ) continuity, constrained by a second-order Lagrange multiplier tensor. Coupled with a realistic constitutive law for triaxial softening damage, such as microplane model M7, the known limitations of the classical Crack Band Model were eliminated. Here, we show that the slCBM closely reproduces the size effect revealed by the gap test at various crack-parallel stresses. To describe it, we present an approximate corrective formula, although a strong loading-path dependence limits its applicability. Except for the rare case of zero crack-parallel stresses, the fracture predictions of the line crack models (linear elastic fracture mechanics, phase-field, extended finite element method (XFEM), cohesive crack models) can be as much as 100% in error. We argue that the localization limiter concept must be extended by including the resistance to material rotation gradients. We also show that, without this resistance, the existing strain-gradient damage theories may predict a wrong fracture pattern and have, for Mode II and III fractures, a load capacity error as much as 55%. Finally, we argue that the crack-parallel stress effect must occur in all materials, ranging from concrete to atomistically sharp cracks in crystals.

Science & Technology - Other Topics↗
DOE OSTI2020

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↗
DOE OSTI2022

Generalized quasiharmonic approximation via space group irreducible derivatives

The quasiharmonic approximation (QHA) is the simplest nontrivial approximation for interacting phonons under constant pressure, bringing the effects of anharmonicity into temperature-dependent observables. Nonetheless, the QHA is often implemented with additional approximations due to the complexity of computing phonons under arbitrary strains, and the generalized QHA, which employs constant stress boundary conditions, has not been completely developed. In this work we formulate the generalized QHA, providing a practical algorithm for computing the strain state and other observables as a function of temperature and true stress. We circumvent the complexity of computing phonons under arbitrary strains by employing irreducible second-order displacement derivatives of the Born-Oppenheimer potential and their strain dependence, which are efficiently and precisely computed using the lone irreducible derivative approach. We formulate two complementary strain parametrizations: a discretized strain grid interpolation and a Taylor series expansion in symmetrized strain. We illustrate our approach by evaluating the temperature and pressure dependence of select elastic constants and the thermal expansion in thoria (ThO 2 ) using density functional theory with three exchange-correlation functionals. The QHA results are compared to our measurements of the elastic constant tensor using time-domain Brillouin scattering and inelastic neutron scattering. Our irreducible derivative approach simplifies the implementation of the generalized QHA, which will facilitate reproducible, data-driven applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

In-situ digital image correlation and thermal monitoring in directed energy deposition additive manufacturing

As in welding, directed energy deposition (DED) additive manufacturing (AM) generates complex residual stresses and distortions commensurate with the complexity of the scan pattern used for deposition. To date, measuring DED distortions on complex geometries has only been achieved post process, discarding the complex thermomechanical history that leads to that final material state. In this work, surround stereo digital image correlation (DIC) is used to 3D map surfaces and strain tensors in-situ in a powder-blown laser DED system. Infrared thermography is then projected onto these surfaces to record the full thermomechanical history of printed parts. DIC presents a unique challenge to DED AM, as no part exists at the beginning of deposition, which (a) prevents application of an appropriate speckle pattern and (b) denies the user a zero strain reference frame. Solutions to these problems are proposed and their limitations explored herein. In sum, this work presents a relatively low-cost solution to monitoring and optimizing the unique temporal artifacts induced by complex scan strategies that was previously unobtainable.

42 ENGINEERING↗
DOE OSTI2024

Stress inside the pion in holographic light-front QCD

In this work, we propose a method to compute the gravitational form factor D ( Q 2 ) in holographic QCD by exploiting the remarkable correspondence between semiclassical light-front QCD and semiclassical field theories in wrapped spacetime in five dimensions. The use of light-front holography bridges physics at large Q 2 as attained in light-front QCD and physics at small Q 2 where the coupling to scalar and tensor fields, e.g. glueballs, are dominant. As an application, we compute the D -term for the pion and compare the results with recent lattice simulations. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

The use of digital thread for reconstruction of local fiber orientation in a compression molded pin bracket via deep learning

A deep convolutional neural network (DCNN) was used for microstructure reconstruction using artificial intelligence (MR-AI) by predicting local average fiber orientation distributions (FOD) in a 3D prepreg platelet molded composite (PPMC) pin bracket. To train the MR-AI model, surface strain fields from residual stresses simulated in PPMC plates were used as the input to the DCNN. A training dataset included PPMC plates with various degrees of global fiber alignment, based on the information obtained from high-fidelity flow simulation of a pin bracket. Further, the MR-AI model was then deployed to analyze FOD in the 3D pin bracket by conducting thermo-elastic residual stress analysis. Initially, the MR-AI model was established entirely on the synthetic simulation data. Then, a μCT scan of a physically molded pin bracket was used to create a finite element model that provided data for additional validation of the DCNN model. For the μCT scan finite element pin bracket the MR-AI model predicted the distribution of fiber orientation tensor components with MAE of 0.10 indicating a global prediction error of 10%. For the flow simulated pin bracket, the MR-AI model predicted the distribution of fiber orientation tensor components with a global prediction error of 11%. The MR-AI model showed the ability to predict regions of varying alignment in the base and flange of the pin bracket. The proposed MR-AI methodology allows for rapid prediction of FOD in geometrically complex parts and offers a promising path to detecting unique fiber orientation states in molded components.

42 ENGINEERING↗
DOE OSTI2022

Fabric-enriched continuum breakage mechanics (F-CBM )

A fabric-enriched continuum breakage mechanics (F-CBM) framework is developed to examine the relation between elastic and inelastic anisotropy in granular materials subjected to high pressure. For this purpose, emphasis is given to the influence of the anisotropic fabric on the plastic cap regime – that is, a macroscopic signature of the inelasticity of granular materials widely dependent on the particle crushability. Anisotropic behaviour prior to yielding is introduced through a symmetric second-order fabric tensor embedded in the expression of the elastic energy potential. It is shown that, thanks to its energy-based formulation, the F-CBM framework provides a platform to predict the rotation and distortion of the yield cap of granular materials as an outcome of anisotropic strain energy storage prior to grain crushing. Parametric analyses indicate that both the degree of anisotropy and the orientation of the fabric have an impact on the energy release resulting from grain crushing and contribute to the resulting stress–strain behaviour. It is shown that, despite the use of very few parameters, the proposed framework can accurately predict the anisotropic yielding and inelastic stress–strain response solely on the basis of the directional properties inferred from measurements of elastic anisotropy.

Engineering↗
DOE OSTI2019

Mesoscale Design of Magnetoelectric Heterostructures and Nanocomposites

Multiferroic materials have a coexistence of at least two ferroic orders. The coupling interaction between the different order parameters can produce new phenomena, such as the magnetoelectric (ME) effect. The ME response is the appearance of an electric polarization P upon applying a magnetic field H, and/or vice versa. Magnetoelectricity has been observed as an intrinsic effect in some natural materials at low temperature; however, such single phase materials suffer from extremely weak ME exchange. Better alternatives are two-phase composites consisting of magnetostrictive and piezoelectric phases, which have been shown to have large ME coefficients. Said composites exploit a unique product ME tensor property that depends upon the individual strictions of the piezoelectric and magnetostrictive phases, the elastic stiffness and integrity of the interphase interfaces, the phase distribution and the dimensionality of its connectivity, and any elastic constraint; amongst other things. This research program between D. Viehland PI and S.Priya at Virginia Tech investigated the mesoscale physics and chemistry of such two-phase ME heterostructures. The results of the proposed research are relevant to new generations of multifunctional materials involving heterogeneous components. How the product ME property be optimized? The answers lies with the interphase-interfaces, where the atomic layers coordinate stress transfer from one phase to another. It was the primary goal of this research program to discover such fundamental relationship for two phase ME heterostructures and nanocomposites. This research program integrated all the components required to fully investigate the ME phenomenon: starting from synthesis of nanomaterials with diversified interfacial characteristics; to full electrical, magnetic, and structural characterizations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Computational micromechanics model based failure criteria for chopped carbon fiber sheet molding compound composites

Chopped carbon fiber sheet molding compound has a great potential in lightweight automotive, marine, and aerospace applications. One of the most challenging tasks is to predict the failure strength of the material due to its anisotropy and heterogeneity, as well as complex stress states in real-world working conditions. In this work, a novel constitutive model of carbon fiber chip is proposed to capture the pre- and post-failure behaviors under different loading modes. On this basis, we propose a new computational micromechanics model, which is calibrated and validated by uniaxial tensile, compressive, and in-plane shear experiments. Furthermore, a set of microstructures representative volume element (RVE) models under complex loading conditions are reconstructed to understand the relationship between the microstructure characteristics and the failure envelopes. Finally, several modified versions of classical failure criteria are proposed for anisotropic materials with consideration of the fiber orientation tensor. The modified Tsai-Wu failure criterion, which shows the best accuracy among all failure criteria, is highlighted in the comparative study.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Determination of single crystal thermal expansion in Uranium-6wt%Niobium shape memory alloy using in-situ diffraction and modeling of textured polycrystalline samples

In-situ neutron diffraction and ex-situ dilatometry experiments were performed on a chemically banded, quenched uranium-6wt% niobium shape memory alloy to study the impact of deformation-(detwinning-) induced texture on the evolution of thermally induced strains (and associated stresses). Thermal heating and cooling cycles between 5K and 473K were performed in-situ with neutron diffraction, and the lattice strain evolution is reported for the observed monoclinic α" phase. Comparisons between the measured diffraction strains and macroscopic dilatometry experiments reveal relationships between micro- and macro-level thermal expansion. Softening of the texture during heating suggests that twin boundary motion can accommodate the large internal thermal strains which are approximately 2% greater than that typically observed during the 175K heating interval used to age more randomly oriented polycrystalline material. Assuming weak constraint of neighboring grains in randomly textured polycrystals, the single crystal thermal expansion tensor is extracted from measurements of lattice strains over the range 5K to 473K. Predictions of polycrystalline thermal expansion, based upon this single crystal thermal expansion tensor, are shown to compare favorably with bulk thermal expansion observations of the as quenched microstructure. However, such a lower-bound estimate is insufficient to explain all aspects of the behavior of the textured material, where the matrix is not isotropic. In conclusion, it is hypothesized that relaxation processes which occur within the quenched microstructure during heating are responsible for the distinct thermal expansion behavior observed during the initial heating cycle compared to cooling and subsequent cycling.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The weakly nonlinear response and nonaffine interpretation of the Johnson–Segalman/Gordon–Schowalter model

We derive new analytical solutions for the non-affine Johnson-Segalman/Gordon-Schowalter (JS/GS) constitutive equation with a general relaxation kernel in medium-amplitude oscillatory shear (MAOS) deformation. The results show time-strain separable (TSS) nonlinearity, therefore providing new physically-meaningful interpretation to the heuristic TSS nonlinear parameter in MAOS. The upper-convected, lower-convected, and corotational Maxwell models are all subsets of the results presented here. The model assumes that the microscale elements causing stress in the material slip compared to the continuum deformation. We introduce a visualization of the non-affine deformation field that acts on stress-generating elements to reinforce the hysical interpretation of the JS/GS class of models. Lastly, a case study is presented where previously published results, from fitting TSS models to MAOS data, can be re-interpreted based on the concept of non-affine motion of the JS/GS framework.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Viscoelastic relaxation and topological fluctuations in glass-forming liquids

In this study, a method for characterizing the topological fluctuations in liquids is proposed. This approach exploits the concept of the weighted gyration tensor of a collection of particles and permits the definition of a local configurational unit (LCU). The first principal axis of the gyration tensor serves as the director of the LCU, which can be tracked and analyzed by molecular dynamics simulations. Analysis of moderately supercooled Kob–Andersen mixtures suggests that orientational relaxation of the LCU closely follows viscoelastic relaxation and exhibits a two-stage behavior. The slow relaxing component of the LCU corresponds to the structural, Maxwellian mechanical relaxation. Additionally, it is found that the mean curvature of the LCUs is approximately zero at the Maxwell relaxation time with the Gaussian curvature being negative. This observation implies that structural relaxation occurs when the configurationally stable and destabilized regions interpenetrate each other in a bicontinuous manner. Finally, the mean and Gaussian curvatures of the LCUs can serve as reduced variables for the shear stress correlation, providing a compelling proof of the close connection between viscoelastic relaxation and topological fluctuations in glass-forming liquids.

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

A High-Throughput Solver for Marginalized Graph Kernels on GPU

Here, we present the design and optimization of a solver for efficient and high-throughput computation of the marginalized graph kernel on General Purpose GPUs. The graph kernel is computed using the conjugate gradient method to solve a generalized Laplacian of the tensor product between a pair of graphs. To cope with the large gap between the instruction throughput and the memory bandwidth of the GPUs, our solver forms the graph tensor product on-the-fly without storing it in memory. This is achieved by using threads in a warp cooperatively to stream the adjacency and edge label matrices of individual graphs by small square matrix blocks called tiles, which are then staged in registers and the shared memory for later reuse. Warps across a thread block can further share tiles via the shared memory to increase data reuse. We exploit the sparsity of the graphs hierarchically by storing only non-empty tiles using a coordinate format and nonzero elements within each tile using bitmaps. We propose a new partition-based reordering algorithm for aggregating nonzero elements of the graphs into fewer but denser tiles to further exploit sparsity. We carry out extensive theoretical analyses on the graph tensor product primitives for tiles of various density and evaluate their performance on synthetic and real-world datasets. Our solver delivers three to four orders of magnitude speedup over existing CPU-based solvers such as GraKeL and GraphKernels. The capability of the solver enables kernel-based learning tasks at unprecedented scales.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Machine-Learning-Based Multiscale Methods for 3D Modelling of Granular Materials by Incorporating History-Dependent State Variables

Over the past decades, the prevalence of machine learning (ML) methods has made the development of ML-based constitutive models for granular materials undoubtedly a popular subject. Numerous studies have been made to feature the loading path or history-dependent stress-strain response of granular media using neural networks. In this work, a novel finite element method (FEM)–ML multiscale approach was developed by incorporating internal variables to improve the simulation accuracy of 3D history-dependent granular materials for the first time. To this end, a surrogate constitutive model based on the single-step-based multi-layer perceptron (MLP) neural network was used to replace representative volume element (RVE) simulations conducted by the discrete element method (DEM) in the multiscale FEM–DEM approach. Although the prediction principle of the MLP aligns with the FEM algorithm, artificially added internal variables are required to differentiate the loading history. To address this issue, history variables associated with the Frobenius norm are proposed to be fed into the MLP coupled with the strain tensor to extract the history-dependent behaviour of granular assemblies. The developed FEM–ML approach was demonstrated in 3D conventional triaxial compression (CTC) simulations. Compared to the multiscale FEM–DEM approach, the proposed FEM–ML method exhibits a significantly improved computational efficiency.

granular materials↗