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.
Forced-convection heat transfer analysis for liquid metal flow in rectangular channels with heat source within fluid stream
By measurement and analysis, the relationship between the distribution of the outflow of acoustic energy over the jet boundary and the far-field intensity is considered. The physical quantity used is the gradient of the pressure evaluated on a geometrical plane at the smallest possible radial distance from the jet axis, but outside the vortical region, in the area where the homogeneous wave equation is reasonably well satisfied. The numerical and experimental procedures involved have been checked out by using a known source. Results indicate that the acoustic power output per unit length of the jet, in the region from which the sound emanates, peaks at approximately 9 diameters downstream. The acoustic emission for a jet Strouhal number of about 0.3 exceeds the emission for all other Strouhal numbers nearly everywhere along the measurement plane. However, the far-field peak intensity distribution obtained from the contribution of each station was found to depend on the spatial extent of the region where sound emanates from the jet, which, in turn, depends more on the far-field angle than on the Strouhal number.
Reference is made to the measurements of a cosmic ray antiproton flux at a few hundred MeV reported by Buffington et al. (1981), noting that one of the final background processes to be removed by the data analysis in that study was helium-induced events which satisfied the criteria for topology and timing. The response in the third scintillator S3 was used to identify and remove these events. For the top two scintillators S1 and S2, pulse size information was lost during the data-taking. A method is reported here for the partial retrieval of pulse size information for the scintillator S2. This is possible because a portion of this signal was subtracted from the Cerenkov response before trigger discrimination and data recording to remove scintillation from the Cerenkov response. For separating protons from more highly charged particles, the method is considered sufficient. It is pointed out that the sample of events identified as antiprotons, for which the method can be applied, has the expected unit charge in scintillator S2.
Modern numerical techniques employing properties of flux Jacobian matrices are extended to general, equilibrium gas laws. Generalizations of the Beam-Warming scheme, Steger-Warming and van Leer flux-vector splittings, and Roe's approximate Riemann solver are presented for three-dimensional, time-varying grids. The approximations inherent in previous generalizations are discussed.
This paper examines the behavior of flux and slope limiters on non-uniform grids in multiple dimensions. Many slope limiters in standard use do not preserve linear solutions on irregular grids impacting both accuracy and convergence. We rewrite some well-known limiters to highlight their underlying symmetry, and use this form to examine the proper - ties of both traditional and novel limiter formulations on non-uniform meshes. A consistent method of handling stretched meshes is developed which is both linearity preserving for arbitrary mesh stretchings and reduces to common limiters on uniform meshes. In multiple dimensions we analyze the monotonicity region of the gradient vector and show that the multidimensional limiting problem may be cast as the solution of a linear programming problem. For some special cases we present a new directional limiting formulation that preserves linear solutions in multiple dimensions on irregular grids. Computational results using model problems and complex three-dimensional examples are presented, demonstrating accuracy, monotonicity and robustness.
This article presents three-dimensional B field solutions for the cylindrical Halbach array in an axial orientation. This arrangement has applications in the design of axial motors and passive axial magnetic bearings and couplers. The analytical model described here assumes ideal magnets with fixed and uniform magnetization. The field component functions are expressed as sums of 2-D definite integrals that are easily computed by a number of mathematical analysis software packages. The analysis is verified with sample calculations and the results are compared to equivalent results from traditional finite-element analysis (FEA). The field solutions are then approximated for use in flux linkage and induced EMF calculations in nearby stator windings by expressing the field variance with angular displacement as pure sinusoidal function whose amplitude depends on radial and axial position. The primary advantage of numerical implementation of the analytical approach presented in the article is that it lends itself more readily to parametric analysis and design tradeoffs than traditional FEA models.
The surface temperature distributions due to thermocapillary convections in a thin liquid layer with heat fluxes imposed on the free surface were investigated. The nondimensional analysis predicts that, when convection is important, the characteristics length scale in the flow direction L, and the characteristic temperature difference delta T sub o can be represented by L and delta T sub o approx. (A2Ma)/1/4 delta T sub R, respectively, where L sub R and delta sub R are the reference scales used in the conduction dominant situations with A denoting the aspect ratio and Ma the Marangoni number. Having L and delta sub o defined, the global surface temperature gradient delta sub o/L, the global thermocapillary driving force, and other interesting features can be determined. Numerical calculations involving a Gaussian heat flux distribution are presented to justify these two relations.
The surface temperature distributions due to thermocapillary convections in a thin liquid layer with heat fluxes imposed on the free surface were investigated. The nondimensional analysis predicts that, when convection is important, the characteristics length scale in the flow direction L, and the characteristic temperature difference delta T sub o can be represented by L and delta T sub o approx. (A2Ma)/1/4 delta T sub R, respectively, where L sub R and delta sub R are the reference scales used in the conduction dominant situations with A denoting the aspect ratio and Ma the Marangoni number. Having L and delta sub o defined, the global surface temperature gradient delta sub o/L, the global thermocapillary driving force, and other interesting features can be determined. Numerical calculations involving a Gaussian heat flux distribution are presented to justify these two relations.
The relation between rotational modulation of the ultraviolet solar irradiance and variations in atmospheric ozone has been investigated using Fourier transform harmonic analysis and cross-correlations. Ozone variations with the same period and phase as 13.5 day or 27-day solar flux variations occur at tropical and subtropical latitudes over a range of pressure levels centered about 3 mbar. The solar-forced oscillation is stronger in the summer hemisphere; as temperature-related variations would be stronger in winter. Changes in solar irradiance over the 11-year cycle can be estimated by scaling rotational modulation. Using this estimate and the ozone-sun relation obtained for rotational modulation yields solar cycle changes of 3.5 percent in 3 percent mixing ratio comparable to that predicted from halocarbons and 0.7 mbar in total ozone.
It is known from computer calculations that if a force-free magnetic-field configuration is stressed progressively by footpoint displacements, the configuration expands and approaches the open configuration with the same surface flux distribution, and, in the process, the energy of the field increases progressively. Analysis of a simple model of force-free fields of cylindrical symmetry leads to simple asymptotic expressions for the extent and energy of such a configuration. The analysis is carried through for both spherical and planar source surfaces. According to this model, the field evolves in a well-behaved manner with no indication of instability or loss of equilibrium.
We present a broadband (approx. 0.5 - 79 keV) spectral and temporal analysis of multiple NuSTAR observations combined with archival Suzaku and Chandra data of NGC4945, the brightest extragalactic source at 100 keV. We observe hard X-ray (> 10 keV) flux and spectral variability, with flux variations of a factor 2 on timescales of 20 ksec. A variable primary continuum dominates the high energy spectrum (> 10 keV) in all the states, while the reflected/scattered flux which dominates at E< 10 keV stays approximately constant. From modelling the complex reflection/transmission spectrum we derive a Compton depth along the line of sight of Thomson approx.2.9, and a global covering factor for the circumnuclear gas of approx. 0.15. This agrees with the constraints derived from the high energy variability, which implies that most of the high energy flux is transmitted, rather that Compton-scattered. This demonstrates the effectiveness of spectral analysis in constraining the geometric properties of the circumnuclear gas, and validates similar methods used for analyzing the spectra of other bright, Compton-thick AGN. The lower limits on the e-folding energy are between 200 − 300 keV, consistent with previous BeppoSAX, Suzaku and Swift BAT observations. The accretion rate, estimated from the X-ray luminosity and assuming a bolometric correction typical of type 2 AGN, is in the range approx. 0.1 − 0.3 lambda(sub Edd) depending on the flux state. The substantial observed X-ray luminosity variability of NGC4945 implies that large errors can arise from using single-epoch X-ray data to derive L/L(sub Edd) values for obscured AGNs.
Vertically integrated back and forward trajectories for the 300-200, 700-500 and surface-800 hPa levels are calculated using Pretoria as point of origin for the Southern Africa Fire-Atmosphere Research Initiative (SAFARI) period September-October 1992. The transport fields are then combined to show both horizontal and vertical transport of air to and from Pretoria at the different levels. Air transport patterns in the vertical are linked to the occurrence of absolutely stable layers which are also evident in the 16 ozonesonde profiles recorded at Pretoria during SAFARI. The coherence of the stratification based on dynamical and ozone analysis permits the use of mean ozone profiles with air volume fluxes to interpret the ozone in terms of photochemistry and transport within stable layers. Extensive recirculation across the meridional plane at Pretoria implies that advection of ozone is slow and that photochemistry is responsible for the observed vertical structure over central southern Africa in September and October 1992. Requisite ozone formation rates are supported by model analysis of ozone and ozone precursors measured from SAFARI and Transport and Atmospheric Research Chemistry near the Equator-Atlantic aircraft.
Analysis has been performed for MAVEN mission. Due to the elliptical orbit, large pressure variations in orbit will be experienced, there is a need to understand how internal pressures change and the flux of gas from vents could potentially bias instrument measurements. Goal of this analysis is to predict the effect that atmospheric gases trapped and vented from spacecraft volumes could have on instrument measurements.
The mixing and heat transfer phenomena in small rocket flow fields with fuel film cooling is not well understood. An instrumented, water-cooled chamber with a gaseous hydrogen/gaseous oxygen injector was used to gather steady-state inner and outer wall temperature profiles. The chamber was tested at 414 kPa (60 psia) chamber pressure, from mixture ratios of 3.41 to 8.36. Sixty percent of the fuel was used for film cooling. These temperature profiles were used as boundary conditions in a finite element analysis program, MSC/NASTRAN, to calculate the local radial and axial heat fluxes in the chamber wall. The normal heat fluxes were then calculated and used as a diagnostic of the rocket's flow field behavior. The normal heat fluxes determined were on the order of 1.0 to 3.0 MW/meters squared (0.6 to 1.8 Btu/sec-inches squared). In the cases where mixture ratio was 5 or above, there was a sharp local heat flux maximum in the barrel section of the chamber. This local maximum seems to indicate a reduction or breakdown of the fuel film cooling layer, possibly due to increased mixing in the shear layer between the film and core flows. However, the flow was thought to be completely laminar, as the throat Reynolds numbers were below 50,000 for all the cases. The increased mixing in the shear layer in the higher mixture ratio cases appeared not to be due to the transition of the flow from laminar to turbulent, but rather due to increased reactions between the hydrogen film and oxidizer-rich core flows.
The behavior of in situ measurements of surface fluxes obtained during FIFE 1987 is examined by using correlative and spectral techniques in order to assess the significance of fluctuations on various time scales, from subdiurnal up to synoptic, intraseasonal, and annual scales. The objectives of this analysis are: (1) to determine which temporal scales have a significant impact on areal averaged fluxes and (2) to design a procedure for filtering an extended flux time series that preserves the basic diurnal features and longer time scales while removing high frequency noise that cannot be attributed to site-induced variation. These objectives are accomplished through the use of a two-dimensional cross-time Fourier transform, which serves to separate processes inherently related to diurnal and subdiurnal variability from those which impact flux variations on the longer time scales. A filtering procedure is desirable before the measurements are utilized as input with an experimental biosphere model, to insure that model based intercomparisons at multiple sites are uncontaminated by input variance not related to true site behavior. Analysis of the spectral decomposition indicates that subdiurnal time scales having periods shorter than 6 hours have little site-to-site consistency and therefore little impact on areal integrated fluxes.
This study examines oceanic surface heat fluxes in the Norwegian, Greenland, and Barents seas using the gridded Navy Fleet Numerical Oceanography Central surface analysis and the First GARP Global Experiment (FGGE) IIc cloudiness data bases. Monthly and annual means of net and turbulent heat fluxes are computed for the FGGE year 1979. The FGGE IIb data base consisting of individual observations provides particularly good data coverage in this region for a comparison with the gridded Navy winds and air temperatures. The standard errors of estimate between the Navy and FGGE IIb winds and air temperatures are 3.6 m/s and 2.5 C, respectively. The computations for the latent and sensible heat fluxes are based on bulk formulas with the same constant heat exchange coefficient of 0.0015. The results show extremely strong wintertime heat fluxes in the northern Greenland Sea and especially in the Barents Sea in contrast to previous studies.
Considerable progress has been made in recent years toward advancing the state-of-the-art in solid mechanics boundary element technology. In the present work, much of this new technology is applied in the development of a general-purpose boundary element method (BEM) for steady-state heat conduction. In particular, the BEM implementation involves the use of higher-order conforming elements, self-adaptive integration and multi-region capability. Two- and three-dimensional, as well as axisymmetric analysis, are incorporated within a unified framework. In addition, techniques are introduced for the calculation of boundary flux, and for the inclusion of thermal resistance across interfaces. As a final extension, an efficient formulation is developed for the analysis of solid three-dimensional bodies with embedded holes. For this last class of problems, the new BEM formulation is particularly attractive, since use of the alternatives (i.e. finite element or finite difference methods) is not practical. A number of detailed examples illustrate the suitability and robustness of the present approach for steady-state heat conduction.
A magnetic flux rope model is developed and used to determine whether the principal axis analysis (PAA) of magnetometer signatures from a single satellite pass is sufficient to obtain the magnetic topology of plasmoids. The model is also used to determine if plasmoid observations are best explained by the flux rope, closed loop, or large-amplitude wave picture. It was found that the principal axis directions is highly dependent on the satellite trajectory through the structure and, therefore, the PAA of magnetometer data from a single satellite pass is insufficient to differentiate between magnetic closed loop and flux rope models. Results also indicate that the flux rope model of plasmoid formation is well suited to unify the observations of various magnetic structures observed by ISEE 3.