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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Design and performance of the University of Michigan 6.6-inch hypersonic wind tunnel

The tunnel described has several design features intended to maintain laminar flow in the boundary layer of its nozzle. Measurements show that transition to turbulence in the nozzle wall boundary layer begins at the throat and is sensitive to surface roughness, heat transfer rate, and longitudinal radius of curvature. The observed dependence of transition on heat transfer rate is the reverse of that predicted by stability theory for infinitesimal disturbances. Tests include boundary layer surveys of a contoured nozzle and a conical nozzle with four interchangeable throats.

Amick, J. L.↗
DOE OSTI2024

Phase Pattern Formation in Grain Boundaries

Here, in this Letter, we derive conditions that predict the existence of two-phase periodic-pattern grain boundary structures that are stable against coarsening. While previous research has established that elastic effects can lead to phase pattern formation on crystal surfaces, the possibility of stable grain boundary structures composed of alternating grain boundary phases has not been previously analyzed. Our theory identifies the specific combination of grain boundary and materials properties that enable the emergence of patterned grain boundary states and shows that the dislocation content of grain boundary phase junctions, absent in surface phenomena, weakens the stability of the patterned structures. The predictions of the theory are tested using a model copper grain boundary that exhibits multiple phases and two-phase pattern formation. We discuss how, similarly to surfaces, elastic effects associated with grain boundary phase junctions have profound implications for how grain boundary phases transform.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Numerical algorithms for computations of feedback laws arising in control of flexible systems

Several continuous models will be examined, which describe flexible structures with boundary or point control/observation. Issues related to the computation of feedback laws are examined (particularly stabilizing feedbacks) with sensors and actuators located either on the boundary or at specific point locations of the structure. One of the main difficulties is due to the great sensitivity of the system (hyperbolic systems with unbounded control actions), with respect to perturbations caused either by uncertainty of the model or by the errors introduced in implementing numerical algorithms. Thus, special care must be taken in the choice of the appropriate numerical schemes which eventually lead to implementable finite dimensional solutions. Finite dimensional algorithms are constructed on a basis of a priority analysis of the properties of the original, continuous (infinite diversional) systems with the following criteria in mind: (1) convergence and stability of the algorithms and (2) robustness (reasonable insensitivity with respect to the unknown parameters of the systems). Examples with mixed finite element methods and spectral methods are provided.

Lasiecka, Irena↗
DOE OSTI2025

Analysis of the weighted shifted boundary method for the Poisson and Stokes problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. Accuracy is maintained by properly shifting the location and values of the boundary conditions. This avoids integration over cut cells and the associated implementation issues. Recently, the Weighted SBM (WSBM) was proposed for the Navier-Stokes equations with free surfaces and the Stokes flow with moving boundaries. The attribute “weighted” in the name WSBM stems from the fact that its variational form is weighted with the elemental volume fraction of active fluid. The motivation for the development of the WSBM was the preservation of the volume of active fluid to a higher degree of accuracy, which in turn resulted in improved stability and robustness characteristics in moving-boundary, time-dependent simulations. In this article, we present the numerical analysis of the WSBM formulations for the Poisson and Stokes problems. We give mathematical conditions under which the bilinear forms defining the discrete variational formulations are uniformly coercive (Poisson problem) or inf-sup stable (Stokes problem). By these results, stability and optimal convergence is proven in the natural norm; L2-error estimates can also be derived.

Approximate domain boundaries↗
DOE OSTI2022

Behavior and mechanisms of Doppler wind lidar error in varying stability regimes

Abstract. Wind lidars are widespread and important tools in atmospheric observations. An intrinsic part of lidar measurement error is due to atmospheric variability in the remote-sensing scan volume. This study describes and quantifies the distribution of measurement error due to turbulence in varying atmospheric stability. While the lidar error model is general, we demonstrate the approach using large ensembles of virtual WindCube V2 lidar performing a profiling Doppler-beam-swinging scan in quasi-stationary large-eddy simulations (LESs) of convective and stable boundary layers. Error trends vary with the stability regime, time averaging of results, and observation height. A systematic analysis of the observation error explains dominant mechanisms and supports the findings of the empirical results. Treating the error under a random variable framework allows for informed predictions about the effect of different configurations or conditions on lidar performance. Convective conditions are most prone to large errors (up to 1.5 m s−1 in 1 Hz wind speed in strong convection), driven by the large vertical velocity variances in convective conditions and the high elevation angle of the scanning beams (62∘). Range-gate weighting induces a negative bias into the horizontal wind speeds near the surface shear layer (−0.2 m s−1 in the stable test case). Errors in the horizontal wind speed and direction computed from the wind components are sensitive to the background wind speed but have negligible dependence on the relative orientation of the instrument. Especially during low winds and in the presence of large errors in the horizontal velocity estimates, the reported wind speed is subject to a systematic positive bias (up to 0.4 m s−1 in 1 Hz measurements in strong convection). Vector time-averaged measurements can improve the behavior of the error distributions (reducing the 10 min wind speed error standard deviation to <0.3 m s−1 and the bias to <0.1 m s−1 in strong convection) with a predictable effectiveness related to the number of decorrelated samples in the time window. Hybrid schemes weighting the 10 min scalar- and vector-averaged lidar measurements are shown to be effective at reducing the wind speed biases compared to cup measurements in most of the simulated conditions, with time averages longer than 10 min recommended for best use in some unstable conditions. The approach in decomposing the error mechanisms with the help of the LES flow field could be extended to more complex measurement scenarios and scans.

17 WIND ENERGY↗
NASA NTRSDate not supplied

Laminar to Turbulence Transition in Boundary Layers due to Tripping Devices

Direct numerical simulations (DNS)and linear stability analysis were performed to identify the mechanisms by which two-and three-dimensional trips promote early transition in subsonic boundary layers over a flat plate. Investigations were conducted for two boundary layer transition trip configurations: one is a two-dimensional rectangular-shaped rod and the other is a three-dimensional zigzag tape. Reynolds numbers based on the trip heights are 683 and 300 for the two-and three-dimensional trips, respectively. Simulations with the two-dimensional trip showed that the flow separates at the top of the trip and forms a long recirculation zone downstream. However, the flow remains steady in the simulation without any external disturbances. Linear stability analysis showed that the inflectional profile in the recirculation zone is strongly unstable to frequencies in the range of 100-300 Hz. Simulations with the external disturbances forced at the frequency of 200 Hz showed that disturbances grow strongly inside the recirculation zone as observed in the Klebanoff1 experiment.

Transition↗
NASA NTRS1990

The inviscid stability of supersonic flow past a sharp cone

The laminar boundary layer which forms on a sharp cone in a supersonic freestream, where lateral curvature plays a key role in the physics of the problem is considered. This flow is then analyzed from the point of view of linear, temporal, inviscid stability. The basic, non-axisymmetric disturbance equations are derived for general flows of this class, and a so called triply generalized inflexion condition is found for the existence of subsonic neutral modes of instability. This condition is analogous to the well-known generalized inflexion condition found in planar flows, although in the present case the condition depends on both axial and aximuthal wavenumbers. Extensive numerical results are presented for the stability problem at a freestream Mach number of 3.8, for a range of streamwise locations. These results reveal that a new mode of instability may occur, peculiar to flows of this type involving curvature. Additionally, asymptotic analyses valid close to the tip of the cone, far downstream of the cone are presented, and these give a partial (asymptotic) description of this additional mode of instability.

Duck, Peter W.↗
DOE OSTI2024

Uncovering the Network Modifier for Highly Disordered Amorphous Li‐Garnet Glass‐Ceramics

Abstract Highly disordered amorphous Li 7 La 3 Zr 2 O 12 (aLLZO) is a promising class of electrolyte separators and protective layers for hybrid or all‐solid‐state batteries due to its grain‐boundary‐free nature and wide electrochemical stability window. Unlike low‐entropy ionic glasses such as Li x PO y N z (LiPON), these medium‐entropy non‐Zachariasen aLLZO phases offer a higher number of stable structure arrangements over a wide range of tunable synthesis temperatures, providing the potential to tune the LBU‐Li + transport relation. It is revealed that lanthanum is the active “network modifier” for this new class of highly disordered Li + conductors, whereas zirconium and lithium serve as “network formers”. Specifically, within the solubility limit of La in aLLZO, increasing the La concentration can result in longer bond distances between the first nearest neighbors of Zr─O and La─O within the same local building unit (LBU) and the second nearest neighbors of Zr─La across two adjacent network‐former and network‐modifier LBUs, suggesting a more disordered medium‐ and long‐range order structure in LLZO. These findings open new avenues for future designs of amorphous Li + electrolytes and the selection of network‐modifier dopants. Moreover, the wide yet relatively low synthesis temperatures of these glass‐ceramics make them attractive candidates for low‐cost and more sustainable hybrid‐ or all‐solid‐state batteries for energy storage.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Crack‐Free Single‐Crystalline LiNiO 2 for High Energy Density All‐Solid‐State Batteries

Abstract Single‐crystalline layered oxide (LiNi 1‐ x ‐ y Mn x Co y O 2 ) cathodes have been found to exhibit exceptional electrochemical properties when coupled with various inorganic solid electrolytes (ISEs) in all‐solid‐state batteries (ASSBs). Their advantages stem from the robust morphological integrity with the absence of grain boundaries and the high electrochemical oxidative stability. Here, ASSBs featuring single‐crystalline LiNiO 2 (LNO) with the highest Ni content are reported, offering a high theoretical specific capacity of 275 mAh g ‐1 alongside a high average discharge voltage (3.7 V vs Li + /Li). Through a careful investigation, it is demonstrated that micron‐sized single‐crystalline LNO (µSC‐LNO) composite cathodes with a halide ISE exhibit a high initial discharge capacity of 205 mAh g ‐1 with an outstanding cycle performance over 200 cycles in room‐temperature ASSBs. The significance of engineering parameters is emphasized, such as particle size and specific density, in promoting a homogeneous and fast Li + transport within the composite cathodes. Furthermore, the formation of undesirable interphase between the halide ISE in the cathode and sulfide ISE separator is elucidated, which may be a critical factor impeding long‐term cyclability of ASSBs. This work provides insights into the design of composite cathodes for high‐energy‐density ASSBs.

Chemistry↗
DOE OSTI2021

Considering computational speed vs. accuracy: Choosing appropriate mesoscale RVE boundary conditions

Modeling a material’s microstructure using continuum theories allows for inspection of the relationship between coarse scale and fine scale behaviors. Computational limits generally require selection of a sub-volume from a bulk sample in order to directly model the microstructure. Boundary conditions are applied to the sub-volume to mimic the excluded bulk material. Appropriate selection of boundary conditions helps effectively determine the appropriate spatial scale required of the sub-volume. Applicable boundary conditions include direct displacement, periodic, and uniform traction. While direct displacement and periodic boundary conditions are commonly used, uniform traction boundary conditions have seen limited use due to rigid body stability issues in simulations of compression or shear deformation. A new application of uniform traction boundary conditions was developed through linear constraint equations, similar to approaches employed by direct displacement and periodic boundary conditions, to quench rigid body motions with minimal interference of the relative deformation of the model. These boundary conditions were tested by compressing several synthetically generated periodic microstructures using the finite element method. Evaluating the effective stiffness along the compression axis, the direct displacement boundary condition produced the stiffest response, whereas the uniform traction boundary condition produced the most compliant. Periodic boundary conditions produced the same response for all volumes analyzed and both the direct displacement and uniform traction boundary conditions trended toward the periodic response as the domain volume increased. Computational performance was also evaluated for each boundary condition using implicit and explicit solvers. Direct displacement boundary conditions presented the lowest computational cost of all of the boundary conditions followed by periodic then uniform traction. The computational expense of periodic and uniform traction boundary conditions limited the viable spatial scale and mesh resolutions able to be simulated. Selection of appropriate boundary conditions for specific uses need to be a balance between allowable computational expense and accuracy of the method. Techniques for evaluating which boundary conditions to use are discussed.

42 ENGINEERING↗
DOE OSTI2025

X-ray Diffraction Studies of Single-Crystal Materials for Broad Battery Applications

Single-crystal materials have attracted growing interest in battery research due to their well-defined crystallographic orientation, absence of grain boundaries, and enhanced mechanical and electrochemical stability. This Review provides a comprehensive overview of recent advances in the synthesis, structural evolution, and performance optimization of single-crystal electrodes and solid electrolytes. Particular focus is placed on the application of advanced X-ray diffraction (XRD) techniques, including operando synchrotron diffraction, reciprocal space mapping, and Bragg coherent diffraction imaging, which have enabled in-depth investigations of lattice strain, cation disorder, phase transitions, and defect formation. Representative case studies across Ni-rich layered oxides, spinel-type cathodes, and garnet-based electrolytes are examined to highlight the structural features unique to single crystals. Additionally, the synergistic integration of XRD with machine learning, tomography, and spectroscopy is discussed as a powerful direction for real-time analysis and predictive modeling. Furthermore, these insights provide critical guidance for the rational design of high-performance single-crystal materials in lithium, sodium, and solid-state battery systems.

25 ENERGY STORAGE↗
DOE OSTI2024

Heuristic predictions of RMP configurations for ELM suppression in ITER burning plasmas and their impact on divertor performance

Abstract A subspace of resonant magnetic perturbation (RMP) configurations for edge localized mode (ELM) suppression is predicted for H-mode burning plasmas at 15 MA current and 5.3 T magnetic field in ITER. Perturbations to the core plasma can be reduced by a factor of 2 for equivalent edge stability proxies, while the perturbed plasma boundary geometry remains mostly resilient. The striated domain of perturbed field lines connecting from the main plasma (normalized poloidal flux < 1 ) to the divertor targets is found to be significantly larger than the expected heat load width in the absence of RMPs. This facilitates heat load spreading with peak values at an acceptable level below 10 MW m − 2 on the outer target already at moderate gas fueling and low Ne seeding for additional radiative dissipation of the 100 MW of power into the scrape-off layer (SOL). On the inner target, however, re-attachment is predicted away from the equilibrium strike point due to increased upstream heat flux, higher downstream temperature and less efficient impurity radiation.

3D plasma boundary modeling↗
NASA NTRS1975

A cubic spline approximation for problems in fluid mechanics

A cubic spline approximation is presented which is suited for many fluid-mechanics problems. This procedure provides a high degree of accuracy, even with a nonuniform mesh, and leads to an accurate treatment of derivative boundary conditions. The truncation errors and stability limitations of several implicit and explicit integration schemes are presented. For two-dimensional flows, a spline-alternating-direction-implicit method is evaluated. The spline procedure is assessed, and results are presented for the one-dimensional nonlinear Burgers' equation, as well as the two-dimensional diffusion equation and the vorticity-stream function system describing the viscous flow in a driven cavity. Comparisons are made with analytic solutions for the first two problems and with finite-difference calculations for the cavity flow.

Rubin, S. G.↗
NASA NTRS1975

Viscous flow solutions with a cubic spline approximation

A cubic spline approximation is used for the solution of several problems in fluid mechanics. This procedure provides a high degree of accuracy even with a nonuniform mesh, and leads to a more accurate treatment of derivative boundary conditions. The truncation errors and stability limitations of several typical integration schemes are presented. For two-dimensional flows a spline-alternating-direction-implicit (SADI) method is evaluated. The spline procedure is assessed and results are presented for the one-dimensional nonlinear Burgers' equation, as well as the two-dimensional diffusion equation and the vorticity-stream function system describing the viscous flow in a driven cavity. Comparisons are made with analytic solutions for the first two problems and with finite-difference calculations for the cavity flow.

Rubin, S. G.↗
NASA NTRS1980

Numerical studies of new stellarator concepts

A three-dimensional computer code has been developed to study the magnetohydrodynamic equilibrium and stability of a diffuse or sharp boundary plasma in toroidal geometry. It is shown how equilibria with net toroidal current identically zero can be determined and how growth rates of instabilities can be calculated. Applications are made to an l = 2, 3 stellarator configuration that offers the possibility of achieving a critical value as high as 10% for the plasma parameter

Bauer, F.↗
NASA NTRS1982

Finite-difference solution of the compressible stability eigenvalue problem

A compressible stability analysis computer code is developed. The code uses a matrix finite difference method for local eigenvalue solution when a good guess for the eigenvalue is available and is significantly more computationally efficient than the commonly used initial value approach. The local eigenvalue search procedure also results in eigenfunctions and, at little extra work, group velocities. A globally convergent eigenvalue procedure is also developed which may be used when no guess for the eigenvalue is available. The global problem is formulated in such a way that no unstable spurious modes appear so that the method is suitable for use in a black box stability code. Sample stability calculations are presented for the boundary layer profiles of a Laminar Flow Control (LFC) swept wing.

Malik, M. R.↗