SEARCH · Engineering Papers

Results for “coupling constant”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,585 records · Page 88

All sources NASA NTRS DOE OSTI

DOE OSTI2024

First-principles investigation of elastic, vibrational, and thermodynamic properties of kagome metals CsM 3 Te 5 (M = Ti, Zr, Hf)

Kagome metals are a unique class of quantum materials characterized by their distinct atomic lattice arrangement, featuring interlocking triangles and expansive hexagonal voids. These lattice structures impart exotic properties, including superconductivity, interaction-driven topological many-body phenomena, and magnetism, among others. The kagome metal CsM 3 ⁢Te 5 (where M = Ti, Zr, or Hf) exhibits both superconductivity and nontrivial topological electronic properties, offering a promising platform for exploring topological superconductivity. This study employs first-principles density functional theory calculations to systematically analyze the elastic, mechanical, vibrational, thermodynamic, and electronic properties of CsM 3 ⁢Te 5 (M = Ti, Zr, Hf). Our calculations reveal that the studied compounds—CsTi 3 ⁢Te 5 , CsZr 3 ⁢Te 5 , and CsHf 3 ⁢Te 5 —are ductile metals with elastic properties akin to the hexagonal Bi and Sb, with average elastic constants, including a bulk modulus of 27 GPa, a shear modulus of 11 GPa, and Young's modulus of 29 GPa. We observe peculiar dispersionless, flat, phonon branches in the vibrational spectra of these metals. Additionally, we thoroughly analyze the symmetries of the zone-center phonon eigenvectors and predict vibrational fingerprints of the Raman- and infrared-active phonon modes. The analysis of thermodynamic properties reveals the Einstein temperature for CsTi 3 ⁢Te 5 , CsZr 3 ⁢Te 5 , and CsHf 3 ⁢Te 5 to be 66, 54, and 53 K, respectively. Our orbital-decomposed electronic structure calculations reveal significant in-plane steric interactions and multiple Dirac band crossings near the Fermi level. We further investigate the role of spin-orbit coupling effect on the studied properties. Furthermore, this theoretical investigation sheds light on the intriguing quantum behavior of kagome metals.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Using Strain in the Vadose Zone During Sinusoidal-Rate Well Tests to Characterize an Aquifer

Well tests with sinusoidal variation in pumping rate are appealing because of the ability to filter competing signals and they also can be set up with no net removal of water, which is an advantage at contaminated sites. Current methods of interpreting sinusoidal-rate well tests analyze the phase lag and amplitude of pressure signals in the aquifer. This approach is restricted by the availability of monitoring wells. An alternative approach could be to measure the stain in the vadose zone instead of pressure in the aquifer. The goal of this research is to evaluate the feasibility of using signals from in-situ sensors in the vadose zone to characterize sinusoidal-rate well tests. Well tests with a variety of periods and amplitudes were performed using a well completed in an unconfined saprolite aquifer with a water table at 9m depth at our test site near Clemson, South Carolina. Vertical strain was measured every meter between 1 to 7 m using a distributed optical fiber sensing technique called coherence length gated microwave photonics interferometry (CMPI). A novel packaging system was designed to couple the optical fiber to the saprolite. The results show that strain measured in the vadose zone during sinusoidal-rate pumping varies periodically with amplitudes from 0.1 to 3 microstrain, depending on the amplitude of the sinusoidal pumping rate. Hydraulic heads at similar radial distances were also periodic with changing amplitudes ranging up to ±0.2 meters. The strain in the vadose zone responds to pressure in the underlying aquifer—the strain is tensile when the pressure is falling, and compressive when the pressure is rising. The phase lag between the pumping rate and the strain was 3 min with a pumping period of 15 minutes, and this gives a hydraulic diffusivity of 0.06 m2/s. This is the same diffusivity obtained by analyzing a constant rate pumping test, and it is slightly larger than the diffusivity (0.044 m2/s) obtained by analyzing the phase lag of the pressure in a nearby monitoring well. Strain measurements in the vadose zone could be made with much less cost and effort than installing a deep monitoring well. These results suggest that strain measurements in the vadose zone could be used to supplement monitoring wells, thereby increasing the resolution, and reducing costs of analyzing well tests.

well test, sinusoidal, periodic pumping, strain↗

DOE OSTI2024

Liquid Piston with Spray Cooling Near-Isothermal Compressor

The goal of this project was to prototype and characterize the performance of a liquid-piston spray-cooled gas compressor. The working principle of the compressor enables optimized high-efficiency operation over a very wide range of operating conditions, unlike conventional compressors that are optimized for a narrow range of operating conditions. The compressor technology is suitable for many applications, such as gas pipeline transport, gas storage, and commercial and residential heat pumps. Both physical testing and computational fluid dynamics (CFD) modeling of the processes using the Oak Ridge National Laboratory high-performance computing center were completed. The experimental and CFD studies focused on a near-isothermal liquid-piston compressor (LPC) that uses propylene glycol to compress CO 2 . The first prototype demonstrated isothermal operation during several sequentially executed cycles of CO 2 compression and raised the temperature of the compressed CO 2 by only 2 K, compared with approximately 6 K when the gas was compressed non-isothermally. Isothermal operation was demonstrated at CO 2 flow rates of up to 2 L/min. The second prototype was designed with two compression chambers to allow continuous flow of high-pressure CO 2 . However, the design of the valve train to direct flow between the compression chambers was not sufficient to allow demonstration of CO 2 compression. Numerical simulations of the LPC in which the compression chamber was filled with propylene glycol injected from the bottom inlet were performed using large eddy simulation (LES) with the wall-adapting local eddy-viscosity subgrid-scale model coupled with the multiphase volume of fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfers within the compression chamber. In this effort, the effects of boundary conditions applied to the LES-VOF calculations (i.e., no wall, an adiabatic wall, and a wall with a heat flux subscribed) on the overall pressure and temperature of the CO 2 gas as well as the transient evolution of flow and heat transfer within the compression chamber were investigated and are discussed in this report. The LES calculation with no wall showed no dynamical flow patterns, and the volume-averaged temperature of CO 2 increased from 305 to 392.7 K, whereas LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO 2 temperatures. The results of the LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO 2 gas layers and circulation flow patterns within the chamber that contributed to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor. Though an industrial partnership for commercializing the compressor was not secured, the technology was attractive for an industrial partner to use in two research proposals in response to US Department of Energy funding opportunity announcements.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1993

Accurate Quartic Force Fields and Vibrational Frequencies for HCN and HNC

The quartic force fields of HCN and HNC are determined using atomic natural orbital one-particle basis sets of spdf/spd and spdfg/spdf quality in conjunction with the CCSD(T) electron correlation method (singles and doubles coupled-cluster theory plus a perturbation estimate of the effects of connected triple excitations). The HCN force field is in good agreement with a recent experimentally derived force field and also with the force field recently computed by Wong and Bacskay. On the basis of the good agreement obtained for HCN, it is argued that the ab initio quartic force field for HNC is superior to a prior force field derived from experiment. The harmonic frequencies of HNC are predicted to be 3822 +/- 10, 472 +/- 5, and 2051 +/- 10 cm(exp -1) for omega(sub 1), omega(sub 2), and omega(sub 3), respectively; the experimentally derived values are above these values and fall outside the estimated uncertainties. Using the quartic force field, spectroscopic constants are predicted for HNC based on a vibrational second-order perturbation theory analysis. It is also asserted that the gas-phase fundamental nu(sub 3) for HNC is slightly lower than the matrix isolation value. The range of validity of the quartic force fields is investigated by comparison of variational vibrational energies computed with the quartic force fields to those obtained from our recently reported global HCN/HNC potential energy surface and also to experimental data.

Lee, Timothy J.↗

NASA NTRS2012

Accurate ab initio Quartic Force Fields of Cyclic and Bent HC2N Isomers

Highly correlated ab initio quartic force field (QFFs) are used to calculate the equilibrium structures and predict the spectroscopic parameters of three HC2N isomers. Specifically, the ground state quasilinear triplet and the lowest cyclic and bent singlet isomers are included in the present study. Extensive treatment of correlation effects were included using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T). Dunning s correlation-consistent basis sets cc-pVXZ, X=3,4,5, were used, and a three-point formula for extrapolation to the one-particle basis set limit was used. Core-correlation and scalar relativistic corrections were also included to yield highly accurate QFFs. The QFFs were used together with second-order perturbation theory (with proper treatment of Fermi resonances) and variational methods to solve the nuclear Schr dinger equation. The quasilinear nature of the triplet isomer is problematic, and it is concluded that a QFF is not adequate to describe properly all of the fundamental vibrational frequencies and spectroscopic constants (though some constants not dependent on the bending motion are well reproduced by perturbation theory). On the other hand, this procedure (a QFF together with either perturbation theory or variational methods) leads to highly accurate fundamental vibrational frequencies and spectroscopic constants for the cyclic and bent singlet isomers of HC2N. All three isomers possess significant dipole moments, 3.05D, 3.06D, and 1.71D, for the quasilinear triplet, the cyclic singlet, and the bent singlet isomers, respectively. It is concluded that the spectroscopic constants determined for the cyclic and bent singlet isomers are the most accurate available, and it is hoped that these will be useful in the interpretation of high-resolution astronomical observations or laboratory experiments.

Inostroza, Natalia↗

NASA NTRS1993

Ultra low-noise charge coupled device

Special purpose CCD designed for ultra low-noise imaging and spectroscopy applications that require subelectron read noise floors, wherein a non-destructive output circuit operating near its 1/f noise regime is clocked in a special manner to read a single pixel multiple times. Off-chip electronics average the multiple values, reducing the random noise by the square-root of the number of samples taken. Noise floors below 0.5 electrons rms are possible in this manner. In a preferred embodiment of the invention, a three-phase CCD horizontal register is used to bring a pixel charge packet to an input gate adjacent a floating gate amplifier. The charge is then repeatedly clocked back and forth between the input gate and the floating gate. Each time the charge is injected into the potential well of the floating gate, it is sensed non-destructively. The floating gate amplifier is provided with a reference voltage of a fixed value and a pre-charge gate for resetting the amplifier between charge samples to a constant gain. After the charge is repeatedly sampled a selected number of times, it is transferred by means of output gates, back into the horizontal register, where it is clocked in a conventional manner to a diffusion MOSFET amplifier. It can then be either sampled (destructively) one more time or otherwise discarded.

Janesick, James R.↗

DOE OSTI2023

Electrochemically induced fracture in LLZO: How the interplay between flaw density and electrostatic potential affects operability

Fracture and short circuit in the Li 7 La 3 Zr 2 O 12 (LLZO) solid electrolyte are two key issues that prevent its adoption in battery cells. In this paper, we utilize phase-field simulations that couple electrochemistry and fracture to evaluate the maximum electric potential that LLZO electrolytes can support as a function of crack density. In the case of a single crack, we find that the applied potential at the onset of crack propagation exhibits inverse square root scaling with respect to crack length, analogous to classical fracture mechanics. Here, we further find that the short-circuit potential scales linearly with crack length. In the realistic case where the solid electrolyte contains multiple cracks, we reveal that failure fits the Weibull model. The failure distributions shift to favor failure at lower overpotentials as areal crack density increases. Furthermore, when flawless interfacial buffers are placed between the applied potential and the bulk of the electrolyte, failure is mitigated. When constant currents are applied, current focuses in near-surface flaws, leading to crack propagation and short circuit. We find that buffered samples sustain larger currents without reaching unstable overpotentials and without failing. Our findings suggest several mitigation strategies for improving the ability of LLZO to support larger currents and improve operability.

25 ENERGY STORAGE↗

DOE OSTI2023

Late Holocene Cliff Retreat in Del Mar, CA, Revealed From Shore Platform 10 Be Concentrations and Numerical Modeling

Rocky coast cliff retreat presents a hazard to coastal communities and infrastructure that is potentially amplified under rising sea level conditions, among other factors. Unfortunately, constraints on retreat rates are typically limited to those derived from imagery and maps spanning the last ~100 years. Here, we use a newly developed coupled model of shore platform profile evolution and cosmogenic radionuclide production that considers the influence of relative sea level (RSL) rise, weathering, material resistance, and wave height decay to model cliff retreat over millennial timescales and its potential drivers in Del Mar, California. We demonstrate the ability to use topographic and bathymetric measurements from a narrow shore platform along with a limited data set of nine cosmogenic 10Be concentrations extending ~125 m from the cliff base to obtain modeled cliff retreat rates that steadily range from 5.0 to 12.5 cm yr -1 over the last two millennia until 100 years before present. These rates are consistent with modern retreat rates of about 2–19 cm yr -1 here. RSL rise in Southern California remained relatively constant during the late Holocene, potentially explaining the relatively stable modeled cliff retreat rate over this time. We also explore the relative influence of weathering, material resistance, and wave erosion efficacy and find that both weathering and wave-driven erosion are necessary to replicate the measured data at this location, with the latter exerting a stronger control on model acceptance, suggesting that waves may provide a possible mechanism by which RSL rise may influence coastal cliff erosion in southern California.

58 GEOSCIENCES↗

NASA NTRS1994

Sensitivity of lag-damping correlations to structural and aerodynamic approximations of isolated experimental rotors in forward flight

The predictions of regressive lag-mode damping levels are correlated with the database of an isolated, soft-inplane, three-blade rotor operated untrimmed. The database was generated at the Army Aeroflightdynamics Directorate at Ames. The correlation covers a broad range of data, from near-zero thrust conditions in hover to high-thrust and highly stalled conditions in forward flight with advance ratio as high as 0.55 and shaft angle as high as 20 degrees. In the experimental rotor, the airfoil or blade portion has essentially uniform mass and stiffness distributions, but the root flexure has highly nonuniform mass and stiffness distributions. Accordingly, the structural approximations refer to four models of root-flexure-blade assembly. They range from a rigid flap-lag model to three elastic flap-lag-torsion models, which differ in modeling the root flexure. The three models of root-flexure are: three root springs in which the bending-torsion couplings are fully accounted for; a finite-length beam element with some average mass and stiffness distributions such that the fundamental frequencies match those of the experimental model; and accurate modal representation in which the actual mass and stiffness distributions of the experimental root-flexure-blade assembly are used in calculating the nonrotating mode shapes. The four models of root-flexure-blade assembly are referred to as the rigid flap-lag model, spring model, modified model and modal model. For each of these four models of the root-flexure-blade assembly, the predictions are based on the following five aerodynamic theories: ear theory, which accounts for large angle-of-attack and reverse-flow effects on lift, and has constant drag and pitching moment; quasisteady stall theory, which includes quasisteady stall lift, drag and pitching moment characteristics of the airfoil section, dynamics stall theory, which uses the ONERA dynamic stall models of lift, drag and pitching moment; dynamic wake theory, which is based on a finite-state three-dimensional wake model and includes all wake effects including both shed and trailing vorticity; and dynamics and wake theory, which combines both dynamic stall theory and dynamic wake theory and is a relatively complete aerodynamic representation.

Gaonkar, G. H.↗

NASA NTRSDate not supplied

Parametric Mechanism Design Through Numerical Optimization and Physics Simulation

Design-Build-Test approaches for developing spaceflight hardware are prohibitively time and cost intensive and often lead to suboptimal mechanism designs. Approaches that couple machine learning and high-fidelity physics simulation could eliminate the need for hardware prototyping and dramatically accelerate the engineering design cycle, ultimately reducing cost. This work presents a modular NASA-developed toolchain to optimize hardware mechanisms in a virtual environment using numerical optimization and multi-body physics simulation. The toolchain enables multi-objective optimization, generates parametric CAD files that can be further post-processed by an end user, and can be expanded to optimize full systems and non-mechanical parameters such as feedback control variables. We demonstrate the toolchain through an independently verifiable design problem that optimizes wheel radius to achieve a desired linear velocity in a rigid-body physics environment when the wheel rotates at a constant angular speed, and then post-process the parametric CAD file of the optimal design generated by the tool before ultimately manufacturing it via 3D printing. We end with a discussion of how the toolchain can incorporate other analysis tools, including finite element analysis, computational fluid dynamics, and granular media simulations.

Optimization↗

NASA NTRSDate not supplied

An Optimization-Based Toolchain for Parametric Mechanism Design

Optimization↗

DOE OSTI2025

Deciphering Reaction Mechanisms of Molecular Proton Reduction Catalysts with Cyclic Voltammetry: Kinetic vs Thermodynamic Control

The kinetics and thermodynamics of elementary reaction steps involved in the catalytic reduction of protons to hydrogen define the reaction landscape for catalysis. The mechanisms can differ in the order of the elementary proton transfer, electron transfer, and bond-forming steps and can be further differentiated by the sites at which protons and electrons localize. Access to fully elucidated mechanistic, kinetic, and thermochemical details of molecular catalysts is crucial to facilitate the development of new catalysts that operate with optimal efficiency, selectivity, and durability. The mechanism by which a catalyst operates, as well as the kinetics and thermodynamics associated with the individual steps, can often be accessed through electroanalytical studies. Here, this Account details the application of cyclic voltammetry to interrogate reaction mechanisms and quantify the kinetics and thermodynamics of elementary reaction steps for a series of molecular catalysts that mediate electrochemical proton reduction. I distinguish the limiting scenarios wherein a catalyst operates under kinetic control vs thermodynamic control, with a focus on detecting how cyclic voltammetry features shift with proton source strength and concentration, as well as scan rate. For systems that operate under kinetic control, catalytic currents are observed at, or slightly positive toward, the formal potential for the redox process that triggers catalysis. Under thermodynamic control, catalytic responses shift as a function of the proton source pKa and effective pH of the solution. After drawing this distinction, we introduce the appropriate voltammetry experiments and accompanying analytical expressions for extracting key metrics from the data. To illustrate analytical strategies to quantify elementary reaction steps of catalysts operating under kinetic control, I describe our studies of proton reduction catalysts Co(dmgBF 2 ) 2 (CH 3 CN) 2 (dmgBF 2 = difluoroboryl-dimethylglyoxime) and [Ni(P 2 Ph N 2 Ph ) 2 ] 2+ (P 2 Ph N 2 Ph = 1,5-phenyl-3,7-phenyl-1,5-diaza-3,7-diphosphacyclooctane). Here, peak shift analysis, foot-of-the-wave analysis, and plateau current analysis are applied to data sets wherein voltammetric response are recorded as a function of catalyst concentration, proton source concentration, proton source strength, and scan rate to quantify rate constants for elementary proton transfer and bond-forming steps in a catalytic cycle. Further, the case study of [Ni(P 2 Ph N 2 Ph ) 2 ] 2+ illustrates how complementary spectroscopic methods can bolster the mechanistic assignment. Collectively, these two studies showcase how detailed mechanistic studies inform on rate-limiting elementary steps in catalysis and other key processes underpinning catalysis. Second, I present analytical strategies to interrogate catalysts operating under thermodynamic control, centered on the case study of [Ni II (P 2 Ph N 2 Bn ) 2 ] 2+ (P 2 Ph N 2 Bn = 1,5-dibenzyl-3,7-diphenyl-1,5-diaza-3,7-diphosphacyclooctane). Here, the application of nonaqueous Pourbaix theory to extract thermodynamic information is introduced, and the construction of a coupled Pourbaix diagram is detailed. This study identifies ligand-based protonation as the key process that places catalysis under thermodynamic control and influences the reaction mechanism. Together, the work detailed in this Account showcases the utility of electroanalytical methods to disentangle complex reaction mechanisms and extract key thermochemical and kinetic parameters for elementary steps of catalysis. Through detailed presentation of the key analytical expressions that underpin these analyses, this Account seeks to facilitate the adoption of cyclic voltammetry by the community to fully extract kinetic, thermochemical, and mechanistic information on electrochemical small-molecule activation.

catalysts↗

DOE OSTI2026

Organic Cation Influence on Organic–Inorganic Thermal Equilibration within 2D Metal Halide Perovskites

Two-dimensional inorganic–organic lead halide perovskites exhibit tunable optoelectronic properties that are dictated by alternating layers of metal halide octahedra and organic cations. In such hybrid materials, vibrational coupling and thermalization between the low mass, insulating, organic cation spacers, and high mass inorganic octahedra remains poorly understood. Here, in this study, using femtosecond infrared-pump electronic-probe (IPEP) spectroscopy, we investigate the role of organic cation identity regarding the kinetics of vibrational energy exchange and sublattice mechanical coupling. Linear aliphatic cation-containing 2D perovskites with 4, 6, or 8 carbons (butylammonium = BA, hexylammonium = HA, and octylammonium = OA) were produced as thin films and mid-infrared pump pulses and then selectively excited organic cation stretch vibrations. Thermal energy introduction was evaluated via visible-wavelength optical probing that conveys inorganic octahedra response, owing to changes in electronic absorption. Multiple distinct spectral shifts appear upon excitation of the organic cations. Following an initial optical Stark shift caused by pump–probe temporal overlap, a redshift occurs within the first 10 ps that we attribute to compression of the octahedra by the expanded organic layers. Subsequently, a blueshift of the bandgap occurs commensurate with equilibration of both sublattices at an elevated temperature with a time constant of 21.8 to 31.6 ps depending on the spacer identity, with thermal energy transfer slowing by nearly 50% for the longest linker. By examining the effect that altered cation identity has on vibrational energy exchange, this work begins to offer routes to tune nonequilibrium response as well as provides insight into the fundamental impacts of organic spacers on thermal energy exchange in perovskite materials.

Peifer, Shoshanna E. [Northwestern Univ., Evanston↗

NASA NTRS2007

Effect of Electron Beam Freeform Fabrication (EBF3) Processing Parameters on Composition of Ti-6-4

The Electron Beam Freeform Fabrication (EBF3) process developed at NASA Langley Research Center was evaluated using a design of experiments approach to determine the effect of processing parameters on the composition and geometry of Ti-6-4 deposits. The effects of three processing parameters: beam power, translation speed, and wire feed rate, were investigated by varying one while keeping the remaining parameters constant. A three-factorial, three-level, fully balanced mutually orthogonal array (L27) design of experiments approach was used to examine the effects of low, medium, and high settings for the processing parameters on the chemistry, geometry, and quality of the resulting deposits. Single bead high deposits were fabricated and evaluated for 27 experimental conditions. Loss of aluminum in Ti-6-4 was observed in EBF3 processing due to selective vaporization of the aluminum from the sustained molten pool in the vacuum environment; therefore, the chemistries of the deposits were measured and compared with the composition of the initial wire and base plate to determine if the loss of aluminum could be minimized through careful selection of processing parameters. The influence of processing parameters and coupling between these parameters on bulk composition, measured by Direct Current Plasma (DCP), local microchemistries determined by Wavelength Dispersive Spectrometry (WDS), and deposit geometry will also be discussed.

Lach, Cynthia L.↗

NASA NTRS2006

Recent Results from Epitaxial Growth on Step Free 4H-SiC Mesas

This paper updates recent progress made in growth, characterization, and understanding of high quality homoepitaxial and heteroepitaxial films grown on step-free 4H-SiC mesas. First, we report initial achievement of step-free 4H-SiC surfaces with carbon-face surface polarity. Next, we will describe further observations of how step-free 4H-SiC thin lateral cantilever evolution is significantly impacted by crystal faceting behavior that imposes non-uniform film thickness on cantilever undersides. Finally, recent investigations of in-plane lattice constant mismatch strain relief mechanisms observed for heteroepitaxial growth of 3C-SiC as well as 2H-AlN/GaN heterofilms on step-free 4H-SiC mesas will be reviewed. In both cases, the complete elimination of atomic heterointerface steps on the mesa structure enables uniquely well-ordered misfit dislocation arrays to form near the heterointerfaces with remarkable lack of dislocations threading vertically into the heteroepilayers. In the case of 3C-SiC heterofilms, it has been proposed that dislocation half-loops nucleate at mesa edges and glide laterally along the step-free 3C/4H interfaces. In contrast, 3C-SiC and 2H-AlN/GaN heterofilms grown on 4H-SiC mesas with steps exhibit highly disordered interface misfit dislocation structure coupled with 100X greater density of dislocations threading through the thickness of the heteroepilayers. These results indicate that the presence of steps at the heteroepitaxial interface (i.e., on the initial heteroepitaxial nucleation surface) plays a highly important role in the defect structure, quality, and relaxation mechanisms of single-crystal heteroepitaxial films.

Neudeck, Philip G.↗

DOE OSTI2021

Error Analysis on Numerical Integration Algorithms in a Hypoelasticity Framework

This report determines local truncation errors for common stress integration algorithms used in explicit finite element codes with hypoelastic material models. The hypoelastic integration algorithms in question utilize an operator splitting procedure in a rotation neutralized configuration, where the stress response is determined from de- coupling the total deformation into rotational and strain dependent components. This document analyzes the error in evolving the stress given a one-step time increment Δt and compares the errors associated with both the rotational and strain components of the operator splitting method. A slight modification to a traditional algorithm is proposed and studied, where the rate of deformation is appropriately rotated from the midstep configuration at t n+1/2 to the end step configuration at t n+1 before the constitutive evaluation. The proposed modification either completely eliminates the error associated with the rotation rate or is of the same order of magnitude as the original algorithm for the three test cases considered in this report. These cases consist of an unaxial stretch with a constant true strain rate with a rigid body rotation, an uniaxial stretch with a constant engineering strain rate with a rigid body rotation, and a simple shear deformation. All three cases are compared to a closed form solution, and in almost every test case the alternative algorithm yields the most accurate one-step local truncation error.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Fixed Bit Rotary Drilling Failure Criteria Effect on Drilling Vibration

Polycrystalline diamond compact (PDC) bit use has been expanded to geothermal drilling applications due to their improved performance in drilling compared to roller cone bits. Multiple vibration models have been developed yet and there are limited use of field data to study the ability of numerical vibration model to predict recorded vibration in the field. The objective of this paper is to evaluate the capabilities of two different rock-bit interaction models with two different ROP formulations for predicting the actual behavior seen in the field. The first model (Model-I) considers the rock-bit interaction, incorporating the effect of subsurface formation characteristics. The second model (Model-II) is a PDC bit model that uses calculated threshold weight-on-cutter (WOCt) to study the interaction of the bit with the rock surface. The models, which couple the drill stem axial and torsional motions, are solved numerically using the fourth-order Runge-Kutta method and compared with recorded downhole drilling data from the Chocolate Mountain Aerial Gunnery Range (CMAGR) geothermal well located in California. An 8–½" PDC bit with 8 blades was modeled and vibration simulation was performed to investigate stick-slip vibrations at 1400 ft. depth. Simulation results using field data indicate that Model-II encounters a low range of stick-slip severity (SSV) as compared to Model-I, which is in the moderate range. Comparison with the field value for SSV shows Model-II is the best representative of the field case for this depth. The sensitivity analysis shows that results are sensitive to model constants making prediction of SSV from models uncertain.

15 GEOTHERMAL ENERGY↗

DOE OSTI2020

Machine-learned interatomic potentials by active learning: amorphous and liquid hafnium dioxide

We propose an active learning scheme for automatically sampling a minimum number of uncorrelated configurations for fitting the Gaussian Approximation Potential (GAP). Our active learning scheme consists of an unsupervised machine learning (ML) scheme coupled with a Bayesian optimization technique that evaluates the GAP model. We apply this scheme to a Hafnium dioxide (HfO 2 ) dataset generated from a "melt-quench" ab initio molecular dynamics (AIMD) protocol. Our results show that the active learning scheme, with no prior knowledge of the dataset, is able to extract a configuration that reaches the required energy fit tolerance. Further, molecular dynamics (MD) simulations performed using this active learned GAP model on 6144 atom systems of amorphous and liquid state elucidate the structural properties of HfO 2 with near ab initio precision and quench rates (i.e., 1.0 K/ps) not accessible via AIMD. The melt and amorphous X-ray structural factors generated from our simulation are in good agreement with experiment. In addition, the calculated diffusion constants are in good agreement with previous ab initio studies.

36 MATERIALS SCIENCE↗