Search NASA⌕ Search

SEARCH · Search NASA

Results for “stabilization method”

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,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Semi-discrete approximations to nonlinear systems of conservation laws; consistency and L(infinity)-stability imply convergence

A convergence theory for semi-discrete approximations to nonlinear systems of conservation laws is developed. It is shown, by a series of scalar counter-examples, that consistency with the conservation law alone does not guarantee convergence. Instead, a notion of consistency which takes into account both the conservation law and its augmenting entropy condition is introduced. In this context it is concluded that consistency and L(infinity)-stability guarantee for a relevant class of admissible entropy functions, that their entropy production rate belongs to a compact subset of H(loc)sup -1 (x,t). One can now use compensated compactness arguments in order to turn this conclusion into a convergence proof. The current state of the art for these arguments includes the scalar and a wide class of 2 x 2 systems of conservation laws. The general framework of the vanishing viscosity method is studied as an effective way to meet the consistency and L(infinity)-stability requirements. How this method is utilized to enforce consistency and stability for scalar conservation laws is shown. In this context we prove, under the appropriate assumptions, the convergence of finite difference approximations (e.g., the high resolution TVD and UNO methods), finite element approximations (e.g., the Streamline-Diffusion methods) and spectral and pseudospectral approximations (e.g., the Spectral Viscosity methods).

Tadmor, Eitan↗
NASA NTRS2021

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗
NASA NTRS1971

Film stability in a vertical rotating tube with a core-gas flow.

The linear hydrodynamic stability of a thin-liquid layer flowing along the inside wall of a vertical tube rotating about its axis in the presence of a core-gas flow is examined. The stability problem is formulated under the conditions that the liquid film is thin, the density and viscosity ratios of gas to liquid are small and the relative (axial) pressure gradient in the gas is of the same order as gravity. The resulting eigenvalue problem is first solved by a perturbation method appropriate to axisymmetric long-wave disturbances. The damped nature (to within the thin-film and other approximations made) of the nonaxisymmetric and short-wave disturbances is noted. In view of the limitations on a truncated perturbation solution when the disturbance wavenumber is not small, an initial value method using digital computer is presented. Stability characteristics of neutral, growing, and damped modes are presented showing the influences of rotation, surface tension, and the core-gas flow. Energy balance in a neutral mode is also illustrated.

Sarma, G. S. R.↗
NASA NTRS1991

Nonlinear rotordynamics analysis

Effective analysis tools were developed for predicting the nonlinear rotordynamic behavior of the Space Shuttle Main Engine (SSME) turbopumps under steady and transient operating conditions. Using these methods, preliminary parametric studies were conducted on both generic and actual HPOTP (high pressure oxygen turbopump) models. In particular, a novel modified harmonic balance/alternating Fourier transform (HB/AFT) method was developed and used to conduct a preliminary study of the effects of fluid, bearing and seal forces on the unbalanced response of a multi-disk rotor in the presence of bearing clearances. The method makes it possible to determine periodic, sub-, super-synchronous and chaotic responses of a rotor system. The method also yields information about the stability of the obtained response, thus allowing bifurcation analyses. This provides a more effective capability for predicting the response under transient conditions by searching in proximity of resonance peaks. Preliminary results were also obtained for the nonlinear transient response of an actual HPOTP model using an efficient, newly developed numerical method based on convolution integration. Currently, the HB/AFT is being extended for determining the aperiodic response of nonlinear systems. Initial results show the method to be promising.

Noah, Sherif T.↗
DOE OSTI2022

Stability Analysis of Coupled Advection-Diffusion Models with Bulk Interface Condition

Numerical stability is of critical importance in general circulation models because it affects the design of algorithms, time to solution, and computational costs associated with the simulations, which are very expensive in practice. In this paper we extend the stability analysis for ocean-atmosphere coupling proposed in [Zhang et al., J. Sci. Comput. 84, 44(2020)] to a more realistic model that includes horizontal advection. We analyze various time-stepping strategies. We find that advection has a stabilizing effect in scenarios common to climate models when bulk interface condition and explicit flux coupling are used. We also show that our method can be used to study the stability impact of advection for other interface conditions such as Dirichlet-Neumann conditions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Development of the Active Correlation Method: Theoretical-Methodological Aspect

The Dubna Gas-Filled Recoil Separator (DGFRS) of the Flerov Laboratory of Nuclear Reactions (FLNR), JINR, is the most advanced facility for the synthesis and study of new superheavy nuclei. In the recent years, new elements with Z = 114 to 118 (Fl, Mc, Lv, Ts, Og) have been successfully synthesized. The DGFRS detection system and a unique method of active correlations for background suppression have played a significant role in these discoveries. Theoretical-methodological aspects of further development of the active correlation method are considered, especially in view of the upcoming commissioning of the new FLNR high-intensity DC-280 cyclotron for acceleration of heavy ions and the new gas-filled recoil separator. A numerical model of the edge effects between the neighboring strips on the p – n junction side of the DSSSD detector is presented. The corresponding empirical examples are given. A more flexible real-time algorithm is considered as a possible substitute for the current version with the rigidly set parameters. Since stability of the calibration parameters is strongly required for applying the method, the radiation stability factor is also briefly considered.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Impact of Humidity and Temperature on the Stability of the Optical Properties and Structure of MAPbI3, MA0.7FA0.3PbI3 and (FAPbI3)0.95(MAPbBr3)0.05 Perovskite Thin Films

In situ real-time spectroscopic ellipsometry (RTSE) measurements have been conducted on MAPbI3, MA0.7FA0.3PbI3, and (FAPbI3)0.95(MAPbBr3)0.05 perovskite thin films when exposed to different levels of relative humidity at given temperatures over time. Analysis of RTSE measurements track changes in the complex dielectric function spectra and structure, which indicate variations in stability influenced by the underlying material, preparation method, and perovskite composition. MAPbI3 and MA0.7FA0.3PbI3 films deposited on commercial fluorine-doped tin oxide coated glass are more stable than corresponding films deposited on soda lime glass directly. (FAPbI3)0.95(MAPbBr3)0.05 films on soda lime glass showed improved stability over the other compositions regardless of the substrate, and this is attributed to the preparation method as well as the final composition.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Intermetallic PtCo Catalysts with Enhanced Performance and Stability [Slides]

Pt-based catalysts need systematic modifications to be adapted to different applications (LDV or HDV). L1 0 -CoPt particles have been demonstrated as highly active and durable cathode materials for fuel cells. Colloidal synthesis can generate catalysts with highly uniform particles size but the removal of surfactants is challenging. Acid leaching warrants a Pt shell with certain thickness, which mitigates the Co leaching. Controlling the size of L1 0 -CoPt particles is a robust method to alter their catalytic behaviours and tailor the function of catalysts. Modifying the carbon features (e.g. porosity and graphitization) is an effective method to regulate the performance and stability of catalysts to adapt the catalysts to different applications.

30 DIRECT ENERGY CONVERSION↗
NASA NTRS1980

Experimental determination of pure rotary stability derivatives using curved and rolling flow wind tunnel

The technique of using a curved and rolling flow wind tunnel to extract pure rotary stability derivatives is presented. Descriptions of the curved flow and the rolling flow test sections of the Virginia Tech Stability Wind Tunnel are given including methods for obtaining the proper velocity profiles and correcting the data acquired. Results of testing current fighter configurations in this facility are presented with particular attention given to comparing pure rotary derivatives with combined rotary and unsteady derivatives obtained by standard oscillation tests. Also the effect of curved and rolling flow on lateral static stability derivatives is examined.

Lutze, F. H.↗
NASA NTRS2019

High-Order Finite-Difference Nonlinear Filter Methods for Subsonic Turbulence Simulation with Stochastic Forcing

Numerical stability of high-order filter schemes developed by Yee & Sjogreen is tested on three-dimensional turbulence simulations with stochastic forcing and their performance is compared with that of TVD and WENO schemes. The best­ performing filter method employs an eighth-order central base scheme with the Kennedy & Gruber skew-symmetric splitting of the inviscid flux derivative, a wavelet-based local flow sensor, a nonlinear filter utilizing the dissipative portion of seventh-order \VENO scheme, and an explicit third - or fourth-order Runge­-Kutta time integration. We show that the filter scheme is more computational]y efficient and provides a wider spectral bandwidth compared to the seventh-order WENO scheme. The method also demonstrates robust long-time integration for moderately compressible turbulence. In contrast, the fifth- and seventh­ order WENO schemes show non-trivial evolution of the velocity and density power spectra. over a. few dozen dynamical times, where both TVD and filter schemes recover a so lid statistically stationary turbulent state

Kritsuk, Alexei G.↗
DOE OSTI2023

Additive Polynomial Time Integrators, Part I: Framework and Fully Implicit-Explicit Collocation Methods

In this paper we generalize the polynomial time integration framework to additively partitioned initial value problems. The framework we present is general and enables the construction of many new families of additive integrators with arbitrary order-of-accuracy and varying degree of implicitness. In this first work, we focus on a new class of implicit-explicit polynomial block methods that are based on fully implicit Runge–Kutta methods with Radau nodes and that possess high stage order. Here, we show that the new fully implicit-explicit (FIMEX) integrators have improved stability compared to existing IMEX Runge–Kutta methods, while also being more computationally efficient due to recent developments in preconditioning techniques for solving the associated systems of nonlinear equations. For PDEs on periodic domains where the implicit component is trivial to invert, we will show how parallelization of the right-hand side evaluations can be exploited to obtain significant speedup compared to existing serial IMEX Runge–Kutta methods. For parallel (in space) finite element discretizations, the new methods can achieve orders of magnitude better accuracy than existing IMEX Runge–Kutta methods and/or achieve a given accuracy several times times faster in terms of computational runtime.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1972

Some methods for determining the profile functions of inhomogeneous media

The problem of determining the electric permittivity of an inhomogeneous dielectric medium is considered. Two analytical methods and one numerical method for solving this problem are presented. All these methods utilize the reflection coefficient measured as a function of the angular spectrum variable. This allows measurements at a fixed frequency, and these methods, therefore, can be used when the medium being studied is dispensive. A few examples are provided and some comments are made regarding the stability of solutions obtained using the various methods.

Mittra, R.↗
NASA NTRS1999

An Evaluation of Oceanographic Optical Instruments and Deployment Methodologies

The primary objective of the Sea-viewing, Wide Field-of-view Sensor (SeaWiFS) Project is to produce water- leaving radiances with an uncertainty of 5% in clear-water regions and chlorophyll a concentrations within +/- 35% over the range of 0.05-50 mg/cu m. Any global mission, like SeaWiFS, requires validation data be submitted from a wide variety of investigators which places a significant challenge on quantifying the total uncertainty associated with the in situ measurements, because each investigator follows slightly different practices when it comes to implementing all of the steps associated with collecting field data, even those with a prescribed set of protocols. This study uses data from multiple cruises to quantify the uncertainties associated with implementing data collection procedures while utilizing differing in-water optical instruments and deployment methods. A comprehensive approach is undertaken and includes: (1) the use of a portable light source and in-water intercomparisons to monitor the stability of the field radiometers, (2) alternative methods for acquiring reference measurements, and (3) different techniques for making in-water profiles. The only system to meet the 5% radiometric objective of the SeaWiFS Project was a free-fall profiler using (relatively inexpensive) modular components, although a more sophisticated (and comparatively expensive) profiler using integral components was very close and only 1% higher. A relatively inexpensive system deployed with a winch and crane was also close, but the ship shadow contamination problem increased the total uncertainty to approximately 6.5%.

Hooker, Stanford B.↗
NASA NTRS2022

Transition Analysis for the CRM-NLF Wind Tunnel Configuration using Transport Equation Models and Linear Stability Correlations

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework rely upon transition correlations that were derived from a limited number of low-speed experiments. Furthermore, these models often account for only a subset of the relevant transition mechanisms and/or cannot accurately predict the sensitivity of those mechanisms to the changes in significant flow parameters. A preceding investigation had targeted the assessment of the transport-equation-based transition models in NASA's OVERFLOW 2.3b solver, namely, the amplification factor transport (AFT-2017b) equation model coupled with the Spalart-Allmaras RANS model and the Langtry-Menter transition models (LM2009 without crossflow effects and LM2015 including the modeling of crossflow transition) implemented with Menter’s shear-stress transport equation (SST2003) RANS model. Comparisons with recent measurements at transonic freestream conditions on the Common Research Model with Natural Laminar Flow (CRM-NLF) reinforced our earlier finding that all three of the above models significantly underpredict the reported extent of the laminar flow region over the entire span of the wing, regardless of the dominant instability mechanism(s) underlying the onset of the transition process. The underprediction of the laminar flow extent was attributed to the failure of the above models in accounting for the stabilizing effect of compressibility on the amplification of Tollmien-Schlichting instabilities. Based on previous linear stability studies related to compressibility effects, the present work proposes modifications to the two classes of transition models that reduce to the original form of each model at low subsonic speeds and do not require any nonlocal flow information or additional transport equation(s). The modifications are shown to significantly improve the predicted laminar extent of the flow and compare well against the data from the CRM-NLF experiment. Additionally, a previous assessment of transition prediction based on the dual, nonparallel N -factor method in conjunction with linear parabolized stability equations (PSE) is extended to additional angles of attack to provide the first comprehensive assessment of transition models based on nonparallel disturbance amplification over the CRM-NLF. In general, the transition criterion based on the dual, nonparallel N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at R eMAC = 15 million for all angles of attack investigated herein and provides additional validation of the improved predictions from the compressibility-corrected transition models.

CFD modeling↗
DOE OSTI2019

Kinetics of Early Decomposition Stages in Diluted bcc Fe–Сu–Ni–Al Alloy: MC+MD Simulation

A combined approach including the Monte Carlo and molecular-dynamics simulation, decomposition kinetics and segregation formation in the multicomponent low-alloy Fe–1.5Cu–2Ni–1.5Al (at %) alloy is studied. It is shown that the formation of Cu nanoparticles surface-enriched with Ni and Al (coprecipitation mode) includes several stages: (i) the formation of clusters consisting of several Cu atoms, (ii) their enrichment with Ni and Al atoms, and (iii) redistribution of Ni and Al atoms with the formation of a surface layer providing stabilization of Cu nanoparticles. Observed structural features of segregations and their stability in Fe–Cu–Ni–Al alloys is discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1991

A comparison of numerical methods for the Rayleigh equation in unbounded domains

A second-order finite difference and two spectral methods, including a Chebyshev tau and a Chebyshev collocation method were implemented to determine the linear hydrodynamic stability of an unbounded shear flow. The velocity profile of the basic flow in the stability analysis mimicks that of a one-stream free mixing layer. Local and global eigenvalue solution methods are used to determine individual eigenvalues and the eigenvalue spectrum, respectively. The calculated eigenvalue spectrum includes a discrete mode, a continuous spectrum associated with the equation singularity and a continuous spectrum associated with the domain unboundedness. The efficiency and the accuracy of these discretization methods in the prediction of the eigensolutions of the discrete mode were evaluated by comparison with a conventional shooting procedure. Their capabilities in mapping out the continuous eigenvalue spectra are also discussed.

Liou, W. W.↗
NASA NTRS2003

A Discrete Constraint for Entropy Conservation and Sound Waves in Cloud-Resolving Modeling

Ideal cloud-resolving models contain little-accumulative errors. When their domain is so large that synoptic large-scale circulations are accommodated, they can be used for the simulation of the interaction between convective clouds and the large-scale circulations. This paper sets up a framework for the models, using moist entropy as a prognostic variable and employing conservative numerical schemes. The models possess no accumulative errors of thermodynamic variables when they comply with a discrete constraint on entropy conservation and sound waves. Alternatively speaking, the discrete constraint is related to the correct representation of the large-scale convergence and advection of moist entropy. Since air density is involved in entropy conservation and sound waves, the challenge is how to compute sound waves efficiently under the constraint. To address the challenge, a compensation method is introduced on the basis of a reference isothermal atmosphere whose governing equations are solved analytically. Stability analysis and numerical experiments show that the method allows the models to integrate efficiently with a large time step.

Zeng, Xi-Ping↗