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 829 records · Page 46

Robustness in the presence of mixed parametric uncertainty and unmodeled dynamics

Continuing the development of the structured singular value approach to robust control design, the authors investigate the problem of computing the structured singular value (mu) in the case of mixed real parametric and complex uncertainty. The problem is shown to be equivalent to a smooth constrained finite-dimensional optimization problem. In view of the fact that the functional to be maximized may have several local extrema, an upper bound on mu whose computation is numerically tractable is established; this leads to a sufficient condition of robust stability and performance. A historical perspective on the development of the mu theory is included.

Fan, Michael K. H.↗

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative, manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. These results are highlighted. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable. Another one of our goals is to examine the stability of our predicted flows. Because time-dependent three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulations). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, Robert W.↗

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program was to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, with the funding from this NASA program, we concentrated on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we developed unique global models that allowed us to study the coupling between the different regions. These results are highlighted in the next section. Another important aspect of our NASA Theory Program concerned the effect that localized 'structure' had on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkland current patterns) or time variations in these input due to storms and substorms. Also, some of the plasma flows that we predicted with our macroscopic models could be unstable, and another one of our goals was to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulations). Therefore, another goal of our NASA Theory Program was to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This could involve a detailed comparison of kinetic, semi-kinetic, and hydrodynamic predictions for a given polar wind scenario or it could involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations provides insight into when the various models can be used with confidence.

Schunk, Robert W.↗

A multi-layer discrete-ordinate method for vector radiative transfer in a vertically-inhomogeneous, emitting and scattering atmosphere. I - Theory. II - Application

A theory is developed for discretizing the vector integro-differential radiative transfer equation including both solar and thermal radiation. A complete solution and boundary equations are obtained using the discrete-ordinate method. An efficient numerical procedure is presented for calculating the phase matrix and achieving computational stability. With natural light used as a beam source, the Stokes parameters from the model proposed here are compared with the analytical solutions of Chandrasekhar (1960) for a Rayleigh scattering atmosphere. The model is then applied to microwave frequencies with a thermal source, and the brightness temperatures are compared with those from Stamnes'(1988) radiative transfer model.

Weng, Fuzhong↗

Liquid rocket spray combustion stability analysis

A computational approach to the analysis of spray combustion stability in liquid rocket combustors is proposed which is based on the unsteady quasi-two-dimensional Euler equations with interphase source terms derived from a Lagrangian treatment of the combusting spray. Based on a preliminary evaluation, the computational methodology presented here is a promising research tool and a potential design/development aid for investigating the stability characteristics of liquid rocket engines. The method is characterized by low numerical noise; the Lagrangian treatment of the spray offers improved flexibility for the direct modeling of spray combustion.

Litchford, Ron J.↗

The flow of plasma in the solar terrestrial environment

The overall goal of our NASA Theory Program is to study the coupling, time delays, and feedback mechanisms between the various regions of the solar-terrestrial system in a self-consistent, quantitative manner. To accomplish this goal, it will eventually be necessary to have time-dependent macroscopic models of the different regions of the solar-terrestrial system and we are continually working toward this goal. However, our immediate emphasis is on the near-earth plasma environment, including the ionosphere, the plasmasphere, and the polar wind. In this area, we have developed unique global models that allow us to study the coupling between the different regions. Another important aspect of our NASA Theory Program concerns the effect that localized structure has on the macroscopic flow in the ionosphere, plasmasphere, thermosphere, and polar wind. The localized structure can be created by structured magnetospheric inputs (i.e., structured plasma convection, particle precipitation or Birkeland current patterns) or time variations in these inputs due to storms and substorms. Also, some of the plasma flows that we predict with our macroscopic models may be unstable, and another one of our goals is to examine the stability of our predicted flows. Because time-dependent, three-dimensional numerical models of the solar-terrestrial environment generally require extensive computer resources, they are usually based on relatively simple mathematical formulations (i.e., simple MHD or hydrodynamic formulation). Therefore, another long-range goal of our NASA Theory Program is to study the conditions under which various mathematical formulations can be applied to specific solar-terrestrial regions. This may involve a detailed comparison of kinetic, semikinetic, and hydrodynamic predictions for a given polar wind scenario or it may involve the comparison of a small-scale particle-in-cell (PIC) simulation of a plasma expansion event with a similar macroscopic expansion event. The different mathematical formulations have different strengths and weaknesses and a careful comparison of model predictions for similar geophysical situations will provide insight into when the various models can be used with confidence.

Schunk, R. W.↗

Average-cost based robust structural control

A method is presented for the synthesis of robust controllers for linear time invariant structural systems with parameterized uncertainty. The method involves minimizing quantities related to the quadratic cost (H2-norm) averaged over a set of systems described by real parameters such as natural frequencies and modal residues. Bounded average cost is shown to imply stability over the set of systems. Approximations for the exact average are derived and proposed as cost functionals. The properties of these approximate average cost functionals are established. The exact average and approximate average cost functionals are used to derive dynamic controllers which can provide stability robustness. The robustness properties of these controllers are demonstrated in illustrative numerical examples and tested in a simple SISO experiment on the MIT multi-point alignment testbed.

Hagood, Nesbitt W.↗

Predicting the stability of a compressible periodic parallel jet flow

It is known that mixing enhancement in compressible free shear layer flows with high convective Mach numbers is difficult. One design strategy to get around this is to use multiple nozzles. Extrapolating this design concept in a one dimensional manner, one arrives at an array of parallel rectangular nozzles where the smaller dimension is omega and the longer dimension, b, is taken to be infinite. In this paper, the feasibility of predicting the stability of this type of compressible periodic parallel jet flow is discussed. The problem is treated using Floquet-Bloch theory. Numerical solutions to this eigenvalue problem are presented. For the case presented, the interjet spacing, s, was selected so that s/omega =2.23. Typical plots of the eigenvalue and stability curves are presented. Results obtained for a range of convective Mach numbers from 3 to 5 show growth rates omega(sub i)=kc(sub i)/2 range from 0.25 to 0.29. These results indicate that coherent two-dimensional structures can occur without difficulty in multiple parallel periodic jet nozzles and that shear layer mixing should occur with this type of nozzle design.

Miles, Jeffrey H.↗

Studying Turbulence Using Numerical Simulation Databases - X Proceedings of the 2004 Summer Program

This Proceedings volume contains 32 papers that span a wide range of topics that reflect the ubiquity of turbulence. The papers have been divided into six groups: 1) Solar Simulations; 2) Magnetohydrodynamics (MHD); 3) Large Eddy Simulation (LES) and Numerical Simulations; 4) Reynolds Averaged Navier Stokes (RANS) Modeling and Simulations; 5) Stability and Acoustics; 6) Combustion and Multi-Phase Flow.

Moin, Parviz↗

NASA-JSC Protocol for the Characterization of Single Wall Carbon Nanotube Material Quality

It is well known that the raw as well as purified single wall carbon nanotube (SWCNT) material always contain certain amount of impurities of varying composition (mostly metal catalyst and non-tubular carbon). Particular purification method also creates defects and/or functional groups in the SWCNT material and therefore affects the its dispersability in solvents (important to subsequent application development). A number of analytical characterization tools have been used successfully in the past years to assess various properties of nanotube materials, but lack of standards makes it difficult to compare these measurements across the board. In this work we report the protocol developed at NASA-JSC which standardizes measurements using TEM, SEM, TGA, Raman and UV-Vis-NIR absorption techniques. Numerical measures are established for parameters such as metal content, homogeneity, thermal stability and dispersability, to allow easy comparison of SWCNT materials. We will also report on the recent progress in quantitative measurement of non-tubular carbon impurities and a possible purity standard for SWCNT materials.

Arepalli, Sivaram↗

Free Body Dynamics of a Spinning Cylinder with Planar Restraint-(a.k.a. Barrel of Fun)

The dynamic motion of a cylinder is analyzed based on rotation about its center of mass and is restrained by a plane normal to the axis passing through its center of mass at an angle. The first part of this work presented an analysis of the stability of the motion. In the current report, the governing equations are numerically integrated in time and the steady state is obtained as a limit of the transient numerical solution. The calculated data are compared with observed behaviors.

Moraru, Laurentiu↗

Color Changing Hydrogen Sensors

During the Space Shuttle Program, one of the most hazardous operation that occurred was the loading of liquid hydrogen (LH2) during fueling operations of the spacecraft. Due to hydrogen's low explosive limit, any amount leaked could lead to catastrophic event. Hydrogen's chemical properties make it ideal as a rocket fuel; however, the fuel is deemed unsafe for most commercial use because of the inability to easily detect the gas leaking. The increased use of hydrogen over traditional fossil fuels would reduce greenhouse gases and America's dependency on foreign oil. Therefore a technology that would improve safety at NASA and in the commercial sector while creating a new economic sector would have a huge impact to NASA's mission. The Chemochromic Detector for sensing hydrogen gas leakage is a color-changing detector that is useful in any application where it is important to know not only the presence but also the location of the hydrogen gas leak. This technology utilizes a chemochromicpigment and polymer matrix that can be molded or spun into rigid or pliable shapes useable in variable temperature environments including atmospheres of inert gas, hydrogen gas, or mixtures of gases. A change in color of the detector material indicates where gaseous hydrogen leaks are occurring. The irreversible sensor has a dramatic color change from beige to dark grey and remains dark grey after exposure. A reversible pigment changes from white to blue in the presence of hydrogen and reverts back to white in the presence of oxygen. Both versions of the sensor's pigments were comprised of a mixture of a metal oxide substrate and a hydro-chromic compound (i.e., the compound that changed color in the presence of hydrogen) and immediately notified the operator of the presence of low levels of hydrogen. The detector can be used in a variety of formats including paint, tape, caulking, injection molded parts, textiles and fabrics, composites, and films. This technology brings numerous benefits over the traditional hydrogen sensors: The technology has excellent temperature stability (4K to 373 K), it can be used in cryogenic fluid applications, it is easy to apply and remove; it requires no power to operate; it has a quick response time; the leak points can be detected visually or electronically; it is nonhazardous, thus environmentally friendly; it can be reversible or irreversible; it does not require on-site monitoring; has a long shelf life; the detector is very durable; and the technology is inexpensive to manufacture.

Chemochromic sensor↗

Use of an Inertia Sphere to Damp the Angular Motions of Spinning Space Vehicles

The general moment equations for a spin-stabilized vehicle with an inertia-reaction angular rate damper were considered, and it was noted that simplification would result if the damper had a spherical inertia distribution. A control system incorporating such a damper was postulated. The resulting equations were linearized, and conditions for stability were obtained from an analysis of the cubic characteristic equation. Two numerical examples were included.

Suddath, Jerrold H.↗

Dynamic stability of unidirectional fiber-reinforced viscoelastic composite plates

This paper deals with a dynamic stability analysis of unidirectional fiber-reinforced composite viscoelastic plates subjected to compressive edge loads. The integrodifferential equations governing the stability problem are obtained by using, in conjunction with a Boltzmann hereditary constitutive law for a three-dimensional viscoelastic medium, a higher-order shear deformation theory of orthotropic plates. Such a theory incorporates transverse shear deformation, transverse normal stress, and rotatory inertia effects. The solution of the stability problem as considered within this paper concerns the determination of the critical in-plane edge loads yielding the asymptotic instability. Numerical applications, based on material properties derived within the framework of Aboudi's micromechanical model, are presented and pertinent conclusions concerning the nature of the loss of stability and the influence of various parameters are outlined.

Chandiramani, N. K.↗

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

Numerical simulations of weakly magnetized and strongly magnetized relativistic jets embedded in a weakly magnetized and strongly magnetized stationary or weakly relativistic (v = c/2) sheath have been performed. A magnetic field parallel to the flow is used in these simulations performed by the new GRMHD numerical code RAISHIN used in its RMHD configuration. In the numerical simulations the Lorentz factor gamma = 2.5 jet is precessed to break the initial equilibrium configuration. In the simulations sound speeds are less than or equal to c/the square root of 3 in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The Alfven wave speed is less than or equal to 0.07 c in the weakly magnetized simulations and less than or equal to 0.56 c in the strongly magnetized simulations. The results of the numerical simulations are compared to theoretical predictions from a normal mode analysis of the linearized relativistic magnetohydrodynamic (RMHD) equations capable of describing a uniform axially magnetized cylindrical relativistic jet embedded in a uniform axially magnetized relativistically moving sheath. The theoretical dispersion relation allows investigation of effects associated with maximum possible sound speeds, Alfven wave speeds near light speed and relativistic sheath speeds. The prediction of increased stability of the weakly magnetized system resulting from c/2 sheath speeds and the stabilization of the strongly magnetized system resulting from c/2 sheath speeds is verified by the numerical simulation results.

Mizuno, Yosuke↗

Stability analysis of spectral methods for hyperbolic initial-boundary value systems

A constant coefficient hyperbolic system in one space variable, with zero initial data is discussed. Dissipative boundary conditions are imposed at the two points x = + or - 1. This problem is discretized by a spectral approximation in space. Sufficient conditions under which the spectral numerical solution is stable are demonstrated - moreover, these conditions have to be checked only for scalar equations. The stability theorems take the form of explicit bounds for the norm of the solution in terms of the boundary data. The dependence of these bounds on N, the number of points in the domain (or equivalently the degree of the polynomials involved), is investigated for a class of standard spectral methods, including Chebyshev and Legendre collocations.

Gottlieb, D.↗

Receptivity of Hypersonic Boundary Layers to Distributed Roughness and Acoustic Disturbances

Boundary-layer receptivity and stability of Mach 6 flows over smooth and rough seven-degree half-angle sharp-tipped cones are numerically investigated. The receptivity of the boundary layer to slow acoustic disturbances, fast acoustic disturbances, and vortical disturbances is considered. The effects of three-dimensional isolated roughness on the receptivity and stability are also simulated. The results for the smooth cone show that the instability waves are generated in the leading edge region and that the boundary layer is much more receptive to slow acoustic waves than to the fast acoustic waves. Vortical disturbances also generate unstable second modes, however the receptivity coefficients are smaller than that of the slow acoustic wave. Distributed roughness elements located near the nose region decreased the receptivity of the second mode generated by the slow acoustic wave by a small amount. Roughness elements distributed across the continuous spectrum increased the receptivity of the second mode generated by the slow and fast acoustic waves and the vorticity wave. The largest increase occurred for the vorticity wave. Roughness elements distributed across the synchronization point did not change the receptivity of the second modes generated by the acoustic waves. The receptivity of the second mode generated by the vorticity wave increased in this case, but the increase is lower than that occurred with the roughness elements located across the continuous spectrum. The simulations with an isolated roughness element showed that the second mode waves generated by the acoustic disturbances are not influenced by the small roughness element. Due to the interaction, a three-dimensional wave is generated. However, the amplitude is orders of magnitude smaller than the two-dimensional wave.

Balakumar, P.↗

Mixing and Transport in the Solar Nebula

Boss & Vanhala (2000, 2001) prepared reviews of triggered collapse and injection models, using Prudence Foster's finite differences code at very high spatial resolution (440 x 1440 cells) to demonstrate the convergence of the R-T fingers in triggered injection models. A two dimensional hydrodynamical calculation with unprecedentedly high spatial resolution (960 x 2880 zones, or almost 3 million grid points) demonstrated that it suitable shock front can both trigger the collapse of an otherwise stable presolar cloud, and inject shock front particles into the collapsing cloud through the formation of what become Rayleigh-Taylor fingers of compressed fluid layers falling into the gravitational potential well of the growing protostar. These calculations suggest that heterogeneity derived from these R-T fingers will persist down to the scale of their injection onto the surface of the solar nebula. Haghighipour developed a numerical code capable of calculating the orbital evolution of dust grains of varied sizes in a gaseous nebula, subject to Epstein and Stokes drag as well as the self-gravity of the disk. In collaboration with the PI and George W. Wetherill, Haghighipour has been involved in development of a new idea on the possibility of rapid formation of ice giant planets via the disk instability mechanism. Haghighipour studied the stability of a five-body system consisting of the Sun and four protoplanets by numerically integrating their equations of motions. Using Levison and Duncan s SWIFT integrator, Haghighipour showed that, depending on the orbital parameters of the bodies, such a system can be stable for 0.1-10 Myr. Time periods of 1 Myr or more are long enough to be consistent with the time scale proposed for the formation of giant planets by the disk instability mechanism and the photoevaporation of the gaseous envelopes of the outermost protoplanets by a nearby OB star, resulting in the formation of ice giant planets. The PI has used his three dimensional models of marginally gravitationally unstable disks to study the preservation of isotopic heterogeneity in evolving protoplanetary disks. Such heterogeneity might arise from the infall onto the disk s surface of solids processed in the X-wind region of the disk, or derived from stellar nucleosynthesis and injected by R-T fingers. The technique used consists of solving a color equation, identical to the gas continuity equation, which follows the time evolution in three space dimensions of an arbitrarily placed initial color field, i.e., a dye inserted the disk. The models show that significant concentrations of color could persist for time periods of about a thousand years or more, even in the most dynamically active region of such a disk. Such a time period might be long enough for solids to coagulate and grow to significant sizes while retaining the isotopic signature of their birth region in the nebula.

Boss, Alan P.↗