In this paper, the author derives an
The Cable equation is one of the most fundamental equations for modeling neuronal dynamics. In this article, we consider a high order compact finite difference numerical solution for the fractional Cable equation, which is a generalization of the classical Cable equation by taking into account the anomalous diffusion in the movement of the ions in neuronal system. The resulting finite difference scheme is unconditionally stable and converges with the convergence order of
- NSF-PAR ID:
- 10060870
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Numerical Methods for Partial Differential Equations
- Volume:
- 34
- Issue:
- 6
- ISSN:
- 0749-159X
- Page Range / eLocation ID:
- p. 2237-2266
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
‐superconvergence for the piecewise linear Ritz‐Galerkin finite element approximations for the second‐order elliptic equation equipped with Dirichlet boundary conditions. This superconvergence error estimate is established between the finite element solution and the usual Lagrange nodal point interpolation of the exact solution, and thus the superconvergence at the nodal points of each element. The result is based on a condition for the finite element partition characterized by the coefficient tensor and the usual shape functions on each element, called ‐equilateral assumption in this paper. Several examples are presented for the coefficient tensor and finite element triangulations which satisfy the conditions necessary for superconvergence. Some numerical experiments are conducted to confirm this new theory of superconvergence. -
more » « less
Abstract We consider the mapping properties of the integral operator arising in nonlocal slender body theory (SBT) for the model geometry of a straight, periodic filament. It is well known that the classical singular SBT integral operator suffers from high wavenumber instabilities, making it unsuitable for approximating the
slender body inverse problem , where the fiber velocity is prescribed and the integral operator must be inverted to find the force density along the fiber. Regularizations of the integral operator must therefore be used instead. Here, we consider two regularization methods: spectral truncation and the‐regularization of Tornberg and Shelley (2004). We compare the mapping properties of these approximations to the underlying partial differential equation (PDE) solution, which for the inverse problem is simply the Stokes Dirichlet problem with data constrained to be constant on cross sections. For the straight‐but‐periodic fiber with constant radius , we explicitly calculate the spectrum of the operator mapping fiber velocity to force for both the PDE and the approximations. We prove that the spectrum of the original SBT operator agrees closely with the PDE operator at low wavenumbers but differs at high frequencies, allowing us to define a truncated approximation with a wavenumber cutoff . For both the truncated and ‐regularized approximations, we obtain rigorous ‐based convergence to the PDE solution as : A fiber velocity with regularity gives convergence, while a fiber velocity with at least regularity yields convergence. Moreover, we determine the dependence of the ‐regularized error estimate on the regularization parameter . -
more » « less
Abstract We investigate in this work two different types of instabilities that set limits on the rotation rates of neutron (compact) stars. The first one is that caused by rotation at the Kepler frequency, at which mass shedding at the star's equator sets in. The second limit is set by instabilities driven by the growth of gravitational radiation‐reaction (GRR) driven
‐modes of order , which are moderated by shear and bulk viscosity. The calculations are performed for two relativistic models for the nuclear equation of state, DD2 and ACB4. The latter accounts for a phase transition that gives rise to the existence of so‐called mass‐twin compact stars. Our results confirm that the stable rotation periods of cold neutron stars are determined by the modes and that these modes are excited at rotation periods between 1 and 1.4 ms (20–30% above the Kepler periods of these stars). The situation is reversed in hot neutron stars where bulk viscosity damps the GRR modes, pushing the excitation period of the ‐mode instability to values below the Kepler period. For cold mass‐twin compact stars, we find that the instability sets in at rotation periods between 0.8 and 1 ms (25–30% below the Kepler period). This feature may allow one to distinguish conventional neutron stars from their possibly existing mass‐twin counterparts observationally, provided the ‐mode instability, which is expected to compete with the ‐mode instability, sets the limit on stable rotation of compact stars. -
more » « less
Abstract We consider a higher‐dimensional version of the Benjamin‐Ono (HBO) equation in the 2D setting:
, which is ‐critical, and investigate properties of solutions both analytically and numerically. For a generalized equation (fractional 2D gKdV) after deriving the Pohozaev identities, we obtain nonexistence conditions for solitary wave solutions, then prove uniform bounds in the energy space or conditional global existence, and investigate the radiation region, a specific wedge in the negative ‐direction. We then introduce our numerical approach in a general context, and apply it to obtain the ground state solution in the 2D critical HBO equation, then show that its mass is a threshold for global versus finite time existing solutions, which is typical in the focusing (mass‐)critical dispersive equations. We also observe that globally existing solutions tend to disperse completely into the radiation in this nonlocal equation. The blow‐up solutions travel in the positive ‐direction with the rescaled ground state profile while also radiating dispersive oscillations into the radiative wedge. We conclude with examples of different interactions of two solitary wave solutions, including weak and strong interactions. -
more » « less
Key points Extreme aviation is accompanied by ever‐present risks of hypobaric hypoxia and decompression sickness. Neuroprotection against those hazards is conferred through fractional inspired oxygen (
) concentrations of 60–100% (hyperoxia). Hyperoxia reduces global cerebral perfusion (gCBF), increases reactive oxygen species within the brain and leads to cell death within the hippocampus. However, an understanding of hyperoxia's effect on cortical activity and concomitant levels of cognitive performance is lacking. This limits our understanding of whether hyperoxia could lower the brain's threshold of tolerance to physiological stressors inherent to extreme aviation, such as high gravitational forces.
This study aimed to quantify the impact of hyperoxia upon global cerebral perfusion (gCBF), cognitive performance and cortical electroencephalography (EEG).
Hyperoxia evoked a rapid reduction in gCBF, yet cognitive performance and vigilance were enhanced. EEG measurements revealed enhanced alpha power, suggesting less desynchrony, within the cortical temporal regions.
Collectively, this work suggests hyperoxia‐induced brain hypoperfusion is accompanied by enhanced cognitive processing and cortical arousal.
Abstract Extreme aviators continually inspire hyperoxic gas to mitigate risk of hypoxia and decompression injury. This neuroprotection carries a physiological cost: reduced cerebral perfusion (CBF). As reduced CBF may increase vulnerability to ever‐present physiological challenges during extreme aviation, we defined the magnitude and duration of hyperoxia‐induced changes in CBF, cortical electrical activity and cognition in 30 healthy males and females. Magnetic resonance imaging with pulsed arterial spin labelling provided serial measurements of global CBF (gCBF), first during exposure to 21% inspired oxygen (
) followed by a 30‐min exposure to 100% . High‐density EEG facilitated characterization of cortical activity during assessment of cognitive performance, also measured during exposure to 21% and 100% . Acid‐base physiology was measured with arterial blood gases. We found that exposure to 100% reduced gCBF to 63% of baseline values across all participants. Cognitive performance testing at 21% was accompanied by increased theta and beta power with decreased alpha power across multiple cortical areas. During cognitive testing at 100% , alpha activity was less desynchronized within the temporal regions than at 21% . The collective hyperoxia‐induced changes in gCBF, cognitive performance and EEG were similar across observed partial pressures of arterial oxygen ( ), which ranged between 276–548 mmHg, and partial pressures of arterial carbon dioxide ( ), which ranged between 34–50 mmHg. Sex did not influence gCBF response to 100% . Our findings suggest hyperoxia‐induced reductions in gCBF evoke enhanced levels of cortical arousal and cognitive processing, similar to those occurring during a perceived threat.
