A non-oscillatory balanced scheme for an idealized tropical climate model

We use the non-oscillatory balanced numerical scheme developed in Part I to track the dynamics of a dry highly nonlinear barotropic/baroclinic coupled solitary wave, as introduced by Biello and Majda (2004), and of the moisture fronts of Frierson et al. (2004) in the presence of dry gravity waves, a barotropic trade wind, and the beta effect. It is demonstrated that, for the barotropic/baroclinic solitary wave, except for a little numerical dissipation, the scheme utilized here preserves total energy despite the strong interactions and exchange of energy between the baroclinic and barotropic components of the flow. After a short transient period where the numerical solution stays close to the asymptotic predictions, the flow develops small scale eddies and ultimately becomes highly turbulent. It is found here that the interaction of a dry gravity wave with a moisture front can either result in a reflection of a fast moistening front or the pure extinction of the precipitation. The barotropic trade wind stretches the precipitation patches and increases the lifetime of the moisture fronts which decay naturally by the effects of dissipation through precipitation while the Coriolis effect makes the moving precipitation patches disappear and appear at other times and places.     

Development and assessment of a variable-order non-oscillatory scheme for convection term discretization

A new scheme for convection term discretization is developed, called VONOS (variable-order non-oscillatory scheme). The development of the scheme is based on the behaviour of well-known non-oscillatory schemes in the pure convection of a step profile test case. The new scheme is a combination of the QUICK and BSOU (bounded second-order upwind) schemes. These two schemes do not have the same formal order of accuracy and for that reason the formal order of accuracy of the new scheme is variable. The scheme is conservative, bounded and accurate. The performance of the new scheme was assessed in three test cases. The results showed that it is more accurate than currently used higher-order schemes, so it can be used in a general purpose algorithm in order to save computational time for the same level of accuracy. © 1998 John Wiley & Sons, Ltd.     

Experiments with three-dimensional riblets as an idealized model of shark skin

The skin of fast sharks exhibits a rather intriguing three-dimensional rib pattern. Therefore, the question arises whether or not such three-dimensional riblet surfaces may produce an equivalent or even higher drag reduction than straight two-dimensional riblets. Previously, the latter have been shown to reduce turbulent wall shear stress by up to 10%. Hence, the drag reduction by three-dimensional riblet surfaces is investigated experimentally. Our idealized 3D-surface consists of sharp-edged fin-shaped elements arranged in an interlocking array. The turbulent wall shear stress on this surface is measured using direct force balances. In a first attempt, wind tunnel experiments with about 365,000 tiny fin elements per test surface have been carried out. Due to the complexity of the surface manufacturing process, a comprehensive parametric study was not possible. These initial wind tunnel data, however, hinted at an appreciable drag reduction. Subsequently, in order to have a better judgement on the potential of these 3D-surfaces, oil channel experiments are carried out. In our new oil channel, the geometrical dimensions of the fins can be magnified 10 times in size as compared to the initial wind tunnel experiments, i.e., from typically 0.5 mm to 5 mm. For these latter oil channel experiments, novel test plates with variable fin configuration have been manufactured, with 1,920–4,000 fins. This enhanced variability permits measurements with a comparatively large parameter range. As a result of our measurements, it can be concluded, that 3D-riblet surfaces do indeed produce an appreciable drag reduction. We found as much as 7.3% decreased turbulent shear stress, as compared to a smooth reference plate. However, in direct comparison with 2D riblets, the performance of 3D-riblets is still inferior by about 1.7%. On the other hand, it appears conceivable, with a careful design of the fin shape (possibly supported by theory), that this inferiority in performance might be reduced. Nevertheless, at present, it seems to be rather unlikely, that 3D-riblets can significantly outperform 2D-riblets. Finally, one interesting finding remains to be mentioned: The optimum drag reduction for short 3D-riblets occurs at a lower rib height than for longer 3D-riblets or for infinitely long 2D-riblets. The same observation had been made previously on shark scales of different species with differing rib lengths, but no explanation could be given. Received: 1 March 1999/Accepted: 16 July 1999     

Flow behavior and deposition study of pollen-shape carrier particles in an idealized inhalation path model

Pollen-shape (spiked sphere) hydroxyapatite (HA) particles for drug carrier application are studied. The particle shape and size effect on flow characteristics and deposition are assessed. The pollen-shape HA particles are synthesized to have comparable size as typical carrier particles with mean diameter of 30-50 m and effective density less than 0.3 g/cm 3 . The flow behaviors of HA and commonly used lactose (LA) carrier particles are characterized by the Carr's compressibility index (CI). The HA particle...     

Assessment of an unstructured exponential scheme discrete ordinates radiation model for non-gray media

For radiative transfer in complex geometries, Sakami and his co-workers have developed a discrete ordinates method (DOM) exponential scheme for unstructured meshes which was mainly applied to gray media. The present study investigates the application of the unstructured exponential scheme to a wider range of non-gray scenarios found in fire and combustion applications, with the goal to implement it in an in-house Computational Fluid Dynamics (CFD) code for fire simulations. The original unstructured gray exponential scheme is adapted to non-gray applications by employing a statistical narrow-band/correlated-k (SNB-CK) gas model and meshes generated using the authors’ own mesh generator. Different non-gray scenarios involving spectral gas absorption by H₂O and CO₂ are investigated and a comparative analysis is carried out between heat flux and radiative source terms predicted and literature data based on ray-tracing and Monte Carlo methods. The maximum discrepancies for total radiative heat flux do not typically exceed 5%.     

A mixed finite element scheme for viscoelastic flows with XPP model

A mixed finite element formulation for viscoelastic flows is derived in this paper, in which the FIC （finite incremental calculus） pressure stabilization process and the DEVSS （discrete elastic viscous stress splitting） method using the Crank-Nicolson-based split are introduced within a general framework of the iterative version of the fractional step algorithm. The SU （streamline-upwind） method is particularly chosen to tackle the convective terms in constitutive equations of viscoelastic flows. Thanks to the proposed scheme the finite elements with equal low-order interpolation approximations for stress-velocity-pressure variables can be successfully used even for viscoelastic flows with high Weissenberg numbers. The XPP （extended Pom-Pom） constitutive model for describing viscoelastic behaviors is particularly integrated into the proposed scheme. The numerical results for the 4：1 sudden contraction flow problem demonstrate prominent stability, accuracy and convergence rate of the proposed scheme in both pressure and stress distributions over the flow domain within a wide range of the Weissenberg number, particularly the capability in reproducing the results, which can be used to explain the ＂die swell＂ phenomenon observed in the polymer injection molding process.     

A semi-implicit lagrangian scheme for 3D shallow water flow with a two-layer turbulence model

A semi-implicit Lagrangian finite difference scheme for 3D shallow water flow has been developed to include an eddy viscosity model for turbulent mixing in the vertical direction. The α-co-ordinate system for the vertical direction has been introduced to give accurate definition of bed and surface boundary conditions. The simple two-layer mixing length model for rough surfaces is used with the standard assumption that the shear stress across the wall region at a given horizontal location is constant. The bed condition is thus defined only by its roughness height (avoiding the need for a friction formula relating to depth-averaged flow, e.g. Chezy, used previously). The method is shown to be efficient and stable with an explicit Lagrangian formulation for convective terms and terms for surface elevation and vertical mixing handled implicitly. The method is applied to current flow around a circular island with gently sloping sides which produce periodic recirculation zones (vortex shedding). Comparisons are made with experimental measurements of velocity using laser Doppler anemometry (time histories at specific points) and surface particle-tracking velocimetry.     

A shock-capturing scheme for body-fitted meshes

A finite difference scheme based on flux difference splitting is presented for the solution of the two-dimensional Euler equations of gas dynamics in a generalized co-ordinate system. The scheme is based on numerical characteristic decomposition and solves locally linearized Riemann problems using upwind differencing. The decomposition is for a generalized co-ordinate system and a convex equation of state. This ensures good shock-capturing properties when incorporated with operator splitting and the advantage of using body-fitted co-ordinates. The resulting scheme is applied to supersonic flow of real air' past a circular cylinder.     

一种小型静电陀螺仪三位一体系统方案

在比较了国外两种小型静电陀螺仪方案的基础上，研究了一种小型静电陀螺仪质量不平衡调制信号读取、数字式支承和电场恒速三位一体工作系统方案。介绍了三位一体系统的工作原理，讨论了系统各组成模块的功能，对制约系统性能的关键因素作了相关分析。分析表明，在抵消变压器线圈阻抗匹配的条件下，可以忽略加力和信号读取之间的相互影响。     

A novel scheme for nonlinear displacement-dependent dampers

In this paper, a novel scheme for nonlinear displacement-dependent (NDD) damper is introduced. The damper is attached to a simple mass-spring-damper vibration system. The vibration system equipped with a NDD damper is mathematically modeled and the nonlinear governing differential equation of the system is derived. To obtain the displacement of the system, the approximate analytical solution of the governing equation is elaborated using the multiple scales method. The advised approximate analytical algorithm is performed for several case studies and is also verified by the numerical fourth-order Runge?CKutta method. In addition, the performance of the NDD damper is analyzed and compared with the performance of the traditional linear damper. It is found that the proposed NDD damper scheme along with the multiple scales method is not only feasible for vibration reduction but also yields satisfactory response performance rather than the existing traditional linear damper.     

A pressure-smoothing scheme for incompressible flow problems

A pressure-smoothing scheme for Stokes and Navier–Stokes flows of Newtonian fluids and for Stokes flow of Maxwell fluids is described. The stress deviator obtained from the calculated velocity field is substituted into the governing equilibrium equation. The resulting equation is then solved to obtain a new, smoothed pressure by a least square finite element method.     

A split operator scheme for ocean wave simulation

A coupled discrete spectral model was developed for the prediction of ocean waves by solving the energy conservation equation of the two-dimensional wave spectrum. The model includes the dispersion correction terms in the governing equation to account for the dispersive effect due to the frequency-dependent velocities of waves. A split operator scheme is used to deal with the numerical problems arising from different terms of the governing equation. The advection terms are solved by the proven accurate minimax characteristics method to avoid excessive numerical damping or oscillations. The dispersion correction terms are solved by central differencing. The source and sink terms are solved by a quasi-second-order explicit scheme with limitation on energy growth per time step to allow the use of a large time step. The model was verified by ideal test cases and wave-hindcasting studies under typhoon conditions in the South China Sea near Hong Kong.     

A two-step godunov-type scheme for the euler equations

A second-order Godunov-type scheme for the Euler equations in conservation form is derived. The method is based on the ENO formulation proposed by Harten et al. The fundamental difference lies in the use of a two-step scheme to compute the time evolution. The scheme is TVD in the linear scalar case, and gives oscillation-free solutions when dealing with nonlinear hyperbolic systems. The admissible time step is twice that of classical Godunovtype schemes. This feature makes it computationally cheaper than one-step schemes, while requiring the same computer storage.

---

Sommario Viene data una nuova estensione al secondo ordine del metodo di Godunov per la soluzione delle equazioni di Eulero in forma conservativa. Il metodo é basato sulla formulazione ENO proposta da Harten et al. La differenza fondamentale consiste nel calcolo dell'evoluzione temporale, ottenuta mediante uno schema a due passi. Questo consente l'uso di un passo di integrazione nel tempo doppio rispetto agli altri schemi alla Godunov ad un solo passo. Il metodo proposto risulta quindi piú efficiente e puó inoltre essere implementato senza alcun aumento dell'occupazione di memoria. Viene dimostrato che lo schema é TVD nel caso lineare, e che fornisce soluzioni prive di oscillazioni spurie nel caso di sistemi non-lineari.

     

A finite volume scheme for two-phase compressible flows

In gas-particle two-phase flows, when the concentration of the disperesed phase is low, certain assumptions may be made which simplify considerably the equations one has to solve. The gas and particle flows are then linked only via the interaction terms. One may therefore uncouple the full system of equations into two subsystems: one for the gas phase, whose homogeneous part reduces to the Euler equations; and a second system for the particle motion, whose homogeneous part is a degenerate hyperbolic system. The equations governing the gas phase flow may be solved using a high-resolution scheme, while the equations describing the motion of the dispersed phase are treated by a donor-cell method using the solution of a particular Riemann problem. Coupling is then achieved via the right-hand-side terms. To illustrate the capabilities of the proposed method, results are presented for a case specially chosen from among the most difficult to handle, since it involves certain geometrical difficulties, the treatment of regions in which particles are absent and the capturing of particle fronts.