Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “power flow equation”

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 91 records · Page 5

A three-dimensional axisymmetric calculation procedure for turbulent flows in a radial vaneless diffuser

An analytical model is proposed to calculate the three-dimensional axisymmetric turbulent flowfield in a radial vneless diffuser. The model assumes that the radial and tangential boundary layer profiles are approximated by power law profiles. Then, using the integrated radial and tangential momentum and continuity equations for the boundary layer and corresponding inviscid equations for the core flow, there results six ordinary differential equation in six unknowns which are easily solved using a Runge-Kutta technique. A model is also proposed for fully developed flow. The results using this technique were compared with the results from a three-dimensional viscous, axisymmetric duct code and with experimental data and good quantitative agreement was obtained.

Schumann, L. F.↗

Solar Supergranulation Revealed as a Superposition of Traveling Waves

40 years ago two new solar phenomena were described: supergranulation and the five-minute solar oscillations. While the oscillations have since been explained and exploited to determine the properties of the solar interior, the supergranulation has remained unexplained. The supergranules, appearing as convective-like cellular patterns of horizontal outward flow with a characteristic diameter of 30 Mm and an apparent lifetime of 1 day, have puzzling properties, including their apparent superrotation and the minute temperature variations over the cells. Using a 60-day sequence of data from the MDI (Michelson-Doppler Imager) instrument onboard the SOHO (Solar and Heliospheric Observatory) spacecraft, we show that the supergranulation pattern is formed by a superposition of traveling waves with periods of 5-10 days. The wave power is anisotropic with excess power in the direction of rotation and toward the equator, leading to spurious rotation rates and north-south flows as derived from correlation analyses. These newly discovered waves could play an important role in maintaining differential rotation in the upper convection zone by transporting angular momentum towards the equator.

Gizon, L.↗

Estimating unsteady aerodynamic forces on a cascade in a three-dimensional turbulence field

An analytical method has been developed to estimate the unsteady aerodynamic forces caused by flow field turbulence on a wind tunnel turning vane cascade system (vane set). This method approximates dynamic lift and drag by linearly perturbing the appropriate steady state force equations, assuming that the dynamic loads are due only to free stream turbulence and that this turbulence is homogeneous, isotropic, and Gaussian. Correlation and unsteady aerodynamic effects are also incorporated into the analytical model. Using these assumptions, equations relating dynamic lift and drag to flow turbulence, mean velocity, and vane set geometry are derived. From these equations, estimates for the power spectra and rms (root mean squared value, delta) loading of both lift and drag can be determined.

Norman, T.↗

Convergence acceleration of iterative solutions of Euler equations for transonic flow computations

In this paper, two acceleration techniques for Euler calculations are investigated. The first technique is an extrapolation procedure based on the Power Method; it is applicable when the iterative matrix has dominant eigenvalues. Both real and complex conjugate roots are allowed. The second technique is a generalization of the Minimal Residual Method, where the extrapolation step consists of a weighted combination of the corrections at different iteration levels and the weights are chosen to minimize the Lz norm of the residual. Numerical results, using Jameson's Runge-Kutta Multigrid Code, are presented. The extra computational work to apply either technique is negligible and the extra storage is not a problem on current supercomputers.

Hafez, M.↗

Estimating unsteady aerodynamic forces on a cascade in a three-dimensional turbulence field

An analytical method has been developed to estimate tne unsteady aerodynamic forces caused by flow field turbulence on a wind tunnel turning vane cascade system (vane set). This method approximates dynamic lift and drag by linearly perturbing the appropriate steady state force equations, assuming that the dynamic loads are due only to free stream turbulence and that this turbulence is homogeneous, isotropic, and Gaussian. Correlation and unsteady aerodynamic effects are also incorporated into the analytical model. Using these assumptions, equations relating dynamic lift and drag to flow turbulence, mean velocity, and vane set geometry are derived. From these equations, estimates for the power spectra and rms (root mean squared value, delta) loading of both lift and drag can be determined.

Norman, T. R.↗

Second-order small-disturbance solutions for hypersonic flow over power-law bodies

Similarity solutions were found which give the adiabatic flow of an ideal gas about two-dimensional and axisymmetric power-law bodies at infinite Mach number to second order in the body slenderness parameter. The flow variables were expressed as a sum of zero-order and perturbation similarity functions for which the axial variations in the flow equations separated out. The resulting similarity equations were integrated numerically. The solutions, which are universal functions, are presented in graphic and tabular form. To avoid a singularity in the calculations, the results are limited to body power-law exponents greater than about 0.85 for the two-dimensional case and 0.75 for the axisymmetric case. Because of the entropy layer induced by the nose bluntness (for power-law bodies other than cones and wedges), only the pressure function is valid at the body surface. The similarity results give excellent agreement with the exact solutions for inviscid flow over wedges and cones having half-angles up to about 20 deg. They give good agreement with experimental shock-wave shapes and surface-pressure distributions for 3/4-power axisymmetric bodies, considering that Mach number and boundary-layer displacement effects are not included in the theory.

Townsend, J. C.↗

Advanced numerics for multi-dimensional fluid flow calculations

In recent years, there has been a growing interest in the development and use of mathematical models for the simulation of fluid flow, heat transfer and combustion processes in engineering equipment. The equations representing the multi-dimensional transport of mass, momenta and species are numerically solved by finite-difference or finite-element techniques. However despite the multiude of differencing schemes and solution algorithms, and the advancement of computing power, the calculation of multi-dimensional flows, especially three-dimensional flows, remains a mammoth task. The following discussion is concerned with the author's recent work on the construction of accurate discretization schemes for the partial derivatives, and the efficient solution of the set of nonlinear algebraic equations resulting after discretization. The present work has been jointly supported by the Ramjet Engine Division of the Wright Patterson Air Force Base, Ohio, and the NASA Lewis Research Center.

Vanka, S. P.↗

Aerodynamics of powered missile separation from a wing

A 3D dynamic 'chimera' algorithm that solves the thin-layer Navier-Stokes equations over multiple moving bodies was modified to numerically simulate the aerodynamics, missile dynamics, and missile plume of a finless missile separating from a wing in transonic flow. A powered missile separation case was considered to examine the influence of the missile and plume on the wing. The wing and missile is at a two degree angle of attack. The computational results show the details of the flow field.

Shanks, S. P.↗

The large Reynolds number - Asymptotic theory of turbulent boundary layers.

A self-consistent, asymptotic expansion of the one-point, mean turbulent equations of motion is obtained. Results such as the velocity defect law and the law of the wall evolve in a relatively rigorous manner, and a systematic ordering of the mean velocity boundary layer equations and their interaction with the main stream flow are obtained. The analysis is extended to the turbulent energy equation and to a treatment of the small scale equilibrium range of Kolmogoroff; in velocity correlation space the two-thirds power law is obtained. Thus, the two well-known 'laws' of turbulent flow are imbedded in an analysis which provides a great deal of other information.

Mellor, G. L.↗

Impact of Reordering on the LU Factorization Performance of Bordered Block-Diagonal Sparse Matrix

Power engineers rely on computer-based simulation tools to assess grid performance and ensure security. At the core of these tools are solvers for sparse linear equations. When transformed into a bordered block-diagonal (BBD) structure, part of the sparse linear equation solving can be parallelized. This work focuses on using the Schur-complement-based method for LU factorization on BBD matrices, specifically, Jacobian matrices from large-scale systems. Our findings show that the natural ordering method outperforms the default ordering method in computational performance for each block of the BBD matrix. This observation is validated using synthetic 25k-bus and 70k-bus cases, showing a speedup of up to 38% when using natural ordering without permutation. Additionally, the impact of the number of partitions is studied, and the result shows that computational performance improves with more, smaller partitions in the BBD matrices.

BBD matrix↗

Integral method for the calculation of three-dimensional, laminar and turbulent boundary layers

The method for turbulent flows is a further development of an existing method; profile families with two parameters and a lag entrainment method replace the simple entrainment method and power profiles with one parameter. The method for laminar flows is a new development. Moment of momentum equations were used for the solution of the problem, the profile families were derived from similar solutions of boundary layer equations. Laminar and turbulent flows at the wings were calculated. The influence of wing tapering on the boundary layer development was shown. The turbulent boundary layer for a revolution ellipsoid is calculated for 0 deg and 10 deg incidence angles.

Stock, H. W.↗

Computational Algorithms for Unit Commitment with AC Power Flows (Final Report)

Security-constrained unit commitment (SCUC) is a key component in power system operations. When AC power flow constraints are considered in the SCUC model (AC-SCUC), the problem becomes extremely difficult due to its discrete and non-convex nature, as described in “Grid Optimization Competition Challenge 3 Problem Formulation (GOCC)”. There are four main challenges: (i) Discrete decisions regarding unit online/offline status and start-up/shut-down procedures for every single unit. The number of discrete decision variables increases considerably when a system integrates multiple generators; (ii) Configuration-based combined-cycle formulations, and multi-commodity models that include ramping products, spin/non-spin products, and regulation up/down products. The combined-cycle units introduce additional discrete decision variables and auxiliary service products further complicate the model by connecting multi-commodity products’ continuous and discrete variables; (iii) SCUC models with AC power flow constraints are far more complex due to massive bilinear terms in the large-scale nonlinear power balance equations. The nonlinear power balance equations are further complicated by the discrete step control variables of shunts; (iv) N − 1 contingency analysis. The size of the model increases linearly with the number of contingencies considered, greatly increasing the size of the optimization model. Accordingly, there is an emergent need to develop a robust algorithm capable of deriving a high-quality solution in a short time and passing through contingency tests simultaneously. In this project, we explore innovative techniques to address this challenging problem by integrating advanced polyhedral theory, approximation methods, relaxation strategies, decomposition techniques, and parallel computing. Each technique approaches the problem from a different perspective, leveraging its specific strengths to tackle distinct challenges. Each individual method has demonstrated its effectiveness in the PI’s previous research. Their integration is expected to significantly reduce the computational time required to solve the proposed complex problem. Successful completion of this project has the potential to transform the industry by enhancing optimization solvers capable of handling large-scale day-ahead energy market clearing models within strict time constraints, while incorporating AC power flow constraints. This advancement will lead to reduced overall generation costs and, consequently, increased social welfare.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Response of panels to turbulence-induced, surface-pressure fluctuations and resulting acoustic radiation to the flow field.

The response of finite rectangular panels, excited by turbulent boundary layers and the resulting acoustical radiation into the flow field is studied. The method of normal modes and the technique of spectral analysis are used to compute the power spectra of the panel amplitude and strain response. Computed results are compared with the corresponding response measurements for panels exposed to attached and separated flows at supersonic speeds and attached flow at subsonic speeds. The thin plate and the convected wave equations are coupled through the boundary condition on the panel to give a formula for the radiation pressure power spectra. A formula for the acoustical damping ratio of the panel is obtained through the principle of virtual work.-

Chyu, W. J.↗

The free compressible viscous vortex

The present study investigates the effects of compressibility on free (unsteady) viscous heat-conducting vortices. Analytical solutions are found in the limit of large but finite Reynolds number and small but finite Mach number. It is shown that the spreading of the vortex causes a radial flow. This flow is given by the solution of an ordinary differential equation, which gives the dependence of the radial velocity on the tangential velocity, density, and temperature profiles of the vortex. Estimates of the radial velocity found by solving this equation are found to be in good agreement with numerical solutions of the full equations. The equations for the viscous evolution are expanded in powers of Mach number to obtain detailed analytical solutions. It is shown that swirling axisymmetric compressible flows generate negative radial velocities far from the vortex core owing to viscous effects, regardless of the initial distributions of vorticity, density, and entropy.

Colonius, Tim↗

Optimization of suppression for two-element treatment liners for turbomachinery exhaust ducts

Sound wave propagation in a soft-walled rectangular duct with steady uniform flow was investigated at exhaust conditions, incorporating the solution equations for sound wave propagation in a rectangular duct with multiple longitudinal wall treatment segments. Modal analysis was employed to find the solution equations and to study the effectiveness of a uniform and of a two-sectional liner in attenuating sound power in a treated rectangular duct without flow (M = 0) and with uniform flow of Mach 0.3. Two-segment liners were shown to increase the attenuation of sound as compared to a uniform liner. The predicted sound attenuation was compared with measured laboratory results for an optimized two-segment suppressor. Good correlation was obtained between the measured and predicted suppressions when practical variations in the modal content and impedance were taken into account. Two parametric studies were also completed.

Motsinger, R. E.↗

Master State Distributed Estimator (masde)

The MSE technology provides a distributed state estimation scheme for power utilities by making it a real time dynamic system, versus the static snapshot method used today. This will allow utilities to verify power grid readings and identify false data on the communication network. The prevailing method of estimation takes all the data from the system into a single load flow equation at the utilities command center. This becomes a large algorithm that is very time consuming to solve, providing the engineers only snap shots of the system. The MFD distributes an algorithm to devises that are already installed at various locations on the power grid. This allows for much faster computer times because the algorithms contain magnitudes less data. When each distributed estimator completes its calculation, it sends the results back to the master device to check against all other distributed estimators.

Reen, DylanW.↗

Second-order small disturbance theory for hypersonic flow over power-law bodies

A mathematical method for determining the flow field about power-law bodies in hypersonic flow conditions is developed. The second-order solutions, which reflect the effects of the second-order terms in the equations, are obtained by applying the method of small perturbations in terms of body slenderness parameter to the zeroth-order solutions. The method is applied by writing each flow variable as the sum of a zeroth-order and a perturbation function, each multiplied by the axial variable raised to a power. The similarity solutions are developed for infinite Mach number. All results obtained are for no flow through the body surface (as a boundary condition), but the derivation indicates that small amounts of blowing or suction through the wall can be accommodated.

Townsend, J. C.↗