Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical stability”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,009 records · Page 56

A theoretical framework for convergence and continuous dependence of estimates in inverse problems for distributed parameter systems

Numerical techniques for parameter identification in distributed-parameter systems are developed analytically. A general convergence and stability framework (for continuous dependence on observations) is derived for first-order systems on the basis of (1) a weak formulation in terms of sesquilinear forms and (2) the resolvent convergence form of the Trotter-Kato approximation. The extension of this framework to second-order systems is considered.

Banks, H. T.↗

Application of precomputed control laws in a reconfigurable aircraft flight control system

A self-repairing flight control system concept in which the control law is reconfigured after actuator and/or control surface damage to preserve stability and pilot command tracking is described. A key feature of the controller is reconfigurable multivariable feedback. The feedback gains are designed off-line and scheduled as a function of the aircraft control impairment status so that reconfiguration is performed simply by updating the gain schedule after detection of an impairment. A novel aspect of the gain schedule design procedure is that the schedule is calculated using a linear quadratic optimization-based simultaneous stabilization algorithm in which the scheduled gain is constrained to stabilize a collection of plant models representing the aircraft in various control failure modes. A description and numerical evaluation of a controller design for a model of a statically unstable high-performance aircraft are given.

Moerder, Daniel D.↗

Scattering from coated structures and antenna pattern control using impedance surfaces, part A/B

The scattering from coated, conducting structures, specifically the coated dihedral corner reflector configuration and the coated strip/plate configuration is examined. The formulation uses impedance-wedge Uniform Theory of Diffraction scattering coefficients to calculate the diffracted fields. A finite-thickness coating is approximated using the impedance boundary condition to arrive at an equivalent impedance for the coating. The formulation of the impedance wedge coefficients is outlined. Far-field, perfectly conducting approximations are discussed. Problems with the present dihedral corner reflector model for certain angles of incidence and observation are discussed along with a potentially rectifying modification. Also, the capacity to measure the electromagnetic properties of lossy materials was developed. The effects of using multiple material coatings on the radiation pattern of the horn antenna were studied. Numerous computations were devoted toward the inverse problem of synthesizing desired radiation patterns using the impedance surfaces. Stabilizing the equivalent sheet impedance using the linear control condition was attempted, and it was found to be a very difficult task.

Balanis, Constantine A.↗

Comparison between pressure gradient method and MAC method on high Re calculation

A cavity flow driven by shear and buoyancy forces is used as a test problem in the application of a nonstaggered pressure gradient (PG) method in solving the two-dimensional incompressible Navier-Stokes equations. Twelve finite differencing schemes are used to solve the cavity flow problem. The schemes consist of various combinations of grid arrangements, upwinding treatments, and conservativeness of convection terms. An artificial source term is introduced, and the solutions are compared with those obtained by the conventional marker-and-cell (MAC) method. The comparisons favor the PG method. Numerical results obtained by the twelve schemes are compared with exact solutions in order to assess the stability and accuracy of each scheme.

Tan, C.-H.↗

On the nonlinear stability of a high-speed, axisymmetric boundary layer

The stability of a high-speed, axisymmetric boundary layer is investigated using secondary instability theory and direct numerical simulation. Parametric studies based on the temporal secondary instability theory identify subharmonic secondary instability as a likely path to transition on a cylinder at Mach 4.5. The theoretical predictions are validated by direct numerical simulation at temporally-evolving primary and secondary disturbances in an axisymmetric boundary-layer flow. At small amplitudes of the secondary disturbance, predicted growth rates agree to several significant digits with values obtained from the spectrally-accurate solution of the compressible Navier-Stokes equations. Qualitative agreement persists to large amplitudes of the secondary disturbance. Moderate transverse curvature is shown to significantly affect the growth rate of axisymmetric second mode disturbances, the likely candidates of primary instability. The influence of curvature on secondary instability is largely indirect but most probably significant, through modulation of the primary disturbance amplitude. Subharmonic secondary instability is shown to be predominantly inviscid in nature, and to account for spikes in the Reynolds stress components at or near the critical layer.

Pruett, C. David↗

Buoyancy-driven instabilities and the nonlinear breakup of a sheared magnetic layer

Motivated by problems concerning the storage and subsequent escape of the solar magnetic field, a study has been made of how a magnetic layer embedded in a convectively stable atmosphere evolves due to axisymmetric instabilities driven by magnetic buoyancy. The initial equilibrium consists of a toroidal field sheared by a weaker poloidal component. The linear stability problem is investigated for both ideal and resistive MHD, and the nonlinear evolution is followed by numerical integration of the equations of motion. In all cases, the instability is greatly affected by the distribution and strength of the poloidal field. In particular, both the horizontal and vertical scales of the motions are controlled by the location of the surface on which the poloidal field vanishes: the resonant surface. In the nonlinear regime, a resonant surface close to the interface between the magnetized and field-free fluid leads to the localization of the instability, so that only a fraction of the magnetic region is disrupted by the motions. By contrast, a deeply seated resonant surface leads to the complete disruption of the layer and to the formation of large, helical magnetic fragments whose identity is preserved for the entire simulation.

Cattaneo, Fausto↗

Computational fluid dynamic analysis of liquid rocket combustion instability

The paper presents a computational analysis of liquid rocket combustion instability. Consideration is given to both a fully nonlinear unsteady calculation as well as a new CFD-based linearized stability analysis. An analytical solution for the linear stability problem in a constant area combustion chamber with uniform mean flow is developed to verify the numerical analyses.

Venkateswaran, Sankaran↗

Variation of crack-opening stresses in three-dimensions - Finite thickness plate

A 3D elastic-plastic finite-element analysis is conducted to study crack-growth behavior of thin and thick center-cracked specimens under constant-amplitude loading conditions. The numerical analysis and the specimen configuration and loading are described for both the thin and thick conditions. Stabilized crack-opening stresses of interior and exterior regions are given as are the closure and opening profiles of the crack-surface plane after the tenth cycle. The effect of thickness is discussed with respect to the crack-opening stress levels and the plastic zones of the interior and exterior regions. A load-reduced-displacement technique allows the calculation of the crack-opening stresses at three locations on the crack surface plane. The constraint effect related to thickness gives a lower stabilized crack-opening stress level for the thick specimens.

Chermahini, R. G.↗

On the nonlinear stability of a high-speed, axisymmetric boundary layer

The stability of a high-speed, axisymmetric boundary layer is investigated using secondary instability theory and direct numerical simulation. Parametric studies based on the temporal secondary instability theory identify subharmonic secondary instability as a likely path to transition on a cylinder at Mach 4.5. The theoretical predictions are validated by direct numerical simulation at temporally-evolving primary and secondary disturbances in an axisymmetric boundary-layer flow. At small amplitudes of the secondary disturbance, predicted growth rates agree to several significant digits with values obtained from the spectrally-accurate solution of the compressible Navier-Stokes equations. Qualitative agreement persists to large amplitudes of the secondary disturbance. Moderate transverse curvature is shown to significantly affect the growth rate of axisymmetric second mode disturbances, the likely candidates of primary instability. The influence of curvature on secondary instability is largely indirect but most probably significant, through modulation of the primary disturbance amplitude. Subharamonic secondary instability is shown to be predominantly inviscid in nature, and to account for spikes in the Reynolds stress components at or near the critical layer.

Pruett, C. D.↗

CO2 laser preionisation

The final report for work done during the reporting period of January 25, 1990 to January 24, 1991 is presented. A literature survey was conducted to identify the required parameters for effective preionization in TEA CO2 lasers and the methods and techniques for characterizing preionizers are reviewed. A numerical model of the LP-140 cavity was used to determine the cause of the transverse mode stability improvement obtained when the cavity was lengthened. The measurement of the voltage and current discharge pulses on the LP-140 were obtained and their subsequent analysis resulted in an explanation for the low efficiency of the laser. An assortment of items relating to the development of high-voltage power supplies is also provided. A program for analyzing the frequency chirp data files obtained with the HP time and frequency analyzer is included. A program to calculate the theoretical LIMP chirp is also included and a comparison between experiment and theory is made. A program for calculating the CO2 linewidth and its dependence on gas composition and pressure is presented. The program also calculates the number of axial modes under the FWHM of the line for a given resonator length. A graphical plot of the results is plotted.

Spiers, Gary D.↗

A finite element study of incompressible flows past oscillating cylinders and aerofoils

We present our numerical results for certain unsteady flows past oscillating cylinders and aerofoils. The computations are based on the stabilized space-time finite element formulation. The implicit equation systems resulting from the space-time finite element discretizations are solved using iterative solution techniques. One of the problems studied is flow past a cylinder which is forced to oscillate in the horizontal direction. In this case we observe a change from an unsymmetric mode of vortex shedding to a symmetric one. An extensive study was carried out for the case in which a cylinder is mounted on lightly damped springs and allowed to oscillate in the vertical direction. In this case the motion of the cylinder needs to be determined as part of the solution, and under certain conditions this motion changes the vortex-shedding pattern of the flow field significantly. This nonlinear fluid-structure interaction exhibits certain interesting behavior such as 'lock-in' and 'hysteresis', which are in good agreement with the laboratory experiments carried out by other researchers in the past. Preliminary results for flow past a pitching aerofoil are also presented.

Mittal, S.↗

Comparison of Artificial Compressibility Methods

Various artificial compressibility methods for calculating three-dimensional, steady and unsteady, laminar and turbulent, incompressible Navier-Stokes equations are compared in this work. Each method is described in detail along with appropriate physical and numerical boundary conditions. Analysis of well-posedness and numerical solutions to test problems for each method are provided. A comparison based on convergence behavior, accuracy, stability and robustness is used to establish the relative positive and negative characteristics of each method.

Kiris, Cetin↗

Numerical Mean Element Orbital Analysis with Morbiter

The Morbiter software numerically averages an osculating orbit s equations of motion (EOM) to arrive at the mean orbit s EOMs, which are then numerically propagated to obtain the long-term orbital ephemerides. The long-term evolution characteristics, and stability, of an orbit are best characterized using a mean element propagation of the perturbed, two-body variational equations of motion. The average process eliminates short period terms, leaving only secular and long period effects. Doing this avoids the Fourier series expansions and truncations required by the traditional analytic methods.

Ely, Todd A.↗

Lattice Boltzmann Method for 3-D Flows with Curved Boundary

In this work, we investigate two issues that are important to computational efficiency and reliability in fluid dynamics applications of the lattice, Boltzmann equation (LBE): (1) Computational stability and accuracy of different lattice Boltzmann models and (2) the treatment of the boundary conditions on curved solid boundaries and their 3-D implementations. Three athermal 3-D LBE models (D3QI5, D3Ql9, and D3Q27) are studied and compared in terms of efficiency, accuracy, and robustness. The boundary treatment recently developed by Filippova and Hanel and Met et al. in 2-D is extended to and implemented for 3-D. The convergence, stability, and computational efficiency of the 3-D LBE models with the boundary treatment for curved boundaries were tested in simulations of four 3-D flows: (1) Fully developed flows in a square duct, (2) flow in a 3-D lid-driven cavity, (3) fully developed flows in a circular pipe, and (4) a uniform flow over a sphere. We found that while the fifteen-velocity 3-D (D3Ql5) model is more prone to numerical instability and the D3Q27 is more computationally intensive, the 63Q19 model provides a balance between computational reliability and efficiency. Through numerical simulations, we demonstrated that the boundary treatment for 3-D arbitrary curved geometry has second-order accuracy and possesses satisfactory stability characteristics.

Mei, Renwei↗

Robust control design with real parameter uncertainty using absolute stability theory

The purpose of this thesis is to investigate an extension of mu theory for robust control design by considering systems with linear and nonlinear real parameter uncertainties. In the process, explicit connections are made between mixed mu and absolute stability theory. In particular, it is shown that the upper bounds for mixed mu are a generalization of results from absolute stability theory. Both state space and frequency domain criteria are developed for several nonlinearities and stability multipliers using the wealth of literature on absolute stability theory and the concepts of supply rates and storage functions. The state space conditions are expressed in terms of Riccati equations and parameter-dependent Lyapunov functions. For controller synthesis, these stability conditions are used to form an overbound of the H2 performance objective. A geometric interpretation of the equivalent frequency domain criteria in terms of off-axis circles clarifies the important role of the multiplier and shows that both the magnitude and phase of the uncertainty are considered. A numerical algorithm is developed to design robust controllers that minimize the bound on an H2 cost functional and satisfy an analysis test based on the Popov stability multiplier. The controller and multiplier coefficients are optimized simultaneously, which avoids the iteration and curve-fitting procedures required by the D-K procedure of mu synthesis. Several benchmark problems and experiments on the Middeck Active Control Experiment at M.I.T. demonstrate that these controllers achieve good robust performance and guaranteed stability bounds.

How, Jonathan P.↗

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

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

Palmer, Grant↗

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

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

Cardelino, Beatriz H.↗

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

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

Cardelino, Beatriz H.↗