格子Boltzman方法数值模拟微通道内不互溶液液两相流动

利用基于场介子格子Boltzman(LB)方法对T形微通道内油-水不互溶两相流动进行了数值模拟.通过静止液滴算例改进该两相流LB模型且验证了模型的有效性,成功数值模拟了T形微通道内油水两相流的五种不同流型.数值模拟结果和实验数据定性和定量吻合很好,表明本文数值方法的有效性.此外,采用改进后的LB两相流模型模拟了不同界面张力和接触角条件下的液滴形成过程,讨论了这两个变量对流型和液滴形状的影响.结果表明LB方法是研究微通道内不互溶两相流的一种有效的数值工具,能够对各种涉及不互溶多相流微流体设备的设计提供指导.     

Numerical simulation of immiscible liquid-liquid flow in microchannels using lattice Boltzmann method

Immiscible kerosene-water two-phase flows in microchannels connected by a T-junction were numerically studied by a Lattice Boltzmann (LB) method based on field mediators.The two-phase flow lattice Boltzmann model was first validated and improved by several test cases of a still droplet.The five distinct flow regimes of the kerosene-water system,previously identified in the experiments from Zhao et al.,were reproduced.The quantitative and qualitative agreement between the simulations and the experimental dat...     

Numerical simulation of a three-layer porous heat exchanger considering lattice Boltzmann method simulation of fluid flow

Journal of Thermal Analysis and Calorimetry - The present study investigates the thermal characteristics of a proposed porous heat exchanger (PHE). This heat exchanger consists of three sections,...     

Pore-scale modeling of rarefied gas flow in fractal micro-porous media,using lattice Boltzmann method (LBM)

Due to the widespread use of rarefied gas flow in micro-porous media in industrial and engineering problems, a pore-scale modeling of rarefied gas flow through two micro-porous media with fractal geometries is presented, using lattice Boltzmann method. For this purpose, square- and circular-based Sierpinski carpets with fractal geometries are selected due to their inherent behavior for real porous media. Diffusive reflection slip model is used and developed for these porous media through this study. With this respect, the planar Poiseuille flow is selected as a benchmark and validated with the literature. The effect of Knudsen number (Kn) on the permeability is investigated and compared in each geometry. It is shown that as Knudsen number increases, the permeability will increase due to the gas slippage effect on the solid blocks. In addition, it is observed that the permeability is more sensitive to the gaseous flow behavior at the slip and beginning of transition flow regimes. At last, the permeability relationship with Knudsen number is presented with a higher coefficient of determination for both fractal geometries, showing that this relation is logarithmic.

     

Investigation of nanofluids on heat transfer enhancement in a louvered microchannel with lattice Boltzmann method

Numerical studies of laminar forced convective heat transfer and fluid flow in a 2D louvered microchannel with Al₂O₃/water nanofluids are performed by the lattice Boltzmann method (LBM). Eight louvers are arranged in tandem within the single-pass microchannel. The Reynolds number based on channel hydraulic diameter and bulk mean velocity ranges from 100 to 400, where the Al₂O₃ fraction varies from 0 to 4%. A double distribution function approach is adopted for modeling fluid flow and heat transfer. Code validations are performed by comparing the streamwise Nusselt number (Nu) profiles and Fanning friction factors of the present LBM and those of the analytical solutions. Good agreements are obtained. Simulated results show that the louver microstructure can disturb the core flow and guide coolant toward the heated walls, thus enhancing the heat transfer significantly. Furthermore, the addition of nanoparticles in microchannels can also augment the heat transfer, but it creates an unnoticeable pressure loss. With both the louver microstructure and nanofluid, a maximum overall Nu enhancement of 7.06 is found relative to that of the fully developed smooth channel.

     

Investigation of pumping mechanism for non‐Newtonian blood flow with AC electrothermal forces in a microchannel by hybrid boundary element method and immersed boundary‐lattice Boltzmann method

Efficient pumping of blood flow in a microfluidic device is essential for rapid detection of bacterial bloodstream infections (BSI) using alternating current (AC) electrokinetics. Compared with AC electro‐osmosis (ACEO) phenomenon, the advantage of AC electrothermal (ACET) mechanism is its capability of pumping biofluids with high electrical conductivities at a relatively high AC voltage frequency. In the current work, the microfluidic pumping of non‐Newtonian blood flow using ACET forces is investigated in detail by modeling its multi‐physics process with hybrid boundary element method (BEM) and immersed boundary‐lattice Boltzmann method (IB‐LBM). The Carreau–Yasuda model is used to simulate the realistic rheological behavior of blood flow. The ACET pumping efficiency of blood flow is studied in terms of different AC voltage magnitudes and frequencies, thermal boundary conditions of electrodes, electrode configurations, channel height, and the channel length per electrode pair. Besides, the effect of rheological behavior on the blood flow velocity is theoretically analyzed by comparing with the Newtonian fluid flow using scaling law analysis under the same physical conditions. The results indicate that the rheological behavior of blood flow and its frequency‐dependent dielectric property make the pumping phenomenon of blood flow different from that of the common Newtonian aqueous solutions. It is also demonstrated that using a thermally insulated electrode could enhance the pumping efficiency dramatically. Besides, the results conclude that increasing the AC voltage magnitude is a more economical pumping approach than adding the number of electrodes with the same energy consumption when the Joule heating effect is acceptable.     

Inclined surface slip flow of nanoparticles with subject to mixed convection phenomenon: Fractional calculus applications

The interest of researchers towards the nanofluids is noticed in recent years due to leading applications in thermal systems and industrial framework. Referring to such motivations, current study explores the role of velocity slip effects for the mixed convection flow of nanofluid endorsed due to inclined surface. The Casson base fluid model for which the thermal impact needs to be improved. The analysis is observed when the role of velocity slip is important. The modeling of unsteady free convective flow problem yields partial differential system. The Atangana-Baleanu (AB) and Caputo-Fabrizio (CF) fractional operators are implemented in order to simulates the computation of problem. The graphical presentations are prepared in order to check the physical dynamic of parameters.     

Using artificial neural network to optimize the flow and natural heat transfer of a magnetic nanofluid in a square enclosure with a fin on its vertical wall: a lattice Boltzmann simulation

Journal of Thermal Analysis and Calorimetry - The heat transfer of a nanofluid in a square enclosure is numerically simulated using FORTRAN software in this paper, by considering the radiation...     

The Ehrenfest method with quantum corrections to simulate the relaxation of molecules in solution: equilibrium and dynamics

The use of the Ehrenfest method to simulate the relaxation of molecules in solution is explored. Using the cyanide ion dissolved in water as a test model, the independent trajectory (IT) and the bundle of trajectories (BT) approximations are shown to provide very different results for the time evolution of the vibrational populations of the solute. None of these approximations reproduce the Boltzmann equilibrium vibrational populations accurately. A modification of the Ehrenfest method based on the use of quantum correction factors is thus proposed to solve this problem. The simulations carried out using the modified Ehrenfest method provide IT and BT relaxation times which are closer to each other and which agree quite well with previous hybrid perturbative results.     

Optimization and effect of wall conduction on natural convection in a cavity with constant temperature heat source: Using lattice Boltzmann method and neural network algorithm

Journal of Thermal Analysis and Calorimetry - In this paper, the natural convection (NC) of the Al2O3–H2O nanofluids (NFs) in a cavity with a heat source in its center is numerically...     

An impedance analyzer method to simulate the oscillating characteristic of a series piezoelectric sensor in oscillators with zero or non-zero phases

An impedance analyzer method is employed to simulate the oscillation frequency of a series piezoelectric quartz crystal (SPQC) in electrolyte or non-electrolyte solutions. The influence of the oscillator phase on the oscillation frequency and response sensitivity are theoretically derived and experimentally verified. In non-electrolyte liquids, the oscillator phase has little effect on both the oscillation frequency and the response to the permittivity. But in electrolyte solutions, the oscillator phase has a significant influence on the oscillation frequency and the response sensitivity to the conductivity. Depending on the oscillator phase, the oscillation frequency of the SPQC may increase, be maintained or decrease with increasing conductivity in low or high conductive solutions. The dependence of the oscillation frequency of the SPQC on the supply voltage is explained. As an example of the applications, the SPQC is applied to the determination of the critical micelle concentration of ionic surfactants in aqueous solutions.     

Conformational sampling with Poisson–Boltzmann forces and a stochastic dynamics/Monte Carlo method: Application to alanine dipeptide

We apply a combination of stochastic dynamics and Monte Carlo methods (MC/SD) to alanine dipeptide, with solvation forces derived from a Poisson–Boltzmann model supplemented with apolar terms. Our purpose is to study the effects of the model parameters, such as the friction constant and the size of the electrostatic finite difference grid, on the rate of conformational sampling and on the accuracy of the resulting free energy map. For dialanine, a converged Ramachandran map is produced in significantly less time than what is required by stochastic dynamics or Monte Carlo alone. MC/SD is also shown to be faster, per timestep, than explicit methods. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1750–1759, 1997     

A study on the effect of magnetic field and the sinusoidal boundary condition on free convective heat transfer of non-Newtonian power-law fluid in a square enclosure with two constant-temperature obstacles using lattice Boltzmann method

Journal of Thermal Analysis and Calorimetry - In this paper, the flow pattern and thermal characteristics of free convection of a Newtonian magnetohydrodynamic fluid flow inside a square enclosure...     

Determination of verapamil hydrochloride with 12-tungstophosphoric acid by resonance Rayleigh scattering method coupled to flow injection system

A flow injection analysis (FIA) method coupled to resonance Rayleigh scattering (RRS) detection for the determination of verapamil hydrochloride (VP) was proposed. In pH 1.0 acidic medium, 12-tungstophosphoric acid (TP) reacted with VP to form an ion-association complex, which resulted in a significant enhancement of RRS intensity. The maximum scattering peak was located at 293 nm. RRS intensity was proportional to the concentration of VP in the range of 0.017-13.0 μg mL⁻¹, and the detection limit (3σ) was 5.1 ng mL⁻¹. The proposed method exhibited satisfactory reproducibility with a relative standard derivation (R.S.D.) of 2.1% for 11 successive determinations of 5.0 μg mL⁻¹ VP. Therefore, a novel method for the determination of VP by FIA-RRS was developed. The optimum reaction conditions and the parameters of the FIA operation such as flow rate, injection volume, reactor length, and so on had been optimized in this paper. The present method had been applied to the determination of VP in serum samples and pharmaceuticals with satisfactory results. The maximal sample throughput in the optimized system was 80 h⁻¹.     

Novel oxidation reaction of prochlorperazine with bromate in the presence of synergistic activators and its application to trace determination by flow injection/spectrophotometric method

Valve-based comprehensive two-dimensional gas chromatography (GC × GC) is one of the most compact, robust, and inexpensive GC × GC instrument designs. The major drawback of a valve-based modulation configuration lies in diminished detection sensitivity. This loss in sensitivity is because under typical operating conditions the fraction of the first column (i.e., column 1) effluent transferred to the second column (i.e., column 2) is likely to be ∼5-10%. To address this loss in sensitivity, we report the development of a unique total-transfer (i.e., 100%) valve-based GC × GC, without adding complexity to the instrumentation. The new instrument design relies upon simply blocking one of the appropriate ports of the high-speed six-port diaphragm valve that is used as the modulator between columns 1 and 2. The modulation period and difference in head pressure between columns 1 and 2 are found to be the two primary variables that are controlled to provide good detection sensitivity and 100% mass transfer from column 1 to column 2. The detection sensitivity is better with a longer the modulation period. A limit of detection of 0.03 ng/μl was obtained for octane. This sensitive GC × GC configuration is also shown to provide acceptable separation peak capacity, with good separations achieved for real complex samples: gasoline and Eucalyptus oil, where compounds were spread out over much of the two-dimensional separation space. In principle, this total-transfer, valve-based GC × GC is more portable and less expensive than currently available GC × GC instrumentation.     

Absorption, concentration and determination of trace amounts of air pollutants by flow injection method coupled with a chromatomembrane cell system: application to nitrogen dioxide determination

The concentration distribution of an analyte in a chromatomembrane cell (CMC) was examined by using various air samples of different air pollutant (NO₂) concentrations and volumes, and the results obtained could be explained by a proposed principle of the concentration distribution of the analyte in the CMC. This principle was for the first time proved experimentally in the present study. On-line preconcentration and continuous determination of the air pollutant (NO₂) in air samples were realized by coupling a three-hole CMC with a flow injection analysis (FIA) system, where a triethanolamine (TEA) aqueous solution (2 g l⁻¹) was used as an absorbing solution for NO₂ in the air samples. A calibration method with standard nitrite aqueous solutions was developed for the determination of NO₂ in the air samples. Concentrations of NO₂ in indoor air and its diluted air samples were determined by the proposed CMC/FIA method. The volume of air sample necessary for the measurement was decreased to only 5 ml. The measuring time for one sample was about 5–6 min even when a 20 ml air sample was used.     

An online method combining a Gasbench II with continuous flow isotope ratio mass spectrometry to determine the content and isotopic compositions of minor amounts of carbonate in silicate rocks

An online method using continuous flow isotope ratio mass spectrometry (CF‐IRMS) interfaced with a Gasbench II device was established to analyze carbon and oxygen isotopic compositions and to estimate the content of minor amounts of carbonate in silicate rocks. The mixtures of standard materials and high‐purity quartz are firstly used to calibrate different quantities of carbonate in silicates. The results suggest that the accuracy and precision of the online analysis are both better than those obtained using an offline method. There is a positive correlation between the carbonate weight and the Mass44 ion beam intensity (or peak area). When the weight of carbonate in the mixtures is greater than 70 µg (equal to ～1800 mV Mass44 ion beam intensity), the δ¹³C and δ¹⁸O values of samples usually have accuracy and precision of ±0.1‰ and ±0.2‰ (1σ), respectively. If the weight is less than 70 µg, some limitations (e.g., not perfectly linear) are encountered that significantly reduce the accuracy and precision. The measured δ¹⁸O values are systematically lower than the true values by ?0.3 to ?0.7‰; the lower the carbonate content, the lower the measured δ¹⁸O value. For samples with lower carbonate content, the required phosphoric acid doses are higher and more oxygen isotope exchanges with the water in the phosphoric acid. To guarantee accurate results with high precision, multiple analyses of in‐house standards and an artificial MERCK sample with δ¹³C values from ?35.58 to 1.61‰ and δ¹⁸O from 6.04 to 18.96‰ were analyzed simultaneously with the unknown sample. This enables correction of the measured raw data for the natural sample based on multiple‐point normalization. The results indicate that the method can be successfully applied to a range of natural rocks. Copyright © 2010 John Wiley & Sons, Ltd.     

Develop dissipative particle dynamics method to study the fluid flow and heat transfer of Ar and O2 flows in the micro- and nanochannels with precise atomic arrangement versus molecular dynamics approach

Fluid atomic behavior is an important factor for industrial applications. Computer simulations based on simple models predict Poiseuille flow for these atomic structures with the presence of external force. In this work, we describe the dynamical properties of Ar and O₂ flows with precise atomic arrangement via dissipative particle dynamics (DPD) and molecular dynamics (MD) simulation approaches. In these methods, each model is represented by using Large-scale Atomic/Molecular Massively Parallel Simulator package. Simulation results show that maximum rate for velocity of Ar flow in platinum and copper microchannels is 0.100 (unit less)/0.091 Å ps^?1 and 0.121 (unit less)/0.105 Å ps^?1 by using DPD/MD approach. This atomic parameter changes to 0.111 (unit less)/0.102 Å ps^?1 and 0.125 (unit less)/0.108 Å ps^?1 for O₂ fluid with mentioned approaches. By decreasing the microchannel size, the maximum rate of velocity reaches to 0.101 (unit less)/0.099 Å ps^?1 and maximum temperature rate decreases to 485 (unit less)/440 K with DPD/MD approaches. These calculated parameters can be used in industrial application designing for some processes such as heat transfer in structures. It was seen that the developed DPD approach was able to simulate the fluid flow and heat transfer of various types of fluids at micro- and nanoscales with suitable accuracy versus MD.