Solid and Hollow Gold Nanostructures for Nanomedicine: Comparison of Photothermal Properties

The photothermal properties of solid and hollow gold nanostructures represented by colloidal solutions of spherical nanoparticles, nanoshells, and nanocages upon irradiation with a 100 mW 808 nm continuous-wave laser for the first time were experimentally compared under identical optical density and nanoparticle concentration conditions. Accompanying computer modeling of light absorption by the studied gold nanostructures revealed the general parameters influencing the photothermal efficiency, which is of significance for nanomedical applications. The spectral position of localized plasmonic excitations of the studied nanostructures ranged from 518 nm for solid gold nanoparticles to 718 nm for gold nanocages, which provided a possibility to observe a direct influence of the wavelength proximity between the localized surface plasmon resonance and laser line on the heat generation capability of the nanostructures. As a result, the best photothermal efficiency was registered for gold nanocages, which proves them as an efficient photothermal treatment agent and a possible candidate to build a nanocarrier platform for drug delivery with a controlled release. Light absorption modeling demonstrated an existence of optimal wall thickness for gold nanoshells that should lead to the maximum photothermal efficiency when irradiated with 808 nm light, which varied from about 0.1 to 0.4 in units of external nanoshell radius with an increase of the wall porosity. Additionally, computer modeling results show that increased wall porosity should lead to enhanced photothermal efficiency of polydisperse colloidal solutions of hollow gold nanostructures.     

In vitro determination of the pressure-diameter relationship and velocity profiles by ultrasonic technics. In vivo application

A good knowledge of arterial flow mechanics and of the phenomena associated with fluid-boundary interactions is necessary for the determination of some fundamental parameters such as velocity, pressure and pressure-diameter relationship during a cardiac cycle. Ultrasonic techniques were developed on a test bench and directly applied to animals without major modification. On such a test bench allowing a good simulation of physiological type flows, velocity field and pressure-diameter relationship were determined. In vivo application of these techniques allowed a systematic analysis of velocity profiles in the rabbit abdominal aorta and a precise approach of rheological properties of the vascular wall.     

Experimental and numerical study of pulsatile flows through stenosis: wall shear stress analysis.

Different shapes of pulsatile flows through a model of stenosis are experimentally and numerically modeled to validate both methods and to determine the wall shear stress temporal evolution downstream from the stenosis. Two-dimensional velocity measurements are performed in a 75% severity stenosis using a pulsed Doppler ultrasonic velocimeter. Finite element package is employed for the transient numerical simulations. Polynomial method, based on the experimental velocity values, is proposed to determine the wall shear stress temporal evolution. There is a good agreement between the numerical and experimental results. The wall shear stress temporal analysis shows oscillating wall shear stress values during the cycle with high wall shear stress values at the throat of about 120 dyn/cm2, and low values downstream from the stenosis of about - 2.5 dyn/cm2. The key result of the study is that the presence of the stenosis leads the artery to work in a direction which is opposite to the direction of a healthy artery.     

Molecular dynamics simulations of ion separation in nano-channel water flows using an electric field

ABSTRACT

Removal of undesired substances from water is a field of investigation recently focused at the nanoscale. Towards this direction, molecular dynamics simulations are conducted in this paper to investigate unwanted ion removal in nanochannel flows. The simulation method incorporates a Poiseuille-like water/ion flow system at the nanoscale where an electric field, of various magnitudes in the range of E?=?0.25–1.5?V/Å, is applied perpendicular to the flow, leading anions and cations close to the wall regions, similar to the Capacitive De-Ionization method. The time needed for ions to reach equilibrium, i.e. to flow in the region near the walls while pure water flows in the channel interior, is t?=?1.3?ns when E?=?1.5?V/Å and t?=?4.0?ns when E?=?0.25?V/Å, showing a dependency on the value of the electric field. Calculations on density, velocity, and temperature values report on fluid properties to be used in the proposed desalination configuration and could act as a basis to guide novel technological applications and extend to higher scales.     

Molecular dynamics simulations of energy accommodation coefficients for gas flows in nano-channels

The energy accommodation coefficient (EAC) used in thermal boundary condition in micro- and nano-gas flows is reported to be always less than unity and greatly influenced by wall characteristics, for example, the wall temperature. A statistical EAC definition was described to calculate the EAC for thermal conduction in argon gas between two smooth platinum plates from two-dimensional non-equilibrium molecular dynamics simulations. The non-equilibrium EAC at the upper wall was calculated for different upper wall temperatures and a fixed bottom wall temperature. The equilibrium EAC at each temperature can then be extrapolated from a series of non-equilibrium EACs as the temperature difference between the two walls approaches zero. The analyses of the effects of wall temperature for various Knudsen number on non-equilibrium and equilibrium EACs show that, for a given lower bottom wall temperature, the non-equilibrium EAC at a high temperature wall increases with an increase in the wall temperature. For a given wall temperature difference, the non-equilibrium EAC increases with the increase in the wall temperatures. The equilibrium EAC also becomes larger at higher temperatures.     

Examination of nanoflow in rectangular slits

This paper presents simulations of 3D nanoscale flow in rectangular channel with molecular dynamics simulation method. Rectangular cross section is a frequently encountered geometric shape for nanoscale flow problems. For a given cross sectional height h, we change the width w of the rectangular cross section and analyze the influence of w/h on the flow characteristics. The distributions of density, temperature, boundary slip and flow velocity inside the rectangular cross section are investigated in detail. Liquid argon material and Lennard-Jones potential are used in the simulations. The simulation results are also compared with Navier–Stokes solutions for rectangular channel flows.     

Wetting of Planar Surfaces by a Gay-Berne Liquid Crystal

Molecular simulation of the wetting of an unstructured attractive wall by Gay-Berne liquid crystals are reported. Simulations are performed in the grand canonical ensemble on a wide pore at constant temperatures of T *=0.53 and 0.56, corresponding to temperatures below and above the nematic-isotropic-vapor triple point. Close to the coexistence chemical potential, a thick liquid film wets the solid surface. The film is composed of stratified layers of molecules parallel to the solid surface, which follow to a nematic domain at the lower temperature and an isotropic one at the higher temperature. In both cases, the film is in equilibrium with the corresponding vapor phase. Close to the liquid-vapor interface there is a manifest tendency for the molecules to orient themselves parallel to the interface. The adsorption on the wall varies continuously with the thermodynamic parameters considered and no evidence of a first order prewetting transition is observed.     

Molecular dynamics simulation of fluid properties by the streamwise oscillation of the solid wall

In this paper, the hydrodynamics of streamwise wall oscillations on a Couette flow are studied using the molecular dynamics method. Firstly, based on the two-dimensional Couette flow model which is made up of copper wall and argon fluid, the characters of the fluid near the streamwise oscillation wall are simulated under the condition of varied oscillating parameters. By scrupulous data processing, some significative results such as the velocity distribution, the density distribution, the potential energy curves of the flow field and the frictional force of the wall are obtained. Secondly, the mechanism how the wall oscillation brings about change to the frictional drag at liquid–solid interface is investigated. And the results indicate that the frictional drag can be reduced significantly by applying appropriate streamwise oscillation to the solid wall. The drag reduction rate mainly depends on the oscillation parameters. In addition, the decrease in the fluid’s density near the wall is another important reason behind the frictional drag reduction.     

Numerical and experimental models of post-operative realistic flows in stenosed coronary bypasses.

By means of both experimental and finite element methods, we simulated three-dimensional unsteady flows through coronary bypass anastomosis. The host artery includes a stenosis shape located at two different distances of grafting. The inflow rates are issued from in vivo measurements in patients who had undergone coronary bypass surgery a few days before. We provide a comparison between experimental and numerical velocity profiles coupled with the numerical analysis of spatial and temporal wall shear stress evolution. The interaction between the graft and coronary flows has been demonstrated. The phase inflow difference can partly be responsible for specific flow phenomena: jet deflection towards a preferential wall or feedback phenomenon that causes the flapping of the post-stenotic jet during the cardiac cycle. In conclusion, we showed the sensitivity of these typical flows to distance of grafting, inflows waveforms but also to their phase difference.     

In vitro study of modifications of blood flow patterns induced by a bifurcation

Velocity profiles obtained with atheromatous and normal bifurcation castings in the presence of various types of flows are proposed. In the atheromatous bifurcation, with steady flow, we observe radial positive or negative velocities at distance to the wall smaller than 1 mm, which may be attributed to small local eddy motions. The maximum of velocity needs a larger distance from the apex than in the case of the "normal" bifurcation to be again located on the axes. With periodical flows, the effects are strongly damped. The wall velocity gradients on several geometries of tubings are investigated to separate the effects of the local rugosity of the wall from those incidental to the geometry of the bifurcation. The alterations caused by the atheroma do not seem to be induced by local modifications or rugosity, but by slow modifications of the local diameter. As a consequence, the variations of the velocity gradient caused by atheroma in the total bifurcation, are more likely due to distance effects of the geometry itself than to local effects of rugosity.     

Unsteady and three-dimensional simulation of blood flow in the human aortic arch

A three-dimensional and pulsatile blood flow in a human aortic arch and its three major branches has been studied numerically for a peak Reynolds number of 2500 and a frequency (or Womersley) parameter of 10. The simulation geometry was derived from the three-dimensional reconstruction of a series of two-dimensional slices obtained in vivo using CAT scan imaging on a human aorta. The numerical simulations were obtained using a projection method, and a finite-volume formulation of the Navier-Stokes equations was used on a system of overset grids. Our results demonstrate that the primary flow velocity is skewed towards the inner aortic wall in the ascending aorta, but this skewness shifts to the outer wall in the descending thoracic aorta. Within the arch branches, the flow velocities were skewed to the distal walls with flow reversal along the proximal walls. Extensive secondary flow motion was observed in the aorta, and the structure of these secondary flows was influenced considerably by the presence of the branches. Within the aorta, wall shear stresses were highly dynamic, but were generally high along the outer wall in the vicinity of the branches and low along the inner wall, particularly in the descending thoracic aorta. Within the branches, the shear stresses were considerably higher along the distal walls than along the proximal walls. Wall pressure was low along the inner aortic wall and high around the branches and along the outer wall in the ascending thoracic aorta. Comparison of our numerical results with the localization of early atherosclerotic lesions broadly suggests preferential development of these lesions in regions of extrema (either maxima or minima) in wall shear stress and pressure.     

Survey of the Influence of the Width of Urban Branch Roads on the Meeting of Two-Way Vehicle Flows

Branch roads, which are densely distributed in cities, allow for the flow of local traffic and provide connections between the city and outlying areas. Branch roads are typically narrow, and two-way traffic flows on branch roads are thus affected when vehicles traveling in opposite directions meet. This study investigates the changes in the velocities of vehicles when they meet on two-way branch roads. Various widths of branch roads were selected, and their influence on traffic flows was investigated via a video survey. The results show that, depending on the average vehicle velocity, branch roads require different widths to prevent a large decrease in velocity when vehicles meet. When the velocity on a branch road is not high (e.g., the average velocity without meeting is approximately 6 m/s), appropriately increasing the road width will notably increase the meeting velocity. However, when the velocity is high (e.g., the average velocity without meeting is greater than 10 m/s), there is a large decrease in velocity when meeting even if the road surface is wide (6.5 m). This study provides a basis for selecting the width of urban branch roads and the simulation of bidirectional traffic on such roads.     

Numerical investigation of the non-Newtonian blood flow in a bifurcation model with a non-planar branch

The non-Newtonian fluid flow in a bifurcation model with a non-planar daughter branch is investigated by using finite element method to solve the three-dimensional Navier–Stokes equations coupled with a non-Newtonian constitutive model, in which the shear thinning behavior of the blood fluid is incorporated by the Carreau–Yasuda model. The objective of this study is to investigate the influence of the non-Newtonian property of fluid as well as of curvature and out-of-plane geometry in the non-planar daughter vessel on wall shear stress (WSS) and flow phenomena. In the non-planar daughter vessel, the flows are typified by the skewing of the velocity profile towards the outer wall, creating a relatively low WSS at the inner wall. In the downstream of the bifurcation, the velocity profiles are shifted towards the flow divider. The low WSS is found at the inner walls of the curvature and the lateral walls of the bifurcation. Secondary flow patterns that swirl fluid from the inner wall of curvature to the outer wall in the middle of the vessel are also well documented for the curved and bifurcating vessels. The numerical results for the non-Newtonian fluid and the Newtonian fluid with original Reynolds number and the corresponding rescaled Reynolds number are presented. Significant difference between the non-Newtonian flow and the Newtonian flow is revealed; however, reasonable agreement between the non-Newtonian flow and the rescaled Newtonian flow is found. Results of this study support the view that the non-planarity of blood vessels and the non-Newtonian properties of blood are an important factor in hemodynamics and may play a significant role in vascular biology and pathophysiology.     

A two-phase model for water and heat transfer within an intermittently-mixed solid-state fermentation bioreactor with forced aeration

A two-phase dynamic model is developed that describes heat and mass transfer in intermittently-mixed solid-state fermentation bioreactors. The model predicts that in the regions of the bed near the air inlet there can be significant differences in the air and solid temperatures, while in the remainder of the bed the gas and solid phases are much closer to equilibrium, although there can be differences in water activity of around 0.05. The increase in the temperature of the gas as it flows through the bed means that it is impossible to prevent the bed from drying out, even if saturated air is used at the air inlet. The substrate can dry to water activities that severely limit growth, unless the bed is intermittently mixed, with the addition of water to bring the water activity back to the desired value. Under the conditions assumed for the simulation, which was designed to mimic the growth of Aspergillus niger on corn, two mixing events were necessary, one at 17.4 and the other at 27.9 h. Even though such a strategy can minimize the restriction of growth by water-limitation, temperature-limitation remains a problem due to the rapid heating dynamics. The model is obviously a useful tool that can be used to guide scale-up and to test control strategies. Such a model, describing the non-equilibrium situation between the gas and solid phases, has not previously been proposed for solid-state fermentation bioreactors. Models in the literature that assume gas-solid temperature and moisture equilibrium cannot describe the large temperature differences between the gas and solid phase which occur within the bed near the air inlet.     

Numerical modelling of simulated blood flow in idealized composite arterial coronary grafts: transient flow

In composite arterial coronary grafts (CACGs), transport phenomena and geometry may considerably alter blood flow dynamics. CACGs aim at revascularizing pathological arteries according to the human anatomy. However, the exact mechanisms causing the failure of coronary bypass grafting are not yet well elucidated. In the present study, computational fluid dynamics (CFD) techniques are applied for the simulation of multi-branched CACGs under physiologically realistic inflow waveforms. The numerical solution is obtained by a finite-volume method formulated in non-orthogonal, curvilinear coordinates and a multi-grid approach. The geometrical models, consisting of idealized and rigid vessels, include the typical T- and a rather new pi-graft configuration. The stenotic effect is also investigated by comparing computational results for three different degrees of area constriction, namely 25%, 50% and 75%, as well as the case without stenosis. Different grafting distances and various inflow rate ratios are imposed, to give an insight into haemodynamical alterations of CACGs and to study the process of restenosis. The results focus on the interaction between the grafts and coronary flows in terms of spatial and temporal variations of velocity and wall shear stress (WSS) distribution. Prominent variations among the different geometries, concerning the velocity profiles and secondary flow motion, are shown. Moreover, the residual flow emerging from different degrees of area constriction shows that low and oscillating shear stresses may arise for even moderate stenotic fields.     

Porohyperelastic finite element modeling of abdominal aortic aneurysms

Abdominal aortic aneurysm (AAA) is the gradual weakening and dilation of the infrarenal aorta. This disease is progressive, asymptomatic, and can eventually lead to rupture--a catastrophic event leading to massive internal bleeding and possibly death. The mechanical environment present in AAA is currently thought to be important in disease initiation, progression, and diagnosis. In this study, we utilize porohyperelastic (PHE) finite element models (FEMs) to investigate how such modeling can be used to better understand the local biomechanical environment in AAA. A 3D hypothetical AAA was constructed with a preferential anterior bulge assuming both the intraluminal thrombus (ILT) and the AAA wall act as porous materials. A parametric study was performed to investigate how physiologically meaningful variations in AAA wall and ILT hydraulic permeabilities affect luminal interstitial fluid velocities and wall stresses within an AAA. A corresponding hyperelastic (HE) simulation was also run in order to be able to compare stress values between PHE and HE simulations. The effect of AAA size on local interstitial fluid velocity was also investigated by simulating maximum diameters (5.5 cm, 4.5 cm, and 3.5 cm) at the baseline values of ILT and AAA wall permeability. Finally, a cyclic PHE simulation was utilized to study the variation in local fluid velocities as a result of a physiologic pulsatile blood pressure. While the ILT hydraulic permeability was found to have minimal affect on interstitial velocities, our simulations demonstrated a 28% increase and a 20% decrease in luminal interstitial fluid velocity as a result of a 1 standard deviation increase and decrease in AAA wall hydraulic permeability, respectively. Peak interstitial velocities in all simulations occurred on the luminal surface adjacent to the region of maximum diameter. These values increased with increasing AAA size. PHE simulations resulted in 19.4%, 40.1%, and 81.0% increases in peak maximum principal wall stresses in comparison to HE simulations for maximum diameters of 35 mm, 45 mm, and 55 mm, respectively. The pulsatile AAA PHE FEM demonstrated a complex interstitial fluid velocity field the direction of which alternated in to and out of the luminal layer of the ILT. The biomechanical environment within both the aneurysmal wall and the ILT is involved in AAA pathogenesis and rupture. Assuming these tissues to be porohyperelastic materials may provide additional insight into the complex solid and fluid forces acting on the cells responsible for aneurysmal remodeling and weakening.     

MHD Free Convective Boundary Layer Flow of a Nanofluid past a Flat Vertical Plate with Newtonian Heating Boundary Condition

Steady two dimensional MHD laminar free convective boundary layer flows of an electrically conducting Newtonian nanofluid over a solid stationary vertical plate in a quiescent fluid taking into account the Newtonian heating boundary condition is investigated numerically. A magnetic field can be used to control the motion of an electrically conducting fluid in micro/nano scale systems used for transportation of fluid. The transport equations along with the boundary conditions are first converted into dimensionless form and then using linear group of transformations, the similarity governing equations are developed. The transformed equations are solved numerically using the Runge-Kutta-Fehlberg fourth-fifth order method with shooting technique. The effects of different controlling parameters, namely, Lewis number, Prandtl number, buoyancy ratio, thermophoresis, Brownian motion, magnetic field and Newtonian heating on the flow and heat transfer are investigated. The numerical results for the dimensionless axial velocity, temperature and nanoparticle volume fraction as well as the reduced Nusselt and Sherwood number have been presented graphically and discussed. It is found that the rate of heat and mass transfer increase as Newtonian heating parameter increases. The dimensionless velocity and temperature distributions increase with the increase of Newtonian heating parameter. The results of the reduced heat transfer rate is compared for convective heating boundary condition and found an excellent agreement.     

The effect of compliance on wall shear in casts of a human aortic bifurcation

Rigid and compliant casts of a human aortic bifurcation were subjected to physiologically realistic pulsatile fluid flows. At a number of sites near the wall in the approximate median plane of the bifurcation of these models, fluid velocity was measured with a laser Doppler velocimeter, and wall motion (in the case of the compliant cast) was determined with a Reticon linescan camera. The velocity and wall motion data were combined to estimate the instantaneous shear rates at the cast wall. Analysis showed that at the outer walls the cast compliance reduced shear rates, while at the walls of the flow divider the shear rate was increased.     

Finite element analysis of biological soft tissue surrounded by a deformable membrane that controls transmembrane flow

Background

Many biological soft tissues are hydrated porous hyperelastic materials, which consist of a complex solid skeleton with fine voids and fluid filling these voids. Mechanical interactions between the solid and the fluid in hydrated porous tissues have been analyzed by finite element methods (FEMs) in which the mixture theory was introduced in various ways. Although most of the tissues are surrounded by deformable membranes that control transmembrane flows, the boundaries of the tissues have been treated as rigid and/or freely permeable in these studies. The purpose of this study was to develop a method for the analysis of hydrated porous hyperelastic tissues surrounded by deformable membranes that control transmembrane flows.

Results

For this, we developed a new nonlinear finite element formulation of the mixture theory, where the nodal unknowns were the pore water pressure and solid displacement. This method allows the control of the fluid flow rate across the membrane using Neumann boundary condition. Using the method, we conducted a compression test of the hydrated porous hyperelastic tissue, which was surrounded by a flaccid impermeable membrane, and a part of the top surface of this tissue was pushed by a platen. The simulation results showed a stress relaxation phenomenon, resulting from the interaction between the elastic deformation of the tissue, pore water pressure gradient, and the movement of fluid. The results also showed that the fluid trapped by the impermeable membrane led to the swelling of the tissue around the platen.

Conclusions

These facts suggest that our new method can be effectively used for the analysis of a large deformation of hydrated porous hyperelastic material surrounded by a deformable membrane that controls transmembrane flow, and further investigations may allow more realistic analyses of the biological soft tissues, such as brain edema, brain trauma, the flow of blood and lymph in capillaries and pitting edema.