Eulerian-Lagrangian method for simulation of cloud cavitation

We present a coupled Eulerian-Lagrangian method to simulate cloud cavitation in a compressible liquid. The method is designed to capture the strong, volumetric oscillations of each bubble and the bubble-scattered acoustics. The dynamics of the bubbly mixture is formulated using volume-averaged equations of motion. The continuous phase is discretized on an Eulerian grid and integrated using a high-order, finite-volume weighted essentially non-oscillatory (WENO) scheme, while the gas phase is modeled as spherical, Lagrangian point-bubbles at the sub-grid scale, each of whose radial evolution is tracked by solving the Keller-Miksis equation. The volume of bubbles is mapped onto the Eulerian grid as the void fraction by using a regularization (smearing) kernel. In the most general case, where the bubble distribution is arbitrary, three-dimensional Cartesian grids are used for spatial discretization. In order to reduce the computational cost for problems possessing translational or rotational homogeneities, we spatially average the governing equations along the direction of symmetry and discretize the continuous phase on two-dimensional or axi-symmetric grids, respectively. We specify a regularization kernel that maps the three-dimensional distribution of bubbles onto the field of an averaged two-dimensional or axi-symmetric void fraction. A closure is developed to model the pressure fluctuations at the sub-grid scale as synthetic noise. For the examples considered here, modeling the sub-grid pressure fluctuations as white noise agrees a priori with computed distributions from three-dimensional simulations, and suffices, a posteriori, to accurately reproduce the statistics of the bubble dynamics. The numerical method and its verification are described by considering test cases of the dynamics of a single bubble and cloud cavitaiton induced by ultrasound fields.Comment: 28 pages, 16 figure

Signal-Theoretic Characterization of Waveguide Mesh Geometries for Models of Two--Dimensional Wave Propagation in Elastic Media

Waveguide Meshes are efficient and versatile models of wave propagation along a multidimensional ideal medium. The choice of the mesh geometry affects both the computational cost and the accuracy of simulations. In this paper, we focus on 2D geometries and use multidimensional sampling theory to compare the square, triangular, and hexagonal meshes in terms of sampling efficiency and dispersion error under conditions of critical sampling. The analysis shows that the triangular geometry exhibits the most desirable tradeoff between accuracy and computational cost.Comment: 9 pages, 6 figures, 1 table, to appear on IEEE Transactions on Speech and Audio Processing, vol. 9, no. 2, february 200

A Space-Time Discontinuous Galerkin Trefftz Method for time dependent Maxwell's equations

We consider the discretization of electromagnetic wave propagation problems by a discontinuous Galerkin Method based on Trefftz polynomials. This method fits into an abstract framework for space-time discontinuous Galerkin methods for which we can prove consistency, stability, and energy dissipation without the need to completely specify the approximation spaces in detail. Any method of such a general form results in an implicit time-stepping scheme with some basic stability properties. For the local approximation on each space-time element, we then consider Trefftz polynomials, i.e., the subspace of polynomials that satisfy Maxwell's equations exactly on the respective element. We present an explicit construction of a basis for the local Trefftz spaces in two and three dimensions and summarize some of their basic properties. Using local properties of the Trefftz polynomials, we can establish the well-posedness of the resulting discontinuous Galerkin Trefftz method. Consistency, stability, and energy dissipation then follow immediately from the results about the abstract framework. The method proposed in this paper therefore shares many of the advantages of more standard discontinuous Galerkin methods, while at the same time, it yields a substantial reduction in the number of degrees of freedom and the cost for assembling. These benefits and the spectral convergence of the scheme are demonstrated in numerical tests

Wave modelling - the state of the art

This paper is the product of the wave modelling community and it tries to make a picture of the present situation in this branch of science, exploring the previous and the most recent results and looking ahead towards the solution of the problems we presently face. Both theory and applications are considered. The many faces of the subject imply separate discussions. This is reflected into the single sections, seven of them, each dealing with a specific topic, the whole providing a broad and solid overview of the present state of the art. After an introduction framing the problem and the approach we followed, we deal in sequence with the following subjects: (Section) 2, generation by wind; 3, nonlinear interactions in deep water; 4, white-capping dissipation; 5, nonlinear interactions in shallow water; 6, dissipation at the sea bottom; 7, wave propagation; 8, numerics. The two final sections, 9 and 10, summarize the present situation from a general point of view and try to look at the future developments