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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Program Management for Concurrent University Satellite Programs, Including Propellant Feed System Design Elements

Propulsion options for CubeSats are limited but are necessary for the CubeSat industry to continue future growth. Challenges to CubeSat propulsion include volume/mass constraints, availability of sufficiently small and certified hardware, secondary payload status, safe storage of high energy, and power requirements. A multi-mode (chemical and electric) thruster was developed by at the Missouri University of Science and Technology (Missouri S&T) to enable CubeSat propulsion missions. Two satellite buses, a 3U and 6U, are being developed by the M-SAT satellite research team at Missouri S&T to demonstrate the multi-mode thruster's capabilities. Two key challenges related to these missions are the development of the feed system to support the thruster and management of the two programs' personnel, resources, timelines, and budgets. The feed system was designed to support the unique needs of the thruster including material compatibility with the custom ionic liquid, variable flow rate, high pressurization, and high voltage within the constraints and budget of a student-designed propulsion system. The feed system was also designed to minimize risk as a secondary payload. This resulted in the development of a unique method to pressurize propellant stored in the feed system tubing. Within the expected operating pressure range, the method was experimentally shown to provide sufficient pressure and propellant volume to the thruster to meet mission success criteria. The 3U and 6U CubeSat buses were designed concurrently with complimentary payloads, hardware, objectives, and team structures. Resources such as personnel, lab space, and time had to be carefully balanced between the two teams. With proper management, the two programs have been able to support one another, providing valuable insight and collaboration between the two teams. Lessons learned include experience with design, testing, and assembly of hardware, team training/mentoring and motivation, improved documentation practices, and risk management.

Withrow, Shannah Nichole↗
NASA NTRS2020

Variability of the Solar Wind Flow Asymmetry in the Martian Magnetosheath Observed by MAVEN

We perform the first statistical analysis of effects that different external conditions20have on a solar wind (SW) flow asymmetry observed in the Martian magnetosheath due to mass loading, making use of∼5 years of Mars Atmosphere and Volatile EvolutioN22(MAVEN) observations. We find the difference between the mean magnetosheath SW velocity component along the SW convective electric field direction in regions in the positive and negative hemispheres displays a strong linear correlation with the ratio between the upstream Interplanetary Magnetic Field (IMF) cross-flow component (By) and the SW proton density (nsw). The asymmetry is maximized (∼25 - 35 kms−1) for low nsw(∼1 cm−3) and large IMF By(∼2 nT). These results suggest the SW flow asymmetry variability is due to a force arising from the differential streaming between SW protons and oxygen ions, with boundary conditions partly defined by the pristine SW properties.

N Romanelli↗
NASA NTRSDate not supplied

Aeroacoustic Computations of a Generic Low Boom Concept in Landing Configuration: Part 3 - Aerodynamic Validation and Noise Source Identification

Wind tunnel test data were used to validate the predicted aerodynamic behavior of a 15%-scale version of a generic, low-boom aircraft. The test was conducted in the NASA Langley 14- by 22-Foot Subsonic Tunnel to determine the low-speed aerodynamic characteristics of the model. Measured steady surface pressures and global forces were used to validate predicted aerodynamic results obtained from high-fidelity simulations of the model as installed in the tunnel. Very good agreement between predicted and measured aerodynamic trends was demonstrated, providing the impetus to proceed with companion airframe noise simulations that were conducted in a free-air setting. Computed near-field flow variables acquired on a permeable data surface were used to generate synthetic pressure records for two 800-element phased microphone arrays positioned overhead and to the side of the model. The array data were beamformed to generate noise source localization maps for the clean model and several landing configurations. Primary and secondary airframe sources were identified and their relative strengths were determined. Far-field integrated noise spectra for the full aircraft, as well as individual components, were obtained from the source maps via integration of tailored regions. The analysis showed that noise produced by the landing gear was the dominant contributor to the far-field acoustic signature of the model, followed by flap noise. The effects of permeable data surface end caps, spatial resolution, array orientation, angle of attack, component interaction, velocity scaling, and numerical precision on synthetic far-field spectra were also evaluated.

low boom↗
NASA NTRSDate not supplied

Lessons Learned on the Use of Data Surfaces for Ffowcs Williams-Hawkings Calculations: Airframe Noise Applications

The use of integral solutions to the acoustic analogy-based Ffowcs Williams-Hawkings (FW-H) formulation is the defacto industrial and academic approach for assessing far-field acoustic behavior from airframe noise simulations. The methodology utilizes time-accurate flow variables obtained at arbitrary data surfaces to determine far-field spectra. The solid surface approach, where the time history of pressure is recorded on the surface of the aircraft, is generally preferred. Permeable surfaces, where pressure and velocity are recorded on a user-defined surface enclosing the aircraft, are also used. This paper seeks to share some of the valuable lessons we learned over the past several years on the limitations of the FW-H formulation and best practices for its successful application. We show results for various configurations comprising isolated landing gears, a business jet, a large commercial transport, and a generic low boom supersonic prototype, all at landing conditions. We demonstrate that the solid FW-H approach has significant accuracy limitations for most realistic cases and that the permeable FW-H approach may also fail in specific scenarios, requiring additional costly steps to extract the correct far-field noise spectra. Suggested guidelines for proper FW-H calculations are also given.

Acoustic propagation, Ffowcs Williams-Hawkings Cal↗
NASA NTRS2022

The Harmonic Linearized Navier-Stokes Equations for Transition Prediction in Three-Dimensional Flows

The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the non-local, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, as roughness elements, steps, gaps, or corners. To solve the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 107 to 108. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear system. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE.

Boundary Layer Stability↗
NASA NTRSDate not supplied

Numerical Analysis of Cavitation-Induced Two-Phase Fluid Transients in Cryogenic Feedlines Using Finite Volume Procedure

This article presents a numerical investigation of the cavitation induced fluid transient in cryogenic feedline due to sudden closure of the valve A one dimensional numerical model based on finite volume procedure consisting of nodes and branches using Generalized Fluid System Simulation Program ( has been employed for the prediction of pressure variations The unsteady friction formulation accurately determines wall shear stress, and a homogeneous mixture model is used for two phase flow modeling Numerical results obtained with the present formulation are compared with experimental results and Method of Characteristic ( based results It is observed that low residual in flow variable calculation improves the solution providing numerically stable results GFSSP is able to reproduce the fluid transient peak due to vapor collapse Subsequently, it can retain the shape of pressure wave for two phase fluid transients, but it is dampening faster compared to experiments Further investigation is required with set of experiments to validate the present numerical approach for different experimental data.

Fluid Transient↗
DOE OSTI2025

Techno-Economic Viability of Flexible Dispatch of Unconventional Geothermal Systems

Flexible geothermal operations could boost project returns through the allocation of improved power purchase agreements and/or exploitation of power price arbitrage opportunities. In this study, we investigated the techno-economic feasibility of variable flow rate control and time-of-day pricing in closed-loop geothermal systems. We considered U-shaped multilateral system configurations and modeled a variety of technical system parameters. These designs were simulated using a slender-body theory (SBT) model for transient heat transfer and fluid flow. This subsurface model was integrated into the flexible geothermal economic model (FGEM) tool to evaluate the overall flexible geothermal system techno-economics. Future hourly ambient temperature conditions were based on the Sup3rCC dataset. Published datasets were used for future hourly wholesale electricity prices. We analyzed four operating strategies: 1) baseload operation, 2) seasonal dispatch (high flow rate during summer and nominal flow rate during the rest of the year), 3) net generation maximization by varying flow rate to maximize net power output, and 4) revenue maximization by varying flow rate to maximize revenue. We ran all four scenarios for a multiloop configuration with 12 lateral passes, 7-km vertical depth and 87-km total drilling length. Furthermore, we assumed a 60 degrees C/km geothermal gradient and ambient temperature and wholesale electricity prices for New Mexico as a typical state location. The nominal flow rate was set to 80 kg/s. When considering drilling costs of $1,000/m and a discount rate of 7%, the generation maximization scenario resulted in the lowest levelized cost of electricity (LCOE) of ~$150/MWh. When considering project return on investment (ROI), defined as lifetime net income divided by upfront capital costs, all flexible operation scenarios performed better than the base case scenario. The highest ROI of 80% was obtained with the revenue maximization scenario. With drilling costs of $200/m and a discount rate of 5%, the generation maximization scenario resulted in LCOE of $49/MWh.

flexible geothermal↗
NASA NTRS1987

Time-marching solution of incompressible Navier-Stokes equations for internal flow

Primitive variables with central differencing on a staggered grid are used in the present, factored ADI finite-difference scheme for artificial compressibility method solution of the incompressible Navier-Stokes equations, leading to a close coupling between velocity and pressure that both enhances stability and eliminates the need for artificial damping. Computational efficiency is enhanced through the use of a spatially variable, fixed Courant number-based time-step. The numerical results obtained for a driven cavity at Re of 10,000, with local cell Re as high as 100, exhibits no flow variable spatial oscillations on a 40 x 40 stretched grid solution.

Soh, W. Y.↗
NASA NTRS1991

Choice of velocity variables for complex flow computation

The issue of adopting the velocity components as dependent velocity variables for the Navier-Stokes flow computations is investigated. The viewpoint advocated is that a numerical algorithm should preferably honor both the physical conservation law in differential form and the geometric conservation law in discrete form. With the use of Cartesian velocity vector, the momentum equations in curvilinear coordinates can retain the full conservation-law form and satisfy the physical conservation laws. With the curvilinear velocity components, source terms appear in differential equations and hence the full conservation law form can not be retained. In discrete expressions, algorithms based on the Cartesian components can satisfy the geometric conservation-law form for convection terms but not for viscous terms; those based on the curvilinear components, on the other hand, cannot satisfy the geometric conservation-law form for either convection or viscous terms. Several flow solutions for domain with 90 and 360 degree turnings are presented to illustrate the issues of using the Cartesian velocity components and the staggered grid arrangement.

Shyy, W.↗
NASA NTRS1974

The solution of three-variable duct-flow equations

This paper establishes a numerical method for the solution of three-variable problems and is applied here to rotational flows through ducts of various cross sections. An iterative scheme is developed, the main feature of which is the addition of a duplicate variable to the forward component of velocity. Two forward components of velocity result from integrating two sets of first order ordinary differential equations for the streamline curvatures, in intersecting directions across the duct. Two pseudo-continuity equations are introduced with source/sink terms, whose strengths are dependent on the difference between the forward components of velocity. When convergence is obtained, the two forward components of velocity are identical, the source/sink terms are zero, and the original equations are satisfied. A computer program solves the exact equations and boundary conditions numerically. The method is economical and compares successfully with experiments on bent ducts of circular and rectangular cross section where secondary flows are caused by gradients of total pressure upstream.

Stuart, A. R.↗
NASA NTRS1985

Evaluation of the satellite derived snow cover area - Runoff forecasting models for the inaccessible basins of western Himalayas

In this study, the existing seasonal snow cover area runoff forecasting models of the Indus, Kabul, Sutlej and Chenab basins were evaluated with the concurrent flow correlation model for the period 1975-79. In all the basins under study, correlation of concurrent flow model explained the variability in flow better than by the snow cover area runoff models. Actually, the concurrent flow correlation model explained more than 90 percent of the variability in the flow of these rivers. Compared to this model, the snow cover area runoff models explained less of the variability in flow. In the Himalayan river basins under study and at least for the period under observation, the concurrent flow correlation model provided a set of results with which to compare the estimates from the snow cover area runoff models.

Dey, B.↗
NASA NTRS1973

Propagation of plane waves in flow through a variable area duct

Duct geometry is an important factor that influences the transmission of sound in a duct. Nonconstant area produces variations in steady flow quantities through the duct, which can cause reflections of acoustic disturbances and the creation of standing wave patterns, and which can attenuate and disperse propagating waves. To study the problem, a second-order-accurate numerical method has been developed. The method of characteristics is used to solve the acoustic equations with the assumption of quasi-one-dimensional flow in the duct. Numerical results are presented over a broad range of frequencies. Low-frequency results are compared with approximate solutions for compact nozzles and high-frequency calculations are compared with results from the short-wave theory. The numerical results are also compared with an analytical solution of the acoustic equations for exponential horn sections with zero mean flow to further establish the validity of the numerical method.

King, L. S.↗
DOE OSTI2020

Extracting a mixing parameter from 2D radiographic imaging of variable-density turbulent flow

We extract a suitably averaged fluctuating density from the two-dimensional radiographic image of a flow. The X-ray attenuation is given by the Beer–Lambert law which exponentially damps the incident beam intensity by a factor proportional to the density, opacity and thickness of the target. By making reasonable assumptions for the mean density, opacity and effective thickness of the target flow, we estimate the density fluctuation contribution to the attenuation. The extracted density fluctuations averaged across the thickness of the flow in the direction of the beam may be used to form the density–specific-volume correlation $b$. In a statistical description of variable-density turbulence, $b$ quantifies the degree of mixedness. The ability to extract a measure of mixedness from experimental data would be a powerful tool that could be used in the validation of mix models. The scheme proposed is tested for DNS data computed for variable density buoyancy-driven mixing. We quantify the deficits in the extracted value of $b$ due to target thickness, Atwood number and modeled signal noise. This analysis justifies using the proposed scheme to infer the mix parameter from thin targets at moderate to low Atwood numbers. Furthermore, to illustrate how the scheme might be used in a practical problem, we demonstrate its application to a radiographic image of counter-shear flow obtained from experiments at the National Ignition Facility.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1983

Multi-grid solution of Neumann pressure problem for viscous flows using primitive variables

The multi-grid (MG) technique has been advanced for use with the Neumann boundary-value problem in clustered curvilinear orthogonal coordinates. This comprises an important step in the analysis of viscous flows using the velocity-pressure formulation of the Navier-Stokes equations. With successive over-relaxation (SOR) as the smoothing operator and with suitably formulated restriction and coarse-grid-correction operators, a 4-grid procedure enhances the efficiency of fine-grid solutions of the Neumann problem by a factor of 3 to 5, depending on the problem parameters. Thy influence of the smoothing operator is also examined by employing the alternating-direction implicit and the strongly implicit techniques instead of SOR.

Ghia, U.↗
NASA NTRS1994

Pool and flow boiling in variable and microgravity

As is well known, boiling is an effective mode of heat transfer in that high heat flux levels are possible with relatively small temperature differences. Its optimal application requires that the process be adequately understood. A measure of the understanding of any physical event lies in the ability to predict its behavior in terms of the relevant parameters. Despite many years of research the predictability of boiling is currently possible only for quite specialized circumstances, e.g., the critical heat flux and film boiling for the pool boiling case, and then only with special geometries. Variable gravity down to microgravity provides the opportunity to test this understanding, but possibly more important, by changing the dimensional and time scales involved permits more detailed observations of elements involved in the boiling process, and perhaps discloses phenomena heretofore unknown. The focus here is on nucleate boiling although, as will be demonstrated below, under but certain circumstances in microgravity it can take place concurrently with the dryout process. In the presence of earth gravity or forced convection effects, the latter process is usually referred to as film boiling. However, no vapor film as such forms with pool boiling in microgravity, only dryout. Initial results are presented here for pool boiling in microgravity, and were made possible at such an early date by the availability of the Get-Away-Specials (GAS). Also presented here are some results of ground testing of a flow loop for the study of low velocity boiling, eventually to take place also in microgravity. In the interim, variable buoyancy normal to the heater surface is achieved by rotation of the entire loop relative to earth gravity. Of course, this is at the expense of varying the buoyancy parallel to the heater surface. Two questions which must be resolved early in the study of flow boiling in microgravity are (1) the lower limits of liquid flow velocity where buoyancy effects become significant to the boiling process (2) the effect of lower liquid flow velocities on the Critical Heat Flux when buoyancy is removed. Results of initial efforts in these directions are presented, albeit restricted currently to the ever present earth gravity.

Merte, Herman, Jr.↗
NASA NTRS2004

Performance of Thermal Mass Flow Meters in a Variable Gravitational Environment

The performance of five thermal mass flow meters, MKS Instruments 179A and 258C, Unit Instruments UFM-8100, Sierra Instruments 830L, and Hastings Instruments HFM-200, were tested on the KC-135 Reduced Gravity Aircraft in orthogonal, coparallel, and counterparallel orientations relative to gravity. Data was taken throughout the parabolic trajectory where the g-level varied from 0.01 to 1.8 times normal gravity. Each meter was calibrated in normal gravity in the orthogonal position prior to flight followed by ground testing at seven different flow conditions to establish a baseline operation. During the tests, the actual flow rate was measured independently using choked-flow orifices. Gravitational acceleration and attitude had a unique effect on the performance of each meter. All meters operated within acceptable limits at all gravity levels in the calibrated orthogonal position. However, when operated in other orientations, the deviations from the reference flow became substantial for several of the flow meters. Data analysis indicated that the greatest source of error was the effect of orientation, followed by the gravity level. This work emphasized that when operating thermal flow meters in a variable gravity environment, it is critical to orient the meter in the same direction relative to gravity in which it was calibrated. Unfortunately, there was no test in normal gravity that could predict the performance of a meter in reduced gravity. When operating in reduced gravity, all meters indicated within 5 percent of the full scale reading at all flow conditions and orientations.

Brooker, John E.↗
DOE OSTI2021

The Leading Modes of Asian Summer Monsoon Variability as Pulses of Atmospheric Energy Flow

Monsoon rainfall anomalies are often organized into dynamical modes that are more predictable than rainfall at specific locations. Here a linear response function (LRF) is constructed to extract the neutral modes (NM) of the dynamical atmospheric system and their optimal forcing that govern precipitation variability, providing a dynamical underpinning for the empirically derived, recurrent rainfall patterns. The leading mode (NM1) features a shift in the Atlantic ITCZ, a drying over the Indian Ocean and Indian subcontinent, and a “South-Flood-North Drought (SFND)” pattern over eastern China, a pattern inferable from changes in the atmospheric energy flow forced by the corresponding optimal forcing and further enhanced by the water vapor and cloud feedbacks. Moreover, the slow evolving component of the NM1 in response to CO2 forcing can be skillfully predicted by running the ocean heat uptake through the LRF, implicating a deterministic ocean dynamic source for the NM1-like trend under climate warming.

monsoon mode of variability, EBM, climate change, ↗