@article{MikuckiSchulerDigeletal.2023, author = {Mikucki, Jill Ann and Schuler, C. G. and Digel, Ilya and Kowalski, Julia and Tuttle, M. J. and Chua, Michelle and Davis, R. and Purcell, Alicia and Ghosh, D. and Francke, G. and Feldmann, M. and Espe, C. and Heinen, Dirk and Dachwald, Bernd and Clemens, Joachim and Lyons, W. B. and Tulaczyk, S.}, title = {Field-Based planetary protection operations for melt probes: validation of clean access into the blood falls, antarctica, englacial ecosystem}, series = {Astrobiology}, volume = {23}, journal = {Astrobiology}, number = {11}, publisher = {Liebert}, address = {New York, NY}, issn = {1557-8070 (online)}, doi = {10.1089/ast.2021.0102}, pages = {1165 -- 1178}, year = {2023}, abstract = {Subglacial environments on Earth offer important analogs to Ocean World targets in our solar system. These unique microbial ecosystems remain understudied due to the challenges of access through thick glacial ice (tens to hundreds of meters). Additionally, sub-ice collections must be conducted in a clean manner to ensure sample integrity for downstream microbiological and geochemical analyses. We describe the field-based cleaning of a melt probe that was used to collect brine samples from within a glacier conduit at Blood Falls, Antarctica, for geomicrobiological studies. We used a thermoelectric melting probe called the IceMole that was designed to be minimally invasive in that the logistical requirements in support of drilling operations were small and the probe could be cleaned, even in a remote field setting, so as to minimize potential contamination. In our study, the exterior bioburden on the IceMole was reduced to levels measured in most clean rooms, and below that of the ice surrounding our sampling target. Potential microbial contaminants were identified during the cleaning process; however, very few were detected in the final englacial sample collected with the IceMole and were present in extremely low abundances (∼0.063\% of 16S rRNA gene amplicon sequences). This cleaning protocol can help minimize contamination when working in remote field locations, support microbiological sampling of terrestrial subglacial environments using melting probes, and help inform planetary protection challenges for Ocean World analog mission concepts.}, language = {en} } @article{AkimbekovDigelTastambeketal.2024, author = {Akimbekov, Nuraly S. and Digel, Ilya and Tastambek, Kuanysh T. and Kozhahmetova, Marzhan and Sherelkhan, Dinara K. and Tauanov, Zhandos}, title = {Hydrogenotrophic methanogenesis in coal-bearing environments: Methane production, carbon sequestration, and hydrogen availability}, series = {International Journal of Hydrogen Energy}, volume = {52}, journal = {International Journal of Hydrogen Energy}, number = {Part D}, publisher = {Elsevier}, address = {New York}, issn = {1879-3487 (online)}, doi = {10.1016/j.ijhydene.2023.09.223}, pages = {1264 -- 1277}, year = {2024}, abstract = {Methane is a valuable energy source helping to mitigate the growing energy demand worldwide. However, as a potent greenhouse gas, it has also gained additional attention due to its environmental impacts. The biological production of methane is performed primarily hydrogenotrophically from H2 and CO2 by methanogenic archaea. Hydrogenotrophic methanogenesis also represents a great interest with respect to carbon re-cycling and H2 storage. The most significant carbon source, extremely rich in complex organic matter for microbial degradation and biogenic methane production, is coal. Although interest in enhanced microbial coalbed methane production is continuously increasing globally, limited knowledge exists regarding the exact origins of the coalbed methane and the associated microbial communities, including hydrogenotrophic methanogens. Here, we give an overview of hydrogenotrophic methanogens in coal beds and related environments in terms of their energy production mechanisms, unique metabolic pathways, and associated ecological functions.}, language = {en} } @article{ZhantlessovaSavitskayaKistaubayevaetal.2024, author = {Zhantlessova, Sirina and Savitskaya, Irina and Kistaubayeva, Aida and Ignatova, Ludmila and Talipova, Aizhan and Pogrebnjak, Alexander and Digel, Ilya}, title = {Correction: Zhantlessova et al. advanced "Green" prebiotic composite of bacterial cellulose/pullulan based on synthetic biology-powered microbial coculture strategy. Polymers 2022, 14, 3224}, series = {Polymers}, volume = {16}, journal = {Polymers}, number = {13}, publisher = {MDPI}, address = {Basel}, issn = {2073-4360}, doi = {10.3390/polym16131802}, pages = {2 Seiten}, year = {2024}, language = {en} } @article{KleefeldPieronek2020, author = {Kleefeld, Andreas and Pieronek, J.}, title = {Elastic transmission eigenvalues and their computation via the method of fundamental solutions}, series = {Applicable Analysis}, volume = {100}, journal = {Applicable Analysis}, number = {16}, publisher = {Taylore \& Francis}, address = {London}, issn = {1563-504X}, doi = {10.1080/00036811.2020.1721473}, pages = {3445 -- 3462}, year = {2020}, abstract = {A stabilized version of the fundamental solution method to catch ill-conditioning effects is investigated with focus on the computation of complex-valued elastic interior transmission eigenvalues in two dimensions for homogeneous and isotropic media. Its algorithm can be implemented very shortly and adopts to many similar partial differential equation-based eigenproblems as long as the underlying fundamental solution function can be easily generated. We develop a corroborative approximation analysis which also implicates new basic results for transmission eigenfunctions and present some numerical examples which together prove successful feasibility of our eigenvalue recovery approach.}, language = {en} } @article{AyalaHarrisKleefeld2024, author = {Ayala, Rafael Ceja and Harris, Isaac and Kleefeld, Andreas}, title = {Direct sampling method via Landweber iteration for an absorbing scatterer with a conductive boundary}, series = {Inverse Problems and Imaging}, volume = {18}, journal = {Inverse Problems and Imaging}, number = {3}, publisher = {AIMS}, address = {Springfield}, issn = {1930-8337}, doi = {10.3934/ipi.2023051}, pages = {708 -- 729}, year = {2024}, abstract = {In this paper, we consider the inverse shape problem of recovering isotropic scatterers with a conductive boundary condition. Here, we assume that the measured far-field data is known at a fixed wave number. Motivated by recent work, we study a new direct sampling indicator based on the Landweber iteration and the factorization method. Therefore, we prove the connection between these reconstruction methods. The method studied here falls under the category of qualitative reconstruction methods where an imaging function is used to recover the absorbing scatterer. We prove stability of our new imaging function as well as derive a discrepancy principle for recovering the regularization parameter. The theoretical results are verified with numerical examples to show how the reconstruction performs by the new Landweber direct sampling method.}, language = {en} } @article{ClausnitzerKleefeld2024, author = {Clausnitzer, Julian and Kleefeld, Andreas}, title = {A spectral Galerkin exponential Euler time-stepping scheme for parabolic SPDEs on two-dimensional domains with a C² boundary}, series = {Discrete and Continuous Dynamical Systems - Series B}, volume = {29}, journal = {Discrete and Continuous Dynamical Systems - Series B}, number = {4}, publisher = {AIMS}, address = {Springfield}, issn = {1531-3492}, doi = {10.3934/dcdsb.2023148}, pages = {1624 -- 1651}, year = {2024}, abstract = {We consider the numerical approximation of second-order semi-linear parabolic stochastic partial differential equations interpreted in the mild sense which we solve on general two-dimensional domains with a C² boundary with homogeneous Dirichlet boundary conditions. The equations are driven by Gaussian additive noise, and several Lipschitz-like conditions are imposed on the nonlinear function. We discretize in space with a spectral Galerkin method and in time using an explicit Euler-like scheme. For irregular shapes, the necessary Dirichlet eigenvalues and eigenfunctions are obtained from a boundary integral equation method. This yields a nonlinear eigenvalue problem, which is discretized using a boundary element collocation method and is solved with the Beyn contour integral algorithm. We present an error analysis as well as numerical results on an exemplary asymmetric shape, and point out limitations of the approach.}, language = {en} } @article{HarrisKleefeld2022, author = {Harris, Isaac and Kleefeld, Andreas}, title = {Analysis and computation of the transmission eigenvalues with a conductive boundary condition}, series = {Applicable Analysis}, volume = {101}, journal = {Applicable Analysis}, number = {6}, publisher = {Taylor \& Francis}, address = {London}, issn = {1563-504X}, doi = {10.1080/00036811.2020.1789598}, pages = {1880 -- 1895}, year = {2022}, abstract = {We provide a new analytical and computational study of the transmission eigenvalues with a conductive boundary condition. These eigenvalues are derived from the scalar inverse scattering problem for an inhomogeneous material with a conductive boundary condition. The goal is to study how these eigenvalues depend on the material parameters in order to estimate the refractive index. The analytical questions we study are: deriving Faber-Krahn type lower bounds, the discreteness and limiting behavior of the transmission eigenvalues as the conductivity tends to infinity for a sign changing contrast. We also provide a numerical study of a new boundary integral equation for computing the eigenvalues. Lastly, using the limiting behavior we will numerically estimate the refractive index from the eigenvalues provided the conductivity is sufficiently large but unknown.}, language = {en} } @article{KleefeldZimmermann2022, author = {Kleefeld, Andreas and Zimmermann, M.}, title = {Computing Elastic Interior Transmission Eigenvalues}, series = {Integral Methods in Science and Engineering}, journal = {Integral Methods in Science and Engineering}, editor = {Constanda, Christian and Bodmann, Bardo E.J. and Harris, Paul J.}, publisher = {Birkh{\"a}user}, address = {Cham}, isbn = {978-3-031-07171-3}, doi = {10.1007/978-3-031-07171-3_10}, pages = {139 -- 155}, year = {2022}, abstract = {An alternative method is presented to numerically compute interior elastic transmission eigenvalues for various domains in two dimensions. This is achieved by discretizing the resulting system of boundary integral equations in combination with a nonlinear eigenvalue solver. Numerical results are given to show that this new approach can provide better results than the finite element method when dealing with general domains.}, language = {en} } @article{Kleefeld2021, author = {Kleefeld, Andreas}, title = {The hot spots conjecture can be false: some numerical examples}, series = {Advances in Computational Mathematics}, volume = {47}, journal = {Advances in Computational Mathematics}, publisher = {Springer}, address = {Dordrecht}, issn = {1019-7168}, doi = {10.1007/s10444-021-09911-5}, year = {2021}, abstract = {The hot spots conjecture is only known to be true for special geometries. This paper shows numerically that the hot spots conjecture can fail to be true for easy to construct bounded domains with one hole. The underlying eigenvalue problem for the Laplace equation with Neumann boundary condition is solved with boundary integral equations yielding a non-linear eigenvalue problem. Its discretization via the boundary element collocation method in combination with the algorithm by Beyn yields highly accurate results both for the first non-zero eigenvalue and its corresponding eigenfunction which is due to superconvergence. Additionally, it can be shown numerically that the ratio between the maximal/minimal value inside the domain and its maximal/minimal value on the boundary can be larger than 1 + 10- 3. Finally, numerical examples for easy to construct domains with up to five holes are provided which fail the hot spots conjecture as well.}, language = {en} } @article{MartinVaqueroKleefeld2020, author = {Mart{\´i}n-Vaquero, J. and Kleefeld, Andreas}, title = {Solving nonlinear parabolic PDEs in several dimensions: Parallelized ESERK codes}, series = {Journal of Computational Physics}, journal = {Journal of Computational Physics}, number = {423}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0021-9991}, doi = {10.1016/j.jcp.2020.109771}, year = {2020}, abstract = {There is a very large number of very important situations which can be modeled with nonlinear parabolic partial differential equations (PDEs) in several dimensions. In general, these PDEs can be solved by discretizing in the spatial variables and transforming them into huge systems of ordinary differential equations (ODEs), which are very stiff. Therefore, standard explicit methods require a large number of iterations to solve stiff problems. But implicit schemes are computationally very expensive when solving huge systems of nonlinear ODEs. Several families of Extrapolated Stabilized Explicit Runge-Kutta schemes (ESERK) with different order of accuracy (3 to 6) are derived and analyzed in this work. They are explicit methods, with stability regions extended, along the negative real semi-axis, quadratically with respect to the number of stages s, hence they can be considered to solve stiff problems much faster than traditional explicit schemes. Additionally, they allow the adaptation of the step length easily with a very small cost. Two new families of ESERK schemes (ESERK3 and ESERK6) are derived, and analyzed, in this work. Each family has more than 50 new schemes, with up to 84.000 stages in the case of ESERK6. For the first time, we also parallelized all these new variable step length and variable number of stages algorithms (ESERK3, ESERK4, ESERK5, and ESERK6). These parallelized strategies allow to decrease times significantly, as it is discussed and also shown numerically in two problems. Thus, the new codes provide very good results compared to other well-known ODE solvers. Finally, a new strategy is proposed to increase the efficiency of these schemes, and it is discussed the idea of combining ESERK families in one code, because typically, stiff problems have different zones and according to them and the requested tolerance the optimum order of convergence is different.}, language = {en} }