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,279 records · Page 71

Heat Shield Cavity Parametric Experimental Aeroheating for a Proposed Mars Smart Lander Aeroshell

The proposed Mars Smart Lander is to be attached through its aeroshell to the main spacecraft bus, thereby producing cavities in the heat shield. To study the effects these cavities will have on the heating levels experienced by the heat shield, an experimental aeroheating investigation was performed at the NASA Langley Research Center in the 20-Inch Mach 6 Air Tunnel. The effects of Reynolds number, angle-of-attack, and cavity size and location on aero-heating levels and distributions were determined and are presented. To aid the discussion on the effects of the cavities, laminar, thin-layer Navier-Stokes flow field solutions were post-processed to calculate relevant boundary layer properties such as boundary layer height and momentum thickness, edge Mach number, and streamwise pressure gradient. It was found that the effect of the cavities varies with angle-of-attack, freestream Reynolds number, and cavity size and location. The presence of a cavity raised the downstream heating rates by as much as 325% as a result of boundary layer transition.

Liechty, Derek S.↗

A study on the criticality safety for the direct disposal of used nuclear fuel in Japan: Application of burnup credit to the criticality safety evaluation for the disposal canister - 14049

Japan has recently started the technical development of the direct disposal of used nuclear fuel (UNF) to prepare various disposal options. Unlike the disposal of radioactive waste from UNF reprocessing, criticality safety is important here because of the presence of a certain amount of fissile nuclides in UNF. This paper provides an outline of the research to be addressed in this field and the relevant studies in Japan. It especially includes the first result of a criticality safety evaluation of a disposal canister model adopting burnup credit. In this evaluation, it was assumed that a steel disposal canister loaded with fuel assemblies was surrounded by a clay buffer material, and the fuel assemblies and the disposal canister were intact. Conditions representative of discharged pressurized-water reactor (PWR) fuel were applied in this study - initial U-235 enrichment of 4.5 wt% and burnup of 45 GWd/t. The isotopic composition of the burnt fuel was evaluated using SWAT3.1 code. In the depletion calculation, conservative parameters were applied within the typical PWR operation range. The effective neutron multiplication factor (k{sub eff}) was then calculated using a continuous-energy Monte Carlo code MVP-2.0, adopting the JENDL-4.0 nuclear data library. Uncertainty in k{sub eff} due to depletion calculation errors was evaluated using the correction factors of isotopic composition, which were obtained via multiple comparisons between the calculated and measured isotopic concentrations of irradiated fuels. In addition, the effects of axial burnup profile and horizontal burnup gradient on k{sub eff} were evaluated using the open burnup distribution database. The k{sub eff} including the uncertainties and conservatism was found to be below 0.95 for both actinide-only and actinide and fission product considerations when the fuel assemblies and the disposal canister were assumed to be intact. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Analysis and calculation of macrosegregation in a casting ingot. MPS solidification model. Volume 1: Formulation and analysis

The physical and numerical formulation of a model for the horizontal solidification of a binary alloy is described. It can be applied in an ingot. The major purpose of the model is to calculate macrosegregation in a casting ingot which results from flow of interdendritic liquid during solidification. The flow, driven by solidification contractions and by gravity acting on density gradients in the interdendritic liquid, was modeled as flow through a porous medium. The symbols used are defined. The physical formulation of the problem leading to a set of equations which can be used to obtain: (1) the pressure field; (2) the velocity field: (3) mass flow and (4) solute flow in the solid plus liquid zone during solidification is presented. With these established, the model calculates macrosegregation after solidification is complete. The numerical techniques used to obtain solution on a computational grid are presented. Results, evaluation of the results, and recommendations for future development of the model are given. The macrosegregation and flow field predictions for tin-lead, aluminum-copper, and tin-bismuth alloys are included as well as comparisons of some of the predictions with published predictions or with empirical data.

Maples, A. L.↗

Scrape-off layer flows with pressure gradient scale length ~ ρ p

A heuristic model for the plasma scrape-off width balances magnetic drifts against parallel loss at c s /2, resulting in a SOL width ~ ρ p . T sep is calculated from Spitzer–Härm parallel thermal conduction. This results in a prediction for the power scrape-off width in quantitative agreement both in magnitude and scaling with recent experimental data. To achieve the ~ c s /2 flow assumed in this model and measured experimentally sets requirements on the ratio of upstream to total SOL particle sources, relative to the square-root of the ratio of target to upstream temperature. The Pfisch-Schlüter model for equilibrium flows has been modified to allow near-sonic flows, appropriate for gradient scale lengths of order ρ p , resulting in a new quadrupole radial flow pattern. The strong parallel flows and plasma charging implied by this model suggest a mechanism for H-mode transition, consistent with many observations

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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 a station upstream of the fan, between all blade rows, and downstream of the last row. Conventional boundary layer parameters were calculated from the measured data. 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 blockage factors calculated from the measured data show an increase in blockage across the rotors and a decrease across the stators. For this fan the calculated blockages for the second stage were essentially the same as those for the first stage.

Ball, C. L.↗

The Uncertainty of Biomass Estimates from LiDAR and SAR Across a Boreal Forest Structure Gradient

In this study, weexamined the uncertainty of aboveground live biomass (AGB) estimates based on light detection and ranging (LiDAR) and synthetic aperture radar (SAR) measurements distributed across a low-biomass vegetation structure gradient fromforest to non-forest in boreal-like ecosystems. The conifer-dominant structure gradient was compiled from ground data amassed from multiple field expeditions in central Maine (USA), Aurskog (Norway), and across central Siberia (Russia). Single variable empirical models were built tomodel AGB fromremote sensingmetrics. Using thesemodels, we calculated a rootmean square error (RMSE) and a 95% confidence interval (CI) of the RMSE fromthe difference between the remote sensing AGB predictions and the ground reference AGB estimates within AGB intervals across a 0-100 Mg ha(exp.−1) boreal forest structure gradient. The results show that the error in AGB predictions (RMSE) and the error uncertainty (the CI) from LiDAR and SAR change across a forest gradient. The errors of airborne LiDAR and SAR metrics and spaceborne LiDAR platforms show a general trend of reduced relative errors as AGBmagnitudes increase, particularly from0 to 60 Mg ha (exp.−1). Empirical models relating spaceborne metrics to AGB and estimates of spaceborne LiDAR error uncertainty demonstrate the difficulty of characterizing differences in AGB at the site-level with current spaceborne sensors, particularly below 80 Mg ha (exp.−1) with less than 50-100% error.

Taiga↗

Piezo-optic tensor of crystals from quantum-mechanical calculations

An automated computational strategy is devised for the ab initio determination of the full fourth-rank piezo-optic tensor of crystals belonging to any space group of symmetry. Elastic stiffness and compliance constants are obtained as numerical first derivatives of analytical energy gradients with respect to the strain and photo-elastic constants as numerical derivatives of analytical dielectric tensor components, which are in turn computed through a Coupled-Perturbed-Hartree-Fock/Kohn-Sham approach, with respect to the strain. Both point and translation symmetries are exploited at all steps of the calculation, within the framework of periodic boundary conditions. The scheme is applied to the determination of the full set of ten symmetry-independent piezo-optic constants of calcium tungstate CaWO{sub 4}, which have recently been experimentally reconstructed. Present calculations unambiguously determine the absolute sign (positive) of the π{sub 61} constant, confirm the reliability of 6 out of 10 experimentally determined constants and provide new, more accurate values for the remaining 4 constants.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Formation Mechanism and Nucleation Effect of Ti{sub 2}O{sub 3}–TiN Complex Nucleus at Solidification Front of 18Cr Ferritic Stainless Steel

Formation mechanism of Ti{sub 2}O{sub 3}–TiN complex nucleus and its nucleation effect on δ-Fe at solidification front of Fe-18Cr ferritic stainless steel was investigated through experimental and theoretical methods. According to the equilibrium calculation, Ti{sub 2}O{sub 3} will crystallize in the molten steel before solidification and TiN will not form until the solid phases appear. However, complex nucleus and equiaxed grains exist in the edge of Ti-bearing ingot, where the solidification velocity and temperature gradient are high. Therefore, complex nucleus is indicated to be formed at solidification front and a non-equilibrium solidification model for multi-component alloy was established to investigate the dendrite tip interface phenomenon. The calculation results exhibit an element enrichment at solidification front and explain well with the formation of complex nucleus before solidification. It is also found that there is approximate lattice matching between TiN, and Ti{sub 2}O{sub 3}, δ-Fe on the basis of disregistry theory. These results are consistent with the observed complex nucleus and the increased equiaxed zone ratio (EZR) of Ti-bearing ferritic stainless steel.

36 MATERIALS SCIENCE↗

Efficient calculation of optimum design sensitivity

Two methods for obtaining the sensitivity of an optimized design to changes in problem parameters are presented. The first-order method uses gradient information to estimate the required sensitivity. The second-order method uses second derivatives of the design objective and constraint functions to provide a quadratic approximation to the new design which results from changing the specified parameter. It is shown that the sensitivity of the optimized design may be discontinuous with respect to the new parameter. This discontinuity is accounted for in the present methods. These methods are compared to a previous approach based on the Kuhn-Tucker necessary conditions for optimality. It is shown that the present methods provide the appropriate sensitivity information in an efficient manner. The methods are demonstrated by a structural synthesis example.

Vanderplaats, G. N.↗

Study of effect of PWR cold leg temperature gradient on reactor core condition

Effects of temperature and flow gradients in Westinghouse designed three-loop Pressurized Water Reactor (PWR) cold legs, the piping between the main coolant pump and the reactor vessel, were evaluated using Computational Fluid Dynamics (CFD) code STAR CCM+ and coupled neutronic and thermal-hydraulic (T/H) code system VERA. In the parametric study, several symmetric and asymmetric temperature gradients that were significantly larger than those observed from plant measurements were applied to the cold leg inlets for comparison with the base case without any temperature gradient. A CFD model using the STAR-CCM+ code was developed for a portion of the RCS region between the Reactor Coolant Pump (RCP) and the core inlet based on previously validated modeling approach. The CFD simulation results were processed for the temperature and flow rate distributions at the core inlet as input to the VERA calculations. The VERA code system consists of COBRA-TF (CTF) for thermal-hydraulics, MPACT for reactor physics and neutron transport, and ORIGEN for isotopic depletion. The VERA model was for depletion calculations of a high-burnup loading pattern with the reactor core in pin-by-pin and subchannel resolution. The results of the study indicate that the postulated temperature gradients within the PWR cold legs do not result in any significant changes in the core inlet temperature distributions and the core power distributions during the reactor operation. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Nonisothermal Pluto atmosphere models

The present thermal profile calculation for a Pluto atmosphere model characterized by a high number fraction of CH4 molecules encompasses atmospheric heating by solar UV flux absorption and conductive transport cooling to the surface of Pluto. The stellar occultation curve predicted for an atmosphere of several-microbar surface pressures (which entail the existence of a substantial temperature gradient close to the surface) agrees with observations and implies that the normal and tangential optical depth of the atmosphere is almost negligible. The minimum period for atmospheric methane depletion is calculated to be 30 years.

Hubbard, W. B.↗

Mathematical simulation of heat and mass transfer in a cylindrical channel versus the tangential momentum accommodation coefficient

The process of heat and mass transfer in a long cylindrical channel has been considered in terms of the mirror–diffuse model of the Maxwell boundary condition. The Hazen–Williams equation is used as a basic equation of the process kinetics. A constant temperature gradient is maintained in the channel. The heat and mass fluxes through the cross section of the channel versus the tangential momentum accommodation coefficient have been calculated in a wide range of the Knudsen number. The heat flux profiles have been constructed. A comparison with relevant published data has been carried out.

42 ENGINEERING↗

Enhancing efficiency and power of quantum-dots resonant tunneling thermoelectrics in three-terminal geometry by cooperative effects

We propose a scheme of multilayer thermoelectric engine where one electric current is coupled to two temperature gradients in three-terminal geometry. This is realized by resonant tunneling through quantum dots embedded in two thermal and electrical resisting polymer matrix layers between highly conducting semiconductor layers. There are two thermoelectric effects, one of which is pertaining to inelastic transport processes (if energies of quantum dots in the two layers are different), while the other exists also for elastic transport processes. These two correspond to the transverse and longitudinal thermoelectric effects, respectively, and are associated with different temperature gradients. We show that cooperation between the two thermoelectric effects leads to markedly improved figure of merit and power factor, which is confirmed by numerical calculation using material parameters. Such enhancement is robust against phonon heat conduction and energy level broadening. Therefore, we demonstrated cooperative effect as an additional way to effectively improve performance of thermoelectrics in three-terminal geometry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

An improved guess for the variational calculation of charge-transfer excitations in large systems

Ab initio quantum-chemical methods that perform well for computing the electronic ground state are not straightforwardly transferable to electronically excited states, particularly in large molecular systems. Wave function theory offers high accuracy, but is often prohibitively expensive. Methods based on time-dependent density functional theory (TD-DFT) are crucially sensitive to the chosen exchange-correlation functional (XCF) parameterization, and system-specific tuning protocols were therefore proposed to address the method's robustness. Methods based on the variational relaxation of the excited-state electron density showcased promising results for the calculation of charge-transfer excitations, but the complex shape of the electronic hypersurface makes convergence to a specific excited state much more difficult than for the ground state when standard variational techniques are applied. We address the latter aspect by providing suitable initial guesses, which we obtain by two separate constrained algorithms. Combined with the squared-gradient minimization algorithm for all-electrons relaxation in a freeze-and-release scheme (FRZ-SGM), we demonstrate that orbital-optimized density functional theory (OO-DFT) calculations can reliably converge to the charge-transfer states of interest even for large molecular systems. We test the FRZ-SGM method on a phenothiazine-anthraquinone CT excitation in a supramolecular Pd(II) coordination cage complex as a function of the cage conformation. This compound has been studied experimentally prior to our work. We compare this freeze-and-release scheme to two XCF reparameterizations, which were recently proposed as low-cost TD-DFT-based alternatives to variational methods. Two dye-semiconductor complexes, which were previously investigated in the context of photovoltaic applications, serve as a second example to investigate the convergence and stability of the FRZ-SGM approach. Our results demonstrate that FRZ-SGM provides reliable convergence for charge-transfer excited states and avoids variational collapse to lower-lying electronic states, whereas time-dependent DFT calculations with an adequate tuning procedure for the range-separation parameter provide a computationally efficient initial estimate of the corresponding energies, with a computational cost comparable to that of configuration-interaction singles (CIS) calculations.

Bogo, Nicola↗

The dynamical evolution of the protosolar nebula

Evolutionary models for protostellar nebulae are calculated under the hypothesis that the only source for the turbulent viscosity is thermal convection. The viscous stress is approximated by an 'alpha' model, and the constant alpha is calculated in terms of the properties of turbulent thermal convection. A relatively sensitive dependence of the Rosseland mean opacity on temperature is needed for the vertical temperature gradient of the nebula to become convectively unstable. However, this requires that the vertical optical depth in the nebula is relatively sensitive to the total surface density, and as the disk is depleted of matter by protostellar accretion, the Rosseland mean optical depth drops and the disk must become convectively unstable. This limits the amount of mass that the nebula can accretionally process before the convective turbulence ceases. The resulting disk evolutionary properties are calculated and comparisons with the solar system are made.

Ruden, Steven P.↗

Dynamics of intrinsic axial flows in unsheared, uniform magnetic fields

A simple model for the generation and amplification of intrinsic axial flow in a linear device, controlled shear decorrelation experiment, is proposed. This model proposes and builds upon a novel dynamical symmetry breaking mechanism, using a simple theory of drift wave turbulence in the presence of axial flow shear. This mechanism does not require complex magnetic field structure, such as shear, and thus is also applicable to intrinsic rotation generation in tokamaks at weak or zero magnetic shear, as well as to linear devices. This mechanism is essentially the self-amplification of the mean axial flow profile, i.e., a modulational instability. Hence, the flow development is a form of negative viscosity phenomenon. Unlike conventional mechanisms where the residual stress produces an intrinsic torque, in this dynamical symmetry breaking scheme, the residual stress induces a negative increment to the ambient turbulent viscosity. The axial flow shear is then amplified by this negative viscosity increment. The resulting mean axial flow profile is calculated and discussed by analogy with the problem of turbulent pipe flow. For tokamaks, the negative viscosity is not needed to generate intrinsic rotation. However, toroidal rotation profile gradient is enhanced by the negative increment in turbulent viscosity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multistate hybrid time-dependent density functional theory with surface hopping accurately captures ultrafast thymine photodeactivation

We report an efficient analytical implementation of first-order nonadiabatic derivative couplings between arbitrary Born–Oppenheimer states in the hybrid time-dependent density functional theory (TDDFT) framework using atom-centered basis functions. Our scheme is based on quadratic response theory and includes orbital relaxation terms neglected in previous approaches. Simultaneous computation of multiple derivative couplings and energy gradients enables efficient multistate nonadiabatic molecular dynamics simulations in conjunction with Tully's fewest switches surface hopping (SH) method. We benchmark the thus obtained multistate TDDFT-SH scheme by simulating ultrafast decay of UV-photoexcited thymine, for which accurate gas-phase data from ultrafast spectroscopy experiments are available. The calculations predict a fast 153 fs decay from the bright S 2 to the dark S 1 excited state, followed by a much slower 14 ps S 1 deactivation to the ground state; statistical uncertainties were estimated using bootstrap sampling. These results agree well with the experimentally observed time constants of 100–200 fs and 5–7 ps, respectively, unlike previous multiconfigurational self-consistent field and second-order algebraic diagrammatic construction calculations. Furthermore, our results support the S 1 -trapping hypothesis. For thymine, the computational cost of a single TDDFT-SH time-step including the lowest 3 states, all couplings and gradients, is ~5 times larger than the cost of a single Born–Oppenheimer dynamics time step for the ground state in our implementation. Thus, ps nonadiabatic dynamics simulations using multistate hybrid TDDFT-SH for systems with up to ~100 atoms are possible without drastic approximations on single workstation nodes. Our implementation will be made available through Turbomole.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Boundary-layer linear stability theory

Most fluid flows are turbulent rather than laminar and the reason for this was studied. One of the earliest explanations was that laminar flow is unstable, and the linear instability theory was first developed to explore this possibility. A series of early papers by Rayleigh produced many notable results concerning the instability of inviscid flows, such as the discovery of inflectional instability. Viscosity was commonly thought to act only to stabilize the flow, and flows with convex velocity profiles appeared to be stable. The investigations that led to a viscous theory of boundary layer instability was reported. The earliest application of linear stability theory to transition prediction calculated the amplitude ratio of the most amplified frequency as a function of Reynolds number for a Blasius boundary layer, and found that this quantity had values between five and nine at the observed Ret. The experiment of Schubauer and Skramstad (1947) completely reversed the prevailing option and fully vindicated the Gottingen proponents of the theory. This experiment demonstrated the existence of instability waves in a boundary layer, their connection with transition, and the quantitative description of their behavior by the theory of Tollmien and Schlichting. It is generally accepted that flow parameters such as pressure gradient, suction and heat transfer qualitatively affect transition in the manner predicted by the linear theory, and in particular that a flow predicted to be stable by the theory should remain laminar. The linear theory, in the form of the e9, or N-factor is today in routine use in engineering studies of laminar flow. The stability theory to boundary layers with pressure gradients and suction was applied. The only large body of numerical results for exact boundary layer solutions before the advent of the computer age by calculating the stability characteristics of the Falkner-Skan family of velocity profiles are given. When the digital computer reached a stage of development which permit the direct solution of the primary differential equations, numerical results were obtained from the linear theory during the next 10 years for many different boundary layer flows: three dimensional boundary layers; free convention boundary layers; compressible boundary layers; boundary layers on compliant walls; a recomputation of Falkner-Skan flows; unsteady boundary layers; and heated wall boundary layers.

Leslie M Mack↗