Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical stability”

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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Receptivity of Supersonic Boundary Layers Due To Acoustic Disturbances Over Blunt Cones

Receptivity and stability of supersonic boundary layers over a 5-degree straight cone with a blunt tip are numerically investigated at a free stream Mach number of 3.5 and at a high Reynolds number of 106/inch. Both the steady and unsteady solutions are obtained by solving the full Navier-Stokes equations using the 5th-order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and using third-order total-variation-diminishing (TVD) Runge-Kutta scheme for time integration. The linear stability results showed that bluntness has less stabilizing effects on the stability of boundary layers over cones than on flat plates and wedges. The unsteady simulations of the interaction of plane threedimensional acoustic waves with the cone showed that the modulation of wavelength and the generation of instability waves first occurred near the leading edge in the plane where the constant acoustic phase lines are perpendicular to the cone axis. Further downstream, this instability region spreads in the azimuthal direction from this plane.

Balakumar, P.↗
DOE OSTI2024

Tri‐Metallic Catalyst for Oxygen Evolution Reaction Enables Continuous Operation of Anion Exchange Membrane Electrolyzer at 1A cm −2 for Hundreds of Hours

Abstract Although numerous efforts are made to synthesize active electrocatalysts for green hydrogen production; catalyst stability, and facile synthesis to scale up the production are still challenging. Herein, the production of novel non‐PGM catalysts for the oxygen reduction reaction (OER) in an alkaline aqueous medium is reported, which is based on the synthesis of a trimetallic metal–organic framework (MOF) precursors. Fine‐tuning of the composition of the metal centers (Ni, Co, and Fe) shows a great effect on OER activity after the MOF undergoes dynamic chemical and structural transformations under OER conditions. In situ characterization reveals the origin of OER activity enhancement as metals’ oxidation state increases, inducing compressive mechanical strain on metal centers, enhancing the electronic conductivity through the formation of oxygen vacancies, and stronger metal–oxygen covalency. Catalysts are used in membrane electrode assembly (MEA) setup within an industrial full‐cell anion exchange membrane electrolyzer (AEMEC), showing a stable performance for 550 h without noticeable decay at 750 and 1000 mA cm −2 industrial level current densities.

Chemistry↗
DOE OSTI2021

Stability Quantification for Consensus-Based Power Flow Control between Transmission and Distribution Power Systems

With increasing integration of distributed energy resources (DERs), distribution systems (DS) with DERs are expected to provide proactive grid services. The result is that power references, or known as dispatch signals, required for DS can become faster changing than in legacy power system. Consensus-based integral controls have been proposed for the purpose of coordinating power generations of DERs in DS to match the power references. These existing controls assume the integral control signal to be sufficiently slow or constant, and ignore the potential dynamics of the integral controller when tracking more varying power references. Therefore, in this paper we present an improvement for such controls by deriving the stability condition utilizing generalized Nyquist criterion. The stability condition is quantified by a set of integral gains that guarantee stability in closed-loop system without assuming a constant integral signal. A rule-of-thumb criterion is also derived to instruct the design of the consensus topology that can provide faster convergence rate for the closed-loop system. Here, the stability and convergence improvements developed in this paper are demonstrated and verified through numerical examples.

42 ENGINEERING↗
DOE OSTI2024

NSA Site Science: Use of ARM Observations from Northern Alaska to Evaluate and Improve Prediction Capabilities

The Arctic is warming at a rate nearly double that of the rest of the planet, leading to profound changes in atmospheric, oceanic, and ice processes. The U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) user facility has played a significant role in Arctic research, operating observatories in Alaska's North Slope for over 25 years. These observatories provide a rich, and wide-reaching dataset that offers insight into atmospheric processes in northern Alaska. This report details the results of a nine-year research project (2015–2024) supported by the DOE Atmospheric Systems Research (ASR) program, that leverages data from ARM’s deployment of observing facilities at Utqiaġvik (known as the North Slope of Alaska, or NSA, site) and Oliktok Point, Alaska. The project was conducted in two phases: - Phase 1 (2015–2019): Focused on understanding key atmospheric processes at Oliktok Point, including cloud formation, high-latitude precipitation, aerosol-cloud interactions, and cloud properties. - Phase 2 (2019–2024): Extended the research to the broader North Slope region, using data from both Oliktok Point and Utqiaġvik. Topics explored included surface energy budgets, atmospheric stability, ice nucleation processes, and microphysics in Arctic clouds. The project resulted in numerous research products, including 50 peer-reviewed publications and dissertations, 169 presentations, and 10 data products. These products cover a variety of topics, including: - Cloud Macro- and Microphysical Properties: Arctic clouds play a crucial role in energy transfer, and accurate representation in models is critical. The study explored cloud transitions, ice crystal shapes, and dual-wavelength radar data to understand ice crystal habits and size distributions. - Aerosol Properties and Processes: The team examined aerosol sources in the Arctic, including industrial emissions and natural sources. Observations showed significant spatial gradients in aerosol concentrations due to human activities and wildfire smoke. The influence of aerosols on cloud formation and the surface energy budget was also assessed. - Aerosol-Cloud Interactions: Research revealed that aerosols might suppress cloud ice production, affecting cloud radiative forcing and precipitation. The impact of local industrial emissions on cloud properties was also investigated. - Contextualizing the North Slope of Alaska in the context of the broader Arctic: To understand broader trends, the project evaluated large-scale circulation patterns and the influence of weather systems on the Arctic. Studies indicated that large-scale processes play a significant role in temperature patterns and the timing of snowmelt. - Advancing ARM Observational and Modeling Capabilities: The project developed new radar data products and advanced measurement techniques, including clutter mitigation and drizzle detection. Uncrewed aerial systems (UAS) and tethered balloon systems (TBS) were deployed to gather detailed atmospheric data. Additionally, the project supported 10 early career scientists, providing training and mentorship to undergraduate interns, graduate students, postdoctoral researchers, and early career researchers. These efforts contributed to the advancement of ARM research capabilities and fostered a new generation of scientists skilled in Arctic atmospheric research. Ultimately, this ASR-supported project has provided valuable insights into Arctic atmospheric processes and their broader climate implications. Recommendations for future work include continuing support for long-term observing at Arctic locations to foster additional research, further exploration of aerosol-cloud interactions and the potential impacts of enhanced industrialization of the Arctic, and expanded use of uncrewed systems to gather data in this remote and harsh environment. Additionally, the data products developed by this work, and the data products developed through the ARM infrastructure, leave a treasure-trove of additional information that should be explored for many years to come to gain additional insight into physical processes in the Arctic atmosphere that drive the rapid changes occurring in at high latitudes and their global impact.

58 GEOSCIENCES↗
DOE OSTI2023

AWAKEN Wind Plant Simulation Comparison

A series of numerical simulations of wind farms, using different model fidelities and for different atmospheric stability conditions, were performed as a part of the American WAKE ExperimeNt. The simulations included using FLORIS wake models, a number of microscale AMR-Wind and Nalu-Wind runs, as well as idealized and complex terrain WRF runs. The largest computations used the AMR-Wind LES solver to simulate a 100 km x 100 km domain containing 541 turbines under unstable atmospheric conditions matching previous measurements, while other LES computations focused on sections of the King Plains wind farm. Results of this qualitative comparison illustrate the interactions with wind farms with large-scale ABL structures in the flow, as well as the extent of downstream wake penetration in the flow and blockage effects around wind farms.

17 WIND ENERGY↗
NASA NTRS1972

Polarization of the auroral electrojet.

Consideration of an idealized model of electrojet polarization. Precipitation from the inner edge of the electron plasma sheet creates a density maximum in the auroral-oval ionosphere, which in turn leads to Hall and Pedersen conductance maximums. It is then assumed that a uniform westward convection electric field is imposed on the lower ionosphere before polarization. Field-aligned currents must flow into the ionosphere equatorward and out of the ionosphere poleward of the Hall conductance maximum. As the convection field and ionospheric density increase during the substorm growth phase, the field-aligned current densities should eventually reach an instability threshold beyond which anomalous resistance should produce field-aligned electric fields. The partial blockage of the field-aligned currents produces an equatorward electric field and therefore a partial Cowling conductivity in the lower ionosphere. Rough numerical estimates indicate that the expected field-aligned currents can exceed the stability threshold estimated by Kindel and Kennel (1971), that 1- to 5-kV field-aligned potential drops correspond to significant electrojet enhancement, and that the required energy dissipation of field-aligned currents in the topside ionosphere, a few ergs/per sq cm per sec column, suggests significant topside modification following auroral breakup.

Coroniti, F. V.↗
NASA NTRS1981

Fast plasma heating by anomalous and inertial resistivity effects

Fast plasma heating by anomalous and inertial resistivity effects is described. A small fraction of the plasma contains strong currents that run parallel to the magnetic field and are driven by an exponentiating electric field. The anomalous character of the current dissipation is caused by the excitation of electrostatic ion cyclotron and/or ion acoustic waves. The role of resistivity due to geometrical effects is considered. Through the use of a marginal stability analysis, equations for the average electron and ion temperatures are derived and numerically solved. The evolution of the plasma is described as a path in the drift velocity diagram, in which the drift velocity is plotted as a function of the electron to ion temperature ratio.

Duijveman, A.↗
NASA NTRS1981

Fast plasma heating by anomalous and inertial resistivity effects in the solar atmosphere

A simple model is presented to describe fast plasma heating by anomalous and inertial resistivity effects. It is noted that a small fraction of the plasma contains strong currents that run parallel to the magnetic field and are driven by an exponentiating electric field. The anomalous character of the current dissipation derives from the excitation of electrostatic ion-cyclotron and/or ion-acoustic waves. The possible role of resistivity deriving from geometrical effects ('inertial resistivity') is also considered. Using a marginal stability analysis, equations for the average electron and ion temperatures are derived and numerically solved. No loss mechanisms are taken into account. The evolution of the plasma is described as a path in the drift velocity diagram, where the drift velocity is plotted as a function of the electron to ion temperature ratio.

Duijveman, A.↗
NASA NTRS1982

A variational principle for compressible fluid mechanics. Discussion of the one-dimensional theory

The second law of thermodynamics is used as a variational statement to derive a numerical procedure to satisfy the governing equations of motion. The procedure, based on numerical experimentation, appears to be stable provided the CFL condition is satisfied. This stability is manifested no matter how severe the gradients (compression or expansion) are in the flow field. For reasons of simplicity only one dimensional inviscid compressible unsteady flow is discussed here; however, the concepts and techniques are not restricted to one dimension nor are they restricted to inviscid non-reacting flow. The solution here is explicit in time. Further study is required to determine the impact of the variational principle on implicit algorithms.

Prozan, R. J.↗
NASA NTRS1981

Influence of exit impedance on finite difference solutions of transient acoustic mode propagation in ducts

The cutoff mode instability problem associated with a transient finite difference solution to the wave equation is explained. The steady-state impedance boundary condition is found to produce acoustic reflections during the initial transient, which cause finite instabilities in the cutoff modes. The stability problem is resolved by extending the duct length to prevent transient reflections. Numerical calculations are presented at forcing frequencies above, below, and nearly at the cutoff frequency, and exit impedance models are presented for use in the practical design of turbofan inlets.

Baumeister, K. J.↗
NASA NTRS1983

Growth of Goertler vortices in compressible boundary layers along curved surfaces

An investigation of the growth of the three-dimensional, counter-rotating, longitudinal type vortices is considered in two-dimensional laminar compressible boundary-layer flow. The basic approximation of the disturbance equations that includes the terms due to boundary layer growth is considered and solved numerically. These terms are shown to have large local effects near the neutral stability region. The study shows that the instability of the boundary layer with respect to the three-dimensional vortices sets in at higher Goertler number as Mach number increases. Also the maximum amplitude ratio of the vortices is reduced by about 20 percent as Mach number increases from 0 to 5.

El-Hady, N. M.↗
NASA NTRS1985

The numerical viscosity of entropy stable schemes for systems of conservation laws

Discrete approximations to hyperbolic systems of conservation laws are studied. The amount of numerical viscosity present in such schemes, is quantified and related to their entropy stability by means of comparison. To this end, conservative schemes which are also entropy conservative are constructed. These entropy conservative schemes enjoy second-order accuracy; moreover, they admit a particular interpretation within the finite-element frameworks, and hence can be formulated on various mesh configurations. It is then shown that conservative schemes are entropy stable if and only if they contain more viscosity than the mentioned above entropy conservative ones.

Tadmor, E.↗
NASA NTRS1987

The numerical viscosity of entropy stable schemes for systems of conservation laws. I

Discrete approximations to hyperbolic systems of conservation laws are studied. The amount of numerical viscosity present in such schemes is quantified and related to their entropy stability by means of comparison.To this end, conservative schemes which are also entropy-conservative are constructed. These entropy-conservative schemes enjoy second-order accuracy; moreover, they can be interpreted as piecewise-linear finite-element methods, and hence can be formulated on various mesh configurations. It is then shown that conservative schemes are entropy stable, if and (for three-point schemes) only they contain more viscosity than that present in the above-mentioned entropy-conservative ones.

Tadmor, Eitan↗
NASA NTRS1988

Quiet Short-Haul Research Aircraft - A summary of flight research since 1981

The Quiet Short-Haul Research Aircraft (QSRA), designed for flight investigation into powered-lift terminal area operations, first flew in 1978 and has flown 600 hours since. This report summarizes QSRA research since 1981. Numerous aerodynamic flight experiments have been conducted including research with an advanced concept stability and control augmentation and pilot display system for category III instrument landings. An electromechanical actuator system was flown to assess performance and reliability. A second ground-based test was conducted to further evaluate circulation-control-wing/upper-surface-blowing performance. QSRA technology has been transferred through reports, guest pilot evaluations and airshow participation. QSRA future research thoughts and an extensive report bibliography are also presented.

Riddle, Dennis W.↗
NASA NTRS1988

Stability of swirling gas flows

The stability of inviscid swirling gas flows to small nonaxisymmetric perturbations is considered. For small Brunt-Vaisala frequencies, the problem reduces to the classical Sturm-Liouville form and the oscillation theorem can be applied. The resulting necessary and sufficient stability condition is compared to various criteria in the literature and a limited numerical study of isothermal rigidly rotating Poiseuille flow. For given azimuthal and axial wavenumbers, it is found numerically that the higher inertial modes become unstable for successively lower Rossby numbers and that this sequence of critical values approaches the theoretical value from above.

Hultgren, Lennart S.↗
NASA NTRS1989

Two-measured variable method for wall interference assessment/correction

An iterative method for wall interference assessment and/or correction is presented for transonic flow conditions in wind tunnels equipped with two component velocity measurements on a single interface. The iterative method does not require modeling of the test article and tunnel wall boundary conditions. Analytical proof for the convergence and stability of the iterative method is shown in the subsonic flow regime. The numerical solutions are given for both 2-D and axisymmetrical cases at transonic speeds with the application of global Mach number correction.

Lo, C. F.↗
NASA NTRS1988

Rotordynamic analysis of a bearing tester

The properties of the solutions of a system of four coupled nonlinear differential equations that model the behavior of the rotating shaft of a bearing tester are studied. In particular, it is shown how the bounds for the rotations of these equations can be obtained from bounds for the solutions of the linearized equations. By studying the behavior of the Fourier transforms of the solution, the approach to the stability boundary can also be predicted. These conclusions are verified by means of numerical solutions of the equations, and of power spectrum density (PSD) plots.

Zalik, Richard A.↗
NASA NTRS1991

Computer modeling of pulsed CO2 lasers for lidar applications

The experimental results will enable a comparison of the numerical code output with experimental data. This will ensure verification of the validity of the code. The measurements were made on a modified commercial CO2 laser. Results are listed as following. (1) The pulse shape and energy dependence on gas pressure were measured. (2) The intrapulse frequency chirp due to plasma and laser induced medium perturbation effects were determined. A simple numerical model showed quantitative agreement with these measurements. The pulse to pulse frequency stability was also determined. (3) The dependence was measured of the laser transverse mode stability on cavity length. A simple analysis of this dependence in terms of changes to the equivalent fresnel number and the cavity magnification was performed. (4) An analysis was made of the discharge pulse shape which enabled the low efficiency of the laser to be explained in terms of poor coupling of the electrical energy into the vibrational levels. And (5) the existing laser resonator code was changed to allow it to run on the Cray XMP under the new operating system.

Spiers, Gary D.↗