NONLINEAR MOTION MECHANICS MODEL OF CURVED BLOOD VESSEL

A geometric model of curved blood vessels is established based on some reason-able hypotheses;the nonlinear motion mechanics model of the curved blood vessel is establishedaccording to basic mechanics laws.This model includes much more physiological factors.It cou-ples the interaction of blood flow with mechanical factors such as the displacement,deformation,strain and stress etc.of the curved blood vessel.It is of great importance for investigating thecirculation rules of the cardiovascular system and the nonlinear pulse wave propagation in curvedblood vessels.     

Wave propagation in a non-ideal dusty gas

In this paper we have studied the behavior of wave motion as propagating wavelets and their culmination into shock waves in a non-ideal gas with dust particles. In the absence of non-ideal effect the gas satisfies an equation of state of Mie–Gruneisen type. An expansion wave resulting from the action of receding piston is considered and the solutions to this problem showing effects of dust particles and non-idealness are obtained. The propagation of weak waves is considered and the flow variables in the region bounded by the piston and the characteristic wave front are found out. The expansive action of a receding piston undergoing an abrupt change in velocity is discussed. Cases of central expansion fan and shock fronts are studied and the solutions up to first order in the physical plane are obtained. The effects of non-idealness and dust particles are discussed in each case.     

On the structure of MHD shock waves in a viscous non-ideal gas

The structure of one-dimensional magnetohydrodynamics (MHD) shock waves is studied using the Navier–Stokes equations for the non-ideal gas phase. The exact solutions are obtained for the flow variables (i.e. particle velocity, temperature, pressure and change-in-entropy) within the shock transition region. The equation of state for a non-ideal gas is considered as given by Landau and Lifshitz. The effects of the non-idealness parameter and coefficient of viscosity of the gas are analysed on the flow variables assuming the magnetic field having only constant axial component. The findings confirm that the thickness of MHD shock front increases with decreasing values of the non-idealness parameter.     

Low-frequency axisymmetric waves in blood vessels of constant cross-section: an asymptotic approach

The asymptotic methods of shell theory are used to study the propagation of axisymmetric waves in blood vessels of constant cross-section. The initial equations are simplified using the assumption that the shell radius is small compared with the wave length. We show that the terms corresponding to the shell inertia cannot be omitted if it is required to describe not only the pressure wave but also the longitudinal wave. We study the influence of external fixation on the pressure wave. In this case, we compare the following two models: in the first model, the ambient medium is modelled by elastic and damping elements uniformly distributed over the shell exterior surface and by additional masses; in the second model, the ambient medium is represented by an infinite elastic space with a cylindrical cavity where the vessel is placed. On the interface between the elastic space and the vessel, we pose the full contact conditions. We show that, from the qualitative standpoint, both models lead to the same result: the pressure wave in the first approximation is a wave in the shell whose walls cannot move in the longitudinal direction. We asymptotically integrate the original equations and hence obtain a one-dimensional equation for the bulk blood flow.     

The mathematical model of the radial artery blood pressure waveform through monitoring of the age-related changes

Age-related changes in blood vessels affect the pulse wave propagation. These changes may cause an increase in wave reflection which leads to amplification of pulse pressure. The pulse pressure changes are associated with certain vascular diseases. Here we present a mathematical model of blood pressure for different age groups. The model is based on one-dimensional wave propagation theory and assumes that the pressure waveform is a superposition of forward propagating wave and backward waves from many reflection sites. The model is based on experimental data obtained by direct measurement of radial artery blood flow. The model clearly shows age-related changes in blood pressure waveform. The results of mathematical model correlate with the radial pressure waveform data. The model can be used in cases where it is not possible to measure the pressure due to movement of subject. Application of model to the direct blood flow measurements data allows the real-time pressure waveform monitoring. Furthermore, this approach enables monitoring of changes in pressure waveform due to the effects of medications.     

Outflow boundary conditions for one-dimensional finite element modeling of blood flow and pressure waves in arteries

Flow and pressure waves, originating due to the contraction of the heart, propagate along the deformable vessels and reflect due to tapering, branching, and other discontinuities. The size and complexity of the cardiovascular system necessitate a “multiscale” approach, with “upstream” regions of interest (large arteries) coupled to reduced-order models of “downstream” vessels. Previous efforts to couple upstream and downstream domains have included specifying resistance and impedance outflow boundary conditions for the nonlinear one-dimensional wave propagation equations and iterative coupling between three-dimensional and one-dimensional numerical methods. We have developed a new approach to solve the one-dimensional nonlinear equations of blood flow in elastic vessels utilizing a space-time finite element method with GLS-stabilization for the upstream domain, and a boundary term to couple to the downstream domain. The outflow boundary conditions are derived following an approach analogous to the Dirichlet-to-Neumann (DtN) method. In the downstream domain, we solve simplified zero/one-dimensional equations to derive relationships between pressure and flow accommodating periodic and transient phenomena with a consistent formulation for different boundary condition types. In this paper, we also present a new boundary condition that accommodates transient phenomena based on a Green’s function solution of the linear, damped wave equation in the downstream domain.     

脑Willis环的一维血流动力学及氧输运特性的数值研究

Willis环是大脑侧枝循环的重要组成部分, 研究其血流动力学特性以及氧输运规律对脑缺血疾病的认知和预防有着非常重要的作用. 该文旨在利用一维血流动力学模型模拟整个Willis环的流量变化和压力分布, 并建立动脉内氧输运的一维模型以模拟Willis 环内氧分压的变化规律, 为研究脑组织内血液流动和氧输运打下基础. 首先, 基于弹性圆管内的一维非线性流动方程和状态方程建立血流动力学模型, 在一维对流扩散方程的基础上, 考虑由管腔向壁面的扩散和壁面细胞的新陈代谢消耗推导出氧输运特性方程. 通过 Lax-Wendroff两步法对血流动力学方程进行离散, 而在进行对流扩散方程的离散时, 则运用迎风格式. 通过数值计算得到了正常情况下Willis环各个血管任意位置的流量、压力和氧分压的变化曲线, 正常情况下各个位置的氧分压处于稳定的平衡状态. 最后, 还通过此模型进一步模拟了右侧颈内动脉狭窄对各个血管内流动的影响. 当狭窄程度达到80%时, 中脑动脉的流量和压力会明显下降, 造成其供应区域的血流减少. 同时, Willis环右侧血管内的氧分压会大大降低, 而左侧血管的氧分压会出现上升趋势, 但幅度要小于右侧血管降低的幅度.     

One dimensional steepening of waves in non-ideal relaxing gas

In this paper, we studied the behavior of different modes of wave propagation and breaking of wave front by employing the theory of singular surfaces in a plane and radially symmetric flow of a non-ideal relaxing gas. The one dimensional steepening of waves is considered and the transport equation for the jump discontinuity of velocity gradient is obtained. The effects of relaxation and van der Waals excluded volume of the medium on the jump discontinuity of velocity gradient are analyzed.     

An improved correlation of the pressure drop in stenotic vessels using Lorentz’s reciprocal theorem

A mathematical model of the human cardiovascular system in conjunction with an accurate lumped model for a stenosis can provide better insights into the pressure wave propagation at pathological conditions. In this study, a theoretical relation between pressure drop and flow rate based on Lorentz’s reciprocal theorem is derived, which offers an identity to describe the relevance of the geometry and the convective momentum transport to the drag force. A voxelbased simulator V-FLOW VOF3 D, where the vessel geometry is expressed by using volume of fluid(VOF) functions, is employed to find the flow distribution in an idealized stenosis vessel and the identity was validated numerically. It is revealed from the correlation that the pressure drop of NS flow in a stenosis vessel can be decomposed into a linear term caused by Stokes flow with the same boundary conditions, and two nonlinear terms. Furthermore, the linear term for the pressure drop of Stokes flow can be summarized as a correlation by using a modified equation of lubrication theory, which gives favorable results compared to the numerical ones. The contribution of the nonlinear terms to the pressure drop was analyzed numerically, and it is found that geometric shape and momentum transport are the primary factors for the enhancement of drag force. This work paves a way to simulate the blood flow and pressure propagation under different stenosis conditions by using 1D mathematical model.     

气压对压力容器前壁超高速撞击损伤的影响

充气压力容器在超高速撞击下的典型损伤包括穿孔及其边缘的裂纹失稳破坏,会导致气体泄漏或爆炸,内压对容器前壁损伤的影响仍不明确。以不同内压的球形铝合金充气压力容器为研究对象,开展了球形铝合金弹丸超高速撞击实验和数值模拟计算,分析了内充气体压强对前壁穿孔形貌特征、穿孔直径、孔边环向应力等的影响规律和影响机理,讨论了气体冲击波的传播行为及影响前壁穿孔边缘裂纹失稳破坏的机制。结果表明:前壁穿孔边缘内翻边形貌与内压相关,内压越高,弯折程度越轻;穿孔直径与内充气体压强正相关,但气体对孔径的影响远小于容器壁厚及撞击速度的影响;穿孔边缘使裂纹失稳破坏的环向拉应力不仅受到后壁反射冲击波的影响,也与容器壁内应力波的传播有关,与内压成正比。     

The nonlinear response of Cattaneo-type thermal loading of a laser pulse on a medium using the generalized thermoelastic model

The nonlinear thermoelastic responses of an elastic medium exposed to laser generated shortpulse heating are investigated in this article. The thermal wave propagation of generalized thermoelastic medium under the impact of thermal loading with energy dissipation is the focus of this research. To model the thermal boundary condition(in the form of thermal conduction),generalized Cattaneo model(GCM) is employed. In the reference configuration, a nonlinear coupled Lord-Shulman-type generalized thermoelasticity formulation using finite strain theory(FST) is developed and the temperature dependency of the thermal conductivity is considered to derive the equations. In order to solve the time-dependent and nonlinear equations, Newmark's numerical time integration technique and an updated finite element algorithm is applied and to ensure achieving accurate continuity of the results, the Hermitian elements are used instead of Lagrangian's. The numerical responses for different factors such as input heat flux and nonlinear terms are expressed graphically and their impacts on the system's reaction are discussed in detail.The results of the study are presented for Green–Lindsay model and the findings are compared with Lord-Shulman model especially with regards to heat wave propagation. It is shown that the nature of the laser's thermal shock and its geometry are particularly determinative in the final stage of deformation. The research also concluded that employing FST leads to achieving more accuracy in terms of elastic deformations; however, the thermally nonlinear analysis does not change the results markedly. For this reason, the nonlinear theory of deformation is required in laser related reviews, while it is reasonable to ignore the temperature changes compared to the reference temperature in deriving governing equations.     

非线性变形节理中纵波传播特性的理论研究

基于波前动量守恒理论和位移不连续方法所提出的时域分析新方法,引入岩石非线性法向本构关系,对弹性纵波在岩石非线性节理中的传播特性进行了理论分析。采用节理变形的双曲线模型(BB模型),获得纵波P波斜入射非线性节理的传播波动方程,并通过参数研究分析了在岩石节理中节理非线性系数、节理初始刚度、应力波入射角和入射波幅值等因素对纵波传播规律的影响。结果表明:所推导的应力波传播方程在考虑多种非线性问题时,通过迭代计算即可方便求出透射波和反射波的数值解,避免了复杂的数学运算;当波斜入射节理面时,产生了波型转换,节理变形的非线性对透射波和反射波有较大影响,透射系数和反射系数并非随着非线性参数的变化而单调变化。时域内所推导的波传播方程更有益于波斜入射时非线性参数的广泛研究,为开展该方面的理论研究工作提供了借鉴。     

血液-血管耦合特性与脉搏波传播特性的关系

脉搏波既不可简单地理解为可压缩血液流体中的压力纵波，也不可简单地理解为沿固体血管传播的涨缩位移横波，而是超乎普通想象的流-固耦合和纵波-横波耦合的复杂波。从分析耦合本构关系的新途径出发，本文中提出了一个流-固耦合/纵波-横波耦合的串联模型，可为解读“位数形势”中医脉诊提供更丰富的信息。结果表明，脉搏波耦合系统的等效体积压缩模量K_s以及相应的耦合系统脉搏波传播速度c_s主要依赖于两个无量纲参数:血液-血管模量比K_b(p)/E(p)和薄壁血管径厚比D(p)/h₀，它们因人而异、因人的不同脉搏位置而异。文中定量分析了它们对c_s的影响，显示人体的K_b/E值在103数量级，从而c_s值在100～101 m/s数量级，以适应人体生理生化反应。由临床有创测量，证实脉搏体积横波与脉搏压力纵波是相耦合地以相同速度传播；还显示脉搏波是在其波阵面上具有氧合生化反应的“生物波”。此外，还讨论了“脉压放大”现象与非线性本构关系和与血管分叉处加载增强反射之间的关系，并讨论了Lewis关于重搏波形成的假设。     

Application of a thermodynamically compatible two‐phase flow model to the high‐resolution simulations of compressible gas–magma flow

This paper reports on the application and development of a fully hyperbolic and fully conservative two‐phase flow model for the simulation of gas and magma flow within volcanic processes. The model solves a set of mixture conservation equations for the gas and magma two‐phase flow with velocity non‐equilibrium. In this model, the effect of the relative velocity is introduced by a kinetic constitutive equation with other equations for volume and mass fractions of the gas phase. The model is examined numerically by the widely used finite volume Godunov methods of centered‐type. Using the Riemann problem, we numerically simulate wave propagation and the development of shocks and rarefactions in volcanic eruptions. These simulations are of magma fragmentation type where the relative velocity continues to dominate. A series of test cases whose solution contains features relevant to gas–magma mixtures are conducted. In particular, numerical results indicate that the model implementation predicts key features of the relative velocity within volcanic processes without any mathematical or physical simplifications. Simulation results are sharply and accurately provided without any spurious oscillations in all of the flow variables. The numerical methods and results are also compared with other numerical methods available in the literature. It is found that the provided resolutions are more accurate for the considered test cases. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical investigation of fluid–particle interactions for embolic stroke

Roughly one-third of all strokes are caused by an embolus traveling to a cerebral artery and blocking blood flow in the brain. The objective of this study is to gain a detailed understanding of the dynamics of embolic particles within arteries. Patient computed tomography image is used to construct a three-dimensional model of the carotid bifurcation. An idealized carotid bifurcation model of same vessel diameters was also constructed for comparison. Blood flow velocities and embolic particle trajectories are resolved using a coupled Euler–Lagrange approach. Blood is modeled as a Newtonian fluid, discretized using the finite volume method, with physiologically appropriate inflow and outflow boundary conditions. The embolus trajectory is modeled using Lagrangian particle equations accounting for embolus interaction with blood as well as vessel wall. Both one- and two-way fluid–particle coupling are considered, the latter being implemented using momentum sources augmented to the discretized flow equations. It was observed that for small-to-moderate particle sizes (relative to vessel diameters), the estimated particle distribution ratio—with and without the inclusion of two-way fluid–particle momentum exchange—were found to be similar. The maximum observed differences in distribution ratio with and without the coupling were found to be higher for the idealized bifurcation model. Additionally, the distribution was found to be reasonably matching the volumetric flow distribution for the idealized model, while a notable deviation from volumetric flow was observed in the anatomical model. It was also observed from an analysis of particle path lines that particle interaction with helical flow, characteristic of anatomical vasculature models, could play a prominent role in transport of embolic particle. The results indicate therefore that flow helicity could be an important hemodynamic indicator for analysis of embolus particle transport. Additionally, in the presence of helical flow, and vessel curvature, inclusion of two-way momentum exchange was found to have a secondary effect for transporting small to moderate embolus particles—and one-way coupling could be used as a reasonable approximation, thereby causing substantial savings in computational resources.     

A hybrid approach for nonlinear computational aeroacoustics predictions

In many aeroacoustics applications involving nonlinear waves and obstructions in the far-field, approaches based on the classical acoustic analogy theory or the linearised Euler equations are unable to fully characterise the acoustic field. Therefore, computational aeroacoustics hybrid methods that incorporate nonlinear wave propagation have to be constructed. In this study, a hybrid approach coupling Navier–Stokes equations in the acoustic source region with nonlinear Euler equations in the acoustic propagation region is introduced and tested. The full Navier–Stokes equations are solved in the source region to identify the acoustic sources. The flow variables of interest are then transferred from the source region to the acoustic propagation region, where the full nonlinear Euler equations with source terms are solved. The transition between the two regions is made through a buffer zone where the flow variables are penalised via a source term added to the Euler equations. Tests were conducted on simple acoustic and vorticity disturbances, two-dimensional jets (Mach 0.9 and 2), and a three-dimensional jet (Mach 1.5), impinging on a wall. The method is proven to be effective and accurate in predicting sound pressure levels associated with the propagation of linear and nonlinear waves in the near- and far-field regions.     

Nonlinear dispersion of long internal waves

The propagation of long weakly nonlinear waves in an atmospheric waveguide is considered. A model system of Kadomtsev-Petviashvili equations [1], which describes the propagation of such waves, is derived. In the case of one excited wave mode the system of model equations goes over into the Kadomtsev-Petviashvili equation, in which, however, the variables x and t are interchanged. The reasons for this are clarified. In the two-dimensional case an approximate solution of the model equations is constructed, and steady nonlinear waves and their interaction in a collision are considered. The results of a numerical verification of the stability of the approximate steady solutions and of the solution to the problem of decay of the wave into quasisolitons are given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 151–157, May–June, 1988.     

Numerical simulations of inviscid and viscous free‐surface flows with application to nonlinear water waves

The purpose of the present study is to establish a numerical model appropriate for solving inviscid/viscous free‐surface flows related to nonlinear water wave propagation. The viscous model presented herein is based on the Navier–Stokes equations, and the free‐surface is calculated through an arbitrary Lagrangian–Eulerian streamfunction‐vorticity formulation. The streamfunction field is governed by the Poisson equation, and the vorticity is obtained on the basis of the vorticity transport equation. For computing the inviscid flow the Laplace streamfunction equation is used. These equations together with the respective (appropriate) fully nonlinear free‐surface boundary conditions are solved using a finite difference method. To demonstrate the model feasibility, in the present study we first simulate collision processes of two solitary waves of different amplitudes, and compute the phenomenon of overtaking of such solitary waves. The developed model is subsequently applied to calculate (both inviscid and the viscous) flow field, as induced by passing of a solitary wave over submerged rectangular structures and rigid ripple beds. Our study provides a reasonably good understanding of the behavior of (inviscid/viscous) free‐surface flows, within the framework of streamfunction‐vorticity formulation. The successful simulation of the above‐mentioned test cases seems to suggest that the arbitrary Lagrangian–Eulerian/streamfunction‐vorticity formulation is a potentially powerful approach, capable of effectively solving the fully nonlinear inviscid/viscous free‐surface flow interactions. Copyright © 2009 John Wiley & Sons, Ltd.     

Lung effect on the hemodynamics in pulmonary artery

The present study investigates blood flow in a pulmonary artery. The aim is to gain a better understanding of offset value in vascular circulation through a two‐dimensional analysis of the Navier–Stokes equations. In this study, the hemodynamics in a blood vessel with truncated outlets at which constant pressure is specified is examined. To simplify the analysis, the vessel walls are regarded as being rigid. In quadratic elements, the streamline upwind Petrov–Galerkin finite element model is employed to simulate the incompressible Newtonian blood flow. The adopted finite element model introduces artificial damping terms solely in the streamline direction. With these terms added to the formulation, the discrete system is enhanced while solution accuracy is maintained without deterioration due to numerical diffusion errors. Copyright © 2001 John Wiley & Sons, Ltd.     

On the Interphase Heat Transfer Effect on Nonlinear Wave Propagation in a Gas-Liquid Mixture

The effect of interphase heat transfer on shock wave propagation is investigated. A multiwave nonlinear equation which in the limiting case of the absence of heat transfer decomposes into two classic generalizations of the Boussinesq equations is derived. Quasi-isothermal and quasi-adiabatic propagation regimes for which the heat transfer is fairly intense are considered. For both regimes, nonlinear equations describing the wave propagation are obtained. The equation describing the first regime is investigated in detail. Exact analytic solutions of this equation are given and used to study the shock wave structures and the solitary wave behavior. Formulas for the dependence of the heat transfer rate on the equilibrium-mixture parameters are obtained.