SEARCH · Engineering Papers

Results for “Diffusion rate”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Numerical simulations in combustion

This paper reviews numerical simulations in reacting flows in general and combustion phenomena in particular. It is shown that use of implicit schemes and/or adaptive mesh strategies can improve convergence, stability, and accuracy of the solution. Difficulties increase as turbulence and multidimensions are considered, particularly when finite-rate chemistry governs the given combustion problem. Particular attention is given to the areas of solid-propellant combustion dynamics, turbulent diffusion flames, and spray droplet vaporization.

Chung, T. J.↗

NASA NTRS1989

Direct numerical simulations of transition in a compressible wake

Direct numerical simulations of a compressible wake are performed to determine the sensitivity of large-scale motions to initial phase lags between the most unstable two-dimensional wave and a pair of oblique waves, and to determine the effect of Mach number on the development of three-dimensional motions. The vorticity dynamics of the flow are studied, and it is determined that, depending upon the initial phase lag, the interaction of the streamwise vorticity with the spanwise vortices may result in the formation of closed vortex loops, or in the removal of streamwise vorticity from the saddle region between the spanwise vortices. Due to lower growth rates experienced at high Mach numbers, the roll-up process is delayed, and the resultant nonlinear motions are much more diffuse.

Chen, Jacqueline H.↗

NASA NTRS1992

Upper bounds for convergence rates of vector extrapolation methods on linear systems with initial iterations

The application of the minimal polynomial extrapolation (MPE) and the reduced rank extrapolation (RRE) to a vector sequence obtained by the linear iterative technique x(sub j) + 1 = Ax(sub j) = b,j = 1,2,..., is considered. Both methods produce a two dimensional array of approximations s(sub n,k) to the solution of the system (I - A)x = b. Here, s(sub n,k) is obtained from the vectors x(sub j), n is less than or equal to j is less than or equal to n + k + 1. It was observed in an earlier publication by the first author that the sequence s(sub n,k), k = 1,2,..., for n greater than 0, but fixed, possesses better convergence properties than the sequence s(sub 0,k), k = 1,2,.... A detailed theoretical explanation for this phenomenon is provided in the present work. This explanation is heavily based on approximations by incomplete polynomials. It is demonstrated by numerical examples when the matrix A is sparse that cycling with s(sub n,k) for n greater than 0, but fixed, produces better convergence rates and costs less computationally than cycling with s(sub 0,k). It is also illustrated numerically with a convection-diffusion problem that the former may produce excellent results where the latter may fail completely. As has been shown in an earlier publication, the results produced by s(sub 0,k) are identical to the corresponding results obtained by applying the Arnoldi method or generalized minimal residual scheme (GMRES) to the system (I - A)x = b.

Sidi, Avram↗

NASA NTRS1999

MRI Measurements and Granular Dynamics Simulation of Segregation of Granular Mixture

A counter intuitive axial segregation phenomenon in a rotating horizontal cylinder has recently captured attention of many researchers in different disciplines. There is a growing consensus that the interplay between the particle dynamics and the evolution of the internal structure during the segregation process must be carefully investigated. Magnetic resonance imaging (MRI) has been used to non-invasively obtain much needed dynamic/static information such as velocity and concentration profiles, and it has proven to be capable of depicting the evolution of segregation processes. Segregation in a rotating cylinder involves two processes: the first is to transport small particles in the radial direction to form a radial core, and the second is to transform the radial core into axially segregated bands. Percolation and/or "stopping" have been proposed as mechanisms for the radial segregation. As to mechanisms for axial band formation, much less is known. The difference in the dynamic angle of repose has been proposed to segregate different components in the axial direction. Recently, Hill and Kakalios have reported that particles mix or demix depending upon the competition between diffusion and preferential drift whose order can be determined by the dynamic angle of repose through the adjustment of the rotation rate. We claim that the dynamic angle of repose could be one of the causes, however, it fails to offer reasonable explanations for certain aspects of the axial migration. For example, we always observe that the radial segregation precedes the axial segregation and small particles migrate in the radial direction to form an axially extended radial core. It then transforms into axially segregated bands. By definition, the effects of the dynamic angle of repose are restricted near the free surface where the flowing layer is present. However, during the process of transforming from the radially segregated core to axially segregated bands, small particles located in the deep core region, which is untouched by the flowing layer, also completely disappear. Usually, the dynamics angle of repose are uniquely defined for individual species to characterize particle properties, and the dynamic angle of repose thus defined provides little information for the dynamic angle of repose of the mixture since the concentration ratio and the internal packing structure do not remain the same during the segregation processes. Under microgravity environment, the dynamics angle of repose argument does not hold since there is simply no flowing layer to influence/determine the preferred directions of segregation. We have thus designed an experiment so that the effects of the dynamic angle of repose can be minimized by filling the cylinder almost completely full. Small particles still formed a radial core and also migrated to form axial bands. As ground based experiments we have designed and conducted both 2D and 3D segregation experiments. The 2D experiments are performed using a thin cylinder (the gap between two end caps is about 5 mm) filled with different combinations of particles. The 3D experiments are conducted with a long cylinder of its length and diameter of 27cm and 7cm, respectively. Results of 2D experiments indicate that different mechanisms govern particle motion in regions near and far from the axis of rotation. Results of 3D experiments indicate that a series of collapses of microstructures of particle packing (micro-collapses) may be responsible for the creation of voids for small particles to migrate through in the axial direction. We have successfully eliminated the dynamic angle of repose as a cause for segregation, however, by almost completely filling the cylinder with the particles, we have lost an opportunity to investigate a possibility of particle "mobility" being a cause for segregation which requires a flowing surface but not the difference in the angle of repose. This is currently being investigated.

Nakagawa, M.↗

NASA NTRS2011

Computing Temperatures in Optically Thick Protoplanetary Disks

We worked with a Monte Carlo radiative transfer code to simulate the transfer of energy through protoplanetary disks, where planet formation occurs. The code tracks photons from the star into the disk, through scattering, absorption and re-emission, until they escape to infinity. High optical depths in the disk interior dominate the computation time because it takes the photon packet many interactions to get out of the region. High optical depths also receive few photons and therefore do not have well-estimated temperatures. We applied a modified random walk (MRW) approximation for treating high optical depths and to speed up the Monte Carlo calculations. The MRW is implemented by calculating the average number of interactions the photon packet will undergo in diffusing within a single cell of the spatial grid and then updating the packet position, packet frequencies, and local radiation absorption rate appropriately. The MRW approximation was then tested for accuracy and speed compared to the original code. We determined that MRW provides accurate answers to Monte Carlo Radiative transfer simulations. The speed gained from using MRW is shown to be proportional to the disk mass.

radiative transfer↗

NASA NTRS2016

Comparison of Turbulent Thermal Diffusivity and Scalar Variance Models

In this study, several variable turbulent Prandtl number formulations are examined for boundary layers, pipe flow, and axisymmetric jets. The model formulations include simple algebraic relations between the thermal diffusivity and turbulent viscosity as well as more complex models that solve transport equations for the thermal variance and its dissipation rate. Results are compared with available data for wall heat transfer and profile measurements of mean temperature, the root-mean-square (RMS) fluctuating temperature, turbulent heat flux and turbulent Prandtl number. For wall-bounded problems, the algebraic models are found to best predict the rise in turbulent Prandtl number near the wall as well as the log-layer temperature profile, while the thermal variance models provide a good representation of the RMS temperature fluctuations. In jet flows, the algebraic models provide no benefit over a constant turbulent Prandtl number approach. Application of the thermal variance models finds that some significantly overpredict the temperature variance in the plume and most underpredict the thermal growth rate of the jet. The models yield very similar fluctuating temperature intensities in jets from straight pipes and smooth contraction nozzles, in contrast to data that indicate the latter should have noticeably higher values. For the particular low subsonic heated jet cases examined, changes in the turbulent Prandtl number had no effect on the centerline velocity decay.

computational fluid dynamics↗

NASA NTRS2018

Thermodynamic and Dynamic Aspects of Ice Nucleation

It is known that ice nucleating particles (INP) immersed within supercooled droplets promote the formation of ice. Common theoretical models used to represent this process assume that the immersed particle lowers the work of ice nucleation without significantly affecting the dynamics of water in the vicinity of the particle. This is contrary to evidence showing that immersed surfaces significantly affect the viscosity and diffusivity of vicinal water. To study how this may affect ice formation this work introduces a model linking the ice nucleation rate to the modification of the dynamics and thermodynamics of vicinal water by immersed particles. It is shown that INP that significantly reduce the work of ice nucleation also pose strong limitations to the growth of the nascent ice germs. This leads to the onset of a new ice nucleation regime, called spinodal ice nucleation, where the dynamics of ice germ growth instead of the ice germ size determines the nucleation rate. Nucleation in this regime is characterized by an enhanced sensitivity to particle area and cooling rate. Comparison of the predicted ice nucleation rate against experimental measurements for a diverse set of species relevant to cloud formation suggests that spinodal ice nucleation may be common in nature.

Barahona, Donifan↗

NASA NTRS2021

Hydroxylation of Apollo 17 Soil Sample 78421 by Solar Wind Protons

Hydroxylation by solar wind protons has been simulated in our laboratory on Apollo 17 lunar sample 78421, a very mature regolith sample that is rich with agglutinates (68%). The goal of this study was to determine the rate of hydroxyl formation and their thermal stability by monitoring changes in the SiOH (hydroxyl) stretching band near 3 μm using diffuse reflectance FTIR spectroscopy (DRIFTS). A 2 keV H2+ ion beam was used to simulate proton implantation on 78421 and on a crushed fused silica sample. We find that the OH band does not change unless the samples have been annealed in vacuum prior to irradiation. Qualitatively, the OH bands for the fused silica and 78421 are very different. The OH band for fused silica is centered at 2.74 μm and is relatively sharp ranging from 2.67 - 3.1 μm at full-width-at-half-maximum (FWHM), while the OH band for 78421 is centered at 3.0 μm and ranges from 2.74 - 3.37 μm at FWHM. The increase in wavelength and broadened nature of the OH band in 78421 may be associated with the OH’s proximity to surface defects and/or lattice vacancies. The lack of the H2O bending mode at 6.1 μm indicates that any adsorbed terrestrial H2O is below our detection limit, and therefore the H2O stretching mode at 2.9 μm is not significantly contributing to the broad 3 μm OH band and implies that proton implantation by itself does not lead to water formation. To simulate the maximum dayside temperature on the lunar surface, the lunar sample was heated after proton irradiation. The proton induced OH concentration was reduced by as much as 25% after heating to 400 K (127 °C).

Apollo 17↗

NASA NTRSDate not supplied

The Effect of Reduced Pressure on the Characteristics of Spreading Flames

Flame spread over solid fuels is a canonical problem in fire science, due to its direct implications on material flammability and importance in fire development. In a microgravity environment, such as onboard a spacecraft, flames can behave very differently than on Earth. This is concerning for spaceflight life safety, especially in higher-oxygen environments. Due to the difficulties associated with microgravity testing, low-pressure environments have been proposed as an alternative to approximately replicate the burning behavior of solid fuels observed in reduced gravity conditions because of similar diffusion and heat transfer mechanisms. However, the roles played by gravity and pressure on flame length, standoff distance, and flame spread rate vary with the burning configuration. In concurrent flame spread, the buoyant flow enhances the spread rate by bringing the flame closer to the fuel surface and increasing the heating of the solid fuel. In opposed flame spread, the sample is preheated by the flame ahead of the flame leading edge, which is strongly affected by the surrounding flow field. In this work, we consider flames spreading over thin cotton samples in both downward (opposed) and upward (concurrent) configurations to investigate the effect of pressure (30-100 kPa) on flame characteristics, such as spread rate and standoff distance. A small forced flow is induced upward so that the flames are exposed to a mixed (forced and free) flow. By reducing pressure, flames become less bright, their standoff distance increases, and their spread rates decrease similar to what is observed in low-gravity environments. These results could in help understanding the differences between flames spreading at low pressure and low gravity environments for these similar, yet very different, spreading configurations. They could also provide more information about potential Earth-based flammability testing of materials for spacecraft applications.

fire safety↗

NASA NTRSDate not supplied

The effect of reduced pressure on the characteristics of spreading flames

fire safety↗

NASA NTRSDate not supplied

Improving Ocean Reanalyses and ENSO Forecasts By Assimilation of Rain-Corrected Satellite Sea Surface Salinity Using the GMOA S2s Forecast System

The ENSO phenomenon has a significant global socio-economic impact and has been the key focus for improving coupled ocean-atmosphere forecasts. Assimilation of satellite altimetry and subsurface temperature and salinity from (mostly) Argo help improve the initialization of the thermocline, while satellite SST aids in constraining surface heat-fluxes, leading to improved coupled system sub-seasonal to seasonal forecasts. However, few studies have focused on improving the near-surface density and mixing through satellite sea surface salinity (SSS) assimilation. The few ocean models that assimilate satellite SSS, bias correct to normalize towards the near-surface Argo data for expediency. This assumption is likely inadequate in rainy regions, where buoyant water forms a fresh surface lens. In previous work, we showed that adjusting SSS to bulk salinity (Sb) using the Rain Impact Model (RIM) of Santos-Garcia et al., 2014 improves the near-surface density and mixed layer depth, leading to deeper thermocline and improved NINO3.4 SST forecasts. We now utilize the Soil Moisture and Ocean Salinity rain-corrected (SMOS_RC) SSS, available in SMOS-CATDS products, to represent Sb more accurately at the first model layer (e.g., 5 m). Rather than a diffusivity model as RIM, SMOS_RC uses a statistical correction dependent on Integrated Multi-satellitE Retrievals for GPM (IMERG) rain rates, established on observed SMOS SSS decreases related to Sb in the presence of rain (Supply et al., 2020). For all experiments, all available along-track absolute dynamic topography and in situ observations are assimilated using the LETKF scheme (Penny et al., 2013). One reanalysis additionally assimilates SMOS SSS data as is, and a separate reanalysis assimilates SMOS_RC. We assess the impact on near-surface and subsurface dynamics by validating against observations and explore how SSS assimilation (SMOS vs SMOS_RC) impacts ENSO forecasts using the NASA GMAO Sub-seasonal to Seasonal coupled forecast system (S2S-v3, Molod et al. 2020). We show that improved estimates of density and near-surface mixing led to more accurate coupled air/sea interaction and better ENSO forecasts. The increased SSS, resulting from the removal of the instantaneous rain effect, modifies the ocean state by enhancing mixing and deepening the thermocline.

Veronica Ruiz Xomchuk↗

NASA NTRS1995

Emission of Whistler-mode waves and diffusion of electrons around interplanetary shocks

Whistler-mode wave emissions are frequently observed at and downstream of interplanetary shocks. Using electron distribution functions measured onboard Ulysses in the energy range 1.6 to 862 eV, we calculate the temperature anisotropy and the wave growth rate of the electromagnetic electron cyclotron instability, Results of the calculations are compared to the whistler wave spectra observed simultaneously. For the studied events there is a good correlation between the wave growth rates and the wave spectra. Particularly, upstream of the shock front where no wave emissions are observed, the anisotropy lies below the wave instability threshold, i.e. the critical anisotropy Ac; on the contrary, downstream of the shock, the anisotropy exceeds Ac in some frequency range. Moreover. the tact that the anisotropy is close to Ac in a large frequency range gives prominence to the effect of velocity space diffusion of the electrons by the waves.

Pierre, F.↗

NASA NTRS1981

An automatic light scattering CCN counter

The counter is a static thermal diffusion chamber which has been modified to include an optical system for the determination of droplet concentration by the measurement of scattered light. The determination of concentration is made by measurement of the peak scattered light signal from the cloud of growing droplets which is a function of both the droplet concentration and chamber supersaturation. Because the formation of the peak is related to the rate of growth of the droplets and sedimentation, both of which are determined by supersaturation, the system calibration can be uniquely determined by comparison with an absolute counter such as a static diffusion chamber with a photographic recording system.

Lala, G. G.↗

NASA NTRS2007

NASA Tech Briefs, November 2007

Topics include: Wireless Measurement of Contact and Motion Between Contact Surfaces; Wireless Measurement of Rotation and Displacement Rate; Portable Microleak-Detection System; Free-to-Roll Testing of Airplane Models in Wind Tunnels; Cryogenic Shrouds for Testing Thermal-Insulation Panels; Optoelectronic System Measures Distances to Multiple Targets; Tachometers Derived From a Brushless DC Motor; Algorithm-Based Fault Tolerance for Numerical Subroutines; Computational Support for Technology- Investment Decisions; DSN Resource Scheduling; Distributed Operations Planning; Phase-Oriented Gear Systems; Freeze Tape Casting of Functionally Graded Porous Ceramics; Electrophoretic Deposition on Porous Non- Conductors; Two Devices for Removing Sludge From Bioreactor Wastewater; Portable Unit for Metabolic Analysis; Flash Diffusivity Technique Applied to Individual Fibers; System for Thermal Imaging of Hot Moving Objects; Large Solar-Rejection Filter; Improved Readout Scheme for SQUID-Based Thermometry; Error Rates and Channel Capacities in Multipulse PPM; Two Mathematical Models of Nonlinear Vibrations; Simpler Adaptive Selection of Golomb Power-of- Two Codes; VCO PLL Frequency Synthesizers for Spacecraft Transponders; Wide Tuning Capability for Spacecraft Transponders; Adaptive Deadband Synchronization for a Spacecraft Formation; Analysis of Performance of Stereoscopic-Vision Software; Estimating the Inertia Matrix of a Spacecraft; Spatial Coverage Planning for Exploration Robots; and Increasing the Life of a Xenon-Ion Spacecraft Thruster.

Source record↗

NASA NTRS2007

Runtime and Pressurization Analyses of Propellant Tanks

Multi-element unstructured CFD has been utilized at NASA SSC to carry out analyses of propellant tank systems in different modes of operation. The three regimes of interest at SSC include (a) tank chill down (b) tank pressurization and (c) runtime propellant draw-down and purge. While tank chill down is an important event that is best addressed with long time-scale heat transfer calculations, CFD can play a critical role in the tank pressurization and runtime modes of operation. In these situations, problems with contamination of the propellant by inclusion of the pressurant gas from the ullage causes a deterioration of the quality of the propellant delivered to the test article. CFD can be used to help quantify the mixing and propellant degradation. During tank pressurization under some circumstances, rapid mixing of relatively warm pressurant gas with cryogenic propellant can lead to rapid densification of the gas and loss of pressure in the tank. This phenomenon can cause serious problems during testing because of the resulting decrease in propellant flow rate. With proper physical models implemented, CFD can model the coupling between the propellant and pressurant including heat transfer and phase change effects and accurately capture the complex physics in the evolving flowfields. This holds the promise of allowing the specification of operational conditions and procedures that could minimize the undesirable mixing and heat transfer inherent in propellant tank operation. It should be noted that traditional CFD modeling is inadequate for such simulations because the fluids in the tank are in a range of different sub-critical and supercritical states and elaborate phase change and mixing rules have to be developed to accurately model the interaction between the ullage gas and the propellant. We show a typical run-time simulation of a spherical propellant tank, containing RP-1 in this case, being pressurized with room-temperature nitrogen at 540 R. Nitrogen, shown in blue on the right-hand side of the figures, enters the tank from the diffuser at the top of the figures and impinges on the RP-1, shown in red, while the propellant is being continuously drained at the rate of 1050 lbs/sec through a pipe at the bottom of the tank. The sequence of frames in Figure 1 shows the resultant velocity fields and mixing between nitrogen and RP-1 in a cross-section of the tank at different times. A vortex is seen to form in the incoming nitrogen stream that tends to entrain propellant, mixing it with the pressurant gas. The RP-1 mass fraction contours in Figure 1 are also indicative of the level of mixing and contamination of the propellant. The simulation is used to track the propagation of the pure propellant front as it is drawn toward the exit with the evolution of the mixing processes in the tank. The CFD simulation modeled a total of 10 seconds of run time. As is seen from Figure 1d, after 5.65 seconds the propellant front is nearing the drain pipe, especially near the center of the tank. Behind this pure propellant front is a mixed fluid of compromised quality that would require the test to end when it reaches the exit pipe. Such unsteady simulations provide an estimate of the time that a high-quality propellant supply to the test article can be guaranteed at the modeled mass flow rate. In the final paper, we will discuss simulations of the LOX and propellant tanks at NASA SSC being pressurized by an inert ullage. Detailed comparisons will be made between the CFD simulations and lower order models as well as with test data. Conditions leading to cryo collapse in the tank will also be identified.

Field, Robert E.↗

NASA NTRS1987

Reynolds-stress and dissipation rate budgets in a turbulent channel flow

The budgets for the Reynolds stresses and for the dissipation rate of the turbulence kinetic energy are computed using direct simulation data of a turbulent channel flow. The budget data reveal that all the terms in the budget become important close to the wall. For inhomogeneous pressure boundary conditions, the pressure-strain term is split into a return term, a rapid term, and a Stokes term. The Stokes term is important close to the wall. The rapid and return terms play different roles depending on the component of the term. A split of the velocity pressure-gradient term into a redistributive term and a diffusion term is proposed, which should be simpler to model. The budget data is used to test existing closure models for the pressure-strain term, the dissipation rate, and the transport rate. In general, further work is needed to improve the models.

Mansour, N. N.↗

NASA NTRS1988

Reynolds-stress and dissipation-rate budgets in a turbulent channel flow

Mansour, N. N.↗

NASA NTRS1999

Water Flow Performance of a Superscale Model of the Fastrac Liquid Oxygen Pump

As part of the National Aeronautics and Space Administration's ongoing effort to lower the cost of access to space, the Marshall Space Flight Center has developed a rocket engine with 60,000 pounds of thrust for use on the Reusable Launch Vehicle technology demonstrator slated for launch in 2000. This gas generator cycle engine, known as the Fastrac engine, uses liquid oxygen and RP-1 for propellants and includes single stage liquid oxygen and RP-1 pumps and a single stage supersonic turbine on a common shaft. The turbopump design effort included the first use and application of new suction capability prediction codes and three-dimensional blade generation codes in an attempt to reduce the turbomachinery design and certification costs typically associated with rocket engine development. To verify the pump's predicted cavitation performance, a water flow test of a superscale model of the Fastrac liquid oxygen pump was conducted to experimentally evaluate the liquid oxygen pump's performance at and around the design point. The water flow test article replicated the flow path of the Fastrac liquid oxygen pump in a 1.582x scale model, including scaled seal clearances for correct leakage flow at a model operating speed of 5000 revolutions per minute. Flow entered the 3-blade axial-flow inducer, transitioned to a shrouded, 6-blade radial impeller, and discharged into a vaneless radial diffuser and collection volute. The test article included approximately 50 total and static pressure measurement locations as well as flush-mounted, high frequency pressure transducers for complete mapping of the pressure environment. The primary objectives of the water flow test were to measure the steady-state and dynamic pressure environment of the liquid oxygen pump versus flow coefficient, suction specific speed, and back face leakage flow rate. Results showed excellent correlation between the predicted and experimentally measured pump head rise at low suction specific speeds. Likewise, only small circumferential variations in steady-state impeller exit and radial diffuser pressure distributions were observed from 80% to 120% of the design flow coefficient, matching the computational predictions and confirming that the integrated design approach has minimized any exit volute-induced distortions. The test article exhibited suction performance trends typically observed in inducer designs with virtually constant head rise with decreasing inlet pressure until complete pump head breakdown. Unfortunately, the net positive suction head at 3% head fall-off occurred far below that predicted at all tested flow coefficients, resulting in a negative net positive suction head margin at the design point in water. Additional testing to map the unsteady pressure environment was conducted and interesting rotating phenomena at the inducer inlet were observed. These rotating phenomena's cell numbers, direction, and speed were correlated with pump operating parameters. The impact of the unsteady phenomena and their corresponding energy losses on the unexpectedly poor pump performance is also discussed.

Skelley, Stephen↗