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:This study aimed to evaluate fatigue of three stent designs when various forces are applied and perform a comparative analyses. A computer simulation using the finite element method was performed. In particular, we constructed a three-dimensional finite element model of nitinol stents with three designs (S6: single-woven wire, wire diameter: 0.006 inch; D6: double-woven wire, wire diameter: 0.006 inch, and D7: double-woven wire, wire diameter: 0.007 inch) that are used to treat canine tracheal collapse (TC). The stents were subjected to a 200 mmHg compression force, a pure torsion force in a perpendicular direction, and a bending-torsion force combining perpendicular and axial forces. The von Mises stress was calculated to evaluate the extent of stent displacement, and Goodman diagrams were plotted to compare fatigue life cycles. D7 exhibited a longer fatigue life compared to S6 and D6. Under compression, pure torsion, and bending-torsion forces, displacement was the smallest for D7, followed by D6 and S6. Similarly, the fatigue life was the longest for D7, followed by D6 and S6. S6 showed the greatest displacement when subjected to external forces; among stents designed using the same wire, D6 displayed less displacement than S6, and D7 exhibited superior fatigue life when subjected to varying degrees of force. This study showed that the structural stability and fatigue life of stents could be effectively compared using finite element method D7 has the greatest stability and structural rigidity under cyclic load.
| S-EPMC10882158 | 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:BackgroundThe spinal hybrid elastic (SHE) rod dynamic stabilization system can provide sufficient spine support and less adjacent segment stress. This study aimed to investigate the biomechanical effects after the internal fracture of SHE rods using finite element analysis.MethodsA three-dimensional nonlinear finite element model was developed. The SHE rod comprises an inner nitinol stick (NS) and an outer polycarbonate urethane (PCU) shell (PS). The fracture was set at the caudal third portion of the NS, where the maximum stress occurred. The resultant intervertebral range of motion (ROM), intervertebral disc stress, facet joint contact force, screw stress, NS stress, and PCU stress were analyzed.ResultsWhen compared with the intact spine model, the overall trend was that the ROM, intervertebral disc stress, and facet joint force decreased in the implanted level and increased in the adjacent level. When compared with the Ns-I, the trend in the Ns-F decreased and remained nearly half effect. Except for torsion, the PS stress of the Ns-F increased because of the sharing of NS stress after the NS fracture.ConclusionsThe study concluded the biomechanical effects still afford nearly sufficient spine support and gentle adjacent segment stress after rod fracture in a worst-case scenario of the thinnest PS of the SHE rod system.
| S-EPMC9413940 | 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