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 1,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Integrated modeling of advanced tokamaks for tearing mode avoidance on DIII-D

Understanding the physics of the onset and evolution of tearing modes (TMs) in tokamak plasmas is crucial for high-performance steady-state operations. The onset condition of (m,n) = (2,1) tearing stability, $Δ′$ > $Δ_c^′$> 0, has been studied in DIII-D steady-state hybrid discharges with accurate equilibrium reconstruction and well-measured plasma profiles. Here, m is the poloidal mode number and n is the toroidal mode number. The onset of n = 1 TMs is carefully identified using the Mirnov array to determine in-phase signals for the onset time and using electron cyclotron emission measurements to trace back from the saturated island width in the experiments for both initial island width and its growth rate. The tearing stability index $Δ′$ at TM onset is calculated in multiple ways, using the modified Rutherford equation and using the PEST-III or resistive DCON magnetohydrodynamic codes, and the results are in reasonable agreement with each other. The calculated $Δ′$ is compared with an analytical formula of the tearing stability threshold $Δ_c^′$ to determine the condition for the mode onset, $Δ′$ > $Δ_c^′$. Finally, the local gradient effect of the plasma current and the pressure profiles on the tearing stability is investigated in terms of the difference between $Δ′$ and $Δ_c^′$ through a series of modeled equilibria with the parametric variations. A database of experimental mode instability onsets is compatible with the simulated regions of instability obtained from the parametric variations.

Linear stability analysis↗
DOE OSTI2023

PBE-GGA predicts the B8↔B2 phase boundary of FeO at Earth’s core conditions

FeO is a crucial component of the Earth’s core, and its thermodynamic properties are essential to developing more accurate core models. It is also a notorious correlated insulator in the NaCl-type (B1) phase at ambient conditions. It undergoes two polymorphic transitions at 300 K before it becomes metallic in the NiAs-type (B8) structure at ~100 GPa. Although its phase diagram is not fully mapped, it is well established that the B8 phase transforms to the CsCl-type (B2) phase at core pressures and temperatures. Here, we report a successful ab initio calculation of the B8↔B2 phase boundary in FeO at Earth’s core pressures. We show that fully anharmonic free energies computed with the Perdew–Burke–Ernzerhof-generalized gradient approximation coupled with thermal electronic excitations reproduce the experimental phase boundary within uncertainties at P > 255 GPa, including the largely negative Clapeyron slope of –52 MPa/K. This study validates the applicability of a standard density functional theory functional to FeO under Earth’s core conditions and demonstrates the theoretical framework that enables complex predictive studies of this region.

58 GEOSCIENCES↗
NASA NTRS1973

Reflection and interference of electromagnetic waves in inhomogeneous media

Solutions were obtained of the wave equation for a plane horizontally polarized electro-magnetic wave incident on a semi infinite two dimensional inhomogeneous medium. Two problems were considered: An inhomogeneous half space, and an inhomogeneous layer of arbitrary thickness. Solutions of the wave equation were obtained in terms of Hankel functions with complex arguments. Numerical calculations were made of the reflection coefficient R at the interface of the homogeneous medium. The startling results show that the reflection coefficient for a complex dielectric constant with gradient, can be less than that of the same medium with zero gradient.

Geiger, F. E.↗
NASA NTRS1975

Coherent ion acceleration at a density gradient

The paper describes a model, in which ions are accelerated against the direction of the electric field at a density gradient. A one-dimensional, inhomogeneous, steady-state is constructed, in which cold electrons stream relative to cold ions. The wavenumber and amplitude of the fluctuating potential are calculated, then used to compute the trajectories of representative ions trapped in the potential troughs. Ion acceleration to energies an order of magnitude greater than the applied potential is demonstrated.

Hsieh, S. L.↗
NASA NTRS1975

The vertical distribution of soluble gases in the troposphere

The thermodynamic properties of several water-soluble gases are reviewed to determine the likely effect of the atmospheric water cycle on their vertical profiles. It is found that gaseous HCl, HNO3, and HBr are sufficiently soluble in water to suggest that their vertical profiles in the troposphere have a similar shape to that of water vapor. It is predicted that HCl, HNO3, and HBr exhibit a steep negative gradient with altitude roughly equal to the altitude gradient of water vapor. Therefore, ground-based sources of inorganic chlorine, odd nitrogen, and inorganic bromine compounds are not likely to affect the stratosphere directly in the mean. Calculations also show that while SO2 and NH3 are significantly affected by the atmospheric water cycle, their abundances may not decrease with altitude as rapidly as does that of water vapor.

Stedman, D. H.↗
NASA NTRS1977

Investigation of supersonic turbulent boundary-layer separation on a compression ramp by an integral method

An investigation was made to determine the feasibility of using a boundary layer integral method to study the separation of a turbulent boundary layer on a two dimensional ramp at supersonic speeds. The numerical calculations were made for a free stream Mach number of 3, a Reynolds number of 10 million, and over a ramp angle range from 0 deg to 30 deg. For ramp angles where no flow separation was indicated, theoretical calculations were in reasonable agreement with experimental data except for a somewhat belated rise in pressure. For larger ramp angles, where separation was present, the investigation produced results that were not in agreement with experiment or with results calculated by time dependent Navier-Stokes methods. This apparently was true because no provision had been made for a proper shock boundary layer interaction where strong normal pressure gradients are induced within the boundary layer under the shock independent of surface curvature effects.

Patel, D. K.↗
NASA NTRS1978

Thermal stress analysis of ceramic structures with NASTRAN isoparametric solid elements

The performance of the NASTRAN level 16.0, twenty node, isoparametric bricks (CIHEX2) at thermal loading was studied. A free ceramic plate was modelled using twenty node bricks of varying thicknesses. The thermal loading for this problem was uniform over the surface with an extremely large gradient through the thickness. No mechanical loading was considered. Temperature-dependent mechanical properties were considered in this analysis. The NASTRAN results were compared to one dimensional stress distributions calculated by direct numerical integration.

Lamberson, S. E.↗
NASA NTRS1980

Velocity and concentration distribution in a Stefan diffusion tube

Numerical solutions for the diffusion- and gravity-driven flow in a cylindrical Stefan tube were obtained from the coupled diffusion and Navier-Stokes equations of Peclet numbers 0.3, 1 and 5. Distributions of binary component concentrations and velocities were calculated. The mass average velocity is parabolic in nature, except at high Peclet numbers. The solvent is not stagnant but recirculates, even in the absence of gravity. Radial concentration gradients develop which act convectively destabilizing. Consequences for the deduction of diffusion coefficients from Stefan tube experiments are discussed.

Markham, B. L.↗
NASA NTRS1989

Analysis of notched metal matrix composites under tensile loading

This paper presents techniques based on a three-dimensional finite-element analysis for the analysis of continuous fiber reinforced metal matrix composite. Examples are shown for specific metal matrix composites such as boron/aluminum and silicon-carbide/aluminum. Specimen stress-strain behavior and stress at first fiber failure were predicted for boron/aluminum laminates containing circular holes and crack-like slits. The predictions compared very well with test data for laminates containing 0 deg fibers and reasonably well for (+/-45)2s laminates. Mesh configuration was shown to have an effect on the calculation of stresses local to the notch. The presence of thin interface layers of matrix material had a significant influence on the slit tip stress state, causing sharper stress gradients near the notch. Interface layers reduced the slit-tip fibers stress in a (+/-45)2s silicon-carbide/aluminum laminate but increased them in a (0/90)2s laminate.

Bigelow, C. A.↗
NASA NTRS1992

Long Duration Exposure Facility (LDEF) attitude measurements of the interplanetary dust experiment

The Long Duration Exposure Facility (LDEF) Interplanetary Dust Experiment (IDE) was unique in providing a time history of impacts of micron-sized particles on six orthogonal faces of LDEF during the first year in orbit. The value of this time resolved data depended on and was enhanced by the proper operation of some basic LDEF systems. Thus, the value of the data is greatly enhanced when the location and orientation of LDEF is known for each time of impact. The location and velocity of LDEF as a function of time can be calculated from the 'two-line elements' published by GSFC during the first year of the LDEF mission. The attitude of LDEF was passively stabilized in a gravity-gradient mode and a magnetically anchored viscous damper was used to dissipate roll, pitch, and yaw motions. Finally, the IDE used a standard LDEF Experiment Power and Data System (EPDS) to collect and store data and also to provide a crystal derived clock pulse (1 count every 13.1072 seconds) for all IDE time measurements. All that remained for the IDE was to provide a system to calibrate the clock, eliminating accumulative errors, and also verify the attitude of LDEF. The IDE used solar cells on six orthogonal faces to observe the LDEF sunrise and provide data about the LDEF attitude. The data was recorded by the EPDS about 10 times per day for the first 345 days of the LDEF mission. This data consist of the number of IDE counts since the last LDEF sunrise and the status of the six solar cells (light or dark) at the time of the last IDE count. The EPDS determined the time that data was recorded and includes, with each record, the master EPDS clock counter (1 count every 1.6384 seconds) that provided the range and resolution for time measurements. The IDE solar cells provided data for an excellent clock calibration, meeting their primary purpose, and the time resolved LDEF attitude measurements that can be gleaned from this data are presented.

Kassel, Philip C., Jr.↗
NASA NTRS2004

Aerodynamic Optimization of Rocket Control Surface Geometry Using Cartesian Methods and CAD Geometry

Aerodynamic design is an iterative process involving geometry manipulation and complex computational analysis subject to physical constraints and aerodynamic objectives. A design cycle consists of first establishing the performance of a baseline design, which is usually created with low-fidelity engineering tools, and then progressively optimizing the design to maximize its performance. Optimization techniques have evolved from relying exclusively on designer intuition and insight in traditional trial and error methods, to sophisticated local and global search methods. Recent attempts at automating the search through a large design space with formal optimization methods include both database driven and direct evaluation schemes. Databases are being used in conjunction with surrogate and neural network models as a basis on which to run optimization algorithms. Optimization algorithms are also being driven by the direct evaluation of objectives and constraints using high-fidelity simulations. Surrogate methods use data points obtained from simulations, and possibly gradients evaluated at the data points, to create mathematical approximations of a database. Neural network models work in a similar fashion, using a number of high-fidelity database calculations as training iterations to create a database model. Optimal designs are obtained by coupling an optimization algorithm to the database model. Evaluation of the current best design then gives either a new local optima and/or increases the fidelity of the approximation model for the next iteration. Surrogate methods have also been developed that iterate on the selection of data points to decrease the uncertainty of the approximation model prior to searching for an optimal design. The database approximation models for each of these cases, however, become computationally expensive with increase in dimensionality. Thus the method of using optimization algorithms to search a database model becomes problematic as the number of design variables is increased.

Nelson, Andrea↗
DOE OSTI2014

A two-dimensional (azimuthal-axial) particle-in-cell model of a Hall thruster

We have developed a two-dimensional Particle-In-Cell model in the azimuthal and axial directions of the Hall thruster. A scaling method that consists to work at a lower plasma density to overcome constraints on time-step and grid-spacing is used. Calculations are able to reproduce the breathing mode due to a periodic depletion of neutral atoms without the introduction of a supplementary anomalous mechanism, as in fluid and hybrid models. Results show that during the increase of the discharge current, an electron-cyclotron drift instability (frequency in the range of MHz and wave number on the order of 3000 rad s{sup −1}) is formed in the region of the negative gradient of magnetic field. During the current decrease, an axial electric wave propagates from the channel toward the exhaust (whose frequency is on the order of 400 kHz) leading to a broadening of the ion energy distribution function. A discussion about the influence of the scaling method on the calculation results is also proposed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

End-wall boundary layer measurements in a two-stage fan

Detailed flow measurements made in the casing boundary layer of a two-stage transonic fan are summarized. These measurements were taken at stations upstream of the fan, between all blade rows, and downstream of the last blade row. At the design tip speed of 429 m/sec the fan achieved a peak efficiency of 0.846 at a pressure ratio of 2.471. The boundary layer data were obtained at three weight flows at the design speed: one near choke flow, one near peak efficiency, and one near stall. Conventional boundary layer parameters were calculated from the data measured at each measuring station for each of the three flows. A classical two dimensional casing boundary layer was measured at the fan inlet and extended inward to approximately 15 percent of span. A highly three dimensional boundary layer was measured at the exit of each blade row and extended inward to approximately 10 percent of span. The steep radial gradient of axial velocity noted at the exit of the rotors was reduced substantially as the flow passed through the stators. This reduced gradient is attributed to flow mixing. The amount of flow mixing was reflected in the radial redistribution of total temperature as the flow passed through the stators. The data also show overturning of the tip flow at the stator exits that is consistent with the expected effect of the secondary flow field. The blockage factors calculated from the measured data show an increase in blockage across the rotors and a decrease across the stators.

Ball, C. L.↗
DOE OSTI2014

Experimental and theoretical demonstration of acoustic Bloch oscillations in porous silicon structures

We report the theoretical calculations and the experimental demonstration of acoustic Bloch oscillations and Wannier-Stark ladders in linear tilted multilayer structures based on porous silicon. The considered structures consist of layers with constant porosity alternated by layers with a linear gradient in the parameter η=1/v{sub L}{sup 2} along the growth direction in order to tilt the acoustic band gap. The purpose of this gradient is to mimic the tilted electronic miniband structure of a superlattice semiconductor under an external electric field. In this way, acoustic Wannier-Stark ladders of equidistant modes are formed and they were experimentally confirmed in the transmission spectrum around 1.2 GHz. Their frequency separation defines the period of the acoustic Bloch oscillations. We fabricated three different structures with the same thicknesses but different values in the η parameter to observe the effect on the period of the Bloch oscillations. We measured the acoustic transmission spectra in the frequency domain, and by using the Fourier transform, we obtained the transmission in the time domain. The transmission spectra of the fabricated samples show acoustic Bloch oscillations with periods of 27, 24, and 19 ns. The experimental results are in good agreement with the transfer matrix calculations. The observed phenomenon is the acoustic counterpart of the well known electronic Bloch oscillations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Design of an additively manufactured functionally graded material of 316 stainless steel and Ti-6Al-4V with Ni-20Cr, Cr, and V intermediate compositions

This study presents a method for designing a computationally informed gradient pathway to fabricate a functionally graded material (FGM) with terminal alloys of 316 stainless steel (SS316) and Ti-6Al-4V via directed energy deposition additive manufacturing with powder feedstock. In this work, the grading is accomplished through the introduction of intermediate elements and alloys (Ni-20Cr, Cr, and V) to avoid the brittle Fe-Ti intermetallic phases that form in the direct liquid phase joining of Ti-alloys and stainless steels. Using a combination of equilibrium calculations and Scheil-Gulliver simulations, a compositional pathway was designed to avoid deleterious phases. FGM samples were fabricated and experimentally characterized to determine the viability of the pathway. A change in phases from fcc to bcc was predicted to occur within the Ni-20Cr/Cr gradient region, and this was validated through experimental characterization. No detrimental phases (intermetallic, Laves, or σ phases) formed along the gradient path, demonstrating a successful computationally-informed design and fabrication of an FGM from SS316 to Ti-6Al-4V.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Formation and Evolution Simulations of Saturn, Including Composition Gradients and Helium Immiscibility

In this paper, the formation of Saturn is modeled by detailed numerical simulations according to the core-nucleated accretion scenario. Previous models are enhanced to include the dissolution of accreting planetesimals, composed of water ice, rock, and iron, in the gaseous envelope of the planet, leading to a nonuniform composition with depth. The immiscibility of helium in metallic hydrogen layers is also considered. The calculations start at a mass of 0.5 Earth masses and are extended to the present day. At 4.57 Gyr, the model, proceeding outward, has the following structure: (i) a central core composed of 100% heavy elements and molecules, (ii) a region with a decreasing heavy-element mass fraction, down to a value of 0.1, (iii) a layer of uniform composition with the helium mass fraction Y enhanced over the primordial value, (iv) a helium rain region with a gradient in Y, (v) an outer convective, adiabatic region with uniform composition in which Y is reduced from the primordial value, and (vi) the very outer layers where cloud condensation of the heavy elements occurs. Models of the distribution of heavy elements as a function of radius are compared with those derived to fit the observations of the Cassini mission, with rough qualitative agreement. The helium mass fraction in Saturn’s outer layers is estimated to be around 20%. Models are found that provide good agreement with Saturn’s intrinsic luminosity and radius.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

The relationship between local liquid density and force applied on a tip of atomic force microscope: A theoretical analysis for simple liquids

The density of a liquid is not uniform when placed on a solid. The structured liquid pushes or pulls a probe employed in atomic force microscopy, as demonstrated in a number of experimental studies. In the present study, the relation between the force on a probe and the local density of a liquid is derived based on the statistical mechanics of simple liquids. When the probe is identical to a solvent molecule, the strength of the force is shown to be proportional to the vertical gradient of ln(ρ{sub DS}) with the local liquid's density on a solid surface being ρ{sub DS}. The intrinsic liquid's density on a solid is numerically calculated and compared with the density reconstructed from the force on a probe that is identical or not identical to the solvent molecule.

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

Radiative precursors driven by converging blast waves in noble gases

A detailed study of the radiative precursor that develops ahead of converging blast waves in gas-filled cylindrical liner z-pinch experiments is presented. The experiment is capable of magnetically driving 20 km s{sup −1} blast waves through gases of densities of the order 10{sup −5} g cm{sup −3} (see Burdiak et al. [High Energy Density Phys. 9(1), 52–62 (2013)] for a thorough description). Data were collected for Ne, Ar, and Xe gas-fills. The geometry of the setup allows a determination of the plasma parameters both in the precursor and across the shock, along a nominally uniform line of sight that is perpendicular to the propagation of the shock waves. Radiation from the shock was able to excite NeI, ArII, and XeII/XeIII precursor spectral features. It is shown that the combination of interferometry and optical spectroscopy data is inconsistent with upstream plasmas being in LTE. Specifically, electron density gradients do not correspond to any apparent temperature change in the emission spectra. Experimental data are compared to 1D radiation hydrodynamics HELIOS-CR simulations and to PrismSPECT atomic physics calculations to assist in a physical interpretation of the observations. We show that upstream plasma is likely in the process of being radiatively heated and that the emission from a small percentage of ionised atoms within a cool background plasma dominates the emission spectra. Experiments were carried out on the MAGPIE and COBRA pulsed-power facilities at Imperial College London and Cornell University, respectively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗