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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Developing disturbances in shear

The origin and growth to moderate amplitude of disturbances in shear flow has been traditionally ascribed to the linear modal instability of the flow. Recent work on initial value problems has suggested that nonmodal growth of perturbations may be of equal and perhaps greater importance, at least in cases of rapid development. Examples of robust growth in model problems which support no instabilities and in baroclinic flows with realistic Ekman damping for which the exponential modes have zero or negative growth rates have been shown. Such examples have focused attention on the perturbations which are configured to tap the energy of mean flows. Here a critical examination of these favorably configured perturbations is given, making use of the simple constant free shear barotropic model which allows construction of exact two-dimensional isolated wave packet solutions.

Farrell, Brian↗
NASA NTRS1987

Concepts for radically increasing the numerical convergence rate of the Euler equations

Integral equation and finite difference methods have been developed for solving transonic flow problems using linearized forms of the transonic small disturbance and Euler equations. A key element is the use of a strained coordinate system in which the shock remains fixed. Additional criteria are developed to determine the free parameters in the coordinate straining; these free parameters are functions of the shock location. An integral equation analysis showed that the shock is located by ensuring that no expansion shocks exist in the solution. The expansion shock appears as oscillations in the solution near the sonic line, and the correct shock location is determined by removing these oscillations. A second objective was to study the ability of the Euler equation to model separated flow.

Nixon, David↗
NASA NTRS1992

Runge-Kutta methods combined with compact difference schemes for the unsteady Euler equations

Recent development using compact difference schemes to solve the Navier-Stokes equations show spectral-like accuracy. A study was made of the numerical characteristics of various combinations of the Runge-Kutta (RK) methods and compact difference schemes to calculate the unsteady Euler equations. The accuracy of finite difference schemes is assessed based on the evaluations of dissipative error. The objectives are reducing the numerical damping and, at the same time, preserving numerical stability. While this approach has tremendous success solving steady flows, numerical characteristics of unsteady calculations remain largely unclear. For unsteady flows, in addition to the dissipative errors, phase velocity and harmonic content of the numerical results are of concern. As a result of the discretization procedure, the simulated unsteady flow motions actually propagate in a dispersive numerical medium. Consequently, the dispersion characteristics of the numerical schemes which relate the phase velocity and wave number may greatly impact the numerical accuracy. The aim is to assess the numerical accuracy of the simulated results. To this end, the Fourier analysis is to provide the dispersive correlations of various numerical schemes. First, a detailed investigation of the existing RK methods is carried out. A generalized form of an N-step RK method is derived. With this generalized form, the criteria are derived for the three and four-step RK methods to be third and fourth-order time accurate for the non-linear equations, e.g., flow equations. These criteria are then applied to commonly used RK methods such as Jameson's 3-step and 4-step schemes and Wray's algorithm to identify the accuracy of the methods. For the spatial discretization, compact difference schemes are presented. The schemes are formulated in the operator-type to render themselves suitable for the Fourier analyses. The performance of the numerical methods is shown by numerical examples. These examples are detailed. described. The third case is a two-dimensional simulation of a Lamb vortex in an uniform flow. This calculation provides a realistic assessment of various finite difference schemes in terms of the conservation of the vortex strength and the harmonic content after travelling a substantial distance. The numerical implementation of Giles' non-refelctive equations coupled with the characteristic equations as the boundary condition is discussed in detail. Finally, the single vortex calculation is extended to simulate vortex pairing. For the distance between two vortices less than a threshold value, numerical results show crisp resolution of the vortex merging.

Yu, Sheng-Tao↗
DOE OSTI2019

Hydrodynamic description of long-distance spin transport through noncollinear magnetization states: Role of dispersion, nonlinearity, and damping

Nonlocal compensation of magnetic damping by spin injection has been theoretically shown to display dynamic, noncollinear magnetization states that carry spin currents over micrometer distances. Such states can be generically referred to as dissipative exchange flows (DEFs) because spatially diffusing spin currents are established by the mutual exchange torque exerted by neighboring spins. Analytical studies to date have been limited to the weak spin injection assumption whereby the equation of motion for the magnetization is mapped to hydrodynamic equations describing spin flow and then linearized. In this work, we analytically and numerically study easy-plane ferromagnetic channels subject to spin injection of arbitrary strength at one extremum under a unified hydrodynamic framework. We propose that DEFs generally exhibit a nonlinear profile along the channel accompanied by a nonlinear frequency tunability. At large injection strengths, we fully characterize a magnetization state we call a contact-soliton DEF (CS-DEF) composed of a stationary soliton at the injection site, which smoothly transitions into a DEF and exhibits a negative frequency tunability. The transition between a DEF and a CS-DEF occurs at the maximum precessional frequency and coincides with the Landau criterion: a subsonic to supersonic flow transition. Leveraging the hydraulic-electrical analogy, the current-voltage characteristics of a nonlinear DEF circuit are presented. Micromagnetic simulations of nanowires that include magnetocrystalline anisotropy and nonlocal dipole fields are in qualitative agreement with the analytical results. The magnetization states found here along with their characteristic profile and spectral features provide quantitative guidelines to pursue an experimental demonstration of DEFs in ferromagnetic materials and establish a unified description for long-distance spin transport.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1999

Disturbance Dynamics in Transitional and Turbulent Flows

In order to expand the predictive capability of single-point turbulence closure models to account for the early-stage transition regime, a methodology for the formulation and calibration of model equations for the ensemble-averaged disturbance kinetic energy and energy dissipation rate is presented. First the decay of laminar disturbances and turbulence in mean shear-free flows is studied. In laminar flows, such disturbances are linear superpositions of modes governed by the Orr-Sommerfeld equation. In turbulent flows, disturbances are described through transport equations for representative mean quantities. The link between a description based on a deterministic evolution equation and a probability based mean transport equation is established. Because an uncertainty in initial conditions exists in the laminar as well as the turbulent regime, a probability distribution must be defined even in the laminar case. Using this probability distribution, it is shown that the exponential decay of the linear modes in the laminar regime can be related to a power law decay of both the (ensemble) mean disturbance kinetic energy and the dissipation rate. The evolution of these mean disturbance quantities is then described by transport equations similar to those for the corresponding turbulent decaying flow. Second, homogeneous shear flow, where disturbances can be described by rapid distortion theory (RDT), is studied. The relationship between RDT and linear stability theory is exploited in order to obtain a closed set of modeled equations. The linear disturbance equations are solved directly so that the numerical simulation yields a database from which the closure coefficients in the ensemble-averaged disturbance equations can be determined.

Grosch, Chester E.↗
DOE OSTI2019

Hierarchical Distributed Voltage Regulation in Networked Autonomous Grids

We propose a novel algorithm to solve optimal power flow (OPF) that aims at dispatching controllable distributed energy resources (DERs) for voltage regulation at minimum cost. The proposed algorithm features unprecedented scalability to large distribution networks by utilizing an information structure based on networked autonomous grids (AGs). Specifically, each AG is a subtree of a large distribution network that has a tree topology. The topology and line parameters of each AG are known only to a regional coordinator (RC) that is responsible for communicating with and dispatching the DERs within this AG. The reduced network, where each AG is treated as a node, is managed by a central coordinator (CC), which knows the topology and line parameters of the reduced network only and communicates with all the RCs. We jointly explore this information structure and the structure of the linearized distribution power flow (LinDistFlow) model to derive a hierarchical, distributed implementation of the primal-dual gradient algorithm that solves the OPF. The proposed implementation significantly reduces the computation burden compared to the centrally coordinated implementation of the primal-dual algorithm. Numerical results on a 4,521-node test feeder show that the proposed hierarchical distributed algorithm can achieve an improvement of more than tenfold in the speed of convergence compared to the centrally coordinated primal-dual algorithm.

Zhou, Xinyang↗
DOE OSTI2019

Hierarchical Distributed Voltage Regulation in Networked Autonomous Grids: Preprint

We propose a novel algorithm to solve optimal power flow (OPF) that aims at dispatching controllable distributed energy resources (DERs) for voltage regulation at minimum cost. The proposed algorithm features unprecedented scalability to large distribution networks by utilizing an information structure based on networked autonomous grids (AGs). Specifically, each AG is a subtree of a large distribution network that has a tree topology. The topology and line parameters of each AG are known only to a regional coordinator (RC) that is responsible for communicating with and dispatching the DERs within this AG. The reduced network, where each AG is treated as a node, is managed by a central coordinator (CC), which knows the topology and line parameters of the reduced network only and communicates with all the RCs. We jointly explore this information structure and the structure of the linearized distribution power flow (LinDistFlow) model to derive a hierarchical, distributed implementation of the primal-dual gradient algorithm that solves OPF. The proposed implementation significantly reduces the computation burden compared to the centrally coordinated implementation of the primal-dual algorithm. Numerical results on a 4,521-node test feeder show that the proposed hierarchical distributed algorithm can achieve more than 10-fold improvement in the speed of convergence compared to the centrally coordinated primal-dual algorithm.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

Hierarchical Distributed Voltage Regulation in Networked Autonomous Grids

We propose a novel algorithm to solve optimal power flow (OPF) that aims at dispatching controllable distributed energy resources (DERs) for voltage regulation at minimum cost. The proposed algorithm features unprecedented scalability to large distribution networks by utilizing an information structure based on networked autonomous grids (AGs). Specifically, each AG is a subtree of a large distribution network that has a tree topology. The topology and line parameters of each AG are known only to a regional coordinator (RC) that is responsible for communicating with and dispatching the DERs within this AG. The reduced network, where each AG is treated as a node, is managed by a central coordinator (CC), which knows the topology and line parameters of the reduced network only and communicates with all the RCs. We jointly explore this information structure and the structure of the linearized distribution power flow (LinDistFlow) model to derive a hierarchical, distributed implementation of the primal-dual gradient algorithm that solves the OPF. The proposed implementation significantly reduces the computation burden compared to the centrally coordinated implementation of the primal-dual algorithm. Numerical results on a 4,521-node test feeder show that the proposed hierarchical distributed algorithm can achieve an improvement of more than tenfold in the speed of convergence compared to the centrally coordinated primal-dual algorithm.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1990

Accurate calculation and instability of supersonic wake flows

This study is concerned with the computation and linear stability of a class of laminar compressible wake flows. The emphasis is on correct basic flow profiles that satisfy the steady equations of motion, and to this end the unperturbed state is obtained through numerical integration of the compressible boundary-layer equations. The linear stability of the flow is examined via the Rayleigh equation that describes evolution of inviscid disturbances. Analytical results are given for short- and long-wavelength disturbances and some numerical results of the general eigenvalue problem are also reported.

Papageorgiou, Demetrius T.↗
DOE OSTI2018

Prediction of radon removal efficiency for a flow-through activated charcoal system and radon mitigation characteristics

Activated charcoal was widely used in radon researches and mitigation applications because of good adsorption ability. Here in the paper, a simple model that can predict radon removal efficiency (η) for radon mitigation applications using flow-through activated charcoal system was developed based on the previous works. Experiments were conducted with a flow-through activated charcoal system to validate the model and study the influence of factors on η. The experimental results of η are consistent with the theory to within 4% at the flow rates of 1–3 L/min, but range from 2% to 14% at the flow rates of 4–5 L/min. The equilibrium time of radon linearly decreases with rising flow rate, and η decreases with rising flow rate. With moisture competing with radon for adsorption sites on the activated charcoal, radon reaches the equilibrium concentration more quickly than moisture does. As relative humidity increases, η decreases, but the radon equilibrium concentration increases. The values of η are bigger for those smaller particle sizes with the same type of activated charcoal. The results of this paper are useful for designing activated charcoal systems at indoor environmental conditions and predicting η for radon mitigation applications.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Convex Optimization of Integrated Power-Gas Energy Flow Model With Applications to Probabilistic Energy Flow

Energy flow calculation is a fundamental problem of the integrated power and gas system (IPGS) operation and planning. However, the nonlinear gas flow model introduces major challenges to the energy flow calculation. In this paper, we propose a tractably convex optimization model to solve the energy flow problem in IPGSs. It is demonstrated that the proposed optimization model has the same optimal solution as the original nonlinear steady energy flow model. Also, piecewise linearization is adopted to tightly linearize the nonlinear objective function of the model, which transforms the formulated convex optimization into a linear program one. Thus, the computation complexity of the proposed energy flow model is significantly reduced as compared with the existing methods. In addition, the proposed model can be extended to probabilistic energy flow estimation. Extensive case studies are conducted to validate the effectiveness of the proposed energy flow model using three IPGSs.

42 ENGINEERING↗
NASA NTRS1994

Geometry of thin liquid sheet flows

Incompresible, thin sheet flows have been of research interest for many years. Those studies were mainly concerned with the stability of the flow in a surrounding gas. Squire was the first to carry out a linear, invicid stability analysis of sheet flow in air and compare the results with experiment. Dombrowski and Fraser did an experimental study of the disintegration of sheet flows using several viscous liquids. They also detected the formulation of holes in their sheet flows. Hagerty and Shea carried out an inviscid stability analysis and calculated growth rates with experimental values. They compared their calculated growth rates with experimental values. Taylor studied extensively the stability of thin liquid sheets both theoretically and experimentally. He showed that thin sheets in a vacuum are stable. Brown experimentally investigated thin liquid sheet flows as a method of application of thin films. Clark and Dumbrowski carried out second-order stability analysis for invicid sheet flows. Lin introduced viscosity into the linear stability analysis of thin sheet flows in a vacuum. Mansour and Chigier conducted an experimental study of the breakup of a sheet flow surrounded by high-speed air. Lin et al. did a linear stability analysis that included viscosity and a surrounding gas. Rangel and Sirignano carried out both a linear and nonlinear invisid stability analysis that applies for any density ratio between the sheet liquid and the surrounding gas. Now there is renewed interest in sheet flows because of their possible application as low mass radiating surfaces. The objective of this study is to investigate the fluid dynamics of sheet flows that are of interest for a space radiator system. Analytical expressions that govern the sheet geometry are compared with experimental results. Since a space radiator will operate in a vacuum, the analysis does not include any drag force on the sheet flow.

Chubb, Donald L.↗
NASA NTRS1996

Droplet Deformation in an Extensional Flow: The Role of Surfactant Physical Chemistry

Surfactant-induced Marangoni effects strongly alter the stresses exerted along fluid particle interfaces. In low gravity processes, these stresses can dictate the system behavior. The dependence of Marangoni effects on surfactant physical chemistry is not understood, severely impacting our ability to predict and control fluid particle flows. A droplet in an extensional flow allows the controlled study of stretching and deforming interfaces. The deformations of the drop allow both Marangoni stresses, which resist tangential shear, and Marangoni elasticities, which resist surface dilatation, to develop. This flow presents an ideal model system for studying these effects. Prior surfactant-related work in this flow considered a linear dependence of the surface tension on the surface concentration, valid only at dilute surface concentrations, or a non-linear framework at concentrations sufficiently dilute that the linear approximation was valid. The linear framework becomes inadequate for several reasons. The finite dimensions of surfactant molecules must be taken into account with a model that includes surfaces saturation. Nonideal interactions between adsorbed surfactant molecules alter the partitioning of surfactant between the bulk and the interface, the dynamics of surfactant adsorptive/desorptive exchange, and the sensitivity of the surface tension to adsorbed surfactant. For example, cohesion between hydrocarbon chains favors strong adsorption. Cohesion also slows the rate of desorption from interfaces, and decreases the sensitivity of the surface tension to adsorbed surfactant. Strong cohesive interactions result in first order surface phase changes with a plateau in the surface tension vs surface concentration. Within this surface concentration range, the surface tension is decoupled from surface concentration gradients. We are engaged in the study of the role of surfactant physical chemistry in determining the Marangoni stresses on a drop in an extensional flow in a numerical and experimental program. Using surfactants whose dynamics and equilibrium behavior have been characterized in our laboratory, drop deformation will be studied in ground-based experiment. In an accompanying numerical study, predictive drop deformations will be determined based on the isotherm and equation of state determined in our laboratory. This work will improve our abilities to predict and control all fluid particle flows.

Stebe, Kathleen J.↗
DOE OSTI2022

Floquet stability analyses of stratified oscillating boundary layers on adiabatic slopes

The presence of a no-slip, impermeable, adiabatic, sloped boundary in an otherwise quiescent, stably stratified, Boussinesq flow generates baroclinic vorticity within a diffusive boundary layer. Such conditions are typical of the oscillating boundary layers on adiabatic abyssal slopes, sloped lake bathymetry and sloped coastal bathymetry in the absence of high-wavenumber internal waves, mean flows, far-field turbulence on larger scales, and resonant tidal–bathymetric interaction. We investigate the linear stability of the oscillating flow within non-dimensional parameter space typical of the $M_2$ tide and hydraulically smooth, middle-latitude abyssal slopes through Floquet linear stability analysis. The flow dynamics depends on three non-dimensional variables: the Reynolds number for Stokes’ second problem, the Prandtl number, and a frequency ratio that accounts for the resonance conditions ($C$, criticality) of the buoyant restoring force and the tidal forcing. The Floquet analysis results suggest that oscillating laminar boundary layers on adiabatic abyssal slopes are increasingly unstable as Reynolds number, criticality parameter and/or spanwise disturbance wavenumber are increased. We also show that the two-dimensional Floquet linear instability necessarily generates three-dimensional baroclinic vorticity, which suggests that the evolution of the gravitational instabilities may be nonlinear as $t→∞$.

42 ENGINEERING↗
DOE OSTI2019

CERN Large Hadron Collider system size scan predictions for PbPb, XeXe, ArAr, and OO with relativistic hydrodynamics

In recent years the understanding on the limits of the smallest possible droplet of the Quark Gluon Plasma has been called into question. Experimental results from both the Large Hadron Collider and the Relativistic Heavy Ion Collider have provided hints that the Quark Gluon Plasma may be produced in systems as small as that formed in p Pb or d Au collisions. Yet alternative explanations still exist from correlations arising from quarks and gluons in a color glass condensate picture. In order to resolve these two scenarios, a system size scan has been proposed at the Large Hadron Collider for collisions of ArAr and OO. Here we make predictions for a possible future run of ArAr and OO collisions at the Large Hadron Collider and study the system size dependence of a variety of flow observables. We find that linear response (from the initial conditions to the final flow harmonics) becomes more dominant in smaller systems whereas linear + cubic response can accurately predict multi-particle cumulants for a wide range of centralities in large systems.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2005

Exposure to a Rotating Virtual Environment During Treadmill Locomotion Causes Adaptation in Heading Direction

The goal of the present study was to investigate the adaptive effects of variation in the direction of optic flow, experienced during linear treadmill walking, on modifying locomotor trajectory. Subjects (n = 30) walked on a motorized linear treadmill at 4.0 km/h for 24 minutes while viewing the interior of a 3D virtual scene projected onto a screen 1.5 m in front of them. The virtual scene depicted constant self-motion equivalent to either 1) walking around the perimeter of a room to one s left (Rotating Room group) 2) walking down the center of a hallway (Infinite Hallway group). The scene was static for the first 4 minutes, and then constant rate self-motion was simulated for the remaining 20 minutes. Before and after the treadmill locomotion adaptation period, subjects performed five stepping trials where in each trial they marched in place to the beat of a metronome at 90 steps/min while blindfolded in a quiet room. The subject s final heading direction (deg), final X (for-aft, cm) and final Y (medio-lateral, cm) positions were measured for each trial. During the treadmill locomotion adaptation period subject s 3D torso position was measured. We found that subjects in the Rotating Room group as compared to the Infinite Hallway group: 1) showed significantly greater deviation during post exposure testing in the heading direction and Y position opposite to the direction of optic flow experienced during treadmill walking 2) showed a significant monotonically increasing torso yaw angular rotation bias in the direction of optic flow during the treadmill adaptation exposure period. Subjects in both groups showed greater forward translation (in the +X direction) during the post treadmill stepping task that differed significantly from their pre exposure performance. Subjects in both groups reported no perceptual deviation in position during the stepping tasks. We infer that 3 viewing simulated rotary self-motion during treadmill locomotion causes adaptive modification of sensory-motor integration in the control of position and trajectory during locomotion which functionally reflects adaptive changes in the integration of visual, vestibular, and proprioceptive cues. Such an adaptation in the control of position and heading direction during locomotion due to the congruence of sensory information demonstrates the potential for adaptive transfer between sensorimotor systems and suggests a common neural site for the processing and self-motion perception and concurrent adaptation in motor output. This will result in lack of subjects perception of deviation of position and trajectory during the post treadmill step test while blind folded.

Mulavara, A. P.↗
NASA NTRS2005

Exposure to a Rotating Virtual Environment During Treadmill Locomotion Causes Adaptation in Heading Direction

The goal of the present study was to investigate the adaptive effects of variation in the direction of optic flow, experienced during linear treadmill walking, on modifying locomotor trajectory. Subjects (n = 30) walked on a motorized linear treadmill at 4.0 kilometers per hour for 24 minutes while viewing the interior of a 3D virtual scene projected onto a screen 1.5 in in front of them. The virtual scene depicted constant self-motion equivalent to either 1) walking around the perimeter of a room to one s left (Rotating Room group) 2) walking down the center of a hallway (Infinite Hallway group). The scene was static for the first 4 minutes, and then constant rate self-motion was simulated for the remaining 20 minutes. Before and after the treadmill locomotion adaptation period, subjects performed five stepping trials where in each trial they marched in place to the beat of a metronome at 90 steps/min while blindfolded in a quiet room. The subject's final heading direction (deg), final X (for-aft, cm) and final Y (medio-lateral, cm) positions were measured for each trial. During the treadmill locomotion adaptation period subject's 3D torso position was measured. We found that subjects in the Rotating Room group as compared to the Infinite Hallway group: 1) showed significantly greater deviation during post exposure testing in the heading direction and Y position opposite to the direction of optic flow experienced during treadmill walking 2) showed a significant monotonically increasing torso yaw angular rotation bias in the direction of optic flow during the treadmill adaptation exposure period. Subjects in both groups showed greater forward translation (in the +X direction) during the post treadmill stepping task that differed significantly from their pre exposure performance. Subjects in both groups reported no perceptual deviation in position during the stepping tasks. We infer that viewing simulated rotary self-motion during treadmill locomotion causes adaptive modification of sensory-motor integration in the control of position and trajectory during locomotion which functionally reflects adaptive changes in the integration of visual, vestibular, and proprioceptive cues. Such an adaptation in the control of position and heading direction during locomotion due to the congruence of sensory information demonstrates the potential for adaptive transfer between sensorimotor systems and suggests a common neural site for the processing and self-motion perception and concurrent adaptation in motor output. This will result in lack of subjects perception of deviation of position and trajectory during the post treadmill step test while blind folded.

Ruttley, T↗
DOE OSTI2019

Optimal Voltage Reference for Droop-Based DERs in Distribution Systems

Due to the newly issued standards on distributed energy resources (DERs), the droop-based operation modes (Q-V and P-V droops) become a requirement for all the interconnected DER units at the distribution level. Voltage references of the DER units can be assigned to improve the operational behavior of the distribution system in the smart grid context. This paper presents a novel formulation for optimally setting voltage reference of the droop-based DER units in distribution systems with the objective of minimizing the system ohmic loss. This problem is formulated as a mixed-integer quadratic convex optimization problem in MATLAB and solved by the MATLAB and CPLEX solvers. The IEEE 123-bus system and the EPRI test system Ckt5 are used to validate the proposed formulation. The result obtained from the proposed formulation is compared with the result of the original problem solved by the heuristic method. Based on the numerical study results, system loss can be reduced after the proposed voltage reference assignment. The computational expense of the new formulation is also much lower than that of the heuristic approach.

distributed energy resource integration↗