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.     

New efficient numerical procedures for solving stochastic variational problems with a priori maximum pointwise error estimates

In a previous paper we gave a new formulation and derived the Euler equations and other necessary conditions to solve strong, pathwise, stochastic variational problems with trajectories driven by Brownian motion. Thus, unlike current methods which minimize the control over deterministic functionals (the expected value), we find the control which gives the critical point solution of random functionals of a Brownian path and then, if we choose, find the expected value.This increase in information is balanced by the fact that our methods are anticipative while current methods are not. However, our methods are more directly connected to the theory and meaningful examples of deterministic variational theory and provide better means of solution for free and constrained problems. In addition, examples indicate that there are methods to obtain nonanticipative solutions from our equations although the anticipative optimal cost function has smaller expected value.In this paper we give new, efficient numerical methods to find the solution of these problems in the quadratic case. Of interest is that our numerical solution has a maximal, a priori, pointwise error of O(h3/2) where h is the node size. We believe our results are unique for any theory of stochastic control and that our methods of proof involve new and sophisticated ideas for strong solutions which extend previous deterministic results by the first author where the error was O(h²).We note that, although our solutions are given in terms of stochastic differential equations, we are not using the now standard numerical methods for stochastic differential equations. Instead we find an approximation to the critical point solution of the variational problem using relations derived from setting to zero the directional derivative of the cost functional in the direction of simple test functions.Our results are even more significant than they first appear because we can reformulate stochastic control problems or constrained calculus of variations problems in the unconstrained, stochastic calculus of variations formulation of this paper. This will allow us to find efficient and accurate numerical solutions for general constrained, stochastic optimization problems. This is not yet being done, even in the deterministic case, except by the first author.     

Numerical analysis of time-dependent Boussinesq models

In this paper we analyse numerical models for time-dependent Boussinesq equations. These equations arise when so-called Boussinesq terms are introduced into the shallow water equations. We use the Boussinesq terms proposed by Katapodes and Dingemans. These terms generalize the constant depth terms given by Broer. The shallow water equations are discretized by using fourth-order finite difference formulae for the space derivatives and a fourth-order explicit time integrator. The effect on the stability and accuracy of various discrete Boussinesq terms is investigated. Numerical experiments are presented in the case of a fourth-order Runge-Kutta time integrator.     

Simplified analytic solutions and a novel fast algorithm for Yb-doped double-clad fiber lasers

In this paper, continuous wave Yb³⁺-doped double-clad fiber lasers (DCFLs) with linear-cavity are investigated theoretically and numerically using the rate equations. Under the steady state conditions, the simplified analytic solutions of Yb³⁺-doped DCFLs under considering the scattering loss are deduced in the strongly pump condition. Compared with the known analytic solutions in published literatures, our analytic solutions are more accurate, especially, at higher reflectivity of output mirror. In addition, a fast and stable algorithm based on the Newton-Raphson method is proposed to simulate numerically Yb³⁺-doped DCFLs. The results by simplified analytic solutions are in good agreement with those by the numerical simulation. Moreover, we have performed the optimization of an Yb³⁺-doped DCFL using the simplified analytic solutions and the numerical simulations, respectively.     

Interpolation on spherical geodesic grids: A comparative study

Most operational models in atmospheric physics, meteorology and climatology nowadays adopt spherical geodesic grids and require “ad hoc” developed interpolation procedures. The author does a comparison between chosen representatives of linear, distance-based and cubic interpolation schemes outlining their advantages and drawbacks in this specific application field. Numerical experiments on a standard test problem, while confirming a good performance of linear and distance-based schemes in a single interpolation step, also show their minor accuracy with respect to the cubic scheme in the more realistic simulation of advection of a meteorological field.     

Numerical simulation of soft palate movement and airflow in human upper airway by fluid-structure interaction method

In this paper, the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing. On the basis of the data taken from the spiral computerized tomography images of a healthy person and a patient with Obstructive Sleep Apnea-Hypopnea Syndrome (OSAHS), three-dimensional models of upper airway cavity and soft palate are reconstructed by the method of surface rendering. Numerical simulation is performed for airflow in the upper airway and displacement of soft palate by fluid-structure interaction analysis. The reconstructed three-dimensional models precisely preserve the original configuration of upper airways and soft palate. The results of the pressure and velocity distributions in the airflow field are quantitatively determined, and the displacement of soft palate is presented. Pressure gradients of airway are lower for the healthy person and the airflow distribution is quite uniform in the case of free breathing. However, the OSAHS patient remarkably escalates both the pressure and velocity in the upper airway, and causes higher displacement of the soft palate. The present study is useful in revealing pathogenesis and quantitative mutual relationship between configuration and function of the upper airway as well as in diagnosing diseases related to anatomical structure and function of the upper airway. The project supported by the National Natural Science Foundation of China (10672036, 10472025 and 10421002), the Natural Science Foundation of Liaoning Province (20032109). English text was polished by Yunming Chen.     

Nonlinear modelling and control of a common rail injection system for diesel engines

Modern trends in designing mechatronic systems call for a synergic design of the separated subsystems (mechanic, electronic parts, control modules, etc.) concurring to the overall performance. Following this point of view, this paper presents a control oriented model and a nonlinear control design for a Common Rail injection system. First a model is developed, which is tuned in a virtual simulation environment, representing the injection system in details in a reliable replication of reality. Then a sliding mode control is developed. Both the model of the injection process and of the control law are validated by a virtual detailed simulation environment. The prediction capability of the model and the control efficiency are clearly shown.     

On the Numerical Range of Operator Polynomials

An upper bound for the number of components of the numerical range of matrix polynomials is established. We also establish a necessary condition and a sufficient condition for the connectedness of the numerical range of a quadratic selfadjoint matrix polynomial. The boundary of the numerical range of linear matrix polynomials is also considered.     

The Effect of Boundary Conditions on Gas-phase Synthesised Silver Nanoparticles

We have prepared spherical non-agglomerated silver nanoparticles by an evaporation–condensation–dilution/cooling technique. Silver was evaporated from a crucible in a tubular flow reactor. A porous tube diluter was used to quench the carrier gas at the outlet of the reactor to enhance the formation of small particles and to suppress agglomeration and other particle growth mechanisms. The number size distribution of the prepared particles was measured with a differential mobility analyser–condensation nucleus counter combination and the size and the shape of the particles were analysed with transmission electron microscope. The system was modelled using a sectional aerosol dynamics computer code to estimate the importance of different aerosol processes. In all conditions the particles obtained were non-agglomerated and spherical. The mean particle diameter varied from 4 to 10-nm depending on boundary conditions. From the modelling studies it can be concluded that the nucleation rate is the most important parameter controlling the final particle size.     

Ab initio modelling of crosslinking in polymers. A case of chains with furan rings

A quantum chemical model (Abinitio HF-MO and DFT-B3LYP) of polyfurfuryl alcohol crosslinking is shown. Two pathways were considered: (a) a Diels–Alder (D–A) reaction between a dihydrofuran moiety and a furan ring; (b) the addition reaction of electrophilic specimens on conjugated double bonds. The ability for the formation of D–A adducts (dienes and dienophiles) was investigated by the frontier molecular orbital (FMO) theory. Energy gaps between the FMO's of diene and dienophile structures suggest that the occurrence of the D–A reaction is possible. The results suggest that the carbons of the exo double bonds attached to dihydrofuran rings are the sites most likely to be attacked by electrophilic species, while the C₃ and C₄ atoms display the opposite tendency due to their low HOMO electron densities. Results of a thermochemical approach to these reactions at the MP2 level was in agreement with reactivity modelled by MO predictions.