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 361 records · Page 20

Scaling Analysis of Two–Phase Flow in Fractal Permeability Fields

Fluid mixing in permeable media is essential in many practical applications. The mixing process is a consequence of velocity fluctuations owing to geological heterogeneities and mobility contrast of fluids. Heterogeneities in natural rocks are often spatially correlated, and their properties, such as permeability, may be described using fractal distributions. This work models the fractal characteristics of such permeability fields in which the covariance function is expressed as a power-law function. A generalized scaling relation is derived relating various fractal permeability fields using the magnitude of their fluctuations. Here, this relation reveals the self-similar behavior of two-phase flow in such permeable media. To that end, a recently developed, high-resolution numerical simulator is employed to validate the analytically derived scaling relations. Two flow problems are considered in which flow is governed by 1) a linear, and 2) a nonlinear transport equation. Due to the probabilistic representation of the fractal permeability fields, a sensitivity study is conducted for each flow scenario to determine the number of realizations required for statistical convergence. Scaling analysis is performed using ensemble averages of simulated saturation profiles and their mixing lengths. Results support the validity of the developed scaling relation across the range of investigated flow conditions at intermediate times. The dynamics of linear flow in the asymptotic regime is affected by the correlation structure of heterogeneity. In nonlinear flow, scaling behavior appears to be dominated by the degree of nonlinearity.

58 GEOSCIENCES↗

Diagnostics of seeded RF plasmas: An experimental study related to the gaseous core reactor

Measurements of the temperature profiles in an RF argon plasma were made over magnetic field intensities ranging from 20 amp turns/cm to 80 amp turns/cm. The results were compared with a one-dimensional numerical treatment of the governing equations and with an approximate closed form analytical solution that neglected radiation losses. The average measured temperatures in the plasma compared well with the numerical treatment, though the experimental profile showed less of an off center temperature peak than predicted by theory. This may be a result of the complex turbulent flow pattern present in the experimental torch and not modeled in the numerical treatment. The radiation term cannot be neglected for argon at the power levels investigated. The closed form analytical approximation that neglected radiation led to temperature predictions on the order of 1000 K to 2000 K higher than measured or predicted by the numerical treatment which considered radiation losses.

Thompson, S. D.↗

Euler analysis of the three dimensional flow field of a high-speed propeller: Boundary condition effects

The results of an investigation of the effects of far field boundary conditions on the solution of the three dimensional Euler equations governing the flow field of a high speed single rotation propeller are presented. The results show that the solutions obtained with the nonreflecting boundary conditions are in good agreement with experimental data. The specification of nonreflecting boundary conditions is effective in reducing the dependence of the solution on the location of the far field boundary. Details of the flow field within the blade passage and the tip vortex are presented. The dependence of the computed power coefficient on the blade passage and the tip vortex are presented. The dependence of the computed power coefficient on the blade setting angle is examined.

Nallasamy, M.↗

Application of turbulence modeling to predict surface heat transfer in stagnation flow region of circular cylinder

A theoretical analysis and numerical calculations for the turbulent flow field and for the effect of free-stream turbulence on the surface heat transfer rate of a stagnation flow are presented. The emphasis is on the modeling of turbulence and its augmentation of surface heat transfer rate. The flow field considered is the region near the forward stagnation point of a circular cylinder in a uniform turbulent mean flow. The free stream is steady and incompressible with a Reynolds number of the order of 10 to the 5th power and turbulence intensity of less than 5 percent. For this analysis, the flow field is divided into three regions: (1) a uniform free-stream region where the turbulence is homogeneous and isotropic; (2) an external viscid flow region where the turbulence is distorted by the variation of the mean flow velocity; and, (3) an anisotropic turbulent boundary layer region over the cylinder surface. The turbulence modeling techniques used are the kappa-epsilon two-equation model in the external flow region and the time-averaged turbulence transport equation in the boundary layer region. The turbulence double correlations, the mean velocity, and the mean temperature within the boundary layer are solved numerically from the transport equations. The surface heat transfer rate is calculated as functions of the free-stream turbulence longitudinal microlength scale, the turbulence intensity, and the Reynolds number.

Wang, Chi R.↗

Euler analysis of the three-dimensional flow field of a high-speed propeller - Boundary condition effects

The results of an investigation of the effects of far field boundary conditions on the solution of the three dimensional Euler equations governing the flow field of a high speed single rotation propeller are presented. The results show that the solutions obtained with the nonreflecting boundary conditions are in good agreement with experimental data. The specification of nonreflecting boundary conditions is effective in reducing the dependence of the solution on the location of the far field boundary. Details of the flow field within the blade passage and the tip vortex are presented. The dependence of the computed power coefficient on the blade passage and the tip vortex are presented. The dependence of the computed power coefficient on the blade setting angle is examined.

Nallasamy, M.↗

Acoustic Mode Decomposition in Rectangular Ducts with Sheared Flow

The performance of new acoustic liner concepts are, in general, characterized and assessed in grazing flow rigs early in the development cycle. These test rigs expose an acoustic liner sample, installed on the side wall of the duct, to a grazing flow and incident acoustic field. The process to characterize these liners involves educing the impedance on the wall where the sample is installed and examining the acoustic power attenuation. Standard approaches to computing impedance or power attenuation generally consider only the effects of a 2D shear flow or uniform flow, on the acoustic field. In this study, the objective is to incorporate 3D shear flow effects in the analysis of acoustic mode attenuation in a rectangular duct flow rig. A modal analysis of microphone measurements obtained on the side walls of the duct upstream and downstream of the test section of the rig is developed. A Galerkin projection of the Pridmore-Brown equation is performed with Chebyshev basis functions in order to incorporate the effects of the Mach number profile on the computation of the axial wavenumber of each mode. Measurements of the Mach number profile are obtained in the test rig and used as input to compute the modes. Comparisons made between the sound field computed with traditional convective Helmholtz modes and the new procedure using Pridmore-Brown modes indicate that the computed acoustic field using Pridmore-Brown modes more accurately reconstructs the acoustic signal at each microphone in the array. The mode structure of the lowest-order mode is shown to be significantly impacted by shear flow refraction effects, and higher-order mode structures are also affected at higher frequencies and centerline Mach number. An assessment of the acoustic mode attenuation for two acoustic liner samples demonstrates that the computed mode amplitudes for both the traditional and new approach are in agreement for the lowest-order mode, but discrepancies arise when higher-order modes are the dominant component of the acoustic field.

Duct acoustics↗

Acoustic Mode Decomposition in Rectangular Ducts with Sheared Flow

The performance of new acoustic liner concepts are, in general, characterized and assessed in grazing flow rigs early in the development cycle. These test rigs expose an acoustic liner sample, installed on the side wall of the duct, to a grazing flow and incident acoustic field. The process to characterize these liners involves educing the impedance on the wall where the sample is installed and examining the acoustic power attenuation. Standard approaches to computing impedance or power attenuation generally consider only the effects of a 2D shear flow or uniform flow, on the acoustic field. In this study, the objective is to incorporate 3D shear flow effects in the analysis of acoustic mode attenuation in a rectangular duct flow rig. A modal analysis of microphone measurements obtained on the side walls of the duct upstream and downstream of the test section of the rig is developed. A Galerkin projection of the Pridmore-Brown equation is performed with Chebyshev basis functions in order to incorporate the effects of the Mach number profile on the computation of the axial wavenumber of each mode. Measurements of the Mach number profile are obtained in the test rig and used as input to compute the modes. Comparisons made between the sound field computed with traditional convective Helmholtz modes and the new procedure using Pridmore-Brown modes indicate that the computed acoustic field using Pridmore-Brown modes more accurately reconstruct the acoustic signals at each microphone in the array. The mode structure of the lowest order mode is shown to be significantly impacted by shear flow refraction effects, and higher order mode structures are also effected at higher frequencies and centerline Mach number. An assessment of the acoustic mode attenuation for two acoustic liner samples demonstrates that the computed mode amplitudes for both the traditional and new approach are in agreement for the lowest order mode, but discrepancies arise when higher order modes are the dominant component of the acoustic field.

Duct acoustics↗

Acoustic scattering from locally homogeneous turbulence

An equation is derived for acoustic scattering from the velocity and temperature fields of a drifting blob of air turbulence. The equation is solved for scattering from the mean and turbulent portions of both fields, and the time autocorrelation and power spectral density of the received signal are calculated. The mean velocity and temperature fields, as well as the corresponding turbulence intensity distributions, are allowed to be spatially nonuniform. The turbulent fields are required to be only locally homogeneous and locally stationary. It is shown that the common practice of representing antenna patterns by either truncating the flow or tapering the strength of the flow is, at least in the case of the transmitting pattern, not a bad approximation. Spectral broadening of the receiver signal due to convection of the small scattering eddies by macroeddies and by the mean flow in modeled, and the broadening by macroeddy convection is seen to render negligible the broadening which arises from the drift of the target flow through the antenna beam. The analysis reveals that only for target flows having a high degree of spatial uniformity and/or symmetry does the received positive-frequency power spectral density turn out to be symmetrical about a center frequency.

Ramsey, V. W.↗

Wake shield

Progress on the modeling of the flow field around a wake shield using a recently obtained code based on the Monte Carlo method is discussed. The direct simulation Monte Carlo method is a method for solving the Boltzman Equation using an approximation to the collision integral term. The collision integrand is evaluated for randomly selected values of its arguments and the summation will approach the integral for large enough samples. The collision effects may be modeled for either hard sphere or various power law potentials. The convective side of the Boltzman equation is approximated over a time step using a simple trajectory calculation of molecules as they travel through the domain of interest.

Bannister, Tommy↗

Adaptive Suction and Blowing for Twin-Tail Buffet Control

Adaptive active flow control for twin-tail buffet alleviation is investigated. The concept behind this technique is to place control ports on the tail outer and inner surfaces with flow suction or blowing applied through these ports in order to minimize the pressure difference across the tail. The suction or blowing volume flow rate from each port is proportional to the pressure difference across the tail at this location. A parametric study of the effects of the number and location of these ports on the buffet response is carried out. The computational model consists of a sharp-edged delta wing of aspect ratio one and swept-back flexible twin tail with taper ratio of 0.23. This complex multidisciplinary problem is solved sequentially using three sets of equations for the fluid flow, aeroelastic response and grid deformation, using a dynamic multi-block grid structure. The computational model is pitched at 30 deg angle of attack. The freestream Mach number and Reynolds number are 0.3 and 1.25 million, respectively. The model is investigated for the inboard position of the twin tails, which corresponds to a separation distance between the twin tails of 33% of the wing span. Comparison of the time history and power spectral density responses of the tails for various distributions of the control ports are presented and discussed.

Kandil, Osama A.↗

Water Vapor Effects on Silica-Forming Ceramics

Silica-forming ceramics such as SiC and Si3N4 are proposed for applications in combustion environments. These environments contain water vapor as a product of combustion. Oxidation of silica-formers is more rapid in water vapor than in oxygen. Parabolic oxidation rates increase with the water vapor partial pressure with a power law exponent value close to one. Molecular water vapor is therefore the mobile species in silica. Rapid oxidation rates and large amounts of gases generated during the oxidation reaction in high water vapor pressures may result in bubble formation in the silica and nonprotective scale formation. It is also shown that silica reacts with water vapor to form Si(OH)4(g). Silica volatility has been modeled using a laminar flow boundary layer controlled reaction equation. Silica volatility depends on the partial pressure of water vapor, the total pressure, and the gas velocity. Simultaneous oxidation and volatilization reactions have been modeled with paralinear kinetics.

Opila, E. J.↗

Molecular flow venting of a volume with an outgassing or desorbing source

A venting equation commonly used to describe the transient evolution of pressure within a volume containing an outgassing or offgassing source and a restrictive vent conductance under conditions of molecular flow has been solved analytically. Solutions are found for sources of finite thickness characterized by classical diffusion-limited behavior (proportional to inverse square root of time), as well as responses for thick material sources often observed in testing that are characterized by a more general form of power-law decay, up to inverse time behavior associated with surface desorption. Solutions involve evaluating integrals where both numerators and denominators of the integrands diverge with time, making wide-ranging transient solutions difficult to directly compute numerically. Usually, one can avoid evaluating these integrals by assuming quasistatic conditions at long time scales. A novel approach is used in this work to analytically produce solutions by generating bespoke mathematical functions, some of which solve integrals that have apparently had no previous analytical solution.

Michael Woronowicz↗

Molecular flow venting of a volume with an outgassing or desorbing source

"A venting equation commonly used to describe the transient evolution of pressure within a volume containing an outgassing or offgassing source and a restrictive vent conductance under conditions of molecular flow has been solved analytically. Solutions are found for sources of finite thickness characterized by classical diffusion-limited behavior (proportional to inverse square root of time), as well as responses for thick material sources often observed in testing that are characterized by a more general form of power-law decay, up to inverse time behavior associated with surface desorption. Solutions involve evaluating integrals where both numerators and denominators of the integrands diverge with time, making wide-ranging transient solutions difficult to directly compute numerically. Usually, one can avoid evaluating these integrals by assuming quasistatic conditions at long time scales. A novel approach is used in this work to analytically produce solutions by generating bespoke mathematical functions, some of which solve integrals that have apparently had no previous analytical solution."

Michael Woronowicz↗

Perovskite Sorbent Oxygen Separation Modeling with MFiX

This document chronicles the development and implementation of computational kinetic rate models that capture absorption and desorption characteristics of the National Energy Technology Laboratory (NETL) developed perovskite, Sr 1-x Ca x FeO 3-δ . Two paths to create accurate kinetic rates were followed: (1) an isothermal rate approach where thermogravimetric (TGA) data are recast as oxygen capacities through a pseudo-second order Lagergren equation (He et al., 2009); and (2) a more traditional Arrhenius approach where experimental data are fit with a power law model to derive associate activation energies (Bulfin et al., 2020a). For reference, the mathematics and associate experimental strategies that support these derivations are included in this report. In addition, computational fluid dynamics (CFD) models were developed to utilize both kinetic rate formulations and applied to simulate oxygen uptake and release in small scale scenarios. The program Multiphase Flow with interphase eXchanges (MFiX) was used to create: (1) discrete element method (DEM) simulations of a single tube of granular perovskite experiencing isothermal O 2 -absorption; and desorption and (2) two-fluid-model (TFM) non-isothermal simulations of perovskite O 2 -absorption and desorption tubes that share a wall. Conjugate heat transfer between steel walled tubes and the perovskite bed were managed with user-defined functions. As the project moves to simulating larger scale devices that will require more robust conjugate heat transfer methods, developed kinetic rates and associate methodologies have been recast for use in the ANSYS Fluent CFD program.

36 MATERIALS SCIENCE↗

Study of one-dimensional spectral dynamic equations of the Reynolds stresses in homogeneous anisotropic turbulence: Application to split-spectrum modeling

The CTR numerical data base generated by direct simulation of homogeneous anisotropic turbulence was used to calculate all of the terms in the spectral balance equations for the turbulent Reynolds stresses. The aim in not only to test the main closure assumptions used in the split-spectrum models, but also to try to devise improved hypotheses deduced from the statistical information. Numerical simulations of turbulent flows provide a large amount of data, a thought provoking wealth of information. The main advantage of this type of comparison is that a great variety of flows can be considered, and this is necessary to test closure hypotheses. Moreover various initial conditions can be introduced in the calculation, even if they are not experimentally feasible. All the terms in the spectral equations can be calculated. The limited Reynolds numbers of the simulations and the statistical noise caused by a small sample, particularly at the large scales, causes some difficulty in the interpretation of the results, but the method of approach proved to be a powerful tool for testing and improving spectral closures.

Schiestel, R.↗

QPOs from Random X-ray Bursts around Rotating Black Holes

We continue our earlier studies of quasi-periodic oscillations (QPOs) in the power spectra of accreting, rapidly-rotating black holes that originate from the geometric 'light echoes' of X-ray flares occurring within the black hole ergosphere. Our present work extends our previous treatment to three-dimensional photon emission and orbits to allow for arbitrary latitudes in the positions of the distant observers and the X-ray sources in place of the mainly equatorial positions and photon orbits of the earlier consideration. Following the trajectories of a large number of photons we calculate the response functions of a given geometry and use them to produce model light curves which we subsequently analyze to compute their power spectra and autocorrelation functions. In the case of an optically-thin environment, relevant to advection-dominated accretion flows, we consistently find QPOs at frequencies of order of approximately kHz for stellar-mass black hole candidates while order of approximately mHz for typical active galactic nuclei (approximately equal to 10(exp 7) solar mass) for a wide range of viewing angles (30 degrees to 80 degrees) from X-ray sources predominantly concentrated toward the equator within the ergosphere. As in out previous treatment, here too, the QPO signal is produced by the frame-dragging of the photons by the rapidly-rotating black hole, which results in photon 'bunches' separated by constant time-lags, the result of multiple photon orbits around the hole. Our model predicts for various source/observer configurations the robust presence of a new class of QPOs, which is inevitably generic to curved spacetime structure in rotating black hole systems.

Kukumura, Keigo↗

Nonlinear evolution of a disturbance in an unbounded viscous fluid with uniform shear

The evolution of a disturbance in the presence of a uniform mean velocity gradient is calculated by a power-series solution of the incompressible Navier-Stokes equations. Terms through those in time cubed are retained in the solution. For the initial condition a three-dimensional cosine distribution with two harmonic terms is assumed. The nonlinear interaction of these harmonic terms produces new harmonics which in turn interact. For large velocity gradients the energy of the disturbance can grow with time. The results shed some light on the maintenance or growth of turbulence in a shear flow.

Deissler, R. G.↗

Flight test evaluation of predicted light aircraft drag, performance, and stability

A technique was developed which permits simultaneous extraction of complete lift, drag, and thrust power curves from time histories of a single aircraft maneuver such as a pull up (from V max to V stall) and pushover (to V max for level flight). The technique, which is an extension of nonlinear equations of motion of the parameter identification methods of Iliff and Taylor and includes provisions for internal data compatibility improvement as well, was shown to be capable of correcting random errors in the most sensitive data channel and yielding highly accurate results. Flow charts, listings, sample inputs and outputs for the relevant routines are provided as appendices. This technique was applied to flight data taken on the ATLIT aircraft. Lack of adequate knowledge of the correct full throttle thrust horsepower true airspeed variation and considerable internal data inconsistency made it impossible to apply the trajectory matching features of the technique.

Smetana, F. O.↗