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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Unsteady Aerodynamic Response of a Linear Cascade of Airfoils in Separated Flow

The overall objective of this research program was to investigate methods to modify the leading edge separation region, which could lead to an improvement in aeroelastic stability of advanced airfoil designs. The airfoil section used is representative of current low aspect ratio fan blade tip sections. The experimental potion of this study investigated separated zone boundary layer from removal through suction slots. Suction applied to a cavity in the vicinity of the separation onset point was found to be the most effective location. The computational study looked into the influence of front camber on flutter stability. To assess the influence of the change in airfoil shape on stability the work-per-cycle was evaluated for torsion mode oscillations. It was shown that the front camberline shape can be an important factor for stabilizing the predicted work-per-cycle and reducing the predicted extent of the separation zone. In addition, data analysis procedures are discussed for reducing data acquired in experiments that involve periodic unsteady data. This work was conducted in support of experiments being conducted in the NASA Glenn Research Center Transonic Flutter Cascade. The spectral block averaging method is presented. This method is shown to be able to account for variations in airfoil oscillation frequency that can occur in experiments that force oscillate the airfoils to simulate flutter.

Capece, Vincent R.↗
DOE OSTI2024

SEI Formation and Lithium-Ion Electrodeposition Dynamics in Lithium Metal Batteries via First-Principles Kinetic Monte Carlo Modeling

The stabilization and enhanced performance of lithium metal batteries (LMBs) depend on the formation and evolution of the Solid Electrolyte Interphase (SEI) layer as a critical component for regulating the Li metal electrodeposition processes. This study employs a first-principles kinetic Monte Carlo (kMC) model to simulate the SEI formation and Li + electrodeposition processes on a lithium metal anode, integrating both the electrochemical electrolyte reduction reactions and the diffusion events giving place to the SEI aggregation processes during battery charge and discharge processes. The model replicates the competitive interactions between organic and inorganic SEI components, emphasizing the influence of the cycling regime. Results indicate that grain boundaries within the SEI facilitate faster lithium-ion transport compared to crystalline regions, crucial for improving the performance and stability of LMBs. The findings underscore the importance of dynamic SEI modeling for further development of next-generation high-energy-density batteries.

25 ENERGY STORAGE↗
NASA NTRS1993

A lag model for turbulent boundary layers developing over rough bleed surfaces

Boundary layer mass removal (bleed) through spanwise bands of holes on a surface is used to prevent or control separation and to stabilize the normal shock in supersonic inlets. The addition of a transport equation lag relationship for eddy viscosity to the rough wall algebraic turbulence model of Cebeci and Chang was found to improve agreement between predicted and measured mean velocity distributions downstream of a bleed band. The model was demonstrated for a range of bleed configurations, bleed rates, and local free stream Mach numbers. In addition, the model was applied to the boundary layer development over acoustic lining materials for the inlets and nozzles of commercial aircraft. The model was found to yield accurate results for integral boundary layer properties unless there was a strong adverse pressure gradient.

Lee, J.↗
NASA NTRSDate not supplied

Integrated GNC Modifications and Performance Assessments for the Mars Ascent Vehicle, Spin-Stabilized Upper Stage Configuration

The objective of the Mars Ascent Vehicle (MAV), an element of the Mars Sample Return (MSR) campaign, is to successfully transport Martian surface samples from the Martian surface to a stable Low-Martian Orbit. Complete autonomy is required throughout ascent, and orbital insertion is constrained by tight dispersion boundaries. Given mass requirements to merely transport the MAV onto the Martian surface, reducing the total mass of the MAV reduces stress on all other campaign elements. An unguided, spin-stabilized second stage MAV allows for lower overall vehicle mass, at the cost of reduced GNC capability. This paper presents GNC modifications developed and employed to improve the performance of an unguided, spin-stabilized second stage MAV configuration. These updates include a predictive calculation scheme to solve the spin-stabilized attitude necessary for successful orbital insertion, and a spin-stabilized control law for a preliminary low-rate first stage spin-up. Trades were conducted focusing on first and second stage low-rate and high-rate spin rates (and other factors) in MSFC's MANTIS toolkit through 6-degree-of-freedom simulation and Monte Carlo analysis, and are presented here.

Dane Erickson↗
DOE OSTI2024

Development of a neural network model for peeling–ballooning stability analysis in the KSTAR tokamak pedestals

The neural network model, MISHKA-NN is developed to mitigate the computational burden associated with the linear ideal magnetohydrodynamic (MHD) stability analysis of the pedestal based on the peeling–ballooning (P–B) model. By utilizing both 1D plasma profiles (current density, pressure gradient, and safety factor) and 0D parameters (plasma geometry, total current, and toroidal mode number), the model predicts linear growth rate of edge-localized ideal MHD instability in a given equilibrium state. By enabling the prediction of each instability within a second, the model reduces the time required for plotting a pedestal P–B stability diagram (j - α diagram) from approximately 100 CPU hours to a few CPU minutes. Notably, even with the utilization of parametric pressure and current profiles and plasma boundary shapes for the training dataset, the model shows a satisfactory level of performance in benchmarking the j - α diagram for the reconstructed equilibrium from a KSTAR tokamak experiment. We anticipate the model to serve as a versatile alternative to 2D linear MHD stability codes, alleviating numerical costs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

High-Fidelity, Low-Dissipation/Symmetry-Preserving Numerical Scheme for Solving the Euler Equations with Unstructured, Metric-Based Mesh Adaptation

This work presents an overview of a high-fidelity compressible Euler solver that utilizes the continuous Galerkin (CG) method with added artificial numerical diffusion for stabilization to solve a variety of unsteady and steady benchmark inviscid flow problems. This work shows that discretizing the Euler equations with this CG approach and first order basis functions produces a cost-effective stencil as well as simple well-posed boundary conditions. We show through convergence testing with manufactured solutions that the reduced stencil of CG, combined with the low amount of artificial diffusion required when using the stabilization method outlined in this work, leads to stable and highly accurate results for a variety of unsteady and steady applications. When combined with the adaptive mesh refinement approach used for many of the cases in this work, our results show that the flow solver achieves even more accurate results. A variety of inviscid flow cases are presented in this work, including transient 2D cases with complex shock structures and several steady 3D airfoils sections with a constant span.

Doetsch, Kevin [ORNL] (ORCID:0000000267051705)↗
NASA NTRS1994

Physical description of boundary-layer transition: Experimental evidence

The problems of understanding the origins of turbulent flow and transition to turbulent flow are the most important unsolved problems of fluid mechanics and aerodynamics. It is well known that the stability, transition, and turbulent characteristics of bounded shear layers are fundamentally different from those of free shear layers. Likewise, the stability, transition, and turbulent characteristics of open systems are fundamentally different from those of closed systems. Because of the influence of indigenous disturbances, surface geometry and roughness, sound, heat transfer, and ablation, it is not possible to develop general prediction schemes for transition location and the nature of turbulent structures in boundary-layer flows. At the present time no mathematical model exists that can predict the transition Reynolds number on a flat plate. The recent progress in this area is encouraging, in that a number of distinct transition mechanisms have been found experimentally. The theoretical work finds them to be amplitude and Reynolds-number dependent. The theory remains rather incomplete with regard to predicting transition. Amplitude and spectral characteristics of the disturbances inside the laminar viscous layer strongly influence which type of transition occurs. The major need in this area is to understand how freestream disturbances are entrained into the boundary layer, i.e., to answer the question of receptivity. We refer receptivity to the mechanism(s) that cause freestream disturbances to enter the boundary layer and create the initial amplitudes for unstable waves.

Saric, William S.↗
DOE OSTI2022

Innovative Method for Welding in Generation 3 CSP to Enable Reliable Manufacturing of Solar Receivers to withstand Daily Cycling at Temperatures Above 700°C (Final Technical Report)

Inconel® Alloy 740H® (alloy 740H) was the first age-hardenable nickel-based alloy approved by the ASME Boiler & Pressure Vessel Code for use in pressure-boundary applications. Over the past ~20 years the alloy has been optimized for weldability and high-temperature stability, approved for use in different applications. Development of a supply chain combined with the advantageous properties of the alloy (high-temperature creep strength, oxidation and corrosion resistance, etc.) have resulted in the alloy being applied to new high-temperature power cycle demonstration projects, and of particular interest are applications to concentrating solar power (CSP) to enable higher-efficiency Generation 3 CSP systems and the corresponding supercritical CO 2 (sCO 2 ) power cycle components (heat exchangers, piping, etc.). The high allowable stresses of alloy 740H also make it a desirable material for current Generation 2 CSP solar power receivers to improve cyclic capability and/or reduce receiver height. Recent experiences in demonstration projects utilizing alloy 740H identified cracking issues during welding and fabrication. In this project, a detailed study was done to confirm and clarify the Stress Relaxation Cracking (SRxC) mechanism, also known as stress relief cracking or strain-age cracking (SAC), during post-weld heat-treatment (PWHT). This involved detailed microscopy and advanced characterization to understand the root cause(s) of three failures obtained from industry. Based in-part on these findings, a targeted laboratory based SRxC test method was utilized to evaluate variables such as heat-to-heat variations, strain level, PWHT temperature, and starting material condition on three heats of alloy 740H. Industrial shop welding of cold-worked plates was also conducted. The research showed the following: SRxC was confirmed as the cracking mechanism for all field failures; Stress state (from residual stresses, constraint, deformation, and local stress concentrations) was playing a significant role in field failures and laboratory testing confirmed increasing susceptibility for all heats with increasing strain levels. High levels of microstructural strain were identified at crack initiation locations, in some cases leading to local recrystallization; Precipitate free zones (PFZs) at grain boundaries were found at relaxation cracks and crack initiation locations uniquely associated with SRxC in alloy 740H. Laboratory testing reproduced this microstructural feature which had only previously been reported in long-term creep testing of weldments. Advanced nano-scale characterization confirmed the presence of a moving boundary leading to coarsening of precipitates and PFZs where damage accumulated; The research suggested heat-to-heat variations due to local chemistry and processing may influence SRxC susceptibility, but more work is needed to fully clarify these effects. To disseminate the key learnings from this research to the scientific and engineering communities and alloy 740H end users, multiple technical publications and presentations were made, an industrial alloy 740H users meeting was held, and a new industry guideline specification document which can be directly implemented by end-users of alloy 740H was produced.

14 SOLAR ENERGY↗
DOE OSTI2022

On the closedness and geometry of tensor network state sets

Tensor network states (TNS) are a powerful approach for the study of strongly correlated quantum matter. The curse of dimensionality is addressed by parametrizing the many-body state in terms of a network of partially contracted tensors. These tensors form a substantially reduced set of effective degrees of freedom. In practical algorithms, functionals like energy expectation values or overlaps are optimized over certain sets of TNS. Concerning algorithmic stability, it is important whether the considered sets are closed because, otherwise, the algorithms may approach a boundary point that is outside the TNS set and tensor elements diverge. Here we discuss the closedness and geometries of TNS sets, and we propose regularizations for optimization problems on non-closed TNS sets. We show that sets of matrix product states (MPS) with open boundary conditions, tree tensor network states, and the multiscale entanglement renormalization ansatz are always closed, whereas sets of translation-invariant MPS with periodic boundary conditions (PBC), heterogeneous MPS with PBC, and projected entangled pair states are generally not closed. The latter is done using explicit examples like the W state, states that we call two-domain states, and fine-grained versions thereof.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Combined plasma–coil optimization algorithms

Combined plasma–coil optimization approaches for designing stellarators are discussed and a new method for calculating free-boundary equilibria for multiregion relaxed magnetohydrodynmics (MRxMHD) is proposed. Four distinct categories of stellarator optimization, two of which are novel approaches, are the fixed-boundary optimization, the generalized fixed-boundary optimization, the quasi-free-boundary optimization, and the free-boundary (coil) optimization. These are described using the MRxMHD energy functional, the Biot–Savart integral, the coil-penalty functional and the virtual casing integral and their derivatives. The proposed free-boundary equilibrium calculation differs from existing methods in how the boundary-value problem is posed, and for the new approach it seems that there is not an associated energy minimization principle because a non-symmetric functional arises. We propose to solve the weak formulation of this problem using a spectral-Galerkin method, and this will reduce the free-boundary equilibrium calculation to something comparable to a fixed-boundary calculation. In our discussion of combined plasma–coil optimization algorithms, we emphasize the importance of the stability matrix.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Emergent Ferromagnetism in CaRuO 3 /CaMnO 3 (111)-Oriented Superlattices

The boundary between CaRuO 3 and CaMnO 3 is an ideal test bed for emergent magnetic ground states stabilized through interfacial electron interactions. In this system, nominally antiferromagnetic and paramagnetic materials combine to yield interfacial ferromagnetism in CaMnO 3 due to electron leakage across the interface. In this work, we show that the crystal symmetry at the surface is a critical factor determining the nature of the interfacial interactions. Specifically, by growing CaRuO 3 /CaMnO 3 heterostructures along the (111) instead of the (001) crystallographic axis, we achieve a 3-fold enhancement of the magnetization and involve the CaRuO 3 layers in the ferromagnetism, which now spans both constituent materials. Importantly, the stabilization of a net magnetic moment in CaRuO 3 through strain effects has been long-sought but never consistently achieved, and our observations demonstrate the importance of interface engineering in the development of new functional heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Plateau–Rayleigh instability with a grain boundary twist

We demonstrate using theoretical, computational, and experimental studies a morphological instability, in which a polycrystalline nanorod breaks up at grain boundaries (GBs) into an array of isolated domains. Our theoretical model is used to establish a neutral stability surface demarcating stable and unstable perturbations. It is shown that GBs play a destabilizing role in which the critical wavelength for the instability decreases with the increase in the GB energy. We carry out phase field simulations, which reveal accelerated pinch-off kinetics with the increase in the GB energy and predict temporal evolution of interfacial profiles in quantitative agreement with experimental observations.

Hussein, Omar↗
DOE OSTI2023

Asymptotic vacuum solution at tokamak X-point tip

In the H-mode regime of diverted tokamaks, the presence of strong pressure gradients in the pedestal gives rise to a sizable bootstrap current, together with the Ohmic and Pfirsch–Schlueter currents, close to the separatrix. For such equilibria, the presence of finite current density close to the separatrix requires the reexamination of equilibrium properties. It is almost universally assumed that the two branches of the separatrix (the stable and unstable manifolds) are straight as they cross at the X-point. However, the opposite angles of the plasma-filled segment and vacuum one cannot be equal if the current density does not vanish at the separatrix on the plasma side. We solve this difficulty by chipping off a thin layer of plasma edge so that the sharp corner of the plasma-filled segment becomes a hyperbola. Using the conformal transformation, we found that in the assumption of a hyperbolic boundary, the X point moves beyond the plasma boundary to fall in the vacuum region. An acute angle of the plasma-filled segment leads to an obtuse opposite angle of vacuum segment and vice versa. In the case of an acute angle of the plasma-filled segment, the new X point shifts inside the X point formed by the asymptotes of a hyperbolic boundary; in the case of an obtuse angle of the plasma-filled segment, the new X point shifts outside the X point formed by the asymptotes of a hyperbolic plasma boundary. Furthermore, the results are important for understanding the X point features, which affect the tokamak edge stability and transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1984

Natural laminar flow experiments on modern airplane surfaces

Flight and wind-tunnel natural laminar flow experiments have been conducted on various lifting and nonlifting surfaces of several airplanes at unit Reynolds numbers between 0.63 x 10 to the 6th power/ft and 3.08 x 10 to the 6th power/ft, at Mach numbers from 0.1 to 0.7, and at lifting surface leading-edge sweep angles from 0 deg to 63 deg. The airplanes tested were selected to provide relatively stiff skin conditions, free from significant roughness and waviness, on smooth modern production-type airframes. The observed transition locations typically occurred downstream of the measured or calculated pressure peak locations for the test conditions involved. No discernible effects on transition due to surface waviness were observed on any of the surfaces tested. None of the measured heights of surface waviness exceeded the empirically predicted allowable surface waviness. Experimental results consistent with spanwise contamination criteria were observed. Large changes in flight-measured performance and stability and control resulted from loss of laminar flow by forced transition. Rain effects on the laminar boundary layer caused stick-fixed nose-down pitch-trim changes in two of the airplanes tested. No effect on transition was observed for flight through low-altitude liquid-phase clouds. These observations indicate the importance of fixed-transition tests as a standard flight testing procedure for modern smooth airframes.

Holmes, B. J.↗
NASA NTRS1987

Reynolds stress closure in jet flows using wave models

Numerical methods were developed that will form the computational part of the turbulence closure scheme. A wave model was developed for the two-dimensional shear layer. This configuration is being used as a test case for the closure schemes. Various numerical schemes were examined to give efficient solutions of the Rayleigh equation for this geometry. These include both spectral and finite difference methods. Secondly, numerical methods are under development to solve the non-separable Rayleigh equation. This solution is required for the closure scheme in more complex geometries. A model problem was used to assist in the algorithm development. Two-dimensional spectral methods and a hybrid spectral/finite difference technique were developed. An analytic solution of the Rayleigh equation for a basic elliptic flow was obtained. This will be used to verify the stability codes developed for arbitrary geometries. Other numerical methods for solving the Rayleigh equation based on the boundary element technique were also examined. These solutions are forming the basis of a model for the shock structure in jets of arbitrary geometry.

Morris, P. J.↗
NASA NTRS1986

Reynolds stress closure in jet flows using wave models

A wave model was developed for the two dimensional shear layer. This configuration is being used as a test case for the closure schemes. Numerical methods are under development to solve the nonseparable Rayleigh equation. A model problem is being used to assist in the algorithm development. An analytic solution of the Rayleigh equation for a basic elliptic flow was obtained. This will be used to verify the stability codes developed for arbitrary geometries. Other numerical methods for solving the Rayleigh equation based on the boundary element technique are being examined.

Morris, P. J.↗
NASA NTRS1987

The method of projected characteristics for the evolution of magnetic arches

A numerical method of solving fully nonlinear MHD equation is described. In particular, the formulation based on the newly developed method of projected characteristics (Nakagawa, 1981) suitable to study the evolution of magnetic arches due to motions of their foot-points is presented. The final formulation is given in the form of difference equations; therefore, the analysis of numerical stability is also presented. Further, the most important derivation of physically self-consistent, time-dependent boundary conditions (i.e. the evolving boundary equations) is given in detail, and some results obtained with such boundary equations are reported.

Nakagawa, Y.↗
NASA NTRS1987

Roll-up of vorticity in adverse-pressure-gradient boundary layers

It is shown how the unsteady, nonlinear critical-layer equation determines the evolution of instability waves in a weak adverse-pressure-gradient boundary layer. Numerical solutions show that the nonlinearity halts the growth of these inviscidly unstable waves. The stabilizing effect of nonlinearity, in the present case, can be described as a consequence of either the increase (toward zero) of the phase jump across the critical layer or the roll-up of the critical-layer disturbance vorticity.

Goldstein, M. E.↗