Divergence-free tangential finite element methods for incompressible flows on surfaces.
Ontology highlight
ABSTRACT:
SUBMITTER: Lederer PL
PROVIDER: S-EPMC8611805 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Computational modeling of tissue-scale cardiac electrophysiology requires numerically converged solutions to avoid spurious artifacts. The steep gradients inherent to cardiac action potential propagation necessitate fine spatial scales and therefore a substantial computational burden. The use of high-order interpolation methods has previously been proposed for these simulations due to their theoretical convergence advantage. In this study, we compare the convergence behavior of linear Lagrange, cubic Hermite, and the newly proposed cubic Hermite-style serendipity interpolation methods for finite element simulations of the cardiac monodomain equation. The high-order methods reach converged solutions with fewer degrees of freedom and longer element edge lengths than traditional linear elements. Additionally, we propose a dimensionless number, the cell Thiele modulus, as a more useful metric for determining solution convergence than element size alone. Finally, we use the cell Thiele modulus to examine convergence criteria for obtaining clinically useful activation patterns for applications such as patient-specific modeling where the total activation time is known a priori.
| S-EPMC4525671 | biostudies-literature
Project description:Shape plays an important role in determining the biomechanical response of a structure. Specimen-specific finite element (FE) models have been developed to capture the details of the shape of biological structures and predict their biomechanics. Shape, however, can vary considerably across individuals or change due to aging or disease, and analysis of the sensitivity of specimen-specific models to these variations has proven challenging. An alternative to specimen-specific representation has been to develop generic models with simplified geometries whose shape is relatively easy to parameterize, and can therefore be readily used in sensitivity studies. Despite many successful applications, generic models are limited in that they cannot make predictions for individual specimens. We propose that it is possible to harness the detail available in specimen-specific models while leveraging the power of the parameterization techniques common in generic models. In this work we show that this can be accomplished by using morphing techniques to parameterize the geometry of specimen-specific FE models such that the model shape can be varied in a controlled and systematic way suitable for sensitivity analysis. We demonstrate three morphing techniques by using them on a model of the load-bearing tissues of the posterior pole of the eye. We show that using relatively straightforward procedures these morphing techniques can be combined, which allows the study of factor interactions. Finally, we illustrate that the techniques can be used in other systems by applying them to morph a femur. Morphing techniques provide an exciting new possibility for the analysis of the biomechanical role of shape, independently or in interaction with loading and material properties.
| S-EPMC2813388 | biostudies-literature
Project description:This study investigates the helical secondary flows in the aortic arch using finite element analysis. The relationship between helical flow and the configuration of the aorta in patients of whose three-dimensional images constructed from computed tomography scans was examined. A finite element model of the pressurized root, arch, and supra-aortic vessels was developed to simulate the pattern of helical secondary flows. Calculations indicate that most of the helical secondary flow was formed in the ascending aorta. Angle α between the zero reference point and the aortic ostium (correlation coefficient (r) = -0.851, P = 0.001), the dispersion index of the cross section of the ascending (r = 0.683, P = 0.021) and descending aorta (r = 0.732, P = 0.010), all correlated closely with the presence of helical flow (P < 0.05). Stepwise multiple linear regression analysis confirmed angel α to be independently associated with the helical flow pattern in therein (standardized coefficients = -0.721, P = 0.023). The presence of helical fluid motion based on the atherosclerotic risks of patients, including those associated with diabetes, hypertension, hyperlipidemia, or renal insufficiency, was also evaluated. Numerical simulation of the flow patterns in aortas incorporating the atherosclerotic risks may better explain the mechanism of formation of helical flows and provide insight into causative factors that underlie them.
| S-EPMC4000403 | biostudies-other
Project description:Designing numerical methods for incompressible fluid flow involving moving interfaces, for example, in the computational modeling of bubble dynamics, swimming organisms, or surface waves, presents challenges due to the coupling of interfacial forces with incompressibility constraints. A class of methods, denoted interfacial gauge methods, is introduced for computing solutions to the corresponding incompressible Navier-Stokes equations. These methods use a type of "gauge freedom" to reduce the numerical coupling between fluid velocity, pressure, and interface position, allowing high-order accurate numerical methods to be developed more easily. Making use of an implicit mesh discontinuous Galerkin framework, developed in tandem with this work, high-order results are demonstrated, including surface tension dynamics in which fluid velocity, pressure, and interface geometry are computed with fourth-order spatial accuracy in the maximum norm. Applications are demonstrated with two-phase fluid flow displaying fine-scaled capillary wave dynamics, rigid body fluid-structure interaction, and a fluid-jet free surface flow problem exhibiting vortex shedding induced by a type of Plateau-Rayleigh instability. The developed methods can be generalized to other types of interfacial flow and facilitate precise computation of complex fluid interface phenomena.
| S-EPMC4928983 | biostudies-literature
Project description:Intelligent mobile sensors, such as uninhabited aerial or underwater vehicles, are becoming prevalent in environmental sensing and monitoring applications. These active sensing platforms operate in unsteady fluid flows, including windy urban environments, hurricanes and ocean currents. Often constrained in their actuation capabilities, the dynamics of these mobile sensors depend strongly on the background flow, making their deployment and control particularly challenging. Therefore, efficient trajectory planning with partial knowledge about the background flow is essential for teams of mobile sensors to adaptively sense and monitor their environments. In this work, we investigate the use of finite-horizon model predictive control (MPC) for the energy-efficient trajectory planning of an active mobile sensor in an unsteady fluid flow field. We uncover connections between trajectories optimized over a finite-time horizon and finite-time Lyapunov exponents of the background flow, confirming that energy-efficient trajectories exploit invariant coherent structures in the flow. We demonstrate our findings on the unsteady double gyre vector field, which is a canonical model for chaotic mixing in the ocean. We present an exhaustive search through critical MPC parameters including the prediction horizon, maximum sensor actuation, and relative penalty on the accumulated state error and actuation effort. We find that even relatively short prediction horizons can often yield energy-efficient trajectories. We also explore these connections on a three-dimensional flow and ocean flow data from the Gulf of Mexico. These results are promising for the adaptive planning of energy-efficient trajectories for swarms of mobile sensors in distributed sensing and monitoring.
| S-EPMC8808707 | biostudies-literature
Project description:Linear elastic fracture modeling coupled with empirical material tensile data result in good quantitative agreement with the experimental determination of mode I fracture for both brittle and toughened epoxy nanocomposites. The nanocomposites are comprised of diglycidyl ether of bisphenol A cured with Jeffamine D-230 and some were filled with core-shell rubber nanoparticles of varying concentrations. The quasi-static single-edge notched bending (SENB) test is modeled using both the surface-based cohesive zone (CZS) and extended finite element methods (XFEM) implemented in the Abaqus software. For each material considered, the critical load predicted by the simulated SENB test is used to calculate the mode I fracture toughness. Damage initiates in these models when nodes at the simulated crack tip attain the experimentally measured yield stress. Prediction of fracture processes using a generalized truncated linear traction-separation law (TSL) was significantly improved by considering the case of a linear softening function. There are no adjustable parameters in the XFEM model. The CZS model requires only optimization of the element displacement at the fracture parameter. Thus, these continuum methods describe these materials in mode I fracture with a minimum number of independent parameters.
| S-EPMC8512390 | biostudies-literature
Project description:Based on the Itô’s isometry and the properties of the solution operator defined by the Mittag-Leffler function, this paper gives a detailed numerical analysis of the finite element method for fractional stochastic Navier–Stokes equations driven by white noise. The discretization in space is derived by the finite element method and the time discretization is obtained by the backward Euler scheme. The noise is approximated by using the generalized \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$L_{2}$\end{document}L2-projection operator. Optimal strong convergence error estimates in the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$L_{2}$\end{document}L2-norm are obtained.
| S-EPMC6208620 | biostudies-other
Project description:PurposeThe study presents an averaged anterior eye geometry model combined with a localised material model that is straightforward, appropriate and amenable for implementation in finite element (FE) modelling.MethodsBoth right and left eye profile data of 118 subjects (63 females and 55 males) aged 22-67 years (38.5 ± 7.6) were used to build an averaged geometry model. Parametric representation of the averaged geometry model was achieved through two polynomials dividing the eye into three smoothly connected volumes. This study utilised the collagen microstructure x-ray data of 6 ex-vivo healthy human eyes, 3 right eyes and 3 left eyes in pairs from 3 donors, 1 male and 2 females aged between 60 and 80 years, to build a localised element-specific material model for the eye.ResultsFitting the cornea and the posterior sclera sections to a 5th-order Zernike polynomial resulted in 21 coefficients. The averaged anterior eye geometry model recorded a limbus tangent angle of 37° at a radius of 6.6 mm from the corneal apex. In terms of material models, the difference between the stresses generated in the inflation simulation up to 15 mmHg in the ring-segmented material model and localised element-specific material model were significantly different (p < 0.001) with the ring-segmented material model recording average Von-Mises stress 0.0168 ± 0.0046 MPa and the localised element-specific material model recording average Von-Mises stress 0.0144 ± 0.0025 MPa.ConclusionsThe study illustrates an averaged geometry model of the anterior human eye that is easy to generate through two parametric equations. This model is combined with a localised material model that can be used either parametrically through a Zernike fitted polynomial or non-parametrically as a function of the azimuth angle and the elevation angle of the eye globe. Both averaged geometry and localised material models were built in a way that makes them easy to implement in FE analysis without additional computation cost compared to the limbal discontinuity so-called idealised eye geometry model or ring-segmented material model.
| S-EPMC10123217 | biostudies-literature
Project description:Introduction: Three-dimensional (3D)-printed custom pelvic implants have become a clinically viable option for patients undergoing pelvic cancer surgery with resection of the hip joint. However, increased clinical utilization has also necessitated improved implant durability, especially with regard to the compression screws used to secure the implant to remaining pelvic bone. This study evaluated six different finite element (FE) screw modeling methods for predicting compression screw pullout and fatigue failure in a custom pelvic implant secured to bone using nine compression screws. Methods: Three modeling methods (tied constraints (TIE), bolt load with constant force (BL-CF), and bolt load with constant length (BL-CL)) generated screw axial forces using functionality built into Abaqus FE software; while the remaining three modeling methods (isotropic pseudo-thermal field (ISO), orthotropic pseudo-thermal field (ORT), and equal-and-opposite force field (FOR)) generated screw axial forces using iterative physics-based relationships that can be implemented in any FE software. The ability of all six modeling methods to match specified screw pretension forces and predict screw pullout and fatigue failure was evaluated using an FE model of a custom pelvic implant with total hip replacement. The applied hip contact forces in the FE model were estimated at two locations in a gait cycle. For each of the nine screws in the custom implant FE model, likelihood of screw pullout failure was predicted using maximum screw axial force, while likelihood of screw fatigue failure was predicted using maximum von Mises stress. Results: The three iterative physics-based modeling methods and the non-iterative Abaqus BL-CL method produced nearly identical predictions for likelihood of screw pullout and fatigue failure, while the other two built-in Abaqus modeling methods yielded vastly different predictions. However, the Abaqus BL-CL method required the least computation time, largely because an iterative process was not needed to induce specified screw pretension forces. Of the three iterative methods, FOR required the fewest iterations and thus the least computation time. Discussion: These findings suggest that the BL-CL screw modeling method is the best option when Abaqus is used for predicting screw pullout and fatigue failure in custom pelvis prostheses, while the iterative physics-based FOR method is the best option if FE software other than Abaqus is used.
| S-EPMC11372789 | biostudies-literature
Project description:Optimal tip sonication settings, namely tip position, input power, and pulse durations, are necessary for temperature sensitive procedures like preparation of viable cell extract. In this paper, the optimum tip immersion depth (20-30% height below the liquid surface) is estimated which ensures maximum mixing thereby enhancing thermal dissipation of local cavitation hotspots. A finite element (FE) heat transfer model is presented, validated experimentally with (R2 > 97%) and used to observe the effect of temperature rise on cell extract performance of E. coli BL21 DE3 star strain and estimate the temperature threshold. Relative yields in the top 10% are observed for solution temperatures maintained below 32°C; this reduces below 50% relative yield at temperatures above 47°C. A generalized workflow for direct simulation using the COMSOL code as well as master plots for estimation of sonication parameters (power input and pulse settings) is also presented.
| S-EPMC9161774 | biostudies-literature