Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic mass 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 271 records · Page 15

A physics-based model for frost buildup under turbulent flow using direct numerical simulations

We present a new model for frost buildup under turbulent (and laminar) flow using direct numerical simulations. The physical model consists of two layers, the air and the frost. The air layer is fully resolved and consists of solving for the velocity, temperature, and vapor mass fraction fields. The frost layer thickness is resolved using conservation of mass and energy. Both phases are dynamically coupled using the immersed boundary method. Three-dimensional simulations are conducted in an open-channel configuration. A number of challenges need to be overcome to make these simulations feasible. First, to enforce far-field conditions of zero gradient and prescribed mean temperature and humidity, a source term is added to the energy and transport equations in the flow solver. Second, the mean frost thickness is subtracted after each time step to ensure a constant mean flow thickness and level of turbulence in the numerical domain. Third, a slow-time acceleration approach, which accelerates the frost buildup by a predetermined factor, is employed to bridge the gap between the fast turbulent and slow frost buildup time scales. Finally, a frost densification scheme is used to overcome the difficulties of vertically varying frost properties. The model is validated by comparing the frost thickness and frost thickness buildup rate over a period of one hour from a cooled flat plate experiment. As a result, both quantities compare favorably with experiments.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

System risk quantification and decision making support using functional modeling and dynamic Bayesian network

Risk-informed decision-making requires a probabilistic assessment of the likelihood of success of control action, given the system status. This paper presents a systematic state transition modeling approach integrating dynamic probabilistic risk assessment with a decision-making process using a dynamic Bayesian network (DBN) coupled with functional modeling. A functional model designed with multilevel flow modeling (MFM) technique was used to build a system state structure inferred by energy, mass, and information flow so that one can verify the developed model with respect to system functionality. The MFM model represents the causal relationship among the nodes, which captures the structure of process parameters and control units. Each node may have multiple possible states, and the DBN structured by the MFM model represents the time-domain transitions among the defined states. Furthermore, the MFM-DBN integrated state transition modeling is a white-box approach that allows one to draw the system's risk profile by updating the system states and supports the decisions probabilistically with physical inference. An example of a simple heating system has been used to illustrate this process, including decision-making support based on quantitative risk profile. For demonstrating its applicability to a complex system operational decision making, a case study of station blackout accident scenario leading to the seal loss of coolant accident in a nuclear power plant is presented. The proposed approach effectively provided the risk profile along time for each option so that the operators can make the best decision, which minimizes the plant risk.

42 ENGINEERING↗

Aeroelastic Calculations of Quiet High- Speed Fan Performed

An advanced high-speed fan was recently designed under a cooperative effort between the NASA Glenn Research Center and Honeywell Engines & Systems. The principal design goals were to improve performance and to reduce fan noise at takeoff. Scale models of the Quiet High-Speed Fan were tested for operability, performance, and acoustics. During testing, the fan showed significantly improved noise characteristics, but a self-excited aeroelastic vibration known as flutter was encountered in the operating range. Flutter calculations were carried out for the Quiet High-Speed Fan using a three-dimensional, unsteady aerodynamic, Reynolds-averaged Navier-Stokes turbomachinery code named "TURBO." The TURBO code can accurately model the viscous flow effects that can play an important role in various aeroelastic problems such as flutter with flow separation, flutter at high loading conditions near the stall line (stall flutter), and flutter in the presence of shock and boundary-layer interaction. Initially, calculations were performed with no blade vibrations. These calculations were at a constant rotational speed and a varying mass flow rate. The mass flow rate was varied by changing the backpressure at the exit boundary of the computational domain. These initial steady calculations were followed by aeroelastic calculations in which the blades were prescribed to vibrate harmonically in a natural mode, at a natural frequency, and with a fixed interblade phase angle between adjacent blades. The AE-prep preprocessor was used to interpolate the in-vacuum mode shapes from the structural dynamics mesh onto the computational fluid dynamics mesh and to smoothly propagate the grid deformations from the blade surface to the interior points of the grid. The aeroelastic calculations provided the unsteady aerodynamic forces on the blade surface due to blade vibrations. These forces were vector multiplied with the structural dynamic mode shape to calculate the work done on the blade during one vibration period, then this result was converted to an aerodynamic damping. Flutter occurs when the aerodynamic damping becomes negative, if structural damping is ignored. The results of these aeroelastic calculations are summarized in this plot of aerodynamic damping versus mass flow rate at a constant rotational speed. As the backpressure is increased, the mass flow rate through the fan decreases and the fan operating point moves towards the stall line. The aeroelastic calculations showed that the aerodynamic damping decreases as the stall line is approached, as observed during testing. In addition, the aeroelastic calculations with the TURBO code correctly predicted the aeroelastic parameters: the most unstable vibration mode and interblade phase angle, as observed during testing. The Quiet High-Speed Fan demonstrated significant noise reductions during testing, but flutter imposed limits on its operating range. The accurate calculation of the aeroelastic characteristics using the TURBO code is a significant step toward eliminating flutter from the operating range and toward realizing the benefits of reduced fan noise. The aeroelastic calculations described here were performed under a grant by University of Toledo researchers in collaboration with Glenn's researchers.

Bakhle, Milind A.↗

Mixing and NOx Emission Calculations of Confined Reacting Jet Flows in Cylindrical and Annular Ducts

Rapid mixing of cold lateral jets with hot cross-stream flows in confined configurations is of practical interest in gas turbine combustors as it strongly affects combustor exit temperature quality, and gaseous emissions in for example rich-lean combustion. It is therefore important to further improve our fundamental understanding of the important processes of dilution jet mixing especially when the injected jet mass flow rate exceeds that of the cross-stream. The results reported in this report describe some of the main flow characteristics which develop in the mixing process in a cylindrical duct. A three-dimensional computational fluid dynamics (CFD) code has been used to predict the mixing flow field characteristics and NOx emission in a quench section of a rich-burn/quick-mix/lean-burn (RQL) combustor. Sixty configurations have been analyzed in both circular and annular geometries in a fully reacting environment simulating the operating condition of an actual RQL gas turbine combustion liner. The evaluation matrix was constructed by varying the number of orifices per row and orifice shape. Other parameters such as J (momentum-flux ratio), MR (mass flowrate ratio), DR (density ratio), and mixer sector orifice ACd (effective orifice area) were maintained constant throughout the entire study. The results indicate that the mixing flow field can be correlated with the NOx production if they are referenced with the stoichiometric equivalence ratio value and not the equilibrium value. The mixing flowfields in both circular and annular mixers are different. The penetration of equal jets in both annular and circular geometries is vastly different which significantly affects the performance of the mixing section. In the computational results with the circular mixer, most of the NOx formation occurred behind the orifice starting at the orifice wake region. General trends have been observed in the NOx production as the number of orifices is changed and this appears to be common for all hole configurations and mixer types (circular or annular). The performance of any orifice shape (in producing minimum NOx) appears to be acceptable if the number of orifices can be freely varied in order to attain the optimum jet penetration.

Oechsle, Victor L.↗

Embedding hard physical constraints in neural network coarse-graining of three-dimensional turbulence

In recent years, deep learning approaches have shown much promise in modeling complex systems in the physical sciences. A major challenge in deep learning of partial differential equations is enforcing physical constraints and boundary conditions. In this work, we propose a general framework to directly embed the notion of an incompressible fluid into convolutional neural networks, and apply this to coarse-graining of turbulent flow. These physics-embedded neural networks leverage interpretable strategies from numerical methods and computational fluid dynamics to enforce physical laws and boundary conditions by taking advantage the mathematical properties of the underlying equations. Here, we demonstrate results on three-dimensional fully developed turbulence, showing that this technique drastically improves local conservation of mass, without sacrificing performance according to several other metrics characterizing the fluid flow.

97 MATHEMATICS AND COMPUTING↗

Unsteady Navier-Stokes computations on a wing-body configuration in ramp motions

Unsteady Navier-Stokes computations are conducted for transonic flows over a wing-body configuration undergoing prescribed ramp motions. The ramp motion from 0 to 15 deg includes angles of attack where vortex breakdown is observed experimentally. The vortex breakdown is found to be delayed until after the ramp motion ends. The dynamic effects on the loads are also demonstrated. The moment coefficient is found to be sensitive to the effect of the virtual mass of the fluid. To verify the numerical results, a grid refinement study has been performed for both steady and unsteady flow conditions, by using up to approximately one million grid points. The results show that the coarse grid can give integrated quantities reasonably well, whereas the finer grids give a more detailed flow structure. Comparisons are also made with available steady-state experimental data.

Obayashi, Shigeru↗

Assessment of the Fluid Dynamics Boundary Condition in Ablating or Blowing Flows

Improved models of ablative thermal protection systems have enabled the treatment of materials and fluid behavior in a coupled manner. This paper reports a new approach to modeling the interface between fluid and material, with attention to the conservation of species mass flux and energy on the fluid side of the interface. The general equation is presented and is shown to recover the traditional uncoupled fluid/materials response interface. Including the chemical reaction terms on the CFD side of the interface makes the heat flux exchange independent of the thermodynamic reference state and, therefore, a measurable quantity. Doing so allows the material response solver to take as input the surface heat flux rather than a film coefficient. Removing the film coefficient approximation enables more direct solution of vehicle thermal response but requires consistency in the wall state. The mixing of the shock layer and pyrolysis gas is then computed with finite rate chemistry within the fluid solver. The boundary conditions described have been implemented in the DPLR v4.05.1 code. Char removal is captured using finite rate chemistry in DPLR’s gas surface interaction module. Aspects of coupling these solutions to material response are discussed.

Ablation↗

Assessment of the Fluid Dynamics Boundary Condition in Ablating or Blowing Flows

Improved models of ablative thermal protection systems have enabled the treatment of materials and fluid behavior in a coupled manner. This paper reports a new approach to modeling the interface between fluid and material, with attention to the conservation of species mass flux and energy on the fluid side of the interface. The general equation is presented and is shown to recover the traditional uncoupled fluid/materials response interface. Including the chemical reaction terms on the CFD side of the interface makes the heat flux exchange independent of the thermodynamic reference state and, therefore, a measurable quantity. Doing so allows the material response solver to take as input the surface heat flux rather than a film coefficient. Removing the film coefficient approximation enables more direct solution of vehicle thermal response but requires consistency in the wall state. The mixing of the shock layer and pyrolysis gas is then computed with finite rate chemistry within the fluid solver. The boundary conditions described have been implemented in the DPLR v4.05.1 code. Char removal is captured using finite rate chemistry in DPLR’s gas surface interaction module. Aspects of coupling these solutions to material response are discussed.

Ablation↗

Thermal Protection System Evaluation Using Arc-jet Flows: Flight Simulation or Research Tool?

The arc-jet has been used to evaluate thermal protection systems (TPS) and materials for the past forty years. Systems that have been studied in this environmerd include ablators, active, and passive TPS concepts designed for vehicles entering planetary and Earth atmospheres. The question of whether arc-jet flow can simulate a flight environment or is it a research tool that provides an aero-thermodynamic heating environment to obtain critical material properties will be addressed. Stagnation point tests in arc-jets are commonly used to obtain material properties such as mass loss rates, thermal chemical stability data, optical properties, and surface catalytic efficiency. These properties are required in computational fluid dynamic codes to accurately predict the performance of a TPS during flight. Special facilities have been developed at NASA Ames Research Center to approximate the flow environment over the mid-fuselage and body flap regions of proposed space-planes type vehicles. This paper compares flow environments generated in flight over a vehicle with those created over an arc-jet test articles in terms of scale, chemistry, and fluid dynamic properties. Flight experiments are essential in order to validate the material properties obtained from arc-jet tests and used to predict flight performance of any TPS being considered for use on a vehicle entering the Earth atmosphere at hypersonic speed.

Stewart, David A.↗

Challenges to Reliable Production Nitrogen Doping of Nb for SRF Accelerating Cavities

Over the last several years, alloying of the surface layer of niobium SRF cavities has been demonstrated to beneficially lower the superconducting RF surface resistance. Nitrogen, titanium, and oxygen have all been demonstrated as effective alloying agents, occupying interstitial sites in the niobium lattice within the RF penetration depth and even deeper, when allowed to thermally diffuse into the surface at appropriate temperatures. The use of nitrogen for this function has been often termed ’nitrogen doping’ and is being applied in the LCLS-II and LCLS-II HE projects. We report characterization studies of the distribution of nitrogen into the exposed niobium surface and how such distribution is affected by the quality of the vacuum furnace environment in which the doping takes place, and the complexity of nitride crystal growth on different grain orientations of surface niobium. Using state-of-the-art quantification methods by dynamic secondary ion mass spectrometry (SIMS) depth profiling in niobium, we identify several phenomena involving furnace-sourced contamination. We also highlight a potential issue with N₂ flow constraints from the flange ’caps’ used during heat treatments.

Reece, Charles↗

Analysis of nitrogen condensation in an expanding nozzle flow

Condensation of nitrogen flow in an expanding nozzle flow is analyzed using one-dimensional gas dynamic equations and the equations for nucleation and droplet growth. Effects of variations in the Tolman constant and the mass accommodation factor are discussed as well as the effect of foreign nuclei. Comparisons are made with experimental data obtained from a small, contoured nozzle.

Wang, F. C.↗

Evolution of the Orszag-Tang vortex system in a compressible medium. II - Supersonic flow

A study is presented on the effect of embedded supersonic flows and the resulting emerging shock waves on phenomena associated with MHD turbulence, including reconnection, the formation of current sheets and vortex structures, and the evolution of spatial and temporal correlations among physical variables. A two-dimensional model problem, the Orszag-Tang (1979) vortex system, is chosen, which involves decay from nonrandom initial conditions. The system is doubly periodic, and the initial conditions consist of single-mode solenoidal velocity and magnetic fields, each containing X points and O points. The initial mass density is flat, and the initial pressure fluctuations are incompressible, balancing the local forces for a magnetofluid of unit mass density. Results on the evolution of the local structure of the flow field, the global properties of the system, and spectral correlations are presented. The important dynamical properties and observational consequences of embedded supersonic regions and emerging shocks in the Orszag-Tang model of an MHD system undergoing reconnection are discussed. Conclusions are drawn regarding the effects of local supersonic regions on MHD turbulence.

Picone, J. Michael↗

UV Observations of Prominence Activation and Cool Loop Dynamics

In this paper we investigate the thermal and dynamic properties of dynamic structures in and around a prominence channel observed on the limb on 17 April 2003. Observations were taken with the Solar and Heliospheric Observatory's Solar Ultraviolet Measurements of Emitted Radiation (SOHO/SUMER) in lines formed at temperatures from 80,000 to 1.6 MK. The instrument was pointed to a single location and took a series of 90 s exposures. Two-dimensional context was provided by the Transition Region and Coronal Explorer (TRACE) in the UV and EUV and the Kanzelhohe Solar Observatory in H-alpha. Two dynamic features were studied in depth: an activated prominence and repeated motions in a loop near the prominence. We calculated three-dimensional geometries and trajectories, differential emission measure, and limits on the mass, pressure, average density, and kinetic and thermal energies. These observations provide important tests for models of dynamics in prominences and cool (approx. 10(exp 5) K)loops, which will ultimately lead to a better understanding the mechanism(s) leading to energy and mass flow in these solar features.

Kucera, Therese A.↗

Large eddy simulation of particle dispersion analysis in a pressurized reactive fluid–particle system

Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in power generation. Designing an effective burner plays a vital role in developing new combustion technologies. Because of high pressure, the volumetric fraction of particles in a pressurized oxy-combustor could be higher than that in the conventional combustor, consequently, the particle dynamics in the pressurized vessel need to be evaluated accordingly. Specifically, the present computational work analyzes the pattern of the particle dispersion by extracting the particle mass flow rate and the particle size distribution in the pressurized combustor. A large eddy simulation (LES) with the environment pressure 15-bar and 100 kWth thermal input of fuel particle has been performed. The particle mass flow rates sampled in several cross-section planes strongly fluctuate, indicating the impacts of the flame as well as the gas-phase flow on the particle movement. Animations of the particle-gas flow show that the flow recirculation plays a major role in producing the pulsing phenomenon of the particle mass flow rate.

Li, Lei↗

Forebody and Inlet Design for the HIFiRE 2 Flight Test

A forebody and inlet have been designed for the HIFiRE 2 scramjet flight test. The test will explore the operating, performance, and stability characteristics of a simple hydrocarbon-fueled scramjet combustor as it transitions from dual-mode to scramjet-mode operation and during supersonic combustion at Mach 8+ flight conditions. Requirements for the compression system were derived from inlet starting and combustor inflow requirements as well as physical size constraints. The design process is described. A planar, fixed geometry, mixed compression concept was used to produce laterally uniform flow at the inlet entrance and a conservative amount of internal contraction with respect to inlet starting. A grid sensitivity study was performed so that important flow physics caused by three-dimensional shock boundary layer interactions could be captured with confidence. Results from low Mach number operability studies, nominal trajectory cases, and high dynamic pressure heat load cases are discussed. The forebody and inlet solutions provide information for on-going combustor calculations, mass capture across the trajectory for fuel system design, and surface heating rates for thermal/structural analysis. The design has a one freestream Mach number margin for inlet starting, exceeds the high Mach number combustor entrance pressure requirement, produces high quality flow at the inlet exit for all Mach numbers and vehicle attitudes in the design space, and fits inside the booster shroud.

Ferlemann, Paul G.↗

Optimization-Inspired Pin-Fin Array for Supercritical Carbon Dioxide Recuperator

Additively manufactured heat exchangers are one possible route to cost-effective sCO 2 power cycles. In this paper, experimental results are obtained for two helical pin fin tubes that were designed following parametric optimization of the fin array. Neither the numerical optimization nor the experimental testing have been previously reported in the literature. The two tube designs were, (1) the optimization inspired design, and (2) the optimization-inspired design with a fin diameter increased by a factor of 2. To characterize the print, the designs were scanned using X-ray computed tomography to measure feature sizes and heat transfer area. The optimization was conducted in a commercial, computational fluid dynamics code. The code solved the Reynolds Averaged Navier-Stokes (RANS) and energy equations with turbulence closure provided by the shear stress transport (SST) k-ω model. In the experiments, the Nusselt number augmentation was measured using the Wilson plot technique and the friction factor was determined with mass flow and pressure drop measurements. The experimental testing indicated that the optimization-inspired design had a friction factor that was four times less than the baseline tube design at equal Nusselt number. Additionally, the optimization-inspired design had a 14% improvement in Nusselt number, at equal friction factor, relative to the best performing computational fluid dynamics (CFD) trial points. Tube design (2), with the larger diameter pin fins, had similar performance, within experimental error, as the tube with the smaller diameter pin fins (1). Both tubes achieved overall fin array efficiencies near 1. A performance factor, V/V 0 , equal to the volume of the enhanced heat exchanger divided by the volume of the baseline (no-fins) heat exchanger, is recommended to quantify internal cooling performance. In conclusion, the experimental shell and tube heat exchanger, using the additively manufactured tube, is competitive with printed circuit heat exchangers in its pressure drop class and could be further improved by optimizing a shell-and-tube heat exchanger utilizing this heat transfer enhancement feature.

42 ENGINEERING↗

Optimization of Hydrogen Production from Pyrolysis of Biomass Waste [Poster]

We developed an CFD aided ML-based tool for rapid assessment and optimization of different compositions of biomass in a fluidized bed reactor. First, we use CFD to simulate fluidized bed reactors with known inlet biomass mixtures and obtain corresponding syngas yields. A lumped kinetic mechanisms represents the conversion of cellulose, hemicellulose, and lignin, as well as subsequent cracking of tars into non-condensable gases (H 2 , CO, CO 2 , CH 4 ). We use Bayesian analysis/optimization to obtain the ideal operational and mass flow conditions for hydrogen production.

08 HYDROGEN↗

Nonreactive mixing study of a scramjet swept-strut fuel injector

The results are presented of a cold-mixing investigation performed to supply combustor design information and to determine optimum normal fuel-injector configurations for a general scramjet swept-strut fuel injector. The experimental investigation was made with two swept struts in a closed duct at a Mach number of 4.4 and a nominal ratio of jet mass flow to air mass flow of 0.0295, with helium used to simulate hydrogen fuel. Four injector patterns were evaluated; they represented the range of hole spacing and the ratio of jet dynamic pressure to free-stream dynamic pressure. Helium concentration, pitot pressure, and static pressure in the downstream mixing region were measured to generate the contour plots needed to define the mixing-region flow field and the mixing parameters. Experimental results show that the fuel penetration from the struts was less than the predicted values based on flat-plate data; but the mixing rate was faster and produced a mixing length less than one-half that predicted.

Mcclinton, C. R.↗