The long computation times required for fluid dynamics simulations (CFD) has lead the industry to look for some alternatives to boost high performance computing (HPC). This paper is focused on the acceleration of fluid dynamics simulations for industrial complex configurations using modern graphics cards (GPUs) that exhibits a substantial parallel architecture. First, some general facts about g...

In this paper we treat subjects which are relevant in the context of iterative methods in implicit time integration for compressible flow simulations. We present a novel renumbering technique, some strategies for choosing the time step in the implicit time integration, and a novel implementation of a matrix-free evaluation for matrix-vector products. For the linearized compressible Euler equati...

Low speed gas flows are studied using a Unified Flow Solver (UFS) that automatically selects continuum-fluiddynamic and kinetic solvers for near-equilibrium and kinetic domains. Direct Numerical Solution of the Boltzmann equation is used in the kinetic domains at moderate and high local Knudsen number, while kinetic schemes for Euler or Navier Stokes equations are used in the continuum domains....

We present core components of a hair modeling and dynamics solution for the feature lm industry. Recent research results in hair simulation are exploited as a dynamics model based on solving the Euler-Lagrange equations of motion for a discretized Cosserat curve is implemented in its entirety. Solutions to the dynamics equations are derived and a framework for symbolic integration is outlined. ...

The application of Krylov subspace iterative methods to unsteady three-dimensional Navier–Stokes codes on massively parallel and distributed computing environments is investigated. Previously, the Euler mode of the Navier–Stokes flow solver Transonic Unsteady Rotor Navier–Stokes (TURNS) has been coupled with a Newton–Krylov scheme which uses two Conjugate-Gradient-like (CG) iterative methods. F...

A recent work of Li et al. [Numer. Math. Theor. Meth. Appl., Vol. 1, pp. 92-112(2008)] proposed a finite volume solver to solve 2D steady Euler equations. Although the Venkatakrishnan limiter is used to prevent the non-physical oscillations nearby the shock region, the overshoot or undershoot phenomenon can still be observed. Moreover, the numerical accuracy is degraded by using Venkatakrishnan...

A scheme for the twedimensional Euler equations that uses flow parameters to determine the direction for upwind-differencing is described. This approach reduces the grid-dependence of conventional schemes. Upwinding at the local flow angle and the local pressure-gradient angle is tested. The upwindbiased fluxes are calculated with Roe's approximate Riemann solver. A centered flux is used for th...

We discuss a leading-edge model used in the computation of the run-out length of dry-flowing avalanches. The model has the form of a singular initial value problem for a scalar ordinary differential equation describing the avalanche dynamics. Existence, uniqueness and smoothness properties of the analytical solution are shown. We also prove the existence of a unique root of the solution. Moreov...

Interaction between a shock wave and an inhomogeneity in fluid has complicated behavior, including vortex and turbulence generating, mixing, shock wave scattering and reflection. In the present paper we deal with the numerical simulation of the considered process. The Euler equations of unsteady inviscid compressible three-dimensional flow are used into the four-equation model of multicomponent...

The vibration problems of uniform and nonuniform Euler-Bernoulli beams have been solved analytically or approximately [1-5] for various end conditions. In order to calculate fundamental natural frequencies and related mode shapes, well known variational techniques such as Rayleigh_Ritz and Galerkin methods have been applied in the past. Besides these techniques, some discretized numerical metho...