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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Parameter margins for stabilized conservative multilinear systems

Simple and elegant derivations of recent results concerning the computation of infinity-norm real-parameter margins for stabilized, mass-spring dynamical systems with the masses and the spring constraints as uncertain parameters are presented. The authors introduce the concept of critical frequency and gain for stabilized conservative systems whose uncertain parameters do not necessarily appear multilinearly in the numerator and denominator of the plant transfer function. An approach to parameter margin computation is presented.

Warren, Wayne↗
NASA NTRS1993

The effect of gravitational modulation on convection in vertical Bridgman growth

During vertical directional solidification of a binary alloy at constant velocity, buoyancy-driven solutal convection may occur due to the solute gradient associated with the solidification process. This problem is further complicated if time-periodic forcing is considered, which is relevant to materials processing in a microgravity environment or as a means of dynamic control of flow instabilities. The effect of time-periodic modulation is studied by introducing a gravitational acceleration which is a sinusoidal function of time. The onset of solutal convection is treated by a stability analysis of the linearized governing equations and boundary conditions. Solutions are obtained numerically employing two distinct computational implementations of Floquet theory. Results are presented for materials with large Schmidt number, and an analysis for large frequency and large Schmidt number yields a more complete description of the behavior in this relevant limit.

Murray, B. T.↗
NASA NTRS1992

The effect of gravitational modulation on convection in vertical Bridgman growth

During vertical directional solidification of a binary alloy at constant velocity, buoyancy driven solutal convection may occur due to the solute gradient associated with the solidification process. This problem is further complicated if time periodic forcing is considered, which is relevant to materials processing in a microgravity environment or as a means of dynamic control of flow instabilities. The effect of time periodic modulation is studied by introducing a gravitational acceleration which is a sinusoidal function of time. The onset of solutal convection is treated by a stability analysis of the linearized governing equations and boundary conditions. Solutions are obtained numerically by employing two distinct computational implementations of Floquet theory. Results for materials with large Schmidt number are presented, and an analysis for large frequency and large Schmidt number yields a more complete description of the behavior in this relevant limit.

Murray, Bruce T.↗
NASA NTRS1998

Research on the Problem of High-Precision Deployment for Large-Aperture Space-Based Science Instruments

The present paper summarizes results from an ongoing research program conducted jointly by the University of Colorado and NASA Langley Research Center since 1994. This program has resulted in general guidelines for the design of high-precision deployment mechanisms, and tests of prototype deployable structures incorporating these mechanisms have shown microdynamically stable behavior (i.e., dimensional stability to parts per million). These advancements have resulted from the identification of numerous heretofore unknown microdynamic and micromechanical response phenomena, and the development of new test techniques and instrumentation systems to interrogate these phenomena. In addition, recent tests have begun to interrogate nanomechanical response of materials and joints and have been used to develop an understanding of nonlinear nanodynamic behavior in microdynamically stable structures. The ultimate goal of these efforts is to enable nano-precision active control of micro-precision deployable structures (i.e., active control to a resolution of parts per billion).

Lake, Mark S.↗
NASA NTRS2003

Low Stretch Solid-Fuel Flame Transient Response to a Step Change in Gravity

The effect of a step change in gravity level on the stability of low stretch diffusion flames over a solid fuel is studied both numerically and experimentally. Drop tower experiments have been conducted in NASA Glenn Research Center's 5.2 Zero Gravity Facility. In the experiments burning PMMA cylinders, a dynamic transition is observed when the steadily burning 1g flame is dropped and becomes a 0g flame. To understand the physics behind this dynamic transition, a transient stagnation point model has been developed which includes gas-phase radiation and solid phase coupling to describe this dynamic process. In this paper, the experimental results are compared with the model predictions. Both model and experiment show that the interior of the solid phase does not have time to change significantly in the few seconds of drop time, so the experimental results are pseudo-steady in the gas-phase, but the solid is inherently unsteady over long time scales. The model is also used to examine the importance of fractional heat losses on extinction, which clearly demonstrates that as the feedback from the flame decreases, the importance of the ongoing heat losses becomes greater, and extinction is observed when these losses represent 80% or more of the flame feedback.

Armstrong, J. B.↗
NASA NTRS2012

Linear Estimation of Particle Bulk Parameters from Multi-Wavelength Lidar Measurements

An algorithm for linear estimation of aerosol bulk properties such as particle volume, effective radius and complex refractive index from multiwavelength lidar measurements is presented. The approach uses the fact that the total aerosol concentration can well be approximated as a linear combination of aerosol characteristics measured by multiwavelength lidar. Therefore, the aerosol concentration can be estimated from lidar measurements without the need to derive the size distribution, which entails more sophisticated procedures. The definition of the coefficients required for the linear estimates is based on an expansion of the particle size distribution in terms of the measurement kernels. Once the coefficients are established, the approach permits fast retrieval of aerosol bulk properties when compared with the full regularization technique. In addition, the straightforward estimation of bulk properties stabilizes the inversion making it more resistant to noise in the optical data. Numerical tests demonstrate that for data sets containing three aerosol backscattering and two extinction coefficients (so called 3 + 2 ) the uncertainties in the retrieval of particle volume and surface area are below 45% when input data random uncertainties are below 20 %. Moreover, using linear estimates allows reliable retrievals even when the number of input data is reduced. To evaluate the approach, the results obtained using this technique are compared with those based on the previously developed full inversion scheme that relies on the regularization procedure. Both techniques were applied to the data measured by multiwavelength lidar at NASA/GSFC. The results obtained with both methods using the same observations are in good agreement. At the same time, the high speed of the retrieval using linear estimates makes the method preferable for generating aerosol information from extended lidar observations. To demonstrate the efficiency of the method, an extended time series of observations acquired in Turkey in May 2010 was processed using the linear estimates technique permitting, for what we believe to be the first time, temporal-height distributions of particle parameters.

linear estimation↗
NASA NTRS2021

Positional Dependence of SNPP VIIRS Solar Diffuser BRDF Change Factor: An Empirical Approach

The Earth-observing Visible Infrared Imaging Radiometer Suite (VIIRS) on the Suomi National Polar-orbiting Partnership satellite regularly performs on-orbit radiometric calibration of its reflective solar bands (RSBs), primarily through observations of an onboard sunlit solar diffuser (SD). The on-orbit change of the SD bidirectional reflectance distribution function (BRDF) value, quantified by a numerical factor called the H-factor, is determined by the onboard SD stability monitor. Our previous study showed that the H-factor is solar angle and view direction dependent. In this study, we determine the dependence of the H-factor on the detector SD view footprint location. We fit an empirical model to the NASA Collection 1 SNPP VIIRS Level 1B (L1B) spectral reflectance difference across the detectors in an RSB over uniform Earth scenes of the Libya 4 desert and deep convective clouds. We apply the model predicted SD positional dependent H-factor to calibrate the RSBs. Under this new calibration scheme, the original unreal striping is removed from the homogeneous Libya 4 desert and the deep convective cloud images, as well as the original unreal striping from the Dunhuang desert image. The SD positional dependent H-factor has been used to calculate the SNPP VIIRS RSB radiometric correction factor for the NASA Collection 2.0 SNPP VIIRS L1B products.

SNPP VIIRS↗
NASA NTRS1983

Numerical simulation of electrothermal de-icing systems

Transient simulations of de-icing of composite aircraft components by electrothermal heating have been computed for both one and two-dimensional rectangular geometries. The implicit Crank-Nicolson formulation is used to insure stability of the finite-differenced heat conduction equations and the phase change in the ice layer is simulated using the Enthalpy method. Numerical solutions illustrating de-icer performance for various composite aircraft blades and environmental conditions are presented. Comparisons are made with previous studies and with available experimental data. Initial results using a coordinate mapping technique to describe the actual blade geometry are discussed.

De Witt, K. J.↗
DOE OSTI2021

Stabilization of kink/peeling modes by coupled rotation and ion diamagnetic drift effects in quiescent H-mode plasmas in DIII-D and JT-60U

Magnetohydrodynamic stability at the edge pedestal in several quiescent H-mode (QH-mode) plasmas in DIII-D and JT-60U experiments was analyzed by considering plasma rotation and ion diamagnetic drift effects. It was identified that a kink/peeling mode, which is a prime candidate for a trigger of edge harmonic oscillation in QH-mode, is stabilized by plasma rotation when considering the ion diamagnetic drift simultaneously in both experiments. Furthermore, the stabilizing effect by rotation becomes more effective in case using the rotation profile of the main ion species evaluated by assuming radial force balance. In addition, when inverting the rotation direction, it was found that the kink/peeling mode is more stabilized when considering the rotation of the main ion species, though the mode is less stabilized by the rotation of impurity ion species. The result implies that the kink/peeling mode stability in QH-mode plasmas is sensitive to how the rotation profile is evaluated, but it is shown that a qualitative trend stabilizing the kink/peeling mode by rotation can be reproduced with the poloidal rotation profile of an impurity predicted numerically based on the neoclassical theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1972

Slowly oscillating lifting surfaces at subsonic and supersonic speeds.

The paper presents an unsteady aerodynamic influence coefficient method based on the low-frequency approximation. The influence coefficients are of a type which have been used to compute steady flow about wing-body combinations; therefore, the new method may be extended readily to low-frequency unsteady flow about wing-body combinations. The validity of the method is demonstrated by comparisons with numerical results from conventional, unsteady lifting surface methods. The method is valid for arbitrary wings in supersonic flow and for wings of finite span in subsonic flow. The method, when extended to include wing-body-tail interactions, will have important applications for predicting stability, control, and gust response characteristics of large airplanes. Dynamic stability derivatives and pressure distributions are given for several planforms. The comparison with either analytical or other well established numerical methods shows good agreement.-

Brune, G. W.↗
NASA NTRS1992

Stability of compressible Taylor-Couette flow

The objectives of this paper are to: (1) develop both analytical and numerical tools that can be used to predict the onset of instability and subsequently to simulate the transition process by which the originally laminar flow evolves into a turbulent flow; and (2) conduct the preliminary investigations with the purpose of understanding the mechanisms of the vortical structures of the compressible flow between tow concentric cylinders.

Kao, K.↗
NASA NTRS1981

Numerical approximation of boundary conditions with applications to inviscid equations of gas dynamics

A comprehensive overview of the state of the art of well-posedness and stability analysis of difference approximations for initial boundary value problems of the hyperbolic type is presented. The applicability of recent theoretical development to practical calculations for nonlinear gas dynamics is examined. The one dimensional inviscid gas dynamics equations in conservation law form are selected for numerical experiments. The class of implicit schemes developed from linear multistep methods in ordinary differential equations is chosen and the use of linear extrapolation as an explicit or implicit boundary scheme is emphasized. Specification of boundary data in the primitive variables and computation in terms of the conservative variables in the interior is discussed. Some numerical examples for the quasi-one-dimensional nozzle are given.

Yee, H. C.↗
NASA NTRS1999

Navier-Stokes Analysis of the Flowfield Characteristics of an Ice Contaminated Aircraft Wing

An analytical study was performed as part of the NASA Lewis support of a National Transportation Safety Board (NTSB) aircraft accident investigation. The study was focused on the performance degradation associated with ice contamination on the wing of a commercial turbo-prop-powered aircraft. Based upon the results of an earlier numerical study conducted by the authors, a prominent ridged-ice formation on the subject aircraft wing was selected for detailed flow analysis using 2-dimensional (2-D), as well as, 3-dimensional (3-D) Navier-Stokes computations. This configuration was selected because it caused the largest lift decrease and drag increase among all the ice shapes investigated in the earlier study. A grid sensitivity test was performed to find out the influence of grid spacing on the lift, drag, and associated angle-of-attack for the maximum lift (C(sub lmax)). This study showed that grid resolution is important and a sensitivity analysis is an essential element of the process in order to assure that the final solution is independent of the grid. The 2-D results suggested that a severe stability and control difficulty could have occurred at a slightly higher angle-of-attack (AOA) than the one recorded by the Flight Data Recorder (FDR). This stability and control problem was thought to have resulted from a decreased differential lift on the wings with respect to the normal loading for the configuration. The analysis also indicated that this stability and control problem could have occurred whether or not natural ice shedding took place. Numerical results using an assumed 3-D ice shape showed an increase of the angle at which this phenomena occurred of about 4 degrees. As it occurred with the 2-D case, the trailing edge separation was observed but started only when the AOA was very close to the angle at which the maximum lift occurred.

Chung, J.↗
DOE OSTI2024

DG-IMEX method for a two-moment model for radiation transport in the $\mathscr{O}$($v$/$c$) limit

Here, we consider neutral particle systems described by moments of a phase-space density and propose a realizability-preserving numerical method to evolve a spectral two-moment model for particles interacting with a background fluid moving with nonrelativistic velocities. The system of nonlinear moment equations, with special relativistic corrections to $\mathscr{O}$($v$/$c$), expresses a balance between phase-space advection and collisions and includes velocity-dependent terms that account for spatial advection, Doppler shift, and angular aberration. The model is conservative for the correct $\mathscr{O}$($v$/$c$) Eulerian-frame number density and is consistent, to $\mathscr{O}$($v$/$c$), with Eulerian-frame energy and momentum conservation. This model is closely related to the one promoted by Lowrie et al. and similar to models currently used to study transport phenomena in large-scale simulations of astrophysical environments. The proposed numerical method is designed to preserve moment realizability, which guarantees that the moments correspond to a nonnegative phase-space density. The realizability-preserving scheme consists of the following key components: (i) a strong stability-preserving implicit-explicit (IMEX) time-integration method; (ii) a discontinuous Galerkin (DG) phase-space discretization with carefully constructed numerical uxes; (iii) a realizability-preserving implicit collision update; and(iv) a realizability-enforcing limiter. In time integration, nonlinearity of the moment model necessitates solution of nonlinear equations, which we formulate as fixed-point problems and solve with tailored iterative solvers that preserve moment realizability with guaranteed global convergence. We also analyze the simultaneous Eulerian-frame number and energy conservation properties of the semi-discrete DG scheme and propose a "spectral redistribution" scheme that promotes Eulerian-frame energy conservation. Through numerical experiments, we demonstrate the accuracy and robustness of this DG-IMEX method and investigate its Eulerian-frame energy conservation properties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1987

Aerothermal modeling program. Phase 2, element A: Improved numerical methods for turbulent viscous recirculating flows

The objective of this effort is to develop improved numerical schemes for predicting combustor flow fields. Various candidate numerical schemes were evaluated, and promising schemes were selected for detailed assessment. The criteria for evaluation included accuracy, computational efficiency, stability, and ease of extension to multidimensions. The candidate schemes were assessed against a variety of simple one- and two-dimensional problems. These results led to the selection of the following schemes for further evaluation: flux spline schemes (linear and cubic) and controlled numerical diffusion with internal feedback (CONDIF). The incorporation of the flux spline scheme and direct solution strategy in a computer program for three-dimensional flows is in progress.

Karki, K. C.↗
DOE OSTI2022

Application of linear prolongation to coarse mesh finite difference acceleration in CASMO5

The Coarse-Mesh Finite Difference (CMFD) method has been used for over a decade to accelerate the convergence of the Method of Characteristics (MOC) solution to the two- dimensional particle transport equation in CASMO5. Numerical testing, along with widespread use in production-level calculations, have shown that the current CMFD implementation provides stability and robustness for a wide range of realistic reactor physics problems. However, the recent development of linear prolongation has attracted attention from the community as a way to further improve the performance and stability of CMFD. Two interpolation methods for linear prolongation are presented in this work and implemented into a test version of CASMO5. The performance of the proposed interpolations, referred to as the bilinear and linear directional schemes, is evaluated in terms of runtime relative to the default constant or uniform scaling update. Numerical results indicate that the use of linear prolongation can reduce the transport solver runtime on average by approximately 10% when tested with two hundred randomly selected cases. The new directional linear interpolation, combined with default constant boundary updates, is found to provide the highest reduction in runtime for the cases analyzed. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Provably-Stable Overload Ride-Through Control for Grid-Forming Inverters Using System-Wide Lyapunov Function Analysis

A key challenge associated with a grid-forming (GFM) inverter based resource (IBR) is its behavior during severe grid disturbances: since a GFM inverter regulates voltage in the fast timescale instead of current or power, it may experience a transient overload of current, power and/or energy during a severe grid disturbance. While many promising control strategies for overload ride-through have been proposed over the past two decades, transient stability of the system during and after the transition to an overload ride-through control mode remains difficult to guarantee. In this work, a novel overload ride-through control strategy is proposed for a system of grid-forming inverters that takes both self-protection and system-wide transient stability into account. A proposed system-level supervisory control uses slow communication to pre-emptively assign a set of local ride through control parameters to individual GFM IBR, including a current-limiting virtual reactance, that guarantees that synchronism is still preserved for any set of anticipated grid disturbances. At the core of the supervisory control lies a Lyapunov-function-based routine capable of establishing a strong, albeit conservative, transient stability guarantee for the system. Here, the proposed overload ride-through control strategy is validated via numerical integration of a reduced-order model, as well as through detailed electromagnetic transient (EMT) simulation.

30 DIRECT ENERGY CONVERSION↗