BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Ansgar Jüngel

Ansgar Jüngel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.APmath.NANumerical Analysismath.PRmath.FADistributed, Parallel, and Cluster Computingmath-phmath.MPMathematical SoftwarePerformancephysics.comp-ph

Catalog footprint

What is connected

41works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

41 published item(s)

preprint2026arXiv

A generalized Scharfetter-Gummel scheme for nonlocal cross-diffusion systems

An implicit Euler finite-volume scheme for a nonlocal cross-diffusion system on the multidimensional torus is analyzed. The equations describe the dynamics of population species with repulsive or attractive interactions. The numerical scheme is based on a generalized Scharfetter-Gummel discretization of the nonlocal flux term. For merely integrable kernel functions, the scheme preserves the positivity, total mass, and entropy structure. The existence of a discrete solution and its convergence to a solution to the continuous problem, as the mesh size tends to zero, are shown. A key difficulty is the degeneracy of the generalized Bernoulli function in the Scharfetter-Gummel approximation. This issue is overcome by proving a uniform estimate for the discrete Fisher information, which requires both the Boltzmann and Rao entropy inequalities. Numerical simulations illustrate the features of the scheme in one and two space dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

A convergent entropy-dissipating BDF2 finite-volume scheme for a population cross-diffusion system

A second-order backward differentiation formula (BDF2) finite-volume discretization for a nonlinear cross-diffusion system arising in population dynamics is studied. The numerical scheme preserves the Rao entropy structure and conserves the mass. The existence and uniqueness of discrete solutions and their large-time behavior as well as the convergence of the scheme are proved. The proofs are based on the G-stability of the BDF2 scheme, which provides an inequality for the quadratic Rao entropy and hence suitable a priori estimates. The novelty is the extension of this inequality to the system case. Some numerical experiments in one and two space dimensions underline the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A coupled stochastic differential reaction-diffusion system for angiogenesis

A coupled system of nonlinear mixed-type equations modeling early stages of angiogenesis is analyzed in a bounded domain. The system consists of stochastic differential equations describing the movement of the positions of the tip and stalk endothelial cells, due to chemotaxis, durotaxis, and random motion; ordinary differential equations for the volume fractions of the extracellular fluid, basement membrane, and fibrin matrix; and reaction-diffusion equations for the concentrations of several proteins involved in the angiogenesis process. The drift terms of the stochastic differential equations involve the gradients of the volume fractions and the concentrations, and the diffusivities in the reaction-diffusion equations depend nonlocally on the volume fractions, making the system highly nonlinear. The existence of a unique solution to this system is proved by using fixed-point arguments and Hölder regularity theory. Numerical experiments in two space dimensions illustrate the onset of formation of vessels.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Analysis of a fractional cross-diffusion system for multi-species populations

The global in time existence of weak solutions to a cross-diffusion system with fractional diffusion in the whole space is proved. The equations describe the evolution of multi-species populations in the regime of large-distance interactions; they have been derived in the many-particle limit from moderately interacting particle systems with Lévy noise. The existence proof is based on a three-level approximation scheme, entropy and moment estimates, and a new Aubin-Lions compactness lemma in the whole space.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Existence and weak-strong uniqueness for Maxwell-Stefan-Cahn-Hilliard systems

A Maxwell-Stefan system for fluid mixtures with driving forces depending on Cahn-Hilliard-type chemical potentials is analyzed. The corresponding parabolic cross-diffusion equations contain fourth-order derivatives and are considered in a bounded domain with no-flux boundary conditions. The main difficulty of the analysis is the degeneracy of the diffusion matrix, which is overcome by proving the positive definiteness of the matrix on a subspace and using the Bott--Duffin matrix inverse. The global existence of weak solutions and a weak-strong uniqueness property are shown by a careful combination of (relative) energy and entropy estimates, yielding $H^2(Ω)$ bounds for the densities, which cannot be obtained from the energy or entropy inequalities alone.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Shigesada-Kawasaki-Teramoto cross-diffusion system beyond detailed balance

The existence of global weak solutions to the cross-diffusion model of Shigesada, Kawasaki, and Teramoto for an arbitrary number of species is proved. The model consists of strongly coupled parabolic equations for the population densities in a bounded domain with no-flux boundary conditions, and it describes the dynamics of the segregation of the population species. The diffusion matrix is neither symmetric nor positive semidefinite. A new logarithmic entropy allows for an improved condition on the coefficients of heavily nonsymmetric diffusion matrices, without imposing the detailed-balance condition that is often assumed in the literature. Furthermore, the large-time convergence of the solutions to the constant steady state is proved by using the relative entropy associated to the logarithmic entropy.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Three-species drift-diffusion models for memristors

A system of drift-diffusion equations for the electron, hole, and oxygene vacancy densities in a semiconductor, coupled to the Poisson equation for the electric potential, is analyzed in a bounded domain with mixed Dirichlet-Neumann boundary conditions. This system describes the dynamics of charge carriers in a memristor device. Memristors can be seen as nonlinear resistors with memory, mimicking the conductance response of biological synapses. In the fast-relaxation limit, the system reduces to a drift-diffusion system for the oxygene vacancy density and electric potential, which is often used in neuromorphic applications. The following results are proved: the global existence of weak solutions to the full system in any space dimension; the uniform-in-time boundedness of the solutions to the full system and the fast-relaxation limit in two space dimensions; the global existence and weak-strong uniqueness analysis of the reduced system. Numerical experiments in one space dimension illustrate the behavior of the solutions and reproduce hysteresis effects in the current-voltage characteristics.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Random-batch method for multi-species stochastic interacting particle systems

A random-batch method for multi-species interacting particle systems is proposed, extending the method of S. Jin, L. Li, and J.-G. Liu [J. Comput. Phys. 400 (2020), 108877]. The idea of the algorithmus is to randomly divide, at each time step, the ensemble of particles into small batches and then to evolve the interaction of each particle within the batches until the next time step. This reduces the computational cost by one order of magnitude, while keeping a certain accuracy. It is proved that the $L^2$ error of the error process behaves like the square root of the time step size, uniformly in time, thus providing the convergence of the scheme. The numerical efficiency is tested for some examples, and numerical simulations of the opinion dynamics in a hierarchical company, consisting of workers, managers, and CEOs, are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A minimizing-movements approach to GENERIC systems

We present a new time discretization scheme adapted to the structure of GENERIC systems. The scheme is variational in nature and is based on a conditional incremental minimization. The GENERIC structure of the scheme provides stability and convergence of the scheme. We prove that the scheme can be rigorously implemented in the case of the damped harmonic oscillator. Numerical evidence is collected, illustrating the performance of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Convergence of a finite-volume scheme for a degenerate-singular cross-diffusion system for biofilms

An implicit Euler finite-volume scheme for a cross-diffusion system modeling biofilm growth is analyzed by exploiting its formal gradient-flow structure. The numerical scheme is based on a two-point flux approximation that preserves the entropy structure of the continuous model. Assuming equal diffusivities, the existence of nonnegative and bounded solutions to the scheme and its convergence are proved. Finally, we supplement the study by numerical experiments in one and two space dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Entropy-dissipating finite-difference schemes for nonlinear fourth-order parabolic equations

Structure-preserving finite-difference schemes for general nonlinear fourth-order parabolic equations on the one-dimensional torus are derived. Examples include the thin-film and the Derrida-Lebowitz-Speer-Spohn equations. The schemes conserve the mass and dissipate the entropy. The scheme associated to the logarithmic entropy also preserves the positivity. The idea of the derivation is to reformulate the equations in such a way that the chain rule is avoided. A central finite-difference discretization is then applied to the reformulation. In this way, the same dissipation rates as in the continuous case are recovered. The strategy can be extended to a multi-dimensional thin-film equation. Numerical examples in one and two space dimensions illustrate the dissipation properties.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantum drift-diffusion equations for a two-dimensional electron gas with spin-orbit interaction

Quantum drift-diffusion equations are derived for a two-dimensional electron gas with spin-orbit interaction of Rashba type. The (formal) derivation turns out to be a non-standard application of the usual mathematical tools, such as Wigner transform, Moyal product expansion and Chapman-Enskog expansion. The main peculiarity consists in the fact that a non-vanishing current is already carried by the leading-order term in the Chapman-Enskog expansion. To our knowledge, this is the first example of quantum drift-diffusion equations involving the full spin vector. Indeed, previous models were either quantum bipolar (involving only the spin projection on a given axis) or full spin but semiclassical.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Global martingale solutions for a stochastic population cross-diffusion system

The existence of global nonnegative martingale solutions to a stochastic cross-diffusion system for an arbitrary but finite number of interacting population species is shown. The random influence of the environment is modeled by a multiplicative noise term. The diffusion matrix is generally neither symmetric nor positive definite, but it possesses a quadratic entropy structure. This structure allows us to work in a Hilbert space framework and to apply a stochastic Galerkin method. The existence proof is based on energy-type estimates, the tightness criterion of Brzeźniak and co-workers, and Jakubowski's generalization of the Skorokhod theorem. The nonnegativity is proved by an extension of Stampacchia's truncation method due to Chekroun, Park, and Temam.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Discrete Bochner inequalities via the Bochner-Bakry-Emery approach for Markov chains

Discrete convex Sobolev inequalities and Beckner inequalities are derived for time-continuous Markov chains on finite state spaces. Beckner inequalities interpolate between the modified logarithmic Sobolev inequality and the Poincaré inequality. Their proof is based on the Bakry-Emery approach and on discrete Bochner-type inequalities established by Caputo, Dai Pra, and Posta and recently extended by Fathi and Maas for logarithmic entropies. The abstract result for convex entropies is applied to several Markov chains, including birth-death processes, zero-range processes, Bernoulli-Laplace models, and random transposition models, and to a finite-volume discretization of a one-dimensional Fokker-Planck equation, applying results by Mielke.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Displacement convexity for the entropy in semidiscrete nonlinear Fokker-Planck equations

The displacement $λ$-convexity of a nonstandard entropy with respect to a nonlocal transportation metric in finite state spaces is shown using a gradient flow approach. The constant $λ$ is computed explicitly in terms of a priori estimates of the solution to a finite-difference approximation of a nonlinear Fokker-Planck equation. The key idea is to employ a new mean function, which defines the Onsager operator in the gradient flow formulation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Energy-transport models for spin transport in ferromagnetic semiconductors

Explicit energy-transport equations for the spinorial carrier transport in ferromagnetic semiconductors are calculated from a general spin energy-transport system that was derived by Ben Abdallah and El Hajj from a spinorial Boltzmann equation. The novelty of our approach are the simplifying assumptions leading to explicit models which extend both spin drift-diffusion and semiclassical energy-transport equations. The explicit models allow us to examine the interplay between the spin and charge degrees of freedom. In particular, the monotonicity of the entropy (or free energy) and gradient estimates are shown for these models and the existence of weak solutions to a time-discrete version of one of the models is proved, using novel truncation arguments. Numerical experiments in one-dimensional multilayer structures using a finite-volume discretization illustrate the effect of the temperature and the polarization parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Existence analysis of a single-phase flow mixture model with van der Waals pressure

The transport of single-phase fluid mixtures in porous media is described by cross-diffusion equations for the mass densities. The equations are obtained in a thermodynamic consistent way from mass balance, Darcy's law, and the van der Waals equation of state for mixtures. The model consists of parabolic equations with cross diffusion with a hypocoercive diffusion operator. The global-in-time existence of weak solutions in a bounded domain with equilibrium boundary conditions is proved, extending the boundedness-by-entropy method. Based on the free energy inequality, the large-time convergence of the solution to the constant equilibrium mass density is shown. For the two-species model and specific diffusion matrices, an integral inequality is proved, which reveals a minimum principle for the mass fractions. Without mass diffusion, the two-dimensional pressure is shown to converge exponentially fast to a constant. Numerical examples in one space dimension illustrate this convergence.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Pipelined Iterative Solvers with Kernel Fusion for Graphics Processing Units

We revisit the implementation of iterative solvers on discrete graphics processing units and demonstrate the benefit of implementations using extensive kernel fusion for pipelined formulations over conventional implementations of classical formulations. The proposed implementations with both CUDA and OpenCL are freely available in ViennaCL and are shown to be competitive with or even superior to other solver packages for graphics processing units. Highest performance gains are obtained for small to medium-sized systems, while our implementations are on par with vendor-tuned implementations for very large systems. Our results are especially beneficial for transient problems, where many small to medium-sized systems instead of a single big system need to be solved.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Finite-Volume Scheme for a Spinorial Matrix Drift-Diffusion Model for Semiconductors

An implicit Euler finite-volume scheme for a spinorial matrix drift-diffusion model for semiconductors is analyzed. The model consists of strongly coupled parabolic equations for the electron density matrix or, alternatively, of weakly coupled equations for the charge and spin-vector densities, coupled to the Poisson equation for the elec-tric potential. The equations are solved in a bounded domain with mixed Dirichlet-Neumann boundary conditions. The charge and spin-vector fluxes are approximated by a Scharfetter-Gummel discretization. The main features of the numerical scheme are the preservation of positivity and L $\infty$ bounds and the dissipation of the discrete free energy. The existence of a bounded discrete solution and the monotonicity of the discrete free energy are proved. For undoped semiconductor materials, the numerical scheme is uncon-ditionally stable. The fundamental ideas are reformulations using spin-up and spin-down densities and certain projections of the spin-vector density, free energy estimates, and a discrete Moser iteration. Furthermore, numerical simulations of a simple ferromagnetic-layer field-effect transistor in two space dimensions are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A higher-order gradient flow scheme for a singular one-dimensional diffusion equation

A nonlinear diffusion equation, interpreted as a Wasserstein gradient flow, is numerically solved in one space dimension using a higher-order minimizing movement scheme based on the BDF (backward differentiation formula) discretization. In each time step, the approximation is obtained as the solution of a constrained quadratic minimization problem on a finite-dimensional space consisting of piecewise quadratic basis functions. The numerical scheme conserves the mass and dissipates the $G$-norm of the two-step BDF time approximation. Numerically, also the discrete entropy and variance are decaying. The decay turns out to be exponential in all cases. The corresponding decay rates are computed numerically for various grid numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Analysis of a coupled spin drift-diffusion Maxwell-Landau-Lifshitz system

The existence of global weak solutions to a coupled spin drift-diffusion and Maxwell-Landau-Lifshitz system is proved. The equations are considered in a two-dimensional magnetic layer structure and are supplemented with Dirichlet-Neumann boundary conditions. The spin drift-diffusion model for the charge density and spin density vector is the diffusion limit of a spinorial Boltzmann equation for a vanishing spin polarization constant. The Maxwell-Landau-Lifshitz system consists of the time-dependent Maxwell equations for the electric and magnetic fields and of the Landau-Lifshitz-Gilbert equation for the local magnetization, involving the interaction between magnetization and spin density vector. The existence proof is based on a regularization procedure, $L^2$-type estimates, and Moser-type iterations which yield the boundedness of the charge and spin densities. Furthermore, the free energy is shown to be nonincreasing in time if the magnetization-spin interaction constant in the Landau-Lifshitz equation is sufficiently small.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Analysis of degenerate cross-diffusion population models with volume filling

A class of parabolic cross-diffusion systems modeling the interaction of an arbitrary number of population species is analyzed in a bounded domain with no-flux boundary conditions. The equations are formally derived from a random-walk lattice model in the diffusion limit. Compared to previous results in the literature, the novelty is the combination of general degenerate diffusion and volume-filling effects. Conditions on the nonlinear diffusion coefficients are identified, which yield a formal gradient-flow or entropy structure. This structure allows for the proof of global-in-time existence of bounded weak solutions and the exponential convergence of the solutions to the constant steady state. The existence proof is based on an approximation argument, the entropy inequality, and new nonlinear Aubin-Lions compactness lemmas. The proof of the large-time behavior employs the entropy estimate and convex Sobolev inequalities. Moreover, under simplifiying assumptions on the nonlinearities, the uniqueness of weak solutions is shown by using the $H^{-1}$ method, the $E$-monotonicity technique of Gajewski, and the subadditivity of the Fisher information.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Entropy-dissipating semi-discrete Runge-Kutta schemes for nonlinear diffusion equations

Semi-discrete Runge-Kutta schemes for nonlinear diffusion equations of parabolic type are analyzed. Conditions are determined under which the schemes dissipate the discrete entropy locally. The dissipation property is a consequence of the concavity of the difference of the entropies at two consecutive time steps. The concavity property is shown to be related to the Bakry-Emery approach and the geodesic convexity of the entropy. The abstract conditions are verified for quasilinear parabolic equations (including the porous-medium equation), a linear diffusion system, and the fourth-order quantum diffusion equation. Numerical experiments for various Runge-Kutta finite-difference discretizations of the one-dimensional porous-medium equation show that the entropy-dissipation property is in fact global.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hypocoercivity for a linearized multi-species Boltzmann system

A new coercivity estimate on the spectral gap of the linearized Boltzmann collision operator for multiple species is proved. The assumptions on the collision kernels include hard and Maxwellian potentials under Grad's angular cut-off condition. Two proofs are given: a non-constructive one, based on the decomposition of the collision operator into a compact and a coercive part, and a constructive one, which exploits the "cross-effects" coming from collisions between different species and which yields explicit constants. Furthermore, the essential spectra of the linearized collision operator and the linearized Boltzmann operator are calculated. Based on the spectral-gap estimate, the exponential convergence towards global equilibrium with explicit rate is shown for solutions to the linearized multi-species Boltzmann system on the torus. The convergence is achieved by the interplay between the dissipative collision operator and the conservative transport operator and is proved by using the hypocoercivity method of Mouhot and Neumann.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Qualitative behavior of solutions to cross-diffusion systems from population dynamics

A general class of cross-diffusion systems for two population species in a bounded domain with no-flux boundary conditions and Lotka-Volterra-type source terms is analyzed. Although the diffusion coefficients are assumed to depend linearly on the population densities, the equations are strongly coupled. Generally, the diffusion matrix is neither symmetric nor positive definite. Three main results are proved: the existence of global uniformly bounded weak solutions, their convergence to the constant steady state in the weak competition case, and the uniqueness of weak solutions. The results hold under appropriate conditions on the diffusion parameters which are made explicit and which contain simplified Shigesada-Kawasaki-Teramoto population models as a special case. The proofs are based on entropy methods, which rely on convexity properties of suitable Lyapunov functionals.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A degenerate fourth-order parabolic equation modeling Bose-Einstein condensation. Part I: Local existence of solutions

A degenerate fourth-order parabolic equation modeling condensation phenomena related to Bose-Einstein particles is analyzed. The model is a Fokker-Planck-type approximation of the Boltzmann-Nordheim equation, only keeping the leading order term. It maintains some of the main features of the kinetic model, namely mass and energy conservation and condensation at zero energy. The existence of a local-in-time nonnegative continuous weak solution is proven. If the solution is not global, it blows up with respect to the $L^\infty$ norm in finite time. The proof is based on approximation arguments, interpolation inequalities in weighted Sobolev spaces, and suitable a priori estimates for a weighted gradient $L^2$ norm.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A degenerate fourth-order parabolic equation modeling Bose-Einstein condensation. Part II: Finite-time blow-up

A degenerate fourth-order parabolic equation modeling condensation phenomena related to Bose-Einstein particles is analyzed. The model is a Fokker-Planck-type approximation of the Boltzmann-Nordheim equation, only keeping the leading order term. It maintains some of the main features of the kinetic model, namely mass and energy conservation and condensation at zero energy. The existence of local-in-time weak solutions satisfying a certain entropy inequality is proven. The main result asserts that if a weighted $L^1$ norm of the initial data is sufficiently large and the initial data satisfies some integrability conditions, the solution blows up with respect to the $L^\infty$ norm in finite time. Furthermore, the set of all such blow-up enforcing initial functions is shown to be dense in the set of all admissible initial data. The proofs are based on approximation arguments and interpolation inequalities in weighted Sobolev spaces. By exploiting the entropy inequality, a nonlinear integral inequality is proved which implies the finite-time blow-up property.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Boundedness of weak solutions to cross-diffusion systems from population dynamics

The global-in-time existence of nonnegative bounded weak solutions to a class of cross-diffusion systems for two population species is proved. The diffusivities are assumed to depend linearly on the population densities in such a way that a certain formal gradient-flow structure holds. The main feature of these systems is that the diffusion matrix may be neither symmetric nor positive definite. The key idea of the proof is to employ the boundedness-by-entropy principle which yields at the same time the existence of global weak solutions and their boundedness. In particular, the uniform boundedness of weak solutions to the population model of Shigesada, Kawasaki, and Teramoto in several space dimensions under certain conditions on the diffusivities is shown for the first time.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The boundedness-by-entropy principle for cross-diffusion systems

A novel principle is presented which allows for the proof of bounded weak solutions to a class of physically relevant, strongly coupled parabolic systems exhibiting a formal gradient-flow structure. The main feature of these systems is that the diffusion matrix may be generally neither symmetric nor positive semi-definite. The key idea of the principle is to employ a transformation of variables, determined by the entropy density, which is defined by the gradient-flow formulation. The transformation yields at the same time a positive semi-definite diffusion matrix, suitable gradient estimates as well as lower and/or upper bounds of the solutions. These bounds are a consequence of the transformation of variables and are obtained without the use of a maximum principle. Several classes of cross-diffusion systems are identified which can be solved by this technique. The systems are formally derived from continuous-time random walks on a lattice modeling, for instance, the motion of ions, cells, or fluid particles.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A note on Aubin-Lions-Dubinskii lemmas

Strong compactness results for families of functions in seminormed nonnegative cones in the spirit of the Aubin-Lions-Dubinski\uı lemma are proven, refining some recent results in the literature. The first theorem sharpens slightly a result of Dubinski\uı (1965) for seminormed cones. The second theorem applies to piecewise constant functions in time and sharpens slightly the results of Dreher and Jüngel (2012) and Chen and Liu (2012). An application is given, which is useful in the study of porous-medium or fast-diffusion type equations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Analysis of an incompressible Navier-Stokes-Maxwell-Stefan system

The incompressible Navier-Stokes equations coupled to the Maxwell-Stefan relations for the molar fluxes are analyzed in bounded domains with no-flux boundary conditions. The system models the dynamics of a multicomponent gaseous mixture under isothermal conditions. The global-in-time existence of bounded weak solutions to the strongly coupled model and their exponential decay to the homogeneous steady state are proved. The mathematical difficulties are due to the singular Maxwell-Stefan diffusion matrix, the cross-diffusion terms, and the Navier-Stokes coupling. The key idea of the proof is the use of a new entropy functional and entropy variables, which allows for a proof of positive lower and upper bounds of the mass densities without the use of a maximum principle.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Entropy dissipative one-leg multistep time approximations of nonlinear diffusive equations

New one-leg multistep time discretizations of nonlinear evolution equations are investigated. The main features of the scheme are the preservation of the nonnegativity and the entropy-dissipation structure of the diffusive equations. The key ideas are to combine Dahlquist's G-stability theory with entropy-dissipation methods and to introduce a nonlinear transformation of variables which provides a quadratic structure in the equations. It is shown that G-stability of the one-leg scheme is sufficient to derive discrete entropy dissipation estimates. The general result is applied to a cross-diffusion system from population dynamics and a nonlinear fourth-order quantum diffusion model, for which the existence of semi-discrete weak solutions is proved. Under some assumptions on the operator of the evolution equation, the second-order convergence of solutions is shown. Moreover, some numerical experiments for the population model are presented, which underline the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Entropy-dissipative discretization of nonlinear diffusion equations and discrete Beckner inequalities

The time decay of fully discrete finite-volume approximations of porous-medium and fast-diffusion equations with Neumann or periodic boundary conditions is proved in the entropy sense. The algebraic or exponential decay rates are computed explicitly. In particular, the numerical scheme dissipates all zeroth-order entropies which are dissipated by the continuous equation. The proofs are based on novel continuous and discrete generalized Beckner inequalities. Furthermore, the exponential decay of some first-order entropies is proved in the continuous and discrete case using systematic integration by parts. Numerical experiments in one and two space dimensions illustrate the theoretical results and indicate that some restrictions on the parameters seem to be only technical.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Global existence analysis for degenerate energy-transport models for semiconductors

A class of energy-transport equations without electric field under mixed Dirichlet-Neumann boundary conditions is analyzed. The system of degenerate and strongly coupled parabolic equations for the particle density and temperature arises in semiconductor device theory. The global-in-time existence of weak nonnegative solutions is shown. The proof consists of a variable transformation and a semi-discretization in time such that the discretized system becomes elliptic and semilinear. Positive approximate solutions are obtained by Stampacchia truncation arguments and a new cut-off test function. Nonlogarithmic entropy inequalities yield gradient estimates which allow for the limit of vanishing time step sizes. Exploiting the entropy inequality, the long-time convergence of the weak solutions to the constant steady state is proved. Because of the lack of appropriate convex Sobolev inequalities to estimate the entropy dissipation, only an algebraic decay rate is obtained. Numerical experiments indicate that the decay rate is typically exponential.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Perfectly Matched Layers versus discrete transparent boundary conditions in quantum device simulations

Discrete transparent boundary conditions (DTBC) and the Perfectly Matched Layers (PML) method for the realization of open boundary conditions in quantum device simulations are compared, based on the stationary and time-dependent Schrödinger equation. The comparison includes scattering state, wave packet, and transient scattering state simulations in one and two space dimensions. The Schrödinger equation is discretized by a second-order Crank-Nicolson method in case of DTBC. For the discretization with PML, symmetric second-, fourth, and sixth-order spatial approximations as well as Crank-Nicolson and classical Runge-Kutta time-integration methods are employed. In two space dimensions, a ring-shaped quantum waveguide device is simulated in the stationary and transient regime. As an application, a simulation of the Aharonov-Bohm effect in this device is performed, showing the excitation of bound states localized in the ring region. The numerical experiments show that the results obtained from PML are comparable to those obtained using DTBC, while keeping the high numerical efficiency and flexibility as well as the ease of implementation of the former method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A finite volume scheme for a Keller-Segel model with additional cross-diffusion

A finite volume scheme for the (Patlak-) Keller-Segel model in two space dimensions with an additional cross-diffusion term in the elliptic equation for the chemical signal is analyzed. The main feature of the model is that there exists a new entropy functional yielding gradient estimates for the cell density and chemical concentration. The main features of the numerical scheme are positivity preservation, mass conservation, entropy stability, and - under additional assumptions - entropy dissipation. The existence of a discrete solution and its numerical convergence to the continuous solution is proved. Furthermore, temporal decay rates for convergence of the discrete solution to the homogeneous steady state is shown using a new discrete logarithmic Sobolev inequality. Numerical examples point out that the solutions exhibit intermediate states and that there exist nonhomogeneous stationary solutions with a finite cell density peak at the domain boundary.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Entropy-stable and entropy-dissipative approximations of a fourth-order quantum diffusion equation

Structure-preserving numerical schemes for a nonlinear parabolic fourth-order equation, modeling the electron transport in quantum semiconductors, with periodic boundary conditions are analyzed. First, a two-step backward differentiation formula (BDF) semi-discretization in time is investigated. The scheme preserves the nonnegativity of the solution, is entropy stable and dissipates a modified entropy functional. The existence of a weak semi-discrete solution and, in a particular case, its temporal second-order convergence to the continuous solution is proved. The proofs employ an algebraic relation which implies the G-stability of the two-step BDF. Second, an implicit Euler and q-step BDF discrete variational derivative method are considered. This scheme, which exploits the variational structure of the equation, dissipates the discrete Fisher information (or energy). Numerical experiments show that the discrete (relative) entropies and Fisher information decay even exponentially fast to zero.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Existence analysis for a simplified transient energy-transport model for semiconductors

A simplified transient energy-transport system for semiconductors subject to mixed Dirichlet-Neumann boundary conditions is analyzed. The model is formally derived from the non-isothermal hydrodynamic equations in a particular vanishing momentum relaxation limit. It consists of a drift-diffusion-type equation for the electron density, involving temperature gradients, a nonlinear heat equation for the electron temperature, and the Poisson equation for the electric potential. The global-in-time existence of bounded weak solutions is proved. The proof is based on the Stampacchia truncation method and a careful use of the temperature equation. Under some regularity assumptions on the gradients of the variables, the uniqueness of solutions is shown. Finally, numerical simulations for a ballistic diode in one space dimension illustrate the behavior of the solutions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Existence analysis of Maxwell-Stefan systems for multicomponent mixtures

Maxwell-Stefan systems describing the dynamics of the molar concentrations of a gas mixture with an arbitrary number of components are analyzed in a bounded domain under isobaric, isothermal conditions. The systems consist of mass balance equations and equations for the chemical potentials, depending on the relative velocities, supplemented with initial and homogeneous Neumann boundary conditions. Global-in-time existence of bounded weak solutions to the quasilinear parabolic system and their exponential decay to the homogeneous steady state are proved. The mathematical difficulties are due to the singular Maxwell-Stefan diffusion matrix, the cross-diffusion coupling, and the lack of standard maximum principles. Key ideas of the proofs are the Perron-Frobenius theory for quasi-positive matrices, entropy-dissipation methods, and a new entropy variable formulation allowing for the proof of nonnegative lower and upper bounds for the concentrations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Transient Schrödinger-Poisson Simulations of a High-Frequency Resonant Tunneling Diode Oscillator

Transient simulations of a resonant tunneling diode oscillator are presented. The semiconductor model for the diode consists of a set of time-dependent Schrödinger equations coupled to the Poisson equation for the electric potential. The one-dimensional Schrödinger equations are discretized by the finite-difference Crank-Nicolson scheme using memory-type transparent boundary conditions which model the injection of electrons from the reservoirs. This scheme is unconditionally stable and reflection-free at the boundary. An efficient recursive algorithm due to Arnold, Ehrhardt, and Sofronov is used to implement the transparent boundary conditions, enabling simulations which involve a very large number of time steps. Special care has been taken to provide a discretization of the boundary data which is completely compatible with the underlying finite-difference scheme. The transient regime between two stationary states and the self-oscillatory behavior of an oscillator circuit, containing a resonant tunneling diode, is simulated for the first time.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Cross diffusion and nonlinear diffusion preventing blow up in the Keller-Segel model

A parabolic-parabolic (Patlak-) Keller-Segel model in up to three space dimensions with nonlinear cell diffusion and an additional nonlinear cross-diffusion term is analyzed. The main feature of this model is that there exists a new entropy functional, yielding gradient estimates for the cell density and chemical concentration. For arbitrarily small cross-diffusion coefficients and for suitable exponents of the nonlinear diffusion terms, the global-in-time existence of weak solutions is proved, thus preventing finite-time blow up of the cell density. The global existence result also holds for linear and fast diffusion of the cell density in a certain parameter range in three dimensions. Furthermore, we show $L^\infty$ bounds for the solutions to the parabolic-elliptic system. Sufficient conditions leading to the asymptotic stability of the constant steady state are given for a particular choice of the nonlinear diffusion exponents. Numerical experiments in two and three space dimensions illustrate the theoretical results.

0 citations0 reviewsNo code linkedOpen access