Experimental study of needle–tissue interaction forces: Effect of needle geometries,insertion methods and tissue characteristics

A thorough understanding of needle–tissue interaction mechanics is necessary to optimize needle design, achieve robotically needle steering, and establish surgical simulation system. It is obvious that the interaction is influenced by numerous variable parameters, which are divided into three categories: needle geometries, insertion methods, and tissue characteristics. A series of experiments are performed to explore the effect of influence factors (material samples n=5 for each factor) on the insertion force. Data were collected from different biological tissues and a special tissue-equivalent phantom with similar mechanical properties, using a 1-DOF mechanical testing system instrumented with a 6-DOF force/torque (F/T) sensor. The experimental results indicate that three basic phases (deformation, insertion, and extraction phase) are existent during needle penetration. Needle diameter (0.7–3.2 mm), needle tip (blunt, diamond, conical, and beveled) and bevel angle (10–85°) are turned out to have a great influence on insertion force, so do the insertion velocity (0.5–10 mm/s), drive mode (robot-assisted and hand-held), and the insertion process (interrupted and continuous). Different tissues such as skin, muscle, fat, liver capsule and vessel are proved to generate various force cures, which can contribute to the judgement of the needle position and provide efficient insertion strategy.     

Effect of plaque compositions on fractional flow reserve in a fluid–structure interaction analysis

Biomechanics and Modeling in Mechanobiology - Coronary artery disease involves the reduction of blood flow to the myocardium due to atherosclerotic plaques. The findings of myocardial ischemia may...     

Effect of ConA—receptor interaction on the structure of cell membrane

Recently,the effect of ligand receptor interaction on the membrane structure of liposomes has been studied extensively,However,little is known about how it exists on biological membranes,In this paper,the effect of Concanavalin A(ConA) receptorinteratcion on the structure of cell membranes was studied by Circular DIchrosim(CD) and 31P Nuclear Magnetic Resonance(NMR).CD results of both the purified macrophage membranes and human erythrocyte hgosts(EG) showed that the conformation of membrane proteins changed after ConA binding.For further research,31P-NMR was used to detect the orgainzation of phosp[holipid molecules on macrophage membranes.After ConA binding,the tendercy to form non bilayer structure increased with the amount of ConA.The changes of 31P-NMR spectra of living macrophages might be partly due to the above stated reason too.In addition,ConA-receptor interaction also induced similar results of 31P-NMR spectra in EG.In contrast,wheat germ agglutinin (WGA),another kind of lectin,rarely showed the same influence.     

Systolic fluid–structure interaction model of the congenitally bicuspid aortic valve: assessment of modelling requirements

A transient fluid–structure interaction (FSI) model of a congenitally bicuspid aortic valve has been developed which allows simultaneous calculation of fluid flow and structural deformation. The valve is modelled during the systolic phase (the stage when blood pressure is elevated within the heart to pump blood to the body). The geometry was simplified to represent the bicuspid aortic valve in two dimensions. A congenital bicuspid valve is compared within the aortic root only and within the aortic arch. Symmetric and asymmetric cusps were simulated, along with differences in mechanical properties. A moving arbitrary Lagrange–Euler mesh was used to allow FSI. The FSI model requires blood flow to induce valve opening and induced strains in the region of 10%. It was determined that bicuspid aortic valve simulations required the inclusion of the ascending aorta and aortic arch. The flow patterns developed were sensitive to cusp asymmetry and differences in mechanical properties. Stiffening of the valve amplified peak velocities, and recirculation which developed in the ascending aorta. Model predictions demonstrate the need to take into account the category, including any existing cusp asymmetry, of a congenital bicuspid aortic valve when simulating its fluid flow and mechanics.     

Spatial velocity distributions in pulse-wave propagation based on fluid–structure interaction

In this paper, spatial velocity distributions in pulse-wave propagation based on a fluid–structure interaction model are presented. The investigation is performed using the assumption of laminar flow and a linear-elastic wall. The fluid–structure interaction scheme is constructed using the finite element method. The results show that velocity distributions embody an obvious time delay in an elastic tube model. Further, the fully developed flow is delayed and the velocity values are increased in comparison with a rigid tube model. The increase in the wall thickness makes the time delay between the velocity peaks of different sites smaller while the time delay between the velocity minima is unchanged. Similarly, the time delay between the velocity bottoms is more easily found when decreasing the internal radius. The model gives valid results for spatial velocity distributions, which provide important information for wave propagation.     

Fluid flow in the osteocyte mechanical environment: a fluid–structure interaction approach

Osteocytes are believed to be the primary sensor of mechanical stimuli in bone, which orchestrate osteoblasts and osteoclasts to adapt bone structure and composition to meet physiological loading demands. Experimental studies to quantify the mechanical environment surrounding bone cells are challenging, and as such, computational and theoretical approaches have modelled either the solid or fluid environment of osteocytes to predict how these cells are stimulated in vivo. Osteocytes are an elastic cellular structure that deforms in response to the external fluid flow imposed by mechanical loading. This represents a most challenging multi-physics problem in which fluid and solid domains interact, and as such, no previous study has accounted for this complex behaviour. The objective of this study is to employ fluid–structure interaction (FSI) modelling to investigate the complex mechanical environment of osteocytes in vivo. Fluorescent staining of osteocytes was performed in order to visualise their native environment and develop geometrically accurate models of the osteocyte in vivo. By simulating loading levels representative of vigorous physiological activity ( $3,000\,\upmu \upvarepsilon $ compression and 300 Pa pressure gradient), we predict average interstitial fluid velocities $(\sim 60.5\,\upmu \text{ m/s })$ and average maximum shear stresses $(\sim 11\, \text{ Pa })$ surrounding osteocytes in vivo. Interestingly, these values occur in the canaliculi around the osteocyte cell processes and are within the range of stimuli known to stimulate osteogenic responses by osteoblastic cells in vitro. Significantly our results suggest that the greatest mechanical stimulation of the osteocyte occurs in the cell processes, which, cell culture studies have indicated, is the most mechanosensitive area of the cell. These are the first computational FSI models to simulate the complex multi-physics mechanical environment of osteocyte in vivo and provide a deeper understanding of bone mechanobiology.     

An investigation of dentinal fluid flow in dental pulp during food mastication: simulation of fluid–structure interaction

This study uses fluid–structure interaction (FSI) simulation to investigate the relationship between the dentinal fluid flow in the dental pulp of a tooth and the elastic modulus of masticated food particles and to investigate the effects of chewing rate on fluid flow in the dental pulp. Three-dimensional simulation models of a premolar tooth (enamel, dentine, pulp, periodontal ligament, cortical bone, and cancellous bone) and food particle were created. Food particles with elastic modulus of 2,000 and 10,000 MPa were used, respectively. The external displacement loading $(5\,\upmu \hbox {m})$ was gradually directed to the food particle surface for 1 and 0.1 s, respectively, to simulate the chewing of food particles. The displacement and stress on tooth structure and fluid flow in the dental pulp were selected as evaluation indices. The results show that masticating food with a high elastic modulus results in high stress and deformation in the tooth structure, causing faster dentinal fluid flow in the pulp in comparison with that obtained with soft food. In addition, fast chewing of hard food particles can induce faster fluid flow in the pulp, which may result in dental pain. FSI analysis is shown to be a useful tool for investigating dental biomechanics during food mastication. FSI simulation can be used to predict intrapulpal fluid flow in dental pulp; this information may provide the clinician with important concept in dental biomechanics during food mastication.     

A fast strong coupling algorithm for the partitioned fluid–structure interaction simulation of BMHVs

The numerical simulation of Bileaflet Mechanical Heart Valves (BMHVs) has gained strong interest in the last years, as a design and optimisation tool. In this paper, a strong coupling algorithm for the partitioned fluid–structure interaction simulation of a BMHV is presented. The convergence of the coupling iterations between the flow solver and the leaflet motion solver is accelerated by using the Jacobian with the derivatives of the pressure and viscous moments acting on the leaflets with respect to the leaflet accelerations. This Jacobian is numerically calculated from the coupling iterations. An error analysis is done to derive a criterion for the selection of useable coupling iterations. The algorithm is successfully tested for two 3D cases of a BMHV and a comparison is made with existing coupling schemes. It is observed that the developed coupling scheme outperforms these existing schemes in needed coupling iterations per time step and CPU time.     

A numerical model of cellular blebbing: A volume-conserving,fluid–structure interaction model of the entire cell

In animal cells, blebs are smooth, quasi-hemispherical protrusions of the plasma membrane that form when a section of the membrane detaches from the underlying actin cytoskeleton and is inflated by flowing cytosol. The mechanics behind this common cellular activity are not yet clear. As a first step in the development of a full computational framework, we present a numerical model of overall cell behavior based upon the interaction between a background Newtonian-fluid cytosol and elastic structures modeling the membrane and filaments. The detailed micromechanics of the cytoskeletal network are the subject of future work. Here, the myosin-driven contraction of the actin network is modeled through stressed elastic filaments. Quantitative models of cytoskeletal micromechanics and biochemistry require accurate estimates of local stress and flow conditions. The main contribution of this paper is the development of a computationally efficient fluid–structure interaction model based on operator splitting, to furnish this data. Cytosol volume conservation (as supported by experimental evidence) is enforced through an intermediate energy minimization step. Realistic bleb formation and retraction is observed from this model, offering an alternative formulation to positing complex continuum behavior of the cytoplasm (e.g. poroelastic model of Charras et al., 2008).     

A simple fluid–structure coupling algorithm for the study of the anastomotic mechanics of vascular grafts

Vascular anastomoses constitute a main factor in poor graft performance due to mismatches in distensibility between the host artery and the graft. This work aims at computational fluid–structure investigations of proximal and distal anastomoses of vein grafts and synthetic grafts. Finite element and finite volume models were developed and coupled with a user-defined algorithm. Emphasis was placed on the simplicity of the coupling algorithm. An artery and vein graft showed a larger dilation mismatch than an artery and synthetic graft. The vein graft distended nearly twice as much as the artery while the synthetic graft displayed only approximately half the arterial dilation. For the vein graft, luminal mismatching was aggravated by development of an anastomotic pseudo-stenosis. While this study focused on end-to-end anastomoses as a vehicle for developing the coupling algorithm, it may serve as useful point of departure for further investigations such as other anastomotic configurations, refined modelling of sutures and fully transient behaviour.     

Effect of surface hydroxyls on dimethyl ether synthesis over the γ-Al2O3 in liquid paraffin: a computational study

In a recent paper (Zuo et al., Appl Catal A 408:130–136, 2011), the mechanism of dimethyl ether (DME) synthesis from methanol dehydration over γ-Al₂O₃ (110) was studied using density functional theory (DFT). Using the same method, the effect of surface hydroxyls on γ-Al₂O₃ in liquid paraffin during DME synthesis from methanol dehydration is investigated. It is found that DME is mainly formed from two adsorbed CH₃O groups via methanol dehydrogenation on both dehydrated and hydrated γ-Al₂O₃ in liquid paraffin. No close correlation between catalytic activity and acid intensity was found. Before and after water adsorption at typical catalytic conditions (e.g., 553 K), the reaction rate is not obviously changed on γ-Al₂O₃(100) surface in liquid paraffin, but the reaction rate decreases by about 11 times on the (110) in liquid paraffin. Considering the area of the (110) and (100) surfaces under actual conditions, the catalytic activity decreased mainly because the Al3 sites on the (110) surface gradually become inactive. Catalytic activity decreased mainly due to surface hydrophilicity. The calculated results were consistent with the experiment.

The effect of the entry and re-entry size in the aortic dissection: a two-way fluid–structure interaction simulation

Biomechanics and Modeling in Mechanobiology - Aortic dissection (AD) is one of the most catastrophic cardiovascular diseases. AD occurs when a layer inside the aorta is disrupted and gives rise to...     

Effect of surface-screening parameter of the Yukawa potential model on vapour–liquid phase coexistence and critical-point properties of confined Yukawa fluid

We report the effect of surface-screening parameter of Yukawa potential model on vapour–liquid phase coexistence and critical-point properties of slit–pore-confined Yukawa fluid, using grand canonical transition-matrix Monte Carlo along with the histogram reweighting method. The effect of surface-screening parameter on the vapour–phase coexistence density is insignificant for the studied system. On the other hand, significant effect of surface-screening parameter is observed on liquid phase coexistence density. With increasing surface-screening parameter, liquid phase coexistence density decreases. Critical-point properties have shown monotonic decreasing trends with increase in surface-screening parameter. Moreover, the effect of change of surface-screening parameter is least on critical temperature changes as compared to critical density and critical pressure changes for the studied Yukawa system in this work.     

Effect of zinc salts on the structure–function of actomyosin from pelagic fish

The effect of zinc salts, zinc chloride and zinc sulfate on the structure and ATPase enzyme activity of actomyosin from pelagic fish (Sardinella longiceps) has been investigated. ATPase enzyme activity decreased in the presence of both the zinc salts. The inhibitory effect is present in both pH of 7.0 and 9.0. At concentration of 1 × 10^?3 M of zinc salts, complete inhibition of ATPase enzyme activity is observed. With the increase in temperature from 25 °C to 45 °C the ATPase activity decreased by nearly 80%. The solubility profile of actomyosin in the presence of zinc salts shows a sigmoid pattern as the concentration of both the zinc salts increases. Free SH content of actomyosin decreased with the increase in concentration of zinc salts. Intrinsic fluorescence indicated significant decrease in relative fluorescence intensity of actomyosin. This indicates significant alterations in the structure of actomyosin. Analysis of secondary structure also indicates significant alteration in the α-helical content upon binding of both zinc salts.     

Numerical simulation of fluid–structure interaction in bypassed DeBakey III aortic dissection

It was found that bypass graft alone could achieve great effects in treating aortic dissection. In order to investigate the mechanical mechanism and the haemodynamic validity of the bypassing treatment for DeBakey III aortic dissection, patient-specific models of DeBakey III aortic dissection treated with different bypassing strategies were constructed. One of the bypassing strategies is bypassing between ascending aorta and abdominal aorta, and the other is bypassing between left subclavian artery and abdominal aorta. Numerical simulations under physiological flow conditions based on fluid–structure interaction were performed using finite element method. The results show that blood flow velocity, pressure and vessel wall displacement of false lumen are all reduced after bypassing. This phenomenon indicates that bypassing is an effective surgery for the treatment of DeBakey III aortic dissection. The effectiveness to cure through lumen is better when bypassing between left subclavian artery and abdominal aorta, while the effectiveness to cure blind lumen is better when bypassing between ascending aorta and abdominal aorta.     

Advances in the study of protein–DNA interaction

Protein-DNA interaction plays an important role in many biological processes. The classical methods and the novel technologies advanced have been developed for the interaction of protein-DNA. Recent developments of these methods and research achievements have been reviewed in this paper.     

Effect of glycerol–water binary mixtures on the structure and dynamics of protein solutions

We have performed 20?ns of fully atomistic molecular dynamics simulations of Hen Egg-White Lysozyme in 0, 10, 20, 30, and 100% by weight of glycerol in water to better understand the microscopic physics behind the bioprotection offered by glycerol to naturally occuring biological systems. The solvent exposure of protein surface residues changes when glycerol is introduced. The dynamic behavior of the protein, as quantified by the incoherent intermediate scattering function, shows a nonmonotonic dependence on glycerol content. The fluctuations of the protein residues with respect to each other were found to be similar in all water-containing solvents, but different from the pure glycerol case. The increase in the number of protein–glycerol hydrogen bonds in glycerol–water binary mixtures explains the slowing down of protein dynamics as the glycerol content increases. We also explored the dynamic behavior of the hydration layer. We show that the short length scale dynamics of this layer are insensitive to glycerol concentration. However, the long length scale behavior shows a significant dependence on glycerol content. We also provide insights into the behavior of bound and mobile water molecules.     

Effect of nutrients on the accumulation of glutamyl endopeptidase in the culture liquid ofBacillus intermedius 3–19

The effect of nutrients and growth conditions on the accumulation of glutamyl endopeptidase in the culture liquid ofBacillus intermedius 3–19 was studied. Glucose and other readily metabolizable carbon sources were found to suppress the production of the enzyme, whereas inorganic phosphate and ammonium cations enhanced it. Protein substrates, such as casein, gelatin, and hemoglobin, did not affect enzyme production. Some bivalent cations (Ca²⁺, Mg²⁺, Co²⁺) increased the production of glutamyl endopeptidase, but others (Zn²⁺, Fe²⁺, Cu²⁺) acted in the opposite way. The rate of enzyme accumulation in the culture liquid increased as the growth rate of the bacterium decreased, so that the maximum enzyme activity was observed in the stationary growth phase. Based on the results of this investigation, an optimal medium for the maximum production of glutamyl endopeptidase byB. intermedius 3–19 was elaborated.