A VARIATIONAL EXPECTATION-MAXIMIZATION METHOD FOR THE INVERSE BLACK BODY RADIATION PROBLEM

The inverse black body radiation problem, which is to reconstruct the area temperature distribution from the measurement of power spectrum distribution, is a well-known ill-posed problem. In this paper, a variational expectation-maximization （EM） method is developed and its convergence is studied. Numerical experiments demonstrate that the variational EM method is more efficient and accurate than the traditional methods, including the Tikhonov regularization method, the Landweber method and the conjugate gradient method.

THE COUPLING METHOD OF FINITE ELEMENTS AND BOUNDARY ELEMENTS FOR RADIATION PROBLEM

The paper presents the variational formulation and well posedness of the coupling method offinite elements and boundary elements for radiation problem. The convergence and optimal errorestimate for the approximate solution and numerical experiment are provided.

The enhanced volume source boundary point method for the calculation of acoustic radiation problem

The Volume Source Boundary Point Method (VSBPM) is greatly improved so that it will speed up the VSBPM's solution of the acoustic radiation problem caused by the vibrating body. The fundamental solution provided by Helmholtz equation is enforced in a weighted residual sense over a tetrahedron located on the normal line of the boundary node to replace the coefficient matrices of the system equation. Through the enhanced volume source boundary point analysis of various examples and the sound field of a vibrating rectangular box in a semi-anechoic chamber, it has revealed that the calculating speed of the EVSBPM is more than 10 times faster than that of the VSBPM while it works on the aspects of its calculating precision and stability, adaptation to geometric shape of vibrating body as well as its ability to overcome the non-uniqueness problem.

Numerical Approximation of a Nonlinear 3D Heat Radiation Problem

In this paper,we are concerned with the numerical approximation of a steady-state heat radiation problem with a nonlinear Stefan-Boltzmann boundary condition in IR^(3).We first derive an equivalent minimization problem and then present a finite element analysis to the solution of such a minimization problem.Moreover,we apply the Newton iterative method for solving the nonlinear equation resulting from the minimization problem.A numerical example is given to illustrate theoretical results.

WATER SURFACE WAVE RADIATION GENERATED BY MULTIPLE CYLINDERS OSCILLATING WITH IDENTICAL FREQUENCY

The water surface wave radiation problem caused by multiple cylinders oscillating with identical frequency was solved in frequency domain by the boundary element method using simple Green's function in the inner water region combined with the eigenfunction expansions in the outer water region. The numerical method is suitable to the situation of constant depth of outer regions and complicated boundary conditions of inner region, while the oscillating modes, motion amplitudes and phases of the cylinders may be different from one another. The second order potential and hydrodynamic forces acting on each cylinder were evaluated completely by perturbation method. Compared with the case of single oscillating cylinder, hydrodynamic interference phenomena, such as wave resonance and negative added mass, of the radiation problem due to the oscillatory motions of multiple cylinders are identified which is of engineering importance to the design of moorings and other facilities involving multiple structures.

Simultaneous estimation of aerosol optical constants and size distribution from angular light-scattering measurement signals

The angular light-scattering measurement（ALSM） method combined with an improved artificial bee colony algorithm is introduced to determine aerosol optical constants and aerosol size distribution（ASD） simultaneously. Meanwhile, an optimized selection principle of the ALSM signals based on the sensitivity analysis and principle component analysis（PCA）is proposed to improve the accuracy of the retrieval results. The sensitivity analysis of the ALSM signals to the optical constants or characteristic parameters in the ASD is studied first to find the optimized selection region of measurement angles. Then, the PCA is adopted to select the optimized measurement angles within the optimized selection region obtained by sensitivity analysis. The investigation reveals that, compared with random selection measurement angles, the optimized selection measurement angles can provide more useful measurement information to ensure the retrieval accuracy. Finally,the aerosol optical constants and the ASDs are reconstructed simultaneously. The results show that the retrieval accuracy of refractive indices is better than that of absorption indices, while the characteristic parameters in ASDs have similar retrieval accuracy. Moreover, the retrieval accuracy in studying L-N distribution is a little better than that in studying Gamma distribution for the difference of corresponding correlation coefficient matrixes of the ALSM signals. All the results confirm that the proposed technique is an effective and reliable technique in estimating the aerosol optical constants and ASD simultaneously.

Asymptotic-Preserving Discrete Schemes for Non-Equilibrium Radiation Diffusion Problem in Spherical and Cylindrical Symmetrical Geometries

We study the asymptotic-preserving fully discrete schemes for nonequilibrium radiation diffusion problem in spherical and cylindrical symmetric geometry.The research is based on two-temperature models with Larsen’s flux-limited diffusion operators.Finite volume spatially discrete schemes are developed to circumvent the singularity at the origin and the polar axis and assure local conservation.Asymmetric second order accurate spatial approximation is utilized instead of the traditional first order one for boundary flux-limiters to consummate the schemes with higher order global consistency errors.The harmonic average approach in spherical geometry is analyzed,and its second order accuracy is demonstrated.By formal analysis,we prove these schemes and their corresponding fully discrete schemes with implicitly balanced and linearly implicit time evolutions have first order asymptoticpreserving properties.By designing associated manufactured solutions and reference solutions,we verify the desired performance of the fully discrete schemes with numerical tests,which illustrates quantitatively they are first order asymptotic-preserving and basically second order accurate,hence competent for simulations of both equilibrium and non-equilibrium radiation diffusion problems.

ON VELOCITY POTENTIALS DUE TO PULSATING PRESSURE DISTRIBUTIONS ON THE FREE SURFACE OF INFINITE-DEPTH WATERS AND THE RADIATION CONDITIONS FOR SECOND-ORDER DIFFRACTION PROBLEMS

In the present paper,two-and three-dimensional velocity potentials generated by pulsating pressure distributions of infinite extent on the free surface of infinite-depth waters are strictly derived based on special cases of concentrated pulsating pressure.The far-field asymptotic behaviour of the potentials and the radiation conditions to be satisfied by them are discussed. It is proved in a general sense that the potentials should be composed of a forced wave component,a free wave component and a local disturbance component.The radiation condition of the forced wave component should correspond to the far-field asymptotic behaviour of the pressure distribution,Hence,the formulation of radiation conditions for the second-order diffraction potentials has theoretically become clear,The radiation conditions for two-and three-dimensional problems are explicitly given in the paper.

An application of radiation boundary condition to scattering problem of acoustic waves in fluids

A numerical method of solving acoustic wave scattering pnblem in fluids is described. Radiation boundary condition (RBC) obtained by factorization method of Helmholtz equation is applied to transforming the exterior boundary value problem in unbounded region into one in a finite region. Combined with RBC and scatterer surface boundary condition, Helmholtz equation is solved numerically by the finite difference method. Computational results for sphere and prolate spheroidal scatterers are in excellent agreement with eigenfunction solutions and much better than the results of OSRC method.

CALCULATION OF NONLINEAR FREE－SURFACE FLOWS RESULTING FROM LARGE￣AMPLITUDE FORCED OSCILLATION OF A TWO－DIMENSIONAL BODY

The BEM combined with the time- st6pping scheme has been applied to the numerical calculation of fully nonlinear free surface flows generated by large amplitude forced transverse oscillation of two-dimensional body. Particular attention is paid on the compatibility of free surface and body surface conditions at the intersection point, and moving radiation boundary is adopted. A new calculating formula of the exact force on the body is also presented.The results demonstrate some nonlinear phenomena and indicate the stability and correctness of the numerical simulation.

SPLITTING A CONCAVE DOMAIN TO CONVEX SUBDOMAINS

We examine a steady-state heat radiation problem and its finite element approximation in Rd, d = 2, 3. A nonlinear Stefan-Boltzmann boundary condition is considered. Another nonlinearity is due to the fact that the temperature is always greater or equal than 0[K]. We prove two convergence theorems for piecewise linear finite element solutions.