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

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

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

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

Accuracy of magnetic resonance studies in the detection of chondral and labral lesions in femoroacetabular impingement : systematic review and meta-analysis

Background: Several types of Magnetic resonance imaging (MRI) are commonly used in imaging of femoroacetabular impingement (FAI), however till now there are no clear protocols and recommendations for each type. The aim of this meta-analysis is to detect the accuracy of conventional magnetic resonance imaging (cMRI), direct magnetic resonance arthrography (dMRA) and indirect magnetic resonance arthrography (iMRA) in the diagnosis of chondral and labral lesions in femoroacetabular impingement (FAI). Methods: A literature search was finalized on the 17th of May 2016 to collect all studies identifying the accuracy of cMRI, dMRA and iMRA in diagnosing chondral and labral lesions associated with FAI using surgical results (arthroscopic or open) as a reference test. Pooled sensitivity and specificity with 95% confidence intervals using a random-effects meta-analysis for MRI, dMRA and iMRA were calculated also area under receiver operating characteristic (ROC) curve (AUC) was retrieved whenever possible where AUC is equivocal to diagnostic accuracy. Results: The search yielded 192 publications which were reviewed according inclusion and exclusion criteria then 21 studies fulfilled the eligibility criteria for the qualitative analysis with a total number of 828 cases, lastly 12 studies were included in the quantitative meta-analysis. Meta-analysis showed that as regard labral lesions the pooled sensitivity, specificity and AUC for cMRI were 0.864, 0.833 and 0.88 and for dMRA were 0.91, 0.58 and 0.92. While in chondral lesions the pooled sensitivity, specificity and AUC for cMRI were 0.76, 0.72 and 0.75 and for dMRA were 0.75, 0.79 and 0.83, while for iMRA were sensitivity of 0.722 and specificity of 0.917. Conclusions: The present meta-analysis showed that the diagnostic test accuracy was superior for dMRA when compared with cMRI for detection of labral and chondral lesions. The diagnostic test accuracy was superior for labral lesions when compared with chondral lesions in both cMRI and dMRA. Promising results are obtained concerning iMRA but further studies still needed to fully assess its diagnostic accuracy

The finite volume-complete flux scheme for advection- diffusion-reaction equations The finite volume-complete flux scheme for advection-diffusion-reaction equations by The finite volume-complete flux scheme for advection-diffusion-reaction equations

Abstract We present a new finite volume scheme for the advection-diffusion-reaction equation. The scheme is second order accurate in the grid size, both for dominant diffusion and dominant advection, and has only a three-point coupling in each spatial direction. Our scheme is based on a new integral representation for the flux of the one-dimensional advection-diffusion-reaction equation, which is derived from 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 suitable quadrature rules to the integral representation gives the complete flux scheme. Extensions of the complete flux scheme to two-dimensional and time-dependent problems are derived, containing the cross flux term or the time derivative in the inhomogeneous flux, respectively. The resulting finite volume-complete flux scheme is validated for several test problems

The finite volume-complete flux scheme for one- dimensional advection-diffusion-reaction equations The finite volume-complete flux scheme for one-dimensional advection-diffusion-reaction equations

Abstract 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 suitable quadrature rules to the integral representation gives the complete flux scheme, which is second order accurate, uniformly in the local Peclet numbers. The flux approximation is combined with a finite volume method, and the resulting finite volume-complete flux scheme is validated for several test problems

Discretization and Parallel Iterative Schemes for Advection-Diffusion-Reaction Problems

Conservation laws of advection-diffusion-reaction (ADR) type are ubiquitous in continuum physics. In this paper we outline discretization of these problems and iterative schemes for the resulting linear system. For discretization we use the finite volume method in combination with the complete flux scheme. The numerical flux is the superposition of a homogeneous flux, corresponding to the advection-diffusion operator, and the inhomogeneous flux, taking into account the effect of the source term (ten Thije Boonkkamp and Anthonissen, J Sci Comput 46(1): 47-70, 2011). For a three-dimensional conservation law this results in a 27point coupling for the unknown as well as the source term. Direct solution of the sparse linear systems that arise in 3D ADR problems is not feasible due to fill-in. Iterative solution of such linear systems remains to be the only efficient alternative which requires less memory and shorter time to solution compared to direct solvers. Iterative solvers require a preconditioner to reduce the number of iterations. We present results using several different preconditioning techniques and study their effectiveness