Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow conservation”

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 937 records · Page 52

Transmission-distribution long-term volt-var planning considering reactive power support capability of distributed PV

High penetration of grid-edge, inverter-based photovoltaic (PV) can cause significant voltage fluctuations not only at the distribution but also at the sub-transmission levels due to PV output intermittency. Traditional reactive power planning approaches do not consider intermittency, nor the possibility of coordinating the control of existing and future volt-ampere reactive resources. This paper proposes a reactive power planning tool for sub-transmission systems to mitigate voltage violations and fluctuations caused by high PV penetration and intermittency with a minimum investment cost. The planning tool coordinates with an optimization-based volt-var operational tool for: a) modeling the coordination of all existing var assets in both sub-transmission and distribution systems to reduce the need of new equipment, and b)selecting a set of scenarios with voltage violations, derived from PV intermittency c) testing the final investment decision. The tool obtains an investment need for each intermittency scenario with a proposed optimal power-flow framework with efficient techniques to handle a high number of discrete variables. Two options are provided for final planning decision: i) a conservative direct combination of investment need solutions and ii) a machine learning-based selection of representative investment needs at most time steps. The final investment decision options are verified using a realistic large-scale sub-transmission system and 5-minute PV and load data. The results show a significant voltage performance improvement with a lower investment cost for additional var equipment compared to conventional approaches.

14 SOLAR ENERGY↗

A two-dimensional numerical study of the flow inside the combustion chambers of a motored rotary engine

A numerical study was performed to investigate the unsteady, multidimensional flow inside the combustion chambers of an idealized, two-dimensional, rotary engine under motored conditions. The numerical study was based on the time-dependent, two-dimensional, density-weighted, ensemble-averaged conservation equations of mass, species, momentum, and total energy valid for two-component ideal gas mixtures. The ensemble-averaged conservation equations were closed by a K-epsilon model of turbulence. This K-epsilon model of turbulence was modified to account for some of the effects of compressibility, streamline curvature, low-Reynolds number, and preferential stress dissipation. Numerical solutions to the conservation equations were obtained by the highly efficient implicit-factored method of Beam and Warming. The grid system needed to obtain solutions were generated by an algebraic grid generation technique based on transfinite interpolation. Results of the numerical study are presented in graphical form illustrating the flow patterns during intake, compression, gaseous fuel injection, expansion, and exhaust.

Shih, T. I. P.↗

A two-dimensional numerical study of the flow inside the combustion chamber of a motored rotary engine

A numerical study was performed to investigate the unsteady, multidimensional flow inside the combustion chambers of an idealized, two-dimensional, rotary engine under motored conditions. The numerical study was based on the time-dependent, two-dimensional, density-weighted, ensemble-averaged conservation equations of mass, species, momentum, and total energy valid for two-component ideal gas mixtures. The ensemble-averaged conservation equations were closed by a K-epsilon model of turbulence. This K-epsilon model of turbulence was modified to account for some of the effects of compressibility, streamline curvature, low-Reynolds number, and preferential stress dissipation. Numerical solutions to the conservation equations were obtained by the highly efficient implicit-factored method of Beam and Warming. The grid system needed to obtain solutions were generated by an algebraic grid generation technique based on transfinite interpolation. Results of the numerical study are presented in graphical form illustrating the flow patterns during intake, compression, gaseous fuel injection, expansion, and exhaust.

Shih, T. I-P.↗

A Moving Embedded Boundary Approach for the Compressible Navier-Stokes Equations in a Block-Structured Adaptive Refinement Framework

A computational technique has been developed to perform compressible flow simulations involving moving boundaries using an embedded boundary approach within the block-structured adaptive mesh refinement (SAMR) framework of AMReX [1], [91], [92]. We leverage the SAMR capability to obtain quantitatively accurate results whilst using robust, second-order finite volume schemes. A conservative, unsplit, cut-cell approach is utilized and a ghost-cell approach is developed for computing the flux on the moving, embedded boundary faces. A third-order least-squares formulation has been developed to compute the wall velocity gradients, and was found to significantly improve the performance of the solver in terms of the quantitative comparison of surface quantities such as the skin friction coefficient. Various test cases are performed to validate the method, and compared with analytical, experimental, and other numerical results in literature. Inviscid and viscous test cases are performed that span a wide regime of flow speeds - acoustic (harmonically pulsating sphere), smooth flows (expansion fan created by a receding piston) and flows with shocks (shock-cylinder interaction, shock-wedge interaction, pitching NACA 0012 airfoil and shock-cone interaction). A closed system with moving boundaries - an oscillating piston in a cylinder, showed that the percentage error in mass within the system decreases with refinement, demonstrating that the numerical scheme is conservative with grid refinement, but is not discretely conservative. Viscous test cases involve that of a horizontally moving cylinder at Re = 40, an inline oscillating cylinder at Re = 100, and a transversely oscillating cylinder at Re = 185. The judicious use of adaptive mesh refinement with appropriate refinement criteria to capture the regions of interest leads to well-resolved flow features, and good quantitative comparison is observed with the results available in literature.

adaptive refinement↗

In tube condensation of low global warming potential refrigerants in an axial micro-fin aluminum tube

Environmental regulations have driven the development of refrigerants with low global warming potential (GWP). To design heat exchangers using these new refrigerants, data are needed concerning the heat transfer coefficient and pressure drop in two-phase flow. Another change is the increasing use of aluminum tubes rather than copper tubes to reduce heat exchanger cost. Hence, this study presents an experimental investigation of flow condensation using an expanded axial micro-fin aluminum tube with a fin-tip diameter of 5.96 mm. Here, the experiments included single compounds R-32, R-1234yf, and R-1234ze(E), zeotropic mixtures with low glide (R-454B), and zeotropic mixtures with high-glide (R-454C and R-455A). Experiments were conducted at condensation temperatures ranging from 40 °C to 50 °C, reduced pressures ranging from 0.21 to 0.55, and mass fluxes ranging from 150 to 350 kg/(m 2 s). Data obtained for these refrigerants constitute one of the first reports for high-glide refrigerants using axial micro-fin aluminum tubes. An evaluation of heat transfer degradation of zeotropic mixtures due to mass transfer resistance at the liquid/vapor interface is presented. This information can be used to design heat exchangers for next generation air conditioning and refrigeration systems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Analytical solution for partial heating on the exterior of the pipe with application to measuring fluid flow rate

An analytical model using Green’s functions for partial external heating of a pipe is developed, which results in an exact mathematical solution for the radial and axial temperature distribution in the pipe wall. Partial heating consists of a constant heat flux function imposed over a small section of the exterior of a pipe and for a limited time duration. The solution comprises steady-state and transient parts, and an algebraic identity is used to decrease the number of summation terms in the slowly-converging steady- state part. Intrinsic verification principles are used to verify the solution. As an example application, this transient solution is applied toward the development of a simple, noninvasive method for in-field measurement of the flow rate in pipes. Furthermore, to simulate this application, a pulse of energy is imposed to the wall of the pipe, and the developed mathematical solution is used to find the flow rate of the fluid inside the pipe. An optimal experiment is designed to find the best measurement location and time.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Characterization of the Effect of Wing Surface Instrumentation on UAV Airfoil Performance

Recently proposed flight research at NASA Dryden Flight Research Center (DFRC) has prompted study into the aerodynamic effects of modifications made to the surfaces of laminar airfoils. The research is focused on the high-aspect ratio, laminar-flow type wings commonly found on UAVs and other aircraft with a high endurance requirement. A broad range of instrumentation possibilities, such as structural, pressure, and temperature sensing devices may require the alteration of the airfoil outer mold line as part of the installation process. This study attempts to characterize the effect of installing this additiona1 instrumentation on key airfoil performance factors, such as transition location, lift and drag curves, and stall point. In particular, the general case of an airfoil that is channeled in the spanwise direction is considered, and the impact on key performance characteristics is assessed. Particular attention is focused on exploring the limits of channel depth and low-Reynolds number on performance and stall characteristics. To quantify the effect of increased skin friction due to premature transition caused by protruding or recessed instrumentation, two simplified, conservative scenarios are used to consider two potential sources of diaturbance: A) that leading edge alterations would cause linearly expanding areas (triangles) of turbulent flow on both surfaces of the wing upstream of the natural transition point, and B) that a channel or bump on the upper surface would trip turbulent flow across the whole upper surface upstream of the natural transition point. A potentially more important consideration than the skin friction drag increment is the change in overall airfoil performance due to the installation of instrumentation along most of the wingspan. To quantify this effect, 2D CFD simulations of the flow over a representative mid-span airfoil section were conducted in order to assess the change in lift and drag curves for the airfoil in the presence of disturbed flow due to the installed instrumentation. A discussion as to the impact on high-altitude and low-speed operation of this and similar aircraft is provided.

Ratnayake, Nalin A.↗

A modified Dodge algorithm for the parabolized Navier-Stokes equation and compressible duct flows

A revised version of Dodge's split-velocity method for numerical calculation of compressible duct flow was developed. The revision incorporates balancing of mass flow rates on each marching step in order to maintain front-to-back continuity during the calculation. The (checkerboard) zebra algorithm is applied to solution of the three dimensional continuity equation in conservative form. A second-order A-stable linear multistep method is employed in effecting a marching solution of the parabolized momentum equations. A checkerboard iteration is used to solve the resulting implicit nonlinear systems of finite-difference equations which govern stepwise transition. Qualitive agreement with analytical predictions and experimental results was obtained for some flows with well-known solutions.

Cooke, C. H.↗

A modified Dodge algorithm for the parabolized Navier-Stokes equations and compressible duct flows

A revised version of Dodge's split-velocity method for numerical calculation of compressible duct flow was developed. The revision incorporates balancing of mass flow rates on each marching step in order to maintain front-to-back continuity during the calculation. The (checkerboard) zebra algorithm is applied to solution of the three dimensional continuity equation in conservative form. A second-order A-stable linear multistep method is employed in effecting a marching solution of the parabolized momentum equations. A checkerboard iteration is used to solve the resulting implicit nonlinear systems of finite-difference equations which govern stepwise transition. Qualitative agreement with analytical predictions and experimental results was obtained for some flows with well-known solutions. Previously announced in STAR as N82-16363

Cooke, C. H.↗

Fuel-rich oxidation of gasoline surrogate components in an atmospheric flow reactor

Fuel-rich oxidation of three typical gasoline surrogate components, toluene, isooctane, and n-heptane, was investigated in an atmospheric-pressure flow reactor at mean gas temperatures from 1050 to 1350 K, equivalence ratio of 9.0, and residence times of 0.45 and 1.2 s. Not only polycyclic aromatic hydrocarbons (PAHs) up to 3 ring structure but also small intermediate products from C 1 to C 5 were quantified by a gas chromatograph mass spectrometry coupled with photon ionization and gas chromatograph with flame ionization detector, respectively. The kinetic model recently developed by Lawrence Livermore National Laboratory was revised to reflect the results of many recent investigations. Basically, the updated model could satisfactorily reproduce the experimental mole fractions of many species. The experimental and simulated results showed that PAH mole fractions produced were in the order of toluene, isooctane, and n-heptane. The kinetic analysis using the model was carried out to explore PAH formation pathways, especially focusing on naphthalene, acenaphthylene, and phenanthrene. Through rate of production analysis, it was found that the main formation pathways of many PAHs were affected by the fuels. Although resonantly stabilized radicals, such as benzyl and fulvenallenyl radicals, played a crucial role in the formation pathways of many PAHs in every fuel, they were produced through hydrogen elimination of toluene in toluene fuel, while they were formed from small intermediate products in isooctane and n-heptane fuels. Sensitivity analysis revealed the difference and similarity of the reactions with large positive coefficients according to the fuels studied here. In conclusion, the molecular growth reactions of aromatic species were influential in PAH production in every fuel, whereas the ring formation reactions from small species and the reactions involving toluene had large positive sensitivity coefficients in isooctane/n-heptane and toluene fuels, respectively.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Transport Analysis & Optimization in a MW-Scale CO2 Electrolyzer (Final Report)

As Twelve continues to scale up their CO2 electrolyzers, both in the size of a single cell and in the number of cells used in a stack, thermal management becomes a growing concern, since excess heat can affect reaction yield and accelerate degradation. In this project, we aim to computationally explore how the anode flow fields used in Twelve’s CO2 electrolyzers function as heat exchangers. In particular, using a homogenized model of a CO2 electrolyzer, we first estimate the amount of heat generated in a cell. Then, we develop a computational fluid dynamics (CFD) model of the so-called “flow field”, i.e. a flow manifold, based on Twelve’s CAD drawings, to evaluate how these flow fields perform as a heat exchanger for the generated heat. We explore both a single cell and a 3-cell stack operating in parallel, where heat generated in one cell can now be transferred to another cell. We evaluate how performance is affected when environmental heat losses are taken into account. Finally, we leverage topology optimization to explore the types of design features a computational optimization algorithm would suggest to supplement our intuition. Overall, our work aims to provide design recommendations for CO2 electrolyzer flow fields and provides a foundation for future studies of flow field optimization.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Numerical solutions of reactive fluid flows during postignition transients in hybrid rocket systems.

A computational method has been developed for the study of the post-ignition transients in hybrid rocket systems. The particular system chosen consisted of a gaseous oxidizer flowing within a tube of solid fuel, resulting in heterogeneous combustion. With the appropriate assumptions, two-dimensional, time-dependent conservation equations were derived for the reacting gas phase, and for the solid phase, in a cylindrical coordinate system. These were then programmed for numerical computation, using two implicit finite-difference schemes, the Lax-Wendroff scheme for the gas phase, and the Crank-Nicolson scheme for the solid phase. Appropriate initial and boundary conditions were represented, including heat and mass conservation at the interface between gas and solid. Initially, no attempt was made to relate the recession rate at the surface to the surface temperature, or to include heat transfer by radiation. A simple case was selected for preliminary calculations, with aluminum and oxygen as fuel and oxidizer, and aluminum oxide as the product.

Hung, W. S. Y.↗

Tone noise of three supersonic helical tip speed propellers in a wind tunnel

Three eight bladed supersonic helical tip speed propellers were tested in the NASA Lewis Wind Tunnel at through flow Mach numbers of 0.6, 0.7, 0.75, 0.8, and 0.85. Such propellers are being considered as propulsors for a future energy conservative aircraft. The three propellers of equal performance, incorporate different plan forms and different amounts of sweep in order to investigate their noise generation, which could cause a cabin environment problem. The acoustically designed propeller SR-3 had 45 deg of tip sweep and was significantly quieter at Mach 0.8 cruise than the straight bladed propeller, SR-2. An intermediate 30 deg tip sweep propeller, SR-1m exhibited a noise level between the other two. Enhanced pressure-time traces indicated that SR-2 and SR-1m exhibited shock-like pressure traces at Mach-0.8 cruise while the SR-3 did not. Noise trends with varying helical tip Mach number and blade loading were also observed.

Dittmar, J. H.↗

Computer program for calculating full potential transonic, quasi-three-dimensional flow through a rotating turbomachinery blade row

A fast, reliable computer code is described for calculating the flow field about a cascade of arbitrary two dimensional airfoils. The method approximates the three dimensional flow in a turbomachinery blade row by correcting for stream tube convergence and radius change in the throughflow direction. A fully conservative solution of the full potential equation is combined with the finite volume technique on a body-fitted periodic mesh, with an artificial density imposed in the transonic region to insure stability and the capture of shock waves. The instructions required to set up and use the code are included. The name of the code is QSONIC. A numerical example is also given to illustrate the output of the program.

Farrell, C. A.↗

Analysis of pultrusion processing for long fiber reinforced thermoplastic composite system

Pultrusion is one of the composite processing technology, commonly recognized as a simple and cost-effective means for the manufacturing of fiber-reinforced, resin matrix composite parts with different regular geometries. Previously, because the majority of the pultruded composite parts were made of thermosetting resin matrix, emphasis of the analysis on the process has been on the conservation of energy from various sources, such as heat conduction and the curing kinetics of the resin system. Analysis on the flow aspect of the process was almost absent in the literature for thermosetting process. With the increasing uses of thermoplastic materials, it is desirable to obtain the detailed velocity and pressure profiles inside the pultrusion die. Using a modified Darcy's law for flow through porous media, closed form analytical solutions for the velocity and pressure distributions inside the pultrusion die are obtained for the first time. This enables us to estimate the magnitude of viscous dissipation and it's effects on the pultruded parts. Pulling forces refined in the pultrusion processing are also analyzed. The analytical model derived in this study can be used to advance our knowledge and control of the pultrusion process for fiber reinforced thermoplastic composite parts.

Tso, W.↗

Features and applications of the Groove Analysis Program (GAP)

An IBM Personal Computer (PC) version of the Groove Analysis program (GAP) was developed to predict the steady state heat transport capability of an axially grooved heat pipe for a specified groove geometry and working fluid. In the model, the capillary limit is determined by the numerical solution of the differential equation for momentum conservation with the appropriate boundary conditions. This governing equation accounts for the hydrodynamic losses due to friction in liquid and vapor flows and due to liquid/vapor shear interaction. Back-pumping in both 0-g and 1-g is accounted for in the boundary condition at the condenser end. Slug formation in 0-g and puddle flow in 1-g are also considered in the model. At the user's discretion, the code will perform the analysis for various fluid inventories (undercharge, nominal charge, overcharge, or a fixed fluid charge) and heat pipe elevations. GAP will also calculate the minimum required heat pipe wall thickness for pressure containment at design temperatures that are greater than or lower than the critical temperature of the working fluid. This paper discusses the theory behind the development of the GAP model. It also presents the many useful and powerful capabilities of the model. Furthermore, a correlation of flight test performance data and the predictions using GAP are presented and discussed.

Ku, Jentung↗

Numerical Modeling of Cavitating Venturi: A Flow Control Element of Propulsion System

In a propulsion system, the propellant flow and mixture ratio could be controlled either by variable area flow control valves or by passive flow control elements such as cavitating venturies. Cavitating venturies maintain constant propellant flowrate for fixed inlet conditions (pressure and temperature) and wide range of outlet pressures, thereby maintain constant, engine thrust and mixture ratio. The flowrate through the venturi reaches a constant value and becomes independent of outlet pressure when the pressure at throat becomes equal to vapor pressure. In order to develop a numerical model of propulsion system, it is necessary to model cavitating venturies in propellant feed systems. This paper presents a finite volume model of flow network of a cavitating venturi. The venturi was discretized into a number of control volumes and mass, momentum and energy conservation equations in each control volume are simultaneously solved to calculate one-dimensional pressure, density, and flowrate and temperature distribution. The numerical model predicts cavitations at the throat when outlet pressure was gradually reduced. Once cavitation starts, with further reduction of downstream pressure, no change in flowrate is found. The numerical predictions have been compared with test data and empirical equation based on Bernoulli's equation.

Majumdar, Alok↗

Freshwater Requirements of Large-Scale Bioenergy Plantations for Limiting Global Warming to 1.5°C

Limiting mean global warming to well below 2 °C will probably require substantial negative emissions (NEs) within the 21st century. To achieve these, bioenergy plantations with subsequent carbon capture and storage (BECCS) may have to be implemented at a large scale. Irrigation of these plantations might be necessary to increase the yield, which is likely to put further pressure on already stressed freshwater systems. Conversely, the potential of bioenergy plantations (BPs) dedicated to achieving NEs through CO2 assimilation may be limited in regions with low freshwater availability. This paper provides a first-order quantification of the biophysical potentials of BECCS as a negative emission technology contribution to reaching the 1.5 °C warming target, as constrained by associated water availabilities and requirements. Using a global biosphere model, we analyze the availability of freshwater for irrigation of BPs designed to meet the projected NEs to fulfill the 1.5 °C target, spatially explicitly on areas not reserved for ecosystem conservation or agriculture. We take account of the simultaneous water demands for agriculture, industries, and households and also account for environmental flow requirements (EFRs) needed to safeguard aquatic ecosystems. Furthermore, we assess to what extent different forms of improved water management on the suggested BPs and on cropland may help to reduce the freshwater abstractions. Results indicate that global water withdrawals for irrigation of BPs range between ~400 and ~3000 km(exp 3) yr(exp -1), depending on the scenario and the conversion efficiency of the carbon capture and storage process. Consideration of EFRs reduces the NE potential significantly, but can partly be compensated for by improved on-field water management.

Stenzel, Fabian↗