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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Model free simulations of a high speed reacting mixing layer

The effects of compressibility, chemical reaction exothermicity and non-equilibrium chemical modeling in a combusting plane mixing layer were investigated by means of two-dimensional model free numerical simulations. It was shown that increased compressibility generally had a stabilizing effect, resulting in reduced mixing and chemical reaction conversion rate. The appearance of 'eddy shocklets' in the flow was observed at high convective Mach numbers. Reaction exothermicity was found to enhance mixing at the initial stages of the layer's growth, but had a stabilizing effect at later times. Calculations were performed for a constant rate chemical rate kinetics model and an Arrhenius type kinetics prototype. The Arrhenius model was found to cause a greater temperature increase due to reaction than the constant kinetics model. This had the same stabilizing effect as increasing the exothermicity of the reaction. Localized flame quenching was also observed when the Zeldovich number was relatively large.

Steinberger, Craig J.↗
NASA NTRS1992

Materials processing in a centrifuge - Numerical modeling of macrogravity effects

The fluid mechanics associated with crystal growth processes on a centrifuge is investigated. A simple scaling analysis is used to examine the relative magnitudes of the forces acting on the system and good agreement is obtained with previous studies. A two-dimensional model of crystal growth on a centrifuge is proposed and calculations are undertaken to help in understanding the fundamental transport processes within the crystal growth cell. Results from three-dimensional calculations of actual centrifuge-based crystal growth systems are presented both for the thermodynamically stable and unstable configurations. The calculations show the existence of flow bifurcations in certain configurations but not in all instances. The numerical simulations also show that the centrifugal force is the dominant stabilizing force on fluid convection in the stable configuration. The stabilizing influence of the Coriolis force is found to be only secondary in nature. No significant impact of gravity gradient is found in the calculations. Simulations of unstable configurations show that the Coriolis force has a stabilizing influence on fluid motion by delaying the onset of unsteady convection. Detailed flow and thermal field characteristics are presented for all the different cases that are simulated.

Ramachandran, N.↗
NASA NTRS1955

Determination of lateral-stability derivatives and transfer-function coefficients from frequency-response data for lateral motions

A method is presented for determining the lateral-stability derivatives, transfer-function coefficients, and the modes for lateral motion from frequency-response data for a rigid aircraft. The method is based on the application of the vector technique to the equations of lateral motion, so that the three equations of lateral motion can be separated into six equations. The method of least squares is then applied to the data for each of these equations to yield the coefficients of the equations of lateral motion from which the lateral-stability derivatives and lateral transfer-function coefficients are computed. Two numerical examples are given to demonstrate the use of the method.

Donegan, James J↗
NASA NTRS1994

The extensional rheology of non-Newtonian materials

It has been proposed to measure the extensional viscosity function of a non-Newtonian polymer solution in a reduced gravity environment as part of the Advanced Fluid Module. In ground-based extensional measurements, the no-sip boundary condition at solid-fluid interfaces always result in appreciable shear gradients in the test fluid; however the removal of gravitational body forces permits controlled extensional deformation of containerless test samples and the first unambiguous measurements of this kind. Imperative to successful implementation of this experiment is the generation and subsequent deformation of a stable cylindrical column of test fluid. A study of the generation and deformation of liquid bridges demonstrates that Newtonian liquid bridges undergo capillary breakup as anticipated when stretched beyond a critical aspect ratio; non-Newtonian liquid bridges, however, are stabilized by the strain-hardening phenomenon exhibited by these materials. Numerical simulations of Newtonian breakup are compared with experimental results, and show that previous ground-based attempts at measuring the extensional viscosity of Newtonian fluids are of limited accuracy.

Spiegelberg, Stephen H.↗
NASA NTRS1995

The initial-value problem for viscous channel flows

Plane viscous channel flows are perturbed and the ensuing initial-value problems are investigated in detail. Unlike traditional methods where traveling wave normal modes are assumed for solution, this works offers a means whereby completely arbitrary initial input can be specified without having to resort to eigenfunction expansions. The full temporal behavior, including both early time transients and the long time asymptotics, can be determined for any initial disturbance. Effects of three-dimensionality can be assessed. The bases for the analysis are: (a) linearization of the governing equations; (b) Fourier decomposition in the spanwise and streamwise directions of the flow; and (c) direct numerical integration of the resulting partial differential equations. All of the stability data that are known for such flows can be reproduced. Also, the optimal initial condition can be determined in a straight forward manner and such optimal conditions clearly reflect transient growth data that is easily determined by a rational choice of a basis for the initial conditions. Although there can be significant transient growth for subcritical values of the Reynolds number using this approach it does not appear possible that arbitrary initial conditions will lead to the exceptionally large transient amplitudes that have been determined by optimization of normal modes. The approach is general and can be applied to other classes of problems where only a finite discrete spectrum exists, such as the boundary layer for example.

Criminale, W. O.↗
NASA NTRS2009

Orbiter Boundary Layer Transition Stability Modeling at Flight Entry Conditions

State of the art boundary layer stability modeling capabilities are increasingly seeing application to entry flight vehicles. With the advent of user friendly and robust implementations of two-dimensional chemical nonequilibrium stability modeling with the STABL/PSE-CHEM software, the need for flight data to calibrate such analyses capabilities becomes more critical. Recent efforts to perform entry flight testing with the Orbiter geometry related to entry aerothermodynamics and boundary layer transition is allowing for a heightened focus on the Orbiter configuration. A significant advancement in the state of the art can likely be achieved by establishing a basis of understanding for the occurrence of boundary layer transition on the Orbiter due to discrete protruding gap fillers and the nominal distributed roughness of the actual thermal protection system. Recent success in demonstrating centerline two-dimensional stability modeling on the centerline of the Orbiter at flight entry conditions provides a starting point for additional investigations. The more detailed paper will include smooth Orbiter configuration boundary layer stability results for several typical orbiter entry conditions. In addition, the numerical modeling approach for establishing the mean laminar flow will be reviewed and the method for determining boundary layer disturbance growth will be overviewed. In addition, if actual Orbiter TPS surface data obtained via digital surface scans become available, it may be possible to investigate the effects of an as-flown flight configuration on boundary layer transition compared to a smooth CAD reference.

Bartkowicz, Matt↗
NASA NTRS2023

2023 Astrophotonics Roadmap: Pathways to Realizing Multi-Functional Integrated Astrophotonic Instruments

Photonic technologies offer numerous functionalities that can be used to realize astrophotonic instruments. The most spectacular example to date is the ESO Gravity instrument at the Very Large Telescope in Chile that combines the light-gathering power of four 8 m telescopes through a complex photonic interferometer. Fully integrated astrophotonic devices stand to offer critical advantages for instrument development, including extreme miniaturization when operating at the diffraction-limit, as well as integration, superior thermal and mechanical stabilization owing to the small footprint, and high replicability offering significant cost savings. Numerous astrophotonic technologies have been developed to address shortcomings of conventional instruments to date, including for example the development of photonic lanterns to convert from multimode inputs to single mode outputs, complex aperiodic fiber Bragg gratings to filter OH emission from the atmosphere, complex beam combiners to enable long baseline interferometry with for example, ESO Gravity, and laser frequency combs for high precision spectral calibration of spectrometers. Despite these successes, the facility implementation of photonic solutions in astronomical instrumentation is currently limited because of (1) low throughputs from coupling to fibers, coupling fibers to chips, propagation and bend losses, device losses, etc, (2) difficulties with scaling to large channel count devices needed for large bandwidths and high resolutions, and (3) efficient integration of photonics with detectors, to name a few. In this roadmap, we identify 24 key areas that need further development. We outline the challenges and advances needed across those areas covering design tools, simulation capabilities, fabrication processes, the need for entirely new components, integration and hybridization and the characterization of devices. To realize these advances the astrophotonics community will have to work cooperatively with industrial partners who have more advanced manufacturing capabilities. With the advances described herein, multi-functional integrated instruments will be realized leading to novel observing capabilities for both ground and space based platforms, enabling new scientific studies and discoveries.

astrophotonics↗
NASA NTRS2015

A Numerical and Experimental Study of Coflow Laminar Diffusion Flames: Effects of Gravity and Inlet Velocity

In this work, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting behavior of laminar coflow methane-air diffusion flames was investigated both computationally and experimentally. A series of flames measured in the Structure and Liftoff in Combustion Experiment (SLICE) was assessed numerically under microgravity and normal gravity conditions with the fuel stream CH4 mole fraction ranging from 0.4 to 1.0. Computationally, the MC-Smooth vorticity-velocity formulation of the governing equations was employed to describe the reactive gaseous mixture; the soot evolution process was considered as a classical aerosol dynamics problem and was represented by the sectional aerosol equations. Since each flame is axisymmetric, a two-dimensional computational domain was employed, where the grid on the axisymmetric domain was a nonuniform tensor product mesh. The governing equations and boundary conditions were discretized on the mesh by a nine-point finite difference stencil, with the convective terms approximated by a monotonic upwind scheme and all other derivatives approximated by centered differences. The resulting set of fully coupled, strongly nonlinear equations was solved simultaneously using a damped, modified Newton's method and a nested Bi-CGSTAB linear algebra solver. Experimentally, the flame shape, size, lift-off height, and soot temperature were determined by flame emission images recorded by a digital camera, and the soot volume fraction was quantified through an absolute light calibration using a thermocouple. For a broad spectrum of flames in microgravity and normal gravity, the computed and measured flame quantities (e.g., temperature profile, flame shape, lift-off height, and soot volume fraction) were first compared to assess the accuracy of the numerical model. After its validity was established, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting tendency of laminar coflow methane-air diffusion flames was explored further by examining quantities derived from the computational results.

microgravity↗
NASA NTRS1990

Staggered solution procedures for multibody dynamics simulation

The numerical solution procedure for multibody dynamics (MBD) systems is termed a staggered MBD solution procedure that solves the generalized coordinates in a separate module from that for the constraint force. This requires a reformulation of the constraint conditions so that the constraint forces can also be integrated in time. A major advantage of such a partitioned solution procedure is that additional analysis capabilities such as active controller and design optimization modules can be easily interfaced without embedding them into a monolithic program. After introducing the basic equations of motion for MBD system in the second section, Section 3 briefly reviews some constraint handling techniques and introduces the staggered stabilized technique for the solution of the constraint forces as independent variables. The numerical direct time integration of the equations of motion is described in Section 4. As accurate damping treatment is important for the dynamics of space structures, we have employed the central difference method and the mid-point form of the trapezoidal rule since they engender no numerical damping. This is in contrast to the current practice in dynamic simulations of ground vehicles by employing a set of backward difference formulas. First, the equations of motion are partitioned according to the translational and the rotational coordinates. This sets the stage for an efficient treatment of the rotational motions via the singularity-free Euler parameters. The resulting partitioned equations of motion are then integrated via a two-stage explicit stabilized algorithm for updating both the translational coordinates and angular velocities. Once the angular velocities are obtained, the angular orientations are updated via the mid-point implicit formula employing the Euler parameters. When the two algorithms, namely, the two-stage explicit algorithm for the generalized coordinates and the implicit staggered procedure for the constraint Lagrange multipliers, are brought together in a staggered manner, they constitute a staggered explicit-implicit procedure which is summarized in Section 5. Section 6 presents some example problems and discussions concerning several salient features of the staggered MBD solution procedure are offered in Section 7.

Park, K. C.↗
NASA NTRS2008

3D RMHD 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 flow using the RAISHIN code. In the numerical simulations a jet with Lorentz factor gamma=2.5 is processed 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 sheath flow. 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 NTRS1984

Orbital stability of the unseen solar companion linked to periodic extinction events

Evidence from three-dimensional numerical modelling is presented that only cometary orbits with a limited range in inclination with respect to the galactic plane are formally stable for the length of time required to cause periodic extinction events. The calculations were done using Cowell's method employing a fourth-order Runge-Kutta integration scheme in an inertial reference frame in orbit about the Galaxy. Tidal perturbations in the radial direction due to the Galaxy and the Coriolis forces are included. The vertical component of the gravitational field of the galactic disk is superimposed on these forces. The results indicate that orbits for Nemesis that are inclined at more than 30 deg to the galactic plane are not allowed and suggests that the search for Nemesis should be concentrated toward the plane of the Galaxy. Perturbations by passing stars or molecular clouds may make even the low-inclination orbits unstable.

Torbett, M. V.↗
NASA NTRS1990

Ballooning mode instability at the plasmapause

The ballooning mode instability, which can excite hydromagnetic waves at the plasmapause, is studied in the presence of azimuthal plasma flows induced during geomagnetically disturbed periods. A general sufficient criterion for the ballooning mode stability is derived, for a change in the potential energy greater than or equal to 0, which involves the integration over an entire field line. A local stability analysis at the equatorial plasmapause region shows that the ballooning modes could be spontaneously generated via instability under at least two conditions: one is similar to the usual interchange condition, and the second to the quasi-interchange modes. Both of these local instability conditions can be derived from the general stability criterion. Finally an exact solution for the equilibrium state with flow is derived analytically, and the change in the potential energy is computed numerically. It is found that, in the cases studied, the flow does not spontaneously excite the ballooning modes; it only further stabilizes (or destabilizes) the ballooning spectrum if originally the system is stable (or unstable). The analysis would be useful for the interpretation of some of the low-frequency modes observed at the ground and near the equatorial plasmapause.

Lakhina, G. S.↗
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.↗