Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “relaxation”

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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Cardiac troponin T N-domain variant destabilizes the actin interface resulting in disturbed myofilament function

Missense variant Ile79Asn in human cardiac troponin T (cTnT-I79N) has been associated with hypertrophic cardiomyopathy and sudden cardiac arrest in juveniles. cTnT-I79N is located in the cTnT N-terminal (TnT1) loop region and is known for its pathological and prognostic relevance. A recent structural study revealed that I79 is part of a hydrophobic interface between the TnT1 loop and actin, which stabilizes the relaxed (OFF) state of the cardiac thin filament. Given the importance of understanding the role of TnT1 loop region in Ca 2+ regulation of the cardiac thin filament along with the underlying mechanisms of cTnT-I79N-linked pathogenesis, we investigated the effects of cTnT-I79N on cardiac myofilament function. Transgenic I79N (Tg-I79N) muscle bundles displayed increased myofilament Ca 2+ sensitivity, smaller myofilament lattice spacing, and slower crossbridge kinetics. These findings can be attributed to destabilization of the cardiac thin filament’s relaxed state resulting in an increased number of crossbridges during Ca 2+ activation. Additionally, in the low Ca 2+ -relaxed state (pCa8), we showed that more myosin heads are in the disordered-relaxed state (DRX) that are more likely to interact with actin in cTnT-I79N muscle bundles. Dysregulation of the myosin super-relaxed state (SRX) and the SRX/DRX equilibrium in cTnT-I79N muscle bundles likely result in increased mobility of myosin heads at pCa8, enhanced actomyosin interactions as evidenced by increased active force at low Ca 2+ , and increased sinusoidal stiffness. These findings point to a mechanism whereby cTnT-I79N weakens the interaction of the TnT1 loop with the actin filament, which in turn destabilizes the relaxed state of the cardiac thin filament.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Ultrafast dynamics in the high-symmetry and in the charge density wave phase of 2 H –NbSe 2

Here we investigate carrier and collective mode dynamics in 2H-NbSe 2 using time-resolved optical pump-probe spectroscopy and compare the results with first-principles calculations. Broadband ultrafast reflectivity studies of 2H-NbSe 2 in a wide temperature interval covering the normal, charge density wave (CDW) and superconducting phase were performed. Spectral features observed in the transient reflectivity experiment were associated with specific optical transitions obtained from band structure calculations. Displacive excitation of coherent phonons showed CDW-associated coherent oscillations of the soft phonon mode across the whole spectral range. Temperature evolution of this coherent phonon mode in the low-excitation linear regime shows softening of the mode down to the CDW transition temperature T CDW with subsequent hardening below T CDW . The global fit of the broadband probe data reveals four different relaxation times associated with characteristic electron-electron, electron-phonon, and phonon-phonon relaxation processes. From first-principles calculations of electron-phonon coupling we associate the few picosecond electron-phonon relaxation time $\tau_2$ with a specific group of phonons with frequencies around 20 meV. On the other hand, the anomalously long relaxation time of $\tau_3 \sim$ to 100 ps is associated with anharmonicity-driven phonon-phonon scattering. All relaxation processes result from anomalies near the second order CDW phase transition that are reflected in the temperature dependencies of the characteristic relaxation times and amplitudes of coherent oscillations. At the highest fluences we observe electronic melting of the CDW and disappearance of the mode hardening below T CDW .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Effect of quasiperiodic and random noise on many-body dynamical decoupling protocols

Symmetries (and their spontaneous rupturing) can be used to protect and engender novel quantum phases and lead to interesting collective phenomena. In the work of K. Agarwal and I. Martin, the authors described a general dynamical decoupling (polyfractal) protocol that can be used to engineer multiple discrete symmetries in many-body systems. This work expands on the former by studying the effect of quasiperiodic and random noise on such a dynamical scheme. We find generally that relaxation of engineered symmetry generators proceeds by (i) an initial relaxation on microscopic timescales to a prethermal plateau whose height is independent of noise, (ii) a linear relaxation regime with a noise-dependent rate, followed by (iii) a slow logarithmic relaxation regime that is only present for quasiperiodic noise. We glean the essential features of these regimes via scaling collapses and show that they can be generally explained by the spectral properties of the various noise waveforms considered. In particular, the quasiperiodic noise is characterized by highly time-dependent spectrum with a noise floor that mimics white noise, and peaks that grow sharper with time. We argue that both the noise floor and peaks contribute to the initial linear-in-time relaxation while the logarithmic regime is initiated when the peaks become sufficiently well resolved and cease to contribute to further relaxation. Here, we provide numerical evidence to justify these findings.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Galvanostatic Intermittent Titration Technique Reinvented: Part I. A Critical Review

The galvanostatic intermittent titration technique (GITT), introduced in 1977 by Weppner and Huggins, provided a readily accessible means to measuring the chemical diffusion coefficient of electrochemical electrode materials. The method continues to be widely used today, but the reported diffusivity values are highly inconsistent, ranging as much as four orders of magnitude for some Li layered oxide compositions. Even qualitative trends of diffusivity are inconsistent, suggesting significant flaws in the implementation of the method. Other variants of the GITT method also suffer from similar inconsistency problems. Here we identify numerous sources of significant error including composition-dependent reaction overpotentials, mathematical flaws in the relaxation analysis methods, finite-size and non-planar geometry effects, inter-particle inhomogeneity issues, early transient effects, and surface area uncertainties. We propose a simple relaxation analysis scheme using the time variable t relax + τ - t relax , where t relax is relaxation time and τ is the galvanostatic pulse duration. We also propose to use dense diffusion-limited samples to isolate the bulk-diffusion process in the time domain. Chemical diffusivity can be extracted much more reliably with this improved implementation of the GITT method.

25 ENERGY STORAGE↗
DOE OSTI2024

A Spectacular Viscoelastic Model Calibration for 828/D230/Alox Generated from Legacy Sandia Data

The nonlinear viscoelastic Spectacular model is calibrated to the thermo-mechanical behavior of 828/D230/Alox with an alox volume fraction of 20 %. Legacy experimental data from Sandia’s polymer properties database (PPD) is used to calibrate the model. Based on known densities of the epoxy 828/D230 and the alox filler, the alox volume fractions listed on the PPD were likely reported incorrectly. The alox volume fractions are recalculated here. Using the recalculated alox volume fractions, the PPD contains experimental data for 828/D230/Alox with alox volume fractions of 16 %, 24 %, and 33 %, so the thermo-mechanical behavior at 20 % alox volume fraction is estimated by interpolating between the bounding cases of of 16 % and 24 %. Because the Spectacularmodel can be fairly challenging to calibrate, the calibration procedure is described in detail. Several of the calibration steps involve inverse parameter identification, where an experiment is simulated and parameters are iteratively updated until the model response matches the experimental data. As the PPD does not fully describe all experimental procedures, the experimental simulations use assumed thermal and mechanical loading rates that are typical for the viscoelastic characterization of epoxies. Spectacular uses four independent relaxation functions related to volumetric (ƒ 1 ), shear (ƒ 2 ), thermal strain (ƒ 3 ), and thermal relaxations (ƒ 4 ). The previous SPEC model form, also known as the universal_polymer model, uses two independent relaxation functions related to volumetric and thermal relaxation (ƒ ν = ƒ 1 = ƒ 3 = ƒ 4 ) and shear relaxation (ƒ s = ƒ 2 ). The two constitutive choices are briefly evaluated here, where it is found that the four relaxation function approach of Spectacular was better suited for fitting the coefficient of thermal expansion during both heating and cooling.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Initial Testing of an In Situ Load Retention Aging Vessel

A thermal aging vessel instrumented with load cells was fabricated. The primary function of the vessel is to continuously monitor the in situ load retention of up to three compressed polymer coupons undergoing thermally accelerated aging under nitrogen. A secondary function is to enable gas sampling of the vessel headspace during thermal aging. Heating of the vessel is achieved using a custom heater jacket. To improve upon our conventional aging study methods which require periodic interruption of aging to perform load testing in an Instron machine at room temperature, this technology aims to automate/facilitate data acquisition/analysis, improve data quality, and enable uninterrupted compression of the polymer which represents the service condition. As an example case to assess functionality of the in situ vessel, the load retention of a siloxane elastomer material additively manufactured by direct-ink-writing (DIW) was measured at three different isothermal aging temperatures for ~1 month. Initial compression of the coupons while near the aging temperature was achieved by temporarily opening the heated vessel to access the interior chamber and manually tightening four nuts to drive the heated compression plate down onto the heated coupons. Initial testing demonstrated achievement of the primary load retention monitoring function. Unfortunately, the vessel leaked which prevented gas sampling; an active purge was used to maintain a nitrogen atmosphere. Welded or otherwise sealed joints, which could be implemented in a future design, would likely eliminate leak paths. To apply time-temperature superposition (TTS), a technique used to provide long-term prediction of the load retention from short-term isothermal data, the load retention needed to be calculated relative to the load at an estimated “equilibrium” time, after most of the transient viscoelastic physical relaxation occurred. The peak load immediately after compression could not be used as the load retention basis for two reasons: (1) age-related changes must be isolated from non-age-related physical relaxation before applying TTS and (2) the manual mechanism used to compress the specimens at the aging temperature was neither smooth nor repeatable which affected the peak load value. To better understand the effect of the mode of initial compression on the measured load, and possibly better estimate “equilibrium” physical relaxation times, systematic stress relaxation experiments were performed using an Instron machine with a thermal chamber. At a given temperature, the DIW polymer was compressed to a fixed strain in either a stepped or continuous manner at two different rates, then held at that strain for 24 hrs. The results indicated that, at a given temperature, the different stress relaxation curves appeared to converge to the same curve at some “equilibrium” time when the non-age-related physical relaxation was mostly complete. Though this observation suggests that the discontinuous manual compression employed by the vessel is feasible, a compression mechanism that is rapid, smooth, and repeatable would enhance its use.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Computer analysis of effects of altering jet fuel properties on refinery costs and yields

This study was undertaken to evaluate the adequacy of future U.S. jet fuel supplies, the potential for large increases in the cost of jet fuel, and to what extent a relaxation in jet fuel properties would remedy these potential problems. The results of the study indicate that refiners should be able to meet jet fuel output requirements in all regions of the country within the current Jet A specifications during the 1990-2010 period. The results also indicate that it will be more difficult to meet Jet A specifications on the West Coast, because the feedstock quality is worse and the required jet fuel yield (jet fuel/crude refined) is higher than in the East. The results show that jet fuel production costs could be reduced by relaxing fuel properties. Potential cost savings in the East (PADDs I-IV) through property relaxation were found to be about 1.3 cents/liter (5 cents/gallon) in January 1, 1981 dollars between 1990 and 2010. However, the savings from property relaxation were all obtained within the range of current Jet A specifications, so there is no financial incentive to relax Jet A fuel specifications in the East. In the West (PADD V) the potential cost savings from lowering fuel quality were considerably greater than in the East. Cost savings from 2.7 to 3.7 cents/liter (10-14 cents/gallon) were found. In contrast to the East, on the West Coast a significant part of the savings was obtained through relaxation of the current Jet A fuel specifications.

Breton, T.↗
NASA NTRS1988

Vibrational-vibrational coupling in air at low humidities

Calculations of sound absorption in air are traditionally based on the assumption that molecular relaxations in N2 and O2 are independent. In binary mixtures of these two gases, however, they are not independent; rather, molecular relaxation is known to be controlled by a very strong vibrational-vibrational (V-V) coupling, which influences both the relaxation frequencies and the relaxation strengths. This article shows that small concentrations of the air constituents CO2 and H2O, which themselves possess a strong V-V coupling to N2 and O2, serve to decouple the N2 and O2 relaxations. To characterize the N2-O2 coupling a coupling strength is derived which depends upon the constituent concentrations and the related reaction rate constants. It is found that the molecular relaxations associated with N2 and O2 in air experience a gradual transition from strong to weak coupling as the humidity increases beyond approximately 0.001 mole percent.

Zuckerwar, Allan J.↗
NASA NTRS1991

The improved robustness of multigrid elliptic solvers based on multiple semicoarsened grids

Multigrid convergence rates degenerate on problems with stretched grids or anisotropic operators, unless one uses line or plane relaxation. For 3-D problems, only plane relaxation suffices, in general. While line and plane relaxation algorithms are efficient on sequential machines, they are quite awkward and inefficient on parallel machines. A new multigrid algorithm is presented based on the use of multiple coarse grids, that eliminates the need for line or plane relaxation in anisotropic problems. This algorithm was developed and the standard multigrid theory was extended to establish rapid convergence for this class of algorithms. The new algorithm uses only point relaxation, allowing easy and efficient parallel implementation, yet achieves robustness and convergence rates comparable to line and plane relaxation multigrid algorithms. The algorithm described is a variant of Mulder's multigrid algorithm for hyperbolic problems. The latter uses multiple coarse grids to achieve robustness, but is unsuitable for elliptic problems, since its V-cycle convergence rate goes to one as the number of levels increases. The new algorithm combines the contributions from the multiple coarse grid via a local switch, based on the strength of the discrete operator in each coordinate direction.

Naik, Naomi H.↗
NASA NTRS2002

Structural Fluctuations and Thermophysical Properties of Molten II-VI Compounds

The objectives of the project are to conduct ground-based experimental and theoretical research on the structural fluctuations and thermophysical properties of molten II-VI compounds to enhance the basic understanding of the existing flight experiments in microgravity materials science programs as well as to study the fundamental heterophase fluctuation phenomena in these melts by: 1) conducting neutron scattering analysis and measuring quantitatively the relevant thermophysical properties of the II-VI melts (such as viscosity, electrical conductivity, thermal diffusivity and density) as well as the relaxation characteristics of these properties to advance the understanding of the structural properties and the relaxation phenomena in these melts and 2) performing theoretical analyses on the melt systems to interpret the experimental results. All the facilities required for the experimental measurements have been procured, installed and tested. Thermal diffusivity of molten tellurium has been measured by a laser flash method in the temperature range of 500 C to 900 C. The measured diffusivity as a function of temperature agrees fairly well with published data. However, a relaxation phenomenon, which shows a slow drift of the measured thermal conductivity toward the equilibrium value after cooling of the melt, was observed for the first time. An apparatus based on the transient torque induced by a rotating magnetic field has been developed to determine the viscosity and electrical conductivity of semiconducting liquids. Viscosity measurements on molten tellurium showed a similar relaxation behavior to the measured diffusivity. The density and volume expansion coefficients for pure Te and HgTe melts were measured as a function of temperature using a pycnometric method. A density maximum was found for both melts but no relaxation behavior was observed. Neutron scattering experiments were performed on the HgTe and HgZnTe melts and the results on pair distribution showed better resolution than previously reported. A simple kinetic theory, which contains the formation reactions of Te polymerization that increases the molecular weight, was proposed to interpret the relaxation behavior of HgZnTe viscosity data.

Su, Ching-Hua↗
NASA NTRS2006

Application of p-Multigrid to Discontinuous Galerkin Formulations of the Poisson Equation

We investigate p-multigrid as a solution method for several different discontinuous Galerkin (DG) formulations of the Poisson equation. Different combinations of relaxation schemes and basis sets have been combined with the DG formulations to find the best performing combination. The damping factors of the schemes have been determined using Fourier analysis for both one and two-dimensional problems. One important finding is that when using DG formulations, the standard approach of forming the coarse p matrices separately for each level of multigrid is often unstable. To ensure stability the coarse p matrices must be constructed from the fine grid matrices using algebraic multigrid techniques. Of the relaxation schemes, we find that the combination of Jacobi relaxation with the spectral element basis is fairly effective. The results using this combination are p sensitive in both one and two dimensions, but reasonable convergence rates can still be achieved for moderate values of p and isotropic meshes. A competitive alternative is a block Gauss-Seidel relaxation. This actually out performs a more expensive line relaxation when the mesh is isotropic. When the mesh becomes highly anisotropic, the implicit line method and the Gauss-Seidel implicit line method are the only effective schemes. Adding the Gauss-Seidel terms to the implicit line method gives a significant improvement over the line relaxation method.

Helenbrook, B. T.↗
NASA NTRS2016

Cryovolcanic Emplacement of Domes on Europa

Here we explore the hypothesis that certain domes on Europa may have been produced by the extrusion of viscous cryolavas. A new mathematical method for the emplacement and relaxation of viscous lava domes is presented and applied to putative cryovolcanic domes on Europa. A similarity solution approach is applied to the governing equation for fluid flow in a cylindrical geometry, and dome relaxation is explored assuming a volume of cryolava has been rapidly emplaced onto the surface. Nonphysical sin- gularities inherent in previous models for dome relaxation have been eliminated, and cryolava cooling is represented by a time-variable viscosity. We find that at the onset of relaxation, bulk kinematic viscosities may lie in the range between 10(exp 3) and 10(exp 6) sq m/s, while the actual fluid lava viscosity may be much lower. Plausible relaxation times to form the domes, which are linked to bulk cryolava rheology, are found to range from 3.6 days to 7.5 years. We find that cooling of the cryolava, while dominated by conduction through an icy skin, should not prevent fluids from advancing and relaxing to form domes within the timescales considered. Determining the range of emplacement conditions for putative cryolava domes will shed light on Europa's resurfacing history. In addition, the rheologies and compositions of erupted cryolavas have implications for subsurface cryomagma ascent and local surface stress conditions on Europa.

volcanism↗
NASA NTRS2016

New Approaches to Inferences for Steep-Sided Domes on Venus

New mathematical approaches for the relaxation and emplacement of viscous lava domes are presented and applied to steep-sided domes on Venus. A similarity solution approach is applied to the governing equation for fluid flow in a cylindrical geometry for two distinct scenarios. In the first scenario, dome relaxation is explored assuming a constant volume of fluid (i.e. lava) has been rapidly emplaced onto the surface. Cooling of lava is represented by a time-variable viscosity and singularities inherent in previous models for dome relaxation have been eliminated. At the onset of relaxation, bulk dynamic viscosities lie in the range between 10(exp 10) - 10(exp 16) Pa s, consistent with basaltic-andesite to rhyolitic compositions. Plausible relaxation times range from 5 to 5000 years, depending on initial lava viscosity. The first scenario, however, is only valid during the final stages of dome relaxation and does not consider the time taken for lava to be extruded onto the surface. In the second scenario, emplacement and growth of a steep-sided dome is considered when the volume of lava on the surface increases over time (i.e. time-variable volume approach). The volumetric flow rate may depend on an arbitrary power of the dome thickness, thus embracing Newtonian as well as other rheologies for describing terrestrial and planetary mass flows. The approach can be used to distinguish between basic flow rate models for fluid emplacement. The formalism results in radial expansion of a dome proportional to t(sup 1/2), consistent with the diffusive nature of the governing equation. The flow at the front is shown to thicken as the front advances for a constant rate of lava supply. Emplacement times are intimately correlated with the bulk rheology. Comparison of the theoretical profiles with the shape of a typical dome on Venus indicates that a Newtonian bulk rheology is most appropriate, consistent with prior studies. However, results here suggest a bulk dynamic viscosity of 10(exp 12) - 10(exp 13) Pa s and emplacement times of approximately 2-16 years. Both scenarios investigated give emplacement times significantly less than prior estimates and compositions consistent with basaltic andesite.

steep-sided domes↗
DOE OSTI2022

Photoluminescence of Cis-Polyacetylene Semiconductor Material

Photoluminescence (PL) is one of the key experimental characterizations of optoelectronic materials, including conjugated polymers (CPs). In this study, a simplified model of an undoped cis-polyacetylene (cis-PA) oligomer was selected and used to explain the mechanism of photoluminescence (PL) of the CPs. Using a combination of the ab initio electronic structure and a time-dependent density matrix methodology, the photo-induced time-dependent excited state dynamics were computed. We explored the phonon-induced relaxation of the photoexcited state for a single oligomer of cis-PA. Here, the dissipative Redfield equation of the motion was used to compute the dissipative excited state dynamics of electronic degrees of freedom. This equation used the nonadiabatic couplings as parameters. The computed excited state dynamics showed that the relaxation rate of the electron is faster than the relaxation rate of the hole. The dissipative excited-state dynamics were combined with radiative recombination channels to predict the PL spectrum. The simulated results showed that the absorption and emission spectra both have a similar transition. The main result is that the computed PL spectrum demonstrates two mechanisms of light emission originating from (i) the inter-band transitions, corresponding to the same range of transition energies as the absorption spectrum and (ii) intra-band transitions not available in the absorption spectra. However, the dissipative Redfield equation of the motion was used to compute the electronic degrees of freedom of the nonadiabatic couplings, which helped to process the time propagation of the excited dynamic state. This excited dynamic state shows that the relaxation rate of the electron is faster than the relaxation rate of the hole, which can be used for improving organic semiconductor materials for photovoltaic and LED applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

McCormick envelopes in mixed-integer PDE-constrained optimization

McCormick envelopes are a standard tool for deriving convex relaxations of optimization problems that involve polynomial terms. Such McCormick relaxations provide lower bounds, for example, in branch-and-bound procedures for mixed-integer nonlinear programs but have not gained much attention in PDE-constrained optimization so far. This lack of attention may be due to the distributed nature of such problems, which on the one hand leads to infinitely many linear constraints (generally state constraints that may be difficult to handle) in addition to the state equation for a pointwise formulation of the McCormick envelopes and renders bound-tightening procedures that successively improve the resulting convex relaxations computationally intractable. We analyze McCormick envelopes for a model problem class that is governed by a semilinear PDE involving a bilinearity and integrality constraints. We approximate the nonlinearity and in turn the McCormick envelopes by averaging the involved terms over the cells of a partition of the computational domain on which the PDE is defined. This yields convex relaxations that underestimate the original problem up to an a priori error estimate that depends on the mesh size of the discretization. These approximate McCormick relaxations can be improved by means of an optimization-based bound-tightening procedure. We show that their minimizers converge to minimizers to a limit problem with a pointwise formulation of the McCormick envelopes when driving the mesh size to zero. We provide a computational example, for which we certify all of our imposed assumptions. The results point to both the potential of the methodology and the gaps in the research that need to be closed. Our methodology provides a framework first for obtaining pointwise underestimators for nonconvexities and second for approximating them with finitely many linear inequalities in an infinite-dimensional setting.

Approximations and Expansions↗
DOE OSTI2022

Molecular simulations and NMR reveal how lipid fluctuations affect membrane mechanics

Lipid bilayers form the main matrix of functional cell membranes, and their dynamics underlie a host of physical and biological processes. Here we show that elastic membrane properties and collective molecular dynamics (MD) are related by the mean-square amplitudes (order parameters) and relaxation rates (correlation times) of lipid acyl chain motions. We performed all-atom MD simulations of liquid-crystalline bilayers that allow direct comparison with carbon-hydrogen (CH) bond relaxations measured with NMR spectroscopy. Previous computational and theoretical approaches have assumed isotropic relaxation, which yields inaccurate description of lipid chain dynamics and incorrect data interpretation. Instead, the new framework includes a fixed bilayer normal (director axis) and restricted anisotropic motion of the CH bonds in accord with their segmental order parameters, enabling robust validation of lipid force fields. Simulated spectral densities of thermally excited CH bond fluctuations exhibited well-defined spin-lattice (Zeeman) relaxations analogous to those in NMR measurements. Their frequency signature could be fit to a simple power-law function, indicative of nematic-like collective dynamics. Moreover, calculated relaxation rates scaled as the squared order parameters yielding an apparent K C modulus for bilayer bending. Our results show a strong correlation with K C values obtained from solid-state NMR studies of bilayers without and with cholesterol as validated by neutron spin-echo measurements of membrane elasticity. The simulations uncover a critical role of interleaflet coupling in membrane mechanics and thus provide important insights into molecular sites of emerging elastic properties within lipid bilayers.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Effect of plasticity on the experimental characterization of the creep property of a cement paste: Spherical vs. conical microindentation

Microindentation is today a powerful tool for rapidly characterizing the elastic modulus and the logarithmic basic creep of a cement paste. However, a theoretical work proved that a microindentation test with a sharp indenter shape involves stress concentration and delayed plasticity which can cause an overestimation of the estimated creep property. This work aims at experimentally investigating the effect of the shape indenter on the identification of the visco-elastic property of a cement paste by microindentation. First, several microindentation grids were carried out on a cement paste with both a Berkovich indenter and a spherical indenter at different penetration depths or load levels for load control (creep) tests or displacement control (relaxation) tests, respectively. As for a visco-elastic material the duality principle between creep and relaxation holds, a numerical method was developed for predicting the relaxation response from the assumed logarithmic creep compliance. The present results show that only spherical microindentation satisfy the duality principle for which creep and relaxation response can be predicted by a unique logarithmic creep compliance function. Finally, spherical microindentation creep tests at low load as well as spherical microindentation relaxation tests provide a greater indentation modulus M and contact creep modulus C{sub v} as likely less affected by possible plasticity effects.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A machine learning approach for increased throughput of density functional theory substitutional alloy studies

In this study, a machine learning-based technique is developed to reduce the computational cost required to explore large design spaces of substitutional alloys. The first advancement is based on a neural network (NN) approach to predict the initial position of minority and majority ions prior to DFT relaxation. The second advancement is to allow the NN to predict the total energy for every possibility minority ion position and select the most stable configuration in the absence of relaxing each trial minority configuration. A bismuth oxide materials system, (Bi x La y Yb z ) 2 MoO 6 , was used as the model system to demonstrate the developed methods and quantify the resulting computational speedup. Compared to a brute force method that requires the calculation of every permutation of minority configuration and subsequent DFT relaxation, a 1.3× speedup was realized if the NN predicted the initial configuration of ions prior to relaxation. Implementation of the second advancement allowed the NN to predict the total energy for all possible trial configurations and downselect the most stable configurations prior to relaxation, resulting in a speedup of approximately 37×. Validation was done by comparing position and energy between the NN and DFT predictions. A maximum position vector mean squared error (MSE) of 1.6 × 10 -2 and a maximum energy MSE of 2.3 × 10 -7 was predicted for the worst case configuration. This method demonstrates a significant computational speedup, which has the potential for even greater computational savings for larger compositional design spaces.

36 MATERIALS SCIENCE↗