Lattice Boltzmann method for flow and heat transfer in microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Lattice Boltzmann Method for Flow and Heat Transfer in Microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Analysis of conical-cylinder shells fundamental characteristics for structural health diagnostics

This paper describes the formulation for the free vibration of joined conical-cylindrical shells with uniform thickness using the transfer of influence coefficient for identification of structural characteristics. These characteristics are importance for structural health monitoring to develop model. This method was developed based on successive transmission of dynamic influence coefficients, which were defined as the relationships between the displacement and the force vectors at arbitrary nodal circles of the system. The two edges of the shell having arbitrary boundary conditions are supported by several elastic springs with meridional/axial, circumferential, radial and rotational stiffness, respectively. The governing equations of vibration of a conical shell, including a cylindrical shell, are written as a coupled set of first order differential equations by using the transfer matrix of the shell. Once the transfer matrix of a single component has been determined, the entire structure matrix is obtained by the product of each component matrix and the joining matrix. The natural frequencies and the modes of vibration were calculated numerically for joined conical-cylindrical shells. The validity of the present method is demonstrated through simple numerical examples, and through comparison with the results of previous researchers.

Information theory and one-dimensional disordered conductors

A novel method is proposed to treat the problem of the random resistance of a strictly one-dimensional conductor with static disorder. It is suggested, for the probability distribution of the transfer matrix of the conductor, the distribution of maximum information-entropy, constrained by the following physical requirements: 1) flux conservation, 2) time-reversal invariance and 3) scaling, with the length of the conductor, of the two lowest cumulants of ζ, where = sh2ζ. The preliminary results discussed in the text are in qualitative agreement with those obtained by sophisticated microscopic theories.

Information theory and resistance fluctuations in one-dimensional disordered conductors

A novel method is proposed to treat the problem of the random resistance of a strictly one-dimensional conductor with static disorder. For the probability distribution of the transfer matrix R of the conductor we propose a distribution of maximum information entropy, constrained by the following physical requirements: (1) flux conservation, (2) time-reversal invariance, and (3) scaling with the length of the conductor of the two lowest cumulants of ω, where R=exp(iω→⋅Jbhat). The preliminary results discussed in the text are in qualitative agreement with those obtained by sophisticated microscopic theories.

Novel signal demodulation technique to estimate the amount of chirp in fiber Bragg gratings

A novel detection technique to estimate the amount of chirp in fiber Bragg gratings (FBGs) is proposed. This method is based on the fact that reflectivity at central wavelength of FBG reflection changes with strain/temperature gradient (linear chirp) applied to the same. Transfer matrix approach was used to vary different grating parameters (length, strength and apodization) to optimize variation of reflectivity with linear chirp. Analysis is done for different sets of `FBG length-refractive index strength' combinations for which reflectivity vary linearly with linear chirp over a decent measurement range. This article acts as a guideline to choose appropriate grating parameters in designing sensing apparatus based on change in reflectivity at central wavelength of FBG reflection.

Transverse plane wave analysis of short elliptical chamber mufflers: An analytical approach

Short elliptical chamber mufflers are used often in the modern day automotive exhaust systems. The acoustic analysis of such short chamber mufflers is facilitated by considering a transverse plane wave propagation model along the major axis up to the low frequency limit. The one dimensional differential equation governing the transverse plane wave propagation in such short chambers is solved using the segmentation approaches which are inherently numerical schemes, wherein the transfer matrix relating the upstream state variables to the downstream variables is obtained. Analytical solution of the transverse plane wave model used to analyze such short chambers has not been reported in the literature so far. This present work is thus an attempt to fill up this lacuna, whereby Frobenius solution of the differential equation governing the transverse plane wave propagation is obtained. By taking a sufficient number of terms of the infinite series, an approximate analytical solution so obtained shows good convergence up to about 1300 Hz and also covers most of the range of muffler dimensions used in practice. The transmission loss (TL) performance of the muffler configurations computed by this analytical approach agrees excellently with that computed by the Matrizant approach used earlier by the authors, thereby offering a faster and more elegant alternate method to analyze short elliptical muffler configurations. (C) 2010 Elsevier Ltd. All rights reserved.

Noise reduction with perforated three-duct muffler components

This paper describes the authors’ distributed parameter approach for derivation of closed-form expressions for the four-pole parameters of the perforated three-duct muffler components. In this method, three simultaneous second-order partial differential equations are first reduced to a set of six first-order ordinary differential equations. These equations are then uncoupled by means of a modal matrix. The resulting 6 × 6 matrix is reduced to the 2 × 2 transfer matrix using the relevant boundary conditions. This is combined with transfer matrices of other elements (upstream and downstream of this perforated element) to predict muffler performance like noise reduction, which is also measured. The correlation between experimental and theoretical values of noise reduction is shown to be satisfactory.

Study of Low Temperature Magnetic Properties of a Single Chain Magnet with Alternate Isotropic and Non-collinear Anisotropic Units

Here we study thermodynamic properties of an important class of single-chain magnets (SCMs), where alternate units are isotropic and anisotropic with anisotropy axes being non-collinear. This class of SCMs shows slow relaxation at low temperatures which results from the interplay of two different relaxation mechanisms, namely dynamical and thermal. Here anisotropy is assumed to be large and negative, as a result, anisotropic units behave like canted spins at low temperatures; but even then simple Ising-type model does not capture the essential physics of the system due to quantum mechanical nature of the isotropic units. We here show how statistical behavior of this class of SCMs can be studied using a transfer matrix (TM) method. We also, for the first time, discuss in detail how weak inter-chain interactions can be treated by a TM method. The finite size effect is also discussed which becomes important for low temperature dynamics. At the end of this paper, we apply this technique to study a real helical chain magnet.

长周期莫尔光栅的理论研究

长周期光纤光栅作为一种损耗型的带阻滤波器，在光纤通信和光纤传感器领域有着广泛的应用。结合长周期光纤光栅的传输特性，利用传输矩阵法对多种长周期莫尔光栅的透射谱进行了理论模拟。分析结果表明，莫尔光栅技术在制作高性能的相移长周期光栅和变迹长周期光栅上有着很好的应用前景。

Ultra-narrow bandwidth resonant reflection grating filters using the second diffracted orders

In this paper, we design resonant reflection grating filters employing the second diffracted orders as the leaky modes, then analyze the bandwidth of the reflection peak and the electric field distributions inside the wavegude under resonance. The numeric calculation confirms that ultra-narrow resonant reflection peaks can be observed in these structures. At the same time, strong electric field enhancement appears under resonance. It provides a new approach to diversify the resonant reflection filters and may open a new way to the realization of ultra-narrow bandwidth filters. (C) 2008 Elsevier B.V. All rights reserved.

DTADH and quantum critical phenomena caused by anisotropy and external magnetic field for spin-1/2 Heisenberg diamond chains

Using the transfer matrix renormalization group (TMRG) method, we study the connection between the first derivative of the thermal average of driving-term Hamiltonian (DTADH) and the trace of quantum critical behaviors at finite temperatures. Connecting with the exact diagonalization method, we give the phase diagrams and analyze the properties of each phase for both the ferromagnetic and anti-ferromagnetic frustrated J(3) anisotropy diamond chain models. The finite-temperature scaling behaviors near the critical regions are also investigated. Further, we show the critical behaviors driven by external magnetic field, analyze the formation of the 1/3 magnetic plateau and the influence of different interactions on those critical points for both the ferrimagnetic and anti-ferromagnetic distorted diamond chains.

Transmission properties of electron in quantum rings

We investigated the transmission probability of a single electron transmission through a quantum ring device based on the single-band effective mass approximation method and transfer matrix theory. The time-dependent Schrodinger equation is applied on a Gaussian wave packet passing through the quantum ring system. The electron tunneling resonance peaks split when the electron transmits through a double quantum ring. The splitting energy increases as the distance between the two quantum rings decreases. We studied the tunneling time through the single electron transmission quantum ring from the temporal evolution of the Gaussian wave packet. The electron probability density is sensitive to the thickness of the barrier between the two quantum rings. (C) 2008 American Institute of Physics.

Thermodynamic properties of Cu-3(CO3)(2)(OH)(2) and the spin-1/2 distorted diamond Heisenberg chain

The thermodynamic properties of the spin-1/2 diamond quantum Heisenberg chain model have been investigated by means of the transfer matrix renormalization group (TMRG) method. Considering different crystal structures, by changing the interactions among different spins and the external magnetic fields, we first investigate the magnetic susceptibility, magnetization, and specific heat of the distorted diamond chain as a model of ferrimagnetic spin systems. The susceptibility and the specific heat show different features for different ferromagnetic (F) and antiferromagnetic (AF) interactions and different magnetic fields. A 1/3 magnetization plateau is observed at low temperature in a magnetization curve. Then, we discuss the theoretical mechanism of the double-peak structure of the magnetic susceptibility and the three-peak structure of the specific heat of the compound Cu-3(CO3)(2)(OH)(2), on which an elegant measurement was performed by Kikuchi [Phys. Rev. Lett. 94, 227201 (2005)]. Our computed results are consistent with the main characteristics of the experimental data. Meanwhile, we find that the double-peak structure of susceptibility can be found in several different kinds of spin interactions in the diamond chain. Moreover, a three-peak behavior is observed in the TMRG results of magnetic susceptibility. In addition, we perform calculations relevant for some experiments and explain the characteristics of these materials. (c) 2007 American Institute of Physics.

Investigation of mode characteristics for microdisk resonators by S-matrix and three-dimensional finite-difference time-domain technique

The mode characteristics of a three-dimensional (3D) microdisk with a vertical refractive index distribution of n(2)/3.4/n(2) are investigated by the S-matrix method and 3D finite-difference time-domain (FDTD) technique. For the microdisk with a thickness of 0.2 mu m. and a radius of 1 mu m, the mode wavelengths and quality factors for the HE7,1 mode obtained by 3D FDTD simulation and the S-matrix method are in good agreement as n(2) increases from 1.0 to 2.6. But the Q factor obtained by the 3D FDTD rapidly decreases from 1.12 X 10(4) to 379 as n2 increases from 2.65 to 2.8 owing to the vertical radiation losses, which cannot be predicted by the proposed S-matrix method. The comparisons also show that quality factors obtained from the analytical solution of two-dimensional microdisks under the effective index approximation are five to seven times smaller than those of the 3D FDTD as n(2) = 1 and R = 1 mu m. (c) 2006 Optical Society of America.