Hybrid finite element and Brownian dynamics method for diffusion-controlled reactions.
Ontology highlight
ABSTRACT:
SUBMITTER: Bauler P
PROVIDER: S-EPMC3350536 | biostudies-other | 2012 Apr
REPOSITORIES: biostudies-other
Similar Datasets
Project description:Reaction-diffusion systems are used to represent many biological and physical phenomena. They model the random motion of particles (diffusion) and interactions between them (reactions). Such systems can be modelled at multiple scales with varying degrees of accuracy and computational efficiency. When representing genuinely multiscale phenomena, fine-scale models can be prohibitively expensive, whereas coarser models, although cheaper, often lack sufficient detail to accurately represent the phenomenon at hand. Spatial hybrid methods couple two or more of these representations in order to improve efficiency without compromising accuracy. In this paper, we present a novel spatial hybrid method, which we call the auxiliary region method (ARM), which couples PDE- and Brownian-based representations of reaction-diffusion systems. Numerical PDE solutions on one side of an interface are coupled to Brownian-based dynamics on the other side using compartment-based 'auxiliary regions'. We demonstrate that the hybrid method is able to simulate reaction-diffusion dynamics for a number of different test problems with high accuracy. Furthermore, we undertake error analysis on the ARM which demonstrates that it is robust to changes in the free parameters in the model, where previous coupling algorithms are not. In particular, we envisage that the method will be applicable for a wide range of spatial multi-scales problems including filopodial dynamics, intracellular signalling, embryogenesis and travelling wave phenomena.
| S-EPMC6124063 | biostudies-literature
Project description:To relate the subcellular molecular events to organ level physiology in heart, we have developed a three-dimensional finite-element-based simulation program incorporating the cellular mechanisms of excitation-contraction coupling and its propagation, and simulated the fluid-structure interaction involved in the contraction and relaxation of the human left ventricle. The FitzHugh-Nagumo model and four-state model representing the cross-bridge kinetics were adopted for cellular model. Both ventricular wall and blood in the cavity were modeled by finite element mesh. An arbitrary Lagrangian Eulerian finite element method with automatic mesh updating has been formulated for large domain changes, and a strong coupling strategy has been taken. Using electrical analog of pulmonary circulation and left atrium as a preload and the windkessel model as an afterload, dynamics of ventricular filling as well as ejection was simulated. We successfully reproduced the biphasic filling flow consisting of early rapid filling and atrial contraction similar to that reported in clinical observation. Furthermore, fluid-structure analysis enabled us to analyze the wave propagation velocity of filling flow. This simulator can be a powerful tool for establishing a link between molecular abnormality and the clinical disorder at the macroscopic level.
| S-EPMC1304609 | biostudies-literature
Project description:The interiors of all living cells are highly crowded with macromolecules, which differs considerably the thermodynamics and kinetics of biological reactions between in vivo and in vitro. For example, the diffusion of green fluorescent protein (GFP) in E. coli is ~10-fold slower than in dilute conditions. In this study, we performed Brownian dynamics (BD) simulations of rigid macromolecules in a crowded environment mimicking the cytosol of E. coli to study the motions of macromolecules. The simulation systems contained 35 70S ribosomes, 750 glycolytic enzymes, 75 GFPs, and 392 tRNAs in a 100 nm × 100 nm × 100 nm simulation box, where the macromolecules were represented by rigid-objects of one bead per amino acid or four beads per nucleotide models. Diffusion tensors of these molecules in dilute solutions were estimated by using a hydrodynamic theory to take into account the diffusion anisotropy of arbitrary shaped objects in the BD simulations. BD simulations of the system where each macromolecule is represented by its Stokes radius were also performed for comparison. Excluded volume effects greatly reduce the mobility of molecules in crowded environments for both molecular-shaped and equivalent sphere systems. Additionally, there were no significant differences in the reduction of diffusivity over the entire range of molecular size between two systems. However, the reduction in diffusion of GFP in these systems was still 4-5 times larger than for the in vivo experiment. We will discuss other plausible factors that might cause the large reduction in diffusion in vivo.
| S-EPMC4295211 | biostudies-literature
Project description:The mechanics of biological fluids is an important topic in biomechanics, often requiring the use of computational tools to analyze problems with realistic geometries and material properties. This study describes the formulation and implementation of a finite element framework for computational fluid dynamics (CFD) in FEBio, a free software designed to meet the computational needs of the biomechanics and biophysics communities. This formulation models nearly incompressible flow with a compressible isothermal formulation that uses a physically realistic value for the fluid bulk modulus. It employs fluid velocity and dilatation as essential variables: The virtual work integral enforces the balance of linear momentum and the kinematic constraint between fluid velocity and dilatation, while fluid density varies with dilatation as prescribed by the axiom of mass balance. Using this approach, equal-order interpolations may be used for both essential variables over each element, contrary to traditional mixed formulations that must explicitly satisfy the inf-sup condition. The formulation accommodates Newtonian and non-Newtonian viscous responses as well as inviscid fluids. The efficiency of numerical solutions is enhanced using Broyden's quasi-Newton method. The results of finite element simulations were verified using well-documented benchmark problems as well as comparisons with other free and commercial codes. These analyses demonstrated that the novel formulation introduced in FEBio could successfully reproduce the results of other codes. The analogy between this CFD formulation and standard finite element formulations for solid mechanics makes it suitable for future extension to fluid-structure interactions (FSIs).
| S-EPMC5816258 | biostudies-literature
Project description:This article describes a numerical solution of the steady-state Poisson-Boltzmann-Smoluchowski (PBS) and Poisson-Nernst-Planck (PNP) equations to study diffusion in biomolecular systems. Specifically, finite element methods have been developed to calculate electrostatic interactions and ligand binding rate constants for large biomolecules. The resulting software has been validated and applied to the wild-type and several mutated avian influenza neurominidase crystal structures. The calculated rates show very good agreement with recent experimental studies. Furthermore, these finite element methods require significantly fewer computational resources than existing particle-based Brownian dynamics methods and are robust for complicated geometries. The key finding of biological importance is that the electrostatic steering plays the important role in the drug binding process of the neurominidase.
| S-EPMC3107071 | biostudies-literature
Project description:The lack of a universal simulation method for triboelectric nanogenerator (TENG) makes the device design and optimization difficult before experiment, which protracts the research and development process and hinders the landing of practical TENG applications. The existing electrostatic induction models for TENGs have limitations in simulating TENGs with complex geometries and their dynamic behaviors under practical movements due to the topology change issues. Here, a dynamic finite element method (FEM) model is proposed. The introduction of air buffer layers and the moving mesh method eliminates the topology change issues during practical movement and allows simulation of dynamic and time-varying behaviors of TENGs with complex 2D/3D geometries. Systematic investigations are carried out to optimize the air buffer thickness and mesh densities, and the optimized results show excellent consistency with the experimental data and results based on other existing methods. It also shows that a 3D disk-type rotating TENG can be simulated using the model, clearly demonstrating the capability and superiority of the dynamic FEM model. Moreover, the dynamic FEM model is used to optimize the shape of the tribo-material, which is used as a preliminary example to demonstrate the possibility of designing a TENG-based sensor.
| S-EPMC7506595 | biostudies-literature
Project description:We report Brownian dynamics simulations of tracer diffusion in regularly crosslinked polymer networks in order to elucidate the transport of a tracer particle in polymer networks. The average mesh size of homogeneous polymer networks is varied by assuming different degrees of crosslinking or swelling, and the size of a tracer particle is comparable to the average mesh size. Simulation results show subdiffusion of a tracer particle at intermediate time scales and normal diffusion at long times. In particular, the duration of subdiffusion is significantly prolonged as the average mesh size decreases with increasing degree of crosslinking, for which long-time diffusion occurs via the hopping processes of a tracer particle after undergoing rattling motions within a cage of the network mesh for an extended period of time. On the other hand, the cage dynamics and hopping process are less pronounced as the mesh size decreases with increasing polymer volume fractions. The interpretation is provided in terms of fluctuations in network mesh size: at higher polymer volume fractions, the network fluctuations are large enough to allow for collective, structural changes of network meshes, so that a tracer particle can escape from the cage, whereas, at lower volume fractions, the fluctuations are so small that a tracer particle remains trapped within the cage for a significant period of time before making infrequent jumps out of the cage. This work suggests that fluctuation in mesh size, as well as average mesh size itself, plays an important role in determining the dynamics of molecules and nanoparticles that are embedded in tightly meshed polymer networks.
| S-EPMC7569880 | biostudies-literature
Project description:Molecular dynamics simulations were performed for the prediction of the finite-size effects of Maxwell-Stefan diffusion coefficients of molecular mixtures and a wide variety of binary Lennard-Jones systems. A strong dependency of computed diffusivities on the system size was observed. Computed diffusivities were found to increase with the number of molecules. We propose a correction for the extrapolation of Maxwell-Stefan diffusion coefficients to the thermodynamic limit, based on the study by Yeh and Hummer ( J. Phys. Chem. B , 2004 , 108 , 15873 - 15879 ). The proposed correction is a function of the viscosity of the system, the size of the simulation box, and the thermodynamic factor, which is a measure for the nonideality of the mixture. Verification is carried out for more than 200 distinct binary Lennard-Jones systems, as well as 9 binary systems of methanol, water, ethanol, acetone, methylamine, and carbon tetrachloride. Significant deviations between finite-size Maxwell-Stefan diffusivities and the corresponding diffusivities at the thermodynamic limit were found for mixtures close to demixing. In these cases, the finite-size correction can be even larger than the simulated (finite-size) Maxwell-Stefan diffusivity. Our results show that considering these finite-size effects is crucial and that the suggested correction allows for reliable computations.
| S-EPMC5943679 | biostudies-literature
Project description:The high demand for biodegradable implants in bone fracture fixations has dramatically increased the use of polymers for biomedical applications as well. However, the replacement of stainless steel and titanium screws by biodegradable materials represents one of the most critical aspects of biomechanics. In this study, the mechanical behavior of polycaprolactone (PCL) in tension and compression is examined. Driven by the advanced technology of computational mechanics, the fixation of the posterior malleolus fracture has been designed and analyzed. The core idea depicts the static analysis of screws made of PCL fixed in the ankle joint. The focus of the study is on this bio-absorbable, polymer-based material performance under constant compression. Parametric analysis is employed for the optimization of the PCL scaffold. Future studies will focus on the experimental verification of the numerical analysis results.
| S-EPMC7151589 | biostudies-literature
Project description:Considering necessary fundamental and structural changes in the production and manufacturing industries to fulfill the industry 4.0 paradigm, the proposal of new ideas and frameworks for operations management of production and manufacturing system is inevitable. This research focuses on traditional methods proposed for storage assignment problem and struggles for new methods and definitions for industry 4.0 based storage assignment concepts. At the first step, the paper proposes a new definition of storage assignment and layout problem for fulfilling storage mechanism agility in terms of automated store and retrieval process (AS/RS) in modern inventories. Then considering the shortcomings of traditional algorithms for storage assignment problem, the paper contributes a new algorithm called SAO/FEM (storage assignment optimization technique), inspired from mechanical engineering discipline for analysis and optimization of storage assignment problem. The proposed new algorithm about stress distribution analogy, and the help of the Finite Element Method and minimum total potential energy theory, proposes a new model for storage assignment optimization. The efficiency of the proposed algorithm in terms of calculation time and the best answer investigated through numerical examples. The article has developed an application for SAO/FEM algorithm as a value creation module and applied new optimized storage positioning in the warehouses.
| S-EPMC7959603 | biostudies-literature