Synchronous exploration for the control of real-time external noise suppression in a three-dimensional subdomain

A mathematical (difference) model is proposed for a real-time active shielding device that shields an acoustic field in a given subdomain from the influence of sound sources located in an additional subdomain. An algorithm for computing the current control ensuring a prescribed process is based on information produced by the author’s technique of synchronous weak noise exploration. This information can be measured in real time. Active control problems for nonstationary solutions of linear difference equations in a three-dimensional domain consisting of two subdomains are studied using the difference potential method. The shape of the domain and the boundary conditions may depend on time, while the coefficients may depend on time and spatial coordinates. If the difference problem is a mathematical model of sound propagation, the goal of control is to change the acoustic field in the given subdomains, for example, to shield the acoustic field in one subdomain from the undesirable influence (noise) of sources located in the other subdomain.     

Mathematical modelling of noise around a Y-shaped windy road

In this study, we investigated the sound radiation from vehicles moving on a road with a Y-shaped fork in windy conditions. First, we consider the problem for two point sound sources moving in opposite directions along a road with an arbitrary bend. The background of the road is modelled as an elastic half-space. The solution to this problem is obtained using the integral Fourier transforms over space variables and time, and the stationary phase method. Next, the solution to the general case for many sources moving along a Y-shaped road fork is obtained as the superposition of the partial solutions. We performed a numerical analysis of the traffic noise characteristics for a Y-shaped fork on a city road in the town of Lodz, Poland as an example.     

Non-stationary problem of active sound control in bounded domains

The present paper deals with the non-stationary problem of active shielding of a domain from undesirable external sources of noise. Active shielding is achieved by constructing additional (secondary) sources in such a way that the total contribution of all sources leads to the desirable effect. The problem is formulated as an inverse source problem with the secondary sources positioned outside the domain to be shielded. Along with the undesirable field (noise) to be shielded the presence of a desirable component (“friendly” sound) is accepted in the analysis. The constructed solution of the problem requires only the knowledge of the total field (noise) on the perimeter of the shielded domain. Some important aspects of the problem are addressed in the paper for the first time, such as the non-stationary formulation of the problem, existence of the resonance regimes and sensitivity of the solution to the input errors. The obtained solution is applicable to aeroacoustics problems.     

Near-field fluctuations and far-field noise of a three-element airfoil system by a discrete vortex method

A circulation-based discrete vortex method is used on a three-element airfoil system. Kutta conditions and Kelvin’s circulation theorem are additional conditions required for this method to determine the circulation distributions on each element and to determine vortex shedding. Discrete shed vortices are introduced near the four sharp edges to represent the sharp-edge vortex shedding caused by unsteady flow separation. The computational procedure warrants neutrally stable solutions of the self-sustained fluctuating flowfield that can provide broad-band spectral information for far-field noise predictions. The near-field vortex method directly calculates the parameters used in an asymptotic formula for far-field sound computation that attributes the noise sources to vortex interactions among the shed vortices and the surface circulations of the three-element airfoil system. The far-field noise characters are then analyzed and compared to the experimental data in the literature.     

Nonlinear problem of active sound control

We consider the problem of active sound control, in which some domain is protected from the field generated outside. The active shielding is realized via the implementation of additional sources in such a way that the total contribution of all sources leads to the wanted effect. Mathematically the problem is reduced to seeking the source terms satisfying some a priori described requirements and belongs to the class of inverse source problems. From the application standpoint, this problem can be closely related to active noise shielding and active vibration. It is important that along with unwanted field to be shielded a wanted field is accepted in the analysis. The solution to the problem requires only the knowledge of the total field at the perimeter of the shielded domain. For the first time the active shielding sources are obtained for the nonlinear statement of the problem. It is obtained via the theory of potentials, and the solution is represented in the form of a simple layer. For this purpose, the theory of Calderón-Ryaben’kii potentials is first extended to nonlinear formulations. In the solution, we also take into account the feedback of the secondary sources on the input data.     

Lavrentiev regularization + Ritz approximation = uniform finite element error estimates for differential equations with rough coefficients

We consider a parametric family of boundary value problems for a diffusion equation with a diffusion coefficient equal to a small constant in a subdomain. Such problems are not uniformly well-posed when the constant gets small. However, in a series of papers, Bakhvalov and Knyazev have suggested a natural splitting of the problem into two well-posed problems. Using this idea, we prove a uniform finite element error estimate for our model problem in the standard parameter-independent Sobolev norm. We also study uniform regularity of the transmission problem, needed for approximation. A traditional finite element method with only one additional assumption, namely, that the boundary of the subdomain with the small coefficient does not cut any finite element, is considered.

One interpretation of our main theorem is in terms of regularization. Our FEM problem can be viewed as resulting from a Lavrentiev regularization and a Ritz-Galerkin approximation of a symmetric ill-posed problem. Our error estimate can then be used to find an optimal regularization parameter together with the optimal dimension of the approximation subspace.

     

On the impact of thermal sources on noise generation

Sound emission is nowadays considered as a major environmental issue. The sound emission is generated, amongst other sources, due to an increasing amount of traffic and transport, i.e. road transport, rail and air traffic. Here, sound emission presents a significant risk to public health and a major cause of stress, especially in industrial countries. In this framework the present work is dedicated to the topic of active noise control with the target of noise reduction. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A spatial domain decomposition method for parabolic optimal control problems

We present a non-overlapping spatial domain decomposition method for the solution of linear–quadratic parabolic optimal control problems. The spatial domain is decomposed into non-overlapping subdomains. The original parabolic optimal control problem is decomposed into smaller problems posed on space–time cylinder subdomains with auxiliary state and adjoint variables imposed as Dirichlet boundary conditions on the space–time interface boundary. The subdomain problems are coupled through Robin transmission conditions. This leads to a Schur complement equation in which the unknowns are the auxiliary state adjoint variables on the space-time interface boundary. The Schur complement operator is the sum of space–time subdomain Schur complement operators. The application of these subdomain Schur complement operators is equivalent to the solution of an subdomain parabolic optimal control problem. The subdomain Schur complement operators are shown to be invertible and the application of their inverses is equivalent to the solution of a related subdomain parabolic optimal control problem. We introduce a new family of Neumann–Neumann type preconditioners for the Schur complement system including several different coarse grid corrections. We compare the numerical performance of our preconditioners with an alternative approach recently introduced by Benamou.     

A transient eddy current problem on a moving domain. Numerical analysis

The aim of this paper is to introduce and analyze a numerical method to solve a transient eddy current problem which arises from the modeling of electromagnetic forming in the axisymmetric case. The resulting problem is degenerate parabolic with the time derivative acting on a moving subdomain. This paper is the sequel of Bermúdez et al. (SIAM J. Math. Anal. 45, 3629–3650, 2013), where a weak formulation of this problem was proved to be well posed and additional regularity of the solution was also established. In the present paper, we propose a finite element method in space combined with a backward Euler time scheme for its numerical solution. We obtain error estimates and report numerical results which allow us to assess the performance of the proposed method.     

Effective optimization-based approach for designing active noise control systems in enclosures

An efficient constrained optimization-based prototype program OPTANC is developed to expedite the optimum design of active noise control systems in enclosures. The boundary element method is used to model the sound field of enclosures in which the walls provide complex impedance and point noise sources may exist at arbitrary locations. A sequential quadratic programming algorithm is selected as the optimizer for the deisgn because of its accuracy, efficiency, and reliability. The program is coded in C with portability on micro, mini, and mainframe computers, and is also modularized for future expansion. Simulations show that the software can effectively and efficiently produce the optimal locations and sound strengths of the control sources for active noise control problems.     

A PRACTICAL PARALLEL DIFFERENCE SCHEME FOR PARABOLIC EQUATIONS

A practical parallel difference scheme for parabolic equations is constructed as follows: to decompose the domain Ω into some overlapping subdomains, take flux of the last time layer as Neumann boundary conditions for the time layer on inner boundary points of subdomains, solve it with the fully implicit scheme on each subdomain, then take correspondent values of its neighbor subdomains as its values for inner boundary points of each subdomain and mean of its neighbor subdomain and itself at overlapping points. The scheme is unconditionally convergent. Though its truncation error is O(τ h), the convergent order for the solution can be improved to O(τ h2).     

On the convergence rate of a parallel nonoverlapping domain decomposition method

In recent years,a nonoverlapping domain decomposition iterative procedure,which is based on using Robin-type boundary conditions as information transmission conditions on the subdomain interfaces,has been developed and analyzed.It is known that the convergence rate of this method is 1-O(h),where h is mesh size.In this paper,the convergence rate is improved to be 1-O(h1/2 H-1/2)sometime by choosing suitable parameter,where H is the subdomain size.Counter examples are constructed to show that our convergence estimates are sharp,which means that the convergence rate cannot be better than 1-O(h1/2H-1/2)in a certain case no matter how parameter is chosen.     

Mathematical simulation of acoustic wave refraction near a caustic

Issues related to the computation of wave fields in an acoustic medium near caustics are considered. A boundary condition on a caustic is established, and the Green’s function of a boundary value problem for the general case of a varying speed of sound is constructed. For this purpose, an auxiliary Goursat problem is considered and a system of its particular solutions is constructed using hyper-geometric functions. A Volterra integral equation for the Green’s function is obtained, and an algorithm for its expansion with respect to smoothness is described. A finite difference scheme approximating the solution of the differential problem with an unbounded coefficient is proposed. Numerical results are presented.     

Statistical analysis of stochastic resonance with ergodic diffusion noise

A subthreshold signal is transmitted through a channel and may be detected when some noise--with known structure and proportional to some level--is added to the data. There is an optimal noise level, called stochastic resonance that corresponds to the highest Fisher information in the problem of estimation of the signal. As noise we consider an ergodic diffusion process and the asymptotic is considered as time goes to infinity. We propose consistent estimators of the subthreshold signal and we solve further a problem of hypotheses testing. We also discuss evidence of stochastic resonance for both estimation and hypotheses testing problems via examples.     

Train sound radiation

The problem of sound radiation from a train is considered. This object has been simulated as sets of point sources uniformly distributed in the domain of a moving lengthened rectangle. The solution of the problem is obtained using integral Fourier transforms over the space coordinates and time. The integrals are calculated by using the stationary phase method. Numerical analysis is carried out for acoustic pressure and sound intensity.     

Homogenization of high frequency nonlinear acoustics equations

High frequency nonlinear acoustics equations proposed in [1] are considered. The acoustic characteristics of the medium rapidly oscillate. This model describes the propagation of pulses of sound shocks generated by supersonic planes, blast waves in atmosphere and ocean, continuous radiation of sound sources. An asymptotic solutions of high frequency nonlinear acoustics equations is constructed (when the characteristic size of inhomogeneity is small with respect to the height of the layer where the problem is posed     

Singular value decomposition in an inverse source problem

An inverse source problem for the wave equation with additional information measured on some parts of the boundary is considered. The degree of ill-posedness of the inverse problem is investigated. A numerical algorithm based on the SVD of a discrete inverse problem is constructed and tested.     

Active sound control in composite regions

In active sound control, noise shielding of a target region is achieved via additional sources (called controls) situated at the perimeter of the region. The sources protect the target region by adjusting the acoustic field near the boundary of the region. In the present paper a numerical model of active sound control based on surface potentials in 3D bounded composite regions is numerically studied. In the composite region setup, it is required that the regions be shielded from noise while allowing admissible sound that is generated in the shielded regions to be preserved. The admissible sound is usually required to propagate freely inside the protected regions or in a (selective) predetermined pattern. The adjusting approach used here does not require any knowledge of the sound sources or the properties of the propagation medium in order to obtain the controls. Moreover, the approach differs sharply from some other approaches where the detailed knowledge of the sound sources and the propagation medium is required. For the first time, numerical test cases involving both free communication and predetermined communication pattern between the regions in three dimensions are considered. In all test cases, these regions are effectively shielded from the noise while any present admissible sound is preserved. In addition, selective propagation of the admissible sound between the regions is enforced. The effect of the number of controls on their operation is also studied. Whether admissible sound is present or not, the level of noise cancellation decreases linearly as fewer controls are used. In addition to the increase in size of the interference zone, the controls become individually distinguishable.     

On the linearized theory op flow around bodies. Method of force sources

The method of force sources is proposed for solving linear problems related to the interaction between rigid bodies, and fluids, or gases. Method is based on the introduction of perturbation force sources into equation of motion of fluid media. Boundary conditions at the rigid body surface make it possible to reduce the problem of hydrodynamic reactions to an integral equation defining the function of force sources. Method is illustrated by the solution of three simple problems in the field of acoustics, and of viscous, and compressible media flow around bodies.

In the linearized theory of flow around rigid bodies, as well as in acoustics, an important part of the sound wave generation analysis concerns the determination of hydrodynamic reactions of the medium on moving, pulsating, or oscillating bodies. Such reactions make themselves felt as constant, or variable mechanical forces, such as drag and lift, or in the case of sound wave emitters, as the wave resistance. Various methods had been proposed for the computation of such forces, as for example, in the monographs [1 to 6].

Here, a different approach to the problem of determination of surface forces exerted by liquids and gases on the rigid body is proposed. By resorting to the formalism of the generalized functions it is possible to introduce into the equations of motion of fluid media a perturbation source in the form volume density of forces exercised by the body on the gas. The distribution of surface tension entering into the expression of this force is selected in such a manner as to satisfy boundary conditions at the body surface. It becomes possible with the use of this device to reduce the problem of determination of forces acting on the body surface to the solution of certain Integral equations. The proposed method is in all respects completely analogous to the well-known method of sources and sinks [1 to 1]. Both methods reduce the problem of interaction between body and gas to the solution of Integral equations. The method of sources and sinks, however, leads to an integral equation which describes the distribution of fictitious sources and sinks in the volume of the body having the density of the medium, while the method of force sources yields an integral equation which directly defines the distribution of mechanical forces over the surface of the body (*).

We may note that the method of force sources had to a certain extent been already used in papers [6 and 7] for the determination of sound radiation by means of point-force sources.     

A two-level domain decomposition method for the iterative solution of high frequency exterior Helmholtz problems

Summary. We present a Lagrange multiplier based two-level domain decomposition method for solving iteratively large-scale systems of equations arising from the finite element discretization of high-frequency exterior Helmholtz problems. The proposed method is essentially an extension of the regularized FETI (Finite Element Tearing and Interconnecting) method to indefinite problems. Its two key ingredients are the regularization of each subdomain matrix by a complex interface lumped mass matrix, and the preconditioning of the interface problem by an auxiliary coarse problem constructed to enforce at each iteration the orthogonality of the residual to a set of carefully chosen planar waves. We show numerically that the proposed method is scalable with respect to the mesh size, the subdomain size, and the wavenumber. We report performance results for a submarine application that highlight the efficiency of the proposed method for the solution of high frequency acoustic scattering problems discretized by finite elements. Received March 17, 1998 / Revised version received June 7, 1999 / Published online January 27, 2000