Turbulent viscosity optimized by data assimilation

International audienceAs an alternative approach to classical turbulence modelling using a first or second order closure, the data assimilation method of optimal control is applied to estimate a time and space-dependent turbulent viscosity in a three-dimensional oceanic circulation model. The optimal control method, described for a 3-D primitive equation model, involves the minimization of a cost function that quantifies the discrepancies between the simulations and the observations. An iterative algorithm is obtained via the adjoint model resolution. In a first experiment, a k + L model is used to simulate the one-dimensional development of inertial oscillations resulting from a wind stress at the sea surface and with the presence of a halocline. These results are used as synthetic observations to be assimilated. The turbulent viscosity is then recovered without the k + L closure, even with sparse and noisy observations. The problems of controllability and of the dimensions of the control are then discussed. A second experiment consists of a two-dimensional schematic simulation. A 2-D turbulent viscosity field is estimated from data on the initial and final states of a coastal upwelling event

Particle exchange and residence times in the North Western Mediterranean

The effects of the hydrodynamic processes on the distribution of passive drifters in the Gulf of Lions (GoL) have been investigated using a Lagrangian approach coupled with a 3D circulation model (Symphonie). We consider passive drifters, for which transport processes are determined solely by the 3D flow fields, which are in turn primarily forced by the North Mediterranean Current (NMC) and by the Rhˆone fresh-water inputs. The model reproduces 600 3D Lagrangian trajectories of particles released along the coastal area of the GoL during the winter period (January-February). The GoL has been divided into four sectors, each corresponding to a zone playing a strategic role in the hydrodynamics of the study area. The macroscopic characteristics of the transport on the shelf zone are analyzed in terms of total concentration and residence times of the cluster released in the basin. Particle distributions are strongly related to the mesoscale and sub-mesoscale hydrodynamic structures on the shelf and to the offshore circulation associated with the NMC. Two crucial areas are identified: a dispersive zone, corresponding to the central part of the continental shelf, and a wide offshore zone, representing an area of both aggregation and transition

Measurement of the permanent electric dipole moment of the neutron

We present the result of an experiment to measure the electric dipole moment EDM) of the neutron at the Paul Scherrer Institute using Ramsey's method of separated oscillating magnetic fields with ultracold neutrons (UCN). Our measurement stands in the long history of EDM experiments probing physics violating time reversal invariance. The salient features of this experiment were the use of a Hg-199 co-magnetometer and an array of optically pumped cesium vapor magnetometers to cancel and correct for magnetic field changes. The statistical analysis was performed on blinded datasets by two separate groups while the estimation of systematic effects profited from an unprecedented knowledge of the magnetic field. The measured value of the neutron EDM is d_{\rm n} = (0.0\pm1.1_{\rm stat}\pm0.2_{\rmsys})\times10^{-26}e\,{\rm cm}

The n2EDM experiment at the Paul Scherrer Institute

We present the new spectrometer for the neutron electric dipole moment (nEDM) search at the Paul Scherrer Institute (PSI), called n2EDM. The setup is at room temperature in vacuum using ultracold neutrons. n2EDM features a large UCN double storage chamber design with neutron transport adapted to the PSI UCN source. The design builds on experience gained from the previous apparatus operated at PSI until 2017. An order of magnitude increase in sensitivity is calculated for the new baseline setup based on scalable results from the previous apparatus, and the UCN source performance achieved in 2016

Magnetic-field uniformity in neutron electric-dipole-moment experiments

Magnetic field uniformity is of the utmost importance in experiments to measure the electric dipole moment of the neutron. A general parametrization of the magnetic field in terms of harmonic polynomial modes is proposed, going beyond the linear-gradients approximation. We review the main undesirable effects of non-uniformities: depolarization of ultracold neutrons, and Larmor frequency shifts of neutrons and mercury atoms. The theoretical predictions for these effects were verified by dedicated measurements with the single-chamber nEDM apparatus installed at the Paul Scherrer Institute

Optically pumped Cs magnetometers enabling a high-sensitivity search for the neutron electric dipole moment

An array of 16 laser-pumped scalar Cs magnetometers was part of the neutron electric dipole moment (nEDM) experiment taking data at the Paul Scherrer Institute in 2015 and 2016. It was deployed to measure the gradients of the experiment's magnetic field and to monitor their temporal evolution. The originality of the array lies in its compact design, in which a single near-infrared diode laser drives all magnetometers that are located in a high-vacuum chamber, with a selection of the sensors mounted on a high-voltage electrode. We describe details of the Cs sensors' construction and modes of operation, emphasizing the accuracy and sensitivity of the magnetic-field readout. We present two applications of the magnetometer array directly beneficial to the nEDM experiment: (i) the implementation of a strategy to correct for the drift of the vertical magnetic-field gradient and (ii) a procedure to homogenize the magnetic field. The first reduces the uncertainty of the nEDM result. The second enables transverse neutron spin relaxation times exceeding 1500 s, improving the statistical sensitivity of the nEDM experiment by about 35% and effectively increasing the rate of nEDM data taking by a factor of 1.8

Data blinding for the nEDM experiment at PSI

Psychological bias towards, or away from, prior measurements or theory predictions is an intrinsic threat to any data analysis. While various methods can be used to try to avoid such a bias, e.g. actively avoiding looking at the result, only data blinding is a traceable and trustworthy method that can circumvent the bias and convince a public audience that there is not even an accidental psychological bias. Data blinding is nowadays a standard practice in particle physics, but it is particularly difficult for experiments searching for the neutron electric dipole moment (nEDM), as several cross measurements, in particular of the magnetic field, create a self-consistent network into which it is hard to inject a false signal. We present an algorithm that modifies the data without influencing the experiment. Results of an automated analysis of the data are used to change the recorded spin state of a few neutrons within each measurement cycle. The flexible algorithm may be applied twice (or more) to the data, thus providing the option of sequentially applying various blinding offsets for separate analysis steps with independent teams. The subtle manner in which the data are modified allows one subsequently to adjust the algorithm and to produce a re-blinded data set without revealing the initial blinding offset. The method was designed for the 2015/2016 measurement campaign of the nEDM experiment at the Paul Scherrer Institute. However, it can be re-used with minor modification for the follow-up experiment n2EDM, and may be suitable for comparable projects elsewhere

The design of the n2EDM experiment: nEDM Collaboration

We present the design of a next-generation experiment, n2EDM, currently under construction at the ultracold neutron source at the Paul Scherrer Institute (PSI) with the aim of carrying out a high-precision search for an electric dipole moment of the neutron. The project builds on experience gained with the previous apparatus operated at PSI until 2017, and is expected to deliver an order of magnitude better sensitivity with provision for further substantial improvements. An overview is of the experimental method and setup is given, the sensitivity requirements for the apparatus are derived, and its technical design is described

Turbulent viscosity optimized by data assimilation

As an alternative approach to classical turbulence modelling using a first or second order closure, the data assimilation method of optimal control is applied to estimate a time and space-dependent turbulent viscosity in a three-dimensional oceanic circulation model. The optimal control method, described for a 3-D primitive equation model, involves the minimization of a cost function that quantifies the discrepancies between the simulations and the observations. An iterative algorithm is obtained via the adjoint model resolution. In a first experiment, a k + L model is used to simulate the one-dimensional development of inertial oscillations resulting from a wind stress at the sea surface and with the presence of a halocline. These results are used as synthetic observations to be assimilated. The turbulent viscosity is then recovered without the k + L closure, even with sparse and noisy observations. The problems of controllability and of the dimensions of the control are then discussed. A second experiment consists of a two-dimensional schematic simulation. A 2-D turbulent viscosity field is estimated from data on the initial and final states of a coastal upwelling event.Key words. Oceanography: general (numerical modelling) · Oceanography: physical (turbulence · diffusion · and mixing processes

Geostatistical analysis of three dimensional current patterns in coastal oceanography: application to the Gulf of Lions (nw mediterranean sea)

International audienc