We consider the numerical solution of a fourth‐order total variation flow problem representing surface relaxation below the roughening temperature. Based on a regularization and scaling of the nonlinear fourth‐order parabolic equation, we perform an implicit discretization in time and a C0Interior Penalty Discontinuous Galerkin (C0IPDG) discretization in space. The C0IPDG approximation can be derived from a mixed formulation involving numerical flux functions where an appropriate choice of the flux functions allows to eliminate the discrete dual variable. The fully discrete problem can be interpreted as a parameter dependent nonlinear system with the discrete time as a parameter. It is solved by a predictor corrector continuation strategy featuring an adaptive choice of the time step sizes. A documentation of numerical results is provided illustrating the performance of the C0IPDG method and the predictor corrector continuation strategy. The existence and uniqueness of a solution of the C0IPDG method will be shown in the second part of this paper.
We prove the existence and uniqueness of a solution of a C0Interior Penalty Discontinuous Galerkin (C0IPDG) method for the numerical solution of a fourth‐order total variation flow problem that has been developed in part I of the paper. The proof relies on a nonlinear version of the Lax‐Milgram Lemma. It requires to establish that the nonlinear operator associated with the C0IPDG approximation is Lipschitz continuous and strongly monotone on bounded sets of the underlying finite element space.
more » « less- NSF-PAR ID:
- 10086677
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Numerical Methods for Partial Differential Equations
- Volume:
- 35
- Issue:
- 4
- ISSN:
- 0749-159X
- Page Range / eLocation ID:
- p. 1477-1496
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
Beginning from the shallow water equations (SWEs), a nonlinear self-similar analytic solution is derived for barotropic flow over varying topography. We study conditions relevant to the ocean slope where the flow is dominated by Earth's rotation and topography. The solution is found to extend the topographic β-plume solution of Kuehl (2014) in two ways. (1) The solution is valid for intensifying jets. (2) The influence of nonlinear advection is included. The SWEs are scaled to the case of a topographically controlled jet, and then solved by introducing a similarity variable, η = cxnxyny. The nonlinear solution, valid for topographies h = h0 − αxy3, takes the form of the Lambert W-function for pseudo velocity. The linear solution, valid for topographies h = h0 − αxy−γ, takes the form of the error function for transport. Kuehl's results considered the case −1 ≤ γ < 1 which admits expanding jets, while the new result considers the case γ < −1 which admits intensifying jets and a nonlinear case with γ = −3.more » « less
-
more » « less
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 of
O (k α +h 4)u (x ,t ) ∈C α ([0,T ];C 6(Ω))0 < , which improves previous regularity assumptions that requireα < 1C 2[0,T ]u (x ,t ) ∈C 1([0,T ];C 6(Ω))O (k ) -
more » « less
Abstract Octreotide acetate, the active pharmaceutical ingredient in the long‐acting release (LAR) drug product Sandostatin®, is a cyclic octapeptide that mimics the naturally occurring somatostatin peptide hormone. Modern NMR can be a robust analytical method to identify and quantify octreotide molecules. Previous1H chemical shift assignments were mostly performed in organic solvents, and no assignments for heteronuclear13C,15N, and aromatic1H nuclei are available. Here, using state‐of‐the‐art 1D and 2D homo‐ and heteronuclear NMR experiments, octreotide was fully assigned, including water exchangeable amide protons, in aqueous buffer except for13CO and15NH of F1,15NH of C2, and15N
ζ Hζ of K5 that were not observed because of water exchange or conformational exchange. The solution NMR spectra were then directly compared with 1D1H/13C/15N solid‐state NMR (SSNMR) spectra showing the potential applicability of13C/15N SSNMR for octreotide drug product characterization. -
more » « less
Abstract The T1of a hyperpolarized site in solution is a key parameter that determines the time‐window in which its NMR signals are observable. For13C sites adjacent to protons,1H‐decoupling has been shown to increase the hyperpolarized signal resolution and SNR. Additionally, polarization transfer to protons has shown utility in increasing the sensitivity of detection. However,1H‐decoupling could lead to a change in the decay rate of a hyperpolarized13C site. Here we tested this possible effect in a case where the protons are directly bound to an sp3hyperpolarized13C site (using [1,2‐13C2]choline) and1H‐decoupling was applied continuously throughout the hyperpolarized decay measurement. We found that1H‐decoupling did not lead to any significant changes in the13C polarization decay time but did result in the expected collapse of J‐coupling and produced sharper signals. This result suggested that1H‐decoupling did not affect the decay rate of hyperpolarized sp3 13C sites. The deuterium‐substitution approach (using [1,1,2,2‐D4,1‐13C]choline) showed a dramatic prolongation of T1. Upper bounds on the T1of all investigated sites were calculated.
