Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stability boundary”

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,333 records · Page 74

Robust, High-Performing Maize–Perovskite-Based Solar Cells with Improved Stability

Herein, we focus on improving the long-term chemical and thermomechan-ical stability of perovskite solar cells (PSCs), two major challenges currently limiting their commercial deployment. Our strategy incorporates a long-chain starch polymer into the perovskite precursor. The starch polymer confers multiple beneficial effects by forming hydrogen bonds with the methylammonium iodide precursor, templating perovskite growth that results in a compact and homogeneous film deposited in a simple one-step coating (antisolvent-free). The inclusion of starch in the methylammonium lead iodide films strongly improves their thermomechanical and environmental stability while maintaining a high photovoltaic performance. The fracture energy (G c ) of the film is increased to above 5 J/m 2 by creating a nanocomposite that provides intrinsic reinforcement at grain boundaries. Additionally, improved optoelectronic properties achieved with the starch polymer enable good photostability of the active layer and enhanced resistance to thermal cycling.

14 SOLAR ENERGY↗

Numerical tests of a modified full implicit continuous Eulerian (FICE) scheme with projected normal characteristic boundary conditions for MHD flows

A numerical method has been developed based on a modified full implicit continuous Eulerian (FICE) scheme and projected normal characteristic boundary conditions for simulating MHD flows which undergo a long process of evolution. An astrophysical flow is chosen for illustration of this procedure, and numerical tests are made to verify the computational stability and physically realistic solution. Three computational tests have been accomplished; they are tests of solving methods, characteristic boundary condition, and time steps. The tests show that the program from the modified FICE scheme with proper boundary conditions and time steps can be made numerically stable for a time long enough to obtain physically plausible solutions.

Wu, S. T.↗

High temperature stabilization of ultrafine grain tungsten alloys through synergistic compositional complexities

Thermally-stabilized fine-grained microstructures in tungsten provide a pathway to harnessing enhanced properties such as a reduced ductile-to-brittle transition temperature and improved strength while mitigating the adverse effects of grain growth and recrystallization. Here, in this study, we employ a material design strategy that relies on grain boundary segregation in the nanocrystalline state driven by alloy thermodynamics balanced with impurity scavenging through in situ formation of kinetically-stabilizing metal carbides. A W-Ti-Cr alloy is designed through lattice Monte Carlo methods and subsequently synthesized in a single-phase nanocrystalline state, which upon annealing, evolves into an ultrafine grained microstructure containing chromium grain boundary segregation collectively with a dispersed TiX (X=C,O) phase through the reaction of titanium with carbon and oxygen impurities. In situ synchrotron X-ray diffraction experiments demonstrate that increased Ti promotes stabilization across a larger temperature range but with diminishing returns above 10 at.% Ti. Long-term stability was confirmed through the retention of the ultrafine grained microstructure for a total of 8 days (192 h) at 1300 °C without grain/carbide growth and/or recrystallization. Our results demonstrate that, through strategic tailoring of composition and microstructure, one can harness the benefits of both thermodynamic and kinetic stabilization mechanisms, opening pathways for future alloy formulations that expand the window of stability.

36 MATERIALS SCIENCE↗

Calculation of three-dimensional compressible laminar and turbulent boundary layers. An implicit finite-difference procedure for solving the three-dimensional compressible laminar, transitional, and turbulent boundary-layer equations

An implicit finite-difference procedure is presented for solving the compressible three-dimensional boundary-layer equations. The method is second-order accurate, unconditionally stable (conditional stability for reverse cross flow), and efficient from the viewpoint of computer storage and processing time. The Reynolds stress terms are modeled by (1) a single-layer mixing length model and (2) a two-layer eddy viscosity model. These models, although simple in concept, accurately predicted the equilibrium turbulent flow for the conditions considered. Numerical results are compared with experimental wall and profile data for a cone at an angle of attack larger than the cone semiapex angle. These comparisons clearly indicate that the numerical procedure and turbulence models accurately predict the experimental data with as few as 21 nodal points in the plane normal to the wall boundary.

Harris, J. E.↗

Developing reliable machine learning interatomic potential for Fe–Cr–Ni austenitic alloys

Gaining atomistic understanding of mechanical behavior of heat-resistant structural materials such as Fe–Cr–Ni-based alloys requires an approach with an accuracy close to density functional theory (DFT) that considers the intrinsic properties of the bulk lattice and important defects such as stacking faults, grain boundaries, and surfaces. This work aims to develop reliable machine learning interatomic potential (MLIAP) at cross-scale for Fe–Cr–Ni ternary alloys with a focus on the face-centered-cubic (fcc) solid solution structure. Leveraging the advantages of moment tensor potentials, which typically necessitate a relatively small training dataset and enable rapid calculations using the large-scale atomic/molecular massively parallel simulator package, we ensure the stability and accuracy of the trained potentials. Important defects such as stacking faults, grain boundaries, and surfaces for wide-range compositions are investigated. Structural, thermal, elastic, and defect properties are determined from molecular dynamics simulations comprising several thousand atoms, generated via canonical Monte Carlo simulations guided by the trained potential. The trained potential allows efficient atomic simulations of structural, thermal, and mechanical properties of fcc Fe–Cr–Ni solid solution alloys as a function of composition and temperature. Therefore, the MLIAP approach represents a major advancement from DFT calculations that are limited to small simulation sizes and traditional molecular dynamics simulations using relatively low accuracy potentials. Furthermore, this work outlines a practical foundation for further investigating the structural evolution and mechanical behavior of austenitic stainless steel and nickel-based alloys in a wide array of applications in extreme environments.

Crystal structure↗

Entropy Splitting and Numerical Dissipation

A rigorous stability estimate for arbitrary order of accuracy of spatial central difference schemes for initial-boundary value problems of nonlinear symmetrizable systems of hyperbolic conservation laws was established recently by Olsson and Oliger (1994) and Olsson (1995) and was applied to the two-dimensional compressible Euler equations for a perfect gas by Gerritsen and Olsson (1996) and Gerritsen (1996). The basic building block in developing the stability estimate is a generalized energy approach based on a special splitting of the flux derivative via a convex entropy function and certain homogeneous properties. Due to some of the unique properties of the compressible Euler equations for a perfect gas, the splitting resulted in the sum of a conservative portion and a non-conservative portion of the flux derivative. hereafter referred to as the "Entropy Splitting." There are several potential desirable attributes and side benefits of the entropy splitting for the compressible Euler equations that were not fully explored in Gerritsen and Olsson. The paper has several objectives. The first is to investigate the choice of the arbitrary parameter that determines the amount of splitting and its dependence on the type of physics of current interest to computational fluid dynamics. The second is to investigate in what manner the splitting affects the nonlinear stability of the central schemes for long time integrations of unsteady flows such as in nonlinear aeroacoustics and turbulence dynamics. If numerical dissipation indeed is needed to stabilize the central scheme, can the splitting help minimize the numerical dissipation compared to its un-split cousin? Extensive numerical study on the vortex preservation capability of the splitting in conjunction with central schemes for long time integrations will be presented. The third is to study the effect of the non-conservative proportion of splitting in obtaining the correct shock location for high speed complex shock-turbulence interactions. The fourth is to determine if this method can be extended to other physical equations of state and other evolutionary equation sets. If numerical dissipation is needed, the Yee, Sandham, and Djomehri (1999) numerical dissipation is employed. The Yee et al. schemes fit in the Olsson and Oliger framework.

Yee, H. C.↗

Density limit in peeling-limited pedestals at and above the Greenwald value in DIII-D high poloidal beta plasmas

High pressure, peeling limited pedestals with pedestal normalized beta β N,ped >2 and pedestal top density n ped at or above the Greenwald density nG have been achieved in DIII-D high poloidal beta plasmas, with high global normalized beta β N >3 and energy confinement H 98 ~1.2-1.7. Higher β N allows higher pedestal density above the Greenwald value and higher pedestal pressure, even with a low injected torque. MHD modeling confirms that the experimental profiles lie near the peeling-mode unstable boundary with high normalized pressure gradient βMHD and high edge current density. Experimental analysis and stability calculations indicate that the high poloidal beta with strong Shafranov shift, high βMHD and weak/negative magnetic shear improves the pedestal stability by decoupling the peeling and ballooning modes and stabilizing the ballooning modes, thus facilitating access to the second stable region of peeling-ballooning mode. The access to the second stable peeling-ballooning stability region opens Super-H-like channels without extremely strong shaping or strong torque injection. The high-pressure peeling pedestal allows the pedestal density to go beyond the Greenwald limit with strong ExB shear maintained: pedestal pressure increases with pedestal density even when n ped >n G , until reaching the ideal MHD instability boundary, where giant ELMs occur. The giant ELMs are dominated by a strong n=1 component and cause a large reduction of the edge pressure, but a negligible change of the core pressure, consistent with kink/peeling-mode induced instability. The pedestal recovers from the collapse and typically sustains a high baseline density, around the Greenwald limit, during the whole discharge duration. Experiments also found that internal transport barriers and n ped ~0.9nG, peeling limited pedestals could be simultaneously achieved in high β N plasmas, while an internal feedback between ITB strength and pedestal performance is found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Receptivity of Hypersonic Boundary Layers to Acoustic and Vortical Disturbances (Invited)

Boundary-layer receptivity to two-dimensional acoustic and vortical disturbances for hypersonic flows over two-dimensional and axi-symmetric geometries were numerically investigated. The role of bluntness, wall cooling, and pressure gradients on the receptivity and stability were analyzed and compared with the sharp nose cases. It was found that for flows over sharp nose geometries in adiabatic wall conditions the instability waves are generated in the leading-edge region and that the boundary layer is much more receptive to slow acoustic waves as compared to the fast waves. The computations confirmed the stabilizing effect of nose bluntness and the role of the entropy layer in the delay of boundary layer transition. The receptivity coefficients in flows over blunt bodies are orders of magnitude smaller than that for the sharp cone cases. Wall cooling stabilizes the first mode strongly and destabilizes the second mode. However, the receptivity coefficients are also much smaller compared to the adiabatic case. The adverse pressure gradients increased the unstable second mode regions.

Balakumar, P.↗

DHARMA - Discriminant hyperplane abstracting residuals minimization algorithm for separating clusters with fuzzy boundaries

Learning of discriminant hyperplanes in imperfectly supervised or unsupervised training sample sets with unreliably labeled samples along the fuzzy joint boundaries between sample clusters is discussed, with the discriminant hyperplane designed to be a least-squares fit to the unreliably labeled data points. (Samples along the fuzzy boundary jump back and forth from one cluster to the other in recursive cluster stabilization and are considered unreliably labeled.) Minimization of the distances of these unreliably labeled samples from the hyperplanes does not sacrifice the ability to discriminate between classes represented by reliably labeled subsets of samples. An equivalent unconstrained linear inequality problem is formulated and algorithms for its solution are indicated. Landsat earth sensing data were used in confirming the validity and computational feasibility of the approach, which should be useful in deriving discriminant hyperplanes separating clusters with fuzzy boundaries, given supervised training sample sets with unreliably labeled boundary samples.

Dasarathy, B. V.↗

Intermittency Model for Coupled CFD-Stability Transition Analysis in Hypersonic Flow

Accurate prediction of aerothermodynamic loads on hypersonic vehicles requires precise modeling of surface quantities across the boundary layer transition zone. The peak heating loads and total heat transfer are determined by parameters such as the transition zone length and the magnitude of potential overshoots in heat flux and skin friction beyond their respective values in fully turbulent flows. While previous studies on CFD integrated transition modeling using linear stability correlations have shown promise in modeling these features, they did not develop an intermittency model for hypersonic flows. This paper presents a data-driven approach to develop a model of this type by correlating experimental transition data with the relevant flow parameters. The resulting model demonstrates significant improvements over previous low-speed models in terms of predicting heat transfer distributions during the transition process associated with first and second mode instabilities in axisymmetric high-speed flows. Separate correlations for flight and ground test conditions are developed, and the potential to combine these correlations is discussed.

Transition↗

Computational Study of Low-Speed Boundary-Layer Laminar-Turbulent Transition over a 30º Swept Plate in a Channel

Reliable boundary-layer laminar-turbulent transition predictions are critical for the design of aircraft that use natural or controlled laminar flow to achieve low drag and better fuel efficiency. During cruise flight through low turbulence, the three-dimensional boundary layer over swept surfaces can transition prematurely due to stationary crossflow instabilities induced by surface roughness. Parabolized stability equations analyses have been successfully used to explain the different phases of the transition process under well controlled perturbations such as steps or discrete roughness elements. However, methodologies are still lacking to model the boundary layer receptivity to distributed random roughness, and the following non-linear instability stages leading to transition. For engineering applications, several semi-empirical methods based on linear amplification factors or local helicity have found growing adoption, but their reliability is limited by their level of simplification and the sparsity of data they are based on. To address this gap, a new study of the effect of distributed random surface roughness on crossflow transition was started, using a 30◦ swept plate as a test vehicle. The plate, which is mounted in the 2 Foot × 3 Foot Low-Speed Boundary-Layer Channel at the NASA Langley Research Center, has been tested with micro-scale roughness of varied geometry near the leading edge. In order to analyze the flow in the test section and prepare Direct Numerical Simulations with resolved roughness, Reynolds-Averaged Navier-Stokes (RANS) simulations of the flow over the plate installed in the wind tunnel test section were performed using the Langtry-Menter transition model with stationary crossflow extension. This paper presents the results of these simulations, with comparison to the experiment.

ARMD↗

Subharmonic Route to Boundary-Layer Transition - Critical Layer Nonlinearity

The linear and nonlinear dynamics of a triad of initially linear stability waves comprising a single plane wave at fundamental frequency and two symmetric oblique waves with half the frequency and streamwise wave number of the plane wave are presented. Analysis is performed for the initial nonlinear development of the waves where the order of the oblique waves' amplitude is equal to or less than that of the plane wave. Results show that the fundamental basically follows the linear theory, while the subharmonic follows an exponential-of-an-exponential growth.

Mankbadi, Reda R.↗

Interface stability in a slowly rotating low-gravity tank Theory

The equilibrium configuration of a bubble in a rotating liquid confined by flat axial boundaries (baffles) is found. The maximum baffle spacing assuring bubble confinement is bounded from above by the natural length of a bubble in an infinite medium under the same conditions. Effects of nonzero contact angle are minimal. The problem of dynamic stability is posed. It can be solved in the limit of rapid rotation, for which the bubble is a long cylinder. Instability is to axisymmetric perturbations; nonaxisymmetric perturbations are stable. The stability criterion agrees with earlier results.

Gans, R. F.↗

Nonlinear development and secondary instability of large-amplitude Goertler vortices in hypersonic boundary layers

The weakly nonlinear development of large-amplitude Goertler vortices in a small neighborhood of the neutral stability position given by the linear theory is considered. It is shown that if the basic state and the wall curvature satisfy a certain condition, the Goertler vortices will grow until they become so large as to drive a mean flow correction as large as the basic state. Attention is then given to the further downstream evolution of these large-amplitude vortices, and it is shown that the vortices have a triple layer structure which consists of a region of vortex activity bounded by two transition layers over which the amplitude of the harmonic part of the vortex decays to zero exponentially. The positions of the two transition layers are found to be governed by a free boundary problem, which is solved for a number of curvature distribution cases. The secondary instability of the two transition layers is investigated with respect to traveling waves which are pi/2 out of phase in the spanwise direction with the steady Goertler vortices.

Fu, Y.↗

Strain-relaxed tetragonal MAPbI 3 results in efficient mesoporous solar cells

In mesoporous organo-halide perovskite solar cells, the coexistence of tetragonal phase and cubic phase in the perovskite layer is found to result in strained crystal structure at the mesoporous TiO 2 (mp-TiO 2 )/perovskite heterointerface. This results in structural defects pinholes, grain boundaries, and interfacial porosity which is detrimental towards photovoltaic performance and device stability. The strain at the interface originates from the lattice mismatch between the mp-TiO 2 nanoparticles (NPs) and perovskite. Here, a transformative approach is demonstrated to realize relaxed and high purity tetragonal phase MAPbI 3 (RP) near the mp-TiO 2 interface region. This approach involves inserting a ~2 nm lattice matched buffer layer of cubic CsPbBr 3 between mp-TiO 2 NPs and MAPbI 3 , which serves as a template for epitaxial growth for top MAPbI 3 . As a result, the solar cell with relaxed MAPbI3 shows a power conversion efficiency (PCE) of 22.12% with significantly enhanced environmental, light, and thermal stability.

14 SOLAR ENERGY↗

The effects of curvature and viscosity on baroclinic instability: A two-layer model

A linear stability analysis of a baroclinic zonal current contained between two parallel rigid boundaries is presented. Curvature is included by performing the analysis on a beta b-plane and viscosity by allowing for the effects of Ekman layers on the rigid boundaries. A two-layer model is used. This calculation was carried out to assist in the design of a spherical model of the general circulation of the earth's atmosphere for Spacelab. In the low-gravity environment on an orbiting vehicle, a dominant radial dielectric body force, analogous to planetary gravity, can be achieved over a volume of liquid held between two concentric spheres. The results show the Eady short wavelength cutoff, and long wavelength cutoffs due to Ekman damping and curvature.

Fowlis, W. W.↗

Numerical method for gas dynamics combining characteristic and conservation concepts

An efficient implicit numerical method that solves the compressible Navier-Stokes equations in arbitrary curvilinear coordinates by the finite-volume technique is presented. An intrinsically dissipative difference scheme and a fully implicit treatment of boundary conditions, based on characteristic and conservation concepts, are used to improve stability and accuracy. Efficiency is achieved by using a diagonal form of the implicit algorithm and spatially varying time-steps. Comparisons of various schemes and methods are presented for one- and two-dimensional flows, including transonic separated flow past a thick circular-arc airfoil in a channel. The new method is equal to or better than a version of MacCormack's hybrid method in accuracy and it converges to a steady state up to an order of magnitude faster.

Coakley, T. J.↗

Stability of zonal flows on Jupiter

Gierasch et al. (1979) have shown that the question of the stability of zonal flows on Jupiter is influenced strongly by the lower boundary condition. In the case of Jupiter, the troposphere rests above a deep layer of fluid which is probably extremely close to adiabatic. The present study extends the earlier investigations to a wider range of conditions, taking into account the nongeostrophic limit and the presence of horizontal shear. The major difference between the behavior of the model derived in the study and the behavior of a model with a rigid lower boundary, corresponding to the terrestrial conditions, is that baroclinic instability by the Eady mechanism is suppressed. The influence of beta (Charney-type instability) operates only at the top of the troposphere, not at its base. The suppression of the Eady modes seems to apply all the way down to the nongeostrophic limit where the Rossby number is order unity.

Conrath, B. J.↗