Conformal Mapping on Rough Boundaries II: Applications to bi-harmonic problems

We use a conformal mapping method introduced in a companion paper to study the properties of bi-harmonic fields in the vicinity of rough boundaries. We focus our analysis on two different situations where such bi-harmonic problems are encountered: a Stokes flow near a rough wall and the stress distribution on the rough interface of a material in uni-axial tension. We perform a complete numerical solution of these two-dimensional problems for any univalued rough surfaces. We present results for sinusoidal and self-affine surface whose slope can locally reach 2.5. Beyond the numerical solution we present perturbative solutions of these problems. We show in particular that at first order in roughness amplitude, the surface stress of a material in uni-axial tension can be directly obtained from the Hilbert transform of the local slope. In case of self-affine surfaces, we show that the stress distribution presents, for large stresses, a power law tail whose exponent continuously depends on the roughness amplitude

Conformal Mapping on Rough Boundaries I: Applications to harmonic problems

The aim of this study is to analyze the properties of harmonic fields in the vicinity of rough boundaries where either a constant potential or a zero flux is imposed, while a constant field is prescribed at an infinite distance from this boundary. We introduce a conformal mapping technique that is tailored to this problem in two dimensions. An efficient algorithm is introduced to compute the conformal map for arbitrarily chosen boundaries. Harmonic fields can then simply be read from the conformal map. We discuss applications to "equivalent" smooth interfaces. We study the correlations between the topography and the field at the surface. Finally we apply the conformal map to the computation of inhomogeneous harmonic fields such as the derivation of Green function for localized flux on the surface of a rough boundary

A weakly stable algorithm for general Toeplitz systems

We show that a fast algorithm for the QR factorization of a Toeplitz or Hankel matrix A is weakly stable in the sense that R^T.R is close to A^T.A. Thus, when the algorithm is used to solve the semi-normal equations R^T.Rx = A^Tb, we obtain a weakly stable method for the solution of a nonsingular Toeplitz or Hankel linear system Ax = b. The algorithm also applies to the solution of the full-rank Toeplitz or Hankel least squares problem.Comment: 17 pages. An old Technical Report with postscript added. For further details, see http://wwwmaths.anu.edu.au/~brent/pub/pub143.htm

Adult Raphe-Specific Deletion of Lmx1b Leads to Central Serotonin Deficiency

The transcription factor Lmx1b is essential for the differentiation and survival of central serotonergic (5-HTergic) neurons during embryonic development. However, the role of Lmx1b in adult 5-HTergic neurons is unknown. We used an inducible Cre-LoxP system to selectively inactivate Lmx1b expression in the raphe nuclei of adult mice. Pet1-CreERT2 mice were generated and crossed with Lmx1bflox/flox mice to obtain Pet1-CreERT2; Lmx1bflox/flox mice (which termed as Lmx1b iCKO). After administration of tamoxifen, the level of 5-HT in the brain of Lmx1b iCKO mice was reduced to 60% of that in control mice, and the expression of tryptophan hydroxylase 2 (Tph2), serotonin transporter (Sert) and vesicular monoamine transporter 2 (Vmat2) was greatly down-regulated. On the other hand, the expression of dopamine and norepinephrine as well as aromatic L-amino acid decarboxylase (Aadc) and Pet1 was unchanged. Our results reveal that Lmx1b is required for the biosynthesis of 5-HT in adult mouse brain, and it may be involved in maintaining normal functions of central 5-HTergic neurons by regulating the expression of Tph2, Sert and Vmat2

A bootstrap method for sum-of-poles approximations

A bootstrap method is presented for finding efficient sum-of-poles approximations of causal functions. The method is based on a recursive application of the nonlinear least squares optimization scheme developed in (Alpert et al. in SIAM J. Numer. Anal. 37:1138–1164, 2000), followed by the balanced truncation method for model reduction in computational control theory as a final optimization step. The method is expected to be useful for a fairly large class of causal functions encountered in engineering and applied physics. The performance of the method and its application to computational physics are illustrated via several numerical examples

An Integrative Model for Soil Biogeochemistry and Methane Processes: I. Model Structure and Sensitivity Analysis

Environmental changes are anticipated to generate substantial impacts on carbon cycling in peatlands, affecting terrestrial-climate feedbacks. Understanding how peatland methane (CH4) fluxes respond to these changing environments is critical for predicting the magnitude of feedbacks from peatlands to global climate change. To improve predictions of CH4 fluxes in response to changes such as elevated atmospheric CO2 concentrations and warming, it is essential for Earth system models to include increased realism to simulate CH4 processes in a more mechanistic way. To address this need, we incorporated a new microbial-functional group-based CH4 module into the Energy Exascale Earth System land model (ELM) and tested it with multiple observational data sets at an ombrotrophic peatland bog in northern Minnesota. The model is able to simulate observed land surface CH4 fluxes and fundamental mechanisms contributing to these throughout the soil profile. The model reproduced the observed vertical distributions of dissolved organic carbon and acetate concentrations. The seasonality of acetoclastic and hydrogenotrophic methanogenesis—two key processes for CH4 production—and CH4 concentration along the soil profile were accurately simulated. Meanwhile, the model estimated that plant-mediated transport, diffusion, and ebullition contributed to ∼23.5%, 15.0%, and 61.5% of CH4 transport, respectively. A parameter sensitivity analysis showed that CH4 substrate and CH4 production were the most critical mechanisms regulating temporal patterns of surface CH4 fluxes both under ambient conditions and warming treatments. This knowledge will be used to improve Earth system model predictions of these high-carbon ecosystems from plot to regional scales

The translation, validity and reliability of the German version of the Fremantle Back Awareness Questionnaire

Background: The Fremantle Back Awareness Questionnaire (FreBAQ) claims to assess disrupted self-perception of the back. The aim of this study was to develop a German version of the Fre-BAQ (FreBAQ-G) and assess its test-retest reliability, its known-groups validity and its convergent validity with another purported measure of back perception. Methods: The FreBaQ-G was translated following international guidelines for the transcultural adaptation of questionnaires. Thirty-five patients with non-specific CLBP and 48 healthy participants were recruited. Assessor one administered the FreBAQ-G to each patient with CLBP on two separate days to quantify intra-observer reliability. Assessor two administered the FreBaQ-G to each patient on day 1. The scores were compared to those obtained by assessor one on day 1 to assess inter-observer reliability. Known-groups validity was quantified by comparing the FreBAQ-G score between patients and healthy controls. To assess convergent validity, patient\u27s FreBAQ-G scores were correlated to their two-point discrimination (TPD) scores. Results: Intra- and Inter-observer reliability were both moderate with ICC3.1 = 0.88 (95%CI: 0.77 to 0.94) and 0.89 (95%CI: 0.79 to 0.94), respectively. Intra- and inter-observer limits of agreement (LoA) were 6.2 (95%CI: 5.0±8.1) and 6.0 (4.8±7.8), respectively. The adjusted mean difference between patients and controls was 5.4 (95%CI: 3.0 to 7.8, p\u3c0.01). Patient\u27s FreBAQ-G scores were not associated with TPD thresholds (Pearson\u27s r = -0.05, p = 0.79). Conclusions: The FreBAQ-G demonstrated a degree of reliability and known-groups validity. Interpretation of patient level data should be performed with caution because the LoA were substantial. It did not demonstrate convergent validity against TPD. Floor effects of some items of the FreBAQ-G may have influenced the validity and reliability results. The clinimetric properties of the FreBAQ-G require further investigation as a simple measure of disrupted self-perception of the back before firm recommendations on its use can be made

Order-of-magnitude speedup for steady states and traveling waves via Stokes preconditioning in Channelflow and Openpipeflow

Steady states and traveling waves play a fundamental role in understanding hydrodynamic problems. Even when unstable, these states provide the bifurcation-theoretic explanation for the origin of the observed states. In turbulent wall-bounded shear flows, these states have been hypothesized to be saddle points organizing the trajectories within a chaotic attractor. These states must be computed with Newton's method or one of its generalizations, since time-integration cannot converge to unstable equilibria. The bottleneck is the solution of linear systems involving the Jacobian of the Navier-Stokes or Boussinesq equations. Originally such computations were carried out by constructing and directly inverting the Jacobian, but this is unfeasible for the matrices arising from three-dimensional hydrodynamic configurations in large domains. A popular method is to seek states that are invariant under numerical time integration. Surprisingly, equilibria may also be found by seeking flows that are invariant under a single very large Backwards-Euler Forwards-Euler timestep. We show that this method, called Stokes preconditioning, is 10 to 50 times faster at computing steady states in plane Couette flow and traveling waves in pipe flow. Moreover, it can be carried out using Channelflow (by Gibson) and Openpipeflow (by Willis) without any changes to these popular spectral codes. We explain the convergence rate as a function of the integration period and Reynolds number by computing the full spectra of the operators corresponding to the Jacobians of both methods.Comment: in Computational Modelling of Bifurcations and Instabilities in Fluid Dynamics, ed. Alexander Gelfgat (Springer, 2018

On a global differential geometric approach to the rational mechanics of deformable media

In the past the rational mechanics of deformable media was largely concerned with materials governed by linear constitutive equations. In recent years, the theory has expanded considerably towards covering materials for which the constitutive equations are inherently nonlinear, and/or whose mechanical properties resemble in some respects those of a fluid and in others those of a solid. In the present article we formulate a satisfactory global mathematical theory of moving deformable media, which includes all these aspects