Input and benchmarking data for flow simulations in discrete fracture networks.
Ontology highlight
ABSTRACT:
SUBMITTER: Fumagalli A
PROVIDER: S-EPMC6231035 | biostudies-literature | 2018 Dec
REPOSITORIES: biostudies-literature
Publications
Input and benchmarking data for flow simulations in discrete fracture networks.
Fumagalli Alessio A Keilegavlen Eirik E Scialò Stefano S
Data in brief 20181026
This article reports and describes the data related to the paper "Conforming, non-conforming and non-matching discretization couplings in discrete fracture network simulations" (Fumagalli et al., 2019). The data provided include a set of geometrical input data of Discrete Fracture Networks (DFNs) and a set of simulation results. The geometrical data describe the geometry of fracture networks of increasing complexity. These data also include the geometry of a DFN extruded from a real fracture out ...[more]
Similar Datasets
Project description:Since the seminal work of Cundall and Strack (1979), the Discrete Element Method (DEM) has now become accepted as a key tool amongst researchers exploring the fundamental behavior of granular materials. Along with a sustained increase in the number of publications documenting use of DEM in research, intensive development of new open-source and commercial DEM codes has taken place in the last decades. The credibility of these software packages depends on their capacity to replicate physical observations and to reproduce theoretical expressions. Researchers often calibrate DEM codes against laboratory data to gain confidence about their predictions, however, theoretical verifications at the macro and particle levels are often omitted or not explicitly documented or acknowledged. The validation of DEM codes against theoretical expressions is fundamental to guarantee reproducibility and generality of the software, and to avoid bias in more complex simulations. In this article, a dataset providing numerical simulation data along with input files is presented. The dataset relates to a series of theoretical validation approaches, previously documented in the literature, that were here applied to verify the open-source DEM code LAMMPS. The ability of LAMMPS to capture the macroscopic behaviour of granular packages is evaluated by shearing a face-center-cubic (FCC) array of monosized spheres. The calculation of particle translational/rotational motions and forces/torques is checked by considering a clump rolling down an inclined plane. Additionally, the stress-strain behavior of Toyoura sand under "drained" and "undrained" shearing is characterized by a series of LAMMPS outputs. The dataset collected from these simulations can be employed by users to benchmark new or existing DEM codes. Both the LAMMPS input scripts and the simulation results for all the cases are available in a public repository.
| S-EPMC10964069 | biostudies-literature
Project description:We show the error in water-limited yields simulated by crop models which is associated with spatially aggregated soil and climate input data. Crop simulations at large scales (regional, national, continental) frequently use input data of low resolution. Therefore, climate and soil data are often generated via averaging and sampling by area majority. This may bias simulated yields at large scales, varying largely across models. Thus, we evaluated the error associated with spatially aggregated soil and climate data for 14 crop models. Yields of winter wheat and silage maize were simulated under water-limited production conditions. We calculated this error from crop yields simulated at spatial resolutions from 1 to 100 km for the state of North Rhine-Westphalia, Germany. Most models showed yields biased by <15% when aggregating only soil data. The relative mean absolute error (rMAE) of most models using aggregated soil data was in the range or larger than the inter-annual or inter-model variability in yields. This error increased further when both climate and soil data were aggregated. Distinct error patterns indicate that the rMAE may be estimated from few soil variables. Illustrating the range of these aggregation effects across models, this study is a first step towards an ex-ante assessment of aggregation errors in large-scale simulations.
| S-EPMC4824533 | biostudies-literature
Project description:Long-time molecular dynamics (MD) simulations are now able to fold small proteins reversibly to their native structures [Lindorff-Larsen K, Piana S, Dror RO, Shaw DE (2011) Science 334(6055):517-520]. These results indicate that modern force fields can reproduce the energy surface near the native structure. To test how well the force fields recapitulate the other regions of the energy surface, MD trajectories for a variant of protein G are compared with data from site-resolved hydrogen exchange (HX) and other biophysical measurements. Because HX monitors the breaking of individual H-bonds, this experimental technique identifies the stability and H-bond content of excited states, thus enabling quantitative comparison with the simulations. Contrary to experimental findings of a cooperative, all-or-none unfolding process, the simulated denatured state ensemble, on average, is highly collapsed with some transient or persistent native 2° structure. The MD trajectories of this protein G variant and other small proteins exhibit excessive intramolecular H-bonding even for the most expanded conformations, suggesting that the force fields require improvements in describing H-bonding and backbone hydration. Moreover, these comparisons provide a general protocol for validating the ability of simulations to accurately capture rare structural fluctuations.
| S-EPMC4234613 | biostudies-literature
Project description:Mammalian cells adapt their functional state in response to external signals in form of ligands that bind receptors on the cell-surface. Mechanistically, this involves signal-processing through a complex network of molecular interactions that govern transcription factor activity patterns. Computer simulations of the information flow through this network could help predict cellular responses in health and disease. Here we develop a recurrent neural network framework constrained by prior knowledge of the signaling network with ligand-concentrations as input and transcription factor -activity as output. Applied to synthetic data, it predicts unseen test-data (Pearson correlation r=0.98) and the effects of gene knockouts (r=0.8). We stimulate macrophages with 59 different ligands, with and without addition of lipopolysaccharide, and collect transcriptomics data. The framework predicts this data under cross-validation (r=0.8) and knockout simulations suggest a role for RIPK1 in modulating the lipopolysaccharide response. This work demonstrates the feasibility of genome-scale simulations of intracellular signaling.
2022-05-09 | GSE202515 | GEO
Project description:Performed FloChIP for PU.1, RUNX3 and CEBPB to demonstrate feasibility of TF ChIP-seq on multiple factors.
2019-12-08 | E-MTAB-8553 | biostudies-arrayexpress
Project description:A general method for facilitating the interpretation of computer simulations of protein folding with minimally frustrated energy landscapes is detailed and applied to a designed ankyrin repeat protein (4ANK). In the method, groups of residues are assigned to foldons and these foldons are used to map the conformational space of the protein onto a set of discrete macrobasins. The free energies of the individual macrobasins are then calculated, informing practical kinetic analysis. Two simple assumptions about the universality of the rate for downhill transitions between macrobasins and the natural local connectivity between macrobasins lead to a scheme for predicting overall folding and unfolding rates, generating chevron plots under varying thermodynamic conditions, and inferring dominant kinetic folding pathways. To illustrate the approach, free energies of macrobasins were calculated from biased simulations of a non-additive structure-based model using two structurally motivated foldon definitions at the full and half ankyrin repeat resolutions. The calculated chevrons have features consistent with those measured in stopped flow chemical denaturation experiments. The dominant inferred folding pathway has an "inside-out", nucleation-propagation like character.
| S-EPMC3520928 | biostudies-literature
Project description:BackgroundHistorical, social and economic reasons can lead to major differences in the allocation of health system resources and research funding. These differences might endanger the progress in diagnostic and therapeutic approaches of socio-economic important diseases. The present study aimed to assess different benchmarking approaches that might be used to analyse these disproportions. Research in two categories was analysed for various output parameters and compared to input parameters. Germany was used as a high income model country. For the areas of cardiovascular and respiratory medicine density equalizing mapping procedures visualized major geographical differences in both input and output markers.ResultsAn imbalance in the state financial input was present with 36 cardiovascular versus 8 respiratory medicine state-financed full clinical university departments at the C4/W3 salary level. The imbalance in financial input is paralleled by an imbalance in overall quantitative output figures: The 36 cardiology chairs published 2708 articles in comparison to 453 articles published by the 8 respiratory medicine chairs in the period between 2002 and 2006. This is a ratio of 75.2 articles per cardiology chair and 56.63 articles per respiratory medicine chair. A similar trend is also present in the qualitative measures. Here, the 2708 cardiology publications were cited 48337 times (7290 times for respiratory medicine) which is an average citation of 17.85 per publication vs. 16.09 for respiratory medicine. The average number of citations per cardiology chair was 1342.69 in contrast to 911.25 citations per respiratory medicine chair. Further comparison of the contribution of the 16 different German states revealed major geographical differences concerning numbers of chairs, published items, total number of citations and average citations.ConclusionDespite similar significances of cardiovascular and respiratory diseases for the global burden of disease, large input and output imbalances have been revealed in the present study which point to a need for changes in funding policies. The present study supplies data that could be used for decision making in the field of health systems funding.
| S-EPMC2533656 | biostudies-literature
Project description:The density and configurational changes of crystal dislocations during plastic deformation influence the mechanical properties of materials. These influences have become clearest in nanoscale experiments, in terms of strength, hardness and work hardening size effects in small volumes. The mechanical characterization of a model crystal may be cast as an inverse problem of deducing the defect population characteristics (density, correlations) in small volumes from the mechanical behavior. In this work, we demonstrate how a deep residual network can be used to deduce the dislocation characteristics of a sample of interest using only its surface strain profiles at small deformations, and then statistically predict the mechanical response of size-affected samples at larger deformations. As a testbed of our approach, we utilize high-throughput discrete dislocation simulations for systems of widths that range from nano- to micro- meters. We show that the proposed deep learning model significantly outperforms a traditional machine learning model, as well as accurately produces statistical predictions of the size effects in samples of various widths. By visualizing the filters in convolutional layers and saliency maps, we find that the proposed model is able to learn the significant features of sample strain profiles.
| S-EPMC7237459 | biostudies-literature
Project description:The lumenal pH of an organelle is one of its defining characteristics and central to its biological function. Experiments have elucidated many of the key pH regulatory elements and how they vary from compartment-to-compartment, and continuum mathematical models have played an important role in understanding how these elements (proton pumps, counter-ion fluxes, membrane potential, buffering capacity, etc.) work together to achieve specific pH setpoints. While continuum models have proven successful in describing ion regulation at the cellular length scale, it is unknown if they are valid at the subcellular level where volumes are small, ion numbers may fluctuate wildly, and biochemical heterogeneity is large. Here, we create a discrete, stochastic (DS) model of vesicular acidification to answer this question. We used this simplified model to analyze pH measurements of isolated vesicles containing single proton pumps and compared these results to solutions from a continuum, ordinary differential equations (ODE)-based model. Both models predict similar parameter estimates for the mean proton pumping rate, membrane permeability, etc., but, as expected, the ODE model fails to report on the fluctuations in the system. The stochastic model predicts that pH fluctuations decrease during acidification, but noise analysis of single-vesicle data confirms our finding that the experimental noise is dominated by the fluorescent dye, and it reveals no insight into the true noise in the proton fluctuations. Finally, we again use the reduced DS model explore the acidification of large, lysosome-like vesicles to determine how stochastic elements, such as variations in proton-pump copy number and cycling between on and off states, impact the pH setpoint and fluctuations around this setpoint.
| S-EPMC6946529 | biostudies-literature
Project description:High-throughput genome-wide chromatin conformation capture assay (Hi-C) is routinely used to profile long-range genomic interactions and three-dimensional organization of genomes. A key application of Hi-C is the comparative analysis of genomic interactions across different time points, cellular conditions, or multiple stimuli. While operating characteristics of methods for Hi-C data processing such as normalization, pairwise interaction and higher-order organization detection have been relatively well studied, properties of methods for differential chromatin interaction detection are less investigated. We have recently developed FreeHi-C to enable data-driven non-parametric simulations from Hi-C experiments. Here, we extend FreeHi-C with a user/data-driven spike-in module to facilitate comparisons of differential chromatin interaction detection methods where the ground truth differential chromatin interactions are known under a wide variety of settings. We use FreeHi-C to benchmark four differential chromatin interaction detection methods, namely HiCcompare, multiHiCcompare, diffHic, and Selfish, using three comparative analysis settings with different sequencing depths and spike-in proportions. This comparison reveals distinguished performances in terms of the standard metrics such as the false discovery rate control, detection power, significance order, precision-recall curve, and receiver operating characteristic curve as well as overall genomic properties of the types of differential chromatin interactions detectable by each method. Furthermore, it highlights the lack of power for all methods in small replication settings.
| S-EPMC7855360 | biostudies-literature