Multiobjective muffler shape optimization with hybrid acoustics modeling

This paper considers the combined use of a hybrid numerical method for the modeling of acoustic mufflers and a genetic algorithm for multiobjective optimization. The hybrid numerical method provides accurate modeling of sound propagation in uniform waveguides with non-uniform obstructions. It is based on coupling a wave based modal solution in the uniform sections of the waveguide to a finite element solution in the non-uniform component. Finite element method provides flexible modeling of complicated geometries, varying material parameters, and boundary conditions, while the wave based solution leads to accurate treatment of non-reflecting boundaries and straightforward computation of the transmission loss (TL) of the muffler. The goal of optimization is to maximize TL at multiple frequency ranges simultaneously by adjusting chosen shape parameters of the muffler. This task is formulated as a multiobjective optimization problem with the objectives depending on the solution of the simulation model. NSGA-II genetic algorithm is used for solving the multiobjective optimization problem. Genetic algorithms can be easily combined with different simulation methods, and they are not sensitive to the smoothness properties of the objective functions. Numerical experiments demonstrate the accuracy and feasibility of the model-based optimization method in muffler design.     

Cloaks,editors, and bubbles: applications of spacetime transformation theory

Spacetime or ‘event’ cloaking was recently introduced as a concept, and the theoretical design for such a cloak was presented for illumination by electromagnetic waves [McCall et al., J. Opt. 2011]. Here it is described how event cloaks can be designed for simple wave systems, using either an approximate ‘speed cloak’ method, or an exact full‐wave one. Further, details of many of the implications of spacetime transformation devices are discussed, including their (usually) directional nature, spacetime distortions (as opposed to cloaks), and how leaky cloaks manifest themselves. More exotic concepts are also addressed, in particular concepts that follow naturally on from considerations of simple spacetime transformation devices, such as spacetime modeling and causality editors. A proposal for implementing an interrupt‐without‐interrupt concept is described. Finally, the design for a time‐dependent ‘bubbleverse’ is presented, based on temporally modulated Maxwell's Fisheye transformation device (T‐device) in a flat background spacetime.     

Augmented Riemann solvers for the shallow water equations over variable topography with steady states and inundation

We present a class of augmented approximate Riemann solvers for the shallow water equations in the presence of a variable bottom surface. These belong to the class of simple approximate solvers that use a set of propagating jump discontinuities, or waves, to approximate the true Riemann solution. Typically, a simple solver for a system of m conservation laws uses m such discontinuities. We present a four wave solver for use with the the shallow water equations—a system of two equations in one dimension. The solver is based on a decomposition of an augmented solution vector—the depth, momentum as well as momentum flux and bottom surface. By decomposing these four variables into four waves the solver is endowed with several desirable properties simultaneously. This solver is well-balanced: it maintains a large class of steady states by the use of a properly defined steady state wave—a stationary jump discontinuity in the Riemann solution that acts as a source term. The form of this wave is introduced and described in detail. The solver also maintains depth non-negativity and extends naturally to Riemann problems with an initial dry state. These are important properties for applications with steady states and inundation, such as tsunami and flood modeling. Implementing the solver with LeVeque’s wave propagation algorithm [R.J. LeVeque, Wave propagation algorithms for multi-dimensional hyperbolic systems, J. Comput. Phys. 131 (1997) 327–335] is also described. Several numerical simulations are shown, including a test problem for tsunami modeling.     

Simulating Free Surface Flows with SPH

The SPH (smoothed particle hydrodynamics) method is extended to deal with free surface incompressible flows. The method is easy to use, and examples will be given of its application to a breaking dam, a bore, the simulation of a wave maker, and the propagation of waves towards a beach. Arbitrary moving boundaries can be included by modelling the boundaries by particles which repel the fluid particles. The method is explicit, and the time steps are therefore much shorter than required by other less flexible methods, but it is robust and easy to program.     

A Design of Flattop-beam Shaper at Millimeter Wavelengths

Beam shaping technique for generating flattop beam has long been used in optics, and there are many methods to design flattop-beam shapers in optics. Flattop beam is also important in some applications at millimeter wavelengths such as in the quasi-optical power divider, however, there are few design method for the flattop-beam shaper used at millimeter wave frequencies and the design method of optics are not suitable to the beam shaper at millimeter wave frequencies, so a design method for the flattop-beam shaper at millimeter wave frequencies is presented in this paper. Several flattop-beam shapers in the form of diffractive element have been designed and simulation results are presented.     

On the use of Gegenbauer reconstructions for shock wave propagation modeling

In therapeutic ultrasound, the presence of shock waves can be significant due to the use of high intensity beams, as well as due to shock formation during inertial cavitation. Although modeling of such strongly nonlinear waves can be carried out using spectral methods, such calculations are typically considered impractical, since accurate calculations often require hundreds or even thousands of harmonics to be considered, leading to prohibitive computational times. Instead, time-domain algorithms which generally utilize Godunov-type finite-difference schemes are commonly used. Although these time domain methods can accurately model steep shock wave fronts, unlike spectral methods they are inherently unsuitable for modeling realistic tissue dispersion relations. Motivated by the need for a more general model, the use of Gegenbauer reconstructions as a postprocess tool to resolve the band-limitations of the spectral methods are investigated. The present work focuses on eliminating the Gibbs phenomenon when representing a steep wave front using a limited number of harmonics. Both plane wave and axisymmetric 2D transducer problems will be presented to characterize the proposed method.     

CIVA: an expertise platform for simulation and processing NDT data

Ultrasonic modeling and simulation are more and more widely used by the different actors of industrial NDT. The applications are numerous and show a great variety: help for diagnosis, data reconstruction, performance demonstration, probe design and inspection parameters settling, virtual testing etc. The CEA (the French Atomic Energy Commission) is strongly involved in this evolution with the development of the CIVA expertise platform which gathers in the same software advanced processing and modeling tools. In the aim of fulfilling requirements of an intensive use the choice has been made to mainly adopt semi-analytical approximated methods. The wave propagation modeling is based on an integral formulation of the radiated field and applies the so-called pencil method. The modeling of beam-defect interaction and echoes formation mechanisms apply approximated theories such as Kirchhoff approximation or GTD. Over the years and with successive versions of the software, this approach is enriched by adaptations and improvements of the existing models or by new models, in order to extend the field of applicability of the simulation.     

Elastoacoustic model with uncertain mechanical properties for ultrasonic wave velocity prediction: application to cortical bone evaluation

The axial transmission technique can measure the longitudinal wave velocity of an immersed solid. An elementary model of the technique is developed with a set of source and receivers placed in a semi-infinite fluid coupled at a plane interface with a semi-infinite solid. The acoustic fluid is homogeneous. The solid is homogeneous, isotropic, and linearly elastic. The work is focused on the prediction of the measured velocity (apparent velocity) when the solid is considered to have random material properties. The probability density functions of the random variables modeling each mechanical parameter of the solid are derived following the maximum entropy principle. Specific attention is paid to the modeling of Poisson's ratio so that the second-order moments of the velocities remain finite. The stochastic solver is based on a Monte Carlo numerical simulation and uses an exact semianalytic expression of the acoustic response derived with the Cagniard-de Hoop method. Results are presented for a solid with the material properties of cortical bone. The estimated mean values and confidence regions of the apparent velocity are presented for various dispersion levels of the random parameters. A sensibility analysis with respect to the source and receivers locations is presented.     

Real-time synthesis of clarinet-like instruments using digital impedance models

A real-time synthesis model of wind instruments sounds, based upon a classical physical model, is presented. The physical model describes the nonlinear coupling between the resonator and the excitor through the Bernoulli equation. While most synthesis methods use wave variables and their sampled equivalent in order to describe the resonator of the instrument, the synthesis model presented here uses sampled versions of the physical variables all along the synthesis process, and hence constitutes a straightforward digital transposition of each part of the physical model. Moreover, the resolution scheme of the problem (i.e., the synthesis algorithm) is explicit and all the parameters of the algorithm are expressed analytically as functions of the physical and the control parameters.     

基于Matlab Simulink的物理实验——简谐振动仿真研究

以简谐振动和阻尼振动仿真为例介绍了用Simulink对物理实验进行仿真建模的方法。给出了位移、速度等振动曲线;可对振动过程中的动能、势能以及机械能进行监测。体现出Simulink仿真物理实验的优越性。     

大尺寸物光波面彩色数字全息高质量重建研究

为解决检测面尺寸较大时CCD难以得到高质量数字全息图的问题,本文利用负透镜设计光学系统让CCD接收来自物体的缩小虚像,以球面波为参考波,使用单色CCD近距离得到三种色光照射下的大尺寸彩色物体的数字全息图,然后采用可控放大率波面重建算法得到同一尺寸的数字全息重建像,合成彩色数字全息重建像.同时,使用两种消零级方法去除零级干扰,提高重建像质量,一种方法利用空间光调制器相移技术在参考光中加入一次任意相移,记录两幅数字全息图,消除重建零级像|另一种方法使用“无干扰全息图”消除重建零级像及共轭像.本文讨论结果可为大物体彩色数字全息及多波长数字全息检测应用提供有益的参考.     

Full wave modeling of therapeutic ultrasound: efficient time-domain implementation of the frequency power-law attenuation

For the simulation of therapeutic ultrasound applications, a method including frequency-dependent attenuation effects directly in the time domain is highly desirable. This paper describes an efficient numerical time-domain implementation of the power-law attenuation model presented by Szabo [Szabo, J. Acoust. Soc. Am. 96, 491-500 (1994)]. Simulations of therapeutic ultrasound applications are feasible in conjunction with a previously presented finite differences time-domain (FDTD) algorithm for nonlinear ultrasound propagation [Ginter et al., J. Acoust. Soc. Am. 111, 2049-2059 (2002)]. Szabo implemented the empirical frequency power-law attenuation using a causal convolutional operator directly in the time-domain equation. Though a variety of time-domain models has been published in recent years, no efficient numerical implementation has been presented so far for frequency power-law attenuation models. Solving a convolutional integral with standard time-domain techniques requires enormous computational effort and therefore often limits the application of such models to 1D problems. In contrast, the presented method is based on a recursive algorithm and requires only three time levels and a few auxiliary data to approximate the convolutional integral with high accuracy. The simulation results are validated by comparison with analytical solutions and measurements.     

Modeling of acoustic wave propagation in time-domain using the discontinuous Galerkin method – A comparison with measurements

Simulations of acoustic wave propagation in time-domain are presented. In the simulations, the discontinuous Galerkin method for spatial derivatives and the low-storage Runge–Kutta approach for time derivatives are used. Three different simulation cases are studied. First, the directivity of loudspeaker is simulated. In the second case, acoustic wave propagation in free space is studied using a short pulse. In the last case, acoustic wave scattering from a metallic cylinder is simulated. All simulation results are compared with measurement results. The measurements for the acoustic wave scattering from the metallic cylinder are made in 2D planes using an automated measurement system. Comparison between the simulation and measurement results are made both temporally and spatially and a good agreement between the simulation and measurement results is found. The results suggest that the discontinuous Galerkin method coupled with the low-storage Runge–Kutta approach is a viable tool for modeling acoustic wave propagation in the time-domain.     

Method for eliminating zero-order diffraction in lensless Fourier transform digital holography

A method to induce phase shifting in lensless Fourier of digital holography system is presented. In this method, by computer simulation and theoretical analysis on the technology to eliminate the influence of zero-order diffraction in force, it can be found that, a reference light induced in a random non-2π integral number of phase shifting shooting to get the second hologram, and wave reconstruction can be got by the difference value image of the two holograms. And, the method of different digital holography record system with different phase shifting should be used. In this paper, theoretical analyses have been done in detail to discuss the problems that exist in the unsuitable phase-shifting methods. Furthermore, some experiments have been done to prove the reliability of this method. This method can significantly improve the image quality and give better resolution of the reconstructed image.     

基于等阻抗差分段法的变阻抗线电路模拟

新型Z箍缩驱动源系统采用变阻抗线结构形式的径向阻抗变换器进行能量传递和变换,为了获得最大的能量效率,需要应用电路模拟方法研究变阻抗线传递功率的特性。采用了等阻抗差分段法对变阻抗线进行了建模计算,与以往的等长度分段方法比较,在相同的分段数量下,等阻抗差分段法在阻抗变化的连续性上优于等长度分段法。给出了具有指数形式的变阻抗线的算例,分别采用两种分段法进行了计算,对计算结果的比较表明:等阻抗差分段法可以使计算更快地收敛到真实解。     

螺旋线型三导体脉冲形成线快波抵消方法

针对螺旋线型三导体脉冲形成线的快波振荡问题,提出一种阻抗分两段阶梯变化的螺旋线结构,分析了此类螺旋线快波反射与慢波反射相互抵消快波的条件,给出两段变阻抗螺旋线慢波系数之间的关系式以及匝数密度之比计算方法,通过仿真计算证明该方法能有效消除快波反射对螺旋线输出脉冲的影响。     

合成孔径数字全息的分析模拟及多参考光合成孔径数字全息

分析了合成孔径数字全息的理论和实现方法,提出了多参考光合成孔径数字全息方法及其应用前景,以双狭缝为例,模拟了单幅及单参考光、多参考光两种合成孔径数字全息术的记录和再现方法,对三种方法再现像的分辨率进行了比较.     

Digital simulation of random processes and its applications

Efficient methods are presented for digital simulation of a general homogeneous process (multidimensional or multivariate or multivariate-multidimensional) as a series of cosine functions with weighted amplitudes, almost evenly spaced frequencies, and random phase angles. The approach is also extended to the simulation of a general non-homogeneous oscillatory process characterized by an evolutionary power spectrum. Generalized forces involved in the modal analysis of linear or non-linear structures can be efficiently simulated as a multivariate process using the cross-spectral density matrix computed from the spectral density function of the multidimensional excitation process. Possible applications include simulation of (i) wind-induced ocean wave elevation, (ii) spatial random variation of material properties, (iii) the fluctuating part of atmospheric wind velocities and (iv) random surface roughness of highways and airport runways.     

Efficient equalization of multi-exciter distributed mode loudspeakers

An efficient digital equalization method is applied successfully to the problem of spectral equalization of multi-exciter distributed mode loudspeakers (DML). It is based on a chain of second-order sections of infinite impulse response parametric filters with very low computational cost. The method compensates for the measured multi-exciter DML response in order to achieve a desired frequency response. The sound radiation of these flat loudspeakers is a complex superposition of excited modes that vary strongly with frequency. Therefore, the characteristic multi-exciter DML spectrum is very irregular and is equalized with the method presented here for a natural, uncolored response. In multichannel systems, such as wave field synthesis (WFS), the use of efficient filters to equalize a large amount of drivers is an advantageous approach. The equalization process has been applied to two multi-exciter DML prototypes, comprising three and five exciters per panel. Both panel and exciter equalization have been addressed, which consequences on the filtered responses are discussed. Finally, some subjective assessments are carried out to optimize the order of the filter while maintaining the perceived quality of the equalization.     

同轴相移数字全息中相移角的选取及相移误差的消除

基于两步同轴相移数字全息,首先从理论上分析了记录时不同相移角的选取及相移误差对再现像的影响,并给出了一种利用再现像所有抽样点的强度偏差之和作为评价标准,通过逐步改变理论设定相移角值来寻找实际相移角的相移误差消除新方法;其次对二步相移数字全息中记录时参考光波最佳相移角的选取作了计算机模拟,发现只有将参考光波的相移角选择在一定范围内,再现像的噪声较小;最后利用计算机模拟了相移误差消除,验证了所提方法的可行性.