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

Constraints on the IBEX Ribbon’s Origin from Its Evolution over a Solar Cycle

Abstract In 2009, the Interstellar Boundary Explorer (IBEX) discovered a narrow “ribbon” of energetic neutral atom emissions across the sky with properties correlated with the solar wind latitudinal structure and the interstellar magnetic field draped around the heliosphere. It is widely believed that the ribbon is formed from the escape of heliospheric ENAs into the local interstellar medium and their eventual return as secondary ENAs. However, there is no consensus on the rate of pitch angle scattering of these PUIs before they become secondary ENAs. We test two opposing limits of scattering rates (“weak” versus “strong”) by solving a time-dependent model of the ribbon that evolves with the solar cycle, and we compare them to IBEX observations over 2009–2019. First, we find that both models qualitatively reproduce the evolution of IBEX fluxes for most of the data set, with a few exceptions, although the strong (or “spatial retention”) scattering model greatly underestimates the observed fluxes. Regardless, time dependence of fluxes cannot distinguish these models. Second, the ribbon’s geometric properties, i.e., its center and radius, are significantly different between the models. The spatial retention model reproduces the observed ribbon centers as a function of energy and time slightly better than the weak scattering model, and the spatial retention model reproduces the observed ribbon radius over energy and time almost perfectly, whereas the weak scattering model compares poorly. Our analysis favors the spatial retention mechanism as the source of the IBEX ribbon, but it requires modification to increase the flux of ENAs observed at 1 au.

Zirnstein, E. J. (ORCID:0000000172400618)↗
NASA NTRS1982

Large Prandtl number finite-amplitude thermal convection with Maxwell viscoelasticity

It has long been known that the earth behaves viscoelastically. Viscoelasticity may be of importance in two aspects of mantle convection, including time-dependent behavior and local storage of recoverable work. The present investigation makes use of thermal convection in a box as a prototype of mantle flow. It is demonstrated that recoverable work can be important to the local mechanical energy balance in the descending lithosphere. It is shown that, even when assuming large viscoelastic parameters, an inherent time-dependence of viscoelastic convection appears only in local exchanges of mechanical energy. There is no strong exchange between buoyant potential energy and recoverable strain energy in the Rayleigh number range investigated. The investigation is mainly concerned with viscoelastic effects occurring on a buoyant time scale. It is found that viscoelastic effects have a negligible influence on the long term thermal energetics of mantle convection.

Ivins, E. R.↗
DOE OSTI2020

Quantum Spectral Methods for Differential Equations

Recently developed quantum algorithms address computational challenges in numerical analysis by performing linear algebra in Hilbert space. Such algorithms can produce a quantum state proportional to the solution of a d-dimensional system of linear equations or linear differential equations with complexity poly(logd). While several of these algorithms approximate the solution to within ϵ with complexity poly(log(1/ϵ)), no such algorithm was previously known for differential equations with time-dependent coefficients. In this work, we develop a quantum algorithm for linear ordinary differential equations based on so-called spectral methods, an alternative to finite difference methods that approximates the solution globally. Using this approach, we give a quantum algorithm for time-dependent initial and boundary value problems with complexity poly(logd, log(1/ϵ)).

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1996

Aerodynamic Characteristics of an Aerospace Vehicle During a Subsonic Pitch-Over Maneuver

Time-dependent CFD has been used to predict aerospace vehicle aerodynamics during a subsonic rotation maneuver. The inviscid 3D3U code is employed to solve the 3-D unsteady flow field using an unstructured grid of tetrahedra. As this application represents a challenge to time-dependent CFD, observations concerning spatial and temporal resolution are included. It is shown that even for a benign rotation rate, unsteady aerodynamic effects are significant during the maneuver. Possibly more significant, however, the rotation maneuver creates ow asymmetries leading to yawing moment, rolling moment, and side force which are not present in the quasi-steady case. A series of steady solutions at discrete points in the maneuver are also computed for comparison with wind tunnel measurements and as a means of quantifying unsteady effects.

Kleb, William L.↗
NASA NTRS1989

Simultaneous stochastic inversion for geomagnetic main field and secular variation. II - 1820-1980

With the aim of producing readable time-dependent maps of the geomagnetic field at the core-mantle boundary, the method of simultaneous stochastic inversion for the geomagnetic main field and secular variation, described by Bloxham (1987), was applied to survey data from the period 1820-1980 to yield two time-dependent geomagnetic-field models, one for the period 1900-1980 and the other for 1820-1900. Particular consideration was given to the effect of crustal fields on observations. It was found that the existing methods of accounting for these fields as sources of random noise are inadequate in two circumstances: (1) when sequences of measurements are made at one particular site, and (2) for measurements made at satellite altitude. The present model shows many of the features in the earth's magnetic field at the core-mantle boundary described by Bloxham and Gubbins (1985) and supports many of their earlier conclusions.

Bloxham, Jeremy↗
NASA NTRS2000

Growth of Matrix Cracks During Intermediate Temperature Stress Rupture of a SiC/SiC Composite in Air

The crack density of woven Hi-Nicalon(sup TM) (Nippon Carbon, Japan) fiber, BN interphase, melt-infiltrated SiC matrix composites was determined for specimens subjected to tensile stress rupture at 815 C. A significant amount of matrix cracking occurs due to the growth of fiber-bridged microcracks even at stresses below the run-out condition. This increased cracking corresponded to time dependent strain accumulation and acoustic emission activity during the constant load test. However, the portion of the rupture specimens subjected to cooler temperatures (< 600 C than the hot section had significantly lower crack densities compared to the hotter regions. From the acoustic emission and time dependent strain data it can be inferred that most of the matrix crack growth occurred within the first few hours of the tensile rupture experiment. The crack growth was attributed to an interphase recession mechanism that is enhanced by the presence of a thin carbon layer between the fiber and the matrix as a result of the composite fabrication process. One important consequence of matrix crack growth at the lower stresses is poor retained strength at room temperature for specimens that did not fail.

Morscher, Gregory N.↗
DOE OSTI2024

Nonlinear control of the minimum safety factor in tokamaks by optimal allocation of spatially moving electron cyclotron current drive

The minimum value of the safety factor profile is related to the magnetohydrodynamic (MHD) stability of the plasma confined in a tokamak. Therefore, active control of the minimum safety factor may mitigate MHD instabilities that can degrade or even terminate plasma confinement. Typically, in most tokamak scenarios, the minimum safety factor evolves spatially with time, i.e., the location at which the safety factor achieves the minimum value changes with time. In addition to the inherent nonlinearities in the minimum safety factor evolution, its spatial variation makes the control design challenging. In particular, complexity in control design may arise from the need for time-dependent nonlinear models that account for spatial variation of the minimum safety factor. Furthermore, the minimum safety factor may drift to locations where the actuator authority is low. The problem of minimum safety factor control with target location tracking and moving electron cyclotron current drive (ECCD) is addressed in this work. A nonlinear time-dependent model that incorporates the spatial variation of the minimum safety factor is presented. A nonlinear controller based on optimal feedback linearization is developed to track a target minimum safety factor. The proposed controller treats the ECCD position as a controllable variable. In other words, the controller prescribes the ECCD position (in addition to the non-inductive powers) in real time based on an optimal criterion that is defined a priori. This work also presents the steps necessary to integrate the minimum safety factor controller with a total energy controller to achieve multiple control objectives simultaneously. In conclusion, the proposed integrated control algorithm is tested using nonlinear simulations in the Control Oriented Transport SIMulator (COTSIM) for a DIII-D tokamak scenario.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Production of neutron-rich heavy nuclei in deep-inelastic 208 Pb + 208 Pb collisions within the stochastic mean-field theory

In deep-inelastic collisions of heavy nuclei, reaction products with a wide range of mass and charge are produced. Such collisions been considered as a possible way to produce superheavy nuclei, as an alternative to fusion reactions. To provide reliable theoretical predictions, it is desired to develop microscopic approaches that correctly and accurately describe nucleon transfer processes in dissipative collisions of heavy nuclei. The purpose of the present work is (1) to investigate the mechanism of nucleon transfers in dissipative collisions of two heavy nuclei, and (2) to explore possible pathways to produce neutron-rich heavy nuclei, through detailed theoretical analyses of fluctuations and correlations in nucleon transfers in 208 Pb + 208 Pb reactions. Three-dimensional time-dependent Hartree-Fock (TDHF) calculations are performed for the collisions of 208 Pb + 208 Pb at 𝐸 c.m. = 832, 936, and 1040 MeV, using the Skyrme SLy4d energy density functional. To calculate fluctuations and correlations in nucleon transfers, we employ the stochastic mean-field (SMF) theory, and the results are compared with another theoretical framework currently available, the time-dependent random phase approximation (TDRPA). Primary and secondary production cross sections are calculated with the SMF theory combined with a statistical model, GEMINI ++ . Using information of nucleon flow across a neck of colliding nuclei in TDHF calculations, we solve quantal diffusion equations for fluctuations and correlations in nucleon transfers based on the SMF theory. From the SMF calculations, we obtain the time evolution of diffusion coefficients as well as fluctuations and correlations in nucleon transfers for a range of initial orbital angular momenta. We compare the results of the SMF calculations with those of TDRPA, showing that TDRPA tends to predict substantially larger fluctuations and correlations in strongly damped collisions of heavy nuclei, which exhibit complex initial angular momentum dependence, while the SMF results provide almost constant (stable) values. Using the obtained fluctuations and correlations, we calculate primary and secondary production cross sections for the 208 Pb + 208 Pb collisions. From the results, we find that both lighter and heavier reaction products as compared to 208 Pb are produced for a wide region in the 𝑁−𝑍 plane as primary products, thanks to the quantal diffusion mechanism in the dissipative collisions. However, we show that cross sections for production of heavy nuclei with 𝑍 ≳ 90 or 𝑁 ≳ 135 are washed out due to secondary particle evaporation and/or fission processes. On the other hand, we find that there remain sizable cross sections for production of neutron-rich nuclei along 𝑁 = 126 with 𝑍< 82, even after secondary disintegration processes. We demonstrate that the secondary production cross sections depend weakly on incident energies, but lower (higher) energy is slightly preferred for production of nuclei with smaller (larger) atomic numbers as compared to 𝑍 = 82. Here, based on the microscopic SMF calculations, it has been shown that deep-inelastic collisions of heavy nuclei, such as 208 Pb + 208 Pb examined in this study, can be a promising means to produce neutron-rich heavy nuclei along 𝑁=126. Discrepancies between the SMF and TDRPA approaches are left unsolved for future investigations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The evolution of HCO + in molecular clouds using a novel chemical post-processing algorithm

Modelling the chemistry of molecular clouds is critical to accurately simulating their evolution. To reduce computational cost, 3D simulations generally restrict their chemistry to species with strong heating and cooling effects. Time-dependent information about the evolution of other species is therefore often neglected. We address this gap by post-processing tracer particles in the SILCC-Zoom molecular cloud simulations. Using a chemical network of 39 species and 301 reactions (including freeze-out of CO and H 2 O) and a novel algorithm to reconstruct a density grid from sparse tracer particle data, we produce time-dependent density distributions for various species. We focus upon the evolution of HCO + , which is a critical formation reactant of CO but is not typically modelled on the fly. We find that ∼ 90 per cent of the HCO + content of the cold molecular gas forms in situ around n HCO + ∼ 10 3 –10 4 cm −3 , over a time-scale of approximately 1 Myr. The remaining ∼ 10 per cent forms at high extinction sites, with minimal turbulent mixing out into the less dense gas. We further show that the dominant HCO + formation pathway is dependent on the visual extinction, with the reaction H 3 + + CO contributing 90 per cent of the total HCO + production above A V, 3D = 3. We produce the very first maps of the HCO + column density, N(HCO + ), and show that it reaches values as high as 10 15 cm −2 . We find that 50 per cent of the HCO + mass is located within AV ∼ 10–30 in a density range of 10 3.5 –10 4.5 cm −3 . Our maps of N(HCO + ) are shown to be in good agreement with recent observations of the W49A star-forming region.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1976

Diffraction of a shock wave by a compression corner. I - Regular reflection

The unsteady, two-dimensional flowfield resulting from the interaction of a moving planar shock wave with a compression corner is determined using a second-order, discontinuity-fitting, finite-difference approach. The time-dependent Euler equations are transformed to normalize the distance between the body and peripheral shock and to include the existing self-similar property of the flow. The resulting set of partial differential equations in conservation-law form is then solved in a time-dependent fashion using MacCormack's scheme. The vortical singularity, which lies on the body surface, and the single reflected shock are both treated as discontinuities in the numerical procedure. The results of the numerical simulation compare quite favorably with existing experimental interferograms and yield better flowfield resolution than previous first-order, shock-capturing, numerical solutions.

Kutler, P.↗
NASA NTRS1983

Initial-value problems for Rossby waves in a shear flow with critical level

The time-dependent evolution of sheared Rossby waves starting from an initial disturbance is studied for the simple case in which the shear is uniform. The uniform-shear assumption allows explicit solutions to be obtained which are useful in addressing the issue of the long-time asymptotic approach to normal modes and in assessing the relative importance of viscosity, nonlinearity and time-dependence in the evolution of Rossby waves in the presence of critical layers.

Tung, K. K.↗
NASA NTRS1995

Effect of Gravity on the Mammalian Cell Deformation

The effect of human cell immersed in culture liquid under a micro-gravity environment has been investigated. The study is based on the numerical simulation of the configuration of human cell affected by the time dependent variation of gravity acceleration ranging from 10(exp -3) to 2 g(sub o) (g(sub o) = 9.81 m/s(exp 2)) in 15 seconds. Both the free floating cell and the cell contacted to the upper and lower inclined walls imposed by the time-dependent reduced gravity acceleration are considered in this study. The results show that the cell configuration changes from spherical to horizontally elongated ellipsoid for both the free floating cell and the cell sitting on the lower inclined wall while the cell configuration varies from spherical to vertically elongated ellipsoid for the cell hanging to the upper inclined wall when the gravity acceleration increases. Experimental observations, carried out of human cells exposed to the variation of gravity levels, show that the results of experimental observations agree exactly with the theoretical model computation described in this paper. These results sre significant for humans exposed to the micro-gravity environment.

Hung, R. J.↗
NASA NTRS2002

Local Discontinuous Galerkin Methods for Partial Differential Equations with Higher Order Derivatives

In this paper we review the existing and develop new continuous Galerkin methods for solving time dependent partial differential equations with higher order derivatives in one and multiple space dimensions. We review local discontinuous Galerkin methods for convection diffusion equations involving second derivatives and for KdV type equations involving third derivatives. We then develop new local discontinuous Galerkin methods for the time dependent bi-harmonic type equations involving fourth derivatives, and partial differential equations involving fifth derivatives. For these new methods we present correct interface numerical fluxes and prove L(exp 2) stability for general nonlinear problems. Preliminary numerical examples are shown to illustrate these methods. Finally, we present new results on a post-processing technique, originally designed for methods with good negative-order error estimates, on the local discontinuous Galerkin methods applied to equations with higher derivatives. Numerical experiments show that this technique works as well for the new higher derivative cases, in effectively doubling the rate of convergence with negligible additional computational cost, for linear as well as some nonlinear problems, with a local uniform mesh.

Yan, Jue↗
DOE OSTI2026

Evaluation of the beam-induced depolarization of the HJET target at the EIC

The Polarized Atomic Hydrogen Gas Jet Target (HJET) has played a central role in the absolute calibration of proton beam polarization at RHIC and is foreseen as a key component of the hadron polarimetry program at the future Electron–Ion Collider (EIC). The substantially higher beam current, reduced bunch spacing, and shorter bunch length planned for EIC operation motivate a careful reassessment of possible beam-induced depolarization of the jet target. In this paper, the depolarization of ground-state hydrogen atoms caused by the time-dependent magnetic field of the circulating polarized proton beam is quantitatively evaluated. The hydrogen atom is treated as a four-level hyperfine system in a holding magnetic field, and transitions driven by harmonic components of the bunch-induced magnetic field are analyzed using time-dependent quantum-mechanical evolution along atomic trajectories. Numerical tracking of hydrogen atoms through the beam region is performed using nominal EIC beam parameters. It is shown that, for a holding field of 120 mT (as used at RHIC), the resulting depolarization of the jet target at the EIC is negligibly small, ≲ 0.01 %, and well below the level relevant for EIC polarization accuracy requirements. The stability of this result with respect to plausible variations of the EIC proton beam parameters is also evaluated. In addition, possible effects under alternative experimental conditions are also examined.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Long-term survival and second malignant tumor prediction in pediatric, adolescent, and young adult cancer survivors using Random Survival Forests: a SEER analysis

Abstract Survival and second malignancy prediction models can aid clinical decision making. Most commonly, survival analysis studies are performed using traditional proportional hazards models, which require strong assumptions and can lead to biased estimates if violated. Therefore, this study aims to implement an alternative, machine learning (ML) model for survival analysis: Random Survival Forest (RSF). In this study, RSFs were built using the U.S. Surveillance Epidemiology and End Results to (1) predict 30-year survival in pediatric, adolescent, and young adult cancer survivors; and (2) predict risk and site of a second tumor within 30 years of the first tumor diagnosis in these age groups. The final RSF model for pediatric, adolescent, and young adult survival has an average Concordance index (C-index) of 92.9%, 94.2%, and 94.4% and average time-dependent area under the receiver operating characteristic curve (AUC) at 30-years since first diagnosis of 90.8%, 93.6%, 96.1% respectively. The final RSF model for pediatric, adolescent, and young adult second malignancy has an average C-index of 86.8%, 85.2%, and 88.6% and average time-dependent AUC at 30-years since first diagnosis of 76.5%, 88.1%, and 99.0% respectively. This study suggests the robustness and potential clinical value of ML models to alleviate physician burden by quickly identifying highest risk individuals.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1979

Global magnetohydrodynamic simulation of the two-dimensional magnetosphere

The time-dependent magnetohydrodynamic interaction of the solar wind with a two-dimensional dipole magnetic field has been simulated using a novel Lagrangian particle type of MHD code that can treat local low density or vacuum regions without numerical instability. This enables one to simulate the time-dependent magnetic tail. When the solar wind field is southward, a magnetic field line topology consistent with Dungey's model emerges in steady state. The tail, however, is short, and the x-points are only slightly shifted from their vacuum locations, because of strong numerical resistivity. Different configurations resulting from different relative orientations of the solar wind magnetic field and dipole axis are also presented. While the magnetic field is relatively steady, the density and flow in the magnetosheath are turbulent, as are the bow shock and magnetopause; the Kelvin-Helmholtz instability may account for these phenomena. We also model a 'substorm' as the passage of a rotational discontinuity in the solar wind over the dipole. Both 90 and 180 deg shifts to a southward solar wind field cause a violent readjustment of the magnetic tail which eventually settles down to the Dungey configuration.

Leboeuf, J. N.↗
NASA NTRS1982

Coronal closed structures. IV - Hydrodynamical stability and response to heating perturbations

The response of magnetically confined atmospheres to perturbations in the temperature and density distribution, and the local heating rate by means of a one-dimensional time-dependent hydrodynamical code, which incorporates the full energy, momentum and mass conservation equations is studied. These studies extend the linear instability analysis of Habbal and Rosner (1979) into the finite-amplitude regime, and generalize the confined atmosphere models of Serio et al., to the time-dependent domain. The results show that closed coronal atmospheres are stable against finite-amplitude perturbations if the chromospheric response is taken into account; and observed correlated increases in coronal density and temperature can only be achieved under quiescent conditions by increasing the heat deposition rate relatively more in the chromosphere than in the corona.

Peres, G.↗
NASA NTRS1986

Bifurcations in unsteady aerodynamics

Nonlinear algebraic functional expansions are used to create a form for the unsteady aerodynamic response that is consistent with solutions of the time dependent Navier-Stokes equations. An enumeration of means of invalidating Frechet differentiability of the aerodynamic response, one of which is aerodynamic bifurcation, is proposed as a way of classifying steady and unsteady aerodynamic phenomena that are important in flight dynamics applications. Accomodating bifurcation phenomena involving time dependent equilibrium states within a mathematical model of the aerodynamic response raises an issue of memory effects that becomes more important with each successive bifurcation.

Tobak, M.↗