非牛顿流体内流传热研究

在本文中,研究了注入轴对称模腔非牛顿流体非定常流动.本文的第二部份研究了上随体Maxwell流体管内热流动.对于注入模腔流动.其本构方程采用幂律流体模型方程.为了避免在表现粘度中温度关系引起的非线性.引进了一特征粘度的概念.描述本力学过程的基本方程是,本构方程、定常状态的运动方程、非定常能量方程及连续方程.该方程组在空间是二维问题,在数学上是三维问题.采用分裂差分格式求得本方程组的数值解答.分裂法曾成功应用于求解牛顿流体问题.在本文中,首次将分裂法成功地应用解决非牛顿流体流动问题.对于圆管内热流,给出了差分格式,使基本方程组化为一个三对角方程组.其结果,给出了不同时刻的模腔内二维温度分布.     

第五届北京高中数学知识应用竞赛初赛试题及参考答案

1.用宽度为 0 .3m的带子缠绕包扎圆柱型管道 ,管道长 30 m,截面周长为 0 .5 m .(1)如果用带子全部包住管道 ,最少要用多长带子 .请你给出计算这个最小长度的公式 ,并且依此计算出所需长度数值 .(2 )现有一条长度为 5 1m的带子 ,想将这条带子全部用于缠绕包扎这个管道 ,可以使带子的接缝处重叠褡接 .请你给出用这条带子缠绕包扎这个管道的方案 .(计算结果精确到 0 .0 0 1m)解　 (1)记圆管截面周长为 C,圆管长度为 L ,带子的宽度为 W,带子的长度为 M.则使用带子缠绕包扎管道时 ,如果带子全部包住管道且互不重叠 ,它们应该满足如下关系WM - …     

基于FFT的区域变换对数均匀跳扩散模型期权定价

本文考虑具有区域变换跳跃幅度服从对数均匀分布的跳扩散模型的期权定价问题.本文给出了这样模型的期权定价方法和计算过程,当中采用了FFT(快速傅里叶变换法),最后给出了数值计算结果.     

发展流动与形成流动

本文着重研究振荡状态下的管流问题.从基本的Navier-Stokes方程出发,进行线性化之后采用虚宗量Bessel函数的方法求解,得出了一组描述圆管振荡状态下发展流动的速度分布、压力分布公式.它们较之Atabek等人推导的公式更为简明和易于计算.并且,在相同的条件下简化为形成流动的表达式后,两者的公式完全一致.数值计算还表明,本文公式的理论结果无论是和Atabek等人的理论计算结果比较,还是和他们的实验结果比较,都是相当一致的.     

全附体ONRT船模在波浪中自航的数值模拟

采用基于重叠网格技术的CFD方法数值研究了全附体ONRT船模在迎浪工况中自航的水动力特性.文中数值计算采用自主开发的面向船舶与海洋工程的CFD求解器naoe-FOAM-SJTU.自航计算中船体运动及螺旋桨转动等通过重叠网格技术完成,波浪环境则采用求解器中的三维数值造波和消波模块实现.计算中自航船模的螺旋桨转速通过静水自航数值计算得出,波浪工况计算采用东京2015 CFD会议中标准算例进行设置.数值计算结果,如船体运动、实时航速变化等,与试验数据进行了对比分析.此外,给出了数值预报的推力和扭矩系数,并且通过详细的流场信息来分析和解释了船模在波浪中自航过程中的水动力变化情况.数值预报结果同试验值吻合较好,说明采用当前结合重叠网格技术和CFD的数值方法可以很好地预报波浪中自航问题.     

全附体ONRT船模在波浪中自航的数值模拟（英文）

采用基于重叠网格技术的CFD方法数值研究了全附体ONRT船模在迎浪工况中自航的水动力特性.文中数值计算采用自主开发的面向船舶与海洋工程的CFD求解器naoe-FOAM-SJTU.自航计算中船体运动及螺旋桨转动等通过重叠网格技术完成,波浪环境则采用求解器中的三维数值造波和消波模块实现.计算中自航船模的螺旋桨转速通过静水自航数值计算得出,波浪工况计算采用东京2015 CFD会议中标准算例进行设置.数值计算结果,如船体运动、实时航速变化等,与试验数据进行了对比分析.此外,给出了数值预报的推力和扭矩系数,并且通过详细的流场信息来分析和解释了船模在波浪中自航过程中的水动力变化情况.数值预报结果同试验值吻合较好,说明采用当前结合重叠网格技术和CFD的数值方法可以很好地预报波浪中自航问题.     

半无穷长圆管内的低雷诺数入口流

1969年Lew及Fung[1]计算了圆管内的低雷诸数入口流.1982年Dagan等人[2]得到了有限长圆柱形孔道内蠕动流的级数解.[1]中所得的数值解实质上代表有限长圆管内的低雷诺数入口流,因为一般解中的富氏积分已用富氏级数代替.本文直接计算富氏积分,更精确地求出了真正的半无穷长圆管内Stokes入口流的速度分布,压力分布以及流函数,与此对应的入口段长度为圆管半径的1.2倍,接近于Lew及Fung得到的结果1.3倍.此外,本文还研究了配置法的收敛性,证明了此法在入口流问题中具有很好的收敛性,因此可以在其他类似的问题中采用.     

圆管层流与湍流进口段效应修正系数的研究

本文对于光滑圆管层流与湍流进口段长度以及压力损失和流量进口段效应修正系数.提出统一的计算方法通过具体算例,提供了层流与湍流进口段效应压力损失和流量修正系数的理论计算公式.并且给出计及湍流进口段效应的流量计算方法.理论计算结果与实验数据比较表明,本文提供的公式.是简便而又可靠的.     

计算Riemann-Liouville分数阶积分和导数的数值算法

研究计算Riemann-Liouville (RL)分数阶积分和导数的数值算法.首先,分析了RL分数阶积分和导数的定义式,由于定义式中包含一个积分瑕点,使RL分数阶积分和导数难于计算.然后,给出了一种去掉积分瑕点的方法,在此基础上设计出计算RL分数阶积分和导数的数值算法,并证明了此数值算法具有一阶精度.最后,给出了计算实例,计算结果说明提出的算法是有效的.     

等离子体反应器流场的数值模拟

本文采用颗粒轨道模型对等离子体反应器进行数值模拟.数值格式用SIMPLER法.湍流模型用简单的次网格模型.计算结果给出速度场、温度场及颗粒轨道,对工程设计和工艺过程参数控制有参考价值.     

Numerical Simulations Of Three Dimensional Glottis Flows

A detailed numerical investigation of the transient three dimensional flow downstream of the human glottis is performed. For the present work a threefoldly scaled up realistic static model of the vocal folds and a transient inflow function, which is obtained from experimental data is applied to the CFD model. In the numerical simulations transient flow features like forming ring vortexes and axis switching are resolved which are also known for elliptical and rectangular jets. Similar structures can also be identified in the results of the time–resolved PIV measurements. The influence of the ventricular folds and the geometrical variation of the glottis contacts are also considered in the numerical study. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Effect of fiber orientation and cross section of composite tubes on their energy absorption ability in axial dynamic loading

The crushing behavior of composite tubes in axial impact loading is investigated. Tubes of circular and rectangular cross section are simulated using an LS-DYNA software. The effect of fiber orientation on the energy absorbed in laminated composite tubes is also studied. The results obtained show that rectangular tubes absorb less energy than circular ones, and their maximum crushing load is also lower. The composite tubes with a [+θ/ -θ] lay-up configuration absorb a minimum amount of energy at θ = 15°. The simulation results for a rectangular composite tube with a [+30/–30] lay-up configuration are compared with available experimental data. Cylindrical composite tubes fabricated from woven glass/polyester composites with different lay-ups were also tested using a drop-weight impact tester, and very good agreement between experimental and numerical results is achieved.     

Eigenfrequencies of a tube bundle immersed in a fluid

In this paper we study a simplified version of a mathematical model that describes the eigenfrequencies and eigenmotions of a coupled system consisting of a set of tubes (or a tube bundle) immersed in an incompressible perfect fluid. The fluid is assumed to be contained in a rectangular cavity, and the tubes are assumed to be identical, and periodically distributed in the cavity. The mathematical model that governs this physical problem is an elliptic differential eigenvalue problem consisting of the Laplace equation with a nonlocal boundary condition on the holes, and a homogeneous Neumann boundary condition on the boundary of the cavity. In the simplified model that we study in this paper, the Neumann condition is replaced by a periodic boundary condition. Our goal in studying this simple version is to derive some basic properties of the problem that could serve as a guide to envisage similar properties for the original model. In practical situations, this kind of problem needs to be solved for tube bundles containing a very large number of tubes. Then the numerical analysis of these problems is in practice very expensive. Several approaches to overcome this difficulty have been proposed in the last years using homogenization techniques. Alternatively, we propose in this paper an approach that consists in obtaining an explicit decomposition of the problem into a finite family of subproblems, which can be easily solved numerically. Our study is based on a generalized notion of periodic function, and on a decomposition theorem for periodic functions that we introduce in the paper. Our results rely on the theory of almost periodic functions, and they provide a simple numerical method for obtaining approximations of all the eigenvalues of the problem for any number of tubes in the cavity. We also discuss a numerical example.     

ON THE REASONING EXPRESSIONS AND THE NATURE OF DEDUCTION RULES

Reasoning expressions are those which express the reaasoning procedure by means ofonly deduction rules and the initial formulas(axioms or assumptions)without the helpof any intermediate results.They express the procedure systematically,completely andconcisely.The deduction rules are mappings from formulas(premises)to formula(conclusion).The elementary rules are certain propositional connectives(but notnecessarily truth functions)while the higher rules are certain quantifiers.Besides,thedetachment rule is an inverse of the connective implication,and is itself the kernel ofdeduction method;while another inverse of implication(i.e.the suggestion rule)is thekernel of induction method.     

A corrected smoothed particle hydrodynamics method for solving transient viscoelastic fluid flows

In this work, a corrected smoothed particle hydrodynamics (CSPH) method is proposed and extended to the numerical simulation of transient viscoelastic fluid flows due to that its approximation accuracy in solving the Navier–Stokes equations is higher than that of the smoothed particle hydrodynamics (SPH) method, especially near the boundary of the domain. The CSPH approach comes with the idea of combining the SPH approximation for the interior particles with the modified smoothed particle hydrodynamics (MSPH) method for the exterior particles, this is because that the later method has higher accuracy than the SPH method although it also needs more computational cost. In order to show the validity of CSPH method to simulate unsteady viscoelastic flows problems, the planar shear flow problems, including transient Poiseuille, Couette flow and transient combined Poiseuille and Couette flow for the Oldroyd-B fluid are solved and compared with the analytical and SPH results. Subsequently, the general viscoelastic fluid based on the eXtended Pom–Pom (XPP) model is numerically investigated and the viscoelastic free surface phenomena of impacting drop are simulated by the CSPH for its extended application and the purpose of illustrating the ability of the proposed method. The numerical results are presented and compared with available solutions, which shows a very good agreement. All the numerical results show the higher accuracy and better stability of the CSPH than the SPH, especially for larger Weissenberg numbers.     

一种基于一维元胞自动机的灾害后果推演模型

突发事件发生后,其灾害后果随时间不断演变。了解灾害后果的发展趋势,有助于快速应急响应,为此提出了一种灾害后果推演模型。首先从承灾体角度,考虑灾害后果的共性,描述了突发事件灾害后果的演化机理;然后提出基于一维元胞自动机的灾害后果推演模型,完成了对元胞、元胞空间、邻居和规则的设置。最后,通过两个算例说明了本文提出模型的有效性。该模型适用于多种突发事件,推演的结果能为应急决策提供支持。     

Continuation and stability analysis for Bloch waves of the Gross-Pitaevskii equation

We describe a two-parameter continuation algorithm for computing Bloch waves of Bose-Einstein condensates (BEC) in optical lattices which is governed by the Gross-Pitaevskii equation (GPE). The Fourier collocation method and fourth-order Adini’s elements with penalty are used to discretize the GPE. We propose two different approaches so that the two-parameter continuation algorithm can be modified to compute closed tubes at the four corners of the Bloch band. We also study linear stability analysis for the GPE. We show that all the discrete steady-state solutions are numerically neutrally stable. Numerical results show that the four edges of the Bloch waves are surrounded by closed loops if the coefficient of the cubic nonlinear term is greater than that of the periodic potential. Moreover, closed tubes at the four corners of the Bloch band are obtained. The numerical results display superfluidity of BEC.     

Optimization on the cold drawing process of 6063 aluminium tubes

To obtain the T832 material behaviour, 6063-T4 tubes need to be cold drawn before artificial aging to enhance their mechanical properties. Since a high area reduction is required, multiple drawing passes are generally necessary to avoid fracture of tubes during the drawing process. These several steps need a time consuming process which could be reduced by optimization. The optimization of the tools geometry has been done through a finite element model with the objective of reducing the maximum stress level to draw tubes in one pass. Based on the optimum design, a new tool was built and acceptable accordance was observed between experiments and numerical results.     

瞬态热传导的奇异边界法及其MATLAB实现

基于动力学问题时间依赖基本解的奇异边界法是一种无网格边界配点法.该方法引入源点强度因子的概念从而避免了基本解的源点奇异性,具有数学简单、编程容易、精度高等优点.将该方法用于瞬态热传导问题的数值模拟,运用MATLAB实现该问题的数值研究,并创建相应的MATLAB工具箱.针对二维和三维瞬态热传导问题,进行了基于反插值技术和经验公式的奇异边界法MATLAB算例实现.针对支撑圆坯低温瞬态温度场的模拟结果表明,瞬态热传导奇异边界法的MATLAB工具箱具有简单、方便、精确可靠的优点.研究成果有助于发展瞬态热传导的奇异边界法,并为瞬态热传导问题的数值分析和仿真提供了一种简单高效的模拟工具.     

On the Green’s function for the Helmholtz operator in an impedance circular cylindrical waveguide

This paper addresses the problem of finding a series representation for the Green’s function of the Helmholtz operator in an infinite circular cylindrical waveguide with impedance boundary condition. Resorting to the Fourier transform, complex analysis techniques and the limiting absorption principle (when the undamped case is analyzed), a detailed deduction of the Green’s function is performed, generalizing the results available in the literature for the case of a complex impedance parameter. Procedures to obtain numerical values of the Green’s function are also developed in this article.