A Unified Framework for Modeling Continuum and Rarefied Gas Flows.
Ontology highlight
ABSTRACT:
SUBMITTER: Xiao H
PROVIDER: S-EPMC5638930 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:While the rheology of non-brownian suspensions in the dilute regime is well understood, their behavior in the dense limit remains mystifying. As the packing fraction of particles increases, particle motion becomes more collective, leading to a growing length scale and scaling properties in the rheology as the material approaches the jamming transition. There is no accepted microscopic description of this phenomenon. However, in recent years it has been understood that the elasticity of simple amorphous solids is governed by a critical point, the unjamming transition where the pressure vanishes, and where elastic properties display scaling and a diverging length scale. The correspondence between these two transitions is at present unclear. Here we show that for a simple model of dense flow, which we argue captures the essential physics near the jamming threshold, a formal analogy can be made between the rheology of the flow and the elasticity of simple networks. This analogy leads to a new conceptual framework to relate microscopic structure to rheology. It enables us to define and compute numerically normal modes and a density of states. We find striking similarities between the density of states in flow, and that of amorphous solids near unjamming: both display a plateau above some frequency scale ω(∗) ∼ |z(c) - z|, where z is the coordination of the network of particle in contact, z(c) = 2D where D is the spatial dimension. However, a spectacular difference appears: the density of states in flow displays a single mode at another frequency scale ω(min) ≪ ω(∗) governing the divergence of the viscosity.
| S-EPMC3324009 | biostudies-literature
Project description:Granular media take on great importance in industry and geophysics, posing a severe challenge to materials science. Their response properties elude known soft rheological models, even when the yield-stress discontinuity is blurred by vibro-fluidization. Here we propose a broad rheological scenario where average stress sums up a frictional contribution, generalizing conventional ?(I)-rheology, and a kinetic collisional term dominating at fast fluidization. Our conjecture fairly describes a wide series of experiments in a vibrofluidized vane setup, whose phenomenology includes velocity weakening, shear thinning, a discontinuous thinning transition, and gaseous shear thickening. The employed setup gives access to dynamic fluctuations, which exhibit a broad range of timescales. In the slow dense regime the frequency of cage-opening increases with stress and enhances, with respect to ?(I)-rheology, the decrease of viscosity. Diffusivity is exponential in the shear stress in both thinning and thickening regimes, with a huge growth near the transition.
| S-EPMC5141475 | biostudies-literature
Project description:Nucleosomes form the fundamental building blocks of eukaryotic chromatin, and previous attempts to understand the principles governing their genome-wide distribution have spurred much interest and debate in biology. In particular, the precise role of DNA sequence in shaping local chromatin structure has been controversial. This paper rigorously quantifies the contribution of hitherto-debated sequence features-including G+C content, 10.5?bp periodicity, and poly(dA:dT) tracts-to three distinct aspects of genome-wide nucleosome landscape: occupancy, translational positioning and rotational positioning. Our computational framework simultaneously learns nucleosome number and nucleosome-positioning energy from genome-wide nucleosome maps. In contrast to other previous studies, our model can predict both in vitro and in vivo nucleosome maps in Saccharomyces cerevisiae. We find that although G+C content is the primary determinant of MNase-derived nucleosome occupancy, MNase digestion biases may substantially influence this GC dependence. By contrast, poly(dA:dT) tracts are seen to deter nucleosome formation, regardless of the experimental method used. We further show that the 10.5?bp nucleotide periodicity facilitates rotational but not translational positioning. Applying our method to in vivo nucleosome maps demonstrates that, for a subset of genes, the regularly-spaced nucleosome arrays observed around transcription start sites can be partially recapitulated by DNA sequence alone. Finally, in vivo nucleosome occupancy derived from MNase-seq experiments around transcription termination sites can be mostly explained by the genomic sequence. Implications of these results and potential extensions of the proposed computational framework are discussed.
| S-EPMC6170202 | biostudies-literature
Project description:Mechanics plays a key role in the development of higher organisms. However, understanding this relationship is complicated by the difficulty of modeling the link between local forces generated at the subcellular level and deformations observed at the tissue and whole-embryo levels. Here we propose an approach first developed for lipid bilayers and cell membranes, in which force-generation by cytoskeletal elements enters a continuum mechanics formulation for the full system in the form of local changes in preferred curvature. This allows us to express and solve the system using only tissue strains. Locations of preferred curvature are simply related to products of gene expression. A solution, in that context, means relaxing the system's mechanical energy to yield global morphogenetic predictions that accommodate a tendency toward the local preferred curvature, without a need to explicitly model force-generation mechanisms at the molecular level. Our computational framework, which we call SPHARM-MECH, extends a 3D spherical harmonics parameterization known as SPHARM to combine this level of abstraction with a sparse shape representation. The integration of these two principles allows computer simulations to be performed in three dimensions on highly complex shapes, gene expression patterns, and mechanical constraints. We demonstrate our approach by modeling mesoderm invagination in the fruit-fly embryo, where local forces generated by the acto-myosin meshwork in the region of the future mesoderm lead to formation of a ventral tissue fold. The process is accompanied by substantial changes in cell shape and long-range cell movements. Applying SPHARM-MECH to whole-embryo live imaging data acquired with light-sheet microscopy reveals significant correlation between calculated and observed tissue movements. Our analysis predicts the observed cell shape anisotropy on the ventral side of the embryo and suggests an active mechanical role of mesoderm invagination in supporting the onset of germ-band extension.
| S-EPMC5985003 | biostudies-literature
Project description:Dense granular materials display a complicated set of flow properties, which differentiate them from ordinary fluids. Despite their ubiquity, no model has been developed that captures or predicts the complexities of granular flow, posing an obstacle in industrial and geophysical applications. Here we propose a 3D constitutive model for well-developed, dense granular flows aimed at filling this need. The key ingredient of the theory is a grain-size-dependent nonlocal rheology--inspired by efforts for emulsions--in which flow at a point is affected by the local stress as well as the flow in neighboring material. The microscopic physical basis for this approach borrows from recent principles in soft glassy rheology. The size-dependence is captured using a single material parameter, and the resulting model is able to quantitatively describe dense granular flows in an array of different geometries. Of particular importance, it passes the stringent test of capturing all aspects of the highly nontrivial flows observed in split-bottom cells--a geometry that has resisted modeling efforts for nearly a decade. A key benefit of the model is its simple-to-implement and highly predictive final form, as needed for many real-world applications.
| S-EPMC3637726 | biostudies-literature
Project description:Mathematical models are powerful tools to investigate, simulate, and evaluate potential interventions for infectious diseases dynamics. Much effort has focused on the Susceptible-Infected-Recovered (SIR)-type compartment models. These models consider host populations and measure change of each compartment. In this study, we propose an alternative patch dynamic modeling framework from pathogens' perspective. Each patch, the basic module of this modeling framework, has four standard mechanisms of pathogen population size change: birth (replication), death, inflow, and outflow. This framework naturally distinguishes between-host transmission process (inflow and outflow) and within-host infection process (replication) during the entire transmission-infection cycle. We demonstrate that the SIR-type model is actually a special cross-sectional and discretized case of our patch dynamics model in pathogens' viewpoint. In addition, this patch dynamics modeling framework is also an agent-based model from hosts' perspective by incorporating individual host's specific traits. We provide an operational standard to formulate this modular-designed patch dynamics model. Model parameterization is feasible with a wide range of sources, including genomics data, surveillance data, electronic health record, and from other emerging technologies such as multiomics. We then provide two proof-of-concept case studies to tackle some of the existing challenges of SIR-type models: sexually transmitted disease and healthcare acquired infections. This patch dynamics modeling framework not only provides theoretical explanations to known phenomena, but also generates novel insights of disease dynamics from a more holistic viewpoint. It is also able to simulate and handle more complicated scenarios across biological scales such as the current COVID-19 pandemic.
| S-EPMC7561140 | biostudies-literature
Project description:We present the Yale Stocks and Flows Database (YSTAFDB), which comprises most of the material stocks and flows (STAF) data generated at the Center for Industrial Ecology at Yale University since the early 2000s. These data describe material cycles, criticality, and recycling in terms of 62 elements and various engineering materials, e.g., steel, on spatial scales and timeframes ranging from cities to global and the 1800s to ca. 2013. YSTAFDB integrates this diverse collection of STAF data, previously scattered across various non-uniformly formatted electronic files, into a single data structure and file format. Here, we discuss this data structure as well as the usage and formatting of data records in YSTAFDB. YSTAFDB contains 100,000+ data records that are all situated in their systems contexts, with additional metadata included as available. YSTAFDB offers a comprehensive basis upon which STAF data can be accumulated, integrated, and exchanged, and thereby improves their accessibility. Therefore, YSTAFDB facilitates deeper understanding of sustainable materials use and management, which are key goals of contemporary sustainability science.
| S-EPMC6555823 | biostudies-other
Project description:Much development has been directed toward improving the performance and automation of spike sorting. This continuous development, while essential, has contributed to an over-saturation of new, incompatible tools that hinders rigorous benchmarking and complicates reproducible analysis. To address these limitations, we developed SpikeInterface, a Python framework designed to unify preexisting spike sorting technologies into a single codebase and to facilitate straightforward comparison and adoption of different approaches. With a few lines of code, researchers can reproducibly run, compare, and benchmark most modern spike sorting algorithms; pre-process, post-process, and visualize extracellular datasets; validate, curate, and export sorting outputs; and more. In this paper, we provide an overview of SpikeInterface and, with applications to real and simulated datasets, demonstrate how it can be utilized to reduce the burden of manual curation and to more comprehensively benchmark automated spike sorters.
| S-EPMC7704107 | biostudies-literature
Project description:In this paper we develop a continuum theory of clustering in ensembles of self-propelled inelastically colliding rods with applications to collective dynamics of common gliding bacteria Myxococcus Xanthus. A multiphase hydrodynamic model that couples densities of oriented and isotropic phases is described. This model is used for the analysis of an instability that leads to spontaneous formation of directionally moving dense clusters within initially dilute isotropic "gas" of myxobacteria. Numerical simulations of this model confirm the existence of stationary dense moving clusters and also elucidate the properties of their collisions. The results are shown to be in a qualitative agreement with experiments.
| S-EPMC3663047 | biostudies-literature
Project description:BackgroundGene and protein interaction experiments provide unique opportunities to study the molecular wiring of a cell. Integrating high-throughput functional genomics data with this information can help identifying networks associated with complex diseases and phenotypes.ResultsHere we introduce an integrated statistical framework to test network properties of single and multiple genesets under different interaction models. We implemented this framework as an open-source software, called Python Geneset Network Analysis (PyGNA). Our software is designed for easy integration into existing analysis pipelines and to generate high quality figures and reports. We also developed PyGNA to take advantage of multi-core systems to generate calibrated null distributions on large datasets. We then present the results of extensive benchmarking of the tests implemented in PyGNA and a use case inspired by RNA sequencing data analysis, showing how PyGNA can be easily integrated to study biological networks. PyGNA is available at http://github.com/stracquadaniolab/pygna and can be easily installed using the PyPi or Anaconda package managers, and Docker.ConclusionsWe present a tool for network-aware geneset analysis. PyGNA can either be readily used and easily integrated into existing high-performance data analysis pipelines or as a Python package to implement new tests and analyses. With the increasing availability of population-scale omic data, PyGNA provides a viable approach for large scale geneset network analysis.
| S-EPMC7579948 | biostudies-literature