Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical stabilization”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Optimization of cascade blade mistuning. II - Global optimum and numerical optimization

The values of the mistuning which yield the most stable eigenvectors are analytically determined, using the simplified equations of motion which were developed in Part I of this work. It is shown that random mistunings, if large enough, may lead to the maximal stability, whereas the alternate mistunings cannot. The problem of obtaining maximum stability for minimal mistuning is formulated, based on numerical optimization techniques. Several local minima are obtained using different starting mistuning vectors. The starting vectors which lead to the global minimum are identified. It is analytically shown that all minima appear in multiplicities which are equal to the number of compressor blades. The effect of mistuning on the flutter speed is studied using both an optimum mistuning vector and an alternate mistuning vector. Effects of mistunings in elastic axis locations are shown to have a negligible effect on the eigenvalues. Finally, it is shown that any general two-dimensional bending-torsion system can be reduced to an equivalent uncoupled torsional system.

Nissim, E.↗
NASA NTRS1982

Application of matrix singular value properties for evaluating gain and phase margins of multiloop systems

A stability margin evaluation method in terms of simultaneous gain and phase changes in all loops of a multiloop system is presented. A universal gain-phase margin evaluation diagram is constructed by generalizing an existing method using matrix singular value properties. Using this diagram and computing the minimum singular value of the system return difference matrix over the operating frequency range, regions of guaranteed stability margins can be obtained. Singular values are computed for a wing flutter suppression and a drone lateral attitude control problem. The numerical results indicate that this method predicts quite conservative stability margins. In the second example if the eigenvalue magnitude is used instead of the singular value, as a measure of nearness to singularity, more realistic stability margins are obtained. However, this relaxed measure generally cannot guarantee global stability.

Mukhopadhyay, V.↗
DOE OSTI2022

Numerical investigation of the influence of shear and thermal stratification on the wind turbine tip‐vortex stability

Summary The interaction between wind turbine wakes and atmospheric turbulence is characterised by complex dynamics. In this study, two major components of the atmospheric boundary layer dynamics have been isolated, namely, the mean velocity profile shear and the thermal stratification, to examine their impact on the near‐wake development by undertaking a series of highly resolved large‐eddy simulations. Subsequently, instantaneous flow fields are extracted from the simulations and used to conduct Fourier analysis and proper orthogonal decomposition (POD) and compute the mean kinetic energy fluxes by different POD modes to better understand the tip‐vortex instability mechanisms. Our findings indicate that shear can significantly affect the breakup of the wind turbine tip‐vortices and the shape and stable length of the wake, whereas thermal stratification seems to only have limited contribution to the spatial development of the near‐wake field. Finally, our analysis shows that the applied perturbation frequency determines the tip‐vortex breakup location as it controls the onset of the mutual inductance instability.

17 WIND ENERGY↗
DOE OSTI2014

High-Order Numerical Methods for the Simulation of Linear and Nonlinear Waves: High-Frequency Radiation and Dynamic Stability (Final Report)

Over the past four years the Principal Investigator (PI) David Nicholls has worked on several projects in connection with award DE-SC0001549. Of the greatest import has been the continued supervision of five Ph.D. students (Robyn Canning, Travis McBride, Andrew Sward, Zheng Fang, and Venu Tammali). Canning and McBride defended their theses and graduated in May 2012, while Sward defended his thesis and graduated in May 2013. Both Fang and Tammali plan to defend their theses within the year and graduate in May 2015. Fang is now a very experienced graduate researcher with one paper accepted for publication and another in preparation. Tammali is nearly to the point of writing a paper and will work this summer as an intern at Argonne National Laboratory in the Mathematics and Computer Science Division under the supervision of Paul Fischer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Optimal damping profile ratios for stabilization of perfectly matched layers in general anisotropic media

Conventional perfectly matched layers (PML) can be unstable for certain kinds of anisotropic media. Multi-axial PML removes such instability using nonzero damping coe cients in the directions tangential with the PML interface. While using non-zero damping pro le ratios can stabilize PML, it is important to obtain the smallest possible damping pro le ratios to minimize arti cial re ections caused by these non-zero ratios, particularly for 3D general anisotropic media. Using the eigenvectors of the PML system matrix, we develop a straightforward and e cient numerical algorithm to determine the optimal damping pro le ratios to stabilize PML in 2D and 3D general anisotropic media. Numerical examples show that our algorithm provides optimal damping pro le ratios to ensure the stability of PML and complex-frequency-shifted PML for elastic-wave modeling in 2D and 3D general anisotropic media.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

Convergence acceleration of implicit schemes in the presence of high aspect ratio grid cells

The performance of Navier-Stokes codes are influenced by several phenomena. For example, the robustness of the code may be compromised by the lack of grid resolution, by a need for more precise initial conditions or because all or part of the flowfield lies outside the flow regime in which the algorithm converges efficiently. A primary example of the latter effect is the presence of extended low Mach number and/or low Reynolds number regions which cause convergence deterioration of time marching algorithms. Recent research into this problem by several workers including the present authors has largely negated this difficulty through the introduction of time-derivative preconditioning. In the present paper, we employ the preconditioned algorithm to address convergence difficulties arising from sensitivity to grid stretching and high aspect ratio grid cells. Strong grid stretching is particularly characteristic of turbulent flow calculations where the grid must be refined very tightly in the dimension normal to the wall, without a similar refinement in the tangential direction. High aspect ratio grid cells also arise in problems that involve high aspect ratio domains such as combustor coolant channels. In both situations, the high aspect ratio cells can lead to extreme deterioration in convergence. It is the purpose of the present paper to address the reasons for this adverse response to grid stretching and to suggest methods for enhancing convergence under such circumstances. Numerical algorithms typically possess a maximum allowable or optimum value for the time step size, expressed in non-dimensional terms as a CFL number or vonNeumann number (VNN). In the presence of high aspect ratio cells, the smallest dimension of the grid cell controls the time step size causing it to be extremely small, which in turn results in the deterioration of convergence behavior. For explicit schemes, this time step limitation cannot be exceeded without violating stability restrictions of the scheme. On the other hand, for implicit schemes, which are typically unconditionally stable, there appears to be room for improvement through careful tailoring of the time step definition based on results of linear stability analyses. In the present paper, we focus on the central-differenced alternating direction implicit (ADI) scheme. The understanding garnered from this analyses can then be applied to other implicit schemes. In order to systematically study the effects of aspect ratio and the methods of mitigating the associated problems, we use a two pronged approach. We use stability analyses as a tool for predicting numerical convergence behavior and numerical experiments on simple model problems to verify predicted trends. Based on these analyses, we determine that efficient convergence may be obtained at all aspect ratios by getting a combination of things right. Primary among these are the proper definition of the time step size, proper selection of viscous preconditioner and the precise treatment of boundary conditions. These algorithmic improvements are then applied to a variety of test cases to demonstrate uniform convergence at all aspect ratios.

Buelow, B. E. O.↗
DOE OSTI2019

Extended virtual element method for the Laplace problem with singularities and discontinuities

In this paper, we propose the extended virtual element method (X-VEM) to treat singularities and crack discontinuities that arise in the Laplace problem. The virtual element method (VEM) is a stabilized Galerkin formulation on arbitrary polytopal meshes, wherein the basis functions are implicit (virtual)—they are not known explicitly nor do they need to be computed within the problem domain. In this work, suitable projection operators are used to decompose the bilinear form on each element into two parts: a consistent term that reproduces the first-order polynomial space and a correction term that ensures stability. A similar approach is pursued in the X-VEM with a few notable extensions. To capture singularities and discontinuities in the discrete space, we augment the standard virtual element space with an additional contribution that consists of the product of virtual nodal basis (partition-of-unity) functions with enrichment functions. For discontinuities, basis functions are discontinuous across the crack and for singularities a weakly singular enrichment function that satisfies the Laplace equation is chosen. For the Laplace problem with a singularity, we devise an extended projector that maps functions that lie in the extended virtual element space onto linear polynomials and the enrichment function, whereas for the discontinuous problem, the consistent element stiffness matrix has a block-structure that is readily computed. An adaptive homogeneous numerical integration method is used to accurately and efficiently (no element-partitioning is required) compute integrals with integrands that are weakly singular. Once the element projection matrix is computed, the same steps as in the standard VEM are followed to compute the element stabilization matrix. Numerical experiments are performed on quadrilateral and polygonal (convex and nonconvex elements) meshes for the problem of an $\textit{L}$-shaped domain with a corner singularity and the problem of a cracked membrane under mode III loading, and results are presented that affirm the sound accuracy and demonstrate the optimal rates of convergence in the $\textit{L}^2$ norm and energy of the proposed method.

42 ENGINEERING↗
DOE OSTI2022

Statics and dynamics of skyrmions interacting with disorder and nanostructures

Magnetic skyrmions are topologically stable nanoscale particlelike objects that were discovered in 2009. Since that time, intense research interest in the field has led to the identification of numerous compounds that support skyrmions over a range of conditions spanning from cryogenic to room temperatures. Skyrmions can be set into motion under various types of driving, and the combination of their size, stability, and dynamics makes them ideal candidates for numerous applications. At the same time, skyrmions represent a new class of system in which the energy scales of the skyrmion-skyrmion interactions, sample disorder, temperature, and drive can compete. A growing body of work indicates that the static and dynamic states of skyrmions can be influenced strongly by pinning or disorder in the sample; thus, an understanding of such effects is essential for the eventual use of skyrmions in applications. The current state of knowledge regarding individual skyrmions and skyrmion assemblies interacting with quenched disorder or pinning is reviewed. Additionally, the microscopic mechanisms for skyrmion pinning, including the repulsive and attractive interactions that can arise from impurities, grain boundaries, or nanostructures, are outlined. This is followed by descriptions of depinning phenomena, sliding states over disorder, the effect of pinning on the skyrmion Hall angle, the competition between thermal and pinning effects, the control of skyrmion motion using ordered potential landscapes such as one- or two-dimensional periodic asymmetric substrates, the creation of skyrmion diodes, and skyrmion ratchet effects. Highlighted are the distinctions arising from internal modes and the strong gyrotropic or Magnus forces that cause the dynamical states of skyrmions to differ from those of other systems with pinning, such as vortices in type-II superconductors, charge density waves, or colloidal particles. Throughout this review future directions and open questions related to the pinning and dynamics in skyrmion systems are also discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2016

The stability of H/V spectral ratios from noise measurements in Armutlu Peninsula (Turkey)

The horizontal to vertical spectral ratio (H/V) method has been successfully using in order to estimate the fundamental resonance frequency of the sedimentary cover, its thickness and amplification factor since at least 3 decades. There are numerous studies have been carried out on the stability of the H/V spectral ratios. Almost all studies showed that fundamental frequency is stable even measurements are repeated at different times. From this point of view, the results will show us an approach whether the stations are suitable for accurate estimate of earthquake studies and engineering purposes or not. Also we want to see if any effects of the amplification factor changing on the seismograms for Armutlu Seismic Network (ARNET) even though seismic stations are established far away from cultural noise and located on hard rock sites. It has been selected one hour recorded data of all stations during the most stationary times. The amplification and resonant frequency variations of H/V ratio were calculated to investigate temporal stability in time. There is a total harmony in fundamental frequencies values and H/V spectral ratio values in time-lagged periods. Some stations shows secondary minor peaks in high frequency band due to a shallow formation effect or cultural noises around. In the east side of the area ILYS station shows amplitude peak in lower fundamental frequency band from expected. This could compose a high amplification in lower frequencies and so that yield less reliable results in local earthquakes studies. By the experimental results from ambient noise analysis, it could be worked up for relocation of one station.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1988

Stability and transition in supersonic boundary layers

The three-dimensional time-dependent, compressible Navier-Stokes equations are numerically solved by a Fourier-Chebyshev collocation algorithm to study the stability of a Mach 4.5 flow over a flat plate. Several nonlinear direct simulations suggest the existence of a secondary instability which might provide a possible route to transition. Pertinent differences in the energy content of the various Fourier modes between this instability and the more common incompressible K-type instabilities are pointed out.

Erlebacher, Gordon↗
NASA NTRS1986

Boundary layer transition

The boundary layer stability, its active control by sound and surface heating and the effect of curvature are studied numerically and experimentally for subsonic flow. In addition, the experimental and flight test data are correlated using the stability theory for supersonic Mach numbers. Active transition fixing and feedback control of boundary layer by sound interactions are experimentally investigated at low speed over an airfoil. Numerical simulation of active control by surface heating and cooling in air shows that by appropriate phase adjustment a reduction in the level of perturbation can be obtained. This simulation is based on the solution of two-dimensional compressible Navier-Stokes equations for a flat plate. Goertler vortices are studied experimentally on an airfoil in the Low Turbulence Pressure Tunnel (LTPT). The flow pattern was visualized using the sublimating chemical technique and data were obtained using a three component laser velocimeter. The effect of curvature on swept leading-edge stability on a cylinder was numerically studied. The results suggest that transition is dominated by traveling disturbance waves and that the waves with the greatest total amplification has an amplitude ratio of e sup 11. Experimental data from the quiet supersonic tunnel and flight tests are analyzed using linear compressible stability theory.

Maestrello, L.↗
DOE OSTI2021

Influence of Protection Layers on Thermal Stability of Nitride Thin Films

Nitrides are commonly studied for their properties relevant to numerous coating applications, and have been recently used to synthesize newly predicted materials in thin film form. However, thermodynamic stability of such nitride materials is difficult to experimentally evaluate, due to limited thermal budget of thin films. Here we show that thermal stability of nitride films can be extended using protection layers. Specifically, we find that zirconium nitride (ZrN) thin films are stable up to at least 1200°C annealing temperature in N2 atmosphere, if aluminum nitride (AlN) diffusion barriers and capping layers are used. X-ray diffraction (XRD) measurements show the expected thermodynamically stable rocksalt structure of ZrN, and scanning electron microscopy (SEM) confirms compact microstructure of these ZrN films. Without such protection layers, transient bixbyite Zr2ON2 phase and high-temperature ZrSi2 phase are observed after annealing by XRD, SEM, and auger electron spectroscopy (AES), due to side reaction with native surface oxide (ZrOx) and substrate (Si) respectively. These results demonstrate that protection layers are beneficial to increase the thermal budget of nitride thin films, for evaluating thermodynamic stability of new nitride materials, for capturing transient metastable phases during crystallization, and for using functional nitride coatings in extreme environments.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Efficient solution of the Wigner–Liouville equation using a spectral decomposition of the force field

The Wigner–Liouville equation is reformulated using a spectral decomposition of the classical force field instead of the potential energy. The latter is shown to simplify the Wigner–Liouville kernel both conceptually and numerically as the spectral force Wigner–Liouville equation avoids the numerical evaluation of the highly oscillatory Wigner kernel which is nonlocal in both position and momentum. The quantum mechanical evolution is instead governed by a term local in space and non-local in momentum, where the non-locality in momentum has only a limited range. An interpretation of the time evolution in terms of two processes is presented; a classical evolution under the influence of the averaged driving field, and a probability-preserving quantum-mechanical generation and annihilation term. Using the inherent stability and reduced complexity, a direct deterministic numerical implementation using Chebyshev and Fourier pseudo-spectral methods is detailed. For the purpose of illustration, we present results for the time-evolution of a one-dimensional resonant tunneling diode driven out of equilibrium.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1982

Instability of a solidifying binary mixture

An analysis is performed on the stability of a solidifying binary mixture due to surface tension variation of the free liquid surface. The basic state solution is obtained numerically as a nonstationary function of time. Due to the time dependence of the basic state, the stability analysis is of the global type which utilizes a variational technique. Also due to the fact that the basic state is a complex function of both space and time, the stability analysis is performed through numerical means.

Antar, B. N.↗
NASA NTRS1986

On compactness of admissible parameter sets: Convergence and stability in inverse problems for distributed parameter systems

A series of numerical examples is reported and several algorithms compared for estimation of coefficients in differential equation models. Unconstrained, constrained and Tikhonov regularization methods are tested for their behavior with regard to both convergence (of approximation methods for the states and parameters) and stability (continuity of the estimates with respect to perturbations in the data or observed states).

Banks, H. T.↗
DOE OSTI2016

From stochastic processes to numerical methods: A new scheme for solving reaction subdiffusion fractional partial differential equations

We have introduced a new explicit numerical method, based on a discrete stochastic process, for solving a class of fractional partial differential equations that model reaction subdiffusion. The scheme is derived from the master equations for the evolution of the probability density of a sum of discrete time random walks. We show that the diffusion limit of the master equations recovers the fractional partial differential equation of interest. This limiting procedure guarantees the consistency of the numerical scheme. The positivity of the solution and stability results are simply obtained, provided that the underlying process is well posed. We also show that the method can be applied to standard reaction–diffusion equations. This work highlights the broader applicability of using discrete stochastic processes to provide numerical schemes for partial differential equations, including fractional partial differential equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

ML-Based Power System Stability Assessment Considering Network Topology Changes: WECC 20,000+ Bus System Case Study

Modern power grids are fast-changing and thus require real-time monitoring and online stability assessment. With the rapid development of machine learning (ML) techniques, using data-driven models to provide fast and accurate estimations of power system stability marginal information, such as frequency nadir for frequency stability and critical clearing time (CCT) for transient stability, have become possible. However, despite the numerous research on ML-based methods for frequency nadir and CCT prediction, there is limited work on the impact of different network topology changes. Furthermore, most previous studies only focused on small or synthetic systems, and there is a lack of research on actual large power system models. In this paper, the above issues are addressed by studying the actual U.S. Western Electricity Coordinating Council (WECC) system model with more than 20,000 buses. Massive simulations are conducted in PowerWorld Simulator to study the impact of various topology change scenarios on both frequency stability and transient stability. System operating information is extracted from the success dispatch cases of various network topologies to generate a comprehensive dataset for ML-based models. Two ML methods, random forest (RF) and multilayer perceptron (MLP) neural network, are trained and tested for both frequency nadir prediction and CCT prediction. Test results have proven the models are capable of online stability assessment for large power networks such as the WECC system with sufficient accuracy.

critical clearing time↗
DOE OSTI2015

Stabilization of ring dark solitons in Bose-Einstein condensates

Earlier work has shown that ring dark solitons in two-dimensional Bose-Einstein condensates are generically unstable. In this work, we propose a way of stabilizing the ring dark soliton via a radial Gaussian external potential. We investigate the existence and stability of the ring dark soliton upon variations of the chemical potential and also of the strength of the radial potential. Numerical results show that the ring dark soliton can be stabilized in a suitable interval of external potential strengths and chemical potentials. Furthermore, we also explore different proposed particle pictures considering the ring as a moving particle and find, where appropriate, results in very good qualitative and also reasonable quantitative agreement with the numerical findings.

74 ATOMIC AND MOLECULAR PHYSICS↗