Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “perturbation 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 865 records · Page 48

Subsonic panel method for designing wing surfaces from pressure distribution

An iterative method has been developed for designing wing section contours corresponding to a prescribed subcritical distribution of pressure. The calculations are initialized by using a surface panel method to analyze a baseline wing or wing-fuselage configuration. A first-order expansion to the baseline panel method equations is then used to calculate a matrix containing the partial derivative of potential at each control point with respect to each unknown geometry parameter. In every iteration cycle, the matrix is used both to calculate the geometry perturbation and to analyze the perturbed geometry. The distribution of potential on the perturbed geometry is established by simple linear extrapolation from the baseline solution. The extrapolated potential is converted to pressure by Bernoulli's equation. Not only is the accuracy of the approach good for very large perturbations, but the computing cost of each complete iteration cycle is substantially less than one analysis solution by a conventional panel method.

Bristow, D. R.↗

Dynamic analysis of turbulent annular seals based on Hirs' lubrication equation

Expressions are derived which define dynamic coefficients for high-pressure annular seals typical of neck-ring and interstage seals employed in multistage centrifugal pumps. Completely developed turbulent flow is assumed in both the circumferential and axial directions, and is modeled in this analysis by Hirs' turbulent lubrication equations. Linear zeroth and first-order 'short-bearing' perturbation solutions are developed by an expansion in the eccentricity ratio. The influence of inlet swirl is accounted for in the development of the circumferential flow field. Comparisons are made between the stiffness, damping, and inertia coefficients derived herein based on Hirs' model and previously published results based on other models. Finally, numerical results are presented for interstage seals in the Space Shuttle Main Engine High Pressure Fuel Turbopump and a water pump.

Childs, D. W.↗

Structural Tailoring of SSME Blades (vanes)

The engine blade design optimization program STAEBL (Structural Tailoring of Engine Blades) is available at the NASA Lewis computer facility. The analysis capabilities of this program were extended to typical loading conditions for SSME turbopump blades including thermal and pressure loading. Input files for representative SSME blade designs were developed and sample optimization studies for these blades completed. The structural tailoring program combines a general optimization package and a finite element blade analysis package. The analysis package's capabilities include natural frequency, maximum stress, and forced response computation, and fatigue life and flutter analysis. Optimization is performed using the feasible directions method. The current design is modified by perturbing the design variables so that the design constraints are satisfied while the objective function, such as blade weight, is reduced at the maximum rate. The program's geometric design variables include blade thickness distribution, thickness to chord ratios, and root chord.

Rubinstein, R.↗

Airborne measurements of stratospheric constituents over Antarctica in the austral spring 1987. I - Method and ozone observations

A Fourier transform spectrometer was flown aboard a DC-8 on 10 flights over Antarctica during August and September, 1987, as part of the Airborne Antarctic Ozone Experiment (AAOE). Observing the sun at infrared wavelengths, it was possible to determine the integrated column amount above the flight altitude for ozone and a number of other chemical species that are believed to be important in the perturbed chemistry of the 'ozone hole'. The paper describes the method, the observations, the data analysis procedure, and the ozone results. During the observation period, ozone developed a steep gradient near the edge of the polar vortex; deep within the vortex, the average ozone column decreased by about 1.6 percent per day during September.

Mankin, William G.↗

Numerical modeling of passive microwave O2 observations over precipitation

Radiometric observations at 118 GHz are compared to a planar-stratified numerical radiative transfer model, in order to study the microwave scattering and emission properties of precipitation cells. Agreement between observed and computed brightness perturbations is shown to be within 10 percent over the radiometrically opaque, mature regions of the cell. Transparent-channel cell top reflectivities of up to 50 percent in the convective core region and 6 percent in the anvil region are determined, using perturbed and nonscattering weighting functions. A parameterized rain cell model is used in a content analysis of 118 GHz precipitation cell spectra. The dominant 118-GHz spectral modes are found to contain useful data on the cell top altitude; they are applicable for detection of transparent anvil regions.

Gasiewski, A. J.↗

Interface stability of a fluid column subject to periodic accelerations normal to the longitudinal axis of the column

This project investigates the interface stability of an inviscid fluid column in the presence of a periodic acceleration field which has a component normal to the longitudinal axis of the column. A ramification of this configuration is that perturbations cannot be considered axisymmetric. The column is taken to be infinite in length. Floquet analysis will be utilized in the stability investigation.

Lyell, M. J.↗

High resolution numerical simulation of the linearized Euler equations in conservation law form

A linearized Euler solver based on a high resolution numerical scheme is presented. The approach is to linearize the flux vector as opposed to carrying through the complete linearization analysis with the dependent variable vector written as a sum of the mean and the perturbed flow. This allows the linearized equations to be maintained in conservation law form. The linearized equations are used to compute unsteady flows in turbomachinery blade rows arising due to blade vibrations. Numerical solutions are compared to theoretical results (where available) and to numerical solutions of the nonlinear Euler equations.

Sreenivas, Kidambi↗

Fully Threaded Tree for Adaptive Refinement Fluid Dynamics Simulations

A fully threaded tree (FTT) for adaptive refinement of regular meshes is described. By using a tree threaded at all levels, tree traversals for finding nearest neighbors are avoided. All operations on a tree including tree modifications are O(N), where N is a number of cells, and are performed in parallel. An efficient implementation of the tree is described that requires 2N words of memory. A filtering algorithm for removing high frequency noise during mesh refinement is described. A FTT can be used in various numerical applications. In this paper, it is applied to the integration of the Euler equations of fluid dynamics. An adaptive mesh time stepping algorithm is described in which different time steps are used at different l evels of the tree. Time stepping and mesh refinement are interleaved to avoid extensive buffer layers of fine mesh which were otherwise required ahead of moving shocks. Test examples are presented, and the FTT performance is evaluated. The three dimensional simulation of the interaction of a shock wave and a spherical bubble is carried out that shows the development of azimuthal perturbations on the bubble surface.

FINITE ELEMENT ANALYSIS↗

Engine System Loads Analysis Compared to Hot-Fire Data

Early implementation of structural dynamics finite element analyses for calculation of design loads is considered common design practice for high volume manufacturing industries such as automotive and aeronautical industries. However with the rarity of rocket engine development programs starts, these tools are relatively new to the design of rocket engines. In the NASA MC-1 engine program, the focus was to reduce the cost-to-weight ratio. The techniques for structural dynamics analysis practices, were tailored in this program to meet both production and structural design goals. Perturbation of rocket engine design parameters resulted in a number of MC-1 load cycles necessary to characterize the impact due to mass and stiffness changes. Evolution of loads and load extraction methodologies, parametric considerations and a discussion of load path sensitivities are important during the design and integration of a new engine system. During the final stages of development, it is important to verify the results of an engine system model to determine the validity of the results. During the final stages of the MC-1 program, hot-fire test results were obtained and compared to the structural design loads calculated by the engine system model. These comparisons are presented in this paper.

Frady, Gregory P.↗

Integrated Aerodynamic/Structural/Dynamic Analyses of Aircraft with Large Shape Changes

The conceptual and preliminary design processes for aircraft with large shape changes are generally difficult and time-consuming, and the processes are often customized for a specific shape change concept to streamline the vehicle design effort. Accordingly, several existing reports show excellent results of assessing a particular shape change concept or perturbations of a concept. The goal of the current effort was to develop a multidisciplinary analysis tool and process that would enable an aircraft designer to assess several very different morphing concepts early in the design phase and yet obtain second-order performance results so that design decisions can be made with better confidence. The approach uses an efficient parametric model formulation that allows automatic model generation for systems undergoing radical shape changes as a function of aerodynamic parameters, geometry parameters, and shape change parameters. In contrast to other more self-contained approaches, the approach utilizes off-the-shelf analysis modules to reduce development time and to make it accessible to many users. Because the analysis is loosely coupled, discipline modules like a multibody code can be easily swapped for other modules with similar capabilities. One of the advantages of this loosely coupled system is the ability to use the medium-to high-fidelity tools early in the design stages when the information can significantly influence and improve overall vehicle design. Data transfer among the analysis modules are based on an accurate and automated general purpose data transfer tool. In general, setup time for the integrated system presented in this paper is 2-4 days for simple shape change concepts and 1-2 weeks for more mechanically complicated concepts. Some of the key elements briefly described in the paper include parametric model development, aerodynamic database generation, multibody analysis, and the required software modules as well as examples for a telescoping wing, a folding wing, and a bat-like wing.

Samareh, Jamshid A.↗

Sloshing Gas in the Core of the Most Luminous Galaxy Cluster RXJ1347.5-1145

We present new constraints on the merger history of the most X-ray luminous cluster of galaxies, RXJ1347.5-1145, based its unique multiwavelength morphology. Our X-ray analysis confirms the core gas is undergoing "sloshing" resulting from a prior, large scale, gravitational perturbation. In combination with extensive multiwavelength observations, the sloshing gas points to the primary and secondary clusters having had at least two prior strong gravitational interactions. The evidence supports a model in which the secondary subcluster with mass M=4.8+/-2.4 x 10(exp 14) Stellar Mass has previously (> or approx.=0.6 Gyr ago) passed by the primary cluster, and has now returned for a subsequent crossing where the subcluster's gas has been completely stripped from its dark matter halo. RXJ1347 is a prime example of how core gas sloshing may be used to constrain the merger histories of galaxy clusters through multiwavelength analyses.

Markevitch, Maxim↗

Sloshing Gas in the Core of the Most Luminous Galaxy Cluster RXJ1347.5-1145

We present new constraints on the merger history of the most X-ray luminous cluster of galaxies, RXJ1347.5-1145, based on its unique multiwavelength morphology. Our X-ray analysis confirms the core gas is undergoing "sloshing" resulting from a prior, large scale, gravitational perturbation. In combination with extensive multiwavelength observations, the sloshing gas points to the primary and secondary clusters having had at least two prior strong gravitational interactions. The evidence supports a model in which the secondary subcluster with mass M=4.8+/-2.4x10(exp 14) solar Mass has previously (> or approx.0.6 Gyr ago) passed by the primary cluster, and has now returned for a subsequent crossing where the subcluster's gas has been completely stripped from its dark matter halo. RXJ1347 is a prime example of how core gas sloshing may be used to constrain the merger histories of galaxy clusters through multiwavelength analyses.

Johnson, Ryan E.↗

Singular Perturbation-Based Large-Signal Order Reduction of Microgrids for Stability and Accuracy Synthesis With Control

The increasing penetration of distributed energy resources (DERs) highlights the growing importance of microgrids (MGs) in enhancing power system reliability. Employing electromagnetic transient (EMT) analysis in MGs becomes crucial for controlling the rapid transients. However, this requires an accurate but high-order model of power electronics and their underlying control loops, complexifying the stability analysis from the viewpoint of a higher control level. To overcome these challenges, this paper proposes a large-signal order reduction (LSOR) method for MGs with considerations of external control inputs and the detailed dynamics of underlying control levels based on singular perturbation theory (SPT). Specially, we innovatively proposed and strictly proved a general stability and accuracy assessment theorem that allows us to analyze the dynamic stability of a full-order nonlinear system by only leveraging our derived reduced-order model (ROM) and boundary layer model (BLM). Furthermore, this theorem furnishes a set of conditions that determine the accuracy of the developed ROM. Lastly, by embedding such a theorem into the SPT, we propose a novel LSOR approach with guaranteed accuracy and stability analysis equivalence. Case studies are conducted on MG systems to show the effectiveness of the proposed approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Resilient Cooperative Voltage Control for Distribution Network with High Penetration Distributed Energy Resources

Here, this paper considers the problem of designing a resilient distributed voltage control algorithm for distribution systems with high penetration of distributed energy resources in the presence of an unknown cyber-attack. The purpose of the attack is to force the system to violate the operating voltage limit by intercepting its communication channels and inserting exogenous signals to perturb and/or modify the information being exchanged. We first review the cooperative voltage control proposed in our previous work and provide a new stability analysis for it. Next, we present a resilient cooperative voltage control algorithm by introducing a virtual system interconnected with the original system such that the voltage can be maintained within the operational limit under unknown attacks. The resiliency of the proposed algorithm is demonstrated via simulations on the IEEE 8500-node system when subjected to an attack which consists of corrupting the data being exchanged in the communication network between two generation units.

cooperative control↗

Efficient perturbation-tracking method for directly probing the spectral phonon properties from molecular dynamics simulations

Existing methods for directly extracting the spectral phonon properties from molecular dynamics (MD) simulations, like the normal mode analysis (NMA) and spectral energy density analysis, all require a very long simulation time to produce reliable results with good convergence. So far, these methods are mainly applied in studies using small systems and with empirical potentials, as the heavy computational load has greatly hindered their further applications. Here we propose a perturbation-tracking (PT) method for directly probing the mode-wise phonon anharmonic frequencies and lifetimes. We show that results obtained from our method are in excellent agreement with those from the conventional NMA approach, using Si as the model material system. Comparing with the NMA approach, the PT method offers a greater accuracy and significant improvement of efficiency. It takes an average of two orders of magnitude and up to three orders of magnitude less simulation time to obtain the same lifetime result of a phonon mode with intermediate to high accuracy. Meanwhile, our method preserves all the dynamics of probed phonon mode from a particular state, which means it is capable of studying the transient thermal transport processes in a nonequilibrium system. Besides the exceptional efficiency, our method also comes with freedom to choose to probe only those modes of interest. This makes it ideal for use with large systems and in computationally demanding applications, such as ab initio MD simulations. Moreover, the PT method we propose here is very straightforward and easy to implement.

74 ATOMIC AND MOLECULAR PHYSICS↗

Once-Through Steam Generator Model Analysis Using Python and Advanced Optimization Tools (Summer Internship Report)

This study focuses on the parametric analysis of design parameters for a once-through steam generator (OTSG) model, using python and advanced optimization tools to facilitate applications such as the flowing autoclave steam generator (FASG) test cases. Building on previous research involving another OTSG with a different design, this project aims to enhance our understanding of how steam generators (SGs) behave and how their outputs are influenced by changes in design. The reason for this design change is to allow for more precise modeling and optimization of SG performance, to provide a comparative analysis between the two designs, and to set up the model for integration with the FASG test case. The OTSG python-model is a mathematical representation (including fluid flow and heat transfer equations/models/correlations) of a steam-generating unit in a pressurized water reactor-type small modular reactor system. Design studies involve changing the model’s input design parameters to observe the resulting effects on the output of the system. By using advanced optimization tools, such as the Risk Analysis Virtual Environment (RAVEN) developed at Idaho National Laboratory, detailed design parametric studies and model optimization were performed. Six input parameters—pressure, temperature, and mass flow rate for the inlet of the primary-side (hot fluid) and secondary-side (cold fluid), respectively, of the SG—were randomly perturbed via RAVEN’s Monte Carlo Sampler module, using uniform distributions (i.e., ±1%, ±5% and ±10% relative changes) for 600 samples. The analysis provides valuable insights into SG optimization and can be used for sensor placement optimization to effectively monitor and obtain experimental data in other tests.

20 FOSSIL-FUELED POWER PLANTS↗

Boundary layer structure over the ocean observed by LEANDRE 1 during a tramontane event

A new airborne backscatter lidar, has been developed by CNRS (Service d'Aeronomie, (SA) Laboratoire de Meteorologie Dynamique (LMD) and the Institut des Sciences de 1'Univers) in the frame of the LEANDRE research program. It has been qualified on board the ARAT in autumn 1989 and spring 1990 and was involved in its first cooperative campaign during PYREX in October and November 1990. During this campaign, lidar observations of the perturbations induced on tropospheric flow and boundary layer structure were performed, and results are presented. A large number of experiments were performed, for synoptic situation description (meteorological radiosoundings, constant level balloons) and local flow analysis (aircrafts, radars, sodars). For the first time in such an experiment, a lidar has been flown on a research aircraft to perform altitude resolved observations of these perturbations, and we will present here results obtained for deflected flow structure. In the presence of a synoptic northerly flow, part of it is deflected to the east by the Pyrenees, and accelerated over the Mediterranean by the mountain surroundings. In this case, a low level wind is generated (the Tramontane) bringing cold and dry air over the Mediterranean Sea. As the sea is still at a warm temperature in November (around 17 degrees), an Internal Marine Boundary Layer rapidly grows over the first tens of kilometers and stabilizes at about 1 km depth, corresponding to an altitude just below the Lifting Condensation Level. The whole Marine Atmospheric Boundary Layer (MABL) is characterized by highly turbulent motions bringing large humid particles from the surface up to its top. The lidar signal due to scattering by these particles is then representative of the turbulent kinetic energy in this layer.

Flamant, C.↗

Evaluation of residual gas fraction estimation methods for cycle-to-cycle combustion variability analysis and modeling

Cycle-to-cycle combustion variability in spark-ignition engines during normal operation is mainly caused by random perturbations of the in-cylinder conditions such as the flow velocity field, homogeneity of the air-fuel distribution, spark energy discharge, and turbulence intensity of the flame front. Such perturbations translate into the variability of the energy released observed at the end of the combustion process. During normal operating conditions, the cycle-to-cycle variability (CCV) of the energy release behaves as random uncorrelated noise. However, during diluted combustion, in either the form of exhaust gas recirculation (EGR) or excess air (lean operation), the CCV tends to increase as dilution increases. Moreover, when the ignition limit is reached at high dilution levels, the combustion CCV is exacerbated by sporadic occurrences of incomplete combustion events, and the uncorrelation assumption no longer holds. The low or null energy released by partial burns and misfires has an impact on the following combustion event due to the residual gas that carries burned and unburned gases, which contributes to the deterministic coupling between engine cycles. Many residual gas fraction estimation methods, however, only address the nominal case where complete combustion occurs and combustion events are uncorrelated. Here we evaluate the efficacy of such methods on capturing the effects of partial burns and misfires on the residual gas estimate for high-EGR operation. The advantages and disadvantages of each method are discussed based on their ability to generate cycle-to-cycle estimates. Finally, a comparison between the different estimation techniques is presented based on their usefulness for control-oriented modeling.

42 ENGINEERING↗