Search NASA⌕ Search

SEARCH · Search NASA

Results for “Higher order 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 127 records · Page 7

Target correlation effects on neutron-nucleus total, absorption, and abrasion cross sections

Second order optical model solutions to the elastic scattering amplitude were used to evaluate total, absorption, and abrasion cross sections for neutron nucleus scattering. Improved agreement with experimental data for total and absorption cross sections is found when compared with first order (coherent approximation) solutions, especially below several hundred MeV. At higher energies, the first and second order solutions are similar. There are also large differences in abrasion cross section calculations; these differences indicate a crucial role for cluster knockout in the abrasion step.

Cucinotta, Francis A.↗

Eddy mixing coefficient upper limit derived from the photochemical balance of O2

This work is based on the study of the photochemical balance of molecular oxygen in the martian atmosphere by using a one-dimensional model of photochemical reactions involving species derived from CO2 and H2O. The model is basically similar to one used previously for the study of the regulation of CO on a global scale, but the chemical rates are taken from another source. In the present scheme, the regulation of molecular oxygen is studied over timescales of the order of its photochemical lifetime (approximately equals 30 yr), which is much shorter than typical escape timescales. Thus, the escape fluxes are fixed to the values given by 3 and 4. We examine the calculated equilibrium abundances of O2 for given thermal, eddy diffusion coefficients and H2O profiles. The thermal profile is taken from in the lower atmosphere. At higher levels, in order to include the diurnal and seasonal thermal profile variability, we have also used the IRTM data. In order to study the influence of both temperature and pressure profiles on the O2 mixing ratio, we have made several tests corresponding to different martian seasons. The results show that the influence of pressure and temperature is quantitatively weak compared to the one of K and of the water vapor density (H2O). Thus, in the following we have fixed the pressure at the surface to a value of 7 mbar and we have used unique standard thermal profile corresponding to a profile roughly averaged over the year, the season, and the day: T equal 205 K at 0 km altitude, 175 K at 25 km, and 145 K at 50 km.

Rosenqvist, J.↗

Computation of hypersonic flow through a narrow expansion slot

The compressible Navier-Stokes equations are numerically solved for hypersonic flow over a three-dimensional ramp with a narrow expansion slot. In a two-dimensional test case, it is shown that a higher order scheme is needed to avoid too much numerical dissipation and hence too much total-pressure loss. The calculation demonstrates the role of viscosity in the expansion process of hypersonic flow through the narrow slot. Cases with various wall temperatures and slot widths are studied. Calculations show that wall cooling reduces the thickness of boundary layer, and hence increases the flow expansion substantially. In the lower portion of the slot, inviscidly, the flow is dominated by a highly expanded low density fluid, and viscously, by a viscous layer. As a direct consequence, as soon as the wedge angle is large enough, the mass flux through the slot is almost constant.

Hung, Ching-Mao↗

Radiation Heat Transfer Between Diffuse-Gray Surfaces Using Higher Order Finite Elements

This paper presents recent work on developing methods for analyzing radiation heat transfer between diffuse-gray surfaces using p-version finite elements. The work was motivated by a thermal analysis of a High Speed Civil Transport (HSCT) wing structure which showed the importance of radiation heat transfer throughout the structure. The analysis also showed that refining the finite element mesh to accurately capture the temperature distribution on the internal structure led to very large meshes with unacceptably long execution times. Traditional methods for calculating surface-to-surface radiation are based on assumptions that are not appropriate for p-version finite elements. Two methods for determining internal radiation heat transfer are developed for one and two-dimensional p-version finite elements. In the first method, higher-order elements are divided into a number of sub-elements. Traditional methods are used to determine radiation heat flux along each sub-element and then mapped back to the parent element. In the second method, the radiation heat transfer equations are numerically integrated over the higher-order element. Comparisons with analytical solutions show that the integration scheme is generally more accurate than the sub-element method. Comparison to results from traditional finite elements shows that significant reduction in the number of elements in the mesh is possible using higher-order (p-version) finite elements.

Gould, Dana C.↗

The neutron flux in space following a polar cap neutron event on November 15, 1960

A neutron detector flown on an Atlas pod during the 15 November 1960 solar event measured a larger neutron flux than another detector flown on a similar vehicle during a quiet solar period. The excess neutron flux at low latitudes obtained by comparing these two measurements is about an order of magnitude higher than expected from the neutron production in the atmosphere by solar protons as deduced from Lingenfelter's calculations and Webber's model for the 15 November event. It is suggested that most of the discrepancy can be accounted for by considering neutron production by energetic solar alpha particles and lowered geomagnetic cutoffs.

NEUTRON FLUX↗

An evaluation of a coupled microstructural approach for the analysis of functionally graded composites via the finite-element method

A comparison is presented between the predictions of the finite-element analysis and a recently developed higher-order theory for functionally graded materials subjected to a thorough-thickness temperature gradient. In contrast to existing micromechanical theories that utilize classical (i.e., uncoupled) homogenization schemes to calculate micro-level and macro-level stress and displacement fields in materials with uniform or nonuniform fiber spacing (i.e., functionally graded materials), the new theory explicitly couples the microstructural details with the macrostructure of the composite. Previous thermo-elastic analysis has demonstrated that such coupling is necessary when: the temperature gradient is large with respect to the dimension of the reinforcement; the characteristic dimension of the reinforcement is large relative to the global dimensions of the composite and the number of reinforcing fibers or inclusions is small. In these circumstances, the standard micromechanical analyses based on the concept of the representative volume element used to determine average composite properties produce questionable results. The comparison between the predictions of the finite-element method and the higher-order theory presented herein establish the theory's accuracy in predicting thermal and stress fields within composites with a finite number of fibers in the thickness direction subjected to a thorough-thickness thermal gradient.

Pindera, Marek-Jerzy↗

Study of diamond film growth and properties

The objective was to study diamond film growth and its properties in order to enhance the laser damage threshold of substrate materials. Calculations were performed to evaluate laser induced thermal stress parameter, R(sub T) of diamond. It is found that diamond has several orders of magnitude higher in value for R(sub T) compared to other materials. Thus, the laser induced damage threshold (LIDT) of diamond is much higher. Diamond films were grown using a microwave plasma enhanced chemical vapor deposition (MPECVD) system at various conditions of gas composition, pressure, temperature, and substrate materials. A 0.5 percent CH4 in H2 at 20 torr were ideal conditions for growing of high quality diamond films on substrates maintained at 900 C. The diamond films were polycrystalline which were characterized by scanning electron microscopy (SEM) and Raman scattering spectroscopy. The top surface of the growing film is always rough due to the facets of polycrystalline film while the back surface of the film replicates the substrate surface. An analytical model based on two dimensional periodic heat flow was developed to calculate the effective in-plane (face parallel) diffusivity of a two layer system. The effective diffusivity of diamond/silicon samples was measured using a laser pulse technique. The thermal conductivity of the films was measured to be 13.5 W/cm K, which is better than that of a type Ia natural diamond. Laser induced damage experiments were performed on bare Si substrates, diamond film coated Si, and diamond film windows. Significant improvements in the LIDT were obtained for diamond film coated Si compared to the bare Si.

Albin, Sacharial↗

Absorbed Dose and Dose Equivalent Calculations for Modeling Effective Dose

While in orbit, Astronauts are exposed to a much higher dose of ionizing radiation than when on the ground. It is important to model how shielding designs on spacecraft reduce radiation effective dose pre-flight, and determine whether or not a danger to humans is presented. However, in order to calculate effective dose, dose equivalent calculations are needed. Dose equivalent takes into account an absorbed dose of radiation and the biological effectiveness of ionizing radiation. This is important in preventing long-term, stochastic radiation effects in humans spending time in space. Monte carlo simulations run with the particle transport code FLUKA, give absorbed and equivalent dose data for relevant shielding. The shielding geometry used in the dose calculations is a layered slab design, consisting of aluminum, polyethylene, and water. Water is used to simulate the soft tissues that compose the human body. The results obtained will provide information on how the shielding performs with many thicknesses of each material in the slab. This allows them to be directly applicable to modern spacecraft shielding geometries.

Welton, Andrew↗

Temperature dependence of the formation of sulfate aerosols in the stratosphere

Classical nucleation theory is used in calculations of the temperature dependence of the characteristics and nucleation rates of sulfate aerosols in the binary H2SO4-H2O vapor mixture, in order to assess the influence of temperature on the formation of sulfate aerosols in the stratosphere, and to explore the possibility of new particle formation through homogeneous nucleation processes at regions where temperature is as low as -75 C, rather than the often-assumed -50 or -55 C. Calculation results indicate that the number of particles formed at a lower temperature is larger by several orders of magnitude than at higher temperatures, when water and sulfuric acid vapor concentrations are kept constant, and that large quantities of ultrafine particles which cannot be detected by conventional methods may exist at low-temperature stratospheric regions.

Yue, G. K.↗

Fluxes and distances of gamma-ray bursts

Gamma-gamma pair production and annihilation in gamma-ray burst source regions are calculated. These processes do not produce effective collimation of the flux above pair production threshold, but when repeated approximately 1/alpha times, higher order processes soften the spectrum and permit escape of intense fluxes.

Carrigan, B. J.↗

Boundary-layer equations in generalized curvilinear coordinates

A set of higher-order boundary-layer equations is derived valid for three-dimensional compressible flows. The equations are written in a generalized curvilinear coordinate system, in which the surface coordinates are nonorthogonal; the third axis is restricted to be normal to the surface. Also, higher-order viscous terms which are retained depend on the surface curvature of the body. Thus, the equations are suitable for the calculation of the boundary layer about arbitrary vehicles. As a starting point, the Navier-Stokes equations are derived in a tensorian notation. Then by means of an order-of-magnitude analysis, the boundary-layer equations are developed. To provide an interface between the analytical partial differentiation notation and the compact tensor notation, a brief review of the most essential theorems of the tensor analysis related to the equations of the fluid dynamics is given. Many useful quantities, such as the contravariant and the covariant metrics and the physical velocity components, are written in both notations.

Panaras, Argyris G.↗

Higher-order viscous shock-layer solutions for high altitude flows

A higher-order viscous shock-layer (HVSL) method has been developed and is used to obtain physically consistent results under varying degrees of low-density conditions for perfect-gas and nonequilibrium flows past long slender bodies. This method is implemented with higher-order body and shock slip conditions. Detailed comparisons with the direct simulation Monte Carlo (DSMC) method and Navier-Stokes (NS) calculations establish HVSL as a computationally efficient code for the low-density calculations. Further, it is shown that the deficiencies in the standard viscous shock-layer (VSL) method when predicting low-density flows lie in the governing equations and not with the slip boundary conditions as considered by earlier researchers.

Gupta, Roop N.↗

Morphology of nitric acid and water ice films

Ice films have been used to simulate stratospheric cloud surfaces in order to obtain laboratory data on solubilities and heterogeneous reaction rates. In the present study, environmental scanning electron microscopy (ESEM) is used to study thin films of both water ice and nitric acid ice near the composition of the trihydrate. The ices are formed by vapor deposition onto aluminum or borosilicate-glass substrates cooled to about 200 K. Micrographs are recorded during the deposition process and during subsequent annealing at higher temperatures. The results show that the ice films are composed of loosely consolidated granules, which range from about 1 to 20 microns in size at temperatures between 197 and 235 K. Cubic water ice is sometimes observed at 200 K, which converts to the hexagonal form at slightly higher temperatures. The loose packing of the granules confirms the high porosities of these films obtained from separate bulk porosity measurements. Average surface areas calculated from the observed granule sizes range from about 0.2 to 1 sq m/g and agree with surface areas obtained by gas-adsorption (BET) analysis of annealed ice films. For unannealed films, the BET areas are about an order of magnitude higher than the ESEM results, implying that the unannealed ices contain microporosity which is lost during the annealing process.

Keyser, Leon F.↗

Calculation of turbulent boundary layers with heat transfer and pressure gradient utilizing a compressibility transformation. Part 1: Summary report

The analysis uses a compressibility transformation and utilizes higher order closure rules to complete the transformation. By requiring that the momentum equations in differential form be satisfied at the wall and at the sublayer edge, correspondence rules are obtained which relate the variable property (VP) flow to a constant property (CP) flow in which mass transfer and pressure gradient occur simultaneously. A new CP formulation is developed and numerical results for a variety of cases are presented. Comparisons with earlier forms of the transformation and with experiment are included. For the zero pressure gradient case some differences between the various predictions are observed. For the several pressure gradient cases examined, the results are found to be essentially identical to those given by first order closure rules; i.e., by a form of transformation which relates the VP flow to a CP flow with pressure gradient but zero mass transfer.

Economos, C.↗

On the surface-to-bulk mode conversion of Rayleigh waves.

Surface-to-bulk wave conversion phenomena occurring at a discontinuity characterized by a surface contour deformation are shown to be usable as a means for tapping Rayleigh waves in a nonpiezoelectric solid. A boundary perturbation technique is used in the treatment of the mode conversion problem. A systematic procedure is presented for calculating not only the first-order scattered waves, which include the reflected surface wave and the converted bulk wave, but also the higher order terms.

Chang, C.-P.↗

Aerocapture and aeroassisted orbital transfer - A high performance leverage space technology

Single pass aerocapture and aeroassisted orbital transfer designs are discussed for modifying spacecraft orbits using planetary atmospheres. Aerodynamic drag can be used to shift a hyperbolic flyby trajectory to a planetary orbit by having the spacecraft altering its in-plane L/D ratio to maintain constant drag until the desired velocity is attained. The spacecraft would then move out of the atmosphere, sometimes with a small propulsive assist. Sample missions to Uranus, Jupiter, and Titan are mentioned, with mission scenarios commencing with a Shuttle launch. Biconic-shaped spacecraft are cited as having the best L/D for insertion accuracy, a low hypersonic ballistic coefficient, and volumeric efficiency. Aeroassisted orbital transfer involves atmospheric braking in order to shed velocity on a return trip from higher orbits. No decision has yet been made as to the spacecraft shape, but calculations indicate that aeroassist braking before return to LEO for rendezvous with the Shuttle has the potential for the lowest cost/mass to GEO of any delivery system currently under consideration.

Cruz, M. I.↗

The direct simulation of high-speed mixing-layers without and with chemical heat release

A direct numerical simulation of high speed reacting and non-reacting flows for H2-air systems is presented. The calculations are made for a convective Mach number of 0.38 with hyperbolic tangent initial profile and finite rate chemical reactions. A higher-order numerical method is used in time accurate mode to time advance the solution to a statistical steady state. About 600 time slices of all the variables are then stored for statistical analysis. It is shown that most of the problems of high-speed combustion with air are characterized by relatively weak heat release. The present study shows that: (1) the convective speed is reduced by heat release by about 10 percent at this convective Mach number M(sub c) = 0.38; (2) the variation of the mean and rms fluctuation of temperature can be explained on the basis of temperature fluctuation between the flame temperature and the ambient; (3) the growth rate with heat release is reduced by 7 percent; and (4) the entrainment is reduced by 25 percent with heat release. These differences are small in comparison with incompressible flow dynamics, and are argued to be due to the reduced importance of heat release in comparison with the large enthalpy gradients resulting from the large-scale vortex dynamics. It is finally suggested that the problems of reduced mixing in high-speed flows are not severely complicated by heat release.

Sekar, B.↗

Postbuckling of long thick plates in compression including higher order transverse shearing effects

Aluminum plates of various length/width ratios loaded in compression are examined analytically in terms of buckling and postbuckling for comparison with classical theory. The plates are considered to be subjected to longitudinal compressive displacement with no edge stress and are assumed to be simply supported at the side edges. The average longitudinal direct stresses are computed to plot characteristic curves by means of the following theories: 3D flexibility, higher-order shear deformation, first-order shear deformation, and classical von Karman. The 3D flexibility approach yields the lowest results with more unknowns than the higher-order shear-deformation method. The 3D flexibility approach is considered to be the most accurate, and it is shown that the calculated resultants and displacements tend to vary when consideration is given to the effects of transverse shear.

Stein, Manuel↗