Search NASA⌕ Search

SEARCH · Search NASA

Results for “PHASE VELOCITY”

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,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Characteristics of Non-Premixed Turbulent Flames in Microgravity

The overall objectives of this research are: (1) to obtain and analyze experimental data on flame images, and the spatial and temporal distributions of temperature, radiation, velocity and gas-phase species in microgravity turbulent gas-jet diffusion flames; and (2) to utilize these data to validate and refine the existing predictive capabilities. Work on this project commenced in June 1996. The first investigations on turbulent gas-jet diffusion flames in microgravity were initiated by Bahadori and co-workers in 1991. These studies have shown that significant differences exist in the transition processes in normal-gravity and microgravity flames, and that the turbulent flames in microgravity behave very differently as compared to their buoyancy-dominated normal-gravity counterparts. For example, in the transition regime while the visible flame height, for given fuel and nozzle size, in normal gravity decreases, the height of the microgravity flame increases. In the fully developed turbulent regime, the normal-gravity flame height is independent of injection velocity, whereas the microgravity flame height continues to increase, although at a lower rate than in the laminar and transitional regimes. Other differences between the normal-gravity and microgravity flames arise in the jet shear-layer instability characteristics, extent of the transitional regime and blow-off limit characteristics.

Hegde, Uday↗
NASA NTRS2011

Further Constraints and Uncertainties on the Deep Seismic Structure of the Moon

The Apollo Passive Seismic Experiment (APSE) consisted of four 3-component seismometers deployed between 1969 and 1972, that continuously recorded lunar ground motion until late 1977. The APSE data provide a unique opportunity for investigating the interior of a planet other than Earth, generating the most direct constraints on the elastic structure, and hence the thermal and compositional evolution of the Moon. Owing to the lack of far side moonquakes, past seismic models of the lunar interior were unable to constrain the lowermost 500 km of the interior. Recently, array methodologies aimed at detecting deep lunar seismic reflections found evidence for a lunar core, providing an elastic model of the deepest lunar interior consistent with geodetic parameters. Here we study the uncertainties in these models associated with the double array stacking of deep moonquakes for imaging deep reflectors in the Moon. We investigate the dependency of the array stacking results on a suite of parameters, including amplitude normalization assumptions, polarization filters, assumed velocity structure, and seismic phases that interfere with our desired target phases. These efforts are facilitated by the generation of synthetic seismograms at high frequencies (approx. 1Hz), allowing us to directly study the trade-offs between different parameters. We also investigate expected amplitudes of deep reflections relative to direct P and S arrivals, including predictions from arbitrarily oriented focal mechanisms in our synthetics. Results from separate versus combined station stacking help to establish the robustness of stacks. Synthetics for every path geometry of data were processed identically to that done with data. Different experiments were aimed at examining various processing assumptions, such as adding random noise to synthetics and mixing 3 components to some degree. The principal stacked energy peaks put forth in recent work persist, but their amplitude (which maps into reflector impedance contrast) and timing (which maps into reflector depth) depend on factors that are not well constrained -- most notably, the velocity structure of the overlying lunar interior. Thus, while evidence for the lunar core remains strong, the depths of imaged reflectors have associated uncertainties that will require new seismic data and observations to constrain. These results strongly advocate further investigations on the Moon to better resolve the interior (e.g., Selene missions), for the Moon apparently has a rich history of construction and evolution that is inextricably tied to that of Earth.

Lin, Pei-Ying Patty↗
DOE OSTI2014

Note: Phase retrieval method for analyzing single-phase displacement interferometry data

We present a phase retrieval method (PRM) for analyzing single-phase displacement interferometry measurements on rapidly changing velocity histories, including photon Doppler velocimetry (PDV). PRM identifies the peaks and valleys as well as zero-crossing points in a PDV time series, performs normalization and extracts point-by-point phase and thus velocity information. PRM does not require a wide time window as in sliding window Fourier transformation, and thus improves the effective temporal resolution. This method is implemented in analyzing PDV data obtained from gas gun experiments, and validated against simultaneous measurements with velocity interferometer system for any reflector.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

Characteristics of transverse waves in chromospheric mottles

Using data obtained by the high temporal and spatial resolution Rapid Oscillations in the Solar Atmosphere instrument on the Dunn Solar Telescope, we investigate at an unprecedented level of detail transverse oscillations in chromospheric fine structures near the solar disk center. The oscillations are interpreted in terms of propagating and standing magnetohydrodynamic kink waves. Wave characteristics including the maximum transverse velocity amplitude and the phase speed are measured as a function of distance along the structure's length. Solar magnetoseismology is applied to these measured parameters to obtain diagnostic information on key plasma parameters (e.g., magnetic field, density, temperature, flow speed) of these localized waveguides. The magnetic field strength of the mottle along the ∼2 Mm length is found to decrease by a factor of 12, while the local plasma density scale height is ∼280 ± 80 km.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Acceleration of electrons by a laser pulse at its output onto an optical surface of the vacuum – transparent medium interface. Laser synchrotron

We consider the electron dynamics in the field of an electromagnetic wave produced at the vacuum – transparent medium interface upon reflection from the boundary, close to total internal reflection. The propagation velocity of a constant phase of the electromagnetic wave along the interface can vary from c/n to infinity (c is the speed of light in vacuum, and n is the refractive index of the medium at the interface). In this case, there emerge regions of positive and negative phases of the field with wavelengths, approximately equal to half the wavelength of the original laser beam, which can propagate at a speed close to that of light in vacuum. If a beam of relativistic electrons propagates along the surface, they can gain energy and accelerate, as well as radiate. With closed trajectories of electron motion, a laser synchrotron will be implemented as a result of many acceleration cycles. (extreme laser radiation: physics and fundamental applications)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Experimental and first-principles studies on the elastic properties of α-hafnium metal under pressure

Compressional and shear wave velocities of the α phase of hafnium have been measured up to 10.4 GPa at room temperature using ultrasonic interferometry in a multi-anvil apparatus. A finite strain equation of state analysis yielded K s0 = 110.4 (5) GPa, G 0 = 54.7(5) GPa,K s0 ' = 3.7 and G 0 ' = 0.6 for the elastic bulk and shear moduli and their pressure derivatives at ambient conditions. Complementary to the experimental data, the single crystal elastic constants, elastic anisotropy and the unit cell axial ratio c/a of α-hafnium at high pressures were investigated by Density Functional Theory (DFT) based first principles calculations. A c/a value of 1.605 is predicted for α-Hf at 40 GPa, which is in excellent agreement with previous experimental results. The low-pressure derivative of the shear modulus observed in our experimental data up to 10 GPa was found to originate from the elastic constant C44 which exhibits negligible pressure dependence within the current experimental pressure range. At higher pressures (>10 GPa), C 44 was predicted to soften and the shear wave velocity ν S trended to decrease with pressure, which can be interpreted as a precursor to the α-ω transition similar to that observed in other group IV elements (titanium and zirconium). Here, the acoustic velocities, bulk and shear moduli, and the acoustic Debye temperature (θ D = 240.1 K) determined from the current experiments were all compared well with those predicted by our theoretical DFT calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Flow of immiscible ferrofluids in a planar gap in a rotating magnetic field

Analytical solutions are obtained for the steady, fully developed flow of two layers of immiscible ferrofluids of different thicknesses between two parallel plates. Interfacial linear and internal angular momentum balance relations are derived for the case when there is a ferrofluid-ferrofluid interface to obtain the translational and spin velocity profiles in the gap. As expected for the limit of low applied field amplitude, the magnitude of the translational velocity is directly proportional to the frequency of the applied magnetic field and to the square of the magnetic field amplitude. Expressions for the velocity profiles are obtained for the zero spin viscosity and non-zero spin viscosity cases and the effect of applied pressure gradient on the flows is studied. The spin velocity in both ferrofluid phases is in the direction of the rotating magnetic field, except for cases of extreme applied pressure gradients for which the fluid vorticity opposes the spin. We find that for the case of non-zero spin viscosity, flow reversals are predicted using representative ferrofluid property values and field conditions. The unique predictions of the solution with non-zero spin viscosity could be used to experimentally test the existence of couple stresses in ferrofluids and the validity of previously reported values of the so-called spin viscosity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1993

Current Status in Cavitation Modeling

Cavitation is a common problem for many engineering devices in which the main working fluid is in liquid state. In turbomachinery applications, cavitation generally occurs on the inlet side of pumps. The deleterious effects of cavitation include: lowered performance, load asymmetry, erosion and pitting of blade surfaces, vibration and noise, and reduction of the overall machine life. Cavitation models in use today range from rather crude approximations to sophisticated bubble dynamics models. Details about bubble inception, growth and collapse are relevant to the prediction of blade erosion, but are not necessary to predict the performance of pumps. An engineering model of cavitation is proposed to predict the extent of cavitation and performance. The vapor volume fraction is used as an indicator variable to quantify cavitation. A two-phase flow approach is employed with the assumption of the thermal equilibrium between liquid and vapor. At present velocity slip between the two phases is selected. Preliminary analyses of 2D flows shows qualitatively correct results.

Singhal, Ashok K.↗
DOE OSTI2015

In-flight observations of low-mode ρR asymmetries in NIF implosions

Charged-particle spectroscopy is used to assess implosion symmetry in ignition-scale indirect-drive implosions for the first time. Surrogate D{sup 3}He gas-filled implosions at the National Ignition Facility produce energetic protons via D+{sup 3}He fusion that are used to measure the implosion areal density (ρR) at the shock-bang time. By using protons produced several hundred ps before the main compression bang, the implosion is diagnosed in-flight at a convergence ratio of 3–5 just prior to peak velocity. This isolates acceleration-phase asymmetry growth. For many surrogate implosions, proton spectrometers placed at the north pole and equator reveal significant asymmetries with amplitudes routinely ≳10%, which are interpreted as ℓ=2 Legendre modes. With significant expected growth by stagnation, it is likely that these asymmetries would degrade the final implosion performance. X-ray self-emission images at stagnation show asymmetries that are positively correlated with the observed in-flight asymmetries and comparable in magnitude, contradicting growth models; this suggests that the hot-spot shape does not reflect the stagnated shell shape or that significant residual kinetic energy exists at stagnation. More prolate implosions are observed when the laser drive is sustained (“no-coast”), implying a significant time-dependent asymmetry in peak drive.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

The 1992/93 eclipse of 31 Cygni

Extensive new ultraviolet and optical spectra of an atmospheric eclipse define the physical properties throughout the wind and chromosphere of 31 Cyg. These data require mass loss of approximately 3 x 10(exp -8) solar masses/yr in a wind that may merge smoothly into the chromosphere. Considerations of how energy is injected into the wind, however, suggest that the chromosphere and wind are separate structures. Most, if not all, of the velocity structure in metallic lines, which we have heretofore simulated with Doppler widths in the range 15-25 km/s, results from differential expansion of the atmosphere. Electron densities in the inner R(sub star) of the chromosphere are in the range 1.5 x 10(exp 9) - 2 x 10(exp 8)/cu cm, which implies clumping of the gas. The ionization of oxygen and nitrogen is consistent with clumping by factors of 3-30 in the outer chromosphere, roughly the amount required to give enough gas pressure to support the chromosphere. Chromospheric gas in 31 Cyg becomes hotter with increasing height, thus with decreasing optical depth, in a way that seems similar for all the zeta Aur binaries. Excitation temperature for Fe II in 1992/93 rises from about 5000 K at the deepest points sampled to about 12,500 K high in the wind. Strengths of violet Balmer lines give an excitation temperature for hydrogen of 6200-6500 K above a radial mass column density of approximately 0.01 g/sq cm. This amount of excitation implies that Lyman-alpha is thermalized beneath about 4 x 10(exp -3) g/sq cm. The outer atmosphere was symmetrical to within a factor of 2 in 1992/93, although it was clearly variable at this level, and it had similar mass column densities as in 1982. One manifestation of the variability was a flow toward the B star at phases 0.013-0.022 spanning velocities 45-100 km/s. Several lines of evidence point to a complicated and variable ionization in the wind: At large distance from the K star, measured mass column densities are less by up to a factor of 3 than required by a smoothly flowing wind. Also, much of the gas beyond r = 350 solar radii = 1.75 R(sub star) has a negative radial velocity. Radial velocities of the shell lines imply the outer atmosphere rotates in the direction of orbital motion, possibly through the deflection of wind flow lines in this direction.

Eaton, Joel A.↗
DOE OSTI2016

The roll-up and merging of coherent structures in shallow mixing layers

The current study seeks a fundamental explanation to the development of two-dimensional coherent structures (2DCSs) in shallow mixing layers. A nonlinear numerical model based on the depth-averaged shallow water equations is used to investigate the temporal evolution of shallow mixing layers, where the mapping from temporal to spatial results is made using the velocity at the center of the mixing layers. The flow is periodic in the streamwise direction. Transmissive boundary conditions are used in the cross-stream boundaries to prevent reflections. Numerical results are compared to linear stability analysis, mean-field theory, and secondary stability analysis. Results suggest that the onset and development of 2DCS in shallow mixing layers are the result of a sequence of instabilities governed by linear theory, mean-field theory, and secondary stability theory. The linear instability of the shearing velocity gradient gives the onset of 2DCS. When the perturbations reach a certain amplitude, the flow field of the perturbations changes from a wavy shape to a vortical (2DCS) structure because of nonlinearity. The development of the vertical 2DCS does not appear to follow weakly nonlinear theory; instead, it follows mean-field theory. After the formation of 2DCS, separate 2DCSs merge to form larger 2DCS. In this way, 2DCSs grow and shallow mixing layers develop and grow in scale. The merging of 2DCS in shallow mixing layers is shown to be caused by the secondary instability of the 2DCS. Eventually 2DCSs are dissipated by bed friction. The sequence of instabilities can cause the upscaling of the turbulent kinetic energy in shallow mixing layers.

42 ENGINEERING↗
DOE OSTI2018

Interparameter trade-off quantification and reduction in isotropic-elastic full-waveform inversion: synthetic experiments and Hussar land data set application

The problem of inverting for multiple physical parameters in the subsurface using seismic full-waveform inversion (FWI) is complicated by interparameter trade-off arising from inherent ambiguities between different physical parameters. Parameter resolution is often characterized using scattering radiation patterns, but these neglect some important aspects of interparameter trade-off. More general analysis and mitigation of interparameter trade-off in isotropic-elastic FWI is possible through judiciously chosen multiparameter Hessian matrix-vector products. We show that products of multiparameter Hessian off-diagonal blocks with model perturbation vectors, referred to as interparameter contamination kernels , are central to the approach. We then apply the multiparameter Hessian to various vectors designed to provide information regarding the strengths and characteristics of interparameter contamination, both locally and within the whole volume. With numerical experiments, we observe that S-wave velocity perturbations introduce strong contaminations into density and phase-reversed contaminations into P-wave velocity, but themselves experience only limited contaminations from other parameters. Based on these findings, we introduce a novel strategy to mitigate the influence of interparameter trade-off with approximate contamination kernels. Furthermore, we recommend that the local spatial and interparameter trade-off of the inverted models be quantified using extended multiparameter point spread functions (EMPSFs) obtained with pre-conditioned conjugate-gradient algorithm. Compared to traditional point spread functions, the EMPSFs appear to provide more accurate measurements for resolution analysis, by de-blurring the estimations, scaling magnitudes and mitigating interparameter contamination. Approximate eigenvalue volumes constructed with stochastic probing approach are proposed to evaluate the resolution of the inverted models within the whole model. With a synthetic Marmousi model example and a land seismic field data set from Hussar, Alberta, Canada, we confirm that the new inversion strategy suppresses the interparameter contamination effectively and provides more reliable density estimations in isotropic-elastic FWI as compared to standard simultaneous inversion approach.

58 GEOSCIENCES↗
NASA NTRSDate not supplied

Coherence Analysis of the Space Launch System using Unsteady Pressure Sensitive Paint

Transonic buffet during atmospheric ascent is a major source of unsteady loading on launch vehicles that, in the past, has led to structural failures. Thus, determining accurate buffet forcing functions (BFFs) to properly predict the vehicle response to buffet is of vital importance in launch vehicle design. The state of the art for obtaining the BFFs relies on wind-tunnel tests where the fluctuating pressures are measured by pressure transducers (PTs) at sparse locations on a rigid, geometrically-scaled buffet model. To compute the BFFs, the outer mold line (OML) of the vehicle is mapped onto contiguous panels centered at the location of the PTs and the measured fluctuating pressures are integrated over the panels’ areas. To mitigate conservatism due to the assumption that the measured buffet pressures act in phase across each panel, coherence factors are applied, effectively reducing the integration areas and, therefore, the buffet forces. For coherence factors to provide the proper level of BFF attenuation, accurate knowledge of the spatial and temporal correlation of the buffet pressures is paramount. Unfortunately, even with hundreds of unsteady pressure transducers instrumenting the models, compromises must be made on the spatial resolution of the PTs. Typically, the PT layout aims at resolving the pressure correlation along the longitudinal axis of the vehicle, assuming that coherent structures propagate mostly in the longitudinal direction. As a result, the distribution of azimuthal/cross-stream correlation is not well known and its impact on the estimated BFFs is often neglected. To fill this and other knowledge gaps, extremely high spatial-resolution uPSP data were collected for three different production-design configurations of the Space Launch System in the NASA Ames Research Center Unitary Plan Wind Tunnel 11-Foot Transonic Wind Tunnel. Specifically, two cargo configurations, the Block 1 and Block 1B, and one crew configuration, the Block 1B crew, were surveyed at resolutions ranging from 600,000 to over 1 million uPSP measurement locations. To shed light on the temporal and spectral behavior of are presented. Several OML regions and flow features of interest are investigated, from theexpansion/shock on the Orion Multi-Purpose Crew Vehicle, to the terminal shock environment on the core stage, and the Strouhal shedding behind the boosters forward attach. The sensitivity of these environments to the vehicle attitude is examined. Furthermore, for selected panels,coherence factors that accurately capture the azimuthal coherence distribution of the buffetpressures are derived and their impact on the estimated BFFs is discussed. Finally, distributions of the local convection velocity and cross spectrum phase are presented.

buffet↗
DOE OSTI2014

A second-order accurate immersed boundary-lattice Boltzmann method for particle-laden flows

A new immersed boundary-lattice Boltzmann method (IB-LBM) is presented for fully resolved simulations of incompressible viscous flows laden with rigid particles. The immersed boundary method (IBM) recently developed by Breugem (2012) [19] is adopted in the present method, development including the retraction technique, the multi-direct forcing method and the direct account of the inertia of the fluid contained within the particles. The present IB-LBM is, however, formulated with further improvement with the implementation of the high-order Runge–Kutta schemes in the coupled fluid–particle interaction. The major challenge to implement high-order Runge–Kutta schemes in the LBM is that the flow information such as density and velocity cannot be directly obtained at a fractional time step from the LBM since the LBM only provides the flow information at an integer time step. This challenge can be, however, overcome as given in the present IB-LBM by extrapolating the flow field around particles from the known flow field at the previous integer time step. The newly calculated fluid–particle interactions from the previous fractional time steps of the current integer time step are also accounted for in the extrapolation. The IB-LBM with high-order Runge–Kutta schemes developed in this study is validated by several benchmark applications. It is demonstrated, for the first time, that the IB-LBM has the capacity to resolve the translational and rotational motion of particles with the second-order accuracy. The optimal retraction distances for spheres and tubes that help the method achieve the second-order accuracy are found to be around 0.30 and −0.47 times of the lattice spacing, respectively. Simulations of the Stokes flow through a simple cubic lattice of rotational spheres indicate that the lift force produced by the Magnus effect can be very significant in view of the magnitude of the drag force when the practical rotating speed of the spheres is encountered. This finding may lead to more comprehensive studies of the effect of the particle rotation on fluid–solid drag laws. It is also demonstrated that, when the third-order or the fourth-order Runge–Kutta scheme is used, the numerical stability of the present IB-LBM is better than that of all methods in the literature, including the previous IB-LBMs and also the methods with the combination of the IBM and the traditional incompressible Navier–Stokes solver. - Highlights: • The IBM is embedded in the LBM using Runge–Kutta time schemes. • The effectiveness of the present IB-LBM is validated by benchmark applications. • For the first time, the IB-LBM achieves the second-order accuracy. • The numerical stability of the present IB-LBM is better than previous methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Coarsening of complex microstructures following spinodal decomposition

Coarsening plays a pivotal role in materials engineering, but our understanding of the dynamics of coarsening in morphologically complex systems is still limited. In this paper, we examine the correlations between the interfacial velocity and interfacial morphologies, and then predict the evolution of mean curvature based on the correlations. Three simulated structures with varying volume fractions, two bicontinuous and one nonbicontinuous, are generated using the Cahn-Hilliard equation. We find general correlations between interfacial velocity and mean curvature, as well as between interfacial velocity and the surface Laplacian of the mean curvature. Furthermore, we find that the probability of finding a patch of interface with a given normal velocity and the same local principal curvatures is described well by a Gaussian distribution, independent of the principal curvature values and the volume fractions of the structures. We also find that average interfacial velocity is described by a polynomial of the mean curvature and the net curvature. Based on this finding, we develop a semi-analytical approach to predicting the rate of change of the mean curvature, which determines the morphological evolution of complex microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effects of oxide surface chemistry on diffusioosmosis

Diffusioosmosis is the movement of fluid induced by gradients in solute concentration. Recent studies suggest that in low permeability rocks, it may be the dominant mode of reactant transport and thus control rates of diagenesis, which cannot be adequately explained by pressure-driven flow alone. In this paper we investigate how the equilibrium between an oxide mineral and the surrounding fluid phase influences the diffusioosmotic velocity. We have developed a theory for the case of a thin double layer, where the Debye length is smaller than the characteristic pore size. Several factors contribute to the total fluid velocity: the chemical structure of the mineral surface, the electrolyte type and concentration gradient, and the solution pH. Individual factors can act in concert or in opposition, leading to widely varying magnitudes and directions of the velocity. The numerical results are within the range of the limited experimental data. Our results highlight how surface charging and surface complexation impact the flow, and how they depend on pH.

Diagenesis↗
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↗