Source author record

Yantao Yang

Yantao Yang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2022arXiv

Flow induced dissolution of femtoliter surface droplet arrays

The dissolution of liquid nanodroplets is a crucial step in many applied processes, such as separation and dispersion in food industry, crystal formation of pharmaceutical products, concentrating and analysis in medical diagnosis, and drug delivery in aerosols. In this work, using both experiments and numerical simulations, we \textit{quantitatively} study the dissolution dynamics of femtoliter surface droplets in a highly ordered array under a uniform flow. Our results show that the dissolution of femoliter droplets strongly depends on their spatial positions relative to the flow direction, drop-to-drop spacing in the array, and the imposed flow rate. In some particular case, the droplet at the edge of the array can dissolve about 30% faster than the ones located near the centre. The dissolution rate of the droplet increases by 60% as the inter-droplet spacing is increased from 2.5 $μ$m to 20 $μ$m. Moreover, the droplets close to the front of flow commence to shrink earlier than those droplets in the center of the array. The average dissolution rate is faster for faster flow. As a result, the dissolution time $T_{i}$ decreases with the Reynolds number Re of the flow as $T_{i}\propto Re^{-3/4}$. The experimental results are in good agreement with numerical simulations where the advection-diffusion equation for the concentration field is solved and the concentration gradient on the surface of the drop is computed. The findings suggest potential approaches to manipulate nanodroplet sizes in droplet arrays simply by dissolution controlled by an external flow. The obtained droplets with varying curvatures may serve as templates for generating multifocal microlens in one array.

preprint2021arXiv

Layering and vertical transport in sheared double diffusive convection in the diffusive regime

A sequence of two and three-dimensional simulations is conducted for the double diffusive convection (DDC) flows in the diffusive regime subjected to an imposed shear. The flow is confined between two horizontal plates which are maintained at different constant temperature, salinity, and different velocity, thus setting up a shear across the flow. The lower plate is fixed at higher temperature and salinity, while the overall (unperturbed) density gradient is statically stable. For a wide range of control parameters, and for sufficiently strong perturbation of the conductive initial state, we find that staircase-like structures spontaneously develop, with relatively well-mixed layers separated by sharp interfaces of enhanced scalar gradient. Such staircases appear to be robust even in the presence of strong shear over very long times, although we typically observe early time coarsening of the number of observed layers. For the same set of control parameters, different asymptotic layered states, with markedly different vertical scalar fluxes, can arise for different initial perturbation structures. The imposed shear does significantly spatio-temporally modify the vertical transport of the various scalars. The flux ratio (i.e., the ratio between the density fluxes due to the total (convective and diffusive) salt flux and the total heat flux) is found, at steady state, to be essentially equal to the square root of the ratio of the salt diffusivity to the thermal diffusivity, consistently with the physical model originally proposed by Linden and Shirtcliffe (1978) and the variational arguments presented by Stern (1982) for unsheared double diffusive convection.

preprint2020arXiv

Subcritical behaviour in double diffusive convection within the diffusive regime

We conduct two- and three-dimensional simulations for double diffusive convection in the diffusive regime, where the fluid flow is driven by a destabilizing temperature gradient and stabilized by a stably stratified salinity gradient. We study how the heat flux, Reynolds number, and flow structures change with the density ratio $Λ$, which is the ratio of the buoyancy force induced by the salinity gradient to that by the temperature gradient. When $Λ$ increases from zero, the flow first behaves similarly as in pure Rayleigh-Bénard (RB) convection, both with respect to flow structure and to heat transport. The linear stability analysis of Baines & Gill (J. Fluid Mech., vol. 37, 1969, pp. 289-306) had estimated the critical density ratio $Λ_c$, above which the flow becomes stable. However, here we show that by using a large-scale circulation as initial condition (rather than the linear profiles assumed in the linear stability analysis), DDC in the diffusive regime can exhibit subcritical behaviour when $Λ> Λ_c$, i.e., coexistence of states at the same control parameters. Even though the density ratio becomes thousands times that of the critical value $Λ_c$, there is still convection with strongly enhanced heat transfer properties compared to the pure conduction case. We reveal the corresponding flow structures and find an unstably-stratified region sandwiched between two stably-stratified layers. Our results demonstrate the importance of the initial condition for DDC in the diffusive regime, especially in the situation of a large density ratio, which occurs in high-latitude ocean regions.

preprint2020arXiv

What rotation rate maximizes heat transport in rotating Rayleigh-Bénard convection with Prandtl number larger than one?

The heat transfer and flow structure in rotating Rayleigh-Bénard convection are strongly influenced by the Rayleigh ($Ra$), Prandtl ($Pr$), and Rossby ($Ro$) number. For $Pr\gtrsim 1$ and intermediate rotation rates, the heat transfer is increased compared to the non-rotating case. We find that the regime of increased heat transfer is subdivided into a low and a high $Ra$ number regime. For $Ra\lesssim 5\times10^8$ the heat transfer at a given $Ra$ and $Pr$ is highest at an optimal rotation rate, at which the thickness of the viscous and thermal boundary layer is about equal. From the scaling relations of the thermal and viscous boundary layer thicknesses, we derive that the optimal rotation rate scales as $1/Ro_\mathrm{opt} \approx 0.12 Pr^{1/2}Ra^{1/6}$. In the low $Ra$ regime the heat transfer is similar in a periodic domain and cylindrical cells with different aspect ratios, i.e.\ the ratio of diameter to height. This is consistent with the view that the vertically aligned vortices are the dominant flow structure. For $Ra\gtrsim 5\times10^8$ the above scaling for the optimal rotation rate does not hold anymore. It turns out that in the high $Ra$ regime, the flow structures at the optimal rotation rate are very different than for lower $Ra$. Surprisingly, the heat transfer in the high $Ra$ regime differs significantly for a periodic domain and cylindrical cells with different aspect ratios, which originates from the sidewall boundary layer dynamics and the corresponding secondary circulation.

preprint2016arXiv

Scaling laws and flow structures of double diffusive convection in the finger regime

Direct numerical simulations are conducted for double diffusive convection (DDC) bounded by two parallel plates, with fluid properties similar to the values of seawater. The DDC flow is driven by an unstable salinity difference and stabilized at the same time by a temperature difference. For these conditions the flow can be in the finger regime. We develop scaling laws for three key response parameters of the system: The non-dimensional salinity flux $Nu_S$ mainly depends on the salinity Rayleigh number $Ra_S$, which measures the strength of the salinity difference, and exhibits a very weak dependence on the density ratio $Λ$, which is the ratio of the buoyancy forces induced by two scalar differences. The non-dimensional flow velocity $Re$ and the non-dimensional heat flux $Nu_T$ are dependent on both $Ra_S$ and $Λ$. However, the rescaled Reynolds number $ReΛ^{α^{\rm eff}_u}$ and the rescaled convective heat flux $(Nu_T-1)Λ^{α^{\rm eff}_T}$ depend only on $Ra_S$. The two exponents are dependent on the fluid properties and are determined from the numerical results. Moreover, the behaviors of $Nu_S$ and $ReΛ^{α^{\rm eff}_u}$ agree with the predictions of the Grossmann-Lohse theory which was originally developed for the Rayleigh-Bénard flow. The non-dimensional salt-finger width and the thickness of the velocity boundary layers, after being rescaled by $Λ^{α^{\rm eff}_u/2}$, collapse and obey a similar power-law scaling relation with $Ra_S$. When $Ra_S$ is large enough, salt fingers do not extend from one plate to the other and horizontal zonal flows emerge in the bulk region. We then show that the current scaling strategy can be successfully applied to the experimental results of a heat-copper-ion system~(Hage and Tilgner, Phys. Fluids, 22, 076603, 2010).

preprint2016arXiv

Shock responses of nanoporous aluminum by molecular dynamics simulations

We present systematic investigations on the shock responses of nanoporous aluminum (np-Al) by nonequilibrium molecular dynamics simulations. The dislocation nucleation sites are found to concentrate in low latitude region near the equator of the spherical void surfaces. We propose a continuum wave reflection theory and a resolved shear stress model to explain the distribution of dislocation nucleation sites. The simulations reveals two mechanisms of void collapse: the plasticity mechanism and the internal jetting mechanism. The plasticity mechanism, which leads to transverse collapse of voids, prevails under relatively weaker shocks; while the internal jetting mechanism, which leads to longitudinal filling of the void vacuum, plays more significant role as the shock intensity increases. In addition, an abnormal thermodynamic phenomenon (i.e., arising of temperature with pressure dropping) in shocked np-Al is discovered. This phenomenon is incompatible with the conventional Rankine-Hugoniot theory, and is explained by the nonequilibrium processes involved in void collapse. The influences of void collapse on spall fracture of np-Al is studied. Under the same loading velocity, the spall strength of np-Al is found to be lower than that of single-crystal Al; but the spall resistance is higher in np-Al than in single-crystal Al. This is explained by the combined influences of thermal dissipation and stress attenuation during shock wave propagation in np-Al.

preprint2016arXiv

Vertically bounded double diffusive convection in the fingering regime: comparing no-slip vs free-slip boundary conditions

Vertically bounded fingering double diffusive convection (DDC) is numerically investigated, focusing on the influences of different velocity boundary conditions, i.e. the no-slip condition which is inevitable in the lab-scale experimental research, and the free-slip condition which is an approximation for the interfaces in many natural environments, such as the oceans. For both boundary conditions the flow is dominated by fingers and the global responses follow the same scaling laws, with enhanced prefactors for the free-slip cases. Therefore, the laboratory experiments with the no-slip boundaries serve as a good model for the finger layers in the ocean. Moreover, in the free-slip case although the tangential shear stress is eliminated at the boundaries, the local dissipation rate in the near-wall region may exceed the value found in the no-slip cases, which is caused by the stronger vertical motions of fingers and sheet structures near the free-slip boundaries. This counter intuitive result might be relevant for properly estimating and modelling the mixing and entrainment phenomena at free-surfaces and interfaces.

preprint2015arXiv

A multiple--resolution strategy for Direct Numerical Simulation of scalar turbulence

In this paper a numerical procedure to simulate low diffusivity scalar turbulence is presented. The method consists of using a grid for the advected scalar with a higher spatial resolutions than that of the momentum. The latter usually requires a less refined mesh and integrating both fields on a single grid tailored to the most demanding variable, produces an unnecessary computational overhead. A multiple resolution approach is used also in the time integration in order to maintain the stability of the scalars on the finer grid. The method is the more advantageous the less diffusive the scalar is with respect to momentum, therefore it is particularly well suited for large Prandtl or Schmidt number flows. However, even in the case of equal diffusivities the present procedure gives CPU time and memory occupation savings. The reason is that the absence of the pressure term in the scalar equation leads to much steeper gradients in the scalar field as compared to the velocity field.

preprint2015arXiv

From convection rolls to finger convection in double-diffusive turbulence

Double diffusive convection (DDC), which is the buoyancy driven flow with fluid density depending on two scalar components, is ubiquitous in many natural and engineering enviroments. Of great interests are scalers transfer rate and flow structures. Here we systematically investigate DDC flow between two horizontal plates, driven by an unstable salinity gradient and stabilized by a temperature gradient. Counterintuitively, when increasing the stabilizing temperature gradient, the salinity flux first increases, even though the velocity monotonically decreases, before it finally breaks down to the purely diffusive value. The enhanced salinity transport is traced back to a transition in the overall flow pattern, namely from large scale convection rolls to well-organised vertically-oriented salt fingers. We also show and explain that the unifying theory of thermal convection originally developed by Grossmann and Lohse for Rayleigh-Bénard convection can be directly applied to DDC flow for a wide range of control parameters (Lewis number and density ratio), including those which cover the common values relevant for ocean flows.

preprint2015arXiv

Inertial waves and mean velocity profiles in a rotating pipe and a circular annulus with axial flow

In this paper we solve the inviscid inertial wave solutions in a circular pipe or annulus rotating constantly about its axis with moderate angular speed. The solutions are constructed by the so-called helical wave functions. We reveal that the mean velocity profiles must satisfy certain conditions to accommodate the inertial waves at the bulk region away from boundary. These conditions require the axial and azimuthal components of the mean velocity take the shapes of the zeroth and first order Bessel functions of the first kind, respectively. The theory is then verified by data obtained from direct numerical simulations for both rotating pipe and circular annulus, and excellent agreement is found between theory and numerical results. Large scale vortex clusters are found in the bulk region where the mean velocity profiles match the theoretical predictions. The success of the theory in rotating pipe, circular annulus, and streamwise rotating channel suggests that such inertial waves are quite common in wall bounded flow with background rotation.

preprint2015arXiv

Salinity transfer in bounded double diffusive convection

The double diffusive convection between two parallel plates is numerically studied for a series of parameters. The flow is driven by the salinity difference and stabilized by the thermal field. Our simulations are directly compared to experiments by Hage and Tilgner (\emph{Phys. Fluids} 22, 076603 (2010)) for several sets of parameters and reasonable agreement is found. This in particular holds for the salinity flux and its dependence on the salinity Rayleigh number. Salt fingers are present in all simulations and extend through the entire height. The thermal Rayleigh number seems to have minor influence on salinity flux but affects the Reynolds number and the morphology of the flow. Next to the numerical calculation, we apply the Grossmann-Lohse theory for Rayleigh-Bénard flow to the current problem without introducing any new coefficients. The theory successfully predicts the salinity flux both with respect to the scaling and even with respect to the absolute value for the numerical and experimental results.