A finite element approach combining a reduced‐order system,Padé approximants,and an adaptive frequency windowing for fast multi‐frequency solution of poro‐acoustic problems

In this work, a solution strategy is investigated for the resolution of multi‐frequency structural‐acoustic problems including 3D modeling of poroelastic materials. The finite element method is used, together with a combination of a modal‐based reduction of the poroelastic domain and a Padé‐based reconstruction approach. It thus takes advantage of the reduced‐size of the problem while further improving the computational efficiency by limiting the number of frequency resolutions of the full‐sized problem. An adaptive procedure is proposed for the discretization of the frequency range into frequency intervals of reconstructed solution. The validation is presented on a 3D poro‐acoustic example. Copyright © 2013 John Wiley & Sons, Ltd.     

Computed tomography sound velocity reconstruction using incomplete data

An approach based on limited-angle transmission tomography for reconstruction of the sound velocity distribution in the breast is proposed. The imaging setup is similar to that of x-ray mammography. With this setup, the time-of-flight data are acquired by a linear array positioned at the top of the compressed breast that both transmits and receives, and a metal plate is placed at the bottom as a reflector. The setup allows acoustic data acquisition for simultaneous B-mode image formation and the tomographic sound velocity reconstruction. In order to improve the sound velocity estimation accuracy, a new reconstruction algorithm based on a convex programming formulation has been developed. Extensive simulations for both imaging and time-of-flight data based on a 5-MHz linear array were performed on tissues with different geometries and acoustic parameters. Results show that the sound velocity error was generally 1-3 m/s, with a maximum of 5.8 m/s. The radii of the objects under investigation varied from 2 to 6 mm, and all of them were detected successfully. Thus, the proposed approach has been shown to be both feasible and accurate. The approach can be used to complement conventional B-mode imaging to further enhance the detection of breast cancer.     

Topology optimization of acoustic–structure interaction problems using a mixed finite element formulation

The paper presents a gradient‐based topology optimization formulation that allows to solve acoustic–structure (vibro‐acoustic) interaction problems without explicit boundary interface representation. In acoustic–structure interaction problems, the pressure and displacement fields are governed by Helmholtz equation and the elasticity equation, respectively. Normally, the two separate fields are coupled by surface‐coupling integrals, however, such a formulation does not allow for free material re‐distribution in connection with topology optimization schemes since the boundaries are not explicitly given during the optimization process. In this paper we circumvent the explicit boundary representation by using a mixed finite element formulation with displacements and pressure as primary variables (a u /p‐formulation). The Helmholtz equation is obtained as a special case of the mixed formulation for the elastic shear modulus equating to zero. Hence, by spatial variation of the mass density, shear and bulk moduli we are able to solve the coupled problem by the mixed formulation. Using this modelling approach, the topology optimization procedure is simply implemented as a standard density approach. Several two‐dimensional acoustic–structure problems are optimized in order to verify the proposed method. Copyright © 2006 John Wiley & Sons, Ltd.     

层析扫描声学显微镜相关算法及其比较

为了获得物体内部清晰的三维图像,提高激光扫描声学显微镜的分辨能力,在激光扫描声学显微镜原理的基础上,阐述了层析扫描声学显微镜的结构和工作方式。结合层析技术并且根据超声波在介质中的传播规律介绍了三种层析算法:频域内的插值算法;往返传播算法;滤波反向传播算法。最后通过理论和实验比较,结果表明:往返传播算法通过对滤波函数进行改进,在有限孔径的前提下,应用于层析扫描声学显微镜时比其它两种算法可以获得更好的重建图像效果。     

Short-Range Image-Based Method for the Inspection of Strong Scatterers Using Microwaves

In this paper, a new inverse-scattering-based reconstruction method is presented. The aim of the approach is to "invert" wavefleld samples, which are collected by experimental tomographic systems working at radio frequencies and microwaves. Imaging systems operating in this band require the solution of a nonlinear and highly ill-posed inverse problem. The need to regularize the inverse problem is addressed here by considering an efficient inexact Newton method that is able to inspect strong scatterers. In this paper, the mathematical formulation of the approach is detailed and discussed. Moreover, the results of several numerical simulations concerning the reconstructions of dielectric structures in noisy environments and in several applicative scenarios are reported.     

A coupled finite-element, boundary-integral method for simulating ultrasonic flowmeters

Today's most popular technology of ultrasonic flow measurement is based on the transit-time principle. In this paper, a numerical simulation technique applicable to the analysis of transit-time flowmeters is presented. A flowmeter represents a large simulation problem that also requires computation of acoustic fields in moving media. For this purpose, a novel boundary integral method, the Helmholtz integral-ray tracing method (HIRM), is derived and validated. HIRM is applicable to acoustic radiation problems in arbitrary mean flows at low Mach numbers and significantly reduces the memory demands in comparison with the finite-element method (FEM). It relies on an approximate free-space Green's function which makes use of the ray tracing technique. For simulation of practical acoustic devices, a hybrid simulation scheme consisting of FEM and HIRM is proposed. The coupling of FEM and HIRM is facilitated by means of absorbing boundaries in combination with a new, reflection-free, acoustic-source formulation. Using the coupled FEM-HIRM scheme, a full three-dimensional (3-D) simulation of a complete transit-time flowmeter is performed for the first time. The obtained simulation results are in good agreement with measurements both at zero flow and under flow conditions.     

A generalized two point ellipsoidal anisotropic ray tracing for converted waves

Rocks can be anisotropic due to a variety of reasons. When estimating rock velocities from seismic data, failure to introduce anisotropy into earth models could generate distortions in the final images that can have enormous economic impact. To estimate anisotropic earth velocities by tomographic methods, it is necessary to trace rays or to solve the wave equation in models where anisotropy has been properly considered. Thus, in this work we present a 3-D generalized ellipsoidal travel time formulation that allow us to trace rays in an anisotropic medium. We propose to trace rays in anisotropic media by solving a set of nonlinear optimization problems, where the group velocities for P and S wave propagation modes are 3-D ellipsoidal approximations that have been recently obtained. Moreover, we prove that this 3-D ellipsoidal anisotropic ray tracing formulation is a convex nonlinear optimization problem, and therefore any solution of the problem is a global minimum. Each optimization problem is solved by the global spectral gradient method, which requires first order information and has low computation and low storage requirements. Our approach for tracing rays in anisotropic media is a generalization in the sense that handles titled axis of symmetry and, close to the axis of symmetry, it is an accurate formulation for 2-D transversely isotropic media and 3-D orthorhombic media, depending on the input parameters. Moreover, this formulation gives the exact ray trajectories in 2-D and 3-D homogeneous isotropic media. The simplicity of the formulation and the low computational cost of the optimization method allow us to present a variety of numerical results that illustrate the behavior and computational advantages of the approach, and the difficulties when working in anisotropic media. Partially supported by Fonacit project UCV-97-003769     

声学数值计算的分区光滑径向点插值无网格法

针对标准的有限元法分析声学问题时由于数值色散导致高波数计算结果不可靠问题,将分区光滑径向点插值法(cell-based smoothed radial point interpolation method, CS-RPIM)应用到二维声学分析中,推导了分区光滑径向点插值法分析二维声学问题的原理公式。该方法将问题域划分为三角形背景单元,每个单元进一步分成若干个光滑域,对每个光滑域进行声压梯度光滑处理,运用光滑Galerkin弱形式构造系统方程,并按有限元中方法施加必要的边界条件。CS-RPIM提供了合适的模型硬度,能有效降低色散效应,提高计算精度。对管道和二维轿车声学问题的数值分析结果表明,与标准有限元法相比,CS-RPIM具有更高的精度和准确度,在高波数计算时这种优势特别明显。     

A symmetric weak form of Biot's equations based on redundant variables representing the fluid,using a Helmholtz decomposition of the fluid displacement vector field

A novel symmetric weak formulation of Biot's equations for linear acoustic wave propagation in layered poroelastic media is presented. The primary variables used are the frame displacement, the acoustic pore pressure, the scalar potential and the vector potential obtained from a Helmholtz decomposition of the fluid displacement. Also a symmetric weak form based on the frame displacement, the pore pressure and the fluid displacement is obtained as an intermediate result. hp finite element simulations of a double leaf partition based on this new weak formulation is verified against simulation results from the classical frame displacement, fluid displacement formulation and a frame displacement pore pressure formulation. All three formulations simulated, displays the same rate of convergence with respect to finite element bases polynomial degree. The novel formulation also extends a previously published frame displacement, pore pressure, scalar fluid displacement potential formulation with an implicit irrotational fluid displacement assumption to a full representation of Biot's equations. Copyright © 2010 John Wiley & Sons, Ltd.     

Elasto–acoustic and acoustic–acoustic coupling on non‐matching grids

Flexible discretization techniques for the approximative solution of coupled wave propagation problems are investigated, focussing on aero–acoustic and elasto–acoustic coupling. In particular, the advantages of using non‐matching grids are presented, when one subregion has to be resolved by a substantially finer grid than the other subregion. For the elasto–acoustic coupling, the problem formulation remains essentially the same as for the matching situation, while for the aero–acoustic coupling, the formulation is enhanced with Lagrange multipliers within the framework of mortar finite element methods. Several numerical examples are presented to demonstrate the flexibility and applicability of the approach. Copyright © 2006 John Wiley & Sons, Ltd.     

Three-dimensional tomographic reconstruction of mesospheric airglow structures using two-station ground-based image measurements

A new methodology is presented to create two-dimensional (2D) and three-dimensional (3D) tomographic reconstructions of mesospheric airglow layer structure using two-station all-sky image measurements. A fanning technique is presented that produces a series of cross-sectional 2D reconstructions, which are combined to create a 3D mapping of the airglow volume. The imaging configuration is discussed and the inherent challenges of using limited-angle data in tomographic reconstructions have been analyzed using artificially generated imaging objects. An iterative reconstruction method, the partially constrained algebraic reconstruction technique (PCART), was used in conjunction with a priori information of the airglow emission profile to constrain the height of the imaged region, thereby reducing the indeterminacy of the inverse problem. Synthetic projection data were acquired from the imaging objects and the forward problem to validate the tomographic method and to demonstrate the ability of this technique to accurately reconstruct information using only two ground-based sites. Reconstructions of the OH airglow layer were created using data recorded by all-sky CCD cameras located at Bear Lake Observatory, Utah, and at Star Valley, Wyoming, with an optimal site separation of ~100 km. The ability to extend powerful 2D and 3D tomographic methods to two-station ground-based measurements offers obvious practical advantages for new measurement programs. The importance and applications of mesospheric tomographic reconstructions in airglow studies, as well as the need for future measurements and continued development of techniques of this type, are discussed.     

Formulation of a thin shell finite element with continuity C 0 and convected material frame notion

 This paper presents a formulation for a new family of thin shell finite elements. The element is formulated by using a convected material frame notion which offers an interesting framework to take into account large transformations. Bending behaviour is calculated from the Love–Kirchhoff assumptions and from a finite difference technique between adjacent elements. We therefore called this element SFE for semi-finite-element. This method allows us to keep C ⁰ continuity without introducing other variables than the 3 classical displacements, which reduces computational time. In this paper, a full formulation of this element is described more precisely. It takes into account the coupling effect between both membrane and bending behaviour. Various sample solutions that illustrate the effectiveness of the element in linear and nonlinear analysis are presented, with some sheet metal forming examples. Received 10 January 2001     

A fast BE‐FE coupling scheme for partly immersed bodies

Fluid–structure coupled problems are investigated to predict the vibro‐acoustic behavior of submerged bodies. The finite element method is applied for the structural part, whereas the boundary element method is used for the fluid domain. The focus of this paper is on partly immersed bodies. The fluid problem is favorably modeled by a half‐space formulation. This way, the Dirichlet boundary condition on the free fluid surface is incorporated by a half‐space fundamental solution. A fast multipole implementation is presented for the half‐space problem. In case of a high density of the fluid, the forces due to the acoustic pressure, which act on the structure, cannot be neglected. Thus, a strong coupling scheme is applied. An iterative solver is used to handle the coupled system. The efficiency of the proposed approach is discussed using a realistic model problem. Copyright © 2009 John Wiley & Sons, Ltd.     

An Improved Method for Ray Tracing Through Curved Inhomogeneities in Composite Materials

Ultrasonic NDE of anisotropic materials, due to its inherent complexity, demands a stringent approach to determine wave propagation path in presence of inhomogeneities. In addition to developing inspection procedures, the ray path information needs to be integrated with the tomographic reconstruction algorithm of composite materials. The present study proposes a method to identify inhomogeneity boundaries using edge detection approach and to implement the same in two-dimensional ray tracing model. In general, inhomogeneity boundaries can be detected on the basis of first-order gradient of a field value (such as Young's modulus). The boundaries can be represented in terms of location and orientation. This information is included in the algorithm for more accurate ray tracing. Examples show the usefulness of the method to take into effect poor penetration of ultrasound energy in certain regions of inhomogeneous materials.     

Isogeometric FEM‐BEM coupled structural‐acoustic analysis of shells using subdivision surfaces

We introduce a coupled finite and boundary element formulation for acoustic scattering analysis over thin‐shell structures. A triangular Loop subdivision surface discretisation is used for both geometry and analysis fields. The Kirchhoff‐Love shell equation is discretised with the finite element method and the Helmholtz equation for the acoustic field with the boundary element method. The use of the boundary element formulation allows the elegant handling of infinite domains and precludes the need for volumetric meshing. In the present work, the subdivision control meshes for the shell displacements and the acoustic pressures have the same resolution. The corresponding smooth subdivision basis functions have the C¹ continuity property required for the Kirchhoff‐Love formulation and are highly efficient for the acoustic field computations. We verify the proposed isogeometric formulation through a closed‐form solution of acoustic scattering over a thin‐shell sphere. Furthermore, we demonstrate the ability of the proposed approach to handle complex geometries with arbitrary topology that provides an integrated isogeometric design and analysis workflow for coupled structural‐acoustic analysis of shells.     

BOUNDARY ELEMENT–FINITE ELEMENT COUPLED EIGENANALYSIS OF FLUID–STRUCTURE SYSTEMS

The eigenanalysis of acoustical cavities with flexible structure boundaries, such as a fluid-filled container or an automobile cabin enclosure, is considered. An algebraic eigenvalue problem formulation for the fluid–structure problem is presented by combining the acoustic fluid boundary element eigenvalue analysis method and the structural finite elements. For many practical eigenproblems, use of finite elements to discretize the fluid domain leads to large stiffness and mass matrices. Since the acoustic boundary element discretization requires putting nodes only on the wetted surface of the structure, the size of the eigenproblem is reduced considerably, thus reducing the eigenvalue extraction effort. Futhermore, unlike in ordinary cases, the finite element discretization of pressure–displacement based fluid–structure problem gives rise to unsymmetric matrices. Therefore, the fact that the boundary element formulation produces unsymmetric matrices does not introduce additional difficulties here compared to the finite element case in the choice of an eigenvalue extraction procedure. Examples are included to demonstrate the fluid–structure eigenanalysis using boundary elements for the fluid domain and finite elements for the structure.     

Optical phase contrast measurement of ultrasonic fields

This report describes an optical phase contrast imaging technique for the measurement of wide bandwidth ultrasound fields in water. In this method, a collimated optical wavefront (λ_l = 810 nm) impinges on a wide bandwidth ultrasound pulse. The method requires that refractive index perturbations induced by the ultrasound field be sufficiently small. Specifically, on exit from the acoustic field, the phase of the optical wavefront must be proportional to the ray sum of local density taken in the direction of propagation of the incident optical wave. A similar restriction is placed on the dimensions of the ultrasound pulse. Repeated measurement of this phase as the ultrasound field is rotated through 180° about an axis normal to the direction of propagation of the incident optical wave generates the Radon transform of the ultrasonically induced refractive index perturbation. Standard tomographic reconstruction techniques are used to reconstruct the full three-dimensional refractive index perturbation. A simple two-lens imaging system and an optical signal processing element from phase contrast microscopy provide a method of directly measuring an affine function of the desired optical phase for small optical phase shifts. The piezo- and elasto-optic coefficients (the first partial derivatives of refractive index with respect to density and pressure) relate refractive index to density and pressure via a linear model. The optical measurement method described in this paper provides a direct, quantitative measurement of the piezo- and elasto-optic coefficients (from the density or pressure fields)     

Optimal wave-front reconstruction strategies for multiconjugate adaptive optics

We propose an optimal approach for the phase reconstruction in a large field of view (FOV) for multiconjugate adaptive optics. This optimal approach is based on a minimum-mean-square-error estimator that minimizes the mean residual phase variance in the FOV of interest. It accounts for the C2n profile in order to optimally estimate the correction wave front to be applied to each deformable mirror (DM). This optimal approach also accounts for the fact that the number of DMs will always be smaller than the number of turbulent layers, since the C2n profile is a continuous function of the altitude h. Links between this optimal approach and a tomographic reconstruction of the turbulence volume are established. In particular, it is shown that the optimal approach consists of a full tomographic reconstruction of the turbulence volume followed by a projection onto the DMs accounting for the considered FOV of interest. The case where the turbulent layers are assumed to match the mirror positions [model-approximation (MA) approach], which might be a crude approximation, is also considered for comparison. This MA approach will rely on the notion of equivalent turbulent layers. A comparison between the optimal and MA approaches is proposed. It is shown that the optimal approach provides very good performance even with a small number of DMs (typically, one or two). For instance, good Strehl ratios (greater than 20%) are obtained for a 4-m telescope on a 150-arc sec x 150-arc sec FOV by using only three guide stars and two DMs.     

Perfectly matched layers for frequency-domain integral equation acoustic scattering problems

Simulations of acoustic wavefields in inhomogeneous media are always performed on finite numerical domains. If contrasts actually extend over the domain boundaries of the numerical volume, unwanted, non-physical reflections from the boundaries will occur. One technique to suppress these reflections is to attenuate them in a locally reflectionless absorbing boundary layer enclosing the spatial computational domain, a perfectly matched layer (PML). This technique is commonly applied in time-domain simulation methods like finite element methods or finite-difference time-domain, but has not been applied to the integral equation method. In this paper, a PML formulation for the three-dimensional frequency-domain integral-equation-based acoustic scattering problem is derived. Three-dimensional acoustic scattering configurations are used to test the PML formulation. The results demonstrate that strong attenuation (a factor of 200 in amplitude) of the scattered pressure field is achieved for thin layers with a thickness of less than a wavelength, and that the PMLs themselves are virtually reflectionless. In addition, it is shown that the integral equation method, both with and without PMLs, accurately reproduces pressure fields by comparing the obtained results with analytical solutions.     

Solution of FE‐BE coupled eigenvalue problems for the prediction of the vibro‐acoustic behavior of ship‐like structures

To predict the vibro‐acoustic behavior of structures, both a structural problem and an acoustic problem have to be solved. For thin structures immersed in water, a strong interaction between the structural domain and fluid domain occurs. This significantly alters the resonance frequencies. In this work, the structure is modeled by the finite element method. The exterior acoustic problem is solved by a fast boundary element method employing hierarchical matrices. An FE‐BE formulation is presented, which allows the solution of the coupled eigenvalue problem and thus the prediction of the coupled eigenfrequencies and mode shapes. It is based on a Schur complement formulation of the FE‐BE system yielding a generalized eigenvalue problem. A Krylov–Schur solver is applied for its efficient solution. Hereby, the compressibility of the fluid is neglected. The coupled eigensolution is then used for a model reduction strategy allowing fast frequency sweep calculations. The efficiency of the proposed formulations is investigated with respect to memory consumption, accuracy, and computation time. Copyright © 2011 John Wiley & Sons, Ltd.