Separation of viscothermal losses and scattering in ultrasonic characterization of porous media

This paper presents a method for separating viscothermal and scattering losses in ultrasonic characterization of porous media. This method is based on variations of the static pressure of the saturating fluid. Experimental results were already presented in previous papers and the losses separation was verified experimentally. The aim of this paper is to present an analytic justification of this losses separation in the case of this experimental method and to show that it is possible to estimate acoustic parameters without the knowledge of scattering characteristics. The standard scattering length is used to renormalize speed and transmission through the porous medium, described as an equivalent fluid. Under certain assumptions corresponding to a weak scattering regime, it is shown how viscothermal and scattering losses can be separated easily without knowing scattering characteristics. Application of this model is presented in the case of weak scattering in a polyurethane foam and in the limit case of stronger scattering in a glass beads sample.     

Scattering of Millimeter Waves from Acid Rain

A theoretical investigation of scattering characteristics of millimeter waves from acid rain is conducted using the method of vector spherical function expansion. The particles of acid rain are considered as the double-layered spheroid. The method is used to derive an analytical solution to the problem of scattering of millimeter waves by a coated dielectric spheroid. The normalized bistatic cross section versus scattering angle for the acid raindrop is calculated at frequencies of 35 GHz and 94 GHz. This is the first part of our study scattering characteristics of millimeter waves from acid rain.     

Quasi-equilibrium closure hierarchies for the Boltzmann equation

In this paper, explicit method of constructing approximations (the triangle entropy method) is developed for nonequilibrium problems. This method enables one to treat any complicated nonlinear functionals that fit best the physics of a problem (such as, for example, rates of processes) as new independent variables.

The work of the method is demonstrated on the Boltzmann's-type kinetics. New macroscopic variables are introduced (moments of the Boltzmann collision integral, or scattering rates). They are treated as independent variables rather than as infinite moment series. This approach gives the complete account of rates of scattering processes. Transport equations for scattering rates are obtained (the second hydrodynamic chain), similar to the usual moment chain (the first hydrodynamic chain). Various examples of the closure of the first, of the second, and of the mixed hydrodynamic chains are considered for the hard sphere model. It is shown, in particular, that the complete account of scattering processes leads to a renormalization of transport coefficients.

The method gives the explicit solution for the closure problem, provides thermodynamic properties of reduced models, and can be applied to any kinetic equation with a thermodynamic Lyapunov function.     

拉曼光谱测定散射介质的温度

本文发展了测量散射体温度的拉曼光谱方法,选用Ｓｒ（ＮＯ３）２粉晶做散射介质,观察到７３８ｃｍ＾－１,１０５７ｃｍ＾－１两个拉曼散射带,记录它们的Ｓｔｏｋｅｓ、Ａｎｔｉ－Ｓｔｏｋｅｓ分量的光谱强度,得到硝酸锶粉晶散射体的温度。实验结果显示,这种非接触光谱测温法具有重要的实用价值。     

A Monte Carlo approach for the calculation of polarized light: application to an incident narrow beam

Monte Carlo approaches to compute multiple scattering of polarized light are examined. A Backward Monte Carlo (BMC) method is developed to solve the Stokes vector of the multiple scattered light for an inhomogeneous scattering medium with boundaries. A generalized form of the BMC method in vector notation is proposed. This method can determine the scattered light with sufficient accuracy in both intensity and polarization compared to the same calculation using the doubling-adding method for a plane parallel medium.For application to a narrow incident beam and an inhomogeneous medium, a modified BMC method is developed, borrowing a concept from the Forward Monte Carlo (FMC) method for the first scattering events. Furthermore, a modification of the total scattering matrix, i.e., the combination of the derived scattering matrix with its time inverse, is discussed. This BMC method can be used successfully for model calculations of lidar and other laser measurements of polarized light.     

In-plane bistatic backward scattering from seabottom with randomly inhomogeneous sediment and rough interface

1 Introduction Backward scattering of sound due to sediment is the main source of shallow waterreverberation. In order to predict the reverberation or detect sediment properties frommeasured reverberation data, a reasonable in-plane bistatic backward scattering (BBS)model is essential. The scattering can be caused by the roughness of water-sediment in-terface or by inhomogeneities within the volume of sediment. A great deal of researchhas been done on sediment backscattering, most of which h…     

Modeling scattering from azimuthally symmetric bathymetric features using wavefield superposition

In this paper, an approach for modeling the scattering from azimuthally symmetric bathymetric features is described. These features are useful models for small mounds and indentations on the seafloor at high frequencies and seamounts, shoals, and basins at low frequencies. A bathymetric feature can be considered as a compact closed region, with the same sound speed and density as one of the surrounding media. Using this approach, a number of numerical methods appropriate for a partially buried target or facet problem can be applied. This paper considers the use of wavefield superposition and because of the azimuthal symmetry, the three-dimensional solution to the scattering problem can be expressed as a Fourier sum of solutions to a set of two-dimensional scattering problems. In the case where the surrounding two half spaces have only a density contrast, a semianalytic coupled mode solution is derived. This provides a benchmark solution to scattering from a class of penetrable hemispherical bosses or indentations. The details and problems of the numerical implementation of the wavefield superposition method are described. Example computations using the method for a simple scattering feature on a seabed are presented for a wide band of frequencies.     

Open rotor tone scattering

This paper describes a method for calculating the scattering of open rotor tones by an adjacent surface such as the aircraft fuselage or wing. This is of practical importance because shielding by these structures could be used to significantly reduce community noise levels. In the first part of this paper the propeller is modelled as a rotating point source and the effect of propeller translation is neglected. The scattering of sound from this rotating point source by either a rigid circular cylinder or a rigid half-plane is calculated and the reflection/shielding effect is clearly demonstrated. In the second part of this paper, the point source/stationary medium model is extended to include the effect of distributed sources and a moving medium. The extended model is then applied to calculate the effect of scattering of rotor-alone tones by a wing (which is modelled as an infinite half-plane). The effect of wing sweep and the application of the Kutta-condition at the trailing edge is investigated. An approximate method for quickly calculating open rotor tone scattering is also developed and validated.     

Radiative properties of densely packed spheres in semitransparent media: A new geometric optics approach

This contribution presents a new Ray-tracing method for calculating effective radiative properties of densely packed spheres in non-absorbing or semitransparent host medium. The method is restricted to the geometric optic objects and neglects the wave effects. The effective radiative properties such as the absorption and scattering coefficients, and phase function are retrieved from the calculation of mean-free paths of scattering and absorption, and the angular scattering probability of radiation propagating in the dispersed medium. The model accounts for the two geometric effects called here as non-point scattering and ray transportation effects. The successful comparison of the current model with data of radiative properties and transmittances of particle beds in a non-absorbing medium reported in the literature confirm its suitability. It is shown that: (i) for opaque or absorbing particles (not systematically opaque), the non-point scattering is the dominant geometric effects whereas both non-point scattering and ray transportation effects occur for weakly absorbing and transparent particles. In the later cases, these two geometric effects oppose and may cancel out. This may explain why the Independent scattering theory works well for packed of quasi-transparent particles; (ii) the non-point scattering and ray transportation effects can be captured through the scattering and absorption coefficients while using the classical form of phase function. This enables using the standard radiative transfer equation (RTE); (iii) the surrounding medium absorption can be accounted for without any homogenization rule. It contributes to increasing the effective absorption coefficient of the composite medium as expected but, at the same time, it reduces the particle extinction; and (iv) the current transfer calculation predicts remarkably the results of direct Monte Carlo (MC) simulation. This study tends therefore to confirm that the RTE can be applied to densely packed media by using effective radiative properties.     

粒子模拟中激励源的一种等效设置方法

 基于计算电磁学中对强迫激励源消除虚假反射的算法分析,提出了用等效电流和等效磁流在FDTD公式中引入电场激励源和磁场激励源的方法。从粒子模拟方法的基本方程和迭代公式出发,分析了激励源的引入过程,推导出激励源所等效的电流项和磁流项表达式,实现了新的激励源设置方法,并进行了数值验证和结果讨论。研究表明：这类等效模型与标准FDTD公式能紧密结合,引入非常方便;不必专门设置一个附加的散射场区来处理散射场的计算,大大节省了计算时间和计算内存,比常规总场/散射场体系方法的效率高20%以上,对粒子模拟这类耗时的计算较为适用。通过对2维柱坐标系下相对论速调管放大器(RKA)的模拟,证明了此类激励源设置方法的实用性。     

Fast and accurate radiance calculations using truncation approximation for anisotropic scattering phase functions

The phase function for solar light scattering by large particles such as cloud droplets is strongly anisotropic due to very strong peaking in the forward direction. This creates numerical difficulties when attempting to calculate accurate reflected and transmitted radiances, which are important for remote sensing of atmospheric and surface properties. A popular approach uses the delta function to approximate the forward-scattering peak in a fraction of energy and a limited number of polynomial terms or a geometrically truncated function for the remaining fraction (so-called truncation approximations). This article compares and discusses several methods for fast and accurate calculations using truncation approximations. When using a single truncation approximation for all scattering orders, large biases appear in directions near the solar and anti-solar points. As shown here, high accuracy can be obtained using different truncation approximations depending on the order of scattering. Of particular importance is the use of phase functions close to the exact phase functions for the first few orders of scattering. Applying the method in combination with the Monte Carlo (MC) method, in which the truncation fraction for a scattering order depends on the scattering angle at the previous scattering event, obtains accurate radiance calculations under almost all geometrical and optical conditions, including in directions near the solar point. Because the method also reduces computational noise due to the MC sampling of radiance, it is useful for fast and accurate radiance calculations for cloudy atmospheres.     

周期结构表面散射系数测试及修正方法*

无规入射散射系数是描述结构声反射特性的重要参数,在周期结构等典型扩散体的声学设计和室内声场模拟中具有重要的作用。本文给出了混响室法测试无规入射散射系数的原理,对正弦型周期结构进行了实验测试,并针对测试中存在的缺陷,提出了一种考虑空气间隙吸声的修正方法,经验证具有较好的修正效果。     

Regularization of boundary integral equations in the problems of wave field diffraction by curved boundaries

A regularization of the exact Fredholm integral equations for the field or its derivative on a scattering surface is proposed. This approach allows one to calculate the scattering or diffraction of pulsed wave fields by curved surfaces of arbitrary geometry. Mathematically, the method is based on the replacement of the exact Fredholm integral equations by their truncated analogs, in which the contributions of the geometrically shadowed regions are cancelled. This approach has a clear physical meaning and provides stable solutions even when the direct numerical solution of mathematically exact initial integral equations leads to unstable results. The method is mathematically substantiated and tested using the problem of plane-wave scattering by a cylinder as an example.     

Theoretical analysis and measurements on near specular scattering

In scattering measurement it is important and also difficult to separate instrument scattering from the scattering of a sample. Particularly at angles near the specular reflection, it is believed that the small angle limit for scattering measurement is determined not only by detector proximity to the specular beam, but also by the sample scattering. When the scattering angles are smaller than a certain limit, the total scattering combining the instrument scattering and the sample scattering is lower than the instrument scattering. This prevents us from discriminating the instrument scattering and the sample scattering at angles smaller than this limit. In this paper a data processing method is reported that enables us to distinguish the instrument scattering at angles smaller than this limit. Based on this method, we can measure the scattering at angles restricted only by the performance of the instrument.     

Finite difference computation of acoustic scattering by small surface inhomogeneities and discontinuities

The use of finite difference schemes to compute the scattering of acoustic waves by surfaces made up of different materials with sharp surface discontinuities at the joints would, invariably, result in the generations of spurious reflected waves of numerical origin. Spurious scattered waves are produced even if a high-order scheme capable of resolving and supporting the propagation of the incident wave is used. This problem is of practical importance in jet engine duct acoustic computation. In this work, the basic reason for the generation of spurious numerical waves is first examined. It is known that when the governing partial differential equations of acoustics are discretized, one should only use the long waves of the computational scheme to represent or simulate the physical waves. The short waves of the computational scheme have entirely different propagation characteristics. They are the spurious numerical waves. A method by which high wave number components (short waves) in the wave scattering process is intentionally removed so as to minimize the scattering of spurious numerical waves is proposed. This method is implemented in several examples from computational aeroacoustics to illustrate its effectiveness, accuracy and efficiency. This method is also employed to compute the scattering of acoustic waves by scatterers, such as rigid wall acoustic liner splices, with width smaller than the computational mesh size. Good results are obtained when comparing with computed results using much smaller mesh size. The method is further extended for applications to computations of acoustic wave reflection and scattering by very small surface inhomogeneities with simple geometries.     

Three-dimensional visualization of objects in scattering medium using integral imaging and spectral analysis

This paper presents a three-dimensional visualization method of 3D objects in a scattering medium. The proposed method employs integral imaging and spectral analysis to improve the visual quality of 3D images. The images observed from 3D objects in the scattering medium such as turbid water suffer from image degradation due to scattering. The main reason is that the observed image signal is very weak compared with the scattering signal. Common image enhancement techniques including histogram equalization and contrast enhancement works improperly to overcome the problem. Thus, integral imaging that enables to integrate the weak signals from multiple images was discussed to improve image quality. In this paper, we apply spectral analysis to an integral imaging system such as the computational integral imaging reconstruction. Also, we introduce a signal model with a visibility parameter to analyze the scattering signal. The proposed method based on spectral analysis efficiently estimates the original signal and it is applied to elemental images. The visibility-enhanced elemental images are then used to reconstruct 3D images using a computational integral imaging reconstruction algorithm. To evaluate the proposed method, we perform the optical experiments for 3D objects in turbid water. The experimental results indicate that the proposed method outperforms the existing methods.     

Simple method for obtaining electron scattering phase shifts from energies of an atom in a cavity

We present a simple method for obtaining elastic scattering phase shifts and cross sections from precise ab initio many-body perturbation theory energies of atoms in variable cavities. This method does not require calculations of wave functions of continuum states, can be generalized to many atoms and ions, and is extremely convenient because existing codes developed for energy calculations can be used without modification. The high precision of the method and close agreement with experiment are illustrated on examples of e-Ar and e-Kr scattering. Correlations as well as relativistic corrections are systematically considered.     

Diffuse surface calibration method to improve accuracy and dynamic range of aerosol elastic light scattering measurements

A new method to calibrate detectors for elastic light scattering (ELS) measurement based on diffuse scattering from a Lambertian surface is presented. The method produces a calibration signal that is approximately seven orders of magnitude larger than a propane gas Rayleigh scattering calibration. The method also allows for calibration of detectors such as photodiodes, which are not sensitive enough to detect Rayleigh scattering for calibration but possess characteristics desirable for the measurement of soot ELS. Since the method is only suitable for backward scattering calibrations, transfer of calibration data from a backward- to a forward-oriented detector is accomplished with a secondary laser and integrating sphere. In demonstration experiments, calibration constants for photomultiplier tube (PMT) detectors obtained using both Rayleigh scattering and diffuse surface scattering agreed within experimental uncertainties as did measurements of in-flame scattering coefficients obtained with PMTs and photodiodes. However, achievable uncertainties with the diffuse-surface calibration approach were significantly reduced. More importantly, by enabling the use of photodiode detectors in ELS measurements, the new method facilitates operation at higher photon fluxes resulting in improved signal-to-noise ratios, reduced influence of photon shot noise, and the ability to achieve higher dynamic range in transient measurements.     

Reconstruction of object location for diffuse fluorescence tomography on the basis of hybrid models of light scattering in biotissues

We propose a new method for reconstructing the spatial distribution of fluorophore in a highly scattering object from its images obtained by the method of diffuse fluorescence tomography. This method is intended for diagnostics of the fluorophore-marked tumors and is based on the algebraic-reconstruction principle combined with a new theoretical model and simulation of light propagation in randomly scattering media, such as biotissues, by the Monte Carlo method. The model experiments show that for 18-mm thick objects, the developed method allows one to determine location of the geometric center of a fluorescent inhomogeneity and its transverse and longitudinal dimensions with accuracies of up to 0.5 and 1.5 mm, respectively.