skip to main content

Title: On an inverse potential problem for a fractional reaction–diffusion equation

This paper considers an inverse problem for a reaction diffusion equation from overposed final time data. Specifically, we assume that the reaction termis known but modified by a space-dependent coefficientto obtain. Thus the strength of the reaction can vary with location. The inverse problem is to recover this coefficient. Our technique is to use iterative Newton-type methods although we also use and analyse higher order schemes of Halley type. We show that such schemes are well defined and prove convergence results. Our assumption about the diffusion process is also more general in that we will extend the traditional parabolic equation paradigm to include the subdiffusion case based on non-local fractional order operators in time. The final section of the paper shows numerical reconstructions based on the above methods and compares our methodology to previous work based on the linear model withas well as to the nonlinear case. We also show the interdependence between effective reconstruction ofqand the coupling between the value of the final time of measurement and the subdiffusion parameter.

more » « less
Author(s) / Creator(s):
Publisher / Repository:
IOP Publishing
Date Published:
Journal Name:
Inverse Problems
Page Range / eLocation ID:
Article No. 065004
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Both thin (55μm) composite and thick (350μm) all active material battery porous electrodes were prepared for estimating the diffusion coefficient of Li+(DLi+)in tellurium (Te) during electrochemical lithiation. Galvanostatic intermittent titration technique (GITT), cyclic voltammetry (CV), and electrochemical impedance spectroscopy (EIS) were applied to quantify the chemical lithium solid-state diffusion coefficient within the Te active material in the electrodes. Multiple methods of GITT and EIS were assessed. For the composite Te electrodes, theDLi+was on the order of 10−11cm2s−1from both CV and GITT methods, but 10−9cm2s−1from EIS. For the thick tellurium electrodes, both GITT and EIS resulted in lithium diffusion coefficient estimates in the range of 10−11–10−12cm2s−1. The general trend across all methods that quantified the diffusion coefficient as a function of lithiation of tellurium was that theDLi+decreased rapidly when the Te material was initially lithiated. TheDLi+at the phase transition voltage plateau (∼1.7 V, vs Li/Li+, where both Te and Li2Te were expected) had the lowestDLi+,while theDLi+both before and after the plateau was generally higher. Among all the electrochemical measurements ofDLi+,the modified GITT method with modelling the relaxation region resulted in relatively low scatter in the data, provided values as a function of lithiation, and was well suited to thick electrodes with a flat discharge plateau as was the case herein.

    more » « less
  2. Abstract

    Gravitational-wave (GW) detections of binary black hole (BH) mergers have begun to sample the cosmic BH mass distribution. The evolution of single stellar cores predicts a gap in the BH mass distribution due to pair-instability supernovae (PISNe). Determining the upper and lower edges of the BH mass gap can be useful for interpreting GW detections of merging BHs. We useMESAto evolve single, nonrotating, massive helium cores with a metallicity ofZ= 10−5, until they either collapse to form a BH or explode as a PISN, without leaving a compact remnant. We calculate the boundaries of the lower BH mass gap for S-factors in the range S(300 keV) = (77,203) keV b, corresponding to the ±3σuncertainty in our high-resolution tabulated12C(α,γ)16O reaction rate probability distribution function. We extensively test temporal and spatial resolutions for resolving the theoretical peak of the BH mass spectrum across the BH mass gap. We explore the convergence with respect to convective mixing and nuclear burning, finding that significant time resolution is needed to achieve convergence. We also test adopting a minimum diffusion coefficient to help lower-resolution models reach convergence. We establish a new lower edge of the upper mass gap asMlower6014+32Mfrom the ±3σuncertainty in the12C(α,γ)16O rate. We explore the effect of a larger 3αrate on the lower edge of the upper mass gap, findingMlower6918+34M. We compare our results with BHs reported in the Gravitational-Wave Transient Catalog.

    more » « less
  3. Abstract

    While it is well known that cosmic rays (CRs) can gain energy from turbulence via second-order Fermi acceleration, how this energy transfer affects the turbulent cascade remains largely unexplored. Here, we show that damping and steepening of the compressive turbulent power spectrum are expected once the damping timetdampρv2/ĖCRECR1becomes comparable to the turbulent cascade time. Magnetohydrodynamic simulations of stirred compressive turbulence in a gas-CR fluid with diffusive CR transport show clear imprints of CR-induced damping, saturating atĖCRϵ˜, whereϵ˜is the turbulent energy input rate. In that case, almost all of the energy in large-scale motions is absorbed by CRs and does not cascade down to grid scale. Through a Hodge–Helmholtz decomposition, we confirm that purely compressive forcing can generate significant solenoidal motions, and we find preferential CR damping of the compressive component in simulations with diffusion and streaming, rendering small-scale turbulence largely solenoidal, with implications for thermal instability and proposed resonant scattering ofE≳ 300 GeV CRs by fast modes. When CR transport is streaming dominated, CRs also damp large-scale motions, with kinetic energy reduced by up to 1 order of magnitude in realisticECREgscenarios, but turbulence (with a reduced amplitude) still cascades down to small scales with the same power spectrum. Such large-scale damping implies that turbulent velocities obtained from the observed velocity dispersion may significantly underestimate turbulent forcing rates, i.e.,ϵ˜ρv3/L.

    more » « less
  4. Abstract

    The genericity of Arnold diffusion in the analytic category is an open problem. In this paper, we study this problem in the followinga prioriunstable Hamiltonian system with a time-periodic perturbationHε(p,q,I,φ,t)=h(I)+i=1n±12pi2+Vi(qi)+εH1(p,q,I,φ,t),where(p,q)Rn×Tn,(I,φ)Rd×Tdwithn,d⩾ 1,Viare Morse potentials, andɛis a small non-zero parameter. The unperturbed Hamiltonian is not necessarily convex, and the induced inner dynamics does not need to satisfy a twist condition. Using geometric methods we prove that Arnold diffusion occurs for generic analytic perturbationsH1. Indeed, the set of admissibleH1isCωdense andC3open (a fortiori,Cωopen). Our perturbative technique for the genericity is valid in theCktopology for allk∈ [3, ∞) ∪ {∞,ω}.

    more » « less
  5. Abstract

    We combine our dynamical modeling black-hole mass measurements from the Lick AGN Monitoring Project 2016 sample with measured cross-correlation time lags and line widths to recover individual scale factors,f, used in traditional reverberation-mapping analyses. We extend our sample by including prior results from Code for AGN Reverberation and Modeling of Emission Lines (caramel) studies that have utilized our methods. Aiming to improve the precision of black-hole mass estimates, as well as uncover any regularities in the behavior of the broad-line region (BLR), we search for correlations betweenfand other AGN/BLR parameters. We find (i) evidence for a correlation between the virial coefficientlog10(fmean,σ)and black-hole mass, (ii) marginal evidence for a similar correlation betweenlog10(frms,σ)and black-hole mass, (iii) marginal evidence for an anticorrelation of BLR disk thickness withlog10(fmean,FWHM)andlog10(frms,FWHM), and (iv) marginal evidence for an anticorrelation of inclination angle withlog10(fmean,FWHM),log10(frms,σ), andlog10(fmean,σ). Last, we find marginal evidence for a correlation between line-profile shape, when using the root-mean-square spectrum,log10(FWHM/σ)rms, and the virial coefficient,log10(frms,σ), and investigate how BLR properties might be related to line-profile shape usingcaramelmodels.

    more » « less