Search NASA⌕ Search

SEARCH · Search NASA

Results for “gradient calculation”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS1977

Effects of finite plasma beta on the lower-hybrid-drift instability

The local dispersion relation for the lower-hybrid-drift (LHD) instability is derived and analyzed, taking into account the finite-beta effects associated with transverse electromagnetic perturbations as well as with resonant and nonresonant electron-orbit modifications due to magnetic-field gradients. The influence of finite-beta effects on the LHD instability is calculated in a fully self-consistent manner for arbitrary values of electron-ion temperature ratio, local beta, cross-field ExB velocity/ion thermal speed ratio, and other plasma parameters. Stability properties are investigated analytically for the case of cold electrons, and the local dispersion relation is solved numerically in the parameter regime of most interest for high-density plasma pinches. The results show that for all parameter regimes studied, the net effect of finite plasma beta is to reduce the maximum growth rate of the LHD instability, although the details can vary, depending on the plasma parameters. Except in the limit where the electron/ion temperature ratio tends to zero, it is found that there is a critical value of plasma beta above which the LHD instability is completely stabilized.

Davidson, R. C.↗
NASA NTRS1983

Vegetation classification based on Advanced Very High Resolution Radiometer /AVHRR/ satellite imagery

Data from the NOAA-6 spacecraft Advanced Very High Resolution Radiometer (AVHRR) were tested for effectiveness for vegetation classification. Vegetation, climatological, and meteorological data were gathered for three days over 12 locations, and the normalized differences between the AVHRR bands 1 and 2 were determined. A vegetative greenness index was compared with a hydrologic factor and vegetation characteristics as measured by ground truth. A multivariate vegetation gradient model was formulated, incorporating AVHRR and climatological data. The hydrologic factor was calculated in terms of the precipitation, evaporation, maximum and minimum temperatures, and the hydrologic capacity. The observations were taken over Texas, which has a wide range of climates. A high correlation was found in the vegetation-HF index. The AVHRR data are concluded to be an effective tool for analysis of vegetation/climate relationships.

Norwine, J.↗
NASA NTRS1986

Cellular and dendritic growth in a binary melt - A marginal stability approach

A simple model for the constrained growth of an array of cells or dendrites in a binary alloy in the presence of an imposed positive temperature gradient in the liquid is proposed, with the dendritic or cell tip radius calculated using the marginal stability criterion of Langer and Muller-Krumbhaar (1977). This approach, an approach adopting the ad hoc assumption of minimum undercooling at the cell or dendrite tip, and an approach based on the stability criterion of Trivedi (1980) all predict tip radii to within 30 percent of each other, and yield a simple relationship between the tip radius and the growth conditions. Good agreement is found between predictions and data obtained in a succinonitrile-acetone system, and under the present experimental conditions, the dendritic tip stability parameter value is found to be twice that obtained previously, possibly due to a transition in morphology from a cellular structure with just a few side branches, to a more fully developed dendritic structure.

Laxmanan, V.↗
NASA NTRS1987

Meridional transport of magnetic flux in the solar wind between 1 and 10 AU - A theoretical analysis

Observations suggesting that the mean solar wind azimuthal field strength B(theta) near the ecliptic plane falls off more rapidly with heliocentric distance than would be expected in a classic Parker expansion is reexamined from a theoretical perspective using a three-dimensional MHD nonlinear numerical model for steady, corotating flow. For realistic solar wind parameters, it is found that a purely axisymmetric expansion can produce sizable magnetic flux deficits only when there are substantial meridional gradients in mean flow conditions localized about the ecliptic plane near the sun. Calculations on three-dimensional cororating flows are presented which demonstrate that latitudinal transport of magnetic flux by stream interactions may be an important consideration in generating the deficits in mean B(theta).

Pizzo, V. J.↗
NASA NTRS1988

All resistojet control of the NASA dual keel Space Station

This paper describes a control system that uses small (0.1-N) hydrogen-fueled, high-expansion-ratio resistojets to control a Space Station's position and attitude. The number and size of the resistojets is calculated based on requirements for stationkeeping, compensation for steady-state gravity-gradient and aerodynamic torques, momentum loading due to mobile remote manipulator motion, Shuttle Orbiter docking, and translational evasive maneuvering. The all-resistojet control system performs continuous orbit maintenance while keeping the linear accelerations on the spacecraft lower than they would be without any control system operating.

Paluszek, M. A.↗
NASA NTRS1993

A numerical study of confined turbulent jets

A numerical investigation is reported of turbulent incompressible jets confined in two ducts, one cylindrical and the other conical with a 5 degree divergence. In each case, three Craya-Curtet numbers are considered which correspond, respectively, to flow situations with no moderate and strong recirculation. Turbulence closure is achieved by using the k-epsilon model and a recently proposed realizable Reynolds stress algebraic equation model that relates the Reynolds stresses explicitly to the quadratic terms of the mean velocity gradients and ensures the positiveness of each component of the turbulent kinetic energy. Calculations are carried out with a finite-volume procedure using boundary-fitted curvilinear coordinates. A second-order accurate, bounded convection scheme and sufficiently fine grids are used to prevent the solutions from being contaminated by numerical diffusion. The calculated results are compared extensively with the available experimental data. It is shown that the numerical methods presented are capable of capturing the essential flow features observed in the experiments and that the realizable Reynolds stress algebraic equation model performs much better than the k-epsilon model for this class of flows of great practical importance.

Zhu, J.↗
NASA NTRS1959

Structure of Weak Shock Waves in a Monatomic Gas

The profiles and thicknesses of normal shock waves in argon at Mach numbers of 1.335, 1.454, 1.576, and 1-713 were determined experimentally by means of a free-molecule probe whose equilibrium temperature is related by kinetic theory to the local flow properties and their gradients. Comparisons were made between the experimental shock profiles and the theoretical profiles calculated from the Navier-Stokes equations, the Grad 13-moment equations, and the Burnett equations. New, very accurate numerical integrations of the Burnett equations were obtained for this purpose with results quite different from those found by Zoller, to whom the solution of this problem is frequently attributed. The experimental shock profiles were predicted with approximately equal success by the Navier-Stokes and Burnett theories, while the 13-moment method was definitely less satisfactory. A surprising feature of the theoretical results is the relatively small difference in predictions between the Navier-Stokes and Burnett theories in the present range of shock strengths and the contrastingly large difference between predictions of Burnett and the 13-moment theories. It is concluded that the Navier-Stokes equations are correct for weak shocks and that within the present shock strength range the Burnett equations make no improvement which merits the trouble of solving them. For shocks of noticeably greater strength, say with a shock Mach number of more than 2.5, it remains fundamentally doubtful that any of these theories can be correct.

Sherman, F. S.↗
NASA NTRS1997

Ring Current Modeling in a Realistic Magnetic Field Configuration

A 3-dimensional kinetic model has been developed to study the dynamics of the storm time ring current in a dipole magnetic field. In this paper, the ring current model is extended to include a realistic, time-varying magnetic field model. The magnetic field is expressed as the cross product of the gradients of two Euler potentials and the bounce-averaged particle drifts are calculated in the Euler potential coordinates. A dipolarization event is modeled by collapsing a tail-like magnetosphere to a dipole-like configuration. Our model is able to simulate the sudden enhancements in the ring current ion fluxes and the corresponding ionospheric precipitation during the substorm expansion.

Fok, M.-C.↗
NASA NTRS2006

Magnetic Field Solver

The Magnetic Field Solver computer program calculates the magnetic field generated by a group of collinear, cylindrical axisymmetric electromagnet coils. Given the current flowing in, and the number of turns, axial position, and axial and radial dimensions of each coil, the program calculates matrix coefficients for a finite-difference system of equations that approximates a two-dimensional partial differential equation for the magnetic potential contributed by the coil. The program iteratively solves these finite-difference equations by use of the modified incomplete Cholesky preconditioned-conjugate-gradient method. The total magnetic potential as a function of axial (z) and radial (r) position is then calculated as a sum of the magnetic potentials of the individual coils, using a high-accuracy interpolation scheme. Then the r and z components of the magnetic field as functions of r and z are calculated from the total magnetic potential by use of a high-accuracy finite-difference scheme. Notably, for the finite-difference calculations, the program generates nonuniform two-dimensional computational meshes from nonuniform one-dimensional meshes. Each mesh is generated in such a way as to minimize the numerical error for a benchmark one-dimensional magnetostatic problem.

Ilin, Andrew V.↗
DOE OSTI2014

Ultraviolet laser patterning of porous silicon

This work reports on the fabrication of 1D fringed patterns on nanostructured porous silicon (nanoPS) layers (563, 372, and 290 nm thick). The patterns are fabricated by phase-mask laser interference using single pulses of an UV excimer laser (193 nm, 20 ns pulse duration). The method is a single-step and flexible approach to produce a large variety of patterns formed by alternate regions of almost untransformed nanoPS and regions where its surface has melted and transformed into Si nanoparticles (NPs). The role of laser fluence (5–80 mJ cm{sup −2}), and pattern period (6.3–16 μm) on pattern features and surface structuring are discussed. The results show that the diameter of Si NPs increases with fluence up to a saturation value of 75 nm for a fluence ≈40 mJ cm{sup −2}. In addition, the percentage of transformed to non-transformed region normalized to the pattern period follows similar fluence dependence regardless the period and thus becomes an excellent control parameter. This dependence is fitted within a thermal model that allows for predicting the in-depth profile of the pattern. The model assumes that transformation occurs whenever the laser-induced temperature increase reaches the melting temperature of nanoPS that has been found to be 0.7 of that of crystalline silicon for a porosity of around 79%. The role of thermal gradients across the pattern is discussed in the light of the experimental results and the calculated temperature profiles, and shows that the contribution of lateral thermal flow to melting is not significant for pattern periods ≥6.3 μm.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Simplified Pore Network Model for Analysis of Flow Capacity and Residual Oil Saturation Distribution Based on Computerized Tomography

The pore network (connected tunnels) can be used not only to describe pore geometry inside a rock but also to calculate various fluid flow properties on a micro scale. Based on a “digital core” extracted from micro computerized tomography (μ-CT), we have developed a simplified pore network model. This study simplifies the pore structure as a lattice of cubes. Pore size, pore distribution, tunnel connectivity, fluid viscosity, interfacial tension, and the external driving pressure gradient are considered for each adjacent cube to evaluate the flow capacity of the pore network. We calculated the permeability, effective permeability, and effective sweep efficiency. The results show that fluid flow initially starts in the largest tunnel, which corresponds to the threshold pressure of the rock sample (fluid flow onset pressure). Permeability changes nonlinearly before all the tunnels are involved in the flow. Changes in the microscopic residual oil saturation and the sweep efficiency were visualized based on the three-dimensional pore distribution. Visualization shows that more oil remains in the formation, the heterogeneity is greater, and the oil viscosity is higher.

58 GEOSCIENCES↗
NASA NTRS1959

Charts and Tables for Estimating the Stability of the Compressible Laminar Boundary Layer with Heat Transfer and Arbitrary Pressure Gradient

The minimum critical Reynolds numbers for the similar solutions of the compressible laminar boundary layer computed by Cohen and Reshotko and also for the Falkner and Skan solutions as recomputed by Smith have been calculated by Lin's rapid approximate method for two-dimensional disturbances. These results enable the stability of the compressible laminar boundary layer with heat transfer and pressure gradient to be easily estimated after the behavior of the boundary layer has been computed by the approximate method of Cohen and Reshotko. The previously reported unusual result (NACA Technical Note 4037) that a highly cooled stagnation point flow is more unstable than a highly cooled flat-plate flow is again encountered. Moreover, this result is found to be part of the more general result that a favorable pressure gradient is destabilizing for very cool walls when the Mach number is less than that for complete stability. The minimum critical Reynolds numbers for these wall temperature ratios are, however, all larger than any value of the laminar-boundary-layer Reynolds number likely to be encountered. For Mach numbers greater than those for which complete stability occurs a favorable pressure gradient is stabilizing, even for very cool walls.

Tetervin, Neal↗
DOE OSTI2016

Theoretical study of NMR, infrared and Raman spectra on triple-decker phthalocyanines

Electronic structures and magnetic properties of multi-decker phthalocyanines were studied by theoretical calculation. Electronic structures, excited processes at multi-states, isotropic chemical shifts of {sup 13}C, {sup 14}N and {sup 1}H-nuclear magnetic resonance (NMR), principle V-tensor in electronic field gradient (EFG) tensor and asymmetry parameters (η), vibration mode in infrared (IR) and Raman spectra of triple-decker phthalocyanines were calculated by density functional theory (DFT) and time-dependent DFT using B3LYP as basis function. Electron density distribution was delocalized on the phthalocyanine rings with electron static potential. Considerable separation of chemical shifts in {sup 13}C, {sup 14}N and {sup 1}H-NMR was originated from nuclear spin interaction between nitrogen and carbon atoms, nuclear quadrupole interaction based on EFG and η of central metal under crystal field. Calculated optical absorption at multi-excited process was derived from overlapping π-orbital on the phthalocyanine rings. The vibration modes in IR and Raman spectra were based on in-plane deformation and stretching vibrations of metal-ligand coordination bond on the deformed structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Sensitivity study of hydrogen Soret transport in yttrium Hydride-Based nuclear fuel

Yttrium hydride is an excellent solid neutron moderator material for high temperature nuclear reactor applications due to its high hydrogen density and exceptional hydride stability at high temperatures. Despite these attractive characteristics, the details of how hydrogen behaves within yttrium hydride while temperature gradients exist are still not well understood. The evolution of the hydrogen composition profile resulting from a temperature gradient requires knowledge of hydrogen’s heat of transport, a critical parameter that has not yet been measured for this material. In this work, we perform hydride redistribution, hydrogen dissociation, and hydrogen leakage calculations while varying the Soret heat of transport of hydrogen in yttrium hydride to elucidate the sensitivity of hydride stability under temperature gradients to this parameter. This study analyzes hydride stability of a hypothetical uranium-yttrium hydride nuclear fuel design during operation of a high temperature liquid metal-cooled nuclear reactor. Assuming U-YH x could be fabricated in a physically stabilized manner, this fuel system can likely maintain hydride stability while operating at very high power densities and temperatures. We find that even though the hydrogen dissociation pressure in the gas gap does vary by several percent as the heat of transport temperature parameter is varied, the hydrogen content in the U-YH x fuel meat is relatively insensitive to this parameter over the course of a high burnup fuel cycle; this is due to yttrium hydride’s excellent hydrogen retention under the high temperature conditions considered here. Here, this suggests that hydride stability analyses are insensitive to the value of the Soret heat of transport in U-YH x under steady state liquid metal-cooled reactor conditions. However, the susceptibility to internal gas overpressurization-induced stress-rupture of the cladding during a high temperature transient is more sensitive to this parameter due to the non-linear dependence of hydrogen gas dissociation pressure vs. composition and temperature.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1986

Analysis of aerodynamic coefficients using gradient data: Spanwise turbulence effects on airplane response

The influence of spanwise turbulence on airplane dynamic behavior is determined. Calculations are based on data collected from test flights with the NASA B-57 research aircraft. The approach is to first compute aerodynamic forces and moments due to a spanwise distribution of angle of attack and airspeed. Secondly, these quantities are incorporated into the equations of motion. Simulation of flights done with the effects of spanwise turbulence included are compared to simulations without any spanwise turbulence. The findings of the study are that the moments developed by turbulence along the span are significant and that more realistic flight simulation can be achieved by including the spanwise turbulence terms.

Ringnes, E. A.↗
NASA NTRS1995

Radio wave scattering in the outer heliosphere

Current models for the 2-3 kHz emissions observed by the Voyager spacecraft in the outer heliosphere involve 2f(p) radiation generated near the termination shock or the heliopause. Radio wave scattering by solar wind density irregularities strongly affects observed sources of f(p) and 2f(p) emission in the inner heliosphere and the characteristics of astrophysical sources. In particular, the angular size, brightness temperature, and time variability of the source are strongly affected by scattering, thereby having major implications for the inferred size, energy budget, time variability, location, and nature of the source if scattering is ignored. This paper addresses whether scattering is important for interpreting the Voyager 2-3 kHz emissions. Quantitative calculations (with and without diffraction) are performed for the angular broadening of an outer heliospheric source as a function of path length, radiation frequency relative to f(p) and the spectrum of density irregularities. The effects of scattering in both the solar wind and the heliosheath are considered. Predictions for radial gradients in the source's apparent angular size and in the source's modulation index are presented. The calculations are compared with observations and the results discussed. First estimates suggest that scattering plausibly dominates the observed source size. The observed trend in modulation index with heliocentric distance is consistent with scattering being important and the source being in the outer heliosphere. Additional arguments for scattering being important are summarized.

Cairns, Iver H.↗
NASA NTRS1991

Numerical study of high-area-ratio H2/O2 rocket nozzles

Calculations were made for high-area-ratio H2/O2 rocket nozzles by using the RPLUS code which employs an implicit finite volume, lower-upper symmetric successive over-relaxation scheme to solve the Navier-Stokes equations and the species equations in a coupled manner. The combustion processes of hydrogen and oxygen are modeled by an 8-species and 18-step reaction mechanism, and the turbulence is simulated by the Baldwin-Lomax turbulence model with a pressure gradient correction on the van Driest's damping constant in the inner-layer model. The calculated results for two different rocket nozzles are presented as contours, profiles, and performance values. The predicted specific impulse and thrust agree well with the experimental data and those from the industry standard code. It is demonstrated that the RPLUS code can be used for analysis of rocket nozzles.

Kim, S. C.↗
DOE OSTI2014

Reduced model prediction of electron temperature profiles in microtearing-dominated National Spherical Torus eXperiment plasmas

A representative H-mode discharge from the National Spherical Torus eXperiment is studied in detail to utilize it as a basis for a time-evolving prediction of the electron temperature profile using an appropriate reduced transport model. The time evolution of characteristic plasma variables such as β{sub e}, ν{sub e}{sup ∗}, the MHD α parameter, and the gradient scale lengths of T{sub e}, T{sub i}, and n{sub e} were examined as a prelude to performing linear gyrokinetic calculations to determine the fastest growing micro instability at various times and locations throughout the discharge. The inferences from the parameter evolutions and the linear stability calculations were consistent. Early in the discharge, when β{sub e} and ν{sub e}{sup ∗} were relatively low, ballooning parity modes were dominant. As time progressed and both β{sub e} and ν{sub e}{sup ∗} increased, microtearing became the dominant low-k{sub θ} mode, especially in the outer half of the plasma. There are instances in time and radius, however, where other modes, at higher-k{sub θ}, may, in addition to microtearing, be important for driving electron transport. Given these results, the Rebut-Lallia-Watkins (RLW) electron thermal diffusivity model, which is based on microtearing-induced transport, was used to predict the time-evolving electron temperature across most of the profile. The results indicate that RLW does a good job of predicting T{sub e} for times and locations where microtearing was determined to be important, but not as well when microtearing was predicted to be stable or subdominant.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗