skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Quantification of mixing in vesicle suspensions using numerical simulations in two dimensions

ORCiD logo [1];  [2];  [3]
  1. Department of Mechanical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA
  2. Department of Scientific Computing, Florida State University, Tallahassee, Florida 32306, USA
  3. Department of Mechanical Engineering, The University of Texas at Austin, Austin, Texas 78712, USA, Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas 78712, USA
Publication Date:
Sponsoring Org.:
OSTI Identifier:
Grant/Contract Number:
SC0010518 and DE-SC0009286
Resource Type:
Journal Article: Publisher's Accepted Manuscript
Journal Name:
Physics of Fluids
Additional Journal Information:
Journal Volume: 29; Journal Issue: 2; Related Information: CHORUS Timestamp: 2018-02-14 12:33:20; Journal ID: ISSN 1070-6631
American Institute of Physics
Country of Publication:
United States

Citation Formats

Kabacaoğlu, G., Quaife, B., and Biros, G.. Quantification of mixing in vesicle suspensions using numerical simulations in two dimensions. United States: N. p., 2017. Web. doi:10.1063/1.4975154.
Kabacaoğlu, G., Quaife, B., & Biros, G.. Quantification of mixing in vesicle suspensions using numerical simulations in two dimensions. United States. doi:10.1063/1.4975154.
Kabacaoğlu, G., Quaife, B., and Biros, G.. Wed . "Quantification of mixing in vesicle suspensions using numerical simulations in two dimensions". United States. doi:10.1063/1.4975154.
title = {Quantification of mixing in vesicle suspensions using numerical simulations in two dimensions},
author = {Kabacaoğlu, G. and Quaife, B. and Biros, G.},
abstractNote = {},
doi = {10.1063/1.4975154},
journal = {Physics of Fluids},
number = 2,
volume = 29,
place = {United States},
year = {Wed Feb 01 00:00:00 EST 2017},
month = {Wed Feb 01 00:00:00 EST 2017}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1063/1.4975154

Citation Metrics:
Cited by: 1work
Citation information provided by
Web of Science

Save / Share:
  • A stochastic approach based on generalized polynomial chaos (gPC) is used to quantify the error in large-eddy simulation (LES) of a spatially evolving mixing layer flow and its sensitivity to different simulation parameters, viz., the grid stretching in the streamwise and lateral directions and the subgrid-scale (SGS) Smagorinsky model constant (C{sub S}). The error is evaluated with respect to the results of a highly resolved LES and for different quantities of interest, namely, the mean streamwise velocity, the momentum thickness, and the shear stress. A typical feature of the considered spatially evolving flow is the progressive transition from a laminarmore » regime, highly dependent on the inlet conditions, to a fully developed turbulent one. Therefore, the computational domain is divided in two different zones (inlet dependent and fully turbulent) and the gPC error analysis is carried out for these two zones separately. An optimization of the parameters is also carried out for both these zones. For all the considered quantities, the results point out that the error is mainly governed by the value of the C{sub S} constant. At the end of the inlet-dependent zone, a strong coupling between the normal stretching ratio and the C{sub S} value is observed. The error sensitivity to the parameter values is significantly larger in the inlet-dependent upstream region; however, low-error values can be obtained in this region for all the considered physical quantities by an ad hoc tuning of the parameters. Conversely, in the turbulent regime the error is globally lower and less sensitive to the parameter variations, but it is more difficult to find a set of parameter values leading to optimal results for all the analyzed physical quantities. A similar analysis is also carried out for the dynamic Smagorinsky model, by varying the grid stretching ratios. Comparing the databases generated with the different subgrid-scale models, it is possible to observe that the error cost function computed for the streamwise velocity and for the momentum thickness is not significantly sensitive to the used SGS closure. Conversely, the prediction of the shear stress is much more accurate when using a dynamic subgrid-scale model and the variance of the error is lower in magnitude.« less
  • A recently developed model for simulating the dynamic behavior of silicate materials is applied to the loading and unloading properties of granite. Four time-resolved wave profile measurements on granite are presented and used to supplement Hugoniot data to constrain the model. This model is based on a dynamic process-dependent description of phase mixing. In the solid--solid mixed-phase region, the loading states are determined from a simple two-parameter constraint equation which relates the mass fraction of the high pressure phase {lambda} to the Gibbs free-energy difference between the phases {Delta}{ital G}. On release, the reverse phase transition is modeled as amore » two-stage transition with each stage described by a two-parameter equation with the same form as the loading equation. In all four of the simulations, the loading behavior of the shock-wave experiment is well represented. For three of the experiments the release behavior is accurately simulated up to the point at which the measured interface speed reaches its first minimum. For the remaining experiment the simulation matches only the earliest release behavior and then rapidly deviates from the measured data.« less
  • The Euler similarity criteria for laboratory experiments and time-dependent mixing transition are important concepts introduced recently for application to prediction and analysis of astrophysical phenomena. However, Euler scaling by itself provides no information on the distinctive spectral range of high Reynolds number turbulent flows found in astrophysics situations. On the other hand, time-dependent mixing transition gives no indication on whether a flow that just passed the mixing transition is sufficient to capture all of the significant dynamics of the complete astrophysical spectral range. In this paper, a new approach, based on additional insight gained from review of Navier-Stokes turbulence theory,more » is developed. It allows for revelations about the distinctive spectral scale dynamics associated with high Reynolds number astrophysical flows. From this perspective, the energy-containing range of the turbulent flow measured in a laboratory setting must not be unintentionally contaminated in such a way that the interactive influences of this spectral scale range in the corresponding astrophysical situation cannot be faithfully represented. In this paper, the concept of a minimum state is introduced as the lowest Reynolds number turbulent flow that a time-dependent mixing transition must achieve to fulfill this objective. Later in the paper, the Reynolds number of the minimum state is determined as 1.6x10{sup 5}. The temporal criterion for the minimum state is also obtained. The efforts here can be viewed as a unification and extension of the concepts of both similarity scaling and transient mixing transition concepts. Finally, the implications of our approach in planning future intensive laser experiments or massively parallel numerical simulations are discussed. A systematic procedure is outlined so that as the capabilities of the laser interaction experiments and supporting results from detailed numerical simulations performed in recently advanced supercomputing facilities increase progressively, a strategy is developed so that a progressively increasing range of dynamic structures and their statistical influences on evolving astrophysical flows can be attained in laboratory investigations.« less
  • An analytical model of a two-stroke cycle, reciprocating, compression ignition engine was used to investigate the ignition and combustion characteristics of coal/water slurry fuels. The engine cycle simulation was based on a thermodynamic analysis of the cylinder contents and consisted of models for the injection, ignition, combustion, mixing, heat transfer, work, and scavenging processes. The thermodynamic analysis resulted in a set of first-order nonlinear, ordinary differential equations which were numerically integrated to obtain instantaneous cylinder gas, droplet, and particle conditions. The simulation results were first compared to experimental data from a large, slow-speed (120 rpm) engine using a coal/water slurrymore » fuel.« less
  • We have achieved validation in the form of simulation-experiment agreement for Rayleigh-Taylor turbulent mixing rates (known as {alpha}) over the past decade. The problem was first posed sixty years ago. Recent improvements in our simulation technology allow sufficient precision to distinguish between mixing rates for different experiments. We explain the sensitivity and non-universality of the mixing rate. These playa role in the difficulties experienced by many others in efforts to compare experiment with simulation. We analyze the role of initial conditions, which were not recorded for the classical experiments of Youngs et al. Reconstructed initial conditions with error bars aremore » given. The time evolution of the long and short wave length portions of the instability are analyzed. We show that long wave length perturbations are strong at t = 0, but are quickly overcome by the rapidly growing short wave length perturbations. These conclusions, based solely on experimental data analysis, are in agreement with results from theoretical bubble merger models and numerical simulation studies but disagree with models based on superposition of modes.« less