Alternative boundary conditions for solving optical diffusion equations by a finite difference time domain analysis in three-dimensional scattering medium

Alternative boundary conditions for solving optical diffusion equations in three-dimensional (3-D) scattering medium by a finite difference time domain (FDTD) analysis formulated by the author are proposed. The previous boundary conditions were defined only by fluence rate, which, although essential, is only one factor needed to solve approximated diffusion equations for fluence rate. In this paper, alternative boundary conditions defined both by fluence rate and radiant flux have been proposed for use in the FDTD analysis, which is derived from the two coupled differential equations for fluence rate and radiant flux. It has been become clear that these boundary conditions are almost equivalent to the previous boundary conditions in the FDTD analysis for sufficiently fine grid spacing. For the analysis with coarser grid spacing, the proposed boundary conditions suppress analytical errors, especially in intensity of time-resolved reflectance and transmittance.     

等离子体粒子模拟中的有效电流分配方法

介绍Villasenor和Buneman提出的一种可以自动满足静电修正要求的有效电流分配方法.着重讨论在实际模拟中对出界粒子施加粒子边界条件时这种方法有效性的破坏,针对几种典型的粒子边界条件,提出相应的修正方案.模拟实例表明,采用此方法及修正方案后,得到的电磁场能自动地满足Poisson方程.     

耗散粒子动力学处理复杂固体壁面的一种有效方法

耗散粒子动力学(dissipative particle dynamics,DPD)作为一种介观尺度拉格朗日型粒子方法,已经成功地应用于微纳米流动和生化科技的研究中.复杂固体壁面的处理和壁面边界条件的实施一直是DPD方法发展及应用的一个障碍.提出了处理复杂固体壁面的一种新的方法.复杂固体区域通过冻结随机分布并且达到平衡状态的DPD粒子代表;所冻结的DPD粒子位于临近流动区域的一个截距内;在靠近固体壁面的流动区域中设置流动反弹层,当流动DPD粒子进入此流动层后反弹回流动区域.应用这种固体壁面处理方法对简单流动区域的Poiseuille流动和复杂多孔介质内的流动进行了分析.研究表明,这种新的固体壁面处理方法能够有效模拟复杂固体区域,准确实施壁面边界条件.     

边界扫描测试算法的分析与优化

伴随着电子技术和边界扫描测试技术的飞速发展,新的边界扫描测试算法也在不断涌现;而边界扫描测试的算法,一般都是指在互联模型的基础上,边界扫描测试向量的生成算法;生成合理的测试向量集可以以最短的测试时间来覆盖尽可能多的故障;从对一些常用的边界扫描测试算法进行了粗略的分析,到后来对等权值算法和二进制计数算法进行了详细的分析,通过引入和分析边界扫描测试算法的定理、公式以及推论等,分别提出了等权值优化算法和权值递加算法;与优化前的算法作为比较,等权值优化算法降低了征兆混淆出现的概率,而权值递加算法同时降低了征兆误判率和征兆混淆率;综合分析,新的算法更好的权衡了测试向量集的完备性指标和紧凑型指标。     

Implicit spectrally-accurate method for moving boundary problems using immersed boundary conditions concept

A fully implicit, spectral algorithm for the analysis of moving boundary problem is described. The algorithm is based on the concept of immersed boundary conditions (IBC), i.e., the computational domain is fixed while the time dependent physical domain is submerged inside the computational domain, and is described in the context of the diffusion-type problems. The physical conditions along the edges of the physical domain are treated as internal constraints. The method eliminates the need for adaptive grid generation that follows evolution of the physical domain and provides sharp resolution of the location of the boundary. Various tests confirm the spectral accuracy in space and the first- and second-order accuracy in time. The computational cost advantage of the IBC method as compared with the more traditional algorithm based on the mapping concept is demonstrated.     

A meshless algorithm with moving least square approximations for elliptic Signorini problems

Based on the moving least square(MLS) approximations and the boundary integral equations(BIEs), a meshless algorithm is presented in this paper for elliptic Signorini problems. In the algorithm, a projection operator is used to tackle the nonlinear boundary inequality conditions. The Signorini problem is then reformulated as BIEs and the unknown boundary variables are approximated by the MLS approximations. Accordingly, only a nodal data structure on the boundary of a domain is required. The convergence of the algorithm is proven. Numerical examples are given to show the high convergence rate and high computational efficiency of the presented algorithm.     

光纤模式分析的伽辽金有限元模型

提出了采用二阶吸收边界条件的全矢量平面伽辽金有限元模型,用于分析任意横截面形状和各向异性折射率分布的光纤的传导模式和泄漏模式,能精确求出各模式传输常量的实部和虚部以及模场分布,既不出现伪解,又不漏解.推导了各向异性介质全矢量耦合波动方程及其变分形式,给出了基于结点的二阶三角形单元的离散公式和单元矩阵,成功的将二阶吸收边界条件加入外边界二次线性单元的离散公式.计算表明采用该模型分析光子晶体光纤模式有效折射率与采用多极方法和基于离散函数展开的有限差分法所得结果吻合很好,采用二阶吸收边界条件计算限制损耗比一阶吸收边界条件结果精确.     

共轭梯度法求解双曲传热多宗量反问题

提出双曲传热反问题热物性参数和边界条件多宗量联合反演的一般数值求解模式,考虑了非均质和分布参数的影响,时域上采用时域精细算法进行离散,建立了便于敏度分析的有限元正演模型.由最小二乘原理建立反演模型,应用共轭梯度法进行求解.探讨了时间步长和测量误差对反演结果的影响,并进行了数值验证.     

Hybrid power flow analysis using coupling loss factor of SEA for low-damping system—Part I: Formulation of 1-D and 2-D cases

This paper proposes a hybrid method for the prediction of vibrational and acoustic responses of low-damping system in the medium-to-high frequency ranges by using the power flow analysis (PFA) algorithm and statistical energy analysis (SEA) coupling concepts. The main part of this method is the application of the coupling loss factor (CLF) of SEA to the boundary condition of PFA in reverberant system. First, for hybrid PFA, the hybrid boundary conditions on 1-D and 2-D cases were derived in the general form. To verify the derived boundary conditions, numerical analyses for each case were performed. The hybrid PFA solutions using derived boundary conditions were compared with the classical PFA solutions with various reverberance factors including the effects of the characteristic length, excitation frequency and group velocity besides damping loss factor of the subsystem. Additionally, the hybrid PFA on 3-D case and the hybrid power flow finite element method (PFFEM) for hybrid PFA of built-up structures are described in the other companion paper.     

DSMC方法下给定热流边界条件的实施

本文详细讨论了DSMC方法中流体温度、能量及边界热流的统计方法,发展了一种从边界热流求得与壁面碰撞分子反射速度的方法。该方法被称为逆温度抽样算法(ITS,Inverse Temperature Sampling)方法。在此基础上,本文发展了DSMC方法中壁面处给定热流边界条件的实施方法。计算结果表明: ITS方法能准确抽样反射分子的特征温度,进而求得分子反射速度。基于该方法的给定热流边界条件可以准确求得壁面处温度分布,以及流场内的压力、速度。     

Heat transfer of non-Newtonian liquid flows in convergent-divergent channels with fourth-kind thermal boundary conditions

In the present paper, we consider the heat transfer of non-Newtonian liquid flows moving in convergent- divergent channels with common boundary. The case of flows moving in opposite directions was addressed. At the common boundary, continuity conditions for temperature and heat flow were adopted. A mathematical model, a calculation algorithm, and simulation data are reported. Plots of numerical data characterizing the channel flows are presented. A comparative analysis of heat transfer in smooth and convergent-divergent channels is given.     

An explicit algorithm for imbedding solid boundaries in Cartesian grids for the reactive Euler equations

We present a local and point-wise scheme for imposing reflective boundary conditions to stationary internal boundaries for solving the reactive Euler equations on Cartesian grids. The scheme is presented in two and three dimensions and can run efficiently on parallel machines while still maintaining the same advantages over other methods for enforcing internal boundary conditions. Level sets are used to represent internal solid regions along with a new local node sorting algorithm that decouples internal boundary nodes by establishing their connectivity to other internal boundary nodes. This approach allows us to enforce boundary conditions via a direct procedure, removing the need to solve a coupled system of equations numerically. We examine the accuracy and fidelity of our internal boundary algorithm by simulating flows past various solid boundaries in two and three dimensions, showing good agreement between our numerical results and experimental data.     

Generalized analysis of a muffler with any number of interacting ducts

The muffler elements that use perforated elements for acoustic attenuation are common in practice. In typical commercial mufflers perforated elements are used involving two, three, four or more interacting ducts. Analysis of such configurations involves writing down the basic governing equations of mass continuity, momentum balance, etc., and then elimination of velocity variables to obtain the coupled ordinary differential equations in terms of acoustic pressure variables. Mathematical modelling and the consequent analytical derivation of the transmission loss for these multi-duct acoustical elements become increasingly tedious, as just not the number of ducts, but also their relative arrangement along with the boundary conditions dictate the analysis considerably. In the present paper, authors have proposed a generalization and thus an algebraic algorithm to directly produce the system matrix, eliminating the tedium of writing the basic governing equations and elimination of velocity variables. Also, a convenient approach for applying the boundary conditions is outlined here.     

Exact solution for the fractional cable equation with nonlocal boundary conditions

The fractional cable equation is studied on a bounded space domain. One of the prescribed boundary conditions is of Dirichlet type, the other is of a general form, which includes the case of nonlocal boundary conditions. In real problems nonlocal boundary conditions are prescribed when the data on the boundary can not be measured directly. We apply spectral projection operators to convert the problem to a system of integral equations in any generalized eigenspace. In this way we prove uniqueness of the solution and give an algorithm for constructing the solution in the form of an expansion in terms of the generalized eigenfunctions and three-parameter Mittag-Leffler functions. Explicit representation of the solution is given for the case of double eigenvalues. We consider some examples and as a particular case we recover a recent result. The asymptotic behavior of the solution is also studied.     

基于改进Mallat算法的倍频程分析方法研究

倍频程分析是噪声测量系统中的重要环节,针对传统倍频程分析方法中滤波器相对衰减值测量精度不高和声压级谱测量不准确的缺陷,本文提出了一种基于改进Mallat算法的倍频程分析方法。本文提出的改进Mallat算法:第一,对信号进行FFT运算,删除由于滤波器的非锐截止特性导致的频率混淆分量;第二,利用矩形窗函数与信号相乘,消除由于卷积运算时会导致的边界效应;最后,应用虚拟仪器软件开发平台Lab VIEW搭建了基于改进Mallat算法的倍频程分析系统。仿真对比实验结果验证了改进Mallat倍频程分析方法对噪声信号的相对衰减值测量精度更高,声压级谱显示更加准确。     

Exact boundary conditions for the initial value problem of convex conservation laws

The initial value problem of convex conservation laws, which includes the famous Burgers’ (inviscid) equation, plays an important rule not only in theoretical analysis for conservation laws, but also in numerical computations for various numerical methods. For example, the initial value problem of the Burgers’ equation is one of the most popular benchmarks in testing various numerical methods. But in all the numerical tests the initial data have to be assumed that they are either periodic or having a compact support, so that periodic boundary conditions at the periodic boundaries or two constant boundary conditions at two far apart spatial artificial boundaries can be used in practical computations. In this paper for the initial value problem with any initial data we propose exact boundary conditions at two spatial artificial boundaries, which contain a finite computational domain, by using the Lax’s exact formulas for the convex conservation laws. The well-posedness of the initial-boundary problem is discussed and the finite difference schemes applied to the artificial boundary problems are described. Numerical tests with the proposed artificial boundary conditions are carried out by using the Lax–Friedrichs monotone difference schemes.     

Algorithm for particle packing based on Aristotelian mechanics

A dynamic algorithm is proposed for three-dimensional packing of spherical solid particles. The particles are deposited within a specified region with a fixed rigid boundary. The velocity of each particle is proportional to its weight and forces due to contact of the particle with the boundary and neighbor particles. Dimensional analysis of the equations of particle motion is performed. The average density and coordination number distribution for an equilibrium packing are calculated. The dependence of these characteristics on viscosity, granulometric composition, and representation of initial conditions (numerical analogue of material pouring into a specified volume) is studied.     

Modelling transient heat conduction in solids at multiple length and time scales: A coupled non-equilibrium molecular dynamics/continuum approach

A method for controlling the thermal boundary conditions of non-equilibrium molecular dynamics simulations is presented. The method is simple to implement into a conventional molecular dynamics code and independent of the atomistic model employed. It works by regulating the temperature in a thermostatted boundary region by feedback control to achieve the desired temperature at the edge of an inner region where the true atomistic dynamics are retained. This is necessary to avoid intrinsic boundary effects in non-equilibrium molecular dynamics simulations. Three thermostats are investigated: the global deterministic Nosé–Hoover thermostat and two local stochastic thermostats, Langevin and stadium damping. The latter thermostat is introduced to avoid the adverse reflection of phonons that occurs at an abrupt interface. The method is then extended to allow atomistic/continuum models to be thermally coupled concurrently for the analysis of large steady state and transient heat conduction problems. The effectiveness of the algorithm is demonstrated for the example of heat flow down a three-dimensional atomistic rod of uniform cross-section subjected to a variety of boundary conditions.     

Binary scattering and breakup in the three-nucleon system

We present the further development of the three-particle formalism for differential Faddeev equations. The asymptotic boundary conditions in the hyperspherical adiabatic representation have been constructed. We prove that these conditions are asymptotically equivalent to the standard Merkuriev boundary conditions. With these boundary conditions we have formulated the boundary-value problem for Faddeev equations which has the property that the binary channel and the breakup channel are explicitly orthogonal. The effective numerical scheme for solving the formulated boundary-value problem is given.