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,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

The development of an explicit thermochemical nonequilibrium algorithm and its application to compute three dimensional AFE flowfields

This study presents a three-dimensional explicit, finite-difference, shock-capturing numerical algorithm applied to viscous hypersonic flows in thermochemical nonequilibrium. The algorithm employs a two-temperature physical model. Equations governing the finite-rate chemical reactions are fully-coupled to the gas dynamic equations using a novel coupling technique. The new coupling method maintains stability in the explicit, finite-rate formulation while allowing relatively large global time steps. The code uses flux-vector accuracy. Comparisons with experimental data and other numerical computations verify the accuracy of the present method. The code is used to compute the three-dimensional flowfield over the Aeroassist Flight Experiment (AFE) vehicle at one of its trajectory points.

Palmer, Grant↗
NASA NTRS2004

Advanced Computational Modeling of Vapor Deposition in a High-pressure Reactor

In search of novel approaches to produce new materials for electro-optic technologies, advances have been achieved in the development of computer models for vapor deposition reactors in space. Numerical simulations are invaluable tools for costly and difficult processes, such as those experiments designed for high pressures and microgravity conditions. Indium nitride is a candidate compound for high-speed laser and photo diodes for optical communication system, as well as for semiconductor lasers operating into the blue and ultraviolet regions. But InN and other nitride compounds exhibit large thermal decomposition at its optimum growth temperature. In addition, epitaxy at lower temperatures and subatmospheric pressures incorporates indium droplets into the InN films. However, surface stabilization data indicate that InN could be grown at 900 K in high nitrogen pressures, and microgravity could provide laminar flow conditions. Numerical models for chemical vapor deposition have been developed, coupling complex chemical kinetics with fluid dynamic properties.

Cardelino, Beatriz H.↗
NASA NTRS2004

Advanced Computational Modeling of Vapor Deposition in a High-Pressure Reactor

In search of novel approaches to produce new materials for electro-optic technologies, advances have been achieved in the development of computer models for vapor deposition reactors in space. Numerical simulations are invaluable tools for costly and difficult processes, such as those experiments designed for high pressures and microgravity conditions. Indium nitride is a candidate compound for high-speed laser and photo diodes for optical communication system, as well as for semiconductor lasers operating into the blue and ultraviolet regions. But InN and other nitride compounds exhibit large thermal decomposition at its optimum growth temperature. In addition, epitaxy at lower temperatures and subatmospheric pressures incorporates indium droplets into the InN films. However, surface stabilization data indicate that InN could be grown at 900 K in high nitrogen pressures, and microgravity could provide laminar flow conditions. Numerical models for chemical vapor deposition have been developed, coupling complex chemical kinetics with fluid dynamic properties.

Cardelino, Beatriz H.↗
DOE OSTI2021

Impact of Increased Inverter Penetration on Power System Small-Signal Stability: Preprint

The transformation of the power system to include more distributed energy resources (DER) implies an increase in the number of inverter-based resources deployed on the grid. Envisioning future scenarios, this paper presents a small-signal stability analysis for a power grid comprising synchronous generators and inverter-based resources. Three types of inverter control are considered: grid following, droop-controlled grid forming, and virtual oscillator control grid forming. Although small-signal stability of power systems is a widely studied topic, systematic analysis of mixed machine-inverter systems with detailed control models at various inverter penetration levels are limited. This paper addresses the gap with numerical simulations tailored to the IEEE 39-bus system. Results show the system may become unstable at high inverter penetration level of grid-following inverters, and grid-forming inverter control can potentially improve system stability, thereby enabling very high level of DERs.

41 EE - Solar Energy Technologies Office (EE-4S)↗
NASA NTRS1985

Symposium on Numerical and Physical Aspects of Aerodynamic Flows, 3rd, California State University, Long Beach, CA, January 21-24, 1985, Proceedings

Among the topics discussed are hydrodynamic stability and turbulent transition, separated flow simulations employing finite difference methods, a quasi-simultaneous finite difference approach for strongly interacting flow, the Newton solution of coupled Euler and boundary layer equations, multielement airfoil aerodynamics, vorticity dynamics, the significance of the thin layer Navier-Stokes approximation, a comparison of interactive and Navier-Stokes calculations of separating boundary layer flows, and a Navier-Stokes analysis for the mechanism of shock wave/boundary layer interactions. Also discussed are transonic interactions of unsteady vortical flows, boundary layers on oscillating airfoils, a laser-interferometer skin friction meter, flow over surface-mounted semibluff bodies, computations for spinning and spin-stabilized supersonic projectiles, transition calculations in three-dimensional flows, the effect of projectile boattailing in transonic flow, flow characteristics of a body of revolution at incidence, the influence of the boundary layer on the leading edge flow of a swept wing, and an interactive scheme for three-dimensional transonic flows.

Source record↗
NASA NTRS1985

Effects of stratospheric lapse rate on numerically simulated thunderstorm cloud top structure

The three-dimensional anelastic model of Schlesinger (1975, 1978, 1984) is utilized to study the effects of stratospheric stability on thunderstorm cloud top structure. Three experiments were conducted; in the first the stratosphere is assumed isothermal, in the second test there is a temperature lapse of 3 K/km, and in the third case there is a 3 K/km inversion. Time variations of vertical velocity, mature storm structure, and cloud top are examined. It is observed that the cloud summit height is lowered as stratospheric stability increases, a cloud-top thermal couplet occurs in all cases, the cold region has a U or V shape, and there is a height minimum in the warm regions.

Schlesinger, Robert E.↗
NASA NTRS1989

Simulation of collisional transport processes and the stability of planetary rings

The utility of the phase-space fluid method for the study of planetary ring dynamics is presently demonstrated through the numerical solution of a model kinetic equation for a flattened Keplerian disk. Attention is given to ringlets composed of single-sized particles, as well as to ringlets composed of two different-sized particles; in the latter case, the ringlets evolve in such a way that the lighter particles are confined by the heavier ones. The results obtained indicate that some natural process may sharpen the optical depth profile of edges even without an external forcing mechanism, and that intermediate optical depths are dynamically preferred in some cases.

Brophy, Thomas G.↗
NASA NTRS1990

WKBJ approximation of the stability of a frontal mean state

The baroclinic instability of a frontal mean state is investigated using the WKBJ approximation. The results are compared with numerical calculations performed on the same mean state. Excellent agreement (within 5 percent) is found for jets whose half-width is as a Rossby radius of deformation. For jets 20 percent broader, the agreement is almost perfect.

Lindzen, Richard S.↗
NASA NTRS2003

Structure and Stability of One-Dimensional Detonations in Ethylene-Air Mixtures

The propagation of one-dimensional detonations in ethylene-air mixtures is investigated numerically by solving the one-dimensional Euler equations with detailed finite-rate chemistry. The numerical method is based on a second-order spatially accurate total-variation-diminishing scheme and a point implicit, first-order-accurate, time marching algorithm. The ethylene-air combustion is modeled with a 20-species, 36-step reaction mechanism. A multi-level, dynamically adaptive grid is utilized, in order to resolve the structure of the detonation. Parametric studies over an equivalence ratio range of 0.5 less than phi less than 3 for different initial pressures and degrees of detonation overdrive demonstrate that the detonation is unstable for low degrees of overdrive, but the dynamics of wave propagation varies with fuel-air equivalence ratio. For equivalence ratios less than approximately 1.2 the detonation exhibits a short-period oscillatory mode, characterized by high-frequency, low-amplitude waves. Richer mixtures (phi greater than 1.2) exhibit a low-frequency mode that includes large fluctuations in the detonation wave speed; that is, a galloping propagation mode is established. At high degrees of overdrive, stable detonation wave propagation is obtained. A modified McVey-Toong short-period wave-interaction theory is in excellent agreement with the numerical simulations.

Yungster, S.↗
DOE OSTI2021

Resistive wall mode stability and resonant field amplification in MAST high beta plasma

The n=1 (n is the toroidal mode number) resistive wall mode (RWM) instability and the resonant field amplification (RFA) due to a stable RWM response are numerically investigated for a MAST high-pressure plasma scenario, utilizing the MARS-F (Liu Y.Q. et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu Y.Q. et al 2008 Phys. Plasmas 15 112503) codes, for the purpose of understanding the potential RWM behavior in MAST-U plasmas. Ideal kink stability analysis identifies a target plasma, with parameters similar to that of the reference equilibrium reconstructed from a MAST high-pressure discharge, that accesses the RWM regime. The unstable n=1 RWM for the target equilibrium is subject to strong damping by the plasma toroidal flow and/or the drift kinetic effects from thermal particles. As a result, the mode is found to be stable under the experimental flow conditions, or even without flow stabilization if drift kinetic stabilization is included. The stability prediction is robust against variation of the assumed resistive wall minor radius. Active magneto-hydrodynamic (MHD) spectroscopy modeling, using the magnetic coils designed for controlling the edge localized modes (ELMs) in MAST-U as the antenna, shows strong resonant field amplification due to a stable RWM response in the target plasma. Maximal amplification, of 6-7 times larger than the applied vacuum field within the plasma, is obtained assuming a coil phasing of 120-180 degrees between the upper and lower rows of the ELM control coils. The MHD-kinetic hybrid model predicts generally higher RFA than the fluid model, in particular near the low-field and high-field sides of the plasma boundary surface. Furthermore, the MARS-F/K modeling thus shows that the RWM in MAST-U is subject to strong flow and kinetic stabilization, but can nevertheless be detected by active MHD spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Numerical simulation of two-dimensional heat transfer in composite bodies with application to de-icing of aircraft components

Transient, numerical simulations of the de-icing of composite aircraft components by electrothermal heating were performed for a two dimensional rectangular geometry. The implicit Crank-Nicolson formulation was used to insure stability of the finite-difference heat conduction equations and the phase change in the ice layer was simulated using the Enthalpy method. The Gauss-Seidel point iterative method was used to solve the system of difference equations. Numerical solutions illustrating de-icer performance for various composite aircraft structures and environmental conditions are presented. Comparisons are made with previous studies. The simulation can also be used to solve a variety of other heat conduction problems involving composite bodies.

Chao, D. F. K.↗
DOE OSTI2023

Direct Numerical Simulation of Partial Fuel Stratification Assisted Lean Premixed Combustion for Assessment of Hybrid G-Equation/Well-Stirred Reactor Model

Partial fuel stratification (PFS) is a promising fuel injection strategy to stabilize lean premixed combustion in spark-ignition (SI) engines. PFS creates a locally stratified mixture by injecting a fraction of the fuel, just before spark timing, into the engine cylinder containing homogeneous lean fuel/air mixture. Further, this locally stratified mixture, when ignited, results in complex flame structure and propagation modes similar to partially premixed flames and allows for faster and more stable flame propagation than a homogeneous lean mixture. This study focuses on understanding the detailed flame structures associated with PFS-assisted lean premixed combustion. First, a two-dimensional direct numerical simulation (DNS) is performed using detailed fuel chemistry, experimental pressure trace, and realistic initial conditions mapped from a prior engine large-eddy simulation (LES), replicating practical lean SI operating conditions. DNS results suggest that the conventional triple flame structure is prevalent during the initial stage of flame kernel growth. Both premixed and nonpremixed combustion modes are present with the premixed mode contributing dominantly to the total heat release. Detailed analysis further reveals the effects of flame stretch and fuel pyrolysis on flame displacement speed. Based on the DNS findings, the accuracy of a hybrid G-equation/well-stirred reactor (WSR) combustion model is assessed for the PFS-assisted lean operation in the LES context. It is found that the G-equation model qualitatively captures the premixed branches of the triple flame, while the WSR model predicts the nonpremixed branch of the triple flame. Finally, potential needs for improvements to the hybrid G-equation/WSR modeling approach are discussed.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1990

A computer module used to calculate the horizontal control surface size of a conceptual aircraft design

The creation of a computer module used to calculate the size of the horizontal control surfaces of a conceptual aircraft design is discussed. The control surface size is determined by first calculating the size needed to rotate the aircraft during takeoff, and, second, by determining if the calculated size is large enough to maintain stability of the aircraft throughout any specified mission. The tail size needed to rotate during takeoff is calculated from a summation of forces about the main landing gear of the aircraft. The stability of the aircraft is determined from a summation of forces about the center of gravity during different phases of the aircraft's flight. Included in the horizontal control surface analysis are: downwash effects on an aft tail, upwash effects on a forward canard, and effects due to flight in close proximity to the ground. Comparisons of production aircraft with numerical models show good accuracy for control surface sizing. A modified canard design verified the accuracy of the module for canard configurations. Added to this stability and control module is a subroutine that determines one of the three design variables, for a stable vectored thrust aircraft. These include forward thrust nozzle position, aft thrust nozzle angle, and forward thrust split.

Sandlin, Doral R.↗
DOE OSTI2025

Approximating the particle distribution in rotating and tandem mirror traps

Steady-state distribution functions can be used to calculate stability conditions for modes, radiation energy losses and particle loss rates. Heuristic analytic approximations to these distributions can capture key behaviors of the true distributions such as the relative speeds of different transport processes while possessing computational advantages over their numerical counterparts. In this paper, we motivate and present a closed-form ana- lytic model for a distribution of particles in a centrifugal or tandem mirror. We find that our model outperforms other known models in approximating numerical steady- state simulations outside of a narrow range of low confining potentials. We demonstrate the model’s suitability in the high confining potential regime for applications such as loss-cone stability thresholds, fusion yields and available energy.

Li, G.X.↗
DOE OSTI2022

An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD

The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Additionally, therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in Chacón et al. (2002) to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime (≥ 10 7 ) with the successful resolution of plasmoid instabilities and thin current sheets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2003

Identifying Head-Trunk and Lower Limb Contributions to Gaze Stabilization During Locomotion

The goal of the present study was to determine how the multiple, interdependent full-body sensorimotor subsystems respond to a change in gaze stabilization task constraints during locomotion. Nine subjects performed two gaze stabilization tasks while walking at 6.4 km/hr on a motorized treadmill: 1) focusing on a central point target; 2) reading numeral characters; both presented at 2m in front at the level of their eyes. While subjects performed the tasks we measured: temporal parameters of gait, full body sagittal plane segmental kinematics of the head, trunk, thigh, shank and foot, accelerations along the vertical axis at the head and the shank, and the vertical forces acting on the support surface. We tested the hypothesis that with the increased demands placed on visual acuity during the number recognition task, subjects would modify full-body segmental kinematics in order to reduce perturbations to the head in order to successfully perform the task. We found that while reading numeral characters as - compared to the central point target: 1) compensatory head pitch movement was on average 22% greater despite the fact that the trunk pitch and trunk vertical translation movement control were not significantly changed; 2) coordination patterns between head and trunk as reflected by the peak cross correlation between the head pitch and trunk pitch motion as well as the peak cross correlation between the head pitch and vertical trunk translation motion were not significantly changed; 3) knee joint total movement was on average 11% greater during the period from the heel strike event to the peak knee flexion event in stance phase of the gait cycle; 4) peak acceleration measured at the head was significantly reduced by an average of 13% in four of the six subjects. This was so even when the peak acceleration at the shank and the transmissibility of the shock wave at heel strike (measured by the peak acceleration ratio of the head/shank) remained unchanged. Taken together these results provide further evidence that the full body contributes to gaze stabilization during locomotion, and that its different functional elements can be modified online to contribute to gaze stabilization for different visual task constraints.

Mulavara, Ajitkumar P.↗
DOE OSTI2021

A duality-based coupling of Cosserat crystal plasticity and phase field theories for modeling grain refinement

High-rate deformation processes of metals entail intense grain refinement and special attention needs to be paid to capture the evolution of microstructure. In this article, a new formulation for coupling Cosserat crystal plasticity and phase field is developed. A common approach is to penalize kinematic incompatibility between lattice orientation and displacement-based elastic rotation. However, this can lead to significant solution sensitivity to the penalty parameter, resulting in low accuracy and convergence rates. To address these issues, a duality-based formulation is developed which directly imposes the rotational kinematic compatibility. A weak inf-sup-based skew-symmetric stress projection is introduced to suppress instabilities present in the dual formulation. An additional least squares stabilization is introduced to suppress the spurious lattice rotation with a suitable parameter range derived analytically and validated numerically. The required high-order continuity is attained by the reproducing kernel approximation. It is observed that equal order displacement-rotation-phase field approximations are stable, which allows efficient employment of the same set of shape functions for all independent variables. The proposed formulation is shown to yield superior accuracy and convergence with marginal parameter sensitivity compared to the penalty-based approach and successfully captures the dominant rotational recrystallization mechanism including block dislocation structures and grain boundary migration.

42 ENGINEERING↗
DOE OSTI2021

Structural stability and artificial buckling modes in topology optimization

Abstract This paper demonstrates how a strain energy transition approach can be used to remove artificial buckling modes that often occur in stability constrained topology optimization problems. To simulate the structural response, a nonlinear large deformation hyperelastic simulation is performed, wherein the fundamental load path is traversed using Newton’s method and the critical buckling load levels are estimated by an eigenvalue analysis. The goal of the optimization is to minimize displacement, subject to constraints on the lowest critical buckling loads and maximum volume. The topology optimization problem is regularized via the Helmholtz PDE-filter and the method of moving asymptotes is used to update the design. The stability and sensitivity analyses are outlined in detail. The effectiveness of the energy transition scheme is demonstrated in numerical examples.

Dalklint, Anna (ORCID:0000000346195205)↗