Viscous-slip,thermal-slip,and temperature-jump coefficients based on the linearized Boltzmann equation (and five kinetic models) with the Cercignani–Lampis boundary condition

A polynomial expansion procedure and the ADO (analytical discrete-ordinates) method are used to compute the viscous-slip coefficient, the thermal-slip coefficient, and the temperature-jump coefficient from the linearized Boltzmann equation (LBE) for rigid-sphere interactions and the Cercignani–Lampis (CL) boundary condition. These same quantities are also computed from five kinetic models, with the CL condition, and compared to the LBE result. Equivalent results for the LBE and the kinetic models, all based on the usual Maxwell boundary condition, are also reported.     

The viscous-slip, diffusion-slip, and thermal-creep problems for a binary mixture of rigid spheres described by the linearized Boltzmann equation

An analytical version of the discrete-ordinates method (the ADO method) is used with recently established analytical expressions for the rigid-sphere scattering kernels to develop concise and particularly accurate solutions to the viscous-slip, the diffusion-slip, and the half-space thermal-creep problems for a binary gas mixture described by the linearized Boltzmann equation. In addition to a computation of the viscous-slip, the diffusion-slip, and the thermal-slip coefficients, for the case of Maxwell boundary conditions for each of the two species, the velocity, heat-flow, and shear-stress profiles are established for each species of particles. Numerical results are reported for two binary mixtures (Ne–Ar and He–Xe) with various molar concentrations.     

Dynamical analysis of Mathieu equation with two kinds of van der Pol fractional-order terms

In this paper the dynamics of Mathieu equation with two kinds of van der Pol (VDP) fractional-order terms is investigated. The approximately analytical solution is obtained by the averaging method. The steady-state solution, existence conditions and stability condition for the steady-state solution are presented, and it is found that the two kinds of VDP fractional coefficients and fractional orders remarkably affect the steady-state solution, which is characterized by the additional damping coefficient (ADC) and additional stiffness coefficient (ASC). The comparisons between the analytical and numerical solutions verify the correctness and satisfactory precision of the approximately analytical solution. The presented typical amplitude–frequency curves illustrate the important effects of two kinds of VDP fractional-order terms on system dynamics. The application of two VDP fractional-order terms in vibration control is discussed. At last, the detailed results are summarized and the conclusions are made.     

Primary resonance of fractional-order van der Pol oscillator

In this paper the primary resonance of van der Pol (VDP) oscillator with fractional-order derivative is studied analytically and numerically. At first the approximately analytical solution is obtained by the averaging method, and it is found that the fractional-order derivative could affect the dynamical properties of VDP oscillator, which is characterized by the equivalent damping coefficient and the equivalent stiffness coefficient. Moreover, the amplitude–frequency equation for steady-state solution is established, and the corresponding stability condition is also presented based on Lyapunov theory. Then, the comparisons of several different amplitude–frequency curves by the approximately analytical solution and the numerical integration are fulfilled, and the results certify the correctness and satisfactory precision of the approximately analytical solution. At last, the effects of the two fractional parameters, i.e., the fractional coefficient and the fractional order, on the amplitude–frequency curves are investigated for some typical excitation amplitudes, which are different from the traditional integer-order VDP oscillator.     

Couette flow of a binary mixture of rigid-sphere gases described by the linearized Boltzmann equation

A concise and accurate solution to the problem of plane Couette flow for a binary mixture of rigid-sphere gases described by the linearized Boltzmann equation and general (specular-diffuse) Maxwell boundary conditions for each of the two species of gas particles is developed. An analytical version of the discrete-ordinates method is used to establish the velocity, heat-flow, and shear-stress profiles for both types of particles, as well as the particle-flow and heat-flow rates associated with each of the two species. Accurate numerical results are given for the case of a mixture of helium and argon confined between molybdenum and tantalum plates.     

A NUMERICAL STUDY OF THE FLOW OVER ELLIPSOIDAL OBJECTS INSIDE A CYLINDRICAL TUBE

A numerical scheme is developed to obtain the flow field around one, two and five ellipsoidal objects inside a cylindrical tube. The scheme uses the Galerkin finite element technique and the primitive variable(u–v–p) formulation. The two-dimensional incompressible Navier–Stokes equations are solved numerically by using the direct mixed interpolation method. A Picard iteration scheme is used for the solution of the resulting system of non-linear algebraic equations. The computer code is verified by checking with known analytical solutions for the flow past a sphere. Results for the shear stress distributions along the ellipsoids, forces and drag coefficients are obtained for different geometric ratios and Reynolds numbers. Some of the intermediate computational results on the velocity fields developed are also reported.     

Riemann problem for one-dimensional binary gas enhanced coalbed methane process

With an extended Langmuir isotherm, a Riemann problem for one-dimensional binary gas enhanced coalbed methane (ECBM) process is investigated. A new analytical solution to the Riemann problem, based on the method of characteristics, is developed by introducing a gas selectivity ratio representing the gas relative sorption affinity. The influence of gas selectivity ratio on the enhanced coalbed methane processes is identified.     

Nonlinear fracture of 2D magnetoelectroelastic media: Analytical and numerical solutions

A strip electric–magnetic polarization saturation (SEMPS) model is developed to study the electric and magnetic yielding effects on a crack in magnetoelectroelastic (MEE) media. In this model, the MEE medium is treated as being mechanically brittle, and electrically and magnetically ductile. Analogously to the classic Dugdale model, the electric and magnetic yielding zones in front of the crack are represented for simplicity by two strips. In the electric yielding strip the electric displacement equals the electric displacement saturation and meanwhile in the magnetic yielding zone the magnetic induction equals the magnetic induction saturation. The nonlinear analytical solution of this SEMPS model of crack in an infinite MEE medium is obtained using an integral equation approach. The equivalence between the proposed SEMPS model and the existing strip electric–magnetic breakdown (SEMB) model is demonstrated.To analyze the nonlinear fracture problem in the corresponding finite MEE media, the non-linear hybrid extended displacement discontinuity-fundamental solution (NLHEDD-FS) method is modified, and a multiple iteration approach is adapted to determine the electric and magnetic yielding zones. Comparing with the analytical solution, the applicability and effectiveness of the NLHEDD-FS method is verified. Numerical results based on the SEMPS model for a center crack in infinite and finite MEE strip are presented.     

Simulation of two‐phase flow with moving immersed boundaries

A two‐dimensional multi‐phase model for immiscible binary fluid flow including moving immersed objects is presented. The fluid motion is described by the incompressible Navier–Stokes equation coupled with a phase‐field model based on van der Waals' free energy density and the Cahn–Hilliard equation. A new phase‐field boundary condition was implemented with minimization of the free energy in a direct way, to specifically improve the physical behavior of the contact line dynamics for moving immersed objects. Numerical stability and execution time were significantly improved by the use of the new boundary condition. Convergence toward the analytical solution was demonstrated for equilibrium contact angle, the Lucas–Washburn theory and Stefan's problem. The proposed model may be used for multi‐phase flow problems with moving boundaries of complex geometry, such as the penetration of fluid into a deformable, porous medium. Copyright © 2010 John Wiley & Sons, Ltd.     

Dynamic instability of vibrating carbon nanotubes near small layers of graphite sheets based on nonlocal continuum elasticity

This article presents a new asymptotic method to predict dynamic pull-in instability of nonlocal clamped–clamped carbon nanotubes (CNTs) near graphite sheets. Nonlinear governing equations of carbon nanotubes actuated by an electric field are derived. With due allowance for the van der Waals effects, the pull-in instability and the natural frequency–amplitude relationship are investigated by a powerful analytical method, namely, the parameter expansion method. It is demonstrated that retaining two terms in series expansions is sufficient to produce an acceptable solution. The obtained results from numerical methods verify the strength of the analytical procedure. The qualitative analysis of system dynamics shows that the equilibrium points of the autonomous system include center points and unstable saddle points. The phase portraits of the carbon nanotube actuator exhibit periodic and homoclinic orbits.     

An efficient method to solve the strongly coupled nonlinear differential equations of impact dampers

In this paper, ongoing studies to solve nonlinear differential equations are extended by combining the Newmark-beta integration method and the piecewise linearization approach. The discussed method is illustrated with a practical example. In doing so, the coupled nonlinear differential equations of an impact oscillator, which incorporates the Hertzian contact, are derived. To investigate this problem, an object-oriented computer code, based on the presented method, is written in MATLAB. Furthermore, the discussed problem is solved numerically using the Runge–Kutta commercial code. To verify the calculated results, the contact durations, which are obtained using the discussed methods, are compared with the previous analytical results. In this study, accuracy of solution and the process time (cost) are selected as two main parameters of the solution method. The so-called adequacy factor is presented to combine the two main parameters of solution. Finally, it is shown that in the case of Hertzian contact, the presented method can be more adequate than the Runge–Kutta method.     

Approximate analytical solution in slow-fast system based on modified multi-scale method

A simple, yet accurate modified multi-scale method (MMSM) for an approximately analytical solution in nonlinear oscillators with two time scales under forced harmonic excitation is proposed. This method depends on the classical multi-scale method (MSM) and the method of variation of parameters. Assuming that the forced excitation is a constant, one could easily obtain the approximate analytical solution of the simplified system based on the traditional MSM. Then, this solution for the oscillator under forced harmonic excitation could be established after replacing the harmonic excitation by the constant excitation. To certify the correctness and precision of the proposed analytical method, the van der Pol system with two scales subject to slowly periodic excitation is investigated; this system presents rich dynamical phenomena such as spiking (SP), spiking-quiescence (SP-QS), and quiescence (QS) responses. The approximate analytical expressions of the three types of responses are given by the MMSM, and it can be found that the precision of the new analytical method is higher than that of the classical MSM and better than that of the harmonic balance method (HBM). The results obtained by the present method are considerably better than those obtained by traditional methods, quantitatively and qualitatively, particularly when the excitation frequency is far less than the natural frequency of the system.     

Approximate analytical solution in slow-fast system based on modi?ed multi-scale method

A simple, yet accurate modi?ed multi-scale method(MMSM) for an approximately analytical solution in nonlinear oscillators with two time scales under forced harmonic excitation is proposed. This method depends on the classical multi-scale method(MSM) and the method of variation of parameters. Assuming that the forced excitation is a constant, one could easily obtain the approximate analytical solution of the simpli?ed system based on the traditional MSM. Then, this solution for the oscillator under forced harmonic excitation could be established after replacing the harmonic excitation by the constant excitation. To certify the correctness and precision of the proposed analytical method, the van der Pol system with two scales subject to slowly periodic excitation is investigated; this system presents rich dynamical phenomena such as spiking(SP),spiking-quiescence(SP-QS), and quiescence(QS) responses. The approximate analytical expressions of the three types of responses are given by the MMSM, and it can be found that the precision of the new analytical method is higher than that of the classical MSM and better than that of the harmonic balance method(HBM). The results obtained by the present method are considerably better than those obtained by traditional methods,quantitatively and qualitatively, particularly when the excitation frequency is far less than the natural frequency of the system.     

含有非圆形双孔的无限平板中应力的解析解研究

无限平板中含有任意形状单个孔的问题可以使用复变函数方法获得其应力解析解.对于无限平板中含有两个圆孔或两个椭圆孔的双连通域问题,也可以利用多种方法进行求解,比如双极坐标法、应力函数法、复变函数法以及施瓦茨交替法等.其中复变函数中的保角变换方法是获得应力解析解的一个重要方法.但目前尚未见到用此方法求解无限板中含有一个正方形孔和一个椭圆孔的问题.当板在无穷远处受有均布载荷和孔边作用垂直均布压力时,利用保角变换方法可以求解板中含有两个特定形状孔的问题.该方法将所讨论的区域映射成象平面里的一个圆环,其中最关键的一步是找出相应的映射函数.基于黎曼映射定理,提出了该映射函数一般形式,并利用最优化方法,找到了该问题的具体映射函数,然后通过孔边应力边界条件建立了求解两个解析函数的基本方程,获得了该问题的应力解析解.运用ANSYS有限单元法与结果进行了对比.研究了孔距、椭圆形孔大小和两孔布置方位对边界切向应力的影响,以及不同载荷下两孔中心线上应力分布规律.      

Wet/dry areas method for interfacial (free surface) flow

In this paper, a new numerical method is developed for two‐dimensional interfacial (free surface) flows, based on the control volume method and conservative integral form of the Navier–Stokes equations with a standard staggered grid. The new method deploys two continuity equations, the continuity equation of the mass conservation for better convergence of the implicit scheme and the continuity equation of the volume conservation for the equation of pressure correction. The convection terms (the total momentum flux) on the surfaces of control volume are accurately calculated from the wet area exposed to the water, and the dry area exposed to the air. The numerical results produced by the new numerical method agree very well with the analytical solution, experimental images and experimentally measured velocity. Copyright © 2012 John Wiley & Sons, Ltd.     

Three-dimensional analysis of transient thermal stresses in functionally graded plates

An analytical solution is presented for three-dimensional thermomechanical deformations of a simply supported functionally graded (FG) rectangular plate subjected to time-dependent thermal loads on its top and/or bottom surfaces. Material properties are taken to be analytical functions of the thickness coordinate. The uncoupled quasi-static linear thermoelasticity theory is adopted in which the change in temperature, if any, due to deformations is neglected. A temperature function that identically satisfies thermal boundary conditions at the edges and the Laplace transformation technique are used to reduce equations governing the transient heat conduction to an ordinary differential equation (ODE) in the thickness coordinate which is solved by the power series method. Next, the elasticity problem for the simply supported plate for each instantaneous temperature distribution is analyzed by using displacement functions that identically satisfy boundary conditions at the edges. The resulting coupled ODEs with variable coefficients are also solved by the power series method. The analytical solution is applicable to a plate of arbitrary thickness. Results are given for two-constituent metal-ceramic FG rectangular plates with a power-law through-the-thickness variation of the volume fraction of the constituents. The effective elastic moduli at a point are determined by either the Mori–Tanaka or the self-consistent scheme. The transient temperature, displacements, and thermal stresses at several critical locations are presented for plates subjected to either time-dependent temperature or heat flux prescribed on the top surface. Results are also given for various volume fractions of the two constituents, volume fraction profiles and the two homogenization schemes.     

A viscous vortex particle method for deforming bodies with application to biolocomotion

Bio‐inspired mechanics of locomotion generally consist of the interaction of flexible structures with the surrounding fluid to generate propulsive forces. In this work, we extend, for the first time, the viscous vortex particle method (VVPM) to continuously deforming two‐dimensional bodies. The VVPM is a high‐fidelity Navier–Stokes computational method that captures the fluid motion through evolution of vorticity‐bearing computational particles. The kinematics of the deforming body surface are accounted for via a surface integral in the Biot–Savart velocity. The spurious slip velocity in each time step is removed by computing an equivalent vortex sheet and allowing it to flux to adjacent particles; hence, no‐slip boundary conditions are enforced. Particles of both uniform and variable size are utilized, and their relative merits are considered. The placement of this method in the larger class of immersed boundary methods is explored. Validation of the method is carried out on the problem of a periodically deforming circular cylinder immersed in a stagnant fluid, for which an analytical solution exists when the deformations are small. We show that the computed vorticity and velocity of this motion are both in excellent agreement with the analytical solution. Finally, we explore the fluid dynamics of a simple fish‐like shape undergoing undulatory motion when immersed in a uniform free stream, to demonstrate the application of the method to investigations of biomorphic locomotion. Copyright © 2008 John Wiley & Sons, Ltd.     

Half-space problem of the nonlinear Boltzmann equation for weak evaporation and condensation of a binary mixture of vapors

The half-space problem of evaporation and condensation of a binary mixture of vapors is investigated on the basis of the kinetic theory of gases. Assuming the Mach number of the normal component of the flow is small, a solution of the Boltzmann equation that varies slowly in the scale of the molecular mean-free-path (slowly varying solution) is introduced. Then a fluid-dynamic system that describes the behavior of the slowly varying solution is derived by a systematic asymptotic analysis. The analytical expression of the conditions allowing steady evaporation or condensation is derived from that system. We analyze the qualitative difference between the conditions in the evaporation and condensation cases: four conditions are needed in the former case while only one condition is required in the latter case. The present paper extends a earlier contribution of the first author for the BGK-type model equation [S. Takata, Half-space problem of weak evaporation and condensation of a binary mixture of vapors, in: Capitelli M. (Ed.), Rarefied Gas Dynamics, AIP, New York, 2005, pp. 503–508] to the Boltzmann equation. The extension is achieved by considering the linear stability of the far field in the case of evaporation and the H theorem, the monotonic decrease of the flux of Boltzmann's H function, in the case of condensation.     

Modeling particle sedimentation in viscous fluids with a coupled immersed boundary method and discrete element method

Numerical techniques have increasingly been used to model fluid–particle two-phase flows. Coupling the immersed boundary method (IBM) and discrete element method (DEM) is one promising approach for modeling particulate flows. In this study, IBM was coupled with DEM to improve the reliability and accuracy of IBM for determining the positions of particles during the sedimentation process within viscous fluids. The required ratio of the particle diameter to the grid size (D/dx) was determined by comparing the simulation results with the analytical solution and experimental data. A dynamic mesh refinement model was utilised in the IBM model to refine the computational fluid dynamics grid near the particles. In addition, an optimum coupling interval between the IBM and DEM models was determined based on the experimental results of a single particle sedimentation within silicon oil at a Reynolds number of 1.5. The experimental results and the analytical solution were then utilised to validate the IBM–DEM model at Reynolds numbers of 4.1, 11.6, and 31.9. Finally, the validated model was utilised to investigate the sedimentation process for more than one particle by modeling the drafting-kissing-tumbling process and the Boycott phenomenon. Benchmark tests showed that the IBM–DEM technique preserves the advantages of DEM for tracking a group of particles, while the IBM provides a reliable and accurate approach for modeling the particle–fluid interaction.     

Application of solution structure theorems to Cattaneo–Vernotte heat conduction equation with non-homogeneous boundary conditions

In this study, a non-Fourier heat conduction problem formulated using the Cattaneo–Vernotte (C–V) model with non-homogeneous boundary conditions is solved with the superposition principle in conjunction with solution structure theorems. It is well known that the aforementioned analytical method is not suitable for such a class of thermal problems. However, by performing a functional transformation, the original non-homogeneous partial differential equation governing the physical problem can be cast into a new form such that it consists of a homogeneous part and an additional auxiliary function. As a result, the modified homogeneous governing equation can then be solved with solution structure theorems for temperatures inside a finite planar medium. The methodology provides a convenient, accurate, and efficient solution to the C–V heat conduction equation with non-homogeneous boundary conditions.