A numerical study of MHD generalized Couette flow and heat transfer with variable viscosity and electrical conductivity

The steady flow and heat transfer of an electrically conducting fluid with variable viscosity and electrical conductivity between two parallel plates in the presence of a transverse magnetic field is investigated. It is assumed that the flow is driven by combined action of axial pressure gradient and uniform motion of the upper plate. The governing nonlinear equations of momentum and energy transport are solved numerically using a shooting iteration technique together with a sixth-order Runge-Kutta integration algorithm. Solutions are presented in graphical form and given in terms of fluid velocity, fluid temperature, skin friction and heat transfer rate for various parametric values. Our results reveal that the combined effect of magnetic field, viscosity, exponents of variable properties, various fluid and heat transfer dimensionless quantities and the electrical conductivity variation, have significant impact on the hydromagnetic and electrical properties of the fluid.     

Adomian decomposition method for solving the telegraph equation in charged particle transport

In this paper, the analysis for the telegraph equation in case of isotropic small angle scattering from the Boltzmann transport equation for charged particle is presented. The Adomian decomposition is used to solve the telegraph equation. By means of MAPLE the Adomian polynomials of obtained series (ADM) solution have been calculated. The behaviour of the distribution function are shown graphically. The results reported in this article provide further evidence of the usefulness of Adomain decomposition for obtaining solution of linear and nonlinear problems.     

Beyond Adomian method: The variational iteration method for solving heat-like and wave-like equations with variable coefficients

Recently Adomian method was used to solve various kinds of heat-like and wave-like equations. In this Letter, an alternative approach called the variational iteration method is presented to overcome the demerit of complex calculation of Adomian polynomial. Some examples are given to show the reliability and the efficiency of the variational iteration method.     

Adomian decomposition method and Padé approximants for solving the Blaszak--Marciniak lattice

The Adomian decomposition method （ADM） and Pade approximants are combined to solve the well-known Blaszak-Marciniak lattice, which has rich mathematical structures and many important applications in physics and mathematics. In some cases, the truncated series solution of ADM is adequate only in a small region when the exact solution is not reached. To overcome the drawback, the Pade approximants, which have the advantage in turning the polynomials approximation into a rational function, are applied to the series solution to improve the accuracy and enlarge the convergence domain. By using the ADM-Pade technique, the soliton solutions of the Blaszak-Marciniak lattice are constructed with better accuracy and better convergence than by using the ADM alone. Numerical and figurative illustrations show that it is a promising tool for solving nonlinear problems.     

Analytic treatment of linear and nonlinear Schrödinger equations: A study with homotopy-perturbation and Adomian decomposition methods

In this work, two powerful analytical methods, called homotopy-perturbation method (HPM) and Adomian decomposition method (ADM) are introduced to obtain the exact solutions of linear and nonlinear Schrödinger equations. The main objective is to propose alternative methods of solution, which do not require small parameters and avoid linearization and physically unrealistic assumptions. The results show that these methods are very efficient and convenient and can be applied to a large class of problems. The comparison of the methods shows that although the numerical results of these methods are the same, HPM is much easier, more convenient and efficient than ADM.     

用恒流法自制仪器精确测定液体的黏度

介绍了用恒流法自制仪器精确测定液体的黏度.该自制仪器,能有效地保证底部的密封性、液体水位恒定和可操作性.另外,还增设一个玻璃管连通器,并附加一根标尺和游标,可使液体水位的测量精确度提高.在接毛细管输出液体的量筒处加上光电计时,可使测量液体的流速的精确度提高.再将相关数据送入已编好程序的电子仪器进行计算并显示出黏度.另外该仪器可调换不同管径的毛细管,因此应用本方法能快速准确地测量多种液体的黏度.     

一种新的模拟渗流运动的数值方法

根据格子Boltzmann方法及相关理论，建立了一个新的模拟渗流运动的数值模型，所得模型没有在边界上采取相应平均措施，同时还避免了一些非物理副产品的出现-实例计算数值结果与精确解符合较好，证明模型可靠- 关键词： 渗流 格子Boltzmann方法 数值模型     

Stability of hydromagnetic Dean flow between two arbitrarily spaced concentric circular cylinders in the presence of a uniform axial magnetic field

The effect of a uniform axial magnetic field on the stability of the flow of an incompressible viscous electrically conducting fluid between two arbitrarily spaced concentric circular cylinders driven by a constant azimuthal pressure gradient is studied. The linearized stability equations for steady axisymmetric disturbances form an eigenvalue problem, which are solved by using a classical Runge–Kutta scheme combined with a shooting method, termed unit disturbance method. It is observed that for fixed gap width, the magnetic field has a stabilizing influence on the flow for both perfectly conducting and nonconducting walls. It is also found that for a given value of magnetic parameter, stabilization is more as the gap width increases. Further the electrically nonconducting walls are found to be more destabilizing than the perfectly conducting walls. The critical value of the radii ratio (0<η<1) beyond which the first unstable mode becomes nonaxisymmetric is determined for various values of the magnetic parameter.     

Theoretical and numerical study of a thermal convection problem with temperature-dependent viscosity in an infinite layer

A convection problem with temperature-dependent viscosity in an infinite layer is presented. This problem has important applications in mantle convection. The existence of a stationary bifurcation is proven together with a condition to obtain the critical parameters at which the bifurcation takes place. A numerical strategy has been developed to calculate the critical bifurcation curves and the most unstable modes for a general dependence of viscosity on temperature. An exponential dependence of viscosity on temperature has been considered in the numerical calculations. Comparisons with the classic Rayleigh-Bénard problem with constant viscosity indicate that the critical temperature difference threshold decreases as the exponential rate parameter increases. The vertical velocity of the marginal mode exhibits motion concentrated in the region where viscosity is smaller.     

A hybrid numerical method for interfacial fluid flow with soluble surfactant

We address a significant difficulty in the numerical computation of fluid interfaces with soluble surfactant that occurs in the physically representative limit of large bulk Peclet number Pe. At the high values of Pe in typical fluid-surfactant systems, there is a transition layer near the interface in which the surfactant concentration varies rapidly, and large gradients at the interface must be resolved accurately to evaluate the exchange of surfactant between the interface and bulk flow. We use the slenderness of the layer to develop a fast and accurate ‘hybrid’ numerical method that incorporates a separate, singular perturbation analysis of the dynamics in the transition layer into a full numerical solution of the interfacial free boundary problem. The accuracy and efficiency of the method is assessed by comparison with a more ‘traditional’ numerical approach that uses finite differences on a curvilinear coordinate system exterior to the bubble, without the separate transition layer reduction. The traditional method implemented here features a novel fast calculation of fluid velocity off the interface.     

A low-dissipation and time-accurate method for compressible multi-component flow with variable specific heat ratios

A low-dissipation method for calculating multi-component gas dynamics flows with variable specific heat ratio that is capable of accurately simulating flows which contain both high- and low-Mach number features is proposed. The technique combines features from the double-flux multi-component model, nonlinear error-controlled WENO, adaptive TVD slope limiters, rotated Riemann solvers, and adaptive mesh refinement to obtain a method that is both robust and accurate. Success of the technique is demonstrated using an extensive series of numerical experiments including premixed deflagrations, Chapman–Jouget detonations, re-shocked Richtmyer–Meshkov instability, shock-wave and hydrogen gas column interaction, and multi-dimensional detonations. This technique is relatively straight-forward to implement using an existing compressible Navier–Stokes solver based on Godunov’s method.     

变系数瞬态热传导问题边界元格式的特征正交分解降阶方法

提出了一种基于边界元法求解变系数瞬态热传导问题的特征正交分解(POD)降阶方法,重组并推导出变系数瞬态热传导问题适合降阶的边界元离散积分方程,建立了变系数瞬态热传导问题边界元格式的POD降阶模型,并用常数边界条件下建立的瞬态热传导问题的POD降阶模态,对光滑时变边界条件瞬态热传导问题进行降阶分析.首先,对一个变系数瞬态热传导问题,建立其边界域积分方程,并将域积分转换成边界积分;其次,离散并重组积分方程,获得可用于降阶分析的矩阵形式的时间微分方程组;最后,用POD模态矩阵对该时间微分方程组进行降阶处理,建立降阶模型并对其求解.数值算例验证了本文方法的正确性和有效性.研究表明:1)常数边界条件下建立的低阶POD模态矩阵,能够用来准确预测复杂光滑时变边界条件下的温度场结果;2)低阶模型的建立,解决了边界元法中采用时间差分推进技术求解大型时间微分方程组时求解速度慢、算法稳定性差的问题.     

An analytical discrete-ordinates solution of the S-model kinetic equations for flow in a cylindrical tube

An integro-differential form of the linearized S-model kinetic equations for describing flow in a cylindrical tube is projected in such a way as to yield a pair of coupled transport equations that defines the desired velocity and heat-flow profiles. This system is then solved symbolically to yield a pair of coupled integral equations for the physical quantities required. At this point some transformations are carried out to yield a restatement of the original problem in terms of a “pseudo-problem” defined by plane-geometry variables. An analytical version of the discrete-ordinates method is then used to solve the pseudo-problem, and so, after both MATLAB and FORTRAN versions of the developed algorithm are implemented, results thought to be highly accurate are obtained for the case of diffuse reflection from the walls of a cylindrical tube. In addition to the velocity and heat-flow profiles, for the cases of Poiseuille flow and thermal-creep flow, the velocity slips, the heat-flow profiles evaluated at the wall, the particle-flow rates and the heat-flow rates for these two problems are reported for selected values of the tube radius.     

Regularized variable metric method versus the conjugate gradient method in solution of radiative boundary design problem

In this paper, an appropriate distribution of the heating elements’ strengths in a radiation furnace is estimated using inverse methods so that a pre-specified temperature and heat flux distribution is attained on the design surface. Minimization of the sum of the squares of the error function is performed using the variable metric method (VMM), and the results are compared with those obtained by the conjugate gradient method (CGM) established previously in the literature. It is shown via test cases and a well-founded validation procedure that the VMM, when using a “regularized” estimator, is more accurate and is able to reach at a higher quality final solution as compared to the CGM. The test cases used in this study were two-dimensional furnaces filled with an absorbing, emitting, and scattering gas.     

A new domain decomposition method with overlapping patches for ultrascale simulations: Application to biological flows

We address the failure in scalability of large-scale parallel simulations that are based on (semi-)implicit time-stepping and hence on the solution of linear systems on thousands of processors. We develop a general algorithmic framework based on domain decomposition that removes the scalability limitations and leads to optimal allocation of available computational resources. It is a non-intrusive approach as it does not require modification of existing codes. Specifically, we present here a two-stage domain decomposition method for the Navier–Stokes equations that combines features of discontinuous and continuous Galerkin formulations. At the first stage the domain is subdivided into overlapping patches and within each patch a C⁰ spectral element discretization (second stage) is employed. Solution within each patch is obtained separately by applying an efficient parallel solver. Proper inter-patch boundary conditions are developed to provide solution continuity, while a Multilevel Communicating Interface (MCI) is developed to provide efficient communication between the non-overlapping groups of processors of each patch. The overall strong scaling of the method depends on the number of patches and on the scalability of the standard solver within each patch. This dual path to scalability provides great flexibility in balancing accuracy with parallel efficiency. The accuracy of the method has been evaluated in solutions of steady and unsteady 3D flow problems including blood flow in the human intracranial arterial tree. Benchmarks on BlueGene/P, CRAY XT5 and Sun Constellation Linux Cluster have demonstrated good performance on up to 96,000 cores, solving up to 8.21B degrees of freedom in unsteady flow problem. The proposed method is general and can be potentially used with other discretization methods or in other applications.     

The DRESOR method for radiative heat transfer in a one-dimensional medium with variable refractive index

This paper extends the DRESOR (Distribution of Ratios of Energy Scattered by the medium Or Reflected by the boundary surface) method to radiative transfer in a variable refractive index medium. In this method, the intensity is obtained from the source term along the curved integration paths determined only by the variable refractive index, and the DRESOR values are calculated by the Monte Carlo method in which the propagation of the energy bundles are affected by Snell's law. With given temperatures on the black boundaries of a one-dimensional medium, the temperature distribution inside the medium with a variable scattering property is calculated under the condition of radiative equilibrium. It is shown that the DRESOR method has a good accuracy in the cases studied. For an isotropic-scattering medium with the same optical thickness, the scattering albedo has no effect on the temperature distribution, which can be obtained from the general equations and can be seen as an extension of what exists for a constant refractive index; however, the different refractive index causes obvious changes in the temperatures inside the medium. The effect of anisotropic scattering on the temperature distribution cannot be ignored, although it is still weaker than the effect caused by variation in the refractive index.     

Theoretical modeling of the Casimir force-induced instability in freestanding nanowires with circular cross-section

The Casimir force can induce instability and adhesion in freestanding nanostructures. Previous research efforts in this area have exclusively focused on modeling the instability in structures with planar or rectangular cross-section, while, to the best knowledge of the authors, no attention has been paid to investigate this phenomenon for nanowires with circular cross-section. In this study, effects of the Casimir force on the instability and adhesion of freestanding Cylinder–Plate and Cylinder–Cylinder geometries are investigated, which are commonly encountered in real nanodevices. To compute the Casimir force, two approaches, i.e. the proximity force approximation (PFA) for small separations and Dirichlet asymptotic approximation (scattering theory) for large separations, are considered. A continuum mechanics theory is employed, in conjunction with the Euler-beam model, to obtain constitutive equations of the systems. The governing nonlinear constitutive equations of the nanostructures are solved using two different approaches, i.e. the analytical modified Adomian decomposition (MAD) and the numerical finite difference method (FDM). The detachment length and minimum gap, both of which prevent the Casimir force-induced adhesion, are computed for both configurations.     

Effects of variable properties on MHD heat and mass transfer flow near a stagnation point towards a stretching sheet in a porous medium with thermal radiation

The effect of variable viscosity and thermal conductivity on steady magnetohydrodynamic(MHD) heat and mass transfer flow of viscous and incompressible fluid near a stagnation point towards a permeable stretching sheet embedded in a porous medium are presented,taking into account thermal radiation and internal heat genberation/absorbtion.The stretching velocity and the ambient fluid velocity are assumed to vary linearly with the distance from the stagnation point.The Rosseland approximation is used to describe the radiative heat flux in the energy equation.The governing fundamental equations are first transformed into a system of ordinary differential equations using a scaling group of transformations and are solved numerically by using the fourth-order Rung-Kutta method with the shooting technique.A comparison with previously published work has been carried out and the results are found to be in good agreement.The results are analyzed for the effect of different physical parameters,such as the variable viscosity and thermal conductivity,the ratio of free stream velocity to stretching velocity,the magnetic field,the porosity,the radiation and suction/injection on the flow,and the heat and mass transfer characteristics.The results indicate that the inclusion of variable viscosity and thermal conductivity into the fluids of light and medium molecular weight is able to change the boundary-layer behavior for all values of the velocity ratio parameter λ except for λ = 1.In addition,the imposition of fluid suction increases both the rate of heat and mass transfer,whereas fluid injection shows the opposite effect.     

Variable viscosity effects on the flow of MHD hybrid nanofluid containing dust particles over a needle with Hall current—a Xue model exploration

This study scrutinizes the flow of engine oil-based suspended carbon nanotubes magneto-hydrodynamics (MHD) hybrid nanofluid with dust particles over a thin moving needle following the Xue model. The analysis also incorporates the effects of variable viscosity with Hall current. For heat transfer analysis, the effects of the Cattaneo-Christov theory and heat generation/absorption with thermal slip are integrated into the temperature equation. The Tiwari-Das nanofluid model is used to develop the envisioned mathematical model. Using similarity transformation, the governing equations for the flow are translated into ordinary differential equations. The bvp4c method based on Runge-Kutta is used, along with a shooting approach. Graphs are used to examine and depict the consequences of significant parameters on involved profiles. The results revealed that the temperature of the fluid and boundary layer thickness is diminished as the solid volume fraction is raised. Also, with an enhancement in the variable viscosity parameter, the velocity distribution becomes more pronounced. The results are substantiated by assessing them with an available study.     

一种滑移区气体流动的格子Boltzmann曲边界处理新格式

针对滑移区复杂气-固边界存在速度滑移现象,提出了一种基于格子Boltzmann方法的非平衡态外推与有限差分相结合的曲边界处理新格式.该格式具有可考虑实际物理边界与网格线偏移量的优势,较传统half-way DBB(diffusive bounce-back)格式更能准确反映实际边界情况,同时还可获取壁面处气体宏观量及其法向梯度等信息.采用本文所提曲边界处理格式模拟分析了滑移区气体平直/倾斜微通道Poiseuille流、微圆柱绕流和同心微圆柱面旋转Couette流问题.研究结果表明,采用曲边界处理新格式所得结果与理论值以及文献结果符合良好,适用于滑移区气体流动的复杂边界处理,且比half-way DBB格式具有更高的精度,较修正DBB格式具有更好的适应性.