High order generalized upwind schemes and numerical solution of singular perturbation problems

High even order generalizations of the traditional upwind method are introduced to solve second order ODE-BVPs without recasting the problem as a first order system. Both theoretical analysis and numerical comparison with central difference schemes of the same order show that these new methods may avoid typical oscillations and achieve high accuracy. Singular perturbation problems are taken into account to emphasize the main features of the proposed methods. AMS subject classification (2000)  65L10, 65L12, 65L50     

A new approach based on the surrogating method in the project time compression problems

This paper develops a mathematical model for project time compression problems in CPM/PERT type networks. It is noted this formulation of the problem will be an adequate approximation for solving the time compression problem with any continuous and non-increasing time-cost curve. The kind of this model is Mixed Integer Linear Program (MILP) with zero-one variables, and the Benders' decomposition procedure for analyzing this model has been developed. Then this paper proposes a new approach based on the surrogating method for solving these problems. In addition, the required computer programs have been prepared by the author to execute the algorithm. An illustrative example solved by the new algorithm, and two methods are compared by several numerical examples. Computational experience with these data shows the superiority of the new approach.     

On the solution of Maxwell's equations in polygonal domains

This paper is concerned with the structure of the singular and regular parts of the solution of time‐harmonic Maxwell's equations in polygonal plane domains and their effective numerical treatment. The asymptotic behaviour of the solution near corner points of the domain is studied by means of discrete Fourier transformation and it is proved that the solution of the boundary value problem does not belong locally to H² when the boundary of the domain has non‐acute angles. A splitting of the solution into a regular part belonging to the space H², and an explicitly described singular part is presented. For the numerical treatment of the boundary value problem, we propose a finite element discretization which combines local mesh grading and the singular field methods and derive a priori error estimates that show optimal convergence as known for the classical finite element method for problems with regular solutions. Copyright © 2006 John Wiley & Sons, Ltd.     

微分算子扰动问题的值域研究

本文研究作用在C2周期函数空间上的微分算子u→u″ g(u) ,其中g(u)为连续有界函数.我们将证明上述微分算子的值域限制在周期函数空间的“超曲面”中.     

Effect of plasma treatment on interface property of BCN/GaN structure

Interface properties of BCN/GaN metal-insulator-semiconductor (MIS) structures are investigated by X-ray photoelectron spectroscopy (XPS) and capacitance versus voltage (C-V) characteristics measurements. The BCN/GaN samples are fabricated by in situ process consisting of plasma treatment and deposition of BCN film in the plasma-assisted chemical vapor deposition (PACVD) apparatus. XPS measurement shows that the oxide formation at the BCN/GaN interface is suppressed by nitrogen (N₂) and hydrogen (H₂) plasma treatment. The interface state density is estimated from C-V characteristics measured at 1 MHz using Terman method. The minimum interface state density appears from 0.2 to 0.7 eV below the conduction band edge of GaN. The minimum value of the interface state density is estimated to be 3.0 × 10¹⁰ eV⁻¹ cm⁻² for the BCN/GaN structure with mixed N₂ and H₂ plasma treatment for 25 min. Even after annealing at 430 °C for 10 min, the interface state density as low as 6.0 × 10¹⁰ eV⁻¹ cm⁻² is maintained.     

Microstructure of laser melted Al18B4O33w/2024Al composite

The microstructure of a laser treated Al₁₈B₄O₃₃w/2024Al composite has been investigated using transmission electron microscope (TEM), low-angle (glancing angle) X-ray diffraction (XRD) techniques. Various surface microstructures were observed in the laser treated composite. The Al₁₈B₄O₃₃ whisker on the surface of the composite was decomposed during laser surface melting, various decomposition products were studied in the laser treated composite. Eutectic phases and the precipitation in the matrix of the composite with laser-treated were observed. The main phases detected in the molten zone were aluminum and decomposition products Al₂O₃. The effect of laser treatment on the hardness of the composite was also examined. A surface hardness of 400 Hv was noted.     

放射治疗剂量分布的Monte Carlo模拟

利用Monte Carlo方法模拟放射治疗剂量分布,改写现有的放射治疗软件DPM,并进行剂量计算。与图形处理方便的Matlab软件相结合,给出直观的显示结果。该方法是一个很有前景的研究方向,因为它既可以获得很快的仿真速度,又能得到直观的仿真结果图形。运行结果表明剂量计算能在1min之内完成,计算速度和精度满足要求。该方法对于提高放射治疗水平具有重要的指导意义和应用价值。     

A Bayesian analysis for the seismic data on Taiwan

A Bayesian approach is used to analyze the seismic events with magnitudes at least 4.7 on Taiwan. Following the idea proposed by Ogata (1988,Journal of the American Statistical Association,83, 9–27), an epidemic model for the process of occurrence times given the observed magnitude values is considered, incorporated with gamma prior distributions for the parameters in the model, while the hyper-parameters of the prior are essentially determined by the seismic data in an earlier period. Bayesian inference is made on the conditional intensity function via Markov chain Monte Carlo method. The results yield acceptable accuracies in predicting large earthquake events within short time periods.     

一类仓库货物储存模型经典解的存在唯一性

文章用特征线法给出一类仓库货物储存模型解的递推表达式,并证明其光滑性从而得到了经典解的存在唯一性。     

某库岸滑坡稳定分析及治理措施研究

从地质条件出发,详细分析了滑坡的成因,以大型的岩土工程数值分析软件(ADINA)为研究手段,建立有限元计算模型,用强度折减的方法分析了滑坡体在初始工况、蓄水工况、暴雨工况、水位骤降工况和地震工况下的稳定情况,根据滑坡的失稳情况,在综合考虑的基础上,提出了削坡减载和地表排水相结合的治理措施,计算结果表明,治理措施是可行的,研究结果对设计和施工具有较大的参考价值.