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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,369 records · Page 76

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Boundary Corrections for Kernel Approximation to Differential Operators

The kernel-based approach to operator approximation for partial differential equations has been shown to be unconditionally stable for linear PDEs and numerically exhibit unconditional stability for non-linear PDEs. These methods have the same computational cost as an explicit finite difference scheme but can exhibit order reduction at boundaries. In previous work on periodic domains, order reduction was addressed, yielding high-order accuracy. The issue addressed in this work is the elimination of order reduction of the kernel-based approach for a more general set of boundary conditions. Further, we consider the case of both first and second order operators. To demonstrate the theory, we provide not only the mathematical proofs but also experimental results by applying various boundary conditions to different types of equations. The results agree with the theory, demonstrating a systematic path to high order for kernel-based methods on bounded domains.

97 MATHEMATICS AND COMPUTING↗

NASA NTRSDate not supplied

Modeling Boundary-Layer Transition in Subsonic Flow over a Swept Wing

Predicting the onset of boundary-layer transition is often more accurate using physics-based models that directly compute disturbance growth rather than phenomenological models often implemented into industrial CFD codes. The aim of this ongoing study is to calibrate linear, physics-based computations of transition in subsonic flows over swept wings against a large set of experimental data. Advancing the calibration of linear models of transition contributes to the CFD-Vision-2030 goal of automated boundary-layer transition prediction. This progress report uses the dual N-factor method to model transition over the swept NACA 64-2-015A wing. The flow conditions match selected test conditions from an extensive experimental dataset acquired from the NASA Ames 12-ft Pressure Tunnel. The OVERFLOW 2.4b flow solver is used to obtain laminar basic states based on an infinite-span assumption. Stability analyses are performed on 365 distinct configurations with linear stability theory (LST) and parabolized stability equations (PSE) from the Langley Stability and Transition Analysis Codes (LASTRAC), modeling the growth of Tollmien-Schlichting (TS) and stationary crossflow (SCF) disturbances. From a total of 67 data points for unswept, i.e., TS-dominant configurations, the critical N-factor based on PSE is found to be N_TS = 9. The SCF critical N-factor is found to be near 8 for the highly swept, SCF-dominant configurations. Dual N-factor curves for both LST and PSE computations demonstrate a high level of interaction between TS and SCF. It may be worthwhile to investigate an alternate metric to visualize maximal SCF amplification upstream of the transition location to account for the growth of SCF modes near the leading edge, which is not considered in the conventional applications of the dual N-factor criterion.

boundary-layer transition↗

NASA NTRSDate not supplied

Stability Analysis of the Flow over a Swept Forward-Facing Step using PIV Base Flowsin a Non-Orthogonal Coordinate System

Understanding the flow of a crossflow-vortex-dominated boundary layer over a forward-facing step excrescence is necessary in order to mitigate an early transition scenario on the wings and tail of commercial aircraft. By performing BiGlobal stability analysis on the flowfield measured with high-resolution, stereographic Particle Image Velocimetry (PIV), previous work has shown that a family of unstable disturbances exists in the direct downstream vicinity of a supercritical, forward-facing step. The goal of the present paper is to improve the previously used stability approach on two fronts: 1) the stability problem is formulated in a non-orthogonal coordinate system and 2) the contribution of the out-of-plane base-flow derivatives is accounted for locally. The resulting eigen solutions feature significantly stronger growth rates and eigenfunctions that are localized above regions of reverse flow. The changes in the problem formulation furthermore establish a significant improvement in the comparison of the stability solutions with the results from the Spectral Proper Orthogonal Decomposition (SPOD) of a time-resolved measurement of the perturbation content. A large destabilizing effect by the out-of-plane base-flow derivatives is determined to be justified, despite the assumption that these derivatives are small, because the effect can be reconstructed upon using an eigenvalue-correction formula that assumes the responsible terms are infinitesimally small. The unstable perturbation mechanisms are demonstrated to have a convective nature by assessing the relation between the frequency and the out-of-plane wavenumber, which indicates that their group speed does not approach zero. The last 2 facts, that 1) the perturbations are convective and 2) that the out-of-plane base-flow-derivative terms are small, remove any qualitative suspicions that the perturbation problem is not conducive to parabolization. This opens the path to analyzing this problem with a plane-marching approach.

Boundary- Layer Transition↗

NASA NTRSDate not supplied

CFD Simulation of an Inward-Turning Supersonic Inlet Unstart at Flight Mach Number 1.7

Interest in inward-turning supersonic inlets renewed during NASA’s Commercial Supersonic Technology Project, as a means of reducing disturbances caused by propulsion integration. Inward-turning designs have the potential to greatly reduce the strength of shock waves emanating outward from the nacelle, and therefore their contribution to the overall sonic boom signature. For top-mounted, or over-wing installations, disturbances to aircraft aerodynamics are reduced. The nature of the inward-turning geometry results in a boundary-layer thickening effect, opposite that of the familiar thinning of an outward-turning cone boundary-layer. The application of a small amount of boundary-layer bleed at the throat has been shown to alleviate the negative effects of the boundary-layer on subsonic diffuser performance. Boundary-layer bleed, however, comes at the expense of increased complexity, and external disturbances at bleed exits that would partially defeat the purpose of the inward-turning scheme. Vortex generators were proposed as an alternative to bleed, and an initial CFD study of their effectiveness in this application forms the basis of the present work. A source-term model in the Wind-US CFD code was used to model the effects of an array of vortex generators placed at the entrance to the subsonic diffuser. Baseline results with no vortex generators compared favorably to previously reported results with no bleed at back-pressures near the design point, but stability margin in the present case was limited by a curious non-linear phenomena similar to unstart in mixed-compression inlets. While some improvement was noted with vortex generators, their effectiveness was diminished by the limited stability margin. A time-accurate calculation was subsequently performed to better understand this process. Results of this calculation are presented and instantaneous solutions compared with steady-state, subcritical results previously reported. The stability margin limitation is attributed to the use of a constant-pressure outflow boundary condition, which differed from that used in the prior simulations.

Pre-Compressor Cooling, Water Injection, Hypersoni↗

DOE OSTI2020

On the Relationship Between the Marine Cold Air Outbreak M Parameter and Low-Level Cloud Heights in the Midlatitudes

Focusing on conditions of subsidence when low clouds are present, ground-based observations in both the North Atlantic and the Southern Ocean reveal strong relationships between cloud boundary (base and top heights) and different measures of lower tropospheric instability. The difference in potential temperature between the surface and 800 hPa (a metric called M) provides a stronger relationship than measures of inversion strength such as the lower tropospheric stability and estimated inversion strength. This is because (1) inversion strength itself does not correlate well with cloud boundaries, and (2) M contains information that appears important for cloud boundaries. These include the surface forcing through the use of sea surface rather than near-surface air temperature and an upper level close to the real cloud top. These results expand upon previous work on the importance of M as a predictor of cloud morphology. However, important differences are found in low-cloud conditions for the North Atlantic as compared to the Southern Ocean (for a given value of M): stronger inversions, deeper boundary layers, and much larger sea level pressures. Therefore, the relationship between cloud boundaries and M differs between the two regions. A general circulation model provides similar relationships as observed between M and both cloud top height and temperature but tends to place clouds higher and at colder temperatures than observed for a given M. This might cause issues with the representation of precipitation, cloud cover and radiation in the Southern Ocean.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2017

Lightweight ZERODUR: Validation of Mirror Performance and Mirror Modeling Predictions

Upcoming spaceborne missions, both moderate and large in scale, require extreme dimensional stability while relying both upon established lightweight mirror materials, and also upon accurate modeling methods to predict performance under varying boundary conditions. We describe tests, recently performed at NASA's XRCF chambers and laboratories in Huntsville Alabama, during which a 1.2 m diameter, f/1.2988% lightweighted SCHOTT lightweighted ZERODUR(TradeMark) mirror was tested for thermal stability under static loads in steps down to 230K. Test results are compared to model predictions, based upon recently published data on ZERODUR(TradeMark). In addition to monitoring the mirror surface for thermal perturbations in XRCF Thermal Vacuum tests, static load gravity deformations have been measured and compared to model predictions. Also the Modal Response(dynamic disturbance) was measured and compared to model. We will discuss the fabrication approach and optomechanical design of the ZERODUR(TradeMark) mirror substrate by SCHOTT, its optical preparation for test by Arizona Optical Systems (AOS). Summarize the outcome of NASA's XRCF tests and model validations

technology↗

DOE OSTI2023

On the Interfacial Assembly of Anisotropic Amphiphilic Janus Particles

It has been shown both theoretically and experimentally that amphiphilic Janus particles are the most effective solid surfactants to stabilize interfaces. In most cases, the Janus particles investigated have uniform morphologies with Janus boundaries dividing the particle into halves. However, there are many examples of Janus particles where the hydrophilic and hydrophobic domains are not equally distributed. The effects of this uneven domain distribution on the mechanism and kinetics of Janus particle assembly, and final equilibrium state are not well-understood. Dynamic pendant drop tensiometry offers a means to probe both the equilibrium assembly and the kinetics and mechanism of assembly. Here, in this work, the interfacial kinetics and assembly of spherical anisotropic Janus particles are investigated using dynamic pendant drop tensiometry. Systematic studies quantifying the time-dependent interfacial behavior as a function of Janus particle morphology, chemical composition, particle concentration, and NaOH and HCl concentration are performed. These studies shed light on the assembly mechanism of more complex Janus particle morphologies and highlight their effectiveness as interface stabilizers.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Singular Perturbation-Based Large-Signal Order Reduction of Microgrids for Stability and Accuracy Synthesis With Control

The increasing penetration of distributed energy resources (DERs) highlights the growing importance of microgrids (MGs) in enhancing power system reliability. Employing electromagnetic transient (EMT) analysis in MGs becomes crucial for controlling the rapid transients. However, this requires an accurate but high-order model of power electronics and their underlying control loops, complexifying the stability analysis from the viewpoint of a higher control level. To overcome these challenges, this paper proposes a large-signal order reduction (LSOR) method for MGs with considerations of external control inputs and the detailed dynamics of underlying control levels based on singular perturbation theory (SPT). Specially, we innovatively proposed and strictly proved a general stability and accuracy assessment theorem that allows us to analyze the dynamic stability of a full-order nonlinear system by only leveraging our derived reduced-order model (ROM) and boundary layer model (BLM). Furthermore, this theorem furnishes a set of conditions that determine the accuracy of the developed ROM. Lastly, by embedding such a theorem into the SPT, we propose a novel LSOR approach with guaranteed accuracy and stability analysis equivalence. Case studies are conducted on MG systems to show the effectiveness of the proposed approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗

NASA NTRS1994

The stabilizing effect of compressibility in turbulent shear flow

Direct numerical simulation of turbulent homogeneous shear flow is performed in order to clarify compressibility effects on the turbulence growth in the flow. The two Mach numbers relevant to homogeneous shear flow are the turbulent Mach number M(t) and the gradient Mach number M(g). Two series of simulations are performed where the initial values of M(g) and M(t) are increased separately. The growth rate of turbulent kinetic energy is observed to decrease in both series of simulations. This 'stabilizing' effect of compressibility on the turbulent energy growth rate is observed to be substantially larger in the DNS series where the initial value of M(g) is changed. A systematic companion of the different DNS cues shows that the compressibility effect of reduced turbulent energy growth rate is primarily due to the reduced level of turbulence production and not due to explicit dilatational effects. The reduced turbulence production is not a mean density effect since the mean density remains constant in compressible homogeneous shear flow. The stabilizing effect of compressibility on the turbulence growth is observed to increase with the gradient Mach number M(g) in the homogeneous shear flow DNS. Estimates of M(g) for the mixing and the boundary layer are obtained. These estimates show that the parameter M(g) becomes much larger in the high-speed mixing layer relative to the high-speed boundary layer even though the mean flow Mach numbers are the same in the two flows. Therefore, the inhibition of turbulent energy production and consequent 'stabilizing' effect of compressibility on the turbulence (over and above that due to the mean density variation) is expected to be larger in the mixing layer relative to the boundary layer in agreement with experimental observations.

Sarkar, S.↗

NASA NTRS1985

Energy Stability of Thermocapillary Convection in Models of the Float Zone Process

The energy-stability of thermocapillary convection in models of the float-zone, crystal-growing process was studied. Stability limits, as functions of pertinent parameters, that will identify conditions which will not allow the existence of an undesirable oscillatory flow instability were determined. Such instabilities may occur in the space processing of semiconductor materials. The determination of the stability limits will involve two sets of numerical computations: (1) solution of the nonlinear governing equations together with the appropriate boundary conditions to determine the basic state (in general, velocity, pressure and temperature fields and the displacement of free surfaces and interfaces); and (2) solution of a nonlinear Euler-Lagrange systems for the energy-stability limit. Both computations, while difficult, should be within the scope of available computer capability and available concepts in numerical analysis. Finite-element methods are attractive candidates for the numerical work.

Neitzel, G. P.↗

NASA NTRS2013

Aeromechanics Analysis of a Boundary Layer Ingesting Fan

Boundary layer ingesting propulsion systems have the potential to significantly reduce fuel burn but these systems must overcome the challe nges related to aeromechanics-fan flutter stability and forced response dynamic stresses. High-fidelity computational analysis of the fan a eromechanics is integral to the ongoing effort to design a boundary layer ingesting inlet and fan for fabrication and wind-tunnel test. A t hree-dimensional, time-accurate, Reynolds-averaged Navier Stokes computational fluid dynamics code is used to study aerothermodynamic and a eromechanical behavior of the fan in response to both clean and distorted inflows. The computational aeromechanics analyses performed in th is study show an intermediate design iteration of the fan to be flutter-free at the design conditions analyzed with both clean and distorte d in-flows. Dynamic stresses from forced response have been calculated for the design rotational speed. Additional work is ongoing to expan d the analyses to off-design conditions, and for on-resonance conditions.

Bakhle, Milind A.↗

NASA NTRSDate not supplied

Isolating the Surface Type Influence on Arctic Cloud Properties

In response to anthropogenic climate change, substantial declines in sea ice have been observed. An urgent question is whether and how clouds respond to the changing Arctic surface. Due to the important radiative influence of clouds, the response of clouds to sea ice cover change constitutes a potentially important climate feedback within the Arctic. The sign of this feedback is strongly influenced by the seasonality of the cloud response to sea ice loss. Given the importance of this phenomenon, previous research focused on using observations of inter-annual variability to quantify the cloud-sea ice relationship. However, the reliance on inter-annual variability to assess the covariance between clouds and sea ice make it challenging to control for the influence of large-scale meteorology on the clouds. We have devised an approach that applies an in situ, airborne observational strategy to quantify the influence of surface type on Arctic clouds to satellite data. Our event-based method composites cloud property differences for ocean, marginal ice zone, and sea ice surface types for >5000 MIZ crossing events. We find that cloud properties in non-summer months below ~1.5 km are influenced by surface type such that greater cloud fraction and water content are found over ocean relative to sea ice regions. Results are statistically significantly at the 95% confidence level. During summer, the surface type does not influence the cloud properties. The primary cause of the surface type dependent cloud property differences is the thermodynamic profiles differences that occur between the two surface types: namely, ocean footprints are warmer, moister, have more positive surface turbulent fluxes and are less stable than their sea ice counterparts. We find that the differences in thermodynamic profiles by surface type that correspond to the cloud property differences also explains a significant portion of the variability across events. Our analysis provides evidence that the influence of surface type on cloud properties fundamentally occurs due to the inherent surface temperature differences between the surface types. Significant surface type influences on cloud properties are only found in the presence of surface temperature differences between the surface types and only non-summer months meet this criterion. A conceptual model of this process is as follows, fundamental properties of the ocean and sea ice surface types (e.g., surface albedo, heat capacity/thermal inertia, surface roughness, and thermal conductance) lead to systematic differences in the surface temperature by surface type. In response to this surface temperature difference, atmosphere boundary layer processes (including turbulent mixing and radiation) communicate the surface temperature differences to the lower atmosphere driving the greater stability over sea ice than ocean. These differences in the background thermodynamic state profiles between the surface types yield atmospheric conditions that are more conducive to cloud development over ocean than sea ice. The details of these boundary layer processes are sensitive to the surface type but may not be fundamental to the resulting thermodynamic structure and cloud differences. Thus, our results indicate that to accurately model that cloud response to sea ice loss in climate models the surface type properties that yield the systematic differences in ocean and sea ice temperature (e.g., surface albedo, surface turbulent fluxes parameters, and thermal conductance) must be accurately represented.

Arctic radiation budget↗

DOE OSTI2023

Decahedra and Icosahedra Everywhere: The Anomalous Crystallization of Au and Other Metals at the Nanoscale

Throughout history, objects with fivefold symmetry have been a popular topic of interest for artists, philosophers, and scientists. This may be because fivefold symmetry is very conspicuous in nature. In the case of crystals, fivefold rotational symmetry is mathematically forbidden, and macroscopic crystals exhibiting fivefold symmetry have never been shown to exist. Nevertheless, in the nanoworld, nanoparticles are often found with decahedral and icosahedral shapes that have an overall pseudo–fivefold symmetry. These structures are observed at many length scales, from 1 nm up to 1 µm. In this review, several reasons for the stability of fivefold nanoparticles are discussed. Finally, these include the formation of twin boundaries, surface reconstruction faceting, and other factors.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Uncovering the role of nanoscale precipitates on martensitic transformation and superelasticity

Here we characterize the role of coherent nanoscale B2 Ni 50 Al 50 precipitates on the temperature- and stress-induced martensitic phase transformation in nanocrystalline Ni 63 Al 37 shape memory alloys using multi-million-atoms molecular dynamics (MD) simulations. We studied two types of precipitates: one with single crystal precipitates (SXP) and a second where grain boundaries cut through precipitates (PXP). Simulations reveal that the presence of B2 precipitates stabilizes the cyclic flag-shaped stress-strain response, characteristic of superelasticity, and reduces remnant strain. In contrast, single-phase nanocrystalline Ni 63 Al 37 exhibits degradation of the reverse transformation during cyclic loading and, eventually, incomplete reversible transformation within a few cycles. This is consistent with previous experimental findings of ultra-low fatigue in Ni-Ti-Cu alloys with Ti 2 Cu precipitates. The simulations reveal that the presence of precipitates significantly improves the reversibility of the transformation by acting as elastic zones that partially shield the martensitic transformation and drive the reverse transformation. A detailed analysis of the MD trajectories reveals that the martensitic transformation of the matrix induces ultra-large elastic deformation in some of the B2 precipitates (approximately 12%) to the point of resulting in a martensite-like atomic structure.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Effects of Nonequilibrium Atomic Structure on Ionic Diffusivity in LLZO: A Classical and Machine Learning Molecular Dynamics Study

To improve the performance of electrochemical devices, it is essential to understand the effects of nonequilibrium motifs in solids, such as grain boundaries, amorphous phases, and highly strained regions, on atomic-scale transport and stability. Molecular dynamics simulations are used to explore the combined effect of far-from-equilibrium atomic structures and the choice of interatomic potential on ionic diffusivity predictions for Li 7 La 3 Zr 2 O 12 (LLZO), a promising solid electrolyte for all-solid-state batteries. Amorphization and high strain are considered using both classical Buckingham interatomic potentials and machine learning force fields. Here we find that both crystalline expansion and amorphization tend to slow diffusion, although the different physical encodings in the two potentials impact the properties in different ways. We trace these variations to a combination of structural and transport factors, the contributions of which are deconvoluted computationally. Graph-based analysis reveals that the variations for amorphous LLZO arise from the connectivity of diffusion pathways within the predicted structures, which generally correlates with diffusivity and is notably higher for structures generated by the machine learning force fields. Our study provides additional insight into the relationship between atomic structure and diffusivity in LLZO, while also highlighting the need for care in choosing and validating potentials to simulate far from equilibrium structures.

25 ENERGY STORAGE↗

DOE OSTI2025

Intralattice-bonded phase-engineered ultrahigh-Ni single-crystalline cathodes suppress strain evolution

Single crystallization remains a debated strategy for advancing Ni-rich cathode materials. While it mitigates particle cracking and improves tap density by eliminating particle boundaries, extended diffusion pathways introduce volumetric and lattice distortions, compromising electrochemical and structural stability. These challenges hinder the commercialization of high-Ni single-crystal cathodes, calling for a reassessment of their viability. Here, in this study, we report a structural design: intralattice-bonded phase single-crystal LiNi 0.92 Co 0.03 Mn 0.05 O 2 (IBP-SC92). This architecture maintains structural integrity while shortening diffusion pathways, resulting in almost zero electrochemical degradation during cycling. The robust structure and fast ion transport mitigate lattice strain, as confirmed by multiscale high-resolution diffraction and imaging techniques, preventing intragranular cracks and irreversible phase transitions. As a result, IBP-SC92 shows outstanding cycling stability, with nearly 100% capacity retention after 100 cycles in half cells and 94.5% retention after 1,000 cycles in full cells. This redefined single-crystal cathode represents a significant step towards the industrial adoption of high-energy-density materials.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Volumetric recombination in EMC3-EIRENE: Implementation and first application to the pre-fusion power operation phase in ITER

An upgrade to the 3D plasma boundary model EMC3-EIRENE for detached divertor plasmas is introduced and evaluated. Stabilization of the iterative approximation of a self-consistent plasma–neutral gas solution is required at low divertor plasma temperatures, and this has been achieved by linearization of the energy loss term associated with the electron–neutral gas interaction. An implicit method for volumetric recombination has been implemented, which improves agreement with SOLPS-ITER results for the pre-fusion power operation phase in ITER. Resonant magnetic perturbations for control of edge localized modes are found to result in a partially detached divertor plasma with a significantly broader reduction of particle loads in the traditional strike zone, which is attributed to volumetric recombination. The additional non-axisymmetric far scrape-off layer strike points, on the other hand, remain attached at relatively high temperature so that volumetric recombination is not effective there.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Energizing charged particles by an orbit instability in a slowly rotating magnetic field

In this work, the stability of charged particle motion in a uniform magnetic field with an added spatially uniform transverse rotating magnetic field (RMF) is studied analytically. It is found that the stability diagram of a single particle's orbit depends critically on the chosen boundary conditions. We show that for many boundary conditions and wide regions in the parameter space, RMFs oscillating far below the cyclotron frequency can cause linear instabilities in the motion which break μ invariance and energize particles. Such energization may appear at odds with the adiabatic invariance of μ; however, adiabatic invariance is an asymptotic result and does not preclude such heating by magnetic fields oscillating at slow frequencies. This mechanism may contribute to heating in the edge plasma of field-reversed configurations (FRCs) in rotamak-FRC experiments. Furthermore, these RMF-driven instabilities may significantly enhance azimuthal current drive during the formation of FRCs in such devices.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗