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.

SourceResearch connectionFrom yearThrough year

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

At least 1,225 records · Page 68

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. Part 1: The ODE connection and its implications for algorithm development in computational fluid dynamics

Spurious stable as well as unstable steady state numerical solutions, spurious asymptotic numerical solutions of higher period, and even stable chaotic behavior can occur when finite difference methods are used to solve nonlinear differential equations (DE) numerically. The occurrence of spurious asymptotes is independent of whether the DE possesses a unique steady state or has additional periodic solutions and/or exhibits chaotic phenomena. The form of the nonlinear DEs and the type of numerical schemes are the determining factor. In addition, the occurrence of spurious steady states is not restricted to the time steps that are beyond the linearized stability limit of the scheme. In many instances, it can occur below the linearized stability limit. Therefore, it is essential for practitioners in computational sciences to be knowledgeable about the dynamical behavior of finite difference methods for nonlinear scalar DEs before the actual application of these methods to practical computations. It is also important to change the traditional way of thinking and practices when dealing with genuinely nonlinear problems. In the past, spurious asymptotes were observed in numerical computations but tended to be ignored because they all were assumed to lie beyond the linearized stability limits of the time step parameter delta t. As can be seen from the study, bifurcations to and from spurious asymptotic solutions and transitions to computational instability not only are highly scheme dependent and problem dependent, but also initial data and boundary condition dependent, and not limited to time steps that are beyond the linearized stability limit.

Yee, H. C.↗

NASA NTRS1991

Dynamical approach study of spurious steady-state numerical solutions of nonlinear differential equations. I - The dynamics of time discretization and its implications for algorithm development in computational fluid dynamics

Yee, H. C.↗

DOE OSTI2021

Quantum-Based Molecular Dynamics Simulations Using Tensor Cores

Tensor cores, along with tensor processing units, represent a new form of hardware acceleration specifically designed for deep neural network calculations in artificial intelligence applications. Tensor cores provide extraordinary computational speed and energy efficiency but with the caveat that they were designed for tensor contractions (matrix–matrix multiplications) using only low-precision floating-point operations. Despite this perceived limitation, we demonstrate how tensor cores can be applied with high efficiency to the challenging and numerically sensitive problem of quantum-based Born–Oppenheimer molecular dynamics, which requires highly accurate electronic structure optimizations and conservative force evaluations. The interatomic forces are calculated on-the-fly from an electronic structure that is obtained from a generalized deep neural network, where the computational structure naturally takes advantage of the exceptional processing power of the tensor cores and allows for high performance in excess of 100 Tflops on a single Nvidia A100 GPU. Stable molecular dynamics trajectories are generated using the framework of extended Lagrangian Born–Oppenheimer molecular dynamics, which combines computational efficiency with long-term stability, even when using approximate charge relaxations and force evaluations that are limited in accuracy by the numerically noisy conditions caused by the low-precision tensor core floating-point operations. A canonical ensemble simulation scheme is also presented, where the additional numerical noise in the calculated forces is absorbed into a Langevin-like dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1975

Non-numeric computation for high eccentricity orbits

Geocentric orbits of large eccentricity (e = 0.9 to 0.95) are significantly perturbed in cislunar space by the sun and moon. The time-history of the height of perigee, subsequent to launch, is particularly critical. The determination of 'launch windows' is mostly concerned with preventing the height of perigee from falling below its low initial value before the mission lifetime has elapsed. Between the extremes of high accuracy digital integration of the equations of motion and of using an approximate, but very fast, stability criteria method, this paper is concerned with the developement of a method of intermediate complexity using non-numeric computation. The computer is used as the theory generator to generalize Lidov's theory using six osculating elements. Symbolic integration is completely automatized and the output is a set of condensed formulae well suited for repeated applications in launch window analysis. Examples of applications are given.

Sridharan, R.↗

NASA NTRS2014

Effects of Forward- and Backward-Facing Steps on the Crossflow Receptivity and Stability in Supersonic Boundary Layers

The effects of forward- and backward-facing steps on the receptivity and stability of three-dimensional supersonic boundary layers over a swept wing with a blunt leading edge are numerically investigated for a freestream Mach number of 3 and a sweep angle of 30 degrees. The flow fields are obtained by solving the full Navier-Stokes equations. The evolution of instability waves generated by surface roughness is simulated with and without the forward- and backward-facing steps. The separation bubble lengths are about 5-10 step heights for the forward-facing step and are about 10 for the backward-facing step. The linear stability calculations show very strong instability in the separated region with a large frequency domain. The simulation results show that the presence of backward-facing steps decreases the amplitude of the stationary crossflow vortices with longer spanwise wavelengths by about fifty percent and the presence of forward-facing steps does not modify the amplitudes noticeably across the steps. The waves with the shorter wavelengths grow substantially downstream of the step in agreement with the linear stability prediction.

Balakumar, P.↗

DOE OSTI2020

Finite spatial-grid effects in energy-conserving particle-in-cell algorithms

Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. We also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length. Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. Finally, we also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Stabilization of Ultrasmall Platinum Nanoparticles by Nitrogen-Doped Carbon: Implications for Catalysis and Electrocatalysis

Heterogeneous materials comprising platinum nanoparticles on carbon supports have numerous applications including fuel cell electrodes and heterogeneous catalysts. The effective application of these materials for fuel cells and catalysis will be greatly advanced by the ability to control the oxidation and sintering of the nanoparticles by modifications of the carbon support. One attempt of such control has been doping carbon supports with nitrogen. Here, in this work, a cutting-edge, high-sensitivity, in situ XRD instrument, which allows observation of ultrasmall Pt nanoparticles, has been combined with in situ XPS to provide unprecedented clarity in the characterization of supported Pt nanoparticles in oxidizing and high-temperature environments. On a nitrogen-doped carbon support derived from poly-phenylporphyrin, Pt nanoparticles show increased stability to oxidation and thermal sintering. The enhanced Pt–support interaction arising from the N dopant versus the N-free carbon is manifested by (1) decreased initial Pt particle sizes, (2) small particle size at higher surface densities, (3) increased resistance of Pt nanoparticles to oxidation, (4) increased electron binding energy of Pt0, and (5) increased resistance of Pt nanoparticles to sintering. It is expected that the higher stability of Pt on NC will be manifested in higher activity in fuel cells and high-temperature catalytic reactions.

catalytic reactions↗

NASA NTRS1979

Numerical methods of solving a system of multi-dimensional nonlinear equations of the diffusion type

The principles of conservation and stability of difference schemes achieved using the iteration control method were examined. For the schemes obtained of the predictor-corrector type, the conversion was proved for the control sequences of approximate solutions to the precise solutions in the Sobolev metrics. Algorithms were developed for reducing the differential problem to integral relationships, whose solution methods are known, were designed. The algorithms for the problem solution are classified depending on the non-linearity of the diffusion coefficients, and practical recommendations for their effective use are given.

Agapov, A. V.↗

NASA NTRS1980

Numerical studies of new stellarator concepts

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

Bauer, F.↗

NASA NTRS1990

Numerical simulation of particle-wave interaction in boundary layers

The effects of wall injection and particle motion on the spatial stability of two-dimensional plane channel flow are investigated. For this purpose, an accurate Navier-Stokes solver to simulate the space-time evolution of disturbances in three-dimensional flows has been developed. The code is operational on the NASA Langley CRAY2 and can be ported to any other supercomputer. The code has been tested extensively in tracking the spatial evolution of two-dimensional disturbances in plane channel flow and provided excellent agreement with the linear theory including at the inflow/outflow boundaries. Preliminary calculations have been performed to investigate the effects of stationary and moving sources of vortical disturbances simulating a particle traveling in the flow field. Results suggest that even at very low amplitudes, vortical disturbances act as amplifiers on the Tollmien-Schlichting waves promoting rapid instability. It is also found that slow moving particles are more dangerous than both stationary and fast moving particles for the same disturbance levels.

Biringen, S.↗

NASA NTRS1977

Effects of finite plasma beta on the lower-hybrid-drift instability

The local dispersion relation for the lower-hybrid-drift (LHD) instability is derived and analyzed, taking into account the finite-beta effects associated with transverse electromagnetic perturbations as well as with resonant and nonresonant electron-orbit modifications due to magnetic-field gradients. The influence of finite-beta effects on the LHD instability is calculated in a fully self-consistent manner for arbitrary values of electron-ion temperature ratio, local beta, cross-field ExB velocity/ion thermal speed ratio, and other plasma parameters. Stability properties are investigated analytically for the case of cold electrons, and the local dispersion relation is solved numerically in the parameter regime of most interest for high-density plasma pinches. The results show that for all parameter regimes studied, the net effect of finite plasma beta is to reduce the maximum growth rate of the LHD instability, although the details can vary, depending on the plasma parameters. Except in the limit where the electron/ion temperature ratio tends to zero, it is found that there is a critical value of plasma beta above which the LHD instability is completely stabilized.

Davidson, R. C.↗

NASA NTRS1996

Mixing and Transition Control Studied

Considerable progress in understanding nonlinear phenomena in both unbounded and wallbounded shear flow transition has been made through the use of a combination of high- Reynolds-number asymptotic and numerical methods. The objective of this continuing work is to fully understand the nonlinear dynamics so that ultimately (1) an effective means of mixing and transition control can be developed and (2) the source terms in the aeroacoustic noise problem can be modeled more accurately. Two important aspects of the work are that (1) the disturbances evolve from strictly linear instability waves on weakly nonparallel mean flows so that the proper upstream conditions are applied in the nonlinear or wave-interaction streamwise region and (2) the asymptotic formulations lead to parabolic problems so that the question of proper out-flow boundary conditions--still a research issue for direct numerical simulations of convectively unstable shear flows--does not arise. Composite expansion techniques are used to obtain solutions that account for both mean-flow-evolution and nonlinear effects. A previously derived theory for the amplitude evolution of a two-dimensional instability wave in an incompressible mixing layer (which is in quantitative agreement with available experimental data for the first nonlinear saturation stage for a plane-jet shear layer, a circular-jet shear layer, and a mixing layer behind a splitter plate) have been extended to include a wave-interaction stage with a three-dimensional subharmonic. The ultimate wave interaction effects can either give rise to explosive growth or an equilibrium solution, both of which are intimately associated with the nonlinear self-interaction of the three dimensional component. The extended theory is being evaluated numerically. In contrast to the mixing-layer situation, earlier comparisons of theoretical predictions based on asymptotic methods and experiments in wall-bounded shear-flow transition have been somewhat lacking in one aspect or another. The current work strongly suggests that the main weakness is the underlying asymptotic representation of the linear "part" of the problem and not the explicit modeling of the nonlinear/wave-interaction effects. Consequently, the long-wave-length/high-Reynolds-number asymptotic limit for the Blasius boundary-layer stability problem was reexamined, and a new dispersion relationship for the instability waves that is uniformly valid for both the upper- and lower branch regions to the required order of approximation was obtained. A comparison with numerical results, obtained by solving the Orr-Sommerfeld stability problem, shows that the asymptotic formula provides surprisingly good results, even for values of the frequency parameter usually encountered in experimental investigations. This is particularly evident in the dynamically important upper-branch region, where much of the nonlinear interactions in transition experiments are believed to take place. The result is important in that it can be used to greatly improve the accuracy of weakly nonlinear critical-layer-based theories, and a consistent nonlinear theory is currently under evaluation.

Source record↗

DOE OSTI2025

An adaptive and stability-promoting layerwise training approach for sparse deep neural network architecture

This work presents a two-stage adaptive framework for progressively developing deep neural network (DNN) architectures that generalize well for a given training data set. In the first stage, a layerwise training approach is adopted where a new layer is added each time and trained independently by freezing parameters in the previous layers. We impose desirable structures on the DNN by employing manifold regularization, sparsity regularization, and physics-informed terms. We introduce a ε – δ – stability-promoting concept as a desirable property for a learning algorithm and show that employing manifold regularization yields a ε – δ stability-promoting algorithm. Further, we also derive the necessary conditions for the trainability of a newly added layer and investigate the training saturation problem. In the second stage of the algorithm (post-processing), a sequence of shallow networks is employed to extract information from the residual produced in the first stage, thereby improving the prediction accuracy. Numerical investigations on prototype regression and classification problems demonstrate that the proposed approach can outperform fully connected DNNs of the same size. Moreover, by equipping the physics-informed neural network (PINN) with the proposed adaptive architecture strategy to solve partial differential equations, we numerically show that adaptive PINNs not only are superior to standard PINNs but also produce interpretable hidden layers with provable stability. As a result, we also apply our architecture design strategy to solve inverse problems governed by elliptic partial differential equations.

42 ENGINEERING↗

NASA NTRS1974

Aeroelastic stability of coupled flap-lag motion of hingeless helicopter blades at arbitrary advance ratios

Equations for large amplitude coupled flap-lag motion of a hingeless elastic helicopter blade in forward flight are derived. Only a torsionally rigid blade excited by quasi-steady aerodynamic loads is considered. The effects of reversed flow together with some new terms due to radial flow are included. Using Galerkin's method the spatial dependence is eliminated and the equations are linearized about a suitable equilibrium position. The resulting system of homogeneous periodic equations is solved using multivariable Floquet-Liapunov theory, and the transition matrix at the end of the period is evaluated by two separate methods. Computational efficiency of the two numerical methods is compared. Results illustrating the effects of forward flight and various important blade parameters on the stability boundaries are presented.

Friedmann, P.↗

NASA NTRS1991

Controls for space structures

Assembly and operation of large space structures (LSS) in orbit will require robot-assisted docking and berthing of partially-assembled structures. These operations require new solutions to the problems of controls. This is true because of large transient and persistent disturbances, controller-structure interaction with unmodeled modes, poorly known structure parameters, slow actuator/sensor dynamical behavior, and excitation of nonlinear structure vibrations during control and assembly. For on-orbit assembly, controllers must start with finite element models of LSS and adapt on line to the best operating points, without compromising stability. This is not easy to do, since there are often unmodeled dynamic interactions between the controller and the structure. The indirect adaptive controllers are based on parameter estimation. Due to the large number of modes in LSS, this approach leads to very high-order control schemes with consequent poor stability and performance. In contrast, direct model reference adaptive controllers operate to force the LSS to track the desirable behavior of a chosen model. These schemes produce simple control algorithms which are easy to implement on line. One problem with their use for LSS has been that the model must be the same dimension as the LSS - i.e., quite large. A control theory based on the command generator tracker (CGT) ideas of Sobel, Mabins, Kaufman and Wen, Balas to obtain very low-order models based on adaptive algorithms was developed. Closed-loop stability for both finite element models and distributed parameter models of LSS was proved. In addition, successful numerical simulations on several LSS databases were obtained. An adaptive controller based on our theory was also implemented on a flexible robotic manipulator at Martin Marietta Astronautics. Computation schemes for controller-structure interaction with unmodeled modes, the residual mode filters or RMF, were developed. The RMF theory was modified to compensate slow actuator/sensor dynamics. These new ideas are being applied to LSS simulations to demonstrate the ease with which one can incorporate slow actuator/sensor effects into our design. It was also shown that residual mode filter compensation can be modified for small nonlinearities to produce exponentially stable closed-loop control. A theory for disturbance accommodating controllers based on reduced order models of structures was developed, and stability results for these controllers in closed-loop with large-scale finite element models of structures were obtained.

Balas, Mark↗

DOE OSTI2024

A second-order-in-time, explicit approach addressing the redundancy in the low-Mach, variable-density Navier-Stokes equations

A novel algorithm for explicit temporal discretization of the variable-density, low-Mach Navier-Stokes equations is presented here in this study. Recognizing there is a redundancy between the mass conservation equation, the equation of state, and the transport equation(s) for the scalar(s) which characterize the thermochemical state, and that it destabilizes explicit methods, we demonstrate how to analytically eliminate the redundancy and propose an iterative scheme to solve the resulting transformed scalar equations. The method obtains second-order accuracy in time regardless of the number of iterations, so one can terminate this subproblem once stability is achieved. Hence, flows with larger density ratios can be simulated while still retaining the efficiency, low cost, and parallelizability of an explicit scheme. The temporal discretization algorithm is used within a pseudospectral direct numerical simulation which extends the method of Kim, Moin, and Moser for incompressible flow to the variable-density, low-Mach setting, where we demonstrate stability for density ratios up to ~25.7.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Applications of atomic layer deposition and chemical vapor deposition for perovskite solar cells

Metal halide perovskite solar cells (PSCs) have rapidly evolved over the past decade to become a photovoltaic technology on the cusp of commercialization. In the process, numerous fabrication strategies have been explored with the goal of simultaneously optimizing for device efficiency, stability, and scalability. Atomic layer deposition (ALD) and chemical vapor deposition (CVD) have proven to be effective tools for the fabrication of various components of PSCs. This review article examines the application of ALD and CVD for the deposition and modification of charge transport layers, passivation layers, absorber materials, encapsulants, and electrodes. It outlines the use of these vapor deposition techniques in state-of-the-art, multi-junction solar cell devices, and also contains a discussion of the stability of metal halide perovskite materials under ALD and CVD conditions based on in situ characterization reported in literature. Finally, this article concludes with insights into future ALD and CVD research directions that could be undertaken to further aid the deployment of PSCs in emerging solar photovoltaic markets.

14 SOLAR ENERGY↗

NASA NTRS1963

INVESTIGATION OF PROBLEMS IN THERMAL CONVECTION

Numerical solution of thermal convection equations for velocity and temperature fluctuations as well as the stability analysis and heat transport for free boundary conditions

THERMAL CONVECTION↗