首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 109 毫秒
1.
带源参数的二维热传导反问题的无网格方法   总被引:2,自引:1,他引:1  
程荣军  程玉民 《力学学报》2007,39(6):843-847
利用无网格有限点法求解带源参数的二维热传导反问题,推导了相应的离散方程. 与 其它基于网格的方法相比,有限点法采用移动最小二乘法构造形函数,只需要节点信息,不 需要划分网格,用配点法离散控制方程,可以直接施加边界条件,不需要在区域内部求积分. 用有限点法求解二维热传导反问题具有数值实现简单、计算量小、可以任意布置节点等优点. 最后通过算例验证了该方法的有效性.  相似文献   

2.
在空间域上采用只与结点有关的无网格方法离散,在时间域上采用精细积分方法求 解. 无网格离散过程中,利用伽辽金积分等效弱形式代替微分形式的控制方程,并 用修正变分原理满足位移边界条件,采用移动最小二乘法求解离散的形函数,把形 函数代入等效积分弱形式得到离散的二阶方程;精细积分过程中非齐次项采 用Romberg积分. 同时给出了两种不同边界条件的谐响 应求解的两个数值算例,得到了精确的数值结果.  相似文献   

3.
本文采用形式解法研究三维有限元离散模型中波动的传播规律。文中由波动方程得到它的特征方程,进而给出了三维离散网格中波的频散关系,分析了三维网格中平面波、非匀匀波、寄生振荡及其它可能存在的运动形式的基本特点;详细讨论了平面波的频散、截止频率及偏振漂移等重要的不同于连续介质中的波动传播性质的特点,提示了其中不同于低维网格中 波动的现象。  相似文献   

4.
磁流体方程的数值求解在等离子体物理学、天体物理研究以及流动控制等领域具有重要意义,本文构造了用于求解理想磁流体动力学方程的基于移动网格的熵稳定格式,此方法将Roe型熵稳定格式与自适应移动网格算法结合,空间方向采用熵稳定格式对磁流体动力学方程进行离散,利用变分法构造网格演化方程并通过Gauss-Seidel迭代法对其迭代求解实现网格的自适应分布,在此基础上采用守恒型插值公式实现新旧节点上的量值传递,利用三阶强稳定Runge-Kutta方法将数值解推进到下一时间层。数值实验表明,该算法能有效捕捉解的结构(特别是激波和稀疏波),分辨率高,通用性好,具有强鲁棒性。  相似文献   

5.
两套节点格林元嵌入式离散裂缝模型数值模拟方法   总被引:1,自引:1,他引:0  
程林松  杜旭林  饶翔  曹仁义  贾品 《力学学报》2022,54(10):2892-2903
对于原始嵌入式离散裂缝模型(EDFM), 在计算包含裂缝单元的基质网格内的压力分布时采用了线性分布假设, 这导致了油藏开发早期对非稳态窜流量的计算精度不足. 因此, 本文提出了一种两套节点格林元法的EDFM数值模拟方法. 两套节点格林元法的核心思想是将压力节点与流量节点区分开, 一套压力节点设置在单元顶点, 另一套流量节点设置在网格边的中点, 满足局部物质守恒、具有二阶精度的同时, 可适用于任意网格类型. 本文将两套节点格林元法与EDFM耦合, 采用了非稳态渗流控制方程的边界积分形式推导了基质网格与裂缝网格之间传质量的新格式, 代替了线性分布假设以提高模拟精度; 此外, 修正后的EDFM能适应任意形态的基质网格剖分, 拓展了原始EDFM仅适用于矩形基质网格、难以考虑复杂油藏边界的局限性. 研究表明: 通过对比商业模拟软件tNavigator? LGR模块与原始EDFM, 验证了本文模型具有较高的早期计算精度; 以复杂油藏边界?缝网?SRV分区模型为例, 通过对比SFEM-COMSOL商业模拟软件, 验证了本文模型处理复杂问题的适应性. 本文研究可用于裂缝性油藏开发动态的精确模拟.   相似文献   

6.
磁流体流动在现代工业和科研中有着广泛的应用,但磁流体的流动受到磁场的影响,与一般流体区别较大,需要对其进行深入的研究。磁流体的流动受到流体力学流动方程和麦克斯韦方程的共同影响,其精确解在有限条件下才能得到,因此对磁流体的流动进行数值模拟具有重要的意义。本文采用移动最小二乘法计算形函数,利用无网格局部Petrov-Galerkin(MLPG)法得到控制方程的离散形式,在管壁为任意电导率及任意方向外加磁场的条件下,对方形直管道中定常流动的磁流体进行了数值计算。MLPG法的计算是基于节点的,不需要任何网格或单元,是一种真正的无网格方法。计算结果与Scheriff精确解进行了比较,表明该方法适用于中等以下哈特曼数的磁流体流动计算。  相似文献   

7.
将不规则区域嵌入到规则的矩形区域,在矩形区域上将弹性平面问题的控制方程采用重心Lagrange插值离散,得到控制方程矩阵形式的离散表达式。在边界节点上利用重心插值离散边界条件,规则区域采用置换法施加边界条件,不规则区域采用附加法施加边界条件,得到求解平面弹性问题的过约束线性代数方程组,采用最小二乘法进行求解,得到整个规则区域上的位移数值解。利用重心插值计算得到不规则区域内任意节点的位移值,计算精度可到10-14以上。数值算例验证了所建立方法的有效性和计算精度。  相似文献   

8.
对于已含初始裂纹平台管节点的寿命预测很大程度上依靠应力强度因子的精确值,而复杂载荷条件下的节点应力强度因子的计算尚无参数方程直接确定。本文提出了一种含表面裂纹的K节点的有限元网格产生方法,即把整个K节点划分为几个子区域,每个子区域的网格具有不同类型的单元和不同的密度。这种方法在控制网格密度,尤其是控制沿着裂纹边缘单元的边长比方面有其独特的优越性,当所有子区域的网格自动产生后,容易得到整个结构的有限元模型。同时用J积分和位移外推插值法分别计算了一个K型节点沿着裂纹前缘的应力强度因子值,发现:试验得到的应力强度因子值和提出的模型计算结果非常吻合,证明了所提有限元模型的准确性。  相似文献   

9.
彭林欣  李知闲  项嘉诚  覃霞 《力学学报》2022,54(12):3366-3382
在加肋板无网格模型中, 肋条的位置对各种工况下加肋板受力性能的影响至关重要. 文章基于一阶剪切变形和移动最小二乘法理论提出一种考虑非线性影响的加肋板无网格模型, 并利用遗传算法优化肋条位置. 首先, 采用离散节点分别对平板和肋条进行离散, 得到加肋板的无网格离散模型; 其次, 通过冯·卡门大挠度理论得到非矩形板几何非线性问题的弯曲控制方程; 再次, 通过哈密顿原理得到加肋非矩形板自由振动问题的控制方程; 最后引入遗传算法, 以肋条的位置为设计变量、非矩形加肋板中心点挠度最小或自振频率最大为目标函数, 对肋条位置进行优化. 在考虑了几何非线性影响的肋条位置优化过程中, 肋条位置改变时只需重新计算位移转换矩阵, 避免了网格重构. 本文以全局荷载下单肋条菱形板为例与理论解进行对比, 进行有效性验证. 再以板的中点挠度最小和自振频率最大为优化目标, 对局部荷载作用下不同形状、不同肋条布置方式的加肋板进行优化, 分析方法的收敛性及稳定性.   相似文献   

10.
平面弹性力学问题的离散元法   总被引:1,自引:0,他引:1  
根据离散元的基本原理,基于变形体的理论提出了适用于平面弹性力学问题的界面位移、应变和应力模式,建立了求解平面弹性力学问题的离散元方程和相应的迭代求解方法.通过界面位移可以简洁地将位移和力的边界条件引入离散系统的控制方程,也可以方便地求解节点位移.数值算例表明,与具有相同网格的有限元结果相比,离散元能同时给出精度相对较高的应力解和精度相当的位移解.  相似文献   

11.
The lattice‐BGK method has been extended by introducing additional, free parameters in the original formulation of the lattice‐BGK methods. The relationship between these parameters and the macroscopic moment equations is analysed by Taylor series and Chapman–Enskog expansion. The parameters are determined from the macroscopic moment equations by comparisons with the governing equations to be modelled. Extensions are presented for the Navier–Stokes equations at low Mach numbers in Cartesian or axisymmetric coordinates with constant or variable density, for scalar convection–diffusion equations and for equations of Poisson type. The generalized lattice‐BGK concept is demonstrated by two applications of chemical engineering. These are the computation of chemically reacting flow through an axisymmetric reactor and of the transport and deposition of particles to filters under the action of different forces. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

12.
A perturbation technique is proposed for solution of the generalized equations governing the thermal behaviour of thin metal films described by a hyperbolic two-step model. The generalized equations of this model contains diffusion terms in both the electron and lattice energy equations and assumes that incident laser radiation is absorbed by both the electron gas and solid lattice to account for the thermal behaviour of semiconducting and impure metals. A perturbation technique is utilized to eliminate the coupling between the electron and phonon energy equations when the normalized temperature difference between electrons and phonons is a small quantity, which is true in materials exhibit high coupling factors. Received on 6 January 1999  相似文献   

13.
In this paper, the newly developed lattice Boltzmann flux solver (LBFS) is developed into a version in the rotating frame of reference for simulation of turbomachinery flows. LBFS is a finite volume solver for the solution of macroscopic governing differential equations. Unlike conventional upwind or Godunov‐type flux solvers which are constructed by considering the mathematical properties of Euler equations, it evaluates numerical fluxes at the cell interface by reconstructing local solution of lattice Boltzmann equation (LBE). In other words, the numerical fluxes are physically determined rather than by some mathematical approximation. The LBE is herein expressed in a relative frame of reference in order to correctly recover the macroscopic equations, which is also the basis of LBFS. To solve the LBE, an appropriate lattice Boltzmann model needs to be established in advance. This includes both the determinations of the discrete velocity model and its associated equilibrium distribution functions. Particularly, a simple and effective D1Q4 model is adopted, and the equilibrium distribution functions could be efficiently obtained by using the direct method. The present LBFS is validated by several inviscid and viscous test cases. The numerical results demonstrate that it could be well applied to typical and complex turbomachinery flows with favorable accuracy. It is also shown that LBFS has a delicate dissipation mechanism and is thus free of some artificial fixes, which are often needed in conventional schemes. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

14.
This paper presents a new model of lattice Boltzmann method for full compressible flows. On the basis of multi‐speed model, an extra potential energy distribution function is introduced to recover the full compressible Navier–Stokes equations with a flexible specific‐heat ratio and Prandtl number. The Chapman–Enskog expansion of the kinetic equations is performed, and the two‐dimension‐seventeen‐velocity density equilibrium distribution functions are obtained. The governing equations are discretized using the third order monotone upwind scheme for scalar conservation laws finite volume scheme. The van Albada limiter is used to avoid spurious oscillations. In order to verify the accuracy of this double‐distribution‐function model, the Riemann problems, Couette flows, and flows around a NACA0012 airfoil are simulated. It is found that the proposed lattice Boltzmann model is suitable for compressible flows, even for strong shock wave problem, which has an extremely large pressure ratio, 100,000. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

15.
A possible improvement of a continuum model for diatomic crystals is examined using continuum limit of a discrete diatomic model. For this purpose, various discrete models of diatomic lattice are compared at the linearized and weakly nonlinear levels. The suitable numbering of the atoms in the lattice is found which is better adopted for continualization than the familiar pair numbering introducing two sub-lattices. The coupled governing partial nonlinear differential equations for longitudinal strain and relative distance between the atoms are obtained in the continuum limit that allows us to describe localization of the strains due to the presence of the atoms of two kinds. It is found, that the equations obtained possess two kinds of localized wave solutions, one related to the acoustical branch and the other one related to the optical branch.  相似文献   

16.
A coupled numerical method for the direct numerical simulation of particle–fluid systems is formulated and implemented, resolving an order of magnitude smaller than particle size. The particle motion is described by the time-driven hard-sphere model, while the hydrodynamic equations governing fluid flow are solved by the lattice Boltzmann method (LBM). Particle–fluid coupling is realized by an immersed boundary method (IBM), which considers the effect of boundary on surrounding fluid as a restoring force added to the governing equations of the fluid. The proposed scheme is validated in the classical flow-around-cylinder simulations, and preliminary application of this scheme to fluidization is reported, demonstrating it to be a promising computational strategy for better understanding complex behavior in particle–fluid systems.  相似文献   

17.
In this paper we derive a macroscopic kinetic law for twin boundary motion from a lattice dynamical model. The model is developed for compound and type-1 twins and it is explicitly illustrated for a Cu-Al-Ni shape memory alloy. The governing multiple-well energy is calculated using an effective interatomic potential; a Frenkel-Kontorowa type model is developed for the dynamics at the lattice scale; and a quasi-continuum approximation is used to determine the continuum-scale kinetics. The model predicts that compound twins in the Cu-Al-Ni shape memory alloy are an order of magnitude more mobile than type-1 twins which is consistent with experimental observations.  相似文献   

18.
Lattice models with long-range interactions of power-law type are suggested as a new type of microscopic model for fractional non-local elasticity. Using the transform operation, we map the lattice equations into continuum equation with Riesz derivatives of non-integer orders. The continuum equations that are obtained from the lattice model describe fractional generalization of non-local elasticity models. Particular solutions and correspondent asymptotic of the fractional differential equations for displacement fields are suggested for the static case.  相似文献   

19.
基于电子-声子相互作用的双曲两步热传导模型的超快热弹性理论,计及晶格的热传导效应,利用有限元方法研究了无限大金属薄膜在短暂激光冲击下诱导的位移、应力、应变和温度等物理量的演化特点,与已有文献比较,说明该方法的合理性与正确性.比较了计及电子热爆发力与不计电子热爆发力对位移、应力等物理量的影响,说明计及电子热爆发力的必要性...  相似文献   

20.
Analytical solution for electron and lattice temperature distribution in the solid initially heated by a laser short-pulse is presented. Strained parameters method is introduced when formulating electron and lattice temperature distributions. Laser short pulse heating of gold film is simulated numerically and temperature data at the end of the heating pulse are adopted as initial condition to the governing equations of energy transport for analytical solutions. This enables to solve the governing equations of energy analytically in the cooling period. It is found that electron temperature decays sharply while lattice site temperature increases gradually in the surface regions during the cooling cycle. As the depth from the surface increases change in both temperatures become gradual.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号