Coupled modeling of tumour angiogenesis,tumour growth,and blood perfusion

This paper proposes a more realistic mathematical simulation method to investigate the dynamic process of tumour angio-genesis by fully coupling the vessel growth, tumour growth and associated blood perfusion. The tumour growth and angiogenesis are coupled by the chemical microenvironment and the cell-matrix interaction. The haemodynamic calculation is carried out on the new vasculature, and an estimation of vessel collapse is made according to the wall shear stress criterion. The results are consistent with physiological observations, and further confirm the application of the coupled model feedback mechanism. The model is available to examine the interactions between angiogenesis and tumour growth, to study the change in the dynamic process of chemical environment and the vessel remodeling.     

Numerical simulation of solid tumor angiogenesis with Endostatin treatment： a combined analysis of inhibiting effect of anti-angiogenic factor and micro mechanical environment of extracellular matrix

To investigate the influence of anti-angiogenesis drug Endostatin on solid tumor angiogenesis, a mathematical model of tumor angiogenesis was developed with combined influences of local extra-cellular matrix mechanical environment, and the inhibiting effects of Angiostatin and Endostatin. Simulation results show that Angiostatin and Endostatin can effectively inhibit the process of tumor angiogenesis, and decrease the number of blood vessels in the tumor. The present model could be used as a valid theoretical method in the investigation of anti-angiogenic therapy of tumors.     

Two-dimensional discrete mathematical model of tumor-induced angiogenesis

A 2D discrete mathematical model of a nine-point finite difference scheme is built to simulate tumor-induced angiogenesis. Nine motion directions of an individual endothelial cell and two parent vessels are extended in the present model. The process of tumor-induced angiogenesis is performed by coupling random motility, chemotaxis, and haptotaxis of endothelial cell in different mechanical environments inside and outside the tumor. The results show that nearly realistic tumor microvascular networks with neoplastic pathophysiological characteristics can be generated from the present model. Moreover, the theoretical capillary networks generated in numerical simulations of the discrete model may provide useful information for further clinical research.     

3D numerical simulation of avascular tumour growth: effect of hypoxic micro-environment in host tissue

A three-dimensional （3D） mathematical model of tumour growth at the avascular phase and vessel remodelling in host tissues is proposed with emphasis on the study of the interactions of tumour growth and hypoxic micro-environment in host tissues. The hybrid based model includes the continuum part, such as the distributions of oxygen and vascular endothelial growth factors （VEGFs）, and the discrete part of tumour cells （TCs） and blood vessel networks. The simulation shows the dynamic process of avascular tumour growth from a fewinitial cells to an equilibrium state with varied vessel networks. After a phase of rapidly increasing numbers of the TCs, more and more host vessels collapse due to the stress caused by the growing tumour. In addition, the consumption of oxygen expands with the enlarged tumour region. The study also discusses the effects of certain factors on tumour growth, including the density and configuration of pre- existing vessel networks and the blood oxygen content. The model enables us to examine the relationship between early tumour growth and hypoxic micro-environment in host tissues, which can be useful for further applications, such as tumour metastasis and the initialization of tumour angiogenesis.     

Numerical simulation of tumor-induced angiogenesis influenced by the extra-cellular matrix mechanical environment

To investigate tumor-induced angiogenesis under the influence of the mechanical environments inside and outside the tumor, mathematical model of tumor angiogenesis was developed. In the model, extra-cellular matrix （ECM） was treated as a thin plane. The displacement of ECM is obtained from the force balance equation consisted of the ECs traction, the ECM visco-elastic forces and the exter- nal forces. Simulation results show that a layered capillary network is obtained with a well vascularized region at the periphery of the tumor. The present model can be used as a valid theoretical method in the basic researches in tumorinduced angiogenesis.     

混凝土中化学-热-湿-力耦合过程的数值方法

提出了一个火灾下混凝土中化学-热-湿-力耦合过程分析的两级数学模型. 混凝 土模型化为充满两种非混溶孔隙流体的非饱和变形多孔多相介质. 数学模型基于控制干空 气、湿份及基质溶解物的质量守恒、混凝土介质混合体的动量守恒和焓(能量)守恒的耦合 偏微分方程组. 模型中特别考虑到了高温下的脱盐过程. 构造了一个用于数值模拟 化学-热-湿-力耦合行为的有限元求解过程的混合弱形式. 并且针对其中具有非自伴随算子特性的 双曲线控制方程的空间离散进行了特殊考虑. 数值结果例题显示所发展的数学模型和数值方 法在重现火灾条件下的混凝土中化学-热-湿-力耦合行为的有效性.     

A mathematical model for solid-gas coupled problems on the methane flowing in coal seam

Based on the fundamental theories of the solid deformation and methane flowing, a mathematical model for solid-gas coupled problem on methane flowing is developed. The smoothness and existence of the solution of this model are proved, and the numerical method for solving the present problem is presented.This project was supported by the Shanxi Province Natural Science Foundation.     

混凝土爆破损伤的SPH-FEM耦合法数值模拟

为了提高计算效率以及更好展现爆炸荷载下混凝土破坏过程,采用SPH-FEM耦合法对混凝土爆破成坑进行模拟。首先结合前人给出的C30混凝土Holmquist-Johnson-Cook（HJC）部分本构参数,通过理论推导等方法确定出剩余的参数;然后代入模型中计算,将数值解与实测数据进行对比;最后以峰值压力和峰值加速度作为考察对象,对HJC模型中21个参数敏感性进行分析。结果表明：SPH-FEM耦合法能直观地模拟爆炸荷载作用下爆坑的发展全过程,且能够较好地处理SPH边界问题;基于所给出的C30混凝土HJC本构参数,采用SPH-FEM耦合法对混凝土爆破破坏进行模拟,计算结果与实测数据吻合度高,表明HJC本构参数的确定具有合理性。此外,还发现HJC本构参数对爆破问题结果的敏感度各不相同,指出对峰值压力和峰值加速度均有较大影响的参数在确定的时候需引起足够的重视。

     

Numerical simulation of polyhedral crystal growth based on a mathematical model arising from non–local thermomechanics

In this work we present some results of the numerical simulation of the growth of a crystal from its melt, taking into account faceting. The simulation is based on a numerical solution of a three–dimensional generalized Stefan problem. That problem arises from a non–local thermomechanical theory applied to a continuous system with an interface and embodies ideas from the dislocation theory of crystal growth. In the model, the crystal surface is an isotherm and the growth velocity of a crystal face depends on the velocities of the other faces and on the whole crystal configuration as well as on the temperature gradient. A front fixing formulation of the model is considered. This is a conservative form of the Isotherm Migration Method [6, 7, 8, 9, 10, 11] in spherical coordinates. The numerical solution is based on an explicit finite difference discretization of the resulting non–linear equations. We develop a theoretical analysis of the interface equations that drive the crystal face motion. Numerical results, showing evolution of complex crystals with configuration changing during the growth, are in accord with experimental results. Furthermore, numerical experiments offer useful information on the influence of certain parameters in the model on the growth process. Received: March 21, 1996     

Numerical simulation of avascular tumor growth based on p27 gene regulation

A multi-scale continuous-discrete model based on the effects of the p27 gene control is built to simulate the avascular tumor growth. At the tissue level, the continuous Eulerian model is adopted to determine the distribution of the concentration of oxygen, the extracellular matrix (ECM), and the matrix-degradative enzyme (MDE). At the cellular level, the discrete Lagrangien model is adopted to determine the movement, the proliferation, and the death of single tumor cells (TCs). At the genetic level, whether a cell is committed to mitosis is determined by solving a set of equations modeling the effects of the p27 gene control. The avascular morphological evolution of the solid tumor growth is simulated, including the radius of the solid tumor, the number of the TCs, the oxygen distribution over time, and the inhibiting effect of the up-regulating p27 gene expression.     

Three-dimensional mathematical model and its application in the neighborhood of the Three Gorges Reservoir dam in the Yangtze River

The calculation of flow and sediment transport is one of the most important tasks in river engineering. The task is particularly difficult because a number of complex physical phenomena should be accounted for more realistically in a model with a predictive power. Three-dimensional calculations of river flow and suspended sediment transport are performed in this paper with application in the Three Gorges Reservoir in the Yangtze River. A period of 76 years after the dam is built is simulated and the results are compared with laboratory measurements obtained by Tsinghua University whereby the model is verified and calibrated. Generally speaking, the calculated results agree well with the experiments, demonstrating that the present model can be used for flow and sediment transport prediction in rivers. The project supported by the National Natural Science Foundation of China (50009004)     

Numerical simulation of blood flow and interstitial fluid pressure in solid tumor microcirculation based on tumor-induced angiogenesis

A coupled intravascular–transvascular–interstitial fluid flow model is developed to study the distributions of blood flow and interstitial fluid pressure in solid tumor microcirculation based on a tumor-induced microvascular network. This is generated from a 2D nine-point discrete mathematical model of tumor angiogenesis and contains two parent vessels. Blood flow through the microvascular network and interstitial fluid flow in tumor tissues are performed by the extended Poiseuille’s law and Darcy’s law, respectively, transvascular flow is described by Starling’s law; effects of the vascular permeability and the interstitial hydraulic conductivity are also considered. The simulation results predict the heterogeneous blood supply, interstitial hypertension and low convection on the inside of the tumor, which are consistent with physiological observed facts. These results may provide beneficial information for anti-angiogenesis treatment of tumor and further clinical research. The project supported by the National Natural Science Foundation of China (10372026).     

Numerical simulation of laminar natural convection in a laterally heated vertical cylindrical enclosure: application to crystal growth

Laminar natural convection has been studied in a laterally heated vertical cylindrical enclosure with a free insulated surface and a centrally located constant temperature wall at the top. These conditions are a simplification of the conditions existing in a Czochralski crystal pulling system. The laminar, axisymmetric flow of a Newtonian, constant physical properties fluid under Boussinesq’s approximation has been considered. Governing equations in primitive variable form are solved numerically by control volume method. SIMPLE algorithm due to Patankar has been used for the numerical simulation. The effects of the constant wall heat flux boundary condition at the side wall have been investigated whereas the bottom wall is considered to be insulated. Streamlines and isotherms are presented for various Rayleigh numbers and Prandtl numbers. Heat flux vectors through the melt are plotted for selected cases. The axial velocity and temperature variations at different horizontal sections of the crucible have been presented graphically to explain the transport processes inside the crucible. It has been observed that in case of low Pr and high Ra, flow separation occurs at the vertical wall of the crucible which leads to an oscillatory flow as Ra increases. The investigation has been extended to the oscillatory regime of flow in the zone of supercritical Rayleigh numbers and some unsteady results are also presented. Finally a heat transfer correlation has been developed for steady-state case.     

Cosserat连续体模型中本构参数对应变局部化模拟结果影响的数值分析

基于所发展的压力相关弹塑性Cosserat连续体模型及相应的数值方法,以一维剪切层及二维平板压缩问题为例,数值分析了Cosserat连续体模型中的本构参数Cosserat剪模、软化模量及内部长度参数对应变局部化数值模拟结果的影响.结果表明在一定取值范围内,Cosserat剪模对数值模拟结果几乎没有影响,并给出了具体数值计算时的取值范围;软化模量绝对值越大,后破坏段的荷载-位移曲线越陡,计算得到的剪切带宽度越窄;内部长度参数越大,后破坏段的荷载-位移曲线越平缓,计算得到的剪切带越宽.     

基于Vreman亚格子模型的有限开口空间内热驱动流数值模拟

利用基于Vreman亚格子模型的大涡模拟技术对有开口的单室和双室房间内热驱动流进行了数值模拟,利用函数分析法定量分析了模拟结果的准确性,并与Smagorinsky亚格子模型的模拟结果进行了比较.结果表明,Vreman和Smagorinsky亚格子模型的计算结果均能够满足工程的需求,但Vreman亚格子模型在开口附近区域的温度和U速度计算结果在整体上比Smagorinsky亚格子模型更接近实验值;Vreman亚格子模型未像Sma-gorinsky亚格子模型那样过高地估算壁面附近高温区域的粘性耗散;对于单室房间内热烟气层高度的预测,采用Vreman模型得到的计算结果准确性比Smagorinsky亚格子模型提高近50％.     

基于扰动状态概念与E-B模型的粗粒料力学行为模拟

考虑到基于复杂弹塑性模型与扰动状态概念的本构模型在数值积分时的困难,以Duncan-Chang(E-B)模型描述材料在相对完整状态(RI)下的响应,用临界状态模型描述材料在完全调整状态(FA)下的力学行为。依据常规三轴试验曲线判断材料状态的转变并建议了相应的扰动因子计算方法。利用Abaqus软件实现了相应的程序代码,数值算例模拟了堆石料的三轴试验,结果表明邓肯-张模型与扰动状态概念的结合在模拟应变软化、剪缩与剪胀方面具有良好能力。在此基础上进一步模拟了某混凝土面板堆石坝的施工、蓄水、退水的过程,数值结果表明变形分布符合一般规律,且在定量上与原型观测吻合较好。