中低纬电离层电场理论模式

介绍了一个电离层电场理论模式. 该模式从电离层发电机理论的基本方程出发, 采用地磁偶极坐标系, 推导出电离层电势满足的微分方程作为电场理论模式的出发方程. 模式的主要输入参量为电离层背景中性风(由经验模式HWM93给出)和电导率(由经验模式IRI90, MSISE90分别给出的电子和离子的浓度、温度以及中性大气的浓度、温度计算得出). 采用松弛迭代法求解电势的偏微分方程得出电势, 进而获得中低纬电离层电场、电流随时间、高度和地磁纬度的分布. 模式很好地再现了电离层电势、电场和赤道电急流的基本结构和形态, 可应用于对高层大气和电离层的电动力学过程的研究.     

一类植物病虫害模型及数值模拟

通过考虑某一区域内的植物在喷洒杀虫剂之后分成的三种类别,建立了一个植物病虫害模型.然后给出它的正平衡点,又利用稳定性第一近似方法讨论正平衡点的稳定性.最后,结合千岛湖发生的松毛虫害使用MATLAB进行数据模拟,多角度比较不同的移除系数和树木间的有效接触率分析其对虫害发展的影响,对病虫害的治理具有较好的指导作用.     

等离子体反应器流场的数值模拟

本文采用颗粒轨道模型对等离子体反应器进行数值模拟.数值格式用SIMPLER法.湍流模型用简单的次网格模型.计算结果给出速度场、温度场及颗粒轨道,对工程设计和工艺过程参数控制有参考价值.     

一种基于积分方程的数值微分方法

本文研究了目前一些求解数值微分的方法无法求出端点导数或是求出的端点附近导数不可用的问题.利用构造一类积分方程的方法,将数值微分问题转化为这类积分方程的求解,并用一种加速的迭代正则化方法来求解积分方程. 数值实验结果表明该算法可以有效求出端点的导数,且具有数值稳定、计算简单等优点.     

多种环境外力作用下作物生长系统的动力学模型及过程数值模拟

依据生物力学的理论与方法,综合探讨了多种环境因子(光照、温度、土壤水分与养分等)外力作用下作物生长的动力学机制.采用连续马尔柯夫过程方法建立了作物生长系统的动力学模型及数值模拟方法,研究了生长速率对外力条件变化的反应与生物量饱和值的变化,并以半干旱区作物为例进行了数值模拟计算与分析,为大田管理提供了定量依据.该动力模型与其它植物生长模型相比,其优越性在于从宏观和微观相结合上反映了资源利用的多维性,揭示了资源转化过程.所建作物生长系统模拟方法,具有先进性和可操作性,实际模拟结果与大田观测结果能很好吻合.     

非常规井数值模拟技术研究进展与发展趋势

与常规井数值模拟相比,非常规井数值模拟侧重于刻画垂直、水平井筒内的多相流动,解释由于摩擦、静水力、加速度以及滑移造成的各相之间的速度差及压力损失,模拟井下流量控制设备及复杂现象.从复杂结构井的流量方程、多相井筒管流与油藏渗流耦合模拟、多段井模型、井下流量控制设备模拟、近井区域粗化、井下复杂现象模拟和嵌入式裂缝模型7个方面出发总结了非常规井数值模拟技术研究进展和主要缺陷.基于多段井模型的井筒 油藏耦合数值模拟、扩展井模型与油藏模型耦合模拟、流固耦合数值模拟、对井下监测调控设备的数值模拟以及多段压裂复杂结构井的数值模拟是未来的发展趋势.     

一类边界耦合的扩散系统爆破现象的数值模拟

对一类边界上非平凡耦合的抛物型系统进行了数值模拟．为验证已有的关于是否出现爆破的理论成果,先用固定网格算法针对四种具体的边界条件进行试验,并根据不同的初始数据分组．为了进一步探讨爆破发生的时刻、位置以及爆破速率,再用移动网格算法针对可能出现爆破现象的两组边界条件和初始数据进行试验,并根据不同的监测函数分组．随后对算法的有效性做出说明并分析试验结果．最后对系统的一种特殊情况给出一个算例．     

昆虫振翅飞行的数值模拟

在非惯性参考系下对昆虫振翅扰动的二维非定常流场进行了数值模拟,避免了计算中的移动边界困难,从而缩短了计算时间,模型具有3个自由度,可以模拟任意已知的翅的平面运动．通过模拟相对复杂的自然界昆虫的振翅运动,研究昆虫是如何控制飞行．计算结果表明,有2个参数可能被昆虫用来控制飞行:翅平动和转动间的相位差以及垂直于平均振翅平面方向的横向振幅．     

基于迭代正则化高斯-牛顿法的非线性Urysohn积分方程数值解

非线性Urysohn积分方程在许多领域中都有广泛的应用,但由于该方程具有不适定性的特点,数据的微小扰动可能导致解的巨大变化,给数值求解带来很大困难.为了获得稳定的、准确的数值解,本文利用迭代正则化高斯-牛顿法对此方程进行求解,给出了利用Sigmoid-型函数确定迭代正则化参数的方法.对一类重力测定问题进行了数值模拟,将得到的数值解和相应的精确解作比较.结果表明,本文提出的方法在求解非线性Urysohn积分方程时是可行的也是有效的.     

数值模拟混溶驱动问题的迎风格式及理论分析

本文提出数据模拟混溶驱动问题的一类广义迎风格式，严格证明了其满足离散极值原理和离散质量守恒原理，并获得收敛性定理，对模型问题的试算，得到令人满意的数值结果。     

Theoretical and numerical analysis on multispectral bioluminescence tomography

Recently, molecular imaging has been rapidly developed to studyphysiological and pathological processes in vivo at the cellularand molecular levels. Among molecular imaging modalities, opticalimaging has attracted a major attention for its unique advantages.In this paper, we establish a mathematical framework for multispectralbioluminescence tomography (BLT) that allows simultaneous studiesof multiple optical reporters. We show solution existence, uniquenessand continuous dependence on data as well as the limiting behaviourswhen the regularization parameter approaches zero or when thepenalty parameter approaches infinity. Then, we propose twonumerical schemes for multispectral BLT and derive error estimatesfor the corresponding solutions.     

Temporally regularized direct numerical simulation

Experience with fluid-flow simulation suggests that, in some instances, under-resolved direct numerical simulation (DNS), without a residual-stress model per se but with artificial damping of small scales to account for energy lost in the cascade from resolved to unresolved scales, may be as reliable as simulations based on more complex models of turbulence. One efficient and versatile manner to selectively damp under-resolved spatial scales is by a relaxation regularization, e.g. Stolz and Adams [S. Stolz, N.A. Adams, An approximate deconvolution procedure for large eddy simulation, Phys. Fluids II (1999) 1699-1701]. We consider the analogous approach based on time scales, time filtering and damping of under-resolved temporal features. The paper explores theoretical and practical aspects of temporally damped fluid-flow simulations. We prove existence of solutions to the resulting continuum model. We also establish the effect of the damping of under-resolved temporal features as the energy balance and dissipation and prove that the time fluctuations → 0 in a precise sense. The method is then demonstrated to obtain both steady-state and time-dependent coarse-grid solutions of the Navier-Stokes equations.     

Theoretical study and numerical simulation of textiles

We propose a new approach for developing continuum models fit to describe the mechanical behavior of textiles. We develop a physically motivated model, based on the properties of the yarns, which can predict and simulate the textile behavior. The approach relies on the selection of a suitable topological model for the patch of the textile, coupled with constitutive models for the yarn behavior. The textile structural configuration is related to the deformation through an energy functional, which depends on both the macroscopic deformation and the distribution of internal nodes. We determine the equilibrium positions of these latter, constrained to an assigned macroscopic deformation. As a result, we derive a macroscopic strain energy function, which reflects the possibly nonlinear character of the yarns as well as the anisotropy induced by the microscopic topological pattern. By means of both analytical estimates and numerical experiments, we show that our model is well suited for both academic test cases and real industrial textiles, with particular emphasis on the tricot textile.     

On numerical simulation of tsunami run-up

The improved adaptive grid method is presented for the numerical modelling of propagation of long surface waves in large water areas and their interaction with coasts. (© 2015 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Turbulent diffusion: Direct numerical simulation

Direct numerical simulations (DNS) of a vibrating grid turbulence at two different Reynolds numbers (based on an amplitude and frequency of the grid) Re = 500, 1000 are presented. The evolution of a turbulent/non-turbulent interface was detected and compared to the theoretical results obtained in the paper [1] using Lie group analysis of the governing equations. Good agreement was found between DNS and theory. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Direct numerical simulation of compressible isotropic turbulence

Direct numerical simulation (DNS) of decaying compressible isotropic turbulence at turbulence Mach numbers of M_t = 0.2-0.7 and Taylor Reynolds numbers of 72 and 153 is performed by using the 7th order upwind-biased difference and 8th order center difference schemes. Results show that proper upwind-biased difference schemes can release the limit of“start-up” problem to Mach numbers. Compressibility effects on the statistics of turbulent flow as well as the mechanics of shocklets in compressible turbulence are also studied, and the conclusion is drawn that high Mach number leads to more dissipation. Scaling laws in compressible turbulence are also analyzed. Evidence is obtained that scaling laws and extended self similarity (ESS) hold in the compressible turbulent flow in spite of the presence of shocklets, and compressibility has little effect on scaling exponents.     

稠油开采数值模拟研究

针对稠油油藏使用水平井并且注蒸汽的开采方式,建立了相应的数学模型,模型考虑了油、气、水三个相态及所有组分,和由于稠油流体特性而引起的启动压力,对数学模型使用九点差分进行离散,用全隐式方法将数值模型线性化,给出了模型等效化的处理方法,将计算结果与模拟区块的实际生产数据做动态历史拟合,并提出了开发调整方案,取得了良好效果.该模型与其它方法做了比较,模型所得结果与实际产量的绝对值偏差为1.9％,是所有方法中最小的.模型与求解方法适用于稠油油藏注蒸汽水平并开采数值模拟.     

Dynamical analysis and numerical simulation of bloch chaotic system

In this article, some dynamics of Bloch chaotic system have been studied. Based on Lagrange multiplier method, optimization theory, and the generalized positively definite and radially unbound Lyapunov functions with respect to the parameters of the system, we derive the ultimate bound and a family of mathematical expressions of globally exponentially attractive sets for this system with respect to the parameters of system. The results obtained in this article provides theory basis for chaotic synchronization, chaotic control, Hausdorff dimension, and Lyapunov dimension of chaotic attractors of Bloch chaotic system. © 2016 Wiley Periodicals, Inc. Complexity 21: 201–206, 2016     

Mathematical and numerical aspects of one-dimensional laminar flame simulation

The aim of this paper is to solve several mathematical and numerical questions related to the simulation of stationary and nonstationary premixed flat flames. Most of the results are obtained in the general context of complex chemical and diffusion mechanisms. The main mathematical results concern: (i) thea priori positivity of the mass fractions, and (ii) the sensitivity of the flame speed to the computational domain. The numerical method proposed for solving the stationary problem is a new combination of the pseudo-nonstationary approach, the Newton iterations, and the adaptive gridding. The computation of H₂-O₂-N₂ flames with various initial concentrations (including the chemical extinction zone) shows the efficiency of this method.