Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Pulsating Poiseuille flow of a cement slurry

In this paper, we investigate the pulsatile flow of a cement slurry in a pipe. The constitutive relation for the viscous stress tensor is based on the power-law model, where the shear viscosity not only depends on the shear rate but also the volume fraction of the cement particles. To solve for the volume fraction field, a convection-diffusion equation is used. The dimensionless form of the governing equations and the boundary conditions are solved numerically at different pulsatile cycles. Further, a parametric study is conducted to study the effect of different dimensionless numbers at different cycles on the velocity and the volume fraction profiles.

36 MATERIALS SCIENCE↗
DOE OSTI2023

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

Angle-resolved polarized Raman study of layered b-As x P x-1 alloys: Identification of As-P vibrational modes

In this study, the polarization-resolved Raman spectra of b-As x P x-1 flakes with varying arsenic concentration (x = 0, 0.4, and 0.8) are systematically investigated. A clear polarization dependence is observed in the Raman intensity of all the P-P, As-P, and As-As modes and confirmed by the first-principles calculations. The observed angle-dependence can be related to the structure of the Raman tensors of each peak. In particular, the Raman intensities corresponding to As-P vibration modes are studied in detail. The out of plane As-P mode ($A$$^{1}_{g}$) is clearly identified, but the butterfly-shaped behavior observed in the polar plot around 350 cm -1 is interpreted as the combination of the two in-plane vibrations, one along the armchair and the other along the zigzag direction. This work allows unambiguous assignment of the $A$$^{1}_{g}$, B 2g , and $A$$^{2}_{g}$ modes of the As-P vibrational modes and the importance of the consideration of the superposition of the in-plane B 2g and $A$$^{2}_{g}$ modes. Such study can be utilized to further characterize the anisotropic structural properties and analyze the local distribution of the As and P in b-As x P x-1 alloys which will play key roles in electronic and optical properties. Furthermore, this methodology can be further extended to alloyed two-dimensional material systems with in-plane anisotropy in general.

2D materials↗
DOE OSTI2021

Flow-driven spectral chaos (FSC) method for simulating long-time dynamics of arbitrary-order non-linear stochastic dynamical systems

Uncertainty quantification techniques such as the time-dependent generalized polynomial chaos (TD-gPC) use an adaptive orthogonal basis to better represent the stochastic part of the solution space (aka random function space) in time. However, because the random function space is constructed using tensor products, TD-gPC-based methods are known to suffer from the curse of dimensionality. Here, we introduce a new numerical method called the flow-driven spectral chaos (FSC) which overcomes this curse of dimensionality at the random-function-space level. The proposed method is not only computationally more efficient than existing TD-gPC-based methods but is also far more accurate. The FSC method uses the concept of enriched stochastic flow maps to track the evolution of a finite-dimensional random function space efficiently in time. To transfer the probability information from one random function space to another, two approaches are developed and studied herein. In the first approach, the probability information is transferred in the mean-square sense, whereas in the second approach the transfer is done exactly using a new theorem that was developed for this purpose. The FSC method can quantify uncertainties with high fidelity, especially for the long-time response of stochastic dynamical systems governed by ODEs of arbitrary order. Six representative numerical examples, including a nonlinear problem (the Van-der-Pol oscillator), are presented to demonstrate the performance of the FSC method and corroborate the claims of its superior numerical properties. Finally, a parametric, high-dimensional stochastic problem is used to demonstrate that when the FSC method is used in conjunction with Monte Carlo integration, the curse of dimensionality can be overcome altogether.

(nonlinear) stochastic dynamical systems↗
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

A conservative implicit-PIC scheme for the hybrid kinetic-ion fluid-electron plasma model on curvilinear meshes

We report that the hybrid kinetic-ion fluid-electron plasma model is widely used to study challenging multi-scale problems in space and laboratory plasma physics. Here, a novel conservative scheme for this model employing implicit particle-in-cell techniques is extended to arbitrary coordinate systems via curvilinear maps from logical to physical space. The scheme features a fully non-linear electromagnetic formulation with a multi-rate time advance - including sub-cycling and orbit-averaging for the kinetic ions. By careful choice of compatible particle-based kinetic-ion and mesh-based fluid-electron discretizations in curvilinear coordinates, as well as particle-mesh interpolations and implicit midpoint time advance, the scheme is proven to conserve total energy for arbitrary curvilinear meshes. In the electrostatic limit, the method is also proven to conserve total momentum for arbitrary curvilinear meshes. Although momentum is not conserved for arbitrary curvilinear meshes in the electromagnetic case, it is for an important subset of Cartesian tensor-packed meshes. The scheme and its novel conservation properties are demonstrated for several challenging numerical problems using different curvilinear meshes, including a merging flux-rope simulation for a space weather application, and a helical m = 1 mode simulation for magnetic fusion energy application.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Accelerating high-order mesh optimization using finite element partial assembly on GPUs

In this paper we present a new GPU-oriented mesh optimization method based on high order finite elements. Our approach relies on node movement with fixed topology, through the Target-Matrix Optimization Paradigm (TMOP) and uses a global nonlinear solve over the whole computational mesh, i.e., all mesh nodes are moved together. A key property of the method is that the mesh optimization process is recast in terms of finite element operations, which allows us to utilize recent advances in the field of GPU-accelerated high order finite element algorithms. For example, we reduce data motion by using tensor factorization and matrix-free methods, which have superior performance characteristics compared to traditional full finite element matrix assembly and offer advantages for GPU based HPC hardware. Furthermore, we describe the major mathematical components of the method along with their efficient GPU-oriented implementation. In addition, we propose an easily reproducible mesh optimization test that can serve as a performance benchmark for the mesh optimization community.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Sylvester-preconditioned adaptive-rank implicit time integrators for advection-diffusion equations with variable coefficients

Here, we consider the adaptive-rank integration of multi-dimensional time-dependent advection-diffusion partial differential equations (PDEs) with variable coefficients. We employ a standard finite-difference method for spatial discretization coupled with high-order diagonally implicit Runge-Kutta temporal schemes. The discrete equation is a generalized Sylvester equation (GSE), which we solve with a projection-based adaptive-rank algorithm structured around two key strategies: (i) constructing dimension-wise subspaces using a novel atypical extended Krylov strategy, and (ii) efficiently solving the basis coefficient matrix with a preconditioned GMRES solver. The low-rank decomposition is performed in 2D using SVD and with high-order SVD (HOSVD) in 3D to represent the tensor in a compressed Tucker format. For d-dimensional problems (here, d = 2 or 3), the computational complexity and memory storage of the approach are found numerically to scale as and $\mathscr{O}(Nr^2) + \mathscr{O} (r^{d+1})$ and $\mathscr{O}(Nr) + \mathscr{O} (r^{d})$, respectively, with the one-dimensional resolution and the maximal rank during the Krylov iteration (which we find to be largely independent of on our numerical examples). We present numerical examples that illustrate the advertised properties of the algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Magnetic pair distribution function and half polarized neutron powder diffraction at the HB-2A powder diffractometer

Local magnetic order and anisotropy are often central for understanding fundamental behavior and emergent functional properties in quantum materials and beyond. Advances in neutron powder diffraction experiments and analysis tools now allow for quantitative determination. Here, we demonstrate this here with complementary total neutron scattering and polarized neutron measurements on the HB-2A neutron powder diffractometer at the High Flux Isotope Reactor (HFIR). In recent years, magnetic pair distribution function (mPDF) analysis has emerged as a powerful technique for probing local magnetic spin ordering of magnetic materials. This method can be broadly applied to any magnetic material but is particularly effective for studying systems with short-range magnetic order, such as materials with reduced dimensionality, geometrically frustrated magnets, thermoelectrics, multiferroics, and correlated paramagnets. Magnetic anisotropy often underpins the short-range order adopted. Half-polarized neutron powder diffraction (pNPD) can be used to determine the local susceptibility tensor on the magnetic sites to quantify the magnetic anisotropy. Combining the techniques of mPDF and pNPD can therefore provide valuable insights into local magnetic behavior. A series of measurements optimized for these techniques are presented as exemplar cases focused on frustrated materials where short-range order dominates, these include measurements to ultra-low temperature (<100 mK) not typically accessible for such experiments.

Half polarized neutron scattering↗
DOE OSTI2021

On the implementation of dislocation reactions in continuum dislocation dynamics modeling of mesoscale plasticity

The continuum dislocation dynamics framework for mesoscale plasticity is intended to capture the dislocation density evolution and the deformation of crystals when subjected to mechanical loading. It does so by solving a set of transport equations for dislocations concurrently with crystal mechanics equations, with the latter being cast in the form of an eigenstrain problem. Incorporating dislocation reactions in the dislocation transport equations is essential for making such continuum dislocation dynamics predictive. A formulation is proposed to incorporate dislocation reactions in the transport equations of the vector density-based continuum dislocation dynamics. This formulation aims to rigorously enforce dislocation line continuity using the concept of virtual dislocations that close all dislocation loops involved in cross slip, annihilation, and glissile and sessile junction reactions. The addition of virtual dislocations enables us to accurately enforce the divergence free condition upon the numerical solution of the dislocation transport equations for all slip systems individually. A set of tests were performed to illustrate the accuracy of the formulation and the solution of the transport equations within the vector density-based continuum dislocation dynamics. Comparing the results from these tests with an earlier approach in which the divergence free constraint was enforced on the total dislocation density tensor or the sum of two densities when only cross slip is considered shows that the new approach yields highly accurate results. Bulk simulations were performed for a face centered cubic crystal based on the new formulation and the results were compared with discrete dislocation dynamics predictions of the same. The microstructural features obtained from continuum dislocation dynamics were also analyzed with reference to relevant experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Significant acceleration of solid-state NMR simulations via three-angle powder averaging

The anisotropic frequency shifts imparted onto the NMR resonance frequency depend on the spherical angular coordinates that describe the orientations of the NMR interaction tensors with respect to the applied magnetic field direction. Experiments performed using magic-angle spinning, however, gain a dependence on a third angle: the rotor phase γ. Traditionally, a carousel average is performed to integrate over γ, which leads to a slow convergence of intensities without contributing to the underlying powder patterns. Herein, we show an order of magnitude acceleration in computation time may be obtained by including the γ-averaging into the main powder average to eliminate redundant calculation of resonance frequencies.

magic-angle spinning↗
DOE OSTI2020

A review of the application of 2D isotropic-anisotropic correlation NMR spectroscopy in structural studies of chalcogenide glasses

The application of solid-state high-resolution NMR spectroscopy in the structural investigation of chalcogenide glasses in Ge/As/P/Si-X (X = S, Se,Te) systems has remained challenging even for the spin-1/2 nuclides ( 29 Si, 31 P, 77 Se, 125 Te), owing to their low natural abundance (except for 31 P), slow spin-lattice relaxation rate and large CSA and chemical shift distribution induced line broadening effects. However, most of these deleterious effects can be successfully overcome in two-dimensional (2D) isotropic-anisotropic correlation NMR experiments, especially when performed at high magnetic field and in conjunction with the Car-Purcell-Meiboom-Gill (CPMG) echo train acquisition. Here, we present a short introduction to the basic principles of such experiments and review their applications over the last decade in deciphering various short- and intermediate- range structural characteristics of chalcogenide glasses in S-Se, Se-Te, Ge-Se, As-Se and Si-Se systems as well as in investigating the molecular dynamics in a P-Se supercooled liquid. We anticipate possible future applications of these 2D isotropic-anisotropic correlation NMR experiments, particularly in conjunction with density functional theory-based calculations of NMR chemical shift tensor parameters and additional signal enhancement schemes, in addressing complex structural correlations and distributions in chalcogenide glasses.

2D NMR↗
DOE OSTI2021

Density functional theory calculations of the thermodynamic and kinetic properties of point defects in β -U

Density functional theory (DFT) calculations of the thermodynamic and kinetic properties of point defects in the β phase of uranium are reported. Defect energies and entropies were calculated using 2 x 2 x 2 supercells and the Generalized Gradient Approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) for the exchange-correlation potential. Due to computational cost, calculations of the vibrational properties governing entropies were performed by only displacing atoms within (roughly) the 3rd or 4th nearest neighbor shell of the defect, which implicitly assumes that atoms beyond this distance are unaffected by the defect. Migration barriers were estimated by nudged elastic band (NEB) calculations. The low symmetry of the β-U phase (the unit cell is tetragonal and contains 30 atoms) results in many point defect configurations and even more migration pathways. A connectivity map, starting from the most stable point defects, was developed in order to identify the rate-limiting step controlling the net diffusion rate in each crystallographic direction. The uranium self-diffusivity tensor was calculated by combining the defect formation energies, entropies, migration barriers and attempt frequencies. The fastest diffusion rate was determined to be a vacancy mechanism in the z crystallographic direction. The predicted uranium self-diffusivity for this mechanism agrees well with available experimental data. Finally, the diffusion mechanisms and rates identified in this study will inform fuel performance models of swelling and gas evolution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Crystal mechanics-based thermo-elastic constitutive modeling of orthorhombic uranium using generalized spherical harmonics and first-order bounding theories

In earlier works, a mathematical procedure for invertible microstructure-property linkages was developed using computationally efficient spectral methods for polycrystalline cubic and hexagonal metals. This paper formulates such invertible microstructure–property linkages for orthorhombic polycrystalline metals relying on the generalized spherical harmonics (GSH) spectral basis. The procedure is used to compute property closures of orthorhombic polycrystals. The closures represent the complete set of theoretically possible combinations of effective properties for a selected material. The procedure relies on the first-order bounding theories and considers orientation distribution functions (ODFs) as the main microstructural descriptor influencing homogenized properties. Numerous examples of these closures involving second-rank thermal expansion and fourth-rank elastic stiffness tensorial properties over a broad range of temperatures are presented for α-uranium (α-U). In doing so, certain key properties of these closures are exploited to facilitate their computation with drastically reduced computational effort. Along with the recently developed GSH-based interpolation procedure for ODFs from coarsely spaced experimental measurement grids to finely spaced finite element mesh resolution grids presented in Barrett et al., the developed computationally efficient ODF-effective property linkages are used to establish a crystal mechanics-based simulation framework coupled with the finite element method (FEM). The ODF dependent thermal expansion and elastic stiffness tensors are efficiently calculated at every integration point and used by the FEM to predict the overall distortion of a hemispherical part made of α-U during heating. In conclusion, it is shown that the developed framework can be used to simulate microstructurally heterogeneous components under thermo-mechanical loadings in a computationally efficient manner.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Anisotropic temperature-dependent elastic constants and thermal conductivities of TRISO particle coatings

Tristructural isotropic (TRISO) particles show great promise as a candidate fuel for use in several next-generation high-temperature nuclear reactor designs due to their structural integrity and fuel performance at high temperatures and burnups. Computational codes exist that can simulate TRISO fuel performance characteristics and failure probabilities under extreme conditions which require knowledge of the TRISO coatings’ thermophysical properties. The thermophysical descriptions of the TRISO particle's layers (i.e., buffer, pyrolytic carbon, and silicon carbide) currently used in fuel performance codes, however, assume that many of these properties are constant with respect to temperature or texture. In this paper, we obtain the full elastic stiffness tensors of the carbon and silicon carbide layers, which have transversely isotropic symmetry. Additionally, the calculated elastic properties of each of the coating layers are in remarkable agreement with the current models at room temperature and correct orientations. Additionally, the textured 3C-SiC layer was found to exhibit novel auxetic behavior above 1500 °C. The anisotropic temperature-dependent thermal conductivities of the carbon layers were calculated using acoustical Grüneisen-Debye theory which are in excellent agreement with current models at room temperature and correct orientations. These texture- and temperature-dependent relationships can be incorporated into the thermophysical description of TRISO particles in order to more accurately model fuel performance and failure probabilities under extreme conditions in forthcoming high-fidelity computational simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Measuring residual stress in Ti-6Al-4V with HR-EBSD, using reference patterns from annealed material

In this study, elastic strain and residual stress values in an engineering material were determined by high-resolution EBSD (HR-EBSD) using reference patterns from annealed, recrystallized material. Specifically, rolled Ti-6Al-4V sheet was analyzed using HR EBSD in both an as-annealed and a strained state, the latter established via a static two-point bend fixture. The commercial software CrossCourt4 (CC4) was used to evaluate the relative intragranular strain tensor components of the alpha grains mapped by HR-EBSD. The intragranular reference patterns identified in each grain by the software were then replaced with similarly-oriented patterns from scans of the as-annealed, recrystallized material. The CC4 analysis was then repeated to provide calibrated strain and residual stress values. Elastic strains on the order of 10–3 were detected, with corresponding residual stress values on the order of 1 GPa. Subsequent analyses using high energy synchrotron x-ray diffraction showed good correlation with the elastic strains measured by HR-EBSD using the annealed reference patterns. Additionally, the maximum macroscopic stress analytically predicted for two-point bending was in agreement with the residual stress observed. Good repeatability of the methodology and insight into the state of individual grains were also observed. The findings provide encouraging evidence that accurate intragranular and long-range residual stress measurements of engineering significance can be made with HR-EBSD using a physical reference pattern library.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Grain boundary metastability controls irradiation resistance in nanocrystalline metals

Grain boundaries (GBs) in polycrystalline materials are powerful sinks for irradiation defects. While standard theories assume that a GB’s efficiency as a sink is defined solely by its character before irradiation, recent evidence conclusively shows that the irradiation sink efficiency is a highly dynamic property controlled by the intrinsic metastability of GBs under far-from-equilibrium irradiation conditions. In this paper, we reveal that the denuded (i.e., defect-free) zone, typically the signature of a strong sink, can collapse as irradiation damage accumulates. We propose a radiation damage evolution model that captures this behavior based on the emergence of a series of irradiation defect-enabled metastable GB microstate changes that dynamically alter the ability of the GB to absorb further damage. We show that these microstate changes control further defect absorption and give rise to the formation of a defect network that manifests itself as a net Nye-tensor signal detectable via lattice curvature experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Tunable Second Harmonic Generation in Twisted Bilayer Graphene

Twisted stacking of van der Waals (vdW) materials introduces a new way in bandstructure engineering and has given rise to numerous extraordinary physical phenomena. Despite the absence of second harmonic generation (SHG) in non-gated monolayer graphene, artificially twisted bilayer graphene (tBLG) possesses more possible point-group symmetries including those with broken inversion-symmetry. Here, we report twist-angle-dependent SHG from tBLG, which is the first demonstration of an elemental material with intrinsically tunable nonlinearity. We show that depending on the twist angle, the susceptibility of the dominant chiral tensor component of tBLGs can vary from 0 to 28×10 4 pm 2 /V, which is at the same order of magnitude as on-resonance susceptibility of monolayer MoS 2 . These results shed light on the underlying symmetry of tBLG systems and electronic band-structure near van Hove singularities (vHs). More importantly, they introduce a new degree of freedom, the twisting, in creating efficient second-order nonlinear material form centrosymmetric constituents

36 MATERIALS SCIENCE↗