Kernel Interpolation of Incident Sound Field in Region Including Scattering Objects

A method for estimating the incident sound field inside a region containing scattering objects is proposed. The sound field estimation method has various applications, such as spatial audio capturing and spatial active noise control; however, most existing methods do not take into account the presence of scatterers within the target estimation region. Although several techniques exist that employ knowledge or measurements of the properties of the scattering objects, it is usually difficult to obtain them precisely in advance, and their properties may change during the estimation process. Our proposed method is based on the kernel ridge regression of the incident field, with a separation from the scattering field represented by a spherical wave function expansion, thus eliminating the need for prior modeling or measurements of the scatterers. Moreover, we introduce a weighting matrix to induce smoothness of the scattering field in the angular direction, which alleviates the effect of the truncation order of the expansion coefficients on the estimation accuracy. Experimental results indicate that the proposed method achieves a higher level of estimation accuracy than the kernel ridge regression without separation.Comment: Accepted to IEEE Workshop on Applications of Signal Processing to Audio and Acoustics (WASPAA) 202

Spatial Active Noise Control Method Based On Sound Field Interpolation From Reference Microphone Signals

A spatial active noise control (ANC) method based on the interpolation of a sound field from reference microphone signals is proposed. In most current spatial ANC methods, a sufficient number of error microphones are required to reduce noise over the target region because the sound field is estimated from error microphone signals. However, in practical applications, it is preferable that the number of error microphones is as small as possible to keep a space in the target region for ANC users. We propose to interpolate the sound field from reference microphones, which are normally placed outside the target region, instead of the error microphones. We derive a fixed filter for spatial noise reduction on the basis of the kernel ridge regression for sound field interpolation. Furthermore, to compensate for estimation errors, we combine the proposed fixed filter with multichannel ANC based on a transition of the control filter using the error microphone signals. Numerical experimental results indicate that regional noise can be sufficiently reduced by the proposed methods even when the number of error microphones is particularly small.Comment: Accepted to International Conference on Acoustics, Speech and Signal Processing (ICASSP) 202

Numerical Computation, Data Analysis and Software in Mathematics and Engineering

The present book contains 14 articles that were accepted for publication in the Special Issue “Numerical Computation, Data Analysis and Software in Mathematics and Engineering” of the MDPI journal Mathematics. The topics of these articles include the aspects of the meshless method, numerical simulation, mathematical models, deep learning and data analysis. Meshless methods, such as the improved element-free Galerkin method, the dimension-splitting, interpolating, moving, least-squares method, the dimension-splitting, generalized, interpolating, element-free Galerkin method and the improved interpolating, complex variable, element-free Galerkin method, are presented. Some complicated problems, such as tge cold roll-forming process, ceramsite compound insulation block, crack propagation and heavy-haul railway tunnel with defects, are numerically analyzed. Mathematical models, such as the lattice hydrodynamic model, extended car-following model and smart helmet-based PLS-BPNN error compensation model, are proposed. The use of the deep learning approach to predict the mechanical properties of single-network hydrogel is presented, and data analysis for land leasing is discussed. This book will be interesting and useful for those working in the meshless method, numerical simulation, mathematical model, deep learning and data analysis fields

Inverse problems in thermoacoustics

Thermoacoustics is a branch of fluid mechanics, and is as such governed by the conservation laws of mass, momentum, energy and species. While computational fluid dynamics (CFD) has entered the design process of many applications in fluid mechanics, its success in thermoacoustics is limited by the multi-scale, multi-physics nature of the subject. In his influential monograph from 2006, Prof. Fred Culick writes about the role of CFD in thermoacoustic modeling: The main reason that CFD has otherwise been relatively helpless in this subject is that problems of combustion instabilities involve physical and chemical matters that are still not well understood. Moreover, they exist in practical circumstances which are not readily approximated by models suitable to formulation within CFD. Hence, the methods discussed and developed in this book will likely be useful for a long time to come, in both research and practice. [. . . ] It seems to me that eventually the most effective ways of formulating predictions and theoretical interpretations of combustion instabilities in practice will rest on combining methods of the sort discussed in this book with computational fluid dynamics, the whole confirmed by experimental results. Despite advances in CFD and large-eddy simulation (LES) in particular, unsteady simulations for more than a few selected operating points are computationally infeasible. The ‘methods discussed in this book’ refer to reduced-order models of thermoacoustic oscillations. Whether intentional or not, the last sentence anticipates the advent of data-driven methods, and encapsulates the philosophy behind this work. This work brings together two workhorses of the design process: physics-informed reduced-order models and data from higher-fidelity sources such as simulations and experiments. The three building blocks to all our statistical inference frameworks are: (i) a hierarchical view of reduced-order models consisting of states, parameters and governing equations; (ii) probabilistic formulations with random variables and stochastic processes; and (iii) efficient algorithms from statistical learning theory and machine learning. While leveraging advances in statistical and machine learning, we demonstrate the feasibility of Bayes’ rule as a first principle in physics-informed statistical inference. In particular, we discuss two types of inverse problems in thermoacoustics: (i) implicit reduced-order models representative of nonlinear eigenproblems from linear stability analysis; and (ii) time-dependent reduced-order models used to investigate nonlinear dynamics. The outcomes of statistical inference are improved predictions of the state, estimates of the parameters with uncertainty quantification and an assessment of the reduced-order model itself. This work highlights the role that data can play in the future of combustion modeling for thermoacoustics. It is increasingly impractical to store data, particularly as experiments become automated and numerical simulations become more detailed. Rather than store the data itself, the techniques in this work optimally assimilate the data into the parameters of a physics-informed reduced-order model. With data-driven reduced-order models, rapid prototyping of combustion systems can feed into rapid calibration of their reduced-order models and then into gradient-based design optimization. While it has been shown, e.g. in the context of ignition and extinction, that large-eddy simulations become quantitatively predictive when augmented with data, the reduced-order modeling of flame dynamics in turbulent flows remains challenging. For these challenging situations, this work opens up new possibilities for the development of reduced-order models that adaptively change any time that data from experiments or simulations becomes available.Schlumberger Cambridge International Scholarshi

Computational automation for efficient design of acoustic metamaterials

Acoustic metamaterials (AMMs) are an exciting technology because they are capable of responding to vibrations in ways that are impossible to achieve with conventional materials. However, realization of AMMs requires engineering design to provide a connection between first-principles research and production of parts that perform as expected. Designing AMMs is a challenging endeavor because evaluating designs is costly and manufacturing metamaterials requires precise techniques with small minimum resolutions. To address these challenges, new computational tools are necessary to aid design. This work proposes three tasks that improve the capabilities of design for AMM while being extensible to other engineering design automation tasks. The first task is to develop a design exploration tool that improves the computational efficiency of identifying sets of high-performing designs in a design space that is sparse and comprises mixed discrete/continuous data. The second task is to develop a process for designers to evaluate manufacturability of difficult-to-manufacture parts and drive co-development of manufacturing methods and AMM. In the final task, a machine learning based method is developed to efficiently model AMM with heterogeneous arrangements of their microstructures such that strict homogenization is infeasible. The outcomes from completing these tasks will provide a significant and novel improvement over existing methods of designing AMMs.Mechanical Engineerin

Solar Seismology from Space. a Conference at Snowmass, Colorado

The quality of the ground based observing environment suffers from several degrading factors: diurnal interruptions and thermal variations, atmospheric seeing and transparency fluctuations and adverse weather interruptions are among the chief difficulties. The limited fraction of the solar surface observable from only one vantage point is also a potential limitation to the quality of the data available without going to space. Primary conference goals were to discuss in depth the scientific return from current observations and analyses of solar oscillations, to discuss the instrumental and site requirements for realizing the full potential of the seismic analysis method, and to help bring new workers into the field by collecting and summarizing the key background theory. At the conclusion of the conference there was a clear consensus that ground based observation would not be able to provide data of the quality required to permit a substantial analysis of the solar convection zone dynamics or to permit a full deduction of the solar interior structure

Advances in modeling gas adsorption in porous materials for the characterization applications

The dissertation studies methods for mesoporous materials characterization using adsorption at various levels of scale and complexity. It starts with the topic introduction, necessary notations and definitions, recognized standards, and a literature review. Synthesis of novel materials requires tailoring of the characterization methods and their thorough testing. The second chapter presents a nitrogen adsorption characterization study for silica colloidal crystals (synthetic opals). These materials have cage-like pores in the range of tens of nanometers. The adsorption model can be described within a macroscopic approach, based on the Derjaguin-Broekhoff-de Boer (DBdB) theory of capillary condensation. A kernel of theoretical isotherms is built and applied to the solution of the adsorption integral equation to derive the pore-size distribution from experimental data. The technique is validated with a surface modification of the samples so that it changes the interaction but not the pore size. The second chapter deals with the characterization of three-dimensional ordered mesoporous (3DOm) carbons. Similar to opals, these materials have cage-like mesopores, however, these pores are connected with large windows. These windows affect the adsorption process and calculated pore-size distributions. The grand canonical Monte Carlo simulations with derived solid-fluid potentials, which take into account the 3DOm carbons geometry, confirm the critical role of interconnections, their size, and number, for correct interpretation of adsorption data for the PSD calculations. The fourth chapter discusses a method for the pore size estimation that can serve as an alternative to the adsorption isotherms analysis. It is based on measurements of elastic properties of liquid that can be useful for the pore size estimation. A Vycor glass sample, a disordered mesoporous material with channel-like pores having a characteristic size of ca. 6-8 nm, is considered. The changes in longitudinal and shear moduli from the experimental data and molecular simulations are predicted with a near-quantitative agreement. Then, it follows by their relation of the moduli to the pore size, which is promising for characterization. The last fifth chapter considers a promising Monte Carlo method, the Kinetic Monte Carlo (kMC) algorithm. This method is efficient for the vapor-liquid equilibrium prediction in dense regions. This chapter shows a benchmark with conventional Metropolis et al algorithms as well as a parallelization scheme of the kMC algorithm

Snapshot-Based Methods and Algorithms

An increasing complexity of models used to predict real-world systems leads to the need for algorithms to replace complex models with far simpler ones, while preserving the accuracy of the predictions. This two-volume handbook covers methods as well as applications. This second volume focuses on applications in engineering, biomedical engineering, computational physics and computer science

Nonparametric enrichment in computational and biological representations of distributions

This thesis proposes nonparametric techniques to enhance unsupervised learning methods in computational or biological contexts. Representations of intractable distributions and their relevant statistics are enhanced by nonparametric components trained to handle challenging estimation problems. The first part introduces a generic algorithm for learning generative latent variable models. In contrast to traditional variational learning, no representation for the intractable posterior distributions are computed, making it agnostic to the model structure and the support of latent variables. Kernel ridge regression is used to consistently estimate the gradient for learning. In many unsupervised tasks, this approach outperforms advanced alternatives based on the expectation-maximisation algorithm and variational approximate inference. In the second part, I train a model of data known as the kernel exponential family density. The kernel, used to describe smooth functions, is augmented by a parametric component trained using an efficient meta-learning procedure; meta-learning prevents overfitting as would occur using conventional routines. After training, the contours of the kernel become adaptive to the local geometry of the underlying density. Compared to maximum-likelihood learning, our method better captures the shape of the density, which is the desired quantity in many downstream applications. The final part sees how nonparametric ideas contribute to understanding uncertainty computation in the brain. First, I show that neural networks can learn to represent uncertainty using the distributed distributional code (DDC), a representation similar to the nonparametric kernel mean embedding. I then derive several DDC-based message-passing algorithms, including computations of filtering and real-time smoothing. The latter is a common neural computation embodied in many postdictive phenomena of perception in multiple modalities. The main idea behind these algorithms is least-squares regression, where the training data are simulated from an internal model. The internal model can be concurrently updated to follow the statistics in sensory stimuli, enabling adaptive inference