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 865 records · Page 48

Brief Communication: A Simplified Approach to Transient Convective Droplet Evaporation and Burning

Empirical correlations for evaporation rates from single fuel droplets have existed since the 1930s. These correlations, which will be referred to in this article as Froessling/Ranz-Marshall types of correlations, are appropriate to the special cases of steady-state evaporation in the absence of chemical reaction. In a previous article by the authors, the quasi-steady evaporation and burning processes associated with a fuel drop in a convective environment are examined through a droplet model based on the boundary layer approach. For droplet Reynolds numbers of practical interest, this model produces very reasonable steady state as well as quasi-time-dependent droplet simulations, requiring relatively short computational times and yielding good agreement with the above-mentioned empirical correlations. The steady-state case, however, is usually relevant to practical combustor situations only when the drop has reached a nearly uniform temperature since the heating process of the drop cannot be considered to be quasi-steady. In the present study, the transient heating process of the droplet interior during evaporation and/or burning is taken into account, and thus calculations pertaining to the entire life-time of the droplet are carried out. It is of particular interest here to obtain simplified correlations to describe the transient behavior of evaporating and burning droplets; these may be incorporated with greater ease into spray calculations. Accordingly, we have chosen to use stagnation conditions in the present model in a modification of the Froessling/Ranz-Marshall correlations. These modified correlations, incorporating an effective transfer number, produce a fairly accurate representation of droplet evaporation and burning, while requiring only one tenth the computational effort used in a full boundary layer solution.

Madooglu, K.↗

Development and Implementation of Non-Newtonian Rheology Into the Generalized Fluid System Simulation Program (GFSSP)

One of the goals of this program was to develop the experimental and analytical/computational tools required to predict the flow of non-Newtonian fluids through the various system components of a propulsion system: pipes, valves, pumps etc. To achieve this goal we selected to augment the capabilities of NASA's Generalized Fluid System Simulation Program (GFSSP) software. GFSSP is a general-purpose computer program designed to calculate steady state and transient pressure and flow distributions in a complex fluid network. While the current version of the GFSSP code is able to handle various systems components the implicit assumption in the code is that the fluids in the system are Newtonian. To extend the capability of the code to non-Newtonian fluids, such as silica gelled fuels and oxidizers, modifications to the momentum equations of the code have been performed. We have successfully implemented in GFSSP flow equations for fluids with power law behavior. The implementation of the power law fluid behavior into the GFSSP code depends on knowledge of the two fluid coefficients, n and K. The determination of these parameters for the silica gels used in this program was performed experimentally. The n and K parameters for silica water gels were determined experimentally at CFDRC's Special Projects Laboratory, with a constant shear rate capillary viscometer. Batches of 8:1 (by weight) water-silica gel were mixed using CFDRC s 10-gallon gelled propellant mixer. Prior to testing the gel was allowed to rest in the rheometer tank for at least twelve hours to ensure that the delicate structure of the gel had sufficient time to reform. During the tests silica gel was pressure fed and discharged through stainless steel pipes ranging from 1", to 36", in length and three diameters; 0.0237", 0.032", and 0.047". The data collected in these tests included pressure at tube entrance and volumetric flowrate. From these data the uncorrected shear rate, shear stress, residence time, and viscosity were evaluated using formulae for non-Newtonian, power law fluids. The maximum shear rates (corrected for entrance effects) obtained in the rheometer with the current setup were in the 150,000 to 170,000sec- range. GFSSP simulations were performed with a flow circuit simulating the capillary rheometer and using Power Law gel viscosity coefficients from the experimental data. The agreement between the experimental data and the simulated flow curves was within +/-4% given quality entrance effect data.

DiSalvo, Roberto↗

Recent Studies of the Behavior of the Sun's White-Light Corona Over Time

Predictions of upcoming solar cycles are often related to the nature and dynamics of the Sun's polar magnetic field and its influence on the corona. For the past 30 years we have a more-or-less continuous record of the Sun's white-light corona from groundbased and spacebased coronagraphs. Over that interval, the large scale features of the corona have varied in what we now consider a 'predictable' fashion--complex, showing multiple streamers at all latitudes during solar activity maximum; and a simple dipolar shape aligned with the rotational pole during solar minimum. Over the past three decades the white-light corona appears to be a better indicator of 'true' solar minimum than sunspot number since sunspots disappear for months (even years) at solar minimum. Since almost all predictions of the timing of the next solar maximum depend on the timing of solar minimum, the white-light corona is a potentially important observational discriminator for future predictors. In this contribution we describe recent work quantifying the large-scale appearance of the Sun's corona to correlate it with the sunspot record, especially around solar minimum. These three decades can be expanded with the HAO archive of eclipse photographs which, although sparse compared to the coronagraphic coverage, extends back to 1869. A more extensive understanding of this proxy would give researchers confidence in using the white-light corona as an indicator of solar minimum conditions.

SaintCyr, O. C.↗

Inelastic deformation of metal matrix composites: Plasticity and damage mechanisms, part 2

The inelastic deformation mechanisms for the SiC (SCS-6)/Ti-15-3 system were studied at 538 C (1000 F) using a combination of mechanical measurements and detailed microstructural examinations. The objectives were to evaluate the contributions of plasticity and damage to the overall MMC response, and to compare the room temperature and elevated temperature deformation behaviors. Four different laminates were studied: (0)8, (90)8,(+ or -45)2s, and (0/90)2s, with the primary emphasis on the unidirectional (0)8, and (90)8 systems. The elevated temperature responses were similar to those at room temperature, involving a two-stage elastic-plastic type of response for the (0)8 system, and a characteristic three-stage deformation response for the (90)8 and (+ or -45)2s systems. The primary effects of elevated temperatures included: (1) reduction in the 'yield' and failure strengths; (2) plasticity through diffused slip rather than concentrated planar slip (which occurred at room temperature); and (3) time-dependent deformation. The inelastic deformation mechanism for the (0)8 MMC was dominated by plasticity at both temperatures. For the (90)8 and (+ or -45)2s MMCs, a combination of damage and plasticity contributed to the deformation at both temperatures.

Majumdar, B. S.↗

Size Segregation and Number Density Enhancement of Particles in Accretion Disk Eddies

We investigate the conditions for trapping solid dust particles in eddies and discuss the behavior of particles in a non-laminar protoplanetary accretion disk. We considered particle sizes from small dust grains to larger objects, 10(exp -4) cm less than a(sub p) less than 10(exp 2) cm. Independent of the source of turbulence, one can expect eddies to exist in the gas flow of a accretion disk, in the form of randomly occurring turbulent features or as convective cells. Due to the centrifugal force, solid particles are driven out of an eddy. It will be shown that this process is inhibited by the gravitational force induced by the protostar. Because of the mass dependence of the friction time, a given eddy becomes a trap for particles of a characteristic size and causes a local change in the dust density. Thus, the size distribution of the grains is no longer spatially homogeneous on small scales. Our general estimates do not depend on special turbulence or convection models. We calculate the maximal inhomogeneity due to this process. The strongest effect was observed for mm-sized particles, which can be concentrated by a factor of 100 within only 100 years.

Klahr, H. H.↗

Star Formation in the Galaxy and the Fluctuating UV Radiation Field

We examine the formation of massive stars in the Galaxy, the resultant fluctuating UV radiation field, and the effect of this Field on the star-forming interstellar medium. Following previous researchers such as Habing (1968), we calculate the average interstellar radiation field at the Solar Circle of the Galaxy. However, our new calculations follow more closely the time dependence of the field at any point. We show that there is a significant difference between the mean field and the median field, and that there are substantial fluctuations of the field (on timescales of order 100 million years) at a given point. Far Ultraviolet Radiation (FUV, photon energies of 6 eV - 13.6 eV) has been recognized as the main source of heating of the neutral interstellar gas. Given the pressure of the interstellar medium (ISM) the FUV field determines whether the thermal balance of the neutral gas results in cold (T approximately 50 - 100 K) clouds (CNM), warm (T about 10,000 K) (WNM), for a combination of the two (the two phase ISM) We present results for the time history of the FUV field for points in the local ISM of the Milky Way Galaxy. The presence of this fluctuating heating rate converts CNM to WNM and vice versa. We show how to calculate the average fractions of the gas in the CNM and WNM when the interstellar gas is subject to this fluctuating FUV field. The knowledge of how these fractions depend on the gas properties (i.e. mean density and composition) and on the FUV-sources (i.e. the star formation rate, or the IMF, or the size distribution of associations) is a basic step in building any detailed model of the large scale behavior of the ISM and the mutual relation between the ISM and the SFR.

Hollenbach, David↗

Melt-Crucible Wetting Behavior in Semiconductor Melt Growth Systems

The wetting angles of several semiconductor-substrate combinations that are of practical importance for crystal growth have been measured: Ga-GaSb-Sb on fused quartz; Ge on fused quartz and carbon-based substrates, each with different roughness; Si on fused quartz plates and on plates coated with fused quartz, Si3N4, and BN powders. The Ga-GaSb-Sb system showed no significant dependence of the wetting angle on the composition despite a large composition dependence of the surface tension. For Ge, the effect of the roughness on the angle could initially be seen on both types of substrates, but for the fused quartz substrates an equilibrium value independent of the surface treatment was reached after several hours of contact time. For Si, total wetting was found for Si3N4 powders. A reduction of the angle over time was found for both fused quartz and BN powders, with the BN powder showing the highest angle at 150-120 deg.

Croell, A.↗

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.↗

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↗