Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphase”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Electrolytes in Multiple-Phase Hydrogen Storage Reactions

Multiple-phase hydrogen storage materials such as metal alanates and borohydrides, and destabilized systems offer the possibility of high hydrogen storage capacity with favorable thermodynamics. However, the multiphase nature of these materials intrinsically limits the kinetics due to the required transport of species between phases, which are typically in dry powder form. To address this limitation, the influence of added electrolytes is explored. This approach is motivated by analogy with similar multiphase battery reactions that show reduced kinetic limitations while necessarily containing electrolytes. Previous experimental results showing improved kinetics for MgH2/Sn (using a LiBH4/KBH4 eutectic electrolyte) and NaAlH4 (using a diglyme electrolyte) are further analyzed in terms of this analogy. The results show that the analogy is useful and rate constants are increased. Importantly, the inclusion of an electrolyte also appears to alleviate the continuously decreasing rates with the extent of reaction, which is characteristic of many multiphase hydrides. Instead, reaction rates are approximately constant until near completion. Together, these effects can lead to >10× shorter overall reaction times. In addition, new results are presented for the hydrogenation of MgB2 using Li/K/CsI and Li/K/CsCl eutectic electrolytes, where >60% conversion to Mg(BH4)2 was demonstrated at 350 bar.

08 HYDROGEN↗
DOE OSTI2017

The Presence of Thermally Unstable X-Ray Filaments and the Production of Cold Gas in the NGC 5044 Group

We present the results of a deep Chandra observation of the X-ray bright moderate-cooling flow group NGC 5044 along with the observed correlations between the ionized, atomic, and molecular gas in this system. The Chandra observation shows that the central AGN has undergone two outbursts in the past 10{sup 8} years, based on the presence of two pairs of nearly bipolar X-ray cavities. The molecular gas and dust within the central 2 kpc is aligned with the orientation of the inner pair of bipolar X-ray cavities, suggesting that the most recent AGN outburst had a dynamical impact on the molecular gas. NGC 5044 also hosts many X-ray filaments within the central 8 kpc, but there are no obvious connections between the X-ray and H α filaments and the more extended X-ray cavities that were inflated during the prior AGN outburst. Using the line width of the blended Fe-L line complex as a diagnostic for multiphase gas, we find that the majority of the multiphase thermally unstable gas in NGC 5044 is confined within the X-ray filaments. While the cooling time and entropy of the gas within the X-ray filaments are very similar, not all filaments show evidence of gas cooling or an association with H α emission. We suggest that the various observed properties of the X-ray filaments are suggestive of an evolutionary sequence where thermally unstable gas begins to cool, becomes multiphased, develops H α emitting plasma, and finally produces cold gas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Production of Cool Gas in Thermally Driven Outflows

Galactic outflows commonly contain multiphase gas, and its physical origin requires explanation. Using the Cholla Galactic OutfLow Simulations suite of high-resolution isolated galaxy models, we demonstrate the viability of rapid radiative cooling as a source of fast-moving ( v ~ 1000 km s -1 ), cool (10 4 K) gas observed in absorption-line studies of outflows around some star-forming galaxies. By varying the mass loading and geometry of the simulated winds, we identify a region of parameter space that leads to cool gas in outflows. In particular, when using an analytically motivated central feedback model, we find that cooling flows can be produced with reasonable mass-loading rates (${\dot{M}}_{\mathrm{wind}}/{\dot{M}}_{\mathrm{SFR}}\sim 0.5$), provided that the star formation rate surface density is high. When a more realistic clustered feedback model is applied, destruction of high-density clouds near the disk and interactions between different outflow regions indicate that lower mass-loading rates of the hot gas within the feedback region may still produce multiphase outflows. These results suggest an origin for fast-moving cool gas in outflows that does not rely on directly accelerating cool gas from the interstellar medium. Finally, these cooling flows may additionally provide an explanation for the multiphase gas ubiquitously observed in the halos of star-forming galaxies at low redshift.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Reduced Order Description of Experimental Two-Phase Pipe Flows: Characterization of Flow Structures and Dynamics via Proper Orthogonal Decomposition

Multiphase pipe flow is investigated using proper orthogonal decomposition for tomographic X-ray data, where holdup, cross-sectional phase distributions and phase interface characteristics within the pipe are obtained. Six cases of stratified and mixed flow with water content of 10%, 30% and 80% are investigated to gain insight into effects of velocity and proportion of water on the flow fields. Dispersed and slug flows are separately analyzed to consider the added interface complexity of the flow fields. These regimes are also highly applicable to industry operational flows. Instantaneous and fluctuating phase fractions of the four flow regime are analyzed and reduced order dynamical descriptions are generated. Stratified flow cases display coherent structures that highlight the liquid-liquid interface location while the mixed flow cases show minimal coherence of the eigenmodes. The dispersed flow displays coherent structures for the first few modes near the horizontal center of the pipe, representing the liquid-liquid interface location while the slug flow case shows coherent structures that correspond to the cyclical formation and break up of the slug in the first 5 modes. The low order descriptions of the high water content, stratified flow field indicates that main characteristics can be captured with minimal degrees of freedom. Reconstructions of the dispersed flow and slug flow cases indicate that dominant features are observed in the low order dynamical description utilizing less than 1%of the full order model. POD temporal coefficients a1, a2 and a3 show a high level of interdependence for the slug flowcase. The coefficients also describe the phase fraction holdup as a function of time for both dispersed and slug flow. The second coefficient, a2, and the centerline holdup profile show a mean percent difference below 9% between the two curves. The mathematical description obtained from the decomposition will deepen the understanding of multiphase flow characteristics and is applicable to long distance multiphase transport pipelines, fluidized beds, hydroelectric power and nuclear processes to name a few.

Viggiano, Bianca↗
DOE OSTI2025

Experimental Investigation of Vapor Formation in Liquid CO2 Flow Through a Converging-Diverging Nozzle

Carbon dioxide is an attractive working fluid for many cycles, including for pumped thermal energy storage (PTES). A challenge with some proposed sCO2 PTES cycles is the operation of sCO2 machinery outside the typical bounds of experience, with local phase change from the liquid state being particularly unknown. Presently, there is insufficient data in the literature regarding multiphase CO2 to adequately design a multiphase-tolerant turbine, so generation of foundational data is required. This experimental study investigates the flow characteristics of sub-sonic liquid CO2 undergoing expansion and phase change in a converging-diverging nozzle. The nozzle is instrumented to measure static pressure, unsteady pressure, temperature, and density. The static pressure transducers are located at 27 axial locations to accurately characterize the pressure profile in the nozzle. High-accuracy RTDs are located at the entrance and exit of the nozzle, and three dynamic pressure transducers are strategically located to capture any unsteady phenomena. During testing, values of mass flow and nozzle inlet pressure are swept to vary the pressure drop and fluid properties. The measured total pressure drop in the nozzle is compared to a homogenous model and the Lockhart-Martinelli correlation method, with the latter predicting loss quite closely. The resulting data set is valuable for validating multiphase numerical models in a simple geometry before implementation of these models in turbomachinery design.

25 ENERGY STORAGE↗
DOE OSTI2025

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗
DOE OSTI2012

Open-source MFIX-DEM software for gas-solids flows: Part II — Validation studies

With rapid advancements in computer hardware and numerical algorithms, computational fluid dynamics (CFD) has been increasingly employed as a useful tool for investigating the complex hydrodynamics inherent in multiphase flows. An important step during the development of a CFD model and prior to its application is conducting careful and comprehensive verification and validation studies. Accordingly, efforts to verify and validate the open-source MFIX-DEM software, which can be used for simulating the gas–solids flow using an Eulerian reference frame for the continuum fluid and a Lagrangian discrete framework (Discrete Element Method) for the particles, have been made at the National Energy Technology Laboratory (NETL). In part I of this paper, extensive verification studies were presented and in this part, detailed validation studies of MFIX-DEM are presented. A series of test cases covering a range of gas–solids flow applications were conducted. In particular the numerical results for the random packing of a binary particle mixture, the repose angle of a sandpile formed during a side charge process, velocity, granular temperature, and voidage profiles from a bounded granular shear flow, lateral voidage and velocity profiles from a monodisperse bubbling fluidized bed, lateral velocity profiles from a spouted bed, and the dynamics of segregation of a binary mixture in a bubbling bed were compared with available experimental data, and in some instances with empirical correlations. In addition, sensitivity studies were conducted for various parameters to quantify the error in the numerical simulation.

42 ENGINEERING↗
DOE OSTI2017

Predicting the performance uncertainty of a 1-MW pilot-scale carbon capture system after hierarchical laboratory-scale calibration and validation

A challenging problem in designing pilot-scale carbon capture systems is to predict, with uncertainty, the adsorber performance and capture efficiency under various operating conditions where no direct experimental data exist. Motivated by this challenge, we previously proposed a hierarchical framework in which relevant parameters of physical models were sequentially calibrated from different laboratory-scale carbon capture unit (C2U) experiments. Specifically, three models of increasing complexity were identified based on the fundamental physical and chemical processes of the sorbent-based carbon capture technology. Results from the corresponding laboratory experiments were used to statistically calibrate the physical model parameters while quantifying some of their inherent uncertainty. The parameter distributions obtained from laboratory-scale C2U calibration runs are used in this study to facilitate prediction at a larger scale where no corresponding experimental results are available. In this paper, we first describe the multiphase reactive flow model for a sorbent-based 1-MW carbon capture system then analyze results from an ensemble of simulations with the upscaled model. The simulation results are used to quantify uncertainty regarding the design’s predicted efficiency in carbon capture. In particular, we determine the minimum gas flow rate necessary to achieve 90% capture efficiency with 95% confidence.

bubbling bed↗
DOE OSTI2018

Modelling of Blast Furnace with Respective Chemical Reactions in Coke and Ore Burden Layers

The ironmaking blast furnace (BF) is an efficient chemical reactor for producing liquid iron from solid iron ore, where the solids of coke and iron ore are charged in alternative layers and different chemical reactions occur in the two solid layers as they descend. Such respective reacting burden layers have not been considered explicitly in the previous BF models. In this article, a mathematical model based on multi-fluid theory is developed for describing the multiphase reacting flows considering the respective reacting burden layers. Then, this model is applied to a BF, covering the area from the burden surface at the furnace top to the liquid surface above the hearth, to describe the inner states of a BF in terms of the multiphase flows, temperature distribution and reduction process. The results show that some key important features in the layered burden with respective chemical reactions are captured, including fluctuating iso-lines in terms of gas flow and thermochemical behaviours; particularly the latter cannot be well captured in the previous BF models. The temperature difference between gas–solid phases is found to be larger near the raceway, at the cohesive zone and at the furnace top, and the thermal reserved zone can be identified near the shaft. Three chemical reserve zones of hematite, magnetite and wustite can also be observed near the stockline, in the shaft near the wall and near centre, respectively. Inside each reserve zone, the corresponding ferrous oxides stay constantly high in alternative layers; the overall performance indicators including gas utilization efficiency and reduction degree also stay stable in an alternative-layered structure. This model provides a cost-effective tool to investigate the BF in-furnace process and optimize BF operation.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Automated contact angle estimation for three-dimensional X-ray microtomography data

Multiphase flow in capillary regimes is a fundamental process in a number of geoscience applications. The ability to accurately define wetting characteristics of porous media can have a large impact on numerical models. In this paper, a newly developed automated three-dimensional contact angle algorithm is described and applied to high-resolution X-ray microtomography data from multiphase bead pack experiments with varying wettability characteristics. The algorithm calculates the contact angle by finding the angle between planes fit to each solid/fluid and fluid/fluid interface in the region surrounding each solid/fluid/fluid contact point. Results show that the algorithm is able to reliably compute contact angles using the experimental data. The in situ contact angles are typically larger than flat surface laboratory measurements using the same material. Furthermore, wetting characteristics in mixed-wet systems also change significantly after displacement cycles.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2019

An extension of Darcy’s law incorporating dynamic length scales

We propose a physics-based, macroscale formulation of multiphase porous-media flows that both honors the validity of Darcy’s law in steady or near-steady flows and accommodates the effects of heterogeneities and nonlinearities in unsteady flows. The new formulation recognizes that parameters characterizing the system operate at different length-scales. In particular, the use of Darcy’s law, predicated on the assumption of near-steady flows, requires dynamic length scales, owing to the possibility of rapid fluctuations in rock properties and fluid saturations attributable to heterogeneity and nonlinearity. We accommodate these dynamic length-scales through dynamic spatial averaging. The length and position of the averaging window are characterized by the length of the mixing zone and the direction of propagation of information in the transport process. Here, we validate the proposed formulation by comparing highly accurate,two dimensional numerical solutions against core-scale displacement experiments. The proposed paradigm is consistent with the classical multiphase Darcy formulation, in the sense that the latter also represents the results of an averaging approach.

42 ENGINEERING↗
DOE OSTI2020

Comparison of Surface Tension Generation Methods in Smoothed Particle Hydrodynamics for Dynamic Systems

Developing robust numerical models of dynamic surface tension dominated multiphase systems is an ongoing challenge, especially in scenarios with large density and viscosity ratios. This is critical to the design and understanding of various physical and engineering systems, such as fluidized beds, fuel injectors, and drug delivery schemes. Much of the computational work in surface tension dominated multiphase flows has employed the continuum surface force method (CSF) of Brackbill et al. [1], which recasts surface tension from a surface force to a volumetric force that can be imposed in the vicinity of an interface. The CSF method produces accurate results across a variety of systems, however it relies on the identication of surface normals, which can be unreliable under certain conditions. Alternative methods of simulating surface tension have been proposed. Here the advantages and disadvantages of the CSF method in comparison to a pairwise forces (PF) method proposed by Tartakovsky and Panchenko [2] are explored. The CSF and PF methods are used in a smoothed particle hydrodynamics (SPH) framework to model dynamic systems. Results are compared to existing test cases from the literature and to analytic solutions derived from fundamental normal mode behavior of bubbles and droplets. A physical system for which the PF method is more stable and physically appropriate than the CSF method is identied.

smoothed particle hydrodynamics, multiphase flow↗
DOE OSTI2025

Enhanced MPM framework with multipatch isogeometric analysis for geotechnical applications

Achieving stable stress solutions at large strains using the Material Point Method (MPM) is challenging due to the accumulation of errors associated with geometry discretization, cell-crossing noise, and volumetric locking. Several simplified attempts exist in the literature to mitigate these errors, including higher-order frameworks. However, the stability of the MPM solution in such frameworks has been limited to simple geometries and the single-phase formulation (i.e., neglecting pore fluid). Although never explored, multipatch isogeometric analysis offers desirable qualities to simulate complex geometries while mitigating errors in the MPM. The degree of required high-order spatial integration has also never been investigated to infer a minimum limit for the stability of the stress solution in MPM. This paper presents a general-purpose numerical framework for simulating stable stresses in porous media, capturing both near incompressibility and multiphase interactions. First, the numerical framework is presented considering Non-Uniform Rational B-splines (NURBS) to perform isogeometric analysis (IGA) in MPM. Additionally, a volumetric strain smoothing algorithm is used to alleviate errors associated with volumetric locking. Second, the manifestation of cell-crossing errors is assessed via a series of problems with orders ranging from linear to cubic interpolation functions. Third, the use of NURBS is investigated and verified for problems with circular geometries. Finally, multipatch analysis is deployed to simulate plane strain and 3D penetration in soils, considering nearly incompressible elastoplastic (total stress) analysis and fully-coupled hydro-mechanical (effective stress) analysis. The stability of the solution is also analyzed for different constitutive models. From the results, it can be concluded that the framework using cubic interpolation functions with strain smoothing is the most convenient, presenting stable stress solutions for a broad range of multiphase geotechnical applications.

58 GEOSCIENCES↗
DOE OSTI2022

Transient simulation of biomass combustion in a circulating fluidized bed riser

Interest in circulating fluidized bed (CFB) boilers as a power generation technology has sky-rocketed in recent years because of several advantages this technology offers over conventional boilers, such as increased gas-solid mixing, which results in higher combustion efficiency and the ability to use lower rank fuels. CFB combustors are operated at lower temperatures than conventional thermal power generation combustors, thus reducing NO x emissions, while SO 2 emissions can be conveniently controlled through the addition of Ca-based sulfur sorbents within the combustor. Herin this paper summarizes the modeling effort on a 50 kW th CFB combustor designed, built, and operated at CanmetENERGY in Ottawa, Canada. The numerical model employs the multiphase particle-in-cell (PIC) approach in the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite. The MFiX-PIC model parameters for the simulation are tuned against cold-flow experiments from CanmetENERGY using olivine sand as the inert bed material. It is shown that for the relatively coarse fluid meshes and large parcel sizes necessitated by the scale of the simulation, filter size dependent corrections to the drag law must be incorporated to ensure accuracy of the simulation results. The validated cold flow model is extended to simulate reacting flow with torrefied hardwood as the feedstock and to validate the combustion reaction scheme. The species concentrations at the riser outlet are compared against CanmetENERGY’s experiments and show satisfactory agreement. The simulations demonstrate the ability of MFiX-PIC to accurately capture both the physics and chemistry of a CFB combustor at bench scales, which can be further extended to pilot- and industrial-scale systems.

09 BIOMASS FUELS↗
DOE OSTI2018

Burial-driven methane recycling in marine gas hydrate systems

Natural gas hydrate may be buried with sediments until it is no longer stable at a given pressure and temperature, resulting in conversion of hydrate into free gas. This gas may migrate upward and recycle back into the hydrate stability zone to form hydrate. As of yet, however, no quantitative description of the methane recycling process has been developed using multiphase flow simulations to model burial-driven gas hydrate recycling. In this study, we present a series of 1D multiphase transport simulations to investigate the methane recycling process in detail. By invoking the effects of capillary phenomena on hydrate and gas formation in pores of varying size, we find that a free gas phase can migrate a significant distance above the bulk base of hydrate stability. Since the top of the free gas occurrence is often identified as the base of the hydrate stability zone from seismic data, our results demonstrate that not only could this assumption mischaracterize a hydrate system, but that under recycling conditions the highest hydrate saturations can occur beneath the top of the free gas occurrence. We show that the presence of pore size distributions requires a replacement zone through which hydrate saturations progressively decrease with depth and are replaced with free gas. This replacement zone works to buffer against significant gas buildup that could lead to fracturing of overlying sediments. In conclusion, this work provides a framework for simulating flow and transport of methane within the 3-phase stability zone from a mass conservation perspective.

58 GEOSCIENCES↗
DOE OSTI2020

Development of a CPU/GPU portable software library for Lagrangian–Eulerian simulations of liquid sprays

The Lagrangian–Eulerian method is widely used for simulations of fuel sprays in turbulent combustion because of the advantage of treating the spray droplets as discrete points. One challenge of the Lagrangian–Eulerian method is the intense computational requirement when tracking the large number of Lagrangian particles needed for high fidelity. We have developed a performance-portable library, Grit, to track the Lagrangian particles in parallel on central processing unit (CPU) and graphics processing unit (GPU) accelerated high performance computing (HPC) architectures. Grit is a C++ library which employs Message Passing Interface (MPI) for distributed memory parallelism and Kokkos programming model for on-node shared memory parallelism with performance portability across different architectures of GPUs and multi-core/manycore CPUs. The parallel algorithms, key parallel kernels, and their performances on the pre-exascale supercomputer, Summit, are presented. Grit is coupled with a direct numerical simulation (DNS) solver, S3D, for multiphase simulations. A conservative formulation has been developed and implemented in Grit for phase coupling with thermodynamic consistency. The formulation separates the conservation of mass, momentum and energy from the physical models to prevent accidental violation of conservation laws due to inconsistent models. The formulation also enforces consistent definitions of enthalpies of the fuel for both phases and the latent heat of evaporation. Finally, simulations of turbulent particle-laden flow and the evaporation of dilute turbulent spray jet are performed to verify the software implementation, and to demonstrate the scalability of Grit for large-scale multiphase simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

A mass-conserving lattice Boltzmann method with dynamic grid refinement for immiscible two-phase flows

A mass-conserving lattice Boltzmann method (LBM) for multiphase flows is presented in this paper. The proposed LBM improves a previous model (Lee and Liu, 2010 [21]) in terms of mass conservation, speed-up, and efficiency, and also extends its capabilities for implementation on non-uniform grids. The presented model consists of a phase-field lattice Boltzmann equation (LBE) for tracking the interface between different fluids and a pressure-evolution LBM for recovering the hydrodynamic properties. In addition to the mass conservation property and the simplicity of the algorithm, the advantages of the current phase-field LBE are that it is an order of magnitude faster than the previous interface tracking LBE proposed by Lee and Liu (2010) [21] and it requires less memory resources for data storage. Meanwhile, the pressure-evolution LBM is equipped with a multi-relaxation-time (MRT) collision operator to facilitate attainability of small relaxation rates thereby allowing simulation of multiphase flows at higher Reynolds numbers. Additionally, we reformulate the presented MRT-LBM on nonuniform grids within an adaptive mesh refinement (AMR) framework. Various benchmark studies such as a rising bubble and a falling drop under buoyancy, droplet splashing on a wet surface, and droplet coalescence onto a fluid interface are conducted to examine the accuracy and versatility of the proposed AMR-LBM. The proposed model is further validated by comparing the results with other LB models on uniform grids. A factor of about 20 in savings of computational resources is achieved by using the proposed AMR-LBM. As a more demanding application, the Kelvin–Helmholtz instability (KHI) of a shear-layer flow is investigated for both density-matched and density-stratified binary fluids. The KHI results of the density-matched fluids are shown to be in good agreement with the benchmark AMR results based on the sharp-interface approach. When a density contrast between the two fluids exists, a typical chaotic structure in the flow field is observed at a Reynolds number of 10000, which indicates that the proposed model is a promising tool for direct numerical simulation of two-phase flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗