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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An in-line burn-up algorithm: proposal and preliminary analysis

A new burn-up algorithm is proposed for coupling neutronics and isotopic depletion. The algorithm is analogous to those proposed for enforcing xenon equilibrium and converging neutronics and thermal-hydraulics. During the corrector step of standard predictor-corrector methods, the corrector nuclide density is updated as neutronics converges. This is expected to improve stability while reducing computational expense compared to other implicit algorithms. The stability of the algorithm is investigated by Fourier analysis for a simple case whose results have been successfully compared with numerical experiments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1979

Optical absorption and radiative heat transport in olivine at high temperature

Results are presented of measurements of the optical absorption spectra (300-8000 nm) of olivine as a function of temperature (300-1700 K) under conditions of controlled and known oxygen fugacity within the stability field of the samples. The absorption spectra are used to calculate the temperature-dependent radiative transfer coefficient of olivine and to numerically study the accuracy of the method. The present absorption measurements in olivine under oxidizing conditions known to be within the olivine stability field indicate that the effective radiative conductivity K(R) is lower than that obtained in previous studies under different experimental conditions. The lower value of K(R) makes it more likely that some of the earth's internal heat is removed by convection and less likely that thermal models involving conduction and radiation alone will satisfactorily explain thermal conditions in the earth's mantle.

Shankland, T. J.↗
NASA NTRS2007

3D Relativistic MHD Simulations of Magnetized Spine-Sheath Relativistic Jets

We have performed numerical simulations of weakly and strongly magnetized relativistic jets embedded in a weakly and strongly magnetized stationary or mildly relativistic (0.5c) sheath using the RAISHIN code. In the numerical simulations a jet with Lorentz factor gamma=2.5 is precessed to break the initial equilibrium configuration. Results of the numerical simulations are compared to theoretical predictions from a normal mode-analysis of the linearized RMHD equations describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized moving sheath. The prediction of increased stability of a weakly-magnetized system with mildly relativistic sheath flow to Kelvin-Helmholtz instabilities and the stabilization of a strongly-magnetized system with mildly relativistic sheath flow is confirmed by the numerical simulations.

Mizuno, Yosuke↗
NASA NTRS1986

Introduction to GRASP - General rotorcraft aeromechanical stability program - A modern approach to rotorcraft modeling

The General Rotorcraft Aeromechanical Stability Program (GRASP), which is a hybrid between finite element programs and spacecraft-oriented multibody programs, is described in terms of its design and capabilities. Numerical results from GRASP are presented and compared with the results from an existing, special-purpose coupled rotor/body aeromechanical stability program and with experimental data of Dowell and Traybar (1975 and 1977) for large deflections of an end-loaded cantilevered beam. The agreement is excellent in both cases.

Hodges, D. H.↗
NASA NTRS1985

Outdoor performance testing of thin-film devices

The Advanced Systems Research Group supports the photovoltaic advanced R&D (PV AR&D) project by providing outdoor (global) testing of PV cells, submodules, modules, and arrays. The group also provides in house engineering and analysis to identify and determine how technical issues such as cell/module/system adaptations, long term stability, reliability, economics, materials availability, safety, and environmental impacts affect the development and ultimate use of advanced PV thin film, innovative cell, and material technologies. A major thrust of the research effort is to develop and utilize instrumentation and procedures for monitoring and analyzing PV cells and submodules including outdoor performance and stability testing and life cycle accelerated stress testing. To accomplish the above, the solar energy research institute (SERI) outdoor PV test facility was established in 1982. The group has designed testing systems and analysis procedures for, and has tested, numerous amorphous silicon thin film submodules provided by SERI subcontractors and has performed long term outdoor stability tests on CdS/CuIr Se sub 2 and hydrogen passivated silicon solar cells. A significant contribution from this facility over the past year was the testing of large area amorphous silicon submodules.

Source record↗
NASA NTRS1984

Second- and third-order upwind difference schemes for hyperbolic conservation laws

Second- and third-order two time-level five-point explicit upwind-difference schemes are described for the numerical solution of hyperbolic systems of conservation laws and applied to the Euler equations of inviscid gas dynamics. Nonliner smoothing techniques are used to make the schemes total variation diminishing. In the method both hyperbolicity and conservation properties of the hyperbolic conservation laws are combined in a very natural way by introducing a normalized Jacobian matrix of the hyperbolic system. Entropy satisfying shock transition operators which are consistent with the upwind differencing are locally introduced when transonic shock transition is detected. Schemes thus constructed are suitable for shockcapturing calculations. The stability and the global order of accuracy of the proposed schemes are examined. Numerical experiments for the inviscid Burgers equation and the compressible Euler equations in one and two space dimensions involving various situations of aerodynamic interest are included and compared.

Yang, J. Y.↗
NASA NTRS1987

Development of a generic combustion stability code for liquid propellant rocket engines

The mathematical framework for a combustion stability analysis code is outlined. The goal for the code is to be general enough in problem treatment so that its validity and accuracy extend over a wide range of problem applications and that it lends the convenience for any future model improvement if necessary. An approach for modeling the combustion dynamics is devised to meet both requirements. An open-loop numerical procedure is also developed to mechanistically model various combustion processes for determining the stability parameters.

Fang, J. J.↗
NASA NTRS1990

Counterrotating prop-fan simulations which feature a relative-motion multiblock grid decomposition enabling arbitrary time-steps

Improvements are presented of a computer algorithm developed for the time-accurate flow analysis of rotating machines. The flow model is a finite volume method utilizing a high-resolution approximate Riemann solver for interface flux definitions. The numerical scheme is a block LU implicit iterative-refinement method which possesses apparent unconditional stability. Multiblock composite gridding is used to orderly partition the field into a specified arrangement of blocks exhibiting varying degrees of similarity. Block-block relative motion is achieved using local grid distortion to reduce grid skewness and accommodate arbitrary time step selection. A general high-order numerical scheme is applied to satisfy the geometric conservation law. An even-blade-count counterrotating unducted fan configuration is chosen for a computational study comparing solutions resulting from altering parameters such as time step size and iteration count. The solutions are compared with measured data.

Janus, J. Mark↗
NASA NTRS2000

Improving Tiltrotor Whirl-Mode Stability with Rotor Design Variations

Further increases in tiltrotor speeds are limited by coupled wing/rotor whirl-mode aeroelastic instability. Increased power, thrust, and rotor efficiency are not enough: the whirl-mode stability boundary must also be improved. With current technology, very stiff, thick wings of limited aspect ratio are essential to meet the stability requirements, which severely limits cruise efficiency and maximum speed. Larger and more efficient tiltrotors will need longer and lighter wings, for which whirl-mode flutter is a serious design issue. Numerous approaches to improving the whirl-mode airspeed boundary have been investigated, including tailored stiffness wings, active stability augmentation, variable geometry rotors, highly swept tips, and at one extreme, folding rotors. The research reported herein began with the much simpler approach of adjusting the chordwise positions of the rotor blade aerodynamic center and center of gravity, effected by offsetting the airfoil quarter chord or structural mass with respect to the elastic axis. The research was recently extended to include variations in blade sweep, control system stiffness, and pitch-flap coupling (delta(sub 3)). As an introduction to the subject, and to establish a baseline against which to measure stability improvements, this report will first summarize results. The paper will then discuss more advanced studies of swept blades and control-system modifications.

Acree, C. W., Jr.↗
NASA NTRS1994

Anisotropic Alfven-ballooning modes in Earth's magnetosphere

We have carried out a theoretical analysis of the stability and parallel structure of coupled shear Alfven and slow magnetosonic waves in Earth's inner magnetopause (i.e., at equatorial distances between about five and ten Earth radii) including effects of finite anisotropic Grad-Shafranov equation yields an approximate self-consistent magnetohydrodynamic (MHD) equilibrium. This MHD equilibrium is used in the numerical solution of a set of eigenmode equations which describe the field line eigenfrequency, linear stability, and parallel eigenmode structure. We call these modes anisotropic Alfven-ballooning modes. The main results are: (1) The field line eigenfrequency can be significantly lowered by finite pressure effects. (2) The parallel mode structure of the transverse wave components is fairly insensitive to changes in the plasma pressure, but the compressional magnetic component can become highly peaked near the magnetic equator as a result of increased pressure, especially when P(sub perpendicular to) is greater than P(sub parallel) (here P(sub perpendicular to) and P(sub parallel) are the perpendicular and parallel plasma pressure). (3) For the isotropic (P(sub parallel) = P(sub perpendicular to) = P) case ballooning instability can occur when the ratio of the plasma presure to the magnetic pressure, beta = P/(B squared/8 pi), exceeds a critical value beta(sup B)(sub 0) is approximately equal to 3.5 at the equator. (4) Compared to the isotropic case the critical beta value is lowered by anisotropy, either due to decreased field line bending stabilization when P(sub parallel) is greater than P(sub perpendicular to) or due to increased ballooning-mirror destabilization when P(sub perpendicular to) is greater than P(sub parallel). (5) We use a beta-delta stability diagram to display the regions of instability with respect to the equatorial values of the parameters bar beta and delta, where bar beta = (1/3)(beta(sub parallel) + 2 beta(sub perpendicular to)) is an average beta value and delta = 1 - P(sub parallel)/P(sub perpendicular to) is a measure of the plasma anisotropy. The diagram is divided into regions corresponding to the firehose, mirror and ballooning instabilities. It appears that observed values of the plasma pressure are below the critical value for the isotropic ballooning instability but it may be possible to approach a ballooning-mirror instability when P(sub perpendicular to)/P(sub parallel) is greater than or approximately 2.

Chan, Anthony A.↗
DOE OSTI2023

Reinforcement Learning of Structured Stabilizing Control for Linear Systems With Unknown State Matrix

This paper delves into designing feedback control gains for a continuous-time linear quadratic regulator (LQR) problem that is constrained to certain predefined structure with unknown state matrix. We bring forth the ideas from reinforcement learning (RL) in conjunction with sufficient stability and performance guarantees in order to design these structured gains using the trajectory measurements of states and controls. Here we first formulate a model-based framework using dynamic programming (DP) to embed the structural constraint to the LQR gain computation in the continuous-time setting, and then subsequently, formulate a policy iteration RL algorithm that can alleviate the requirement of known state matrix in conjunction with maintaining the feedback gain structure. The design enables a distributed learning control design which is necessary for many large-scale cyber-physical systems. Theoretical guarantees are provided for stability and convergence of the structured reinforcement learning (SRL) algorithm. We validate our theoretical results with numerical simulations on a multi-agent networked linear time-invariant (LTI) dynamic system.

42 ENGINEERING↗
NASA NTRS2012

Receptivity of Hypersonic Boundary Layers to Distributed Roughness and Acoustic Disturbances

Boundary-layer receptivity and stability of Mach 6 flow over smooth and rough 7 half-angle sharp-tipped cones are numerically investigated. The receptivity of the boundary layer to slow acoustic disturbances, fast acoustic disturbances, and vortical disturbances are considered. The effects of two-dimensional isolated and distributed roughness on the receptivity and stability are also simulated. The results show that the instability waves are generated in the leading edge region and that the boundary layer is much more receptive to slow acoustic waves than to the fast waves. Vortical disturbances also generate unstable second modes, however the receptivity coefficients are smaller than that of the slow acoustic wave. An isolated two-dimensional roughness element of height h/delta =1/4 did not produce any difference in the receptivity or in the stability of the boundary layer. Distributed roughness elements produced a small decrease in the receptivity coefficient and also stabilized the boundary layer by small amounts.

Balakumar, Ponnampalam↗
NASA NTRS1996

The effects of complex chemistry on triple flames

The structure, ignition, and stabilization mechanisms for a methanol (CH3OH)-air triple flame are studied using Direct Numerical Simulations (DNS). The methanol (CH3OH)-air triple flame is found to burn with an asymmetric shape due to the different chemical and transport processes characterizing the mixture. The excess fuel, methanol (CH3OH), on the rich premixed flame branch is replaced by more stable fuels CO and H2, which burn at the diffusion flame. On the lean premixed flame side, a higher concentration of O2 leaks through to the diffusion flame. The general structure of the triple point features the contribution of both differential diffusion of radicals and heat. A mixture fraction-temperature phase plane description of the triple flame structure is proposed to highlight some interesting features in partially premixed combustion. The effects of differential diffusion at the triple point add to the contribution of hydrodynamic effects in the stabilization of the triple flame. Differential diffusion effects are measured using two methods: a direct computation using diffusion velocities and an indirect computation based on the difference between the normalized mixture fractions of C and H. The mixture fraction approach does not clearly identify the effects of differential diffusion, in particular at the curved triple point, because of ambiguities in the contribution of carbon and hydrogen atoms' carrying species.

Echekki, T.↗
DOE OSTI2024

Enhanced Stability and Sensitivity for CA-125 Detection Under Microfluidic Shear Flow Using Polyethylene Glycol-Coated Biosensor

The microfluidic-based point-of-care (POC) diagnostic tool has garnered significant interest in recent years, offering rapid and cost-effective disease detection. There is a growing trend toward integrating microfluidic platforms with biosensors, aligning lab-on-a-chip technologies with POC diagnostic devices. Despite numerous efforts to incorporate biosensors into microfluidic systems, researchers have performed very limited investigations on the stability of biomarker detection when biosensors operate under microfluidic shear flow conditions. Gold nanoparticles (AuNPs) are a widely employed material in capacitive biosensors for antibody immobilization and sensitivity enhancement. However, AuNPs have limitations in providing stable detection of biomarkers within microfluidic shear flow due to their agglomeration nature. This study addresses these limitations by employing 2 kDa polyethylene glycol (PEG) as an intermediate biofunctional layer to immobilize CA-125 antibodies on gold-interdigitated electrodes for the stable and accurate detection of CA-125 antigens. The stabilities and sensitivities of AuNPs and PEG-coated biosensors are evaluated under both static drop and microfluidic shear flow conditions for CA-125 antigen detection. The experimental results demonstrate a capacitive signal response (5660 pF at 10 kHz) 2.2 times higher using the PEG-coated biosensor than the signal (2551 pF at 10 kHz) measured by the AuNP-coated biosensor in the detection of CA-125 antigen–antibody conjugation under static drop conditions, indicating the higher sensitivity of the PEG-coated biosensor. Additionally, the PEG-coated biosensor exhibits better consistency for the CA-125 antigen detection between static drop and microfluidic shear flow conditions (Cp decrease in percentage (ΔCp%↓) = 2.9% at 10 kHz) compared to the electrical signals measured using the AuNP-coated biosensor (ΔCp%↓ = 32.4% at 10 kHz), which suggests that the PEG-coated biosensor demonstrates higher stability for CA-125 antigen detection under microfluidic shear flow conditions. With these significant improvements brought by the PEG-coated biosensor, especially under microfluidic conditions, a substantial hurdle in developing electrical biosensors for POC diagnostic applications has been overcome, expediting further advancements in the field.

36 MATERIALS SCIENCE↗
NASA NTRS2001

Rotor Design Options for Improving V-22 Whirl-Mode Stability

With current tilt-rotor technology, very stiff, thick wings of limited aspect ratio are essential to meet aeroelastic stability requirements, which severely limits cruise efficiency and maximum speed. Larger and more efficient tiltrotors; will need longer and lighter wings, for which whirl-mode flutter is a serious design issue. Numerous approaches to improving the whirl-mode airspeed boundary have been investigated, typically focusing on wing design, active stability augmentation, and variable geometry rotors. The research to be reported applies the much simpler approaches of sweeping the outboard blade sections and extending balance weights at the tip. Improving proprotor whirl-mode stability margins is an ongoing research activity at NASA Ames Research Center. Previous efforts focussed on the XV-15 rotor. This is the first application of these design principles to a full V-22 rotor model in CAMRAD II (a previous publication presented results for a preliminary analytical model). The analysis to be discussed in the proposed paper includes a dual load-path blade model, a compressible-flow aerodynamic model, a representation of swashplate/rotor modal coupling, and other new features. The paper will include discussions of the mechanisms by which blade sweep and tip-mass offsets improve whirl-flutter margins. The V-22 was modeled with CAMRAD II; the rotor model is shown. It is also shows an example swept blade derived from the V-22 rotor system. For this rotor, the primary significance of sweep is the improved whirl-flutter boundary, not the reduced Mach-number effects. An offset tip mass is also shown.

Acree, C. W., Jr.↗
NASA NTRS2000

Effects of Rotor Design Variations on Tiltrotor Whirl-Mode Stability

Further increases in tiltrotor speeds are limited by coupled wing/rotor whirl-mode aeroelastic instability. Increased power, thrust, and efficiency are not enough: the whirl-mode stability boundary must also be improved. With current technology, very stiff, thick wings of limited aspect ratio are essential to meet the stability requirements, which severely limits cruise efficiency and maximum speed. Larger and more efficient tiltrotors will need longer and lighter wings, for which whirl-mode flutter is a serious design issue. Numerous approaches to improving the whirl-mode airspeed boundary have been investigated, typically focusing on wing design, active stability augmentation, and variable geometry rotors. The research reported herein applies the much simpler approach of sweeping the outboard blade sections. The effects-of control-system stiffness were also studied. Reference covers the first phase of the current research and discusses the evolution of the concept.

Acree, C. W., Jr.↗
DOE OSTI2020

Negative Charge Confined Amine Carriers within the Nanowire Network for Stable and Efficient Membrane Carbon Capture

Membrane-based carbon dioxide (CO 2 ) capture has attracted great research interest owing to its potential for higher separation efficiency and lower energy consumption. However, it is still a challenging task to capture CO 2 with membrane from flue gas, especially under moderate-temperature and high-humidity conditions. Here, a stable CO 2 -selective membrane that can operate at temperatures up to 90 °C and under high humidity is reported. The positively charged amine carriers for CO 2 are confined within the negatively charged polymer modified carbon nanotube (CNT) network. In this structure, interconnected CNTs act as the framework for the selective layer and provide numerous nanochannels for gas transport. The negatively charged polymer ensures the carrier stability and further regulates the size of nanochannels in the CNT network. By virtue of carrier-facilitated transport, high CO 2 permeance (up to 3300 gas permeation units) and high CO 2 /N 2 selectivity (400) are achieved under simulated flue gas conditions. Moreover, theoretical calculations verify that the stable separation performance is due to the strong electrostatic interaction between the amine carriers and polymer matrix. The high performance and good stability indicate the great potential of this novel membrane structure for practical application in CO 2 capture.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Operator inference with roll outs for learning reduced models from scarce and low-quality data

Data-driven modeling has become a key building block in computational science and engineering. However, data that are available in science and engineering are typically scarce, often polluted with noise and affected by measurement errors and other perturbations, which makes learning the dynamics of systems challenging. Here, in this work, we propose to combine data-driven modeling via operator inference with the dynamic training via roll outs of neural ordinary differential equations. Operator inference with roll outs inherits interpretability, scalability, and structure preservation of traditional operator inference while leveraging the dynamic training via roll outs over multiple time steps to increase stability and robustness for learning from low-quality and noisy data. Numerical experiments with data describing shallow water waves and surface quasi-geostrophic dynamics demonstrate that operator inference with roll outs provides predictive models from training trajectories even if data are sampled sparsely in time and polluted with noise of up to 10%.

97 MATHEMATICS AND COMPUTING↗