Search Results
Auditory cortex modelled as a dynamical network of oscillators: Understanding event-related fields and their adaptation
Adaptation, the reduction of neuronal responses by repetitive stimulation, is a ubiquitous feature of auditory cortex (AC). It is not clear what causes adaptation, but short-term synaptic depression (STSD) is a potential candidate for the underlying mechanism. We examined this hypothesis via a computational model based on AC anatomy, which includes serially connected core, belt, and parabelt areas. The model replicates the event-related field (ERF) of the magnetoencephalogram as well as ERF adaptation. The model dynamics are described by excitatory and inhibitory state variables of cell populations, with the excitatory connections modulated by STSD. We analysed the system dynamics by linearizing the firing rates and solving the STSD equation using time-scale separation. This allows for characterization of AC dynamics as a superposition of damped harmonic oscillators, so-called normal modes. We show that repetition suppression of the N1m is due to a mixture of causes, with stimulus repetition modifying both the amplitudes and the frequencies of the normal modes. In this view, adaptation results from a complete reorganization of AC dynamics rather than a reduction of activity in discrete sources. Further, both the network structure and the balance between excitation and inhibition contribute significantly to the rate with which AC recovers from adaptation. This lifetime of adaptation is longer in the belt and parabelt than in the core area, despite the time constants of STSD being spatially constant. Finally, we critically evaluate the use of a single exponential function to describe recovery from adaptation.
Weak solutions and weak-strong uniqueness for a thermodynamically consistent phase-field model
In this paper we prove the existence of weak solutions for a thermodynamically consistent phase-field model introduced in [26] in two and three dimensions of space. We use a notion of solution inspired by [18], where the pointwise internal energy balance is replaced by the total energy inequality complemented with a weak form of the entropy inequality. Moreover, we prove existence of local-in-time strong solutions and, finally, we show weak-strong uniqueness of solutions, meaning that every weak solution coincides with a local strong solution emanating from the same initial data, as long as the latter exists.
Periodic Lp estimates by R-boundedness: Applications to the Navier--Stokes equations
General evolution equations in Banach spaces are investigated. Based on an operator-valued version of de Leeuw's transference principle, time-periodic Lp estimates of maximal regularity type are established from R-bounds of the family of solution operators (R-solvers) to the corresponding resolvent problems. With this method, existence of time-periodic solutions to the Navier--Stokes equations is shown for two configurations: in a periodically moving bounded domain and in an exterior domain, subject to prescribed time-periodic forcing and boundary data.
Maximal dissipative solutions for incompressible fluid dynamics
We introduce the new concept of maximal dissipative solutions for the Navier--Stokes and Euler equations and show that these solutions exist and the solution set is closed and convex. The concept of maximal dissipative solutions coincides with the concept of weak solutions as long as the weak solutions inherits enough regularity to be unique. A maximal dissipative solution is defined as the minimizer of a convex functional and we argue that this definition bears several advantages.
Weak-strong uniqueness for energy-reaction-diffusion systems
We establish weak-strong uniqueness and stability properties of renormalised solutions to a class of energy-reaction-diffusion systems, which genuinely feature cross-diffusion effects. The systems considered are motivated by thermodynamically consistent models, and their formal entropy structure allows us to use as a key tool a suitably adjusted relative entropy method. Weak-strong uniqueness is obtained for general entropy-dissipating reactions without growth restrictions, and certain models with a non-integrable diffusive flux. The results also apply to a class of (isoenergetic) reaction-cross-diffusion systems.
Fast reaction limits via $Gamma$-convergence of the flux rate functional
We study the convergence of a sequence of evolution equations for measures supported on the nodes of a graph. The evolution equations themselves can be interpreted as the forward Kolmogorov equations of Markov jump processes, or equivalently as the equations for the concentrations in a network of linear reactions. The jump rates or reaction rates are divided in two classes; `slow' rates are constant, and `fast' rates are scaled as 1/∈, and we prove the convergence in the fast-reaction limit ∈ → 0. We establish a Γ-convergence result for the rate functional in terms of both the concentration at each node and the flux over each edge (the level-2.5 rate function). The limiting system is again described by a functional, and characterizes both fast and slow fluxes in the system. This method of proof has three advantages. First, no condition of detailed balance is required. Secondly, the formulation in terms of concentration and flux leads to a short and simple proof of the Γ-convergence; the price to pay is a more involved compactness proof. Finally, the method of proof deals with approximate solutions, for which the functional is not zero but small, without any changes.
On the existence of weak solutions in the context of multidimensional incompressible fluid dynamics
We define the concept of energy-variational solutions for the Navier--Stokes and Euler equations. This concept is shown to be equivalent to weak solutions with energy conservation. Via a standard Galerkin discretization, we prove the existence of energy-variational solutions and thus weak solutions in any space dimension for the Navier--Stokes equations. In the limit of vanishing viscosity the same assertions are deduced for the incompressible Euler system. Via the selection criterion of maximal dissipation we deduce well-posedness for these equations.
Numerical analysis for nematic electrolytes
We consider a system of nonlinear PDEs modeling nematic electrolytes, and construct a dissipative solution with the help of its implementable, structure-inheriting space-time discretization. Computational studies are performed to study the mutual effects of electric, elastic, and viscous effects onto the molecules in a nematic electrolyte.
High order discretization methods for spatial-dependent epidemic models
In this paper, an SIR model with spatial dependence is studied and results regarding its stability and numerical approximation are presented. We consider a generalization of the original Kermack and McKendrick model in which the size of the populations differs in space. The use of local spatial dependence yields a system of integro-differential equations. The uniqueness and qualitative properties of the continuous model are analyzed. Furthermore, different choices of spatial and temporal discretizations are employed, and step-size restrictions for population conservation, positivity, and monotonicity preservation of the discrete model are investigated. We provide sufficient conditions under which high order numerical schemes preserve the discrete properties of the model. Computational experiments verify the convergence and accuracy of the numerical methods.
Analysis of a quasi-variational contact problem arising in thermoelasticity
We formulate and study two mathematical models of a thermoforming process involving a membrane and a mould as implicit obstacle problems. In particular, the membrane-mould coupling is determined by the thermal displacement of the mould that depends in turn on the membrane through the contact region. The two models considered are a stationary (or elliptic) model and an evolutionary (or quasistatic) one. For the first model, we prove the existence of weak solutions by solving an elliptic quasi-variational inequality coupled to elliptic equations. By exploring the fine properties of the variation of the contact set under non-degenerate data, we give sufficient conditions for the existence of regular solutions, and under certain contraction conditions, also a uniqueness result. We apply these results to a series of semi-discretised problems that arise as approximations of regular solutions for the evolutionary or quasistatic problem. Here, under certain conditions, we are able to prove existence for the evolutionary problem and for a special case, also the uniqueness of time-dependent solutions.