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,351 records · Page 75

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.↗

Utilization of Unmanned Aircraft Systems for Environmental Purposes at the Savannah River Site – 26578

Born in the 1950s, the Savannah River Plant was constructed as a plutonium and tritium production plant. As the Cold War era came and went, a change of direction was signaled when the name of the facility changed from Savannah River Plant to Savannah River Site (1989) – the main mission at Savannah River shifted from “production” to “cleanup.” The Department of Energy is responsible for managing the 310 square miles of land that is the Savannah River Site and the cleanup/remediation activities that occur. The Savannah River National Laboratory utilizes unmanned aircraft systems to assist with some of those environmental monitoring and remediation activities. One function of unmanned aircraft systems at the Savannah River National Laboratory is conducting aerial photography and videography. Various infrastructure evaluations have been performed with unmanned aircraft – facility rooftop inspections, water tower lock verifications, earthen dam surveys, etc. An unmanned aircraft system has been used for progress footage of remediation projects – Dunbarton Bay remediation, Saltstone Disposal Unit construction, etc. Thermal cameras on an unmanned aircraft system are used to help identify cracks/leaks in structures from vantage points not readily accessible by personnel. Virtual tours of different waste units are conducted with unmanned aircraft systems for Site groups and federal/state regulators to save costs and travel time. Additionally, advanced remote sensing technologies are used on unmanned aircraft systems at the Savannah River Site. Light Detection and Ranging uses laser pulses to measure distances to the Earth's surface or other targets to create highly detailed topographic maps that accurately portray terrain (e.g., elevation changes, slopes, and contours). Data collection with a Light Detection and Ranging unit mounted on an unmanned aircraft system is quick and efficient – large-area surveys are conducted in shorter time frames. Because Light Detection and Ranging can penetrate through foliage and vegetation to ground level, it is being used in conjunction with a watershed model to study the General Separations Area of the Site to determine stormwater flow volume and direction, monitor stream stability, and forecast river flow behavior. A hyperspectral camera captures a wide spectrum of electromagnetic wavelengths across numerous narrow bands, beyond what the human eye can see. It provides detailed spectral information about the objects or surfaces in an image, which can be used to identify and analyze materials based on their spectral signatures. A hyperspectral camera mounted on an unmanned aircraft system has been used at the Savannah River Site for non-destructive evaluation of concrete and concrete structures. Unmanned aircraft systems at the Savannah River Site are also used for the application of herbicide on identified rooftops. P-Reactor and R-Reactor are entombed nuclear reactors at the Savannah River Site. Unwanted vegetation is often present on these rooftops – over time roots can penetrate seams, cracks, and joints of the general roof and concrete caps, leading to water ingress and deterioration of the structural components. For time and cost efficiencies as well as in the interest of personnel safety, an unmanned aircraft is used to dispense herbicide on rooftop areas where vegetation is noticed to help alleviate the issues/hazards.

Lorier, Troy [Savannah River National Laboratory (↗

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↗

Omega-Limit Sets and Input-To-State Stability in Power Grids with Switching Equilibria

This paper studies a power transmission system with both conventional generators (CGs) and distributed energy assets (DEAs) providing frequency control. We consider an operating condition with demand aggregating two dynamic components: one that switches between different values on a finite set, and one that varies smoothly over time. Such dynamic operating conditions may result from protection scheme activations, external cyber-attacks, or due to the integration of dynamic loads, such as data centers. Mathematically, the dynamics of the resulting system are captured by a system that switches between a finite number of vector fields - or modes -, with each mode having a distinct equilibrium point induced by the demand aggregation. To analyze the stability properties of the resulting switching system, we leverage tools from hybrid dynamic inclusions and the concept of ..omega.. -limit sets from sets. Specifically, we characterize a compact set that is semi-globally practically asymptotically stable under the assumption that the switching frequency and load variation rate are sufficiently slow. For arbitrarily fast variations of the load, we use a level-set argument with multiple Lyapunov functions to establish input-to-state stability of a larger set and with respect to the rate of change of the loads. The theoretical results are illustrated via numerical simulations on the IEEE 39-bus test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

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.↗

Quantitative proteomics reveals extensive lysine ubiquitination and transcription factor stability states in Arabidopsis

Abstract Protein activity, abundance, and stability can be regulated by post-translational modification including ubiquitination. Ubiquitination is conserved among eukaryotes and plays a central role in modulating cellular function; yet, we lack comprehensive catalogs of proteins that are modified by ubiquitin in plants. In this study, we describe an antibody-based approach to enrich ubiquitinated peptides coupled with isobaric labeling to enable quantification of up to 18-multiplexed samples. This approach identified 17,940 ubiquitinated lysine sites arising from 6,453 proteins from Arabidopsis (Arabidopsis thaliana) primary roots, seedlings, and rosette leaves. Gene ontology analysis indicated that ubiquitinated proteins are associated with numerous biological processes including hormone signaling, plant defense, protein homeostasis, and metabolism. We determined ubiquitinated lysine residues that directly regulate the stability of three transcription factors, CRYPTOCHROME-INTERACTING BASIC-HELIX-LOOP-HELIX 1 (CIB1), CIB1 LIKE PROTEIN 2 (CIL2), and SENSITIVE TO PROTON RHIZOTOXICITY1 (STOP1) using in vivo degradation assays. Furthermore, codon mutation of CIB1 to create a K166R conversion to prevent ubiquitination, via CRISPR/Cas9-derived adenosine base editing, led to an early flowering phenotype and increased expression of FLOWERING LOCUS T (FT). These comprehensive site-level ubiquitinome profiles provide a wealth of data for future functional studies related to modulation of biological processes mediated by this post-translational modification in plants.

Biochemistry & Molecular Biology↗