1,382 research outputs found

    On the use of finite fault solution for tsunami generation problems

    Get PDF
    The present study is devoted to the problem of tsunami wave generation. The main goal of this work is two-fold. First of all, we propose a simple and computationally inexpensive model for the description of the sea bed displacement during an underwater earthquake, based on the finite fault solution for the slip distribution under some assumptions on the dynamics of the rupturing process. Once the bottom motion is reconstructed, we study waves induced on the free surface of the ocean. For this purpose we consider three different models approximating the Euler equations of the water wave theory. Namely, we use the linearized Euler equations (we are in fact solving the Cauchy-Poisson problem), a Boussinesq system and a novel weakly nonlinear model. An intercomparison of these approaches is performed. The developments of the present study are illustrated on the 17 July 2006 Java event, where an underwater earthquake of magnitude 7.7 generated a tsunami that inundated the southern coast of Java.Comment: 31 pages, 10 figures, 3 tables. Other author's papers can be downloaded at http://www.lama.univ-savoie.fr/~dutykh

    Stochastic representation of the Reynolds transport theorem: revisiting large-scale modeling

    Get PDF
    We explore the potential of a formulation of the Navier-Stokes equations incorporating a random description of the small-scale velocity component. This model, established from a version of the Reynolds transport theorem adapted to a stochastic representation of the flow, gives rise to a large-scale description of the flow dynamics in which emerges an anisotropic subgrid tensor, reminiscent to the Reynolds stress tensor, together with a drift correction due to an inhomogeneous turbulence. The corresponding subgrid model, which depends on the small scales velocity variance, generalizes the Boussinesq eddy viscosity assumption. However, it is not anymore obtained from an analogy with molecular dissipation but ensues rigorously from the random modeling of the flow. This principle allows us to propose several subgrid models defined directly on the resolved flow component. We assess and compare numerically those models on a standard Green-Taylor vortex flow at Reynolds 1600. The numerical simulations, carried out with an accurate divergence-free scheme, outperform classical large-eddies formulations and provides a simple demonstration of the pertinence of the proposed large-scale modeling

    Residual-based stabilization of the finite element approximation to the acoustic perturbation equations for low Mach number aeroacoustics

    Get PDF
    This is the peer reviewed version of the following article: [Guasch, O., Sánchez-Martín, P., Pont, A., Baiges, J., and Codina, R. (2016) Residual-based stabilization of the finite element approximation to the acoustic perturbation equations for low Mach number aeroacoustics. Int. J. Numer. Meth. Fluids, 82: 839–857. doi: 10.1002/fld.4243], which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/fld.4243/abstract. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving.The acoustic perturbation equations (APE) are suitable to predict aerodynamic noise in the presence of a non-uniform mean flow. As for any hybrid computational aeroacoustics approach, a first computational fluid dynamics simulation is carried out from which the mean flow characteristics and acoustic sources are obtained. In a second step, the APE are solved to get the acoustic pressure and particle velocity fields. However, resorting to the finite element method (FEM) for that purpose is not straightforward. Whereas mixed finite elements satisfying an appropriate inf–sup compatibility condition can be built in the case of no mean flow, that is, for the standard wave equation in mixed form, these are difficult to implement and their good performance is yet to be checked for more complex wave operators. As a consequence, strong simplifying assumptions are usually considered when solving the APE with FEM. It is possible to avoid them by resorting to stabilized formulations. In this work, a residual-based stabilized FEM is presented for the APE at low Mach numbers, which allows one to deal with the APE convective and reaction terms in its full extent. The key of the approach resides in the design of the matrix of stabilization parameters. The performance of the formulation and the contributions of the different terms in the equations are tested for an acoustic pulse propagating in sheared-solenoidal mean flow, and for the aeolian tone generated by flow past a two-dimensional cylinder.Peer ReviewedPostprint (author's final draft
    corecore