Filtration of aerosol particles by cylindrical fibers within a parallel and staggered array

The World Health Organization (WHO) in 2013 reported that more than seven million unexpected losses every year are credited to air contamination. Because of incredible adaptability and expense viability of fibrous filters, they are broadly used for removing particulates from gasses. The influence of appropriate parameters, e.g., the fiber arrangement, solid volume fraction (SVF or α), fluid flow face velocity (mean inlet velocity), and filter thickness (I _x ), on pressure drop and deposition efficiency are researched. Furthermore, to study the effects of variation of the laminar flow regime and fiber’s cross-sectional shape on the deposition of particles, only a single square fiber has been placed in a channel. By means of finite volume method (FVM), the 2-D motion of 100–1000 nm particles was investigated numerically. The Lagrangian method has been employed and the Saffman’s lift, Drag, and Brownian forces have been considered to affect this motion. Contribution of increasing the Reynolds number to filtration performance increased with smaller fine aerosols to a level of 59.72 %. However, for over 500 nm, the Re = 100 has more efficient results up to 26.97 %. Remarkably, the single square fiber in Re = 200 regime performs similarly to the optimum choice of multi-fibrous filters. It was portrayed the parallel circular multi-fibrous filter with a ratio of horizontal-to-vertical distances between fibers, l/h = 1.143; α = 0.687, I _x  = 116.572, and h/d _f  = 1.0 is the most efficient filter’s structure. The increase in the ratio of vertical distances between fibers-to-fiber’s diameter (h/d _f ) and decrease in SVF or α, results in a drastically decrement of the filtration performance of both parallel and staggered structures. The obtained results have been validated with previous research findings.     

Study on microchannel flows with a sudden contraction–expansion at a wide range of Knudsen number using lattice Boltzmann method

Numerical simulations have been performed on the pressure-driven rarefied flow through channels with a sudden contraction–expansion of 2:1:2 using isothermal two and three-dimensional lattice Boltzmann method (LBM). In the LBM, a Bosanquet-type effective viscosity and a modified second-order slip boundary condition are used to account for the rarefaction effect on gas viscosity to cover the slip and transition flow regimes, that is, a wider range of Knudsen number. Firstly, the in-house LBM code is verified by comparing the computed pressure distribution and flow pattern with experimental ones measured by others. The verified code is then used to study the effects of the outlet Knudsen number Kn _o , driving pressure ratio P _i /P _o , and Reynolds number Re, respectively, varied in the ranges of 0.001–1.0, 1.15–5.0, and 0.02–120, on the pressure distributions and flow patterns as well as to document the differences between continuum and rarefied flows. Results are discussed in terms of the distributions of local pressure, Knudsen number, centerline velocity, and Mach number. The variations of flow patterns and vortex length with Kn _o and Re are also documented. Moreover, a critical Knudsen number is identified to be Kn _oc  = 0.1 below and above which the behaviors of nonlinear pressure profile and velocity distribution and the variations of vortex length with Re upstream and downstream of constriction are different from those of continuum flows.     

Pressure drop of three-phase liquid–liquid–gas slug flow in round microchannels

In this paper we present a model for the calculation of pressure drop of three-phase liquid–liquid–gas slug flow in microcapillaries of a circular cross section. Introduced models consist of terms attributing for frictional and interfacial pressure drop, incorporating the presence of a stagnant thin film at the wall of the channel. Different formulations of the interfacial pressure drop equation were employed, using expressions developed by Bretherton (J Fluid Mech 10:166–188, 1961), Warnier et al. (Microfluid Nanofluid 8:33–45, 2010) or Ratulowski and Chang (Phys Fluids A 1:1642–1655, 1989). Models were validated experimentally using oleic acid–water–nitrogen and heptane–water–nitrogen three-phase flows in round Teflon or Radel R microchannels of 254- and 508-µm nominal inner diameter, for capillary numbers Ca _b between 10^?4 and 4.9 × 10^?1 and Reynolds numbers Re between 0.095 and 300. Best agreement between measured and calculated values of pressure drop, with relative error between ?22 and 19 % or ?20 and 16 %, is reached for Warnier’s or Ratulowski and Chang’s interfacial pressure drop equation, respectively. The results prove that three-phase slug flow pressure drop can be successfully predicted by extending existing two-phase slug flow correlations. Good agreement of Bretherton’s equation was reached only at lower Ca numbers, indicating that an extension of the interfacial pressure drop equation as performed by Warnier et al. (Microfluid Nanofluid 8:33–45, 2010) or Ratulowski and Chang (Phys Fluids A 1:1642–1655, 1989) for higher capillary numbers is necessary. Additionally it was demonstrated that pressure drop increases substantially if dry slug flow occurs or if microchannels with significant surface roughness are employed. Those influences were not accounted for in the models presented.     

Star reduction among minimal length addition chains

An addition chain for a natural number n is a sequence \({1=a_0 < a_1 < \cdots < a_r=n}\) of numbers such that for each 0 < i ≤ r, a _i  = a _j  + a _k for some 0 ≤ k ≤ j < i. The minimal length of an addition chain for n is denoted by ?(n). If j = i ? 1, then step i is called a star step. We show that there is a minimal length addition chain for n such that the last four steps are stars. Then we conjecture that there is a minimal length addition chain for n such that the last \({\lfloor\frac{\ell(n)}{2}\rfloor}\)-steps are stars. We verify that the conjecture is true for all numbers up to 2¹⁸. An application of the result and the conjecture to generate a minimal length addition chain reduce the average CPU time by 23–29% and 38–58% respectively, and memory storage by 16–18% and 26–45% respectively for m-bit numbers with 14 ≤ m ≤ 22.     

An orthogonal analysis method for decoupling the nozzle geometrical parameters of microthrusters

This paper proposes an orthogonal analysis method for decoupling the multiple nozzle geometrical parameters of microthrusters, thus an reconfigured design can be implemented to generate a proper thrust. In this method, the effects of various nozzle geometrical parameters, including throat width W _t , half convergence angle θ _in , half divergence angle θ _out , exit-to-throat section ratio W _e /W _t and throat radius of the curvature R _t /W _t , on the performance of microthrusters are sorted by range analysis. Analysis results show that throat width seriously affects thrust because range value of 67.53 mN is extremely larger than the range value of other geometry parameters. For average specific impulse (ASI), the range value of exit-to-throat section ratio W _e /W _t and half divergence angle θ _out are 4.82 s and 3.72 s, respectively. Half convergence angle with the range value of 0.39 s and throat radius with 0.32 s have less influence on ASI compared with exit-to-throat section ratio and half divergence angle. When increasing the half convergence angle from 10° to 40° and throat radius of the curvature from 3 to 9, average specific impulse initially decreases and then increases. A MEMS solid propellant thruster (MSPT) with the reconfigured geometrical parameters of nozzle is fabricated to verify the feasibility of the proposed method. The thrust of the microthruster can reach 25 mN. Power is estimated to be 0.84 W. This work provides design guideline to reasonably configure geometry parameters of microthruster.     

Negation-Limited Complexity of Parity and Inverters

In negation-limited complexity, one considers circuits with a limited number of NOT gates, being motivated by the gap in our understanding of monotone versus general circuit complexity, and hoping to better understand the power of NOT gates. We give improved lower bounds for the size (the number of AND/OR/NOT) of negation-limited circuits computing Parity and for the size of negation-limited inverters. An inverter is a circuit with inputs x ₁,…,x _n and outputs ¬ x ₁,…,¬ x _n . We show that: (a) for n=2 ^r ?1, circuits computing Parity with r?1 NOT gates have size at least 6n?log?₂(n+1)?O(1), and (b) for n=2 ^r ?1, inverters with r NOT gates have size at least 8n?log?₂(n+1)?O(1). We derive our bounds above by considering the minimum size of a circuit with at most r NOT gates that computes Parity for sorted inputs x ₁≥???≥x _n . For an arbitrary r, we completely determine the minimum size. It is 2n?r?2 for odd n and 2n?r?1 for even n for ?log?₂(n+1)??1≤r≤n/2, and it is ?3n/2??1 for r≥n/2. We also determine the minimum size of an inverter for sorted inputs with at most r NOT gates. It is 4n?3r for ?log?₂(n+1)?≤r≤n. In particular, the negation-limited inverter for sorted inputs due to Fischer, which is a core component in all the known constructions of negation-limited inverters, is shown to have the minimum possible size. Our fairly simple lower bound proofs use gate elimination arguments in a somewhat novel way.     

Generation and distillation of non-Gaussian entanglement from nonclassical photon statistics

With a product state of the form \({{\rho}_{\rm in} = {\rho}_{a} \otimes |0 \rangle_b {_b} \langle 0|}\) as input to a beam splitter, the output two-mode state ρ _out is shown to be negative under partial transpose (NPT) whenever the photon number distribution (PND) statistics { p(n _a ) } associated with the possibly mixed state ρ _a of the input a-mode is antibunched or otherwise nonclassical, i.e., whenever { p(n _a ) } fails to respect any one of an infinite sequence of necessary and sufficient classicality conditions. Negativity under partial transpose turns out to be a necessary and sufficient test for entanglement of ρ _out which is generically non-Gaussian. The output of a PND distribution is further shown to be distillable if any one of an infinite sequence of three term classicality conditions is violated.     

Fabrication of piezoelectric multilayer thin-film actuators

In this study, we fabricated multilayer ceramics (MLCs) composed of multilayered Pb(Zr,Ti)O₃ (PZT) piezoelectric thin films with internal electrodes and evaluated their dielectric and piezoelectric properties. The stack of PZT ferroelectric layers (550 nm) and SrRuO₃ (SRO, 80 nm) electrodes were alternatively deposited on Pt/Ti-coated silicon-on-insulator substrates by radio-frequency magnetron sputtering. The MLCs composed of one, three, and five PZT layers were fabricated by the alternate sputtering deposition of PZT ferroelectric layers and SRO electrodes through the movable shadow mask. The capacitances of MLCs were proportionally increased with the number of PZT layers, while their relative dielectric constants were almost same among the each MLC. The MLCs exhibited symmetric and saturated P–E hysteresis loops similar to the conventional PZT thin films. We estimated that the piezoelectric properties of MLCs by FEM simulation, and confirmed that the effective transverse piezoelectric coefficients (d _31,eff ) increased with the number of PZT layers. The piezoelectric coefficients calculated to be d _31,eff  = ?2964 pC/N at 25 PZT layers, which is much higher than those of conventional single-layer piezoelectric thin films.     

A Priori knowledge based secure payload estimation

Many contemporary steganographic schemes aim to embed fixed-length secret message in the cover while minimizing the stego distortion. However, in some cases, the secret message sender requires to embed a variable-length secret payload within his expected stego security. This kind of problem is named as secure payload estimation (SPE). In this paper, we propose a practical SPE approach for individual cover. The stego security metric we adopt here is the detection error rate of steganalyzer (P _E ). Our method is based on a priori knowledge functions, which are two kinds of functions to be determined before the estimation. The first function is the relation function of detection error rate and stego distortion (P _E ? D function). The second function reflects the relationship between stego distortion and payload rate (D ? α) of the chosen cover. The P _E ? D is the general knowledge, which is calculated from image library. On the other hand, D ? α is for specific cover, which is needed to be determined on site. The estimating procedure is as follows: firstly, the sender solves the distortion D under his expected P _E via P _E ? D, and then calculates the corresponding secure payload α via D ? α of the cover. For on-site operations, the most time-consuming part is calculating D ? α function for cover image, which costs 1 time of STC coding. Besides this, the rest on-site operations are solving single-variable formulas, which can be easily tackled. Our approach is an efficient and practical solution for SPE problem.     

MDS codes in Doob graphs

The Doob graph D(m, n), where m > 0, is a Cartesian product of m copies of the Shrikhande graph and n copies of the complete graph K ₄ on four vertices. The Doob graph D(m, n) is a distance-regular graph with the same parameters as the Hamming graph H(2m + n, 4). We give a characterization of MDS codes in Doob graphs D(m, n) with code distance at least 3. Up to equivalence, there are m ³/36+7m ²/24+11m/12+1?(m mod 2)/8?(m mod 3)/9 MDS codes with code distance 2m + n in D(m, n), two codes with distance 3 in each of D(2, 0) and D(2, 1) and with distance 4 in D(2, 1), and one code with distance 3 in each of D(1, 2) and D(1, 3) and with distance 4 in each of D(1, 3) and D(2, 2).     

H<Subscript>2</Subscript>S nano-jetting through a p–n junction-like graphene/Au nano-injector: a molecular dynamics study

This paper uses fully atomistic molecular dynamics to outline the dynamics of H₂S nano-jetting through a p–n junction-like graphene/Au nano-injector. We examined the effects of nano-injector diameter (d), system temperature (T), and the extrusion velocity (v) of a graphite piston plate on the formation of H₂S nano-jets, system pressure, and the number of molecules (N _m) in the outflow. The combined effects of high critical pressure and small orifice resulted in a larger jet angle, which increased the number of H₂S molecules stuck to the graphene surface at the outlet. Moving the graphite piston plate toward the orifice of the nano-injector increased in the change in momentum and interactive forces between H₂S molecules, resulting in three phases of pressure establishment in the nano-injector: incubation (phase I), steep pressure increase (phase II), and high pressure (phase III). When operated at T ≥ 300 K and v < 27.912 m/s, the proposed nano-jet device is able to produce a well-dispersed spray of H₂S without H₂S molecules sticking to the graphene surface at the outlet. The p–n junction-like Au-doped graphene surface provides an additional energy barrier preventing the transport of electrons from H₂S molecule to the graphene. This inhibits the accumulation of H₂S molecules and subsequent blockages at the exit of the nano-injector. Simulation results demonstrate the potential of using chemiresistive sensing to monitor H₂S flow patterns during nano-jetting. The findings presented in this study could be transformative to the design of nano-injectors for other gases commonly used as biomarkers.     

Demarcating wetting states in textured microchannels under flow conditions by Poiseuille number

In this work, flow friction in microchannels decorated with micropillars was investigated experimentally, with an interest to understand the wetting transition through two simple means: Poiseuille number and scaling laws. Different wetting states were demarcated by qualitatively assessing the behaviour of Poiseuille number (Po = f·Re, where f is friction factor and Re is Reynolds number), which are further corroborated by confocal microscopy-based measurements and numerical simulations. The wetting transition ensued smoothly with an increase in Re, independent of the gas fraction (a ratio of area covered by the liquid–gas interface to the total projected area), for moderate gas fractions, whereas an early breakdown of the Cassie–Baxter state occurred irrespective of Re at high gas fractions. Additionally, the scaling laws were found to correlate well with the underlying state of the flow. Our observations revealed that the liquid–gas interface exhibits a partial slip, contrary to the common notion that it is shear free. It is inferred that an increase in effective flow area leads to a reduction in flow friction in textured microchannels. The present work underlines three important outcomes. The first is the identification of wetting states in flow conditions shown by tracking the Poiseuille number. The second is that the liquid–gas interface is deduced to behave like a partial slip boundary. The third is that a textured microchannel can be worse than an enlarged dimension microchannel.     

On Some Optimization Problems for the Rényi Divergence

We consider the problem of determining the maximum and minimum of the Rényi divergence D_λ(P||Q) and D_λ(Q||P) for two probability distribution P and Q of discrete random variables X and Y provided that the probability distribution P and the parameter α of α-coupling between X and Y are fixed, i.e., provided that Pr{X = Y } = α.     

Effects of topological changes in microchannel geometries on the asymmetric breakup of a droplet

Passive asymmetric breakups of a droplet could be done in many microchannels of various geometries. In order to study the effects of different geometries on the asymmetric breakup of a droplet, four types of asymmetric microchannels with the topological equivalence of geometry are designed, which are T-90, Y-120, Y-150, and I-180 microchannels. A three-dimensional volume of fluid multiphase model is employed to investigate the asymmetric rheological behaviors of a droplet numerically. Three regimes of rheological behaviors as a function of the capillary numbers Ca and the asymmetries As defined by As = (b1 ? b2)/(b1 + b2) (where b1 and b2 are the widths of two asymmetric sidearms) have been observed. A power law model based on three major factors (Ca, As and the initial volume ratio r ₀) is employed to describe the volume ratio of two daughter droplets. The analysis of pressure fields shows that the pressure gradient inside the droplet is one of the major factors causing the droplet translation during its asymmetric breakup. Besides the above similarities among various microchannels, the asymmetric breakup in them also have some slight differences as various geometries have different enhancement or constraint effects on the translation of the droplet and the cutting action of flows. It is disclosed that I-180 microchannel has the smallest critical capillary number, the shortest splitting time, and is hardest to generate satellite droplets.     

Error bounds for complementarity problems with tridiagonal nonlinear functions

In this paper we consider the complementarity problem NCP(f) with f(x) = Mx + φ(x), where M ∈ R ^n×n is a real matrix and φ is a so-called tridiagonal (nonlinear) mapping. This problem occurs, for example, if certain classes of free boundary problems are discretized. We compute error bounds for approximations \({\hat x}\) to a solution x* of the discretized problems. The error bounds are improved by an iterative method and can be made arbitrarily small. The ideas are illustrated by numerical experiments.     

Formation of inverse Chladni patterns in liquids at microscale: roles of acoustic radiation and streaming-induced drag forces

While Chladni patterns in air over vibrating plates at macroscale have been well studied, inverse Chladni patterns in water at microscale have recently been reported. The underlying physics for the focusing of microparticles on the vibrating interface, however, is still unclear. In this paper, we present a quantitative three-dimensional study on the acoustophoretic motion of microparticles on a clamped vibrating circular plate in contact with water with emphasis on the roles of acoustic radiation and streaming-induced drag forces. The numerical simulations show good comparisons with experimental observations and basic theory. While we provide clear demonstrations of three-dimensional particle size-dependent microparticle trajectories in vibrating plate systems, we show that acoustic radiation forces are crucial for the formation of inverse Chladni patterns in liquids on both out-of-plane and in-plane microparticle movements. For out-of-plane microparticle acoustophoresis, out-of-plane acoustic radiation forces are the main driving force in the near-field, which prevent out-of-plane acoustic streaming vortices from dragging particles away from the vibrating interface. For in-plane acoustophoresis on the vibrating interface, acoustic streaming is not the only mechanism that carries microparticles to the vibrating antinodes forming inverse Chladni patterns: In-plane acoustic radiation forces could have a greater contribution. To facilitate the design of lab-on-a-chip devices for a wide range of applications, the effects of many key parameters, including the plate radius R and thickness h and the fluid viscosity μ, on the microparticle acoustophoresis are discussed, which show that the threshold in-plane and out-of-plane particle sizes balanced from the acoustic radiation and streaming-induced drag forces scale linearly with R and \(\sqrt \mu\), but inversely with \(\sqrt h\).     

When LP Is the Cure for Your Matching Woes: Improved Bounds for Stochastic Matchings

Consider a random graph model where each possible edge e is present independently with some probability p _e . Given these probabilities, we want to build a large/heavy matching in the randomly generated graph. However, the only way we can find out whether an edge is present or not is to query it, and if the edge is indeed present in the graph, we are forced to add it to our matching. Further, each vertex i is allowed to be queried at most t _i times. How should we adaptively query the edges to maximize the expected weight of the matching? We consider several matching problems in this general framework (some of which arise in kidney exchanges and online dating, and others arise in modeling online advertisements); we give LP-rounding based constant-factor approximation algorithms for these problems. Our main results are the following:

• We give a 4 approximation for weighted stochastic matching on general graphs, and a 3 approximation on bipartite graphs. This answers an open question from Chen et al. (ICALP’09, LNCS, vol. 5555, pp. 266–278, [2009]).
• We introduce a generalization of the stochastic online matching problem (Feldman et al. in FOCS’09, pp. 117–126, [2009]) that also models preference-uncertainty and timeouts of buyers, and give a constant factor approximation algorithm.

     

Research on axial bearing capacity of rectangular concrete-filled steel tubular columns based on artificial neural networks

Design of rectangular concrete-filled steel tubular (CFT) columns has been a big concern owing to their complex constraint mechanism. Generally, most existing methods are based on simplified mechanical model with limited experimental data, which is not reliable under many conditions, e.g., columns using high strength materials. Artificial neural network (ANN) models have shown the effectiveness to solve complex problems in many areas of civil engineering in recent years. In this paper, ANN models were employed to predict the axial bearing capacity of rectangular CFT columns based on the experimental data. 305 experimental data from articles were collected, and 275 experimental samples were chosen to train the ANN models while 30 experimental samples were used for testing. Based on the comparison among different models, artificial neural network model1 (ANN1) and artificial neural network model2 (ANN2) with a 20-neuron hidden layer were chosen as the fit prediction models. ANN1 has five inputs: the length (D) and width (B) of cross section, the thickness of steel (t), the yield strength of steel (f _y), the cylinder strength of concrete (f′_c). ANN2 has ten inputs: D, B, t, f _y, f′_c, the length to width ratio (D/B), the length to thickness ratio (D/t), the width to thickness ratio (B/t), restraint coefficient (ξ), the steel ratio (α). The axial bearing capacity is the output data for both models.The outputs from ANN1 and ANN2 were verified and compared with those from EC4, ACI, GJB4142 and AISC360-10. The results show that the implemented models have good prediction and generalization capacity. Parametric study was conducted using ANN1 and ANN2 which indicates that effect law of basic parameters of columns on the axial bearing capacity of rectangular CFT columns differs from design codes.The results also provide convincing design reference to rectangular CFT columns.     

Deterministic lateral displacement (DLD) in the high Reynolds number regime: high-throughput and dynamic separation characteristics

Recent progress in the development of biosensors has created a demand for high-throughput sample preparation techniques that can be easily integrated into microfluidic or lab-on-a-chip platforms. One mechanism that may satisfy this demand is deterministic lateral displacement (DLD), which uses hydrodynamic forces to separate particles based on size. Numerous medically relevant cellular organisms, such as circulating tumor cells (10–15 µm) and red blood cells (6–8 µm), can be manipulated using microscale DLD devices. In general, these often-viscous samples require some form of dilution or other treatment prior to microfluidic transport, further increasing the need for high-throughput operation to compensate for the increased sample volume. However, high-throughput DLD devices will require a high flow rate, leading to an increase in Reynolds numbers (Re) much higher than those covered by existing studies for microscale (≤?100 µm) DLD devices. This study characterizes the separation performance for microscale DLD devices in the high-Re regime (10?<?Re?<?60) through numerical simulation and experimental validation. As Re increases, streamlines evolve and microvortices emerge in the wake of the pillars, resulting in a particle trajectory shift within the DLD array. This differs from previous DLD works, in that traditional models only account for streamlines that are characteristic of low-Re flow, with no consideration for the transformation of these streamlines with increasing Re. We have established a trend through numerical modeling, which agrees with our experimental findings, to serve as a guideline for microscale DLD performance in the high-Re regime. Finally, this new phenomenon could be exploited to design passive DLD devices with a dynamic separation range, controlled simply by adjusting the device flow rate.     

Study on TSV isolation liners for a Via Last approach with the use in 3D-WLP for MEMS

This paper discusses approaches for the isolation of deep high aspect ratio through silicon vias (TSV) with respect to a Via Last approach for micro-electro-mechanical systems (MEMS). Selected TSV samples have depths in the range of 170…270 µm and a diameter of 50 µm. The investigations comprise the deposition of different layer stacks by means of subatmospheric and plasma enhanced chemical vapour deposition (PECVD) of tetraethyl orthosilicate; Si(OC₂H₅)₄ (TEOS). Moreover, an etch-back approach and the selective deposition on SiN were also included in the investigations. With respect to the Via Last approach, the contact opening at the TSV bottom by means of a specific spacer-etching method have been addressed within this paper. Step coverage values of up to 74 % were achieved for the best of those approaches. As an alternative to the SiO₂-isolation liners a polymer coating based on the CVD of Parylene F was investigated, which yields even higher step coverage in the range of 80 % at the lower TSV sidewall for a surface film thickness of about 1000 nm. Leakage current measurements were performed and values below 0.1 nA/cm² at 10 kV/cm were determined for the Parylene F films which represents a promising result for the aspired application to Via Last MEMS-TSV.