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.
ODD-K angular flux tape computer program for radiation shielding design and analysis in nuclear rocket systems
Gas temperatures and pressures were measured around the second test position in the H2/O2 rocket engine at NASA-Lewis. Measured gas temperatures generally varied from 1210 to 1390 C. Measured pressures were in good agreement with other studies for throat tubes in a square chamber rocket engine. Heat transfer coefficients were measured at 90 and 180 degrees from the stagnation point and resulted in values of 27.5 and 8.5 kW/sq m C, respectively. A thermal model was developed to predict temperatures in bare and coated tubes and rods. Agreement between measured and predicted temperatures below the surface of a bare Mar-M 246 tube was very good for most of the heat up and cool down period. Predicted temperatures were significantly below measured temperatures for the coated tubes. A thermal model to simulate heat transfer to the leading edge of an HPFTP blade was developed and showed that TBCs can significantly dampen the thermal transient which occurs in the HPFTP during the startup of the SSME.
To achieve the high enthalpy conditions associated with hypersonic flight, many ground test facilities burn fuel in the air upstream of the test chamber. Unfortunately, the products of combustion contaminate the test gas and alter gas properties and the heat fluxes associated with aerodynamic heating. The difference in the heating rates between clean air and a vitiated test medium needs to be understood so that the thermal management system for hypersonic vehicles can be properly designed. This is particularly important for advanced hypersonic vehicle concepts powered by air-breathing propulsion systems that couple cooling requirements, fuel flow rates, and combustor performance by flowing fuel through sub-surface cooling passages to cool engine components and preheat the fuel prior to combustion. An analytical investigation was performed comparing clean air to a gas vitiated with methane/oxygen combustion products to determine if variations in gas properties contributed to changes in predicted heat flux. This investigation started with simple relationships, evolved into writing an engineering-level code, and ended with running a series of CFD cases. It was noted that it is not possible to simultaneously match all of the gas properties between clean and vitiated test gases. A study was then conducted selecting various combinations of freestream properties for a vitiated test gas that matched clean air values to determine which combination of parameters affected the computed heat transfer the least. The best combination of properties to match was the free-stream total sensible enthalpy, dynamic pressure, and either the velocity or Mach number. This combination yielded only a 2% difference in heating. Other combinations showed departures of up to 10% in the heat flux estimate.
This paper describes the hardware and software upgrades, operation, and performance of the high intensity diffractometer for residual stress analysis (HIDRA) instrument, a residual stress mapping neutron diffractometer located at the High Flux Isotope Reactor at Oak Ridge National Laboratory in Oak Ridge Tennessee, USA. Following a major upgrade in 2018, the new instrument has a single 3 He multiwire 2D 30 × 30 cm 2 position sensitive detector, yielding a field of view of 17° 2θ. Further, the increase in the field of view (from 4° 2θ) from the previous model instrument has contributed to the tremendous improvement in the out of plane solid angle such that the 3D count rate could be obtained easily. Accordingly, the hardware, software, Data Acquisition System (DAS), and so on have also been updated. Finally, all these enhanced features of HIDRA have been ably demonstrated by conducting multi directional diffraction measurements in the quenched 750-T74 aluminum, and the evolved and improved strain/stress mappings are presented.
Equations describing the temporal and spatial behavior of the kinematic moisture and heat flux are introduced. In these nonlinear equations, the contribution by diabatic processes to the large-scale flux is composed of two parts. One part is associated with a Rayleigh damping term, while the other arises from temporal and spatial changes in the pressure gradient term. The influence of diabatic processes on the large-scale moisture fluxes depends greatly on the degree of balance between forcing and damping terms in the governing equations. The existence of a near balance requires a reduction in the large-scale horizontal geostrophic wind speed. From a scale analysis of the moisture flux equations it is argued that reductions in the large-scale horizontal wind speed, observed within major cumulus cloud systems, help conserve large-scale moisture fluxes. The deviation of the wind from geostrophic conditions is easily estimated. This wind modification induces secondary vertical circulations that contribute to the convergence, creating or supporting long-lived mesoscale flows. In the tropics the wind modification has an antitriptic relationship. These diagnostic findings suggest possible modifications to the wind field in the application of cumulus parameterization, and may be important in diabatic initialization of numerical weather prediction models.
The present paper extends the analysis by Sakurai, Goossens, and Hollweg (1991) on resonant Alfven waves in nonuniform magnetic flux tubes. It proves that the fundamental conservation law for resonant Alfven waves found in ideal MHD by Sakurai, Goossens, and Hollweg remains valid in dissipative MHD. This guarantees that the jump conditions of Sakurai, Goossens, and Hollweg, that connect the ideal MHD solutions for xi(sub r), and P' across the dissipative layer, are correct. In addition, the present paper replaces the complicated dissipative MHD solutions obtained by Sakurai, Goossens, and Hollweg for xi(sub r), and P' in terms of double integrals of Hankel functions of complex argument of order 1/3 with compact analytical solutions that allow a straight- forward mathematical and physical interpretation. Finally, it presents an analytical dissipative MHD solution for the component of the Lagrangian displacement in the magnetic surfaces perpen- dicular to the magnetic field lines xi(sub perpendicular) which enables us to determine the dominant dynamics of resonant Alfven waves in dissipative MHD.
In this work, we investigate the ability of mesoscale and large-eddy simulation (LES) model configurations to predict the mean wind speed profile within the boundary layer of tropical cyclones (TCs). To this end, we perform idealized simulations of five hypothetical intense storms ranging from categories 1 to 5 on the Saffir–Simpson scale and extract time-averaged quantities near the eyewall region. We compare the model-generated data against mean wind speed profiles compiled from dropsondes launched from reconnaissance aircraft operating in the North Atlantic basin. Our analysis shows that mesoscale- and LES-generated mean wind fields display important differences in the boundary layer, including the magnitude of shear as well as the height where their low-level wind speed maxima are located. In addition, a comparison between the two model configurations with the dropsonde data shows that both modeling approaches are unable to capture the typical structure of mean winds in the lower part of the TC boundary layer (10–500 m), calling into question the use of simulations of near-axisymmetric storms for investigating the wind structure of past events. To better understand these differences, we conduct a momentum-budget analysis and show that modeled turbulent fluxes are underestimated in the mesoscale boundary layer parameterization compared to the LES model. Based on the analysis of the horizontal turbulent fluxes and their potential impact on mean flow quantities, a TC-specific boundary layer parameterization may be needed.
Neutrons have historically been used for a broad range of biological applications employing techniques such as small-angle neutron scattering, neutron spin echo, diffraction, and inelastic scattering. Unlike neutron scattering techniques that obtain information in reciprocal space, attenuation-based neutron imaging measures a signal in real space that is resolved on the order of tens of micrometers. The principle of neutron imaging follows the Beer-Lambert law and is based on the measurement of the bulk neutron attenuation through a sample. Greater attenuation is exhibited by some light elements (most notably, hydrogen), which are major components of biological samples. Contrast agents such as deuterium, gadolinium, or lithium compounds can be used to enhance contrast in a similar fashion as it is done in medical imaging, including techniques such as optical imaging, magnetic resonance imaging, X-ray, and positron emission tomography. For biological systems, neutron radiography and computed tomography have increasingly been used to investigate the complexity of the underground plant root network, its interaction with soils, and the dynamics of water flux in situ. Moreover, efforts to understand contrast details in animal samples, such as soft tissues and bones, have been explored. This manuscript focuses on the advances in neutron bioimaging such as sample preparation, instrumentation, data acquisition strategy, and data analysis using the High Flux Isotope Reactor CG-1D neutron imaging beamline. The aforementioned capabilities will be illustrated using a selection of examples in plant physiology (herbaceous plant/root/soil system) and biomedical applications (rat femur and mouse lung).
The thermal energy storage (TES) system has the capability to efficiently preserve thermal energy directly derived from the energy source, minimizing any conversion losses. Especially latent heat storage offers distinct advantages, including a substantial increase in energy storage density and minimization of temperature fluctuations within the plants. However, the phase change material (PCM) employed in latent heat storage has low thermal conductivity. Consequently, various studies are being conducted to enhance heat transfer. One approach to enhance heat transfer involves utilizing metal foam to maximize the heat transfer area. However, modeling metal foam with its intricate structure is a challenging task in numerical analysis. For this reason, ongoing research focuses on simplifying the modeling of metal foam. Nevertheless, fully encompassing all the characteristics of actual metal foam proves to be a challenging task for the simplified analytical model. The objective of this paper is to interpret the simple lattice metal foam analysis model from the perspective of behavior induced by buoyancy. When comparing the analysis results of solid PCM and liquid PCM with the same thermal conductivity under changes in porosity and gravity direction, we conducted an analysis to discern the trends in effective thermal conductivity that are overestimated due to convection. In the analysis, a constant heat flux of 10 kW and a constant surface boundary condition of 350 K were applied, and a sensitivity study regarding the mesh was conducted. The results indicate that, from the perspective of gravity in the simple lattice model, the solid analysis yields an effective thermal conductivity 29-47% higher compared to the liquid analysis. Additionally, as porosity increases, there is an observed increase of 24-33% in effective thermal conductivity.
An analytical technique for obtaining the time-resolved heat flux of a turbine blade is applied to the case of a TFE 731-2 hp full-stage rotating turbine. In order to obtain the heat flux values from the thin film gage temperature histories, a finite difference procedure is used to solve the heat equation with variable thermal properties. After setting out the data acquisition and analysis procedures, their application is illustrated for three midspan locations on the blade and operation at the design flow function. Results demonstrate that the magnitude of the heat flux fluctuation due to vane-balde interaction is large by comparison to the time-averaged heat flux at all investigated locations; FFT of a portion of the heat flux record illustrates that the dominant frequencies occur at the wake-cutting frequency and its harmonics.
We present the thermal infrared (5-35 micrometer) spectrum of 956 Elisa as measured by the Spitzer Infrared Spectrograph ("IRS"; Houck,1.R. et .11. [20041. Astrophys, 1. SuppL 154, 18-24) together with new ground-based lightcurve data and near-IR spectra. From the visible lightcurve photometry, we determine a rotation period of 16.494 +/- 0.001 h, identify the rotational phase of the Spitzer observations, and estimate the visible absolute magnitude (Hv) at that rotational phase to be 12.58 +/- 0.04. From radiometric analysis of the thermal flux spectrum, we find that at the time of observation 956 Elisa had a projected radius of 5.3 +/- 0.4 km with a visible albedo pv = 0.142+/- 0.022, significantly lower than that of the prototype V-type asteroid, 4 Vesta. (This corresponds to a radius of 5.2 +/- 0.4 km at lightcurve mean.) Analysis with the standard thermal model (STM) results in a sub-solar temperature of 292.3 +/- 2.8 K and beaming parameter eta = 1.16 +/- 0.05. Thermophysical modeling places a lower limit of 20 J m(exp -2)K(exp -1)s(exp -1/2) on the thermal inertia of the asteroid's surface layer (if the surface is very smooth) but more likely values fall between 30 and 150 J m(exp -2)K(exp -1)s(exp -1/2) depending on the sense of rotation. The emissivity spectrum, calculated by dividing the measured thermal flux spectrum by the modeled thermal continuum, exhibits mineralogically interpretable spectral features within the 9-12 micrometer reststrahlen band, the 15-16.5 micrometer Si-O-Si stretching region, and the 16-25 micrometer reststrahlen region that are consistent with pyroxene of diogenitic composition: extant diogenitic pyroxenes fall within the narrow compositional range W0(sub 2+/-1)En(sub 74+/-2)Fs(sub 24+/-1). Spectral deconvolution of the 9-12 micrometer reststrahlen features indicates that up to approximately 20% olivine may also be present, suggesting an olivine-diogenite-like mineralogy. The mid-IR spectrum is inconsistent with non-cumulate eucrite as the major component on the surface of 956 Elisa, although cumulate eucrite material may be present at abundances lower than that of the diogenite component. Analysis of new near-IR spectra of 956 Elisa with the Modified Gaussian Model (MGM; Sunshine, J,M., Pieters, C.M., Pratt, S.F.[1990]. J. Geophys. Res, 95 (May), 6955-6966) results in two pyroxene compositions: 75% magnesian low-Ca pyroxene and 25% high-Ca pyroxene. High-Ca pyroxene is not evident in the mid-IR data, but may belong to a component that is underrepresented in the mid-IR spectrum either because of its spatial distribution on the asteroid or because of its particle size. High-Ca pyroxenes that occur as exsolution lamellae may also be more evident spectrally in the NIR than in the mid-IR. In any case, we find that the mid-IR spectrum of 956 Elisa is dominated by emission from material of diogenite- like composition, which has very rarely been observed among asteroids.
Techniques for extraction boundary layer parameters from measurements of a short-pulse CO2 Doppler lidar are described. The measurements are those collected during the First International Satellites Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE). By continuously operating the lidar for about an hour, stable statistics of the radial velocities can be extracted. Assuming that the turbulence is horizontally homogeneous, the mean wind, its standard deviations, and the momentum fluxes were estimated. Spectral analysis of the radial velocities is also performed from which, by examining the amplitude of the power spectrum at the inertial range, the kinetic energy dissipation was deduced. Finally, using the statistical form of the Navier-Stokes equations, the surface heat flux is derived as the residual balance between the vertical gradient of the third moment of the vertical velocity and the kinetic energy dissipation. Combining many measurements would normally reduce the error provided that, it is unbiased and uncorrelated. The nature of some of the algorithms however, is such that, biased and correlated errors may be generated even though the raw measurements are not. Data processing procedures were developed that eliminate bias and minimize error correlation. Once bias and error correlations are accounted for, the large sample size is shown to reduce the errors substantially. The principal features of the derived turbulence statistics for two case studied are presented.
The near confluence of the successful launch of the Orbiting Carbon Observatory2 on July 2, 2014 and the acceptance on August 20, 2015 by Google, NASA Ames Research Center and USRA of a 1152 qubit D-Wave 2X Quantum Annealing Computer (QAC), offered an exceptional opportunity to explore the potential of this technology to address the scientific prediction of global annual carbon uptake by land surface processes. At UMBC,we have collected and processed 20 months of global Level 2 light CO2 data as well as fluorescence data. In addition we have collected ARM data at 2sites in the US and Ameriflux data at more than 20 stations. J. Dorband has developed and implemented a multi-hidden layer Boltzmann Machine (BM) algorithm on the QAC. Employing the BM, we are calculating CO2 fluxes by training collocated OCO-2 level 2 CO2 data with ARM ground station tower data to infer to infer measured CO2 flux data. We generate CO2 fluxes with a regression analysis using these BM derived weights on the level 2 CO2 data for three Ameriflux sites distinct from the ARM stations. P. Gentine has negotiated for the access of K34 Ameriflux data in the Amazon and is applying a neural net to infer the CO2 fluxes. N. Talik validated the accuracy of the BM performance on the QAC against a restricted BM implementation on the IBM Softlayer Cloud with the Nvidia co-processors utilizing the same data sets. G. Nearing and K. Harrison have extended the GSFC LIS model with the NCAR Noah photosynthetic parameterization and have run a 10 year global prediction of the net ecosystem exchange. C. Pellisier is preparing a BM implementation of the Kalman filter data assimilation of CO2 fluxes. At UMBC, R. Prouty is conducting OSSE experiments with the LISNoah model on the IBM iDataPlex to simulate the impact of CO2 fluxes to improve the prediction of global annual carbon uptake. J. LeMoigne and D. Simpson have developed a neural net image registration system that will be used for MODIS ENVI and will be converted to a BM algorithm implementation on the QAC. The first integer adder has been implemented on the D-Wave 2X by A. Shehab that will perform HAAR wavelets for image compression of MODIS scenes. Finally, based on the next generations of QACs, we are preparing a 5-year performance road map on the scalability of the current QAC algorithms.
Abstract The Escherichia coli genome‐scale metabolic model (GEM) is an exemplar systems biology model for the simulation of cellular metabolism. Experimental validation of model predictions is essential to pinpoint uncertainty and ensure continued development of accurate models. Here, we quantified the accuracy of four subsequent E. coli GEMs using published mutant fitness data across thousands of genes and 25 different carbon sources. This evaluation demonstrated the utility of the area under a precision–recall curve relative to alternative accuracy metrics. An analysis of errors in the latest (iML1515) model identified several vitamins/cofactors that are likely available to mutants despite being absent from the experimental growth medium and highlighted isoenzyme gene‐protein‐reaction mapping as a key source of inaccurate predictions. A machine learning approach further identified metabolic fluxes through hydrogen ion exchange and specific central metabolism branch points as important determinants of model accuracy. This work outlines improved practices for the assessment of GEM accuracy with high‐throughput mutant fitness data and highlights promising areas for future model refinement in E. coli and beyond.
Remote sensing of ionizing radiation has a significant role in waste management, nuclear material management and nonproliferation, and radiation safety. Robotic platforms can surpass the number of tasks that are achieved by humans. With this technique, the operator's radiation exposure can be decreased. Remote sensing allows for the evaluation and monitoring of radiological contamination. Gamma-ray and neutron sensors were integrated onto the robotic platforms. This approach allows for the radiation sensor data to be dynamically tracked and mapped thus enabling further analysis of the radiation flux in temporal and spatial domains. The goal is to complete scheduled tasks while the robot is being irradiated. To achieve this, electronic components must be shielded and radiation hardened. CZT Detector: Cadmium Zinc Telluride (CZT) detector technology has been a promising solution for gamma-ray and x-ray measurements. Detector data is transferred to the Odroid minicomputer that controls and powers the module via the USB. Robot Operating System (ROS) was utilized for data acquisition and data fusion. The Mariscotti method was employed for the spectrum analysis. A function was programmed in ROS for the automatic identification of photopeaks. CLYC Detector: A Cs{sub 2}LiYCl{sub 6}:Ce{sup 3+} (CLYC) detector was used for simultaneous medium-resolution gamma-ray measurements and neutron counting. A 2.54 cm diameter photomultiplier tube (PMT) was equipped with a high voltage supply and a miniature digitizer. Gamma-ray excitation: fast core-to-valence luminescence (CVL) with 1 ns decay constant, and prompt Ce{sup 3+} emission with 50 ns decay constant. Neutron excitation: slow cerium self-trapped excitation (Ce{sup 3+} STE), 1000 ns decay constant. Radiation Source Localization: Maximum Likelihood Estimation (MLE) and gradient-based methods were used to locate the position of a radiation source based on measured radiation intensities. Multi-Particle Transport Code FLUKA: Estimation of radiation damage of the electronic components is important in order to optimize the robot's operational time while it is irradiated. Displacement per atom (DPA) represents the radiation damage in materials exposed to the ionizing radiation. Various shielding layers of different thickness t were analyzed (< 5% statistical error). The model of the controller of the UAS was designed in FLUKA. Conclusion: CZT and CLYC detectors were integrated onto the robotic platforms. Radiation source localization and contour mapping using robotic platforms were studied. Functions for data analysis and fusion were developed in ROS. FLUKA code was utilized to analyze DPA values. Layers of low-density and high-density materials were used to shield the UAS electronics.
The BOUT++ edge plasma transport code is applied to study the effects of neutral recycling, drifts and heating scheme on edge plasma profiles and the divertor heat flux width of two discharges of the experimental advanced superconducting tokamak (EAST) steady-state H-mode plasmas heated by low hybrid wave (LHW) and neutral beam (NB), respectively. Neutral recycling seems to have played an important role in the plasma density profile. The edge plasma density drops dramatically for the case w/o neutral recycling, while it can be sustained for the case w/ neutral recycling. Drifts are found to have significant influences on edge plasma profiles. Both the amplitude and width of the divertor heat flux are found to have increased a lot due to drifts. The simulated heat flux width w/ drifts for the two discharges shows reasonable agreement with the experiments. However, the width from the simulation and experiment for the LHW heated discharge is much larger than that of the NB heated discharge. Comparison with Goldston's drift-based model and more detailed analysis on the heat flux contributions from drifts versus turbulence show that drifts are the dominant factor in edge plasma transport for both LHW and NB heated discharges, while turbulence may have played a more important role in determining the heat flux width in LHW heated discharges than that in NB heated discharges. This may account for the larger heat flux width in the LHW heated discharge. The magnetic topology and the equilibrium change by the LHW power may also be a potential reason for the larger heat flux width in the LHW heated discharges, which still needs more evidence and studies to support it. Finally, additional SOLPS simulation w/o drifts for the two discharges show reasonable agreement with the counterpart from the BOUT++ simulation, suggesting the two are both suitable codes for EAST edge plasma simulation.
Heat transfer between cladding and coolant during transient scenarios remains a critical area of uncertainty in understanding nuclear reactor safety. To advance the understanding of transient and accident scenarios involving critical heat flux (CHF), an in-pile experiment for the Transient Reactor Test (TREAT) facility at Idaho National Laboratory was developed. The experiment named CHF-Static Environment Rodlet Transient Test Apparatus (CHF-SERTTA), consists of a hollow borated stainless steel heater rod submerged in a static water pool which is heated via the (n, a) reaction in boron-10. This paper presents a novel inverse heat transfer method to determine CHF by using the optimization and uncertainty software DAKOTA to calibrate a RELAP5-3D model of CHF-SERTTA to temperature measurements obtained from a thermocouple welded to the surface of the rod.
Preliminary analysis of the momentum flux (or of the Poynting vector) of the classical electromagnetic version of the quantum vacuum consisting of zero-point radiation impinging on accelerated objects as viewed by an inertial observer suggests that the resistance to acceleration attributed to inertia may be a force of opposition originating in the vacuum. This analysis avoids the ad hoc modeling of particle-field interaction dynamics used previously by Haisck Rueda and Puthoff (1994) to derive a similar result. This present approach is not dependent upon what happens at the particle point but on how an external observer assesses the kinematical characteristics of the zero-point radiation impinging on the accelerated object. A relativistic form of the equation of motion results from the present analysis.