高波数Helmholtz 方程的内罚有限元方法

本文考虑二维和三维区域上高波数Helmholtz 散射问题的线性内罚有限元方法. 该散射问题的边界条件取为一阶吸收边界条件. 本文证明了, 如果加罚参数γ-γ_r+iγ_i 的虚部 γ_i 大于零, 那么内罚有限元方法是绝对稳定的, 即对任意k,h,R > 0 都存在唯一解. 这里k 是波数, h 为网格尺寸, R是区域的直径. 进一步地, 如果|γ_r|≤γ_i≤1, 那么存在与k,h,γ,R 无关的常数C₀;C₁;C₂, 使得当k³h²R ≤ C₀ 时, 该方法的H¹ 误差界为(C₁kh + C₂k³h²R)RM(f, g), 当k³h²R > C₀ 且kh 有界时,H¹ 误差界为(C₁kh + C₂/γ_i)RM(f, g), 其中M(f, g) := (‖f‖_L²(Ω) + R^-1/2‖g‖_L²(Γ)) + R^-1|g|_H^1/2(Γ). 另外, 本文还推导了L² 误差估计. 注意到γ = 0 时内罚有限元方法就是经典的有限元方法, 通过取加罚参数为iγ_>i 并令γ_i 趋于0+, 本文还在k³h²R ≤ C₀ 的条件下, 得到了有限元方法的稳定性和误差估计.作者以前的工作只考虑了加罚参数为纯虚数的情形并且没有考虑对R 的依赖关系.     

Mathematical modelling of the flotation deinking process

The overall flotation deinking process can be divided into four basic microprocesses:

• 1.(1) collision or capture of an (ink) particle by an air bubble
• 2.(2) adhesion of an (ink) particle to the air bubble by sliding
• 3.(3) development of a three-phase contact at the air bubble/water/particle interface, and
• 4.(4) bubble/particle stability or instability after an aggregate is formed each of these microprocesses have an associated probability that they will occur successfully in a flotation cell.

In this paper, the associated probabilities of each microprocess are employed in the development of a kinetic- or population balance-type model of the overall flotation process. The overall model contains two kinetic constants: the first, k₁ governs the overall probability of a free ink particle successfully intercepting and adhering to an air bubble; the second, k₂ is a measure of the probability that a bubble/particle aggregate pair will become unstable and split to yield a “new” free ink particle.The solution to the kinetic model is presented in terms of k₁ and k₂, which are themselves functions of system parameters such as bubble and particle physical properties (e.g., diameter, density), fluid properties (e.g., viscosity, surface tension), etc. From this solution, a definition of a theoretical flotation efficiency, as well as other system performance parameters are presented.     

Modeling 2D transient heat conduction problems by the numerical manifold method on Wachspress polygonal elements

Due to the use of dual cover systems, i.e., the mathematical cover and the physical cover, the numerical manifold method (NMM) is able to solve physical problems with boundary-inconsistent meshes. Meanwhile, n-gons (n>4) are very impressive, due to their greater flexibility in discretization, less sensitivity to volumetric and shear locking, and better suitability for complex microstructures simulation. In this paper, the NMM, combined with Wachspress-type hexagonal elements, is developed to solve 2D transient heat conduction problems. Based on the governing equations, the NMM temperature approximation and the modified variational principle, the NMM discrete formulations are deduced. The solution strategy to time-dependent global equations and the spatial integration scheme are presented. The advantages of the proposed approach in both discretization and accuracy are demonstrated through several typical examples with increasing complexity. The extension of polygonal elements in unsteady thermal analysis within the NMM is realized.     

Mathematical Modeling of Czochralski Type Growth Processes for Semiconductor Bulk Single Crystals

This paper deals with the mathematical modeling and simulation of crystal growth processes by the so-called Czochralski method and related methods, which are important industrial processes to grow large bulk single crystals of semiconductor materials such as, e. g., silicon (Si) or gallium arsenide (GaAs) from the melt. In particular, we investigate a recently developed technology in which traveling magnetic fields are applied in order to control the behavior of the turbulent melt flow. Since numerous different physical effects like electromagnetic fields, turbulent melt flows, high temperatures, heat transfer via radiation, etc., play an important role in the process, the corresponding mathematical model leads to an extremely difficult system of initial-boundary value problems for nonlinearly coupled partial differential equations. In this paper, we describe a mathematical model that is under use for the simulation of real-life growth scenarios, and we give an overview of mathematical results and numerical simulations that have been obtained for it in recent years.     

Environmental Effects on the Mechanical Properties of Glass-FRP and Aramid-FRP Rebars

In the paper, the experimental results on the effect of temperature and moisture on the mechanical properties of FRP (Fiber-Reinforced Polymer) reinforcements are presented. FRP rebars made from glass and aramid fibers were subjected to cyclic thermal actions at temperatures ranging between 20 and 70°C, typical of natural hot-climate environments. Tensile tests were also carried out on FRP rebars. The effect of moisture was investigated by cyclic wetting and drying the FRP rebars under laboratory conditions before their testing in tension. Finally, the elastic modulus and tensile strength of the FRP rebars exposed to these cyclic actions were compared with those obtained for unexposed ones, in order to evaluate the mechanical damage caused by environmental conditions.     

Axi symmetric 2D simulation and numerical heat transfer characteristics for the calibration furnace in a rectangular enclosure

This paper presents axi symmetric 2D numerical investigation of the spherical thermocouple calibration furnace in a rectangular enclosure. The focus is on the flow structure inside the Saturn (a hollow spherical cavity), external flow behavior due to annulus block heating and the surface temperature uniformity. Mesh sensitivity analysis is adopted to extract the mesh with minimum number of nodes but with fast convergent finite element solution. The basic strategy here is that temperature perturbation error at a single point instead of a single element contributed to the total perturbation error qualitatively remains the same. Agreement between numerical simulation results and the experiment results is good with a maximum temperature deviation 10 °C for the cavity temperature 400 °C. Finally, standard numerical temperature uncertainty due to variation in thermal conductivity is computed through the sensitivity coefficient using uncertainty analysis.     

二阶椭圆问题基于泡函数的简化的稳定化二阶混合有限元格式

本文研究二阶椭圆方程基于泡函数的稳定化的二阶混合有限元格式,通过消去泡函数导出一种自由度很少的简化的稳定化的二阶混合有限元格式, 误差分析表明消去泡函数的简化格式与带有泡函数的格式具有相同的精度而可以节省6N_p个自由度(其中N_p三角形剖分中的顶点数目).       

Relationships for moments of order statistics from the right-truncated generalized half logistic distribution

In this paper, we establish several recurrence relations satisfied by the single and the product moments for order statistics from the right-truncated generalized half logistic distribution. These relationships may be used in a simple recursive manner in order to compute the single and the product moments of all order statistics for all sample sizes and for any choice of the truncation parameter P. These generalize the corresponding results for the generalized half logistic distribution derived recently by Balakrishnan and Sandhu (1995, J. Statist. Comput. Simulation, 52, 385–398).Earlier went by the name R. A. Sandhu.     

Computing the moments of order statistics from nonidentically distributed phase-type random variables

In this paper, we derive a recurrence relation for the single moments of order statistics (o.s.) arising from n independent nonidentically distributed phase-type (PH) random variables (r.v.’s). This recurrence relation will enable one to compute all single moments of all o.s. in a simple recursive manner.     

A note on “Price discount based on early order commitment in a single manufacturer-multiple retailer supply chain”

In a recent paper by Xie et al. [Xie, J., Zhou, D., Wei, J.C., Zhao, X., 2010. Price discount based on early order commitment in a single manufacturer-multiple retailer supply chain. European Journal of Operational Research 200, 368-376], the authors have studied the early order commitment (EOC) strategy for a decentralized, two-level supply chain consisting of a single manufacturer and multiple retailers. They fail to provide an algorithm to determine the optimal EOC periods to minimize the total supply chain cost. This note proposes a polynomial-time algorithm to find the optimal solutions, and provides a new set of sufficient conditions under which the wholesale price discount scheme coordinates the whole supply chain.     

A C0 interior penalty discontinuous Galerkin method for fourth‐order total variation flow. I: Derivation of the method and numerical results

We consider the numerical solution of a fourth‐order total variation flow problem representing surface relaxation below the roughening temperature. Based on a regularization and scaling of the nonlinear fourth‐order parabolic equation, we perform an implicit discretization in time and a C⁰ Interior Penalty Discontinuous Galerkin (C⁰IPDG) discretization in space. The C⁰IPDG approximation can be derived from a mixed formulation involving numerical flux functions where an appropriate choice of the flux functions allows to eliminate the discrete dual variable. The fully discrete problem can be interpreted as a parameter dependent nonlinear system with the discrete time as a parameter. It is solved by a predictor corrector continuation strategy featuring an adaptive choice of the time step sizes. A documentation of numerical results is provided illustrating the performance of the C⁰IPDG method and the predictor corrector continuation strategy. The existence and uniqueness of a solution of the C⁰IPDG method will be shown in the second part of this paper.     

Direct numerical simulation of single gas bubbles in pure and contaminated liquids

Disperse gas bubbles play an important role in many industrial applications. Knowing the rising velocity, the interfacial area, or the critical size for break–up or coalescence in different systems can be crucial for the process design. Hence, knowing the fundamental behaviour of a single bubble appears mandatory for the examination of bubble swarms and for the Euler–Lagrange or Euler–Euler modelling of disperse systems. In the present work a level–set–based volume–tracking method is implemented into the CFD–code OpenFOAM to follow the free interface of a single bubble. The volume–tracking method is coupled with a transport model for surfactants on the interface, including adsorption and desorption processes. The dependency of surface tension on the local surfactant concentration on the interface is modelled by a non–linear (Langmuir) equation of state. Marangoni forces, resulting from surface tension gradients, are included. The rise of a single air bubble (i) in pure water and (ii) in the presence of surfactants of different strengths is simulated. The results show good agreement with available correlations from literature. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A finite difference approach for the initial-boundary value problem of the fractional Klein-Kramers equation in phase space

Considering the features of the fractional Klein-Kramers equation (FKKE) in phase space, only the unilateral boundary condition in position direction is needed, which is different from the bilateral boundary conditions in [Cartling B., Kinetics of activated processes from nonstationary solutions of the Fokker-Planck equation for a bistable potential, J. Chem. Phys., 1987, 87(5), 2638–2648] and [Deng W., Li C., Finite difference methods and their physical constrains for the fractional Klein-Kramers equation, Numer. Methods Partial Differential Equations, 2011, 27(6), 1561–1583]. In the paper, a finite difference scheme is constructed, where temporal fractional derivatives are approximated using L1 discretization. The advantages of the scheme are: for every temporal level it can be dealt with from one side to the other one in position direction, and for any fixed position only a tri-diagonal system of linear algebraic equations needs to be solved. The computational amount reduces compared with the ADI scheme in [Cartling B., Kinetics of activated processes from nonstationary solutions of the Fokker-Planck equation for a bistable potential, J. Chem. Phys., 1987, 87(5), 2638–2648] and the five-point scheme in [Deng W., Li C., Finite difference methods and their physical constrains for the fractional Klein-Kramers equation, Numer. Methods Partial Differential Equations, 2011, 27(6), 1561–1583]. The stability and convergence are proved and two examples are included to show the accuracy and effectiveness of the method.     

Regular Score Tests of Independence in Multivariate Extreme Values

The score tests of independence in multivariate extreme values derived by Tawn (Tawn, J.A., “Bivariate extreme value theory: models and estimation,” Biometrika 75, 397–415, 1988) and Ledford and Tawn (Ledford, A.W. and Tawn, J.A., “Statistics for near independence in multivariate extreme values,” Biometrika 83, 169–187, 1996) have non-regular properties that arise due to violations of the usual regularity conditions of maximum likelihood. Two distinct types of regularity violation are encountered in each of their likelihood frameworks: independence within the underlying model corresponding to a boundary point of the parameter space and the score function having an infinite second moment. For applications, the second form of regularity violation has the more important consequences, as it results in score statistics with non-standard normalisation and poor rates of convergence. The corresponding tests are difficult to use in practical situations because their asymptotic properties are unrepresentative of their behaviour for the sample sizes typical of applications, and extensive simulations may be needed in order to evaluate adequately their null distribution. Overcoming this difficulty is the primary focus of this paper. We propose a modification to the likelihood based approaches used by Tawn (Tawn, J.A., “Bivariate extreme value theory: models and estimation,” Biometrika 75, 397–415, 1988) and Ledford and Tawn (Ledford, A.W. and Tawn, J.A., “Statistics for near independence in multivariate extreme values,” Biometrika 83, 169–187, 1996) that provides asymptotically normal score tests of independence with regular normalisation and rapid convergence. The resulting tests are straightforward to implement and are beneficial in practical situations with realistic amounts of data. AMS 2000 Subject Classification Primary—60G70 Secondary—62H15     

水下爆炸气泡与复杂弹塑性结构的相互作用研究

计及结构的弹塑性,将边界元法(BEM)与有限元法(FEM)耦合提出了气泡与弹塑性结构耦合动力学计算方法,并开发了全套的三维水下气泡分析程序（UBA）,计算值与实验值之间误差在10%以内．以水面舰船为例,将三维计算程序工程化．并分析了水下爆炸气泡载荷作用下船体的弹塑性响应,从船体结构典型单元上的应力时历曲线可以看出,在气泡坍塌时出现应力峰值,证实了气泡坍塌压力及射流引起的压力对舰船等结构造成严重毁伤．从气泡与舰船的相互作用中可以看出,舰船低阶垂向振型被激起,在气泡作用下呈鞭状运动,同时舰船随着气泡的膨胀和收缩作升沉运动,通过本文的分析得到了适合于工程应用的规律及结论．     

Spatial Velocity and Temperature Measurements in Non‐Isothermal Liquid Flows

A surface tension driven flow in the liquid vicinity of an air bubble on a heated wall is studied experimentally. The liquid flow caused by the temperature gradient along the surface of the bubble is termed thermocapillary convection. The surface tension force and the buoyancy force oppose one another. The measurement technique is the 3D particle tracking velocity and thermometry, 3D PTV/T, using thermochromic liquid crystals and digital image processing. The paper describes the method in some detail and presents quantitative results for different Marangoni numbers. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Robustness in deterministic multi-objective linear programming with respect to the relative interior and angle deviation

This paper deals with the robustness issue in deterministic multi-objective linear programming from two new standpoints. It is shown that a robustness notion recently reported in the literature is equivalent to strict efficiency. Corresponding to an efficient solution, a new quantity, robustness order (RO) is defined with respect to the interiority order of the cost matrix in the binding cone. A linear programming problem is provided to calculate the RO of a given efficient solution. The second part of the paper is devoted to investigating the robustness with respect to the eligible angle deviation of the cost matrix in the binding cone. Theoretical results are given to obtain the maximum eligible angle deviation. Finally, the relationship between two above-mentioned robustness standpoints is established. To have a better geometrical view, we prove the results for single-objective LP problems at first, and then we extend them to the multi-objective case. In addition to the theoretical results, some clarifying examples are given.     

Stability analysis of frame structures by bubble-function method

In the stability analysis of frame structures, the results by conventional finite element method (FEM) in which one member is taken as one element are sometimes unavailable. This paper took a new basic function system with bubble functions as the shape function of a bar element to develop a bubble function finite element method (BFEM), in which the bending and the geometric stiffness matrices were derived from the principle of virtual work. Bubble functions are finite element modes that are located entirely within a single element and are zero on boundaries of the element, but are nonzero at the other points. BFEM is as concise as conventional bar FEM but has better accuracy, and is adaptable to the buckling analysis of all kinds of frame structures. The use of bubble functions significantly improves the convergence of finite element analysis, and efficiently reduces the computation cost for the buckling analysis of frame structures. Numerical results show that using bubble functions in finite element for the stability analysis of structures is very efficient, especially for high-rise and large-scale frame structures.     

Optimal bandwidth selection for conditional efficiency measures: A data-driven approach

In productivity analysis an important issue is to detect how external (environmental) factors, exogenous to the production process and not under the control of the producer, might influence the production process and the resulting efficiency of the firms. Most of the traditional approaches proposed in the literature have serious drawbacks. An alternative approach is to describe the production process as being conditioned by a given value of the environmental variables (Cazals, C., Florens, J.P., Simar, L., 2002. Nonparametric Frontier estimation: A robust approach. Journal of Econometrics 106, 1–25; Daraio, C., Simar, L., 2005. Introducing environmental variables in nonparametric Frontier models: A probabilistic approach. Journal of Productivity Analysis 24(1), 93–121). This defines conditional efficiency measures where the production set in the input ×$\times$ output space may depend on the value of the external variables. The statistical properties of nonparametric estimators of these conditional measures are now established (Jeong, S.O., Park, B.U., Simar, L., 2008. Nonparametric conditional efficiency measures: Asymptotic properties. Annals of Operations Research doi: 10.1007/s10479-008-0359-5). These involve the estimation of a nonstandard conditional distribution function which requires the specification of a smoothing parameter (a bandwidth). So far, only the asymptotic optimal order of this bandwidth has been established. This is of little interest for the practitioner. In this paper we fill this gap and we propose a data-driven technique for selecting this parameter in practice. The approach, based on a Least Squares Cross Validation procedure (LSCV), provides an optimal bandwidth that minimizes an appropriate (weighted) integrated Squared Error (ISE). The method is carefully described and exemplified with some simulated data with univariate and multivariate environmental factors. An application on real data (performances of Mutual Funds) illustrates how this new optimal method of bandwidth selection works in practice.     

Convergence Analysis of a Numerical Scheme for the Porous Medium Equation by an Energetic Variational Approach

The porous medium equation (PME)is a typical nonlinear degenerate parabolic equation. We have studied numerical methods for PME by an energetic variational approach in [C. Duan et al., J. Comput. Phys., 385 (2019), pp. 13–32], where the trajectory equation can be obtained and two numerical schemes have been developed based on different dissipative energy laws. It is also proved that the nonlinear scheme, based on $f$log $f$ as the total energy form of the dissipative law, is uniquely solvable on an admissible convex set and preserves the corresponding discrete dissipation law. Moreover, under certain smoothness assumption, we have also obtained the second order convergence in space and the first order convergence in time for the scheme. In this paper, we provide a rigorous proof of the error estimate by a careful higher order asymptotic expansion and two step error estimates. The latter technique contains a rough estimate to control the highly nonlinear term in a discrete $W$^1,∞norm and a refined estimate is applied to derive the optimal error order.