Singular perturbation for the weakly nonlinear reaction diffusion equation with boundary perturbation

In this paper, a class of nonlinear singularly perturbed initial boundary value problems for reaction diffusion equations with boundary perturbation are considered under suitable conditions. Firstly, by dint of the regular perturbation method, the outer solution of the original problem is obtained. Secondly, by using the stretched variable and the expansion theory of power series the initial layer of the solution is constructed. And then, by using the theory of differential inequalities, the asymptotic behavior of the solution for the initial boundary value problems is studied. Finally, using some relational inequalities the existence and uniqueness of solution for the original problem and the uniformly valid asymptotic estimation are discussed.     

Asymptotic analysis of perforated plates and membranes. Part 1: Static problems for small holes

Static problems for the elastic plates and rods periodically perforated by small holes of different shapes are solved using the asymptotic approach based on the combination of the asymptotic technique and the multi-scale homogenization method. Using the asymptotic homogenization method the original boundary-value problem is reduced to the combination of two types of problems. First one is a recurrent system of unit cell problems with the conditions of periodic continuation. And the second problem is a homogenized boundary-value problem for the entire domain, characterized by the constant effective coefficients obtained from the solution of the unit cell problems. The combination of the perturbation method and the technique of successive approximations is applied for the solution of the unit cell problems. Taking into the account small size of holes the method of perturbation of the shape of the boundary and the Schwarz alternating method are used. The problems of torsion of a rod with perforated cross-section; deflection of the perforated membrane loaded by a normal load; and bending of perforated plates with circular and square holes are considered consecutively. The error of the applied asymptotic techniques is estimated and the high accuracy of the obtained solutions is demonstrated.     

A CLASS OF NONLINEAR NONLOCAL SINGULARLY PERTURBED PROBLEMS FOR REACTION DIFFUSION EQUATIONS WITH BOUNDARY PERTURBATION

A class of nonlinear nonlocal for singularly perturbed Robin initial boundary value problems for reaction diffusion equations with boundary perturbation is considered. Under suitable conditions, first, the outer solution of the original problem was obtained. Secondly, using the stretched variable, the composing expansion method and the expanding theory of power series the initial layer was constructed. Finally, using the theory of differential inequalities the asymptotic behavior of solution for the initial boundary value problems was studied, and educing some relational inequalities the existence and uniqueness of solution for the original problem and the uniformly valid asymptotic estimation were discussed.     

磁流体动力学在航空工程中的应用与展望

介绍了磁流体动力学在航空工程中的主要应用方式,主要包括:磁流体冲压组合发动机、磁流体涡轮组合发动机、燃烧室后磁流体发电、表面磁流体发电、磁流体加速风洞、磁流体推力矢量、进气道大尺寸磁流体流动控制、边界层分离流动控制、边界层转捩控制、飞行器头部热流控制等;探讨了磁流体技术在应用中存在的关键科学与技术问题,对导电流体的产生、磁流体实验设备与实验技术、多场耦合机理及数值模拟方法等进行了分析;最后对磁流体技术在航空工程上的应用与发展进行了总结与展望.     

奇异摄动在层流预混火焰理论研究中的应用

奇异摄动被广泛应用于求取力学问题的近似解.一个典型问题就是流体力学中的边界层.郭永怀先生曾发展了适用于平板黏性流动边界层问题的奇异摄动理论.类似于流体力学中的边界层,燃烧研究中的层流预混火焰也可以通过奇异摄动理论进行分析,在燃烧研究中通常称其为大活化能渐近分析.本文介绍了大活化能渐近分析在一维平面预混火焰和球形传播火焰中的应用及相关研究进展.首先介绍了预混火焰结构及其涉及的不同特征尺度,分析了大活化能条件下出现的特征尺度分离,并给出了关于平面预混火焰大活化能渐近分析的详细推导,讨论了热辐射对火焰传播的影响;然后介绍了大活化能渐近分析在点火与球形传播火焰方面的应用,指出了只有能够同时描述点火与球形火焰传播的理论才能准确地预测临界点火条件,并讨论了考虑链式反应的点火与火焰传播理论,分析了热辐射对球形火焰传播的影响,给了关于火焰稳定性理论研究的发展趋势.最后,基于当前研究进展对未来的研究方向进行了展望,其中涉及多步化学反应、低温冷火焰、复杂流动、辐射重吸收等.      

Elastic Green’s Function in Anisotropic Bimaterials Considering Interfacial Elasticity

The two-dimensional elastic Green’s function is calculated for a general anisotropic elastic bimaterial containing a line dislocation and a concentrated force while accounting for the interfacial structure by means of a generalized interfacial elasticity paradigm. The introduction of the interface elasticity model gives rise to boundary conditions that are effectively equivalent to those of a weakly bounded interface. The equations of elastic equilibrium are solved by complex variable techniques and the method of analytical continuation. The solution is decomposed into the sum of the Green’s function corresponding to the perfectly bonded interface and a perturbation term corresponding to the complex coupling nature between the interface structure and a line dislocation/concentrated force. Such construct can be implemented into the boundary integral equations and the boundary element method for analysis of nano-layered structures and epitaxial systems where the interface structure plays an important role.     

Flow of viscoelastic fluids through a porous channel—I

The flow of viscoelastic fluids through a porous channel with one impermeable wall is computed. The flow is characterized by a boundary value problem in which the order of the differential equation exceeds the number of boundary conditions. Three solutions are developed: (i) an exact numerical solution, (ii) a perturbation solution for small R, the cross-flow Reynold's number and (iii) an asymptotic solution for large R. The results from exact numerical integration reveal that the solutions for a non-Newtonian fluid are possible only up to a critical value of the viscoelastic fluid parameter, which decreases with an increase in R. It is further demonstrated that the perturbation solution gives acceptable results only if the viscoelastic fluid parameter is also small. Two more related problems are considered: fluid dynamics of a long porous slider, and injection of fluid through one side of a long vertical porous channel. For both the problems, exact numerical and other solutions are derived and appropriate conclusions drawn.     

移动接触下求解二维闭合裂纹的双重迭代边界元法

使用子域边界元法对受移动接触弹性体作用下的二维闭合裂纹问题进行了数值计算。由于两弹性体的接触界面和裂纹表面的接触范围的大小和接触状态事先是未知的 ,对此 ,在两个接触表面同时采用迭代的方法进行了求解。在裂纹的每个裂尖上都采用了四分之一的奇异单元以保证裂尖位移场和应力场奇异性的满足。用我们编制的二维裂纹问题程序对一些中心裂纹问题进行了计算 ,计算结果与经典断裂力学的理论值比较吻合。在无摩擦的条件下 ,对一些具有不同角度且受移动接触弹性体作用下的闭合裂纹问题进行了数值计算 ,得到了一些耦合作用下的应力强度因子的计算结果     

ANTI-PLANE ANALYSIS FOR ELLIPTICAL INCLUSION IN MAGNETOELECTROELASTIC MATERIALS

This paper considers the multi-field coupling in magneroelectroelastic composite materials consisting of the inclusion and the matrix are magnetoelectroelastic materials. The mechanical,electric and magnetic fields around an elliptical cylinder inclusion are formulated by complex potentials. Inside the inclusion,the strain,electric and magnetic fields are found to be uniform and vary with the shape of the ellipse. When the inclusion is reduced to a crack,along the interface,the strain,electric field strength and magnetic field strength equal the corresponding remote ones,which can be used as the boundary condition. Special cases,such as a rigid and permeable inclusion,a soft and impermeable inclusion,a line inclusion and a crack problem are discussed in detail.     

输气管道壁面涂料减阻机理的实验研究

用IFA-300热线风速仪以高于对应最小湍流时间尺度的分辨率精细测量了风洞中不同壁面涂料的管道湍流边界层不同法向位置流向速度分量的时间序列信号,利用湍流边界层近壁区域对数律平均速度剖面与壁面摩擦速度、流体黏性系数等内尺度物理量的关系和壁面摩擦速度与壁面摩擦切应力的关系,在准确测量湍流边界层近壁区域对数律平均速度剖面的基础上,间接测量湍流边界层的壁面摩擦阻力．对不同壁面涂料的壁湍流脉动速度信号用子波分析进行多尺度分解,用子波系数的瞬时强度因子和平坦因子检测管道湍流边界层中的多尺度相干结构,提取不同尺度相干结构的条件相位平均波形,对比研究输气管道壁面涂料的减阻机理．     

A multi-field incremental variational framework for gradient-extended standard dissipative solids

The paper presents a constitutive framework for solids with dissipative micro-structures based on compact variational statements. It develops incremental minimization and saddle point principles for a class of gradient-type dissipative materials which incorporate micro-structural fields (micro-displacements, order parameters, or generalized internal variables), whose gradients enter the energy storage and dissipation functions. In contrast to classical local continuum approaches to inelastic solids based on locally evolving internal variables, these global micro-structural fields are governed by additional balance equations including micro-structural boundary conditions. They describe changes of the substructure of the material which evolve relatively to the material as a whole. Typical examples are theories of phase field evolution, gradient damage, or strain gradient plasticity. Such models incorporate non-local effects based on length scales, which reflect properties of the material micro-structure. We outline a unified framework for the broad class of first-order gradient-type standard dissipative solids. Particular emphasis is put on alternative multi-field representations, where both the microstructural variable itself as well as its dual driving force are present. These three-field settings are suitable for models with threshold- or yield-functions formulated in the space of the driving forces. It is shown that the coupled macro- and micro-balances follow in a natural way as the Euler equations of minimization and saddle point principles, which are based on properly defined incremental potentials. These multi-field potential functionals are outlined in both a continuous rate formulation and a time-space-discrete incremental setting. The inherent symmetry of the proposed multi-field formulations is an attractive feature with regard to their numerical implementation. The unified character of the framework is demonstrated by a spectrum of model problems, which covers phase field models and formulations of gradient damage and plasticity.     

考虑转动梯度的界面问题研究

在考虑转动梯度效应的基础上,对界面附近的应力状态进行了研究,首先应用摄动法求解了基于偶应力理论的准轴对称问题,并在此基础上分析了界面问题和边界约束效应,结果表明,在剪应力作用下,在界面附近和固定边界附近存在一组边界效应解,相对于经典的弹性理论结果,它对剪应力的可观的修正。     

Effect of Internal Heat Generation on the Onset of Marangoni Convection in a Fluid Layer Overlying a Layer of an Anisotropic Porous Medium

Linear stability analysis has been performed to investigate the effect of internal heat generation on the criterion for the onset of Marangoni convection in a two-layer system comprising an incompressible fluid-saturated anisotropic porous layer over which lies a layer of the same fluid. The upper non-deformable free surface and the lower rigid surface are assumed to be insulated to temperature perturbations. The fluid flow in the porous layer is governed by the modified Darcy equation and the Beavers–Joseph empirical slip condition is employed at the interface between the two layers. The resulting eigenvalue problem is solved exactly. Besides, analytical expression for the critical Marangoni number is also obtained by using regular perturbation technique with wave number as a perturbation parameter. The effect of internal heating in the porous layer alone exhibits more stabilizing effect on the system compared to its presence in both fluid and porous layers and the system is least stable if the internal heating is in fluid layer alone. It is found that an increase in the value of mechanical anisotropy parameter is to hasten the onset of Marangoni convection while an opposite trend is noticed with increasing thermal anisotropy parameter. Besides, the possibilities of controlling (suppress or augment) Marangoni convection is discussed in detail.     

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

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

Laminar flow through a channel with uniformly porous walls of different permeability

Summary The steady laminar flow of a viscous incompressible fluid through a two-dimensional channel, having fluid sucked or injected with different velocities through its uniformly porous parallel walls is considered. A solution for small suction Reynolds number has been given by the authors in a previous paper. The purpose of this paper is to present a solution valid for large Reynolds numbers for the cases of (i) suction at both walls, and (ii) suction at one wall and injection at the other. A technique of matching outer and inner expansions is used to obtain an asymptotic solution for both of these cases. Further a perturbation solution for the case of suction at one wall and injection at the other is obtained by choosing the difference between two wall velocities as the perturbation parameter. Both asymptotic and perturbation solutions are confirmed by exact numerical solutions. As expected, the resulting solutions show the presence of the usual suction boundary layers in both types of flow considered in this paper.     

Stress-strain state in the vicinity of a crack tip under mixed loading

A method is proposed to calculate the eigenvalues of the class of nonlinear eigenvalue problems resulting from the problem of determining the stress-strain state in the vicinity of a crack tip in power-law materials over the entire range of mixed modes of deformation, from the opening mode to pure shear. The proposed approach was used to found eigenvalues of the problem that differ from the well-known eigenvalue corresponding to the Hutchinson-Rice-Rosengren solution. The resulting asymptotic form of the stress field is a self-similar intermediate asymptotic solution of the problem of a crack in a damaged medium under mixed loading. Using the new asymptotic form of the stress field and introducing a self-similar variable, we obtained an asymptotic solution of the problem of a crack in a damaged medium and constructed the regions of dispersed material near the crack.     

On the solution for steady state and time dependent flows in certain channels with small wall curvature

An asymptotic scheme is presented for the solution of the steady state and time dependent stream functions for flows in symmetric curved walled channels. In this scheme a class of non-linear Jeffery-Hamel solutions appear at O(1), and thus provide the first approximation to the steady state stream function. This class of Jeffery-Hamel solutions are evaluated by using a simple perturbation about Poiseuille flow. The classic Orr-Sommerfeld eigenproblem appears at O(1) in the asymptotic development of the time dependent stream function, but here there is a slow streamwise dependence. This eigenvalue problem, for a complex wave number, is solved using an algorithm which automatically provides an initial guess which is then used to iterate to the correct eigenvalue. Higher order terms in the asymptotic development, for both the steady state and time dependent stream functions, are evaluated to provide a solution for the total stream function.     

数智流体力学的发展及油气渗流领域应用

大数据及人工智能技术的崛起推动了数智流体力学的快速发展.数智流体力学是将流体力学、大数据和人工智能相结合,以流体力学场景需求为导向,形成以“数”为基础,以“智”为核心,以算力为支撑的新研究范式.核心内涵是要以数据驱动为主,融合物理信息、专家经验等先验知识,利用智能化手段构建“数据+物理”双驱动的数智模型,解决场景需求问题.数智流体力学在建模灵活性、运算效率、计算精度方面具有十分明显的优势,其应用潜力已经在多尺度流动、多场耦合以及流场建模等方面得到验证.数智流体力学研究范式包括数据治理和智能算法构建,其中数据治理工作尤为重要,治理后的数据质量是智能算法能否发挥其价值的关键.智能算法中“数据+物理”协同驱动主要存在四种引入机制,分别是基于输入数据的嵌入机制、基于模型架构的嵌入机制、基于损失函数的嵌入机制和基于模型优化的嵌入机制.以油气领域应用为例,介绍了数智流体力学在储层物性参数预测、压裂效果评价以及注采参数优化等方面的一系列研究进展.数智流体力学是流体力学未来的重要发展方向之一,以场景需求为导向、深度融合物理信息等先验知识的新一代智能理论与方法是数智流体力学发展的必然趋势,能够从崭新的角...     

边界元特征值分析及Burton-Miller法探究BEM eigenvalue analysis and the investigation of the Burton-Miller method

运用围道积分方法将边界元非线性特征值问题转化为规模很小的广义特征值问题,从而构造出一种边界元特征值分析方法。数值算例验证了该方法的求解精度。针对外声场问题,通过对常规、法向导数和Burton‐Miller边界积分方程的虚假特征频率的计算和比较,揭示了Burton‐Miller法规避虚假特征频率的本质,并对其中的叠加常数的最优取值给出了一种新的解释。