-
Cross-Chip Partial Reconfiguration for the Initialisation of Modular and Scalable Heterogeneous Systems
Authors:
Marvin Fuchs,
Hendrik Krause,
Timo Muscheid,
Lukas Scheller,
Luis E. Ardila-Perez,
Oliver Sander
Abstract:
The almost unlimited possibilities to customize the logic in an FPGA are one of the main reasons for the versatility of these devices. Partial reconfiguration exploits this capability even further by allowing to replace logic in predefined FPGA regions at runtime. This is especially relevant in heterogeneous SoCs, combining FPGA fabric with conventional processors on a single die. Tight integratio…
▽ More
The almost unlimited possibilities to customize the logic in an FPGA are one of the main reasons for the versatility of these devices. Partial reconfiguration exploits this capability even further by allowing to replace logic in predefined FPGA regions at runtime. This is especially relevant in heterogeneous SoCs, combining FPGA fabric with conventional processors on a single die. Tight integration and supporting frameworks like the FPGA subsystem in Linux facilitate use, for example, to dynamically load custom hardware accelerators. Although this example is one of the most common use cases for partial reconfiguration, the possible applications go far beyond. We propose to use partial reconfiguration in combination with the AXI C2C cross-chip bus to extend the resources of heterogeneous MPSoC and RFSoC devices by connecting peripheral FPGAs. With AXI C2C it is easily possible to link the programmable logic of the individual devices, but partial reconfiguration on peripheral FPGAs utilising the same channel is not officially supported. By using an AXI ICAP controller in combination with custom Linux drivers, we show that it is possible to enable the PS of the heterogeneous SoC to perform partial reconfiguration on peripheral FPGAs, and thus to seamlessly access and manage the entire multi-device system. As a result, software and FPGA firmware updates can be applied to the entire system at runtime, and peripheral FPGAs can be added and removed during operation.
△ Less
Submitted 16 August, 2024;
originally announced August 2024.
-
Optimal Demodulation Domain for Microwave SQUID Multiplexers in Presence of Readout System Noise
Authors:
M. E. García Redondo,
N. A. Müller,
J. M. Salum,
L. P. Ferreyro,
J. D. Bonilla-Neira,
J. M. Geria,
J. J. Bonaparte,
T. Muscheid,
R. Gartmann,
A. Almela,
M. R. Hampel,
A. E. Fuster,
L. E. Ardila-Perez,
M. Wegner,
M. Platino,
O. Sander,
S. Kempf,
M. Weber
Abstract:
The Microwave SQUID Multiplexer (μMUX) is the device of choice for the readout of a large number of Low-Temperature Detectors in a wide variety of experiments within the fields of astronomy and particle physics. While it offers large multiplexing factors, the system noise performance is highly dependent on the cold and warm-readout electronic systems used to read it out, as well as the demodulatio…
▽ More
The Microwave SQUID Multiplexer (μMUX) is the device of choice for the readout of a large number of Low-Temperature Detectors in a wide variety of experiments within the fields of astronomy and particle physics. While it offers large multiplexing factors, the system noise performance is highly dependent on the cold and warm-readout electronic systems used to read it out, as well as the demodulation domain and parameters chosen. In order to understand the impact of the readout systems in the overall detection system noise performance, first we extended the available μMUX simulation frameworks including additive and multiplicative noise sources in the probing tones (i.e. phase and amplitude noise), along with the capability of demodulating the scientific data, either in resonator's phase or scattering amplitude. Then, considering the additive noise as a dominant noise source, the optimum readout parameters to achieve minimum system noise were found for both open-loop and flux-ramp demodulation schemes in the aforementioned domains. Later, we evaluated the system noise sensitivity to multiplicative noise sources under the optimum readout parameters. Finally, as a case study, we evaluated the optimal demodulation domain and expected system noise level for a typical Software-Defined Radio (SDR) readout system. This work leads to an improved system performance prediction and noise engineering based on the available readout electronics and selected demodulation domain.
△ Less
Submitted 19 August, 2024; v1 submitted 10 June, 2024;
originally announced June 2024.
-
ZynqMP-based board-management mezzanines for Serenity ATCA-blades
Authors:
T. Mehner,
L. E. Ardila-Perez,
M. N. Balzer,
O. Sander,
D. Tcherniakhovski,
M. Schleicher,
M. Fuchs,
G. Fedi,
G. Gimas,
G. M. Iles,
M. Pesaresi,
A. W. Rose,
T. Schuh
Abstract:
In the context of the CMS Phase-2 tracker back-end processing system, two mezzanines based on the Zynq Ultrascale+ Multi-Processor System-on-Chip (MPSoC) device have been developed to serve as centralized slow control and board management solution for the Serenity-family \textcolor{black}{Advanced Telecommunications Computing Architecture (ATCA)} blades.
This paper presents the developments of t…
▽ More
In the context of the CMS Phase-2 tracker back-end processing system, two mezzanines based on the Zynq Ultrascale+ Multi-Processor System-on-Chip (MPSoC) device have been developed to serve as centralized slow control and board management solution for the Serenity-family \textcolor{black}{Advanced Telecommunications Computing Architecture (ATCA)} blades.
This paper presents the developments of the MPSoC mezzanines to execute the Intelligent Platform Management Controller (IPMC) software in the real-time capable processors of the MPSoC. In coordination with the Shelf Manager, once full-power is enabled, a CentOS-based Linux distribution is executed in the application processors of the MPSoC, on which EMPButler and the Serenity Management Shell (SMASH) are running.
△ Less
Submitted 4 April, 2024;
originally announced April 2024.
-
Full-Scale Readout Electronics for the ECHo-100k Experiment
Authors:
T. Muscheid,
R. Gartmann,
N. Karcher,
F. Schuderer,
M. Neidig,
L. E. Ardila-Perez,
M. Balzer,
S. Kempf,
O. Sander
Abstract:
Recent advances in the development of cryogenic particle detectors such as magnetic microcalorimeters (MMCs) allow the fabrication of sensor arrays with an increasing number of pixels. Since these detectors must be operated at the lowest temperatures, the readout of large detector arrays is still quite challenging. This is especially true for the ECHo experiment, which presently aims to simultaneo…
▽ More
Recent advances in the development of cryogenic particle detectors such as magnetic microcalorimeters (MMCs) allow the fabrication of sensor arrays with an increasing number of pixels. Since these detectors must be operated at the lowest temperatures, the readout of large detector arrays is still quite challenging. This is especially true for the ECHo experiment, which presently aims to simultaneously run 6,000 two-pixel detectors to investigate the electron neutrino mass. For this reason, we developed a readout system based on a microwave SQUID multiplexer ($μ$MUX) that is operated by a custom software-defined radio (SDR) at room-temperature. The SDR readout electronics consist of three distinct hardware units: a data processing board with a Xilinx ZynqUS+ MPSoC; a converter board that features DACs, ADCs, and a coherent clock distribution network; and a radio frequency front-end board to translate the signals between the baseband and the microwave domains. Here, we describe the characteristics of the full-scale SDR system. First, the generated frequency comb for driving the $μ$MUX was evaluated. Subsequently, by operating the SDR in direct loopback, the crosstalk of the individual channels after frequency demultiplexing was investigated. Finally, the system was used with a 16-channel $μ$MUX to evaluate the linearity of the SDR, and the noise contributed to the overall readout setup.
△ Less
Submitted 3 April, 2024;
originally announced April 2024.
-
High-Order Integration on regular triangulated manifolds reaches Super-Algebraic Approximation Rates through Cubical Re-parameterizations
Authors:
Gentian Zavalani,
Oliver Sander,
Michael Hecht
Abstract:
We present a novel methodology for deriving high-order quadrature rules (HOSQ) designed for the integration of scalar functions over regular embedded manifolds. To construct the rules, we introduce square-squeezing--a homeomorphic multilinear hypercube-simplex transformation--reparametrizing an initial flat triangulation of the manifold to a hypercube mesh. By employing square-squeezing, we approx…
▽ More
We present a novel methodology for deriving high-order quadrature rules (HOSQ) designed for the integration of scalar functions over regular embedded manifolds. To construct the rules, we introduce square-squeezing--a homeomorphic multilinear hypercube-simplex transformation--reparametrizing an initial flat triangulation of the manifold to a hypercube mesh. By employing square-squeezing, we approximate the integrand and the volume element for each hypercube domain of the reparameterized mesh through interpolation in Chebyshev-Lobatto grids. This strategy circumvents the Runge phenomenon, replacing the initial integral with a closed-form expression that can be precisely computed by high-order quadratures. We prove novel bounds of the integration error in terms of the $r^\text{th}$-order total variation of the integrand and the surface parameterization, predicting high algebraic approximation rates that scale solely with the interpolation degree and not, as is common, with the average simplex size. For smooth integrals whose total variation is constantly bounded with increasing $r$, the estimates prove the integration error to decrease even exponentially, while mesh refinements are limited to achieve algebraic rates. The resulting approximation power is demonstrated in several numerical experiments, particularly showcasing $p$-refinements to overcome the limitations of $h$-refinements for highly varying smooth integrals.
△ Less
Submitted 23 November, 2023;
originally announced November 2023.
-
RFSoC Gen3-Based Software-Defined Radio Characterization for the Readout System of Low-Temperature Bolometers
Authors:
M. E. García Redondo,
T. Muscheid,
R. Gartmann,
J. M. Salum,
L. P. Ferreyro,
N. A. Müller,
J. D. Bonilla-Neira,
J. M. Geria,
J. J. Bonaparte,
A. Almela,
L. E. Ardila-Perez,
M. R. Hampel,
A. E. Fuster,
M. Platino,
O. Sander,
M. Weber,
A. Etchegoyen
Abstract:
This work reports the performance evaluation of an SDR readout system based on the latest generation (Gen3) of the AMD's Radio Frequency System-on-Chip (RFSoC) processing platform, which integrates a full-stack processing system and a powerful FPGA with up to 32 high-speed and high-resolution 14-bit Digital-to-Analog Converters (DACs) and Analog-to-Digital Converters (ADCs). The proposed readout s…
▽ More
This work reports the performance evaluation of an SDR readout system based on the latest generation (Gen3) of the AMD's Radio Frequency System-on-Chip (RFSoC) processing platform, which integrates a full-stack processing system and a powerful FPGA with up to 32 high-speed and high-resolution 14-bit Digital-to-Analog Converters (DACs) and Analog-to-Digital Converters (ADCs). The proposed readout system uses a previously developed multi-band, double-conversion IQ RF-mixing board targeting a multiplexing factor of approximately 1,000 bolometers in a bandwidth between 4 and 8 GHz, in line with state-of-the-art microwave SQUID multiplexers ($μ$MUX). The characterization of the system was performed in two stages, under the conditions typically imposed by the multiplexer and the cold readout circuit. First, in transmission, showing that noise and spurious levels of the generated tones are close to the values imposed by the cold readout. Second, in RF loopback, presenting noise values better than -100 dBc/Hz totally in agreement with the state-of-the-art readout systems. It was demonstrated that the RFSoC Gen3 device is a suitable enabling technology for the next generation of superconducting detector readout systems, reducing system complexity, increasing system integration, and achieving these goals without performance degradation.
△ Less
Submitted 8 May, 2024; v1 submitted 6 November, 2023;
originally announced November 2023.
-
Lessons learned while developing the Serenity-S1 ATCA card
Authors:
T. Mehner,
L. E. Ardila-Perez,
M. Balzer,
G. Fedi,
M. Fuchs,
A. Howard,
G. Iles,
M. Loutit,
S. Mansbridge,
F. Palla,
D. Parker,
M. Pesaresi,
A. Rose,
M. Saleh,
O. Sander,
M. Schleicher,
C. Strohman,
D. Tcherniakhovski,
T. Williams,
J. Zhao
Abstract:
The Serenity-S1 is a Xilinx Virtex Ultrascale+ based Advanced Telecommunications Computing Architecture (ATCA) processing blade that has been optimised for production. It incorporates many developments from the Serenity-A and Serenity-Z prototype cards and, where possible, adopts solutions being used across CERN. It also uses many new parts because commonly used parts have disappeared from the mar…
▽ More
The Serenity-S1 is a Xilinx Virtex Ultrascale+ based Advanced Telecommunications Computing Architecture (ATCA) processing blade that has been optimised for production. It incorporates many developments from the Serenity-A and Serenity-Z prototype cards and, where possible, adopts solutions being used across CERN. It also uses many new parts because commonly used parts have disappeared from the market during the semiconductor crisis, with only some returning.
Improvements to simplify manufacture, the performance of new components, some of the more difficult aspects of procurement, the performance of production-grade Samtec 25\,Gb/s optical firefly parts, and issues with the rack cooling infrastructure are discussed.
△ Less
Submitted 14 December, 2023; v1 submitted 3 November, 2023;
originally announced November 2023.
-
A geometrically nonlinear Cosserat shell model for orientable and non-orientable surfaces: Discretization with geometric finite elements
Authors:
Lisa Julia Nebel,
Oliver Sander,
Mircea Bîrsan,
Patrizio Neff
Abstract:
We investigate discretizations of a geometrically nonlinear elastic Cosserat shell with nonplanar reference configuration originally introduced by Bîrsan, Ghiba, Martin, and Neff in 2019. The shell model includes curvature terms up to order 5 in the shell thickness, which are crucial to reliably simulate high-curvature deformations such as near-folds or creases. The original model is generalized t…
▽ More
We investigate discretizations of a geometrically nonlinear elastic Cosserat shell with nonplanar reference configuration originally introduced by Bîrsan, Ghiba, Martin, and Neff in 2019. The shell model includes curvature terms up to order 5 in the shell thickness, which are crucial to reliably simulate high-curvature deformations such as near-folds or creases. The original model is generalized to shells that are not homeomorphic to a subset of $\mathbb{R}^2$. For this, we replace the originally planar parameter domain by an abstract two-dimensional manifold, and verify that the hyperelastic shell energy and three-dimensional reconstruction are invariant under changes of the local coordinate systems. This general approach allows to determine the elastic response for even non-orientable surfaces like the Möbius strip and the Klein bottle. We discretize the model with a geometric finite element method and, using that geometric finite elements are $H^1$-conforming, prove that the discrete shell model has a solution. Numerical tests then show the general performance and versatility of the model and discretization method.
△ Less
Submitted 1 August, 2023; v1 submitted 22 April, 2023;
originally announced April 2023.
-
An Implementation of a Channelizer based on a Goertzel Filter Bank for the Read-Out of Cryogenic Sensors
Authors:
L. P. Ferreyro,
M. García Redondo,
M. R. Hampel,
A. Almela,
A. Fuster,
J. Salum,
J. M. Geria,
J. Bonaparte,
J. Bonilla-Neira,
N. Müller,
N. Karcher,
O. Sander,
M. Platino,
A. Etchegoyen
Abstract:
In this work we present an application of the Goertzel Filter for the channelization of multi-tonal signals, typically used for the read-out of cryogenic sensors which are multiplexed in the frequency domain (FDM), by means of Microwave Superconducting Quantum Interference Device (SQUID) Multiplexer ($μ$MUX). We demonstrate how implementing a bank of many of these filters, can be used to perform a…
▽ More
In this work we present an application of the Goertzel Filter for the channelization of multi-tonal signals, typically used for the read-out of cryogenic sensors which are multiplexed in the frequency domain (FDM), by means of Microwave Superconducting Quantum Interference Device (SQUID) Multiplexer ($μ$MUX). We demonstrate how implementing a bank of many of these filters, can be used to perform a channelization of the multi-tonal input signal to retrieve the data added by the sensors. We show how this approach can be implemented in a resource-efficient manner in a Field Programmable Gate Array (FPGA) within the state-of-the-art, which allows great scalability for reading thousands of sensors; as is required by Radio Telescopes in Cosmic Microwave Background Radiation (CMB) surveys using cryogenic bolometers, particles detection like Neutrino mass estimation using cryogenic calorimeters or Quantum Computing.
△ Less
Submitted 2 March, 2023;
originally announced March 2023.
-
DTS-100G -- A versatile heterogeneous MPSoC board for cryogenic sensor readout
Authors:
Timo Muscheid,
Artur Boebel,
Nick Karcher,
Tomas Vanat,
Luis Ardila-Perez,
Igor Cheviakov,
Michael Schleicher,
Manfred Zimmer,
Matthias Balzer,
Oliver Sander
Abstract:
Heterogeneous devices such as the Multi-Processor System-on-Chip (MPSoC) from Xilinx are extremely valuable in custom instrumentation systems. This contribution presents the joint development of a heterogeneous MPSoC board called DTS-100G by DESY and KIT. The board is built around a Xilinx Zynq Ultrascale+ chip offering all available high-speed transceivers using QSFP28, 28 Gbps FireFly, FMC, and…
▽ More
Heterogeneous devices such as the Multi-Processor System-on-Chip (MPSoC) from Xilinx are extremely valuable in custom instrumentation systems. This contribution presents the joint development of a heterogeneous MPSoC board called DTS-100G by DESY and KIT. The board is built around a Xilinx Zynq Ultrascale+ chip offering all available high-speed transceivers using QSFP28, 28 Gbps FireFly, FMC, and FMC+ interfaces. The board is not designed for a particular application, but can be used as a generic DAQ platform for a variety of physics experiments. The DTS-100G board was successfully developed, built, and commissioned. ECHo-100k is the first experiment which will employ the board. This contribution shows the system architecture and explains how the DTS-100G board is a crucial component in the DAQ chain.
△ Less
Submitted 25 January, 2023;
originally announced January 2023.
-
Split Boot -- True Network-Based Booting on Heterogeneous MPSoCs
Authors:
Marvin Fuchs,
Luis E. Ardila-Perez,
Torben Mehner,
Oliver Sander
Abstract:
In the context of the High-Luminosity (HL) upgrade of the LHC, many custom ATCA electronics boards are being designed containing heterogeneous System-on-Chip (SoC) devices, more specifically the Xilinx Zynq UltraScale+ (ZUS+) family. While the application varies greatly, these devices are regularly used for performing board management tasks, making them a fundamental element in the correct operati…
▽ More
In the context of the High-Luminosity (HL) upgrade of the LHC, many custom ATCA electronics boards are being designed containing heterogeneous System-on-Chip (SoC) devices, more specifically the Xilinx Zynq UltraScale+ (ZUS+) family. While the application varies greatly, these devices are regularly used for performing board management tasks, making them a fundamental element in the correct operation of the board. The large number of hundreds of SoC devices creates significant challenges in their firmware deployment, maintenance, and accessibility. Even though U-Boot on ZUS+ devices supports network boot through the Preboot Execution Environment (PXE), the standard ZUS+ boot process contains application-specific information at earlier boot steps, particularly within the First Stage Bootloader (FSBL). This prevents the initialization of several devices from a universal image. Inspired by the PXE boot process on desktop PCs, this paper describes split boot, a novel boot method tailored to the specific needs of the ZUS+. All application-specific configuration is moved to a network storage device and automatically fetched during the boot process. We considered the entire process, from firmware and software development to binary distribution in a large-scale system. The developed method nicely integrates with the standard Xilinx development toolset workflow.
△ Less
Submitted 25 January, 2023; v1 submitted 13 January, 2023;
originally announced January 2023.
-
Simultaneous MMC readout using a tailored μMUX based readout system
Authors:
Daniel Richter,
Mathias Wegner,
Felix Ahrens,
Christian Enss,
Nick Karcher,
Oliver Sander,
Constantin Schuster,
Marc Weber,
Thomas Wolber,
Sebastian Kempf
Abstract:
Magnetic microcalorimeters (MMCs) are cryogenic, energy-dispersive single-particle detectors providing excellent energy resolution, intrinsically fast signal rise time, quantum efficiency close to 100\%, large dynamic range as well as almost ideal linear response. One of the remaining challenges to be overcome to ultimately allow for the utilization of large-scale MMC based detector arrays with th…
▽ More
Magnetic microcalorimeters (MMCs) are cryogenic, energy-dispersive single-particle detectors providing excellent energy resolution, intrinsically fast signal rise time, quantum efficiency close to 100\%, large dynamic range as well as almost ideal linear response. One of the remaining challenges to be overcome to ultimately allow for the utilization of large-scale MMC based detector arrays with thousands to millions of individual pixels is the realization of a SQUID based multiplexing technique particularly tailored for MMC readout. Within this context, we report on the first truly multiplexed readout of an MMC based detector array using a frequency-division multiplexing approach realized by a custom microwave SQUID multiplexer based readout system.
△ Less
Submitted 27 March, 2023; v1 submitted 14 November, 2022;
originally announced November 2022.
-
Benchmarking Numerical Algorithms for Harmonic Maps into the Sphere
Authors:
Sören Bartels,
Klaus Böhnlein,
Christian Palus,
Oliver Sander
Abstract:
We numerically benchmark methods for computing harmonic maps into the unit sphere, with particular focus on harmonic maps with singularities. For the discretization we compare two different approaches, both based on Lagrange finite elements. While the first method enforces the unit-length constraint only at the Lagrange nodes, the other one adds a pointwise projection to fulfill the constraint eve…
▽ More
We numerically benchmark methods for computing harmonic maps into the unit sphere, with particular focus on harmonic maps with singularities. For the discretization we compare two different approaches, both based on Lagrange finite elements. While the first method enforces the unit-length constraint only at the Lagrange nodes, the other one adds a pointwise projection to fulfill the constraint everywhere. For the solution of the resulting algebraic problems we compare a nonconforming gradient flow with a Riemannian trust-region method. Both are energy-decreasing and can be shown to converge globally to a stationary point of the Dirichlet energy. We observe that while the nonconforming and the conforming discretizations both show similar behavior, the second-order trust-region method needs less iterations than the solver based on gradient flow.
△ Less
Submitted 27 September, 2022;
originally announced September 2022.
-
Representative volume element approximations in elastoplastic spring networks
Authors:
Sabine Haberland,
Patrick Jaap,
Stefan Neukamm,
Oliver Sander,
Mario Varga
Abstract:
We study the large-scale behavior of a small-strain lattice model for a network composed of elastoplastic springs with random material properties. We formulate the model as an evolutionary rate independent system (ERIS). In an earlier work we derived a homogenized continuum model, which has the form of linearized elastoplasticity, as evolutionary $Γ$-limit as the lattice parameter tends to zero. I…
▽ More
We study the large-scale behavior of a small-strain lattice model for a network composed of elastoplastic springs with random material properties. We formulate the model as an evolutionary rate independent system (ERIS). In an earlier work we derived a homogenized continuum model, which has the form of linearized elastoplasticity, as evolutionary $Γ$-limit as the lattice parameter tends to zero. In the present paper we introduce a periodic representative volume element approximation (RVE) for the homogenized system. As a main result we prove convergence of the RVE approximation as the size of the RVE tends to infinity. We also show that the hysteretic stress-strain relation of the effective system can be described with help of a generalized Prandt-Ishlinskii operator, and we prove convergence of a periodic RVE approximation for that operator. We combine the RVE approximation with a numerical scheme for rate-independent systems and obtain a computational scheme that we use to numerically investigate the homogenized system in the specific case when the original network is given by a two-dimensional lattice model. We simulate the response of the system to cyclic and uniaxial, monotonic loading, and numerically investigate the convergence rate of the periodic RVE approximation. In particular, our simulations show that the RVE error decays with the same rate as the RVE error in the static case of linear elasticity.
△ Less
Submitted 26 September, 2022;
originally announced September 2022.
-
Two-level system hyperpolarization using a quantum Szilard engine
Authors:
Martin Spiecker,
Patrick Paluch,
Nicolas Gosling,
Niv Drucker,
Shlomi Matityahu,
Daria Gusenkova,
Simon Günzler,
Dennis Rieger,
Ivan Takmakov,
Francesco Valenti,
Patrick Winkel,
Richard Gebauer,
Oliver Sander,
Gianluigi Catelani,
Alexander Shnirman,
Alexey V. Ustinov,
Wolfgang Wernsdorfer,
Yonatan Cohen,
Ioan M. Pop
Abstract:
The innate complexity of solid state physics exposes superconducting quantum circuits to interactions with uncontrolled degrees of freedom degrading their coherence. By using a simple stabilization sequence we show that a superconducting fluxonium qubit is coupled to a two-level system (TLS) environment of unknown origin, with a relatively long energy relaxation time exceeding $50\,\text{ms}$. Imp…
▽ More
The innate complexity of solid state physics exposes superconducting quantum circuits to interactions with uncontrolled degrees of freedom degrading their coherence. By using a simple stabilization sequence we show that a superconducting fluxonium qubit is coupled to a two-level system (TLS) environment of unknown origin, with a relatively long energy relaxation time exceeding $50\,\text{ms}$. Implementing a quantum Szilard engine with an active feedback control loop allows us to decide whether the qubit heats or cools its TLS environment. The TLSs can be cooled down resulting in a four times lower qubit population, or they can be heated to manifest themselves as a negative temperature environment corresponding to a qubit population of $\sim 80\,\%$. We show that the TLSs and the qubit are each other's dominant loss mechanism and that the qubit relaxation is independent of the TLS populations. Understanding and mitigating TLS environments is therefore not only crucial to improve qubit lifetimes but also to avoid non-Markovian qubit dynamics.
△ Less
Submitted 22 May, 2024; v1 submitted 1 April, 2022;
originally announced April 2022.
-
A homogenized bending theory for prestrained plates
Authors:
Klaus Böhnlein,
Stefan Neukamm,
David Padilla-Garza,
Oliver Sander
Abstract:
The presence of prestrain can have a tremendous effect on the mechanical behavior of slender structures. Prestrained elastic plates show spontaneous bending in equilibrium -- a property that makes such objects relevant for the fabrication of active and functional materials. In this paper we study microheterogeneous, prestrained plates that feature nonflat equilibrium shapes. Our goal is to underst…
▽ More
The presence of prestrain can have a tremendous effect on the mechanical behavior of slender structures. Prestrained elastic plates show spontaneous bending in equilibrium -- a property that makes such objects relevant for the fabrication of active and functional materials. In this paper we study microheterogeneous, prestrained plates that feature nonflat equilibrium shapes. Our goal is to understand the relation between the properties of the prestrained microstructure and the global shape of the plate in mechanical equilibrium. To this end, we consider a three-dimensional, nonlinear elasticity model that describes a periodic material that occupies a domain with small thickness. We consider a spatially periodic prestrain described in the form of a multiplicative decomposition of the deformation gradient. By simultaneous homogenization and dimension reduction, we rigorously derive an effective plate model as a Γ-limit for vanishing thickness and period. That limit has the form of a nonlinear bending energy with an emergent spontaneous curvature term. The homogenized properties of the bending model (bending stiffness and spontaneous curvature) are characterized by corrector problems. For a model composite -- a prestrained laminate composed of isotropic materials -- we investigate the dependence of the homogenized properties on the parameters of the model composite. Secondly, we investigate the relation between the parameters of the model composite and the set of shapes with minimal bending energy. Our study reveals a rather complex dependence of these shapes on the composite parameters.
△ Less
Submitted 21 March, 2022;
originally announced March 2022.
-
Slowing down light in a qubit metamaterial
Authors:
Jan David Brehm,
Richard Gebauer,
Alexander Stehli,
Alexander N. Poddubny,
Oliver Sander,
Hannes Rotzinger,
Alexey V. Ustinov
Abstract:
The rapid progress in quantum information processing leads to a rising demand for devices to control the propagation of electromagnetic wave pulses and to ultimately realize a universal and efficient quantum memory. While in recent years significant progress has been made to realize slow light and quantum memories with atoms at optical frequencies, superconducting circuits in the microwave domain…
▽ More
The rapid progress in quantum information processing leads to a rising demand for devices to control the propagation of electromagnetic wave pulses and to ultimately realize a universal and efficient quantum memory. While in recent years significant progress has been made to realize slow light and quantum memories with atoms at optical frequencies, superconducting circuits in the microwave domain still lack such devices. Here, we demonstrate slowing down electromagnetic waves in a superconducting metamaterial composed of eight qubits coupled to a common waveguide, forming a waveguide quantum electrodynamics system. We analyze two complementary approaches, one relying on dressed states of the Autler-Townes splitting, and the other based on a tailored dispersion profile using the qubits tunability. Our time-resolved experiments show reduced group velocities of down to a factor of about 1500 smaller than in vacuum. Depending on the method used, the speed of light can be controlled with an additional microwave tone or an effective qubit detuning. Our findings demonstrate high flexibility of superconducting circuits to realize custom band structures and open the door to microwave dispersion engineering in the quantum regime.
△ Less
Submitted 14 February, 2022;
originally announced February 2022.
-
Numerical approaches for investigating quasiconvexity in the context of Morrey's conjecture
Authors:
Jendrik Voss,
Robert J. Martin,
Oliver Sander,
Siddhant Kumar,
Dennis M. Kochmann,
Patrizio Neff
Abstract:
Deciding whether a given function is quasiconvex is generally a difficult task. Here, we discuss a number of numerical approaches that can be used in the search for a counterexample to the quasiconvexity of a given function $W$. We will demonstrate these methods using the planar isotropic rank-one convex function \[
W_{\rm magic}^+(F)=\frac{λ_{\rm max}}{λ_{\rm min}}-\log\frac{λ_{\rm max}}{λ_{\rm…
▽ More
Deciding whether a given function is quasiconvex is generally a difficult task. Here, we discuss a number of numerical approaches that can be used in the search for a counterexample to the quasiconvexity of a given function $W$. We will demonstrate these methods using the planar isotropic rank-one convex function \[
W_{\rm magic}^+(F)=\frac{λ_{\rm max}}{λ_{\rm min}}-\log\frac{λ_{\rm max}}{λ_{\rm min}}+\log\det F=\frac{λ_{\rm max}}{λ_{\rm min}}+2\logλ_{\rm min}\,, \] where $λ_{\rm max}\geqλ_{\rm min}$ are the singular values of $F$, as our main example. In a previous contribution, we have shown that quasiconvexity of this function would imply quasiconvexity for all rank-one convex isotropic planar energies $W:\operatorname{GL}^+(2)\rightarrow\mathbb{R}$ with an additive volumetric-isochoric split of the form \[
W(F)=W_{\rm iso}(F)+W_{\rm vol}(\det F)=\widetilde W_{\rm iso}\bigg(\frac{F}{\sqrt{\det F}}\bigg)+W_{\rm vol}(\det F) \] with a concave volumetric part. This example is therefore of particular interest with regard to Morrey's open question whether or not rank-one convexity implies quasiconvexity in the planar case.
△ Less
Submitted 17 January, 2022;
originally announced January 2022.
-
Online Demodulation and Trigger for Flux-ramp Modulated SQUID Signals
Authors:
N. Karcher,
T. Muscheid,
T. Wolber,
D. Richter,
C. Enss,
S. Kempf,
O. Sander
Abstract:
Due to the periodic characteristics of SQUIDs, a suitable linearization technique is required for SQUID-based readout. Flux-ramp modulation is a common linearization technique and is typically applied for the readout of a microwave-SQUID-multiplexer as well as since recently also for dc-SQUIDs. Flux-ramp modulation requires another stage in the signal processing chain to demodulate the SQUID outpu…
▽ More
Due to the periodic characteristics of SQUIDs, a suitable linearization technique is required for SQUID-based readout. Flux-ramp modulation is a common linearization technique and is typically applied for the readout of a microwave-SQUID-multiplexer as well as since recently also for dc-SQUIDs. Flux-ramp modulation requires another stage in the signal processing chain to demodulate the SQUID output signal before further processing. For cryogenic microcalorimenters, these events are given by fast exponentially rising and slowly exponentially decaying pulses which shall be detected by a trigger engine and recorded by a storage logic. Since the data rate can be decreased significantly by demodulation and event detection, it is desirable to do both steps on the deployed fast FPGA logic during measurement before passing the data to a general-purpose processor.
In this contribution, we show the implementation of efficient multi-channel flux-ramp demodulation computed at run-time on a SoC-FPGA. Furthermore, a concept and implementation for an online trigger and buffer mechanism with its theoretical trigger loss rates depending on buffer size is presented. Both FPGA modules can be operated with up to 500 MHz clock frequency and can efficiently process 32 channels. Correct functionality and data reduction capability of the modules are demonstrated in measurements utilizing magnetic microcalorimeter irradiated with an Iron-55 source for event generation and read out by a microwave SQUID multiplexer.
△ Less
Submitted 22 October, 2021;
originally announced October 2021.
-
Conforming Finite Elements for $H(\text{sym}\,\text{Curl})$ and $H(\text{dev}\,\text{sym}\,\text{Curl})$
Authors:
Oliver Sander
Abstract:
We construct conforming finite elements for the spaces $H(\text{sym}\,\text{Curl})$ and $H(\text{dev}\,\text{sym}\,\text{Curl})$. Those are spaces of matrix-valued functions with symmetric or deviatoric-symmetric $\text{Curl}$ in a Lebesgue space, and they appear in various models of nonstandard solid mechanics. The finite elements are not $H(\text{Curl})$-conforming. We show the construction, pro…
▽ More
We construct conforming finite elements for the spaces $H(\text{sym}\,\text{Curl})$ and $H(\text{dev}\,\text{sym}\,\text{Curl})$. Those are spaces of matrix-valued functions with symmetric or deviatoric-symmetric $\text{Curl}$ in a Lebesgue space, and they appear in various models of nonstandard solid mechanics. The finite elements are not $H(\text{Curl})$-conforming. We show the construction, prove conformity and unisolvence, and point out optimal approximation error bounds.
△ Less
Submitted 18 June, 2021; v1 submitted 26 April, 2021;
originally announced April 2021.
-
The Parallel Full Approximation Scheme in Space and Time for a Parabolic Finite Element Problem
Authors:
Oliver Sander,
Ruth Schöbel,
Robert Speck
Abstract:
The parallel full approximation scheme in space and time (PFASST) is a parallel-in-time integrator that allows to integrate multiple time-steps simultaneously. It has been shown to extend scaling limits of spatial parallelization strategies when coupled with finite differences, spectral discretizations, or particle methods. In this paper we show how to use PFASST together with a finite element dis…
▽ More
The parallel full approximation scheme in space and time (PFASST) is a parallel-in-time integrator that allows to integrate multiple time-steps simultaneously. It has been shown to extend scaling limits of spatial parallelization strategies when coupled with finite differences, spectral discretizations, or particle methods. In this paper we show how to use PFASST together with a finite element discretization in space. While seemingly straightforward, the appearance of the mass matrix and the need to restrict iterates as well as residuals in space makes this task slightly more intricate. We derive the PFASST algorithm with mass matrices and appropriate prolongation and restriction operators and show numerically that PFASST can, after some initial iterations, gain two orders of accuracy per iteration.
△ Less
Submitted 17 March, 2021;
originally announced March 2021.
-
Flux ramp modulation based MHz frequency-division dc-SQUID multiplexer
Authors:
Daniel Richter,
Ludwig Hoibl,
Thomas Wolber,
Nick Karcher,
Andreas Fleischmann,
Christian Enss,
Marc Weber,
Oliver Sander,
Sebastian Kempf
Abstract:
We present a MHz frequency-division dc-SQUID multiplexer that is based on flux ramp modulation and a series array of $N$ identical current-sensing dc-SQUIDs with tightly coupled input coil. By running a periodic, sawtooth-shaped current signal through an additional modulation coil being tightly, but non-uniformly coupled to the individual SQUIDs, the voltage drop across the array changes according…
▽ More
We present a MHz frequency-division dc-SQUID multiplexer that is based on flux ramp modulation and a series array of $N$ identical current-sensing dc-SQUIDs with tightly coupled input coil. By running a periodic, sawtooth-shaped current signal through an additional modulation coil being tightly, but non-uniformly coupled to the individual SQUIDs, the voltage drop across the array changes according to the superposition of the flux-to-voltage characteristics of the individual SQUIDs within each cycle of the modulation signal. In this mode of operation, an input signal injected in the input coil of one of the SQUIDs and being quasi-static within a time frame adds a constant flux offset and leads to a phase shift of the associated SQUID characteristics. The latter is inherently proportional to the input signal and can be inferred by channelizing and down-converting the sampled array output voltage. Using a prototype multiplexer as well as a self-developed high-speed readout electronics for real-time phase determination, we demonstrate the simultaneous readout of four signal sources with MHz bandwidth per channel.
△ Less
Submitted 4 March, 2021; v1 submitted 16 January, 2021;
originally announced January 2021.
-
OpenIPMC: a free and open source Intelligent Platform Management Controller
Authors:
Luigi Calligaris,
André Cascadan,
Luis E. Ardila-Perez,
Bruno Casu,
Alison França da Costa,
Ailton Akira Shinoda,
Lucas Arruda Ramalho,
Oliver Sander
Abstract:
OpenIPMC is a free and open source firmware designed to operate as an Intelligent Platform Management Controller (IPMC). An IPMC is a fundamental component of electronic boards conformant to the Advanced Telecommunications Computing Architecture (ATCA) standard, currently being adopted by a number of high energy physics experiments, and is responsible for monitoring the health parameters of the bo…
▽ More
OpenIPMC is a free and open source firmware designed to operate as an Intelligent Platform Management Controller (IPMC). An IPMC is a fundamental component of electronic boards conformant to the Advanced Telecommunications Computing Architecture (ATCA) standard, currently being adopted by a number of high energy physics experiments, and is responsible for monitoring the health parameters of the board, managing its power states, and providing board control, debug and recovery functions to remote clients. OpenIPMC is based on the FreeRTOS real-time operating system and is designed to be architecture-independent, allowing it to be built for a variety of different Microcontrollers. Having a fully free and open source code is an innovative aspect for this kind of firmware, allowing full customization by the user. In this work we present the features and structure of OpenIPMC and its example implementations on Xilinx Zynq UltraScale+ (ZynqUS+), Espressif ESP32, and ST Microelectronics STM32 architectures.
△ Less
Submitted 1 June, 2021; v1 submitted 2 November, 2020;
originally announced November 2020.
-
Versatile Configuration and Control Framework for Real Time Data Acquisition Systems
Authors:
Nick Karcher,
Richard Gebauer,
Robin Bauknecht,
Rainer Illichmann,
Oliver Sander
Abstract:
Modern physics experiments often utilize FPGA-based systems for real-time data acquisition. Integrated analog electronics demand for complex calibration routines. Furthermore, versatile configuration and control of the whole system is a key requirement. Beside low-level register interface to the FPGA, also access to I$^2$C and SPI buses is often needed to configure the complete system. Calibration…
▽ More
Modern physics experiments often utilize FPGA-based systems for real-time data acquisition. Integrated analog electronics demand for complex calibration routines. Furthermore, versatile configuration and control of the whole system is a key requirement. Beside low-level register interface to the FPGA, also access to I$^2$C and SPI buses is often needed to configure the complete system. Calibration through an FPGA is inflexible and yields a complex hardware implementation. On the contrary, calibration through a remote system is possible but considerably slower due to repetitive network accesses. By using SoC-FPGA solutions with a microprocessor, more sophisticated configuration and calibration solutions, as well as standard remote access protocols, can be efficiently integrated in software.
Based on Xilinx Zynq US+ SoC-FPGAs, we implemented a versatile control framework. This software framework offers a convenient access to the hardware and a flexible abstraction via remote-procedure calls (RPCs). Based on the open source RPC library gRPC, functionality with low-latent control flow, complex algorithms, data conversions and processing, as well as configuration via external buses can be provided to a client via Ethernet. Furthermore, client interfaces for various programming languages can be generated automatically which eases collaboration among different working groups and integration into existing software. This contribution presents the framework as well as benchmarks regarding latency and data throughput.
△ Less
Submitted 30 October, 2020;
originally announced November 2020.
-
Quantum non-demolition dispersive readout of a superconducting artificial atom using large photon numbers
Authors:
Daria Gusenkova,
Martin Spiecker,
Richard Gebauer,
Madita Willsch,
Francesco Valenti,
Nick Karcher,
Lukas Grünhaupt,
Ivan Takmakov,
Patrick Winkel,
Dennis Rieger,
Alexey V. Ustinov,
Nicolas Roch,
Wolfgang Wernsdorfer,
Kristel Michielsen,
Oliver Sander,
Ioan M. Pop
Abstract:
Reading out the state of superconducting artificial atoms typically relies on dispersive coupling to a readout resonator. For a given system noise temperature, increasing the circulating photon number $\bar{n}$ in the resonator enables a shorter measurement time and is therefore expected to reduce readout errors caused by spontaneous atom transitions. However, increasing $\bar{n}$ is generally obs…
▽ More
Reading out the state of superconducting artificial atoms typically relies on dispersive coupling to a readout resonator. For a given system noise temperature, increasing the circulating photon number $\bar{n}$ in the resonator enables a shorter measurement time and is therefore expected to reduce readout errors caused by spontaneous atom transitions. However, increasing $\bar{n}$ is generally observed to also increase these transition rates. Here we present a fluxonium artificial atom in which we measure an overall flat dependence of the transition rates between its first two states as a function of $\bar{n}$, up to $\bar{n}\approx200$. Despite the fact that we observe the expected decrease of the dispersive shift with increasing readout power, the signal-to-noise ratio continuously improves with increasing $\bar{n}$. Even without the use of a parametric amplifier, at $\bar{n}=74$, we measure fidelities of 99% and 93% for feedback-assisted ground and excited state preparation, respectively.
△ Less
Submitted 20 November, 2020; v1 submitted 30 September, 2020;
originally announced September 2020.
-
Truncated Nonsmooth Newton Multigrid for phase-field brittle-fracture problems, with analysis
Authors:
Carsten Gräser,
Daniel Kienle,
Oliver Sander
Abstract:
We propose the Truncated Nonsmooth Newton Multigrid Method (TNNMG) as a solver for the spatial problems of the small-strain brittle-fracture phase-field equations. TNNMG is a nonsmooth multigrid method that can solve biconvex, block-separably nonsmooth minimization problems with linear time complexity. It exploits the variational structure inherent in the problem, and handles the pointwise irrever…
▽ More
We propose the Truncated Nonsmooth Newton Multigrid Method (TNNMG) as a solver for the spatial problems of the small-strain brittle-fracture phase-field equations. TNNMG is a nonsmooth multigrid method that can solve biconvex, block-separably nonsmooth minimization problems with linear time complexity. It exploits the variational structure inherent in the problem, and handles the pointwise irreversibility constraint on the damage variable directly, without regularization or the introduction of a local history field. In the paper we introduce the method and show how it can be applied to several established models of phase-field brittle fracture. We then prove convergence of the solver to a solution of the nonsmooth Euler--Lagrange equations of the spatial problem for any load and initial iterate. On the way, we show several crucial convexity and regularity properties of the models considered here. Numerical comparisons to an operator-splitting algorithm show a considerable speed increase, without loss of robustness.
△ Less
Submitted 10 March, 2023; v1 submitted 23 July, 2020;
originally announced July 2020.
-
The effective mitigation of greenhouse gas emissions from rice paddies without compromising yield by early-season drainage
Authors:
S. F. U. Islam,
J. W. van Groenigen,
L. S. Jensen,
B. O. Sander,
A. de Neergaard
Abstract:
Global rice production systems face two opposing challenges: the need to increase production to accommodate the world's growing population while simultaneously reducing greenhouse gas (GHG) emissions. Adaptations to drainage regimes are one of the most promising options for methane mitigation in rice production. Whereas several studies have focused on mid-season drainage (MD) to mitigate GHG emiss…
▽ More
Global rice production systems face two opposing challenges: the need to increase production to accommodate the world's growing population while simultaneously reducing greenhouse gas (GHG) emissions. Adaptations to drainage regimes are one of the most promising options for methane mitigation in rice production. Whereas several studies have focused on mid-season drainage (MD) to mitigate GHG emissions, early-season drainage (ED) varying in timing and duration has not been extensively studied. However, such ED periods could potentially be very effective since initial available Carbon levels (and thereby the potential for methanogenesis) can be very high in paddy systems with rice straw incorporation. This study tested the effectiveness of seven drainage regimes varying in their timing and duration (combinations of ED and MD) to mitigate methane (CH4) and nitrous oxide (N2O) emissions in a 101-day growth chamber experiment. Emissions were considerably reduced by early-season drainage compared to both conventional continuous flooding (CF) and the MD drainage regime. The results suggest that ED + MD drainage may have the potential to reduce methane emissions and yield-scaled global warming potential by 85 to 90% compared to CF and by 75 to 77% compared to MD only. A combination of (short or long) ED drainage and one MD drainage episode was found to be the most effective in mitigating methane emissions without negatively affecting yield. In particular, compared with CF, the long early-season drainage treatments LE + SM and LE + LM significantly decreased yield scaled global warming potential by 85% and 87% respectively. This was associated with carbon being stabilised early in the season, thereby reducing available carbon for methanogenesis. Overall nitrous oxide emissions were small and not significantly affected by ED.
△ Less
Submitted 4 June, 2020;
originally announced June 2020.
-
Reducing greenhouse gas emissions and grain arsenic and lead levels without compromising yield in organically produced rice
Authors:
S. F. U. Islam,
A. de Neergaard,
B. O. Sander,
L. S. Jensen,
R. Wassmann,
J. W. van Groenigen
Abstract:
Flooded rice production is crucial to global food security, but there are associated environmental concerns. In particular, it is a significant source of methane and nitrous oxide emissions and a large consumer of water resources, while arsenic, cadmium and lead levels in the grain are a serious health concern. There is also a tendency to use more organic fertilisers to close nutrient cycles, posi…
▽ More
Flooded rice production is crucial to global food security, but there are associated environmental concerns. In particular, it is a significant source of methane and nitrous oxide emissions and a large consumer of water resources, while arsenic, cadmium and lead levels in the grain are a serious health concern. There is also a tendency to use more organic fertilisers to close nutrient cycles, posing a threat of even higher greenhouse gas emissions and grain arsenic levels. It has been shown that alternate wetting and drying (AWD) water management reduces both water use and greenhouse gas emissions, but success at maintaining yields varies. This study tested the effect of early AWD (e-AWD) versus continuous flooding (CF) water management practices on grain yields, greenhouse gas emissions and grain arsenic, cadmium and lead levels in a split plot field experiment with organic fertilisers under organic management. The e-AWD water regime showed no difference in yield for the organic treatments. Yields significantly increased by 5 to 16 percent in the combination treatments. Root biomass and root length increased in the e-AWD treatments up to 72 and 41 percent, respectively. The e-AWD water regime reduced seasonal methane emissions by 71 to 85 percent for organic treatments and by 51 to 76 percent for combination treatments; this was linked to a 15 to 47 percent reduction in dissolved organic carbon. Nitrous oxide emissions increased by 23 to 305 percent but accounted for less than 20 percent of global warming potential. Area and yield scaled global warming potentials were reduced by 67 to 83 percent. The e-AWD regime altered soil redox potentials, resulting in a reduction in grain arsenic and lead concentrations of up to 66 % and 73 % respectively. Grain cadmium levels were also reduced up to 33 % in organic treatments.
△ Less
Submitted 4 June, 2020;
originally announced June 2020.
-
Reducing the impact of radioactivity on quantum circuits in a deep-underground facility
Authors:
Laura Cardani,
Francesco Valenti,
Nicola Casali,
Gianluigi Catelani,
Thibault Charpentier,
Massimiliano Clemenza,
Ivan Colantoni,
Angelo Cruciani,
Luca Gironi,
Lukas Grünhaupt,
Daria Gusenkova,
Fabio Henriques,
Marc Lagoin,
Maria Martinez,
Giorgio Pettinari,
Claudia Rusconi,
Oliver Sander,
Alexey V. Ustinov,
Marc Weber,
Wolfgang Wernsdorfer,
Marco Vignati,
Stefano Pirro,
Ioan M. Pop
Abstract:
As quantum coherence times of superconducting circuits have increased from nanoseconds to hundreds of microseconds, they are currently one of the leading platforms for quantum information processing. However, coherence needs to further improve by orders of magnitude to reduce the prohibitive hardware overhead of current error correction schemes. Reaching this goal hinges on reducing the density of…
▽ More
As quantum coherence times of superconducting circuits have increased from nanoseconds to hundreds of microseconds, they are currently one of the leading platforms for quantum information processing. However, coherence needs to further improve by orders of magnitude to reduce the prohibitive hardware overhead of current error correction schemes. Reaching this goal hinges on reducing the density of broken Cooper pairs, so-called quasiparticles. Here, we show that environmental radioactivity is a significant source of nonequilibrium quasiparticles. Moreover, ionizing radiation introduces time-correlated quasiparticle bursts in resonators on the same chip, further complicating quantum error correction. Operating in a deep-underground lead-shielded cryostat decreases the quasiparticle burst rate by a factor fifty and reduces dissipation up to a factor four, showcasing the importance of radiation abatement in future solid-state quantum hardware.
△ Less
Submitted 5 May, 2020;
originally announced May 2020.
-
Accelerating complex control schemes on a heterogeneous MPSoC platform for quantum computing
Authors:
Richard Gebauer,
Nick Karcher,
Jonas Hurst,
Marc Weber,
Oliver Sander
Abstract:
Control and readout of superconducting quantum bits (qubits) require microwave pulses with gigahertz frequencies and nanosecond precision. To generate and analyze these microwave pulses, we developed a versatile FPGA-based electronics platform. While basic functionality is directly handled within the FPGA, guaranteeing highest accuracy on the nanosecond timescale, more complex control schemes rend…
▽ More
Control and readout of superconducting quantum bits (qubits) require microwave pulses with gigahertz frequencies and nanosecond precision. To generate and analyze these microwave pulses, we developed a versatile FPGA-based electronics platform. While basic functionality is directly handled within the FPGA, guaranteeing highest accuracy on the nanosecond timescale, more complex control schemes render impractical to implement in hardware.
To provide deterministic timing and low latency with high flexibility, we developed the Taskrunner framework. It enables the execution of complex control schemes, so-called user tasks, on the real-time processing unit (RPU) of a heterogeneous Multiprocessor System-on-Chip (MPSoC). These user tasks are specified conveniently using standard C language and are compiled automatically by the MPSoC platform when loaded onto the RPU. We present the architecture of the Taskrunner framework as well as timing benchmarks and discuss applications in the field of quantum computing.
△ Less
Submitted 16 April, 2020;
originally announced April 2020.
-
State preparation of a fluxonium qubit with feedback from a custom FPGA-based platform
Authors:
Richard Gebauer,
Nick Karcher,
Daria Gusenkova,
Martin Spiecker,
Lukas Grünhaupt,
Ivan Takmakov,
Patrick Winkel,
Luca Planat,
Nicolas Roch,
Wolfgang Wernsdorfer,
Alexey V. Ustinov,
Marc Weber,
Martin Weides,
Ioan M. Pop,
Oliver Sander
Abstract:
We developed a versatile integrated control and readout instrument for experiments with superconducting quantum bits (qubits), based on a field-programmable gate array (FPGA) platform. Using this platform, we perform measurement-based, closed-loop feedback operations with $428 \, \mathrm{ns}$ platform latency. The feedback capability is instrumental in realizing active reset initialization of the…
▽ More
We developed a versatile integrated control and readout instrument for experiments with superconducting quantum bits (qubits), based on a field-programmable gate array (FPGA) platform. Using this platform, we perform measurement-based, closed-loop feedback operations with $428 \, \mathrm{ns}$ platform latency. The feedback capability is instrumental in realizing active reset initialization of the qubit into the ground state in a time much shorter than its energy relaxation time $T_1$. We show experimental results demonstrating reset of a fluxonium qubit with $99.4\,\%$ fidelity, using a readout-and-drive pulse sequence approximately $1.5 \, \mathrm{μs}$ long. Compared to passive ground state initialization through thermalization, with the time constant given by $T_1 = ~ 80 \, \mathrm{μs}$, the use of the FPGA-based platform allows us to improve both the fidelity and the time of the qubit initialization by an order of magnitude.
△ Less
Submitted 14 December, 2019;
originally announced December 2019.
-
The DUNE Framework: Basic Concepts and Recent Developments
Authors:
Peter Bastian,
Markus Blatt,
Andreas Dedner,
Nils-Arne Dreier,
Christian Engwer,
René Fritze,
Carsten Gräser,
Christoph Grüninger,
Dominic Kempf,
Robert Klöfkorn,
Mario Ohlberger,
Oliver Sander
Abstract:
This paper presents the basic concepts and the module structure of the Distributed and Unified Numerics Environment and reflects on recent developments and general changes that happened since the release of the first Dune version in 2007 and the main papers describing that state [1, 2]. This discussion is accompanied with a description of various advanced features, such as coupling of domains and…
▽ More
This paper presents the basic concepts and the module structure of the Distributed and Unified Numerics Environment and reflects on recent developments and general changes that happened since the release of the first Dune version in 2007 and the main papers describing that state [1, 2]. This discussion is accompanied with a description of various advanced features, such as coupling of domains and cut cells, grid modifications such as adaptation and moving domains, high order discretizations and node level performance, non-smooth multigrid methods, and multiscale methods. A brief discussion on current and future development directions of the framework concludes the paper.
△ Less
Submitted 22 June, 2020; v1 submitted 30 September, 2019;
originally announced September 2019.
-
Phonon traps reduce the quasiparticle density in superconducting circuits
Authors:
Fabio Henriques,
Francesco Valenti,
Thibault Charpentier,
Marc Lagoin,
Clement Gouriou,
Maria Martínez,
Laura Cardani,
Marco Vignati,
Lukas Grünhaupt,
Daria Gusenkova,
Julian Ferrero,
Sebastian T. Skacel,
Wolfgang Wernsdorfer,
Alexey V. Ustinov,
Gianluigi Catelani,
Oliver Sander,
Ioan M. Pop
Abstract:
Out of equilibrium quasiparticles (QPs) are one of the main sources of decoherence in superconducting quantum circuits, and are particularly detrimental in devices with high kinetic inductance, such as high impedance resonators, qubits, and detectors. Despite significant progress in the understanding of QP dynamics, pinpointing their origin and decreasing their density remain outstanding tasks. Th…
▽ More
Out of equilibrium quasiparticles (QPs) are one of the main sources of decoherence in superconducting quantum circuits, and are particularly detrimental in devices with high kinetic inductance, such as high impedance resonators, qubits, and detectors. Despite significant progress in the understanding of QP dynamics, pinpointing their origin and decreasing their density remain outstanding tasks. The cyclic process of recombination and generation of QPs implies the exchange of phonons between the superconducting thin film and the underlying substrate. Reducing the number of substrate phonons with frequencies exceeding the spectral gap of the superconductor should result in a reduction of QPs. Indeed, we demonstrate that surrounding high impedance resonators made of granular aluminum (grAl) with lower gapped thin film aluminum islands increases the internal quality factors of the resonators in the single photon regime, suppresses the noise, and reduces the rate of observed QP bursts. The aluminum islands are positioned far enough from the resonators to be electromagnetically decoupled, thus not changing the resonator frequency, nor the loading. We therefore attribute the improvements observed in grAl resonators to phonon trapping at frequencies close to the spectral gap of aluminum, well below the grAl gap.
△ Less
Submitted 20 August, 2019; v1 submitted 12 August, 2019;
originally announced August 2019.
-
The quasiconvex envelope of conformally invariant planar energy functions in isotropic hyperelasticity
Authors:
Robert J. Martin,
Jendrik Voss,
Ionel-Dumitrel Ghiba,
Oliver Sander,
Patrizio Neff
Abstract:
We consider conformally invariant energies $W$ on the group $\operatorname{GL}^+(2)$ of $2\times2$-matrices with positive determinant, i.e. $W\colon\operatorname{GL}^+(2)\to\mathbb{R}$ such that \[W(AFB) = W(F) \qquad\text{for all }\; A,B\in\{aR\in\operatorname{GL}^+(2) \,|\, a\in(0,\infty)\,,\; R\in\operatorname{SO}(2)\}\,,\] where $\operatorname{SO}(2)$ denotes the special orthogonal group, and…
▽ More
We consider conformally invariant energies $W$ on the group $\operatorname{GL}^+(2)$ of $2\times2$-matrices with positive determinant, i.e. $W\colon\operatorname{GL}^+(2)\to\mathbb{R}$ such that \[W(AFB) = W(F) \qquad\text{for all }\; A,B\in\{aR\in\operatorname{GL}^+(2) \,|\, a\in(0,\infty)\,,\; R\in\operatorname{SO}(2)\}\,,\] where $\operatorname{SO}(2)$ denotes the special orthogonal group, and provide an explicit formula for the (notoriously difficult to compute) quasiconvex envelope of these functions. Our results, which are based on the representation $W(F)=h(\frac{λ_1}{λ_2})$ of $W$ in terms of the singular values $λ_1,λ_2$ of $F$, are applied to a number of example energies in order to demonstrate the convenience of the eigenvalue-based expression compared to the more common representation in terms of the distortion $\mathbb{K}:=\frac12\frac{\lVert F\rVert^2}{\det F}$. Special cases of our results can be obtained from earlier works by Astala et al. and Yan.
△ Less
Submitted 31 December, 2018;
originally announced January 2019.
-
Software-defined Radio Readout System for the ECHo Experiment
Authors:
Oliver Sander,
Nick Karcher,
Oliver Krömer,
Marc Weber,
Sebastian Kempf,
Mathias Wegner,
Christian Enss
Abstract:
Metallic magnetic calorimeters (MMCs) are cryogenic detectors that offer an excellent energy resolution, a signal rise time of <100 ns, a high dynamic range and almost optimal linearity. MMCs are of high interest for many experiments. One of them, the ECHo experiment, requires the utilization of large MMC detector arrays. The readout of such MMC arrays is a challenging task, which can be tackled u…
▽ More
Metallic magnetic calorimeters (MMCs) are cryogenic detectors that offer an excellent energy resolution, a signal rise time of <100 ns, a high dynamic range and almost optimal linearity. MMCs are of high interest for many experiments. One of them, the ECHo experiment, requires the utilization of large MMC detector arrays. The readout of such MMC arrays is a challenging task, which can be tackled using software-defined radios (SDRs). Though SDR is a well-known approach in communications engineering, a dedicated implementation for frequency division multiplexed readout of MMCs is new and one of the technological key elements of the ECHo project. ECHo will be the first experiment to use microwave SQUID multiplexed MMC detectors and therefore pioneering the hardware, firmware and software development in this domain. This paper presents the detailed concepts and current status of the development.
△ Less
Submitted 24 June, 2018;
originally announced June 2018.
-
Function space bases in the dune-functions module
Authors:
Christian Engwer,
Carsten Gräser,
Steffen Müthing,
Oliver Sander
Abstract:
The dune-functions Dune module provides interfaces for functions and function space bases. It forms one abstraction level above grids, shape functions, and linear algebra, and provides infrastructure for full discretization frameworks like dune-pdelab and dune-fem. This document describes the function space bases provided by dune-functions. These are based on an abstract description of bases for p…
▽ More
The dune-functions Dune module provides interfaces for functions and function space bases. It forms one abstraction level above grids, shape functions, and linear algebra, and provides infrastructure for full discretization frameworks like dune-pdelab and dune-fem. This document describes the function space bases provided by dune-functions. These are based on an abstract description of bases for product spaces as trees of simpler bases. From this description, many different numberings of degrees of freedom by multi-indices can be derived in a natural way. We describe the abstract concepts, document the programmer interface, and give a complete example program that solves the stationary Stokes equation using Taylor-Hood elements.
△ Less
Submitted 25 June, 2018;
originally announced June 2018.
-
Projection-Based Finite Elements for Nonlinear Function Spaces
Authors:
Philipp Grohs,
Hanne Hardering,
Oliver Sander,
Markus Sprecher
Abstract:
We introduce a novel type of approximation spaces for functions with values in a nonlinear manifold. The discrete functions are constructed by piecewise polynomial interpolation in a Euclidean embedding space, and then projecting pointwise onto the manifold. We show optimal interpolation error bounds with respect to Lebesgue and Sobolev norms. Additionally, we show similar bounds for the test func…
▽ More
We introduce a novel type of approximation spaces for functions with values in a nonlinear manifold. The discrete functions are constructed by piecewise polynomial interpolation in a Euclidean embedding space, and then projecting pointwise onto the manifold. We show optimal interpolation error bounds with respect to Lebesgue and Sobolev norms. Additionally, we show similar bounds for the test functions, i.e., variations of discrete functions. Combining these results with a nonlinear Céa lemma, we prove optimal $L^2$ and $H^1$ discretization error bounds for harmonic maps from a planar domain into a smooth manifold. All these error bounds are also verified numerically.
△ Less
Submitted 17 March, 2018;
originally announced March 2018.
-
Truncated Nonsmooth Newton Multigrid Methods for Block-Separable Minimization Problems
Authors:
Carsten Gräser,
Oliver Sander
Abstract:
The Truncated Nonsmooth Newton Multigrid (TNNMG) method is a robust and efficient solution method for a wide range of block-separable convex minimization problems, typically stemming from discretizations of nonlinear and nonsmooth partial differential equations. This paper proves global convergence of the method under weak conditions both on the objective functional, and on the local inexact subpr…
▽ More
The Truncated Nonsmooth Newton Multigrid (TNNMG) method is a robust and efficient solution method for a wide range of block-separable convex minimization problems, typically stemming from discretizations of nonlinear and nonsmooth partial differential equations. This paper proves global convergence of the method under weak conditions both on the objective functional, and on the local inexact subproblem solvers that are part of the method. It also discusses a range of algorithmic choices that allows to customize the algorithm for many specific problems. Numerical examples are deliberately omitted, because many such examples have already been published elsewhere.
△ Less
Submitted 25 October, 2017; v1 submitted 14 September, 2017;
originally announced September 2017.
-
A globally convergent filter-trust-region method for large deformation contact problems
Authors:
Jonathan Youett,
Oliver Sander,
Ralf Kornhuber
Abstract:
We present a globally convergent method for the solution of frictionless large deformation contact problems for hyperelastic materials. The discretisation uses the mortar method which is known to be more stable than node-to-segment approaches. The resulting non-convex constrained minimisation problems are solved using a filter-trust-region scheme, and we prove global convergence towards first-orde…
▽ More
We present a globally convergent method for the solution of frictionless large deformation contact problems for hyperelastic materials. The discretisation uses the mortar method which is known to be more stable than node-to-segment approaches. The resulting non-convex constrained minimisation problems are solved using a filter-trust-region scheme, and we prove global convergence towards first-order optimal points. The constrained Newton problems are solved robustly and efficiently using a Truncated Non-smooth Newton Multigrid (TNNMG) method with a Monotone Multigrid (MMG) linear correction step. For this we introduce a cheap basis transformation that decouples the contact constraints. Numerical experiments confirm the stability and efficiency of our approach.
△ Less
Submitted 4 August, 2017;
originally announced August 2017.
-
Solving primal plasticity increment problems in the time of a single predictor-corrector iteration
Authors:
Oliver Sander
Abstract:
The Truncated Nonsmooth Newton Multigrid (TNNMG) method is a well-established method for the solution of strictly convex block-separably nondifferentiable minimization problems. It achieves multigrid-like performance even for non-smooth nonlinear problems, while at the same time being globally convergent and without employing penalty parameters. We show that the algorithm can be applied to the pri…
▽ More
The Truncated Nonsmooth Newton Multigrid (TNNMG) method is a well-established method for the solution of strictly convex block-separably nondifferentiable minimization problems. It achieves multigrid-like performance even for non-smooth nonlinear problems, while at the same time being globally convergent and without employing penalty parameters. We show that the algorithm can be applied to the primal problem of classical linear elastoplasticity with hardening. Numerical experiments show that the method is considerably faster than classical predictor-corrector methods. Indeed, solving an entire increment problem with TNNMG takes less time than a single predictor-corrector iteration for the same problem. Since the algorithm does not rely on differentiability of the objective functional, nonsmooth yield functions like the Tresca yield function can be easily incorporated. The method is closely related to a predictor-corrector scheme with a consistent tangent predictor and line search. We explain the algorithm, prove global convergence, and show its efficiency using a standard benchmark from the literature.
△ Less
Submitted 18 October, 2017; v1 submitted 12 July, 2017;
originally announced July 2017.
-
Simulating Wear On Total Knee Replacements
Authors:
Ansgar Burchardt,
Christian Abicht,
Oliver Sander
Abstract:
Wear on total knee replacements (TKRs) is an important criterion for their performance characteristics. Numerical simulations of such wear has seen increasing attention over the last years. They have the potential to be much faster and less expensive than the in vitro tests in use today. While it is unlikely that in silico tests will replace actual physical tests in the foreseeable future, a judic…
▽ More
Wear on total knee replacements (TKRs) is an important criterion for their performance characteristics. Numerical simulations of such wear has seen increasing attention over the last years. They have the potential to be much faster and less expensive than the in vitro tests in use today. While it is unlikely that in silico tests will replace actual physical tests in the foreseeable future, a judicious combination of both approaches can help making both implant design and pre-clinical testing quicker and more cost-effective.
The challenge today for the design of simulation methods is to obtain results that convey quantitative information, and to do so quickly and reliably. This involves the choice of mathematical models as well as the numerical tools used to solve them. The correctness of the choice can only be validated by comparing with experimental results.
In this paper we present finite element simulations of the wear in TKRs during the gait cycle standardized in the ISO 14243-1 document, used for compliance testing in several countries. As the ISO 14243-1 standard is precisely defined and publicly available, it can serve as an excellent benchmark for comparison of wear simulation methods. Our novel contact algorithm works without Lagrange multipliers and penalty methods, achieving unparalleled stability and efficiency. We compare our simulation results with the experimental data from physical tests using two different actual TKRs, each test being performed three times. We can closely predict the total mass loss due to wear after five million gait cycles. We also observe a good match between the wear patterns seen in experiments and our simulation results.
△ Less
Submitted 3 April, 2017;
originally announced April 2017.
-
Test Function Spaces for Geometric Finite Elements
Authors:
Oliver Sander
Abstract:
We construct test function spaces for geometric finite elements. Geometric finite elements (GFE) are generalizations of Lagrangian finite elements to situations where the unknown function maps into a nonlinear space. Test functions for such spaces arise as variations of GFE functions wherever the GFE function space has a local manifold structure. For any given GFE function $u_h$, the test function…
▽ More
We construct test function spaces for geometric finite elements. Geometric finite elements (GFE) are generalizations of Lagrangian finite elements to situations where the unknown function maps into a nonlinear space. Test functions for such spaces arise as variations of GFE functions wherever the GFE function space has a local manifold structure. For any given GFE function $u_h$, the test functions form a linear space that depends on $u_h$. They generalize Jacobi fields in the same way that the GFE interpolation functions generalize geodesic curves. Having test function spaces allows to extend the GFE method to boundary value problems that do not have a minimization formulation.
△ Less
Submitted 25 July, 2016;
originally announced July 2016.
-
Challenges in QCD matter physics - The Compressed Baryonic Matter experiment at FAIR
Authors:
CBM Collaboration,
T. Ablyazimov,
A. Abuhoza,
R. P. Adak,
M. Adamczyk,
K. Agarwal,
M. M. Aggarwal,
Z. Ahammed,
F. Ahmad,
N. Ahmad,
S. Ahmad,
A. Akindinov,
P. Akishin,
E. Akishina,
T. Akishina,
V. Akishina,
A. Akram,
M. Al-Turany,
I. Alekseev,
E. Alexandrov,
I. Alexandrov,
S. Amar-Youcef,
M. Anđelić,
O. Andreeva,
C. Andrei
, et al. (563 additional authors not shown)
Abstract:
Substantial experimental and theoretical efforts worldwide are devoted to explore the phase diagram of strongly interacting matter. At LHC and top RHIC energies, QCD matter is studied at very high temperatures and nearly vanishing net-baryon densities. There is evidence that a Quark-Gluon-Plasma (QGP) was created at experiments at RHIC and LHC. The transition from the QGP back to the hadron gas is…
▽ More
Substantial experimental and theoretical efforts worldwide are devoted to explore the phase diagram of strongly interacting matter. At LHC and top RHIC energies, QCD matter is studied at very high temperatures and nearly vanishing net-baryon densities. There is evidence that a Quark-Gluon-Plasma (QGP) was created at experiments at RHIC and LHC. The transition from the QGP back to the hadron gas is found to be a smooth cross over. For larger net-baryon densities and lower temperatures, it is expected that the QCD phase diagram exhibits a rich structure, such as a first-order phase transition between hadronic and partonic matter which terminates in a critical point, or exotic phases like quarkyonic matter. The discovery of these landmarks would be a breakthrough in our understanding of the strong interaction and is therefore in the focus of various high-energy heavy-ion research programs. The Compressed Baryonic Matter (CBM) experiment at FAIR will play a unique role in the exploration of the QCD phase diagram in the region of high net-baryon densities, because it is designed to run at unprecedented interaction rates. High-rate operation is the key prerequisite for high-precision measurements of multi-differential observables and of rare diagnostic probes which are sensitive to the dense phase of the nuclear fireball. The goal of the CBM experiment at SIS100 (sqrt(s_NN) = 2.7 - 4.9 GeV) is to discover fundamental properties of QCD matter: the phase structure at large baryon-chemical potentials (mu_B > 500 MeV), effects of chiral symmetry, and the equation-of-state at high density as it is expected to occur in the core of neutron stars. In this article, we review the motivation for and the physics programme of CBM, including activities before the start of data taking in 2022, in the context of the worldwide efforts to explore high-density QCD matter.
△ Less
Submitted 29 March, 2017; v1 submitted 6 July, 2016;
originally announced July 2016.
-
Simulation of Deformation and Flow in Fractured, Poroelastic Materials
Authors:
Katja K. Hanowski,
Oliver Sander
Abstract:
We introduce a coupled system of PDEs for the modeling of the fluid-fluid and fluid-solid interaction in a poroelastic material with a single static fracture. The fluid flow in the fracture is modeled by a lower-dimensional Darcy equation, which interacts with the surrounding rock matrix and the fluid it contains. We explicitly allow the fracture to end within the domain, and the fracture width is…
▽ More
We introduce a coupled system of PDEs for the modeling of the fluid-fluid and fluid-solid interaction in a poroelastic material with a single static fracture. The fluid flow in the fracture is modeled by a lower-dimensional Darcy equation, which interacts with the surrounding rock matrix and the fluid it contains. We explicitly allow the fracture to end within the domain, and the fracture width is an unknown of the problem. The resulting weak problem is nonlinear, elliptic and symmetric, and can be given the structure of a fixed-point problem. We show that the coupled fluid-fluid problem has a solution in a specially crafted Sobolev space, even though the fracture width cannot be bounded away from zero near the crack tip.
For numerical simulations, we combine XFEM discretizations for the rock matrix deformation and pore pressure with a standard lower-dimensional finite element method for the fracture flow problem. The resulting coupled discrete system consists of linear subdomain problems coupled by nonlinear coupling conditions. We solve the coupled system with a substructuring solver and observe very fast convergence. We also observe optimal mesh dependence of the discretization errors even in the presence of crack tips.
△ Less
Submitted 12 August, 2016; v1 submitted 18 June, 2016;
originally announced June 2016.
-
The interface for functions in the dune-functions module
Authors:
Christian Engwer,
Carsten Gräser,
Steffen Müthing,
Oliver Sander
Abstract:
The dune-functions dune module introduces a new programmer interface for discrete and non-discrete functions. Unlike the previous interfaces considered in the existing dune modules, it is based on overloading operator(), and returning values by-value. This makes user code much more readable, and allows the incorporation of newer C++ features such as lambda expressions. Run-time polymorphism is imp…
▽ More
The dune-functions dune module introduces a new programmer interface for discrete and non-discrete functions. Unlike the previous interfaces considered in the existing dune modules, it is based on overloading operator(), and returning values by-value. This makes user code much more readable, and allows the incorporation of newer C++ features such as lambda expressions. Run-time polymorphism is implemented not by inheritance, but by type erasure, generalizing the ideas of the std::function class from the C++11 standard library. We describe the new interface, show its possibilities, and measure the performance impact of type erasure and return-by-value.
△ Less
Submitted 18 December, 2015;
originally announced December 2015.
-
The Dune FoamGrid implementation for surface and network grids
Authors:
Oliver Sander,
Timo Koch,
Natalie Schröder,
Bernd Flemisch
Abstract:
We present FoamGrid, a new implementation of the DUNE grid interface. FoamGrid implements one- and two-dimensional grids in a physical space of arbitrary dimension, which allows for grids for curved domains. Even more, the grids are not expected to have a manifold structure, i.e., more than two elements can share a common facet. This makes FoamGrid the grid data structure of choice for simulating…
▽ More
We present FoamGrid, a new implementation of the DUNE grid interface. FoamGrid implements one- and two-dimensional grids in a physical space of arbitrary dimension, which allows for grids for curved domains. Even more, the grids are not expected to have a manifold structure, i.e., more than two elements can share a common facet. This makes FoamGrid the grid data structure of choice for simulating structures such as foams, discrete fracture networks, or network flow problems. FoamGrid implements adaptive non-conforming refinement with element parametrizations. As an additional feature it allows removal and addition of elements in an existing grid, which makes FoamGrid suitable for network growth problems. We show how to use FoamGrid, with particular attention to the extensions of the grid interface needed to handle non-manifold topology and grid growth. Three numerical examples demonstrate the possibilities offered by FoamGrid.
△ Less
Submitted 11 November, 2015;
originally announced November 2015.
-
Numerical Treatment of a Geometrically Nonlinear Planar Cosserat Shell Model
Authors:
Oliver Sander,
Patrizio Neff,
Mircea Bîrsan
Abstract:
We present a new way to discretize a geometrically nonlinear elastic planar Cosserat shell. The kinematical model is similar to the general 6-parameter resultant shell model with drilling rotations. The discretization uses geodesic finite elements, which leads to an objective discrete model which naturally allows arbitrarily large rotations. Finite elements of any approximation order can be constr…
▽ More
We present a new way to discretize a geometrically nonlinear elastic planar Cosserat shell. The kinematical model is similar to the general 6-parameter resultant shell model with drilling rotations. The discretization uses geodesic finite elements, which leads to an objective discrete model which naturally allows arbitrarily large rotations. Finite elements of any approximation order can be constructed. The resulting algebraic problem is a minimization problem posed on a nonlinear finite-dimensional Riemannian manifold. We solve this problem using a Riemannian trust-region method, which is a generalization of Newton's method that converges globally without intermediate loading steps. We present the continuous model and the discretization, discuss the properties of the discrete model, and show several numerical examples, including wrinkles of thin elastic sheets in shear.
△ Less
Submitted 28 January, 2015; v1 submitted 11 December, 2014;
originally announced December 2014.
-
Unsaturated subsurface flow with surface water and nonlinear in- and outflow conditions
Authors:
Heiko Berninger,
Mario Ohlberger,
Oliver Sander,
Kathrin Smetana
Abstract:
We analytically and numerically analyze groundwater flow in a homogeneous soil described by the Richards equation, coupled to surface water represented by a set of ordinary differential equations (ODE's) on parts of the domain boundary, and with nonlinear outflow conditions of Signorini's type. The coupling of the partial differential equation (PDE) and the ODE's is given by nonlinear Robin bounda…
▽ More
We analytically and numerically analyze groundwater flow in a homogeneous soil described by the Richards equation, coupled to surface water represented by a set of ordinary differential equations (ODE's) on parts of the domain boundary, and with nonlinear outflow conditions of Signorini's type. The coupling of the partial differential equation (PDE) and the ODE's is given by nonlinear Robin boundary conditions. This article provides two major new contributions regarding these infiltration conditions. First, an existence result for the continuous coupled problem is established with the help of a regularization technique. Second, we analyze and validate a solver-friendly discretization of the coupled problem based on an implicit-explicit time discretization and on finite elements in space. The discretized PDE leads to convex spatial minimization problems which can be solved efficiently by monotone multigrid. Numerical experiments are provided using the DUNE numerics framework.
△ Less
Submitted 11 January, 2013;
originally announced January 2013.