Search NASA⌕ Search

SEARCH · Search NASA

Results for “Cauchy”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Pressure dependence of Poisson’s ratio of glassy Baltic amber studied by Brillouin scattering spectroscopy

Amber is a typical natural glass with very long aging time. Elastic properties of amber at high pressures have been studied using Brillouin scattering and a diamond anvil cell. Both longitudinal and transverse acoustic modes have been observed up to 12 GPa. The pressure dependences of longitudinal, shear, Young’s, and bulk moduli, compressibility, and Poisson’s ratio were determined. The longitudinal, shear, Young’s, bulk moduli show the remarkable increase, and compressibility shows a marked decrease with increasing pressure. However, it is found that the pressure dependence of Poisson’s ratio is very small. The mechanism of this small pressure dependence was discussed. Finally, the Cauchy type relation between longitudinal and shear moduli was examined. Its coefficients show the deviation from the Cauchy relation owing to the coexistence of different intermolecular interactions in amber.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2020

Fifth-degree elastic energy for predictive continuum stress-strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

Levitas, Valery↗
DOE OSTI2022

Planar, longitudinal, compressive waves in solids: Thermodynamics and uniaxial strain restrictions

A new tri-axial pressure-based constitutive expression has been found using Cauchy's stress tensor. This stress state emphasizes pressure and shear stress. The description is a pressure plus an effective shear stress allowing for a constitutive law based on atomic solid-state phase changes in crystalline cells due to pressure plus shear-based dislocation motion commonly associated with plasticity. Pressure has a new role in the material's constitutive response as it is separated from plasticity. The thermo-mechanical system describes third-order Gibbs’ expressions without specific volume restrictions placed upon the material. Isothermally, the ratio of heat to shear work in elastic copper is shown to approach zero at a very low temperature and become larger than one as temperature approaches melting. Wave compression models investigated are elastic and plastic: in fully elastic materials, the planar wave is restricted by Poisson's effect although plastic shear changes this constraint. Plastic deformation, dominated by dissipative shear stresses in uniaxial strain, heats the material while excluding phase changes from hydrostatic pressures. The material properties per se across Hugoniot shocks are described with entropy concepts. Shock waves are exceedingly complex since the constitutive laws are linked at extreme temperatures, pressures, and shear stresses. Isothermal, isentropic, isochoric, and iso-shear conditions are used throughout with Jacobian algebra.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Sensitivity analysis of a layered piezoelectric system using ZFEM

The complex variable finite element method (ZFEM) is a numerical technique which aims to find the partial derivatives of the independent variables with respect to variation in dependent parameters declared in the physics. This is done by combining the complex Taylor series expansion within the weak formulation of the governing equation in a coupled system of linear equations forming a complex valued block matrix given by the Cauchy–Riemann matrix representation. In this work, two-dimensional linear first-order elements have been implemented in ZFEM to predict the design derivatives of the mechanical displacement field and the voltage potential field for a layered piezoelectric system in a steady-state study with Dirichlet boundary condition applied at the top and bottom edges of the geometry. This approach allows the standard FEM solution to quantify the sensitivity of the mechanical displacement and voltage potential fields with respect to small variations in the material properties through the information obtained from the computation of the derivatives. The domain is formed by a layered body with PZT-4 and PZT-5 stacked together. For result verification, the numerical solution obtained with ZFEM was compared to results from a commercial FEM package and the solution from the imaginary part was compared to the exact solution of a well-known benchmark problem. In conclusion, comparison of the results showed good agreement for both the real and imaginary parts of the solution and the largest sensitivities were found in PZT-5 specifically in C 13 , C 33 , and ε 33 .

42 ENGINEERING↗
DOE OSTI2026

Coarse-grained fixed-point tensor networks and holographic reflected entropy in 3D gravity

We use the framework of fixed-point BCFT tensor networks to present a microscopic CFT derivation of the correspondence between reflected entropy (RE) and entanglement wedge cross section (EW) in AdS 3 /CFT 2 , for both bipartite and multipartite settings. These fixed-point tensor networks, obtained by triangulating Euclidean CFT path integrals, allow us to explicitly construct the canonical purification via cutting-and-gluing CFT path integrals. Employing modular flow in the large-c limit, we demonstrate that these intrinsic CFT manipulations reproduce bulk geometric prescriptions, without assuming the AdS/CFT dictionary. The emergence of bulk geometry is traced to coarse-graining over heavy states in the large-c limit. Universal coarse-grained BCFT data for compact 2D CFTs, through the relation to Liouville theory with ZZ boundary conditions, yields hyperbolic geometry on the Cauchy slice. The corresponding averaged replica partition functions reproduce all candidate EWs, arising from different averaging patterns, with the dominant one providing the correct RE and EW. In this way, many heuristic tensor-network intuitions in toy models are made precise and established directly from intrinsic CFT data.

AdS-CFT correspondence↗
DOE OSTI2025

Wrong signs are alright

It has been shown that some Lorentz-invariant quantum field theories, such as those with higher-dimensional operators with negative coefficients, lead to superluminality on some classical backgrounds. While superluminality by itself is not logically inconsistent, these theories also predict the formation of closed time-like curves at the classical level, starting from initial conditions without such curves. This leads to the formation of a Cauchy Horizon which prevents a complete description of the time evolution of such systems. Inspired by the chronology protection arguments of General Relativity, we show that quantum mechanical effects from low energy quanta strongly backreact on such configurations, exciting unknown short-distance degrees of freedom and invalidating the classical predictions. Thus, there is no obvious low-energy obstruction to the existence of these operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Covariant bit threads

We derive several new reformulations of the Hubeny-Rangamani-Takayanagi covariant holographic entanglement entropy formula. These include: (1) a minimax formula, which involves finding a maximal-area achronal surface on a timelike hypersurface homologous to D(A) (the boundary causal domain of the region A whose entropy we are calculating) and minimizing over the hypersurface; (2) a max V-flow formula, in which we maximize the flux through D(A) of a divergenceless bulk 1-form V subject to an upper bound on its norm that is non-local in time; and (3) a min U-flow formula, in which we minimize the flux over a bulk Cauchy slice of a divergenceless timelike 1-form U subject to a lower bound on its norm that is non-local in space. The two flow formulas define convex programs and are related to each other by Lagrange duality. For each program, the optimal configurations dynamically find the HRT surface and the entanglement wedges of A and its complement. The V-flow formula is the covariant version of the Freedman-Headrick bit thread reformulation of the Ryu-Takayanagi formula. We also introduce a measure-theoretic concept of a “thread distribution”, and explain how Riemannian flows, V-flows, and U-flows can be expressed in terms of thread distributions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Simple holography in general spacetimes

The simple or “outermost” wedge in AdS is the portion of the entanglement wedge that can be reconstructed with sub-exponential effort from CFT data. Here we furnish a definition in arbitrary spacetimes: given an input wedge a analogous to a CFT boundary region, the simple wedge z(a) is the largest wedge accessible by a “zigzag,” a certain sequence of antinormal lightsheets. We show that z(a) is a throat, and that it is contained in every other throat. This implies that z(a) is unique; that it is contained in the generalized entanglement wedge; and that it reduces to the AdS prescription as a special case. The zigzag explicitly constructs a preferred Cauchy slice that renders the simple wedge accessible from a; thus it adds a novel structure even in AdS. So far, no spacelike construction is known to reproduce these results, even in time-symmetric settings. This may have implications for the modeling of holographic encoding by tensor networks.

AdS-CFT Correspondence↗
DOE OSTI2022

Numerical algorithms for water waves with background flow over obstacles and topography

Abstract We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

Ambrose, David M.↗
DOE OSTI2025

Efficient proximal subproblem solvers for a nonsmooth trust-region method

In [R. J. Baraldi and D. P. Kouri, Mathematical Programming, (2022), pp. 1-40], we introduced an inexact trust-region algorithm for minimizing the sum of a smooth nonconvex and nonsmooth convex function. The principle expense of this method is in computing a trial iterate that satisfies the so-called fraction of Cauchy decrease condition—a bound that ensures the trial iterate produces sufficient decrease of the subproblem model. In this paper, we expound on various proximal trust-region subproblem solvers that generalize traditional trust-region methods for smooth unconstrained and convex-constrained problems. We introduce a simplified spectral proximal gradient solver, a truncated nonlinear conjugate gradient solver, and a dogleg method. Finally, we compare algorithm performance on examples from data science and PDE-constrained optimization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Multi-objective materials bayesian optimization with active learning of design constraints: Design of ductile refractory multi-principal-element alloys

We report Bayesian Optimization (BO) has emerged as a powerful framework to efficiently explore and exploit materials design spaces. To date, most BO approaches to materials design have focused on the materials discovery problem as if it were a single expensive-to-query ‘black box’ in which the target is to optimize a single objective (i.e., material property or performance metric). Also, such approaches tend to be constraint agnostic. Here, we present a novel multi-information BO framework capable of actively learning materials design as a multiple objectives and constraints problem. We demonstrate this framework by optimally exploring a Refractory Multi-Principal-Element Alloy (MPEA) space, here specifically, the system Mo-Nb-Ti-V-W. The MPEAs are explored to optimize two density-functional theory (DFT) derived ductility indicators (Pugh’s Ratio and Cauchy pressure) while learning design constraints relevant to the manufacturing of high-temperature gas-turbine components. Alloys in the BO Pareto-front are analyzed using DFT to gain an insight into fundamental atomic and electronic underpinning for their superior performance, as evaluated within this framework.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Lagrange-Remap strategy for multi-material fluid-solid simulations using compressive limiters

In the present work, the Lagrange-Remap strategy proposed in [1] is extended to multi-material fluid-solid simulations. Both hypo-elastic and hyper-elastic material models are considered to describe the mechanical behavior of the solids. In practice, the deviatoric stress tensor (for hypo-elastic materials) and the left Cauchy-Green tensor (for isotropic hyper-elastic materials) are remapped, while the use of compressive limiters effectively reduces numerical diffusion during the remapping step. The simplicity of this diffuse interface approach is emphasized in the context of multi-material fluid-solid simulations. A series of Lagrange-Remap test cases, involving both solids and fluids, are conducted and compared with reference Lagrangian simulations, demonstrating the robustness and accuracy of the overall numerical strategy.

Compressive limiters↗
DOE OSTI2021

Anisotropic elastoplastic response of double-porosity media

Here we present a continuum framework to simulate fluid flow through anisotropic elastoplastic media with double porosity. Two effective stress measures σ' and δ" emerge from the thermodynamic formulation, which are energy-conjugate to the elastic and plastic components of strain, respectively. Both effective stress measures can be expressed as a combination of the total Cauchy stress σ and the average pore pressure $\bar{p}$ in the two pore scales. In the effective stress for elasticity, $\bar{p}$ is scaled with a rank-2 Biot tensor, whereas the effective stress for plasticity follows the Terzaghi form in which $\bar{p}$ is scaled by the Kronecker delta. The Biot tensor and storage coefficients are derived as functions of elasticity parameters and porosities. A mixed finite element formulation is introduced to discretize the domain and solve initial boundary-value problems. A stabilization scheme is employed on equal-order interpolation for both displacement and pressure fields throughout the entire range of drainage responses. Numerical simulations reproduce the hydromechanical response of Opalinus shale in one-dimensional consolidation tests throughout the range of primary and secondary consolidation under different external loads. Numerical simulations of the consolidation of a rectangular domain subjected to a strip load demonstrate the efficacy of the proposed stabilization scheme, as well as illustrate the impacts of stress history, mass transfer, and different pore systems on the hydromechanical response.

42 ENGINEERING↗
DOE OSTI2023

A fast Fourier transform-based solver for elastic micropolar composites

This work presents a spectral micromechanical formulation for obtaining the full-field and homogenized response of elastic micropolar composites. The algorithm relies on a coupled set of convolution integral equations for the micropolar strains, where periodic Green’s operators associated with a linear homogeneous reference medium are convolved with functions of the Cauchy and couple stress fields that encode the material’s heterogeneity, as well as any potential material nonlinearity. Such convolution integral equations take an algebraic form in the reciprocal Fourier space that can be solved iteratively. In this vein, the fast Fourier transform (FFT) algorithm is leveraged to accelerate the numerical solution, resulting in a mesh-free formulation in which the periodic unit cell representing the heterogeneous material can be discretized by a regular grid of pixels in two dimensions (or voxels in three dimensions). For verification, the numerical solutions obtained with the micropolar FFT solver are compared with analytical solutions for a matrix with a dilute circular inclusion subjected to plane strain loading. The developed computational framework is then used to study length-scale effects and effective (micropolar) moduli of composites with various topological configurations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

MechElastic: A Python library for analysis of mechanical and elastic properties of bulk and 2D materials

We report the MechElastic Python package evaluates the mechanical and elastic properties of bulk and 2D materials using the elastic coefficient matrix ( C ij ) obtained from any ab-initio density-functional theory (DFT) code. The current version of this package reads the output of VASP, ABINIT, and Quantum Espresso codes (but it can be easily generalized to any other DFT code) and performs the appropriate post-processing of elastic constants as per the requirement of the user. This program can also detect the input structure's crystal symmetry and test the mechanical stability of all crystal classes using the Born-Huang criteria. Various useful material-specific properties such as elastic moduli, longitudinal and transverse elastic wave velocities, Debye temperature, elastic anisotropy, 2D layer modulus, hardness, Pugh's ratio, Cauchy's pressure, Kleinman's parameter, and Lame's coefficients, can be estimated using this program. Another existing feature of this program is to employ the ELATE package (2016) [29] and plot the spatial variation of several elastic properties such as Poisson's ratio, linear compressibility, shear modulus, and Young's modulus in three dimensions. Further, the MechElastic package can plot the equation of state (EOS) curves for energy and pressure for a variety of EOS models such as Murnaghan, Birch, Birch-Murnaghan, and Vinet, by reading the inputted energy/pressure versus volume data obtained via numerical calculations or experiments. This package is particularly useful for the high-throughput analysis of elastic and mechanical properties of materials.

2D materials↗
DOE OSTI2022

Elastic behavior of binary and ternary refractory multi-principal-element alloys

Elastic properties are essential to the mechanical performance of materials, and therefore, can be tuned to design high-performance materials. In this study, the elastic constants of the equiatomic ternary refractory multi-principal-element alloys, NbTiV, and MoNbV, were investigated, using in-situ neutron diffraction and first-principles calculations. The experimentally measured and theoretically predicted elastic constants show a good agreement. The alloying effect of adding Ti and Mo into NbV base alloy on the elastic constants is studied. Particularly, adding Ti elements into the NbV alloy results in the increase of the Zener anisotropy ratio from 0.59 to 0.99, leading to the formation of the elastically isotropic NbTiV alloy, while the addition of Mo decreases the anisotropy ratio to 0.52. Pugh’s ratio (B/G), Cauchy pressure (C 12 -C 44 ), and Poisson’s ratio (ν) are used to predict the brittle/ductile nature of the studied alloys, which is consistent with the mechanical results. The present work provides valuable insights into the design of ductile and strong refractory high-entropy alloys by tuning the elastic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A finite-strain rate- and pressure-dependent constitutive framework for analyzing shock compression behavior of cemented tungsten carbides to 100 GPa

In the present study a thermodynamically-consistent finite-strain rate-and-pressure-dependent constitutive framework is implemented to analyze the shock-compression behavior of cemented tungsten carbides to 100 GPa. Central to this framework is the use of logarithmic strain with a set of invariant basis that allow the Cauchy stress tensor to be expressed as a sum of three response terms that are mutually orthogonal, thus permitting a complete separation of the deviatoric and volumetric (pressure) response. An overstress viscoplasticity model that includes strain and strain rate hardening along with thermal softening is used to represent the deviatoric response, while a complete Mie-Grüneisen equation of state (EoS) is used to obtain the pressure response. Using this formulation, the shock-induced compression behavior of cemented tungsten carbide - obtained from planar plate impact experiments using a 30 mm powder gun to peak stresses of up to ~100 GPa - is analyzed to better understand the structure of the measured shock wave profiles and the associated in-material shock quantities. Of particular interest is the evolution of material inelasticity and strength, and temperature in the tungsten carbide samples during the shock compression process.

Cemented tungsten carbide↗
DOE OSTI2023

Limits of Detection for EXAFS Characterization of Heterogeneous Single-Atom Catalysts

Single-atom catalysts (SACs), consisting of individual metal atoms dispersed on a support, attract attention due to their unique reactivity, efficient use of precious metals, and precise chemical tunability. Characterization of the metal species is crucial to substantiate structure–function relationships. Authors often use—and referees often require—X-ray absorption spectroscopy (XAS) data to prove the absence of clustered metal (or metal oxide) structures after pre-treatment and under in situ or operando conditions. However, there has been no critical assessment of the limitations of XAS in substantiating such conclusive statements, which is particularly important given the potential outsized influence of minority catalyst structures in dictating catalytic activity. Here, in this article, we quantitatively assess the detection limits of XAS to identify metal (or metal oxide) clusters in samples containing predominantly single atoms by modeling the extended X-ray absorption fine structure (EXAFS) of mixtures of structures. We identified that a significant fraction of clusters can coexist with SAC active sites (e.g., ~10% metallic Pt or ~40% oxidized Pt clusters in Pt/CeO 2 SACs), while eluding detection via EXAFS with any statistical significance. To generalize these conclusions, a descriptor-based screening of bulk metal oxides using a continuous Cauchy wavelet transform was proposed that suggests certain materials for which differentiating atomically dispersed metal species and metal oxide clusters would be infeasible by EXAFS (e.g., ReO x ). Based on this analysis, we suggest best practices for the study of SACs using EXAFS and provide recommendations to ensure that conclusions do not outpace the evidence used to support them. In this rapidly expanding research area, rigorous characterization will lead to greater understanding of the behavior of SACs and ultimately improved catalytic materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗