SEARCH · Search NASA
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.
Azimuthal sampling of electron flux on spinning satellite provide direction and magnitude of energetic electron space gradients at trapping magnetospheric boundaries
A new approach is presented for analyzing compressor tip clearance flow. The basic idea is that the clearance velocity field can be (approximately) decomposed into independent throughflow and crossflow, since chordwise pressure gradients are much smaller than normal pressure gradients in the clearance region. As in the slender body approximation in external aerodynamics, this description implies that the three-dimensional steady clearance flow can be viewed as a two-dimensional, unsteady flow. Using this approach, a similarity scaling for the crossflow in the clearance region is developed and a generalized description of the clearance vortex is derived. Calculations based on the similarity scaling agree well with a wide range of experimental data in regard to flow features such as crossflow velocity field, static pressure field, and tip clearance vortex trajectory.
The structure and the thermophysical properties of d-UZr 2 are investigated using 0 K density functional theory and ab initio molecular dynamics (AIMD). Modeling the true paramagnetic state of this intermetallic compound has been challenging using first-principles calculations. For the first time, we find that the generalized gradient approximation method without applying an on-site Coulomb interaction term (Hubbard U) can result in a ground state that is antiferromagnetic (AFM). We believe that this weak AFM ground state is the closest to the real paramagnetic state. We found that structure optimization at finite temperatures using AIMD is necessary to achieve this ground state instead of the non-magnetic and ferromagnetic states previously reported in the literature. Our findings indicate that applying the Hubbard U on uranium f-orbitals in this metallic system is unnecessary and not recommended, as it leads to a large overestimation of the volume and introduces an unphysical strong spin polarization. Finally, our approach results in atomic volume, thermal expansion, and heat capacities that have strong agreement with experiments.
MoNbTaVW, a pioneering refractory high-entropy alloy (RHEA), is renowned for its exceptional strength at elevated temperatures. Like most RHEAs, it solidifies into a dendritic microstructure with steep concentration gradients, necessitating heat treatment for equilibration. Steep concentration gradients typically complicate mechanical property analysis and modeling, key challenges for high-throughput alloy discovery. This study examines the effects of high-temperature homogenization (1800 °C for 8 hours) on the microstructure and mechanical properties of MoNbTaVW across temperatures up to 1200 °C. Postmortem microstructural analysis, coupled with chemical, crystallographic, and mechanical assessments, revealed that yield strength was surprisingly insensitive to homogenization, with less than a 3 pct difference between as-cast and heat-treated materials. Nanoindentation confirmed minimal nanohardness changes across the dendritic structure at room temperature. However, homogenization significantly enhanced high-temperature work hardening, producing higher peak compressive strength. Both as-cast and homogenized MoNbTaVW exhibited room-temperature strengths of 1400 MPa, exceeding previously reported values. These findings demonstrate that high-temperature treatment enables microstructural homogenization without compromising strength, a unique behavior among RHEAs. In conclusion, while not universally applicable, this insight highlights the importance of understanding microstructural development in experimental alloys and offers a pathway to design alloys that achieve homogeneous-like properties without heat treatment.
The diffusion of gases confined in nanoporous materials underpins membrane and adsorption-based gas separations, yet relatively few measurements of diffusion coefficients in the promising class of materials, metal–organic frameworks (MOFs), have been reported to date. Recently we reported self-diffusion coefficients for 13 CO 2 in the MOF Zn 2 (dobpdc) (dobpdc 4– = 4,4'-dioxidobiphenyl-3,3'-dicarboxylate) which has one-dimensional channels with a diameter of approximately 2 nm [Forse, A. C.; J. Am. Chem. Soc. 2018, 140, 1663-1673]. By analyzing the evolution of the residual 13 C chemical shift anisotropy line shape at different gradient strengths, we obtained self-diffusion coefficients both along (D ∥ ) and between (D ⊥ ) the one-dimensional MOF channels. The observation of nonzero D ⊥ was unexpected based on the single crystal X-ray diffraction structure and flexible lattice molecular dynamics simulations, and we proposed that structural defects may be responsible for self-diffusion between the MOF channels. In this paper we revisit this analysis and show that homogeneous line broadening must be taken into account to obtain accurate values for D ⊥ . In the presence of homogeneous line broadening, intensity at a particular NMR frequency represents signal from crystals with a range of orientations relative to the applied magnetic field and magnetic gradient field. To quantify these effects, we perform spectral simulations that take into account homogeneous broadening and allow improved D ⊥ values to be obtained. Finally, our new analysis best supports nonzero D ⊥ at all studied dosing pressures and shows that our previous analysis overestimated D ⊥ .
The primary focus of this paper is to investigate how a rotor's unsteady tip clearance flow structure changes in a low speed one and half stage axial compressor when the rotor tip gap size is increased from 0.5 mm (0.49% of rotor tip blade chord, 2% of blade span) to 2.4 mm (2.34% chord, 4% span) at the design condition are investigated. The changes in unsteady tip clearance flow with the 0.62 % tip gap as the flow rate is reduced to near stall condition are also investigated. A Large Eddy Simulation (LES) is applied to calculate the unsteady flow field at these three flow conditions. Detailed Stereoscopic PIV (SPIV) measurements of the current flow fields were also performed at the Johns Hopkins University in a refractive index-matched test facility which renders the compressor blades and casing optically transparent. With this setup, the unsteady velocity field in the entire flow domain, including the flow inside the tip gap, can be measured. Unsteady tip clearance flow fields from LES are compared with the PIV measurements and both LES and PIV results are used to study changes in tip clearance flow structures. The current study shows that the tip clearance vortex is not a single structure as traditionally perceived. The tip clearance vortex is formed by multiple interlaced vorticities. Therefore, the tip clearance vortex is inherently unsteady. The multiple interlaced vortices never roll up to form a single structure. When phased-averaged, the tip clearance vortex appears as a single structure. When flow rate is reduced with the same tip gap, the tip clearance vortex rolls further upstream and the tip clearance vortex moves further radially inward and away from the suction side of the blade. When the tip gap size is increased at the design flow condition, the overall tip clearance vortex becomes stronger and it stays closer to the blade suction side and the vortex core extends all the way to the exit of the blade passage. Measured and calculated unsteady flow fields inside the tip gap agree fairly well. Instantaneous velocity vectors inside the tip gap from both the PIV and LES do show flow separation and reattachment at the entrance of tip gap as some earlier studies suggested. This area at the entrance of tip gap flow (the pressure side of the blade) is confined very close to the rotor tip section. With a small tip gap (0.5mm), the gap flow looks like a simple two-dimensional channel flow with larger velocity near the casing for both flow rates. A small area with a sharp velocity gradient is observed just above the rotor tip. This strong shear layer is turned radially inward when it collides with the incoming flow and forms the core structure of the tip clearance vortex. When tip gap size is increased to 2.4 mm at the design operation, the radial profile of the tip gap flow changes drastically. With the large tip gap, the gap flow looks like a two-dimensional channel flow only near the casing. Near the rotor top section, a bigger region with very large shear and reversed flow is observed.
The primary focus of this paper is to investigate how a rotor's unsteady tip clearance flow structure changes in a low speed one and half stage axial compressor when the rotor tip gap size is increased from 0.5 mm (0.49% of rotor tip blade chord, 2% of blade span) to 2.4 mm (2.34% chord, 4% span) at the design condition are investigated. The changes in unsteady tip clearance flow with the 0.62 % tip gap as the flow rate is reduced to near stall condition are also investigated. A Large Eddy Simulation (LES) is applied to calculate the unsteady flow field at these three flow conditions. Detailed Stereoscopic PIV (SPIV) measurements of the current flow fields were also performed at the Johns Hopkins University in a refractive index-matched test facility which renders the compressor blades and casing optically transparent. With this setup, the unsteady velocity field in the entire flow domain, including the flow inside the tip gap, can be measured. Unsteady tip clearance flow fields from LES are compared with the PIV measurements and both LES and PIV results are used to study changes in tip clearance flow structures. The current study shows that the tip clearance vortex is not a single structure as traditionally perceived. The tip clearance vortex is formed by multiple interlaced vorticities. Therefore, the tip clearance vortex is inherently unsteady. The multiple interlaced vortices never roll up to form a single structure. When phased-averaged, the tip clearance vortex appears as a single structure. When flow rate is reduced with the same tip gap, the tip clearance vortex rolls further upstream and the tip clearance vortex moves further radially inward and away from the suction side of the blade. When the tip gap size is increased at the design flow condition, the overall tip clearance vortex becomes stronger and it stays closer to the blade suction side and the vortex core extends all the way to the exit of the blade passage. Measured and calculated unsteady flow fields inside the tip gap agree fairly well. Instantaneous velocity vectors inside the tip gap from both the PIV and LES do show flow separation and reattachment at the entrance of tip gap as some earlier studies suggested. This area at the entrance of tip gap flow (the pressure side of the blade) is confined very close to the rotor tip section. With a small tip gap (0.5mm), the gap flow looks like a simple two-dimensional channel flow with larger velocity near the casing for both flow rates. A small area with a sharp velocity gradient is observed just above the rotor tip. This strong shear layer is turned radially inward when it collides with the incoming flow and forms the core structure of the tip clearance vortex. When tip gap size is increased to 2.4 mm at the design operation, the radial profile of the tip gap flow changes drastically. With the large tip gap, the gap flow looks like a two-dimensional channel flow only near the casing. Near the rotor top section, a bigger region with very large shear and reversed flow is observed.
The electronic band structure of correlated Ca 3 Ru 2 O 7 featuring an antiferromagnetic (AFM) as well as a structural transition has been determined theoretically at high temperatures, which has led to the understanding of the remarkable properties of Ca 3 Ru 2 O 7 , such as the bulk spin-valve effects. However, its band structure and Fermi surface (FS) below the structural transition have not been resolved, even though a FS consisting of electron pockets was found experimentally. Here we report magnetoelectrical transport and thermoelectric measurements with the electric current and temperature gradient directed along the a and b axes, respectively, of an untwined single crystal of Ca 3 Ru 2 O 7 . The thermopowers obtained along the two crystal axes were found to show opposite signs at low temperatures, demonstrating the presence of both electron and hole pockets on the FS. In addition, how the FS evolves across T* = 30K at which a distinct transition from coherent to incoherent behavior occurs was also inferred: the Hall and Nernst coefficient results suggest a temperature- and momentum-dependent partial gap opening in Ca 3 Ru 2 O 7 below the structural transition with a possible Lifshitz transition occurring at T*. In conclusion, the experimental demonstration of a correlated semimetal ground state in Ca 3 Ru 2 O 7 calls for further theoretical studies of this remarkable material.
During the Flash process, the cross section of a plain-carbon or a low-alloy steel is austenitized through rapid heating and transformed on rapid cooling to a predominantly martensite + bainite structure with small amounts of retained austenite. Unlike conventional heat treating, homogeneity is intentionally avoided during Flash processing of steels. The Flash process assembly consists of a pair of rolls that transfer the steel sheets through the heating and cooling stage of the thermal cycle. The initial microstructure of the steel consists of ferrite (body-centered cubic iron) + carbide ((Fe,X)mCn) mixture. The heating rate through the peak temperature is a function of temperature and reaches a peak of about 300-400°C/s and the cooling rate has a maximum value of 3,000-4,000°C/s. The on-heating phase transformations include carbide dissolution, austenite (face-centered cubic iron) nucleation and growth, and diffusion of carbon and other substitutional elements in the steel. The on-cooling phase transformations include formation of martensite (body-centered tetragonal phase containing supersaturated solute) and bainite (ferrite plates with or without fine carbides). In this project, the focus is on Fe-C-Cr steels that are currently Flash processed for armor applications. The modeling effort proposed here will help optimize the Flash thermal cycle for these low alloy steels to achieve the target performance, which is an ongoing effort at SFP Works. A significant feature of Flash processed Fe-C-Cr steels is the presence of scatter in the through-thickness in the sheet. The variability in hardness results from a variability in the bainite + martensite microstructure that is sensitive to the local chemical concentration of C and Cr. Such a chemical inhomogeneity is intentionally obtained in the Flash process. Although such a microstructural gradient is presumably responsible for the exceptional properties of the Flash processed steel, it is very important to quantify the gradients as a function of Flash variabilities in processing parameters and the input microstructure. Understanding the mechanistic pathway that leads to microstructural gradients could be ground-breaking and instrumental for achieving better process control and optimized microstructural state to meet application-specific strength-ductility requirements. Since the final microstructure depends on setting up precise solute concentration gradients through a rapid heating process, and transforming these regions into various phases, it is important to understand how small changes in steel chemistry, input microstructure (carbide size and distribution), and process variables (Flash thermal cycle) will impact the solute concentration gradients.
Abstract not provided.
The quality of simulated hypersonic stagnation region heating with tetrahedral meshes is investigated by using an updated three-dimensional, upwind reconstruction algorithm for the inviscid flux vector. An earlier implementation of this algorithm provided improved symmetry characteristics on tetrahedral grids compared to conventional reconstruction methods. The original formulation however displayed quantitative differences in heating and shear that were as large as 25% compared to a benchmark, structured-grid solution. The primary cause of this discrepancy is found to be an inherent inconsistency in the formulation of the flux limiter. The inconsistency is removed by employing a Green-Gauss formulation of primitive gradients at nodes to replace the previous Gram-Schmidt algorithm. Current results are now in good agreement with benchmark solutions for two challenge problems: (1) hypersonic flow over a three-dimensional cylindrical section with special attention to the uniformity of the solution in the spanwise direction and (2) hypersonic flow over a three-dimensional sphere. The tetrahedral cells used in the simulation are derived from a structured grid where cell faces are bisected across the diagonal resulting in a consistent pattern of diagonals running in a biased direction across the otherwise symmetric domain. This grid is known to accentuate problems in both shock capturing and stagnation region heating encountered with conventional, quasi-one-dimensional inviscid flux reconstruction algorithms. Therefore the test problems provide a sensitive indicator for algorithmic effects on heating. Additional simulations on a sharp, double cone and the shuttle orbiter are then presented to demonstrate the capabilities of the new algorithm on more geometrically complex flows with tetrahedral grids. These results provide the first indication that pure tetrahedral elements utilizing the updated, three-dimensional, upwind reconstruction algorithm may be used for the simulation of heating and shear in hypersonic flows in upwind, finite volume formulations.
The "4-day wave" is an eastward moving quasi-nondispersive feature with period near 4 days occurring near the winter polar stratopause. This paper presents evidence of the 4-day feature in Microwave Limb Sounder (MLS) temperature, geopotential height, and ozone data from the late southern winters of 1992 and 1993. Space-time spectral analyses reveal a double-peaked temperature structure consisting of one peak near the stratopause and another in the lower mesosphere, with an out-of-phase relationship between the two peaks. This double- peaked structure is reminiscent of recent three-dimensional barotropic/baroclinic instability model predictions and is observed here for the first time. The height variation of the 4-day ozone signal is shown to compare well with a linear advective-photochemical tracer model. Negative regions of quasigeostrophic potential vorticity (PV) gradient and positive Eliassen-Palm flux divergence are shown to occur, consistent with instability dynamics playing a role in wave forcing. Spectral analyses of PV derived from MLS geopotential height fields reveal a 4-day signal peaking near the polar stratopause. The three-dimensional structure of the 4-day wave resembles the potential vorticity "charge" concept, wherein a PV anomaly in the atmosphere (analogous to an electrical charge in a dielectric material) induces a geopotential field, a vertically oriented temperature dipole, and circulation about the vertical axis.
We report first-principles calculations of the elastic properties, electronic structure and magnetic behavior performed over the Ba{sub 2}NiMoO{sub 6} double perovskite. Calculations are carried out through the full-potential linear augmented plane-wave method within the framework of the Density Functional Theory (DFT) with exchange and correlation effects in the Generalized Gradient and Local Density Approximations, including spin polarization. The elastic properties calculated are bulk modulus (B), the elastic constants (C{sub 11}, C{sub 12} and C{sub 44}), the Zener anisotropy factor (A), the isotropic shear modulus (G), the Young modulus (Y) and the Poisson ratio (υ). Structural parameters, total energies and cohesive properties of the perovskite are studied by means of minimization of internal parameters with the Murnaghan equation, where the structural parameters are in good agreement with experimental data. Furthermore, we have explored different antiferromagnetic configurations in order to describe the magnetic ground state of this compound. The pressure and temperature dependence of specific heat, thermal expansion coefficient, Debye temperature and Grüneisen parameter were calculated by DFT from the state equation using the quasi-harmonic model of Debye. A specific heat behavior C{sub V} ≈ C{sub P} was found at temperatures below T = 400 K, with Dulong–Petit limit values, which is higher than those, reported for simple perovskites.
Loss of photosynthetic area can affect soil microbial communities by altering the availability of fixed carbon. We used denaturing gradient gel electrophoresis (DGGE) and Biolog filamentous-fungus plates to determine the effects of artificial defoliation of pines in a mixed pine-spruce forest on the composition of the fungal community in a forest soil. As measured by DGGE, two fungal species were affected significantly by the defoliation of pines (P < 0.001); the frequency of members of the ectomycorrhizal fungus genus Cenococcum decreased significantly, while the frequency of organisms of an unidentified soil fungus increased. The decrease in the amount of Cenococcum organisms may have occurred because of the formation of extensive hyphal networks by species of this genus, which require more of the carbon fixed by their host, or because this fungus is dependent upon quantitative differences in spruce root exudates. The defoliation of pines did not affect the overall composition of the soil fungal community or fungal-species richness (number of species per core). Biolog filamentous-fungus plate assays indicated a significant increase (P < 0.001) in the number of carbon substrates utilized by the soil fungi and the rate at which these substrates were used, which could indicate an increase in fungal-species richness. Thus, either small changes in the soil fungal community give rise to significant increases in physiological capabilities or PCR bias limits the reliability of the DGGE results. These data indicate that combined genetic and physiological assessments of the soil fungal community are needed to accurately assess the effect of disturbance on indigenous microbial systems.
The quasiparticle band structure can be properly described by Hedin's GW approximation (GW), at a high computational cost. For band gaps, semilocal density functionals up to the generalized gradient approximation (GGA) level cannot compete with the accuracy of hybrid-based approximations or GW. Meta-GGA density functionals with a strong dependence on the kinetic energy density ingredient can potentially give wider band gaps compared with GGAs. The recent TASK meta-GGA density functional from Aschebrock and Kümmel [ Phys. Rev. Research 1 , 033082 (2019) ], is constructed with an enhanced nonlocality in the generalized Kohn-Sham scheme and therefore harbors great opportunities for band gap prediction. Although this approximation was found to yield excellent band gaps of bulk solids, this accuracy cannot be straightforwardly transferred to low-dimensional materials. Additionally, the reduced screening of these materials results in larger band gaps compared with their bulk counterparts, as an additional barrier to overcome. In this paper we demonstrate how the alteration of this functional affects the band gaps of monolayers and nanoribbons and present accurate band gaps competing with the revised Heyd-Scuseria-Ernzerhof (HSE06) approximation. In order to achieve this goal, we have modified the TASK functional (a) by changing the tight upper bound for one- or two-electron systems ( h X 0 ) from 1.174 to 1.29 and (b) by changing the limit of the interpolation function f X ( α → ∞ ) of the TASK functional that interpolates the exchange enhancement factor F X ( s , α ) from α = 0 to 1. The resulting modified TASK (mTASK) was tested for various materials from three dimensions to two dimensions to one dimension (nanoribbons) and was compared with the results of the higher-level hybrid functional HSE06 or with the G 0 W 0 approximation within many-body perturbation theory. We find that mTASK systematically improves the band gaps and band structures of two-dimensional (2D) and 1D systems, without significantly affecting the accuracy of the original TASK for the bulk 3D materials, when compared with the Perdew-Burke-Ernzerhof (PBE) GGA and the strongly constrained and appropriately normed (SCAN) meta-GGA. We further demonstrate the applicability of mTASK by assessing the band structures of transition metal dichalcogenide nanoribbons with respect to various bending curvatures.
Stencil loops are a common motif in computations including convolutional neural networks, structured-mesh solvers for partial differential equations, and image processing. Stencil loops are easy to parallelise, and their fast execution is aided by compilers, libraries, and domain-specific languages. Reverse-mode automatic differentiation, also known as algorithmic differentiation, autodiff, adjoint differentiation, or back-propagation, is sometimes used to obtain gradients of programs that contain stencil loops. Unfortunately, conventional automatic differentiation results in a memory access pattern that is not stencil-like and not easily parallelisable. In this paper we present a novel combination of automatic differentiation and loop transformations that preserves the structure and memory access pattern of stencil loops, while computing fully consistent derivatives. The generated loops can be parallelised and optimised for performance in the same way and using the same tools as the original computation. We have implemented this new technique in the Python tool PerforAD, which we release with this paper along with test cases derived from seismic imaging and computational fluid dynamics applications.
COMETBOARDS (Comparative Evaluation Testbed of Optimization and Analysis Routines for Design of Structures) is a design optimization test bed that can evaluate the performance of several different optimization algorithms. A few of these optimization algorithms are the sequence of unconstrained minimization techniques (SUMT), sequential linear programming (SLP) and the sequential quadratic programming techniques (SQP). A genetic algorithm (GA) is a search technique that is based on the principles of natural selection or "survival of the fittest". Instead of using gradient information, the GA uses the objective function directly in the search. The GA searches the solution space by maintaining a population of potential solutions. Then, using evolving operations such as recombination, mutation and selection, the GA creates successive generations of solutions that will evolve and take on the positive characteristics of their parents and thus gradually approach optimal or near-optimal solutions. By using the objective function directly in the search, genetic algorithms can be effectively applied in non-convex, highly nonlinear, complex problems. The genetic algorithm is not guaranteed to find the global optimum, but it is less likely to get trapped at a local optimum than traditional gradient-based search methods when the objective function is not smooth and generally well behaved. The purpose of this research is to assist in the integration of genetic algorithm (GA) into COMETBOARDS. COMETBOARDS cast the design of structures as a constrained nonlinear optimization problem. One method used to solve constrained optimization problem with a GA to convert the constrained optimization problem into an unconstrained optimization problem by developing a penalty function that penalizes infeasible solutions. There have been several suggested penalty function in the literature each with there own strengths and weaknesses. A statistical analysis of some suggested penalty functions is performed in this study. Also, a response surface approach to robust design is used to develop a new penalty function approach. This new penalty function approach is then compared with the other existing penalty functions.
The Low Power Atmospheric Compensation Experiment (LACE) is scheduled for launch in late 1989 into a 556 km altitude circular orbit of 43 deg inclination. The LACE flight dynamics experiment is an experiment secondary to the primary LACE mission. The purpose of the experiment is to provide on-orbit systems identification of the LACE spacecraft. The structure of the LACE spacecraft is of special interest to the CSI community. It incorporates 3 deployable/retractable booms of maximum length 45.72 m (150 ft) mounted on a rectangular parallelepiped bus of mass 1,200 kg. The zenith directed gravity gradient boom is mounted on the top of the bus; the retroreflector boom is mounted forward and deployed along the velocity vector, the balance boom is mounted and pointed aft. Attitude stabilization is accomplished by means of gravity gradient torques and by a momentum wheel. The LACE flight dynamics experiment is designed to measure modal frequencies, damping ratios, and oscillation amplitudes of the LACE spacecraft, as well as the vibration intensity generated by boom deployments and retractions. It is anticipated that this experiment will provide an opportunity for improvements in the accuracy of computer simulations of flexible structures and multibody dynamics.