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,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Extreme Elasticity Anisotropy in Molecular Glasses

Abstract Glasses are generally assumed to be isotropic and there are no literature reports of elastic anisotropy for molecular glasses. However, as glasses formed by physical vapor deposition can be structurally anisotropic, it is of interest to investigate the elastic anisotropy in these materials. Micro‐Brillouin light spectroscopy is used in several experimental geometries to determine the elastic stiffness tensors of three glass films of itraconazole vapor‐deposited at substrate temperatures ( T sub ) of 330, 315, and 290 K, respectively. Significant elastic anisotropy is observed and, in these glasses, the elastic anisotropy shows a strong correlation with the molecular orientation. The out‐of‐plane and in‐plane Young's moduli of the high T sub (330 K) sample, which features a predominantly vertical molecular orientation, exhibit a high anisotropy ratio of 2.2. The observed elastic anisotropy is much larger than those previously observed in liquid crystals and even many crystalline solids.

Cang, Yu↗
DOE OSTI2020

Exploration of Electrochemical Reactions at Organic–Inorganic Halide Perovskite Interfaces via Machine Learning in In Situ Time-of-Flight Secondary Ion Mass Spectrometry

The instability of hybrid organic–inorganic perovskite (HOIP) devices is one of the significant challenges preventing commercialization. Exploring these phenomena is severely limited by the complexity of the intrinsic electrochemistry of HOIPs, the presence of multiple volatile and mobile ionic species, and the possible role of environmentally induced reactions at surfaces and triple-phase junctions. Here, in situ studies of the electrochemistry of methylammonium lead bromide perovskite with the Au electrode interface are reported via light- and voltage-dependent time-of-flight secondary ion mass spectrometry (ToF-SIMS) imaging of lateral perovskite heterostructures. While ToF-SIMS allows for the visualization of the chemical composition along the surface and its evolution with light and electrical bias, the interpretation of the multidimensional data obtained is often limited due to strong correlations between chemical signatures and the need to track multiple peaks at once. Here, a machine learning workflow combining the Hough transform and non-negative matrix factorization and no-negative tensor decomposition is developed to avoid this limitation and extract salient features of associated chemical changes and to separate the light- and voltage-dependent dynamics. Finally, combining these in situ characterizations and the machine learning workflow provides comprehensive information on the chemical nature of moving species, ion accumulation, and interfacial electrochemical reactions in HOIP devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dual Optical Hyperbolicity of PdCoO 2 and PdCrO 2 Delafossite Single Crystals

Hyperbolic materials exhibit a very peculiar optical anisotropy with simultaneously different signs of the dielectric tensor components. This anisotropy allows the propagation of exotic surface-wave excitations like hyperbolic phonons and plasmon polaritons. While hyperbolic materials hold promise for applications in subwavelength photonics and enhanced light-matter interactions, their natural occurrence is limited to a few materials, often accompanied by significant dielectric losses and limited hyperbolic spectral bandwidth. Focusing on PdCoO 2 and PdCrO 2 delafossite transition-metal oxides, in this paper unique dual hyperbolic regimes are demonstrated:: one localized around a phonon absorption in the mid-infrared spectral region, and the other extending into the visible range. Both hyperbolic regimes show exceptional properties including low dissipation and high hyperbolic quality factors. These results pave the way for innovative applications of delafossite layered metals in subwavelength photonics, imaging, and sensing.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Giant Modulation of Refractive Index from Picoscale Atomic Displacements

Abstract It is shown that structural disorder—in the form of anisotropic, picoscale atomic displacements—modulates the refractive index tensor and results in the giant optical anisotropy observed in BaTiS 3 , a quasi‐1D hexagonal chalcogenide. Single‐crystal X‐ray diffraction studies reveal the presence of antipolar displacements of Ti atoms within adjacent TiS 6 chains along the c ‐axis, and threefold degenerate Ti displacements in the a – b plane. 47/49 Ti solid‐state NMR provides additional evidence for those Ti displacements in the form of a three‐horned NMR lineshape resulting from a low symmetry local environment around Ti atoms. Scanning transmission electron microscopy is used to directly observe the globally disordered Ti a–b plane displacements and find them to be ordered locally over a few unit cells. First‐principles calculations show that the Ti a – b plane displacements selectively reduce the refractive index along the ab ‐plane, while having minimal impact on the refractive index along the chain direction, thus resulting in a giant enhancement in the optical anisotropy. By showing a strong connection between structural disorder with picoscale displacements and the optical response in BaTiS 3 , this study opens a pathway for designing optical materials with high refractive index and functionalities such as large optical anisotropy and nonlinearity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Tunable, Homoepitaxial Hyperbolic Metamaterials Enabled by High Mobility CdO

We report propagating light exhibits hyperbolicity in strongly anisotropic materials where the principal components of the dielectric tensor are opposite in sign. While hyperbolicity occurs naturally in anisotropic polar dielectrics, wherein optical phonons along orthogonal crystal axes are nondegenerate, such optical anisotropy can also be engineered in hyperbolic metamaterials (HMMs): thin film superlattices of alternating dielectric and metallic layers. Contrasted with the severely limited tunability of natural hyperbolic materials, the hyperbolic behavior of HMMs can be tailored significantly both through superlattice design and material selection. However, so far HMMs have suffered from high optical losses, hindering their performance. In this report, broadly tunable (λ = 2–5 µm) Type I and II hyperbolic modes with low losses (quality (Q)-factors up to 19.7) are observed through attenuated total reflectance measurements of monolithic, homoepitaxial superlattices of high- and low-doped cadmium oxide (CdO). Further, the low losses offered by CdO enable the first demonstration of real-space imaging of hyperbolic plasmon polaritons in nanoresonators by scattering-type scanning near-field optical microscopy—previously only possible for hyperbolic phonon polariton materials. Tunable, low-loss CdO HMMs promise designability for applications such as on-chip photonics, super-resolution imaging (hyperlensing), enhanced emission, novel emitter designs, and possibly quantum nanophotonic and time variant metasurfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Nonlinear Dispersion Relation and Out‐of‐Plane Second Harmonic Generation in MoSSe and WSSe Janus Monolayers

Janus transition metal dichalcogenides are an emerging class of atomically thin materials with engineered broken mirror symmetry that gives rise to long‐lived dipolar excitons, Rashba splitting, and topologically protected solitons. They hold great promise as a versatile nonlinear optical platform due to their broadband harmonic generation tunability, ease of integration on photonic structures, and nonlinearities beyond the basal crystal plane. Here, second and third harmonic generation in MoSSe and WSSe Janus monolayers is studied. Polarization‐resolved spectroscopy is used to map the full second‐order susceptibility tensor of MoSSe, including its out‐of‐plane components. In addition, the effective third‐order susceptibility and the second‐order nonlinear dispersion close to exciton resonances for both MoSSe and WSSe are measured at room and cryogenic temperatures. This work sets a bedrock for understanding the nonlinear optical properties of Janus transition metal dichalcogenides and probing their use in the next‐generation on‐chip multifaceted photonic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of Stoichiometry on the Structure and Polarization of BaTiO 3

Barium titanate (BaTiO 3 ) is a material of interest for photonic device applications due to its strong optical non-linearity. However, BaTiO 3 -based devices have not found widespread adoption, in part due to the challenges associated with synthesizing high quality thin-films. Here, high-resolution scanning transmission electron microscope (STEM) imaging is used to investigate the atomic structure of both on- and off-stoichiometric BaTiO 3 synthesized by molecular beam epitaxy (MBE). Here, this investigation reveals an asymmetry in the way the BaTiO 3 atomic lattice accommodates off-stoichiometry growth and unveils features beyond what is expected from diffraction or surface characterization techniques. Excess titanium incorporates into the BaTiO 3 lattice to form pervasive defects despite titanium-rich films having a low surface roughness and high-quality appearance in diffraction. Excess barium forms a rough, water-soluble surface layer but does not significantly impact the quality of the BaTiO 3 lattice. STEM is used to map titanium atom displacement in real-space. The average displacement distance is 30–60 pm in the strained thin-films, higher than the <20 pm displacement in bulk BaTiO 3 . Additionally, the titanium atom displacement direction deviates from the c-axis of the unit cell, which may have implications for the material's electro-optic tensor and thus for electro-optic device design.

Cavanagh, Ashley E. [Harvard Univ., Cambridge, MA ↗
DOE OSTI2023

Large Enhancements in Optical and Piezoelectric Properties in Ferroelectric Zn 1- x Mg x O Thin Films through Engineering Electronic and Ionic Anharmonicities

Multifunctionality as a paradigm requires materials exhibiting multiple superior properties. Integrating second-order optical nonlinearity and large bandgap with piezoelectricity can, for example, enable broadband, strain-tunable photonics. Though very different phenomena at distinct frequencies, both second-order optical nonlinearity and piezoelectricity are third-rank polar tensors present only in acentric crystal structures. However, simultaneously enhancing both phenomena is highly challenging since it involves competing effects with tradeoffs. Recently, a large switchable ferroelectric polarization of ≈80 μC cm –2 was reported in Zn 1-x Mg x O films. Here, ferroelectric Zn 1-x Mg x O is demonstrated to be a platform that hosts simultaneously a 30% increase in the electronic bandgap, a 50% enhancement in the second harmonic generation (SHG) coefficients, and a near 200% improvement in the piezoelectric coefficients over pure ZnO. These enhancements are shown to be due to a 400% increase in the electronic anharmonicity and a ≈200% decrease in the ionic anharmonicity with Mg substitution. Precisely controllable periodic ferroelectric domain gratings are demonstrated down to 800 nm domain width, enabling ultraviolet quasi-phase-matched optical harmonic generation as well as domain-engineered piezoelectric devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Solid‐state nuclear magnetic resonance of spin‐9/2 nuclei 115 In and 209 Bi in functional inorganic complex oxides

Abstract Indium and bismuth are technologically important elements, in particular as oxides for optoelectronic applications. 115 In and 209 Bi are both I = 9/2 nuclei with high natural abundances and moderately high frequencies but large nuclear electric quadrupole moments. Leveraging the quadrupolar interaction as a measure of local symmetry and polyhedral distortions for these nuclei could provide powerful insights on a range of applied materials. However, the absence of reported nuclear magnetic resonance (NMR) parameters on these nuclei, particularly in oxides, hinders their use by the broader materials community. In this contribution, solid‐state 115 In and 209 Bi NMR of three recently discovered quaternary bismuth or indium oxides are reported, supported by density functional theory calculations, numerical simulations, diffraction and additional multinuclear ( 27 Al, 69,71 Ga, and 121 Sb) solid‐state NMR measurements. The compounds LiIn 2 SbO 6 , BiAlTeO 6 , and BiGaTeO 6 are measured without special equipment at 9.4 T, demonstrating that wideline techniques such as the QCPMG pulse sequence and frequency‐stepped acquisition can enable straightforward extraction of quadrupolar tensor information in I = 9/2 115 In and 209 Bi even in sites with large quadrupolar coupling constants. Relationships are described between the NMR observables and local site symmetry. These are amongst the first reports of the NMR parameters of 115 In, 121 Sb, and 209 Bi in oxides.

Griffith, Kent J.↗
DOE OSTI2024

Novel pore size-controlled, susceptibility matched, 3D-printed MRI phantoms

We report the design concept and fabrication of MRI phantoms, containing blocks of aligned microcapillaires that can be stacked into larger arrays to construct diameter distribution phantoms or fractured, to create a “powder-averaged” emulsion of randomly oriented blocks for vetting or calibrating advanced MRI methods, that is, diffusion tensor imaging, AxCaliber MRI, MAP-MRI, and multiple pulsed field gradient or double diffusion-encoded microstructure imaging methods. The goal was to create a susceptibility-matched microscopically anisotropic but macroscopically isotropic phantom with a ground truth diameter that could be used to vet advanced diffusion methods for diameter determination in fibrous tissues. Two-photon polymerization, a novel three-dimensional printing method is used to fabricate blocks of capillaries. Double diffusion encoding methods were employed and analyzed to estimate the expected MRI diameter. Susceptibility-matched microcapillary blocks or modules that can be assembled into large-scale MRI phantoms have been fabricated and measured using advanced diffusion methods, resulting in microscopic anisotropy and random orientation. This phantom can vet and calibrate various advanced MRI methods and multiple pulsed field gradient or diffusion-encoded microstructure imaging methods. We demonstrated that two double diffusion encoding methods underestimated the ground truth diameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A Poromechanical Framework for Internal Interactions Induced by Solid Inclusions

The framework of poromechanics is generalized to simulate the multiscale behavior of porous media subjected to internal loadings stemming from the growth of solid inclusions. This generalization is designed to enable the study of anisotropic internal stress generation from solid growth within the pores, while recovering isotropic fluid-induced loading as a particular case. For this purpose, a mathematical strategy to define constitutive tensors in a thermodynamically consistent form is proposed, thus offering new opportunities for determining the poromechanical properties of a porous solid through advanced experimentation or micromechanical models. The framework is specialized by means of established elastic solutions for single pore–matrix interaction, as well as through homogenization schemes considering the interaction among congruent pores. In particular, the second Eshelby solution and the Tanaka–Mori–Benveniste homogenization scheme are used to derive a microporoelastic model. At an elemental scale, the model is tested under mixed control conditions by replicating different scenarios of geomaterial testing. In addition, the model characteristics are outlined with reference to inelastic microscopic loadings replicating chemo-mechanical forcing, such as expansive crystal formation. Through a series of parametric analyses, it is shown that the microstructure of the pores significantly influences the properties of porous media. Most notably, it is shown that the effects of a solid forming within the pores depend in a highly nonlinear fashion on the constitutive characteristics of the inhomogeneities and can therefore not be readily quantified or predicted without models capturing the diverse multiscale interactions among pores, inhomogeneities, and matrix.

chemomechanics↗
DOE OSTI2020

Quantum Fourier transform revisited

Summary The fast Fourier transform (FFT) is one of the most successful numerical algorithms of the 20th century and has found numerous applications in many branches of computational science and engineering. The FFT algorithm can be derived from a particular matrix decomposition of the discrete Fourier transform (DFT) matrix. In this paper, we show that the quantum Fourier transform (QFT) can be derived by further decomposing the diagonal factors of the FFT matrix decomposition into products of matrices with Kronecker product structure. We analyze the implication of this Kronecker product structure on the discrete Fourier transform of rank‐1 tensors on a classical computer. We also explain why such a structure can take advantage of an important quantum computer feature that enables the QFT algorithm to attain an exponential speedup on a quantum computer over the FFT algorithm on a classical computer. Further, the connection between the matrix decomposition of the DFT matrix and a quantum circuit is made. We also discuss a natural extension of a radix‐2 QFT decomposition to a radix‐ d QFT decomposition. No prior knowledge of quantum computing is required to understand what is presented in this paper. Yet, we believe this paper may help readers to gain some rudimentary understanding of the nature of quantum computing from a matrix computation point of view.

Camps, Daan↗
DOE OSTI2020

Matrix-free preconditioning for high-order H (curl) discretizations

The greater arithmetic intensity of high-order finite element discretizations makes them attractive for implementation on next-generation hardware, but assembly of high-order finite element operators as matrices is prohibitively expensive. As a result, the development of general algebraic solvers for such operators has been an open research challenge. Fast matrix-free application of high-order operators has received significant attention in the literature in the context of Poisson-type problems, but preconditioners and solvers for inverting more general operators are not very well-developed. In this paper, we consider the problem of preconditioning a definite Maxwell operator at high polynomial order without assembling a matrix. We show that given efficient preconditioners for high-order H 1 finite element problems on the same mesh, efficient H(curl) preconditioners can be constructed in an auxiliary space framework. We demonstrate the resulting preconditioners in a practical setting with tensor-product basis functions on an unstructured mesh of quadrilaterals. Overall, our approach uses a sparsified H 1 solver constructed on a low-order mesh of the nodal points of the underlying high-order space, and we show that the resulting H(curl) preconditioner is effective at very high polynomial orders for two-dimensional model problems with complicated geometry, varying piecewise constant coefficients, and curved elements. The resulting preconditioner scales with nearly optimal O(p d+1 ) floating point operation count and optimal O(p d ) memory transfer requirements, outperforming existing Maxwell preconditioners in the high-order regime.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A Scalable Reduced‐Order Model for the Steady Navier–Stokes Equations

Scaling up new scientific technologies from laboratory to industry often involves demonstrating performance on a larger scale. Computer simulations can accelerate design and predictions in the deployment process, though traditional numerical methods are computationally intractable even for intermediate pilot plant scales. Recently, the component reduced order modeling method has been developed to tackle this challenge by combining projection reduced order modeling and discontinuous Galerkin domain decomposition. However, while many scientific or engineering applications involve nonlinear physics, this method has only been demonstrated for various linear systems. In this work, the component reduced order modeling method is extended to steady Navier–Stokes flow, with application to general nonlinear physics in view. The large‐scale, global domain is decomposed into a combination of small‐scale unit component. Linear subspaces for flow velocity and pressure are identified via proper orthogonal decomposition over sample snapshots collected from each small‐scale unit component. Velocity bases are augmented with a pressure supremizer to satisfy the inf–sup condition for stable pressure prediction. Two different nonlinear reduced order modeling methods are employed and compared for efficient evaluation of nonlinear advection: A third‐order tensor projection operator and the empirical quadrature procedure. The proposed method is demonstrated on the flow over arrays of five different unit objects, achieving a 23‐fold speedup with less than 4% relative error in domains up to 256 times larger than the unit components. Furthermore, a numerical experiment with the pressure supremizer strongly indicates the need for a supremizer for stable pressure prediction. A comparison between the tensorial approach and the empirical quadrature procedure revealed a slight advantage of the empirical quadrature procedure. The framework is compared with an alternating Schwarz‐based reduced‐order approach, demonstrating improved efficiency and robustness for the DG‐based global solver while retaining flexibility for sub‐scale iterative solvers. The method is further extended to a coupled advection–diffusion and Navier–Stokes system, illustrating its applicability to multi‐physics problems and its potential for more general, inter‐coupled nonlinear systems.

42 ENGINEERING↗
DOE OSTI2024

Data‐driven variational method for discrepancy modeling: Dynamics with small‐strain nonlinear elasticity and viscoelasticity

Abstract The effective inclusion of a priori knowledge when embedding known data in physics‐based models of dynamical systems can ensure that the reconstructed model respects physical principles, while simultaneously improving the accuracy of the solution in the previously unseen regions of state space. This paper presents a physics‐constrained data‐driven discrepancy modeling method that variationally embeds known data in the modeling framework. The hierarchical structure of the method yields fine scale variational equations that facilitate the derivation of residuals which are comprised of the first‐principles theory and sensor‐based data from the dynamical system. The embedding of the sensor data via residual terms leads to discrepancy‐informed closure models that yield a method which is driven not only by boundary and initial conditions, but also by measurements that are taken at only a few observation points in the target system. Specifically, the data‐embedding term serves as residual‐based least‐squares loss function, thus retaining variational consistency. Another important relation arises from the interpretation of the stabilization tensor as a kernel function, thereby incorporating a priori knowledge of the problem and adding computational intelligence to the modeling framework. Numerical test cases show that when known data is taken into account, the data driven variational (DDV) method can correctly predict the system response in the presence of several types of discrepancies. Specifically, the damped solution and correct energy time histories are recovered by including known data in the undamped situation. Morlet wavelet analyses reveal that the surrogate problem with embedded data recovers the fundamental frequency band of the target system. The enhanced stability and accuracy of the DDV method is manifested via reconstructed displacement and velocity fields that yield time histories of strain and kinetic energies which match the target systems. The proposed DDV method also serves as a procedure for restoring eigenvalues and eigenvectors of a deficient dynamical system when known data is taken into account, as shown in the numerical test cases presented here.

Masud, Arif↗
DOE OSTI2024

Atomic Structure and Dynamics of Unusual and Wide-Gap Phase-Change Chalcogenides: A GeTe 2 Case

Brain-inspired computing, reconfigurable optical metamaterials, photonic tensor cores, and many other advanced applications require next-generation phase-change materials (PCMs) with better energy efficiency and a wider thermal and spectral range for reliable operations. Germanium ditelluride (GeTe 2 ), with higher thermal stability and a larger bandgap compared to current benchmark PCMs, appears promising for THz metasurfaces and the controlled crystallization of atomically thin 2D materials. Using high-energy X-Ray diffraction supported by first-principles simulation, the atomic structure in semiconducting pulsed laser deposition films and metallic high-temperature liquids is investigated. The results suggest that the structural and chemical metastability of GeTe 2 , leading to disproportionation into GeTe and Te, is related to high internal pressure during a semiconductor–metal transition, presumably occurring in the supercooled melt. Similar phenomena are expected for canonical GeS 2 and GeSe 2 under high temperatures and pressures.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Investigation of Main Bearing Fatigue Estimate Sensitivity to Synthetic Turbulence Models Using a Novel Drivetrain Model Implemented in OpenFAST

ABSTRACT A coupled medium‐fidelity drivetrain model is developed and implemented in OpenFAST for a 10‐MW land‐based reference turbine. The implementation is verified against a fully coupled multibody wind turbine model, including a detailed drivetrain. The new model can simultaneously and accurately estimate main bearing loads and represent elastic bending of the drivetrain. It has low computational cost and is useful for early design phases, sensitivity analyses and complex systems like wind farms (where computational expense must be expended elsewhere). Here, the model is implemented for a monopile offshore wind turbine and used to investigate the sensitivity of main bearing basic rating life to different synthetic turbulence models. Large‐eddy simulations (LES) targeting stable, neutral, and unstable atmospheric conditions at below‐, near‐ and above‐rated wind speeds are used as a reference. The turbulence models recommended by the International Electrotechnical Commission, the Mann spectral tensor model, and the Kaimal spectral model with exponential coherence are fitted to the LES data. Additionally, a constrained turbulence generator, PyConTurb (short for Python Constrained Turbulence ), based on LES data, is applied in the aero‐hydro‐servo‐elastic simulations. Taking PyConTurb as the baseline, the Kaimal model significantly underestimates fatigue of the downwind main bearing, with between 10% and 40% less damage. The Mann model also underestimates the downwind main bearing fatigue by up to 30%. The upwind main bearing damage is driven by mean loads, and differences between models are less significant, although the trends are similar. Reasons for these discrepancies are investigated and attributed to differences in spatial and temporal variations among the turbulence models.

17 WIND ENERGY↗
DOE OSTI2020

An Improved Fracture Mechanics-Informed Multiscale Thermomechanical Damage Model for Ceramic Matrix Composites

This paper extends recent work done by the authors in modeling length scale-dependent damage behavior of ceramic matrix composites (CMCs) to include effects of local anisotropy introduced by matrix cracking. This model captures scale-dependent damage initiation and propagation behavior of the brittle matrix by employing internal state variable (ISV) theory within a multiscale modeling framework to obtain damaged matrix stress/strain constitutive relationships at each length scale. The damage ISV captures the effects of matrix cracking and growth by using fracture mechanics and the self-consistent scheme to determine the reduced stiffness of the cracked matrix. Matrix cracks, which activate when stress intensity factors near manufacturing induced cavities exceed the fracture toughness of the material, are assumed to be transversely isotropic in the plane of the crack, and matrix anisotropy occurs when the damaged stiffness tensor is rotated from the crack plane to the global axes. The crack progression and temporal evolution of the damage ISV are governed by fracture mechanics and crack growth kinetics. The model effectively captures first matrix cracking, which is the first significant deviation from linear elasticity. The nonlinear predictive capabilities of the material model are demonstrated for monolithic silicon carbide (SiC) and a 2D woven five-harness satin (5HS) carbon fiber SiC matrix (C/SiC) CMC.

Skinner, Travis↗