Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Sample Calculations”

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 73 records · Page 4

The effect of solid film lubricants on the stability of rolling element bearings

A technique has been developed for modeling solid films in rolling element bearings. For a given bearing geometry an effective viscosity can be calculated for a solid film as a function of the film shear modulus and ball-race friction coefficient. The calculated effective viscosity can subsequently be used as an input to a numerical model of cage motion and stability. Results from a sample calculation of effective viscosity and prediction of cage stability for a turbopump bearing are presented for films of molybdenum disulfide (MoS2) and polytetrafluoroethylene ()PTFE).

Kannel, J. W.↗

Calculation method using Clarkson integration for the physical dose at the center of the spread-out Bragg peak in carbon-ion radiotherapy

Purpose: In broad-beam carbon-ion radiotherapy performed using the heavy-ion medical accelerator in Chiba, the number of monitor units is determined by measuring the physical dose at the center of the spread-out Bragg peak (SOBP) for the treatment beam. The total measurement time increases as the number of treatment beams increases, which hinders the treatment of an increased number of patients. Hence, Kusano et al.[Jpn. J. Med. Phys. 23(Suppl. 2), 65-68 (2003)] proposed a method to calculate the physical dose at the center of the SOBP for a treatment beam. Based on a recent study, the authors here propose a more accurate calculation method.Methods: The authors measured the physical dose at the center of the SOBP while varying the circular field size and range-shifter thickness. The authors obtained the physical dose at the center of the SOBP for an irregularly shaped beam using Clarkson integration based on these measurements.Results: The difference between the calculated and measured physical doses at the center of the SOBP varied with a change in the central angle of the sector segment. The differences between the calculated and measured physical doses at the center of the SOBP were within {+-}1% for all irregularly shaped beams that were used to validate the calculation method.Conclusions: The accuracy of the proposed method depends on both the number of angular intervals used for Clarkson integration and the fineness of the basic data used for calculations: sampling numbers for the field size and thickness of the range shifter. If those parameters are properly chosen, the authors can obtain a calculated monitor unit number with high accuracy sufficient for clinical applications.

60 APPLIED LIFE SCIENCES↗

Modeling soot growth in Comp-B using a detonation history

A chemical kinetics procedure for soot formation calculation is summarized. By importing Lagrangian hydrodynamic calculations, species history is consistently accounted for. Detonation product composition is assumed in region with the burnt high explosive (HE). Soot formation is calculated using the density and energy obtained in the hydrodynamics calculation of HE detonation. Ideal gas equation-of-state (EOS) with temperature-dependent specific heat is used in the HE product region. Sample calculations show reasonably-expected temperature and soot formation.

42 ENGINEERING↗

Calculation of stability derivatives for slowly oscillating bodies of revolution at Mach 1.0

A parabolic method for steady transonic flow is extended to bodies of revolution oscillating in a sonic flow field. A Laplace transform technique is employed to derive the dipole solution, and the Adams-Sears iterative technique is used in the stability derivative calculation. A computer program is developed to perform the stability derivative calculation for the slowly oscillating cone and parabolic ogive. Inputs for the program are body geometry thickness ratio, acceleration constant, and pitch axis location. Sample calculations were performed for the parabolic ogive and circular cone and results are compared with those obtained by using other techniques and the available experimental data for circular cones.

Ruo, S. Y.↗

Fuselage design for a specified Mach-sliced area distribution

A procedure for designing a fuselage having a prescribed effective area distribution computed from -90 deg Mach slices is described. This type of calculation is an essential tool in designing a complete configuration with an effective area distribution that corresponds to a desired sonic boom signature shape. Sample calculations are given for M=2 and M=3 designs. The examples include fuselages constrained to have circular cross sections and fuselages having cross sections of arbitrary shape. It is found that, for a prescribed effective area distribution having sharp variations, the iterative procedure converges to a smoothed approximation to the prescribed distribution. For a smooth prescribed area distribution, the solution is not unique.

Barger, Raymond L.↗

Core saturation in a moving medium

A numerical technique for solving the line transfer equation of a two-level atom in static equilibrium is presented. Complete redistribution of emitted photons is assumed, as is saturation at the line core. Emission intensity is calculated either by a generalized Eddington-Barber relation, a first-order differential equation for the specific intensity, or by a formal transfer integral. Sample calculations are performed of the line transfer equation in a semi-infinite atmosphere with a constant Planck function of the collision parameter and for the Mg II resonance line in a model solar atmosphere experiencing shocks. Attention is focused on the line wings in the latter problem. The first order differential equation approach yields the best intensity values and temperature structure.

Kalkofen, W.↗

A Lanczos eigenvalue method on a parallel computer

Eigenvalue analyses of complex structures is a computationally intensive task which can benefit significantly from new and impending parallel computers. This study reports on a parallel computer implementation of the Lanczos method for free vibration analysis. The approach used here subdivides the major Lanczos calculation tasks into subtasks and introduces parallelism down to the subtask levels such as matrix decomposition and forward/backward substitution. The method was implemented on a commercial parallel computer and results were obtained for a long flexible space structure. While parallel computing efficiency for the Lanczos method was good for a moderate number of processors for the test problem, the greatest reduction in time was realized for the decomposition of the stiffness matrix, a calculation which took 70 percent of the time in the sequential program and which took 25 percent of the time on eight processors. For a sample calculation of the twenty lowest frequencies of a 486 degree of freedom problem, the total sequential computing time was reduced by almost a factor of ten using 16 processors.

Bostic, Susan W.↗

A Lanczos eigenvalue method on a parallel computer

Eigenvalue analyses of complex structures is a computationally intensive task which can benefit significantly from new and impending parallel computers. This study reports on a parallel computer implementation of the Lanczos method for free vibration analysis. The approach used here subdivides the major Lanczos calculation tasks into subtasks and introduces parallelism down to the subtask levels such as matrix decomposition and forward/backward substitution. The method was implemented on a commercial parallel computer and results were obtained for a long flexible space structure. While parallel computing efficiency is problem and computer dependent, the efficiency for the Lanczos method was good for a moderate number of processors for the test problem. The greatest reduction in time was realized for the decomposition of the stiffness matrix, a calculation which took 70 percent of the time in the sequential program and which took 25 percent of the time on eight processors. For a sample calculation of the twenty lowest frequencies of a 486 degree of freedom problem, the total sequential computing time was reduced by almost a factor of ten using 16 processors.

Bostic, Susan W.↗

Thermal-infrared remote sensing and Kirchhoff's Law. 1: Laboratory measurements

By the end of this century the Earth Observing System (EOS) will provide worldwide, thermal infrared, multispectral images of the Earth, presenting geologists with a new kind of remote sensing data for interpretation. Thus it has become essential to understand the spectral emittance behavior of terrestrial surface materials. Perhaps the most fundamental question to be answereed is the extent to which such materials follow Kirchhoff's law (epsilon = 1 -R) under laboratory and field conditions, especially when a sample displays a thermal gradient. We present the first rigorous quantitative comparison of directional and hemispherical reflectance and directional emittance of rock and soil samples in the laboratory, with thermal gradients induced by heating them from below and allowing them to radiate to a colder background. The results show that only an extemeley low density sample composed of fine particles sifted into a 'fairy castle' structure displays a thermal gradient steep enough within the infrared skin depth to cause significant (6%) departure from Kirchhoff's law. There is no detectable effect on the more normal terrestrial samples, such as soils and rocks measured in the laboratory, even when semitransparent coatings are involved. Thus both emittance and reflectance measurements can be used to calculate sample emissivity for most terrestrial surface materieals. However, the effect on Kirchhoffian behavior of different field environments, which may induce a steeper thermal gradient in particulate samples, has yet to be determined, and some low-density surface materials like newly fallen snow, frost, and efflorescent salts on playas have yet to be measured in emittance.

Salisbury, John W.↗

Application of a convolutional neural network in permeability prediction: A case study in the Jacksonburg-Stringtown oil field, West Virginia, USA

Permeability is a critical parameter for understanding subsurface fluid flow behavior, managing reservoirs, enhancing hydrocarbon recovery, and sequestering carbon dioxide. In general, permeability is measured in the laboratory based on subsurface core samples, calculated from well logs or estimated from well tests. However, laboratory measurements and well tests are expensive, time-consuming, and usually limited to a few core samples or wells in a hydrocarbon field or carbon storage site. Machine-learning techniques are good options for generating a rapid, robust, and cost-effective permeability prediction model because of their strengths to recognize the potential interrelationships between input and output variables. Convolutional neural networks (CNN), as a good pattern recognition algorithm, are widely used in image processing, natural language processing, and speech recognition, but are rarely used with regression problems and even less often in reservoir characterization. We have developed a CNN regression model to estimate the permeability in the Jacksonburg-Stringtown oil field, West Virginia, which is a potential carbon storage site and enhanced oil recovery operations field. We also evaluate the concept of the geologic feature image, which is converted from geophysical well logs. Five variables, including two commonly available conventional well logs (the gamma rays [GRs] and bulk density) and three well-log-derived variables (the slopes of the GR and bulk density curves, and shale content), are used to generate a geologic feature image. The CNN treats the geologic feature image as the input and the permeability as the desired output. In addition, the permeability predicted using traditional backpropagation artificial neural networks, which are optimized by genetic algorithms and particle swarm optimization, is compared with the permeability estimated using our CNN. Our results indicate that the CNN regression model provides more accurate permeability predictions than the traditional neural network.

Geochemistry & Geophysics↗

Semiempirical fragmentation models on galactic cosmic ray transport calculations with hydrogen target

Nuclear fragmentation cross sections of Silberberg and Tsao that are more accurate for a hydrogen target were implemented in the data base to replace those of Rudstam for a galactic cosmic ray transport code (HZETRN). Sample calculations were made for the transported galactic cosmic ray flux through a liquid hydrogen shield at solar minimum condition to determine the effect of such a change. The transported flux based on the Silberberg-Tsao semiempirical formalism contains fewer high-LET (linear energy transfer) components but more low-LET components than the results based on Rudstam's formalism: and this disparity deepens as the shield thickness increases. A comparison of the results obtained from using both energy-dependent and energy-independent cross sections of Silberberg and Tsao indicates that the energy-independent assumption results in an underestimation of high-LET flux above 100 keV/micron by approximately 40 percent for a 15-g/cm(sup 2) thickness of liquid hydrogen. Similar results were obtained in a previous study when both energy-dependent and energy-independent cross sections of Rudstam were considered. Nonetheless, the present study found that an energy-independent calculation would be best accomplished by using Rudstam's cross sections as done in the past for various engineering applications.

Shinn, Judy L.↗

A time dependent difference theory for sound propagation in ducts with flow

A time dependent numerical solution of the linearized continuity and momentum equation was developed for sound propagation in a two dimensional straight hard or soft wall duct with a sheared mean flow. The time dependent governing acoustic difference equations and boundary conditions were developed along with a numerical determination of the maximum stable time increments. A harmonic noise source radiating into a quiescent duct was analyzed. This explicit iteration method then calculated stepwise in real time to obtain the transient as well as the steady state solution of the acoustic field. Example calculations were presented for sound propagation in hard and soft wall ducts, with no flow and plug flow. Although the problem with sheared flow was formulated and programmed, sample calculations were not examined. The time dependent finite difference analysis was found to be superior to the steady state finite difference and finite element techniques because of shorter solution times and the elimination of large matrix storage requirements.

Baumeister, K. J.↗

A method for calculating heat transfer in the laminar flow region of bodies

This report has been prepared to provide a practical method for determining the chordwise distribution of the rate of heat transfer from the surface of a wing or body of revolution to air. The method is limited in use to the determination of heat transfer from the forward section of such bodies when the flow is laminar. A comparison of the calculated average heat-transfer coefficient for the nose section of the wing of a Lockheed 12-A airplane with that experimentally determined shows a satisfactory agreement. A sample calculation is appended.

Allen, H Julian↗

2012 Accomplishments - Tritium Aging Studies on Stainless Steels

This report summarizes the research and development accomplishments during FY12 for the tritium effects on materials program. The tritium effects on materials program is designed to measure the long-term effects of tritium and its radioactive decay product, helium-3, on the structural properties of forged stainless steels which are used as the materials of construction for tritium reservoirs. The FY12 R&D accomplishments include: (1) Fabricated and Thermally-Charged 150 Forged Stainless Steel Samples with Tritium for Future Aging Studies; (2) Developed an Experimental Plan for Measuring Cracking Thresholds of Tritium-Charged-and-Aged Steels in High Pressure Hydrogen Gas; (3) Calculated Sample Tritium Contents For Laboratory Inventory Requirements and Environmental Release Estimates; (4) Published report on “Cracking Thresholds and Fracture Toughness Properties of Tritium-Charged-and-Aged Stainless Steels”; and, (5) Published report on “The Effects of Hydrogen, Tritium, and Heat Treatment on the Deformation and Fracture Toughness Properties of Stainless Steels”. These accomplishments are highlighted here and references given to additional reports for more detailed information.

36 MATERIALS SCIENCE↗

Melting behavior and phase relations of lunar samples

Cooling rate studies of 12002 were conducted and the results interpreted in terms of the crystallization history of this rock and certain other picritic Apollo 12 samples. Calculations of liquid densities and viscosities during crystallization, crystal settling velocities, and heat loss by the parent rock body are discussed, as are petrographic studies of other Apollo 12 samples. The process of magmatic differentiation that must have accompanied the early melting and chemical fractionation of the moon's outer layers was investigated. The source of regions of both high- and low-titanium mare basalts were also studied.

Hays, J. F.↗

A time dependent difference theory for sound propagation in ducts with flow

A time dependent numerical solution of the linearized continuity and momentum equation is developed for sound propagation in a two-dimensional straight hard or soft wall duct with a sheared mean flow. The time dependent governing acoustic-difference equations and boundary conditions are developed along with a numerical determination of the maximum stable time increments. The analysis begins with a harmonic noise source radiating into a quiescent duct. This explicit iteration method then calculates stepwise in real time to obtain the transient as well as the 'steady' state solution of the acoustic field. Example calculations are presented for sound propagation in hard and soft wall ducts, with no flow and with plug flow. Although the problem with sheared flow has been formulated and programmed, sample calculations have not yet been examined. So far, the time dependent finite difference analysis has been found to be superior to the steady state finite difference and finite element techniques because of shorter solution times and the elimination of large matrix storage requirements.

Baumeister, K. J.↗

Mass Spectrometric Studies of Oxides

Current studies at NASA Glenn on oxide thermodynamics are discussed. Previous studies on the vaporization of B2O3 in reducing atmospheres led to inconsistent studies when B was used as a reductant. It is shown that liquid B2O3 does not wet B and a clear phase separation was noted in the Knudsen cell. This problem was solved by using FeB and Fe2B to supply a different and constant activity of B. The thermodynamic data thus derived are compared to quantum chemical composite calculations. A major problem in high temperature mass spectrometry is the determination of accurate ionization cross sections, particularly for molecules. The method of Deutsch and Mark shows promise and some sample calculations are discussed. Finally current studies on the thermodynamics of rare earth silicates are discussed. Here the problems are obtaining a measurable signal from SiO2 vaporization and non-equilibrium vaporization. The use of a Ta reducing agent provides a stronger signal, which is related to silica activity. The Whitman-Motzfeld relation adapted to KEMS measurements is applied to obtain equilibrium pressures.

Jacobson, Nathan S.↗

Viscous/inviscid interaction procedure for high-amplitude oscillating airfoils

A coupling procedure between an inviscid potential theoretical panel method and an unsteady interactive boundary layer calculation procedure has been developed for calculating unsteady airloads on oscillating airfoils in the high-incidence regime. A coordinate transformation into a time-dependent surface coordinate has been incorporated into the boundary layer code to avoid a quasi-steady treatment of the viscous flow equations in the leading edge region. The upper and lower sides of the airfoil are calculated with the complete set of unsteady equations. The interactive region of the boundary layer calculation can be straightforwardly extended over the complete upper surface of the airfoil while accounting for interaction areas in the region of strong adverse pressure gradients with incipient separation. Comparisons of sample calculations with experimental data show reasonably good agreement even for unsteady drag coefficients.

Geissler, W.↗