SEARCH · Engineering Papers
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.
A 2D numerical model for a turbulent surface jet issuing into a finite-depth channel is presented. The model uses an unstructured, finite element procedure. The calculations confirm the experimental observation that scale similarity may be apparent in the mean velocity profile even in the presence of a large encroaching recirculation cell.
An analysis of linear wave modes associated with supersonic jets confined inside an acoustically lined rectangular duct is presented. Mathematical formulations are given for the vortex-sheet model and continuous mean flow model of the jet flow profiles. Detailed dispersion relations of these waves in a two-dimensional confined jet as well as an unconfined free jet are computed. Effects of the confining duct and the liners on the jet instability and acoustic waves are studied numerically. It is found that the effect of the liners is to attenuate waves that have supersonic phase velocities relative to the ambient flow. Numerical results also show that the growth rates of the instability waves could be reduced significantly by the use of liners. In addition, it is found that the upstream propagating neutral waves of an unconfined jet could become attenuated when the jet is confined.
The production rates of spallogenic radionuclides and stable isotopes in the L-chondrite Knyahinya were investigated using the MCNPX code. Numerous cosmogenic nuclides had been measured in many Knyahinya samples. The pre-atmospheric size and sample locations of Knyahinya are well known, thus Knyahinya is a good test case for cosmogenic-nuclide production-rate calculations. Our calculated profiles were compared to the measurements to determine effective proton fluxes.
This report presents the use of a time-marching three-dimensional compressible Navier-Stokes equation numerical solver with a one-equation turbulence model to simulate the flow fields developed along concave wall surfaces without and with a downstream extension flat wall surface. The 3-D Navier- Stokes numerical solver came from the NASA Glenn-HT code. The one-equation turbulence model was derived from the Spalart and Allmaras model. The computational approach was first calibrated with the computations of the velocity and Reynolds shear stress profiles of a steady flat plate boundary layer flow. The computational approach was then used to simulate developing boundary layer flows along concave wall surfaces without and with a downstream extension wall. The author investigated the computational results of surface friction factors, near surface velocity components, near wall temperatures, and a turbulent shear stress component in terms of turbulence modeling, computational mesh configurations, inlet turbulence level, and time iteration step. The computational results were compared with existing measurements of skin friction factors, velocity components, and shear stresses of the developing boundary layer flows. With a fine computational mesh and a one-equation model, the computational approach could predict accurately the skin friction factors, near surface velocity and temperature, and shear stress within the flows. The computed velocity components and shear stresses also showed the vortices effect on the velocity variations over a concave wall. The computed eddy viscosities at the near wall locations were also compared with the results from a two equation turbulence modeling technique. The inlet turbulence length scale was found to have little effect on the eddy viscosities at locations near the concave wall surface. The eddy viscosities, from the one-equation and two-equation modeling, were comparable at most stream-wise stations. The present one-equation turbulence model is an effective approach for turbulence modeling in the near solid wall surface region of flow over a concave wall.
The vast illuminating power of the Global Positioning System (GPS), which transformed space geodesy in the 199Os, is now serving to probe the earth's fluid envelope in unique ways. Three distinct techniques have emerged: ground-based sensing of the integrated atmospheric moisture; space-based profiling of atmospheric refractivity, pressure, temperature, moisture, and other properties by active limb sounding; and surface (ocean and ice) altimetry and scatterometry with reflected signals detected from space. Ground-based GPS moisture sensing is already in provisional use for numerical weather prediction. Limb sounding, while less mature, offers a bevy of attractions, including high accuracy, stability, and vertical resolution; all-weather operation; and exceptionally low cost. GPS bistatic radar, r 'reflectometry,' is the least advanced but shows promise for a number of niche applications.
We investigated the anthropogenic and volcanic contributions to sulfate aerosol in the stratosphere through modeling and analysis of satellite data. We use a global model GOCART to simulate SO2 and sulfate aerosol in the period of 2000 to 2010, during which numerous volcanic eruptions occurred although nothing like the magnitudes of El Chichon or Pinatubo. We compared the model results with the column SO2 data from OMI and stratospheric SO2 data from MLS instrument on Aura satellite and the aerosol vertical profiles from the SCIAMACHY instrument on Envisat and the CALIOP instrument on CALIPSO satellites. Finally, we assessed the relative contributions of volcanic aerosols vs. anthropogenic aerosols to the observed decadal stratospheric aerosol trends.
Multi-year sets of ground-based sun-photometer measurements conducted at 12 Arctic sites and 9 Antarctic sites were examined to determine daily mean values of aerosol optical thickness tau(lambda) at visible and near-infrared wavelengths, from which best-fit values of Ångström's exponent alpha were calculated. Analyzing these data, the monthly mean values of tau(0.50 micrometers) and alpha and the relative frequency histograms of the daily mean values of both parameters were determined for winter-spring and summer-autumn in the Arctic and for austral summer in Antarctica. The Arctic and Antarctic covariance plots of the seasonal median values of alpha versus tau(0.50 micrometers) showed: (i) a considerable increase in tau(0.50 micrometers) for the Arctic aerosol from summer to winter-spring, without marked changes in alpha; and (ii) a marked increase in tau(0.50 micrometer) passing from the Antarctic Plateau to coastal sites, whereas alpha decreased considerably due to the larger fraction of sea-salt aerosol. Good agreement was found when comparing ground-based sun-photometer measurements of tau(lambda) and alpha at Arctic and Antarctic coastal sites with Microtops measurements conducted during numerous AERONET/MAN cruises from 2006 to 2013 in three Arctic Ocean sectors and in coastal and off-shore regions of the Southern Atlantic, Pacific, and Indian Oceans, and the Antarctic Peninsula. Lidar measurements were also examined to characterize vertical profiles of the aerosol backscattering coefficient measured throughout the year at Ny-Ålesund. Satellite-based MODIS, MISR, and AATSR retrievals of tau(lambda) over large parts of the oceanic polar regions during spring and summer were in close agreement with ship-borne and coastal ground-based sun-photometer measurements. An overview of the chemical composition of mode particles is also presented, based on in-situ measurements at Arctic and Antarctic sites. Fourteen log-normal aerosol number size-distributions were defined to represent the average features of nuclei, accumulation and coarse mode particles for Arctic haze, summer background aerosol, Asian dust and boreal forest fire smoke, and for various background austral summer aerosol types at coastal and high-altitude Antarctic sites. The main columnar aerosol optical characteristics were determined for all 14 particle modes, based on in-situ measurements of the scattering and absorption coefficients. Diurnally averaged direct aerosol-induced radiative forcing and efficiency were calculated for a set of multimodal aerosol extinction models, using various Bidirectional Reflectance Distribution Function models over vegetation-covered, oceanic and snow-covered surfaces. These gave a reliable measure of the pronounced effects of aerosols on the radiation balance of the surface-atmosphere system over polar regions.
Preliminary evaluation of the segmented ion thruster (SIT) ion beam characteristics was performed by measuring the two dimensional ion current density distribution at four axial locations downstream from the thruster. Current density profiles were measured with the concurrent operation of multiple ion sources. Preliminary assessments regarding symmetry of the plume, characterization of plume development as a function of axial position, and the interactions of multiple simultaneously operating ion sources were conducted. A numerical integration scheme was applied to the acquired data to determine the total beam current. A qualitative analysis of the downstream development of the ion beam was conducted. Most of the data acquired so far has been on somewhat of a trial basis and extensive testing is scheduled to resume in August 1994.
The paper presents a three-dimensional CFD simulation for predicting the liquid nitrogen chilldown process through a SS-304 tube. A dispersed mixture model has been employed to predict the two-phase flow behavior whereas the Lee phase change model has been used to predict the boiling process. Cryogenic chilldown process starts with a pure vapor region and proceeds to film, transition and nucleate boiling regimes before reaching the single phase liquid flow and complete chilldown of the transfer line. The present numerical model predicted the film boiling regime accurately and the slope of the chilldown curve deviates slightly from the experimental curve towards the transition and nucleate boiling regime, which needs further investigation. The contours of liquid volume fraction and temperature profiles has been presented to understand the inverted annular film flows in the film boiling regime in detail. In future, the present CFD model will be extended to a wide range of chilldown conditions and different cryogenic fluids to develop a robust methodology.
Advances in micro-fabrication processes have generated tremendous interests in miniaturizing chemical and biomedical analyses into integrated microsystems (Lab-on-Chip devices). To successfully design and operate the micro fluidics system, it is essential to understand the fundamental fluid flow phenomena when channel sizes are shrink to micron or even nano dimensions. One important phenomenon is the electro kinetic effect in micro/nano channels due to the existence of the electrical double layer (EDL) near a solid-liquid interface. Not only EDL is responsible for electro-osmosis pumping when an electric field parallel to the surface is imposed, EDL also causes extra flow resistance (the electro-viscous effect) and flow anomaly (such as early transition from laminar to turbulent flow) observed in pressure-driven microchannel flows. Modeling and simulation of electro-kinetic effects on micro flows poses significant numerical challenge due to the fact that the sizes of the double layer (10 nm up to microns) are very thin compared to channel width (can be up to 100 s of m). Since the typical thickness of the double layer is extremely small compared to the channel width, it would be computationally very costly to capture the velocity profile inside the double layer by placing sufficient number of grid cells in the layer to resolve the velocity changes, especially in complex, 3-d geometries. Existing approaches using "slip" wall velocity and augmented double layer are difficult to use when the flow geometry is complicated, e.g. flow in a T-junction, X-junction, etc. In order to overcome the difficulties arising from those two approaches, we have developed a sub-grid integration method to properly account for the physics of the double layer. The integration approach can be used on simple or complicated flow geometries. Resolution of the double layer is not needed in this approach, and the effects of the double layer can be accounted for at the same time. With this approach, the numeric grid size can be much larger than the thickness of double layer. Presented in this report are a description of the approach, methodology for implementation and several validation simulations for micro flows.
The problem of microwave thermal emission from a half-space random medium is solved. We consider a laminar structure which has a nonuniform temperature distribution in the vertical direction. A radiative transfer approach is applied. For constant absorption and scattering coefficients, the brightness temperature is determined by a simple closed-form formula. Physical interpretations and numerical results are illustrated and discussed for the various cases.
Natural convection heat transfer in a two-dimensional differentially heated semiannular enclosure is studied. The enclosure is isothermally heated and cooled at the inner and outer walls, respectively. A commercial software based on the SIMPLER algorithm was used to simulate the velocity and temperature profiles. Various parameters that affect the momentum and heat transfer processes were examined. These parameters include the Rayleigh number, Prandtl number, radius ratio, and the angle of inclination. A flow regime extending from conduction-dominated to convection-dominated flow was examined. The computed results of heat transfer are presented as a function of flow parameter and geometric factors. It is found that the heat transfer rate attains a minimum when the enclosure is tilted about +50 deg with respect to the gravitational direction.
Global maps of the monthly mean net upward longwave radiation flux at the ocean surface were obtained for April, July, October 1985 and January 1986. These maps were produced by blending information obtained from a combination of general circulation model cloud radiative forcing fields, the top of the atmosphere cloud radiative forcing from ERBE and TOVS profiles and sea surface temperature on ISCCP C1 tapes. The fields are compatible with known meteorological regimes of atmospheric water vapor content and cloudiness. There is a vast area of high net upward longwave radiation flux (greater than 80/sq Wm) in the eastern Pacific Ocean throughout most of the year. Areas of low net upward longwave radiation flux ((less than 40/sq Wm) are the tropical convective regions and extra tropical regions that tend to have persistent low cloud cover.The technique used relies on General Circulation Model simulations and so is subject to some of the uncertainties associated with the model. However, all input information regarding temperature, moisture, and cloud cover is from satellite data having near global coverage. This feature of the procedure alone warrants its consideration for further use in compiling global maps of longwave radiation.
An iteration difference method for the calculation of the incompressible laminar boundary layer is described. The method uses Prandtl's boundary layer equation and the boundary conditions directly and permits the attainment of an arbitrary accuracy. The method has been tested successfully in the continuation of the Blasius profile on the flat plate, on the circular cylinder investigated by Heimenz and on an elliptical cylinder of fineness ratio 1:4. The method makes possible the testing of previously developed methods, all of which contain important assumptions.
Abstract Drying of moist porous media can be very energy inefficient. For example, in the pulp and paper industry, paper drying consumes more than two-thirds of the total energy used in paper machines. Novel drying technologies can decrease the energy used for drying and lessen the manufacturing processes' carbon footprint. Developing next-generation drying technologies to dry moist porous media may require an understanding of removing moisture from a fully saturated porous material with excess water. This paper provides a fundamental understanding of heat and mass transfer in a fully saturated porous medium with excess water. This is relevant, for example, in drying tissue as well as pulp or paper for the purpose of thermal insulation where pressing is preferred to be avoided to overcome the reduction in the sheet thickness. For this purpose, a theoretical drying model is developed where the porous medium corresponds to paper and is assumed to be sandwiched between two excess-water layers (bottom and top). The conjugate model consists of energy and mass conservation equations for each layer. The model is validated with corresponding experimental data. In the model, the thickness of each water layer is calculated as a function of drying time based on local temperature and total moisture content. The numerical model is transient and one-dimensional in space (i.e., in the thickness direction). This paper demonstrates the governing equations, boundary conditions, and results when the saturated porous medium with water layers is heated from one side. Moisture and temperature profiles are estimated in the thickness direction of the porous medium as it dries.
Upon thermal aging, nitroplasticizer (NP) can thermally degrade and release numerous degradation products depending on aging conditions. In this report, we summarize our recent findings, a library of key degradation products in the NP thermal aging is provided. The library contains 25 identified species with their corresponding MS/MS spectrum and intensity profiles, which are derived from our 44-month thermal aging experiment of NP conducted in air, nitrogen, and water at room temperature (RT), 38, 45, 55, and 64°C. To provide a comprehensive yet simplified overview of NP degradation, we summarized the proposed mechanisms of HONO elimination and decomposition, acid generation, PBNA nitrosation and nitration, acid-catalyzed hydrolyses of NP and polynitrated PBNAs, as well as the formation of NP degradants. As the last part of a series of publications and reports on the investigation of NP degradation mechanisms, this work further unfolds complex chemistry that occurs in the entire aging process of NP.
Fractional calculus is a relatively old yet emerging field of mathematics with the widest range of engineering and biomedical applications. Despite being an incredibly powerful tool, it, however, requires promotion in the engineering community. Rheology is undoubtedly one of the fields where fractional calculus has become an integral part of cutting-edge research. There exists extensive literature on the theoretical, experimental, and numerical treatment of various fractional viscoelastic flows in constraint geometries. However, the general theoretical approach that unites several most commonly used models is missing. Here we present exact analytical solutions for fractional viscoelastic flow in a circular pipe. We find velocity profiles and shear stresses for fractional Maxwell, Kelvin–Voigt, Zener, Poynting–Thomson, and Burgers models. The dynamics of these quantities are studied with respect to normalized pipe radius, fractional orders, and elastic moduli ratio. Three different types of behavior are identified: monotonic increase, resonant, and aperiodic oscillations. The models developed are applicable in the widest material range and allow for the alteration of the balance between viscous and elastic properties of the materials.
This report describes the activities during the fourth six month period of the investigation of acoustic propagation in the atmosphere with a realistic lapse temperature profile. A significant error was detected since the previous semi-annual report and has been corrected in both the plane wave and point source solutions. This report then describes both of these problems in some detail along with presenting some numerical results from the model. Work will begin this summer on the model of propagation in an inversion.