特殊坐标系中特殊非牛顿流边界层方程的相似解

给出了在一个特殊坐标系中三阶流体的二维定常运动方程组．该坐标系中由无粘流体的势流确定，即以环绕任意物体的非粘性流动的流线为Ф-坐标，速度势线为ψ-坐标，构成正交曲线坐标系．结果表明，边界层方程与浸没在流体中的物体的形状无关．第一次近似假定第二梯度项与粘性项和第三梯度项相比，可以忽略不计．第二梯度项的存在，将防碍第三梯度流相似解的比例变换的导出．利用李群方法计算了边界层方程的无穷小生成元．将边界层方程组变换为常微分方程组．利用Runge-Kutta法结合打靶技术求解了该非线性微分方程组的数值解．     

二阶流体通过径向延伸平面时滑移、黏性耗散、焦耳热对MHD流动的影响

研究了二阶导电的非Newton流体,在一个可径向放射状延伸,并伴有部分滑动表面上的流动及其热交换.部分滑移用一个无量纲的滑移因子控制,其取值范围从0(全黏着)到无穷大(全滑移).使用适当的相似变换,把待求的非线性偏微分方程转化为常微分方程.讨论了边界条件的不足,在无需增加任何边界条件下,使用有效的数值格式,求解所得到的微分方程.部分滑移、磁场交互参数以及二阶流体的参数对速度场和温度场的综合分析发现,滑移量的增加,流体的动力边界层和热边界层增厚.因为当滑移量的增加,允许更多的流体通过该平面,表面摩擦因数的数值下降,并在更高的滑移参数下,摩擦因数趋于0,即流体无黏性地通过.还研究了磁场对速度场和温度场的重要影响.     

应用逐次Taylor级数线性化法在多孔伸展界面上求解高对流Maxwell流体的磁流体动力学流动

研究不可压缩的高对流Maxwell(UCM)流体,在多孔伸展界面上作磁流体动力学(MHD)的边界层流动.利用相似变换将控制的偏微分方程,变换为非线性常微分方程.采用逐次Taylor级数线性化方法(STLM)求解该非线性问题.对所显现的参数完成速度分量的计算,介绍了表面摩擦因数的数值结果,并分析了问题所显现参数的变化.     

边界元方法的抽象误差估计及其应用

§1.引言 边界元方法是近二十年来发展的一种求解偏微分方程的数值方法,其基本思想是:先利用Green公式或位势将区域上的偏微分方程转化成边界上的积分方程,此时偏微分方程的解由边界积分方程的解表出;然后数值求解边界积分方程,进而求得偏微分方程的近     

一阶双曲型偏微分方程的模糊边界控制

对于一类半线性的双曲型偏微分方程的模糊边界控制问题,通过模糊控制方法,将半线性的偏微分方程系统精确表示为T-S模糊偏微分方程模型.因为控制器仅仅分布于边界上,所以基于T-S模糊偏微分方程模型而设计的模糊边界控制器将更容易执行,并且能够保证闭环系统指数稳定.然后利用Lyapunov方法将给出的闭环系统指数稳定的充分条件转化为求解线性不等式的问题.最后,通过仿真实例说明了模糊边界控制的有效性.     

Ladyzhenskaya模型的非协调有限元逼近

§1.引言 Navier-Stokes方程是描述粘性不可压缩流体运动的偏微分方程,它是研究这类非线性问题的很好的数学模型。但是,当速度梯度较大时,方程的整体解是否唯一可解,这个问题尚未得到解决。为此,对于不可压缩粘性流体的定常情形,Ladyzhenskaya提出利用下面的模型来代替Navier-Stokes方程。 设Ω是R~n(n=2或3)中的有界区域,边界?ΩLipshitz连续,u是流体速度,p是     

外部绕流的Navier-Stokes方程的边界层方程和维数分裂法

采用基于物体表面二维曲面的半测地坐标系(S-coordinate)建立了一个新的外部绕流边界层方程(boundary layer equations,BLE).BLE是一个关于物体的未知法向粘性应力张量和压力的非线性偏微分方程,其解的存在性得到了证明.此外,通过在二维流形上应用若干个2D-3C偏微分方程组来近似Navier-Stokes方程,获得了三维Navier-Stokes方程的维数分裂法.最后,对球和椭球的外部绕流问题给出了算例.     

基于简化的应变梯度理论下Kirchhoff板模型边值问题的提法及其应用北大核心CSCD

考虑应变梯度和速度梯度的影响,建立薄板控制微分方程及给出其边值问题的提法,修正了前人给出的薄板角点条件.采用Levy法,给出受分布力作用下简支板的挠度及自由振动频率的解析解.通过与文献中分子动力学数据对比,验证了该文模型的有效性并提出校核材料参数的一种方法.研究结果表明,增大弹性地基和应变梯度参数可以有效提高板的等效刚度,而速度梯度参数则相反.该文提出的板的边值问题为研究薄板在复杂支撑边界及外荷载等条件提供了理论依据.同时,有望为其有限元法、有限差分法和基于能量原理的Galerkin法等数值方法提供理论依据.     

混合模型下具有动态违约边界的债券定价

在混合模型下,研究了具有动态违约边界的公司债券定价问题.首先利用风险中性定价原理建立此定价问题的数学模型.然后,应用函数代换技巧和偏微分方程镜像法给出模型的显式解.最后,通过一个算例分析动态违约边界对公司债券价格的影响.结果表明:通过调整违约边界的相关参数值,可以得到不同形状的债券价格曲线,进而控制风险或得到更高的债券收益率.     

考虑感应磁场影响时,伸展表面上的MHD驻点流动及其热传递

考虑感应磁场的影响,研究不可压缩粘性流体在伸展表面上,作稳定磁流体动力学(MHD)的驻点流动.通过相似变换,将非线性的偏微分方程,变换成为常微分方程.用打靶法数值地求解变换后的边界层方程,得到不同的磁场参数和Prandtl数Pr时的数值解.对a/c>1和a/c<1两种情况(其中a和c均为正值),讨论感应磁场参数对表面摩擦因数、局部Nusselt数、速度和温度的影响,绘出变化曲线并给予讨论.     

Shape optimization of a body in the Oseen flow

This article is concerned with a numerical simulation of shape optimization of the Oseen flow around a solid body. The shape gradient for shape optimization problem in a viscous incompressible flow is computed by the velocity method. The flow is governed by the Oseen equations with mixed boundary conditions containing the pressure. The structure of continuous shape gradient of the cost functional is derived by using the differentiability of a minimax formulation involving a Lagrange functional with a function space parametrization technique. A gradient type algorithm is applied to the shape optimization problem. Numerical examples show that our theory is useful for practical purpose and the proposed algorithm is feasible. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Perturbed integral equation method on determining unknown geometry in fluid flow

We consider an inverse problem arising in fluid flow. An algorithm to find the shape of a body in uniform flow is proposed when the tangential velocity on its boundary is given a priori. The fluid flow is assumed to be inviscid, incompressible and irrotational.The essential idea to develop our algorithm is the boundary modification process toward the solution shape with the help of the perturbed integral equations. The perturbed integral equations are derived from the boundary perturbation. We also give examples exhibiting the reliability for our proposed algorithm.     

Shape inverse problem for Stokes–Brinkmann equations

In this paper, we propose an imaging technique for the detection of porous inclusions in a stationary flow governed by Stokes–Brinkmann equations. We introduce the velocity method to perform the shape deformation, and derive the structure of shape gradient for the cost functional based on the continuous adjoint method and the function space parametrization technique. Moreover, we present a gradient-type algorithm to the shape inverse problem. The numerical results demonstrate the proposed algorithm is feasible and effective for the quite high Reynolds numbers problems.     

Heat transfer in a Walter’s liquid B fluid over an impermeable stretching sheet with non-uniform heat source/sink and elastic deformation

In this present article an analysis is carried out to study the boundary layer flow behavior and heat transfer characteristics in Walter’s liquid B fluid flow. The stretching sheet is assumed to be impermeable, the effects of viscous dissipation, non-uniform heat source/sink in the presence and in the absence of elastic deformation (which was escaped from attention of researchers while formulating the viscoelastic boundary layer flow problems)on heat transfer are addressed. The basic boundary layer equations for momentum and heat transfer, which are non-linear partial differential equations, are converted into non-linear ordinary differential equations by means of similarity transformation. Analytical solutions are obtained for the resulting boundary value problems. The effects of viscous dissipation, Prandtl number, Eckert number and non-uniform heat source/sink on heat transfer (in the presence and in the absence of elastic deformation) are shown in several plots and discussed. Analytical expressions for the wall frictional drag coefficient, non-dimensional wall temperature gradient and non-dimensional wall temperature are obtained and are tabulated for various values of the governing parameters. The present study reveals that, the presence of work done by deformation in the energy equation yields an augment in the fluid’s temperature.     

A Mixed Analytical/Numerical Method for Velocity and Heat Transfer of Laminar Power-Law Fluids

This paper presents a relatively simple numerical method to investigate the flow and heat transfer of laminar power-law fluids over a semi-infinite plate in the presence of viscous dissipation and anisotropy radiation. On one hand, unlike most classical works, the effects of power-law viscosity on velocity and temperature fields are taken into account when both the dynamic viscosity and the thermal diffusivity vary as a power-law function. On the other hand, boundary layer equations are derived by Taylor expansion, and a mixed analytical/numerical method (a pseudo-similarity method) is proposed to effectively solve the boundary layer equations. This method has been justified by comparing its results with those of the original governing equations obtained by a finite element method. These results agree very well especially when the Reynolds number is large. We also observe that the robustness and accuracy of the algorithm are better when thermal boundary layer is thinner than velocity boundary layer.     

Hele Shaw Convection with Imposed Shear Flows: Boundary-Layer Formulation

The nonlinear convection forced by the boundaries of a Hele Shaw cell to align perpendicular to an imposed shear flow was analytically investigated by the boundary-layer method. The imposed shear flow may be a Couette flow that extends throughout the convecting layer or flow confined to a boundary, depending on the geometry of the Hele Shaw cell. This study examined the case in which the imposed shear flow has a boundary-layer structure and its interaction with the convecting interior. Analytical solutions for both the boundary layer and interior were obtained. The study revealed the following.For large aspect ratio A , the interaction of the imposed shear flow and convection is confined to the boundary layer. The boundary layer is a viscous rather than a thermal layer. The results showed that the range of validity of the Hele Shaw equations used in the literature is of order 1/ A ². For an asymptotically large aspect ratio A up to order 1/ A ², the velocity in the y -direction must be zero. The velocity in the x -direction and the z -direction has a parabolic dependence on y , but the temperature perturbation does not depend on y . These results may have implication for convection in porous media.     

Drag minimization for Navier–Stokes Flow

This paper investigates the drag minimization in a two‐dimensional flow which is governed by a nonhomogeneous Navier–Stokes equations. Two approaches are utilized to derive shape gradient of the cost functional. The first one is to use the shape derivative of the fluid state and its associated adjoint state; the second one is to utilize the differentiability of a minimax formulation involving a Lagrange functional with a function space parametrization technique. Finally, a gradient type algorithm is effectively formulated and implemented for the mentioned drag minimization problem. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A fundamental solution method for three-dimensional viscous flow problems with obstacles in a periodic array

In this paper, we propose a fundamental solution method for three-dimensional viscous flow problems with obstacles in a periodic array. Our problem is mathematically a boundary value problem of the Stokes equation with periodic boundary conditions, to which it is difficult to give a good approximation by the ordinary fundamental solution method. Our method gives an approximate solution by a linear combination of the periodic fundamental solutions. In addition, we can compute the drag forces on the obstacles by using the data obtained in our method. Numerical examples for the problems of flows past spheres show the effectiveness of our method.     

Development of a high order discretization scheme for solving fully nonlinear magnetohydrodynamic equations

We have developed fully fourth order accurate compact finite difference discretization scheme for the Navier-Stokes equations coupled with Maxwell''s equations. The implementation is done in cylindrical polar geometry. Due to the full-MHD modeling of physical flow, the modeled equations are fully nonlinear coupled hydrodynamic equations which are again coupled with Maxwells equations. In our computations, we have accounted for the induced magnetic field in the flow of an electrically conducting fluid in an external magnetic field. The code is tested against available experimental and theoretical data where applicable. It is observed that a smaller grid of $64 \times 64$ is sufficient for weakly nonlinear problems and higher grids up to $512 \times 512$ are needed as the degree of nonlinearities grow in the modeled equation. In the absence of magnetic field, a discontinuity of total drag coefficient and separation length is noted for $Re=73$ which is in agreement with literature. When the magnetic Reynolds number $Rm<1$ separation length decreases linearly with strength of magnetic field on a log-log scale whereas if $Rm>1$, it decreases nonlinearly, at a much faster rate. Thermal boundary layer thickness decreases as the strength of magnetic field increases and it forces the thermal convection to take place in a laminar structure as observed from thermal contour lines. Finally, using divided differences, we establish that the accuracy of the proposed numerical scheme is in fact fourth order.     

复杂边界旋转Navier-Stokes方程的几何方法及二度并行算法

本文提出了一种求解复杂边界旋转Navier-Stokes方程的微分几何方法及其二度并行算法.此方法可用于求解透平机械内部叶片间流动和飞行器外部绕流等复杂流动问题.假设流动区域可以用一系列光滑曲面■_k,k=1,2,…,K分割为一系列子区域(称作流层),通过应用微分几何的方法,三维N-S算子可以分解为两类算子之和:建立在曲面■_k切空间上"膜算子"和曲面■_k法线方向的"挠曲算子",将挠曲算子应用欧拉中心差商来逼近,由此得到建立在■_k上的"2D-3C"N-S方程.求解2D-3C N-S方程并且反复迭代直到收敛.我们得到"二度并行算法",它是2D-3C N-S方程并行算法与k方向的同时并行.这个算法的优点在于,(1)可以改进由于复杂边界造成的不规则三维网格引起的逼近解的精度;(2)为克服边界层的数值效应,在边界层内可以构造很密的流层,形成三维多尺度的网格,是一个很好的边界层算法;(3)这个方法不同于经典的区域分解算法,这里的每个子区域只需要求解一个"2D-3C"N-S方程,而经典区域分解方法要在每个子区域上求解三维问题.