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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Motion and decay of vortex rings submerged in a rotational flow

The interaction between vortex rings of finite strength and an axisymmetric rotational background flow is studied by a singular perturbation method, because it is difficult to use a finite-difference method to analyze the viscous decay in the small core of a vortex ring. The analysis is carried out by combining a composite solution of a vortex ring and an unsteady Euler solution for the background rotational flow. Using the method of averaging, a numerical scheme is developed to obtain an Euler solution in which the grid and time-step sizes depend solely on the length and velocity scales of the background flow. Numerical results are presented to illustrate the interaction between the trajectories and decay rates of the vortex rings and the background rotational flow.

Ishii, K.↗
NASA NTRS1987

Dynamic characteristics of a vibrating beam with periodic variation in bending stiffness

A detailed dynamic analysis is performed of a vibrating beam with bending stiffness periodic in the spatial coordinate. The effects of system parameters on beam response are explored with a perturbation expansion technique. It is found that periodic stiffness acts to modulate the modal displacements from the characteristic shape of a simple sine wave. The results are verified by a finite element solution and through experimental testing.

Townsend, John S.↗
NASA NTRS1988

Optimal reentry guidance for aeroassisted orbit transfer vehicles

A three-state model is presented for analyzing the problem of optimal changes in heading with minimum energy loss for a hypersonic gliding vehicle. A further model order reduction to a single state model is examined using singular perturbation theory. The optimal solution for the reduced problem defines an optimal altitude profile dependent on the current energy of the vehicle. A separate boundary-layer analysis is used to account for altitude and flight path angle dynamics, and to obtain lift and bank angle control solutions. By considering alternative approximations to solve the boundary-layer problem, three guidance laws are obtained, each having a feedback form. The guidance laws are evaluated for a hypothetical vehicle, and compared to an optimal solution obtained using a multiple shooting algorithm.

Calise, Anthony J.↗
NASA NTRS1954

Theoretical study of the transonic lift of a double-wedge profile with detached bow wave

A theoretical study is described of the aerodynamic characteristics at small angle of attack of a thin, double-wedge profile in the range of supersonic flight speed in which the bow wave is detached. The analysis is carried out within the framework of the transonic (nonlinear) small-disturbance theory, and the effects of angle of attack are regarded as a small perturbation on the flow previously calculated at zero angle. The mixed flow about the front half of the profile is calculated by relaxation solution of a suitably defined boundary-value problem for transonic small-disturbance equation in the hodograph plane (i.e., the Tricomi equation). The purely supersonic flow about the rear half is found by an extension of the usual numerical method of characteristics. Analytical results are also obtained, within the framework of the same theory, for the range of speed in which the bow wave is attached and the flow is completely supersonic.

Vincenti, Walter G↗
NASA NTRS1956

Lift hysteresis at stall as an unsteady boundary-layer phenomenon

Analysis of rotating stall of compressor blade rows requires specification of a dynamic lift curve for the airfoil section at or near stall, presumably including the effect of lift hysteresis. Consideration of the magnus lift of a rotating cylinder suggests performing an unsteady boundary-layer calculation to find the movement of the separation points of an airfoil fixed in a stream of variable incidence. The consideration of the shedding of vorticity into the wake should yield an estimate of lift increment proportional to time rate of change of angle of attack. This increment is the amplitude of the hysteresis loop. An approximate analysis is carried out according to the foregoing ideas for a 6:1 elliptic airfoil at the angle of attack for maximum lift. The assumptions of small perturbations from maximum lift are made, permitting neglect of distributed vorticity in the wake. The calculated hysteresis loop is counterclockwise. Finally, a discussion of the forms of hysteresis loops is presented; and, for small reduced frequency of oscillation, it is concluded that the concept of a viscous "time lag" is appropriate only for harmonic variations of angle of attack with time at mean conditions other than maximum lift.

Moore, Franklin K↗
NASA NTRS1994

Airfoil optimization by the one-shot method

An efficient numerical approach for the design of optimal aerodynamic shapes is presented in this paper. The objective of any optimization problem is to find the optimum of a cost function subject to a certain state equation (Governing equation of the flow field) and certain side constraints. As in classical optimal control methods, the present approach introduces a costate variable (Language multiplier) to evaluate the gradient of the cost function. High efficiency in reaching the optimum solution is achieved by using a multigrid technique and updating the shape in a hierarchical manner such that smooth (low-frequency) changes are done separately from high-frequency changes. Thus, the design variables are changed on a grid where their changes produce nonsmooth (high-frequency) perturbations that can be damped efficiently by the multigrid. The cost of solving the optimization problem is approximately two to three times the cost of the equivalent analysis problem.

Kuruvila, G.↗
NASA NTRS1997

Linearized Poststall Aerodynamic and Control Law Models of the X-31A Aircraft and Comparison with Flight Data

The X-31A aircraft has a unique configuration that uses thrust-vector vanes and aerodynamic control effectors to provide an operating envelope to a maximum 70 deg angle of attack, an inherently nonlinear portion of the flight envelope. This report presents linearized versions of the X-31A longitudinal and lateral-directional control systems, with aerodynamic models sufficient to evaluate characteristics in the poststall envelope at 30 deg, 45 deg, and 60 deg angle of attack. The models are presented with detail sufficient to allow the reader to reproduce the linear results or perform independent control studies. Comparisons between the responses of the linear models and flight data are presented in the time and frequency domains to demonstrate the strengths and weaknesses of the ability to predict high-angle-of-attack flight dynamics using linear models. The X-31A six-degree-of-freedom simulation contains a program that calculates linear perturbation models throughout the X-31A flight envelope. The models include aerodynamics and flight control system dynamics that are used for stability, controllability, and handling qualities analysis. The models presented in this report demonstrate the ability to provide reasonable linear representations in the poststall flight regime.

Stoliker, Patrick C.↗
NASA NTRS1995

Airfoil Design and Optimization by the One-Shot Method

An efficient numerical approach for the design of optimal aerodynamic shapes is presented in this paper. The objective of any optimization problem is to find the optimum of a cost function subject to a certain state equation (governing equation of the flow field) and certain side constraints. As in classical optimal control methods, the present approach introduces a costate variable (Lagrange multiplier) to evaluate the gradient of the cost function. High efficiency in reaching the optimum solution is achieved by using a multigrid technique and updating the shape in a hierarchical manner such that smooth (low-frequency) changes are done separately from high-frequency changes. Thus, the design variables are changed on a grid where their changes produce nonsmooth (high-frequency) perturbations that can be damped efficiently by the multigrid. The cost of solving the optimization problem is approximately two to three times the cost of the equivalent analysis problem.

Kuruvila, G.↗
NASA NTRS2005

Integration of an Autopilot for a Micro Air Vehicle

Two autopilots providing autonomous flight capabilities are presented herein. The first is the Pico-Pilot, demonstrated for the 12-inch size class of micro air vehicles. The second is the MicroPilot MP2028(sup g), where its integration into a 36-inch Zagi airframe (tailless, elevons only configuration) is investigated and is the main focus of the report. Analytical methods, which include the use of the Advanced Aircraft Analysis software from DARCorp, were used to determine the stability and control derivatives, which were then validated through wind tunnel experiments. From the aerodynamic data, the linear, perturbed equations of motion from steady-state flight conditions may be cast in terms of these derivatives. Using these linear equations, transfer functions for the control and navigation systems were developed and feedback control laws based on Proportional, Integral, and Derivative (PID) control design were developed to control the aircraft. The PID gains may then be programmed into the autopilot software and uploaded to the microprocessor of the autopilot. The Pico-Pilot system was flight tested and shown to be successful in navigating a 12-inch MAV through a course defined by a number of waypoints with a high degree of accuracy, and in 20 mph winds. The system, though, showed problems with control authority in the roll and pitch motion of the aircraft: causing oscillations in these directions, but the aircraft maintained its heading while following the prescribed course. Flight tests were performed in remote control mode to evaluate handling, adjust trim, and test data logging for the Zagi with integrated MP2028(sup g). Ground testing was performed to test GPS acquisition, data logging, and control response in autonomous mode. Technical difficulties and integration limitations with the autopilot prevented fully autonomous flight from taking place, but the integration methodologies developed for this autopilot are, in general, applicable for unmanned air vehicles within the 36-inch size class or larger that use a PID control based autopilot.

Platanitis, George↗
DOE OSTI2024

Uncertainty Quantification and Sensitivity Analysis of Non-Nuclear Advanced Controls Testbed Reactor Mockup

The research presented in this report describes our progress in applying stochastic methods and uncertainty quantification, parametric study, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated with Griffin, with the goal of developing a reduced order (surrogate) model which can be rapidly sampled while perturbing multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums, which are controlled by a hybrid proportional, integral and derivative (PID) controller, actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time and $\beta$ from a transient Griffin simulation executed with the improved quasi-static method to provide the kinetic parameters) as inputs to functions which control the CD rotation angle. Using a number of stochastic method approaches, we developed a dual purpose training-surrogate model of the NTP system using polynomial regression. The trained model can be rapidly sampled while simultaneously perturbing various input parameters of the model, such as coefficients on the PID control, or temperature (directly affect the neutron cross section). The surrogate model delivers accurate results orders-of-magnitude faster (minutes, not days) than the base model. Once the base model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs and their effect on the output. For example, coefficients used in the PID control system may vary due to some physical interference, or there may be uncertainty in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters and operational boundaries can be determined. The goal of this work is to support development of an advanced control system to operate CDs in a functioning NTP system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Geometry and dynamics of lipid membranes: The Scriven-Love number

The equations governing lipid membrane dynamics in planar, spherical, and cylindrical geometries are presented here. Unperturbed and first-order perturbed equations are determined and nondimensionalized. In membrane systems with a nonzero base flow, perturbed in-plane and out-of-plane quantities are found to vary over different length scales. Here, a new dimensionless number, named the Scriven-Love number, and the well-known Föppl–von Kármán number result from a scaling analysis. The Scriven-Love number compares out-of-plane forces arising from the in-plane, intramembrane viscous stresses to the familiar elastic bending forces, while the Föppl–von Kármán number compares tension to bending forces. Both numbers are calculated in past experimental works, and span a wide range of values in various biological processes across different geometries. In situations with large Scriven-Love and Föppl–von Kármán numbers, the dynamical response of a perturbed membrane is dominated by out-of-plane viscous and surface tension forces—with bending forces playing a negligible role. Calculations of non-negligible Scriven-Love numbers in various biological processes and in vitro experiments show in-plane intramembrane viscous flows cannot generally be ignored when analyzing lipid membrane behavior.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2020

Quantification of the Land Surface and Brown Ocean Influence on Tropical Cyclone Intensification Over Land

An investigation of TC Kelvin in February 2018 over northeast Australia was conducted to understand the mechanisms of the Brown Ocean effect (BOE), and to develop a comprehensive analysis framework for land falling tropical cyclones (TCs) in the process. NASA’s Land Information System (LIS) coupled to the NASA Unified WRF (NU-WRF) system was employed as the numerical model framework for 12 land/soil moisture perturbation experiments. Impacts of soil moisture and surface enthalpy flux conditions on TC Kelvin were investigated by closely evaluating simulated track and intensity, midlevel atmospheric thermodynamic properties, vertical wind shear, total precipitable water (TPW), and surface moisture flux. The results suggest that there were recognized differentiations among the sensitivity simulations as a result of land surface (e.g., soil moisture and texture) conditions. However, the intensification of TC Kelvin over land was more strongly related to atmospheric moisture advection and the diurnal cycle of solar radiation (i.e., radiative cooling) than to overall soil moisture conditions or surface fluxes. The analysis framework employed here for TC Kelvin can serve as a foundation to specifically quantify the factors governing the BOE. It also demonstrates that the BOE is not a binary influence (i.e. all or nothing), but instead operates in a continuum from largely to minimally influential such that it could be utilized to help improve prediction of inland effects for all land falling TCs.

Jinwoong Yoo↗
DOE OSTI2026

Living on the edge: a non-perturbative resolution to the negativity of bulk entropies

Lin, Maldacena, Rozenberg, and Shan (LMRS) presented a new information paradox in black hole physics by noticing that the entanglement and Rényi entropies in a two-sided black hole can become negative when the geometry contains a very large number of matter excitations behind the black hole horizon. While originally this puzzle was presented in the context of BPS two-sided black holes in two-dimensional supergravity, the negativity in fact persists for more general two-sided black holes in the presence of a large number of matter excitations. Since the entanglement and Rényi entropies in ordinary quantum systems cannot be negative, resolving this puzzle is a necessary step towards understanding the quantum mechanical description of black holes. In this paper, we explain how to address the entanglement negativity puzzle, both in the original setting discussed by LMRS and in more general non-supersymmetric settings, by summing over all non-perturbative contributions to the gravitational path integral. We then interpret this result from the point of view of a dual matrix integral, which we use to extend our analysis beyond the regime of validity of the genus re-summation performed in the gravitational path integral. In this regime, positivity is rescued by new saddles of the matrix integral, a one-eigenvalue instanton and a two-eigenvalue instanton. Finally, we formulate a similar puzzle and its resolution using random tensor network techniques.

2D Gravity↗
DOE OSTI2024

Numerical study of implosion instability mitigation in magnetically driven solid liner dynamic screw pinches

Magnetized liner inertial fusion experiments on the Z accelerator suffer from magneto-Rayleigh–Taylor instabilities (MRTI) that compromise integrity of the imploding cylindrical liner, limiting achievable fusion fuel conditions and ultimately reducing magneto-inertial fusion target performance. Dynamic screw pinches (DSP) provide a method to reduce MRTI in-flight via application of magnetic field line tension to the imploding liner outer surface. In contrast with z-pinches that drive implosions with an azimuthal magnetic field, dynamic screw pinches enforce an additional axial drive magnetic field component, making the overall drive magnetic field helical. As the liner implodes, cumulative MRTI development is reduced by dynamically shifting the orientation of the fastest growing instability modes. Three-dimensional magnetohydrodynamic simulations show that the DSP mechanism effectively stabilizes initially solid cylindrical liner implosions driven by Z-scale current pulses, indicating that MRTI mitigation increases with the ratio of axial to azimuthal drive magnetic field components (i.e., the drive field ratio). We also performed a spectral analysis of the simulated imploding density distributions, extracting wavelength and pitch angle of the simulated MRTI structures to study their dynamics during the implosion. Simulations of liners initially perturbed with drive-field-aligned sinusoidal structures indicate that MRTI mitigation in DSP implosions decreases with perturbation wavelength, once again suggestive of magnetic field line tension effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Living on the edge: a non-perturbative resolution to the negativity of bulk entropies

Lin, Maldacena, Rozenberg, and Shan (LMRS) presented a new information paradox in black hole physics by noticing that the entanglement and Rényi entropies in a two-sided black hole can become negative when the geometry contains a very large number of matter excitations behind the black hole horizon. While originally this puzzle was presented in the context of BPS two-sided black holes in two-dimensional supergravity, the negativity in fact persists for more general two-sided black holes in the presence of a large number of matter excitations. Since the entanglement and Rényi entropies in ordinary quantum systems cannot be negative, resolving this puzzle is a necessary step towards understanding the quantum mechanical description of black holes. In this paper, we explain how to address the entanglement negativity puzzle, both in the original setting discussed by LMRS and in more general non-supersymmetric settings, by summing over all non-perturbative contributions to the gravitational path integral. We then interpret this result from the point of view of a dual matrix integral, which we use to extend our analysis beyond the regime of validity of the genus re-summation performed in the gravitational path integral. In this regime, positivity is rescued by new saddles of the matrix integral, a one-eigenvalue instanton and a two-eigenvalue instanton. Finally, we formulate a similar puzzle and its resolution using random tensor network techniques.

FOS: Physical sciences↗
NASA NTRS1975

Stochastic stability and instability of model ecosystems

In this work, we initiate a stability study of multispecies communities in stochastic environment by using Ito's differential equations as community models. By applying the direct method of Liapunov, we obtain sufficient conditions for stability and instability in the mean of the equilibrium populations. The conditions are expressed in terms of the dominant diagonal property of community matrices, which is a suitable mechanism for resolving the central problem of 'complexity vs stability' in model ecosystems. As a by-product of this analysis we exhibit important structural properties of the stochastic density-dependent models, and establish tolerance of community stability to a broad class of nonlinear time-varying perturbations.

Ladde, G. S.↗
NASA NTRS1984

Design sensitivity analysis with Applicon IFAD using the adjoint variable method

A numerical method is presented to implement structural design sensitivity analysis using the versatility and convenience of existing finite element structural analysis program and the theoretical foundation in structural design sensitivity analysis. Conventional design variables, such as thickness and cross-sectional areas, are considered. Structural performance functionals considered include compliance, displacement, and stress. It is shown that calculations can be carried out outside existing finite element codes, using postprocessing data only. That is, design sensitivity analysis software does not have to be imbedded in an existing finite element code. The finite element structural analysis program used in the implementation presented is IFAD. Feasibility of the method is shown through analysis of several problems, including built-up structures. Accurate design sensitivity results are obtained without the uncertainty of numerical accuracy associated with selection of a finite difference perturbation.

Frederick, Marjorie C.↗
NASA NTRS1985

Implementation of design sensitivity analysis with existing finite element codes

A numerical method is presented to implement structural design sensitivity analysis theory, using the versatility and convenience of existing finite element structural analysis programs. Design variables such as thickness and cross-sectional areas of components of individual members and built-up structures are considered. Structural performance functionals considered include displacement and stress. The method is also applicable for eigenvalue problem design sensitivity analysis. It is shown that calculations can be carried out outside existing finite element codes, using postprocessing data only. Thus, design sensitivity analysis software does not have to be embedded in an existing finite element code. Feasability of the method is shown through analysis of several problems, including a built-up structure. Accurate design sensitivity results are obtained without the uncertainty of numerical accuracy associated with selection of finite difference perturbations.

Choi, Kyung K.↗