Search NASA⌕ Search

SEARCH · Search NASA

Results for “Two Phase 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 559 records · Page 31

Spectrally Stabilized Interface Capturing Formulation and Implementation in Nek5000/NekRS

This report documents the formulation of a novel level-set method for incompressible two-phase flows in the continuous Galerkin (CG) high order spectral element framework. The overall method hinges on a novel implementation of the spectral vanishing viscosity (SVV) operator for the stabilization of linear/non-linear hyperbolic problems. The multidimensional SVV convolution kernels, which in essence, have a similar effect as a high pass filter applied to the derivatives, are formulated by exploiting the tensor product form, analogous to the construction of the usual stiffness matrix system. The resulting kernels are directionally decoupled and ensure a linear, symmetric positive definite, elliptic matrix operator. The SVV formulation is demonstrated to provide a robust stabilizing mechanism through challenging linear and non-linear hyperbolic problems, including problems pertinent to the level-set formulation. The two-phase framework conceptualized herein is based on the conservative level-set (CLS) method which represents the interface between the fluids by the 0.5 iso-contour of the smoothed Heaviside function. The CLS method is augmented with a preconditioning procedure for interface normals using the signed distance function which precludes the manifestation of spurious oscillations in the vicinty of the interface. Further, the existing mixed explicit-implicit approach for the solution of Navier-Stokes equations in Nek5000, as described in Tomboulides et al, is augmented with a pressure coefficient splitting approach for the Poisson equation, which greatly accelerated the convergence of pressure solver for two-phase systems with large density ratio. The robustness and accuracy of the overall two-phase method is demonstrated through canonical challenging problems involving high density and viscosity ratios, with and without surface tension. The two-phase formulation is wholly implemented in Nek5000 and the SVV stabilization method is implemented in NekRS, which is the essential precursor to the two-phase framework, undergoing active development.

97 MATHEMATICS AND COMPUTING↗

EWTN: Quantifying Mass Transport to Enable Water Electrolyzer Architectures with Low Flow-Rate Sensitivity

To develop cost-effective and high-performing polymer electrolyte water electrolyzers (PEWEs) for gigawatt-scale applications, researchers have focused on reducing precious metal catalyst loadings and optimizing porous transport layers. However, the performance of PEWEs is also affected by system architecture. Mass transport losses are dependent on localized architecture and material interactions. In-situ measurements, such as current density distribution maps have demonstrated advantages in understanding the intricate characteristics and influence of two-phase flow within PEWEs. This study proposes the parameter of effective water transport number (EWTN) as a quantitative tool to investigate such current density distribution (CDD) measurements for PEWEs. Results show that higher flow-rates have EWTN values of 0.95 and above, indicating no mass transport limitations; while lower flow-rates with large gradients CDD have EWTN values between 0.6-0.8, indicating mass transport limited conditions. The new analysis also identified a correlation between mass transport losses due to bubble accumulation, membrane hydration, and ohmic overpotentials. To address these limitations, an unitized pin-type LGDL/flow-field design was developed, which effectively prevents local gas phase accumulation, resulting in improved mass transport characteristics. The results of this work show reduced flow-rate sensitivity with the pin-type architecture and ∼13% increased performance at 0.24 ml / min / cm 2 .

Electrochemistry↗

Bellows flow-induced vibrations

Results of theoretical and experimental investigations of bellows typical of those found in space shuttle external tanks are presented. Correlation parameters are identified which generalize the alternating stress calculations cited in an earlier study. Alternating stress amplitudes and mean stress levels form the basis of a fatigue analysis incorporating seven ordinate charts for 347 stainless steel, alloy 21 6-9, and Inco 718. A crack propagation model is included with a program for computing bellows fatigue life. Two phase flow and material hardness properties are discussed.

Johnson, J. E.↗

SRM Internal Flow Tests and Computational Fluid Dynamic Analysis: CFD RSRM Full-Scale Analyses - Volume 2

This document presents the full-scale analyses of the CFD RSRM. The RSRM model was developed with a 20 second burn time. The following are presented as part of the full-scale analyses: (1) RSRM embedded inclusion analysis; (2) RSRM igniter nozzle design analysis; (3) Nozzle Joint 4 erosion anomaly; (4) RSRM full motor port slag accumulation analysis; (5) RSRM motor analysis of two-phase flow in the aft segment/submerged nozzle region; (6) Completion of 3-D Analysis of the hot air nozzle manifold; (7) Bates Motor distributed combustion test case; and (8) Three Dimensional Polysulfide Bump Analysis.

Source record↗

New findings and instrumentation from the NASA Lewis Microgravity Facility

The study of fundamental combustion and fluid physics phenomena in microgravity enables a new range of experiments to be performed since buoyancy-type flows are nearly eliminated, normally obscured forces and flows may be isolated, gravitational settling or sedimentation is nearly eliminated, and larger time or length scales in experiments become permissible. Some new phenomena found through ground-based microgravity research at NASA facilities is reviewed along with the instrumentation and drop tower facilities used in the research. Experiments are described in the areas of premixed gas combustion, droplet combustion, particle cloud combustion, liquid reorientation, and two-phase flow patterns. Future plans for instrumentation development are also outlined.

Ross, Howard D.↗

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Flow Boiling Heat Transfer Characteristics of Water for Metal-foam-filled Horizontal Tubes

The goal of this paper is to enhance the in-tube flow boiling heat transfer of water while accounting for potential increase in pressure drop. Metal foams are a class of cellular structure with large surface-area-to-volume ratio and tortuous structure which has shown promising results for various energy conversion and storage applications. The higher heat transfer area and higher nucleation sites density due to the porous media can effectively enhance the heat transfer of water in both single-phase and two-phase operations. This paper presents the local heat transfer and pressure drop measurements of water in a partially metal-foam-filled horizontal copper tube. The ranges of parameters in the experiments are: mass flux from 80 to 200 kg/s-m2, heat flux from 3.5 to 105 kW/m2, and vapor quality from subcooled to 0.3. The results in the metal-foam-filled tube are compared with that of the bare copper tube. The effect of metal foam on the thermal hydraulic characteristics is discussed and analyzed. In addition, the two-phase flow behavior in a transparent tube filled with the same metal foam was visualized and investigated by high-speed imaging system. The experiments were also conducted in the tube without the metal foam for comparison.

Nawaz, Kashif↗

EXPERIMENTAL AND MODELLING ANALYSIS OF A LARGE-SCALE TWO-PHASE LOOP THERMOSYPHON

Liquid pumping requires large quantities of electrical energy, including about 7% of the energy of building heating, ventilation, and air conditioning (HVAC) systems. To replace pumped condenser-cooling tower water loops with a passive alternative system, we implemented a commercial-scale two-phase loop thermosyphon (TPLT). The unit consists of a 13 m riser integrated with a commercially available cooling tower and circulation water heater that simulates heat loads up to 25 kW. In addition to providing passive cooling capabilities, the cooling tower unit is also maintenance free, reliable, and can operate in both dry and wet modes. This study characterizes the performance (temperature difference between the evaporator and condenser) and the two-phase flow behavior of the loop under various refrigerant charges. Sight glasses installed throughout the loop are used to identify the operating flow regimes in the riser and downcomer. Over the range of operating conditions, we identified that there is an optimal refrigerant charge range for a specific heat load at which near-isothermal operation can be achieved. We further developed a model to predict the operating flow rate and gravitational height in the downcomer and compared it with the experimental data. The results show that the model agrees with the experimental data, in particular the threshold gravitational head height at which there will be subcooled liquid at the condenser exit, indicating that it can be used as a predictive tool for identifying the optimal loop charge for a given heat load.

Loop thermosyphon, passive two phase↗

FY23 Report on Water NSTF Testing at Two-Phase Conditions: Off-normal Scenarios

Under support from the Department of Energy (DOE) and the Office of Advanced Reactor Technologies (ART), a large-scale test facility has been constructed at Argonne National Laboratory to generate NQA-1 qualified validation data for passive decay heat removal systems in advanced reactors. The Natural convection Shutdown heat removal Test Facility (NSTF) reflects key features of a ½ scale, water-based, Reactor Cavity Cooling System (RCCS) and is intended to study the behavior, bound performance, and ultimately guide design decisions for passive decay heat removal systems for advanced reactors. In addition to the experimental activities detailed in this report, a supportive computational modeling effort is on-going which has been demonstrated to significantly strengthen the experimental program while also improving accuracy of the computer models. Together these create a mutually beneficial relationship integral to meeting the overall program objective of examining the heat removal performance of the RCCS concept. This report details the experimental activities and upgrades performed during the program’s fifth year of water-based operation. The theme of this year of testing was examining off-normal scenarios. In practice, that meant examining scenarios with blockages or other design basis conditions. From a maintenance and capability standpoint that meant assessing valving options for implementing blockages, assessing drain and refill capabilities, and ensuring in all cases that operators could safely, efficiently, and repeatably bring the facility to an off-normal condition and return it to the normal operating state. In one notable case, no solution existed, nor was a solution commercially available, that met the facility needs for implementing blockages within the two-phase region of the chimney. A valve and feedthrough were custom-designed for NSTF to operate submerged in saturated water while being safely actuated from outside the tank. The fifth year of testing included eight matrix tests consisting of 162 hours of active heating, 9,735 kWh of electrical heating, with seven tests classified as Accepted per NQA-1 and one as Trending. All matrix tests were performed at two-phase flow conditions. The power parametric series from previous years was extended to include another scenario with a decay load equivalent to 1.75 MWt, full scale, further resolving the system flow oscillation response as a function of input power. A depletion test was performed to extend previous work and examine the system’s approach to stagnation and geysering that were not previously observed. This test began at an initial inventory level of 50%, prototypic power of 2.1 MWt, and employed an accelerated drain of 0.75 gpm plus boiloff until stagnation conditions were achieved and multiple geysering events were recorded. Further exploring the geysering phenomenon was a separate effects test examining static boiling in the risers and chimney. This static boiling test was performed at 48 kWe, a higher power than testing in the previous year, and included enhancements to the drain system to remove inventory in a better-controlled and more reliable manner. Two tests were performed, at prototypic conditions of 1.4 MWt and 2.4 MWt, respectively, examining the system response to throttled conditions in the singlephase region of the piping. A ball valve immediately before he riser inlet header was increasingly throttled, progressing the system through a repeatable pattern of two-phase oscillations and stable regimes before ultimately inducing stagnation and geysering. Concluding the year, two tests were performed at baseline conditions of 70% initial inventory level and 2.1 MWt of prototypic decay heat to examine both repeatability and the response to throttling within the two-phase region, at the inventory tank inlet using the custom valve mentioned above. The first such test repeated transient operating conditions and revealed significant sensitivity t

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Influence of Control Jets on Flush Air-data Sensors

Computations are performed to investigate the effect of rocket control motors on flush air-data sensor systems. Such sensors are critical for the control of space vehicles during launch and re-entry, but are prone to interference from rocket motors, hypersonic-flow effects, etc. Computational analyses provide a means for studying these interference effects and exploring opportunities for mitigating them, either through design techniques or through appropriate processing of the sensor outputs. In the present work, the influence of rocket control motors on the nosecone flush air-data sensors of a launch-abort vehicle is studied. Particular attention is paid to the differential effect of various control-jet combinations on surface pressures. The relative effectiveness of inviscid, viscous, turbulent and two-phase-flow approximations in addressing this problem is also investigated.

Woodruff, Stephen↗

The Thermal Response of a Packed Bed Thermal Energy Storage System upon Saturated Steam Injection Using Distributed Temperature Sensing

The effectiveness of a thermal energy storage (TES) system is typically characterized with the help of thermal stratification or temperature gradients along the direction of heat injection, which is typically the flow direction of heat transfer fluid. The steepness of temperature gradients are a direct indicator of the effectiveness or efficiency of the heat storage or dispatch process. The temperature gradient evolution along the packed bed of ceramic particles upon saturated steam injection is presented in this work. Distributed temperature sensing based on optical frequency domain reflectometry was deployed in a packed bed of ceramic particles to capture the thermal front evolution in the axial direction. The physical processes accompanying steam injection in packed beds are complex due to phase change, transitioning two-phase flow, and changes in condensate accumulation. Therefore, the variation of thermal response of the TES system for various steam injection flow rates was experimentally studied using a high-resolution distributed temperature sensing system in a chemically inert alumina particle-packed bed. Distinct zones of different heat transfer modes were observed during the steam injection experiments. A distinct conduction zone, evident from diffuse thermal fronts, was observed at low flow rates, and these thermal gradients became sharper as the flow rate increased. The diffuse thermal fronts in the heat storage media suggest a low exergy efficiency of the TES system, as energy losses started initiating before a significant fraction of the bed was saturated with steam.

25 ENERGY STORAGE↗

Determination of the Residence Time of Food Particles During Aseptic Sterilization

The paper describes a non-invasive method to measure the time an individual particle takes to move through a length of stainless steel pipe. The food product is in two phase flow (liquids and solids) and passes through a pipe with pressures of approximately 60 psig and temperatures of 270-285 F. The proposed problem solution is based on the detection of transitory amplitude and/or phase changes in a microwave transmission path caused by the passage of the particles of interest. The particles are enhanced in some way, as will be discussed later, such that they will provide transitory changes that are distinctive enough not to be mistaken for normal variations in the received signal (caused by the non-homogeneous nature of the medium). Two detectors (transmission paths across the pipe) will be required and place at a known separation. A minimum transit time calculation is made from which the maximum velocity can be determined. This provides the minimum residence time. Also average velocity and statistical variations can be computed so that the amount of 'over-cooking' can be determined.

Carl, J. R.↗

NASA Tech Briefs, September 2003

Topics include: Oxygen-Partial-Pressure Sensor for Aircraft Oxygen Mask; Three-Dimensional Venturi Sensor for Measuring Extreme Winds; Swarms of Micron-Sized Sensors; Monitoring Volcanoes by Use of Air-Dropped Sensor Packages; Capacitive Sensors for Measuring Masses of Cryogenic Fluids; UHF Microstrip Antenna Array for Synthetic- Aperture Radar; Multimode Broad-Band Patch Antennas; 164-GHz MMIC HEMT Frequency Doubler; GPS Position and Heading Circuitry for Ships; Software for Managing Parametric Studies; Software Aids Visualization of Computed Unsteady Flow; Software for Testing Electroactive Structural Components; Advanced Software for Analysis of High-Speed Rolling-Element Bearings; Web Program for Development of GUIs for Cluster Computers; XML-Based Generator of C++ Code for Integration With GUIs; Oxide Protective Coats for Ir/Re Rocket Combustion Chambers; Simplified Waterproofing of Aerogels; Improved Thermal-Insulation Systems for Low Temperatures; Device for Automated Cutting and Transfer of Plant Shoots; Extension of Liouville Formalism to Postinstability Dynamics; Advances in Thrust-Based Emergency Control of an Airplane; Ultrasonic/Sonic Mechanisms for Drilling and Coring; Exercise Device Would Exert Selectable Constant Resistance; Improved Apparatus for Measuring Distance Between Axles; Six Classes of Diffraction-Based Optoelectronic Instruments; Modernizing Fortran 77 Legacy Codes; Active State Model for Autonomous Systems; Shields for Enhanced Protection Against High-Speed Debris; Scaling of Two-Phase Flows to Partial-Earth Gravity; Neutral-Axis Springs for Thin-Wall Integral Boom Hinges.

Source record↗

FY-25 Progress on Computational Modeling of the Water Based NSTF

This report summarizes the system level modeling using RELAP5-3D of the Natural Convection Shutdown Heat Removal Test Facility (NSTF) completed in FY25. This year’s work focuses on a new tank configuration where the inlet of the tank was lowered in elevation by 45”. The stability boundaries of the NSTF are thoroughly studied and stability maps are constructed based on the stability and the oscillation patterns of the system. Five distinct operational modes are identified, namely single-phase liquid, uniform double peak oscillations, uniform sinusoidal oscillations, stable two-phase flow, and non-uniform oscillations. Next, the riser inlet throttling case of experimental test Run-104 is simulated with the RELAP5 model where good agreement is obtained between the model and the experimental data. The simulation also highlights the effects of backflow of water from the tank to the upper region of the chimney. Additionally, the decay heat removal test of Run-99 is simulated with the RELAP5 model. Comparison is carried out between this run and a similar run with the mid-tank inlet of Run-74 performed in FY22. With the lower tank inlet, the RELAP5 model is able to predict the experimental data more accurately than the previous mid tank inlet configuration. The discrepancy in model prediction accuracy highlights the non-symmetrical spatial effects in the tank that would otherwise be more easily captured with higher fidelity models. Lastly, two exploratory studies are conducted to investigate the behaviors of the NSTF when 1) heating is provided to the downcomer and 2) a bypass channel is added between the horizontal chimney section to the downcomer.

42 ENGINEERING↗

Probing multiscale dissolution dynamics in natural rocks through microfluidics and compositional analysis

Mineral dissolution significantly impacts many geological systems. Carbon released by diagenesis, carbon sequestration, and acid injection are examples where geochemical reactions, fluid flow, and solute transport are strongly coupled. The complexity in these systems involves interplay between various mechanisms that operate at timescales ranging from microseconds to years. Current experimental techniques characterize dissolution processes using static images that are acquired with long measurement times and/or low spatial resolution. These limitations prevent direct observation of how dissolution reactions progress within an intact rock with spatially heterogeneous mineralogy and morphology. We utilize microfluidic cells embedded with thin rock samples to visualize dissolution with significant temporal resolution (100 ms) in a large observation window (3 × 3 mm). Here we injected acidic fluid into eight shale samples ranging from 8 to 86 wt % carbonate. The pre- and postreaction microstructures are characterized at the scale of pores (0.1 to 1 µm) and fractures (1 to 1,000 µm). We observe that nonreactive particle exposure, fracture morphology, and loss of rock strength are strongly dependent on both the relative volume of reactive grains and their distribution. Time-resolved images of the rock unveil the spatiotemporal dynamics of dissolution, including two-phase flow effects in real time and illustrate the changes in the fracture interface across the range of compositions. Moreover, the dynamical data provide an approach for characterizing reactivity parameters of natural heterogeneous samples when porous media effects are not negligible. The platform and workflow provide real-time characterization of geochemical reactions and inform various subsurface engineering processes.

58 GEOSCIENCES↗

Multi-Zone Liquid Thrust Chamber Performance Code with Domain Decomposition for Parallel Processing

Computational Fluid Dynamics (CFD) has considerably evolved in the last decade. There are many computer programs that can perform computations on viscous internal or external flows with chemical reactions. CFD has become a commonly used tool in the design and analysis of gas turbines, ramjet combustors, turbo-machinery, inlet ducts, rocket engines, jet interaction, missile, and ramjet nozzles. One of the problems of interest to NASA has always been the performance prediction for rocket and air-breathing engines. Due to the complexity of flow in these engines it is necessary to resolve the flowfield into a fine mesh to capture quantities like turbulence and heat transfer. However, calculation on a high-resolution grid is associated with a prohibitively increasing computational time that can downgrade the value of the CFD for practical engineering calculations. The Liquid Thrust Chamber Performance (LTCP) code was developed for NASA/MSFC (Marshall Space Flight Center) to perform liquid rocket engine performance calculations. This code is a 2D/axisymmetric full Navier-Stokes (NS) solver with fully coupled finite rate chemistry and Eulerian treatment of liquid fuel and/or oxidizer droplets. One of the advantages of this code has been the resemblance of its input file to the JANNAF (Joint Army Navy NASA Air Force Interagency Propulsion Committee) standard TDK code, and its automatic grid generation for JANNAF defined combustion chamber wall geometry. These options minimize the learning effort for TDK users, and make the code a good candidate for performing engineering calculations. Although the LTCP code was developed for liquid rocket engines, it is a general-purpose code and has been used for solving many engineering problems. However, the single zone formulation of the LTCP has limited the code to be applicable to problems with complex geometry. Furthermore, the computational time becomes prohibitively large for high-resolution problems with chemistry, two-equation turbulence model, and two-phase flow. To overcome these limitations, the LTCP code is rewritten to include the multi-zone capability with domain decomposition that makes it suitable for parallel processing, i.e., enabling the code to run every zone or sub-domain on a separate processor. This can reduce the run time by a factor of 6 to 8, depending on the problem.

Homayun K. Navaz↗

The electrospray: Fundamentals and combustion applications

Liquid fuel dispersion in practical systems is typically achieved by spraying the fuel into a polydisperse distribution of droplets evaporating and burning in a turbulent gaseous environment. In view of the nearly unsurmountable difficulties of this two-phase flow, it would be useful to use an experimental arrangement that allow a systematic study of spray evolution and burning in configurations of gradually increasing levels of complexity, starting from laminar sprays to fully turbulent ones. An Electrostatic Spray (ES) of charged droplets lends itself to this type of combustion experiments under well-defined conditions and can be used to synthesize gradually more complex spray environments. In its simplest configuration, a liquid is fed into a small metal tube maintained at several kilovolts relative to a ground electrode few centimeters away. Under the action of the electric field, the liquid meniscus at the outlet of the capillary takes a conical shape, with a thin jet emerging from the cone tip. This jet breaks up farther downstream into a fine spray of charged droplets. Several advantages distinguish the electrospray from alternative atomization techniques: the self-dispersion property of the spray due to coulombic repulsion; the absence of droplet coalescence; the potential control of the trajectories of charged droplets by suitable disposition of electrostatic fields; and the decoupling of atomization, which is strictly electrostatic, from gas flow processes. Furthermore, as recently shown in our laboratory, the electrospray can produce quasi-monodisperse droplets over a very broad size range (1-100 microns). The ultimate objective of this research project is to study the formation and burning of electrosprays of liquid fuels first in laminar regimes and then in turbulent ones. Combustion will eventually be investigated in conditions of three-dimensional droplet-droplet interaction, for which experimental studies have been limited to either qualitative observations in sprays or more quantitative observations on simplified systems consisting of a small number of droplets or droplet arrays. The compactness and potential controllability of this spray generaiton system makes it appealing for studies to be undertaken in the next two years on electrospray combustion in reduced-gravity environments such as those achievable at NASA microgravity test facilities.

Gomez, Alessandro↗

Multigrid Methods for Fully Implicit Oil Reservoir Simulation

In this paper we consider the simultaneous flow of oil and water in reservoir rock. This displacement process is modeled by two basic equations: the material balance or continuity equations and the equation of motion (Darcy's law). For the numerical solution of this system of nonlinear partial differential equations there are two approaches: the fully implicit or simultaneous solution method and the sequential solution method. In the sequential solution method the system of partial differential equations is manipulated to give an elliptic pressure equation and a hyperbolic (or parabolic) saturation equation. In the IMPES approach the pressure equation is first solved, using values for the saturation from the previous time level. Next the saturations are updated by some explicit time stepping method; this implies that the method is only conditionally stable. For the numerical solution of the linear, elliptic pressure equation multigrid methods have become an accepted technique. On the other hand, the fully implicit method is unconditionally stable, but it has the disadvantage that in every time step a large system of nonlinear algebraic equations has to be solved. The most time-consuming part of any fully implicit reservoir simulator is the solution of this large system of equations. Usually this is done by Newton's method. The resulting systems of linear equations are then either solved by a direct method or by some conjugate gradient type method. In this paper we consider the possibility of applying multigrid methods for the iterative solution of the systems of nonlinear equations. There are two ways of using multigrid for this job: either we use a nonlinear multigrid method or we use a linear multigrid method to deal with the linear systems that arise in Newton's method. So far only a few authors have reported on the use of multigrid methods for fully implicit simulations. Two-level FAS algorithm is presented for the black-oil equations, and linear multigrid for two-phase flow problems with strong heterogeneities and anisotropies is studied. Here we consider both possibilities. Moreover we present a novel way for constructing the coarse grid correction operator in linear multigrid algorithms. This approach has the advantage in that it preserves the sparsity pattern of the fine grid matrix and it can be extended to systems of equations in a straightforward manner. We compare the linear and nonlinear multigrid algorithms by means of a numerical experiment.

Molenaar, J.↗