Ancient Chinese algorithm: The Ying Buzu Shu (method of surplus and deficiency) vs Newton iteration method

Air exploratory discussion of an ancient Chinese algorithm, the Ying Buzu Shu, in about 2nd century BC, known as the rule of double false position in the West is given. In addition to pointing out that the rule of double false position is actually a translation version of the ancient Chinese algorithm, a comparison with well-known Newton iteration method is also made. If derivative is introduced, the ancient Chinese algorithm reduces to the Newton method. A modification of the ancient Chinese algorithm is also proposed, and some of applications to nonlinear oscillators are illustrated.

Four-dimensional Newton's constant in brane world with a sloping extra dimension

We show the four-dimensional Newton's constant obtained naturally from five-dimensional brane world with a tinily sloping extra dimension, which is independent of the bulk Weyl tensor. The corresponding universe is stiff fluid dominated when the slope of extra dimension is very small. Otherwise, the universe may be undergoing a self-acceleration at present epoch and have a decelerated phase in very recent past.

喷涂于运动固壁上的模口流变流动

模口膨胀是聚合物加工中的重要现象，通常用流变效应来解释． 射流自由面上温度分布不均匀，必然会产生表面张力梯度驱动的热毛细对流．采用二维非定常有限元方法，数值模拟喷涂于运动固壁上的流变流体的模口膨胀问题．计算中考虑了流变效应和热毛细效应的耦合作用．结果表明，流变效应和热毛细效应两者均使射流截面增大，非Newton流体的流变效应是模口膨胀的主要原因．

A Study of Composite Beam with Shape Memory Alloy Arbitrarily Embedded under Thermal and Mechanical Loadings

The constitutive relations and kinematic assumptions on the composite beam with shape memory alloy (SMA) arbitrarily embedded are discussed and the results related to the different kinematic assumptions are compared. As the approach of mechanics of materials is to study the composite beam with the SMA layer embedded, the kinematic assumption is vital. In this paper, we systematically study the kinematic assumptions influence on the composite beam deflection and vibration characteristics. Based on the different kinematic assumptions, the equations of equilibrium/motion are different. Here three widely used kinematic assumptions are presented and the equations of equilibrium/motion are derived accordingly. As the three kinematic assumptions change from the simple to the complex one, the governing equations evolve from the linear to the nonlinear ones. For the nonlinear equations of equilibrium, the numerical solution is obtained by using Galerkin discretization method and Newton-Rhapson iteration method. The analysis on the numerical difficulty of using Galerkin method on the post-buckling analysis is presented. For the post-buckling analysis, finite element method is applied to avoid the difficulty due to the singularity occurred in Galerkin method. The natural frequencies of the composite beam with the nonlinear governing equation, which are obtained by directly linearizing the equations and locally linearizing the equations around each equilibrium, are compared. The influences of the SMA layer thickness and the shift from neutral axis on the deflection, buckling and post-buckling are also investigated. This paper presents a very general way to treat thermo-mechanical properties of the composite beam with SMA arbitrarily embedded. The governing equations for each kinematic assumption consist of a third order and a fourth order differential equation with a total of seven boundary conditions. Some previous studies on the SMA layer either ignore the thermal constraint effect or implicitly assume that the SMA is symmetrically embedded. The composite beam with the SMA layer asymmetrically embedded is studied here, in which symmetric embedding is a special case. Based on the different kinematic assumptions, the results are different depending on the deflection magnitude because of the nonlinear hardening effect due to the (large) deflection. And this difference is systematically compared for both the deflection and the natural frequencies. For simple kinematic assumption, the governing equations are linear and analytical solution is available. But as the deflection increases to the large magnitude, the simple kinematic assumption does not really reflect the structural deflection and the complex one must be used. During the systematic comparison of computational results due to the different kinematic assumptions, the application range of the simple kinematic assumption is also evaluated. Besides the equilibrium study of the composite laminate with SMA embedded, the buckling, post-buckling, free and forced vibrations of the composite beam with the different configurations are also studied and compared.

Static Study of Cantilever Beam stiction under Electrostatic Force Influence

The model and analysis of the cantilever beam adhesion problem under the action of electrostatic force are given. Owing to the nonlinearity of electrostatic force, the analytical solution for this kind of problem is not available. In this paper, a systematic method of generating polynomials which are the exact beamsolutions of the loads with different distributions is provided. The polynomials are used to approximate the beam displacement due to electrostatic force. The equilibrium equation offers an answer to how the beam deforms but no information about the unstuck length. The derivative of the functional with respect to the unstuck length offers such information. But to compute the functional it is necessary to know the beam deformation. So the problem is iteratively solved until the results are converged. Galerkin and Newton-Raphson methods are used to solve this nonlinear problem. The effects of dielectric layer thickness and electrostatic voltage on the cantilever beamstiction are studied.The method provided in this paper exhibits good convergence. For the adhesion problem of cantilever beam without electrostatic voltage, the analytical solution is available and is also exactly matched by the computational results given by the method presented in this paper.

Energy Conversion in Shape Memory Alloy Heat Engine Part I: Theory

Shape Memory Alloy (SMA) can be easily deformed to a new shape by applying a small external load at low temperature, and then recovers its original configuration upon heating. This unique shape memory phenomenon has inspired many novel designs. SMA based heat engine is one among them. SMA heat engine is an environment-friendly alternative to extract mechanical energy from low-grade energies, for instance, warm wastewater, geothermal energy, solar thermal energy, etc. The aim of this paper is to present an applicable theoretical model for simulation of SMA-based heat engines. First, a micro-mechanical constitutive model is derived for SMAs. The volume fractions of austenite and martensite variants are chosen as internal variables to describe the evolution of microstructure in SMA upon phase transition. Subsequently, the energy equation is derived based on the first thermodynamic law and the previous SMA model. From Fourier’s law of heat conduction and Newton’s law of cooling, both differential and integral forms of energy conversion equation are obtained.

喷涂于运动边界上液体薄膜射流的流变效应

采用非Newton流体的二阶流体模型分析了相对高温的液体熔体薄膜由模口喷出并涂于运动的固体膜上，讨论了由自由面上温度梯度驱动的非Nweton液体薄膜的热毛细流动，考虑热毛细流动的流变效应，分析是工于润滑理论近似和摄动理论近似，得到了液体高度方程和非Newton液体薄膜的热流体力学过程描述，具体求解了弱流变流体效应的情况。

Rheological flow from a die and painting on a moving solid wall

Die swell is an important, phenomenon. in polymer processing, and is explained usually by rheological properties of the fluid. Because of the nonuniform of temperature distribution on the free surface of the liquid jet, the thermo capillary convection driven by surface tension gradient exists. The rheological fluid flowing out of a die and painting on a moving solid wall is studied by the numerical finite element method of a two-dimensional and unsteady model in the present paper, and both the rheological effect of a non-Newtonian fluid and the thermocapillary effect are considered. The results show that both,effects; will enlarge the cross-section of the fluid jet, and the rheological effect of non-Newtonian fluid dominates the process in general.

A unified model for dislocation nucleation, dislocation emission and dislocation free zone

In this paper, a unified model for dislocation nucleation, emission and dislocation free zone is proposed based on the Peierls framework. Three regions are identified ahead of the crack tip. The emitted dislocations, located away from the crack tip in the form of an inverse pileup, define the plastic zone. Between that zone and the cohesive zone immediately ahead of the crack tip, there is a dislocation free zone. With the stress field and the dislocation density field in the cohesive zone and plastic zone being, respectively, expressed in the first and second Chebyshev polynomial series, and the opening and slip displacements in trigonometric series, a set of nonlinear algebraic equations can be obtained and solved with the Newton-Raphson Method. The results of calculations for pure shearing and combined tension and shear loading after dislocation emission are given in detail. An approximate treatment of the dynamic effects of the dislocation emission is also developed in this paper, and the calculation results are in good agreement with those of molecular dynamics simulations.

Dislocation nucleation and emission from crack-tip

The problems of dislocation nucleation and emission from a crack tip are analysed based on Peierls model. The concept adopted here is essentially the same as that proposed by Rice. A slight modification is introduced here to identify the pure linear elastic response of material. A set of new governing equations is developed, which is different from that used by Beltz and Rice. The stress field and the dislocation density field can be expressed as the first and second Chebyshev polynomial series respectively. Then the opening and slip displacements can be expanded as the trigonometric series. The Newton-Raphson Method is used to solve a set of nonlinear algebraic equations. The new governing equations allow us to extend the analyses to the case of dislocation emission. The calculation results for pure shearing, pure tension and combined tension and shear loading are given in detail.

高速剪切层线性稳定性空间模式参数分析

用一类高精度对称紧致差分格式数值离散一阶改型三维可压粘性扰动方程,对导出的非线性离散特征值问题采用二阶修正Newton-Raphson边值迭代局部解法,实现了超声速剪切流的线性空间稳定性分析。基于对流Mach数、Reynolds数、速度比、密度比等参数研究,讨论了压缩性效应、粘性效应、超声速扰动快/慢模态等,结果显示超声速粘性剪切流的弱不稳定性同多种影响因素密切相关。

微重力流体力学

目录

前言
第一章 微重力流体科学概论
一、微重力科学与微重力流体科学
1、微重力环境
2、重力的影响
3、微重力流体科学的发展
二、微重力流体力学概述
1、对流
2、扩散及输运现象
3、液滴和气泡动力学
4、多相流过程
5、残余重力效应
6、其他流体力学问题
三、微重力物理化学概述
1、临界现象
2、燃烧
3、分散体悬浮系统
4、晶体生长的物理化学问题
四、微重力流体科学的研究途径
1、微重力研究的一般途径
2、微重力实验手段
参考文献
第二章 基本方程组和流体运动特性
一、引言
二、连续性方程和迁移方程
三、动量方程
1、流体的粘性——Reynolds应力
2、动量守恒定律
3、Navier-Stokes方程
四、能量方程
1、总能量方程
2、动能方程
3、内能方程
4、粘性耗散函数
5、Fourier定律及另外二种形式的能量方程
6、不可压流体的导热方程
五、Newton流体的运动方程组及定解条件
1、基本方程组和适定性
2、定解条件
六、Boussinesq近似及适用范围
七、相似律和无量纲参数
1、利用Buckinghan〓定理导出相似参数
2、微重力流体力学的有关物理量和无量纲参数
参考文献
第三章 毛细现象以及界面的平衡和稳定
一、引言
二、表面张力的物理描述
三、液体射流的表面不稳定
1、基本方程组和基态
2、小扰动的线性化方程
3、本征值方程及其解
四、等温条件下液桥的平衡位型和稳定
1、表面张力作用下的平衡条件
2、毛细稳定性
3、旋转稳定性
第五章 液桥的流体动力学稳定理论
1、基本假设和液桥的平衡条件
2、稳定问题的数学提法
3、液桥的Liapunov稳定理论
4、特殊情形（Ω〓=μ=0）以及纯半波不稳定（n=1，m=1）
5、小扰动方程的变分方程
6、小Weber数和大Reynolds数情形的不稳定发展率
7、液桥微重力实验的结果的分析
8、讨论和结论
参考文献
第四章 对流和扩散
一、Pearson对流
1、自由面不变形时的小扰动分析
2、自由面可变形情形
3、非线性理论
4、多层不混溶液体系统
二、热毛细对流
1、矩形容器中的热毛细对流
2、柱形液桥的热毛细对流
3、半浮区液桥热毛细对流的数值模拟
4、薄层液体的热毛细对流
三、热毛细振荡对流的实验研究
1、液桥内部的温度振荡
2、热毛细对流的表面振荡
3、综合测量
四、热毛细对流的振荡机理
1、热流体波不稳定性
2、表面波不稳定性
3、有限高度液桥的线性不稳定性
4、三维不定常数值模拟
5、重力的影响
6、一种非稳定性理论
7、关于振荡的激发机制
参考文献
第五章 液滴动力学
一、等温液滴动力学
1、球形液滴的振荡
2、不混溶液体中球形液滴的振荡
3、弱非线性理论
4、实验模拟
二、非等温液滴的Marangoni迁移
1、定常线性化理论（小Reynolds数，小Marangoni数）
2、非线性理论
3、实验结果
三、液滴和气泡的相互作用
1、双气泡的轴对称理论
2、多液滴的轴对称理论
四、旋转液滴的演化序列和分叉理论
1、旋转液滴的演化
2、旋转液滴的Thomson-Tait稳定准则
3、长期稳定性和动力稳定性
4、长期稳定性真实性的实验证明
5、结论
参考文献
第六章 微重力材料流体力学
一、晶体生长过程
二、纯扩散过程
1、一维扩散过程
2、二维扩散过程
3、固-液界面弯曲对径向分凝的影响
三、浮区晶体生长
1、浮区的热毛细对流
2、浮区的熔质毛细对流
3、浮区对流的振荡特征（小Prandtl数对流）
4、耦合过程
四、溶液晶体生长
1、溶液晶体生长的相变界面过程
2、一维纯扩散过程
3、准定常溶液晶体生长过程
4、不定常溶液生长过程
五、气相晶体生长
1、气相晶体生长过程
2、一维模型
3、物理气相输运中的对流效应
4、化学气相沉积（CVD）过程
参考文献

Parallel Fully-Implicit Computation of MHD Acceleration Experiments

A three-dimensional MHD solver is described in the paper. The solver simulates reacting flows with nonequilibrium between translational-rotational, vibrational and electron translational modes. The conservation equations are discretized with implicit time marching and the second-order modified Steger-Warming scheme, and the resulted linear system is solved iteratively with Newton-Krylov-Schwarz method that is implemented by PETSc package. The results of convergence tests are plotted, which show good scalability and convergence around twice faster when compared with the DPLR method. Then five test runs are conducted simulating the experiments done at the NASA Ames MHD channel, and the calculated pressures, temperatures, electrical conductivity, back EMF, load factors and flow accelerations are shown to agree with the experimental data. Our computation shows that the electrical conductivity distribution is not uniform in the powered section of the MHD channel, and that it is important to include Joule heating in order to calculate the correct conductivity and the MHD acceleration.

Fast and effective analysis methods for fiber Raman amplifiers

A series of new single-step methods and their corresponding algorithms with automatic step size adjustment for model equations of fiber Raman amplifiers are proposed and compared in this paper. On the basis of the Newton-Raphson method, multiple shooting algorithms for the two-point boundary value problems involved in solving Raman amplifier propagation equations are constructed. A verified example shows that, compared with the traditional Runge-Kutta methods, the proposed methods can increase the accuracy by more than two orders of magnitude under the same conditions. The simulations for Raman amplifier propagation equations demonstrate that our methods can increase the computing speed by more than 5 times, extend the step size significantly, and improve the stability in comparison with the Dormand-Prince method. The numerical results show that the combination of the multiple shooting algorithms and the proposed methods has the capacity to rapidly and effectively solve the model equations of multipump Raman amplifiers under various conditions such as co-, counter- and bi-directionally pumped schemes, as well as dual-order pumped schemes.