Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time dependent behavior”

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,459 records · Page 81

Topological analysis of computed three-dimensional viscous flow fields

Computed solutions of the time-dependent, Reynolds-averaged Navier-Stokes equations for three dimensional flows having thin shear layers are analyzed using topological concepts. Specific examples include the transonic flow over a body of revolution with conical afterbody at moderate angles of incidence to the free stream. Experimental flow-visualization techniques are simulated graphically to visualize the computed flow. Scalar and vector fluid dynamics properties such as pressure, shear stress, and vorticity on the body surface are presented as topological maps, and their relationship to one another in terms of orientation and singular points is discussed. The extrapolation from these surface topologies toward the understanding of external flow-field behavior is and demonstrated.

Deiwert, G. S.↗

The computation of thermo-chemical nonequilibrium hypersonic flows

Several conceptual designs for vehicles that would fly in the atmosphere at hypersonic speeds have been developed recently. For the proposed flight conditions the air in the shock layer that envelops the body is at a sufficiently high temperature to cause chemical reaction, vibrational excitation, and ionization. However, these processes occur at finite rates which, when coupled with large convection speeds, cause the gas to be removed from thermo-chemical equilibrium. This non-ideal behavior affects the aerothermal loading on the vehicle and has ramifications in its design. A numerical method to solve the equations that describe these types of flows in 2-D was developed. The state of the gas is represented with seven chemical species, a separate vibrational temperature for each diatomic species, an electron translational temperature, and a mass-average translational-rotational temperature for the heavy particles. The equations for this gas model are solved numerically in a fully coupled fashion using an implicit finite volume time-marching technique. Gauss-Seidel line-relaxation is used to reduce the cost of the solution and flux-dependent differencing is employed to maintain stability. The numerical method was tested against several experiments. The calculated bow shock wave detachment on a sphere and two cones was compared to those measured in ground testing facilities. The computed peak electron number density on a sphere-cone was compared to that measured in a flight test. In each case the results from the numerical method were in excellent agreement with experiment. The technique was used to predict the aerothermal loads on an Aeroassisted Orbital Transfer Vehicle including radiative heating. These results indicate that the current physical model of high temperature air is appropriate and that the numerical algorithm is capable of treating this class of flows.

Candler, Graham↗

Laser High-Cycle Thermal Fatigue of Pulse Detonation Engine Combustor Materials Tested

Pulse detonation engines (PDE's) have received increasing attention for future aerospace propulsion applications. Because the PDE is designed for a high-frequency, intermittent detonation combustion process, extremely high gas temperatures and pressures can be realized under the nearly constant-volume combustion environment. The PDE's can potentially achieve higher thermodynamic cycle efficiency and thrust density in comparison to traditional constant-pressure combustion gas turbine engines (ref. 1). However, the development of these engines requires robust design of the engine components that must endure harsh detonation environments. In particular, the detonation combustor chamber, which is designed to sustain and confine the detonation combustion process, will experience high pressure and temperature pulses with very short durations (refs. 2 and 3). Therefore, it is of great importance to evaluate PDE combustor materials and components under simulated engine temperatures and stress conditions in the laboratory. In this study, a high-cycle thermal fatigue test rig was established at the NASA Glenn Research Center using a 1.5-kW CO2 laser. The high-power laser, operating in the pulsed mode, can be controlled at various pulse energy levels and waveform distributions. The enhanced laser pulses can be used to mimic the time-dependent temperature and pressure waves encountered in a pulsed detonation engine. Under the enhanced laser pulse condition, a maximum 7.5-kW peak power with a duration of approximately 0.1 to 0.2 msec (a spike) can be achieved, followed by a plateau region that has about one-fifth of the maximum power level with several milliseconds duration. The laser thermal fatigue rig has also been developed to adopt flat and rotating tubular specimen configurations for the simulated engine tests. More sophisticated laser optic systems can be used to simulate the spatial distributions of the temperature and shock waves in the engine. Pulse laser high-cycle thermal fatigue behavior has been investigated on a flat Haynes 188 alloy specimen, under the test condition of 30-Hz cycle frequency (33-msec pulse period and 10-msec pulse width including a 0.2-msec pulse spike; ref. 4). Temperature distributions were calculated with one-dimensional finite difference models. The calculations show that that the 0.2-msec pulse spike can cause an additional 40 C temperature fluctuation with an interaction depth of 0.08 mm near the specimen surface region. This temperature swing will be superimposed onto the temperature swing of 80 C that is induced by the 10-msec laser pulse near the 0.53-mm-deep surface interaction region.

Zhu, Dong-Ming↗

Spurious Behavior of Shock-Capturing Methods: Problems Containing Stiff Source Terms and Discontinuities

The goal of this paper is to relate numerical dissipations that are inherited in high order shock-capturing schemes with the onset of wrong propagation speed of discontinuities. For pointwise evaluation of the source term, previous studies indicated that the phenomenon of wrong propagation speed of discontinuities is connected with the smearing of the discontinuity caused by the discretization of the advection term. The smearing introduces a nonequilibrium state into the calculation. Thus as soon as a nonequilibrium value is introduced in this manner, the source term turns on and immediately restores equilibrium, while at the same time shifting the discontinuity to a cell boundary. The present study is to show that the degree of wrong propagation speed of discontinuities is highly dependent on the accuracy of the numerical method. The manner in which the smearing of discontinuities is contained by the numerical method and the overall amount of numerical dissipation being employed play major roles. Moreover, employing finite time steps and grid spacings that are below the standard Courant-Friedrich-Levy (CFL) limit on shockcapturing methods for compressible Euler and Navier-Stokes equations containing stiff reacting source terms and discontinuities reveals surprising counter-intuitive results. Unlike non-reacting flows, for stiff reactions with discontinuities, employing a time step and grid spacing that are below the CFL limit (based on the homogeneous part or non-reacting part of the governing equations) does not guarantee a correct solution of the chosen governing equations. Instead, depending on the numerical method, time step and grid spacing, the numerical simulation may lead to (a) the correct solution (within the truncation error of the scheme), (b) a divergent solution, (c) a wrong propagation speed of discontinuities solution or (d) other spurious solutions that are solutions of the discretized counterparts but are not solutions of the governing equations. The present investigation for three very different stiff system cases confirms some of the findings of Lafon & Yee (1996) and LeVeque & Yee (1990) for a model scalar PDE. The findings might shed some light on the reported difficulties in numerical combustion and problems with stiff nonlinear (homogeneous) source terms and discontinuities in general.

Yee, Helen M. C.↗

A New 1D Model for Thermal Mixing and Stratification in Advanced Reactor Transients

Thermal mixing and stratification in large pools and enclosures play a critical role in the safety and performance of pool-type nuclear reactors, particularly during transient scenarios involving significant temperature differences between incoming and bulk coolant. Accurate modeling of these phenomena is essential for predicting system behavior and supporting passive safety features such as natural circulation. Here, this paper presents a new 1D model for thermal mixing and stratification, developed and implemented in the SAM code. The model represents a large pool as 1D coolant jet channels and zero-dimensional bulk pool volumes, enabling the simulation of a wide range of flow configurations, including hot and cold jet interactions, stratified layers, and the influence of complex geometries such as ceilings, free surfaces, and internal obstacles. Heat exchange between jet and pool regions is governed by closure relations calibrated against 3D computational fluid dynamics (CFD) simulations. The model improves upon earlier approaches by incorporating time-dependent jet characteristics and capturing the associated delay effects more accurately. Code-to-code comparisons and validation against experimental data from the Thermal Stratification Test Facility demonstrate the model’s accuracy and flexibility. This work offers two key contributions: (1) an efficient and robust method for simulating thermal mixing and stratification at the system level, eliminating the need for external coupling between system analysis codes and CFD, and (2) a significant enhancement of SAM’s capabilities to analyze thermal stratification phenomena in advanced reactor systems.

SAM↗

Temperature-Dependent Spin-Lattice Relaxation of the Nitrogen-Vacancy Spin Triplet in Diamond

Spin-lattice relaxation within the nitrogen-vacancy (NV) center’s electronic ground-state spin triplet limits its coherence times, and thereby impacts its performance in quantum applications. We report measurements of the relaxation rates on the NV center’s $|m_s = 0 \rangle ↔ |m_s = ± 1 \rangle$ and $|m_s = -1 \rangle ↔ |m_s = + 1\rangle$ transitions as a function of temperature from 9 to 474 K in high-purity samples. We show that the temperature dependencies of the rates are reproduced by an ab initio theory of Raman scattering due to second-order spin-phonon interactions, and we discuss the applicability of the theory to other spin systems. Using a novel analytical model based on these results, we suggest that the high-temperature behavior of NV spin-lattice relaxation is dominated by interactions with two groups of quasilocalized phonons centered at 68.2(17) and 167(12) meV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Long-Term Effects of Humidity on Stainless Steel Pitting in Sea Salt Exposures

Ground 304 stainless steel (SS) samples were exposed to sea salt particles at 35 °C and two relative humidity (RH) levels for durations ranging from 1 week to 2 years. For all exposure times, pit number density and total pit volume at 40% RH were observed to be considerably greater than those at 76% RH. Statistical analysis of distributions of pit populations for both RH conditions showed that pit number density and total pit volume increased rapidly at first but slowed as exposure time increased. Cross-hatched features were observed in the 40% RH pits while ellipsoidal, faceted pits were observed at 76% RH. Optical profilometry indicated that most pits were not hemispherical. X-ray tomography provided evidence of undercutting and fissures. Piecewise curve fitting modeled the 40% RH data closely, predicting that corrosion damage would eventually plateau. However, a similar treatment of the 76% RH data suggested that corrosion damage would continuously increase, which implied that the piecewise power-law fit was limited in its ability to model atmospheric corrosion generally. Based on these observations, the operative mechanisms determining long-term corrosion behavior were hypothesized to be different depending on the RH of exposure.

36 MATERIALS SCIENCE↗

Topological analysis of computed three-dimensional viscous flow fields

Computed solutions of the time-dependent, Reynolds-averaged Navier-Stokes equations for three-dimensional flows having thin shear layers are analyzed, using topological concepts. Specific examples include the transonic flow over a body of revolution with conical afterbody at moderate angles of incidence to the free stream. Experimental flow-visualization techniques are simulated graphically to visualize the computed flow. Scalar and vector fluid dynamic properties, such as pressure, shear stress, and vorticity on the body surface, are presented as topological maps, and their relationship to one another in terms of orientation and singular points is discussed. The extrapolation from these surface topologies toward the understanding of external flow-field behavior is discussed and demonstrated.

Deiwert, G. S.↗

Structure of turbulence at high shear rate

The structure of homogeneous turbulence subject to high shear rate has been investigated by using three-dimensional, time-dependent numerical simulations of the Navier-Stokes equations. This study indicates that high shear rate alone is sufficient for generation of the streaky structures, and that the presence of a solid boundary is not necessary. Evolution of the statistical correlations is examined to determine the effect of high shear rate on the development of anisotropy in turbulence. It is shown that the streamwise fluctuating motions are enhanced so profoundly that a highly anisotropic turbulence state with a 'one-component' velocity field and 'two-component' vorticity field develops asymptotically as total shear increases. Because of high-shear rate, rapid distortion theory predicts remarkably well the anisotropic behavior of the structural quantities.

Lee, Moon Joo↗

Mechanisms of flow control with the unsteady bleed technique

The unsteady bleed technique (a.k.a. internal acoustic forcing) has been shown to be an effective method for control of separation on low Reynolds number airfoils, blunt-end cylinders aligned axially with the flow, cylinders aligned perpendicular to the flow, and forebody geometries at high angles of attack. In many of these investigations, the mechanism for the control has been attributed to enhancement of the shear layer (Kelvin-Helmholtz) instability by the unsteady component of the forcing. However, this is not the only possible mechanism, nor may it be the dominant mechanism under some conditions. In this work it is demonstrated that at least two other mechanisms for flow control are present, and depending on the location and the amplitude of the forcing, these may have significant impact on the flow behavior. Experiments were conducted on a right-circular cylinder with a single unsteady bleed slot aligned along the axis of the cylinder. The effects of forcing frequency, forcing amplitude, and slot location on the azimuthal pressure distribution were studied. The results suggest that a strong vortical structure forms near the unsteady bleed slot when the slot location is upstream of the boundary layer separation point. The structure is unsteady, since it is created by the unsteady forcing. The 'vortex' generates a sizeable pressure spike (C(sub p) = -3.0) in the time-averaged pressure field immediately downstream of the slot. In addition to the pressure spike, the boundary layer separation location moves farther downstream when the forcing is activated. Delay of the separation is believed to be a result of enhancing the Kelvin-Helmholtz instability. When forcing is applied in a quiescent wind tunnel, a weak low-pressure region forms near the slot that is purely the result of the second-order streaming effect.

Williams, D. R.↗

Combustion of Aerosolized Metal Particles in Microgravity

The combustion behavior and interaction effects of multiple metal particles are addressed in this project. The microgravity environment presents a unique opportunity to create an "aerosol" consisting of relatively large, 50-300 micron diameter particles, so that both cloud flame structure and individual particle combustion behavior can be characterized simultaneously. The microgravity experiments are conducted using the 2.2 s NASA Lewis Research Center drop tower. Each test includes aerosolizing metal particles under microgravity, a delay required to produce a steady aerosol, and igniting such an "aerosol" at constant pressure using a hot wire igniter. The flame structure and details of individual particle combustion are visualized using both high speed movie and regular speed video cameras. Combustion products are collected and analyzed after the experiment. A detailed description of the experimental apparatus is given elsewhere. Microgravity experiments with magnesium particle aerosols have been conducted and experiments with zirconium particles are currently in progress. A numerical, time dependent model of the flame propagation in magnesium aerosol is being developed.

Dreizin, E. L.↗

Spectrally Dependent Radiometric Measurement Errors Due to Ccd Serial Pixel-to-Pixel Readout Interference in the Ocean Color Instrument of the NASA PACE Mission

The Ocean Color Instrument on NASA’s PACE mission is a 322 887 nm hyperspectral imager with 1 km x 1 km nadir spatial resolution and 5 nm spectral resolution utilizing charge-coupled devices (CCDs) operating in Time Delay Integration (TDI) mode where each TDI column represents a different wavelength in 0.625 nm increments. After TDI, the charge is moved into serial output pixels and read out. The spatial resolution requires an 8.5 MHz readout rate. This only allows 59 ns for the CCD reset and video to be asserted and settled before sampling. The response exhibits serial pixel-to-pixel readout interference due to the lack of full settling. Each serial pixel value has a dependence on the value of the preceding pixel value. This leads to a spectrally dependent radiometric measurement error of up to 0.3 %. We explain the operation of the detection system, the behavior of the interference, and show the resulting measurement error based on data from both ground testing and on-orbit characterization.

Remote Sensing↗

Spectrally Dependent Radiometric Measurement Errors Due to CCD Serial Pixel-to-Pixel Readout Interference in the Ocean Color Instrument of the NASA PACE Mission

The Ocean Color Instrument on NASA’s PACE mission is a 322 887 nm hyperspectral imager with 1 km x 1 km nadir spatial resolution and 5 nm spectral resolution utilizing charge-coupled devices (CCDs) operating in Time Delay Integration (TDI) mode where each TDI column represents a different wavelength in 0.625 nm increments. After TDI, the charge is moved into serial output pixels and read out. The spatial resolution requires an 8.5 MHz readout rate. This only allows 59 ns for the CCD reset and video to be asserted and settled before sampling. The response exhibits serial pixel-to-pixel readout interference due to the lack of full settling. Each serial pixel value has a dependence on the value of the preceding pixel value. This leads to a spectrally dependent radiometric measurement error of up to 0.3 %. We explain the operation of the detection system, the behavior of the interference, and show the resulting measurement error based on data from both ground testing and on-orbit characterization.

Remote Sensing↗

Assessment of the State of Communication Delay Research in Preparation for Missions Beyond Low Earth Orbit

NASA’s mission-operations paradigm, established during Project Mercury and minimally evolving through the Apollo Program, Space Shuttle Program, and ISS missions, has primarily depended on real-time support from a ground team of experts. This ground team has served as the safety net for crewed spaceflight missions over the past 60 years, managing the combined state of the mission, vehicle, and crew. However, this operational paradigm, which has seen little change in its Human-Systems Integration Architecture (HSIA), will face challenges during long-duration exploration missions beyond low Earth orbit. Lunar missions may experience one-way communication latencies ranging from 3 to 14 seconds, while Mars missions will encounter up to a 44-minute round-trip latency at a maximum distance from Earth. Communication delays negatively impact the behavioral health and performance of individuals and crews operating across the multi-team space-to-ground system. Previous research indicates that increased isolation and the challenges posed by delayed communication lead to heightened stress, frustration, adverse behavioral symptoms, and reduced individual performance. NASA has investigated the nature and effectiveness of managing this shift in complex operations since early 2000, but additional research is necessary to assess the issues associated with communication latencies and identify effective countermeasures. This work aimed to examine the evidence from 20 years of research on comm delays. We systematically reviewed the past 20 years of communication delay literature, focusing on how this research aligns with relevant, high-priority needs in the HSIA, Team, and BMed Risks. This type of panoramic review has not been done since 2013; a current, aggregated picture of what has been studied and how is needed to inform additional studies and mitigation development. Additional perspectives were gained from interviews with the research participants, operations experts, and communication delay researchers. Results from this effort will help characterize the risk posed by communication latency for upcoming Artemis missions and point toward potential mitigations. The final results will be presented. For the literature review, we formed a search term seed set drawing from languages used in two papers representative of the previous state of knowledge and supplemented it with additional search terms. We conducted literature searches on Google Scholar, PubMed, and Web of Science and searched NASA project archives for output from NASA-sponsored research. This search yielded 150 papers, of which 48 were relevant. SME interviews provided perspective on key communication delay issues and concerns and what is being worked on to mitigate those issues. Preliminary results of the literature review have been presented. Final results to be presented include the results of the literature review and findings from the interviews.

literature review↗

Optimization of the light detection system of the ICARUS detector

The ICARUS detector, a key component of the Short Baseline Neutrino (SBN) Program at Fermi National Acelerator Laboratory (FNAL), is a 600-ton Liquid Argon Time Projection Chamber (LArTPC) equipped with a Light Detection System (LDS) that uses 360 Hamamatsu R5912-MOD 8-inch photomultiplier tubes (PMTs), specifically designed to operate under cryogenic conditions ($\sim 87 \ K$). These PMTs feed the trigger signal to the readout, improve the spatial and timing resolution of the events, and contribute to cosmic rays mitigation. During operation at FNAL, a progressive degradation in the PMT gain was observed. We developed an experimental setup to investigate the temperature dependence of PMT performance. Gain measurements were carried out from room temperature to $-70 ^\circ C$ using an environmental chamber. The results show that, while the PMTs exhibit stable performance at room temperature, a significant and irreversible reduction in gain emerges at lower temperatures. Al though $-70 ^\circ C$ remains above the liquid argon temperatures, the trend clearly reveals a gain-sensitive degradation mechanism. A simplified physical model was developed to reproduce and interpret the observed behavior. Based on these findings, a series of mitigation strategies were implemented in the ICARUS detector to preserve PMT performance and ensure reliable operation under cryogenic conditions.

Saia, C. [Catania Astrophys. Observ.]↗

Single-ion properties of the transverse-field Ising model material CoNb 2 O 6

We report CoNb 2 O 6 is one of the few materials that is known to approximate the one-dimensional transverse-field Ising model (1D-TFIM) near its quantum critical point. It has been inferred that Co 2+ acts as a pseudospin 1/2 with anisotropic exchange interactions that are largely Ising-like, enabling the realization of the TFIM. However, the behavior of CoNb 2 O 6 is known to diverge from the ideal TFIM under transverse magnetic fields that are far from the quantum critical point, requiring the consideration of additional anisotropic, bond-dependent (Kitaev-like) terms in the microscopic pseudospin-1/2 Hamiltonian. These terms are expected to be controlled in part by single-ion physics, namely the wave function for the pseudospin-1/2 angular momentum doublet. Here, we present the results of both inelastic neutron scattering measurements and electron paramagnetic resonance spectroscopy on CoNb 2 O 6 , which elucidate the single-ion physics of Co 2+ in CoNb 2 O 6 for the first time. We find that the system is well-described by an intermediate spin-orbit coupled Hamiltonian, and the ground state is a well-isolated Kramers doublet with an anisotropic g tensor. We provide the approximate wave functions for this doublet, which we expect will be useful in theoretical investigations of the anisotropic exchange interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Field-Dependent Ionic Conductivities from Generalized Fluctuation-Dissipation Relations

We derive a relationship for the electric field dependent ionic conductivity in terms of fluctuations of time integrated microscopic variables. We demonstrate this formalism with molecular dynamics simulations of solutions of differing ionic strength with implicit solvent conditions and molten salts. These calculations are aided by a novel nonequilibrium statistical reweighting scheme that allows for the conductivity to be computed as a continuous function of the applied field. In strong electrolytes, we find the fluctuations of the ionic current are Gaussian, and subsequently, the conductivity is constant with applied field. In weaker electrolytes and molten salts, we find the fluctuations of the ionic current are strongly non-Gaussian, and the conductivity increases with applied field. Furthermore, this nonlinear behavior, known phenomenologically for dilute electrolytes as the Onsager-Wien effect, is general and results from the suppression of ionic correlations at large applied fields, as we elucidate through both dynamic and static correlations within nonequilibrium steady states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Superconducting Magnet Qualification Methodology for the Material Plasma Exposure Experiment

The Material Plasma Exposure eXperiment (MPEX), currently under design, is a new linear plasma device to advance the understanding of plasma-material interactions through the generation and delivery of plasmas as they are expected in future fusion reactor divertors. MPEX will be a steady-state device to study high-fluence exposures of plasma-facing materials and components. The requirements for the magnetic field at the target and the heating stages make the application of superconducting coils necessary. Conceptual designs for the superconducting magnets have been developed, and multiple cryostats with warm bore diameters of either 65 cm or 156 cm are envisioned to facilitate their integrated and timely assembly with other systems such as vacuum, water cooling, and RF power. Although design, fabrication, and testing for the magnets as stand-alone units are straightforward, challenges will arise during the integration of the system. Two different field profiles will be used during operation. The magnetic field where the electron cyclotron heating occurs needs to operate at both 1.25 and 2.5 T. It is critical that the magnets all share the same magnetic axis and alignment. The mutual inductance between cryostats will affect the quench behavior of the system. Also, cryostat-to-cryostat forces can be as large as 700 kN, and the magnitude and direction will change depending on which coils are energized. The design of the system must take those characteristics into account along with the quench scenarios. This paper describes the qualification approach that will be used to determine whether stand-alone tests can be used to ensure the success of the integrated system. Fiducials will be used to define the location of the magnetic axis for each cryostat to ensure proper alignment. Quench tests of a single magnet will be performed at a current above the normal operating current to account for additional stored energy from the mutual inductance to adjacent cryostats. Also, a 1018 steel plate will be mounted on either end of a cryostat to simulate the cryostat-to-cryostat forces. Requirements for the size and location of the steel plates are described.

Burkhardt, Earle↗