-
Parallel-in-time solution of hyperbolic PDE systems via characteristic-variable block preconditioning
Authors:
H. De Sterck,
R. D. Falgout,
O. A. Krzysik,
J. B. Schroder
Abstract:
We consider the parallel-in-time solution of hyperbolic partial differential equation (PDE) systems in one spatial dimension, both linear and nonlinear. In the nonlinear setting, the discretized equations are solved with a preconditioned residual iteration based on a global linearization. The linear(ized) equation systems are approximately solved parallel-in-time using a block preconditioner appli…
▽ More
We consider the parallel-in-time solution of hyperbolic partial differential equation (PDE) systems in one spatial dimension, both linear and nonlinear. In the nonlinear setting, the discretized equations are solved with a preconditioned residual iteration based on a global linearization. The linear(ized) equation systems are approximately solved parallel-in-time using a block preconditioner applied in the characteristic variables of the underlying linear(ized) hyperbolic PDE. This change of variables is motivated by the observation that inter-variable coupling for characteristic variables is weak relative to intra-variable coupling, at least locally where spatio-temporal variations in the eigenvectors of the associated flux Jacobian are sufficiently small. For an $\ell$-dimensional system of PDEs, applying the preconditioner consists of solving a sequence of $\ell$ scalar linear(ized)-advection-like problems, each being associated with a different characteristic wave-speed in the underlying linear(ized) PDE. We approximately solve these linear advection problems using multigrid reduction-in-time (MGRIT); however, any other suitable parallel-in-time method could be used. Numerical examples are shown for the (linear) acoustics equations in heterogeneous media, and for the (nonlinear) shallow water equations and Euler equations of gas dynamics with shocks and rarefactions.
△ Less
Submitted 4 July, 2024;
originally announced July 2024.
-
RMF: A Risk Measurement Framework for Machine Learning Models
Authors:
Jan Schröder,
Jakub Breier
Abstract:
Machine learning (ML) models are used in many safety- and security-critical applications nowadays. It is therefore important to measure the security of a system that uses ML as a component. This paper focuses on the field of ML, particularly the security of autonomous vehicles. For this purpose, a technical framework will be described, implemented, and evaluated in a case study. Based on ISO/IEC 2…
▽ More
Machine learning (ML) models are used in many safety- and security-critical applications nowadays. It is therefore important to measure the security of a system that uses ML as a component. This paper focuses on the field of ML, particularly the security of autonomous vehicles. For this purpose, a technical framework will be described, implemented, and evaluated in a case study. Based on ISO/IEC 27004:2016, risk indicators are utilized to measure and evaluate the extent of damage and the effort required by an attacker. It is not possible, however, to determine a single risk value that represents the attacker's effort. Therefore, four different values must be interpreted individually.
△ Less
Submitted 15 June, 2024;
originally announced June 2024.
-
Computing $\vec{\mathcal{S}}$-DAGs and Parity Games
Authors:
Meike Hatzel,
Johannes Schröder
Abstract:
Treewidth on undirected graphs is known to have many algorithmic applications. When considering directed width-measures there are much less results on their deployment for algorithmic results. In 2022 the first author, Rabinovich and Wiederrecht introduced a new directed width measure, $\vec{\mathcal{S}}$-DAG-width, using directed separations and obtained a structural duality for it. In 2012 Berwa…
▽ More
Treewidth on undirected graphs is known to have many algorithmic applications. When considering directed width-measures there are much less results on their deployment for algorithmic results. In 2022 the first author, Rabinovich and Wiederrecht introduced a new directed width measure, $\vec{\mathcal{S}}$-DAG-width, using directed separations and obtained a structural duality for it. In 2012 Berwanger~et~al.~solved Parity Games in polynomial time on digraphs of bounded DAG-width. With generalising this result to digraphs of bounded $\vec{\mathcal{S}}$-DAG-width and also providing an algorithm to compute the $\vec{\mathcal{S}}$-DAG-width of a given digraphs we give first algorithmical results for this new parameter.
△ Less
Submitted 9 May, 2024;
originally announced May 2024.
-
Theoretical Unification of the Fractured Aspects of Information
Authors:
Marcin J. Schroeder
Abstract:
The article has as its main objective the identification of fundamental epistemological obstacles in the study of information related to unnecessary methodological assumptions and the demystification of popular beliefs in the fundamental divisions of the aspects of information that can be understood as Bachelardian rupture of epistemological obstacles. These general considerations are preceded by…
▽ More
The article has as its main objective the identification of fundamental epistemological obstacles in the study of information related to unnecessary methodological assumptions and the demystification of popular beliefs in the fundamental divisions of the aspects of information that can be understood as Bachelardian rupture of epistemological obstacles. These general considerations are preceded by an overview of the motivations for the study of information and the role of the concept of information in the conceptualization of intelligence, complexity, and consciousness justifying the need for a sufficiently general perspective in the study of information, and are followed at the end of the article by a brief exposition of an example of a possible application in the development of the unified theory of information free from unnecessary divisions and claims of superiority of the existing preferences in methodology. The reference to Gaston Bachelard and his ideas of epistemological obstacles and epistemological ruptures seems highly appropriate for the reflection on the development of information study, in particular in the context of obstacles such as the absence of semantics of information, negligence of its structural analysis, separation of its digital and analog forms, and misguided use of mathematics.
△ Less
Submitted 26 February, 2024;
originally announced February 2024.
-
A computational approach to identify the material parameters of the relaxed micromorphic model
Authors:
Mohammad Sarhil,
Lisa Scheunemann,
Peter Lewintan,
Jörg Schröder,
Patrizio Neff
Abstract:
We determine the material parameters in the relaxed micromorphic generalized continuum model for a given periodic microstructure in this work. This is achieved through a least squares fitting of the total energy of the relaxed micromorphic homogeneous continuum to the total energy of the fully-resolved heterogeneous microstructure, governed by classical linear elasticity. The relaxed micromorphic…
▽ More
We determine the material parameters in the relaxed micromorphic generalized continuum model for a given periodic microstructure in this work. This is achieved through a least squares fitting of the total energy of the relaxed micromorphic homogeneous continuum to the total energy of the fully-resolved heterogeneous microstructure, governed by classical linear elasticity. The relaxed micromorphic model is a generalized continuum that utilizes the $\Curl$ of a micro-distortion field instead of its full gradient as in the classical micromorphic theory, leading to several advantages and differences. The most crucial advantage is that it operates between two well-defined scales. These scales are determined by linear elasticity with microscopic and macroscopic elasticity tensors, which respectively bound the stiffness of the relaxed micromorphic continuum from above and below. While the macroscopic elasticity tensor is established a priori through standard periodic first-order homogenization, the microscopic elasticity tensor remains to be determined. Additionally, the characteristic length parameter, associated with curvature measurement, controls the transition between the micro- and macro-scales. Both the microscopic elasticity tensor and the characteristic length parameter are here determined using a computational approach based on the least squares fitting of energies. This process involves the consideration of an adequate number of quadratic deformation modes and different specimen sizes. We conduct a comparative analysis between the least square fitting results of the relaxed micromorphic model, the fitting of a skew-symmetric micro-distortion field (Cosserat-micropolar model), and the fitting of the classical micromorphic model with two different formulations for the curvature...
△ Less
Submitted 28 January, 2024;
originally announced January 2024.
-
Generalized Optimal AMG Convergence Theory for Nonsymmetric and Indefinite Problems
Authors:
Ahsan Ali,
James Brannick,
Karsten Kahl,
Oliver A. Krzysik,
Jacob B. Schroder,
Ben S. Southworth
Abstract:
Algebraic multigrid (AMG) is known to be an effective solver for many sparse symmetric positive definite (SPD) linear systems. For SPD systems, the convergence theory of AMG is well-understood in terms of the $A$-norm but in a nonsymmetric setting such an energy norm is non-existent. For this reason, convergence of AMG for nonsymmetric systems of equations remains an open area of research. Existin…
▽ More
Algebraic multigrid (AMG) is known to be an effective solver for many sparse symmetric positive definite (SPD) linear systems. For SPD systems, the convergence theory of AMG is well-understood in terms of the $A$-norm but in a nonsymmetric setting such an energy norm is non-existent. For this reason, convergence of AMG for nonsymmetric systems of equations remains an open area of research. Existing nonsymmetric AMG algorithms in this setting mostly rely on heuristics motivated by SPD convergence theory. In the SPD setting, the classical form of optimal AMG interpolation provides a useful insight in determining the two grid convergence rate of the method. In this work, we discuss a generalization of the optimal AMG convergence theory targeting nonsymmetric problems by constructing a $2\times 2$ block symmetric indefinite system so that the Petrov-Galerkin AMG process for the nonsymmetric matrix $A$ can be recast as a Galerkin AMG process for a symmetric indefinite system. We show that using this generalization of the optimal interpolation theory, one can obtain the same identity for the two-grid convergence rate as that derived in the SPD setting for optimal interpolation. We also provide supporting numerical results for the convergence result and nonsymmetric advection-diffusion problems.
△ Less
Submitted 20 January, 2024;
originally announced January 2024.
-
Parallel-in-time solution of scalar nonlinear conservation laws
Authors:
H. De Sterck,
R. D. Falgout,
O. A. Krzysik,
J. B. Schroder
Abstract:
We consider the parallel-in-time solution of scalar nonlinear conservation laws in one spatial dimension. The equations are discretized in space with a conservative finite-volume method using weighted essentially non-oscillatory (WENO) reconstructions, and in time with high-order explicit Runge-Kutta methods. The solution of the global, discretized space-time problem is sought via a nonlinear iter…
▽ More
We consider the parallel-in-time solution of scalar nonlinear conservation laws in one spatial dimension. The equations are discretized in space with a conservative finite-volume method using weighted essentially non-oscillatory (WENO) reconstructions, and in time with high-order explicit Runge-Kutta methods. The solution of the global, discretized space-time problem is sought via a nonlinear iteration that uses a novel linearization strategy in cases of non-differentiable equations. Under certain choices of discretization and algorithmic parameters, the nonlinear iteration coincides with Newton's method, although, more generally, it is a preconditioned residual correction scheme. At each nonlinear iteration, the linearized problem takes the form of a certain discretization of a linear conservation law over the space-time domain in question. An approximate parallel-in-time solution of the linearized problem is computed with a single multigrid reduction-in-time (MGRIT) iteration. The MGRIT iteration employs a novel coarse-grid operator that is a modified conservative semi-Lagrangian discretization and generalizes those we have developed previously for non-conservative scalar linear hyperbolic problems. Numerical tests are performed for the inviscid Burgers and Buckley--Leverett equations. For many test problems, the solver converges in just a handful of iterations with convergence rate independent of mesh resolution, including problems with (interacting) shocks and rarefactions.
△ Less
Submitted 10 January, 2024;
originally announced January 2024.
-
Blind Frequency-Domain Equalization Using Vector-Quantized Variational Autoencoders
Authors:
Jinxiang Song,
Vincent Lauinger,
Christian Häger,
Jochen Schröder,
Alexandre Graell i Amat,
Laurent Schmalen,
Henk Wymeersch
Abstract:
We propose a novel frequency-domain blind equalization scheme for coherent optical communications. The method is shown to achieve similar performance to its recently proposed time-domain counterpart with lower computational complexity, while outperforming the commonly used CMA-based equalizers.
We propose a novel frequency-domain blind equalization scheme for coherent optical communications. The method is shown to achieve similar performance to its recently proposed time-domain counterpart with lower computational complexity, while outperforming the commonly used CMA-based equalizers.
△ Less
Submitted 26 December, 2023;
originally announced December 2023.
-
Numerical realization of the Mortensen observer via a Hessian-augmented polynomial approximation of the value function
Authors:
Tobias Breiten,
Karl Kunisch,
Jesper Schröder
Abstract:
Two related numerical schemes for the realization of the Mortensen observer or minimum energy estimator for the state reconstruction of non-linear dynamical systems subject to deterministic disturbances are proposed and compared. Both approaches rely on a polynomial approximation of the value function associated with the energy of the disturbances of the system. Such an approximation is obtained v…
▽ More
Two related numerical schemes for the realization of the Mortensen observer or minimum energy estimator for the state reconstruction of non-linear dynamical systems subject to deterministic disturbances are proposed and compared. Both approaches rely on a polynomial approximation of the value function associated with the energy of the disturbances of the system. Such an approximation is obtained via interpolation considering not only the values but also first and second order derivatives of the value function in a set of sampling points. The scheme is applied to four examples and the results are compared with the well known extended Kalman filter.
△ Less
Submitted 8 November, 2023; v1 submitted 31 October, 2023;
originally announced October 2023.
-
Elasto-plastic residual stress analysis of selective laser sintered porous materials based on 3D-multilayer thermo-structural phase-field simulations
Authors:
Yangyiwei Yang,
Somnath Bharech,
Nick Finger,
Xiandong Zhou,
Joerg Schroeder,
Bai-Xiang Xu
Abstract:
Residual stress and plastic strain in additive manufactured materials can exhibit significant microscopic variation at the powder scale, profoundly influencing the overall properties of printed components. This variation depends on processing parameters and stems from multiple factors, including differences in powder bed morphology, non-uniform thermo-structural profiles, and inter-layer fusion. I…
▽ More
Residual stress and plastic strain in additive manufactured materials can exhibit significant microscopic variation at the powder scale, profoundly influencing the overall properties of printed components. This variation depends on processing parameters and stems from multiple factors, including differences in powder bed morphology, non-uniform thermo-structural profiles, and inter-layer fusion. In this research, we propose a powder-resolved multilayer multiphysics simulation scheme tailored for porous materials through the process of selective laser sintering. This approach seamlessly integrates finite element method (FEM)-based non-isothermal phase-field simulation with thermo-elasto-plastic simulation, incorporating temperature- and phase-dependent material properties. The outcome of this investigation includes a detailed depiction of the mesoscopic evolution of stress and plastic strain within a transient thermo-microstructure, evaluated across a spectrum of beam power and scan speed parameters. Simulation results further reveal the underlying mechanisms. For instance, stress concentration primarily occurs at the necking region of partially melted particles and the junctions between different layers, resulting in the accumulation of plastic strain and residual stress, ultimately leading to structural distortion in the materials. Based on the simulation data, phenomenological relation regarding porosity/densification control by the beam energy input was examined along with the comparison to experimental results. Regression models were also proposed to describe the dependency of the residual stress and the plastic strain on the beam energy input.
△ Less
Submitted 4 April, 2024; v1 submitted 20 October, 2023;
originally announced October 2023.
-
Space-Time Block Preconditioning for Incompressible Resistive Magnetohydrodynamics
Authors:
Federico Danieli,
Ben S. Southworth,
Jacob B. Schroder
Abstract:
This work develops a novel all-at-once space-time preconditioning approach for resistive magnetohydrodynamics (MHD), with a focus on model problems targeting fusion reactor design. We consider parallel-in-time due to the long time domains required to capture the physics of interest, as well as the complexity of the underlying system and thereby computational cost of long-time integration. To ameli…
▽ More
This work develops a novel all-at-once space-time preconditioning approach for resistive magnetohydrodynamics (MHD), with a focus on model problems targeting fusion reactor design. We consider parallel-in-time due to the long time domains required to capture the physics of interest, as well as the complexity of the underlying system and thereby computational cost of long-time integration. To ameliorate this cost by using many processors, we thus develop a novel approach to solving the whole space-time system that is parallelizable in both space and time. We develop a space-time block preconditioning for resistive MHD, following the space-time block preconditioning concept first introduced by Danieli et al. in 2022 for incompressible flow, where an effective preconditioner for classic sequential time-stepping is extended to the space-time setting. The starting point for our derivation is the continuous Schur complement preconditioner by Cyr et al. in 2021, which we proceed to generalise in order to produce, to our knowledge, the first space-time block preconditioning approach for the challenging equations governing incompressible resistive MHD. The numerical results are promising for the model problems of island coalescence and tearing mode, with the overhead computational cost associated with space-time preconditioning versus sequential time-stepping being modest and primarily in the range of 2x-5x, which is low for parallel-in-time schemes in general. Additionally, the scaling results for inner (linear) and outer (nonlinear) iterations are flat in the case of fixed time-step size and only grow very slowly in the case of time-step refinement.
△ Less
Submitted 1 September, 2023;
originally announced September 2023.
-
Advanced methods for analyzing in-situ observations of magnetic reconnection
Authors:
H. Hasegawa,
M. R. Argall,
N. Aunai,
R. Bandyopadhyay,
N. Bessho,
I. J. Cohen,
R. E. Denton,
J. C. Dorelli,
J. Egedal,
S. A. Fuselier,
P. Garnier,
V. Genot,
D. B. Graham,
K. J. Hwang,
Y. V. Khotyaintsev,
D. B. Korovinskiy,
B. Lavraud,
Q. Lenouvel,
T. C. Li,
Y. -H. Liu,
B. Michotte de Welle,
T. K. M. Nakamura,
D. S. Payne,
S. M. Petrinec,
Y. Qi
, et al. (11 additional authors not shown)
Abstract:
There is ample evidence for magnetic reconnection in the solar system, but it is a nontrivial task to visualize, to determine the proper approaches and frames to study, and in turn to elucidate the physical processes at work in reconnection regions from in-situ measurements of plasma particles and electromagnetic fields. Here an overview is given of a variety of single- and multi-spacecraft data a…
▽ More
There is ample evidence for magnetic reconnection in the solar system, but it is a nontrivial task to visualize, to determine the proper approaches and frames to study, and in turn to elucidate the physical processes at work in reconnection regions from in-situ measurements of plasma particles and electromagnetic fields. Here an overview is given of a variety of single- and multi-spacecraft data analysis techniques that are key to revealing the context of in-situ observations of magnetic reconnection in space and for detecting and analyzing the diffusion regions where ions and/or electrons are demagnetized. We focus on recent advances in the era of the Magnetospheric Multiscale mission, which has made electron-scale, multi-point measurements of magnetic reconnection in and around Earth's magnetosphere.
△ Less
Submitted 24 June, 2024; v1 submitted 11 July, 2023;
originally announced July 2023.
-
Real-Time Monitoring of Cable Break in a Live Fiber Network using a Coherent Transceiver Prototype
Authors:
Mikael Mazur,
Lauren Dallachiesa,
Roland Ryf,
Dennis Wallberg,
Erik Börjesson,
Magnus Bergroth,
Börje Josefsson,
Nicolas K. Fontaine,
Haoshuo Chen,
David T. Neilson,
Jochen Schröder,
Per Larsson-Edefors,
Magnus Karlsson
Abstract:
We monitor a 524km live network link using an FPGA-based sensing-capable coherent transceiver prototype during a human-caused cable break. Post-analysis of polarization data reveals minute-level potential warning precursors and baseline-exceeding changes directly preceding the break.
We monitor a 524km live network link using an FPGA-based sensing-capable coherent transceiver prototype during a human-caused cable break. Post-analysis of polarization data reveals minute-level potential warning precursors and baseline-exceeding changes directly preceding the break.
△ Less
Submitted 3 July, 2023;
originally announced July 2023.
-
Constrained Local Approximate Ideal Restriction for Advection-Diffusion Problems
Authors:
Ahsan Ali,
James Brannick,
Karsten Kahl,
Oliver A. Krzysik,
Jacob B. Schroder,
Ben S. Southworth
Abstract:
This paper focuses on developing a reduction-based algebraic multigrid method that is suitable for solving general (non)symmetric linear systems and is naturally robust from pure advection to pure diffusion. Initial motivation comes from a new reduction-based algebraic multigrid (AMG) approach, $\ell$AIR (local approximate ideal restriction), that was developed for solving advection-dominated prob…
▽ More
This paper focuses on developing a reduction-based algebraic multigrid method that is suitable for solving general (non)symmetric linear systems and is naturally robust from pure advection to pure diffusion. Initial motivation comes from a new reduction-based algebraic multigrid (AMG) approach, $\ell$AIR (local approximate ideal restriction), that was developed for solving advection-dominated problems. Though this new solver is very effective in the advection dominated regime, its performance degrades in cases where diffusion becomes dominant. This is consistent with the fact that in general, reduction-based AMG methods tend to suffer from growth in complexity and/or convergence rates as the problem size is increased, especially for diffusion dominated problems in two or three dimensions. Motivated by the success of $\ell$AIR in the advective regime, our aim in this paper is to generalize the AIR framework with the goal of improving the performance of the solver in diffusion dominated regimes. To do so, we propose a novel way to combine mode constraints as used commonly in energy minimization AMG methods with the local approximation of ideal operators used in $\ell$AIR. The resulting constrained $\ell$AIR (C$\ell$AIR) algorithm is able to achieve fast scalable convergence on advective and diffusive problems. In addition, it is able to achieve standard low complexity hierarchies in the diffusive regime through aggressive coarsening, something that has been previously difficult for reduction-based methods.
△ Less
Submitted 14 May, 2024; v1 submitted 1 July, 2023;
originally announced July 2023.
-
Local well-posedness of the Mortensen observer
Authors:
Tobias Breiten,
Jesper Schröder
Abstract:
The analytical background of nonlinear observers based on minimal energy estimation is discussed. It is shown that locally the derivation of the observer equation based on a trajectory with pointwise minimal energy can be done rigorously. The result is obtained by a local sensitivity analysis of the value function based on Pontryagin's maximum principle and the Hamilton-Jacobi-Bellman equation. Th…
▽ More
The analytical background of nonlinear observers based on minimal energy estimation is discussed. It is shown that locally the derivation of the observer equation based on a trajectory with pointwise minimal energy can be done rigorously. The result is obtained by a local sensitivity analysis of the value function based on Pontryagin's maximum principle and the Hamilton-Jacobi-Bellman equation. The consideration of a differential Riccati equation reveals that locally the second derivative of the value function is a positive definite matrix. The local convexity ensures existence of a trajectory minimizing the energy, which is then shown to satisfy the observer equation.
△ Less
Submitted 16 May, 2023;
originally announced May 2023.
-
High-Fidelity Transport of Trapped-Ion Qubits in a Multi-Layer Array
Authors:
Deviprasath Palani,
Florian Hasse,
Philip Kiefer,
Frederick Boeckling,
Jan-Philipp Schroeder,
Ulrich Warring,
Tobias Schaetz
Abstract:
A variety of physical platforms are investigated for quantum control of many particles, and techniques are extended to access multiple dimensions. Here, we present our experimental study of shuttling single Mg$^+$ ions within a scalable trap-array architecture that contains up to thirteen trapping sites in a three-dimensional arrangement. We shuttle ions from a dedicated loading hub to multiple si…
▽ More
A variety of physical platforms are investigated for quantum control of many particles, and techniques are extended to access multiple dimensions. Here, we present our experimental study of shuttling single Mg$^+$ ions within a scalable trap-array architecture that contains up to thirteen trapping sites in a three-dimensional arrangement. We shuttle ions from a dedicated loading hub to multiple sites with a success rate of larger than $0.99999$. In a prototype application, we demonstrate the preservation of the coherence of superposition states of a hyperfine qubit during inter-site shuttling. Our findings highlight the potential of these techniques for use in future large-scale architectures.
△ Less
Submitted 9 May, 2023;
originally announced May 2023.
-
A surrogate model for data-driven magnetic stray field calculations
Authors:
Rainer Niekamp,
Johanna Niemann,
Maximilian Reichel,
Hongbin Zhang,
Jörg Schröder
Abstract:
In this contribution we propose a data-driven surrogate model for the prediction of magnetic stray fields in two-dimensional random micro-heterogeneous materials. Since data driven models require thousands of training data sets, FEM simulations appear to be too time consuming. Hence, a stochastic model based on Brownian motion, which utilizes an efficient evaluation of stochastic transition matric…
▽ More
In this contribution we propose a data-driven surrogate model for the prediction of magnetic stray fields in two-dimensional random micro-heterogeneous materials. Since data driven models require thousands of training data sets, FEM simulations appear to be too time consuming. Hence, a stochastic model based on Brownian motion, which utilizes an efficient evaluation of stochastic transition matrices, is applied for the training data generation. For the encoding of the microstructure and the optimization of the surrogate model, two architectures are compared, i.e. the so-called UResNet model and the Fourier Convolutional neural network (FCNN). Here we analyze two FCNNs, one based on the discrete cosine transformation and one based on the complex-valued discrete Fourier transformation. Finally, we compare the magnetic stray fields for independent microstructures (not used in the training set) with results from the FE$^2$ method, a numerical homogenization scheme, to demonstrate the efficiency of the proposed surrogate model.
△ Less
Submitted 8 April, 2023;
originally announced April 2023.
-
Modeling martensitic transformation in shape memory alloys using multi-phase-field elasticity models based on partial rank-one energy relaxation on pairwise interfaces
Authors:
Mohammad Sarhil,
Oleg Shchyglo,
Dominik Brands,
Jörg Schröder,
Ingo Steinbach
Abstract:
To model the mechanically-driven phase transformations, e.g. martensitic transformation, using the phase-field theory, suitable models are needed for describing the mechanical fields of the individual non-vanishing phase-fields in the interface regions in order to obtain the mechanical driving forces of phase-field motion. Quantitative modeling requires satisfying the interfacial static equilibriu…
▽ More
To model the mechanically-driven phase transformations, e.g. martensitic transformation, using the phase-field theory, suitable models are needed for describing the mechanical fields of the individual non-vanishing phase-fields in the interface regions in order to obtain the mechanical driving forces of phase-field motion. Quantitative modeling requires satisfying the interfacial static equilibrium and kinematic compatibility conditions which have already been achieved in the literature for dual-phase-field materials by using the rank-one relaxation (or convexification) of the energy density. A direct generalization to the multi-phase-field case is not applicable without breaking these conditions partially. To the best of our knowledge, no existing multi-phase-field elasticity model has been able to satisfy the jump conditions between all the locally-active phase-fields on their pairwise normals in triple and higher-order junctions. In this work, we introduce a novel multi-phase-field elasticity model based on the partial rank-one relaxation of the elastic energy density defined on the pairwise interfaces...... (see PDF for the rest of the abstract)
△ Less
Submitted 5 April, 2023;
originally announced April 2023.
-
Compact lithium niobate photonic integrated circuits
Authors:
Yan Gao,
Fuchuan Lei,
Marcello Girardi,
Zhichao Ye,
Raphaël Van Laer,
Victor Torres-Company,
Jochen Schröder
Abstract:
Lithium niobate (LN) is a promising material for future complex photonic-electronic circuits, with wide applications in fields like communications, sensing, quantum optics, and computation. LN took a great stride toward compact photonic integrated circuits (PICs) with the development of partially-etched LN on insulator (LNOI) waveguides. However, integration density is still limited for future hig…
▽ More
Lithium niobate (LN) is a promising material for future complex photonic-electronic circuits, with wide applications in fields like communications, sensing, quantum optics, and computation. LN took a great stride toward compact photonic integrated circuits (PICs) with the development of partially-etched LN on insulator (LNOI) waveguides. However, integration density is still limited for future high-compact PICs due to the partial edge nature of their waveguides. Here, we demonstrate a fully-etched LN PIC platform which, for the first time, simultaneously achieves ultra-low propagation loss and compact circuit size. The tightly-confined fully-etched LN waveguides with smooth sidewalls allow us to bring the bending radius down to 20 $μ$m (corresponds to 1 THz FSR). We have achieved compact high-$Q$ microring resonators with $Q/V$ of 7.1 $\times$ 10$^{4}$ $μ$m$^{-3}$, almost one order of magnitude larger than previous demonstrations. The statistical mean propagation losses of our LN waveguides is 8.5 dB/m (corresponds to mean $Q$-factor of 4.9 $\times$ 10$^{6}$) even with a small bending radius of 40 $μ$m. Our compact and ultra-low-loss LN platform shows great potential in future miniaturized multifunctional integration systems. As complementary evidence to show the utility of our platform, we demonstrate soliton microcombs with an ultra-high repetition rate of 500 GHz in LN.
△ Less
Submitted 2 March, 2023;
originally announced March 2023.
-
Blind Channel Equalization Using Vector-Quantized Variational Autoencoders
Authors:
Jinxiang Song,
Vincent Lauinger,
Yibo Wu,
Christian Häger,
Jochen Schröder,
Alexandre Graell i Amat,
Laurent Schmalen,
Henk Wymeersch
Abstract:
State-of-the-art high-spectral-efficiency communication systems employ high-order modulation formats coupled with high symbol rates to accommodate the ever-growing demand for data rate-hungry applications. However, such systems are more vulnerable to linear and nonlinear transmission impairments, and it is important to mitigate the performance loss via digital signal processing. In this paper, we…
▽ More
State-of-the-art high-spectral-efficiency communication systems employ high-order modulation formats coupled with high symbol rates to accommodate the ever-growing demand for data rate-hungry applications. However, such systems are more vulnerable to linear and nonlinear transmission impairments, and it is important to mitigate the performance loss via digital signal processing. In this paper, we propose a novel machine learning approach for blind channel equalization and estimation using the vector quantized (VQ) \ac{VAE} framework. The proposed approach generalizes the applicability of the conventional \ac{VAE}-based equalizer to nonlinear systems employing high-order modulation formats by introducing a codebook component and an associated novel loss function. We evaluate the performance of the proposed method over a linear additive white Gaussian noise channel with intersymbol interference and two nonlinear scenarios. Simulation results show that the proposed method can achieve similar performance as a data aided equalizer using the \acf{MMSE} criterion, and outperforms the blind\ac{CMA} and the \ac{VAE}-based channel equalizer. Furthermore, we show that for the linear channel, the proposed scheme exhibits better convergence properties than the \ac{MMSE}-based, the \ac{CMA}-based, and the \ac{VAE}-based equalizers in terms of both convergence speed and robustness to variations in training batch size and learning rate.
△ Less
Submitted 22 February, 2023;
originally announced February 2023.
-
End-to-End Learning for VCSEL-based Optical Interconnects: State-of-the-Art, Challenges, and Opportunities
Authors:
Muralikrishnan Srinivasan,
Jinxiang Song,
Alexander Grabowski,
Krzysztof Szczerba,
Holger K. Iversen,
Mikkel N. Schmidt,
Darko Zibar,
Jochen Schröder,
Anders Larsson,
Christian Häger,
Henk Wymeersch
Abstract:
Optical interconnects (OIs) based on vertical-cavity surface-emitting lasers (VCSELs) are the main workhorse within data centers, supercomputers, and even vehicles, providing low-cost, high-rate connectivity. VCSELs must operate under extremely harsh and time-varying conditions, thus requiring adaptive and flexible designs of the communication chain. Such designs can be built based on mathematical…
▽ More
Optical interconnects (OIs) based on vertical-cavity surface-emitting lasers (VCSELs) are the main workhorse within data centers, supercomputers, and even vehicles, providing low-cost, high-rate connectivity. VCSELs must operate under extremely harsh and time-varying conditions, thus requiring adaptive and flexible designs of the communication chain. Such designs can be built based on mathematical models (model-based design) or learned from data (machine learning (ML) based design). Various ML techniques have recently come to the forefront, replacing individual components in the transmitters and receivers with deep neural networks. Beyond such component-wise learning, end-to-end (E2E) autoencoder approaches can reach the ultimate performance through co-optimizing entire parameterized transmitters and receivers. This tutorial paper aims to provide an overview of ML for VCSEL-based OIs, with a focus on E2E approaches, dealing specifically with the unique challenges facing VCSELs, such as the wide temperature variations and complex models.
△ Less
Submitted 25 November, 2022;
originally announced November 2022.
-
Size-effects of metamaterial beams subjected to pure bending: on boundary conditions and parameter identification in the relaxed micromorphic model
Authors:
Mohammad Sarhil,
Lisa Scheunemann,
Jörg Schröder,
Patrizio Neff
Abstract:
In this paper we model the size-effects of metamaterial beams under bending with the aid of the relaxed micromorphic continuum. We analyze first the size-dependent bending stiffness of heterogeneous fully discretized metamaterial beams subjected to pure bending loads. Two equivalent loading schemes are introduced which lead to a constant moment along the beam length with no shear force. The relaxe…
▽ More
In this paper we model the size-effects of metamaterial beams under bending with the aid of the relaxed micromorphic continuum. We analyze first the size-dependent bending stiffness of heterogeneous fully discretized metamaterial beams subjected to pure bending loads. Two equivalent loading schemes are introduced which lead to a constant moment along the beam length with no shear force. The relaxed micromorphic model is employed then to retrieve the size-effects. We present a procedure for the determination of the material parameters of the relaxed micromorphic model based on the fact that the model operates between two well-defined scales. These scales are given by linear elasticity with micro and macro elasticity tensors which bound the relaxed micromorphic continuum from above and below, respectively. The micro elasticity tensor is specified as the maximum possible stiffness that is exhibited by the assumed metamaterial while the macro elasticity tensor is given by standard periodic first-order homogenization. For the identification of the micro elasticity tensor, two different approaches are shown which rely on affine and non-affine Dirichlet boundary conditions of candidate unit cell variants with the possible stiffest response. The consistent coupling condition is shown to allow the model to act on the whole intended range between macro and micro elasticity tensors for both loading cases. We fit the relaxed micromorphic model against the fully resolved metamaterial solution by controlling the curvature magnitude after linking it with the specimen's size. The obtained parameters of the relaxed micromorphic model are tested for two additional loading scenarios.
△ Less
Submitted 28 March, 2023; v1 submitted 31 October, 2022;
originally announced October 2022.
-
Four-periodic infinite staircases for four-dimensional polydisks
Authors:
Caden Farley,
Tara Holm,
Nicki Magill,
Jemma Schroder,
Morgan Weiler,
Zichen Wang,
Elizaveta Zabelina
Abstract:
The ellipsoid embedding function of a symplectic four-manifold measures the amount by which its symplectic form must be scaled in order for it to admit an embedding of an ellipsoid of varying eccentricity. This function generalizes the Gromov width and ball packing numbers. In the one continuous family of symplectic four-manifolds that has been analyzed, one-point blowups of the complex projective…
▽ More
The ellipsoid embedding function of a symplectic four-manifold measures the amount by which its symplectic form must be scaled in order for it to admit an embedding of an ellipsoid of varying eccentricity. This function generalizes the Gromov width and ball packing numbers. In the one continuous family of symplectic four-manifolds that has been analyzed, one-point blowups of the complex projective plane, there is an open dense set of symplectic forms whose ellipsoid embedding functions are completely described by finitely many obstructions, while there is simultaneously a Cantor set of symplectic forms for which an infinite number of obstructions are needed. In the latter case, we say that the embedding function has an infinite staircase. In this paper we identify a new infinite staircase when the target is a four-dimensional polydisk, extending a countable family identified by Usher in 2019. Our work computes the function on infinitely many intervals and thereby indicates a method of proof for a conjecture of Usher.
△ Less
Submitted 28 August, 2023; v1 submitted 26 October, 2022;
originally announced October 2022.
-
LungViT: Ensembling Cascade of Texture Sensitive Hierarchical Vision Transformers for Cross-Volume Chest CT Image-to-Image Translation
Authors:
Muhammad F. A. Chaudhary,
Sarah E. Gerard,
Gary E. Christensen,
Christopher B. Cooper,
Joyce D. Schroeder,
Eric A. Hoffman,
Joseph M. Reinhardt
Abstract:
Chest computed tomography (CT) at inspiration is often complemented by an expiratory CT to identify peripheral airways disease. Additionally, co-registered inspiratory-expiratory volumes can be used to derive various markers of lung function. Expiratory CT scans, however, may not be acquired due to dose or scan time considerations or may be inadequate due to motion or insufficient exhale; leading…
▽ More
Chest computed tomography (CT) at inspiration is often complemented by an expiratory CT to identify peripheral airways disease. Additionally, co-registered inspiratory-expiratory volumes can be used to derive various markers of lung function. Expiratory CT scans, however, may not be acquired due to dose or scan time considerations or may be inadequate due to motion or insufficient exhale; leading to a missed opportunity to evaluate underlying small airways disease. Here, we propose LungViT - a generative adversarial learning approach using hierarchical vision transformers for translating inspiratory CT intensities to corresponding expiratory CT intensities. LungViT addresses several limitations of the traditional generative models including slicewise discontinuities, limited size of generated volumes, and their inability to model texture transfer at volumetric level. We propose a shifted-window hierarchical vision transformer architecture with squeeze-and-excitation decoder blocks for modeling dependencies between features. We also propose a multiview texture similarity distance metric for texture and style transfer in 3D. To incorporate global information into the training process and refine the output of our model, we use ensemble cascading. LungViT is able to generate large 3D volumes of size 320 x 320 x 320. We train and validate our model using a diverse cohort of 1500 subjects with varying disease severity. To assess model generalizability beyond the development set biases, we evaluate our model on an out-of-distribution external validation set of 200 subjects. Clinical validation on internal and external testing sets shows that synthetic volumes could be reliably adopted for deriving clinical endpoints of chronic obstructive pulmonary disease.
△ Less
Submitted 27 August, 2023; v1 submitted 5 October, 2022;
originally announced October 2022.
-
2D reconstruction of magnetotail electron diffusion region measured by MMS
Authors:
J. M. Schroeder,
J. Egedal,
G. Cozzani,
Yu. V. Khotyaintsev,
W. Daughton,
R. E. Denton,
J. L. Burch
Abstract:
Models for collisionless magnetic reconnection in near-Earth space are distinctly characterized as 2D or 3D. In 2D kinetic models, the frozen-in law for the electron fluid is usually broken by laminar dynamics involving structures set by the electron orbit size, while in 3D models the width of the electron diffusion region is broadened by turbulent effects. We present an analysis of in situ spacec…
▽ More
Models for collisionless magnetic reconnection in near-Earth space are distinctly characterized as 2D or 3D. In 2D kinetic models, the frozen-in law for the electron fluid is usually broken by laminar dynamics involving structures set by the electron orbit size, while in 3D models the width of the electron diffusion region is broadened by turbulent effects. We present an analysis of in situ spacecraft observations from the Earth's magnetotail of a fortuitous encounter with an active reconnection region, mapping the observations onto a 2D spatial domain. While the event likely was perturbed by low-frequency 3D dynamics, the structure of the electron diffusion region remains consistent with results from a 2D kinetic simulation. As such, the event represents a unique validation of 2D kinetic, and laminar reconnection models.
△ Less
Submitted 24 September, 2022;
originally announced September 2022.
-
Efficient multigrid reduction-in-time for method-of-lines discretizations of linear advection
Authors:
H. De Sterck,
R. D. Falgout,
O. A. Krzysik,
J. B. Schroder
Abstract:
Parallel-in-time methods for partial differential equations (PDEs) have been the subject of intense development over recent decades, particularly for diffusion-dominated problems. It has been widely reported in the literature, however, that many of these methods perform quite poorly for advection-dominated problems. Here we analyze the particular iterative parallel-in-time algorithm of multigrid r…
▽ More
Parallel-in-time methods for partial differential equations (PDEs) have been the subject of intense development over recent decades, particularly for diffusion-dominated problems. It has been widely reported in the literature, however, that many of these methods perform quite poorly for advection-dominated problems. Here we analyze the particular iterative parallel-in-time algorithm of multigrid reduction-in-time (MGRIT) for discretizations of constant-wave-speed linear advection problems. We focus on common method-of-lines discretizations that employ upwind finite differences in space and Runge-Kutta methods in time. Using a convergence framework we developed in previous work, we prove for a subclass of these discretizations that, if using the standard approach of rediscretizing the fine-grid problem on the coarse grid, robust MGRIT convergence with respect to CFL number and coarsening factor is not possible. This poor convergence and non-robustness is caused, at least in part, by an inadequate coarse-grid correction for smooth Fourier modes known as characteristic components.We propose an alternative coarse-grid that provides a better correction of these modes. This coarse-grid operator is related to previous work and uses a semi-Lagrangian discretization combined with an implicitly treated truncation error correction. Theory and numerical experiments show the coarse-grid operator yields fast MGRIT convergence for many of the method-of-lines discretizations considered, including for both implicit and explicit discretizations of high order. Parallel results demonstrate substantial speed-up over sequential time-stepping.
△ Less
Submitted 20 March, 2023; v1 submitted 14 September, 2022;
originally announced September 2022.
-
Multigrid Reduction in Time for Chaotic Dynamical Systems
Authors:
David A. Vargas,
Robert D. Falgout,
Stefanie Günther,
Jacob B. Schroder
Abstract:
As CPU clock speeds have stagnated and high performance computers continue to have ever higher core counts, increased parallelism is needed to take advantage of these new architectures. Traditional serial time-marching schemes can be a significant bottleneck, as many types of simulations require large numbers of time-steps which must be computed sequentially. Parallel in Time schemes, such as the…
▽ More
As CPU clock speeds have stagnated and high performance computers continue to have ever higher core counts, increased parallelism is needed to take advantage of these new architectures. Traditional serial time-marching schemes can be a significant bottleneck, as many types of simulations require large numbers of time-steps which must be computed sequentially. Parallel in Time schemes, such as the Multigrid Reduction in Time (MGRIT) method, remedy this by parallelizing across time-steps, and have shown promising results for parabolic problems. However, chaotic problems have proved more difficult, since chaotic initial value problems (IVPs) are inherently ill-conditioned. MGRIT relies on a hierarchy of successively coarser time-grids to iteratively correct the solution on the finest time-grid, but due to the nature of chaotic systems, small inaccuracies on the coarser levels can be greatly magnified and lead to poor coarse-grid corrections. Here we introduce a modified MGRIT algorithm based on an existing quadratically converging nonlinear extension to the multigrid Full Approximation Scheme (FAS), as well as a novel time-coarsening scheme. Together, these approaches better capture long-term chaotic behavior on coarse-grids and greatly improve convergence of MGRIT for chaotic IVPs. Further, we introduce a novel low memory variant of the algorithm for solving chaotic PDEs with MGRIT which not only solves the IVP, but also provides estimates for the unstable Lyapunov vectors of the system. We provide supporting numerical results for the Lorenz system and demonstrate parallel speedup for the chaotic Kuramoto- Sivashinsky partial differential equation over a significantly longer time-domain than in previous works.
△ Less
Submitted 26 August, 2022;
originally announced August 2022.
-
Parallel Energy-Minimization Prolongation for Algebraic Multigrid
Authors:
Carlo Janna,
Andrea Franceschini,
Jacob B. Schroder,
Luke Olson
Abstract:
Algebraic multigrid (AMG) is one of the most widely used solution techniques for linear systems of equations arising from discretized partial differential equations. The popularity of AMG stems from its potential to solve linear systems in almost linear time, that is with an O(n) complexity, where n is the problem size. This capability is crucial at the present, where the increasing availability o…
▽ More
Algebraic multigrid (AMG) is one of the most widely used solution techniques for linear systems of equations arising from discretized partial differential equations. The popularity of AMG stems from its potential to solve linear systems in almost linear time, that is with an O(n) complexity, where n is the problem size. This capability is crucial at the present, where the increasing availability of massive HPC platforms pushes for the solution of very large problems. The key for a rapidly converging AMG method is a good interplay between the smoother and the coarse-grid correction, which in turn requires the use of an effective prolongation. From a theoretical viewpoint, the prolongation must accurately represent near kernel components and, at the same time, be bounded in the energy norm. For challenging problems, however, ensuring both these requirements is not easy and is exactly the goal of this work. We propose a constrained minimization procedure aimed at reducing prolongation energy while preserving the near kernel components in the span of interpolation. The proposed algorithm is based on previous energy minimization approaches utilizing a preconditioned restricted conjugate gradients method, but has new features and a specific focus on parallel performance and implementation. It is shown that the resulting solver, when used for large real-world problems from various application fields, exhibits excellent convergence rates and scalability and outperforms at least some more traditional AMG approaches.
△ Less
Submitted 5 August, 2022;
originally announced August 2022.
-
Localization supervision of chest x-ray classifiers using label-specific eye-tracking annotation
Authors:
Ricardo Bigolin Lanfredi,
Joyce D. Schroeder,
Tolga Tasdizen
Abstract:
Convolutional neural networks (CNNs) have been successfully applied to chest x-ray (CXR) images. Moreover, annotated bounding boxes have been shown to improve the interpretability of a CNN in terms of localizing abnormalities. However, only a few relatively small CXR datasets containing bounding boxes are available, and collecting them is very costly. Opportunely, eye-tracking (ET) data can be col…
▽ More
Convolutional neural networks (CNNs) have been successfully applied to chest x-ray (CXR) images. Moreover, annotated bounding boxes have been shown to improve the interpretability of a CNN in terms of localizing abnormalities. However, only a few relatively small CXR datasets containing bounding boxes are available, and collecting them is very costly. Opportunely, eye-tracking (ET) data can be collected in a non-intrusive way during the clinical workflow of a radiologist. We use ET data recorded from radiologists while dictating CXR reports to train CNNs. We extract snippets from the ET data by associating them with the dictation of keywords and use them to supervise the localization of specific abnormalities. We show that this method improves a model's interpretability without impacting its image-level classification.
△ Less
Submitted 14 December, 2022; v1 submitted 20 July, 2022;
originally announced July 2022.
-
Benchmarking and Interpreting End-to-end Learning of MIMO and Multi-User Communication
Authors:
Jinxiang Song,
Christian Häger,
Jochen Schröder,
Timothy J. O'Shea,
Erik Agrell,
Henk Wymeersch
Abstract:
End-to-end autoencoder (AE) learning has the potential of exceeding the performance of human-engineered transceivers and encoding schemes, without a priori knowledge of communication-theoretic principles. In this work, we aim to understand to what extent and for which scenarios this claim holds true when comparing with fair benchmarks. Our particular focus is on memoryless multiple-input multiple-…
▽ More
End-to-end autoencoder (AE) learning has the potential of exceeding the performance of human-engineered transceivers and encoding schemes, without a priori knowledge of communication-theoretic principles. In this work, we aim to understand to what extent and for which scenarios this claim holds true when comparing with fair benchmarks. Our particular focus is on memoryless multiple-input multiple-output (MIMO) and multi-user (MU) systems. Four case studies are considered: two point-to-point (closed-loop and open-loop MIMO) and two MU scenarios (MIMO broadcast and interference channels). For the point-to-point scenarios, we explain some of the performance gains observed in prior work through the selection of improved baseline schemes that include geometric shaping as well as bit and power allocation. For the MIMO broadcast channel, we demonstrate the feasibility of a novel AE method with centralized learning and decentralized execution. Interestingly, the learned scheme performs close to nonlinear vector-perturbation precoding and significantly outperforms conventional zero-forcing. Lastly, we highlight potential pitfalls when interpreting learned communication schemes. In particular, we show that the AE for the considered interference channel learns to avoid interference, albeit in a rotated reference frame. After de-rotating the learned signal constellation of each user, the resulting scheme corresponds to conventional time sharing with geometric shaping.
△ Less
Submitted 15 March, 2022;
originally announced March 2022.
-
Power-efficient soliton microcombs
Authors:
Óskar B. Helgason,
Marcello Girardi,
Zhichao Ye,
Fuchuan Lei,
Jochen Schröder,
Victor Torres Company
Abstract:
Laser frequency combs are enabling some of the most exciting scientific endeavours in the 21st century, ranging from the development of optical clocks to the calibration of the astronomical spectrographs used for searching Earth-like exoplanets. Today, dissipative Kerr solitons generated in microresonators offer the prospect of attaining frequency combs in miniaturized systems by capitalizing on a…
▽ More
Laser frequency combs are enabling some of the most exciting scientific endeavours in the 21st century, ranging from the development of optical clocks to the calibration of the astronomical spectrographs used for searching Earth-like exoplanets. Today, dissipative Kerr solitons generated in microresonators offer the prospect of attaining frequency combs in miniaturized systems by capitalizing on advances in photonic integration. Most of the applications based on soliton microcombs rely on tuning a continuous-wave laser into a longitudinal mode of a microresonator whose dimensions are engineered to display anomalous dispersion at the pump laser frequency. In this configuration, however, nonlinear physics precludes from attaining dissipative Kerr solitons with high power conversion efficiency, with typical comb powers amounting to ~1% of the available laser power. Here, we demonstrate that this fundamental limitation can be overcome by inducing a controllable frequency shift to a selected cavity resonance. Experimentally, we realize this shift using two linearly coupled anomalous-dispersion microresonators (a photonic molecule), resulting in a coherent dissipative Kerr soliton with a conversion efficiency exceeding 50% and excellent line spacing stability. We describe the physical soliton dynamics in this configuration, and discover the system displays unusual characteristics, such as the possibility to backwards initiate solitons and stable operation with a blue detuned pump laser. By optimizing the microcomb power available on chip, these results facilitate the practical implementation of a scalable integrated photonic architecture for energy-efficient applications.
△ Less
Submitted 18 February, 2022;
originally announced February 2022.
-
Toward Parallel in Time for Chaotic Dynamical Systems
Authors:
David A. Vargas,
Robert D. Falgout,
Stefanie Günther,
Jacob B. Schroder
Abstract:
As CPU clock speeds have stagnated, and high performance computers continue to have ever higher core counts, increased parallelism is needed to take advantage of these new architectures. Traditional serial time-marching schemes are a significant bottleneck, as many types of simulations require large numbers of time-steps which must be computed sequentially. Parallel in Time schemes, such as the Mu…
▽ More
As CPU clock speeds have stagnated, and high performance computers continue to have ever higher core counts, increased parallelism is needed to take advantage of these new architectures. Traditional serial time-marching schemes are a significant bottleneck, as many types of simulations require large numbers of time-steps which must be computed sequentially. Parallel in Time schemes, such as the Multigrid Reduction in Time (MGRIT) method, remedy this by parallelizing across time-steps, and have shown promising results for parabolic problems. However, chaotic problems have proved more difficult, since chaotic initial value problems are inherently ill-conditioned. MGRIT relies on a hierarchy of successively coarser time-grids to iteratively correct the solution on the finest time-grid, but due to the nature of chaotic systems, subtle inaccuracies on the coarser levels can lead to poor coarse-grid corrections. Here we propose a modification to nonlinear FAS multigrid, as well as a novel time-coarsening scheme, which together better capture long term behavior on coarse grids and greatly improve convergence of MGRIT for chaotic initial value problems. We provide supporting numerical results for the Lorenz system model problem.
△ Less
Submitted 25 January, 2022;
originally announced January 2022.
-
Comparing radiologists' gaze and saliency maps generated by interpretability methods for chest x-rays
Authors:
Ricardo Bigolin Lanfredi,
Ambuj Arora,
Trafton Drew,
Joyce D. Schroeder,
Tolga Tasdizen
Abstract:
The interpretability of medical image analysis models is considered a key research field. We use a dataset of eye-tracking data from five radiologists to compare the outputs of interpretability methods and the heatmaps representing where radiologists looked. We conduct a class-independent analysis of the saliency maps generated by two methods selected from the literature: Grad-CAM and attention ma…
▽ More
The interpretability of medical image analysis models is considered a key research field. We use a dataset of eye-tracking data from five radiologists to compare the outputs of interpretability methods and the heatmaps representing where radiologists looked. We conduct a class-independent analysis of the saliency maps generated by two methods selected from the literature: Grad-CAM and attention maps from an attention-gated model. For the comparison, we use shuffled metrics, which avoid biases from fixation locations. We achieve scores comparable to an interobserver baseline in one shuffled metric, highlighting the potential of saliency maps from Grad-CAM to mimic a radiologist's attention over an image. We also divide the dataset into subsets to evaluate in which cases similarities are higher.
△ Less
Submitted 19 April, 2023; v1 submitted 22 December, 2021;
originally announced December 2021.
-
Lagrange and $H(\operatorname{curl},{\cal B})$ based Finite Element formulations for the relaxed micromorphic model
Authors:
Jörg Schröder,
Mohammad Sarhil,
Lisa Scheunemann,
Patrizio Neff
Abstract:
Modeling the unusual mechanical properties of metamaterials is a challenging topic for the mechanics community and enriched continuum theories are promising computational tools for such materials. The so-called relaxed micromorphic model has shown many advantages in this field. In this contribution, we present the significant aspects related to the relaxed micromorphic model realization with the f…
▽ More
Modeling the unusual mechanical properties of metamaterials is a challenging topic for the mechanics community and enriched continuum theories are promising computational tools for such materials. The so-called relaxed micromorphic model has shown many advantages in this field. In this contribution, we present the significant aspects related to the relaxed micromorphic model realization with the finite element method. The variational problem is derived and different FEM-formulations for the two-dimensional case are presented. These are a nodal standard formulation $H^1({\cal B}) \times H^1({\cal B})$ and a nodal-edge formulation $H^1({\cal B}) \times H(\operatorname{curl}, {\cal B})$, where the latter employs the Nédélec space. However, the implementation of higher-order Nédélec elements is not trivial and requires some technicalities which are demonstrated. We discuss the convergence behavior of Lagrange-type and tangential-conforming finite element discretizations. Moreover, we analyze the characteristic length effect on the different components of the model and reveal how the size-effect property is captured via this characteristic length.
△ Less
Submitted 1 December, 2021;
originally announced December 2021.
-
Model-Based End-to-End Learning for WDM Systems With Transceiver Hardware Impairments
Authors:
Jinxiang Song,
Christian Häger,
Jochen Schröder,
Alexandre Graell i Amat,
Henk Wymeersch
Abstract:
We propose an AE-based transceiver for a WDM system impaired by hardware imperfections. We design our AE following the architecture of conventional communication systems. This enables to initialize the AE-based transceiver to have similar performance to its conventional counterpart prior to training and improves the training convergence rate. We first train the AE in a single-channel system, and s…
▽ More
We propose an AE-based transceiver for a WDM system impaired by hardware imperfections. We design our AE following the architecture of conventional communication systems. This enables to initialize the AE-based transceiver to have similar performance to its conventional counterpart prior to training and improves the training convergence rate. We first train the AE in a single-channel system, and show that it achieves performance improvements by putting energy outside the desired bandwidth, and therefore cannot be used for a WDM system. We then train the AE in a WDM setup. Simulation results show that the proposed AE significantly outperforms the conventional approach. More specifically, it increases the spectral efficiency of the considered system by reducing the guard band by 37\% and 50\% for a root-raised-cosine filter-based matched filter with 10\% and 1\% roll-off, respectively. An ablation study indicates that the performance gain can be ascribed to the optimization of the symbol mapper, the pulse-shaping filter, and the symbol demapper. Finally, we use reinforcement learning to learn the pulse-shaping filter under the assumption that the channel model is unknown. Simulation results show that the reinforcement-learning-based algorithm achieves similar performance to the standard supervised end-to-end learning approach assuming perfect channel knowledge.
△ Less
Submitted 29 November, 2021;
originally announced November 2021.
-
On the 3-colorability of triangle-free and fork-free graphs
Authors:
Joshua Schroeder,
Zhiyu Wang,
Xingxing Yu
Abstract:
A graph $G$ is said to satisfy the Vizing bound if $χ(G)\leq ω(G)+1$, where $χ(G)$ and $ω(G)$ denote the chromatic number and clique number of $G$, respectively. It was conjectured by Randerath in 1998 that if $G$ is a triangle-free and fork-free graph, where the fork (also known as trident) is obtained from $K_{1,4}$ by subdividing two edges, then $G$ satisfies the Vizing bound. In this paper, we…
▽ More
A graph $G$ is said to satisfy the Vizing bound if $χ(G)\leq ω(G)+1$, where $χ(G)$ and $ω(G)$ denote the chromatic number and clique number of $G$, respectively. It was conjectured by Randerath in 1998 that if $G$ is a triangle-free and fork-free graph, where the fork (also known as trident) is obtained from $K_{1,4}$ by subdividing two edges, then $G$ satisfies the Vizing bound. In this paper, we confirm this conjecture.
△ Less
Submitted 19 November, 2021;
originally announced November 2021.
-
Single volume lung biomechanics from chest computed tomography using a mode preserving generative adversarial network
Authors:
Muhammad F. A. Chaudhary,
Sarah E. Gerard,
Di Wang,
Gary E. Christensen,
Christopher B. Cooper,
Joyce D. Schroeder,
Eric A. Hoffman,
Joseph M. Reinhardt
Abstract:
Local tissue expansion of the lungs is typically derived by registering computed tomography (CT) scans acquired at multiple lung volumes. However, acquiring multiple scans incurs increased radiation dose, time, and cost, and may not be possible in many cases, thus restricting the applicability of registration-based biomechanics. We propose a generative adversarial learning approach for estimating…
▽ More
Local tissue expansion of the lungs is typically derived by registering computed tomography (CT) scans acquired at multiple lung volumes. However, acquiring multiple scans incurs increased radiation dose, time, and cost, and may not be possible in many cases, thus restricting the applicability of registration-based biomechanics. We propose a generative adversarial learning approach for estimating local tissue expansion directly from a single CT scan. The proposed framework was trained and evaluated on 2500 subjects from the SPIROMICS cohort. Once trained, the framework can be used as a registration-free method for predicting local tissue expansion. We evaluated model performance across varying degrees of disease severity and compared its performance with two image-to-image translation frameworks - UNet and Pix2Pix. Our model achieved an overall PSNR of 18.95 decibels, SSIM of 0.840, and Spearman's correlation of 0.61 at a high spatial resolution of 1 mm3.
△ Less
Submitted 15 October, 2021;
originally announced October 2021.
-
REFLACX, a dataset of reports and eye-tracking data for localization of abnormalities in chest x-rays
Authors:
Ricardo Bigolin Lanfredi,
Mingyuan Zhang,
William F. Auffermann,
Jessica Chan,
Phuong-Anh T. Duong,
Vivek Srikumar,
Trafton Drew,
Joyce D. Schroeder,
Tolga Tasdizen
Abstract:
Deep learning has shown recent success in classifying anomalies in chest x-rays, but datasets are still small compared to natural image datasets. Supervision of abnormality localization has been shown to improve trained models, partially compensating for dataset sizes. However, explicitly labeling these anomalies requires an expert and is very time-consuming. We propose a potentially scalable meth…
▽ More
Deep learning has shown recent success in classifying anomalies in chest x-rays, but datasets are still small compared to natural image datasets. Supervision of abnormality localization has been shown to improve trained models, partially compensating for dataset sizes. However, explicitly labeling these anomalies requires an expert and is very time-consuming. We propose a potentially scalable method for collecting implicit localization data using an eye tracker to capture gaze locations and a microphone to capture a dictation of a report, imitating the setup of a reading room. The resulting REFLACX (Reports and Eye-Tracking Data for Localization of Abnormalities in Chest X-rays) dataset was labeled across five radiologists and contains 3,032 synchronized sets of eye-tracking data and timestamped report transcriptions for 2,616 chest x-rays from the MIMIC-CXR dataset. We also provide auxiliary annotations, including bounding boxes around lungs and heart and validation labels consisting of ellipses localizing abnormalities and image-level labels. Furthermore, a small subset of the data contains readings from all radiologists, allowing for the calculation of inter-rater scores.
△ Less
Submitted 28 June, 2022; v1 submitted 29 September, 2021;
originally announced September 2021.
-
Stability of superthermal strahl electrons in the solar wind
Authors:
Jack M. Schroeder,
Stanislav Boldyrev,
Patrick Astfalk
Abstract:
We present a kinetic stability analysis of the solar wind electron distribution function consisting of the Maxwellian core and the magnetic-field aligned strahl, a superthermal electron beam propagating away from the sun. We use an electron strahl distribution function obtained as a solution of a weakly collisional drift-kinetic equation, representative of a strahl affected by Coulomb collisions b…
▽ More
We present a kinetic stability analysis of the solar wind electron distribution function consisting of the Maxwellian core and the magnetic-field aligned strahl, a superthermal electron beam propagating away from the sun. We use an electron strahl distribution function obtained as a solution of a weakly collisional drift-kinetic equation, representative of a strahl affected by Coulomb collisions but unadulterated by possible broadening from turbulence. This distribution function is essentially non-Maxwellian and varies with the heliospheric distance. The stability analysis is performed with the Vlasov-Maxwell linear solver LEOPARD. We find that depending on the heliospheric distance, the core-strahl electron distribution becomes unstable with respect to sunward-propagating kinetic-Alfvén, magnetosonic, and whistler modes, in a broad range of propagation angles. The wavenumbers of the unstable modes are close to the ion inertial scales, and the radial distances at which the instabilities first appear are on the order of 1 AU. However, we have not detected any instabilities driven by resonant wave interactions with the superthermal strahl electrons. Instead, the observed instabilities are triggered by a relative drift between the electron and ion cores necessary to maintain zero electric current in the solar wind frame (ion frame). Contrary to strahl distributions modeled by shifted Maxwellians, the electron strahl obtained as a solution of the kinetic equation is stable. Our results are consistent with the previous studies based on a more restricted solution for the electron strahl.
△ Less
Submitted 13 October, 2021; v1 submitted 16 July, 2021;
originally announced July 2021.
-
Score refinement for confidence-based 3D multi-object tracking
Authors:
Nuri Benbarka,
Jona Schröder,
Andreas Zell
Abstract:
Multi-object tracking is a critical component in autonomous navigation, as it provides valuable information for decision-making. Many researchers tackled the 3D multi-object tracking task by filtering out the frame-by-frame 3D detections; however, their focus was mainly on finding useful features or proper matching metrics. Our work focuses on a neglected part of the tracking system: score refinem…
▽ More
Multi-object tracking is a critical component in autonomous navigation, as it provides valuable information for decision-making. Many researchers tackled the 3D multi-object tracking task by filtering out the frame-by-frame 3D detections; however, their focus was mainly on finding useful features or proper matching metrics. Our work focuses on a neglected part of the tracking system: score refinement and tracklet termination. We show that manipulating the scores depending on time consistency while terminating the tracklets depending on the tracklet score improves tracking results. We do this by increasing the matched tracklets' score with score update functions and decreasing the unmatched tracklets' score. Compared to count-based methods, our method consistently produces better AMOTA and MOTA scores when utilizing various detectors and filtering algorithms on different datasets. The improvements in AMOTA score went up to 1.83 and 2.96 in MOTA. We also used our method as a late-fusion ensembling method, and it performed better than voting-based ensemble methods by a solid margin. It achieved an AMOTA score of 67.6 on nuScenes test evaluation, which is comparable to other state-of-the-art trackers. Code is publicly available at: \url{https://github.com/cogsys-tuebingen/CBMOT}.
△ Less
Submitted 9 July, 2021;
originally announced July 2021.
-
Weighted Relaxation for Multigrid Reduction in Time
Authors:
Masumi Sugiyama,
Jacob B. Schroder,
Ben S. Southworth,
Stephanie Friedhoff
Abstract:
Based on current trends in computer architectures, faster compute speeds must come from increased parallelism rather than increased clock speeds, which are currently stagnate. This situation has created the well-known bottleneck for sequential time-integration, where each individual time-value (i.e., time-step) is computed sequentially. One approach to alleviate this and achieve parallelism in tim…
▽ More
Based on current trends in computer architectures, faster compute speeds must come from increased parallelism rather than increased clock speeds, which are currently stagnate. This situation has created the well-known bottleneck for sequential time-integration, where each individual time-value (i.e., time-step) is computed sequentially. One approach to alleviate this and achieve parallelism in time is with multigrid. In this work, we consider multigrid-reduction-in-time (MGRIT), a multilevel method applied to the time dimension that computes multiple time-steps in parallel. Like all multigrid methods, MGRIT relies on the complementary relationship between relaxation on a fine-grid and a correction from the coarse grid to solve the problem. All current MGRIT implementations are based on unweighted-Jacobi relaxation; here we introduce the concept of weighted relaxation to MGRIT. We derive new convergence bounds for weighted relaxation, and use this analysis to guide the selection of relaxation weights. Numerical results then demonstrate that non-unitary relaxation weights consistently yield faster convergence rates and lower iteration counts for MGRIT when compared with unweighted relaxation. In most cases, weighted relaxation yields a 10%-20% saving in iterations. For A-stable integration schemes, results also illustrate that under-relaxation can restore convergence in some cases where unweighted relaxation is not convergent.
△ Less
Submitted 5 July, 2021;
originally announced July 2021.
-
Polynomial $χ$-binding functions for $t$-broom-free graphs
Authors:
Xiaonan Liu,
Joshua Schroeder,
Zhiyu Wang,
Xingxing Yu
Abstract:
For any positive integer $t$, a \emph{$t$-broom} is a graph obtained from $K_{1,t+1}$ by subdividing an edge once. In this paper, we show that, for graphs $G$ without induced $t$-brooms, we have $χ(G) = o(ω(G)^{t+1})$, where $χ(G)$ and $ω(G)$ are the chromatic number and clique number of $G$, respectively. When $t=2$, this answers a question of Schiermeyer and Randerath. Moreover, for $t=2$, we st…
▽ More
For any positive integer $t$, a \emph{$t$-broom} is a graph obtained from $K_{1,t+1}$ by subdividing an edge once. In this paper, we show that, for graphs $G$ without induced $t$-brooms, we have $χ(G) = o(ω(G)^{t+1})$, where $χ(G)$ and $ω(G)$ are the chromatic number and clique number of $G$, respectively. When $t=2$, this answers a question of Schiermeyer and Randerath. Moreover, for $t=2$, we strengthen the bound on $χ(G)$ to $7ω(G)^2$, confirming a conjecture of Sivaraman. For $t\geq 3$ and \{$t$-broom, $K_{t,t}$\}-free graphs, we improve the bound to $o(ω^{t})$.
△ Less
Submitted 29 November, 2022; v1 submitted 16 June, 2021;
originally announced June 2021.
-
Over-the-fiber Digital Predistortion Using Reinforcement Learning
Authors:
Jinxiang Song,
Zonglong He,
Christian Häger,
Magnus Karlsson,
Alexandre Graell i Amat,
Henk Wymeersch,
Jochen Schröder
Abstract:
We demonstrate, for the first time, experimental over-the-fiber training of transmitter neural networks (NNs) using reinforcement learning. Optical back-to-back training of a novel NN-based digital predistorter outperforms arcsine-based predistortion with up to 60\% bit-error-rate reduction.
We demonstrate, for the first time, experimental over-the-fiber training of transmitter neural networks (NNs) using reinforcement learning. Optical back-to-back training of a novel NN-based digital predistorter outperforms arcsine-based predistortion with up to 60\% bit-error-rate reduction.
△ Less
Submitted 9 June, 2021;
originally announced June 2021.
-
Time-periodic steady-state solution of fluid-structure interaction and cardiac flow problems through multigrid-reduction-in-time
Authors:
Andreas Hessenthaler,
Robert D. Falgout,
Jacob B. Schroder,
Adelaide de Vecchi,
David Nordsletten,
Oliver Röhrle
Abstract:
In this paper, a time-periodic MGRIT algorithm is proposed as a means to reduce the time-to-solution of numerical algorithms by exploiting the time periodicity inherent to many applications in science and engineering. The time-periodic MGRIT algorithm is applied to a variety of linear and nonlinear single- and multiphysics problems that are periodic-in-time. It is demonstrated that the proposed pa…
▽ More
In this paper, a time-periodic MGRIT algorithm is proposed as a means to reduce the time-to-solution of numerical algorithms by exploiting the time periodicity inherent to many applications in science and engineering. The time-periodic MGRIT algorithm is applied to a variety of linear and nonlinear single- and multiphysics problems that are periodic-in-time. It is demonstrated that the proposed parallel-in-time algorithm can obtain the same time-periodic steady-state solution as sequential time-stepping. It is shown that the required number of MGRIT iterations can be estimated a priori and that the new MGRIT variant can significantly and consistently reduce the time-to-solution compared to sequential time-stepping, irrespective of the number of dimensions, linear or nonlinear PDE models, single-physics or coupled problems and the employed computing resources. The numerical experiments demonstrate that the time-periodic MGRIT algorithm enables a greater level of parallelism yielding faster turnaround, and thus, facilitating more complex and more realistic problems to be solved.
△ Less
Submitted 14 November, 2021; v1 submitted 1 May, 2021;
originally announced May 2021.
-
Power Efficient Communication for Low Signal to Noise Ratio Optical Links
Authors:
Ravikiran Kakarla,
Mikael Mazur,
Jochen Schröder,
Peter A. Andrekson
Abstract:
Receiver sensitivity is a particularly important metric in optical communication links operating at low signal-to-noise ratios (SNRs), for example in deep-space communication, since it directly limits the maximum achievable reach and data rate. Pulse position modulation (PPM) with direct detection photon-counting detectors are the most power-efficient solution known, however, the sensitivity gain…
▽ More
Receiver sensitivity is a particularly important metric in optical communication links operating at low signal-to-noise ratios (SNRs), for example in deep-space communication, since it directly limits the maximum achievable reach and data rate. Pulse position modulation (PPM) with direct detection photon-counting detectors are the most power-efficient solution known, however, the sensitivity gain comes at the expense of reduced spectral efficiency. We show that quadrature phase-shift keying (QPSK) modulation with a phase-sensitive ultralow noise pre-amplified coherent receiver outperforms other well-known power-efficient multi-dimensional coherent modulation formats, while simultaneously having higher spectral efficiency. It also results in better sensitivity than PPM for orders up to 64 with ideal direct detection using photon-counting receivers. This is because of the bit error rate characteristics favoring the QPSK format when forward error correction with a large overhead is considered.
△ Less
Submitted 20 April, 2021;
originally announced April 2021.
-
End-to-end Autoencoder for Superchannel Transceivers with Hardware Impairment
Authors:
Jinxiang Song,
Christian Häger,
Jochen Schröder,
Alexandre Graell i Amat,
Henk Wymeersch
Abstract:
We propose an end-to-end learning-based approach for superchannel systems impaired by non-ideal hardware component. Our system achieves up to 60% SER reduction and up to 50% guard band reduction compared with the considered baseline scheme.
We propose an end-to-end learning-based approach for superchannel systems impaired by non-ideal hardware component. Our system achieves up to 60% SER reduction and up to 50% guard band reduction compared with the considered baseline scheme.
△ Less
Submitted 29 March, 2021;
originally announced March 2021.
-
Widely Tunable, Low Linewidth, and High Power Laser Source using an Electro-Optic Comb and Injection-Locked Slave Laser Array
Authors:
J Connor Skehan,
Corentin Naveau,
Jochen Schroder,
Peter Andrekson
Abstract:
We propose a simple approach to implement a tunable, high power and narrow linewidth laser source based on a series of highly coherent tones from an electro-optic frequency comb and a set of 3 DFB slave lasers. We experimentally demonstrate approximately 1.25 THz (10 nm) of tuning within the C-Band centered at 192.9 THz (1555 nm). The output power is approximately 100 mW (20 dBm), with a side band…
▽ More
We propose a simple approach to implement a tunable, high power and narrow linewidth laser source based on a series of highly coherent tones from an electro-optic frequency comb and a set of 3 DFB slave lasers. We experimentally demonstrate approximately 1.25 THz (10 nm) of tuning within the C-Band centered at 192.9 THz (1555 nm). The output power is approximately 100 mW (20 dBm), with a side band suppression ratio greater than 55 dB, and a linewidth below 400 Hz across the full range of tunability. This approach is scalable and may be extended to cover a significantly broader optical spectral range.
△ Less
Submitted 3 March, 2021;
originally announced March 2021.
-
Parallel Velocity Mixing Yielding Enhanced Electron Heating During Magnetic Pumping
Authors:
J. Egedal,
J. Schroeder,
E. Lichko
Abstract:
Magnetic wave perturbations are observed in the solar wind and in the vicinity of Earth's bow shock. For such environments, recent work on magnetic pumping with electrons trapped in the magnetic perturbations have demonstrated the possibility of efficient energization of superthermal electrons. Here we also analyze the energization of such energetic electrons for which the transit time through the…
▽ More
Magnetic wave perturbations are observed in the solar wind and in the vicinity of Earth's bow shock. For such environments, recent work on magnetic pumping with electrons trapped in the magnetic perturbations have demonstrated the possibility of efficient energization of superthermal electrons. Here we also analyze the energization of such energetic electrons for which the transit time through the system is short compared to time scales associated with the magnetic field evolution. In particular, considering an idealized magnetic configuration we show how trapping/detrapping of energetic magnetized electrons can cause effective parallel velocity diffusion. This parallel diffusion, combined with naturally occurring mechanisms known to cause pitch angle scattering, such as Whistler waves, produces enhanced heating rates for magnetic pumping. We find that at low pitch angle scattering rates the combined mechanism enhances the heating beyond the predictions of the recent theory for magnetic pumping with trapped electrons.
△ Less
Submitted 11 February, 2021;
originally announced February 2021.
-
Demonstration of a uniform, high-pressure, high-temperature gas cell with a dual frequency comb absorption spectrometer
Authors:
Ryan K. Cole,
Anthony D. Draper,
Paul J. Schroeder,
Cameron M. Casby,
Amanda S. Makowiecki,
Sean C. Coburn,
Julie E. Steinbrenner,
Nazanin Hoghooghi,
Gregory B. Rieker
Abstract:
Accurate absorption models for gases at high pressure and temperature support advanced optical combustion diagnostics and aid in the study of harsh planetary atmospheres. Developing and validating absorption models for these applications requires recreating the extreme temperature and pressure conditions of these environments in static, uniform, well-known conditions in the laboratory. Here, we pr…
▽ More
Accurate absorption models for gases at high pressure and temperature support advanced optical combustion diagnostics and aid in the study of harsh planetary atmospheres. Developing and validating absorption models for these applications requires recreating the extreme temperature and pressure conditions of these environments in static, uniform, well-known conditions in the laboratory. Here, we present the design of a new gas cell to enable reference-quality absorption spectroscopy at high pressure and temperature. The design centers on a carefully controlled quartz sample cell housed at the core of a pressurized ceramic furnace. The half-meter sample cell is relatively long compared to past high- pressure and temperature absorption cells, and is surrounded by a molybdenum heat spreader that enables high temperature uniformity over the full length of the absorbing gas. We measure the temperature distribution of the sample gas using in situ thermocouples, and fully characterize the temperature uniformity across a full matrix of conditions up to 1000 K and 50 bar. The results demonstrate that the new design enables highly uniform and precisely known conditions across the full absorbing path length. Uniquely, we test the new gas cell with a broadband, high-resolution dual frequency comb spectrometer that enables highly resolved absorption spectroscopy across a wide range of temperature and pressure conditions. With this system, we measure the spectrum of CO$_2$ between 6800 and 7000 cm$^{-1}$ at pressures between 0.2 and 20 bar, and temperatures up to 1000 K. The measurements reveal discrepancies from spectra predicted by the HITRAN2016 database with a Voigt line shape at both low- and high-pressure conditions. These results motivate future work to expand absorption models and databases to accurately model high-pressure and -temperature spectra in combustion and planetary science research.
△ Less
Submitted 8 February, 2021;
originally announced February 2021.
-
The impact of magnetic fields on momentum transport and saturation of shear-flow instability by stable modes
Authors:
A. E. Fraser,
P. W. Terry,
E. G. Zweibel,
M. J. Pueschel,
J. M. Schroeder
Abstract:
The Kelvin-Helmholtz (KH) instability of a shear layer with an initially-uniform magnetic field in the direction of flow is studied in the framework of 2D incompressible magnetohydrodynamics with finite resistivity and viscosity using direct numerical simulations. The shear layer evolves freely, with no external forcing, and thus broadens in time as turbulent stresses transport momentum across it.…
▽ More
The Kelvin-Helmholtz (KH) instability of a shear layer with an initially-uniform magnetic field in the direction of flow is studied in the framework of 2D incompressible magnetohydrodynamics with finite resistivity and viscosity using direct numerical simulations. The shear layer evolves freely, with no external forcing, and thus broadens in time as turbulent stresses transport momentum across it. As with KH-unstable flows in hydrodynamics, the instability here features a conjugate stable mode for every unstable mode in the absence of dissipation. Stable modes are shown to transport momentum up its gradient, shrinking the layer width whenever they exceed unstable modes in amplitude. In simulations with weak magnetic fields, the linear instability is minimally affected by the magnetic field, but enhanced small-scale fluctuations relative to the hydrodynamic case are observed. These enhanced fluctuations coincide with increased energy dissipation and faster layer broadening, with these features more pronounced in simulations with stronger fields. These trends result from the magnetic field reducing the effects of stable modes relative to the transfer of energy to small scales. As field strength increases, stable modes become less excited and thus transport less momentum against its gradient. Furthermore, the energy that would otherwise transfer back to the driving shear due to stable modes is instead allowed to cascade to small scales, where it is lost to dissipation. Approximations of the turbulent state in terms of a reduced set of modes are explored. While the Reynolds stress is well-described using just two modes per wavenumber at large scales, the Maxwell stress is not.
△ Less
Submitted 21 October, 2020;
originally announced October 2020.