Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear stability analysis”

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,099 records · Page 61

Transformation of two and three-dimensional regions by elliptic systems

Finite difference methods for composite grids were analyzed. It was observed that linear interpolation between grids would suffice only where low order accuracy was required. In the context of fluid flow, this would be in regions where the flow was essentially free stream. Higher order interpolation schemes were also investigated. The well known quadratic and cubic interpolating polynomials would increase the formal accuracy of the overall numerical algorithm. However, it can also be shown that the stability of the algorithm may be adversely affected. Further numerical results are needed in order to assess the nature of this instability induced by the interpolation procedure. Finally, error analysis and the order of difference expressions on general curvilinear coordinates are discussed.

Thompson, J. F.↗

Study of Microphonic Effects on the C100 Cryomodule for High Energy Electron Beam Accelerators

The Continuous Electron Beam Accelerator Facility (CEBAF) at Thomas Jefferson National Laboratory (JLab) is a particle accelerator which can accelerate an electron beam to relativistic speeds and apply the beam onto target samples. The C100 superconducting radio frequency (SRF) cavity is the primary accelerating structure of the C100 cryomodule, one of the many cryomodules which compose the CEBAF linear accelerator. SRF cavities are particularly sensitive to internal and external vibrations that can result in a phenomenon called microphonics which degrade the operational stability of a cryomodule. The purpose of this thesis is to investigate the significance of mechanical disturbances on the electromagnetic resonant frequency of a C100 SRF cavity. Knowledge of the mechanical resonance of the cavities and cryomodule sheds light into how these disturbances are most easily realized as deformation which causes radio frequency (RF) detuning. Three studies were conducted: the development and hammer test calibration of a Finite Element Analysis (FEA) model of a C100 cavity, the development and hammer test calibration of an FEA model of a C100 cavity string, and the hammer test of the C100-10R cryomodule at the Cryomodule Test Facility (CMTF). The cavity FEA model was found to accurately predict two modes found in two real cavities in a simply supported configuration. The cavity string FEA model leveraged the calibrated cavity FEA model but was not found to accurately predict the modal behavior of a real cavity string. Even so, the modal behavior of the cavity string inside the C100-10R cryomodule was captured during a hammer test while it was partially assembled. Finally, the C100-10R cryomodule was placed in the CMTF to study RF detuning. The RF detuning spectra during hammer hits and background noise was captured. The results of the hammer testing indicate two strong peaks at low frequencies (9-10 Hz and 22-23 Hz). These two frequencies were found to be nearly coincident to four instances of mechanical resonance found during the hammer testing done on the partially-assembled C100-10R. Because of this, these two modes are believed to contribute to RF detuning of the cryomodule. This test event also included the testing of the effectiveness of a configuration of BNNT canisters designed to act as dampers. While these tests show promising results, the lurking variables render these tests somewhat inconclusive.

Hull, Caleb James↗

Theoretical analysis of the saturation phase of the 1/1 energetic-ion-driven resistive interchange mode

The aim of the present study is to analyze the saturation regime of the energetic-ion-driven resistive interchange mode (EIC) in the LHD plasma. In this work, a set of nonlinear simulations are performed by the FAR3d code that uses a reduced MHD model for the thermal plasma coupled with a gyrofluid model for the energetic particle (EP) species. The hellically trapped EP component is introduced through a modification of the averaged drift velocity operator to include their precessional drift. The nonlinear simulation results show similar 1/1 EIC saturation phases with respect to the experimental observations, reproducing the enhancement of the n/m = 1/1 resistive interchange modes (RIC) amplitude and width as the EP β increases, the EP β threshold for the 1/1 EIC excitation, the further destabilization of the 1/1 EIC as the population of the helically trapped EP increases and the triggering of burst events. The frequency of the 1/1 EIC calculated during the burst event is 9.4 kHz and the 2/2 and 3/3 overtones are destabilized, consistent with the frequency range and the complex mode structure measured in the experiment. In addition, the simulation shows the inward propagation of the 1/1 EIC due to the nonlinear destabilization of the 3/4 and 2/3 energetic particle modes, leading to the partial overlapping between resonances during the burst event. Finally, the analysis of the 1/1 EIC stabilization phase shows the excitation of the 1/1 RIC as soon as the flattening induced by the 1/1 EIC in the pressure profile vanishes, leading to the retrieval of the pressure gradient at the plasma periphery and the overcoming of the RIC stability limit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interim Evaluation of State-of-the-Art Controllers for Low-Power Stirling Generators for Space Applications

Stirling convertors use a linear alternator to extract power from the piston motion; this requires a coil and a magnet (Faraday’s law). The coil introduces unwanted phase shift to degrade power output as well as introduce instability; i.e., low power factor. Tuning capacitor is added to increase power output and stability; widely accepted that tuning capacitor value will affect both the power output and stability.

Free-Piston Stirling Convertor Controller↗

Equilibria, stability and bifurcations of rotating columns of fluid subjected to planar disturbances

Long gyrostatically rotating drops bonded by surface tension are amenable to conventional bifurcation analysis and newer, computer-aided analytical methods, and therefore are useful prototypes of three-dimensional drops. A study is conducted by setting aside instability to Rayleigh's axisymmetric mode and investigating the effects of translationally symmetric (planar) disturbances. The disadvantage of employing single-coordinate representation of drop shapes close to break-up is brought out. It is shown that a family of symmetric two-lobed shapes bifurcates from the main family of perfectly cylindrical shapes when the rotation rate reaches a critical value, in accord with the linearized hydrodynamic analysis of Hocking.

Benner, R. E., Jr.↗

Leapfrog Time-Stepping for Hermite Methods

We introduce here Hermite-leapfrog methods for first order linear wave systems. The new Hermite-leapfrog methods pair leapfrog time-stepping with the Hermite methods of Goodrich and co-authors et al. (Math Comput 75(254):595–630, 2006). The new schemes stagger field variables in both time and space and are high-order accurate for equations with smooth solutions and coefficients. We provide a detailed description of the method and demonstrate that the method conserves variable quantities. Higher dimensional versions of the method are constructed via tensor products. Numerical evidence and rigorous analysis in one space dimension establish stability and high-order convergence. Experiments demonstrating efficient implementations on a graphics processing unit are also presented.

97 MATHEMATICS AND COMPUTING↗

Substructure analysis using NICE/SPAR and applications of force to linear and nonlinear structures

Parallel computing studies are presented for a variety of structural analysis problems. Included are the substructure planar analysis of rectangular panels with and without a hole, the static analysis of space mast, using NICE/SPAR and FORCE, and substructure analysis of plane rigid-jointed frames using FORCE. The computations are carried out on the Flex/32 MultiComputer using one to eighteen processors. The NICE/SPAR runstream samples are documented for the panel problem. For the substructure analysis of plane frames, a computer program is developed to demonstrate the effectiveness of a substructuring technique when FORCE is enforced. Ongoing research activities for an elasto-plastic stability analysis problem using FORCE, and stability analysis of the focus problem using NICE/SPAR are briefly summarized. Speedup curves for the panel, the mast, and the frame problems provide a basic understanding of the effectiveness of parallel computing procedures utilized or developed, within the domain of the parameters considered. Although the speedup curves obtained exhibit various levels of computational efficiency, they clearly demonstrate the excellent promise which parallel computing holds for the structural analysis problem. Source code is given for the elasto-plastic stability problem and the FORCE program.

Razzaq, Zia↗

A Multiplexed Quantitative Analysis of Germline Single Amino Acid Variants by Targeted Proteomics in Nondepleted Human Plasma

Single amino acid variants (SAAVs) in protein sequences are often a direct result of single-nucleotide polymorphisms (SNPs). Certain germline SAAVs have shown biological relevance in different disease conditions but lack precise quantification in circulation, which could hinder functional investigations and progress in biomarker development. Here, we have developed a multiplexed liquid chromatography-selected reaction monitoring (LC-SRM) assay that monitors 5 wild-type and variant peptide pairs (Complement Factor B: CFB-R32Q/R32W, Clusterin: CLU-N317H, Fetuin B: FETUB-K360R, and Kininogen: KNG1-L212P) in nondepleted human plasma. The assay was optimized for imprecision, linearity, stability, and calibration assessments with CVs of under 20%. The wild-type and variant peptide pairs were characterized in a set of healthy individual plasma samples. These target identifications were also validated by SNP genotyping with more than 99% accuracy. For all protein targets, we observed significantly lower concentrations of WT species in the presence variant peptides. In CFB, the concentration of R32Q was significantly lower than its counterpart R32W variant and WT species. Furthermore, our results distinguished phenotypes of homozygosity and heterozygosity of the SAAV presence through direct concentration level characterization. These findings provide some insights into how SAAVs affect quantitative assessments of target peptides. The assay demonstrates a platform for proteogenomic analyses with potential applications in both research and clinical settings.

genetics↗

Applications of singular value analysis and partial-step algorithm for nonlinear orbit determination

An adaptive method in which cruise and nonlinear orbit determination problems can be solved using a single program is presented. It involves singular value decomposition augmented with an extended partial step algorithm. The extended partial step algorithm constrains the size of the correction to the spacecraft state and other solve-for parameters. The correction is controlled by an a priori covariance and a user-supplied bounds parameter. The extended partial step method is an extension of the update portion of the singular value decomposition algorithm. It thus preserves the numerical stability of the singular value decomposition method, while extending the region over which it converges. In linear cases, this method reduces to the singular value decomposition algorithm with the full rank solution. Two examples are presented to illustrate the method's utility.

Ryne, Mark S.↗

Nonlinear system controller design based on domain of attaction: An application to CELSS analysis and control

Nonlinear system controller design based on the domain of attraction is presented. This is particularly suited to investigating Closed Ecological Life Support Systems (CELSS) models. In particular, the dynamic consequences of changes in the waste storage capacity and system mass, and how information is used for control in CELSS models are examined. The models' high dimensionality and nonlinear state equations make them difficult to analyze by any other technique. The domain of attraction is the region in initial conditions that tend toward an attractor and it is delineated by randomly selecting initial conditions from the region of state space being investigated. Error analysis is done by repeating the domain simulations with independent samples. A refinement of this region is the domain of performance which is the region of initial conditions meeting a performance criteria. In nonlinear systems, local stability does not insure stability over a larger region. The domain of attraction marks out this stability region; hence, it can be considered a measure of a nonlinear system's ability to recovery from state perturbations. Considering random perturbations, the minimum radius of the domain is a measure of the magnitude of perturbations for which recovery is guaranteed. Design of both linear and nonlinear controllers are shown. Three CELSS models, with 9 to 30 state variable, are presented. Measures of the domain of attraction are used to show the global behavior of these models under a variety of design and controller scenarios.

Babcock, P. S., IV↗

Numerical prediction of turbulent flame stability in premixed/prevaporized (HSCT) combustors

A numerical analysis of combustion instabilities that induce flashback in a lean, premixed, prevaporized dump combustor is performed. KIVA-II, a finite volume CFD code for the modeling of transient, multidimensional, chemically reactive flows, serves as the principal analytical tool. The experiment of Proctor and T'ien is used as a reference for developing the computational model. An experimentally derived combustion instability mechanism is presented on the basis of the observations of Proctor and T'ien and other investigators of instabilities in low speed (M less than 0.1) dump combustors. The analysis comprises two independent procedures that begin from a calculated stable flame: The first is a linear increase of the equivalence ratio and the second is the linear decrease of the inflow velocity. The objective is to observe changes in the aerothermochemical features of the flow field prior to flashback. It was found that only the linear increase of the equivalence ratio elicits a calculated flashback result. Though this result did not exhibit large scale coherent vortices in the turbulent shear layer coincident with a flame flickering mode as was observed experimentally, there were interesting acoustic effects which were resolved quite well in the calculation. A discussion of the k-e turbulence model used by KIVA-II is prompted by the absence of combustion instabilities in the model as the inflow velocity is linearly decreased. Finally, recommendations are made for further numerical analysis that may improve correlation with experimentally observed combustion instabilities.

Winowich, Nicholas S.↗

EXPERIMENTAL STUDIES OF NONLINEAR INTEGRABLE OPTICS

State-of-the-art accelerators at energy and intensity frontiers require increasingly bright and powerful particle beams. In conventional linear lattices, intense beams suffer from collective instabilities, resulting in beam losses and maximum beam intensity limits. This thesis presents experimental studies of a novel lattice design concept, the nonlinear integrable optics (NIO), aimed at enhancing beam stability limits with little to no beam dynamics degradation. Single-particle beam dynamics measurements of two NIO devices, the quasi-integrable octupole system and the fully integrable Danilov-Nagaitsev system, were carried out at the purpose-built Integrable Optics Test Accelerator (IOTA) at Fermilab. Their simulation, hardware design, and commissioning process are presented. Extensive model and analysis algorithm development and benchmarking is described. Electron beam data from two scientific runs is analyzed, yielding frequency and phase space dynamics consisten t with m odels. These results demonstrate viability and advantages of the NIO design, providing the groundwork for proton studies in the strong space-charge regime and future integrable accelerators.

Kuklev, Nikita↗

The AUREX cell: a versatile operando electrochemical cell for studying catalytic materials using X-ray diffraction, total scattering and X-ray absorption spectroscopy under working conditions

Understanding the structure–property relationship in electrocatalysts under working conditions is crucial for the rational design of novel and improved catalytic materials. This paper presents the Aarhus University reactor for electrochemical studies using X-rays (AUREX) operando electrocatalytic flow cell, designed as an easy-to-use versatile setup with a minimal background contribution and a uniform flow field to limit concentration polarization and handle gas formation. The cell has been employed to measure operando total scattering, diffraction and absorption spectroscopy as well as simultaneous combinations thereof on a commercial silver electrocatalyst for proof of concept. This combination of operando techniques allows for monitoring of the short-, medium- and long-range structure under working conditions, including an applied potential, liquid electrolyte and local reaction environment. The structural transformations of the Ag electrocatalyst are monitored with non-negative matrix factorization, linear combination analysis, the Pearson correlation coefficient matrix, and refinements in both real and reciprocal space. Upon application of an oxidative potential in an Ar-saturated aqueous 0.1 M KHCO 3 /K 2 CO 3 electrolyte, the face-centered cubic (f.c.c.) Ag gradually transforms first to a trigonal Ag 2 CO 3 phase, followed by the formation of a monoclinic Ag 2 CO 3 phase. A reducing potential immediately reverts the structure to the Ag (f.c.c.) phase. Following the electrochemical-reaction-induced phase transitions is of fundamental interest and necessary for understanding and improving the stability of electrocatalysts, and the operando cell proves a versatile setup for probing this. In addition, it is demonstrated that, when studying electrochemical reactions, a high energy or short exposure time is needed to circumvent beam-induced effects.

Frank, Sara (ORCID:0000000163218363)↗

Implementing a Laser Stabilization System for Trapping Ca+ Ions: an Internship Reflection

At Lawrence Livermore National Laboratory, I contributed to a project developing 3D printed micro ion traps for quantum computing. I designed, implemented, and assessed a laser stabilization system that locked lasers to the frequencies required for calibrating our High Finesse WS8-10 wavelength meter and for laser cooling and trapping of Ca+ ions. I also programmed a Python interface for hardware communication, data collection, and statistical analysis. Additionally, I optimized and aligned laser beam paths, and I implemented a closed digital feedback loop using Proportional, Integral, and Derivative (PID) control parameters. I analyzed both the long-term and short-term behavior of our locked lasers and adjusted PID parameters to enhance performance. Furthermore, I used COMSOL to simulate the capacitance of a linear Paul trap design and predict our trap’s performance. The procedures I developed for the interface, analysis, and simulations will continue to support the ion trapping experiment after my appointment. I strengthened my skills in data analysis, Python coding, and optical alignment for laser systems. My confidence as a researcher grew, particularly in communicating my research. This experience taught me the importance of careful planning and consideration in research and solidified my desire to continue exploring novel quantum technology as an undergraduate

42 ENGINEERING↗

Linear and nonlinear dynamic analysis by boundary element method

An advanced implementation of the direct boundary element method (BEM) applicable to free-vibration, periodic (steady-state) vibration and linear and nonlinear transient dynamic problems involving two and three-dimensional isotropic solids of arbitrary shape is presented. Interior, exterior, and half-space problems can all be solved by the present formulation. For the free-vibration analysis, a new real variable BEM formulation is presented which solves the free-vibration problem in the form of algebraic equations (formed from the static kernels) and needs only surface discretization. In the area of time-domain transient analysis, the BEM is well suited because it gives an implicit formulation. Although the integral formulations are elegant, because of the complexity of the formulation it has never been implemented in exact form. In the present work, linear and nonlinear time domain transient analysis for three-dimensional solids has been implemented in a general and complete manner. The formulation and implementation of the nonlinear, transient, dynamic analysis presented here is the first ever in the field of boundary element analysis. Almost all the existing formulation of BEM in dynamics use the constant variation of the variables in space and time which is very unrealistic for engineering problems and, in some cases, it leads to unacceptably inaccurate results. In the present work, linear and quadratic isoparametric boundary elements are used for discretization of geometry and functional variations in space. In addition, higher order variations in time are used. These methods of analysis are applicable to piecewise-homogeneous materials, such that not only problems of the layered media and the soil-structure interaction can be analyzed but also a large problem can be solved by the usual sub-structuring technique. The analyses have been incorporated in a versatile, general-purpose computer program. Some numerical problems are solved and, through comparisons with available analytical and numerical results, the stability and high accuracy of these dynamic analysis techniques are established.

Ahmad, Shahid↗

Flexible aircraft dynamic modeling for dynamic analysis and control synthesis

The linearization and simplification of a nonlinear, literal model for flexible aircraft is highlighted. Areas of model fidelity that are critical if the model is to be used for control system synthesis are developed and several simplification techniques that can deliver the necessary model fidelity are discussed. These techniques include both numerical and analytical approaches. An analytical approach, based on first-order sensitivity theory is shown to lead not only to excellent numerical results, but also to closed-form analytical expressions for key system dynamic properties such as the pole/zero factors of the vehicle transfer-function matrix. The analytical results are expressed in terms of vehicle mass properties, vibrational characteristics, and rigid-body and aeroelastic stability derivatives, thus leading to the underlying causes for critical dynamic characteristics.

Schmidt, David K.↗

An Entropy-Based Approach to Nonlinear Stability

Many numerical methods used in computational fluid dynamics (CFD) incorporate an artificial dissipation term to suppress spurious oscillations and control nonlinear instabilities. The same effect can be accomplished by using upwind techniques, sometimes augmented with limiters to form Total Variation Diminishing (TVD) schemes. An analysis based on numerical satisfaction of the second law of thermodynamics allows many such methods to be compared and improved upon. A nonlinear stability proof is given for discrete scalar equations arising from a conservation law. Solutions to such equations are bounded in the L sub 2 norm if the second law of thermodynamics is satisfied in a global sense over a periodic domain. It is conjectured that an analogous statement is true for discrete equations arising from systems of conservation laws. Analysis and numerical experiments suggest that a more restrictive condition, a positive entropy production rate in each cell, is sufficient to exclude unphysical phenomena such as oscillations and expansion shocks. Construction of schemes which satisfy this condition is demonstrated for linear and nonlinear wave equations and for the one-dimensional Euler equations.

Merriam, Marshal L.↗

Prediction of Stable and Active (Oxy-Hydro) Oxide Nanoislands on Noble-Metal Supports for Electrochemical Oxygen Reduction Reaction

Developing cost-effective oxygen electrocatalysts with high activity and stability is key to their commercialization. However, economical earth-abundant catalysts based on first-row transition-metal oxides suffer from low electrochemical stability, which is difficult to improve without compromising their activity. In this paper, using density functional theory calculations, we demonstrate that noble-metal supports lead to bifunctional enhancement of both the stability and the oxygen reduction reaction (ORR) activity of metal (oxy-hydro) oxide nanoislands. We observe a significant stabilization of supported nanoislands beyond the intrinsic stability limits of bulk phases, which originates from a favorable lattice mismatch and reductive charge transfer from oxophilic supports. Additionally, we discover that interfacial active sites (located between the nanoisland and the support) reinforce the binding strength of reaction intermediates, hence boosting ORR activity. Considering that both stability and activity lead to discovery of CoOOH|Pt, NiOOH|Ag, and FeO 2 |Ag as viable systems for alkaline ORR, we then use a multivariant linear regression method to identify elementary descriptors for efficient screening of promising cost-effective nanoisland|support catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗