我国航空航天结构强度事业的开拓者——黄玉珊

黄玉珊先生是我国著名的力学家、航空航天科学家, 新中国航空高等教育的奠基人. 他不满14岁考入大学, 不满23岁时师从国际著名力学大师铁摩辛柯获博士学位, 随即义无反顾地回到战火连绵的祖国, 受聘中央大学教授. 新中国成立后, 面对航空工业初建时期技术基础薄弱、人才储备不足、教育体系不完善的困境, 他呕心沥血, 鞠躬尽瘁, 建立新学科, 开创新局面, 毕生致力于祖国航空航天教育和科技事业. 他被认为是我国旧飞机疲劳定寿和延寿、损伤容限评定和新飞机损伤容限设计最早的创始人之一, 也是我国飞机自激励振动研究和航天结构环境强度事业的开拓者. 本文通过回顾黄玉珊先生在固体力学、航空航天结构强度领域的主要学术成果和卓越贡献, 缅怀其坚定不移的报国情怀、严谨求实的治学态度、勇于开拓的创新精神、高瞻远瞩的学术视野以及求真务实的教育思想.      

物理增强的流场深度学习建模与模拟方法

流体运动理论上可用Navier?Stokes方程描述, 但由于对流项带来的非线性, 仅在少数情况可求得方程解析解. 对于复杂工程流动问题, 数值模拟难以高效精准计算高雷诺数流场, 实验或现场测量难以获得流场丰富细节. 近年来, 人工智能技术快速发展, 深度学习等数据驱动技术可利用灵活网络结构, 借助高效优化算法, 获得对高维、非线性问题的强大逼近能力, 为研究流体力学计算方法带来新机遇. 有别于传统图像识别、自然语言处理等典型人工智能任务, 深度学习模型预测的流场需满足流体物理规律, 如Navier?Stokes方程、典型能谱等. 近期, 物理增强的流场深度学习建模与模拟方法快速发展, 正逐渐成为流体力学全新研究范式: 根据流体物理规律选取网络输入特征或设计网络架构的方法称为物理启发的深度学习方法, 直接将流体物理规律显式融入网络损失函数或网络架构的方法称为物理融合的深度学习方法. 研究内容涵盖流体力学降阶模型、流动控制方程求解领域.      

浅谈流体力学发展——从Batchelor百年诞辰纪念大会谈起

G.K. Batchelor是20世纪国际流体力学大师,在均匀湍流理论和低雷诺数微流体力学领域做出了开创性的贡献,他对流动追求物理和定量性理解的思想影响了近百年流体力学的发展,是流体力学顶级期刊Journal of Fluid Mechanics的创刊人,也是剑桥大学应用数学与理论物理系(DAMTP)的创建者,培养并影响了一大批在流体动力学理论、实验流体力学、湍流及稳定性、环境流体力学、多相流体力学、磁流体力学、微纳米尺度流体动力学等诸多领域建树卓越的学者.本文以G.K. Batchelor诞辰100周年纪念活动为契机,简要回顾了流体力学近300年的发展历程.概述了流体力学发展历经的以数学和物理为基础建立理论框架的经典流体力学、以应用需求为导向促使自身跨越发展的近代流体力学和以学科融合为特点外延丰富的现代流体力学三个重要阶段;以师承关系、代表性学者及其主要学术贡献为线索,总结了现代流体力学四大学派的形成及其近百年的传承沿革;以历史和发展的视角浅谈当代流体力学发展的动力和趋势,并以风沙环境力学为例,简述流体力学为分支学科发展提供的支撑和引领作用,分支学科的需求为流体力学内生发展提供驱动力...     

On the Navier–Stokes equations simulation of the head-on collision between two surface solitary waves

The scope of this Note is to show the results obtained for simulating the two-dimensional head-on collision of two solitary waves by solving the Navier–Stokes equations in air and water. The work is dedicated to the numerical investigation of the hydrodynamics associated to this highly nonlinear flow configuration, the first numerical results being analyzed. The original numerical model is proved to be efficient and accurate in predicting the main features described in experiments found in the literature. This Note also outlines the interest of this configuration to be considered as a test-case for numerical models dedicated to computational fluid mechanics. To cite this article: P. Lubin et al., C. R. Mecanique 333 (2005).     

In memoriam Professor Ko Saito (1941-2004)

Professor Ko Saito of Hiroshima University suffering from hepatic cancer passed away on May 20, 2004 at the age of 63, one year before his planned retirement in March 2005. He graduated from Hiroshima University in 1964, supervised by Professor Yamamura, and started his career in high-temperature chemical reaction and was appointed to be a Research Associate of the Department of Chemistry at Hiroshima University in 1967. He received a Doctor of Science Degree in chemistry in 1973 and spent his study leave from May 1980 to March 1981 at the University of Essex UK and worked with Professor JN Bradley. He was promoted to Associate Professor in 1986 and then to Professor of the Hiroshima University in 1990.He had one hundred paper publications in well-known scientific journals, the Proceedings of International Symposia on Shock Waves (ISSW) and was an author of a Japanese text book on physical chemistry. He proved himself as an experimentalist as well as a theoretician and often told students that progress in chemical reaction studies should be maintained by conducting both experiments and analyses based on quantum chemistry. He taught a physical chemistry course to all first-year students introducing rigorously his approach.He devoted himself to identify and characterize features of molecules suddenly exposed to shock created high temperatures and particularly concentrated in thermal decomposition processes of carboxylic acid groups. These acid groups are found to decompose molecular and quantum chemical analyses supports the observation. This is one of the highlights of his achievements.He participated regularly in International and Domestic Shock Wave Symposia and intensively presented important results he and his group discovered. His latest works are presented in the Proceedings of ISSW24 held in Beijing, July 2004. Professor Saito was not only an educator but also an active chemical shock wave researcher serving continuously for 40 years and immensely contributed to establish the academic reputation of the Department of Chemistry of Hiroshima University. Undergraduate and graduate students and post doctoral fellows he once supervised are today remembering his philosophy in their mind and try to widen the world of shock wave chemistry.Published online: 18 February 2005[/PUBLISHED]     

On the pressure and stress singularities induced by steady flows of incompressible viscous fluids

Design for structural integrity requires an appreciation of where stress singularities can occur in structural configurations. While there is a rich literature devoted to the identification of such singular behavior in solid mechanics, to date there has been relatively little explicit identification of stress singularities caused by fluid flows. In this study, stress and pressure singularities induced by steady flows of viscous incompressible fluids are asymptotically identified. This is done by taking advantage of an earlier result that the Navier-Stokes equations are locally governed by Stokes flow in angular corners. Findings for power singularities are confirmed by developing and using an analogy with solid mechanics. This analogy also facilitates the identification of flow-induced log singularities. Both types of singularity are further confirmed for two global configurations by applying convergence-divergence checks to numerical results. Even though these flow-induced stress singularities are analogous to singularities in solid mechanics, they nonetheless render a number of structural configurations singular that were not previously appreciated as such from identifications within solid mechanics alone.     

A film flow model for analysing gravity-driven,thin wavy fluid films

To analyse the physics underlying gravity-driven runoff of thin wavy films, a film flow model is developed, and is solved with computational fluid dynamics. This model is based on the lubrication theory, and takes into account the gravitational, wall shear and surface tension forces. A key characteristic of the model is that it assumes only one computational cell over the film height, which enables studying film flow on larger computational domains. A main aim of this study is to perform a detailed validation of the numerical model. The film flow model is validated against several experiments of gravity-driven, thin fluid films on smooth surfaces. The time-averaged film thickness and the fluid speed profiles predicted by the model show very good agreement with experimental results. Similarly, the film flow model is able to predict the wave speeds with sufficient accuracy. The energy spectra of the waves, where higher frequency waves are present in film flows at higher Reynolds numbers, show an exponentially decaying trend at these high frequencies. The model performs better than the Nusselt equation for film flows, which under-predicts the time-averaged film thickness and over-predicts the time-averaged fluid speeds, even for flows at low Reynolds numbers. The film flow model is compared qualitatively for fingering behaviour. This model also allows to investigate film flows on large surfaces, which can be rough, curved and of complex geometrical shape.     

高雷诺数流动的基本控制方程体系和扩散跑物化Navier-Stokes方程组的意义和用途

在计算机发达的时代, 高雷诺($Re$)数绕流计算中有无必要使用简化NS方程组, 本文讨论这个问题. 主要内容如下: (1)高$Re$数绕流包含3种基本流动: 所有方向对流占优流动、所有方向对流扩散竞争流动和部分方向对流占优部分方向对流扩散竞争流动(简称干扰剪切流动), 3个基本流动的特征彼此不同且在流场中所占领域大小彼此相差悬殊, NS方程区域很小,它们的最简单控制方程组Euler、Navier-Stokes (NS)和扩散抛物化(DP) NS方程组的数学性质彼此不同, 因此利用Euler-DPNS-NS方程组体系分析计算高$Re$数绕流流动就是一个合乎逻辑的选择, 该法与利用单一NS方程组的常用方法可以彼此检验和补充. (2)流体之间以及流体与外界的动量、能量和质量交换, 流态从层流到湍流的演化主要发生在干扰剪切流动中, 干扰剪切流及其最简单控制方程------DPNS方程组具有基础意义; DPNS方程组笔者在1967年已提出. (3)诸简化NS方程组: DPNS、抛物化(P)NS、薄层(TL)NS、黏性层(VL)NS方程组的发展、相互关系, 它们的历史贡献和今后的用途; 它们的数学性质均为扩散抛物型, 但它们包含的黏性项彼此有所不同; 从流体力学角度来看, 它们中只有DPNS方程组能够准确描述干扰剪切流动. 提出把诸简化NS方程组统一为DPNS方程组的建议. (4)干扰剪切流------DPNS方程组与无干扰剪切流------边界层方程组之间的关系以及进一步研究干扰剪切流的意义.      

流体力学问题的三次样条配置法

本文给出了三次样条配置法在流体力学问题数值解中的应用以及在这一领域的新进展。给出了流体力学方程中主要的样条函数关系和解算步骤。所有情形都是便于反演的三对角形矩阵。简要评述了SADI方法和样条方法在每一坐标方向的分步计算方法、截断误差和稳定性。给出了处理混合边界条件的一般公式。最后简要讨论了样条近似引起的数值弥散和耗散。      

Numerical solutions of incompressible Euler and Navier-Stokes equations by efficient discrete singular convolution method

An efficient discrete singular convolution (DSC) method is introduced to the numerical solutions of incompressible Euler and Navier-Stokes equations with periodic boundary conditions. Two numerical tests of two-dimensional Navier-Stokes equations with periodic boundary conditions and Euler equations for doubly periodic shear layer flows are carried out by using the DSC method for spatial derivatives and fourth-order Runge-Kutta method for time advancement, respectively. The computational results show that the DSC method is efficient and robust for solving the problems of incompressible flows, and has the potential of being extended to numerically solve much broader problems in fluid dynamics. The project supported by the National Natural Science Foundation of China (No.19902010).     

UNSTEADY CONFINED VISCOUS FLOWS WITH OSCILLATING WALLS AND MULTIPLE SEPARATION REGIONS OVER A DOWNSTREAM-FACING STEP

This paper presents computational solutions for unsteady viscous flows in channels with a downstream-facing step, followed by an oscillating floor. These solutions of the unsteady Navier–Stokes equations are obtained with a time-integration method using artificial compressibility in a fixed computational domain, which is obtained via a time-dependent coordinate transformation from the fluid domain with moving boundaries. The computational method is first validated for steady viscous flows past a downstream-facing step by comparison with previous numerical solutions and experimental results. This method is then used to obtain solutions for unsteady viscous flows with multiple separation regions over a downstream-facing step with oscillating walls, for which there are no previously known solutions. Thus, the present results may be used as benchmark solutions for the unsteady viscous flows with multiple separation regions between fixed and oscillating walls.     

Computational intellectual analytical theory of computational analytical approach to rotating flow of non-Newtonian fluid

Ａｎｅｗｆｉｅｌｄｏｆｍｏｄｅｒｎｃｏｍｐｕｔｅｒｓｃｉｅｎｃｅ—ａｒｔｉｆｉｃｉａｌｉｎｔｅｌｌｉｇｅｎｃｅｉｓｄｅｖｅｌｏｐｅｄｒａｐｉｄｌｙ．Ｔｈｅｓｙｍｂｏｌｉｃｍａｎｉｐｕｌａｔｉｏｎｉｓｔｈｅｆｒｏｎｔｄｉｒｅｃｔｉｏｎｉｎｔｈｅａｒｔｉｆｉｃｉａｌｉｎｔｅｌｌｉｇｅｎｃｅ．Ａｓｅｒｉｅｓｏｆｔｈｅｃｏｍｐｕｔｅｒｓｏｆｔｗａｒｅｉｓｄｅｖｅｌｏｐｅｄｆｏｒｔｈｅｃｏｍｐｕｔａｔｉｏｎａｌｍａｎｉｐｕｌａｔｉｏｎ,ｓｕｃｈａｓＭａｃｓｙｍａ,Ｍａｐｌｅ,Ｍａｔｈｅｍａｔｉｃａ,ｔｏｃ…     

Unified one‐fluid formulation for incompressible flexible solids and multiphase flows: Application to hydrodynamics using the immersed structural potential method (ISPM)

In this paper, we present a two‐dimensional computational framework for the simulation of fluid‐structure interaction problems involving incompressible flexible solids and multiphase flows, further extending the application range of classical immersed computational approaches to the context of hydrodynamics. The proposed method aims to overcome shortcomings such as the restriction of having to deal with similar density ratios among different phases or the restriction to solve single‐phase flows. First, a variation of classical immersed techniques, pioneered with the immersed boundary method (IBM), is presented by rearranging the governing equations, which define the behaviour of the multiple physics involved. The formulation is compatible with the “one‐fluid” formulation for two‐phase flows and can deal with large density ratios with the help of an anisotropic Poisson solver. Second, immersed deformable structures and fluid phases are modelled in an identical manner except for the computation of the deviatoric stresses. The numerical technique followed in this paper builds upon the immersed structural potential method developed by the authors, by adding a level set–based method for the capturing of the fluid‐fluid interfaces and an interface Lagrangian‐based meshless technique for the tracking of the fluid‐structure interface. The spatial discretisation is based on the standard marker‐and‐cell method used in conjunction with a fractional step approach for the pressure/velocity decoupling, a second‐order time integrator, and a fixed‐point iterative scheme. The paper presents a wide d range of two‐dimensional applications involving multiphase flows interacting with immersed deformable solids, including benchmarking against both experimental and alternative numerical schemes.     

Incorporation of polymer diffusivity and migration into constitutive equations

Given a general one-particle constitutive equation for the stress tensor, we discuss how to incorporate the additional effects of polymer diffusivity and migration into that constitutive equation within the framework of continuum mechanics. For the example of an upper-convected Maxwell model representing the polymer contribution to the stress tensor of a dilute polymer solution, we describe i) how to modify the constitutive equation for the stress tensor to include diffusion and migration effects, ii) how to formulate a balance equation for the polymer mass density in order to describe the nonhomogeneous composition of the polymer solution resulting from migration, and iii) how to close the extended set of coupled equations by means of further constitutive equations for the migration velocity and the diffusion tensor. In order to guarantee the material objectivity for all equations, we formulate them in the body tensor formulation of continuum mechanics (and then translate them into Cartesian space). The proposed equations are compared to results of a recent kinetic-theory approach.Dedicated to Professor Arthur S. Lodge on the occasion of his 70th birthday and his retirement from the University of Wisconsin.     

Multiple-relaxation-time lattice Boltzmann method for study of two-lid-driven cavity flow solution multiplicity

As a fundamental subject in fluid mechanics, sophisticated cavity flow patterns due to the movement of multi-lids have been routinely analyzed by the computational fluid dynamics community. Unlike those reported computational studies that were conducted using more conventional numerical methods, this paper features employing the multiple-relaxation-time (MRT) lattice Boltzmann method (LBM) to numerically investigate the two-dimensional cavity flows generated by the movements of two adjacent lids. The obtained MRT-LBM results reveal a number of important bifurcation flow features, such as the symmetry and steadiness of cavity flows at low Reynolds numbers, the multiplicity of stable cavity flow patterns when the Reynolds number exceeds its first critical value, as well as the periodicity of the cavity flow after the second critical Reynolds number is reached. Detailed flow characteristics are reported that include the critical Reynolds numbers, the locations of the vortex centers, and the values of stream function at the vortex centers. Through systematic comparison against the simulation results obtained elsewhere by using the lattice Bhatnagar–Gross–Krook model and other numerical schemes, not only does the MRT-LBM approach exhibit fairly satisfactory accuracy, but also demonstrates its remarkable flexibility that renders the adjustment of its multiple relaxation factors fully manageable and, thus, particularly accommodates the need of effectively investigating the multiplicity of flow patterns with complex behaviors.     

数值摄动算法及其CFD格式

作者提出的数值摄动算法把流体动力学效应耦合进NS方程组和对流扩散（CD）方程离散的数学基本格式（MBS）,特别是耦合进最简单的MBS即一阶迎风和二阶中心格式之中,由此构建成一系列新格式,称呼方便和强调耦合流体动力学起见,称它们为流体力学基本格式（FMBS）。构建FMBS的主要步骤是把MBS中的通量摄动重构为步长的幂级数,利用空间分裂和导出的高阶流体动力学关系式,把结点变量展开成Taylor级数,通过消除重构格式修正微分方程的截断误差诸项求出幂级数的待定系数,由此获得非线性FMBS。FMBS的公式是MBS与 （及 ）之简单多项式的乘积, 和 分别是网格Reynolds数和网格CFL数。FMBS和MBS使用相同结点,简单性彼此相当,但FMBS精度高稳定范围大,例如FMBS包含了许多绝对稳定和绝对正型、高阶迎风和中心有限差分（FD）格式和有限体积（FV）格式,这些格式对网格Reynolds数的任意值均为不振荡格式。可见对不振荡CFD格式的构建,数值摄动算法提供了不同于调节数值耗散等常见的人为构建方法,而利用流体力学自身关系以及把迎风机制通过上、下游摄动重构引入中心MBS的解析构建方法,FMBS除了直接应用于流体计算外;对于通过调节数值耗散、色散和数值群速度特性重构高分辨率格式的研究,最简单FMBS提供了比最简单MBS更精确、但同样简单的基础和起步格式。FMBS用于计算不可压缩流,可压缩流,液滴萃取传质,微通道两相流等,均获得良好数值结果或与已有Benchmark解一致的数值结果。已有文献称数值摄动算法为新型高精度格式和高的算法和高的格式;本文FMBS比数值摄动格式的称呼可更好反映FMBS的物理内容。文中也讨论了值得进一步研究的一些课题,该法亦可用于其它一些数学物理方程（例如,简化Boltzmann方程、磁流体方程、KdV-Burgers方程等）MBS耦合物理动力学效应的重构。      

Computational modeling of biomagnetic micropolar blood flow and heat transfer in a two-dimensional non-Darcian porous medium

We study theoretically and computationally the incompressible, non-conducting, micropolar, biomagnetic (blood) flow and heat transfer through a two-dimensional square porous medium in an (x,y) coordinate system, bound by impermeable walls. The magnetic field acting on the fluid is generated by an electrical current flowing normal to the x–y plane, at a distance l beneath the base side of the square. The flow regime is affected by the magnetization B ₀ and a linear relation is used to define the relationship between magnetization and magnetic field intensity. The steady governing equations for x-direction translational (linear) momentum, y-direction translational (linear) momentum, angular momentum (micro-rotation) and energy (heat) conservation are presented. The energy equation incorporates a special term designating the thermal power per unit volume due to the magnetocaloric effect. The governing equations are non-dimensionalized into a dimensionless (ξ,η) coordinate system using a set of similarity transformations. The resulting two point boundary value problem is shown to be represented by five dependent non-dimensional variables, f _ξ  (velocity), f _η (velocity), g (micro-rotation), E (magnetic field intensity) and θ (temperature) with appropriate boundary conditions at the walls. The thermophysical parameters controlling the flow are the micropolar parameter (R), biomagnetic parameter (N _H ), Darcy number (Da), Forchheimer (Fs), magnetic field strength parameter (Mn), Eckert number (Ec) and Prandtl number (Pr). Numerical solutions are obtained using the finite element method and also the finite difference method for Ec=2.476×10⁻⁶ and Prandtl number Pr=20, which represent realistic biomagnetic hemodynamic and heat transfer scenarios. Temperatures are shown to be considerably increased with Mn values but depressed by a rise in biomagnetic parameter (N _H ) and also a rise in micropolarity (R). Translational velocity components are found to decrease substantially with micropolarity (R), a trend consistent with Newtonian blood flows. Micro-rotation values are shown to increase considerably with a rise in R values but are reduced with a rise in biomagnetic parameter (N _H ). Both translational velocities are boosted with a rise in Darcy number as is micro-rotation. Forchheimer number is also shown to decrease translational velocities but increase micro-rotation. Excellent agreement is demonstrated between both numerical solutions. The mathematical model finds applications in blood flow control devices, hemodynamics in porous biomaterials and also biomagnetic flows in highly perfused skeletal tissue. Dedicated to Professor Y.C. Fung (1919-), Emeritus Professor of Biomechanics, Bioengineering Department, University of California at San Diego, USA for his seminal contributions to biomechanics and physiological fluid mechanics over four decades and his excellent encouragement to the authors, in particular OAB, with computational biofluid dynamics research.     

Towards a constitutive law for the unsteady contact stress in granular media

In order to build a unified modelling for granular media by means of Eulerian averaged equations, it is necessary to study two contributions in the effective Cauchy stress tensor: the first one concerns solid and fluid matter, including contact and collisions between grains; the second one focuses on the random movements of grains and fluid, similar to Reynolds stress for turbulent flows. It is shown that the point of view of piecewise continuous media already used for two phase flows allows one to derive a constitutive equation for the first contribution, under the form of a partial differential equation. Similarly as for the Reynolds stress in turbulent flows, this equation cannot be written only in terms of averaged quantities without adequate approximations. The structure of the closed equation is discussed with respect to irreversible thermodynamics, and in connection with some already existing models. It is emphasised that numerical simulations by the discrete elements method can be used in order to validate these approximations, through numerical experiments statistically considered. Finally an extension of this approach to the second contribution of the effective Cauchy stress tensor is briefly discussed, showing how the modelling of both contributions have to be coupled.      

毕生献中华──纪念林同骥先生诞辰80周年

林同骥先生是我国著名流体力学家．在稀薄气体力学、高超声速、跨声速空气动力学和不可压缩流体动力学等广泛的领域中都有重要的研究成果,尤其是在再入飞行器防热、再人物体的表面烧蚀和不可压缩振荡流等方面,为我国航天事业和海洋工程的发展作出了贡献．林同骥先生在１９７９～１９９０年期间先后任《力学学报》副主编、主编．为《力学学报》在改革开放新形势下的发展和进人国际交流大网络倾注了大量心血．在纪念林同骥先生诞辰８０周年之际,我们向读者介绍这位毕生献中华的科学家,并刊出他的学生们的一些论文,以表达怀念与崇敬之情．     

Practical comparison of differential viscoelastic constitutive equations in finite element analysis of planar 4:1 contraction flow

Finite element modeling of planar 4:1 contraction flow (isothermal incompressible and creeping) around a sharp entrance corner is performed for favored differential constitutive equations such as the Maxwell, Leonov, Giesekus, FENE-P, Larson, White-Metzner models and the Phan Thien-Tanner model of exponential and linear types. We have implemented the discrete elastic viscous stress splitting and streamline upwinding algorithms in the basic computational scheme in order to augment stability at high flow rate. For each constitutive model, we have obtained the upper limit of the Deborah number under which numerical convergence is guaranteed. All the computational results are analyzed according to consequences of mathematical analyses for constitutive equations from the viewpoint of stability. It is verified that in general the constitutive equations proven globally stable yield convergent numerical solutions for higher Deborah number flows. Therefore one can get solutions for relatively high Deborah number flows when the Leonov, the Phan Thien-Tanner, or the Giesekus constitutive equation is employed as the viscoelastic field equation. The close relationship of numerical convergence with mathematical stability of the model equations is also clearly demonstrated.