Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable 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 1,081 records · Page 60

Aerodynamic effects

The tethered satellite concept provides an ideal platform for the study of the interaction of the atmosphere with satellites of various shapes and surfaces under a wide range of flow conditions. From experiments which would measure the drag, lift, and torque acting on the tethered satellite, important information could be obtained which would have application to satellite lifetime prediction, determination of properties of the upper atmosphere, and scientific information on the interaction of high speed molecules with surfaces (the gas surface interaction). These experiments using the tethered satellite concept are described and would measure the following variables: angle of attack, surface roughness, and flow properties.

Karr, G. R.↗

The Onset of Nucleate Boiling Prediction in a Mini Channel with the Eulerian Multiphase Flow

The prediction of boiling incipience is a critical issue for a reactor. Vapor lock causes operating instability and undesirable temperature rise, especially for small channels. This study uses the Eulerian multiphase flow boiling model to predict the onset of nucleate boiling in a millimeter scale rectangular channel. Two wall boiling models, namely the RPI (Rensselaer Polytechnic Institute) and the non-equilibrium sub cooled model are applied. Computational fluid dynamics is employed for this study to investigate the impact of operating conditions like liquid flow rate, operating pressure, heat flux, velocity profile, turbulence model, and inlet liquid temperature. Here, the effect of fundamental boiling parameters on the surface wall temperature and vapor volume fraction is also studied. The parameters under consideration are the bubble departure diameter, bubble departure frequency, nucleate site density, quenching time period, and interface heat transfer coefficients. A series of CFD calculations is conducted by varying the considered variables systematically in a wide range of flow conditions covering laminar, transition, and turbulent flows. The sub-component heat fluxes like the liquid convective, the vapor convective, the evaporative, and the quenching terms are monitored to infer the boiling dynamics in the transition region from the single-phase flow to the nucleate boiling zone. It is found that the CFD approach to detecting the boiling incipience point agrees reasonably with available experimental data. However, its limitations, like inaccuracy for lower flow rates and early transition to boiling flow, are also noticed due to the nature of the considered RPI model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A generalized hyperbolic marching technique for three-dimensional supersonic flow with shocks

The current paper introduces a generalized, nonorthogonal coordinate system for application to supersonic flows with such large angles between streamlines of the initial Cauchy data and the marching direction that the velocity component in the marching direction is actually subsonic. Finite-difference operators are developed to solve the three-dimensional, steady, inviscid conservation equations of fluid flow referenced to this general frame. Moreover, a new method of aligning the difference mesh to the bow shock wave is presented which eliminates both the necessity for differencing the free-stream flow properties and the spurious fluctuations that often arise in conservative variables when they are differenced across the discontinuity. Results are given for the case of flow past a conical surface.

Rizzi, A. W.↗

Detection of Flooding Responses at the River Basin Scale Enhanced by Land use Change

The Georges Creek watershed (area 187.5 sq km) in western Maryland (United States) has experienced land use changes (>17% of area) associated with surface mining of coal. The adjacent Savage River watershed (area 127.2 sq km) is unmined. Moments of flood frequency distributions indicated that climatic variability affected both watersheds similarly. Normalizing annual maximum flows by antecedent streamflow and causative precipitation helped identify trends in flooding response. Analysis of contemporary storm events using Next Generation Weather Radar (NEXRAD) stage III precipitation data showed that Georges Creek floods are characterized by higher peak runoff and a shorter centroid lag than Savage River floods, likely attributable to differences in current land use. Interestingly, Georges Creek produces only two thirds of the storm-flow volume as Savage River, apparently because of infiltration into abandoned deep mine workings and an associated transbasin diversion constructed circa 1900. Empirical trend analysis is thus complicated by both hydroclimatic variability and the legacy of deep mining in the basin.

McCormick, Brian C.↗

Optimizing Flight Departure Delay and Route Selection Under En Route Convective Weather

This paper presents a linear Integer Programming model for managing air traffic flow in the United States. The decision variables in the model are departure delays and predeparture reroutes of aircraft whose trajectories are predicted to cross weather-impacted regions of the National Airspace System. The model assigns delays to a set of flights while ensuring their trajectories are free of any conflicts with weather. In a deterministic setting, there is no airborne holding due to unexpected weather incursion in a flight s path. The model is applied to solve a large-scale traffic flow management problem with realistic weather data and flight schedules. Experimental results indicate that allowing rerouting can reduce departure delays by nearly 57%, but it is associated with an increase in total airborne time due to longer routes flown by aircraft. The computation times to solve this problem were significantly lower than those reported in the earlier studies.

Mukherjee, Avijit↗

Penalty function finite element analysis of steady viscous incompressible flow in rotating coordinates

Finite element methods for incompressible viscous flow in turbomachines have not been presented in the literature previously. This paper develops a penalty function primitive variable method including Coriolis and centrifugal force terms for steady flow in a rotating coordinate system. Simplex elements are used with the result of solution times comparable to equivalent finite different solutions. Example cases considered are Couette flow, Poiseuile flow, flow over a step and flow in a rotating channel. Both laminar and turbulent flows are discussed. The accuracy of computed solutions compares well with theoretical solutions and experimenal measurements.

Rosen, M. C.↗

Combined Numerical/Analytical Perturbation Solutions of the Navier-Stokes Equations for Aerodynamic Ejector/Mixer Nozzle Flows

In spite of rapid advances in both scalar and parallel computational tools, the large number of variables involved in both design and inverse problems make the use of sophisticated fluid flow models impractical, With this restriction, it is concluded that an important family of methods for mathematical/computational development are reduced or approximate fluid flow models. In this study a combined perturbation/numerical modeling methodology is developed which provides a rigorously derived family of solutions. The mathematical model is computationally more efficient than classical boundary layer but provides important two-dimensional information not available using quasi-1-d approaches. An additional strength of the current methodology is its ability to locally predict static pressure fields in a manner analogous to more sophisticated parabolized Navier Stokes (PNS) formulations. To resolve singular behavior, the model utilizes classical analytical solution techniques. Hence, analytical methods have been combined with efficient numerical methods to yield an efficient hybrid fluid flow model. In particular, the main objective of this research has been to develop a system of analytical and numerical ejector/mixer nozzle models, which require minimal empirical input. A computer code, DREA Differential Reduced Ejector/mixer Analysis has been developed with the ability to run sufficiently fast so that it may be used either as a subroutine or called by an design optimization routine. Models are of direct use to the High Speed Civil Transport Program (a joint government/industry project seeking to develop an economically.viable U.S. commercial supersonic transport vehicle) and are currently being adopted by both NASA and industry. Experimental validation of these models is provided by comparison to results obtained from open literature and Limited Exclusive Right Distribution (LERD) sources, as well as dedicated experiments performed at Texas A&M. These experiments have been performed using a hydraulic/gas flow analog. Results of comparisons of DREA computations with experimental data, which include entrainment, thrust, and local profile information, are overall good. Computational time studies indicate that DREA provides considerably more information at a lower computational cost than contemporary ejector nozzle design models. Finally. physical limitations of the method, deviations from experimental data, potential improvements and alternative formulations are described. This report represents closure to the NASA Graduate Researchers Program. Versions of the DREA code and a user's guide may be obtained from the NASA Lewis Research Center.

DeChant, Lawrence Justin↗

Variational energy principle for compressible, baroclinic flow. 2: Free-energy form of Hamilton's principle

The first and second variations are calculated for the irreducible form of Hamilton's Principle that involves the minimum number of dependent variables necessary to describe the kinetmatics and thermodynamics of inviscid, compressible, baroclinic flow in a specified gravitational field. The form of the second variation shows that, in the neighborhood of a stationary point that corresponds to physically stable flow, the action integral is a complex saddle surface in parameter space. There exists a form of Hamilton's Principle for which a direct solution of a flow problem is possible. This second form is related to the first by a Friedrichs transformation of the thermodynamic variables. This introduces an extra dependent variable, but the first and second variations are shown to have direct physical significance, namely they are equal to the free energy of fluctuations about the equilibrium flow that satisfies the equations of motion. If this equilibrium flow is physically stable, and if a very weak second order integral constraint on the correlation between the fluctuations of otherwise independent variables is satisfied, then the second variation of the action integral for this free energy form of Hamilton's Principle is positive-definite, so the action integral is a minimum, and can serve as the basis for a direct trail and error solution. The second order integral constraint states that the unavailable energy must be maximum at equilibrium, i.e. the fluctuations must be so correlated as to produce a second order decrease in the total unavailable energy.

Schmid, L. A.↗

Computations of Internal and External Axisymmetric Nozzle Aerodynamics at Transonic Speeds

Computational Fluid Dynamics (CFD) analyses of axisymmetric circular-arc boattail nozzles have been completed in support of NASA's Next Generation Launch Technology Program to investigate the effects of high-speed nozzle geometries on the nozzle internal flow and the surrounding boattail regions. These computations span the very difficult transonic flight regime, with shock-induced separations and strong adverse pressure gradients. External afterbody and internal nozzle pressure distributions computed with the Wind code are compared with experimental data. A range of turbulence models were examined in Wind, including an Explicit Algebraic Stress model (EASM). Computations on two nozzle geometries have been completed at freestream Mach numbers ranging from 0.6 to 0.9, driven by nozzle pressure ratios (NPR) ranging from 2.9 to 5. Results obtained on converging-only geometry indicate reasonable agreement to experimental data, with the EASM and Shear Stress Transport (SST) turbulence models providing the best agreement. Calculations completed on a converging-diverging geometry involving large-scale internal flow separation did not converge to a true steady-state solution when run with variable timestepping (steady-state). Calculations obtained using constant timestepping (time-accurate) indicate less variations in flow properties compared with steady-state solutions. This failure to converge to a steady-state solution was found to be the result of difficulties in using variable time-stepping with large-scale separations present in the flow. Nevertheless, time-averaged boattail surface pressure coefficient and internal nozzle pressures show fairly good agreement with experimental data. The SST turbulence model demonstrates the best over-all agreement with experimental data.

Dalbello, Teryn↗

Lithospheric extension near Lake Mead, Nevada - A model for ductile flow in the lower crust

Small variations in gravity anomalies and topographic elevation observed in areas that have undergone highly variable amounts of upper crustal thinning can be satisfactorily explained by ductile flow of lower crustal material under the proper conditions. The boundary between the unextended Colorado Plateau and a strongly extended domain in the Basin and Range Province in the Lake Mead (Nevada) region is examined. Finite element modeling of Newtonian flow and power law creep shows that flow over the length scale of the eastern Basin and Range (500) km or more) corresponding to upper crustal extension by a factor of 1.4-3 over 10 million years requires effective viscosities less than 10 to the 18th - 10 to the 20th Pa s for ductile channels 10-25 km thick. Modeling suggests that these effective viscosities may be sustained by lower crustal material deforming at laboratory-derived power law creep rates. The longer-scale flow may require elevated crustal temperatures (more than 700 C), depending on the composition and material properties assumed. Under the boundary conditions assumed in this study the linear viscous flow models yield a satisfactory approximation to deformation by power law creep. This work suggests that flow in the lower crust may be a viable mechanism for producing small variations in total crustal thickness between strongly extended and less extended regions, and thereby explaining the relative uniformity in gravity and topography between such regions.

Kruse, Sarah↗

Intensity variations in plasma flow at the dawn magnetopause

Observations of plasma flows in the region of the dawn magnetopause obtained by the outbound Voyager 1 spacecraft, at a velocity of 11 km/sec, are discussed. Magnetic field and ion data obtained for the period surrounding four magnetopause crossings are presented which reveal energetic anti-sunward flowing ions outside the boundary with a time variability on the order of 400 millisec. These particle intensity variations, observed to vary with frequency, are most likely associated with the particle energization process or with the leakage of magnetospheric protons. The ion flows are considered to have originated sunward of the dawn meridian and were observed to penetrate approximately an ion gyroradius inside the dawn magnetopause.

Lanzerotti, L. J.↗

Applications of variational principles in computing rotational flows

Ecer and Akay (1983) have developed a variational formulation of rotational flow for Euler equations. The present paper provides a summary of these developments. The considered variational formulation provides a transformation of a type considered by Clebsch (1859). In this transformation, a new set of variables replaces the more commonly used primitive variables u(i), rho and p. Here, u(i) denotes the velocity components, while rho is the density, and p the pressure. The employed transformation produces a natural uncoupling of the equations when written in a quasi-linear form. After obtaining the governing equations in terms of the 'Clebsch variables', a solution scheme developed for calculating steady flows is discussed. Attention is given to numerical solutions of Euler equations based on the derived variational principles, and a study of inviscid, separated flows is conducted.

Ecer, A.↗

Development of an integrated BEM approach for hot fluid structure interaction

The development of a boundary element formulation for the study of hot fluid-structure interaction in earth-to-orbit engine hot section components is described. The initial primary thrust of the program to date was directed quite naturally toward the examination of fluid flow, since boundary element methods for fluids are at a much less developed state. This required the development of integral formulations for both the solid and fluid, and some preliminary infrastructural enhancements to a boundary element code to permit coupling of the fluid-structure problem. Boundary element formulations are implemented in two dimensions for both the solid and the fluid. The solid is modeled as an uncoupled thermoelastic medium under plane strain conditions, while several formulations are investigated for the fluid. For example, both vorticity and primitive variable approaches are implemented for viscous, incompressible flow, and a compressible version is developed. All of the above boundary element implementations are incorporated in a general purpose two-dimensional code. Thus, problems involving intricate geometry, multiple generic modeling regions, and arbitrary boundary conditions are all supported.

Dargush, Gary F.↗

Viscous flow computations using a second-order upwind differencing scheme

In the present computations of a wide range of fluid flow problems by means of the primitive variables-incorporating Navier-Stokes equations, a mixed second-order upwinding scheme approximates the convective terms of the transport equations and the scheme's accuracy is verified for convection-dominated high Re number flow problems. An adaptive dissipation scheme is used as a monotonic supersonic shock flow capture mechanism. Many benchmark fluid flow problems, including the compressible and incompressible, laminar and turbulent, over a wide range of M and Re numbers, are presently studied to verify the accuracy and robustness of this numerical method.

Chen, Y. S.↗

Quantifying contributions of natural variability and anthropogenic forcings on increased fire weather risk over the western United States

Previous studies have identified a recent increase in wildfire activity in the western United States (WUS). However, the extent to which this trend is due to weather pattern changes dominated by natural variability versus anthropogenic warming has been unclear. Using an ensemble constructed flow analogue approach, we have employed observations to estimate vapor pressure deficit (VPD), the leading meteorological variable that controls wildfires, associated with different atmospheric circulation patterns. Our results show that for the period 1979 to 2020, variation in the atmospheric circulation explains, on average, only 32% of the observed VPD trend of 0.48 ± 0.25 hPa/decade (95% CI) over the WUS during the warm season (May to September). The remaining 68% of the upward VPD trend is likely due to anthropogenic warming. The ensemble simulations of climate models participating in the sixth phase of the Coupled Model Intercomparison Project suggest that anthropogenic forcing explains an even larger fraction of the observed VPD trend (88%) for the same period and region. These models and observational estimates likely provide a lower and an upper bound on the true impact of anthropogenic warming on the VPD trend over the WUS. During August 2020, when the August Complex “Gigafire” occurred in the WUS, anthropogenic warming likely explains 50% of the unprecedented high VPD anomalies.

54 ENVIRONMENTAL SCIENCES↗

Testing- and Model- Based Optimization of Coal-fired Primary Heater Design for Indirect Supercritical CO 2 Power Cycles (Final Scientific and Technical Report)

The overall objective of this project was to perform the R&D necessary to mitigate the risk associated with the design of a primary heat exchanger for a solid-fired combustion system coupled with an indirect-fired closed-loop Brayton Cycle utilizing supercritical CO 2 . The key technological hurdle was the coupling of a solid-fuel firing system with the primary heater, which poses a singular challenge, which is the management of burner performance and operational conditions in a way to manage heat exchanger tube metal temperatures and temperature ramp rates in the absence of fluid phase change on the inside of the tubes. We designed and built the first ever pseudo power system employing a simple recuperated supercritical CO 2 closed-loop Brayton Cycle coupled to a solid-fuel fired system. Advanced coupled CFD and process modeling were used to design the primary heat exchanger (PHX), which consisted of both radiative and convective sections, to limit tube metal temperatures resulting from the heat release profile of the solid fuel flame near the radiative tubes. The heat exchanger was designed to produce finished CO 2 temperatures of 600 °C a pressure of 20.7 MPa and CO 2 flow of 5.5 kg/s. The constructed PHX was capable of 1.2 MWth heat uptake. During design of the PHX, the modeling showed that most variables influencing flame shape (burner stoichiometric ratio and register velocities and swirl) were not suitable to manage heat flux to the metal surfaces. This is because they substantially increased adiabatic flame temperature through the influence of localized stoichiometric ratio. Excess air and firing rate were the two most powerful variables that could be used to control tube surface temperatures. The coupled system was operated for a total of 407 hours, with the longest continuous run of 248 hours. For 62% of the operational time, the unit was unmanned and in automatic control. The fuels used for the testing included natural gas, two Utah Bituminous coals, woody biomass, and bagasse. During the testing we were able to verify the 1.2 MWth heat uptake and we operated at a finished CO 2 temperature of 607 °C and a pressure of 20.3 MPa simultaneously. The real-time corrosion rate of the Super 304H tube CO 2 surface in the region of the radiative section of the PHX were measured, at an approximate temperature of 550 °C. The two key variables related to corrosion rate are the pressure and flow rate of the CO 2 . A technoeconomic analysis was performed at a scale of 120 MWE. The updated analysis showed that the efficiency of an sCO 2 power producing plant will be related to the pressure drop of the PHX.

01 COAL, LIGNITE, AND PEAT↗

A Nonstationary and Non-Gaussian Moving Average Model for Solar Irradiance

Historically, power has flowed from large power plants to customers. Increasing penetration of distributed energy resources such as solar power from rooftop photovoltaic has made the distribution network a two-way-street with power being generated at the customer level. The incorporation of renewables introduces additional uncertainty and variability into the power grid. Distribution network operation studies are being adapted to include renewables; however, such studies require high quality solar irradiance data that adequately reflect realistic meteorological variability. Data from satellite-based products are spatially complete, but temporally coarse, whereas solar irradiances exhibit high frequency variation at very fine timescales. We propose a new stochastic method for temporally downscaling global horizontal irradiance (GHI) to 1 min resolution, but we do not consider the spatial aspect due to limited availability of the in situ irradiance measurements. Solar irradiance's first and second-order structures vary diurnally and seasonally, and our model adapts to such nonstationarity. Empirical irradiance data exhibits highly non-Gaussian behavior; we develop a nonstationary and non-Gaussian moving average model that is shown to capture realistic solar variability at multiple timescales. We also propose a new estimation scheme based on Cholesky factors of empirical autocovariance matrices, bypassing difficult and inaccessible likelihood-based approaches. The model is demonstrated for a case study of three locations that are located in diverse climates through the United States. The model is compared against competitors from the literature and is shown to provide better uncertainty and variability quantification on testing data.

Cholesky factor↗