Source author record

Julien Salomon

Julien Salomon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2022arXiv

Analysis of the Blade Element Momentum Theory

The blade element momentum (BEM) theory introduced by Lock et al. and formulated in its modern form by Glauert provides a framework to model the aerodynamic interaction between a turbine and a fluid flow. This theory is used either to estimate turbine efficiency or as a design aid. However, a lack of mathematical interpretation limits the understanding of some of its issues. The aim of this paper is to propose an analysis of BEM equations. Our approach is based on a reformulation of Glauert's model which enables us to identify criteria to guarantee the existence of solutions, analyze the convergence of usual and new (and more efficient) solution algorithms, and study turbine design procedures. The mathematical analysis is completed by numerical experiments.

preprint2022arXiv

Optimal periodic resource allocation in reactive dynamical systems: Application to microalgal production

In this article, we focus on a periodic resource allocation problem applied to a dynamical system which comes from a biological system. More precisely, we consider a system with $N$ resources and $N$ activities, each activity use the allocated resource to evolve up to a given time $T > 0$ where a control (represented by a given permutation) will be applied on the system to reallocate the resources. The goal is to find the optimal control strategies which optimize the cost or the benefit of the system. This problem can be illustrated by an industrial biological application, namely, the optimization of a mixing strategy to enhance the growth rate in a microalgal raceway system. A mixing device, such as a paddle wheel, is considered to control the rearrangement of the depth of the algae cultures, hence the light perceived at each lap. We prove that if the dynamics of the system is periodic, then the period corresponds to one reallocation whatever the order of the involved permutation matrix is. A nonlinear optimization problem for one reallocation process is then introduced. Since $N!$ permutations need to be tested in the general case, it can be numerically solved only for a limited number of $N$. To overcome this difficulty, we introduce a second optimization problem which provides a suboptimal solution of the initial problem, but whose solution can be determined explicitly. A sufficient condition to characterize cases where the two problems have the same solution is given. Some numerical experiments are performed to assess the benefit of optimal strategies in various settings.

preprint2020arXiv

Optimization of Bathymetry for Long Waves with Small Amplitude

This paper deals with bathymetry-oriented optimization in the case of long waves with small amplitude. Under these two assumptions, the free-surface incompressible Navier-Stokes system can be written as a wave equation where the bathymetry appears as a parameter in the spatial operator. Looking then for time-harmonic fields and writing the bottom topography as a perturbation of a flat bottom, we end up with a heterogeneous Helmholtz equation with impedance boundary condition. In this way, we study some PDE-constrained optimization problem for a Helmholtz equation in heterogeneous media whose coefficients are only bounded with bounded variation. We provide necessary condition for a general cost function to have at least one optimal solution. We also prove the convergence of a finite element approximation of the solution to the considered Helmholtz equation as well as the convergence of discrete optimum toward the continuous ones. We end this paper with some numerical experiments to illustrate the theoretical results and show that some of their assumptions could actually be removed.

preprint2020arXiv

PARAOPT: A parareal algorithm for optimality systems

The time parallel solution of optimality systems arising in PDE constraint optimization could be achieved by simply applying any time parallel algorithm, such as Parareal, to solve the forward and backward evolution problems arising in the optimization loop. We propose here a different strategy by devising directly a new time parallel algorithm, which we call ParaOpt, for the coupled forward and backward non-linear partial differential equations. ParaOpt is inspired by the Parareal algorithm for evolution equations, and thus is automatically a two-level method. We provide a detailed convergence analysis for the case of linear parabolic PDE constraints. We illustrate the performance of ParaOpt with numerical experiments both for linear and nonlinear optimality systems.

preprint2016arXiv

A fully efficient time-parallelized quantum optimal control algorithm

We present a time-parallelization method that enables to accelerate the computation of quantum optimal control algorithms. We show that this approach is approximately fully efficient when based on a gradient method as optimization solver: the computational time is approximately divided by the number of available processors. The control of spin systems, molecular orientation and Bose-Einstein condensates are used as illustrative examples to highlight the wide range of application of this numerical scheme.

preprint2015arXiv

A discrete-pulse optimal control algorithm with an application to spin systems

This article is aimed at extending the framework of optimal control techniques to the situation where the control field values are restricted to a finite set. We propose a generalization of the standard GRAPE algorithm suited to this constraint. We test the validity and the efficiency of this approach for the inversion of an inhomogeneous ensemble of spin systems with different offset frequencies. It is shown that a remarkable efficiency can be achieved even for a very limited number of discrete values. Some applications in Nuclear Magnetic Resonance are discussed.

preprint2014arXiv

Parareal in time 3D numerical solver for the LWR Benchmark neutron diffusion transient model

We present a parareal in time algorithm for the simulation of neutron diffusion transient model. The method is made efficient by means of a coarse solver defined with large time steps and steady control rods model. Using finite element for the space discretization, our implementation provides a good scalability of the algorithm. Numerical results show the efficiency of the parareal method on large light water reactor transient model corresponding to the Langenbuch-Maurer-Werner (LMW) benchmark [1].

preprint2014arXiv

Reduced basis methods for pricing options with the Black-Scholes and Heston model

In this paper, we present a reduced basis method for pricing European and American options based on the Black-Scholes and Heston model. To tackle each model numerically, we formulate the problem in terms of a time dependent variational equality or inequality. We apply a suitable reduced basis approach for both types of options. The characteristic ingredients used in the method are a combined POD-Greedy and Angle-Greedy procedure for the construction of the primal and dual reduced spaces. Analytically, we prove the reproduction property of the reduced scheme and derive a posteriori error estimators. Numerical examples are provided, illustrating the approximation quality and convergence of our approach for the different option pricing models. Also, we investigate the reliability and effectivity of the error estimators.

preprint2012arXiv

A Reduced Basis Method for the Simulation of American Options

We present a reduced basis method for the simulation of American option pricing. To tackle this model numerically, we formulate the problem in terms of a time dependent variational inequality. Characteristic ingredients are a POD-greedy and an angle-greedy procedure for the construction of the primal and dual reduced spaces. Numerical examples are provided, illustrating the approximation quality and convergence of our approach.

preprint2012arXiv

Control through operators for quantum chemistry

We consider the problem of operator identification in quantum control. The free Hamiltonian and the dipole moment are searched such that a given target state is reached at a given time. A local existence result is obtained. As a by-product, our works reveals necessary conditions on the laser field to make the identification feasible. In the last part of this work, some algorithms are proposed to compute effectively these operators.

preprint2012arXiv

Parareal in time intermediate targets methods for optimal control problem

In this paper, we present a method that enables solving in parallel the Euler-Lagrange system associated with the optimal control of a parabolic equation. Our approach is based on an iterative update of a sequence of intermediate targets that gives rise to independent sub-problems that can be solved in parallel. This method can be coupled with the parareal in time algorithm. Numerical experiments show the efficiency of our method.

preprint2011arXiv

An intermediate targets method for time parallelization in optimal control

In this paper, we present a method that enables to solve in parallel the Euler-Lagrange system associated with the optimal control of a parabolic equation. Our approach is based on an iterative update of a sequence of intermediate targets and gives rise independent sub-problems that can be solved in parallel. Numerical experiments show the efficiency of our method.

preprint2011arXiv

Minimum-weight perfect matching for non-intrinsic distances on the line

Consider a real line equipped with a (not necessarily intrinsic) distance. We deal with the minimum-weight perfect matching problem for a complete graph whose points are located on the line and whose edges have weights equal to distances along the line. This problem is closely related to one-dimensional Monge-Kantorovich trasnport optimization. The main result of the present note is a "bottom-up" recursion relation for weights of partial minimum-weight matchings.

preprint2010arXiv

A monotonic method for solving nonlinear optimal control problems

Initially introduced in the framework of quantum control, the so-called "monotonic algorithms" have demonstrated excellent numerical performance when dealing with bilinear optimal control problems. This paper presents a unified formulation that can be applied to more nonlinear settings compatible with the hypothesis detailed below. In this framework, we show that the well-posedness of the general algorithm is related to a nonlinear evolution equation. We prove the existence of the solution to this equation and give important properties of the optimal control functional. Finally we show how the algorithm works for selected models from the literature and compare it with the gradient algorithm.

preprint2010arXiv

A smoothing monotonic convergent optimal control algorithm for NMR pulse sequence design

The past decade has demonstrated increasing interests in using optimal control based methods within coherent quantum controllable systems. The versatility of such methods has been demonstrated with particular elegance within nuclear magnetic resonance (NMR) where natural separation between coherent and dissipative spin dynamics processes has enabled coherent quantum control over long periods of time to shape the experiment to almost ideal adoption to the spin system and external manipulations. This has led to new design principles as well as powerful new experimental methods within magnetic resonance imaging, liquid-state and solid-state NMR spectroscopy. For this development to continue and expand, it is crucially important to constantly improve the underlying numerical algorithms to provide numerical solutions which are optimally compatible with implementation on current instrumentation and at same time are numerically stable and offer fast monotonic convergence towards the target. Addressing such aims, we here present a smoothing monotonically convergent algorithm for pulse sequence design in magnetic resonance which with improved optimization stability lead to smooth pulse sequence easier to implement experimentally and potentially understand within the analytical framework of modern NMR spectroscopy.

preprint2010arXiv

Analysis of the Toolkit method for the time-dependant Schrödinger equation

The goal of this paper is to provide an analysis of the "toolkit" method used in the numerical approximation of the time-dependent Schrödinger equation. The "toolkit" method is based on precomputation of elementary propagators and was seen to be very efficient in the optimal control framework. Our analysis shows that this method provides better results than the second order Strang operator splitting. In addition, we present two improvements of the method in the limit of low and large intensity control fields.

preprint2010arXiv

Fast transport optimization for Monge costs on the circle

Consider the problem of optimally matching two measures on the circle, or equivalently two periodic measures on the real line, and suppose the cost of matching two points satisfies the Monge condition. We introduce a notion of locally optimal transport plan, motivated by the weak KAM (Aubry-Mather) theory, and show that all locally optimal transport plans are conjugate to shifts and that the cost of a locally optimal transport plan is a convex function of a shift parameter. This theory is applied to a transportation problem arising in image processing: for two sets of point masses on the circle, both of which have the same total mass, find an optimal transport plan with respect to a given cost function satisfying the Monge condition. In the circular case the sorting strategy fails to provide a unique candidate solution and a naive approach requires a quadratic number of operations. For the case of $N$ real-valued point masses we present an O(N |log epsilon|) algorithm that approximates the optimal cost within epsilon; when all masses are integer multiples of 1/M, the algorithm gives an exact solution in O(N log M) operations.

preprint2010arXiv

Local matching indicators for transport with concave costs

In this note, we introduce a class of indicators that enable to compute efficiently optimal transport plans associated to arbitrary distributions of $N$ demands and $N$ supplies in $\mathbf{R}$ in the case where the cost function is concave. The computational cost of these indicators is small and independent of $N$. A hierarchical use of them enables to obtain an efficient algorithm.