2,388 research outputs found

    Solving Nonlinear Parabolic Equations by a Strongly Implicit Finite-Difference Scheme

    Full text link
    We discuss the numerical solution of nonlinear parabolic partial differential equations, exhibiting finite speed of propagation, via a strongly implicit finite-difference scheme with formal truncation error O[(Δx)2+(Δt)2]\mathcal{O}\left[(\Delta x)^2 + (\Delta t)^2 \right]. Our application of interest is the spreading of viscous gravity currents in the study of which these type of differential equations arise. Viscous gravity currents are low Reynolds number (viscous forces dominate inertial forces) flow phenomena in which a dense, viscous fluid displaces a lighter (usually immiscible) fluid. The fluids may be confined by the sidewalls of a channel or propagate in an unconfined two-dimensional (or axisymmetric three-dimensional) geometry. Under the lubrication approximation, the mathematical description of the spreading of these fluids reduces to solving the so-called thin-film equation for the current's shape h(x,t)h(x,t). To solve such nonlinear parabolic equations we propose a finite-difference scheme based on the Crank--Nicolson idea. We implement the scheme for problems involving a single spatial coordinate (i.e., two-dimensional, axisymmetric or spherically-symmetric three-dimensional currents) on an equispaced but staggered grid. We benchmark the scheme against analytical solutions and highlight its strong numerical stability by specifically considering the spreading of non-Newtonian power-law fluids in a variable-width confined channel-like geometry (a "Hele-Shaw cell") subject to a given mass conservation/balance constraint. We show that this constraint can be implemented by re-expressing it as nonlinear flux boundary conditions on the domain's endpoints. Then, we show numerically that the scheme achieves its full second-order accuracy in space and time. We also highlight through numerical simulations how the proposed scheme accurately respects the mass conservation/balance constraint.Comment: 36 pages, 9 figures, Springer book class; v2 includes improvements and corrections; to appear as a contribution in "Applied Wave Mathematics II

    Moderate exercise and chronic stress produce counteractive effects on different areas of the brain by acting through various neurotransmitter receptor subtypes: A hypothesis

    Get PDF
    BACKGROUND: Regular, "moderate", physical exercise is an established non-pharmacological form of treatment for depressive disorders. Brain lateralization has a significant role in the progress of depression. External stimuli such as various stressors or exercise influence the higher functions of the brain (cognition and affect). These effects often do not follow a linear course. Therefore, nonlinear dynamics seem best suited for modeling many of the phenomena, and putative global pathways in the brain, attributable to such external influences. HYPOTHESIS: The general hypothesis presented here considers only the nonlinear aspects of the effects produced by "moderate" exercise and "chronic" stressors, but does not preclude the possibility of linear responses. In reality, both linear and nonlinear mechanisms may be involved in the final outcomes. The well-known neurotransmitters serotonin (5-HT), dopamine (D) and norepinephrine (NE) all have various receptor subtypes. The article hypothesizes that 'Stress' increases the activity/concentration of some particular subtypes of receptors (designated nt(s)) for each of the known (and unknown) neurotransmitters in the right anterior (RA) and left posterior (LP) regions (cortical and subcortical) of the brain, and has the converse effects on a different set of receptor subtypes (designated nt(h)). In contrast, 'Exercise' increases nt(h )activity/concentration and/or reduces nt(s )activity/concentration in the LA and RP areas of the brain. These effects may be initiated by the activation of Brain Derived Neurotrophic Factor (BDNF) (among others) in exercise and its suppression in stress. CONCLUSION: On the basis of this hypothesis, a better understanding of brain neurodynamics might be achieved by considering the oscillations caused by single neurotransmitters acting on their different receptor subtypes, and the temporal pattern of recruitment of these subtypes. Further, appropriately designed and planned experiments will not only corroborate such theoretical models, but also shed more light on the underlying brain dynamics

    Network adaptation improves temporal representation of naturalistic stimuli in drosophila eye: II Mechanisms

    Get PDF
    Retinal networks must adapt constantly to best present the ever changing visual world to the brain. Here we test the hypothesis that adaptation is a result of different mechanisms at several synaptic connections within the network. In a companion paper (Part I), we showed that adaptation in the photoreceptors (R1-R6) and large monopolar cells (LMC) of the Drosophila eye improves sensitivity to under-represented signals in seconds by enhancing both the amplitude and frequency distribution of LMCs' voltage responses to repeated naturalistic contrast series. In this paper, we show that such adaptation needs both the light-mediated conductance and feedback-mediated synaptic conductance. A faulty feedforward pathway in histamine receptor mutant flies speeds up the LMC output, mimicking extreme light adaptation. A faulty feedback pathway from L2 LMCs to photoreceptors slows down the LMC output, mimicking dark adaptation. These results underline the importance of network adaptation for efficient coding, and as a mechanism for selectively regulating the size and speed of signals in neurons. We suggest that concert action of many different mechanisms and neural connections are responsible for adaptation to visual stimuli. Further, our results demonstrate the need for detailed circuit reconstructions like that of the Drosophila lamina, to understand how networks process information

    Binary orbits as the driver of γ-ray emission and mass ejection in classical novae

    Get PDF
    Classical novae are the most common astrophysical thermonuclear explosions, occurring on the surfaces of white dwarf stars accreting gas from companions in binary star systems. Novae typically expel �10,000 solar masses of material at velocities exceeding 1,000 km/s. However, the mechanism of mass ejection in novae is poorly understood, and could be dominated by the impulsive flash of the thermonuclear runaway, prolonged optically thick winds, or binary interaction with the nova envelope. Classical novae are now routinely detected in GeV gamma-rays, suggesting that relativistic particles are accelerated by strong shocks in nova ejecta. Here we present high-resolution imaging of the gamma-ray-emitting nova V959 Mon at radio wavelengths, showing that its ejecta were shaped by binary motion: some gas was expelled rapidly along the poles as a wind from the white dwarf, while denser material drifted out along the equatorial plane, propelled by orbital motion. At the interface between the equatorial and polar regions, we observe synchrotron emission indicative of shocks and relativistic particle acceleration, thereby pinpointing the location of gamma-ray production. Binary shaping of the nova ejecta and associated internal shocks are expected to be widespread among novae, explaining why many novae are gamma-ray emitters

    Rare germline variants in DNA repair genes and the angiogenesis pathway predispose prostate cancer patients to develop metastatic disease

    Get PDF
    Background Prostate cancer (PrCa) demonstrates a heterogeneous clinical presentation ranging from largely indolent to lethal. We sought to identify a signature of rare inherited variants that distinguishes between these two extreme phenotypes. Methods We sequenced germline whole exomes from 139 aggressive (metastatic, age of diagnosis < 60) and 141 non-aggressive (low clinical grade, age of diagnosis ≥60) PrCa cases. We conducted rare variant association analyses at gene and gene set levels using SKAT and Bayesian risk index techniques. GO term enrichment analysis was performed for genes with the highest differential burden of rare disruptive variants. Results Protein truncating variants (PTVs) in specific DNA repair genes were significantly overrepresented among patients with the aggressive phenotype, with BRCA2, ATM and NBN the most frequently mutated genes. Differential burden of rare variants was identified between metastatic and non-aggressive cases for several genes implicated in angiogenesis, conferring both deleterious and protective effects. Conclusions Inherited PTVs in several DNA repair genes distinguish aggressive from non-aggressive PrCa cases. Furthermore, inherited variants in genes with roles in angiogenesis may be potential predictors for risk of metastases. If validated in a larger dataset, these findings have potential for future clinical application
    corecore