SEARCH · Search NASA

Results for “Stiff problems”

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 703 records · Page 39

All sources NASA NTRS DOE OSTI

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 OSTI2018

Resonance broadened quasi-linear (RBQ) model for fast ion distribution relaxation due to Alfvénic eigenmodes

The burning plasma performance is limited by the confinement of the super-alfvénic fusion products such as alpha particles and the auxiliary heating ions capable of exciting the Alfvénic eigenmodes (AEs) (Gorelenkov et al 2014 Nucl. Fusion 54 125001). Here in this work the effect of AEs on fast ions is formulated within the quasi-linear (QL) theory generalized for this problem recently (Duarte 2017 PhD Thesis University of São Paulo, Brazil). The generalization involves the resonance line broadened interaction of energetic particles (EP) with AEs supplemented by the diffusion coefficients depending on EP position in the velocity space. A new resonance broadened QL code (or RBQ1D) based on this formulation allowing for EP diffusion in radial direction is built and presented in details. In RBQ1D applications we reduce the wave particle interaction (WPI) dynamics to 1D case when the particle kinetic energy is nearly constant. The diffusion equation for EP distribution evolution is then solved simultaneously for all particles along the angular momentum direction. We make initial applications of the RBQ1D to a DIII-D plasma with elevated q-profile where the beam ions show stiff transport properties (Collins et al (The DIII-D Team) 2016 Phys. Rev. Lett. 116 095001). In conclusion, AE driven fast ion profile relaxation is studied for validations of the QL approach in realistic conditions of beam ion driven instabilities in DIII-D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

Influencing Factors in Atomic Force Microscopy Based Mechanical Characterization of Biological Cells

AFM-based single cell force spectroscopy is increasingly being used to understand many biological health problems like malaria, cancer etc. A reliable diagnostics needs, accurate measurement and a clear understanding of the influencing factors, which may otherwise sour the measured biomechanics. Results from any successful experimentation should be repeatable and error free. For this, a deeper understanding of the sources of uncertainties, which may affect the results is necessary. In order to assure the accuracy of evaluated properties and to avoid misinterpretation of the experimental data, we have categorized the common causes of uncertainties in an AFM force spectroscopy, based on their sources of origin and discussed possible remedies to them. The present work discusses the assumptions involved in AFM-based biomechanical studies (assumptions in contact model, data analysis, instrument calibration etc.) and their implications in overall estimations of mechanical biomarkers like stiffness, and adhesive strength. Advantages and disadvantages of simultaneous measurement of stiffness and adhesiveness from single force-indentation data have also been discussed.

60 APPLIED LIFE SCIENCES↗

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↗

DOE OSTI2016

Applications of Transient Grating Spectroscopy to Radiation Materials Science

Radiation damage is an inherently multi-scale phenomenon, with much of the key material evolution occurring at the mesoscale. Current techniques for the direct interrogation of material property changes due to irradiation either function with very small volumes, such as in-situ TEM, cannot cover the full range of defect populations produced, such as in positron annihilation spectroscopy, or requires destruction of the test material, such as in standard or micromechanical testing. In addition, the ultimate quantities of interest to radiation materials scientists are the changes in key material properties with irradiation. These include stiffness, ductility, and thermal conductivity for structural materials, and optical properties such as index of refraction or changes in transmittance spectra for optically active materials. The latter of these are often correlated directly to stress- or thermally-induced changes, themselves a function of the specific defect populations responsible for their change. The rapid study of such a mesoscale phenomenon requires an inherently mesoscale technique, where multiple quantities of interest can be accurately measured in very short times. Determining the sensitivity of any technique with known certainly is of critical importance to qualifying a new technique to measure in-situ material property changes during irradiation. In this paper, we present the new applicability of an existing technique, known as transient grating (TG) spectroscopy, to rapidly and accurately study key problems in radiation materials science. Results quantifying the absolute sensitivity of the technique are presented, alongside molecular dynamics (MD) simulations of the TG spectroscopy process to gauge direct agreement between theory, simulation, and experiment. An example problem of interest is demonstrated as an illustration of what is possible, showing elastic and thermal changes in irradiated silicon carbide (SiC) films. Finally, future work to quantify the effects of experimentally-simulated radiation defects is proposed, in support of building a database to use TG spectroscopy to deconvolve irradiated material property changes in-situ in terms of the defect populations responsible for their evolution. The application of the TG spectroscopy technique to in-situ radiation materials science has the potential to answer long-standing questions in the field, relating to ion beam rastering, the onset of void swelling, and direct measurement of radiation-created point defect concentrations. Before it can do so, the technique must carefully be analyzed and tuned to increase its sensitivity, the certainty of its results, and the ability to predict radiation-induced changes by simulations must be developed. Making such measurements in-situ with the aid of predictive MD simulations has the potential to elucidate mechanisms of radiation damage accumulation, as well as more rapidly qualify radiation-resistant alloys for service in nuclear systems. Work is underway to build this database of experimentally-simulated radiation-produced defects, which will be combined and compared with actual irradiated materials for eventual in-situ deconvolution of the TG signal into the defect populations responsible for irradiation-induced microstructural change. (authors)

36 MATERIALS SCIENCE↗

DOE OSTI2016

Bandwidth increasing mechanism by introducing a curve fixture to the cantilever generator

A nonlinear wideband generator architecture by clamping the cantilever beam generator with a curve fixture is proposed. Devices with different nonlinear stiffness can be obtained by properly choosing the fixture curve according to the design requirements. Three available generator types are presented and discussed for polynomial curves. Experimental investigations show that the proposed mechanism effectively extends the operation bandwidth with good power performance. Especially, the simplicity and easy feasibility allow the mechanism to be widely applied for vibration generators in different scales and environments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1985

Model-size reduction technique for the analysis of symmetric anisotropic structures

A two-step computational procedure is presented for reducing the size of the analysis model for an anisotropic symmetric structure to that of the corresponding orthotropic structure. The key elements of the procedure are: (1) decomposition of the stiffness matrix into the sum of an orthotropic and nonorthotropic (anisotropic) parts; and (2) successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate few global approximation vectors (or modes). Then the amplitudes of these modes are computed by using the Rayleigh-Ritz technique. The global approximation vectors are selected to be the solution corresponding to zero nonorthotropic matrix and its various-order derivatives with respect to an anisotropic tracing parameter (identifying the nonorthotropic material coefficients). The size of the analysis model used in generating the global approximation vectors is identical to that of the corresponding orthotropic structure. The effectiveness of the proposed technique is demonstrated by means of numerical examples and its potential for solving other quasi-symmetric problems is discussed.

Noor, A. K.↗

NASA NTRS2006

Convergence Acceleration for Multistage Time-Stepping Schemes

The convergence of a Runge-Kutta (RK) scheme with multigrid is accelerated by preconditioning with a fully implicit operator. With the extended stability of the Runge-Kutta scheme, CFL numbers as high as 1000 could be used. The implicit preconditioner addresses the stiffness in the discrete equations associated with stretched meshes. Numerical dissipation operators (based on the Roe scheme, a matrix formulation, and the CUSP scheme) as well as the number of RK stages are considered in evaluating the RK/implicit scheme. Both the numerical and computational efficiency of the scheme with the different dissipation operators are discussed. The RK/implicit scheme is used to solve the two-dimensional (2-D) and three-dimensional (3-D) compressible, Reynolds-averaged Navier-Stokes equations. In two dimensions, turbulent flows over an airfoil at subsonic and transonic conditions are computed. The effects of mesh cell aspect ratio on convergence are investigated for Reynolds numbers between 5.7 x 10(exp 6) and 100.0 x 10(exp 6). Results are also obtained for a transonic wing flow. For both 2-D and 3-D problems, the computational time of a well-tuned standard RK scheme is reduced at least a factor of four.

Swanson, R. C.↗

NASA NTRS1988

Numerical modeling of multidimensional flow in seals and bearings used in rotating machinery

The rotordynamic behavior of turbomachinery is critically dependent on fluid dynamic rotor forces developed by various types of seals and bearings. The occurrence of self-excited vibrations often depends on the rotor speed and load. Misalignment and rotor wobbling motion associated with differential clearance were often attributed to stability problems. In general, the rotative character of the flowfield is a complex three dimensional system with secondary flow patterns that significantly alter the average fluid circumferential velocity. A multidimensional, nonorthogonal, body-fitted-grid fluid flow model is presented that describes the fluid dynamic forces and the secondary flow pattern development in seals and bearings. Several numerical experiments were carried out to demonstrate the characteristics of this complex flowfield. Analyses were performed by solving a conservation form of the three dimensional Navier-Stokes equations transformed to those for a rotating observer and using the general-purpose computer code PHOENICS with the assumptions that the rotor orbit is circular and that static eccentricity is zero. These assumptions have enabled a precise steady-state analysis to be used. Fluid injection from ports near the seal or bearing center increased fluid-film direct dynamic stiffness and, in some cases, significantly increased quadrature dynamic stiffness. Injection angle and velocity could be used for active rotordynamic control; for example, injection, when compared with no injection, increased direct dynamic stiffness, which is an important factor for hydrostatic bearings.

Hendricks, R. C.↗

DOE OSTI2015

The materials used in bone tissue engineering

Bone tissue engineering looking for an alternative solution to the problem of skeletal injuries. The method is based on the creation of tissue engineered bone tissue equivalent with stem cells, osteogenic factors, and scaffolds - the carriers of these cells. For production of tissue engineered bone equivalent is advisable to create scaffolds similar in composition to natural extracellular matrix of the bone. This will provide optimal conditions for the cells, and produce favorable physico-mechanical properties of the final construction. This review article gives an analysis of the most promising materials for the manufacture of cell scaffolds. Biodegradable synthetic polymers are the basis for the scaffold, but it alone cannot provide adequate physical and mechanical properties of the construction, and favorable conditions for the cells. Addition of natural polymers improves the strength characteristics and bioactivity of constructions. Of the inorganic compounds, to create cell scaffolds the most widely used calcium phosphates, which give the structure adequate stiffness and significantly increase its osteoinductive capacity. Signaling molecules do not affect the physico-mechanical properties of the scaffold, but beneficial effect is on the processes of adhesion, proliferation and differentiation of cells. Biodegradation of the materials will help to fulfill the main task of bone tissue engineering - the ability to replace synthetic construct by natural tissues that will restore the original anatomical integrity of the bone.

60 APPLIED LIFE SCIENCES↗