An official website of the United States government
Abstract We begin with a treatment of the Caputo time‐fractional diffusion equation, by using the Laplace transform, to obtain a Volterra integro‐differential equation. We derive and utilize a numerical scheme that is derived in parallel to the L1‐method for the time variable and a standard fourth‐order approximation in the spatial variable. The main method derived in this article has a rate of convergence ofO(kα + h4)foru(x,t) ∈ Cα([0,T];C6(Ω)),0 < α < 1, which improves previous regularity assumptions that requireC2[0,T]regularity in the time variable. We also present a novel alternative method for a first‐order approximation in time, under a regularity assumption ofu(x,t) ∈ C1([0,T];C6(Ω)), while exhibiting order of convergence slightly more thanO(k)in time. This allows for a much wider class of functions to be analyzed which was previously not possible under the L1‐method. We present numerical examples demonstrating these results and discuss future improvements and implications by using these techniques.
more »
« less
The gradient descent method from the perspective of fractional calculus
Motivated by gradient methods in optimization theory, we give methods based onψ‐fractional derivatives of orderαin order to solve unconstrained optimization problems. The convergence of these methods is analyzed in detail. This paper also presents an Adams–Bashforth–Moulton (ABM) method for the estimation of solutions to equations involvingψ‐fractional derivatives. Numerical examples using the ABM method show that the fractional orderαand weightψare tunable parameters, which can be helpful for improving the performance of gradient descent methods.
more »
« less
- Award ID(s):
- 2027976
- PAR ID:
- 10452785
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Mathematical Methods in the Applied Sciences
- Volume:
- 44
- Issue:
- 7
- ISSN:
- 0170-4214
- Page Range / eLocation ID:
- p. 5520-5547
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT Ternary block copolymer (BCP)‐homopolymer (HP) blends offer a simple method for tuning nanostructure sizes to meet application‐specific demands. Comprehensive dissipative particle dynamic (DPD) simulations were performed to study the impact of polymer interactions, molecular weight, and HP volume fraction (φHP) on symmetric ternary blend morphological stability and domain spacing. DPD reproduces key features of the experimental phase diagram, including lamellar domain swelling with increasingφHP, the formation of an asymmetric bicontinuous microemulsion at a critical HP concentration , and macrophase separation with further HP addition. Simulation results matched experimental values for and lamellar swelling as a function of HP to BCP chain length ratio,α = NHP/NBCP. Structural analysis of blends with fixedφHPbut varyingαconfirmed that ternary blends follow the wet/dry brush model of domain swelling with the miscibility of HPs and BCPs depending onα. Longer HPs concentrate in the center of domains, boosting their swelling efficiencies compared to shorter chains. These results advance our understanding of BCP‐HP blend phase behavior and demonstrate the value of DPD for studying polymeric blends. © 2019 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys.2019,57, 794–803more » « less
-
In the present work, we explored the dynamics of single kinks, kink–anti-kink pairs and bound states in the prototypical fractional Klein–Gordon example of the sine-Gordon equation. In particular, we modified the order β of the temporal derivative to that of a Caputo fractional type and found that, for 1<β<2, this imposes a dissipative dynamical behavior on the coherent structures. We also examined the variation of a fractional Riesz order α on the spatial derivative. Here, depending on whether this order was below or above the harmonic value α=2, we found, respectively, monotonically attracting kinks, or non-monotonic and potentially attracting or repelling kinks, with a saddle equilibrium separating the two. Finally, we also explored the interplay of the two derivatives, when both Caputo temporal and Riesz spatial derivatives are involved.more » « less
-
Abstract Shifting baselines in the Arctic atmosphere‐sea ice‐ocean system have significant potential to alter biogeochemical cycling and ecosystem dynamics. In particular, the impact of increased open water duration on lower trophic level productivity and biological CO2sequestration is poorly understood. Using high‐resolution observations of surface seawater dissolved O2/Ar andpCO2collected in the Pacific Arctic in October 2011 and 2012, we evaluate spatial variability in biological metabolic status (autotrophy vs heterotrophy) as constrained by O2/Ar saturation (∆O2/Ar) as well as the relationship between net biological production and the sea‐air gradient ofpCO2(∆pCO2). We find a robust relationship between∆pCO2and∆O2/Ar(correlation coefficient of −0.74 and −0.61 for 2011 and 2012, respectively), which suggests that biological production in the late open water season is an important determinant of the air‐sea CO2gradient at a timeframe of maximal ocean uptake for CO2in this region. Patchiness in biological production as indicated by∆O2/Arsuggests spatially variable nutrient supply mechanisms supporting late season growth amidst a generally strongly stratified and nutrient‐limited condition.more » « less
-
Abstract Potassium (K) informs on the radiogenic heat production, atmospheric composition, and volatile element depletion of the Earth and other planetary systems. Constraints on the abundance of K in the Earth, Moon, and other rocky bodies have historically hinged on K/U values measured in planetary materials, particularly comparisons of the continental crust and mid‐ocean ridge basalts (MORBs), for developing compositional models of the bulk silicate Earth (BSE). However, a consensus on the most representative K/U value for global MORB remains elusive despite numerous studies. Here, we statistically analyze a critical compilation of MORB data to determine the K/U value of the MORB source. Covariations in the log‐normal abundances of K and U establish that K is 3–7 times less incompatible than U during melting and/or crystallization processes, enabling inverse modeling to infer the K/U of the MORB source region. These comprehensive data have a mean K/U for global MORB = 13,900 ± 200 (2σm;n = 4,646), and define a MORB source region with a K/U between 14,000 and 15,500, depending on the modeled melting regime. However, this range represents strictly a lower limit due to the undefined role of fractional crystallization in these samples and challenges preserving the signatures of depleted components in the MORB mantle source. This MORB source model, when combined with recent metadata analyses of ocean island basalt (OIB) and continental crust, suggests that the BSE has a K/U value >12,100 and contains >260 × 10−6 kg/kg K, resulting in a global production of∼3.5 TW of radiogenic heat today and 1.5 × 1017 kg of40Ar over the lifetime of the planet.more » « less
