Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flux optimization”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

RNS Applications for Interacting Sub- and Supersonic Flows

A solution based grid adaptation method that combines elements of the multigrid method for solution acceleration and the domain decomposition philosophy for grid optimization is described. Unlike other solution based adaptive gridding schemes, wherein the overhead of recomputing the grid and re-evaluating the solution on the adapted grid leads to higher computational costs compared to a non-adapted calculation, the present methodology reduces the computational time required to obtain the solution. The computational effort involved in the present calculation is significantly lower than a non-adapted calculation that utilizes the multigrid method purely as a convergence acceleration tool. In addition to convergence acceleration, the multigrid framework provides a mechanism of information transfer from regions wherein grid refinement is specified to unrefined coarse grid regions. The basis for domain decomposition in the current procedure is the variation in grid refinement requirements for each coordinate direction in different portions of the flow field. The method is demonstrated herein on an efficient set of governing equations termed the reduced Navier Stokes equations, applied in conjunction with a set of physical boundary conditions. The governing equations are discretized through a pressure based flux splitting procedure that is uniformly applicable from incompressible to supersonic Mach numbers.

Rubin, Stanley G.↗
NASA NTRS1999

High-Area-Ratio Rocket Nozzle at High Combustion Chamber Pressure: Experimental and Analytical Validation

Experimental data were obtained on an optimally contoured nozzle with an area ratio of 1025:1 and on a truncated version of this nozzle with an area ratio of 440:1. The nozzles were tested with gaseous hydrogen and liquid oxygen propellants at combustion chamber pressures of 1800 to 2400 psia and mixture ratios of 3.89 to 6.15. This report compares the experimental performance, heat transfer, and boundary layer total pressure measurements with theoretical predictions of the current Joint Army, Navy, NASA, Air Force (JANNAF) developed methodology. This methodology makes use of the Two-Dimensional Kinetics (TDK) nozzle performance code. Comparisons of the TDK-predicted performance to experimentally attained thrust performance indicated that both the vacuum thrust coefficient and the vacuum specific impulse values were approximately 2.0-percent higher than the turbulent prediction for the 1025:1 configurations, and approximately 0.25-percent higher than the turbulent prediction for the 440:1 configuration. Nozzle wall temperatures were measured on the outside of a thin-walled heat sink nozzle during the test fittings. Nozzle heat fluxes were calculated front the time histories of these temperatures and compared with predictions made with the TDK code. The heat flux values were overpredicted for all cases. The results range from nearly 100 percent at an area ratio of 50 to only approximately 3 percent at an area ratio of 975. Values of the integral of the heat flux as a function of nozzle surface area were also calculated. Comparisons of the experiment with analyses of the heat flux and the heat rate per axial length also show that the experimental values were lower than the predicted value. Three boundary layer rakes mounted on the nozzle exit were used for boundary layer measurements. This arrangement allowed total pressure measurements to be obtained at 14 different distances from the nozzle wall. A comparison of boundary layer total pressure profiles and analytical predictions show good agreement for the first 0.5 in. from the nozzle wall; but the further into the core flow that measurements were taken, the more that TDK overpredicted the boundary layer thickness.

Jankovsky, Robert S.↗
NASA NTRS1993

Growth of zinc selenide single crystals by physical vapor transport in microgravity

The goals of this research were the optimization of growth parameters for large (20 mm diameter and length) zinc selenide single crystals with low structural defect density, and the development of a 3-D numerical model for the transport rates to be expected in physical vapor transport under a given set of thermal and geometrical boundary conditions, in order to provide guidance for an advantageous conduct of the growth experiments. In the crystal growth studies, it was decided to exclusively apply the Effusive Ampoule PVT technique (EAPVT) to the growth of ZnSe. In this technique, the accumulation of transport-limiting gaseous components at the growing crystal is suppressed by continuous effusion to vacuum of part of the vapor contents. This is achieved through calibrated leaks in one of the ground joints of the ampoule. Regarding the PVT transport rates, a 3-D spectral code was modified. After introduction of the proper boundary conditions and subroutines for the composition-dependent transport properties, the code reproduced the experimentally determined transport rates for the two cases with strongest convective flux contributions to within the experimental and numerical error.

Rosenberger, Franz↗
NASA NTRS2001

An Implicit Characteristic Based Method for Electromagnetics

An implicit characteristic-based approach for numerical solution of Maxwell's time-dependent curl equations in flux conservative form is introduced. This method combines a characteristic based finite difference spatial approximation with an implicit lower-upper approximate factorization (LU/AF) time integration scheme. This approach is advantageous for three-dimensional applications because the characteristic differencing enables a two-factor approximate factorization that retains its unconditional stability in three space dimensions, and it does not require solution of tridiagonal systems. Results are given both for a Fourier analysis of stability, damping and dispersion properties, and for one-dimensional model problems involving propagation and scattering for free space and dielectric materials using both uniform and nonuniform grids. The explicit Finite Difference Time Domain Method (FDTD) algorithm is used as a convenient reference algorithm for comparison. The one-dimensional results indicate that for low frequency problems on a highly resolved uniform or nonuniform grid, this LU/AF algorithm can produce accurate solutions at Courant numbers significantly greater than one, with a corresponding improvement in efficiency for simulating a given period of time. This approach appears promising for development of dispersion optimized LU/AF schemes for three dimensional applications.

Beggs, John H.↗
NASA NTRSDate not supplied

Computing the effective elasticity of anisotropic porous media from X-ray computed micro-tomography images

Development and optimization of composite materials designed for thermal protection of NASA’s spacecraft requires the understanding of their physical response to high-enthalpy environments. To predict their macro-scale properties and behavior, high-fidelity 3D simulations are performed at the micro-scale on realistic representations of these composites. The digital micro-structures are generated either synthetically or through X-ray micro-tomography reconstructions. One of the main challenges in the prediction of heatshield material structural response is the computation of the effective elasticity of the fibrous composite, as well as the understanding of the deformation and stresses generated at the micro-scale. These are driven by the fiber layout within the micro-structure and the distribution of the infused matrix. In this effort, the micro-mechanical linear elastic behavior of fibrous ablators is modeled through the use of a numerical method based on the Multi-Point Stress Approximation (MPSA) finite volume scheme. The MPSA, a generalization of the more commonly used Multi-Point Flux Approximation (MPFA), was discussed in a previous presentation at the 15th USCCN and in reference. To predict the behavior of fibrous and woven architectures, algorithms that compute the local fiber orientation are used. The implementation of the MPSA was verified using analytical solutions, engineering test cases and compared against legacy Finite Element Analysis (FEA) software. The stress analysis models were then applied to real geometries used by NASA in thermal protection systems such as fibrous preforms and woven materials and the results were compared to experimental data.

Federico Semeraro↗
NASA NTRS2006

Effects of Vegetation and of Heat and Vapor Fluxes from Soil on Snowpack Evolution and Radiobrightness

The radiobrightness of a snowpack is strongly linked to the snow moisture content profile, to the point that the only operational inversion algorithms require dry snow. Forward dynamic models do not include the effects of freezing and thawing of the soil beneath the snowpack and the effect of vegetation within the snow or above the snow. To get a more realistic description of the evolution of the snowpack, we reported an addition to the Snow-Soil-Vegetation-Atmosphere- Transfer (SSVAT) model, wherein we coupled soil processes of the Land Surface Process (LSP) model with the snow model SNTHERM. In the near future we will be adding a radiobrightness prediction based on the modeled moisture, temperature and snow grain size profiles. The initial investigations with this SSVAT for a late winter and early spring snow pack indicate that soil processes warm the snowpack and the soil. Vapor diffusion needs to be considered whenever the ground is thawed. In the early spring, heat flow from the ground into a snow and a strong temperature gradient across the snow lead to thermal convection. The buried vegetation can be ignored for a late winter snow pack. The warmer surface snow temperature will affect radiobrightness since it is most sensitive to snow surface characteristics. Comparison to data shows that SSVAT provides a more realistic representation of the temperature and moisture profiles in the snowpack and its underlying soil than SNTHERM. The radiobrightness module will be optimized for the prediction of brightness when the snow is moist. The liquid water content of snow causes considerable absorption compared to dry snow, and so longer wavelengths are likely to be most revealing as to the state of a moist snowpack. For volumetric moisture contents below about 7% (the pendular regime), the water forms rings around the contact points between snow grains. Electrostatic modeling of these pendular rings shows that the absorption of these rings is significantly higher than a sphere of the same volume. The first implementation of the radiobrightness module will therefore be a simple radiative transfer model without scattering.

soil↗
NASA NTRS2000

HTMT-class Latency Tolerant Parallel Architecture for Petaflops Scale Computation

Computational Aero Sciences and other numeric intensive computation disciplines demand computing throughputs substantially greater than the Teraflops scale systems only now becoming available. The related fields of fluids, structures, thermal, combustion, and dynamic controls are among the interdisciplinary areas that in combination with sufficient resolution and advanced adaptive techniques may force performance requirements towards Petaflops. This will be especially true for compute intensive models such as Navier-Stokes are or when such system models are only part of a larger design optimization computation involving many design points. Yet recent experience with conventional MPP configurations comprising commodity processing and memory components has shown that larger scale frequently results in higher programming difficulty and lower system efficiency. While important advances in system software and algorithms techniques have had some impact on efficiency and programmability for certain classes of problems, in general it is unlikely that software alone will resolve the challenges to higher scalability. As in the past, future generations of high-end computers may require a combination of hardware architecture and system software advances to enable efficient operation at a Petaflops level. The NASA led HTMT project has engaged the talents of a broad interdisciplinary team to develop a new strategy in high-end system architecture to deliver petaflops scale computing in the 2004/5 timeframe. The Hybrid-Technology, MultiThreaded parallel computer architecture incorporates several advanced technologies in combination with an innovative dynamic adaptive scheduling mechanism to provide unprecedented performance and efficiency within practical constraints of cost, complexity, and power consumption. The emerging superconductor Rapid Single Flux Quantum electronics can operate at 100 GHz (the record is 770 GHz) and one percent of the power required by convention semiconductor logic. Wave Division Multiplexing optical communications can approach a peak per fiber bandwidth of 1 Tbps and the new Data Vortex network topology employing this technology can connect tens of thousands of ports providing a bi-section bandwidth on the order of a Petabyte per second with latencies well below 100 nanoseconds, even under heavy loads. Processor-in-Memory (PIM) technology combines logic and memory on the same chip exposing the internal bandwidth of the memory row buffers at low latency. And holographic storage photorefractive storage technologies provide high-density memory with access a thousand times faster than conventional disk technologies. Together these technologies enable a new class of shared memory system architecture with a peak performance in the range of a Petaflops but size and power requirements comparable to today's largest Teraflops scale systems. To achieve high-sustained performance, HTMT combines an advanced multithreading processor architecture with a memory-driven coarse-grained latency management strategy called "percolation", yielding high efficiency while reducing the much of the parallel programming burden. This paper will present the basic system architecture characteristics made possible through this series of advanced technologies and then give a detailed description of the new percolation approach to runtime latency management.

Sterling, Thomas↗
NASA NTRS1995

Effect of the space environment on materials flown on the EURECA/TICCE-HVI experiment

The primary benefit of accurately quantifying and characterizing the space environmental effects on materials is longer instrument and spacecraft life. Knowledge of the limits of materials allows the designer to optimize the spacecraft design so that the required life is achieved. Materials such as radiator coatings that have excellent durability result in the design of smaller radiators than a radiator coated with a lower durability coating. This may reduce the weight of the spacecraft due to a more optimum design. Another benefit of characterizing materials is the quantification of outgassing properties. Spacecraft which have ultraviolet or visible sensor payloads are susceptible to contamination by outgassed volatile materials. Materials with known outgassing characteristics can be restricted in these spacecraft. Finally, good data on material characteristics improves the ability of analytical models to predict material performance. A flight experiment was conducted on the European Space Agency's European Retrievable Carrier (EuReCa) as part of the Timeband Capture Cell Experiment (TICCE). Our main objective was to gather additional data on the dust and debris environments, with the focus on understanding growth as a function of size (mass) for hypervelocity particles 1E-06 cm and larger. In addition to enumerating particle impacts, hypervelocity particles were to be captured and returned intact. Measurements were performed post-flight to determine the flux density, diameters, and subsequent effects on various optical, thermal control and structural materials. In addition to these principal measurements, the experiment also provided a structure and sample holders for the exposure of passive material samples to the space environment, e.g., the effects of thermal cycling, atomic oxygen, etc. Preliminary results are presented, including the techniques used for intact capture of particles.

Maag, Carl R.↗
NASA NTRS2000

Laser Brazing of High Temperature Braze Alloy

The Space Shuttle Main Engine (SSME) consists of 1080 conical tubes, which are furnace brazed themselves, manifolds, and surrounding structural jacket making almost four miles of braze joints. Subsequent furnace braze cycles are performed due to localized braze voids between the coolant tubes. SSME nozzle experiences extremely high heat flux (180 mW/sq m) during hot fire. Braze voids between coolant tubes may result in hot combustion gas escape causing jacket bulges. The nozzle can be disqualified for flight or result in mission failure if the braze voids exceed the limits. Localized braze processes were considered to eliminate braze voids, however, damage to the parent materials often prohibited use of such process. Being the only manned flight reusable rocket engine, it has stringent requirement on the braze process. Poor braze quality or damage to the parent materials limits the nozzle service life. The objective of this study was to develop a laser brazing process to provide quality, localized braze joints without adverse affect on the parent materials. Gold (Au-Cu-Ni-Pd-Mn) based high temperature braze alloys were used in both powder and wire form. Thin section iron base superalloy A286 tube was used as substrate materials. Different Laser Systems including CO2 (10.6 micrometers, 1kW), ND:YAG (1.06 micrometers, 4kW). and direct diode laser (808nm. 150W) were investigated for brazing process. The laser process variables including wavelength. laser power, travel speed and angle of inclination were optimized according to bead geometry and braze alloy wetting at minimum heat input level, The properties of laser brazing were compared to that of furnace brazing. Microhardness profiles were used for braze joint property comparison between laser and furnace brazing. The cooling rate of laser brazing was compared to furnace brazing based on secondary dendritic arm spacing, Both optical and Scanning Electron Microscope (SEM) were used to evaluate the microstructures of the braze materials and tube substrate. Metallography of the laser braze joint was compared to the furnace braze. SEM Energy Disperse X-Ray Spectra (EDX) and back scattered imaging were used to analyze braze alloy segregation. Although all of the laser systems, CO2, ND:YAG, and direct diode laser produced good braze joint, the direct diode laser was selected for its system simplicity, compactness and portability. Excellent laser and braze alloy coupling is observed with powder alloy compared to braze alloy wire. Good wetting is found with different gold based braze alloys. The laser brazing process can be optimized so that the adverse affect on the parent materials can be eliminated. Metallography of the laser braze joint has shown that quality braze joint was produced with laser brazing process. Penetration of the laser braze to the substrate is at neglectable level. Zero penetration is observed. Microstructure examinations shown that no observable changes of the microstructure (grain structure and precipitation) in the HAZ area between laser braze and furnace braze. Wide gaps can be laser brazed with single pass for up to 0.024 inches. Finer dendritic structure is observed in laser brazing compared with equiaxial and coarser grain of the furnace brazing microstructure. Greater segregation is also found in the furnace braze. Higher hardness of the laser braze joint comparing to furnace braze is observed due to the fast cooling rate and Finer microstructure in the laser brazing. Laser braze joint properties meet or exceed the furnace joint properties. Direct diode laser for thin section tube brazing with high temperature braze alloys have been successfully demonstrated. The laser's high energy density and precise control has shown significant advantages in reducing process heat input to the substrates and provide high quality braze joints comparing to other localized braze process such as torch, TIG, and MPTA processes. Significant cost savings can be realized particularly with localized braze comparing to a full furnace braze cycle.

Gao, Y. P.↗
NASA NTRS2008

Evaluations of Risks from the Lunar and Mars Radiation Environments

Protecting astronauts from the space radiation environments requires accurate projections of radiation in future space missions. Characterization of the ionizing radiation environment is challenging because the interplanetary plasma and radiation fields are modulated by solar disturbances and the radiation doses received by astronauts in interplanetary space are likewise influenced. The galactic cosmic radiation (GCR) flux for the next solar cycle was estimated as a function of interplanetary deceleration potential, which has been derived from GCR flux and Climax neutron monitor rate measurements over the last 4 decades. For the chaotic nature of solar particle event (SPE) occurrence, the mean frequency of SPE at any given proton fluence threshold during a defined mission duration was obtained from a Poisson process model using proton fluence measurements of SPEs during the past 5 solar cycles (19-23). Analytic energy spectra of 34 historically large SPEs were constructed over broad energy ranges extending to GeV. Using an integrated space radiation model (which includes the transport codes HZETRN [1] and BRYNTRN [2], and the quantum nuclear interaction model QMSFRG[3]), the propagation and interaction properties of the energetic nucleons through various media were predicted. Risk assessment from GCR and SPE was evaluated at the specific organs inside a typical spacecraft using CAM [4] model. The representative risk level at each event size and their standard deviation were obtained from the analysis of 34 SPEs. Risks from different event sizes and their frequency of occurrences in a specified mission period were evaluated for the concern of acute health effects especially during extra-vehicular activities (EVA). The results will be useful for the development of an integrated strategy of optimizing radiation protection on the lunar and Mars missions. Keywords: Space Radiation Environments; Galactic Cosmic Radiation; Solar Particle Event; Radiation Risk; Risk Analysis; Radiation Protection.

Kim, Myung-Hee↗
NASA NTRS2012

Can Terrestrial Microbes Grow on Mars?

The theme for AbSciCon 2012 is "Exploring Life: Past and Present, Near and Far." The conference will address our current understanding of life - from processes at the molecular level to those which operate at planetary scales. Studying these aspects of life on Earth provides an essential platform from which to examine the potential for life on other worlds, both within our solar system and beyond. Mars exhibits a variety of extreme environments characterized by high UV and ionizing radiation flux, low pressure anoxic atmosphere, scarce or absent liquid water, extreme low temperatures, etc. The ability of terrestrial microorganisms to survive and adapt to the Mars environment has profound implications for astrobiology, planetary protection, and Mars life detection missions. At the NASA Ames Synthetic Biology Initiative, we believe that synthetic biology has the potential to revolutionize human space exploration. As such, the initiative is dedicated to applying the tools and techniques of synthetic biology to space exploration and astrobiology. Biological solutions will be invaluable for space exploration because they are not resource intensive, and they are versatile and self-renewing. An understanding of how to work with DNA in an unfavorable environment is paramount to utilizing biological tools on space missions. Furthermore, the ability to adjust life to the parameters of Mars is vital both to discovering what life on Mars might look like, and to using biological tools under such conditions. As a first step, we need an energy-efficient, low cost means of transporting, storing, and protecting genomic DNA, DNA parts, and whole microbial strains. Our goal is to develop and demonstrate viable and superior alternatives to standard DNA storage methods, which can be optimized to the conditions of space exploration, using synthetic biology as a tool. This includes protocols and kit designs for easy and repeatable DNA and strain recovery from protective storage conditions. We are constructing newly engineered genetic parts for different valuable host organisms, designed to increased long-term survival and functional retention. These methods should be applied for DNA and strain storage and transportation. In parallel, we seek inspiration from natural organisms that have developed means for survival in extreme environmental conditions. We are utilizing novel techniques for analysis of lipid biomarkers in the Antarctic Dry Valleys in order to identify resident microbes in the Antarctic soil and permafrost, as well as biomarker fossils of organisms that survived in the valleys in ages past. Through the identification of these life forms, we hope to understand and draw on new biological tools and strategies for synthetic biological applications on Mars.

Rothschild, Lynn↗
NASA NTRS2023

Direct Measurement of the Spectral Structure of Cosmic-Ray Electrons+Positrons in the Tev Region With Calet on the International Space Station

Detailed measurements of the spectral structure of cosmic-ray electrons and positrons from 10.6 GeV to 7.5 TeV are presented from over 7 years of observations with the CALorimetric Electron Telescope (CALET) on the International Space Station. The instrument, consisting of a charge detector, an imaging calorimeter, and a total absorption calorimeter with a total depth of 30 radiation lengths at normal incidence and a fine shower imaging capability, is optimized to measure the all-electron spectrum well into the TeV region. Because of the excellent energy resolution (a few percent above 10 GeV) and the outstanding e/p separation (10 5 ), CALET provides optimal performance for a detailed search of structures in the energy spectrum. The analysis uses data up to the end of 2022, and the statistics of observed electron candidates has increased more than 3 times since the last publication in 2018. By adopting an updated boosted decision tree analysis, a sufficient proton rejection power up to 7.5 TeV is achieved, with a residual proton contamination less than 10%. The observed energy spectrum becomes gradually harder in the lower energy region from around 30 GeV, consistently with AMS-02, but from 300 to 600 GeV it is considerably softer than the spectra measured by DAMPE and Fermi-LAT. At high energies, the spectrum presents a sharp break around 1 TeV, with a spectral index change from −3.15 to −3.91, and a broken power law fitting the data in the energy range from 30 GeV to 4.8 TeV better than a single power law with 6.9 sigma significance, which is compatible with the DAMPE results. The break is consistent with the expected effects of radiation loss during the propagation from distant sources (except the highest energy bin). We have fitted the spectrum with a model consistent with the positron flux measured by AMS-02 below 1 TeV and interpreted the electron+positron spectrum with possible contributions from pulsars and nearby sources. Above 4.8 TeV, a possible contribution from known nearby supernova remnants, including Vela, is addressed by an event-by-event analysis providing a higher proton-rejection power than a purely statistical analysis.

O. Adriani↗
NASA NTRS1989

Minimum length Pb/SCIN detector for efficient cosmic ray identification

A study was made of the performance of a minimal length cosmic ray shower detector that would be light enough for space flight and would provide efficient identification of positions and protons. Cosmic ray positions are mainly produced in the decay chain of: Pion yields Muon yields Positron and they provide a measure of the matter density traversed by primary protons. Present positron flux measurements are consistent with the Leaky Box and Halo models for sources of cosmic rays. Abundant protons in the space environment are a significant source of background that would wash out the positron signal. Protons and positrons produced very distictive showers of particles when they enter matter; many studies have been published on their behavior on large calorimeter detectors. The challenge is to determine the minimal material necessary (minimal calorimeter depth) for positive particles identification. The primary instrument for the investigation is the Monte Carlo code GEANT, a library of programs from CERN that can be used to model experimental geometry, detector responses and particle interaction processes. The use of the Monte Carlo approach is crucial since statistical fluctuations in shower shape are significant. Studies conducted during the 1988 summer program showed that straightforward approaches to the problem achieved 85 to 90 percent correct identification, but left a residue of 10 to 15 percent misidentified particles. This percentage improved to a few percent when multiple shower-cut criteria were applied to the data. This summer, the same study was extended to employ several physical and statistical methods of identifying response of the calorimeter and the efficiency of the optimal shower cuts to off-normal incidence particle was determined.

Snyder, H. David↗
NASA NTRS2013

Influence of Computational Drop Representation in LES of a Droplet-Laden Mixing Layer

Multiphase turbulent flows are encountered in many practical applications including turbine engines or natural phenomena involving particle dispersion. Numerical computations of multiphase turbulent flows are important because they provide a cheaper alternative to performing experiments during an engine design process or because they can provide predictions of pollutant dispersion, etc. Two-phase flows contain millions and sometimes billions of particles. For flows with volumetrically dilute particle loading, the most accurate method of numerically simulating the flow is based on direct numerical simulation (DNS) of the governing equations in which all scales of the flow including the small scales that are responsible for the overwhelming amount of dissipation are resolved. DNS, however, requires high computational cost and cannot be used in engineering design applications where iterations among several design conditions are necessary. Because of high computational cost, numerical simulations of such flows cannot track all these drops. The objective of this work is to quantify the influence of the number of computational drops and grid spacing on the accuracy of predicted flow statistics, and to possibly identify the minimum number, or, if not possible, the optimal number of computational drops that provide minimal error in flow prediction. For this purpose, several Large Eddy Simulation (LES) of a mixing layer with evaporating drops have been performed by using coarse, medium, and fine grid spacings and computational drops, rather than physical drops. To define computational drops, an integer NR is introduced that represents the ratio of the number of existing physical drops to the desired number of computational drops; for example, if NR=8, this means that a computational drop represents 8 physical drops in the flow field. The desired number of computational drops is determined by the available computational resources; the larger NR is, the less computationally intensive is the simulation. A set of first order and second order flow statistics, and of drop statistics are extracted from LES predictions and are compared to results obtained by filtering a DNS database. First order statistics such as Favre averaged stream-wise velocity, Favre averaged vapor mass fraction, and the drop stream-wise velocity, are predicted accurately independent of the number of computational drops and grid spacing. Second order flow statistics depend both on the number of computational drops and on grid spacing. The scalar variance and turbulent vapor flux are predicted accurately by the fine mesh LES only when NR is less than 32, and by the coarse mesh LES reasonably accurately for all NR values. This is attributed to the fact that when the grid spacing is coarsened, the number of drops in a computational cell must not be significantly lower than that in the DNS.

Bellan, Josette↗
NASA NTRS1998

Seasonal Cycles of Meridional Overturning and Heat Transport of the Indian Ocean

A general circulation model of the Indian Ocean is fitted to monthly averaged climatological temperatures, salinities, and surface fluxes using the adjoint method. Interannual variability is minimized by penalizing the temporal drift from one seasonal cycle to another during a two-year integration. The resultant meridional overturning and heat transport display large seasonal variations, with maximum amplitudes of 18 and 22 (x 10(exp 6) cubic m/s) for the overturning and 1.8 and 1.4 (x 10(exp 15) W) for heat transport near 10 S and 10 N, respectively. A dynamical decomposition of the overturning and heat transport shows that the time-varying Ekman How plus its barotropic compensation can explain a large part of the seasonal variations in overturning and heat transport. The maximum variations at 10 deg N and 10 deg S are associated with monsoon reversal over the northern Indian Ocean and changes of the easterlies over the southern Indian Ocean. An external mode with variable topography has a moderate contribution where the Somali Current and the corresponding gyre reverse direction seasonally. Contribution front vertical shear (thermal wind and ageostrophic shear) is dominant near the southern boundary and large near the Somali Current latitudes. The dominant balance in the zonally integrated heat budget is between heat storage change and heat transport convergence except south of 15 S. Optimization with seasonal forcings improves estimates of sea surface temperatures, but the annual average overturning and heat transport are very similar to previous results with annual mean forcings. The annual average heat transport consists of roughly equal contributions from time-mean and time-varying fields of meridional velocities and temperatures in the northern Indian Ocean. indicating a significant rectification to the heat transport due to the time-varying fields. The time-mean and time-varying contributions are primarily due to the overturning and horizontal gyre, respectively. Inclusion of TOPEX data enhances the seasonal cycles of the estimated overturning and heat transport in the central Indian Ocean significantly and improves the estimated equatorial zonal flows but leads 10 unrealistic estimates of the velocity structure near the Indonesian Throughflow region, most likely owing to the deficiencies in the lateral boundary conditions.

Lee, Tong↗
NASA NTRS2016

Basin-Scale Assessment of the Land Surface Water Budget in the National Centers for Environmental Prediction Operational and Research NLDAS-2 Systems

The purpose of this study is to evaluate the components of the land surface water budget in the four land surface models (Noah, SAC-Sacramento Soil Moisture Accounting Model, (VIC) Variable Infiltration Capacity Model, and Mosaic) applied in the newly implemented National Centers for Environmental Prediction (NCEP) operational and research versions of the North American Land Data Assimilation System version 2 (NLDAS-2). This work focuses on monthly and annual components of the water budget over 12 National Weather Service (NWS) River Forecast Centers (RFCs). Monthly gridded FLUX Network (FLUXNET) evapotranspiration (ET) from the Max-Planck Institute (MPI) of Germany, U.S. Geological Survey (USGS) total runoff (Q), changes in total water storage (dS/dt, derived as a residual by utilizing MPI ET and USGS Q in the water balance equation), and Gravity Recovery and Climate Experiment (GRACE) observed total water storage anomaly (TWSA) and change (TWSC) are used as reference data sets. Compared to these ET and Q benchmarks, Mosaic and SAC (Noah and VIC) in the operational NLDAS-2 overestimate (underestimate) mean annual reference ET and underestimate (overestimate) mean annual reference Q. The multimodel ensemble mean (MME) is closer to the mean annual reference ET and Q. An anomaly correlation (AC) analysis shows good AC values for simulated monthly mean Q and dS/dt but significantly smaller AC values for simulated ET. Upgraded versions of the models utilized in the research side of NLDAS-2 yield largely improved performance in the simulation of these mean annual and monthly water component diagnostics. These results demonstrate that the three intertwined efforts of improving (1) the scientific understanding of parameterization of land surface processes, (2) the spatial and temporal extent of systematic validation of land surface processes, and (3) the engineering-oriented aspects such as parameter calibration and optimization are key to substantially improving product quality in various land data assimilation systems.

NLDAS-2↗
NASA NTRS2004

Carbon Dioxide Dispersion in the Combustion Integrated Rack Simulated Numerically

When discharged into an International Space Station (ISS) payload rack, a carbon dioxide (CO2) portable fire extinguisher (PFE) must extinguish a fire by decreasing the oxygen in the rack by 50 percent within 60 sec. The length of time needed for this oxygen reduction throughout the rack and the length of time that the CO2 concentration remains high enough to prevent the fire from reigniting is important when determining the effectiveness of the response and postfire procedures. Furthermore, in the absence of gravity, the local flow velocity can make the difference between a fire that spreads rapidly and one that self-extinguishes after ignition. A numerical simulation of the discharge of CO2 from PFE into the Combustion Integrated Rack (CIR) in microgravity was performed to obtain the local velocity and CO2 concentration. The complicated flow field around the PFE nozzle exits was modeled by sources of equivalent mass and momentum flux at a location downstream of the nozzle. The time for the concentration of CO2 to reach a level that would extinguish a fire anywhere in the rack was determined using the Fire Dynamics Simulator (FDS), a computational fluid dynamics code developed by the National Institute of Standards and Technology specifically to evaluate the development of a fire and smoke transport. The simulation shows that CO2, as well as any smoke and combustion gases produced by a fire, would be discharged into the ISS cabin through the resource utility panel at the bottom of the rack. These simulations will be validated by comparing the results with velocity and CO2 concentration measurements obtained during the fire suppression system verification tests conducted on the CIR in March 2003. Once these numerical simulations are validated, portions of the ISS labs and living areas will be modeled to determine the local flow conditions before, during, and after a fire event. These simulations can yield specific information about how long it takes for smoke and combustion gases produced by a fire to reach a detector location, how large the fire would be when the detector alarms, and the behavior of the fire until it has been extinguished. This new capability could then be used to optimize the location of fire detectors and fire-suppression ports as well as to evaluate the effectiveness of fire suppressants and response strategies. Numerical data collected from these simulations could also be used to develop a virtual reality fire event for crew training and fire safety awareness. This work is funded by NASA's Bioastronautics Initiative, which has the objective of ensuring and enhancing the health, safety, and performance of humans in space. As part of this initiative, the Microgravity Combustion Science Branch at the NASA Glenn Research Center is conducting spacecraft fire safety research to significantly improve fire safety on inhabited spacecraft.

Wu, Ming-Shin↗
NASA NTRS1999

Vibration-Induced Droplet Atomization

Thermal management is critical to a number of technologies used in a microgravity environment and in Earth-based systems. Examples include electronic cooling, power generation systems, metal forming and extrusion, and HVAC (heating, venting, and air conditioning) systems. One technique that can deliver the large heat fluxes required for many of these technologies is two-phase heat transfer. This type of heat transfer is seen in the boiling or evaporation of a liquid and in the condensation of a vapor. Such processes provide very large heat fluxes with small temperature differences. Our research program is directed toward the development of a new, two-phase heat transfer cell for use in a microgravity environment. In this paper, we consider the main technology used in this cell, a novel technique for the atomization of a liquid called vibration-induced droplet atomization. In this process, a small liquid droplet is placed on a thin metal diaphragm that is made to vibrate by an attached piezoelectric transducer. The vibration induces capillary waves on the free surface of the droplet that grow in amplitude and then begin to eject small secondary droplets from the wave crests. In some situations, this ejection process develops so rapidly that the entire droplet seems to burst into a small cloud of atomized droplets that move away from the diaphragm at speeds of up to 50 cm/s. By incorporating this process into a heat transfer cell, the active atomization and transport of the small liquid droplets could provide a large heat flux capability for the device. Experimental results are presented that document the behavior of the diaphragm and the droplet during the course of a typical bursting event. In addition, a simple mathematical model is presented that qualitatively reproduces all of the essential features we have seen in a burst event. From these two investigations, we have shown that delayed droplet bursting results when the system passes through a resonance condition. This occurs when the initial acceleration of the diaphragm is higher than the critical acceleration and the driving frequency is larger than the initial resonance frequency of the diaphragm-droplet system. We have incorporated this droplet atomization device into a design for a new heat transfer cell for use in a microgravity environment. The cell is essentially a cylindrical container with a hot surface on one end and a cold surface on the other. The vibrating diaphragm is mounted in the center of the cold surface. Heat transfer occurs through droplet evaporation and condensation on the hot and cold ends of the cell. A prototype of this heat transfer cell has been built and tested. It can operate continuously and provides a modest level of heat transfer, about 20 W/sq cm. Our work during the next few years will be to optimize the design of this cell to see if we can produce a device that has significantly better performance than conventional heat exchangers and heat pipes.

Smith, M. K.↗