Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “continuous 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,477 records · Page 82

Excess Water in Astronaut Helmet During EVA on ISS: Mitigations with Flight Demonstrations

Following a second crew report of excess water inexplicably accumulating in the helmet during EVA-80 on March 23, 2022, NASA initiated an aggressive effort to identify, mitigate, and/or eliminate all sources of the potentially life-threatening water. Our narration highlights demonstrations of microgravity flow expectations using terrestrial scale models, mitigations to dangerous water migration within the helmet, low-g two-phase flow separations for the flow entering the helmet, and an investigation of the nature of liquid carry-over from the EMU condensing heat exchanger source. Fast-to-flight demonstrations of each aspect of the work are carried out during hands-on crew interaction with flight scale hardware on ISS during the 2022-2023 timeframe. The results of the tests are described with a focus on the rarely observed, and thus rarely studied, large length scale air-driven wall-bound droplet and rivulet twophase flows in microgravity. The success of the mitigations and directions for continued work is discussed in summary.

spacesuit↗

BETO 2021 Peer Review - Continuous Enzymatic Hydrolysis Development (CEHD) WBS 2.4.1.101

The Continuous Enzymatic Hydrolysis Development (CEHD) project aims to reduce the cost and commercialization scale-up risks of biorefinery sugar-lignin production through development of a deployable continuous enzymatic hydrolysis (CEH) process. Through the use of external cross flow membrane filtration loops coupled to enzymatic hydrolysis (EH) reactors, pretreated biomass solids and enzymes are retained for reaction while solubilized product sugars are removed in situ, with high extents of conversion achieved through a series of reactor-membrane unit stages. The CEHD project is focused on advancing CEH as a transformational, process-intensified, lower-cost method for producing soluble clarified biomass sugars and insoluble lignin-rich streams than traditional batch enzymatic hydrolysis (BEH). The project's primary objective is to reduce the cost of CEH to be compellingly lower than conventional BEH, 10% lower in year 1 and 20% lower in year 3 (end of project). A related objective is to expand CEH's operating envelope to increase process efficiency. Through more thorough de-risking and demonstration of CEH, in conjunction with building a suite of modeling and optimization tools that allow for more facile rigorous in silico evaluation of novel CEH designs and modalities, this project intends to elevate industry interest in adopting CEH as an improvement over conventional batch processing. This project was merit reviewed in FY20. It's now in its first year of a new 3-year plan spanning FY21-FY23.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Estimating Resolution Lengths of Hybrid Turbulence Models

A two-stage procedure has been devised for estimating the spatial resolution achievable in the simulation of a given flow on a given computational grid by a computational fluid dynamics (CFD) code that incorporates a hybrid model of turbulence. The hybrid models to which this procedure is especially relevant are those of the Reynolds-averaged Navier-Stokes (RANS) and the partial-averaged Navier-Stokes (PANS) approaches. This procedure represents the first step toward adding variable-resolution turbulence-modeling capabilities to CFD codes as part of a continuing effort to increase the accuracy and robustness of CFD simulations of unsteady flows. Some background information is prerequisite to a meaningful summary of the procedure. Among experts in CFD, it is well known that combination of the Reynolds-averaged Navier-Stokes (RANS) approach and eddy-viscosity turbulence models offers limited capability for simulating unsteady and complex flows. The RANS approach includes an assumption that most of the energy in a given flow is modeled through turbulence-transport equations and is resolved in a computational grid used to simulate the flow. RANS also overpredicts eddy viscosity, thereby yielding excessive damping of unsteady motion. The eddy viscosity attains an unphysically large value because of unresolved scales, and suppresses most temporal and spatial fluctuations in the resolved flow field. One approach used to overcome this deficiency is to provide a mechanism for the RANS equations to resolve motion only on the largest scales and to use a hybrid model to represent effects at smaller scales. The RANS approach involves the use of a standard two-equation turbulence model in which the effect of turbulence is summarized by a viscosity that is a function of (1) the time-averaged kinetic- energy density (k) associated with the local fluctuating (turbulent) component of flow and (2) the time-averaged rate of dissipation of the turbulent-kinetic- energy density ( ). In PANS, which was developed to overcome the grid dependency associated with other hybrid turbulence models (including that of RANS), the standard two-equation turbulence model is replaced by another two-equation model in which one solves for the previously unresolved k and , which are now allowed to vary in space and time.

Abdol-Hamid, Khaled S.↗

Grid generation about complex three-dimensional aircraft configurations

The problem of obtaining three dimensional grids with sufficient resolution to resolve all the flow or other physical features of interest is addressed. The generation of a computational grid involves a series of compromises to resolve several conflicting requirements. On one hand, one would like the grid to be fine enough and not too skewed to reduce the numerical errors and to adequately resolve the pertinent physical features of the flow field about the aircraft. On the other hand, the capabilities of present or even future supercomputers are finite and the number of mesh points must be limited to a reasonable number: one which is usually much less than desired for numerical accuracy. One technique to overcome this limitation is the 'zonal' grid approach. In this method, the overall field is subdivided into smaller zones or blocks in each of which an independent grid is generated with enough grid density to resolve the flow features in that zone. The zonal boundaries or interfaces require special boundary conditions such that the conservation properties of the governing equations are observed. Much work was done in 3-D zonal approaches with nonconservative zonal interfaces. A 3-D zonal conservative interfacing method that is efficient and easy to implement was developed during the past year. During the course of the work, it became apparent that it would be much more feasible to do the conservative interfacing with cell-centered finite volume codes instead of the originally planned finite difference codes. Accordingly, the CNS code was converted to finite volume form. This new version of the code is named CNSFV. The original multi-zonal interfacing capability of the CNS code was enhanced by generalizing the procedure to allow for completely arbitrarily shaped zones with no mesh continuity between the zones. While this zoning capability works well for most flow situations, it is, however, still nonconservative. The conservative interface algorithm was also implemented but was not completely validated.

Klopfer, Goetz H.↗

Application of the pressure sensitive paint technique to steady and unsteady flow

Pressure sensitive paint is a newly-developed optical measurement technique with which one can get a continuous pressure distribution in much shorter time and lower cost than a conventional pressure tap measurement. However, most of the current pressure sensitive paint applications are restricted to steady pressure measurement at high speeds because of the small signal-to-noise ratio at low speed and a slow response to pressure changes. In the present study, three phases of work have been completed to extend the application of the pressure sensitive paint technique to low-speed testing and to investigate the applicability of the paint technique to unsteady flow. First the measurement system using a commercially available PtOEP/GP-197 pressure sensitive paint was established and applied to impinging jet measurements. An in-situ calibration using only five pressure tap data points was applied and the results showed good repeatability and good agreement with conventional pressure tap measurements on the whole painted area. The overall measurement accuracy in these experiments was found to be within 0.1 psi. The pressure sensitive paint technique was then applied to low-speed wind tunnel tests using a 60 deg delta wing model with leading edge blowing slots. The technical problems encountered in low-speed testing were resolved by using a high grade CCD camera and applying corrections to improve the measurement accuracy. Even at 35 m/s, the paint data not only agreed well with conventional pressure tap measurements but also clearly showed the suction region generated by the leading edge vortices. The vortex breakdown was also detected at alpha=30 deg. It was found that a pressure difference of 0.2 psi was required for a quantitative pressure measurement in this experiment and that temperature control or a parallel temperature measurement is necessary if thermal uniformity does not hold on the model. Finally, the pressure sensitive paint was applied to a periodically changing pressure field with a 12.8s time period. A simple first-order pole model was applied to deal with the phase lag of the paint. The unsteady pressure estimated from the time-changing pressure sensitive paint data agreed well with the pressure transducer data in regions of higher pressure and showed the possibility of extending the technique to unsteady pressure measurements. However, the model still needs further refinement based on the physics of the oxygen diffusion into the paint layer and the oxygen quenching on the paint luminescence.

Shimbo, Y.↗

Numerical Calibration of Mass Flow Plug for Inlet Testing

A simple control volume model has been developed to calculate the discharge coefficient through a mass flow plug (MFP) and validated with a calibration experiment. The maximum error of the model within the operating region of the MFP is 0.54%. The control volume analysis developed work is comprised of a sequence of flow calculations through the MFP. The model uses the MFP geometry and operating pressure and temperature to couple continuity, momentum, energy, an equation of state, and wall shear. The discharge coefficient calculation also includes the effects of boundary layer growth, including the reduction in cross-sectional flow area as characterized by the boundary layer displacement thickness. The last calculation in the sequence uses an integral method to calculate the growth of the boundary layer, from which the displacement thickness is then determined. The result of these successive calculations is an accurate one-dimension model of the velocity, pressure, and temperature through the MFP. For comparison, a computational fluid dynamic (CFD) calibration is shown, which when compared to the presented numerical model, had a lower accuracy with a maximum error of 1.35% in addition to being slower by a factor of 100."

Flow Measurement↗

Joule Thomson refrigerator

A bi-directional Joule Thomson refrigerator is described, which is of simple construction at the cold end of the refrigerator. Compressed gas flowing in either direction through the Joule Thomson expander valve and becoming liquid, is captured in a container in direct continuous contact with the heat load. The Joule Thomson valve is responsive to the temperature of the working fluid near the valve, to vary the flow resistance through the valve so as to maintain a generally constant flow mass between the time that the refrigerator is first turned on and the fluid is warm, and the time when the refrigerator is near its coldest temperature and the fluid is cold. The valve is operated by differences in thermal coefficients of expansion of materials to squeeze and release a small tube which acts as the expander valve.

Chung K Chan↗

Analysis of Hydrology and Interim Measure Performance in the Chromium Plume at Los Alamos National Laboratory - 20498

Hexavalent chromium, Cr(VI), is present in the regional aquifer beneath Los Alamos National Laboratory (LANL) at concentrations greater than the New Mexico groundwater standard (50 ppb). The source of the Cr is blowdown from a power plant where potassium dichromate was used as a corrosion inhibitor from 1956 to 1972. Blowdown comingled with effluent from another source and migrated approximately 3 miles down a canyon, primarily as surface flow, until reaching an alluvial infiltration zone. The infiltrating water moved vertically through a thick (approximately 900 ft), complex stratigraphic sequence of unsaturated tuffs, basalt, and unconsolidated alluvial deposits. Horizontal transport also took place within spatially limited perched groundwater horizons in the vadose zone (VZ). Breakthrough at the regional aquifer water table is thought to have occurred in multiple hydraulic windows or 'drip points' located approximately 3 miles from the initial release site at the power plant. The regional aquifer is highly heterogeneous and anisotropic, and displays unconfined behavior in an upper zone and confined behavior at depth. The Cr(VI) contamination has so far remained confined to the upper approximately 25 m of the aquifer. The ambient hydraulic gradient in the regional aquifer is toward the east/southeast, but the water table is extremely flat in the Cr plume area compared to surrounding areas. There is some local variability in flow direction, potentially related to effects of seasonal pumping from nearby Los Alamos county water-supply wells, which draw from the regional aquifer several hundred feet below the water table. An interim measure (IM) is currently operating at the site to achieve hydraulic control of plume migration using a series of extraction wells, ion-exchange treatment, and injection along the downgradient portion of the plume. The migration of Cr contamination in the aquifer has been previously evaluated using a calibrated probabilistic numerical model, which has also been used to design the spatial configuration of extraction and injection wells for the IM, and to identify locations for monitoring wells. The model is built using the Finite Element Heat and Mass (FEHM) transfer code, and the Model Analysis and Decision Support (MADS) software is used for analyses. To date, the portion of the IM system that has operated in a continuous manner has been very effective at hydraulic control. Modifications to the flow field associated with the influence of IM pumping and injection are evaluated here, with modeling constraints to be provided by recent observations of tracers injected into both monitoring wells and IM injection wells, as well as opportunistic observations of geochemical signatures of treated ion-exchange effluent injected into IM injection wells. The tracer and geochemical data provide insights into modified flow directions within the Cr plume area that cannot be otherwise obtained. Results of the modeling work are being used to inform adaptive management of IM operations and will be a valuable tool for development and refinement of long-term remediation strategies. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Two-phase measurements of a spray in the wake of a bluff body

The dynamics of spray drop interaction with the turbulent, recirculating wake of a flat disk bluff body were investigated using a phase Doppler particle analyzer to determine drop size and velocity and the gas-phase velocity. Detailed measurements obtained included spray drop size, axial and radial velocity, angle of trajectory, and size-velocity correlations. The gas-phase velocity was determined from seeding of the two-phase flow. Results showed dramatic differences in drop behavior for various size classes when interacting with the turbulent flow field. Small drops were quickly entrained and recirculated, while initially, the larger drops continued in the general direction of the spray cone. Further downstream, significant numbers of large drops recirculated, generating a bifurcated size-velocity correlation. These lateral convections and streamwise accelerations and decelerations strongly influenced the number density along with size and velocity distributions. The complex interaction of the spray with the turbulent air-flow points out the need for spatially-resolved measurements that determine drop behavior for individual size classes, rather than characterizing a spray only via simple integral quantities such as the Sauter mean diameter.

Rudoff, R. C.↗

Finite element methods of analysis for high speed viscous flows

Over the past three years a finite element based procedure for the solution of high speed viscous compressible flows was developed. The approach followed was to compute steady state solutions via a false transient, using an explicit time stepping scheme, and to attempt to improve the solution quality by incorporating adaptive mesh procedures. The main thrust of the work was to continue on the extension of the approach to the solution of some realistic compressible viscous flows. When flows at high Reynolds number are investigated, is soon becomes apparent that explicit techniwues have to be supplemented if they are to deal effectively with the large variations in element size and aspect ratio which characterize the computational grids necessary for adequate resolutin of the primary flow features. For this reason, the Taylor-Galerkin solution algorithm was rewritten in an explicit/implicit form. Solutions were computed for the problems of a flow past a flat plate, M-3, Re-100; shock/boundary layer interaction, M-2, Re-296000; flow over a compression corner, M-11.68, Re246000; and unifrom flow past a circular cylinder, M-6.34, Re-39770. A summary of the results is included and demonstrates the numerical performance of the scheme.

Source record↗

Finite-Difference Solutions of the Alternate Turbopump Development High-Pressure Oxidizer Turbopump Pump-End Ball-Bearing Cavity Flows

These analyses were undertaken to aid in the understanding of flow phenomena in the Alternate Turbopump Development (ATD) High-pressure Oxidizer Turbopump (HPOTP) Pump-end ball bearing (PEBB) cavities and their roles in turbopump vibration initiation and bearing distress. This effort was being performed to provide timely support to the program in a decision as to whether or not the program should be continued. In the first case, it was determined that a change in bearing through flow had no significant effect on axial preload. This was a follow-on to a previous study which had resulted in a redesign of the bearing exit cavity which virtually eliminated bearing axial loading. In the second case, a three-dimensional analysis of the inner-race-guided cage configuration was performed so as to determine the pressure distribution on the outer race when the shaft is 0.0002 inches off-center. The results indicate that there is virtually no circumferential pressure difference caused by the offset to contribute to bearing tilt. In the third case, axisymmetric analyses were performed on an outer-race guided cage configuration to determine the magnitude of tangential flow entering the bearing. The removed-shoulder case was analyzed as was the static diverter case. A third analysis where the preload spring was shielded by a sheet of metal for the baseline case was also performed. It was determined that the swirl entering the bearing was acceptable and the project decided to use the outer-race-guided cage configuration. In the fourth case, more bearing configurations were analyzed. These analyses included thermal modeling so as to determine the added benefit of injecting colder fluid directly onto the bearing inner-race contact area. The results of these analyses contributed to a programmatic decision to include coolant injection in the design.

Benjamin, Theodore G.↗

Remobilizing the Interfaces of Thermocapillary Driven Bubbles Retarded by the Adsorption of a Surfactant Impurity on the Bubble Surface

Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is applied to the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow which propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillarity to direct the movement of bubbles in space is the fact that surfactant impurities which are unavoidably present in the continuous phase can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature tension gradient, and decreases to near zero the thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have been retarded by the adsorption onto the bubble surface of a surfactant impurity, Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition the rapid kinetic exchange and high bulk concentration maintain a saturated surface with a uniform surface concentrations. This prevents retarding surface tension gradients and keeps the velocity high. In our first report last year, we detailed experimental results which verified the theory of remobilization in ground based experiments in which the steady velocity of rising bubbles was measured in a continuous phase consisting of a glycerol/water mixture containing a polyethylene glycol surfactant C12E6 (CH3(CH2)11(OCH2CH2)6OH). In our report this year, we detail our efforts to describe theoretically the remobilization observed. We construct a model in which a bubble rises steadily by buoyancy in a continuous (Newtonian) viscous fluid containing surfactant with a uniform far field bulk concentration. We account for the effects of inertia as well as viscosity in the flow in the continuous phase caused by the bubble motion (order one Reynolds number), and we assume that the bubble shape remains spherical (viscous and inertial forces are smaller than capillary forces, i e. small Weber and capillary numbers). The surfactant distribution is calculated by solving the mass transfer equations including convection and diffusion in the bulk, and finite kinetic exchange the bulk and the surface. Convective effects dominate diffusive mass transfer in the bulk of the liquid (high Peclet numbers) except in a thin boundary layer near the surface. A finite volume method is used to numerically solve the hydrodynamic and mass transfer equations on a staggered grid which accounts specifically for the thin boundary layer. We present the results of the nondimensional drag as a function of the bulk concentration of surfactant for different rates of kinetic exchange, from which we develop criteria for the concentration necessary to develop a prescribed degree of remobilization. The criteria compare favorably with the experimental results.

Palaparthi, Ravi↗

Transport phenomena of reactive fluid flow in heterogeneous combustion processes.

A previously developed computer program was used to model two transient hybrid combustion processes involving tubes of solid Plexiglas. In the first study, representing combustion of a hybrid rocket, the oxidizing gas was oxygen, and calculations were continued sufficiently long to obtain steady-state values. Systematic variations were made in reaction rate constant, mass flow rate, and pressure, alternatively using constant and temperature dependent regression rate models for the fuel surface. Consistent results were obtained, as is evidenced by the values for the mass function of the reaction product and the flame temperature, for which plots are supplied. In the second study, fire initiation in a duct was studied, with an air mixture as the oxidizing gas. It was demonstrated that a satisfactory flame spread mechanism could be reproduced on the computer. In both of the above applications, the general, transient, two-dimensional conservation equations were represented, together with chemical reactions, solid-fuel interface conditions, and heat conduction in the solid fuel.

Hung, W. S. Y.↗

A method for calculating unsteady turbulent boundary layers in two- and three-dimensional flows.

The governing equations for an unsteady turbulent boundary layer on a swept infinite cylinder, composed of a continuity equation, a pair of momentum equations and a pair of rate equations for the shear stress, based on the turbulent kinetic-energy equation, are solved numerically. These rate equations for the shear-stress express a balance between the convection, production, dissipation and diffusion of the turbulent shear stress components. Due to the physical model of turbulence employed, this system of equations is hyperbolic. An explicit, second-order accurate, conditionally stable difference scheme for this system, together with appropriate boundary and initial conditions, is formulated and developed into a computer program. Calculations for oscillating free stream flows with no pressure gradient show significant unsteady effects on the turbulent boundary layer.

Singleton, R. E.↗

A continuous physiological data collector

COP-DAC system utilizes oxygen and carbon dioxide analyzers, gas-flow meter, gas breathe-through system, analog computer, and data storage system to provide actual rather than average measurements of physiological and metabolic functions.

Bush, J. C.↗

A Methodology for Actuating RCS Jets in a Continuous, Time-Accurate CFD Simulation

Computational fluid dynamics simulations of Mars entry vehicle reaction control systems have been used to gain a better understanding of the wide range of flow phenomena encountered during various phases of flight. One aspect of this is simulating the system during controlled flight. Current practices model the system in a continuously actuated state, but this methodology is unable to accurately model the dynamic behavior of the system during EDL. One method to improve this is to create a coupled CFD-RBD-Control simulation. The first step to this is establishing a method of accurately actuating an RCS jet during a time accurate CFD simulation. This paper presents a method to achieve this actuation. There are three phases to the methodology: establishing the off state, creating a ramping function that will initialize the jet flow, and establishing the on state. A test case of a jet in a supersonic crossflow is implemented to establish the success of the methodology in being able to actuate a jet in the midst of an ongoing CFD simulation. It shows significant improvement in boundary condition responsiveness for linear and exponential ramping functions in comparison to the step function. A preliminary demonstration of the methodology is presented for the MSL vehicle which shows significant advantages for the implementation of a linear ramping function over a step function.

Alexandra Hickey↗

A Methodology for Actuating RCS Jets in a Continuous, Time-Accurate CFD Simulation

Computational fluid dynamics simulations of Mars entry vehicle reaction control systems have been used to gain a better understanding of the wide range of flow phenomena encountered during various phases of flight. One aspect of this is simulating the system during controlled flight. Current practices model the system in a continuously actuated state, but this methodology is unable to accurately model the dynamic behavior of the system during EDL. One method to improve this is to create a coupled CFD-RBD-Control simulation. The first step to this is establishing a method of accurately actuating an RCS jet during a time accurate CFD simulation. This paper presents a method to achieve this actuation. There are three phases to the methodology: establishing the off state, creating a ramping function that will initialize the jet flow, and establishing the on state. A test case of a jet in a supersonic crossflow is implemented to establish the success of the methodology in being able to actuate a jet in the midst of an ongoing CFD simulation. It shows significant improvement in boundary condition responsiveness for linear and exponential ramping functions in comparison to the step function. A preliminary demonstration of the methodology is presented for the MSL vehicle which shows significant advantages for the implementation of a linear ramping function over a step function.

Alexandra Hickey↗

LSA Large Area Silicon Sheet Task Continuous Czochralski Process Development

A commercial Czochralski crystal growing furnace was converted to a continuous growth facility by installation of a small, in-situ premelter with attendant silicon storage and transport mechanisms. Using a vertical, cylindrical graphite heater containing a small fused quartz test tube linear from which the molten silicon flowed out the bottom, approximately 83 cm of nominal 5 cm diamter crystal was grown with continuous melt addition furnished by the test tube premelter. High perfection crystal was not obtained, however, due primarily to particulate contamination of the melt. A major contributor to the particulate problem was severe silicon oxide buildup on the premelter which would ultimately drop into the primary melt. Elimination of this oxide buildup will require extensive study and experimentation and the ultimate success of continuous Czochralski depends on a successful solution to this problem. Economically, the continuous Czochralski meets near-term cost goals for silicon sheet material.

Rea, S. N.↗