瞬态载荷下L4/L5椎间盘内流固耦合效应

中医正脊治疗通过对腰椎施加瞬态拉伸和旋转来治疗腰椎间盘退变, 本文采用考虑流固耦合效应的数值模拟研究其生物力学机制. 通过实验测量和文献调研, 确定了合理的拉伸和旋转的载荷参数; 发展了使用人体断层扫描图像结合解剖学数据建立详细腰椎几何模型的方法; 将松质骨、终板、椎间盘考虑为多孔弹性介质, 其他组织考虑为线性弹性介质, 进而建立了考虑生物组织中流固耦合效应的物理模型; 通过数值模拟得到了不同瞬态载荷及其组合作用下椎间盘内应力?应变与流体流动的变化规律. 研究发现, 瞬态载荷通过改变L4/L5椎间盘基质应力和髓核内外压力梯度, 在髓核中产生流体流动; 拉伸加载引起流体先流出髓核、再流入髓核, 产生含水量变化; 顺时针旋转加载在髓核左右产生相反的流动, 髓核右侧的含水量变化较左侧大. 本研究所采用的方法为流动过程相关的人体椎间盘退变病理生理机制研究提供了新的方法, 为中医正脊研究提供科学化思路, 也为相关的力学-生物学耦合研究和髓核再生的基础研究提供了一个切入点.      

A Rigorous Derivation of the Coupling of a Kinetic Equation and Burgers’ Equation

In this article, we investigate the kinetic/fluid coupling on a toy model, which we obtain rigorously from a hydrodynamical limit. The idea is that at the level of the full kinetic model, the coupling is obvious. We then investigate the coupling obtained when passing into the limit. We show that, especially in presence of a shock stuck on the interface, the coupling involves a kinetic layer known as the Milne problem. Due to this layer, the limit process is quite delicate and some blow-up techniques are needed to ensure its strong convergence.     

The flow field induced by the torsional oscillations of a spherical cell containing a fluid drop

The flow field induced by the torsional oscillations of a spherical cell containing a fluid drop is examined. It has been found that in addition to the oscillating motion of each fluid particle there exist three standing vortices in every quadrant between the drop and the container. The steady streaming into the drop can be directed either clockwise or counterclockwise, depending on the values of the parameter of the fluids inside and outside the drop. Typical flow fields are shown graphically.     

An objective rotation tensor applied to non-Newtonian fluid mechanics

To develop objective constitutive equations, local frames which translate and rotate with the fluid particle can be used. For example, the corotating frame rotates such that the curl of the velocity calculated in this frame vanishes. From the corotating frame, the Jaumann derivative can be derived. In this paper, a new local frame is developed which causes the cross product of the velocity and acceleration to vanish and is designated as the rigid-rotating frame. The corotating and rigid-rotating frames rotate identically for a rigid-body rotation of the fluid, but rotate differently in flows that contain shearing. This difference in rotation can be used to develop an objective rotation tensor that can be applied to constitutive equations for viscoelastic liquids. The rigid-rotating frame can also be used to develop a rheological time derivative which has been designated the rigid-rotating derivative. These new quantities expand the traditional set of kinematical variables and invariants available for use in constitutive equations. Use of this expanded set of kinematic variables is demonstrated in limiting constitutive equations. Received: 1 March 1999 Accepted: 5 March 1999     

An immersed boundary method for incompressible flows using volume of body function

A simple and effective immersed boundary method using volume of body (VOB) function is implemented on unstructured Cartesian meshes. The flow solver is a second‐order accurate implicit pressure‐correction method for the incompressible Navier–Stokes equations. The domain inside the immersed body is viewed as being occupied by the same fluid as outside with a prescribed divergence‐free velocity field. Under this view a fluid–body interface is similar to a fluid–fluid interface encountered in the volume of fluid (VOF) method for the two‐fluid flow problems. The body can thus be identified by the VOB function similar to the VOF function. In fluid–body interface cells the velocity is obtained by a volume‐averaged mixture of body and fluid velocities. The pressure inside the immersed body satisfies the same pressure Poisson equation as outside. To enhance stability and convergence, multigrid methods are developed to solve the difference equations for both pressure and velocity. Various steady and unsteady flows with stationary and moving bodies are computed to validate and to demonstrate the capability of the current method. Copyright © 2005 John Wiley & Sons, Ltd.     

均布载荷作用下具有光滑中心波纹膜片的非线性分析

采用轴对称旋转壳体的简化Reissner方程，研究了在均布载荷作用下具有光滑中心波纹膜片的非线性弯曲问题。应用格林函数方法，波纹膜片的非线性边值问题化为了非线性积分方程的求解。为了求解积分方程并防止发散，一个插值参数被引入到迭代格式中。计算表明，当载荷很小时，任何插值参数值均能保证迭代的收敛性，取插值参数值接近或等于1获得较快的收敛速度，而当载荷较大时，插值参数值不能取得过大。绘出了波纹膜片的特征曲线，得到的特征曲线可供设计参考。可以断言，当载荷不大时，特征曲线是近似线性的，随着载荷的增大，特征曲线开始向上弯曲，明显偏离线性。本文中提出的解决方法适应于任意轴向截面的波纹壳体。     

Global convergence of inexact Newton methods for transonic flow

In computational fluid dynamics, non-linear differential equations are essential to represent important effects such as shock waves in transonic flow. Discretized versions of these non-linear equations are solved using iterative methods. In this paper an inexact Newton method using the GMRES algorithm of Saad and Schultz is examined in the context of the full potential equation of aerodynamics. In this setting, reliable and efficient convergence of Newton methods is difficult to achieve. A poor initial solution guess often leads to divergence or very slow convergence. This paper examines several possible solutions to these problems, including a standard local damping strategy for Newton's method and two continuation methods, one of which utilizes interpolation from a coarse grid solution to obtain the initial guess on a finer grid. It is shown that the continuation methods can be used to augment the local damping strategy to achieve convergence for difficult transonic flow problems. These include simple wings with shock waves as well as problems involving engine power effects. These latter cases are modelled using the assumption that each exhaust plume is isentropic but has a different total pressure and/or temperature than the freestream.     

Nonlinear bending of corrugated diaphragm with large boundary corrugation under compound load

IntroductionCorrugateddiaphragmisatypeofelasticthinshells .Itsdesignisverycomplicatedbecauseoftoomanyparametersthatinfluenceeachother.Inanumberofinstrumentsmeasuringdisplacements,corrugateddiaphragmissubjectedtoelasticdisplacementthatisatleastthesameorderasitsthickness,sothatitisnecessarytousegeometricalnonlineartheoryofthinshellstoanalyze.Sofarasweknow ,inmostcases,investigatorsdiscussedonlytheproblemofcorrugateddiaphragmwithuniformanddensecorrugationsundertheactionofaunique(uniformlyorconcen…     

裂缝角度对含气多孔介质热流耦合特性影响研究

为研究含裂缝多孔介质内部流体的流动和传热特性,通过热-流耦合传热实验测试对5种裂缝倾角(0°,15°,30°,150°,165°)多孔介质进行了研究.结果表明:裂缝角度y=0°的多孔介质达到稳态需要时间最长,约为720min,传热速率最慢,适用于削弱传热应用;裂缝角度旋转方向对多孔介质内高低温流体的混合能力影响较大,裂...     

Divergent channel flow of a liquid crystalline polymer solution

 The flow of a `model' lyotropic liquid crystal polymer, (hydroxypropyl)cellulose in water, through a rectangular channel with a divergence in the channel width, is studied by in situ light microscopy. Microscopic texture observations are related to measurements of the flow velocity field, in order to characterize the shear and elongational aspects of the flow and to examine the effects of the divergence from a narrow channel to a wide channel. A strong dependence of flow-induced texture on position in the channel is observed and is related to the interplay of shear and elongational strain. The divergence generates both a perpendicular elongational strain due to the widening of the channel, and subsequently an elongational strain along the flow direction due to the change in flow pattern from quasi-radial to unidirectional down the wide channel. Additionally side wall structure is observed to be more complex than a simple strong alignment, displaying a fine birefringent texture. Finally there is a marked dependence of the macroscopic structure on the strain history of the fluid prior to entry into the channel, indicating that very different structures of, for instance, moulded parts, can result from differences in geometry and fluid treatment prior to entry into the mould itself. Received: 12 October 1999/Accepted: 29 October 1999     

含主应力轴旋转的广义塑性位势理论

大量岩土实验与工程实践表明传统塑性位势理论无法合理反映岩土材料的基本变形机制。从塑性位势理论角度来看,当存在主应力轴旋转时,塑性应变增量与应力不共主轴,此时最一般情况下的塑性应变增量须用六个线性无关的塑性势函数来表述,从而提出含主应力轴旋转的广义塑性位势理论一般表达式。通过矩阵分析,将一般应力增量分解成共轴分量与旋转分量之和。在应力增量分解的基础上,提出含主应力轴旋转的广义塑性位势理论的分解表达式     

CONSIDER SAINT-VENANT''''S PRINCIPLE BY MEANS OF CHAIN MODEL

ＩｎｔｒｏｄｕｃｔｉｏｎＡｃｃｏｒｄｉｎｇｔｏＬｏｖｅ[1],Ｓａｉｎｔ＿Ｖｅｎａｎｔ’ｓＰｒｉｎｃｉｐｌｅｃａｎｂｅｅｘｐｒｅｓｓｅｄａｓ:Ｔｈｅｅｆｆｅｃｔｏｆａｓｅｌｆ-ｅｑｕｉｌｉｂｒａｔｅｄｓｙｓｔｅｍｏｆｆｏｒｃｅｓａｃｔｉｎｇｉｎａｓｐｈｅｒｅｗｉｌｌｒｅｓｔｒｉｃｔｉｎａｓｍａｌｌｄｏｍａｉｎｎｅａｒｔｈｅｓｐｈｅｒｅ .Ｔｈｉｓｐｒｉｎｃｉｐｌｅｉｓｔｈｅｂｒｉｄｇｅｂｅｔｗｅｅｎｔｈｅｅｌａｓｔｉｃｉｔｙｔｈｅｏｒｙａｎｄａｐｐｌｉｃａｔｉｏｎ ,ｂｕｔｉｔｈａｓｎｏｔｂｅｅｎｐｒｏｖ…     

A new approach for computing the steady state fluid–structure interaction response of periodic problems

A special type of fluid–structure interaction (FSI) problems are problems with periodic boundary conditions like in turbomachinery. The steady state FSI response of these problems is usually calculated with similar techniques as used for transient FSI analyses. This means that, when the fluid and structure problem are not simultaneously solved with a monolithic approach, the problem is partitioned into a fluid and structural part and that each time step coupling iterations are performed to account for strong interactions between the two sub-domains. This paper shows that a time-partitioned FSI computation can be very inefficient to compute the steady state FSI response of periodic problems. A new approach is introduced in which coupling iterations are performed on periodic level instead of per time step. The convergence behaviour can be significantly improved by implementing existing partitioned solution methods as used for time step coupling (TSC) algorithms in the time periodic coupling (TPC) framework. The new algorithm has been evaluated by comparing the convergence behaviour to TSC algorithms. It is shown that the number of fluid–structure evaluations can be considerably reduced when a TPC algorithm is applied instead of a TSC. One of the most appealing advantages of the TPC approach is that the structural problem can be solved in the frequency domain resulting in a very efficient algorithm for computing steady state FSI responses.     

A convergent method of solution for clamped rectangular plate

The conventional method of solution for clamped rectangular plate lies in releasing the constraints for rotation at the clamped sides and substituting them by equivalent distributed bending moments. This paper gives a general formulation of this method for arbitrary load with proof of the convergence, and discusses the rate of convergence for displacements and bending moments. The method and its convergence can be extended to the case of arbitrary boundary conditions and to plate on elastic foundation.The author expresses his sincere gratitude to Prof. Hu Hai-chang for guidance in this work.     

变角度纤维复合材料层合斜板的颤振分析

假设纤维方向角沿层合板的长度方向线性变化,研究了变角度纤维复合材料层合斜板的颤振．通过坐标变换将斜板变换为正方形板,采用层合板表面连续变化的速度环量来模拟空气对其的作用,速度环量分布利用Cauchy积分公式计算．建立了系统的Lagrange方程并采用Ritz法得到了层合板的自振频率和颤振/不稳定性分离临界速度．通过数值算例验证了本文模型和方法的正确性和收敛性,分析了各个铺层内纤维方向角的变化对自振频率和颤振/不稳定性分离临界速度的影响．研究结果表明,通过纤维的变角度铺设,可有效地提高层合板的基频和颤振/不稳定性分离临界速度．经合理设计的变角度复合材料层合板具有抑制颤振的作用．     

Weakly-imposed Dirichlet boundary conditions for non-Newtonian fluid flow

We propose a new formulation for weakly imposing Dirichlet boundary conditions in non-Newtonian fluid flow. It is based on the Gerstenberger–Wall formulation for Newtonian fluids [1], but extended to non-Newtonian fluids. It uses a stabilization term in the weak form that is independent from the actual fluid model used, except for an adjustable parameter κ, having the physical dimension of a viscosity. The new formulation is tested, combined with an extended finite element method, for the flow past a cylinder between two walls using both a generalized Newtonian and a viscoelastic fluid. It is shown that the convergence is optimal for the generalized Newtonian fluid by comparing with a converged boundary-fitted solution using traditional strong boundary conditions. Also the solution of the viscoelastic fluid compares very well with a traditional solution using a boundary-fitted mesh and strong Dirichlet boundary conditions. For both fluid models we also test various values of the κ parameter and it turns out that a value equal to the zero-shear-viscosity gives good results. But, it is also shown that a wide range of κ values can be chosen without sacrificing accuracy.     

A hybrid scheme based on finite element/volume methods for two immiscible fluid flows

We have successfully extended our implicit hybrid finite element/volume (FE/FV) solver to flows involving two immiscible fluids. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix‐free implicit cell‐centered FV method. The pressure Poisson equation is solved by the node‐based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. This updating strategy can be rigorously proven to be able to eliminate the unphysical pressure boundary layer and is crucial for the correct temporal convergence rate. Our current staggered‐mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centers and the auxiliary variable at vertices. The fluid interface is captured by solving an advection equation for the volume fraction of one of the fluids. The same matrix‐free FV method, as the one used for momentum equations, is used to solve the advection equation. We will focus on the interface sharpening strategy to minimize the smearing of the interface over time. We have developed and implemented a global mass conservation algorithm that enforces the conservation of the mass for each fluid. Copyright © 2009 John Wiley & Sons, Ltd.     

The transient development of the flow in an impulsively rotated annular container

When a fluid-filled container is spun up from rest to a constant angular velocity the fluid responds in such a way that the fluid–container system is ultimately in a state of rigid-body rotation. The fluid can then be said to have traversed a trajectory in phase space from a simple stable equilibrium state of no motion to another stable equilibrium representing full rigid-body rotation. This simple statement belies the fact that during this process the fluid can undergo a series of transitions, from a laminar through a transient turbulent state, before attaining the stable motion that is rigid-body rotation. Using a combination of analytical and computational methods, we focus on the dynamics resulting from an impulsive change in the rotation rate of a fluid-filled annulus, specifically, the impulsive spin-up of a stationary annulus, or the impulsive spin-down of an annulus already in a state of rigid-body rotation. We explore the initial development of the impulsively generated axisymmetric boundary layer, its subsequent instability, and the larger-scale transient features within this class of flows, allowing us to look at the effect these features have on the time it takes for the system to spin up to a steady state, or spin down to rest.     

AN ALGEBRAIC MULTIGRID SOLVER FOR NAVIER–STOKES PROBLEMS IN DISCRETE SECOND–ORDER APPROXIMATION

An algebraic multigrid scheme is presented for solving the discrete Navier–Stokes equations to second-order accuracy using the defect correction method. Solutions for the driven cavity and asymmetric, sudden expansion test problems have been obtained for both structured and unstructured meshes, the resolution and resolution grading being controlled by global and local mesh refinements. The solver is efficient and robust to the extent that, for problems attempted so far, no underrelaxation of variables has been required to ensure convergence. Provided that the computational mesh can resolve the flow field, convergence characteristics are almost mesh-independent. Rates of convergence actually improve with refinement, asymptotically approaching mesh-independent values. For extremely coarse meshes, where dispersive truncation errors would be expected to prevent convergence (or even induce divergence), solutions can still be obtained by using explicit underrelaxation in the iterative cycle.     

Numerical approximation of generalized Newtonian fluids using Powell–Sabin–Heindl elements: I. theoretical estimates

In this paper we consider the numerical approximation of steady and unsteady generalized Newtonian fluid flows using divergence free finite elements generated by the Powell–Sabin–Heindl elements. We derive a priori and a posteriori finite element error estimates and prove convergence of the method of successive approximations for the steady flow case. A priori error estimates of unsteady flows are also considered. These results provide a theoretical foundation and supporting numerical studies are to be provided in Part II. Copyright © 2003 John Wiley & Sons, Ltd.