Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TIME FUNCTION”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Theoretical approach for optical response in electrochemical systems: Application to electrode potential dependence of surface-enhanced Raman scattering

We propose a theoretical approach for optical response in electrochemical systems. The fundamental equation to be solved is based on a time-dependent density functional theory in real-time and real-space in combination with its finite temperature formula treating an electrode potential. Solvation effects are evaluated by a dielectric continuum theory. The approach allows us to treat optical response in electrochemical systems at the atomistic level of theory. We have applied the method to surface-enhanced Raman scattering (SERS) of 4-mercaptopyridine on an Ag electrode surface. It is shown that the SERS intensity has a peak as a function of the electrode potential. Furthermore, the real-space computational approach facilitates visualization of variation of the SERS intensity depending on an electrode potential.

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

Highly efficient and selective pressure-assisted photon-induced polymerization of styrene

The polymerization process of condensed styrene to produce polystyrene as an industrially important polymeric material was investigated using a novel approach by combining external compression with ultraviolet radiation. The reaction evolution was monitored as a function of time and the reaction products were characterized by in situ Fourier transform infrared spectroscopy. By optimizing the loading pressures, we observed highly efficient and selective production of polystyrene of different tacticities. Specifically, at relatively low loading pressures, infrared spectra suggest that styrene monomers transform to amorphous atactic polystyrene (APS) with minor crystalline isotactic polystyrene. In contrast, APS was found to be the sole product when polymerization occurs at relatively higher loading pressures. The time-dependent reaction profiles allow the examination of the polymerization kinetics by analyzing the rate constant and activation volume as a function of pressure. As a result, an optimized pressure condition, which allows a barrierless reaction to proceed, was identified and attributed to the very desirable reaction yield and kinetics. Finally, the photoinitiated reaction mechanism and the growth geometry of the polymer chains were investigated from the energy diagram of styrene and by the topology analysis of the crystal styrene. This study shows strong promise to produce functional polymeric materials in a highly efficient and controlled manner.

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

Time evolution of the electron energy distribution function in pulsed microwave magnetoplasma in H{sub 2}

Time evolution of the Electron Energy Distribution Function (EEDF) is measured in pulsed hydrogen microwave magnetoplasma working at 2.45 GHz. Analysis is performed both in resonance (B = 0.087 T) and off-resonance conditions (B = 0.120 T), at two pressures (0.38 Pa and 0.62 Pa), respectively, and for different incident microwave powers. The important effect of the magnetic field on the electron kinetic is discussed, and a critical analysis of Langmuir probe measurements is given. The Electron Energy Distribution Function is calculated using the Druyvesteyn theory (EEDF) and is corrected using the theory developed by Arslanbekov in the case of magnetized plasma. Three different components are observed in the EEDF, whatever the theory used. They are: (a) a low electron energy component at energy lower than 10 eV, which is ascribed to the electron having inelastic collisions with heavy species (H{sub 2}, H, ions), (b) a high energy component with a mean energy ranging from 10 to 20 eV, which is generally ascribed to the heating of the plasma by the incident microwave power, and (c) a third component observed between the two other ones, mainly at low pressure and in resonance conditions, has been correlated to the electron rotation in the magnetic field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

A Compendium of Wind Statistics and Models for the NASA Space Shuttle and Other Aerospace Vehicle Programs

The wind profile with all of its variations with respect to altitude has been, is now, and will continue to be important for aerospace vehicle design and operations. Wind profile databases and models are used for the vehicle ascent flight design for structural wind loading, flight control systems, performance analysis, and launch operations. This report presents the evolution of wind statistics and wind models from the empirical scalar wind profile model established for the Saturn Program through the development of the vector wind profile model used for the Space Shuttle design to the variations of this wind modeling concept for the X-33 program. Because wind is a vector quantity, the vector wind models use the rigorous mathematical probability properties of the multivariate normal probability distribution. When the vehicle ascent steering commands (ascent guidance) are wind biased to the wind profile measured on the day-of-launch, ascent structural wind loads are reduced and launch probability is increased. This wind load alleviation technique is recommended in the initial phase of vehicle development. The vehicle must fly through the largest load allowable versus altitude to achieve its mission. The Gumbel extreme value probability distribution is used to obtain the probability of exceeding (or not exceeding) the load allowable. The time conditional probability function is derived from the Gumbel bivariate extreme value distribution. This time conditional function is used for calculation of wind loads persistence increments using 3.5-hour Jimsphere wind pairs. These increments are used to protect the commit-to-launch decision. Other topics presented include the Shuttle Shuttle load-response to smoothed wind profiles, a new gust model, and advancements in wind profile measuring systems. From the lessons learned and knowledge gained from past vehicle programs, the development of future launch vehicles can be accelerated. However, new vehicle programs by their very nature will require specialized support for new databases and analyses for wind, atmospheric parameters (pressure, temperature, and density versus altitude), and weather. It is for this reason that project managers are encouraged to collaborate with natural environment specialists early in the conceptual design phase. Such action will give the lead time necessary to meet the natural environment design and operational requirements, and thus, reduce development costs.

Smith, O. E.↗
DOE OSTI2014

Can the ring polymer molecular dynamics method be interpreted as real time quantum dynamics?

The ring polymer molecular dynamics (RPMD) method has gained popularity in recent years as a simple approximation for calculating real time quantum correlation functions in condensed media. However, the extent to which RPMD captures real dynamical quantum effects and why it fails under certain situations have not been clearly understood. Addressing this issue has been difficult in the absence of a genuine justification for the RPMD algorithm starting from the quantum Liouville equation. To this end, a new and exact path integral formalism for the calculation of real time quantum correlation functions is presented in this work, which can serve as a rigorous foundation for the analysis of the RPMD method as well as providing an alternative derivation of the well established centroid molecular dynamics method. The new formalism utilizes the cyclic symmetry of the imaginary time path integral in the most general sense and enables the expression of Kubo-transformed quantum time correlation functions as that of physical observables pre-averaged over the imaginary time path. Upon filtering with a centroid constraint function, the formulation results in the centroid dynamics formalism. Upon filtering with the position representation of the imaginary time path integral, we obtain an exact quantum dynamics formalism involving the same variables as the RPMD method. The analysis of the RPMD approximation based on this approach clarifies that an explicit quantum dynamical justification does not exist for the use of the ring polymer harmonic potential term (imaginary time kinetic energy) as implemented in the RPMD method. It is analyzed why this can cause substantial errors in nonlinear correlation functions of harmonic oscillators. Such errors can be significant for general correlation functions of anharmonic systems. We also demonstrate that the short time accuracy of the exact path integral limit of RPMD is of lower order than those for finite discretization of path. The present quantum dynamics formulation also serves as the basis for developing new quantum dynamical methods that utilize the cyclic nature of the imaginary time path integral.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

THE SPECTRAL SN-GRB CONNECTION: SYSTEMATIC SPECTRAL COMPARISONS BETWEEN TYPE Ic SUPERNOVAE AND BROAD-LINED TYPE Ic SUPERNOVAE WITH AND WITHOUT GAMMA-RAY BURSTS

We present the first systematic investigation of spectral properties of 17 Type Ic Supernovae (SNe Ic), 10 broad-lined SNe Ic (SNe Ic-bl) without observed gamma-ray bursts (GRBs), and 11 SNe Ic-bl with GRBs (SN-GRBs) as a function of time in order to probe their explosion conditions and progenitors. Using a number of novel methods, we analyze a total of 407 spectra, which were drawn from published spectra of individual SNe as well as from the densely time-sampled spectra of Modjaz et al (2014). In order to quantify the diversity of the SN spectra as a function of SN subtype, we construct average spectra of SNe Ic, SNe Ic-bl without GRBs, and SNe Ic-bl with GRBs. We find that SN 1994I is not a typical SN Ic, contrasting the general view, while the spectra of SN 1998bw/GRB 980425 are representative of mean spectra of SNe Ic-bl. We measure the ejecta absorption and width velocities using a new method described here and find that SNe Ic-bl with GRBs, on average, have quantifiably higher absorption velocities, as well as broader line widths than SNe without observed GRBs. In addition, we search for correlations between SN-GRB spectral properties and the energies of their accompanying GRBs. Finally, we show that the absence of clear He lines in optical spectra of SNe Ic-bl, and in particular of SN-GRBs, is not due to them being too smeared-out due to the high velocities present in the ejecta. This implies that the progenitor stars of SN-GRBs are probably free of the He-layer, in addition to being H-free, which puts strong constraints on the stellar evolutionary paths needed to produce such SN-GRB progenitors at the observed low metallicities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Time-dependent orbital-free density functional theory for electronic stopping power: Comparison to the Mermin-Kohn-Sham theory at high temperatures

Electronic stopping power in warm dense matter can affect energy transport and heating in astrophysical processes and internal confinement fusion. For cold condensed matter systems, stopping power can be modeled from first-principles using real-time time-dependent density functional theory (DFT). However, high temperatures (10's to 100's of eV) may be computationally prohibitive for traditional Mermin-Kohn-Sham DFT. New experimental measurements in the warm dense regime motivates the development of first-principles approaches, which can reach these temperatures. Here, we have developed a time-dependent orbital-free density functional theory, which includes a novel nonadiabatic and temperature-dependent kinetic energy density functional, for the simulation of stopping power at any temperature. The approach is nonlinear with respect to the projectile perturbation, includes all ions and electrons, and does not require a priori determination of screened interaction potentials. Finally, our results compare favorably with Kohn-Sham for temperatures in the WDM regime, especially nearing 100 eV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

Ares I-X In-Flight Modal Identification

Operational modal analysis is a procedure that allows the extraction of modal parameters of a structure in its operating environment. It is based on the idealized premise that input to the structure is white noise. In some cases, when free decay responses are corrupted by unmeasured random disturbances, the response data can be processed into cross-correlation functions that approximate free decay responses. Modal parameters can be computed from these functions by time domain identification methods such as the Eigenvalue Realization Algorithm (ERA). The extracted modal parameters have the same characteristics as impulse response functions of the original system. Operational modal analysis is performed on Ares I-X in-flight data. Since the dynamic system is not stationary due to propellant mass loss, modal identification is only possible by analyzing the system as a series of linearized models over short periods of time via a sliding time-window of short time intervals. A time-domain zooming technique was also employed to enhance the modal parameter extraction. Results of this study demonstrate that free-decay time domain modal identification methods can be successfully employed for in-flight launch vehicle modal extraction.

Bartkowicz, Theodore J.↗
DOE OSTI2024

Evaluation of RANS vs. LES simulation of fluid flow through 3 × 3 rod bundle with a simple spacer grid as a precursor to coupled fluid–structure interaction simulations

The research literature on Computational Fluid Dynamics (CFD) of coolant flow through rod bundles with spacer-grids and mixing vanes is replete, ranging from high fidelity Large Eddy Simulation (LES)/Direct Numerical Simulation (DNS) simulations to Reynolds-Averaged Navier–Stokes (RANS) modeled studies. The mixing of flow between subchannels and the pressure drop through the bundle are fundamental quantities useful for comparing and evaluating CFD methods. Less commonly observed and compared are the forces exerted onto the structure by the fluid. The present study seeks to evaluate the use of RANS simulations for predicting the structural response to fluid flow. Wall resolved RANS simulations are benchmarked against LES simulations of fluid flow at a Reynolds number of 15,000 through a 3 × 3 fuel rod bundle with a simple spacer grid. Velocity line-plots are compared showing good agreement between RANS and LES results, ascertaining that the former is capable of capturing the essential time-averaged velocity profile. Additionally, the distribution of forces on the spacer grid and fuel rods are collected as a function of time and space. The RANS methods are evaluated using the frequency and magnitude of the fluctuating forces on various portions of the structure as compared to LES. In conclusion, the power spectral density evaluation of the models reveal underprediction of force amplitude on the rod walls by RANS and also discrepancy in the prediction of high frequency spectra, especially in the immediate vicinity of spacer-grid structure, which may be attributed to the lack of random turbulence fluctuation or insufficient modeling of small-scale eddies in RANS simulation.

FIV↗
DOE OSTI2022

Open-circuit submodule fault diagnosis in MMCs using support vector machines

Series connection of semiconductor submodules (SM) in a modular multilevel converter (MMC) makes the MMC prone to open-circuit (OC) IGBT failures inside SMs. If left undetected, these faults degrade the operation of the MMC and lead to its instability. This article proposes a method to detect, localise, and classify single OC SM faults in an MMC using support vector machines (SVM) trained with data obtained from the capacitor voltage balancing block of the MMC control system. The proposed method relies on data extracted from the sorted capacitor voltage arrays of the upper and lower phase arms. Therefore, it does not require extra measurements and hardware. Additionally, it offers a fixed time for detecting and localising OC SM faults. This method is easy to implement as SVM has a simple decision function. Time-domain simulation case studies are performed on a three-phase nine-level MMC to evaluate the performance of the proposed method.

42 ENGINEERING↗
DOE OSTI2019

Real-time density-matrix coupled-cluster approach for closed and open systems at finite temperature

We extend the coupled-cluster method to correlated quantum dynamics of both closed and open systems at finite temperatures using the thermofield formalism. The approach expresses the time-dependent density matrix in an exponential ansatz and describes time-evolution along the Keldysh path contour. A distinct advantage of the approach is exact trace-preservation as a function of time, ensuring conservation of probability and particle number. Furthermore, the method avoids the computation of correlated bra-states, simplifying the computational implementation. We develop the method in a thermal quasiparticle representation, which allows seamless connection to the projection method and diagrammatic techniques of the traditional coupled-cluster formalism. For comparison, we also apply the thermofield framework to the density-matrix renormalization-group method to obtain reference results for closed and open systems at finite temperature. We test the singles and doubles approximation to the density-matrix coupled-cluster method on the correlated electronic dynamics of the single-impurity Anderson model, demonstrating that the new method successfully captures the correlated dynamics of both closed systems at finite temperature and driven-dissipative open systems. This encouraging performance motivates future applications to nonequilibrium quantum many-body dynamics in realistic systems.

Shushkov, Philip (ORCID:0000000266136117)↗
DOE OSTI2017

Solids Erosion Patterns Developed by Pulse Jet Mixers

Millions of gallons of radioactive waste are stored in underground storage tanks at the Hanford Site in Washington State. This waste will be vitrified at the Waste Treatment and Immobilization Plant that is under construction. Vessels in the pretreatment portion of the plant are being configured for processing waste slurries with challenging physical and rheological properties that range from Newtonian slurries to non-Newtonian sludge. Pulse jet mixing technology has been selected for mobilizing and mixing this waste. In the pulse jet mixing process, slurry is expelled from pulse tube nozzles directed towards the vessel floor. The expelled fluid forms a radial jet that erodes the settled layer of solids. The pulse tubes are configured in a ring or multiple rings and operate concurrently. The expelled fluid and mobilized solids traverse toward the center of the tank. At the tank center the jets from pulse tubes in the ring collide and lift solids upward in a central plume. At the end of the pulse, when the desired fluid volume is expelled from the pulse tube, the applied pressure switches to suction and the pulse tube is refilled. This cycle is used to mobilize and mix the tank contents. An initial step of the process is the erosion of solids from the vessel floor by the radial jets that form on the vessel flow beneath each pulse tube. Experiments have been conducted using simulants to evaluate the ability of the pulse jet mixing system radial jets to combine to develop the central upwell and lift solids into the vessel. These experiments have been conducted at three scales using a range of granular simulants over a range of concentrations. The vessels have elliptical, spherical, or flanged and dished bottoms. Process parameters evaluated include the velocity of fluid expelled from the pulse tube, the duration of the pulse and the duty cycle, the ratio of pulse duration to cycle time. Videos taken from beneath the vessel show the growth of the cleared area from each pulse tube as a function of time. All solids are lifted from the vessel bottom when the system is operating at the critical suspension velocity. The focus of this paper is to compare and contrast erosion patterns developed from different simulants and pulse tube configurations. The cases are evaluated to determine how changes in process parameters affects the PJM ability to mobilize solids from the vessel floor.

erosion↗
DOE OSTI2026

Mathematical Modeling of the Potential and Time Dependence of Ir Dissolution from Hydrous Ir Oxide Oxygen Evolution Catalysts

One of the main degradation mechanisms of hydrous iridium oxide acidic oxygen evolution reaction (OER) catalysts is dissolution and loss into the acidic membrane. While degradation models have been proposed, there is a gap in understanding the potential and time dependence of the iridium dissolution reaction and its mechanistic underpinnings. In this work, Ir dissolution rates measured as a function of time and potential via time-resolved inductively-coupled plasma mass spectrometry (ICP-MS) in aqueous acidic electrolyte are used to establish a mathematical model for Ir dissolution. The mathematical model is generated using proposed formation and dissolution reactions for Ir species. Through comparison with the ICP-MS data and existing information on the potential-dependent Ir phase, we find that the potential and time-dependence of dissolution can be modeled as dissolution of an oxide phase, here represented as IrO 2 , with potential dependent kinetics and formation of a passivating species, a process with a rate-limiting step that is not potential dependent. This understanding of the potential dependence of dissolution and passivation kinetics using aqueous electrolyte half-cell measurements can be used to predict the degradation of Ir oxide in operating energy conversion devices relying on the OER, such as proton-exchange membrane water electrolyzers.

Kariuki, Nancy N. [Argonne National Laboratory (AN↗
DOE OSTI2016

DERIVING THE PROPERTIES OF CORONAL PRESSURE FRONTS IN 3D: APPLICATION TO THE 2012 MAY 17 GROUND LEVEL ENHANCEMENT

We study the link between an expanding coronal shock and the energetic particles measured near Earth during the ground level enhancement of 2012 May 17. We developed a new technique based on multipoint imaging to triangulate the three-dimensional (3D) expansion of the shock forming in the corona. It uses images from three vantage points by mapping the outermost extent of the coronal region perturbed by the pressure front. We derive for the first time the 3D velocity vector and the distribution of Mach numbers, M {sub FM}, of the entire front as a function of time. Our approach uses magnetic field reconstructions of the coronal field, full magnetohydrodynamic simulations and imaging inversion techniques. We find that the highest M {sub FM} values appear near the coronal neutral line within a few minutes of the coronal mass ejection onset; this neutral line is usually associated with the source of the heliospheric current and plasma sheet. We illustrate the variability of the shock speed, shock geometry, and Mach number along different modeled magnetic field lines. Despite the level of uncertainty in deriving the shock Mach numbers, all employed reconstruction techniques show that the release time of GeV particles occurs when the coronal shock becomes super-critical ( M {sub FM} > 3). Combining in situ measurements with heliospheric imagery, we also demonstrate that magnetic connectivity between the accelerator (the coronal shock of 2012 May 17) and the near-Earth environment is established via a magnetic cloud that erupted from the same active region roughly five days earlier.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

A model for the kinetics of a solar-pumped long path laser experiment

A kinetic model for a solar-simulator pumped iodine laser system is developed and compared to an experiment in which the solar simulator output is dispersed over a large active volume (150 cu cm) with low simulator light intensity (approx. 200 solar constants). A trace foreign gas which quenches the upper level is introduced into the model. Furthermore, a constant representing optical absorption of the stimulated emission is introduced, in addition to a constant representing the scattering at each of the mirrors, via the optical cavity time constant. The non-uniform heating of the gas is treated as well as the pressure change as a function of time within the cavity. With these new phenomena introduced into the kinetic model, a best reasonable fit to the experimental data is found by adjusting the reaction rate coefficients within the range of known uncertainty by numerical methods giving a new bound within this range of uncertainty. The experimental parameters modeled are the lasing time, laser pulse energy, and time to laser threshold.

Stock, L. V.↗
NASA NTRS1994

An automated digital data collection and analysis system for the Charpy Impact Tester

The standard Charpy Impact Tester has been modified by the addition of a system of hardware and software to improve the accuracy and consistency of measurements made during specimen fracturing experiments. An optical disc, light source, and detector generate signals that indicate the pendulum position as a function of time. These signals are used by a computer to calculate the velocity and kinetic energy of the pendulum as a function of its position.

Kohne, Glenn S.↗
NASA NTRS2001

Analysis of Ag/FEP Control Blanket Performance from Multiple Satellites

Published data reporting the change in solar absorptance as a function of time on orbit for certain satellites has been re-examined. The optical property change of metallized FEP as a function of total particulate radiation dose received by specific satellites has been estimated. Results tend to follow a common trend and a design curve has been constructed to estimate solar absorptance change as a function of radiation dose for metallized FEP at any altitude. Effects due to apparent contamination on several of the satellites are addressed and taken into account in determining property changes. Data from LDEF, ML-101, NavStar 1-5, NTS-2, the SCATHA experiment, and SOLRAD 11 are used as the basis of the design curve. A similar comparison is made between published atomic oxygen recession data from numerous Space Shuttle flights, a LEO experiment flown by Lockheed, and LDEF results. Variation of atomic oxygen recession rate of FEP as a function of exposure conditions is determined.

Pippin, H. G.↗