Source author record

Nicolae Suciu

Nicolae Suciu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.geo-ph math.NA Numerical Analysis Computational Engineering, Finance, and Science physics.comp-ph physics.flu-dyn

Catalog footprint

What is connected

5works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Global random walk solvers for reactive transport and biodegradation processes in heterogeneous porous media

Flow and multicomponent reactive transport in saturated/unsaturated porous media are modeled by ensembles of computational particles moving on regular lattices according to specific random walk rules. The occupation number of the lattice sites is updated with a global random walk (GRW) procedure which spreads the particles from a lattice site with computational costs comparable to those for a single random walk step in sequential procedures. To cope with the nonlinearity and the degeneracy of the Richards equation the GRW flow solver uses linearization techniques similar to the $L$-scheme developed in finite element/volume approaches. Numerical schemes for reactive transport, coupled with the flow solver via numerical solutions for saturation and water flux, are implemented in splitting procedures. Diffusion-advection steps are solved by GRW algorithms using either biased or unbiased random walk probabilities. Since the number of particles in GRW simulations can be as large as the number of molecules involved in chemical reactions, one avoids the cumbersome problem of rescaling particle densities to approximate concentrations. Reaction steps are therefore formulated in terms of concentrations, as in deterministic approaches. The numerical convergence of the new schemes is demonstrated by comparisons with manufactured analytical solutions. Coupled flow and reactive transport problems of contaminant biodegradation described by the Monod model are further solved and the influence of flow nonlinearity/degeneracy and of the spatial heterogeneity of the medium is investigated numerically.

preprint2020arXiv

Benchmark for numerical solutions of flow in heterogeneous groundwater formations

This article presents numerical investigations on accuracy and convergence properties of several numerical approaches for simulating steady state flows in heterogeneous aquifers. Finite difference, finite element, discontinuous Galerkin, spectral, and random walk methods are tested on one- and two-dimensional benchmark flow problems. Realizations of log-normal hydraulic conductivity fields are generated by Kraichnan algorithms in closed form as finite sums of random periodic modes, which allow direct code verification by comparisons with manufactured reference solutions. The quality of the methods is assessed for increasing number of random modes and for increasing variance of the log-hydraulic conductivity fields with Gaussian and exponential correlation. Experimental orders of convergence are calculated from successive refinements of the grid. The numerical methods are further validated by comparisons between statistical inferences obtained from Monte Carlo ensembles of numerical solutions and theoretical first-order perturbation results. It is found that while for Gaussian correlation of the log-conductivity field all the methods perform well, in the exponential case their accuracy deteriorates and, for large variance and number of modes, the benchmark problems are practically not tractable with reasonably large computing resources, for all the methods considered in this study.

preprint2016arXiv

Building a Bridge from Moments to PDF's: A New Approach to Finding PDF Mixing Models

Probability density function (PDF) methods are a promising alternative to predicting the transport of solutes in groundwater under uncertainty. They make it possible to derive the evolution equations of the mean concentration and the concentration variance, used in moment methods. A mixing model, also known as a dissipation model, is essential for both methods. Finding a satisfactory mixing model is still an open question and due to the rather elaborate PDF methods, a difficult undertaking. Both the PDF equation and the concentration variance equation depend on the same mixing model. This connection is used to find and test an improved mixing model for the much easier to handle concentration variance. Subsequently, this mixing model is transferred to the PDF equation and tested. The newly proposed mixing model yields significantly improved results for both variance modelling and PDF modelling.

preprint2014arXiv

Geophysical implications of a decentered inner core

In a first approximation, the Earth's interior has an isotropic structure with a spherical symmetry. Over the last decades the geophysical observations have revealed, at different spatial scales, the existence of several perturbations from this basic structure. In this paper we discuss the hemispheric perturbations induced to this basic structure if the inner core is displaced from the center of mass of the Earth. Using numerical simulations of the observed hemispheric asymmetry of the seismic waves traveling through the upper inner core, with faster arrival times and higher attenuation in the Eastern Hemisphere, we estimate that the present position of the inner core is shifted by tens of kilometers from the Earth's center eastward in the equatorial plane. If the only forces acting on the inner core were the gravitational forces, then its equilibrium position would be at the Earth's center and the estimated displacement would not be possible. We conjecture that, due to interactions with the flow and the magnetic field inside the outer core, the inner core is in a permanent chaotic motion. To support this hypothesis we analyze more than ten different geophysical phenomena consistent with an inner core motion dominated by time scales from hundreds to thousands of years.

preprint2011arXiv

Seismic hemispheric asymmetry induced by Earth's inner core decentering

In a first approximation the Earth's interior has an isotropic structure with a spherical symmetry. Over the last decades the geophysical observations have revealed, at different spatial scales, the existence of several perturbations from this basic structure. Some of them are situated in the neighborhood of the inner core boundary (ICB). One of the best documented perturbations is the asymmetry at the top of the inner core (ATIC) characterized by faster seismic wave velocity in the eastern hemisphere than in the western hemisphere. All existing explanations are based on a hemispheric variation of the material properties near ICB inside the inner core. Using numerical simulations of the seismic ray propagation, we show that the ATIC can be explained as well by the displacement of the inner core towards east in the equatorial plane tens of kilometers from the Earth's center, without modifying the spherical symmetry in the upper inner core. The hypothesis of a displaced inner core is also sustained by other observed hemispheric asymmetries at the top of the inner core and at the bottom of the outer core. A displaced inner core would have major implications for many mechanical, thermal, and magnetic phenomena in the Earth's interior.