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.     

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     

Monte Carlo Hamiltonian: Linear Potentials

We further study the validity of the Monte Carlo Hamiltonian method. The advantage of the method,in comparison with the standard Monte Carlo Lagrangian approach, is its capability to study the excited states. Weconsider two quantum mechanical models: a symmetric one V(x) = |x|/2; and an asymmetric one V(x) = ∞, forx ＜ 0 and V(x) = x, for x ≥ 0. The results for the spectrum, wave functions and thermodynamical observables are inagreement with the analytical or Runge-Kutta calculations.     

A marker method for the solution of the damped Burgers' equation

A new method for the solution of the damped Burgers' equation is described. The marker method relies on the definition of a convective field associated with the underlying partial differential equation; the information about the approximate solution is associated with the response of an ensemble of markers to this convective field. Some key aspects of the method, such as the selection of the shape function and the initial loading, are discussed in some details. The marker method is applicable to a general class of nonlinear dispersive partial differential equations. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

On accurate boundary conditions for a shape sensitivity equation method

This paper studies the application of the continuous sensitivity equation method (CSEM) for the Navier–Stokes equations in the particular case of shape parameters. Boundary conditions for shape parameters involve flow derivatives at the boundary. Thus, accurate flow gradients are critical to the success of the CSEM. A new approach is presented to extract accurate flow derivatives at the boundary. High order Taylor series expansions are used on layered patches in conjunction with a constrained least‐squares procedure to evaluate accurate first and second derivatives of the flow variables at the boundary, required for Dirichlet and Neumann sensitivity boundary conditions. The flow and sensitivity fields are solved using an adaptive finite‐element method. The proposed methodology is first verified on a problem with a closed form solution obtained by the Method of Manufactured Solutions. The ability of the proposed method to provide accurate sensitivity fields for realistic problems is then demonstrated. The flow and sensitivity fields for a NACA 0012 airfoil are used for fast evaluation of the nearby flow over an airfoil of different thickness (NACA 0015). Copyright © 2005 John Wiley & Sons, Ltd.     

最大似然方法和Bayes方法在结构方程模型分析中的讨论(英文)

本文通过模拟研究,讨论了最大似然方法和Bayes方法在分析结构方程模型中的相似点和不同之处。     

氢原子的固有电偶极矩

用经典力学和双波量子力学计算了氢原子的固有电偶极矩。双波量子理论算得的结果在经典极限下与经典力学的结果一致。普通量子力学对氢原子Stark效应中表现出来的电偶极矩难以做出很好的解释，因为一个波函数描述的是系综而不是单个粒子。经典力学和双波量子力学可描述单个粒子的行为，对永久电偶极矩的计算和解释显得自然而合理。     

Deactivation reactions in the modeled 2,2,6,6‐tetramethyl‐1‐piperidinyloxy‐mediated free‐radical polymerization of styrene: A comparative study with the 2,2,6,6‐tetramethyl‐1‐piperidinyloxy/acrylonitrile system

The competitiveness of the combination and disproportionation reactions between a 1‐phenylpropyl radical, standing for a growing polystyryl macroradical, and a 2,2,6,6‐tetramethyl‐1‐piperidinyloxy (TEMPO) radical in the nitroxide‐mediated free‐radical polymerization of styrene was quantitatively evaluated by the study of the transition geometry and the potential energy profiles for the competing reactions with the use of quantum‐mechanical calculations at the density functional theory (DFT) UB3‐LYP/6‐311+G(3df, 2p)//(unrestricted) Austin Model 1 level of theory. The search for transition geometries resulted in six and two transition structures for the radical combination and disproportionation reactions, respectively. The former transition structures, mainly differing in the out‐of‐plane angle of the N? O bond in the transition structure TEMPO molecule, were correlated with the activation energy, which was determined to be in the range of 8.4–19.4 kcal mol^?1 from a single‐point calculation at the DFT UB3‐LYP/6‐311+G(3df, 2p)//unrestricted Austin Model 1 level. The calculated activation energy for the disproportionation reaction was less favorable by a value of more than 30 kcal mol^?1 in comparison with that for the combination reaction. The approximate barrier difference for the TEMPO addition and disproportionation reaction was slightly smaller for the styrene polymerization system than for the acrylonitrile polymerization system, thus indicating that a β‐proton abstraction through a TEMPO radical from the polymer backbone could diminish control over the radical polymerization of styrene with the nitroxide even more than in the latter system. © 2006 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 45: 232–241, 2007