Modelling of a single drop impact onto liquid film using particle method

The process of single liquid drop impact on thin liquid surface is numerically simulated with moving particle semi‐implicit method. The mathematical model involves gravity, viscosity and surface tension. The model is validated by the simulation of the experimental cases. It is found that the dynamic processes after impact are sensitive to the liquid pool depth and the initial drop velocity. In the cases that the initial drop velocity is low, the drop will be merged with the liquid pool and no big splash is seen. If the initial drop velocity is high enough, the dynamic process depends on the liquid depth. If the liquid film is very thin, a bowl‐shaped thin crown is formed immediately after the impact. The total crown subsequently expands outward and breaks into many tiny droplets. When the thickness of the liquid film increases, the direction of the liquid crown becomes normal to the surface and the crown propagates outward. It is also found that the radius of the crown is described by a square function of time: r_C = [c(t ? t₀)]^0.5. When the liquid film is thick enough, a crown and a deep cavity inside it are formed shortly after the impact. The bottom of the cavity is initially oblate and then the base grows downward to form a sharp corner and subsequently the corner moves downward. Copyright © 2004 John Wiley & Sons, Ltd.     

亚像元动态成像技术中系统调制传递函数的数值分析

亚像元动态成像技术是目前实现遥感器小型化非常有效的新方法,因而研究亚像元动态成像系统的调制传递函数MTF显得非常必要和重要.介绍了亚像元动态成像的基本原理及成像过程,指出了光程差对亚像元动态成像系统MTF的影响.对亚像元动态成像系统的MTF进行了数值分析,并指出了亚像元动态成像系统的成像质量及空间分辨率与一般成像系统相比较的优点.     

A relaxed model and its homogenization for nematic liquid crystals in composite materials

We analyse a model for equilibrium configurations of composite systems of nematic liquid crystal with polymer inclusions, in the presence of an external magnetic field. We assume that the system has a periodic structure, and consider the relaxed problem on the unit length constraint of the nematic director field. The relaxation of the Oseen–Frank energy functional is carried out by including bulk as well as surface energy penalty terms, rendering the problem fully non‐linear. We employ two‐scale convergence methods to obtain effective configurations of the system, as the size of the polymeric inclusions tends to zero. We discuss the minimizers of the effective energies for, both, the constrained as well as the unconstrained models. Copyright © 2004 John Wiley & Sons, Ltd.     

Modelling of an underground waste disposal site by upscaling

In this paper, we study the global behaviour of an underground waste disposal in order to have an accurate upscaled model suitable for the computations involved in safety assessment processes. We start from a detailed model describing the transport of pollutant leaking from a high number of units. Using the method of homogenization, going to the limit, we obtain first a macroscopic model where the sources are now appearing globally. Then we compute a first‐order matched asymptotic expansion and we give the error estimates for this approximation. Copyright © 2004 John Wiley & Sons, Ltd.     

The long‐time behaviour of the thermoconvective flow in a porous medium

For the Boussinesq approximation of the equations of coupled heat and fluid flow in a porous medium we show that the corresponding system of partial differential equations possesses a global attractor. We give lower and upper bounds of the Hausdorff dimension of the attractor depending on a physical parameter of the system, namely the Rayleigh number of the flow. Numerical experiments confirm the theoretical findings and raise new questions on the structure of the solutions of the system. Copyright © 2004 John Wiley & Sons, Ltd.     

晶体微观结构的数值计算

晶体微观结构是晶体材料在特定物理条件下其多个能量极小平衔态在空间形成的某种微尺度的规则分布．几何非线性的连续介质力学理论可以用能量极小化原理来解释晶体微观结构的形成，并用Young测度来刻画平衡态各变体在空间的概率分布．定性的理解与定量地分析和计算晶体材料的微观结构对于发展和改进高级晶体功能材料，如形状记忆合金、铁电体、磁至伸缩材料等，有重要的意义．本文回顾了近年来晶体微观结构数值计算方面的最新进展．介绍了计算晶体微观结构的几种数值方法及有关的数值分析结果。     

Continuation of travelling‐wave solutions of the Navier–Stokes equations

An efficient way of obtaining travelling waves in a periodic fluid system is described and tested. We search for steady states in a reference frame travelling at the wave phase velocity using a first‐order pseudospectral semi‐implicit time scheme adapted to carry out the Newton's iterations. The method is compared to a standard Newton–Raphson solver and is shown to be highly efficient in performing this task, even when high‐resolution grids are used. This method is well suited to three‐dimensional calculations in cylindrical or spherical geometries. Copyright © 2006 John Wiley & Sons, Ltd.     

ΓD‐convergence for a class of quasilinear elliptic equations in thin structures

We study the asymptotic behaviour of the solution of a stationary quasilinear elliptic problem posed in a domain Ω^(ε) of asymptotically degenerating measure, i.e. meas Ω^(ε) → 0 as ε → 0, where ε is the parameter that characterizes the scale of the microstructure. We obtain the convergence of the solution and the homogenized model of the problem is constructed using the notion of convergence in domains of degenerating measure. Proofs are given using the method of local characteristics of the medium Ω^(ε) associated with our problem in a variational form. Copyright © 2005 John Wiley & Sons, Ltd.     

叶栅全三维粘性反问题的数值解

本文发展了一种解叶栅全三维粘性反问题的新的数值方法.基于非正交曲线坐标与相应的非正交速度分量下完全守恒型的Navier－Stokes方程,全三维反问题规定叶片表面的无量纲压力分布反求叶型。计算中叶片表面的边界条件采用一种特殊的方式来处理，即一方面强加给定的压力分布条件,另方面叶面的几何位置在迭代过程中又是可移动的，其移动速度将与Navier—Stokes方程在当地的解联系起来，从而形成一种解定常问题的新的不定常过程.试算证明了本文方法的可行性。     

On parabolic equations with a small parameter

A singularly perturbed parabolic equation with a nonlinear right-hand side of a special form is examined. A numerical analytical study of such equations is performed.