Solution of the hyperbolic mild‐slope equation using the finite volume method

A finite volume solver for the 2D depth‐integrated harmonic hyperbolic formulation of the mild‐slope equation for wave propagation is presented and discussed. The solver is implemented on unstructured triangular meshes and the solution methodology is based upon a Godunov‐type second‐order finite volume scheme, whereby the numerical fluxes are computed using Roe's flux function. The eigensystem of the mild‐slope equations is derived and used for the construction of Roe's matrix. A formulation that updates the unknown variables in time implicitly is presented, which produces a more accurate and reliable scheme than hitherto available. Boundary conditions for different types of boundaries are also derived. The agreement of the computed results with analytical results for a range of wave propagation/transformation problems is very good, and the model is found to be virtually paraxiality‐free. Copyright © 2003 John Wiley & Sons, Ltd.     

Assessment of some iterative methods for non-symmetric linear systems arising in computational fluid dynamics

Various tests have been carried out in order to compare the performances of several methods used to solve the non-symmetric linear systems of equations arising from implicit discretizations of CFD problems, namely the scalar advection-diffusion equation and the compressible Euler equations. The iterative schemes under consideration belong to three families of algorithms: relaxation (Jacobi and Gauss-Seidel), gradient and Newton methods. Two gradient methods have been selected: a Krylov subspace iteration method (GMRES) and a non-symmetric extension of the conjugate gradient method (CGS). Finally, a quasi-Newton method has also been considered (Broyden). The aim of this paper is to provide indications of which appears to be the most adequate method according to the particular circumstances as well as to discuss the implementation aspects of each scheme.     

Superconvergence of discontinuous Galerkin finite element method for the stationary Navier‐Stokes equations

This article focuses on discontinuous Galerkin method for the two‐ or three‐dimensional stationary incompressible Navier‐Stokes equations. The velocity field is approximated by discontinuous locally solenoidal finite element, and the pressure is approximated by the standard conforming finite element. Then, superconvergence of nonconforming finite element approximations is applied by using least‐squares surface fitting for the stationary Navier‐Stokes equations. The method ameliorates the two noticeable disadvantages about the given finite element pair. Finally, the superconvergence result is provided under some regular assumptions. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 421–436, 2007     

Incremental unknowns preconditioning for solving the Helmholtz equation

An efficient preconditioner is developed for solving the Helmholtz problem in both high and low frequency (wavenumber) regimes. The preconditioner is based on hierarchical unknowns on nested grids, known as incremental unknowns (IU). The motivation for the IU preconditioner is provided by an eigenvalue analysis of a simplified Helmholtz problem. The performance of our preconditioner is tested on the iterative solution of two‐dimensional electromagnetic scattering problems. When compared with other well‐known methods, our technique is shown to often provide a better numerical efficacy and, most importantly, to be more robust. Moreover, for the best performance, the number of IU levels used in the preconditioner should be designed for the coarsest grid to have roughly two points per linear wavelength. This result is consistent with the conventional sampling criteria for wave phenomena in contrast with existing IU applications for solving the Laplace/Poisson problem, where the coarsest grid comprises just one interior point. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

曲边介质波导的有限元分析

本文应用曲边有限元方法有效地计算了曲边质波导的色散特性，编制了相应的计算程序,应用本文方法计算得到的椭圆旬质波导和组合椭圆介质波导的主模和高次模的色散曲线与文献[2],[3]的结果相当吻合。     

大型钢结构人行景观天桥的有限元分析

采用有限元方法对大型钢结构人行景观天桥进行了静力学、温度应力、动力学模态和初步的抗震能力分析,并就桥自重作用下的最大竖向位移和实际观测结果进行了比较,两者是一致的.计算结果表明:人行景观天桥在常规、温度极值改变和多遇地震荷载作用等多种工况下,现有的设计方案具有足够的强度、刚度和抗震能力.     

Uniform convergence of Galerkin‐multigrid methods for nonconforming finite elements for second‐order problems with less than full elliptic regularity

In this article we prove uniform convergence estimates for the recently developed Galerkin‐multigrid methods for nonconforming finite elements for second‐order problems with less than full elliptic regularity. These multigrid methods are defined in terms of the “Galerkin approach,” where quadratic forms over coarse grids are constructed using the quadratic form on the finest grid and iterated coarse‐to‐fine intergrid transfer operators. Previously, uniform estimates were obtained for problems with full elliptic regularity, whereas these estimates are derived with less than full elliptic regularity here. Applications to the nonconforming P₁, rotated Q₁, and Wilson finite elements are analyzed. The result applies to the mixed method based on finite elements that are equivalent to these nonconforming elements. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 18: 203–217, 2002; DOI 10.1002/num.10004     

Assessment of thermal expansion coefficient of methylsilsesquioxane thin film

For the methylsilsesquioxane film whose optical birefringence is almost zero, it was recently reported that its vertical thermal expansion coefficient (CTE) was approximately one order of magnitude larger than the lateral CTE. Though the birefringence is not an absolute predictor of anisotropic behavior, the discrepancy in both the CTEs was so remarkable that it was essential to investigate whether the anisotropy was intrinsic property or not. If the effect of Poisson's ratio is considered in the calculation of the vertical CTE and when elastic modulus measured by surface acoustic wave spectroscopy is used in the assessment of the lateral CTE, both the CTEs are coincident with each other. Therefore, it can be concluded that the discrepancy in the CTEs can be attributed to a higher in‐plane polymer chain orientation but it can also arise from the misleadingly assumed modulus and Poisson's ratio. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 3109–3120, 2006