Search NASA⌕ Search

SEARCH · Search NASA

Results for “RELAXATION METHOD”

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 253 records · Page 14

The Electron Thermal Conductivity of Pu and Zr Substituted Gamma-Uranium

Uranium alloys are attractive recycled nuclear fuels because of their high thermal conductivity (k) and fissile density; however, the effects of alloying elements on k remain unclear. Here, the electron thermal conductivity (k_e) of U-Pu-Zr compositions are calculated using density functional theory. The electronic structure is evaluated to understand the effects of plutonium (Pu) and zirconium (Zr) substitution on the k_e of ?-U. Alloys of up to 37.5 at. % Pu and 37.5 at. % Zr are examined. Two methods are applied to calculate k_e; we find that the accuracy of each method depends on the electronic and mass similarities between the solute and solvent atoms. Specifically, when the solute atom is similar in electronic structure and mass, the method that applies the electron relaxation time of ?-U is best, while if the elements are dissimilar, a mixed method that mixes several parameters associated with k_e from each element in the alloy is best. The introduction of all alloying elements decreases k_e; however, in binary compounds, Pu and Zr have different effects. Pu generally flattens the electronic bands but compensates for this deleterious effect by increasing electron density near the Fermi level. Zr flattens the electronic bands more severely without adding electron density near the Fermi level. Therefore, Zr decreases the k_e more than Pu in binary compounds. In ternary compounds, the difference between Pu and Zr is minimal due to the phononic change from the large mass change of Zr substitution, even at 12.5 at. %. Thus, we predict that higher loadings of Pu, and potentially other actinides, can be added to U-Pu-Zr compositions for faster recycling of spent fuel with without sacrificing k. We also note that these k_e calculation methods can be applied to non-fuel alloys that require k_e predictions, such as cladding, heat exchanger, and structural materials.

36 MATERIALS SCIENCE↗

Accurate numerical simulations of open quantum systems using spectral tensor trains

Decoherence between qubits is a major bottleneck in quantum computations. Decoherence results from intrinsic quantum and thermal fluctuations as well as noise in the external fields that perform the measurement and preparation processes. With prescribed colored noise spectra for intrinsic and extrinsic noise, we present a numerical method, Quantum Accelerated Stochastic Propagator Evaluation (Q-ASPEN), to solve the time-dependent noise-averaged reduced density matrix in the presence of intrinsic and extrinsic noise. Q-ASPEN is arbitrarily accurate and can be applied to provide estimates for the resources needed to error-correct quantum computations. We employ spectral tensor trains, which combine the advantages of tensor networks and pseudospectral methods, as a variational ansatz to the quantum relaxation problem and optimize the ansatz using methods typically used to train neural networks. Here, the spectral tensor trains in Q-ASPEN make accurate calculations with tens of quantum levels feasible. We present benchmarks for Q-ASPEN on the spin-boson model in the presence of intrinsic noise and on a quantum chain of up to 32 sites in the presence of extrinsic noise. In our benchmark, the memory cost of Q-ASPEN scales as a low-order polynomial in the size of the system once the number of system states surpasses the number of basis functions used in the spectral expansion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Growing High-Quality InAs Quantum Dots for Infrared Lasers

An improved method of growing high-quality InAs quantum dots embedded in lattice-matched InGaAs quantum wells on InP substrates has been developed. InAs/InGaAs/InP quantum dot semiconductor lasers fabricated by this method are capable of operating at room temperature at wavelengths greater than or equal to 1.8 mm. Previously, InAs quantum dot lasers based on InP substrates have been reported only at low temperature of 77 K at a wavelength of 1.9 micrometers. In the present method, as in the prior method, one utilizes metalorganic vapor phase epitaxy to grow the aforementioned semiconductor structures. The development of the present method was prompted in part by the observation that when InAs quantum dots are deposited on an InGaAs layer, some of the InAs in the InGaAs layer becomes segregated from the layer and contributes to the formation of the InAs quantum dots. As a result, the quantum dots become highly nonuniform; some even exceed a critical thickness, beyond which they relax. In the present method, one covers the InGaAs layer with a thin layer of GaAs before depositing the InAs quantum dots. The purpose and effect of this thin GaAs layer is to suppress the segregation of InAs from the InGaAs layer, thereby enabling the InAs quantum dots to become nearly uniform (see figure). Devices fabricated by this method have shown near-room-temperature performance.

Qiu, Yueming↗

Data-driven equation-free dynamics applied to many-protein complexes: The microtubule tip relaxation

Microtubules (MTs) constitute the largest components of the eukaryotic cytoskeleton and play crucial roles in various cellular processes, including mitosis and intracellular transport. The property allowing MTs to cater to such diverse roles is attributed to dynamic instability, which is coupled to the hydrolysis of GTP (guanosine-5'-triphosphate) to GDP (guanosine-5'-diphosphate) within the β-tubulin monomers. Understanding the equilibrium dynamics and the structural features of both GDP- and GTP-complexed MT tips, especially at an all-atom level, remains challenging for both experimental and computational methods because of their dynamic nature and the prohibitive computational demands of simulating large, many-protein systems. This study employs the “equation-free” multiscale computational method to accelerate the relaxation of all-atom simulations of MT tips toward their putative equilibrium conformation. Using large MT lattice systems (14 protofilaments × 8 heterodimers) comprising ~21-38 million atoms, we applied this multiscale approach to leapfrog through time and nearly double the computational efficiency in realizing relaxed all-atom conformations of GDP- and GTP-complexed MT tips. Commencing from an initial 4 μs unbiased all-atom simulation, we interleave coarse projective “equation-free” jumps with short bursts of all-atom molecular dynamics simulation to realize an additional effective simulation time of 1.875 μs. Our 5.875 μs of effective simulation trajectories for each system expose the subtle yet essential differences in the structures of MT tips as a function of whether β-tubulin monomer is complexed with GDP or GTP, as well as the lateral interactions within the MT tip, offering a refined understanding of features underlying MT dynamic instability. Furthermore, the approach presents a robust and generalizable framework for future explorations of large biomolecular systems at atomic resolution.

Wu, Jiangbo [University of Chicago, IL (United Sta↗

Hole Relaxation Bottlenecks in CdSe/CdTe/CdSe Lateral Heterostructures Lead to Bicolor Emission

Concentric lateral CdSe/CdTe/CdSe heterostructures show bicolor photoluminescence from both a red charge transfer band of the CdSe/CdTe interface and a green fluorescence from CdSe. Further, this work uses visible and near-infrared transient spectroscopy measurements to demonstrate that the deviation from Kasha's rule arises from a hole relaxation bottleneck from CdSe to CdTe. Hole transfer can take up to 1 ns, which permits radiative relaxation of excitons remaining in CdSe. Simulations indicate that the hole relaxation bottleneck arises due to the sparse density of states and poor spatial overlap of hole states at energies near the CdSe band edge. The divergent kinetics of transfer for band edge and hot holes is exploited to vary the ratio of green and red photoluminescence with excitation wavelength, providing another knob to control emission color. These findings support the use of lateral heterojunctions as a method for slowing carrier relaxation in two-dimensional materials.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

An efficient numerical method for solving the incompressible Navier-Stokes equations

This paper describes an efficient numerical method for solving the steady incompressible Navier-Stokes equations. The method is a fully implicit method based on the generalized Galerkin method, and the resulting system of equations is solved in a sweeping mode by iterative line relaxation. Results of the present method are compared with published results for separating and reattaching flows, and parametric studies showing the effects of step size, boundary locations, and Reynolds number are presented. The present method is substantially faster than previously published methods; typical run times range from 10 sec to 1 min of 7600 CPU time; and, based on results obtained to date, it is stable at any Reynolds number.

Murphy, J. D.↗

Application of a multi-level grid method to transonic flow calculations

A multi-level grid method was studied as a possible means of accelerating convergence in relaxation calculations for transonic flows. The method employs a hierarchy of grids, ranging from very coarse to fine. The coarser grids are used to diminish the magnitude of the smooth part of the residuals. The method was applied to the solution of the transonic small disturbance equation for the velocity potential in conservation form. Nonlifting transonic flow past a parabolic arc airfoil is studied with meshes of both constant and variable step size.

South, J. C., Jr.↗

Distributed Relaxation for Conservative Discretizations

A multigrid method is defined as having textbook multigrid efficiency (TME) if the solutions to the governing system of equations are attained in a computational work that is a small (less than 10) multiple of the operation count in one target-grid residual evaluation. The way to achieve this efficiency is the distributed relaxation approach. TME solvers employing distributed relaxation have already been demonstrated for nonconservative formulations of high-Reynolds-number viscous incompressible and subsonic compressible flow regimes. The purpose of this paper is to provide foundations for applications of distributed relaxation to conservative discretizations. A direct correspondence between the primitive variable interpolations for calculating fluxes in conservative finite-volume discretizations and stencils of the discretized derivatives in the nonconservative formulation has been established. Based on this correspondence, one can arrive at a conservative discretization which is very efficiently solved with a nonconservative relaxation scheme and this is demonstrated for conservative discretization of the quasi one-dimensional Euler equations. Formulations for both staggered and collocated grid arrangements are considered and extensions of the general procedure to multiple dimensions are discussed.

Diskin, Boris↗

Application of a multi-level grid method to transonic flow calculations

A multi-level grid method has been studied as a possible means of accelerating convergence in relaxation calculations for transonic flows. The method employs a hierarchy of grids, ranging from very coarse (e.g., 8 x 2 mesh cells) to fine (e.g., 128 x 32); the coarser grids are used to diminish the magnitude of the smooth part of the residuals, hopefully with far less total work than would be required with, say, optimal SLOR iterations on the finest grid. The method was applied to the solution of the transonic small-disturbance equation for the velocity potential in the conservation form. Nonlifting transonic flow past a parabolic-arc airfoil is the example studied, with meshes of both constant and variable step size.

South, J. C., Jr.↗

Hyperbolic/parabolic development for the GIM-STAR code

Flow fields in supersonic inlet configurations were computed using the eliptic GIM code on the STAR computer. Spillage flow under the lower cowl was calculated to be 33% of the incoming stream. The shock/boundary layer interaction on the upper propulsive surface was computed including separation. All shocks produced by the flow system were captured. Linearized block implicit (LBI) schemes were examined to determine their application to the GIM code. Pure explicit methods have stability limitations and fully implicit schemes are inherently inefficient; however, LBI schemes show promise as an effective compromise. A quasiparabolic version of the GIM code was developed using elastical parabolized Navier-Stokes methods combined with quasitime relaxation. This scheme is referred to as quasiparabolic although it applies equally well to hyperbolic supersonic inviscid flows. Second order windward differences are used in the marching coordinate and either explicit or linear block implicit time relaxation can be incorporated.

Spradley, L. W.↗

Effect of Surface Roughness on Hydrodynamic Bearings

A theoretical analysis on the performance of hydrodynamic oil bearings is made considering surface roughness effect. The hydrodynamic as well as asperity contact load is found. The contact pressure was calculated with the assumption that the surface height distribution was Gaussian. The average Reynolds equation of partially lubricated surface was used to calculate hydrodynamic load. An analytical expression for average gap was found and was introduced to modify the average Reynolds equation. The resulting boundary value problem was then solved numerically by finite difference methods using the method of successive over relaxation. The pressure distribution and hydrodynamic load capacity of plane slider and journal bearings were calculated for various design data. The effects of attitude and roughness of surface on the bearing performance were shown. The results are compared with similar available solution of rough surface bearings. It is shown that: (1) the contribution of contact load is not significant; and (2) the hydrodynamic and contact load increase with surface roughness.

Majumdar, B. C.↗

MBE growth technology for high quality strained III-V layers

The III-V films are grown on large automatically perfect terraces of III-V substrates which have a different lattice constant, with temperature and Group III and V arrival rates chosen to give a Group III element stable surface. The growth is pulsed to inhibit Group III metal accumulation of low temperature, and to permit the film to relax to equilibrium. The method of the invention: (1) minimizes starting step density on sample surface; (2) deposits InAs and GaAs using an interrupted growth mode (0.25 to 2 monolayers at a time); (3) maintains the instantaneous surface stoichiometry during growth (As-stable for GaAs, In-stable for InAs); and (4) uses time-resolved RHEED to achieve aspects (1) through (3).

Grunthaner, Frank J.↗

Patch structured surface grid with dynamic curvature clustering

The numerical technique developed by Eiseman (1985) to generate surface grids with lines clustered in areas of greater curvature is described and demonstrated. This approach combines the direct active use of surface curvature with a mean-value relaxation scheme and multigrid methods; it permits an arbitrary number of patches to be organized in a single integer array based on a patch number, an edge number, and a boundary-condition index. Grids for typical CFD applications are presented and briefly characterized.

Wang, YI↗

Viscoelastic characterization of thin-film polymers exposed to low Earth orbit

The materials made available through the Long Duration Exposure Facility (LDEF) satellite provide a set of specimens that can be well characterized and have a known exposure history with reference to atomic oxygen and ultraviolet radiation exposure. Mechanical characteristics measured from control samples and exposed samples provide a data base for predicting the behavior of polymers in low earth orbit. Samples of 1.0 mil thick low density polyethylene were exposed to the low earth orbit environment for a period of six years. These materials were not directly exposed to ram atomic oxygen and offer a unique opportunity for measuring the effect of atomic oxygen and UV radiation on mechanical properties with little concern to the effect of erosion. The viscoelastic characteristics of these materials were measured and compared to the viscoelastic characteristics of control samples. To aid in differentiating the effects of changes in crystallinity resulting from thermal cycling, from the effects of changes in chemical structure resulting from atomic oxygen/UV attack to the polymer, a second set of control specimens, annealed to increase crystallinity, were measured as well. The resulting characterization of these materials will offer insight into the impact of atomic oxygen/UV on the mechanical properties of polymeric materials. The viscoelastic properties measured for the control, annealed, and exposed specimens were the storage and loss modulus as a function of frequency and temperature. From these datum is calculated the viscoelastic master curve derived using the principle of time/temperature superposition. Using the master curve, the relaxation modulus is calculated using the method of Ninomiya and Ferry. The viscoelastic master curve and the stress relaxation modulus provide a direct measure of the changes in the chemical or morphological structure. In addition, the effect of these changes on long-term and short-term mechanical properties is known directly. It should be noted that the dependence on directionality for the polymer films was considered since these films were manufactured by a blown-film process.

Letton, Alan↗

Algorithms for Multiple Fault Diagnosis With Unreliable Tests

In this paper, we consider the problem of constructing optimal and near-optimal multiple fault diagnosis (MFD) in bipartite systems with unreliable (imperfect) tests. It is known that exact computation of conditional probabilities for multiple fault diagnosis is NP-hard. The novel feature of our diagnostic algorithms is the use of Lagrangian relaxation and subgradient optimization methods to provide: (1) near optimal solutions for the MFD problem, and (2) upper bounds for an optimal branch-and-bound algorithm. The proposed method is illustrated using several examples. Computational results indicate that: (1) our algorithm has superior computational performance to the existing algorithms (approximately three orders of magnitude improvement), (2) the near optimal algorithm generates the most likely candidates with a very high accuracy, and (3) our algorithm can find the most likely candidates in systems with as many as 1000 faults.

Shakeri, Mojdeh↗

MBE growth technology for high quality strained III-V layers

III-V films are grown on large automatically perfect terraces of III-V substrates which have a different lattice constant, with temperature and Group II and V arrival rates chosen to give a Group III element stable surface. The growth is pulsed to inhibit Group III metal accumulation to low temperature, and to permit the film to relax to equilibrium. The method of the invention 1) minimizes starting step density on sample surface; 2) deposits InAs and GaAs using an interrupted growth mode (0.25 to 2 mono-layers at a time); 3) maintains the instantaneous surface stoichiometry during growth (As-stable for GaAs, In-stable for InAs); and 4) uses time-resolved RHEED to achieve aspects (1)-14 (3).

Grunthaner, Frank J.↗

Flux-split algorithms for flows with non-equilibrium chemistry and vibrational relaxation

The present consideration of numerical computation methods for gas flows with nonequilibrium chemistry thermodynamics gives attention to an equilibrium model, a general nonequilibrium model, and a simplified model based on vibrational relaxation. Flux-splitting procedures are developed for the fully-coupled inviscid equations encompassing fluid dynamics and both chemical and internal energy-relaxation processes. A fully coupled and implicit large-block structure is presented which embodies novel forms of flux-vector split and flux-difference split algorithms valid for nonequilibrium flow; illustrative high-temperature shock tube and nozzle flow examples are given.

Grossman, B.↗