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 973 records · Page 54

Bulk-Flow Analysis of Hybrid Thrust Bearings for Advanced Cryogenic Turbopumps

A bulk-flow analysis and computer program for prediction of the static load performance and dynamic force coefficients of angled injection, orifice-compensated hydrostatic/hydrodynamic thrust bearings have been completed. The product of the research is an efficient computational tool for the design of high-speed thrust bearings for cryogenic fluid turbopumps. The study addresses the needs of a growing technology that requires of reliable fluid film bearings to provide the maximum operating life with optimum controllable rotordynamic characteristics at the lowest cost. The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a set of bulk-flow mass and momentum conservation and energy transport equations. Mass flow conservation and a simple model for momentum transport within the hydrostatic bearing recesses are also accounted for. The bulk-flow model includes flow turbulence with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing recesses and film lands. The cryogenic fluid properties are obtained from realistic thermophysical equations of state. Turbulent bulk-flow shear parameters are based on Hirs' model with Moody's friction factor equations allowing a simple simulation for machined bearing surface roughness. A perturbation analysis leads to zeroth-order nonlinear equations governing the fluid flow for the thrust bearing operating at a static equilibrium position, and first-order linear equations describing the perturbed fluid flow for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth-order flow field equations renders the bearing flow rate, thrust load, drag torque and power dissipation. Solution to the first-order equations determines the axial stiffness, damping and inertia force coefficients. The computational method uses well established algorithms and generic subprograms available from prior developments. The Fortran9O computer program hydrothrust runs on a Windows 95/NT personal computer. The program, help files and examples are licensed by Texas A&M University Technology License Office. The study of the static and dynamic performance of two hydrostatic/hydrodynamic bearings demonstrates the importance of centrifugal and advection fluid inertia effects for operation at high rotational speeds. The first example considers a conceptual hydrostatic thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The large axial stiffness and damping coefficients of the bearing should provide accurate control and axial positioning of the turbopump and also allow for unshrouded impellers, therefore increasing the overall pump efficiency. The second bearing uses a refrigerant R134a, and its application in oil-free air conditioning compressors is of great technological importance and commercial value. The computed predictions reveal that the LH2 bearing load capacity and flow rate increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial stiffness has a maximum for a recess pressure rati of approx. 0.55. while the axial damping coefficient decreases as the recess pressure ratio increases. The computer results from three flow models are compared. These models are a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the other hand renders the largest damping coefficients and inertia coefficients. The most important findings are related to the reduction of the outflow through the inner radius and the appearance of subambient pressures. The performance of the refrigerant hybrid thrust bearing is evaluated at two operating speeds and pressure drops. The computed results are presented in dimensionless form to evidence consistent trends in the bearing performance characteristics. As the applied axial load increases, the bearing film thickness and flow rate decrease while the recess pressure increases. The axial stiffness coefficient shows a maximum for a certain intermediate load while the damping coefficient steadily increases. The computed results evidence the paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where there is actually an inflow through the bearing inner diameter, accompanied by subambient pressures just downstream of the bearing recess edge. These results are solely due to centrifugal fluid inertia and advection transport effects. Recommendations include the extension of the computer program to handle flexure pivot tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft angular misalignments.

SanAndres, Luis↗

The application of general aerodynamic lifting surface elements to problems in unsteady transonic flow

A study was conducted to investigate the feasibility of using combined subsonic and supersonic linear theory as a means for solving unsteady transonic flow problems in an economical and yet realistic manner. With some modification, existing linear theory methods are combined into a single program and a simple algorithm is derived for determining interference between lifting surface elements of different Mach number. The method is applied to a wide variety of problems for which measured unsteady pressure distributions and Mach number distributions are available. By comparing theory and experiment, the transonic method solutions show a significant improvement over uniform flow solutions. It is concluded that with these refinements the method will provide a means for performing realistic transonic flutter and dynamic response analyses at costs which are compatible with current linear theory based solutions.

Cunningham, A. M., Jr.↗

Growth and decay of disturbances in stratified shear flow in a rotating frame

The stability of stratified plane Couette flow in a rotating frame is investigated for a case in which the gravitational force is parallel to the rotation vector. Partial differential equations describing the behavior of disturbances in the linear regime are derived. Unstratified flow is stable as long as the angular momentum gradient is positive. If the gradient is negative, nonaxisymmetric disturbances grow as a power law in time, if the gradient is sufficiently steep. In flow which is unstable to convection, all perturbations asymptotically grow at the rate given by the Brunt-Vaisala frequency. If heat diffusion is included, all nonaxisymmetric perturbations now eventually decay as t exp -2, even if the flow is unstable to convection. If heat diffusion and viscosity are weak, nonaxisymmetric disturbances in convectively unstable flow will undergo a large transient growth before their eventual decay.

Korycansky, D. G.↗

Study on the well-posedness, convergence and the stability of the semi-implicit upwind numerical solver for the multi-fluid model

The well-posedness, convergence and the stability of the two-fluid code has been studied for a long time. Most of the investigations concern the semi-implicit upwind solution scheme for the six equation two-fluid model such as used in RELAP5 3 or TRACE 21. Since the system code, SPACE 2, adopts one more field, a droplet field, it consists of nine equations (3 mass, 3 momentum and 3 energy balance equations) and thus more involved investigations are necessary to confirm the stability and convergence. For this objective, the old issue of the well-posedness, convergence and the stability is revisited and some general guidelines to develop a well-posed numerical multi-fluid model are derived as follows; (1) Hyperbolicity of the corresponding system of partial differential equations is not a necessary condition for the development of a numerical model for multi-phase flow, but whether or not it is hyperbolic can provide guidance relative to initial conditions, boundary conditions, and expected high frequency behavior of the model. (2) A necessary condition for a well-posed numerical model is stability in the von Neumann sense, i.e. growth factor less than 1.0 for the shortest wave-length, 2{Delta}x. (3) The smallest node size used for convergence studies should be of the order of the characteristic dimension of the average description, i.e. smaller nodes can be used so long as they do not result in unphysical growth of wave-lengths less than the characteristic dimension. The usual mathematical definition of convergence i.e. the behavior of the solution as the node size approaches zero, is not appropriate for the discrete averaged numerical model, since there is diminished physical meaning to behavior at wavelengths less than the characteristic dimension of the average description. Under these guidelines, dispersion analysis and von Neumann stability analysis are performed for the three field multi-fluid, semi-implicit, upwind numerical model to show that the necessary conditions for well-posedness are met. To study the non-linear stability and the convergence, various runs with the torus problem using the SPACE code are utilized to confirm the frequency cascading effects that augment non-linear stability. A robust mechanism for the flow regime change is also a very important factor for developing non-linearly stable and convergent code. A open pipe flow problem is also simulated to investigate the non-linear stability effects in a flow geometry more typical of real safety code applications. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Effect of a forced Couette flow on coupled convective and morphological instabilities during unidirectional solidification

The effect of a forced Couette flow, parallel to a horizontal crystal-melt interface during directional solidification of an alloy of lead containing tin, on the onset of convective and morphological instabilities, is calculated numerically via a linear stability analysis. Such a flow does not affect perturbations with wave vectors perpendicular to the flow. For perturbations with wave vectors parallel to the flow, the onset of morphological instability is somewhat suppressed and thermosolutal convection is greatly suppressed. When instabilities occur, they are oscillatory and correspond to travelling waves. For values of the crystal growth velocity for which mixed morphological and convective modes occur, the presence of a forced flow produces sufficient decoupling to allow otherwise degenerate branches to be identified.

Coriell, S. R.↗

Transition correlation in subsonic flow over a flat plate

A 2D subsonic flow over a flat plate with a freestream Mach number M(infinity) up to 0.8 is considered. The flow can experience continuous uniform suction through the wall and the wall can be heated or cooled continuously with a fixed wall temperature. For a specific combination of M(infinity), suction velocity, and level of heat transfer, the mean flow problem is solved and linear stability calculations are performed to compute the location on the flat plate where the factor representing the integration of growth rates reaches nine. These calculations are repeated for several combinations of flow parameters and the theoretically predicted transition location is presented in the form of a correlation that can account for the effect of wall suction, heat transfer, and Mach number.

Masad, J. A.↗

Mechanically Accelerated Depolymerization of Entangled Linear Polymer Melts

Mechanical forces can enhance the chemical depolymerization of synthetic polymers when shear flow accelerates chain scission. To quantify the extent of mechanically-accelerated scission, the effect of simple shear flow (duration and strength) with low Weissenberg and Deborah numbers was investigated by considering the impact of applied work in both simple shear and shear dominated mixed flows. Hydrogenated polyisoprene was chosen as a model linear, entangled system. The conditions (strain amplitude, frequency, and shearing time) necessary to increase chain scission were assessed in the rubbery melt. Shear flow accelerated chain scission at higher temperatures, suggesting an activated process. Isothermal scission versus work curves were superposed by applying shift factors a T,S , whose Arrhenius-like temperature dependence gave an apparent activation energy for chain scission of ~ 110 kJ/mol, which is likely a combination of the activation energy of viscosity and bond energy. This work provides a base for quantifying the impact of shear on depolymerization of polymer melts and highlight the connection between viscous dissipation and scission chemistry.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Distribution-Agnostic Stochastic Optimal Power Flow for Distribution Grids

This paper outlines a data-driven, distributionally robust approach to solve chance-constrained AC optimal power flow problems in distribution networks. Uncertain forecasts for loads and power generated by photovoltaic (PV) systems are considered, with the goal of minimizing PV curtailment while meeting power flow and voltage regulation constraints. A data- driven approach is utilized to develop a distributionally robust conservative convex approximation of the chance-constraints; particularly, the mean and covariance matrix of the forecast errors are updated online, and leveraged to enforce voltage regulation with predetermined probability via Chebyshev-based bounds. By combining an accurate linear approximation of the AC power flow equations with the distributionally robust chance constraint reformulation, the resulting optimization problem becomes convex and computationally tractable.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Merging high resolution geophysical and geochemical surveys to reduce exploration risk at glass buttes, Oregon

The engineering and studies phase of the Glass Buttes project was aimed at reducing risk during the early stages of geothermal project development. The project’s inclusion of high resolution geophysical and geochemical surveys allowed Ormat to evaluate the value of these surveys both independently and in combination to quantify the most valuable course of action for exploration in an area where structure, permeability, and temperature are the most pressing questions. The sizes of the thermal anomalies at Glass Buttes are unusually large. Over the course of Phase I Ormat acquired high resolution LIDAR data to accurately map fault manifestations at the surface and collected detailed gravity and aeromagnetic surveys to map subsurface structural features. In addition, Ormat collected airborne hyperspectral data to assist with mapping the rock petrology and mineral alteration assemblages along Glass Buttes faults and magnetotelluric (MT) survey to try to better constrain the structures at depth. Direct and indirect identification of alteration assemblages reveal not only the geochemical character and temperature of the causative hydrothermal fluids but can also constrain areas of upflow along specific fault segments. All five datasets were merged along with subsurface lithologies and temperatures to predict the most likely locations for high permeability and hot fluids. The Glass Buttes temperature anomalies include 2 areas, totaling 60 km2 (23 mi2) of measured temperature gradients over 165° C/km (10° F/100ft). The Midnight Point temperature anomaly includes the Strat-1 well with 90°C (194 °F) at 603 m (1981 ft) with a 164 °C/km (10°F/100ft) temperature gradient at bottom hole and the GB-18 well with 71°C (160 °F) at 396 m (1300 ft) with a 182°C/km (11°F/100ft) gradient. The primary area of alteration and elevated temperature occurs near major fault intersections associated with Brothers Fault Zone and Basin and Range systems. Evidence for faulting is observed in each data set as follows. Field observations include fault plane orientations, complicated fault intersections, and hydrothermal alteration apparently pre-dating basalt flows. Geophysical anomalies include large, linear gradients in gravity and aeromagnetic data with magnetic lows possibly associated with alteration. Resistivity low anomalies also appear to have offsets associated with faulting. Hyperspectral and XRF identified alteration and individual volcanic flow units, respectively. When incorporated into a 3D geologic model, the fault intersections near the highest proven temperature and geophysical anomalies provide the first priority targets at Midnight Point. Ormat geologists selected the Midnight Point 52-33 drilling target based on a combination of pre-existing drilling data, geologic field work, geophysical interpretation, and geochemical analysis. Deep temperatures of well 52-33 was lower than anticipated. Temperature gradients in the well mirrored those found in historical drilling, but they decreased below 1500 ft and were isothermal below 2000 ft.

15 GEOTHERMAL ENERGY↗

The effects of suction on the nonlinear stability of the three-dimensional boundary layer above a rotating disc

There exist two types of stationary instability of the flow over a rotating disc corresponding to the upper branch, inviscid mode and the lower branch mode, which has a triple deck structure, of the neutral stability curve. A theoretical study of the linear problem and an account of the weakly nonlinear properties of the lower branch modes have been undertaken by Hall and MacKerrell respectively. Motivated by recent reports of experimental sightings of the lower branch mode and an examination of the role of suction on the linear stability properties of the flow here, the effects are studied of suction on the nonlinear disturbance described by MacKerrell. The additional analysis required in order to incorporate suction is relatively straightforward and enables the derivation of an amplitude equation which describes the evolution of the mode. For each value of the suction, a threshold value of the disturbance amplitude is obtained; modes of size greater than this threshold grow without limit as they develop away from the point of neutral stability.

Bassom, Andrew P.↗

Predicting System Accidents with Model Analysis During Hybrid Simulation

Standard discrete event simulation is commonly used to identify system bottlenecks and starving and blocking conditions in resources and services. The CONFIG hybrid discrete/continuous simulation tool can simulate such conditions in combination with inputs external to the simulation. This provides a means for evaluating the vulnerability to system accidents of a system's design, operating procedures, and control software. System accidents are brought about by complex unexpected interactions among multiple system failures , faulty or misleading sensor data, and inappropriate responses of human operators or software. The flows of resource and product materials play a central role in the hazardous situations that may arise in fluid transport and processing systems. We describe the capabilities of CONFIG for simulation-time linear circuit analysis of fluid flows in the context of model-based hazard analysis. We focus on how CONFIG simulates the static stresses in systems of flow. Unlike other flow-related properties, static stresses (or static potentials) cannot be represented by a set of state equations. The distribution of static stresses is dependent on the specific history of operations performed on a system. We discuss the use of this type of information in hazard analysis of system designs.

Malin, Jane T.↗

A Generalized Copula-Polynomial Chaos Expansion for Probabilistic Power Flow Considering Nonlinear Correlations of PV Injections

This paper develops a generalized Copula-polynomial chaos expansion (PCE) framework for power system probabilistic power flow that can handle both linear and nonlinear correlations of uncertain power injections, such as wind and PVs. A data-driven Copula statistical model is used to capture the correlations of uncertain power injections. This allows us to resort to the Rosenblatt transformation to transform correlated variables into independent ones while preserving the dependence structure. This paves the way of leveraging the PCE for surrogate modeling and uncertainty quantification of power flow results, i.e., achieving the probabilistic distributions of power flows. Simulations carried out on the IEEE 57-bus system show that the proposed framework can get much more accurate results than other alternatives with different linear and nonlinear power injection correlations.

copula↗

Study of ion-gyroscale fluctuations in low-density L-mode plasmas heated by NBI on KSTAR

Broadband density fluctuations with peak frequency ranging from 150 to 400 kHz were measured using a multichannel microwave imaging reflectometer in core region of the low-density L-mode plasmas heated by neutral beam injection on KSTAR. These fluctuations have been studied by comparing the dominant mode scales estimated from the measurement with those predicted from linear gyrokinetic simulation. The measured poloidal wavenumbers are qualitatively comparable to those of the 'fastest growing modes' from simulations, whereas they are larger than those of the 'transport-dominant modes' by about a factor of three. The agreement on wavenumbers between the measurement and linear simulation (for the fastest growing modes) is probably due to sufficiently weak $E \times B$ flow shear compared to the maximum linear growth rate. Meanwhile, the transport-dominant modes seem to be related to the fluctuations in lower frequencies (~80–150 kHz) observed in some of the measurement.

Physics↗

A Hierarchical OPF Algorithm with Improved Gradient Evaluation in Three-Phase Networks

Linear approximation commonly used in solving alternating-current optimal power flow (AC-OPF) simplifies the system models but incurs accumulated voltage errors in large power networks. Such errors will make the primal-dual type gradient algorithms converge to solutions with voltage violation. In this paper, we improve a recent hierarchical OPF algorithm that rested on primal-dual gradients evaluated with a linearized distribution power flow model. Specifically, we propose a more accurate gradient evaluation method based on an unbalanced three-phase nonlinear distribution power flow model to mitigate the errors arising from linearization. The resultant gradients feature a blocked structure that enables our development of an improved hierarchical primal-dual algorithm to solve the OPF problem. Numerical results on the IEEE 123-bus test feeder and a 4,518-node test feeder show that the proposed method can enhance voltage safety at comparable computational efficiency with the linearized algorithm.

approximation algorithms↗

A simple, accurate and efficient algorithm for unsteady transonic flow

The derivation of an algorithm for computing an airfoil's response to small unsteady perturbations is described. The analysis of the far-field boundary conditions for the flow field, and the application of the time-linearization technique to develop an accurate description of the flow field are discussed. The treatment of moving shock waves is examined. The alternating direction implicit scheme of Ballhaus and Steger (1975) is utilized for the discretization of the governing equations and boundary conditions. Numerical examples demonstrating the applicability of the time-linearization method to determine the correlation between airfoil and pressure distribution are provided. The effect of wind tunnel walls on unsteady transonic testing is studied.

Fung, K.-Y.↗

Time integration for diffuse interface models for two-phase flow

We propose a variant of the θ-scheme for diffuse interface models for two-phase flow, together with three new linearization techniques for the surface tension. These involve either additional stabilizing force terms, or a fully implicit coupling of the Navier–Stokes and Cahn–Hilliard equation. In the common case that the equations for interface and flow are coupled explicitly, we find a time step restriction which is very different to other two-phase flow models and in particular is independent of the grid size. We also show that the proposed stabilization techniques can lift this time step restriction. Even more pronounced is the performance of the proposed fully implicit scheme which is stable for arbitrarily large time steps. We demonstrate in a Taylor-flow application that this superior coupling between flow and interface equation can decrease the computation time by several orders of magnitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗