Volume-averaged macroscopic equation for fluid flow in moving porous media

Darcy's law and the Brinkman equation are two main models used for creeping fluid flows inside moving permeable particles. For these two models, the time derivative and the nonlinear convective terms of fluid velocity are neglected in the momentum equation. In this paper, a new momentum equation including these two terms are rigorously derived from the pore-scale microscopic equations by the volume-averaging method, which can reduces to Darcy's law and the Brinkman equation under creeping flow conditions. Using the lattice Boltzmann equation method, the macroscopic equations are solved for the problem of a porous circular cylinder moving along the centerline of a channel. Galilean invariance of the equations are investigated both with the intrinsic phase averaged velocity and the phase averaged velocity. The results demonstrate that the commonly used phase averaged velocity cannot serve as the superficial velocity, while the intrinsic phase averaged velocity should be chosen for porous particulate systems

Seeking for sterile neutrinos with displaced leptons at the LHC

We study the signal of long-lived sterile neutrino at the LHC produced through the decay of the $W$ boson. It decays into charged lepton and jets. The characteristic signature is a hard prompt lepton and a lepton from the displaced decay of the sterile neutrino, which leads to a bundle of displaced tracks with large transverse impact parameter. Different from other studies, we neither reconstruct the displaced vertex nor place requirement on its invariant mass to maintain sensitivity for low sterile neutrino masses. Instead, we focus on the displaced track from the lepton. A difficulty for low mass sterile neutrino study is that the displaced lepton is usually \textit{non-isolated}. Therefore, leptons from heavy flavor quark is the major source of background. We closely follow a search for displaced electron plus muon search at CMS and study their control regions, which is related to our signal regions, in great detail to develop a robust estimation of the background for our signals. After further optimization on the signal limiting the number of jets, low $H_T$ and large lepton displacement $d_0$ to suppress SM background, we reach an exclusion sensitivity of about $10^{-8}$ ($10^{-5}$) for the mixing angle square at 10 (2) GeV sterile neutrino mass respectively. The strategy we propose can cover the light sterile masses complimentary to beam dump and forward detector experiments.Comment: 22 pages, 6 figures, 1 table; v2: matched to Journal version

Simulation of two-phase flows at large density ratios and high Reynolds numbers using a discrete unified gas kinetic scheme

In order to treat immiscible two-phase flows at large density ratios and high Reynolds numbers, a three-dimensional code based on the discrete unified gas kinetic scheme (DUGKS) is developed, incorporating two major improvements. First, the particle distribution functions at cell interfaces are reconstructed using a weighted essentially non-oscillatory scheme. Second, the conservative lower-order Allen-Cahn equation is chosen, instead of the higher-order Cahn-Hilliard equation, to evolve the free-energy based phase field governing the dynamics of two-phase interfaces. Five benchmark problems are simulated to demonstrate the capability of the approach in treating two phase flows at large density ratios and high Reynolds numbers, including three two dimensional problems (a stationary droplet, Rayleigh-Taylor instability, and a droplet splashing on a thin liquid film) and two three-dimensional problems (binary droplets collision and Rayleigh-Taylor instability). All results agree well with the previous numerical and the experimental results. In these simulations, the density ratio and Reynolds number can reach a large value of O(1000). Our improved approach sets the stage for the DUGKS scheme to handle realistic two-phase flow problems

Nondeterminstic ultrafast ground state cooling of a mechanical resonator

We present an ultrafast feasible scheme for ground state cooling of a mechanical resonator via repeated random time-interval measurements on an auxiliary flux qubit. We find that the ground state cooling can be achieved with \emph{several} such measurements. The cooling efficiency hardly depends on the time-intervals between any two consecutive measurements. The scheme is also robust against environmental noises.Comment: 4 pages, 3 figure

Immobilized Talaromyces thermophilus lipase as an efficient catalyst for the production of LML-type structured lipids

LML-type structured lipids are one type of medium- and long-chain triacylglycerols. LML was synthesized using immobilized Talaromyces thermophilus lipase (TTL)-catalyzed interesterification of tricaprylin and ethyl linoleate. The resin AB-8 was chosen, and the lipase/support ratio was determined to be 60 mg/g. Subsequently, the immobilized TTL with strict sn-1,3 regiospecificity was applied to synthesize LML. Under the optimized conditions (60 °C, reaction time 6 h, enzyme loading of 6% of the total weight of substrates, substrate of molar ratio of ethyl linoleate to tricaprylin of 6:1), Triacylglycerols with two long- and one medium-chain FAs (DL-TAG) content as high as 52.86 mol% was obtained. Scale-up reaction further verified the industrial potential of the established process. The final product contained 85.24 mol% DL-TAG of which 97 mol% was LML after purification. The final product obtained with the high LML content would have substantial potential to be used as functional oils

ALP explanation to the muon (g−2)(g-2) and its test at future Tera- ZZ and Higgs factories

Models with an axionlike particle (ALP) can provide an explanation for the discrepancy between experimental measurement of the muon anomalous-magnetic moment $(g-2)_μ$ and the Standard Model prediction. This explanation relies on the couplings of the ALP to the muon and the photon. We also include more general couplings to the electroweak gauge bosons and incorporate them in the calculations up to the 2-loop order. We investigate the existing experimental constraints and find that they do not rule out the ALP model under consideration as a possible explanation for the $(g-2)_μ$ anomaly. At the same time, we find the future Tera-Z and Higgs factories, such as the CEPC and FCC-ee, can completely cover the relevant parameter space through searches with final states $(γγ)γ$, $(μ^+μ^-)γ$, and $(μ^+μ^-)μ^+μ^-$

Kinematic and Dynamic Collision Statistics of Cloud Droplets From High-Resolution Simulations

We study the dynamic and kinematic collision statistics of cloud droplets for a range of flow Taylor microscale Reynolds numbers (up to 500), using a highly scalable hybrid direct numerical simulation approach. Accurate results of radial relative velocity (RRV) and radial distribution function (RDF) at contact have been obtained by taking advantage of their power-law scaling at short separation distances. Three specific but inter-related questions have been addressed in a systematic manner for geometric collisions of same-size droplets (of radius from 10 to 60 μm) in a typical cloud turbulence (dissipation rate at 400 cm2 s-3. Firstly, both deterministic and stochastic forcing schemes were employed to test the sensitivity of the simulation results on the large-scale driving mechanism. We found that, in general, the results are quantitatively similar, with the deterministic forcing giving a slightly larger RDF and collision kernel. This difference, however, is negligible for droplets of radius less than 30 μm. Secondly, we have shown that the dependence of pair statistics on the flow Reynolds number Rλ or larger scale fluid motion is of secondary importance, with a tendency for this effect to saturate at high enough Rλ leading to Rλ-independent results. Both DNS results and theoretical arguments show that the saturation happens at a smaller Rλ for smaller droplets. Finally, since most previous studies of turbulent collision of inertial particles concerned non-sedimenting particles, we have specifically addressed the role of gravity on collision statistics, by simultaneously simulating collision statistics with and without gravity. It is shown that the collision statistics is not affected by gravity when a \u3c ac, where the critical droplet radius ac is found to be around 30 μm for the RRV, and around 20 μm for the RDF. For larger droplets, gravity alters the particle-eddy interaction time and significantly reduces the RRV. The effect of gravity on the RDF is rather complex: gravity reduces the RDF for intermediate-sized droplets but enhances the RDF for larger droplets. In addition, we have studied the scaling exponents of both RDF and RRV, and found that gravity modifies the RDF scaling exponents for both intermediate and large particles, in a manner very similar to the effect of gravity on the RDF at contact. Gravity is shown to cause the scaling exponents for RDF and RRV to level off for large droplets, in contrast to diminishing exponents for non-sedimenting particles

Kinematic and Dynamic Pair Collision Statistics of Sedimenting Inertial Particles Relevant to Warm Rain Initiation

In recent years, direct numerical simulation (DNS) approach has become a reliable tool for studying turbulent collision-coalescence of cloud droplets relevant to warm rain development. It has been shown that small-scale turbulent motion can enhance the collision rate of droplets by either enhancing the relative velocity and collision efficiency or by inertia-induced droplet clustering. A hybrid DNS approach incorporating DNS of air turbulence, disturbance flows due to droplets, and droplet equation of motion has been developed to quantify these effects of air turbulence. Due to the computational complexity of the approach, a major challenge is to increase the range of scales or size of the computation domain so that all scales affecting droplet pair statistics are simulated. Here we discuss our on-going work in this direction by improving the parallel scalability of the code, and by studying the effect of large-scale forcing on pair statistics relevant to turbulent collision. New results at higher grid resolutions show a saturation of pair and collision statistics with increasing flow Reynolds number, for given Kolmogorov scales and small droplet sizes. Furthermore, we examine the orientation dependence of pair statistics which reflects an interesting coupling of gravity and droplet clustering

Effect of hydrodynamic heterogeneity on particle dispersion in a Taylor-Couette flow reactor with variable configurations of inner cylinder

BackgroundEffect of hydrodynamic heterogeneity on particle dispersion in a Taylor-Couette flow (TC) reactor with variable configurations of inner cylinder has been investigated using CFD modelling.MethodsParticle dispersion was tracked based on the Eulerian-Lagrangian approach, where the reactant solution phase was solved in the Eulerian reference frame, while the particle dispersion was calculated by tracking a large number of particles with consideration of the hydrodynamic forces acting on particles and adopting actual particle properties measured from the particle synthesis experiments.Significant FindingsThe simulation reveals that particle dispersion is significantly enhanced by increasing the inner cylinder rotational speed, characterized by particle distribution for both circular inner cylinder Taylor-Couette flow reactor (CTC) and lobed cross-section inner cylinder Taylor-Couette flow reactor (LTC). Particle trajectories or dispersion are influenced by the turbulent Taylor vortices. Particle radial dispersion affects the particle classification by presenting different particle axial velocities in radial direction, while particle axial dispersion can be seen as an indicator for global mixing occurring in the TC reactor, which is enhanced at high rotational speed, especially in the LTC. The calculated dispersion coefficient is found to be similar to the shape of particle size distribution found in the experiments