Stability of Filters for the Navier-Stokes Equation

Data assimilation methodologies are designed to incorporate noisy observations of a physical system into an underlying model in order to infer the properties of the state of the system. Filters refer to a class of data assimilation algorithms designed to update the estimation of the state in a on-line fashion, as data is acquired sequentially. For linear problems subject to Gaussian noise filtering can be performed exactly using the Kalman filter. For nonlinear systems it can be approximated in a systematic way by particle filters. However in high dimensions these particle filtering methods can break down. Hence, for the large nonlinear systems arising in applications such as weather forecasting, various ad hoc filters are used, mostly based on making Gaussian approximations. The purpose of this work is to study the properties of these ad hoc filters, working in the context of the 2D incompressible Navier-Stokes equation. By working in this infinite dimensional setting we provide an analysis which is useful for understanding high dimensional filtering, and is robust to mesh-refinement. We describe theoretical results showing that, in the small observational noise limit, the filters can be tuned to accurately track the signal itself (filter stability), provided the system is observed in a sufficiently large low dimensional space; roughly speaking this space should be large enough to contain the unstable modes of the linearized dynamics. Numerical results are given which illustrate the theory. In a simplified scenario we also derive, and study numerically, a stochastic PDE which determines filter stability in the limit of frequent observations, subject to large observational noise. The positive results herein concerning filter stability complement recent numerical studies which demonstrate that the ad hoc filters perform poorly in reproducing statistical variation about the true signal

Resolving singular forces in cavity flow: Multiscale modeling from atoms to millimeters

A multiscale approach for fluid flow is developed that retains an atomistic description in key regions. The method is applied to a classic problem where all scales contribute: The force on a moving wall bounding a fluid-filled cavity. Continuum equations predict an infinite force due to stress singularities. Following the stress over more than six decades in length in systems with characteristic scales of millimeters and milliseconds allows us to resolve the singularities and determine the force for the first time. The speedup over pure atomistic calculations is more than fourteen orders of magnitude. We find a universal dependence on the macroscopic Reynolds number, and large atomistic effects that depend on wall velocity and interactions.Comment: 4 pages,3 figure

Exploring the Dust Content of Galactic Winds with Herschel. I. NGC 4631

We present a detailed analysis of deep far-infrared observations of the nearby edge-on star-forming galaxy NGC 4631 obtained with the Herschel Space Observatory. Our PACS images at 70 and 160 um show a rich complex of filaments and chimney-like features that extends up to a projected distance of 6 kpc above the plane of the galaxy. The PACS features often match extraplanar Halpha, radio-continuum, and soft X-ray features observed in this galaxy, pointing to a tight disk-halo connection regulated by star formation. On the other hand, the morphology of the colder dust component detected on larger scale in the SPIRE 250, 350, and 500 um data matches the extraplanar H~I streams previously reported in NGC 4631 and suggests a tidal origin. The PACS 70/160 ratios are elevated in the central ~3.0 kpc region above the nucleus of this galaxy (the "superbubble"). A pixel-by-pixel analysis shows that dust in this region has a higher temperature and/or an emissivity with a steeper spectral index (beta > 2) than the dust in the disk, possibly the result of the harsher environment in the superbubble. Star formation in the disk seems energetically insufficient to lift the material out of the disk, unless it was more active in the past or the dust-to-gas ratio in the superbubble region is higher than the Galactic value. Some of the dust in the halo may also have been tidally stripped from nearby companions or lifted from the disk by galaxy interactions.Comment: Accepted for publication in The Astrophysical Journa

Towards precision medicine for pain: diagnostic biomarkers and repurposed drugs

We endeavored to identify objective blood biomarkers for pain, a subjective sensation with a biological basis, using a stepwise discovery, prioritization, validation, and testing in independent cohorts design. We studied psychiatric patients, a high risk group for co-morbid pain disorders and increased perception of pain. For discovery, we used a powerful within-subject longitudinal design. We were successful in identifying blood gene expression biomarkers that were predictive of pain state, and of future emergency department (ED) visits for pain, more so when personalized by gender and diagnosis. MFAP3, which had no prior evidence in the literature for involvement in pain, had the most robust empirical evidence from our discovery and validation steps, and was a strong predictor for pain in the independent cohorts, particularly in females and males with PTSD. Other biomarkers with best overall convergent functional evidence for involvement in pain were GNG7, CNTN1, LY9, CCDC144B, and GBP1. Some of the individual biomarkers identified are targets of existing drugs. Moreover, the biomarker gene expression signatures were used for bioinformatic drug repurposing analyses, yielding leads for possible new drug candidates such as SC-560 (an NSAID), and amoxapine (an antidepressant), as well as natural compounds such as pyridoxine (vitamin B6), cyanocobalamin (vitamin B12), and apigenin (a plant flavonoid). Our work may help mitigate the diagnostic and treatment dilemmas that have contributed to the current opioid epidemic

Apples to apples? Neural correlates of emotion regulation differences between high- and low-risk adolescents

Adolescence has been noted as a period of increased risk taking. The literature on normative neurodevelopment implicates aberrant activation of affective and regulatory regions as key to inhibitory failures. However, many of these studies have not included adolescents engaging in high rates of risky behavior, making generalizations to the most at-risk populations potentially problematic. We conducted a comparative study of nondelinquent community (n = 24, mean age = 15.8 years, 12 female) and delinquent adolescents (n = 24, mean age = 16.2 years, 12 female) who completed a cognitive control task during functional magnetic resonance imaging, where behavioral inhibition was assessed in the presence of appetitive and aversive socioaffective cues. Community adolescents showed poorer behavioral regulation to appetitive relative to aversive cues, whereas the delinquent sample showed the opposite pattern. Recruitment of the inferior frontal gyrus, medial prefrontal cortex, and tempoparietal junction differentiated community and high-risk adolescents, as delinquent adolescents showed significantly greater recruitment when inhibiting their responses in the presence of aversive cues, while the community sample showed greater recruitment when inhibiting their responses in the presence of appetitive cues. Accounting for behavioral history may be key in understanding when adolescents will have regulatory difficulties, highlighting a need for comparative research into normative and nonnormative risk-taking trajectories

Demonstration of superluminal effects in an absorptionless, non-reflective system

We present an experimental and theoretical study of a simple, passive system consisting of a birefringent, two-dimensional photonic crystal and a polarizer in series, and show that superluminal dispersive effects can arise even though no incident radiation is absorbed or reflected. We demonstrate that a vector formulation of the Kramers-Kronig dispersion relations facilitates an understanding of these counter-intuitive effects.Comment: 6 pages, 3 figures, accepted on Physical Review Letter

A new cell primo-culture method for freshwater benthic diatom communities

A new cell primo-culture method was developed for the benthic diatom community isolated from biofilm sampled in rivers. The approach comprised three steps: (1) scraping biofilm from river pebbles, (2) diatom isolation from biofilm, and (3) diatom community culture. With a view to designing a method able to stimulate the growth of diatoms, to limit the development of other microorganisms, and to maintain in culture a community similar to the original natural one, different factors were tested in step 3: cell culture medium (Chu No 10 vs Freshwater “WC” medium modified), cell culture vessel, and time of culture. The results showed that using Chu No 10 medium in an Erlenmeyer flask for cell culture was the optimal method, producing enough biomass for ecotoxicological tests as well as minimising development of other microorganisms. After 96 h of culture, communities differed from the original communities sampled in the two rivers studied. Species tolerant of eutrophic or saprobic conditions were favoured during culture. This method of diatom community culture affords the opportunity to assess, in vitro, the effects of different chemicals or effluents (water samples andindustrial effluents) on diatom communities, as well as on diatom cells, from a wide range of perspectives

A dynamical approach to the spatiotemporal aspects of the Portevin-Le Chatelier effect: Chaos,turbulence and band propagation

Experimental time series obtained from single and poly-crystals subjected to a constant strain rate tests report an intriguing dynamical crossover from a low dimensional chaotic state at medium strain rates to an infinite dimensional power law state of stress drops at high strain rates. We present results of an extensive study of all aspects of the PLC effect within the context a model that reproduces this crossover. A study of the distribution of the Lyapunov exponents as a function of strain rate shows that it changes from a small set of positive exponents in the chaotic regime to a dense set of null exponents in the scaling regime. As the latter feature is similar to the GOY shell model for turbulence, we compare our results with the GOY model. Interestingly, the null exponents in our model themselves obey a power law. The configuration of dislocations is visualized through the slow manifold analysis. This shows that while a large proportion of dislocations are in the pinned state in the chaotic regime, most of them are at the threshold of unpinning in the scaling regime. The model qualitatively reproduces the different types of deformation bands seen in experiments. At high strain rates where propagating bands are seen, the model equations are reduced to the Fisher-Kolmogorov equation for propagative fronts. This shows that the velocity of the bands varies linearly with the strain rate and inversely with the dislocation density, consistent with the known experimental results. Thus, this simple dynamical model captures the complex spatio-temporal features of the PLC effect.Comment: 17 pages, 18 figure