High-Order Leap-Frog Based Discontinuous Galerkin Method for the Time-Domain Maxwell Equations on Non-Conforming Simplicial Meshes

A high-order leap-flog based non-dissipative discontinuous Galerkin timedomain method for solving Maxwell's equations is introduced and analyzed. The proposed method combines a centered approximation for the evaluation of fluxes at the interface between neighboring elements, with a Nth-order leap-frog time scheme. Moreover, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary level hanging nodes. The method is proved to be stable under some CFL-like condition on the time step. The convergence of the semi-discrete approximation to Maxwell's equations is established rigorously and bounds on the global divergence error are provided. Numerical experiments with highorder elements show the potential of the method.     

A compact high‐order HLL solver for nonlinear hyperbolic systems

We present a new finite‐volume method for calculating complex flows on non‐uniform meshes. This method is designed to be highly compact and to accurately capture all discontinuities that may arise within the solution of a nonlinear hyperbolic system. In the first step, we devise a fourth‐degree Hermite polynomial to interpolate the solution. The coefficients defining this polynomial are calculated by using a least‐square method. To introduce monotonicity conditions within the procedure, two constraints are added into the least‐square system. Those constraints are derived by locally matching the high‐order Hermite polynomial with a low‐order TVD polynomial. To emulate these constraints only in regions of discontinuities, data‐depending weights are defined; these weights are based upon normalized indicators of smoothness of the solution and are parameterized by an O(1) quantity. The reconstruction so generated is highly compact and is fifth‐order accurate when the solution is smooth; this reconstruction becomes first order in regions of discontinuities. In the second step, this reconstruction is inserted in an HLL approximate Riemann solver. This solver is designed to correctly capture all discontinuities that may arise into the solution. To this aim, we introduce the contribution of a possible contact discontinuity into the HLL Riemann solver. Thus, a spatially fifth‐order non‐oscillatory method is generated. This method evolves in time the solution and its first derivative. In a one‐dimensional context, a linear spectral analysis and extensive numerical experiments make it possible to assess the robustness and the advantages of the method in computing multi‐scale problems with embedded discontinuities. Copyright © 2008 John Wiley & Sons, Ltd.     

Stokes问题各向异性网格Q2-P1混合元超收敛分析

讨论Stokes问题在各向异性冈格下的Q2-P1混合有限元方法，利用积分恒等式技巧得到了与传统方法相同的超逼近性质，同时基于插值后处理的技巧，构造了速度和压力的一对插值后处理算子，并且前者具有备向异性特征，从而导出了整体超收敛结果．     

K2Fe2B2O7: A transparent nonlinear optical crystal with frustrated magnetism

A new noncentrosymmetric ferroborate crystal, K₂Fe₂B₂O₇, has been grown from high temperature melt. Structure solution from single crystal X-ray diffraction shows that the title compound crystallizes in a trigonal space group P321 with cell dimensions of a=8.7475(12) Å and c=8.5124(17) Å. In the structure, FeO₄ tetrahedron shares its three basal oxygen atoms with BO₃ triangles forming a two-dimensional layer in the ab plane and the layers are connected by the apical Fe-O bonds along the c direction. The crystal is transparent in the visible and near infrared region from 500 to 2000 nm with three pronounced absorption bands ascribed to d-d transitions of tetrahedrally coordinated Fe³⁺ ions. Though, structurally analog to K₂Al₂B₂O₇, the further twisting of the BO₃ groups between adjacent layers reduces its optical nonlinearity to a second-harmonic generation intensity of about 0.4 times that of KDP. Spin-glass behavior is observed at 20 K which is probably due to geometrically magnetic frustration of the triangular Fe net in the ab plane.     

Anti-diffusion method for interface steepening in two-phase incompressible flow

In this paper, we present a method for obtaining sharp interfaces in two-phase incompressible flows by an anti-diffusion correction, that is applicable in a straight-forward fashion for the improvement of two-phase flow solution schemes typically employed in practical applications. The underlying discretization is based on the volume-of-fluid (VOF) interface-capturing method on unstructured meshes. The key idea is to steepen the interface, independently of the underlying volume-fraction transport equation, by solving a diffusion equation with reverse time, i.e. an anti-diffusion equation, after each advection time step of the volume fraction. As the solution of the anti-diffusion equation requires regularization, a limiter based on the directional derivative is developed for calculating the gradient of the volume fraction. This limiter ensures the boundedness of the volume fraction. In order to control the amount of anti-diffusion introduced by the correction algorithm we propose a suitable stopping criterion for interface steepening. The formulation of the limiter and the algorithm for solving the anti-diffusion equation are applicable to 3-dimensional unstructured meshes. Validation computations are performed for passive advection of an interface, for 2-dimensional and 3-dimensional rising-bubbles, and for a rising drop in a periodically constricted channel. The results demonstrate that sharp interfaces can be recovered reliably. They show that the accuracy is similar to or even better than that of level-set methods using comparable discretizations for the flow and the level-set evolution. Also, we observe a good agreement with experimental results for the rising drop where proper interface evolution requires accurate mass conservation.     

Selecting green suppliers based on GSCM practices: Using fuzzy TOPSIS applied to a Brazilian electronics company

Due to an increased awareness and significant environmental pressures from various stakeholders, companies have begun to realize the significance of incorporating green practices into their daily activities. This paper proposes a framework using Fuzzy TOPSIS to select green suppliers for a Brazilian electronics company; our framework is built on the criteria of green supply chain management (GSCM) practices. An empirical analysis is made, and the data are collected from a set of 12 available suppliers. We use a fuzzy TOPSIS approach to rank the suppliers, and the results of the proposed framework are compared with the ranks obtained by both the geometric mean and the graded mean methods of fuzzy TOPSIS methodology. Then a Spearman rank correlation coefficient is used to find the statistical difference between the ranks obtained by the three methods. Finally, a sensitivity analysis has been performed to examine the influence of the preferences given by the decision makers for the chosen GSCM practices on the selection of green suppliers. Results indicate that the four dominant criteria are Commitment of senior management to GSCM; Product designs that reduce, reuse, recycle, or reclaim materials, components, or energy; Compliance with legal environmental requirements and auditing programs; and Product designs that avoid or reduce toxic or hazardous material use.     

A modified Bakhvalov mesh

We consider a few numerical methods for solving a one-dimensional convection–diffusion singularly perturbed problem. More precisely, we introduce a modified Bakvalov mesh generated using some implicitly defined functions. Properties of this mesh and convergence results for several methods on it are given. Numerical results are presented in support of the theoretical considerations.     

比例延迟微分方程的连续有限元法

利用连续有限元法求解比例延迟微分方程,在一致网格下,给出比例延迟微分方程连续有限元解的整体收敛阶,数值实验验证了理论结果的正确性.