SEARCH · Engineering Papers

Results for “stiff problem”

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 631 records · Page 35

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Symbolic construction of the chemical Jacobian of quasi-steady state (QSS) chemistries for Exascale computing platforms

The Quasi-Steady State Approximation (QSSA) can be an effective tool for reducing the size and stiffness of chemical mechanisms for implementation in computational reacting flow solvers. However, for many applications, the resulting model still requires implicit methods for efficient time integration. Here, in this paper, we outline an approach to formulating the QSSA reduction that is coupled with a strategy to generate C++ source code to evaluate the net species production rates, and the chemical Jacobian. The code-generation component employs a symbolic approach enabling a simple and effective strategy to analytically compute the chemical Jacobian. For computational tractability, the symbolic approach needs to be paired with common subexpression elimination which can negatively affect memory usage. Several solutions are outlined and successfully tested on a 3D multipulse ignition problem, thus allowing portable application across chemical model sizes and GPU capabilities. The implementation of the proposed method is available at https://github.com/AMReX-Combustion/PelePhysics under an open-source license.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1978

Solution of the subsynchronous whirl problem in the high-pressure hydrogen turbomachinery of the Space Shuttle Main Engine

Subsynchronous whirl of the high-pressure fuel turbopump limited operation of the Space Shuttle Main Engine for some months in early 1976. The means by which this problem was successfully attacked is of particular interest to the rotor-dynamics community, not only because this machinery was the highest power-to-weight ratio known (77,000 hp...760 pounds), but because of the multiple forcing functions involved and the means, both analytical and experimental, which were utilized in separating variables, pointing toward successful solutions, and evaluating results. The general means of identifying fundamental characteristics, analyzing data, and conducting computer investigations are delineated. The results of analysis and testing are discussed. Since whirl inception occurs at a shaft speed greater than twice the first system critical, this was increased by stiffening the shaft and bearing supports; the system damping and system stiffness was additionally increased by proper interstage seal design to the extent that the instability threshold is now beyond the operating range.

Ek, M. C.↗

NASA NTRS1989

A magnetic damper for first mode vibration reduction in multimass flexible rotors

Many rotating machines such as compressors, turbines and pumps have long thin shafts with resulting vibration problems, and would benefit from additional damping near the center of the shaft. Magnetic dampers have the potential to be employed in these machines because they can operate in the working fluid environment unlike conventional bearings. An experimental test rig is described which was set up with a long thin shaft and several masses to represent a flexible shaft machine. An active magnetic damper was placed in three locations: near the midspan, near one end disk, and close to the bearing. With typical control parameter settings, the midspan location reduced the first mode vibration 82 percent, the disk location reduced it 75 percent and the bearing location attained a 74 percent reduction. Magnetic damper stiffness and damping values used to obtain these reductions were only a few percent of the bearing stiffness and damping values. A theoretical model of both the rotor and the damper was developed and compared to the measured results. The agreement was good.

Kasarda, M. E. F.↗

DOE OSTI2021

Component-wise reduced order model lattice-type structure design

Lattice-type structures can provide a combination of stiffness with light weight that is desirable in a variety of applications. Design optimization of these structures must rely on approximations of the governing physics to render solution of a mathematical model feasible. In this paper, we propose a topology optimization (TO) formulation that approximates the governing physics using component-wise reduced order modeling as introduced in Huynh et al. (2013); Eftang and Patera (2013), which can reduce solution time by multiple orders of magnitude over a full-order finite element model while providing a relative error in the solution of 1%. In addition, the offline training data set from such component-wise models is reusable, allowing its application to many design problems for only the cost of a single offline training phase, and the component-wise method is nearly embarrassingly parallel. We also show how the parameterization chosen in our optimization allows a simplification of the component-wise reduced order model (CWROM) not noted in previous literature, for further speedup of the optimization process. Furthermore, the sensitivity of the compliance with respect to the particular parameterization is derived solely at the component level. In numerical examples, we demonstrate a 1000x speedup over a full-order FEM model with relative error of 1% and show minimum compliance designs for two different cantilever beam examples, one smaller and one larger. Finally, error bounds for the displacement field, compliance, and compliance sensitivity of the CWROM are derived.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1991

On the application of ENO scheme with subcell resolution to conservation laws with stiff source terms

Two approaches are used to extend the essentially non-oscillatory (ENO) schemes to treat conservation laws with stiff source terms. One approach is the application of the Strang time-splitting method. Here the basic ENO scheme and the Harten modification using subcell resolution (SR), ENO/SR scheme, are extended this way. The other approach is a direct method and a modification of the ENO/SR. Here the technique of ENO reconstruction with subcell resolution is used to locate the discontinuity within a cell and the time evolution is then accomplished by solving the differential equation along characteristics locally and advancing in the characteristic direction. This scheme is denoted ENO/SRCD (subcell resolution - characteristic direction). All the schemes are tested on the equation of LeVeque and Yee (NASA-TM-100075, 1988) modeling reacting flow problems. Numerical results show that these schemes handle this intriguing model problem very well, especially with ENO/SRCD which produces perfect resolution at the discontinuity.

Chang, Shih-Hung↗

NASA NTRS2005

Tools for Designing and Analyzing Structures

Structural Design and Analysis Toolset is a collection of approximately 26 Microsoft Excel spreadsheet programs, each of which performs calculations within a different subdiscipline of structural design and analysis. These programs present input and output data in user-friendly, menu-driven formats. Although these programs cannot solve complex cases like those treated by larger finite element codes, these programs do yield quick solutions to numerous common problems more rapidly than the finite element codes, thereby making it possible to quickly perform multiple preliminary analyses - e.g., to establish approximate limits prior to detailed analyses by the larger finite element codes. These programs perform different types of calculations, as follows: 1. determination of geometric properties for a variety of standard structural components; 2. analysis of static, vibrational, and thermal- gradient loads and deflections in certain structures (mostly beams and, in the case of thermal-gradients, mirrors); 3. kinetic energies of fans; 4. detailed analysis of stress and buckling in beams, plates, columns, and a variety of shell structures; and 5. temperature dependent properties of materials, including figures of merit that characterize strength, stiffness, and deformation response to thermal gradients

Luz, Paul L.↗

NASA NTRS1985

Rotor-to-stator rub vibration in centrifugal compressor

One example of excessive vibration encountered during loading of a centrifugal compressor train (H type compressor with HP casing) is discussed. An investigation was made of the effects of the dynamic load on the bearing stiffness and the rotor-bearing system critical speed. The high vibration occurred at a "threshold load," but the machine didn't run smoothly due to rubs even when it had passed through the threshold load. The acquisition and discussion of the data taken in the field as well as a description of the case history which utilizes background information to identify the malfunction conditions is presented. The analysis shows that the failures, including full reverse precession rub and exact one half subharmonic vibration, were caused by the oversize bearings and displacement of the rotor center due to foundation deformation and misalignment between gear shafts, etc. The corrective actions taken to alleviate excessive vibration and the problems which remain to be solved are also presented.

Gao, J. J.↗

NASA NTRS2022

Mars 2020 Perseverance Rover SHERLOC Instrument Isolation System

The NASA Jet Propulsion Laboratory (JPL) successfully landed the Mars 2020 Perseverance Rover at Jezero Crater on the Martian surface. Perseverance’s main mission objective is to cache Martian rock and regolith samples in hermetically sealed tubes to be brought back to Earth by future missions. To achieve this goal the rover is equipped with a 2-meter-long Robotic Arm (RA) which manipulates the Turret to interact with the surface. The Turret is comprised of science instruments and tools that enable surface sampling and science. The backbone of the Turret hardware is the Rotary Percussive Coring Drill. The science instruments, which are directly mounted to the structural housing of this drill, must be able to withstand not only the launch and entry, descent and landing (EDL) loads of the vehicle but also the dynamic percussive environment, induced by the drill throughout surface operations. This paper discusses the technical hardware design, development and testing of a vibration isolation system to protect the Scanning Habitable Environments with Raman and Luminescence for Organics and Chemicals (SHERLOC) Instrument from these dynamic environments. The hardware design aspect of this problem was limited by both size and mass constraints. The Turret is tightly configured and the physical space in which hardware can reside is restricted due to multiple axes of surface interaction points. Another driving requirement is the extreme range of thermal non-operational environment from -135C to +90C. This, in addition to SHERLOC being a heavier instrument than previously integrated on a turret, made the use of isolation systems employed on previous rover missions problematic. The final design incorporates custom wire mesh springs preloaded within titanium hexapod struts equipped with flexure ends. Custom versions of commercially available wire mesh springs were designed and tested. These custom springs differed in dimension, density and stiffness from the off the shelf options. The newly designed springs went through various stages of flight acceptance testing. Detailed physics-based modeling to understand what loads the isolation system would experience was used to define the test conditions for the springs and entire isolation subsystem. The testing program developed around this was extensive to fully characterize the isolation performance over cleanliness levels, temperature and life. Another major aspect of this work was developing a method of replicating the drill’s percussive environment for performance and life testing of Turret mounted hardware. Attaching the instrument to an operating drill for testing was not a viable option. A method of using a Highly Accelerated Life Test (HALT) table was developed and applied across the project as the most representative and tunable option to physically simulate this new environment. This method, along with launch environment random vibration testing, was used to build an entire dynamic qualification and flight hardware test program. The results of the design, modeling, testing and analysis of the SHERLOC Isolation System is described throughout the following paper.

Krafchak, Todd↗

NASA NTRS1994

Numerical Simulations of High-Speed Chemically Reacting Flow

The Essentially NonOscillatory (ENO) shock-capturing scheme for the solution of hyperbolic equations is extended to solve a system of coupled conservation equations governing two-dimensional, time-dependent, compressible chemically reacting flow with full chemistry. The thermodynamic properties of the mixture are modeled accurately, and stiff kinetic terms are separated from the fluid motion by a fractional step algorithm. The methodology is used to study the concept of shock-induced mixing and combustion, a process by which the interaction of a shock wave with a jet of low-density hydrogen fuel enhances mixing through streamwise vorticity generation. Test cases with and without chemical reaction are explored here. Our results indicate that, in the temperature range examined, vorticity generation as well as the distribution of atomic species do not change significantly with the introduction of a chemical reaction and subsequent heat release. The actual diffusion of hydrogen is also relatively unaffected by the reaction process. This suggests that the fluid mechanics of this problem may be successfully decoupled from the combustion processes, and that computation of the mixing problem (without combustion chemistry) can elucidate much of the important physical features of the flow.

Ton, V. T.↗

NASA NTRS1994

Numerical Simulations of High-Speed Chemically Reacting Flow

The essentially nonoscillatory (ENO) shock-capturing scheme for the solution of hyperbolic equations is extended to solve a system of coupled conservation equations governing two-dimensional, time-dependent, compressible chemically reacting flow with full chemistry. The thermodynamic properties of the mixture are modeled accurately, and stiff kinetic terms are separated from the fluid motion by a fractional step algorithm. The methodology is used to study the concept of shock-induced mixing and combustion, a process by which the interaction of a shock wave with a jet of low-density hydrogen fuel enhances mixing through streamwise vorticity generation. Test cases with and without chemical reaction are explored here. Our results indicate that, in the temperature range examined, vorticity generation as well as the distribution of atomic species do not change significantly with the introduction of a chemical reaction and subsequent heat release. The actual diffusion of hydrogen is also relatively unaffected by the reaction process. This suggests that the fluid mechanics of this problem may be successfully decoupled from the combustion processes, and that computation of the mixing problem (without combustion chemistry) can elucidate much of the important physical features of the flow.

Ton, V. T.↗

NASA NTRS2016

Analysis and Characterization of Damage Utilizing an Orthotropic Generalized Composite Material Model Suitable for Use in Impact Problems

The need for accurate material models to simulate the deformation, damage and failure of polymer matrix composites under impact conditions is becoming critical as these materials are gaining increased usage in the aerospace and automotive communities. In order to address a series of issues identified by the aerospace community as being desirable to include in a next generation composite impact model, an orthotropic, macroscopic constitutive model incorporating both plasticity and damage suitable for implementation within the commercial LS-DYNA computer code is being developed. The plasticity model is based on extending the Tsai-Wu composite failure model into a strain hardening-based orthotropic plasticity model with a non-associative flow rule. The evolution of the yield surface is determined based on tabulated stress-strain curves in the various normal and shear directions and is tracked using the effective plastic strain. To compute the evolution of damage, a strain equivalent semi-coupled formulation is used in which a load in one direction results in a stiffness reduction in multiple material coordinate directions. A detailed analysis is carried out to ensure that the strain equivalence assumption is appropriate for the derived plasticity and damage formulations that are employed in the current model. Procedures to develop the appropriate input curves for the damage model are presented and the process required to develop an appropriate characterization test matrix is discussed.

Impact↗

NASA NTRS2016

Analysis and Characterization of Damage Utilizing an Orthotropic Generalized Composite Material Model Suitable for Use in Impact Problems

Impact↗

NASA NTRS1956

The Dynamic Response of a Large Airplane to Continuous Random Atmospheric Disturbances

The statistical approach to the gust-loads problem, which consists in considering flight th1ough turbulent air to be a stationary random process, is extended by including the effect of lateral variations of the instantaneous gust intensity on the aerodynamic forces and on the resultant motions and stresses of rigid and flexible airplanes. By means of some calculations of normal and rolling accelerations, as well as of the root bending moment, it is shown that these effects may be significant for large airplanes.

Aeroplane↗

DOE OSTI2022

High Order Strong Stability Preserving MultiDerivative Implicit and IMEX Runge--Kutta Methods with Asymptotic Preserving Properties

In this article we present a class of high order unconditionally strong stability preserving (SSP) implicit two-derivative Runge--Kutta schemes and SSP implicit-explicit (IMEX) multi-derivative Runge--Kutta schemes where the time-step restriction is independent of the stiff term. The unconditional SSP property for a method of order $p>2$ is unique among SSP methods and depends on a backward-in-time assumption on the derivative of the operator. We show that this backward derivative condition is satisfied in many relevant cases where SSP IMEX schemes are desired. We devise unconditionally SSP implicit Runge--Kutta schemes of order up to $p=4$ and IMEX Runge--Kutta schemes of order up to $p=3$. For the multiderivative IMEX schemes, we also derive and present the order conditions, which have not appeared previously. The unconditional SSP condition ensures that these methods are positivity preserving, and we present sufficient conditions under which such methods are also asymptotic preserving when applied to a range of problems, including a hyperbolic relaxation system, the Broadwell model, and the Bhatnagar--Gross--Krook kinetic equation. We present numerical results to support the theoretical results on a variety of problems.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1980

Thermomechanical phenomena in high speed rubbing

An analytical approach is presented for the modeling of the thermomechanical interactions which occur in high speed sliding situations. These sliding contact problems which are characterized by active and interrelated thermal and mechanical phenomena could be called 'rub energetics' problems. Analytical models were developed to simulate two different rub situations: high energy braking of disk brakes and high speed rubs of gas path seals in turbine engines. The models proved to be particularly useful in predicting the severe temperatures and deformations near hot contact patches on the rubbing surfaces. The size of the hot patches is generally determined by normal load and the properties of the contacting materials. Temperatures at the contact patches can approach the melting point of the materials, especially at high sliding velocities. These high temperatures can lead to large amounts of near-surface deformation and high wear rates. Decreased contact temperatures can result from using materials with increased thermal conductivity and increased heat capacity or choosing mechanical properties (decreased stiffness, yield stress or coefficient of thermal expansion) which give larger hot spot size.

Kennedy, F. E.↗

NASA NTRS1998

Current and Future Critical Issues in Rocket Propulsion Systems

The objective of this research was to tackle several problems that are currently of great importance to NASA. In a liquid rocket engine several complex processes take place that are not thoroughly understood. Droplet evaporation, turbulence, finite rate chemistry, instability, and injection/atomization phenomena are some of the critical issues being encountered in a liquid rocket engine environment. Pulse Detonation Engines (PDE) performance, combustion chamber instability analysis, 60K motor flowfield pattern from hydrocarbon fuel combustion, and 3D flowfield analysis for the Combined Cycle engine were of special interest to NASA. During the summer of 1997, we made an attempt to generate computational results for all of the above problems and shed some light on understanding some of the complex physical phenomena. For this purpose, the Liquid Thrust Chamber Performance (LTCP) code, mainly designed for liquid rocket engine applications, was utilized. The following test cases were considered: (1) Characterization of a detonation wave in a Pulse Detonation Tube; (2) 60K Motor wall temperature studies; (3) Propagation of a pressure pulse in a combustion chamber (under single and two-phase flow conditions); (4) Transonic region flowfield analysis affected by viscous effects; (5) Exploring the viscous differences between a smooth and a corrugated wall; and (6) 3D thrust chamber flowfield analysis of the Combined Cycle engine. It was shown that the LTCP-2D and LTCP-3D codes are capable of solving complex and stiff conservation equations for gaseous and droplet phases in a very robust and efficient manner. These codes can be run on a workstation and personal computers (PC's).

Navaz, Homayun K.↗

DOE OSTI2026

Asymptotic-preserving semi-implicit finite volume scheme for extended magnetohydrodynamics

A Finite Volume (FV) scheme is developed for solving the extended magnetohydrodynamic (XMHD) equations, yielding accurate results in the ideal, resistive, and Hall MHD limits. This is accomplished by first re-writing the XMHD equations such that it allows the algorithm to retain the use of ideal MHD Riemann solvers and the constrained transport method to preserve divergence-free magnetic fields. Incorporation of electron inertia and displacement current introduces additional numerical stiffness which motivates a semi-implicit FV scheme that re-formulates the XMHD model as a relaxation system. The equations are then advanced in time using an explicit 2nd-order Runge–Kutta scheme with operator splitting applied to the implicit source term updates at each sub-stage. For additional numerical stability, a density-dependent slope limiter is implemented to increase flux diffusivity at low density regions where non-ideal effects become significant. The algorithm is subsequently implemented in a scalable adaptive mesh refinement (AMR) framework. As the new algorithm retains many aspects of the ideal MHD formulations, it asymptotes naturally to the ideal MHD limit. Moreover, it shows promising results at the resistive and Hall MHD limits. This is verified against reference test problems for ideal, resistive and Hall MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Cost-constrained adaptive simulations of transient spray combustion in a gas turbine combustor

Predictive high-fidelity simulations of turbulent spray combustion must capture the combined effects of complex chemistry, multiphase evaporating flow and spray-flame interactions to achieve physical accuracy. Finite-rate chemistry (FRC) combined with a realistic chemical mechanism is a combustion model well-suited for this purpose, but has a high computational cost due to the large number and stiffness of transported chemical species. In contrast, flamelet-based models achieve lower cost by transporting a small number of quantities of reduced stiffness, but assumptions regarding local flame topology, boundary conditions and inter-phase coupling limit their physical accuracy. Recently, the Pareto-efficient combustion (PEC) framework was developed to dynamically assign combustion models based on local cost and accuracy metrics in gas-phase reacting flows. In this work, we extend this PEC framework to spray combustion through the rigorous analysis of the multiphase coupling terms in the governing equations. The derivation shows that spray evaporation causes errors in the prediction of species mass fractions for flamelet-based models due to the sensitivity of the local thermo-chemical state to changes in composition caused by fuel vaporization across combustion regimes present in practical spray combustion devices. Sub-model assignment is formulated as a multiple-choice knapsack problem, where computational cost is directly controlled through the fraction of the domain assigned to the FRC sub-model. The extended PEC formulation is applied to the simulation of a realistic rich-quench-lean gas turbine combustor at steady-state conditions, as well as transient operation resulting in lean blow-out (LBO). Analysis of transient simulations during LBO demonstrates the extended PEC formulation’s capacity to dynamically adapt to changing conditions within the combustor. Transient combustor dynamics are shown to approach convergence with limited increases in computational cost, while retaining substantial computational cost reduction compared to monolithic FRC simulations. Through PEC simulations with increasing fractions of the domain assigned to FRC, monolithic flamelet simulations are shown to over-predict flame stability during LBO. Finally, the extended PEC formulation is thus shown to overcome deficiencies of monolithic models by controlling modeling error for multiphase combustion modeling.

33 ADVANCED PROPULSION SYSTEMS↗