An efficient edge based data structure for the compressible Reynolds-averaged Navier–Stokes equations on hybrid unstructured meshes

An efficient edge based data structure has been developed in order to implement an unstructured vertex based finite volume algorithm for the Reynolds-averaged Navier–Stokes equations on hybrid meshes. In the present approach, the data structure is tailored to meet the requirements of the vertex based algorithm by considering data access patterns and cache efficiency. The required data are packed and allocated in a way that they are close to each other in the physical memory. Therefore, the proposed data structure increases cache performance and improves computation time. As a result, the explicit flow solver indicates a significant speed up compared to other open-source solvers in terms of CPU time. A fully implicit version has also been implemented based on the PETSc library in order to improve the robustness of the algorithm. The resulting algebraic equations due to the compressible Navier–Stokes and the one equation Spalart–Allmaras turbulence equations are solved in a monolithic manner using the restricted additive Schwarz preconditioner combined with the FGMRES Krylov subspace algorithm. In order to further improve the computational accuracy, the multiscale metric based anisotropic mesh refinement library PyAMG is used for mesh adaptation. The numerical algorithm is validated for the classical benchmark problems such as the transonic turbulent flow around a supercritical RAE2822 airfoil and DLR-F6 wing-body-nacelle-pylon configuration. The efficiency of the data structure is demonstrated by achieving up to an order of magnitude speed up in CPU times.     

Versatile Near-Infrared Super-Resolution Imaging of Amyloid Fibrils with the Fluorogenic Probe CRANAD-2

CRANAD-2 is a fluorogenic curcumin derivative used for near-infrared detection and imaging in vivo of amyloid aggregates, which are involved in neurodegenerative diseases. We explore the performance of CRANAD-2 in two super-resolution imaging techniques, namely stimulated emission depletion (STED) and single-molecule localization microscopy (SMLM), with markedly different fluorophore requirements. By conveniently adapting the concentration of CRANAD-2, which transiently binds to amyloid fibrils, we show that it performs well in both techniques, achieving a resolution in the range of 45–55 nm. Correlation of SMLM with atomic force microscopy (AFM) validates the resolution of fine features in the reconstructed super-resolved image. The good performance and versatility of CRANAD-2 provides a powerful tool for near-infrared nanoscopic imaging of amyloids in vitro and in vivo.     

A posteriori error analysis of a quadratic finite volume method for nonlinear elliptic problems

In this article, we construct and analyze a residual-based a posteriori error estimator for a quadratic finite volume method (FVM) for solving nonlinear elliptic partial differential equations with homogeneous Dirichlet boundary conditions. We shall prove that the a posteriori error estimator yields the global upper and local lower bounds for the norm error of the FVM. So that the a posteriori error estimator is equivalent to the true error in a certain sense. Numerical experiments are performed to illustrate the theoretical results.     

Concentrated steady vorticities of the Euler equation on 2-d domains and their linear stability

We consider concentrated vorticities for the Euler equation on a smooth domain $\mathrm{\Omega }?{\mathbf{R}}^2$ in the form of

$\omega =\sum _{j=1}^N{\omega }_j{\chi }_{{\mathrm{\Omega }}_j},|{\mathrm{\Omega }}_j|=\pi r_j^2,\underset{{\mathrm{\Omega }}_j}{\int }{\omega }_{j}d\mu ={\mu }_j\ne 0,$

supported on well-separated vortical domains ${\mathrm{\Omega }}_j$, $j=1,\dots ,N$, of small diameters $O(r_j)$. A conformal mapping framework is set up to study this free boundary problem with ${\mathrm{\Omega }}_j$ being part of unknowns. For any given vorticities ${\mu }_1,\dots ,{\mu }_N$ and small $r_1,\dots ,r_N\in {\mathbf{R}}^{+}$, through a perturbation approach, we obtain such piecewise constant steady vortex patches as well as piecewise smooth Lipschitz steady vorticities, both concentrated near non-degenerate critical configurations of the Kirchhoff–Routh Hamiltonian function. When vortex patch evolution is considered as the boundary dynamics of $?{\mathrm{\Omega }}_j$, through an invariant subspace decomposition, it is also proved that the spectral/linear stability of such steady vortex patches is largely determined by that of the 2N-dimensional linearized point vortex dynamics, while the motion is highly oscillatory in the 2N-codim directions corresponding to the vortical domain shapes.     

Maximal Subrings and Covering Numbers of Finite Semisimple Rings

We classify the maximal subrings of the ring of n×n matrices over a finite field, and show that these subrings may be divided into three types. We also describe all of the maximal subrings of a finite semisimple ring, and categorize them into two classes. As an application of these results, we calculate the covering number of a finite semisimple ring.     

A simple three‐wave approximate Riemann solver for the Saint‐Venant‐Exner equations

Erosion and sediments transport processes have a great impact on industrial structures and on water quality. Despite its limitations, the Saint‐Venant‐Exner system is still (and for sure for some years) widely used in industrial codes to model the bedload sediment transport. In practice, its numerical resolution is mostly handled by a splitting technique that allows a weak coupling between hydraulic and morphodynamic distinct softwares but may suffer from important stability issues. In recent works, many authors proposed alternative methods based on a strong coupling that cure this problem but are not so trivial to implement in an industrial context. In this work, we then pursue 2 objectives. First, we propose a very simple scheme based on an approximate Riemann solver, respecting the strong coupling framework, and we demonstrate its stability and accuracy through a number of numerical test cases. However, second, we reinterpret our scheme as a splitting technique and we extend the purpose to propose what should be the minimal coupling that ensures the stability of the global numerical process in industrial codes, at least, when dealing with collocated finite volume method. The resulting splitting method is, up to our knowledge, the only one for which stability properties are fully demonstrated.     

A generalized super AKNS hierarchy associated with orthosymplectic Lie superalgebra OSP(2,2) and its super bi-Hamiltonian structures

For the orthosymplectic Lie superalgebra $\text{◂⋅▸}OSP(2,2)$, we choose a set of basis matrices. A linear combination of those basis matrices presents a spatial spectral matrix. The compatible condition of the spatial part and the corresponding temporal parts of the spectral problem leads to a generalized super AKNS (GSAKNS) hierarchy. By making use of the supertrace identity, the obtained GSAKNS hierarchy can be written as the super bi-Hamiltonian structures.     

A simple,rapid and sensitive method based on modified multiwalled carbon nanotube for preconcentration and determination of lead ions in aqueous media in natural pHs

In this study, multiwalled carbon nanotube (MWCNT) was modified by the pyridine group using a silane agent and characterized by infrared spectroscopy (IR), thermal analysis (TG/DTA), and elemental analysis (CHN) and scanning electron microscopy (SEM). The application of this sorbent was investigated in determination of lead ions in aqueous samples, using flame atomic absorption spectrometry (FAAS). Through this study, different parameters such as pH and sample flow rate on adsorption process and eluent concentration, volume and flow rate were optimized. The limit of detection (LOD), the relative standard deviation and the recovery of the method were 2 ng mL^?1, 1.3% and 99.7%, respectively. Two standard reference materials (NIST 1571 and NIST 1572) were used to verify accuracy of this method. Finally, the sorbent was successfully applied for extraction and determination of low levels of Pb(II) ions in aqueous samples.     

Time-dependent Ginzburg–Landau equations for multi-gap superconductors

We numerically solve the time-dependent Ginzburg–Landau equations for two-gap superconductors using the finite-element technique. The real-time simulation shows that at low magnetic field, the vortices in small-size samples tend to form clusters or other disorder structures. When the sample size is large, stripes appear in the pattern. These results are in good agreement with the previous experimental observations of the intriguing anomalous vortex pattern, providing a reliable theoretical basis for the future applications of multi-gap superconductors.