Macromolecular Transport Through Porous Arterial Walls

Arteries are heterogeneous, composite structures that undergo large cyclic deformations during blood transport. Presence, build-up and consequent rupture of blockages in blood vessels, called atherosclerotic plaques, lead to disruption in the blood flow that can eventually be fatal. Abnormal lipid profile and hypertension are the main risk factors for plaque progression. Treatments span from pharmacological methods, to minimally invasive balloon angioplasty and stent procedures, and finally to surgical alternatives. There is a need to understand arterial disease progression and devise methods to detect, control, treat and manage arterial disease through early intervention. Local delivery through drug eluting stents also provide an attractive option for maintaining vessel integrity and restoring blood flow while releasing controlled amount of drug to reduce and alleviate symptoms. Development of drug eluting stents is hence interesting albeit challenging because it requires an integration of knowledge of mechanical properties with material transport of drug through the arterial wall to produce a desired biochemical effect. Although experimental models are useful in studying such complex multivariate phenomena, numerical models of mass transport in the vessel have proved immensely useful to understand and delineate complex interactions between chemical species, physical parameters and biological variables. The goals of this review are to summarize literature based on studies of mass transport involving low density lipoproteins in the arterial wall. We also discuss numerical models of drug elution from stents in layered and porous arterial walls that provide a unique platform that can be exploited for the design of novel drug eluting stents.

Highly fluorescent GFPm2+-based genome integration-proficient promoter probe vector to study Mycobacterium tuberculosis promoters in infected macrophages

Study of activity of cloned promoters in slow-growing Mycobacterium tuberculosis during long-term growth conditions in vitro or inside macrophages, requires a genome-integration proficient promoter probe vector, which can be stably maintained even without antibiotics, carrying a substrate-independent, easily scorable and highly sensitive reporter gene. In order to meet this requirement, we constructed pAKMN2, which contains mycobacterial codon-optimized gfpm2+ gene, coding for GFPm2+ of highest fluorescence reported till date, mycobacteriophage L5 attP-int sequence for genome integration, and a multiple cloning site. pAKMN2 showed stable integration and expression of GFPm2+ from M. tuberculosis and M. smegmatis genome. Expression of GFPm2+, driven by the cloned minimal promoters of M. tuberculosis cell division gene, ftsZ (MtftsZ), could be detected in the M. tuberculosis/pAKMN2-promoter integrants, growing at exponential phase in defined medium in vitro and inside macrophages. Stable expression from genome-integrated format even without antibiotic, and high sensitivity of detection by flow cytometry and fluorescence imaging, in spite of single copy integration, make pAKMN2 useful for the study of cloned promoters of any mycobacterial species under long-term in vitro growth or stress conditions, or inside macrophages.

Electrochemical Integration of Graphene with Light-Absorbing Copper-Based Thin Films

We present an electrochemical route for the integration of graphene with light-sensitive copper-based alloys used in optoelectronic applications. Graphene grown using chemical vapor deposition (CVD) transferred to glass is found to be a robust substrate on which photoconductive CuxS films of 1-2 mu m thickness can be deposited. The effect of growth parameters on the morphology and photoconductivity of CuxS films is presented. Current-voltage (I-V) characterization and photoconductivity decay experiments are performed with graphene as one contact and silver epoxy as the other.

Numerical analysis of transpiration cooling in carbon dioxide atmosphere at hypersonic Mach numbers

The numerical solutions are obtained for skin friction, heat transfer to the wall and growth of boundary layer along the flat plate by employing two dimensional Navier-Stokes equations governing the hypersonic flow coupled with species continuity equations. Flow fields have been computed along the flat plate in CO2 atmosphere in the presence of transpiration cooling using air and carbon dioxide.

Effect of capillarity and heterogeneity in the numerical modelling of multiphase flow of fluids in unsaturated porous medium

The multiphase flow of fluids in the unsaturated porous medium is considered as a three phase flow of water, NAPL, and air simultaneously in the porous medium. The adaptive solution fully implicit modified sequential method is used for the numerical modelling. The effect of capillarity and heterogeneity effect at the interface between the media is studied and it is observed that the interface criteria has to be taken into account for the correct prediction of NAPL migration especially in heterogeneous media. The modified Newton Raphson method is used for the linearization and Hestines and Steifel Conjugate Gradient method is used as the solver.

Numerical simulation of nano-indentation on rough surfaces

Nano-indentation is a technique used to measure various mechanical properties like hardness, Young's modulus and the adherence of thin films and surface layers. It can be used as a quality control tool for various surface modification techniques like ion-implantation, film deposition processes etc. It is important to characterise the increasing scatter in the data measured at lower penetration depths observed in the nano-indentation, for the technique to be effectively applied. Surface roughness is one of the parameters contributing for the scatter. This paper is aimed at quantifying the nature and the amount of scatter that will be introduced in the measurement due to the roughness of the surface on which the indentation is carried out. For this the surface is simulated using the Weierstrass-Mandelbrot function which gives a self-affine fractal. The contact area of this surface with a conical indenter with a spherical cap at the tip is measured numerically. The indentation process is simulated using the spherical cavity model. This eliminates the indentation size effect observed at the micron and sub-micron scales. It has been observed that there exists a definite penetration depth in relation to the surface roughness beyond which the scatter is reduced such that reliable data could be obtained.

Surface reconstruction from disparate shading: an integration of shape-from-shading and stereopsis

A cooperative integration of stereopsis and shape-from-shading is presented. The integration makes the process of D surface reconstruction better constrained and more reliable. It also obviates the need for surface boundary conditions, and explicit information about the surface albedo and the light source direction, which can now be estimated in an iterative manner

A numerical study of variable coefficient elliptic Cauchy problem via projection method

In this paper, we investigate a numerical method for the solution of an inverse problem of recovering lacking data on some part of the boundary of a domain from the Cauchy data on other part for a variable coefficient elliptic Cauchy problem. In the process, the Cauchy problem is transformed into the problem of solving a compact linear operator equation. As a remedy to the ill-posedness of the problem, we use a projection method which allows regularization solely by discretization. The discretization level plays the role of regularization parameter in the case of projection method. The balancing principle is used for the choice of an appropriate discretization level. Several numerical examples show that the method produces a stable good approximate solution.

Numerical Investigation of Heat Transfer from a Two-Dimensional Sudden Expansion Flow Using Nanofluids

Laminar forced convection heat transfer from two-dimensional sudden expansion flow of different nanofluids is studied numerically. The governing equations are solved using the unsteady stream function-vorticity method. The effect of volume fraction of the nanoparticles and type of nanoparticles on heat transfer is examined and found to have a significant impact. Local and average Nusselt numbers are reported in connection with various nanoparticle, volume fraction, and Reynolds number for expansion ratio 2. The Nusselt number reaches peak values near the reattachment point and reaches asymptotic value in the downstream. Bottom wall eddy and volume fraction shows a significant impact on the average Nusselt number.

A Numerical Scheme for Three-Dimensional Front Propagation and Control of Jordan Mode

As an example of a front propagation, we study the propagation of a three-dimensional nonlinear wavefront into a polytropic gas in a uniform state and at rest. The successive positions and geometry of the wavefront are obtained by solving the conservation form of equations of a weakly nonlinear ray theory. The proposed set of equations forms a weakly hyperbolic system of seven conservation laws with an additional vector constraint, each of whose components is a divergence-free condition. This constraint is an involution for the system of conservation laws, and it is termed a geometric solenoidal constraint. The analysis of a Cauchy problem for the linearized system shows that when this constraint is satisfied initially, the solution does not exhibit any Jordan mode. For the numerical simulation of the conservation laws we employ a high resolution central scheme. The second order accuracy of the scheme is achieved by using MUSCL-type reconstructions and Runge-Kutta time discretizations. A constrained transport-type technique is used to enforce the geometric solenoidal constraint. The results of several numerical experiments are presented, which confirm the efficiency and robustness of the proposed numerical method and the control of the Jordan mode.

Two-dimensional flow past circular cylinders using finite volume lattice Boltzmann formulation

In this article, an extension to the total variation diminishing finite volume formulation of the lattice Boltzmann equation method on unstructured meshes was presented. The quadratic least squares procedure is used for the estimation of first-order and second-order spatial gradients of the particle distribution functions. The distribution functions were extrapolated quadratically to the virtual upwind node. The time integration was performed using the fourth-order RungeKutta procedure. A grid convergence study was performed in order to demonstrate the order of accuracy of the present scheme. The formulation was validated for the benchmark two-dimensional, laminar, and unsteady flow past a single circular cylinder. These computations were then investigated for the low Mach number simulations. Further validation was performed for flow past two circular cylinders arranged in tandem and side-by-side. Results of these simulations were extensively compared with the previous numerical data. Copyright (C) 2011 John Wiley & Sons, Ltd.

Geometrically non-linear dynamics of composite four-bar mechanisms

This work intends to demonstrate the importance of a geometrically nonlinear cross-sectional analysis of certain composite beam-based four-bar mechanisms in predicting system dynamic characteristics. All component bars of the mechanism are made of fiber reinforced laminates and have thin rectangular cross-sections. They could, in general, be pre-twisted and/or possess initial curvature, either by design or by defect. They are linked to each other by means of revolute joints. We restrict ourselves to linear materials with small strains within each elastic body (beam). Each component of the mechanism is modeled as a beam based on geometrically non-linear 3-D elasticity theory. The component problems are thus split into 2-D analyses of reference beam cross-sections and non-linear 1-D analyses along the three beam reference curves. For the thin rectangular cross-sections considered here, the 2-D cross-sectional non-linearity is also overwhelming. This can be perceived from the fact that such sections constitute a limiting case between thin-walled open and closed sections, thus inviting the non-linear phenomena observed in both. The strong elastic couplings of anisotropic composite laminates complicate the model further. However, a powerful mathematical tool called the Variational Asymptotic Method (VAM) not only enables such a dimensional reduction, but also provides asymptotically correct analytical solutions to the non-linear cross-sectional analysis. Such closed-form solutions are used here in conjunction with numerical techniques for the rest of the problem to predict multi-body dynamic responses more quickly and accurately than would otherwise be possible. The analysis methodology can be viewed as a three-step procedure: First, the cross-sectional properties of each bar of the mechanism is determined analytically based on an asymptotic procedure, starting from Classical Laminated Shell Theory (CLST) and taking advantage of its thin strip geometry. Second, the dynamic response of the non-linear, flexible four-bar mechanism is simulated by treating each bar as a 1-D beam, discretized using finite elements, and employing energy-preserving and -decaying time integration schemes for unconditional stability. Finally, local 3-D deformations and stresses in the entire system are recovered, based on the 1-D responses predicted in the previous step. With the model, tools and procedure in place, we identify and investigate a few four-bar mechanism problems where the cross-sectional non-linearities are significant in predicting better and critical system dynamic characteristics. This is carried out by varying stacking sequences (i.e. the arrangement of ply orientations within a laminate) and material properties, and speculating on the dominating diagonal and coupling terms in the closed-form non-linear beam stiffness matrix. A numerical example is presented which illustrates the importance of 2-D cross-sectional non-linearities and the behavior of the system is also observed by using commercial software (I-DEAS + NASTRAN + ADAMS). (C) 2012 Elsevier Ltd. All rights reserved.

Numerical analysis on thermoacoustic prime movers for development of pulse tube cryocoolers

Thermo Acoustic Prime Movers (TAPMs) are being considered as the ideal choice for driving the Pulse Tube Cryocoolers replacing the conventional compressors. The advantages are the absence of moving components and they can be driven by low grade energy as such as fuel, gas, solar energy, waste heat etc. While the development of such TAPMs is in progress in our laboratory, their design and fabrication should be guided by numerical modeling and this may be done by several methods such as solving the energy equation 1], enthalpy flow model 2], CFD 3], etc. We have used CFD technique, since it provides a better insight into the velocity and temperature profiles. The analysis is carried out by varying parameters such as (a) temperature difference across the stack, (b) stack and resonator lengths and (c) different working fluids such as air, nitrogen, argon etc. The theoretical results are compared with the experimental data wherever possible and they are in reasonably good agreement with each other. The analysis indicate that (i) larger temperature difference across the stack leads to increased acoustic amplitude, (ii) longer resonator leads to decrease in frequency with lesser amplitude and (iii) there exists an optimal stack length for the best performance of TAPM. These results are presented here.

A numerical approach to determine the sufficiency of given boundary data sets for uniquely estimating interior elastic properties

We consider an inverse elasticity problem in which forces and displacements are known on the boundary and the material property distribution inside the body is to be found. In other words, we need to estimate the distribution of constitutive properties using the finite boundary data sets. Uniqueness of the solution to this problem is proved in the literature only under certain assumptions for a given complete Dirichlet-to-Neumann map. Another complication in the numerical solution of this problem is that the number of boundary data sets needed to establish uniqueness is not known even under the restricted cases where uniqueness is proved theoretically. In this paper, we present a numerical technique that can assess the sufficiency of given boundary data sets by computing the rank of a sensitivity matrix that arises in the Gauss-Newton method used to solve the problem. Numerical experiments are presented to illustrate the method.

Directional Diffusion Regulator (DDR) for some numerical solvers of hyperbolic conservation laws

A computational tool called ``Directional Diffusion Regulator (DDR)'' is proposed to bring forth real multidimensional physics into the upwind discretization in some numerical schemes of hyperbolic conservation laws. The direction based regulator when used with dimension splitting solvers, is set to moderate the excess multidimensional diffusion and hence cause genuine multidimensional upwinding like effect. The basic idea of this regulator driven method is to retain a full upwind scheme across local discontinuities, with the upwind bias decreasing smoothly to a minimum in the farthest direction. The discontinuous solutions are quantified as gradients and the regulator parameter across a typical finite volume interface or a finite difference interpolation point is formulated based on fractional local maximum gradient in any of the weak solution flow variables (say density, pressure, temperature, Mach number or even wave velocity etc.). DDR is applied to both the non-convective as well as whole unsplit dissipative flux terms of some numerical schemes, mainly of Local Lax-Friedrichs, to solve some benchmark problems describing inviscid compressible flow, shallow water dynamics and magneto-hydrodynamics. The first order solutions consistently improved depending on the extent of grid non-alignment to discontinuities, with the major influence due to regulation of non-convective diffusion. The application is also experimented on schemes such as Roe, Jameson-Schmidt-Turkel and some second order accurate methods. The consistent improvement in accuracy either at moderate or marked levels, for a variety of problems and with increasing grid size, reasonably indicate a scope for DDR as a regular tool to impart genuine multidimensional upwinding effect in a simpler framework. (C) 2012 Elsevier Inc. All rights reserved.