Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLOW RATE”

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 307 records · Page 17

High Efficiency Centrifugal Compressor for Rotorcraft Applications

The report "High Efficiency Centrifugal Compressor for Rotorcraft Applications" documents the work conducted at UTRC under the NRA Contract NNC08CB03C, with cost share 2/3 NASA, and 1/3 UTRC, that has been extended to 4.5 years. The purpose of this effort was to identify key technical barriers to advancing the state-of-the-art of small centrifugal compressor stages; to delineate the measurements required to provide insight into the flow physics of the technical barriers; to design, fabricate, install, and test a state-of-the-art research compressor that is representative of the rear stage of an axial-centrifugal aero-engine; and to acquire detailed aerodynamic performance and research quality data to clarify flow physics and to establish detailed data sets for future application. The design activity centered on meeting the goal set outlined in the NASA solicitation-the design target was to increase efficiency at higher work factor, while also reducing the maximum diameter of the stage. To fit within the existing Small Engine Components Test Facility at NASA Glenn Research Center (GRC) and to facilitate component re-use, certain key design parameters were fixed by UTRC, including impeller tip diameter, impeller rotational speed, and impeller inlet hub and shroud radii. This report describes the design effort of the High Efficiency Centrifugal Compressor stage (HECC) and delineation of measurements, fabrication of the compressor, and the initial tests that were performed. A new High-Efficiency Centrifugal Compressor stage with a very challenging reduction in radius ratio was successfully designed, fabricated and installed at GRC. The testing was successful, with no mechanical problems and the running clearances were achieved without impeller rubs. Overall, measured pressure ratio of 4.68, work factor of 0.81, and at design exit corrected flow rate of 3 lbm/s met the target requirements. Polytropic efficiency of 85.5 percent and stall margin of 7.5 percent were measured at design flow rate and speed. The measured efficiency and stall margin were lower than pre-test CFD predictions by 2.4 percentage points (pt) and 4.5 pt, respectively. Initial impressions from the experimental data indicated that the loss in the efficiency and stall margin can be attributed to a design shortfall in the impeller. However, detailed investigation of experimental data and post-test CFD simulations of higher fidelity than pre-test CFD, and in particular the unsteady CFD simulations and the assessment with a wider range of turbulence models, have indicated that the loss in efficiency is most likely due to the impact of unfavorable unsteady impeller/diffuser interactions induced by diffuser vanes, an impeller/diffuser corrected flow-rate mismatch (and associated incidence levels), and, potentially, flow separation in the radial-to-axial bend. An experimental program with a vaneless diffuser is recommended to evaluate this observation. A subsequent redesign of the diffuser (and the radial-to-axial bend) is also recommended. The diffuser needs to be redesigned to eliminate the mismatching of the impeller and the diffuser, targeting a slightly higher flow capacity. Furthermore, diffuser vanes need to be adjusted to align the incidence angles, to optimize the splitter vane location (both radially and circumferentially), and to minimize the unsteady interactions with the impeller. The radial-to-axial bend needs to be redesigned to eliminate, or at least minimize, the flow separation at the inner wall, and its impact on the flow in the diffuser upstream. Lessons were also learned in terms of CFD methodology and the importance of unsteady CFD simulations for centrifugal compressors was highlighted. Inconsistencies in the implementation of a widely used two-equation turbulence model were identified and corrections are recommended. It was also observed that unsteady simulations for centrifugal compressors require significantly longer integration times than what is current practice in industry.

Gorazd Medic↗

High Efficiency Centrifugal Compressor for Rotorcraft Applications

The report "High Efficiency Centrifugal Compressor for Rotorcraft Applications" documents the work conducted at UTRC under the NRA Contract NNC08CB03C, with cost share 2/3 NASA, and 1/3 UTRC, that has been extended to 4.5 years. The purpose of this effort was to identify key technical barriers to advancing the state-of-the-art of small centrifugal compressor stages; to delineate the measurements required to provide insight into the flow physics of the technical barriers; to design, fabricate, install, and test a state-of-the-art research compressor that is representative of the rear stage of an axial-centrifugal aero-engine; and to acquire detailed aerodynamic performance and research quality data to clarify flow physics and to establish detailed data sets for future application. The design activity centered on meeting the goal set outlined in the NASA solicitation-the design target was to increase efficiency at higher work factor, while also reducing the maximum diameter of the stage. To fit within the existing Small Engine Components Test Facility at NASA Glenn Research Center (GRC) and to facilitate component re-use, certain key design parameters were fixed by UTRC, including impeller tip diameter, impeller rotational speed, and impeller inlet hub and shroud radii. This report describes the design effort of the High Efficiency Centrifugal Compressor stage (HECC) and delineation of measurements, fabrication of the compressor, and the initial tests that were performed. A new High-Efficiency Centrifugal Compressor stage with a very challenging reduction in radius ratio was successfully designed, fabricated and installed at GRC. The testing was successful, with no mechanical problems and the running clearances were achieved without impeller rubs. Overall, measured pressure ratio of 4.68, work factor of 0.81, and at design exit corrected flow rate of 3 lbm/s met the target requirements. Polytropic efficiency of 85.5 percent and stall margin of 7.5 percent were measured at design flow rate and speed. The measured efficiency and stall margin were lower than pre-test CFD predictions by 2.4 percentage points (pt) and 4.5 pt, respectively. Initial impressions from the experimental data indicated that the loss in the efficiency and stall margin can be attributed to a design shortfall in the impeller. However, detailed investigation of experimental data and post-test CFD simulations of higher fidelity than pre-test CFD, and in particular the unsteady CFD simulations and the assessment with a wider range of turbulence models, have indicated that the loss in efficiency is most likely due to the impact of unfavorable unsteady impeller/diffuser interactions induced by diffuser vanes, an impeller/diffuser corrected flow-rate mismatch (and associated incidence levels), and, potentially, flow separation in the radial-to-axial bend. An experimental program with a vaneless diffuser is recommended to evaluate this observation. A subsequent redesign of the diffuser (and the radial-to-axial bend) is also recommended. The diffuser needs to be redesigned to eliminate the mismatching of the impeller and the diffuser, targeting a slightly higher flow capacity. Furthermore, diffuser vanes need to be adjusted to align the incidence angles, to optimize the splitter vane location (both radially and circumferentially), and to minimize the unsteady interactions with the impeller. The radial-to-axial bend needs to be redesigned to eliminate, or at least minimize, the flow separation at the inner wall, and its impact on the flow in the diffuser upstream. Lessons were also learned in terms of CFD methodology and the importance of unsteady CFD simulations for centrifugal compressors was highlighted. Inconsistencies in the implementation of a widely used two-equation turbulence model were identified and corrections are recommended. It was also observed that unsteady simulations for centrifugal compressors require significantly longer integration times than what is current practice in industry.

Centrifugal compressor↗

Acoustic Study of a Sweeping Jet Actuator for Active Flow Control Applications

Active flow control for aircraft has been demonstrated to improve aerodynamic performance but generally introduces additional noise sources. This study presents measurements of the acoustic fields of sweeping jet actuators at full- and half-scale compared to the configuration recently demonstrated in ground- and flight-tests of a full-scale rudder. The effects of scale, flow rate, forward speed, and actuator synchronization state on tone and broadband noise directional power spectral density (PSD) levels are presented. The sweeping jet actuator PSD levels in the audio range increases as the second power of actuator flow rate. A fourth power variation of PSD level with flow rate is observed for a comparable non-sweeping actuator. These measurements will be used to develop noise source models for new aircraft configurations utilizing active flow control.

Active Flow Control Noise↗

Liquid Bismuth Propellant Management System for the Very High Specific Impulse Thruster with Anode Layer

Two prototype bismuth propellant feed systems were constructed and operated in conjunction with a propellant vaporizer. One system provided bismuth to a vaporizer using gas pressurization but did not include a means to measure the flow rate. The second system incorporated an electromagnetic pump to provide fine control of the hydrostatic pressure and a new type of in-line flow sensor that was developed for accurate, real-time measurement of the mass flow rate. High-temperature material compatibility was a driving design requirement for the pump and flow sensor, leading to the selection of Macor for the main body of both components. Posttest inspections of both components revealed no degradation of the material. The gas pressurization system demonstrated continuous pressure control over a range from zero to 200 torr. In separate proof-of-concept experiments, the electromagnetic pump produced a linear pressure rise as a function of current that compared favorably with theoretical pump pressure predictions, producing a pressure rise of 10 kPa at 30 A. Preliminary flow sensor operation indicated a bismuth flow rate of 6 mg/s with an uncertainty of plus or minus 6%. An electronics suite containing a real-time controller was successfully used to control the entire system, simultaneously monitoring all power supplies and performing data acquisition duties.

Polzin, K. A.↗

Performance, Facility Pressure Effects, and Stability Characterization Tests of NASA's Hall Effect Rocket with Magnetic Shielding Thruster

NASAs Hall Effect Rocket with Magnetic Shielding (HERMeS) 12.5 kW Technology Demonstration Unit-1 (TDU-1) Hall thruster has been the subject of extensive technology maturation in preparation for flight system development. Part of the technology maturation effort included experimental evaluation of the TDU-1 thruster with conducting and dielectric front pole cover materials in two different electrical configurations. A graphite front pole cover thruster configuration with the thruster body electrically tied to cathode and an alumina front pole cover thruster configuration with the thruster body floating were evaluated. Both configurations were also evaluated at different facility background pressure conditions to evaluate background pressure effects on thruster operation. Performance characterization tests found that higher thruster performance was attained with the graphite front pole cover configuration with the thruster electrically tied to cathode. A total thrust efficiency of 68 and a total specific impulse of 2,820 s was demonstrated at a discharge voltage of 600 V and a discharge power of 12.5 kW. Thruster stability regimes were characterized with respect to the thruster discharge current oscillations and with maps of the current-voltage-magnetic field (IVB). Analysis of TDU-1 discharge current waveforms found that lower normalized discharge current peak-to-peak and root mean square magnitudes were attained when the thruster was electrically floated with alumina front pole covers. Background pressure effects characterization tests indicated that the thruster performance and stability was mostly invariant to changes in the facility background pressure for vacuum chamber pressure below 110-5 Torr-Xe (for thruster flow rate above 8 mgs). Power spectral density analysis of the discharge current waveform showed that increasing the vacuum chamber background pressure resulted in a higher discharge current dominant frequency. Finally the IVB maps of the TDU-1 thruster taken at elevated magnetic fields indicated that the discharge current became more oscillatory with increased facility background pressure at lower thruster mass flow rates, where thruster operation at higher flow rates resulted in less change to the thrusters IVB characteristics.

Hall Thrusters↗

Investigation of Unsteady Tip Clearance Flow in a Low-Speed One and Half Stage Axial Compressor with LES And PIV

The primary focus of this paper is to investigate how a rotor's unsteady tip clearance flow structure changes in a low speed one and half stage axial compressor when the rotor tip gap size is increased from 0.5 mm (0.49% of rotor tip blade chord, 2% of blade span) to 2.4 mm (2.34% chord, 4% span) at the design condition are investigated. The changes in unsteady tip clearance flow with the 0.62 % tip gap as the flow rate is reduced to near stall condition are also investigated. A Large Eddy Simulation (LES) is applied to calculate the unsteady flow field at these three flow conditions. Detailed Stereoscopic PIV (SPIV) measurements of the current flow fields were also performed at the Johns Hopkins University in a refractive index-matched test facility which renders the compressor blades and casing optically transparent. With this setup, the unsteady velocity field in the entire flow domain, including the flow inside the tip gap, can be measured. Unsteady tip clearance flow fields from LES are compared with the PIV measurements and both LES and PIV results are used to study changes in tip clearance flow structures. The current study shows that the tip clearance vortex is not a single structure as traditionally perceived. The tip clearance vortex is formed by multiple interlaced vorticities. Therefore, the tip clearance vortex is inherently unsteady. The multiple interlaced vortices never roll up to form a single structure. When phased-averaged, the tip clearance vortex appears as a single structure. When flow rate is reduced with the same tip gap, the tip clearance vortex rolls further upstream and the tip clearance vortex moves further radially inward and away from the suction side of the blade. When the tip gap size is increased at the design flow condition, the overall tip clearance vortex becomes stronger and it stays closer to the blade suction side and the vortex core extends all the way to the exit of the blade passage. Measured and calculated unsteady flow fields inside the tip gap agree fairly well. Instantaneous velocity vectors inside the tip gap from both the PIV and LES do show flow separation and reattachment at the entrance of tip gap as some earlier studies suggested. This area at the entrance of tip gap flow (the pressure side of the blade) is confined very close to the rotor tip section. With a small tip gap (0.5mm), the gap flow looks like a simple two-dimensional channel flow with larger velocity near the casing for both flow rates. A small area with a sharp velocity gradient is observed just above the rotor tip. This strong shear layer is turned radially inward when it collides with the incoming flow and forms the core structure of the tip clearance vortex. When tip gap size is increased to 2.4 mm at the design operation, the radial profile of the tip gap flow changes drastically. With the large tip gap, the gap flow looks like a two-dimensional channel flow only near the casing. Near the rotor top section, a bigger region with very large shear and reversed flow is observed.

compressor↗

Investigation of Unsteady Tip Clearance Flow in a Low-Speed One and Half Stage Axial Compressor with LES And PIV

The primary focus of this paper is to investigate how a rotor's unsteady tip clearance flow structure changes in a low speed one and half stage axial compressor when the rotor tip gap size is increased from 0.5 mm (0.49% of rotor tip blade chord, 2% of blade span) to 2.4 mm (2.34% chord, 4% span) at the design condition are investigated. The changes in unsteady tip clearance flow with the 0.62 % tip gap as the flow rate is reduced to near stall condition are also investigated. A Large Eddy Simulation (LES) is applied to calculate the unsteady flow field at these three flow conditions. Detailed Stereoscopic PIV (SPIV) measurements of the current flow fields were also performed at the Johns Hopkins University in a refractive index-matched test facility which renders the compressor blades and casing optically transparent. With this setup, the unsteady velocity field in the entire flow domain, including the flow inside the tip gap, can be measured. Unsteady tip clearance flow fields from LES are compared with the PIV measurements and both LES and PIV results are used to study changes in tip clearance flow structures. The current study shows that the tip clearance vortex is not a single structure as traditionally perceived. The tip clearance vortex is formed by multiple interlaced vorticities. Therefore, the tip clearance vortex is inherently unsteady. The multiple interlaced vortices never roll up to form a single structure. When phased-averaged, the tip clearance vortex appears as a single structure. When flow rate is reduced with the same tip gap, the tip clearance vortex rolls further upstream and the tip clearance vortex moves further radially inward and away from the suction side of the blade. When the tip gap size is increased at the design flow condition, the overall tip clearance vortex becomes stronger and it stays closer to the blade suction side and the vortex core extends all the way to the exit of the blade passage. Measured and calculated unsteady flow fields inside the tip gap agree fairly well. Instantaneous velocity vectors inside the tip gap from both the PIV and LES do show flow separation and reattachment at the entrance of tip gap as some earlier studies suggested. This area at the entrance of tip gap flow (the pressure side of the blade) is confined very close to the rotor tip section. With a small tip gap (0.5mm), the gap flow looks like a simple two-dimensional channel flow with larger velocity near the casing for both flow rates. A small area with a sharp velocity gradient is observed just above the rotor tip. This strong shear layer is turned radially inward when it collides with the incoming flow and forms the core structure of the tip clearance vortex. When tip gap size is increased to 2.4 mm at the design operation, the radial profile of the tip gap flow changes drastically. With the large tip gap, the gap flow looks like a two-dimensional channel flow only near the casing. Near the rotor top section, a bigger region with very large shear and reversed flow is observed.

unsteady↗

Dynamic behavior of liquid sheets

An experiment was conducted to study the aerodynamic instability of liquid sheets issuing from a two-dimensional air-assisted nozzle. Detailed measurements of the frequency of oscillation of the liquid sheet have been made. The measured vibrational frequencies were then correlated with the resulting spray angle. It was shown that the liquid sheet oscillations are dynamically similar to that of hard spring systems. For each air pressure, three distinct modes of breakup are distinguished. At low liquid flow rates, both the sinusoidal and the dilational modes are superimposed on the liquid sheet. With a further increase in liquid flow rate, the liquid sheet oscillations mainly become of the dilational type. It was also shown that the effect of introducing air in the nozzle is similar to the effect of inducing forced vibrations on the nozzle jaws. Thus, for each air flow rate, there is a specific vibration frequency for the nozzle. The frequency of these vibrations is proportional to the air velocity. As the liquid sheet natural frequency approaches that of the nozzle, resonance is established. At resonance, the maximum spray angle is achieved.

Mansour, Adel↗

Experimental Assessment of the Reciprocating Feed System

The primary goal of this project was to design, construct, and test a full scale, high pressure simulated propellant feed system test bed that could evaluate the ability of the Reciprocating Feed System (RFS) to provide essentially constant flow rates and pressures to a rocket engine. The two key issues addressed were the effects of the transition of the drain cycle from tank to tank and the benefits of other hardware such as accumulators to provide a constant pressure flow rate out of the RFS. The test bed provided 500 psi flow at rates of the order of those required for engines in the 20,000 lbf thrust class (e.g., 20 to 40 lb/sec). A control system was developed in conjunction with the test article and automated system operation was achieved. Pre-test planning and acceptance activities such as a documented procedure and hazard analysis were conducted and the operation of the test article was approved by, and conducted in coordination with, appropriate NASA Marshall Space Flight Center personnel under a Space Act Agreement. Tests demonstrated successful control of flow rates and pressures.

Eddleman, David E.↗

Screen Channel Liquid Acquisition Device Outflow Tests in Liquid Hydrogen

This paper presents experimental design and test results of the recently concluded 1-g inverted vertical outflow testing of two 325x2300 full scale liquid acquisition device (LAD) channels in liquid hydrogen (LH2). One of the channels had a perforated plate and internal cooling from a thermodynamic vent system (TVS) to enhance performance. The LADs were mounted in a tank to simulate 1-g outflow over a wide range of LH2 temperatures (20.3 - 24.2 K), pressures (100 - 350 kPa), and flow rates (0.010 - 0.055 kg/s). Results indicate that the breakdown point is dominated by liquid temperature, with a second order dependence on mass flow rate through the LAD. The best performance is always achieved in the coldest liquid states for both channels, consistent with bubble point theory. Higher flow rates cause the standard channel to break down relatively earlier than the TVS cooled channel. Both the internal TVS heat exchanger and subcooling the liquid in the propellant tank are shown to significantly improve LAD performance.

Liquid Acquisition Device↗

Choked flow of fluid nitrogen with emphasis on the thermodynamic critical region

Experimental measurements of critical flow rate and pressure ratio for nitrogen flowing through a nozzle are presented. Data for selected stagnation isotherms from 87.5 to 234 K with pressures to 9.3 MN/m2 are compared to an equilibrium model with real fluid properties and also a nonequilibrium model. Critical flow pressure ratio along an isotherm tends to peak while the flow rate indicates an inflection. The point is closely associated with the transposed critical temperature and represents a change in the fluid structure.

Hendricks, R. C.↗

Choked flow of fluid nitrogen with emphasis on the thermodynamic critical region.

Results of experimental measurements of critical flow rate and pressure ratio for nitrogen flowing through a nozzle. Data for selected stagnation isotherms from 87.5 to 234 K with pressures to 9.3 MN/sq m are compared to an equilibrium model with real fluid properties and also a nonequilibrium model. Critical flow pressure ratio along an isotherm tends to peak while the flow rate indicates an inflection. The point is closely associated with the transposed critical temperature and represents a change in the fluid structure.

Hendricks, R. C.↗

Experiments on transferring helium II with a thermomechanical pump

Porous plugs were tested as thermomechanical pumps for the transfer of helium II from one dewar to another through a transfer line and two valves. SEM images for two different size porous plugs are presented. The present pump design was shown to transfer 80 liters at a maximum flow rate of 90 liters/hr. Temperature, pressure, flow rate, and heater power data are compared with theoretical results.

Mills, G. L.↗

The flow of a thin liquid film on a stationary and rotating disk. I - Experimental analysis and flow visualization

The mean thickness of a thin liquid film of deionized water with a free surface on a stationary and rotating horizontal disk has been measured with a nonobtrusive capacitance technique. The measurements were taken when the rotational speed was 0-300 RPM and the flow rate was 7.0-15.0 LPM. A flow visualization study of the thin film was also performed to determine the characteristics of the waves on the free surface. When the disk was stationary, a circular hydraulic jump was present on the disk. Surface waves were found in the supercritical and subcritical regions at all flow rates studied. When the rotational speed of the disk is low, a standing wave at the edge of the disk was present. As the rotational speed increased, the surface waves changed from the wavy-laminar region to a region in which the waves ran nearly radially across the disk on top of a thin substrate of fluid.

Thomas, S.↗

Loop Heat Pipe Temperature Oscillation Induced by Gravity Assist and Reservoir Heating

The Laser Thermal Control System (LCTS) for the Advanced Topographic Laser Altimeter System (ATLAS) to be installed on NASA's Ice, Cloud, and Land Elevation Satellite (ICESat-2) consists of a constant conductance heat pipe and a loop heat pipe (LHP) with an associated radiator. During the recent thermal vacuum testing of the LTCS where the LHP condenser/radiator was placed in a vertical position above the evaporator and reservoir, it was found that the LHP reservoir control heater power requirement was much higher than the analytical model had predicted. Even with the control heater turned on continuously at its full power, the reservoir could not be maintained at its desired set point temperature. An investigation of the LHP behaviors found that the root cause of the problem was fluid flow and reservoir temperature oscillations, which led to persistent alternate forward and reversed flow along the liquid line and an imbalance between the vapor mass flow rate in the vapor line and liquid mass flow rate in the liquid line. The flow and temperature oscillations were caused by an interaction between gravity and reservoir heating, and were exacerbated by the large thermal mass of the instrument simulator which modulated the net heat load to the evaporator, and the vertical radiator/condenser which induced a variable gravitational pressure head. Furthermore, causes and effects of the contributing factors to flow and temperature oscillations intermingled.

Thermal Control Systems↗

Loop Heat Pipe Temperature Oscillation Induced by Gravity Assist and Reservoir Heating

The Laser Thermal Control System (LCTS) for the Advanced Topographic Laser Altimeter System (ATLAS) to be installed on NASA's Ice, Cloud, and Land Elevation Satellite (ICESat-2) consists of a constant conductance heat pipe and a loop heat pipe (LHP) with an associated radiator. During the recent thermal vacuum testing of the LTCS where the LHP condenser/radiator was placed in a vertical position above the evaporator and reservoir, it was found that the LHP reservoir control heater power requirement was much higher than the analytical model had predicted. Even with the control heater turned on continuously at its full power, the reservoir could not be maintained at its desired set point temperature. An investigation of the LHP behaviors found that the root cause of the problem was fluid flow and reservoir temperature oscillations, which led to persistent alternate forward and reversed flow along the liquid line and an imbalance between the vapor mass flow rate in the vapor line and liquid mass flow rate in the liquid line. The flow and temperature oscillations were caused by an interaction between gravity and reservoir heating, and were exacerbated by the large thermal mass of the instrument simulator which modulated the net heat load to the evaporator, and the vertical radiator/condenser which induced a variable gravitational pressure head. Furthermore, causes and effects of the contributing factors to flow and temperature oscillations intermingled.

0000↗

Investigation of Multiphase Flow in a Packed Bed Reactor Under Microgravity Conditions

In this paper we study the two-phase flow phenomena in a packed bed reactor using an integrated experimental and numerical method. The cylindrical bed is filled with uniformly sized spheres. In the experiment water and air are injected into the bed simultaneously. The pressure distribution along the bed will be measured. The numerical simulation is based on a two-phase flow solver which solves the Navier-Stokes equations on Cartesian grids. A novel coupled level set and moment of fluid method is used to construct the interface. A sequential method is used to position spheres in the cylinder. Preliminary experimental results showed that the tested flow rates resulted in pulse flow. The numerical simulation revealed that air bubbles could merge into larger bubbles and also could break up into smaller bubbles to pass through the pores in the bed. Preliminary results showed that flow passed through regions where the porosity is high. Comparison between the experimental and numerical results in terms of pressure distributions at different flow injection rates will be conducted. Comparison of flow phenomena under terrestrial gravity and microgravity will be made.

computational fluid dynamics↗

Numerical Simulation of Lean Blowout of Alternative Fuels in 7-element Lean Direct Injector

This research presents the result of numerically simulating 7-element swirl-venturi Lean Direct Injector (SV-LDI) lean blowout (LBO) experiments conducted at NASA Glenn Research Center in May of 2019. After simulating a cold flow case to confirm the pressure drop agrees well with the experiment, additional cases with two different fuels (an average jet fuel and a Gevo alcohol-to-jet fuel) from the National Jet Fuels Combustion Program (NJFCP) were computed to numerically determine the LBO condition. The procedure to approach the LBO follows the method used in the experiment where the air mass flow rate is gradually increased while the fuel supply is maintained. Transient history of global heat release rate as a function of air flow rate is presented, as well as temperature contours at different conditions to give a visual representation of the flame state. The Open National Combustion Code (OpenNCC) used in this research adopted reduced HyChem (Hybrid Chemistry) models along with k-LES turbulence model and a Lagrangian spray model that takes into account droplet internal temperature distribution affected by the shear force on the droplet surface. The transport equations and chemical reaction terms are integrated together to enhance conservation of chemical species that are especially important in the near LBO conditions. After showing the computed range of LBO agrees well with the experimental measurements, time averaged solutions of both fuels at their initial condition and at their limiting condition just before LBO are compared in detail to facilitate the understanding of LBO mechanism.

combustion↗