Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time factors”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

The Combined Influence of Molecular Weight and Temperature on the Aging and Viscoelastic Response of a Glassy Thermoplastic Polyimide

The effect of molecular weight on the viscoelastic performance of an advanced polymer (LaRC-SI) was investigated through the use of creep compliance tests. Testing consisted of short-term isothermal creep and recovery with the creep segments performed under constant load. The tests were conducted at three temperatures below the glass transition temperature of five materials of different molecular weight. Through the use of time-aging-time superposition procedures, the material constants, material master curves and aging-related parameters were evaluated at each temperature for a given molecular weight. The time-temperature superposition technique helped to describe the effect of temperature on the timescale of the viscoelastic response of each molecular weight. It was shown that the low molecular weight materials have higher creep compliance and creep rate, and are more sensitive to temperature than the high molecular weight materials. Furthermore, a critical molecular weight transition was observed to occur at a weight-average molecular weight of M (bar) (sub w) 25000 g/mol below which, the temperature sensitivity of the time-temperature superposition shift factor increases rapidly. The short-term creep compliance data were used in association with Struik's effective time theory to predict the long-term creep compliance behavior for the different molecular weights. At long timescales, physical aging serves to significantly decrease the creep compliance and creep rate of all the materials tested.

Nicholson, Lee M.↗
NASA NTRS1980

Numerical methods for large-scale time-dependent partial differential equations

The factors that influence the choice of a numerical solution scheme for large-scale time-dependent partial differential equations are examined with emphasis on the treatment of boundaries. Consideration is given to standard boundary treatments, moving boundaries, and artificial boundaries simulating an infinite domain. Implicit methods and applications to specific problems are discussed with special attention to shocks and steady state solutions.

Turkel, E.↗
DOE OSTI2023

Time-temperature-stress superposition of PBX 9502 compressive creep data for lifetime predictions

The plastic-bonded explosive (PBX) 9502 is a viscoelastic, high-solids loaded polymer bound composite comprised of 95 weight% (wt%) TATB explosive crystals and 5 wt% FK-800 polymer binder. The uniaxial quasi-static mechanical properties have been studied and characterized extensively over 25+ years, including creep. Creep is the strain evolution over time as a result of an asymmetric load typically at temperatures above ambient. Creep measurements to failure can take excessively long periods of time depending on various factors including applied stress, temperature, and loading rate that render a single test to failure impossible. Using the principles of time-temperature-stress superposition (TTSSP), short-term creep tests at various applied stresses and temperatures can be used for long-term creep predictions. The TTSSP model shifts short-term test data using two factors, one for temperature and stress, creating a master curve at a reference test condition. The master curve can be shifted to predict the creep response at any desired test condition. PBX 9502 compression specimens were tested at three temperatures and various applied stresses. This work discusses the TTSSP analysis of PBX 9502 compressive creep data to create a single comprehensive master curve.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Time-Temperature-Stress Superposition of PBX 9502 Compressive Creep Data for Lifetime Predictions

The plastic-bonded explosive (PBX) 9502 is a viscoelastic, high-solids loaded polymer bound composite comprised of 95 weight% (wt%) TATB explosive crystals and 5 wt% FK-800 polymer binder. The uniaxial quasi-static mechanical properties have been studied and characterized extensively over 25+ years, including creep. Creep is the strain evolution over time as a result of an asymmetric load typically at temperatures above ambient. Creep measurements to failure can take excessively long periods of time depending on various factors including applied stress, temperature, and loading rate that render a single test to failure impossible. Using the principles of time-temperature-stress superposition (TTSSP), short-term creep tests at various applied stresses and temperatures can be used for long-term creep predictions. The TTSSP model shifts short-term test data using two factors, one for temperature and stress, creating a master curve at a reference test condition. The master curve can be shifted to predict the creep response at any desired test condition. PBX 9502 compression specimens were tested at three temperatures and various applied stresses. This work discusses the TTSSP analysis of PBX 9502 compressive creep data to create a single comprehensive master curve.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Nonlinear dielectric relaxation of polar liquids

Molecular dynamics of two water models, SPC/E and TIP3P, at a number of temperatures are used to test the Kivelson-Madden equation connecting single-particle and collective dielectric relaxation times through the Kirkwood factor. The relation is confirmed by simulations and used to estimate the nonlinear effect of the electric field on the dielectric relaxation time. We show that the main effect of the field comes through slowing down of the single-particle rotational dynamics and the relative contribution of the field-induced alteration of the Kirkwood factor is insignificant for water. Theories of nonlinear dielectric relaxation need to mostly account for the effect of the field on rotations of a single dipole in a polar liquid.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

A novel potential/viscous flow coupling technique for computing helicopter flow fields

The primary objective of this work was to demonstrate the feasibility of a new potential/viscous flow coupling procedure for reducing computational effort while maintaining solution accuracy. This closed-loop, overlapped velocity-coupling concept has been developed in a new two-dimensional code, ZAP2D (Zonal Aerodynamics Program - 2D), a three-dimensional code for wing analysis, ZAP3D (Zonal Aerodynamics Program - 3D), and a three-dimensional code for isolated helicopter rotors in hover, ZAPR3D (Zonal Aerodynamics Program for Rotors - 3D). Comparisons with large domain ARC3D solutions and with experimental data for a NACA 0012 airfoil have shown that the required domain size can be reduced to a few tenths of a percent chord for the low Mach and low angle of attack cases and to less than 2-5 chords for the high Mach and high angle of attack cases while maintaining solution accuracies to within a few percent. This represents CPU time reductions by a factor of 2-4 compared with ARC2D. The current ZAP3D calculation for a rectangular plan-form wing of aspect ratio 5 with an outer domain radius of about 1.2 chords represents a speed-up in CPU time over the ARC3D large domain calculation by about a factor of 2.5 while maintaining solution accuracies to within a few percent. A ZAPR3D simulation for a two-bladed rotor in hover with a reduced grid domain of about two chord lengths was able to capture the wake effects and compared accurately with the experimental pressure data. Further development is required in order to substantiate the promise of computational improvements due to the ZAPR3D coupling concept.

Summa, J. Michael↗
DOE OSTI2020

A hybrid porous model for full reactor core scale CFD investigation of a prismatic HTGR

Three-dimensional (3-D) Computational Fluid Dynamics (CFD) analysis of a whole nuclear reactor core is a tremendous challenge due to the large geometric volume and complex structures. Here, this research presents a hybrid porous (HP) model to simplify a prismatic High Temperature Gas-cooled Reactor (HTGR) core, so 3-D CFD investigation can be performed on a full reactor core scale. In the HP model, the prototypic small coolant channels in the nuclear fuel blocks are lumped together to form multiple equivalent large coolant channels, and then the porous medium flow model is applied to each of them. Therefore, heat transfer in fuel blocks is computed by a hybrid combination of solid energy and porous flow energy equations. The similarity between the HP model and prototypic model is achieved by deriving the porous flow permeability, inertial resistance factor, and artificial thermophysical properties. Compared with the widely used whole porous (WP) flow model, the HP model preserves more realistic geometric structures, and therefore more accurate physical processes. The General Atomics' Modular High Temperature Gas-cooled Reactor (MHTGR) design was chosen as a prototype to demonstrate the methodology. Simulations were performed using the prototypic CFD model and HP model at steady-state forced circulation, steady-state natural circulation, and transient conditions that correspond to normal operation, extended period of pressurized cool down, and short-term transients after reactor shutdown, respectively. The comparison shows good agreement between the HP model and prototypic model in the maximum fuel temperature, average solid temperature, and helium flow rate, which demonstrates the potential applicability of the HP model for a full reactor core scale simulation in the future. As a benefit, the HP model reduces the mesh quantity by a factor of 50 from a prototypic model. Correspondingly, the computation time was reduced by a factor of at least 30.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Distributed deep reinforcement learning for simulation control

Abstract Several applications in the scientific simulation of physical systems can be formulated as control/optimization problems. The computational models for such systems generally contain hyperparameters, which control solution fidelity and computational expense. The tuning of these parameters is non-trivial and the general approach is to manually ‘spot-check’ for good combinations. This is because optimal hyperparameter configuration search becomes intractable when the parameter space is large and when they may vary dynamically. To address this issue, we present a framework based on deep reinforcement learning (RL) to train a deep neural network agent that controls a model solve by varying parameters dynamically. First, we validate our RL framework for the problem of controlling chaos in chaotic systems by dynamically changing the parameters of the system. Subsequently, we illustrate the capabilities of our framework for accelerating the convergence of a steady-state computational fluid dynamics solver by automatically adjusting the relaxation factors of the discretized Navier–Stokes equations during run-time. The results indicate that the run-time control of the relaxation factors by the learned policy leads to a significant reduction in the number of iterations for convergence compared to the random selection of the relaxation factors. Our results point to potential benefits for learning adaptive hyperparameter learning strategies across different geometries and boundary conditions with implications for reduced computational campaign expenses 4 4 Data and codes available at https://github.com/Romit-Maulik/PAR-RL . .

42 ENGINEERING↗
DOE OSTI2024

Why is My Zero Energy Home Not a Zero Carbon Home?

For years, carbon calculations were done very simply. The method of calculation was to take annual totals of energy consumption and multiply by an average emission factor, either for the grid serving a project or for a larger region (e.g. an EPA eGRID sub region). The level of accuracy of this approximation was reasonably good, although the issue of accuracy was not, to our knowledge, tested. And the data required were minimal – just a year’s worth of bills for each fuel and one lookup factor. But this method assures that a net zero energy home is automatically a net zero carbon home because zero times any possible emission factor is still zero. Starting in the early 2010s, things changed – grids were starting to rely more and more heavily on renewables, and the difference was showing up on aggregate load curves. This was perhaps noticed first in California, where aggressive renewable policies led to significant renewable power generation large enough to affect the overall shape of the diurnal load curve for the Independent Systems Operator.

14 SOLAR ENERGY↗
NASA NTRS1980

Axisymmetric elastodynamic response from normal and radial impact of layered composites with embedded penny-shaped cracks

A method is developed for the dynamic stress analysis of a layered composite containing an embedded penny-shaped crack and subjected to normal and radial impact. Quantitatively, the time-dependent stresses near the crack border can be described by the dynamic stress intensity factors. Their magnitude depends on time, on the material properties of the composite and on the relative size of the crack compared to the composite local geometry. Results obtained show that, for the same material properties and geometry of the composite, the dynamic stress intensity factors for an embedded (penny-shaped) crack reach their peak values within a shorter period of time and with a lower magnitude than the corresponding dynamic stress factors for a through-crack.

Sih, G. C.↗
NASA NTRS1986

Run-time implementation issues for real-time embedded Ada

A motivating factor in the development of Ada as the department of defense standard language was the high cost of embedded system software development. It was with embedded system requirements in mind that many of the features of the language were incorporated. Yet it is the designers of embedded systems that seem to comprise the majority of the Ada community dissatisfied with the language. There are a variety of reasons for this dissatisfaction, but many seem to be related in some way to the Ada run-time support system. Some of the areas in which the inconsistencies were found to have the greatest impact on performance from the standpoint of real-time systems are presented. In particular, a large part of the duties of the tasking supervisor are subject to the design decisions of the implementer. These include scheduling, rendezvous, delay processing, and task activation and termination. Some of the more general issues presented include time and space efficiencies, generic expansions, memory management, pragmas, and tracing features. As validated compilers become available for bare computer targets, it is important for a designer to be aware that, at least for many real-time issues, all validated Ada compilers are not created equal.

Maule, Ruth A.↗
NASA NTRS1979

Normal and radial impact of composites with embedded penny-shaped cracks

A method is developed for the dynamic stress analysis of a layered composite containing an embedded penny-shaped crack and subjected to normal and radial impact. The material properties of the layers are chosen such that the crack lies in a layer of matrix material while the surrounding material possesses the average elastic properties of a two-phase medium consisting of a large number of fibers embedded in the matrix. Quantitatively, the time-dependent stresses near the crack border can be described by the dynamic stress intensity factors. Their magnitude depends on time, on the material properties of the composite and on the relative size of the crack compared to the composite local geometry. Results obtained show that, for the same material properties and geometry of the composite, the dynamic stress intensity factors for an embedded (penny-shaped) crack reach their peak values within a shorter period of time and with a lower magnitude than the corresponding dynamic stress intensity factors for a through-crack.

Sih, G. C.↗
NASA NTRS1999

Dynamics of the Molten Contact Line

The purpose of this program is to develop a basic understanding of how a molten material front spreads over a solid that is below its melting point, arrests, and freezes. Our hope is that the work will contribute toward a scientific knowledge base for certain new applications involving molten droplet deposition, including the "printing" of arbitrary three-dimensional objects by precise deposition of individual molten microdrops that solidify after impact. Little information is available at this time on the capillarity-driven motion and arrest of molten contact line regions. Schiaffino and Sonin investigated the arrest of the contact line of a molten microcrystalline wax spreading over a subcooled solid "target" of the same material. They found that contact line arrest takes place at an apparent liquid contact angle that depends primarily on the Stefan number S=c(T(sub f) -T(sub t)/L based on the temperature difference between the fusion point and the target temperature, and proposed that contact line arrest occurs when the liquid's dynamic contact angle approaches the angle of attack of the solidification front just behind the contact line. They also showed, however, that the conventional continuum equations and boundary conditions have no meaningful solution for this angle. The solidification front angle is determined by the heat flux just behind the contact line, and the heat flux is singular at that point. By comparing experiments with numerical computations, Schiaffino and Sonin estimated that the conventional solidification model must break down within a distance of order 0.1 - 1 microns of the contact line. The physical mechanism for this breakdown is as yet undetermined, and no first-principles theory exists for the contact angle at arrest. Schiaffino and Sonin also presented a framework for understanding how to moderate Weber number molten droplet deposition in terms of similarity laws and experimentation. The study is based on experiments with three molten materials- molten wax on solid wax, water on ice, and mercury on frozen mercury- which between them span a considerable range of the deposition/solidification similarity parameters. Correlations are obtained for the spreading velocity, spreading time scales, the spreading factor (i.e. ratio of deposited drop's final footprint radius and the drop's initial radius), post-spreading liquid oscillation amplitudes and time scales, and bulk solidification time scales. Duthaler carried out an experimental and theoretical investigation of the relationship between the liquid's apparent contact angle and the Capillary number Ca=mu U/sigma based on contact line speed, for molten materials spreading over subcooled solids. This relationship is required for modeling of melt spreading. We have adapted Voinov's methodology to the molten contact line and formulated a theoretical model for the Ca vs. contact angle relationship, based Schiaffino and Sonin#s (1997a,b) wedge-like solidification front model. With the solidification front angle taken from Schiaffino and Sonin, the model is in good agreement with the experimental results for Ca vs. contact angle. Duthaler also extended the experimental investigation of droplet deposition and contact line freezing to more materials, including solder on glass, solder on solder, water on ice, and molten microcrystalline wax on wax. The latter also included tests on inclined targets. Deposition tests have also been done with molten octacosane (C28H58) on various targets. An important objective of our program has been the development of micron-scale sensors for measuring the transient temperature at a point on the substrate surface as a molten contact line moves over it. The expectation is that this temperature history will yield a better understanding of the thermal process in the contact line region. The sensors are of the thermistor type, either 2.5 microns or 1.5 microns square, microfabricated with silicon-based technology on either pure silicon or amorphous silicon dioxide chips. Each chip has 32 sensors distributed on its surface in arrays. The time response is better than 10 ms. At the time of writing, sensor calibration is in progress. Results on thermal transients during contact line passage will be discussed at the conference. While we expect that the data will provide information on the near-contact-line heat transfer process, we also foresee possible problems. First, the spatial resolution of the sensors may be insufficient to resolve the near-contact-line region. Second, the sensors protrude about 0.5 microns above the substrate surface, and may affect the contact line motion. Third, a sensor's temperature history depends on both the heat flux distribution into it from the fusion front and the thermal properties of the substrate below it and the solidified melt between it and the fusion front. The heat flux distribution in the contact line region must therefore be unfolded from computations of the overall system's transient thermal response.

Sonin, Ain A.↗
NASA NTRS2022

SAGE III/ISS Ozone and NO2 Validation Using Diurnal Scaling Factors

We developed a set of solar zenith angle, latitude- and altitude-dependent scaling factors to account for the diurnal variability in ozone (O3) and nitrogen dioxide (NO2) when comparing Stratospheric Aerosol and Gas Experiment (SAGE) III/ISS observations to observations from other times of day. The scaling factors are calculated as a function of solar zenith angle from the four-dimensional output of a global atmospheric chemistry model simulation of 2017–2020 that shows good agreement with observed vertical profiles. Using a global atmospheric chemistry model allows us to account for both chemically and dynamically driven variability. Both year-specific scale factors and a multi-year monthly climatology are available to decrease the uncertainty in inter-instrument comparisons and allow consistent comparisons between observations from different times of day. We describe the variability in the diurnal scale factors as a function of space and time. The quasi-biennial oscillation (QBO) appears to be a contributing factor to interannual variability in the NO2 scaling factors, leading to differences between years that switch sign with altitude. We show that application of these scaling factors improves the comparison between SAGE III/ISS and OSIRIS NO2 and between SAGE III/ISS and OMPS LP, OSIRIS, and ACE-FTS O3 observations. The comparisons between SAGE III/ISS O3 for sunrise or sunset vs. Microwave Limb Sounder (MLS) daytime or nighttime observations are also more consistent when we apply the diurnal scaling factors. There is good agreement between SAGE III/ISS V5.2 ozone and correlative measurements, with differences within 5 % between 20 and 50 km when corrected for diurnal variability. Similarly, the SAGE III/ISS V5.2 NO2 agreement with correlative measurement is mostly within 10 %. While the scale factors were designed for use with SAGE III/ISS observations, they can easily be applied to other observation intercomparisons as well.

Sarah A. Strode↗
DOE OSTI2025

Operator dynamics in Floquet many-body systems

We study operator dynamics in many-body quantum systems, focusing on generic features of systems that are ergodic, spatially extended, and lack conserved densities. Quantum circuits of various types provide simple models for such systems. We focus on Floquet quantum circuits, comparing their behavior with what has been found previously for circuits that are random in time. Floquet circuits, which have discrete time-translation symmetry, represent an intermediate case between circuits that are random in time and lack any symmetry, and systems with a time-independent Hamiltonian and continuous time-translation invariance. By making this comparison, one of our aims is to identify signatures of time-translation symmetry in Floquet operator dynamics. To characterize behavior we examine a variety of quantities in solvable models and numerically: operator autocorrelation functions; the partial spectral form factor; the out-of-time-order correlator (OTOC); and the paths in operator space that make the dominant contributions to the ensemble-averaged autocorrelation functions. Our most striking result is that ensemble-averaged autocorrelation functions show behavior that is distinctively different in Floquet systems compared to systems in which successive time-steps are independent. Specifically, while average autocorrelation functions decay on a microscopic timescale for circuits that are random in time, in Floquet systems they have a late-time tail with a duration that grows parametrically with the size of the operator support. In the simplest models this tail is separated from the initial decay by a minimum, so that the average autocorrelation function has an intermediate-time peak. The existence of these tails provides a way to understand deviations of the spectral form factor from random matrix behavior at times shorter than the Thouless time. In contrast to this feature in autocorrelation functions, we find no new aspects to the behavior of OTOCs for Floquet models compared to random-in-time circuits. We show that this difference between averaged autocorrelation functions and OTOCs can be understood in terms of the paths in operator space that contribute to the two quantities: paths for the former retain a limited support at late times, while paths for the latter are dominated by operator spreading. Published by the American Physical Society 2025

Yoshimura, Takato (ORCID:0000000309159846)↗
DOE OSTI2022

Experimental demonstration of continuous quantum error correction

The storage and processing of quantum information are susceptible to external noise, resulting in computational errors. A powerful method to suppress these effects is quantum error correction. Typically, quantum error correction is executed in discrete rounds, using entangling gates and projective measurement on ancillary qubits to complete each round of error correction. Here we use direct parity measurements to implement a continuous quantum bit-flip correction code in a resource-efficient manner, eliminating entangling gates, ancillary qubits, and their associated errors. An FPGA controller actively corrects errors as they are detected, achieving an average bit-flip detection efficiency of up to 91%. Furthermore, the protocol increases the relaxation time of the protected logical qubit by a factor of 2.7 over the relaxation times of the bare comprising qubits. Our results showcase resource-efficient stabilizer measurements in a multi-qubit architecture and demonstrate how continuous error correction codes can address challenges in realizing a fault-tolerant system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Constraining Type Ia supernova explosions and early flux excesses with the Zwicky Transient Factory

In the new era of time-domain surveys, Type Ia supernovae are being caught sooner after explosion, which has exposed significant variation in their early light curves. Two driving factors for early-time evolution are the distribution of 56 Ni in the ejecta and the presence of flux excesses of various causes. We perform an analysis of the largest young SN Ia sample to date. We compare 115 SN Ia light curves from the Zwicky Transient Facility to the turtls model grid containing light curves of Chandrasekhar mass explosions with a range of 56 Ni masses, 56 Ni distributions, and explosion energies. We find that the majority of our observed light curves are well reproduced by Chandrasekhar mass explosion models with a preference for highly extended 56 Ni distributions. We identify six SNe Ia with an early-time flux excess in our gr-band data (four ‘blue’ and two ‘red’ flux excesses). We find an intrinsic rate of 18 ± 11 percent of early flux excesses in SNe Ia at z < 0.07, based on three detected flux excesses out of 30 (10 percent) observed SNe Ia with a simulated efficiency of 57 percent. This is comparable to rates of flux excesses in the literature but also accounts for detection efficiencies. Two of these events are mostly consistent with circumstellar material interaction, while the other four have longer lifetimes in agreement with companion interaction and 56 Ni-clump models. We find a higher frequency of flux excesses in 91T/99aa-like events (44 ± 13 percent).

79 ASTRONOMY AND ASTROPHYSICS↗