647 research outputs found
A parallel compact-TVD method for compressible fluid dynamics employing shared and distributed-memory paradigms
A novel multi-block compact-TVD finite difference method for the simulation of compressible flows is presented. The method combines distributed and shared-memory paradigms to take advantage of the configuration of modern supercomputers that host many cores per shared-memory node. In our approach a domain decomposition technique is applied to a compact scheme using explicit flux formulas at block interfaces. This method offers great improvement in performance over earlier parallel compact methods that rely on the parallel solution of a linear system. A test case is presented to assess the accuracy and parallel performance of the new method
Grid convergence of high order methods for multiscale complex unsteady viscous compressible flows
Grid convergence of several high order methods for the computation of rapidly developing complex unsteady viscous compressible flows with a wide range of physical scales is studied. The recently developed adaptive numerical dissipation control high order methods referred to as the ACM and wavelet filter schemes are compared with a fifth-order weighted ENO (WENO) scheme. The two 2-D compressible full Navier–Stokes models considered do not possess known analytical and experimental data. Fine grid solutions from a standard second-order TVD scheme and a MUSCL scheme with limiters are used as reference solutions. The first model is a 2-D viscous analog of a shock tube problem which involves complex shock/shear/boundary-layer interactions. The second model is a supersonic reactive flow concerning fuel breakup. The fuel mixing involves circular hydrogen bubbles in air interacting with a planar moving shock wave. Both models contain fine scale structures and are stiff in the sense that even though the unsteadiness of the flows are rapidly developing, extreme grid refinement and time step restrictions are needed to resolve all the flow scales as well as the chemical reaction scales. Our computations were all made on uniform grids, and our conclusions cannot be directly carried over to, for example, curvilinear grids
CABARET in the ocean gyres
Author Posting. © The Author(s), 2009. This is the author's version of the work. It is posted here by permission of Elsevier B.V. for personal use, not for redistribution. The definitive version was published in Ocean Modelling 30 (2009): 155-168, doi:10.1016/j.ocemod.2009.06.009.A new high-resolution Eulerian numerical method is proposed for modelling quasigeostrophic ocean
dynamics in eddying regimes. The method is based on a novel, second-order non-dissipative and lowdispersive
conservative advection scheme called CABARET. The properties of the new method are
compared with those of several high-resolution Eulerian methods for linear advection and gas
dynamics. Then, the CABARET method is applied to the classical model of the double-gyre ocean
circulation and its performance is contrasted against that of the common vorticity-preserving Arakawa
method. In turbulent regimes, the new method permits credible numerical simulations on much coarser
computational grids.Supports from the Royal Society of London and from the Mary Sears Visitor Grant are
acknowledged by SK with gratitude. The work of VG was supported by the Russian Foundation for
Basic Research (RFBR), grant 06-01-00819a. Funding for PB was provided by the NSF grant 0725796
Numerical simulation of conservation laws with moving grid nodes: Application to tsunami wave modelling
In the present article we describe a few simple and efficient finite volume
type schemes on moving grids in one spatial dimension combined with appropriate
predictor-corrector method to achieve higher resolution. The underlying finite
volume scheme is conservative and it is accurate up to the second order in
space. The main novelty consists in the motion of the grid. This new dynamic
aspect can be used to resolve better the areas with large solution gradients or
any other special features. No interpolation procedure is employed, thus
unnecessary solution smearing is avoided, and therefore, our method enjoys
excellent conservation properties. The resulting grid is completely
redistributed according the choice of the so-called monitor function. Several
more or less universal choices of the monitor function are provided. Finally,
the performance of the proposed algorithm is illustrated on several examples
stemming from the simple linear advection to the simulation of complex shallow
water waves. The exact well-balanced property is proven. We believe that the
techniques described in our paper can be beneficially used to model tsunami
wave propagation and run-up.Comment: 46 pages, 7 figures, 7 tables, 94 references. Accepted to
Geosciences. Other author's papers can be downloaded at
http://www.denys-dutykh.com
- …