Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear flow”

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 757 records · Page 42

On the nonlinear interfacial instability of rotating core-annular flow

The interfacial stability of rotating core-annular flows is investigated. The linear and nonlinear effects are considered for the case when the annular region is very thin. Both asymptotic and numerical methods are used to solve the flow in the core and film regions which are coupled by a difference in viscosity and density. The long-term behavior of the fluid-fluid interface is determined by deriving its nonlinear evolution in the form of a modified Kuramoto-Sivashinsky equation. We obtain a generalization of this equation to three dimensions. The flows considered are applicable to a wide array of physical problems where liquid films are used to lubricate higher or lower viscosity core fluids, for which a concentric arrangement is desired. Linearized solutions show that the effects of density and viscosity stratification are crucial to the stability of the interface. Rotation generally destabilizes non-axisymmetric disturbances to the interface, whereas the centripetal forces tend to stabilize flows in which the film contains the heavier fluid. Nonlinear affects allow finite amplitude helically travelling waves to exist when the fluids have different viscosities.

Coward, Aidrian V.↗

Nonlinear Nonmodal Analysis of Hypersonic Flow over Blunt Cones

The linear amplification of modal disturbances that lead to boundary-layer transition in two-dimensional/axisymmetric hypersonic configurations is strongly reduced by the presence of a blunt nosetip, and the mechanisms underlying the observed onset of transition over the cone frustum are currently unknown. Linear nonmodal analysis has shown that both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification for moderate to large nosetip bluntness. The present study extends the previous linear analysis by including the nonlinear effects. Specifically, the perturbation form of the 2D, harmonic Navier-Stokes equations (HNSE) are solved with a fully implicit formulation and the Newton-Raphson method. The increased number of degrees of freedom for the nonlinear system presents difficulties for solution strategies based on direct solution of the linearized system. Such difficulties are overcome by using the GMRES iterative method with a preconditioner corresponding to a simplified Jacobian without the cross derivative terms. The HNSE solver is verified by comparing with nonlinear parabolized stability equation (NPSE) results for the nonlinear evolution of planar waves in an incompressible Blasius boundary layer and in a Mach 6 flow over a blunt cone. Finally, nonlinear nonmodal results are presented for planar traveling disturbances over the blunt cone. The nonmodal analysis demonstrates that entropy-layer disturbances generated close to the nose tip can seed the amplification of higher frequency Mack’s second-mode instabilities further downstream.

boundary layer transition↗

Nonlinear Nonmodal Analysis of Hypersonic Flow over Blunt Cones

The linear amplification of modal disturbances that lead to boundary-layer transition in two-dimensional/axisymmetric hypersonic configurations is strongly reduced by the presence of a blunt nosetip, and the mechanisms underlying the observed onset of transition over the cone frustum are currently unknown. Linear nonmodal analysis has shown that both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification for moderate to large nosetip bluntness. The present study extends the previous linear analysis by including the nonlinear effects. Specifically, the perturbation form of the 2D, harmonic Navier-Stokes equations (HNSE) are solved with a fully implicit formulation and the Newton-Raphson method. The increased number of degrees of freedom for the nonlinear system presents difficulties for solution strategies based on direct solution of the linearized system. Such difficulties are overcome by using the GMRES iterative method with a preconditioner corresponding to a simplified Jacobian without the cross derivative terms. The HNSE solver is verified by comparing with nonlinear parabolized stability equation (NPSE) results for the nonlinear evolution of planar waves in an incompressible Blasius boundary layer and in a Mach 6 flow over a blunt cone. Finally, nonlinear nonmodal results are presented for planar traveling disturbances over the blunt cone. The nonmodal analysis demonstrates that entropy-layer disturbances generated close to the nose tip can seed the amplification of higher frequency Mack’s second-mode instabilities further downstream.

boundary layer transition↗

Residual turbulence from velocity shear stabilized interchange instabilities

The stabilizing effect of velocity shear on the macroscopic, broad bandwidth, ideal interchange instability is studied in linear and nonlinear regimes. A 2D dissipative magnetohydrodynamic (MHD) code is employed to simulate the system. For a given flow shear, V Prime , linear growth rates are shown to be suppressed to below the shear-free level at both the small and large wavelengths. With increasing V Prime , the unstable band in wavenumber-space shrinks so that the peak growth results for modes that correspond to relatively high wavenumbers, on the scale of the density gradient. In the nonlinear turbulent steady state, a similar turbulent spectrum obtains, and the convection cells are roughly circular. In addition, the density fluctuation level and the degree of flattening of the initial inverted density profile are found to decrease as V Prime increases; in fact, unstable modes are almost completely stabilized and the density profile reverts to laminar when V Prime is a few times the classic interchange growth rate. Moreover, the turbulent particle flux diminishes with increasing velocity shear such that all the flux is carried by the classical diffusive flux in the asymptotic limit. The simulations are compared with measurements of magnetic fluctuations from the Maryland Centrifugal Experiment, MCX, which investigated interchange modes in the presence of velocity shear. The experimental spectral data, taken in the plasma edge, are in general agreement with the numerical data obtained in higher viscosity simulations for which the level of viscosity is chosen consistent with MCX Reynolds numbers at the edge. In particular, the residual turbulence in both cases is dominated by elongated convection cells. Finally, concomitant Kelvin-Helmholtz instabilities in the system are also examined. Complete stability to interchanges is obtained only in the parameter space wherein the generalized Rayleigh inflexion theorem is satisfied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Assessing the Optimality of LinDist3Flow for Optimal Tap Selection of Step Voltage Regulators in Unbalanced Distribution Networks

The adoption of distributed energy resources such as photovoltaics (PVs) has increased dramatically during the previous decade. The increased penetration of PVs into distribution networks (DNs) can cause voltage fluctuations that have to be mitigated. One of the key utility assets employed to this end are step-voltage regulators (SVRs). It is desirable to include tap selection of SVRs in optimal power flow (OPF) routines, a task that turns out to be challenging because the resultant OPF problem is nonconvex with added complexities stemming from accurate SVR modeling. While several convex relaxations based on semi-definite programming (SDP) have been presented in the literature for optimal tap selection, SDP based schemes do not scale well and are challenging to implement in large-scale planning or operational frameworks. This paper deals with the optimal tap selection (OPTS) problem for wye-connected SVRs using linear approximations of power flow equations. Specifically, the LinDist3Flow model is adopted and the effective SVR ratio is assumed to be continuous- enabling the formulation of a problem called LinDist3Flow-OPTS, which amounts to a linear program. The scalability and optimality gap of LinDist3Flow-OPTS are evaluated with respect to existing SDP-based and nonlinear programming techniques for optimal tap selection in three standard feeders, namely, the IEEE 13-bus, 123-bus, and 8500-node DNs. For all DNs considered, LinDist3Flow-OPTS achieves an optimality gap of approximately 1% or less while significantly lowering the computational burden.

linear approximations↗

Assessing the Optimality of LinDist3Flow for Optimal Tap Selection of Step Voltage Regulators in Unbalanced Distribution Networks: Preprint

The adoption of distributed energy resources such as photovoltaics (PVs) has increased dramatically during the previous decade. The increased penetration of PVs into distribution networks (DNs) can cause voltage fluctuations that have to be mitigated. One of the key utility assets employed to this end are step-voltage regulators (SVRs). It is desirable to include tap selection of SVRs in optimal power flow (OPF) routines, a task that turns out to be challenging because the resultant OPF problem is nonconvex with added complexities stemming from accurate SVR modeling. While several convex relaxations based on semi-definite programming (SDP) have been presented in the literature for optimal tap selection, SDP based schemes do not scale well and are challenging to implement in large-scale planning or operational frameworks. This paper deals with the optimal tap selection (OPTS) problem for wye-connected SVRs using linear approximations of power flow equations. Specifically, the LinDist3Flow model is adopted and the effective SVR ratio is assumed to be continuous–enabling the formulation of a problem called LinDist3Flow-OPTS, which amounts to a linear program. The scalability and optimality gap of LinDist3Flow-OPTS are evaluated with respect to existing SDP-based and nonlinear programming techniques for optimal tap selection in three standard feeders, namely, the IEEE 13-bus, 123-bus, and 8500-node DNs. For all DNs considered, LinDist3Flow-OPTS achieves an optimality gap of approximately 1% or less while significantly lowering the computational burden.

linear approximations↗

An Analytical Investigation of an Oscillating Wedge in a Supersonic Perfect Gas Flow

Several aspects of the oscillating wedge are investigated to evaluate both the resulting trends for the wedge and methods of analyzing unsteady flows. An existing hypersonic small disturbance theory for an oscillating thin wedge is extended and applied. A perturbation method involving linearization about the known flow is then derived and discussed. Subsequently, a finite difference technique for calculating the complete unsteady flow field of the wedge in motion is presented and discussed in conjunction with some calculated quasi-static nonlinear trends.

Bennett, R. M.↗

An assessment of theoretical models for viscous and transonic flow

Topics discussed include shear flow models, simplified models for treating separation, classical linear theory, a local linearization theory, a transonic linear theory, a transonic nonlinear theory, the experiment of Davis, and the experiment of Tijdeman. It is concluded that shear flow models, which have proven very accurate in taking into account boundary layer effects for panel flutter, are likely to be less so for lifting surface flutter. For many applications in transonic flow, transonic linear theory will be adequate.

Dowell, E. H.↗

Effect of bulges on the stability of boundary layers

The instability of flows around hump and dip imperfections is investigated. The mean flow is calculated using interacting boundary layers, thereby accounting for viscous/inviscid interaction and separation bubbles. Then, the two-dimensional linear stability of this flow is analyzed, and the amplification factors are computed. Results are obtained for several height/width ratios and locations. The theoretical results have been used to correlate the experimental results of Walker and Greening (1942). The observed transition locations are found to correspond to amplification factors varying between 7.4 and 10.0, consistent with previous results for flat plates. The method accounts for both viscous and shear-layer instabilities. Separation is found to increase significantly the amplification factor.

Nayfeh, Ali H.↗

Kinetic simulations of the Kelvin-Helmholtz instability at the magnetopause

Two-dimensional hybrid simulations with particle ions and fluid electrons are used to calculate the kinetic evolution of the Kelvin-Helmholtz instability for a magnetopauselike configuration. The unidirectional magnetic field is essentially transverse to the plasma flow velocity, which is the most unstable case according to linear theory and models the flow dynamics in the subsolar region of the magnetopause for northward interplanetary magnetic field. We recover effects analogous to those found in MHD simulations, including a mode cascade to longer wavelengths. The boundary layer consists of coherent structure and is not well described by a diffusive process. Isolated structures on the order of the ion gyroradius are formed which can cross the boundary in either direction. We describe how the time evolution of these structures represents transport across boundary layers, and we consider the possible connection of these entities to flux transfer events and other structure seen in the low-latitude boundary layer at the Earth's magnetopause as well as to flux ropes commonly observed near the ionopause of Venus. We also discuss the relation of the hybrid calculations to previous MHD simulations and to observations.

Thomas, V. A.↗

On the nonlinear stability of viscous modes within the Rayleigh problem on an infinite flat plate

The stability has been investigated of the unsteady flow past an infinite flat plate when it is moved impulsively from rest, in its own plane. For small times the instantaneous stability of the flow depends on the linearized equations of motion which reduce in this problem to the Orr-Sommerfeld equation. It is known that the flow for certain values of Reynolds number, frequency and wave number is unstable to Tollmien-Schlichting waves, as in the case of the Blasius boundary layer flow past a flat plate. With increase in time, the unstable waves only undergo growth for a finite time interval, and this growth rate is itself a function of time. The influence of finite amplitude effects is studied by solving the full Navier-Stokes equations. It is found that the stability characteristics are markedly changed both by the consideration of the time evolution of the flow, and by the introduction of finite amplitude effects.

Webb, J. C.↗

User Guide for Compressible Flow Toolbox Version 2.1 for Use With MATLAB(Registered Trademark); Version 7

This report provides a user guide for the Compressible Flow Toolbox, a collection of algorithms that solve almost 300 linear and nonlinear classical compressible flow relations. The algorithms, implemented in the popular MATLAB programming language, are useful for analysis of one-dimensional steady flow with constant entropy, friction, heat transfer, or shock discontinuities. The solutions do not include any gas dissociative effects. The toolbox also contains functions for comparing and validating the equation-solving algorithms against solutions previously published in the open literature. The classical equations solved by the Compressible Flow Toolbox are: isentropic-flow equations, Fanno flow equations (pertaining to flow of an ideal gas in a pipe with friction), Rayleigh flow equations (pertaining to frictionless flow of an ideal gas, with heat transfer, in a pipe of constant cross section.), normal-shock equations, oblique-shock equations, and Prandtl-Meyer expansion equations. At the time this report was published, the Compressible Flow Toolbox was available without cost from the NASA Software Repository.

Melcher, Kevin J.↗

Tiltrotor Test Rig Aerodynamic Tares

The Tiltrotor Test Rig (TTR) is a NASA facility designed to test full-scale proprotors in the National Full-Scale Aerodynamics Facility (NFAC). The first wind-tunnel entry included measurements of aerodynamic loads (tares) on the spinner without a rotor. The TTR has an onboard balance for direct measurement of rotor loads, so these tares comprise only spinner loads. Spinner tares also affect shaft loads measured independently of the balance. The results show that in axial flow, spinner drag is a linear function of dynamic pressure as expected, but at other flow angles, spinner loads include substantial lift and drag. Furthermore, there is a pronounced Magnus effect when the spinner is sideways to the flow. There is a sharp change in aerodynamic load at 90-deg yaw angle, indicating separation. Equations were developed to characterize the spinner tares at any arbitrary combination of yaw angle and airspeed.

TTR↗

A scalable mixed-integer decomposition approach for optimal power system restoration

The optimal restoration problem lies at the foundation of the evaluation and improvement of resilience in power systems. In this paper we present a scalable decomposition algorithm, based on the integer L-shaped method, for solving this problem for realistic power systems. The algorithm works by partitioning the problem into a master problem and a slave problem. The master problem optimizes over energization sequences of generators, buses and branches, for the entire restoration horizon, while respecting stylized energization requirements. The slave problem verifies that there exist power flow solutions supporting the energization sequences proposed by the master problem, adding cuts to the master if an infeasible sequence is detected. The power flow verification is carried out using a piece-wise linear approximation of the power flow equations, capable of capturing overvoltage situations common in restoration procedures. We develop specialized binary cuts and hybrid Benders-binary cuts for separating infeasible islands, which are stronger than regular no-good cuts, as well as heuristics for obtaining feasible solutions. We present numerical results for the IEEE 39- and 118-bus systems, and for the Chilean 1500-bus grid, where our algorithm finds near-optimal restoration sequences orders of magnitude faster than the state-of-the-art, demonstrating the effectiveness of the proposed approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Turbulent Spot Pressure Fluctuation Wave Packet Model

Wave packet analysis provides a connection between linear small disturbance theory and subsequent nonlinear turbulent spot flow behavior. The traditional association between linear stability analysis and nonlinear wave form is developed via the method of stationary phase whereby asymptotic (simplified) mean flow solutions are used to estimate dispersion behavior and stationary phase approximation are used to invert the associated Fourier transform. The resulting process typically requires nonlinear algebraic equations inversions that can be best performed numerically, which partially mitigates the value of the approximation as compared to a more complete, e.g. DNS or linear/nonlinear adjoint methods. To obtain a simpler, closed-form analytical result, the complete packet solution is modeled via approximate amplitude (linear convected kinematic wave initial value problem) and local sinusoidal (wave equation) expressions. Significantly, the initial value for the kinematic wave transport expression follows from a separable variable coefficient approximation to the linearized pressure fluctuation Poisson expression. The resulting amplitude solution, while approximate in nature, nonetheless, appears to mimic many of the global features, e.g. transitional flow intermittency and pressure fluctuation magnitude behavior. A low wave number wave packet models also recover meaningful auto-correlation and low frequency spectral behaviors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The role of zonal flows in the saturation of multi-scale gyrokinetic turbulence

The 2D spectrum of the saturated electric potential from gyrokinetic turbulence simulations that include both ion and electron scales (multi-scale) in axisymmetric tokamak geometry is analyzed. The paradigm that the turbulence is saturated when the zonal (axisymmetic) ExB flow shearing rate competes with linear growth is shown to not apply to the electron scale turbulence. Instead, it is the mixing rate by the zonal ExB velocity spectrum with the turbulent distribution function that competes with linear growth. A model of this mechanism is shown to be able to capture the suppression of electron-scale turbulence by ion-scale turbulence and the threshold for the increase in electron scale turbulence when the ion-scale turbulence is reduced. The model computes the strength of the zonal flow velocity and the saturated potential spectrum from the linear growth rate spectrum. The model for the saturated electric potential spectrum is applied to a quasilinear transport model and shown to accurately reproduce the electron and ion energy fluxes of the non-linear gyrokinetic multi-scale simulations. The zonal flow mixing saturation model is also shown to reproduce the non-linear upshift in the critical temperature gradient caused by zonal flows in ion-scale gyrokinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Myocardial metabolism during exposure to carbon monoxide in the conscious dog.

Investigation of the relationship between coronary flow, heart rate, left ventricular function, and myocardial oxygen consumption at increasing levels of carboxyhemoglobin in conscious dogs. The results demonstrate a linear increase in coronary flow and heart rate as the carboxyhemoglobin increases up to 20%. Myocardial oxygen consumption declined during the same period.

Adams, J. D.↗

Fatigue life prediction in bending from axial fatigue information

Bending fatigue in the low cyclic life range differs from axial fatigue due to the plastic flow which alters the linear stress-strain relation normally used to determine the nominal stresses. An approach is presented to take into account the plastic flow in calculating nominal bending stress (S sub bending) based on true surface stress. These functions are derived in closed form for rectangular and circular cross sections. The nominal bending stress and the axial fatigue stress are plotted as a function of life (N sub S) and these curves are shown for several materials of engineering interest.

Manson, S. S.↗