On the efficiency and accuracy of interpolation methods for spectral codes

In this paper a general theory for interpolation methods on a rectangular grid is introduced. By the use of this theory an efficient B-spline based interpolation method for spectral codes is presented. The theory links the order of the interpolation method with its spectral properties. In this way many properties like order of continuity, order of convergence and magnitude of errors can be explained. Furthermore, a fast implementation of the interpolation methods is given. We show that the B-spline based interpolation method has several advantages compared to other methods. First, the order of continuity of the interpolated field is higher than for other methods. Second, only one FFT is needed whereas e.g. Hermite interpolation needs multiple FFTs for computing the derivatives. Third, the interpolation error almost matches the one of Hermite interpolation, a property not reached by other methods investigated.Comment: 19 pages, 5 figure

Design of a freeform two-reflector system to collimate and shape a point source distribution

In this paper we propose a method to compute a freeform reflector system for collimating and shaping a beam from a point source. We construct these reflectors such that the radiant intensity of the source is converted into a desired target. An important generalization in our approach compared to previous research is that the output beam can be in an arbitrary direction. The design problem is approached by using a generalized Monge-Amp\`ere equation. This equation is solved using a least-squares algorithm for non-quadratic cost functions. This algorithm calculates the optical map, from which we can then compute the surfaces. We test our algorithm on two cases. First we consider a uniform source and target distribution. Next, we use the model of a laser diode light source and a ring-shaped target distribution

Design of two-dimensional reflective imaging systems: An approach based on inverse methods

Imaging systems are inherently prone to aberrations. We present an optimization method to design two-dimensional freeform reflectors that minimize aberrations for various parallel ray beams incident on the optical system. We iteratively design reflectors using inverse methods from non-imaging optics and optimize them to obtain a system that produces minimal aberrations. This is done by minimizing a merit function that quantifies aberrations and is dependent on the energy distributions at the source and target of an optical system, which are input parameters essential for inverse freeform design. The proposed method is tested for two configurations: a single-reflector system and a double-reflector system. Classical designs consisting of aspheric elements are well-known for their ability to minimize aberrations. We compare the performance of our freeform optical elements with classical designs. The optimized freeform designs outperform the classical designs in both configurations

A Compact High Order Finite Volume Scheme for Advection-Diffusion-Reaction Equations

We present a new integral representation for the flux of the advection-diffusion-reaction equation, which is based on the solution of a local boundary value problem for the entire equation, including the source term. The flux therefore consists of two parts, corresponding to the homogeneous and particular solution of the boundary value problem. Applying Gauss-Legendre quadrature rules to the integral representation gives the high order finite volume complete flux scheme, which is fourth order accurate for both diffusion dominated and advection dominated flow

A Study of Fluid Interfaces and Moving Contact Lines Using the Lattice Boltzmann Method

AbstractWe study the static and dynamical behavior of the contact line between two fluids and a solid plate by means of the Lattice Boltzmann method (LBM). The different fluid phases and their contact with the plate are simulated by means of standard Shan-Chen models. We investigate different regimes and compare the multicomponent vs. the multiphase LBM models near the contact line. A static interface profile is attained with the multiphase model just by balancing the hydrostatic pressure (due to gravity) with a pressure jump at the bottom. In order to study the same problem with the multicomponent case we propose and validate an idea of a body force acting only on one of the two fluid components. In order to reproduce results matching an infinite bath, boundary conditions at the bath side play a key role. We quantitatively compare open and wall boundary conditions and study their influence on the shape of the meniscus against static and lubrication theory solution

An Iterative Least-Squares Method for the Hyperbolic Monge-Amp\`ere Equation with Transport Boundary Condition

A least-squares method for solving the hyperbolic Monge-Amp\`ere equation with transport boundary condition is introduced. The method relies on an iterative procedure for the gradient of the solution, the so-called mapping. By formulating error functionals for the interior domain, the boundary, both separately and as linear combination, three minimization problems are solved iteratively to compute the mapping. After convergence, a fourth minimization problem, to compute the solution of the Monge-Amp\`ere equation, is solved. The approach is based on a least-squares method for the elliptic Monge-Amp\`ere equation by Prins et al., and is improved upon by the addition of analytical solutions for the minimization on the interior domain and by the introduction of two new boundary methods. Lastly, the iterative method is tested on a variety of examples. It is shown that, when the iterative method converges, second-order global convergence as function of the spatial discretization is obtained.Comment: 30 pages, 24 figure

Analytical modeling for the heat transfer in sheared flows of nanofluids

We developed a model for the enhancement of the heat flux by spherical and elongated nano- particles in sheared laminar flows of nano-fluids. Besides the heat flux carried by the nanoparticles the model accounts for the contribution of their rotation to the heat flux inside and outside the particles. The rotation of the nanoparticles has a twofold effect, it induces a fluid advection around the particle and it strongly influences the statistical distribution of particle orientations. These dynamical effects, which were not included in existing thermal models, are responsible for changing the thermal properties of flowing fluids as compared to quiescent fluids. The proposed model is strongly supported by extensive numerical simulations, demonstrating a potential increase of the heat flux far beyond the Maxwell-Garnet limit for the spherical nanoparticles. The road ahead which should lead towards robust predictive models of heat flux enhancement is discussed.Comment: 14 pages, 10 figures, submitted to PR