Simulation of sound propagation over porous barriers of arbitrary shapes

A time-domain solver using an immersed boundary method is investigated for simulating sound propagation over porous and rigid barriers of arbitrary shapes. In this study, acoustic propagation in the air from an impulse source over the ground is considered as a model problem. The linearized Euler equations are solved for sound propagation in the air and the Zwikker-Kosten equations for propagation in barriers as well as in the ground. In comparison to the analytical solutions, the numerical scheme is validated for the cases of a single rigid barrier with different shapes and for two rigid triangular barriers. Sound propagations around barriers with different porous materials are then simulated and discussed. The results show that the simulation is able to capture the sound propagation behaviors accurately around both rigid and porous barriers

TIME DOMAIN SIMULATION FOR SOUND PROPAGATION OVER VARIOUS OBJECTS AND UNDER VORTICAL BACKGROUND CONDITIONS

Acoustic wave propagations have been studied for a long time with both experimental and numerical methods. Most of the analytical solutions for wave propagations are considered for simple environments such as a homogeneous atmospheres. As a result, the analytical solutions are unable to be applied for complicated environments. Numerical methods have become more and more important in acoustics studies after decades of development. The finite difference time-domain method (FDTD) is one of the most commonly used numerical methods in wave propagation studies. Compared with the other methods, the FDTD method is able to include many aspects of sound wave behaviors such as reflection, refraction, and diffraction in the physical problems. In this thesis, the linearized acoustic Euler equations coupled with the immersed boundary method are applied to investigate the sound wave propagation over complex environments. For the three-dimensional simulations of sound wave propagation in long distance, the moving domain method and parallel computing techniques are applied. Based on these approaches, the computational costs are significantly reduced and the simulation efficiency is greatly improved. When looking into the effects of high subsonic vortical flow, a high order WENO scheme is applied for the simulation. In this way the simulation stability can be achieved and the sound scattering of vortical flow can be studied. Then, the numerical scheme is applied to simulate an ultrasonic plane wave propagating through biological tissue. The linearized Euler acoustic equations coupled with the spatial fractional Laplacian operators are used for numerical simulations. The absorption and attenuation effects of the biological lossy media are successfully observed from the simulation results. Throughout this thesis, the simulation results are compared with either experimental measurements or analytical solutions so that the accuracy of the implemented numerical scheme is validated

TIME DOMAIN SIMULATION FOR SOUND PROPAGATION OVER VARIOUS OBJECTS AND UNDER VORTICAL BACKGROUND CONDITIONS

Acoustic wave propagations have been studied for a long time with both experimental and numerical methods. Most of the analytical solutions for wave propagations are considered for simple environments such as a homogeneous atmospheres. As a result, the analytical solutions are unable to be applied for complicated environments. Numerical methods have become more and more important in acoustics studies after decades of development. The finite difference time-domain method (FDTD) is one of the most commonly used numerical methods in wave propagation studies. Compared with the other methods, the FDTD method is able to include many aspects of sound wave behaviors such as reflection, refraction, and diffraction in the physical problems. In this thesis, the linearized acoustic Euler equations coupled with the immersed boundary method are applied to investigate the sound wave propagation over complex environments. For the three-dimensional simulations of sound wave propagation in long distance, the moving domain method and parallel computing techniques are applied. Based on these approaches, the computational costs are significantly reduced and the simulation efficiency is greatly improved. When looking into the effects of high subsonic vortical flow, a high order WENO scheme is applied for the simulation. In this way the simulation stability can be achieved and the sound scattering of vortical flow can be studied. Then, the numerical scheme is applied to simulate an ultrasonic plane wave propagating through biological tissue. The linearized Euler acoustic equations coupled with the spatial fractional Laplacian operators are used for numerical simulations. The absorption and attenuation effects of the biological lossy media are successfully observed from the simulation results. Throughout this thesis, the simulation results are compared with either experimental measurements or analytical solutions so that the accuracy of the implemented numerical scheme is validated

Time-Domain Simulation of Sound Propagation in Frequency-Dependent Materials

This dissertation investigates sound propagation in frequency-dependent materials. The study provides an improved understanding of how to numerically model the porous impedance materials more accurately under the conditions of complicated geometries. The finite difference time-domain (FDTD) method is implemented on the linearized Euler equation (LEE), along with the immersed boundary (IB) method and other numerical techniques to simulate the acoustic wave propagation in air, water, porous media and biological tissues. When material properties vary in the frequency domain, their time-domain counterpart may contain either convolution operation or fractional derivative operation. Both operations have been studied in this dissertation. Recursive algorithm methods, piece-wise constant recursive methods (PCRC) and piece-wise linear recursive methods (PLRC) are used to numerically solve for convolution operations, and fractional central difference (FCD) methods are used to solve for fractional Laplacians. Both methods show good results in comparison with analytical solutions. A variety of models have been implemented to simulate the acoustic wave propagation inside porous media. The techniques include: the Zwicker and Kosten (ZK) phenomenological model, the Delany and Bazley model, various porosity two-parameter models, the time-domain boundary condition (TDBC) models, and Wilson’s relaxation model (WRX). A new method is also proposed that utilizes the ANSI/ASA-S1.18 measurements to construct a new relaxation function. The new relaxation function can improve the prediction from the TDBC and WRX models significantly. The ZK and WRX models have also been used in predicting the noise reduction of a house. The noise due to transmission and vibration of the wall is modeled as a simple wave transmission through a porous material layer. A curve fitting method is used to match acoustic properties of the wall material. By assembling all the materials together, the over-all acoustic response of a house can be simulated. When acoustic wave propagating in biological tissues, wave propagation equations were previously solved either with convolutions, which consume a large amount of memory, or with pseudo-spectral methods, which cannot handle complicated geometries effectively. The approach described in this study employs FCD method, combined with the IB method for the FDTD simulation. It also works naturally with the IB method which enables a simple Cartesian-type grid mesh to be used to solve problems with complicated geometries. This work also studies acoustic scattering effects caused by 2D or 3D vortices. The LEE is used to investigate sound wave propagation over subsonic vortices. Instead of traditional direct numerical simulation (DNS) methods, the new approach treats vortex flow field as a scattering background flow and solves the acoustic field with the LEE solver. The numerical method uses a high-order WENO scheme to accommodate the highly convective background flow at high Mach numbers. The study focuses on the acoustic field scaling laws scattered by the 2D and 3D vortices

High fidelity fluid-structure turbulence modeling using an immersed-body method

There is an increasing need for turbulence models with fluid-structure interaction (FSI) in many industrial and environmental high Reynolds number flows. Since the complicated structure boundaries move in turbulent flows, it is quite challenging to numerically apply boundary conditions on these moving fluid-structure interfaces. To achieve accurate and reliable results from numerical FSI simulations in turbulent flows, a high fidelity fluid-structure turbulence model is developed using an immersed-body method in this thesis. It does this by coupling a finite element multiphase fluid model and a combined finite-discrete element solid model via a novel thin shell mesh surrounding solid surfaces. The FSI turbulence model presented has four novelties. Firstly, an unsteady Reynolds-averaged Navier-Stokes (URANS) k−ε turbulence model is coupled with an immersed-body method to model FSI by using this thin shell mesh. Secondly, to reduce the computational cost, a log-law wall function is used within this thin shell to resolve the flow near the boundary layer. Thirdly, in order to improve the accuracy of the wall function, a novel shell mesh external-surface intersection approach is introduced to identify sharp solid-fluid interfaces. Fourthly, the model has been extended to simulate highly compressible gas coupled with fracturing solids. This model has been validated by various test cases and results are in good agreement with both experimental and numerical data in published literature. This model is capable to simulate the aerodynamic and hydrodynamic details of fluids and the stress, vibration, deformation and motion of structures simultaneously. Finally, this model has been applied to several industrial applications including a floating structure being moved around by complex hydrodynamic flows involving wave breaking; a blasting engineering simulation with shock waves, fracture propagation, gas-solid interaction and flying fragments; fluid dynamics, flow-induced vibrations, flow-induced fractures of a full-scale vertical axis turbine. Some useful conclusions, e.g. how to model them, how to make them stable and how to predict when they will break, are obtained by this FSI model when applying it to the above applications.Open Acces

DEVELOPMENT OF A COMPUTATIONAL MODEL FOR A SIMULTANEOUS SIMULATION OF INTERNAL FLOW AND SPRAY BREAK-UP OF THE DIESEL INJECTION PROCESS

El proceso de atomización desde una vena o lámina líquida hasta multitud de gotas dispersas en un medio gaseoso ha sido un fenómeno de interés desde hace varias décadas, especialmente en el campo de los motores de combustión interna alternativos. Multitud de estudios experimentales han sido publicados al respecto, pues una buena mezcla de aire-combustible asegura una evaporación y combustión mucho más eficientes, aumentando la potencia del motor y reduciendo la cantidad de contaminantes emitidos. Con el auge de las técnicas computacionales, muchos modelos han sido desarrollados para estudiar este proceso de atomización y mezcla. Uno de los últimos modelos que han aparecido es el llamado ELSA (Eulerian-Lagrangian Spray Atomization), que utiliza un modelo Euleriano para la parte densa del chorro y cambia a un modelo Lagrangiano cuando la concentración de líquido es suficientemente pequeña, aprovechando de esta manera las ventajas de ambos. En el presente trabajo se ha desarrollado un modelo puramente Euleriano para estudiar la influencia de la geometría interna de la tobera de inyección en el proceso de atomización y mezcla. Se ha estudiado únicamente el proceso de inyección diésel. Este modelo permite resolver en un único dominio el flujo interno y el externo, evitando así las comunes simplificaciones y limitaciones de la interpolación entre ambos dominios resueltos por separado. Los resultados actuales son prometedores, el modelo predice con un error aceptable la penetración del chorro, el flujo másico y de cantidad de movimiento, los perfiles de velocidad y concentración, así como otros parámetros característicos del chorro.Martí Gómez-Aldaraví, P. (2014). DEVELOPMENT OF A COMPUTATIONAL MODEL FOR A SIMULTANEOUS SIMULATION OF INTERNAL FLOW AND SPRAY BREAK-UP OF THE DIESEL INJECTION PROCESS [Tesis doctoral]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/43719TESISPremios Extraordinarios de tesis doctorale

Efficient Light and Sound Propagation in Refractive Media with Analytic Ray Curve Tracer

Refractive media is ubiquitous in the natural world, and light and sound propagation in refractive media leads to characteristic visual and acoustic phenomena. Those phenomena are critical for engineering applications to simulate with high accuracy requirements, and they can add to the perceived realism and sense of immersion for training and entertainment applications. Existing methods can be roughly divided into two categories with regard to their handling of propagation in refractive media; first category of methods makes simplifying assumption about the media or entirely excludes the consideration of refraction in order to achieve efficient propagation, while the second category of methods accommodates refraction but remains computationally expensive. In this dissertation, we present algorithms that achieve efficient and scalable propagation simulation of light and sound in refractive media, handling fully general media and scene configurations. Our approaches are based on ray tracing, which traditionally assumes homogeneous media and rectilinear rays. We replace the rectilinear rays with analytic ray curves as tracing primitives, which represent closed-form trajectory solutions based on assumptions of a locally constant media gradient. For general media profiles, the media can be spatially decomposed into explicit or implicit cells, within which the media gradient can be assumed constant, leading to an analytic ray path within that cell. Ray traversal of the media can therefore proceed in segments of ray curves. The first source of speedup comes from the fact that for smooth media, a locally constant media gradient assumption tends to stay valid for a larger area than the assumption of a locally constant media property. The second source of speedup is the constant-cost intersection computation of the analytic ray curves with planar surfaces. The third source of speedup comes from making the size of each cell and therefore each ray curve segment adaptive to the magnitude of media gradient. Interactions with boundary surfaces in the scene can be efficiently handled within this framework in two alternative approaches. For static scenes, boundary surfaces can be embedded into the explicit mesh of tetrahedral cells, and the mesh can be traversed and the embedded surfaces intersected with by the analytic ray curve in a unified manner. For dynamic scenes, implicit cells are used for media traversal, and boundary surface intersections can be handled separately by constructing hierarchical acceleration structures adapted from rectilinear ray tracer. The efficient handling of boundary surfaces is the fourth source of speedup of our propagation path computation. We demonstrate over two orders-of-magnitude performance improvement of our analytic ray tracing algorithms over prior methods for refractive light and sound propagation. We additionally present a complete sound-propagation simulation solution that matches the path computation efficiency achieved by the ray curve tracer. We develop efficient pressure computation algorithm based on analytic evaluations and combine our algorithm with the Gaussian beam for fast acoustic field computation. We validate the accuracy of the simulation results on published benchmarks, and we show the application of our algorithms on complex and general three-dimensional outdoor scenes. Our algorithms enable simulation scenarios that are simply not feasible with existing methods, and they have the potential of being extended and complementing other propagation methods for capability beyond handling refractive media.Doctor of Philosoph