Efficient Computation of Power, Force, and Torque in BEM Scattering Calculations

We present concise, computationally efficient formulas for several quantities of interest -- including absorbed and scattered power, optical force (radiation pressure), and torque -- in scattering calculations performed using the boundary-element method (BEM) [also known as the method of moments (MOM)]. Our formulas compute the quantities of interest \textit{directly} from the BEM surface currents with no need ever to compute the scattered electromagnetic fields. We derive our new formulas and demonstrate their effectiveness by computing power, force, and torque in a number of example geometries. Free, open-source software implementations of our formulas are available for download online

Polarized line formation with J-state interference in the presence of magnetic fields: A heuristic treatment of collisional frequency redistribution

An expression for the partial frequency redistribution (PRD) matrix for line scattering in a two-term atom, which includes the J-state interference between its fine structure line components is derived. The influence of collisions (both elastic and inelastic) and an external magnetic field on the scattering process is taken into account. The lower term is assumed to be unpolarized and infinitely sharp. The linear Zeeman regime in which the Zeeman splitting is much smaller than the fine structure splitting is considered. The inelastic collision rates between the different levels are included in our treatment. We account for the depolarization caused by the collisions coupling the fine structure states of the upper term, but neglect the polarization transfer between the fine structure states. When the fine structure splitting goes to zero, we recover the redistribution matrix that represents the scattering on a two-level atom (which exhibits only m-state interference --- namely the Hanle effect). The way in which the multipolar index of the scattering atom enters into the expression for the redistribution matrix through the collisional branching ratios is discussed. The properties of the redistribution matrix are explored for a single scattering process for an L=0 to 1 to 0 scattering transition with S=1/2 (a hypothetical doublet centered at 5000 A and 5001 A). Further, a method for solving the Hanle radiative transfer equation for a two-term atom in the presence of collisions, PRD, and J-state interference is developed. The Stokes profiles emerging from an isothermal constant property medium are computed.Comment: Accepted for publication in Journal of Quantitative Spectroscopy and Radiative Transfer (JQSRT

Path-tracing Monte Carlo Library for 3D Radiative Transfer in Highly Resolved Cloudy Atmospheres

Interactions between clouds and radiation are at the root of many difficulties in numerically predicting future weather and climate and in retrieving the state of the atmosphere from remote sensing observations. The large range of issues related to these interactions, and in particular to three-dimensional interactions, motivated the development of accurate radiative tools able to compute all types of radiative metrics, from monochromatic, local and directional observables, to integrated energetic quantities. In the continuity of this community effort, we propose here an open-source library for general use in Monte Carlo algorithms. This library is devoted to the acceleration of path-tracing in complex data, typically high-resolution large-domain grounds and clouds. The main algorithmic advances embedded in the library are those related to the construction and traversal of hierarchical grids accelerating the tracing of paths through heterogeneous fields in null-collision (maximum cross-section) algorithms. We show that with these hierarchical grids, the computing time is only weakly sensitivive to the refinement of the volumetric data. The library is tested with a rendering algorithm that produces synthetic images of cloud radiances. Two other examples are given as illustrations, that are respectively used to analyse the transmission of solar radiation under a cloud together with its sensitivity to an optical parameter, and to assess a parametrization of 3D radiative effects of clouds.Comment: Submitted to JAMES, revised and submitted again (this is v2

Thermophysical Phenomena in Metal Additive Manufacturing by Selective Laser Melting: Fundamentals, Modeling, Simulation and Experimentation

Among the many additive manufacturing (AM) processes for metallic materials, selective laser melting (SLM) is arguably the most versatile in terms of its potential to realize complex geometries along with tailored microstructure. However, the complexity of the SLM process, and the need for predictive relation of powder and process parameters to the part properties, demands further development of computational and experimental methods. This review addresses the fundamental physical phenomena of SLM, with a special emphasis on the associated thermal behavior. Simulation and experimental methods are discussed according to three primary categories. First, macroscopic approaches aim to answer questions at the component level and consider for example the determination of residual stresses or dimensional distortion effects prevalent in SLM. Second, mesoscopic approaches focus on the detection of defects such as excessive surface roughness, residual porosity or inclusions that occur at the mesoscopic length scale of individual powder particles. Third, microscopic approaches investigate the metallurgical microstructure evolution resulting from the high temperature gradients and extreme heating and cooling rates induced by the SLM process. Consideration of physical phenomena on all of these three length scales is mandatory to establish the understanding needed to realize high part quality in many applications, and to fully exploit the potential of SLM and related metal AM processes

On the Computation of Power in Volume Integral Equation Formulations

We present simple and stable formulas for computing power (including absorbed/radiated, scattered and extinction power) in current-based volume integral equation formulations. The proposed formulas are given in terms of vector-matrix-vector products of quantities found solely in the associated linear system. In addition to their efficiency, the derived expressions can guarantee the positivity of the computed power. We also discuss the application of Poynting's theorem for the case of sources immersed in dissipative materials. The formulas are validated against results obtained both with analytical and numerical methods for scattering and radiation benchmark cases

Antenna pattern of DUAL detectors of gravitational waves and its exploitation in a network of advanced interferometers

We investigate the directional sensitivity to plane gravitational waves (GWs) of DUAL detectors of cylindrical shape. Calculations make use of the finite element method to simulate the responses to the GW Riemann tensor of a single-mass DUAL (SMD) and of a tapered cylinder (TC) in their wide sensitivity bandwidth. We show that one SMD or a pair of TCs is able to cover both GW polarization amplitudes from almost all incoming directions. We discuss the achievable enhancement in tackling the inverse problem for high frequency [~(2–5) kHz] GWs by adding a TC detector to the future advanced LIGO–VIRGO network