Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Strain calculation”

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,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Pressure tuning of competing interactions on a honeycomb lattice

Exchange interactions are mediated via orbital overlaps across chemical bonds. Thus, modifying the bond angles by physical pressure or strain can tune the relative strength of competing interactions. Here we present a remarkable case of such tuning between the Heisenberg (J) and Kitaev (K) exchange, which respectively establish magnetically ordered and spin liquid phases on a honeycomb lattice. We observe a rapid suppression of the Néel temperature (T N ) with pressure in Ag 3 LiRh 2 O 6 , a spin-1/2 honeycomb lattice with both J and K couplings. Using a combined analysis of x-ray data and first-principles calculations, we find that pressure modifies the bond angles in a way that increases the ∣K/J∣ ratio and thereby suppresses T N . Consistent with this picture, we observe a spontaneous onset of muon spin relaxation (μSR) oscillations below T N at low pressure, whereas in the high pressure phase, oscillations appear only when T < T N /2. Unlike other candidate Kitaev materials, Ag 3 LiRh 2 O 6 is tuned toward a quantum critical point by pressure while avoiding a structural dimerization in the relevant pressure range.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Assessing the source of error in the Thomas–Fermi–von Weizsäcker density functional

We investigate the source of error in the Thomas–Fermi–von Weizsäcker (TFW) density functional relative to Kohn–Sham density functional theory (DFT). In particular, through numerical studies on a range of materials, for a variety of crystal structures subject to strain and atomic displacements, we find that while the ground state electron density in TFW orbital-free DFT is close to the Kohn–Sham density, the corresponding energy deviates significantly from the Kohn–Sham value. Here, we show that these differences are a consequence of the poor representation of the linear response within the TFW approximation for the electronic kinetic energy, confirming conjectures in the literature. In so doing, we find that the energy computed from a non-self-consistent Kohn–Sham calculation using the TFW electronic ground state density is in very good agreement with that obtained from the fully self-consistent Kohn–Sham solution.

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

Compressive strain turns 𝑠 ± - into 𝑑-wave pairing in a one-unit-cell La 3 ⁢Ni 2 ⁢O 7 thin film via substrate-induced hole doping

Motivated by recent reports of ambient-pressure superconductivity in La 3 ⁢Ni 2⁢ O 7 films grown on LaSrAlO 4 , we investigate the superconducting instability in a one-unit-cell (1UC) thin film using ab initio and random-phase approximation techniques. Compared to the high-pressure bulk system, the ratio of interlayer 𝑑 3⁢𝑧 2 −𝑟 2 hopping to intralayer 𝑑 𝑥 2 −𝑦 2 hopping is suppressed in the 1UC thin film, and the crystal-field splitting of the 𝑒 𝑔 orbitals is increased. Here, our calculation indicates that spin-fluctuation-driven pairing correlations are weak for the stoichiometric case at ambient pressure, but increase significantly under hole doping. The leading pairing symmetry is also found to change by hole doping. Specifically, we obtain a leading 𝑑 𝑥 2 −𝑦 2 pairing state at moderate hole doping, followed by a 𝑑 𝑥⁢𝑦 state at higher doping. These states are driven by intraband spin-fluctuation scattering within the 𝛾 hole pocket centered around the 𝑀 point, and arise primarily from states in the Ni layer farther from the substrate. These results strongly suggest that the thin-film superconducting samples are hole-doped and that pairing in this system predominantly arises in the layer, as opposed to the interlayer pairing in the pressurized bulk system.

Zhang, Yang [Oak Ridge National Laboratory (ORNL),↗
NASA NTRS1978

User's guide to computer programs JET 5A and CIVM-JET 5B to calculate the large elastic-plastic dynamically-induced deformations of multilayer partial and/or complete structural rings

These structural ring deflections lie essentially in one plane and, hence, are called two-dimensional (2-d). The structural rings may be complete or partial; the former may be regarded as representing a fragment containment ring while the latter may be viewed as a 2-d fragment-deflector structure. These two types of rings may be either free or supported in various ways (pinned-fixed, locally clamped, elastic-foundation supported, mounting-bracket supported, etc.). The initial geometry of each ring may be circular or arbitrarily curved; uniform-thickness or variable-thickness rings may be analyzed. Strain-hardening and strain-rate effects of initially-isotropic material are taken into account. An approximate analysis utilizing kinetic energy and momentum conservation relations is used to predict the after-impact velocities of each fragment and of the impact-affected region of the ring; this procedure is termed the collision-imparted velocity method (CIVM) and is used in the CIVM-JET 5 B program. This imparted-velocity information is used in conjunction with a finite-element structural response computation code to predict the transient, large-deflection, elastic-plastic responses of the ring. Similarly, the equations of motion of each fragment are solved in small steps in time. Provisions are made in the CIVM-JET 5B code to analyze structural ring response to impact attack by from 1 to 3 fragments, each with its own size, mass, translational velocity components, and rotational velocity. The effects of friction between each fragment and the impacted ring are included.

Wu, R. W. H.↗
NASA NTRS1988

Analysis of crack propagation in roller bearings using the boundary integral equation method - A mixed-mode loading problem

Crack propagation in a rotating inner raceway of a high-speed roller bearing is analyzed using the boundary integral method. The model consists of an edge plate under plane strain condition upon which varying Hertzian stress fields are superimposed. A multidomain boundary integral equation using quadratic elements was written to determine the stress intensity factors KI and KII at the crack tip for various roller positions. The multidomain formulation allows the two faces of the crack to be modeled in two different subregions, making it possible to analyze crack closure when the roller is positioned on or close to the crack line. KI and KII stress intensity factors along any direction were computed. These calculations permit determination of crack growth direction along which the average KI times the alternating KI is maximum.

Ghosn, L. J.↗
DOE OSTI2026

Dirac Fermions and Flat Bands in Phosphorus Carbide Nanotubes: Structural and Quantum Phase Transitions in a Quasi-One-Dimensional Material

Chemically realistic quasi-one-dimensional (1D) materials in which Dirac Fermions and highly degenerate flat bands coexist intrinsically at the Fermi level are exceedingly rare, while representing a highly desirable platform for correlated and topological quantum phenomena. Here, in this work, using specialized symmetry-adapted first-principles calculations we predict a new class of nanomaterials─phosphorus carbide nanotubes (P 2 C 3 NTs)─obtained by rolling monolayer P 2 C 3 , a two-dimensional material shown in a previous letter to host “double Kagome bands”. Both armchair and zigzag P 2 C 3 NTs are stable at room temperature and feature the rare coexistence of Dirac crossings and multiple flat bands at the Fermi level inherited from the underlying honeycomb–Kagome lattice, with the flat bands resilient to elastic deformations. Under large strain, the structure transforms from honeycomb–Kagome to “brick-wall”, accompanied by multiple coupled structural and quantum phase transitions. We also uncover localized edge states, spin splitting from vacancies and dopants, and strain-tunable magnetism. Together, these results establish P 2 C 3 NTs as a chemically specific and mechanically tunable 1D material platform with potential applications in quantum hardware and spintronics.

carbon nanotubes↗
DOE OSTI2022

High‐Pressure and High‐Temperature Single‐Crystal Elasticity of Cr‐Pyrope: Implications for the Density and Seismic Velocity of Subcontinental Lithospheric Mantle

Abstract Single‐crystal X‐ray diffraction and Brillouin spectroscopy experiments were performed on a natural Cr‐pyrope (Prp 71.0 Alm 12.6 Sps 0.7 Grs 3.5 Uvr 12.2 ) at high pressure and high temperature up to 11.0 GPa and 800 K. Fitting the collected data to the third‐order finite strain equation yields bulk modulus ( K S 0 ), shear modulus ( G 0 ), their pressure ((∂ K S /∂ P ) T and (∂ G /∂ P ) T ) and temperature (( ∂K S /∂T ) P and ( ∂G/∂T ) P ) derivatives, K S 0 = 167.7(8) GPa, G 0 = 91.5(5) GPa, (∂ K S /∂ P ) T = 4.3(1), (∂ G /∂ P ) T = 1.4(1), ( ∂K S /∂T ) P = − 0.0175(1) GPa/K and ( ∂G/∂T ) P = − 0.0073(1) GPa/K. Using the obtained results, we examined whether the elastic properties of the Cr‐pyrope can be accurately calculated from those of endmembers including pyrope, almandine, grossular, and uvarovite assuming a linear relationship between elastic properties and composition (end‐member model). The results indicate that the end‐member model provides a sufficient approximation for the elastic properties of Cr‐pyrope in calculating the density and velocity of the subcontinental lithospheric mantle (SCLM). We modeled the densities and velocities of three typical types of SCLM (Archon, Proton, and Tecton) in order to investigate how the variation of chemical composition influences the SCLM. We obtained that the compositional change from the Archon to the Tecton increases the density of the SCLM significantly, which can be an important prerequisite for SCLM delamination. However, the compositional variation only slightly changes the velocity of the SCLM and the change is within the uncertainty of the calculation. Moreover, in comparison to the velocity, ρ / V P and ρ / V S are much more sensitive to the compositional change of the SCLM.

Xu, Jingui↗
DOE OSTI2026

Modulating charge transport via 2 MeV He + irradiation in VO 2

Vanadium dioxide (VO 2 ) is of interest for adaptive electronic applications such as neuromorphic neuristor devices and variable emissivity or tunable thermal control materials, thanks to its key property—a metal–insulator transition (MIT) at 68 °C that is accompanied by a dramatic change in electrical and optical properties. To improve performance in these roles, it is critical to develop approaches to engineer transport properties and the MIT behavior. While many documented techniques exist to modulate the MIT and film resistivities via lattice strain and chemical doping, less is known about the effects of ion irradiation on the intrinsic properties of VO 2 , despite the ability to control the spatial distribution of irradiation beams and the prevalence of high energy ion implantation in the semiconductor industry. The impact of irradiation of different acceleration energies on the responses of VO 2 is of specific interest, as charged particle energy generally impacts both the resulting defect profile and corresponding transport behavior. Here, we demonstrate that 2 MeV He ions at equivalent calculated displacements per atom, in two different types of films, can create remarkable changes to the nature of charge transport in VO 2 , especially in the low-temperature insulating phase. Simulation of resulting changes in electrical conductivity reveals that He ion irradiation offers a strategy to increase both oscillation frequency and the signal transmission. These results provide insights into the intentional design of defect populations to modulate transport for neuromorphic VO 2 devices.

Gurrola, Rebeca M. [Texas A & M Univ., College Sta↗
DOE OSTI2025

Experimental determination of energy-level structures of diamond silicon-vacancy centers in off-axis magnetic fields

Here, we report the development of an experimental approach to efficiently determine the energy level structure of an individual silicon vacancy (SiV) center in a magnetic field along an arbitrary direction. This approach uses two coupling rates (one each for the ground and the excited states) to characterize the combined effects of static strain and dynamic Jahn-Teller coupling and exploits the fact that orbital Zeeman effects vanish when the magnetic field is normal to the SiV axis. With an analytical expression for the energy level structure of the SiV under a transverse magnetic field, the two coupling rates can be directly derived from two measurements: one on the frequency separation between two spin-conserved transitions and the other on the coherent population trapping resonance of the SiV ground spin states. A detailed comparison between the numerical calculation and the experimental result on the dependence of the spin-conserved splitting on both the amplitude and direction of the magnetic field further reveals unequal orbital magnetic coupling for the ground and excited states, indicating that unequal orbital quenching factors are needed for an accurate description of the SiV energy level structure in a magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1978

Global crustal response model

A computer program was developed for calculating the radial displacement due to the body tide as a function of spatial position and time. The positions of the moon and sun were evaluated by means of the Hill-Brown and Newcomb theories, respectively. Green's functions were utilized to simulate the elastic Earth responses due to point loads, a convolution program has been developed to evaluate the elastic Earth responses due to ocean loading by a given ocean tide constituent. These responses include radial and horizontal displacement, tilt and gravity acceleration, and strain tensor components.

Source record↗
DOE OSTI2025

Internally Catalyzed Hydrogen Atom Transfer (I-CHAT)—A New Class of Reactions in Combustion Chemistry

The current paradigm of low-T combustion and autoignition of hydrocarbons is based on the sequential two-step oxygenation of fuel radicals. The key chain-branching occurs when the second oxygenation adduct (OOQOOH) is isomerized releasing an OH radical and a key ketohydroperoxide (KHP) intermediate. The subsequent homolytic dissociation of relatively weak O–O bonds in KHP generates two more radicals in the oxidation chain leading to ignition. Based on the recently introduced intramolecular “catalytic hydrogen atom transfer” mechanism (J. Phys. Chem. 2024, 128, 2169), abbreviated here as I-CHAT, we have identified a novel unimolecular decomposition channel for KHPs to form their classical isomers—enol hydroperoxides (EHP). The uncertainty in the contribution of enols is typically due to the high computed barriers for conventional (“direct”) keto–enol tautomerization. Remarkably, the I-CHAT dramatically reduces such barriers. The novel mechanism can be regarded as an intramolecular version of the intermolecular relay transfer of H-atoms mediated by an external molecule following the general classification of such processes (Catal. Rev.-Sci. Eng. 2014, 56, 403). Here, we present a detailed mechanistic and kinetic analysis of the I-CHAT-facilitated pathways applied to n-hexane, n-heptane, and n-pentane models as prototype molecules for gasoline, diesel, and hybrid rocket fuels. We particularly examined the formation kinetics and subsequent dissociation of the γ-enol-hydroperoxide isomer of the most abundant pentane-derived isomer γ-C5-KHP observed experimentally. To gain molecular-level insight into the I-CHAT catalysis, we have also explored the role of the internal catalyst moieties using truncated models. All applied models demonstrated a significant reduction in the isomerization barriers, primarily due to the decreased ring strain in transition states. In addition, the longer-range and sequential H-migration processes were also identified and illustrated via a combined double keto–enol conversion of heptane-2,6-diketo-4-hydroperoxide as a potential chain-branching model. To assess the possible impact of the I-CHAT channels on global fuel combustion characteristics, we performed a detailed kinetic analysis of the isomerization and decomposition of γ-C5-KHP comparing I-CHAT with key alternative reactions—direct dissociation and Korcek channels. Calculated rate parameters were implemented into a modified version of the n-pentane kinetic model developed earlier using RMG automated model generation tools (ACS Omega, 2023, 8, 4908). Simulations of ignition delay times revealed the significant effect of the new pathways, suggesting an important role of the I-CHAT pathways in the low-T combustion of large alkanes.

Biochemistry & Molecular Biology↗
DOE OSTI2022

First-Principles Study of Ferroelectric and Optical Properties in Derivatives of Thiourea

Organic ferroelectric materials in the form of molecular crystals are promising alternatives to well-known ferroelectric perovskite oxides due to their structural flexibility, tunability, and ease of processing. First-principles density functional theory calculations combined with group theory analyses are employed in this work to discuss polar nature of a class of derivatives of thiourea: CH 4 N 2 A where A = O, S, Se, and Te. Further, optical properties in the UV range along with thermodynamic and dynamic stability of such materials are also evaluated in the study. We find that CH 4 N 2 Te with bandgap energy of ~3.33 eV shows large spontaneous polarization of P = 10.92 μC/cm 2 . The electronic and optical properties become tunable under strain.

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

Aculeaxanthones A–E, new xanthones from the marine-derived fungus Aspergillus aculeatinus WHUF0198

Introduction Dimeric natural products are widespread in plants and microorganisms, which usually have complex structures and exhibit greater bioactivities than their corresponding monomers. In this study, we report five new dimeric tetrahydroxanthones, aculeaxanthones A−E ( 4−8 ), along with the homodimeric tetrahydroxanthone secalonic acid D ( 1 ), chrysoxanthones B and C ( 2 and 3 ), and 4−4’-secalonic acid D ( 9 ), from different fermentation batches of the title fungus. Methods A part of the culture was added to a total of 60 flasks containing 300 ml each of number II fungus liquid medium and culture 4 weeks in a static state at 28˚C. The liquid phase (18 L) and mycelia was separated from the fungal culture by filtering. A crude extract was obtained from the mycelia by ultrasound using acetone. To obtain a dry extract (18 g), the liquid phase combined with the crude extract were further extracted by EtOAc and concentrated in vacuo. The MIC of anaerobic bacteria was examined by a broth microdilution assay. To obtain MICs for aerobic bacteria, the agar dilution streak method recommended in Clinical and Laboratory Standards Institute document (CLSI) M07-A10 was used. Compounds 1−9 was tested against the Bel-7402, A-549 and HCT-116 cell lines according to MTT assay. Results and Discussion The structures of these compounds were elucidated on the base of 1D and 2D NMR and HR-ESIMS data, and the absolute configurations of the new xanthones 4−8 were determined by conformational analysis and time-dependent density functional theory-electronic circular dichroism (TDDFT-ECD) calculations. Compounds 1–9 were tested for cytotoxicity against the Bel-7402, A549, and HCT-116 cancer cell lines. Of the dimeric tetrahydroxanthone derivatives, only compound 6 provided cytotoxicity effect against Bel-7402 cell line (IC50, 1.96 µM). Additionally, antimicrobial activity was evaluated for all dimeric tetrahydroxanthones, including four Gram-positive bacteria including Enterococcus faecium ATCC 19434, Bacillus subtilis 168, Staphylococcus aureus ATCC 25923 and MRSA USA300; four Gram-negative bacteria, including Helicobacter pylori 129, G27, as well as 26,695, and multi drug-resistant strain H. pylori 159, and one Mycobacterium M. smegmatis ATCC 607. However, only compound 1 performed activities against H. pylori G27, H. pylori 26695, H. pylori 129, H. pylori 159, S. aureus USA300, and B. subtilis 168 with MIC values of 4.0, 4.0, 2.0, 2.0, 2.0 and 1.0 μg/mL, respectively.

Wu, Jun↗
NASA NTRS1983

Specimen size and geometry effects on fracture toughness of Al2O3 measured with short rod and short bar chevron-notch specimens

Plane strain fracture toughness measurements were made on Al2O3 using short rod and short bar chevron notch specimens previously calibrated by the authors for their dimensionless stress intensity factor coefficients. The measured toughness varied systematically with variations in specimen size, proportions, and chevron notch angle apparently due to their influence on the amount of crack extension to maximum load (the measurement point). The toughness variations are explained in terms of a suspected rising R curve for the material tested, along with a discussion of an unavoidable imprecision in the calculation of K sub Ic for materials with rising R curves when tested with chevron notch specimens.

Shannon, J. L., Jr.↗
NASA NTRS1982

A time-dependent ice sheet model - Preliminary results

A numerical model of ice sheet flow is developed, and preliminary results are described. This model includes vertical resolution of temperature, stress, and strain rate which represents a considerable improvement over previous vertically averaged ice sheet models. The model follows the flow of ice along a flow line within an ice sheet drainage basin. Longitudinal stresses and basal sliding are included. Basal sliding is dependent on the base shear stress and a specified distribution of basal water pressure. The numerical methods used to solve the coupled set of stress and velocity equations for the static and time-evolutionary cases are discussed. A steady state profile simulating an ice stream is calculated for a particular set of input parameters, and changes in the profile are examined for different choices of parameters. Preliminary studies of response behavior are completed using a simplified ice sheet geometry with a fixed terminus or grounding line. The results of these studies illustrate ice sheet thinning in response to a lowered sea level or to a reduction in the extent of ice rises (or pinning points) within ice shelves.

Bindschadler, R. A.↗
NASA NTRS1990

A real time neural net estimator of fatigue life

A neural net architecture is proposed to estimate, in real-time, the fatigue life of mechanical components, as part of the Intelligent Control System for Reusable Rocket Engines. Arbitrary component loading values were used as input to train a two hidden-layer feedforward neural net to estimate component fatigue damage. The ability of the net to learn, based on a local strain approach, the mapping between load sequence and fatigue damage has been demonstrated for a uniaxial specimen. Because of its demonstrated performance, the neural computation may be extended to complex cases where the loads are biaxial or triaxial, and the geometry of the component is complex (e.g., turbopump blades). The generality of the approach is such that load/damage mappings can be directly extracted from experimental data without requiring any knowledge of the stress/strain profile of the component. In addition, the parallel network architecture allows real-time life calculations even for high frequency vibrations. Owing to its distributed nature, the neural implementation will be robust and reliable, enabling its use in hostile environments such as rocket engines. This neural net estimator of fatigue life is seen as the enabling technology to achieve component life prognosis, and therefore would be an important part of life extending control for reusable rocket engines.

Troudet, T.↗
NASA NTRS1990

A real time neural net estimator of fatigue life

A neural network architecture is proposed to estimate, in real-time, the fatigue life of mechanical components, as part of the intelligent Control System for Reusable Rocket Engines. Arbitrary component loading values were used as input to train a two hidden-layer feedforward neural net to estimate component fatigue damage. The ability of the net to learn, based on a local strain approach, the mapping between load sequence and fatigue damage has been demonstrated for a uniaxial specimen. Because of its demonstrated performance, the neural computation may be extended to complex cases where the loads are biaxial or triaxial, and the geometry of the component is complex (e.g., turbopumps blades). The generality of the approach is such that load/damage mappings can be directly extracted from experimental data without requiring any knowledge of the stress/strain profile of the component. In addition, the parallel network architecture allows real-time life calculations even for high-frequency vibrations. Owing to its distributed nature, the neural implementation will be robust and reliable, enabling its use in hostile environments such as rocket engines.

Troudet, T.↗