Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TIME DEPENDENCE”

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

Prediction of photodynamics of 200 nm excited cyclobutanone with linear response electronic structure and ab initio multiple spawning

Simulations of photochemical reaction dynamics have been a challenge to the theoretical chemistry community for some time. In an effort to determine the predictive character of current approaches, we predict the results of an upcoming ultrafast diffraction experiment on the photodynamics of cyclobutanone after excitation to the lowest lying Rydberg state (S 2 ). A picosecond of nonadiabatic dynamics is described with ab initio multiple spawning. Herein we use both time dependent density functional theory (TDDFT) and equation-of-motion coupled cluster singles and doubles (EOM-CCSD) theory for the underlying electronic structure theory. We find that the lifetime of the S 2 state is more than a picosecond (with both TDDFT and EOM-CCSD). The predicted ultrafast electron diffraction spectrum exhibits numerous structural features, but weak time dependence over the course of the simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamic pair-breaking current, critical superfluid velocity, and nonlinear electromagnetic response of nonequilibrium superconductors

Here, we report numerical calculations of a dynamic pair-breaking current density $J_d$ and a critical superfluid velocity $v_d$ in a nonequilibrium superconductor carrying a uniform, large-amplitude AC current density $J(t) = J_a$sinΩ$\textit{t}$ with Ω well below the gap frequency Ω $\ll Δ_0/ℏ$. The dependencies $J_d$(Ω, $\textit{T}$) and $v_d$(Ω, $\textit{T}$) near the critical temperature $T_c$ were calculated from either the full time-dependent nonequilibrium equations for a dirty $\textit{s}$-wave superconductor or the time-dependent Ginzburg-Landau (TDGL) equations for a gapped superconductor, taking into account the GL relaxation time of the order parameter $τ_{\text{GL}}$ and the inelastic electron-phonon relaxation time of quasiparticles $τ_E$. We show that both approaches give similar frequency dependencies of $J_d$(Ω) and $v_d$(Ω) which gradually increase from their static pair-breaking GL values $J_c$ and $v_c$ at Ω$τ_E \ll$1 to $\sqrt{2}J_c$ and $\sqrt{2}v_c$ at Ω$τ_E \gg$ 1. Here $J_d, v_d$ and a dynamic superheating field at which the Meissner state becomes unstable were calculated in two different regimes of a fixed AC current and a fixed AC superfluid velocity induced by the applied AC magnetic field $H = H_a$ sinΩ$\textit{t}$ in a thin superconducting filament or a type-II superconductor with a large GL parameter. We also calculated a nonlinear electromagnetic response of a nonequilibrium superconducting state, particularly a dynamic kinetic inductance and a dissipative quasiparticle conductivity, taking into account the oscillatory dynamics of superconducting condensate and the kinetics of quasiparticles driven by a strong AC current. It is shown that an AC current density produces multiple harmonics of the electric field, the amplitudes of the higher-order harmonics diminishing as $τ_E$ increases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Rheological transition in mantle convection with a composite temperature-dependent, non-Newtonian and Newtonian rheology

Numerical simulations of mantle convection with a composite temperature-dependent, Newtonian and non-Newtonian creep law have revealed a transition in the dominant creep mechanism with the increasing vigour of convection. Newtonian creep is found to dominate in the low Rayleigh number regime. With sufficiently high effective Rayleigh number, the overall creep mechanism in the convective flow becomes non-Newtonian. The transitional Rayleigh number increases strongly with the activation energy. These results would suggest a scenario that in the early epochs of Earth the flow in the mantle would have been governed by non-Newtonian rheology and would have exhibited both strong spatial and temporal fluctuations. With time the flow mechanism would behave like a Newtonian fluid and would have a different time-dependent character. In time-dependent Newtonian-dominated flows there are still localized features with distinctly non-Newtonian character. Our analysis of the relative contributions to the lateral viscosity field supports the idea that the inference of the nature of lateral viscosity heterogeneities by seismic tomography may be strongly contaminated by the dominant non-Newtonian contributions to the total lateral viscosity field.

Van Den Berg, Arie P.↗

Influence of geochemical reactions on the creep behavior of Mt. Simon sandstone

The capture and subsurface storage of carbon dioxide is a sustainable option that is currently pursued worldwide to mitigate greenhouse gas effect. However, predicting the long-term mechanical integrity of CO 2 underground storage systems remains a challenge. To address that question, it is essential to understand the influence of fluid-rock chemo-mechanical interactions on the long-term and on the time-dependent mechanical properties. In turn, the long-term mechanical response and the time-dependent mechanical behavior can be represented by the creep response. We investigate the impact of CO 2 -induced geochemical reaction on the creep response of Mt. Simon sandstone with a 50–400 micron grain size. We perform static and dynamic flow experiments on Mt. Simon sandstone specimens under geological conditions, at a temperature of T = 50 - 53 °C and for CO 2 pressures of P = 8:62; 17:2 MPa under both static flow and dynamic flow-through conditions. After aging, we employ creep indentation testing, high-resolution SEM-EDS, computer vision, machine learning, and micromechanics modeling to probe changes on the microstructure and mechanical properties. Following both static and dynamic flow-through experiments, we observe a 10–22% decrease in quartz volume fraction and an increase in both the microporosity (7–28%) and nanoporosity (60–65%). Additional CO 2 -induced microstructural changes include an enlargement of pore throats and the formation of channels. These observations point to the presence of CO 2 -induced K-feldspar dissolution and clay dissolution reactions. Here, the macroscopic creep behavior is logarithmic and the macroscopic creep modulus varies depending on the microporosity and the relative quartz and feldspar content. As a result of these geochemical reactions and of the related microstructural changes, a 55-60% decrease in the macroscopic logarithmic creep modulus is predicted.

58 GEOSCIENCES↗

Design Methods, Tools, and Data for Ceramic Solar Receivers Year 1 Continuation Report

This report describes the first year of work on a project to develop the methods, tools, and data required to analyze high temperature ceramic Concentrating Solar Power (CSP) components. This first year focused on developing the methods and data required for a time-independent assessment of potential components, focusing in particular on ceramic solar receivers. The report describes both model development and testing work focused on accomplishing this goal. Our overall conclusion is that high temperature ceramic CSP components are viable and could provide a means to overcome the expected low reliability and short service life for equivalent components constructed from Ni-based superalloys. Based on the results reported here, we recommend the project continue to Phase II which will develop more sophisticated, realistic models for time-dependent failure of ceramics operating in expected CSP component conditions and develop the time-dependent ceramic test data needed to parameterize these models, using commercial SiC as a reference material.

14 SOLAR ENERGY↗

Perturbation theory for a resistivity shunted Josephson element

We present a systematic perturbation theory, extendable in principle to all orders of magnitude, for the solution of the equations of motion of an ideal Josephson element shunted by a resistance and driven by a dc current source and a small time-dependent source. We present second-order results for the case in which the time dependence is that of a single sinusoid, and these results are compared with other numerical and analytical calculations. Near, but not on, the first constant voltage step where the perturbation theory appears divergent, the current-voltage characteristic is calculated by means of a nonperturbative adiabatic procedure. The impedance and responsivity agree with earlier results.

Thompson, E. D.↗

Simulated traveling interplanetary disturbances initiated by various solar phenomena

Earlier studies have shown that interplanetary disturbances observed by spacecraft at large heliocentric distances from the sun can be directly associated with various solar activities such as surges, sprays, eruptive prominences, and flares. A digital computer code to simulate these phenomena is developed. This numerical code solves a set of conservation equations for an adiabatic time-dependent spherically-symmetric fluid flow. The solar events are simulated by perturbations in the dependent variables at the lower boundary of the quiet sun corona, and the time-dependent numerical solution simulates the coronal response. Three physical cases are simulated by using the proposed numerical code: (1) part of a surge, spray and prominence; (2) a single flare; and (3) two flares whose interplanetary consequences interact in space. Results for the first two cases are compared with both ground-based and space probe observations. Fairly good agreement with observations is found. A comparison of the third case with observations is presently under way.

Dryer, M.↗

Viscosity Relaxation in Molten HgZnTe

Rotating cup measurements of the viscosity of the pseudo-binary melt, HgZnTe have shown that the isothermal liquid with zinc mole fraction 0.16 requires tens of hours of equilibration time before a steady viscous state can be achieved. Over this relaxation period, the viscosity at 790 C increases by a factor of two, while the viscosity at 810 C increases by 40%. Noting that the Group VI elements tend to polymerize when molten, we suggest that the viscosity of the melt is enhanced by the slow formation of Te atom chains. To explain the build-up of linear Te n-mers, we propose a scheme, which contains formation reactions with second order kinetics that increase the molecular weight, and decomposition reactions with first order kinetics that inactivate the chains. The resulting rate equations can be solved for the time dependence of each molecular weight fraction. Using these molecular weight fractions, we calculate the time dependence of the average molecular weight. Using the standard semi-empirical relation between polymer average molecular weight and viscosity, we then calculate the viscosity relaxation curve. By curve fitting, we find that the data imply that the rate constant for n-mer formation is much smaller than the rate constant for n-mer deactivation, suggesting that Te atoms only weakly polymerize in molten HgZnTe. The steady state toward which the melt relaxes occurs as the rate of formation of an n-mer becomes exactly balanced by the sum of the rate for its deactivation and the rate for its polymerization to form an (n+1)-mer.

Su, Ching-Hua↗

Aspects of modelling the tectonics of large volcanoes on the terrestrial planets

Analytic solutions for the responses of planetary lithospheres to volcanic loads have been used to model faulting and infer elastic plate thicknesses. Predictions of the distribution of faulting around volcanic loads, based on the application of Anderson's criteria for faulting to the results of the models, do not agree well with observations. Such models do not give the stress state in the load itself, but only suggest a state of horizontal compressive stress there. Further, these models have considered only the effect of an instantaneously emplaced load. They do not address the time evolution of stresses, nor do they consider the effect of a load which grows. A finite element approach allows us to assign elements to the load itself, and thus permits calculation of the stress state and stress history within the edifice. The effects of episodic load growth can also be treated. When these effects are included, models give much better agreement with observations. We use the finite element code TECTON to construct axisymmetric models of volcanoes resting on an elastic lithospheric plate overlying a viscoelastic asthenosphere. We have implemented time-dependent material properties in order to simulate incremental volcano growth. The viscoelastic layer was taken to extend to a sufficient depth so that a rigid lower boundary has no significant influence on the results. The code first calculates elastic deformations and stresses and then determines the time-dependent viscous deformations and stresses. Time in the model scales as the Maxwell time tau(m) in the asthenosphere. We consider a volcano 25 km in height and 200 km in radius on an elastic lithosphere 40 km thick (parameters approximately appropriate to Ascraeus Mons). The volcano consists of three load increments applied at intervals of 1000 tau(m). Contours of maximum deviatoric stress in the fully-grown edifice at the conclusion of flexure (t = 3000 tau(m)) are shown.

Mcgovern, Patrick J.↗

Echo mapping the Balmer-emission region in NGC 3516

We use a maximum-entropy method to show that the transfer function (TF) of the broad-line region (BLR) in the Seyfert-1 galaxy NGC 3516 is time-varying. The TF decreased by nearly a factor of two over a time scale less than half a year during a monitoring campaign in 1990. We conclude that the observed time dependency is likely to be either due to variations of the spectral index of the ionizing continuum spectrum; due to a change in the number or covering fraction of broad-line clouds; due to nonlinear line response; or due to nonstationary anisotropy of the continuum source on time scales of several weeks to months. An extended continuum source could also explain the observed time-dependency. The symmetry of the time lag for variations in the red and blue wings of H-alpha indicates that the BLR kinematics is not dominated by organized radial inflow or outflow. If we assume circular orbits, the observed time lag 14 days at velocity v = 4000 km/s suggests a mass of the central object of approximately 2 x 10(exp 7) solar mass. The H-alpha TF peaks away from zero delay, indicating that the H-alpa BLR is either non-spherical or inwardly emitting.

Wanders, I.↗

Characterization of Fuel-to-Coolant Heat Transfer During Reactivity-Initiated Accidents Using Tightly Coupled Thermal Hydraulics and Fuel Thermomechanics

The reactivity-initiated accident (RIA) is a complex scenario with several tightly interacting physical phenomena. Accurately predicting fuel behavior during these transients is difficult due to limitations in the modeling of fuel-to-coolant heat transfer. Common approaches to simulate RIAs involve standalone calculations using either a fuel performance code or a thermal-hydraulic code. The complex interdependencies of thermal-hydraulic and fuel mechanical behavior suggest that a tight coupling between these codes may provide more accurate predictions of fuel-to-coolant heat transfer and cladding mechanical response. Here, RELAP5-3D and BISON are coupled in this paper to simulate RIAs, and a sensitivity analysis is performed to rank key thermal properties and two-phase heat transfer parameters relevant for fuel-to-coolant heat transfer and cladding failure mechanisms in UO 2 –Zircaloy-4 systems. Gas gap conductance, film boiling heat transfer uncertainty, pulse width, fuel-specific heat capacity, and cladding-specific heat capacity were identified as important parameters. Variations in figures of merit resulting from changes to pulse width and the material thermal properties indicate that time-dependent heat transfer rates are significant for safety-relevant mechanical parameters due to the time dependence of cladding ductility and pellet-cladding mechanical interaction loading. The results suggest that the thermal-hydraulic factors have a nonnegligible influence on the thermomechanical solution and vice versa. Tight coupling of both sets of physics is recommended to improve prediction of fuel behavior during RIAs. Highlights include the following: 1. The RELAP5-3D thermal-hydraulic code and the BISON fuel performance code are tightly coupled for simulation of RIA transients with energy depositions at the Zircaloy-4 cladding failure threshold. 2. Departure from nucleate boiling occurred for all simulated cases. Due to the ductility of fresh fuel, substantial ballooning occurred in most cases. 3. Gas gap conductance, fuel-specific heat capacity, cladding-specific heat capacity, transient pulse width, and film boiling heat transfer were the dominant thermal factors impacting the safety figures of merit at energy depositions.

Critical Heat Flux (CHF)↗

Numerical solutions for laminar and turbulent viscous flow over single and multi-element airfoils using body-fitted coordinate systems

The technique of body-fitted coordinate systems is applied in numerical solutions of the complete time-dependent compressible and incompressible Navier-Stokes equations for laminar flow and to the time-dependent mean turbulent equations closed by modified Kolmogorov hypotheses for turbulent flow. Coordinate lines are automatically concentrated near to the bodies at higher Reynolds number so that accurate resolution of the large gradients near the solid boundaries is achieved. Two-dimensional bodies of arbitrary shapes are treated, the body contour(s) being simply input to the program. The complication of the body shape is thus removed from the problem.

Thompson, J. F.↗

Boundary value problems in magnetohydrodynamics (and fluid dynamics). I - Radiation boundary condition

In view of the difficulties that emerge in numerical simulations of time-dependent phenomena involving MHD waves due to waves' artificial reflection at numerical boundaries, where time-dependent magnetic fields outside those boundaries are unknown, a method previously developed for the determination of radiative boundary conditions is generalized for the treatment of isothermal MHD waves. The illustrative problem thus treated is that of the magnetic braking of a rigidly rotating disk which is embedded in a differentially-rotating fluid. The analytic results obtained are compared with those from numerical simulations with and without the present boundary treatment.

Vanajakshi, T. C.↗

An analysis of CCD images of the coma of Comet Halley

The analysis of selected CCD images of the coma of comet P/Halley is presented. The images were taken using specially designed filters that isolate regions of a comet's spectrum such that only sunlight which has been scattered by the dust in the coma is recorded. The modeling analysis objective is to make use of the skills developed in the development of Monte Carlo particle trajectory models for the distributions of gas species in cometary comae and to use those models as a basis for a new dust coma model. This model will include a self-consistant picture of the time-dependent dusty-gas dynamics of the inner coma and the three-dimensional time-dependent trajectories of the dust particles under the influence of solar gravity and solar radiation pressure in the outer coma. The model is intended to be used as a tool to analyze selected images from the two sets of CCD images with the hope that it will help the understanding of the effects of a number of important processes on the spatial morphology of the observed dust coma. The processes of importance to the observed dust coma include: (1) the dust particle size distribution function; (2) the terminal velocities of various sized dust particles in the inner coma; (3) the radiation scattering properties of dust particles, which are important both in terms of the observe scattered radiation and the radiation pressure acceleration on dust particles; (4) the fragmentation and/or vaporization of dust particles; and (5) the relative importance of CHON and silicate dust particles as they contribute both to the dusty-gasdynamics in the inner coma (that produce the dust particle terminal velocities) and to the observed spatial morphology on the outer dust coma.

Combi, Michael↗

The small amplitude magnetohydrodynamic Riemann problem

The small-amplitude MHD Riemann problem is studied using the Cohen-Kulsrud-Burgers equations. Unlike the coplanar Riemann problem, the evolution of noncoplanar Riemann problems is not self-similar and its flow structures could change in time. But its large-time behavior is very simple and a time-dependent 2 - 3 intermediate shock is always involved for the noncoplanar field rotations. The time-dependent 2 - 3 intermediate shock has a well-defined structure and exists for any degree of field rotation.

Wu, C. C.↗

Disk irradiation and light curves of x ray novae

We study the disk instability and the effect of irradiation on outbursts in the black hole X-ray nova system. In both the optical and soft X-rays, the light curves of several X-ray novae, A0620-00, GH 2000+25, Nova Muscae 1991 (GS 1124-68), and GRO J0422+32, show a main peak, a phase of exponential decline, a secondary maximum or reflare, and a final bump in the late decay followed by a rapid decline. Basic disk thermal limit cycle instabilities can account for the rapid rise and overall decline, but not the reflare and final bump. The rise time of the reflare, about 10 days, is too short to represent a viscous time, so this event is unlikely to be due to increased mass flow from the companion star. We explore the possibility that irradiation by X-rays produced in the inner disk can produce these secondary effects by enhancing the mass flow rate within the disk. Two plausible mechanisms of irradiation of the disk are considered: direct irradiation from the inner hot disk and reflected radiation from a corona or other structure above the disk. Both of these processes will be time dependent in the context of the disk instability model and result in more complex time-dependent behavior of the disk structure. We test both disk instability and mass transfer burst models for the secondary flares in the presence of irradiation.

Kim, S.-W.↗

Integrated Structural/Acoustic Modeling of Heterogeneous Panels

A model for the dynamic response of heterogeneous media is presented. A given medium is discretized into a number of subvolumes, each of which may contain an elastic anisotropic material, void, or fluid, and time-dependent boundary conditions are applied to simulate impact or incident pressure waves. The full time-dependent displacement and stress response throughout the medium is then determined via an explicit solution procedure. The model is applied to simulate the coupled structural/acoustic response of foam core sandwich panels as well as aluminum panels with foam inserts. Emphasis is placed on the acoustic absorption performance of the panels versus weight and the effects of the arrangement of the materials and incident wave frequency.

Bednarcyk, Brett, A.↗

Mathematical Framework Underlying the In Situ Electrochemical Diagnosis of Adsorbed Intermediates Formed during Redox Reactions at Electrode Surfaces

Previously, we have presented an electrochemical technique wherein an electroactive tracer species is employed to probe the rate-limiting factors governing redox reactions at an electrode surface. In this technique, the electrode is first held potentiostatically to facilitate a redox process (step 1), and then the potential is released to open circuit conditions (step 2) so as to monitor the time-dependent re-equilibration of the electrode potential in the presence of the tracer. The time-dependent potential response in step 2 has been shown to contain information about diffusion—limited or desorption—limited steps, enabling in situ probing of the electrochemistry at the electrode surface during step 1. In the present contribution, a theoretical model governing the transient response in step 2 is developed for two limiting cases: diffusion—limited and desorption—limited recovery of the electrode potential. Mathematical modeling shows that, during re-equilibration, the step 2 potential transient corresponding to a case where step 1 involves surface adsorbed species which undergo desorption in step 2 exhibits a much longer time constant than that when re-equilibration occurs under diffusion limitations. The mathematical framework presented herein provides a sound fundamental basis for applying the aforementioned technique to studying adsorption-desorption processes during electrochemistry. Also, technique limitations are presented in light of the modeling findings.

25 ENERGY STORAGE↗