Simulation of the interaction between blood flow and atherosclerotic plaque

It has been well accepted that over 50% of cerebral ischemic events are the result of rupture of vulnerable carotid atheroma and subsequent thrombosis. Such strokes are potentially preventable by carotid interventions. Selection of patients for intervention is currently based on the severity of carotid luminal stenosis. It has been, however, widely accepted that luminal stenosis alone may not be an adequate predictor of risk. To evaluate the effects of degree of luminal stenosis and plaque morphology on plaque stability, we used a coupled nonlinear time-dependent model with flow-plaque interaction simulation to perform flow and stress/strain analysis for stenotic artery with a plaque. The Navier-Stokes equations in the Arbitrary Lagrangian-Eulerian (ALE) formulation were used as the governing equations for the fluid. The Ogden strain energy function was used for both the fibrous cap and the lipid pool. The plaque Principal stresses and flow conditions were calculated for every case when varying the fibrous cap thickness from 0.1 to 2mm and the degree of luminal stenosis from 10% to 90%. Severe stenosis led to high flow velocities and high shear stresses, but a low or even negative pressure at the throat of the stenosis. Higher degree of stenosis and thinner fibrous cap led to larger plaque stresses, and a 50% decrease of fibrous cap thickness resulted in a 200% increase of maximum stress. This model suggests that fibrous cap thickness is critically related to plaque vulnerability and that, even within presence of moderate stenosis, may play an important role in the future risk stratification of those patients when identified in vivo using high resolution MR imaging.

How critical is fibrous cap thickness to carotid plaque stability? A flow-plaque interaction model

Background and Purpose Acute cerebral ischemic events are associated with rupture of vulnerable carotid atheroma and subsequent thrombosis. Factors such as luminal stenosis and fibrous cap thickness have been thought to be important risk factors for plaque rupture. We used a flow-structure interaction model to simulate the interaction between blood flow and atheromatous plaque to evaluate the effect of the degree of luminal stenosis and fibrous cap thickness on plaque vulnerability. Methods A coupled nonlinear time-dependent model with a flow-plaque interaction simulation was used to perform flow and stress/strain analysis in a stenotic carotid artery model. The stress distribution within the plaque and the flow conditions within the vessel were calculated for every case when varying the fibrous cap thickness from 0.1 to 2 mm and the degree of luminal stenosis from 10% to 95%. A rupture stress of 300 kPa was chosen to indicate a high risk of plaque rupture. A 1-sample t test was used to compare plaque stresses with the rupture stress. Results High stress concentrations were found in the plaques in arteries with >70% degree of stenosis. Plaque stresses in arteries with 30% to 70% stenosis increased exponentially as fibrous cap thickness decreased. A decrease of fibrous cap thickness from 0.4 to 0.2 mm resulted in an increase of plaque stress from 141 to 409 kPa in a 40% degree stenotic artery. Conclusions There is an increase in plaque stress in arteries with a thin fibrous cap. The presence of a moderate carotid stenosis (30% to 70%) with a thin fibrous cap indicates a high risk for plaque rupture. Patients in the future may be risk stratified by measuring both fibrous cap thickness and luminal stenosis.

Nonlinear analysis for the pre- and post-yield behaviour of a composite structure with the refined plastic hinge approach

The computational technique of the full ranges of the second-order inelastic behaviour evaluation of steel-concrete composite structure is not always sought forgivingly, and therefore it hinders the development and application of the performance-based design approach for the composite structure. To this end, this paper addresses of the advanced computational technique of the higher-order element with the refined plastic hinges to capture the all-ranges behaviour of an entire steel-concrete composite structure. Moreover, this paper presents the efficient and economical cross-section analysis to evaluate the element section capacity of the non-uniform and arbitrary composite section subjected to the axial and bending interaction. Based on the same single algorithm, it can accurately and effectively evaluate nearly continuous interaction capacity curve from decompression to pure bending technically, which is the important capacity range but highly nonlinear. Hence, this cross-section analysis provides the simple but unique algorithm for the design approach. In summary, the present nonlinear computational technique can simulate both material and geometric nonlinearities of the composite structure in the accurate, efficient and reliable fashion, including partial shear connection and gradual yielding at pre-yield stage, plasticity and strain-hardening effect due to axial and bending interaction at post-yield stage, loading redistribution, second-order P-δ and P-Δ effect, and also the stiffness and strength deterioration. And because of its reliable and accurate behavioural evaluation, the present technique can be extended for the design of the high-strength composite structure and potentially for the fibre-reinforced concrete structure.

Nonlinear saturation of ion-acoustic instability

It is proposed that the wave mediated indirect wave-particle interaction may be responsible for nonlinear saturation of current driven low frequency ion-acoustic turbulence. This process decreases the growth rate and increases the damping rate of the wave. Comparison has been made with some experiments.

Nonlinear saturation of ion-acoustic instability

It is proposed that the wave mediated indirect wave-particle interaction may be responsible for nonlinear saturation of current driven low frequency ion-acoustic turbulence. This process decreases the growth rate and increases the damping rate of the wave. Comparison has been made with some experiments.

Nonlinear aeroelasticity of rotating and flapping wings - a review

The interaction between large deflections, rotation effects and unsteady aerodynamics makes the dynamic analysis of rotating and flapping wing a nonlinear aeroelastic problem. This problem is governed by nonlinear periodic partial differential equations whose solution is needed to calculate the response and loads acting on vehicles using rotary or flapping wings for lift generation. We look at three important problems in this paper. The first problem shows the effect of nonlinear phenomenon coming from piezoelectric actuators used for helicopter vibration control. The second problem looks at the propagation on material uncertainty on the nonlinear response, vibration and aeroelastic stability of a composite helicopter rotor. The third problem considers the use of piezoelectric actuators for generating large motions in a dragonfly inspired flapping wing. These problems provide interesting insights into nonlinear aeroelasticity and show the likelihood of surprising phenomenon which needs to be considered during the design of rotary and flapping wing vehicle

Nonlinear dynamic analysis of dragonfly inspired piezoelectrically driven flapping and pitching wing

The nonlinear equations for coupled elastic flapping-twisting motion of a dragonfly in- spired smart flapping wing are used for a flapping wing actuated from the root by a PZT unimorph in the piezofan configuration. Excitation by the piezoelectric harmonic force generates only the flap bending motion, which in turn, induces the elastic twist motion due to interaction between flexural and torsional vibrations modes. An unsteady aerodynamic model is used to obtain the aerodynamic forces. Numerical simulations are performed using a wing whose size is the same as the dragonfly Sympetrum Frequens wing. It is found that the value of average lift reaches to its maximum when the smart flapping wing is excited at a frequency closer to the natural frequency in torsion. Moreover, consideration of the elastic twisting of flapping wing leads to an increase in the lift force. It is also found that the flapping wing generates sufficient lift to support its own weight and carry a small pay- load. Therefore, the piezoelectrically actuated smart flapping wing based on the geometry of Sympetrum Frequens wing and undergoing flapping-twisting motions may be considered as a potential candidate for use in MAV applications.

Chemical vapor detection using nonlinear electrical properties of carbon nanotube bundles

We demonstrate the electrical transport behavior of carbon nanotubes (CNTs) upon exposure to organic analytes (namely ethanol, benzene, acetone and toluene). The resulting nonlinear current-voltage characteristics revealed a power law dependence of the differential conductivity on the applied bias voltage. Moreover, suppression of differential conductivity at zero bias is found to be dependent on different selective analytes. The power law exponent values have been monitored before, during and after exposure to the chemicals, which revealed a reversible change in the number of electron conducting channels. Therefore, the reduction in the number of conductive paths can be attributed to the interaction of the chemical analyte on the CNT surfaces, which causes a decrease in the differential conductivity of the CNT sample. These results demonstrate chemical selectivity of CNTs due to varying electronic interaction with different chemical analytes.

Nonlinear optical second harmonic generation in ZnS quantum dots and observation on optical properties of ZnS/PMMA nanocomposites

ZnS quantum dots (QDs) of different sizes are synthesized by a simple chemical co-precipitation method at room temperature, by varying pH value of the reaction mixture. Samples are characterized by an X-ray diffractometer, transmission electron microscope, energy-dispersive X-ray analysis, etc. Linear optical properties, including UV-visible absorption and photoluminescence emission characteristics, of as-prepared QDs are measured. Size dependent nonlinear optical property, such as second harmonic generation (SHG) of 1064 nm Nd:YAG laser fundamental radiation in the synthesized ZnS QDs, is reported for the first time, to the best of our knowledge, by using the standard Kurtz-Perry powder method. In not to study the possibility of the synthesized ZnS QDs in different device applications ZnS/PMMA (polymethylmethacrylate) nanocomposites are also synthesized. The presence of weak chemical interaction between the polymer matrix and ZnS QDs is confirmed by Fourier transform infrared spectroscopy. Thermal properties of the nanocomposites are studied by differential scanning calorimetry and thermo-gravimetric analysis techniques, which show that the composites are stable up to similar to 300 degrees C temperature. (C) 2013 Elsevier B.V. All rights reserved.

Evolution induced catastrophe in a nonlinear dynamical model of material failure

In order to study the failure of disordered materials, the ensemble evolution of a nonlinear chain model was examined by using a stochastic slice sampling method. The following results were obtained. (1) Sample-specific behavior, i.e. evolutions are different from sample to sample in some cases under the same macroscopic conditions, is observed for various load-sharing rules except in the globally mean field theory. The evolution according to the cluster load-sharing rule, which reflects the interaction between broken clusters, cannot be predicted by a simple criterion from the initial damage pattern and even then is most complicated. (2) A binary failure probability, its transitional region, where globally stable (GS) modes and evolution-induced catastrophic (EIC) modes coexist, and the corresponding scaling laws are fundamental to the failure. There is a sensitive zone in the vicinity of the boundary between the GS and EIC regions in phase space, where a slight stochastic increment in damage can trigger a radical transition from GS to EIC. (3) The distribution of strength is obtained from the binary failure probability. This, like sample-specificity, originates from a trans-scale sensitivity linking meso-scopic and macroscopic phenomena. (4) Strong fluctuations in stress distribution different from that of GS modes may be assumed as a precursor of evolution-induced catastrophe (EIC).

Probing biomolecular interaction forces using an anharmonic acoustic technique for selective detection of bacterial spores.

Receptor-based detection of pathogens often suffers from non-specific interactions, and as most detection techniques cannot distinguish between affinities of interactions, false positive responses remain a plaguing reality. Here, we report an anharmonic acoustic based method of detection that addresses the inherent weakness of current ligand dependant assays. Spores of Bacillus subtilis (Bacillus anthracis simulant) were immobilized on a thickness-shear mode AT-cut quartz crystal functionalized with anti-spore antibody and the sensor was driven by a pure sinusoidal oscillation at increasing amplitude. Biomolecular interaction forces between the coupled spores and the accelerating surface caused a nonlinear modulation of the acoustic response of the crystal. In particular, the deviation in the third harmonic of the transduced electrical response versus oscillation amplitude of the sensor (signal) was found to be significant. Signals from the specifically-bound spores were clearly distinguishable in shape from those of the physisorbed streptavidin-coated polystyrene microbeads. The analytical model presented here enables estimation of the biomolecular interaction forces from the measured response. Thus, probing biomolecular interaction forces using the described technique can quantitatively detect pathogens and distinguish specific from non-specific interactions, with potential applicability to rapid point-of-care detection. This also serves as a potential tool for rapid force-spectroscopy, affinity-based biomolecular screening and mapping of molecular interaction networks.

Nonlinear Constitutive Behavior Of Ferroelectric Materials

The ferroelectric specimen is considered as an aggregation of many randomly oriented domains. According to this mechanism, a multi-domain mechanical model is developed in this paper. Each domain is represented by one element. The applied stress and electric field are taken to be the stress and electric field in the formula of the driving force of domain switching for each element in the specimen. It means that the macroscopic switching criterion is used for calculating the volume fraction of domain switching for each element. By using the hardening relation between the driving force of domain switching and the volume fraction of domain switching calibrated, the volume fraction of domain switching for each element is calculated. Substituting the stress and electric field and the volume fraction of domain switching into the constitutive equation of ferroelectric material, one can easily get the strain and electric displacement for each element. The macroscopic behavior of the ferroelectric specimen is then directly calculated by volume averaging. Meanwhile, the nonlinear finite element analysis for the ferroelectric specimen is carried out. In the finite element simulation, the volume fraction of domain switching for each element is calculated by using the same method mentioned above. The interaction between different elements is taken into account in the finite element simulation and the local stress and electric field for each element is obtained. The macroscopic behavior of the specimen is then calculated by volume averaging. The computation results involve the electric butterfly shaped curves of axial strain versus the axial electric field and the hysteresis loops of electric displacement versus the electric field for ferroelectric specimens under the uniaxial coupled stress and electric field loading. The present theoretical prediction agrees reasonably with the experimental results.

An Efficient Front-Tracking Method For Fully Nonlinear Interfacial Waves

A fully nonlinear and dispersive model within the framework of potential theory is developed for interfacial (2-layer) waves. To circumvent the difficulties arisen from the moving boundary problem a viable technique based on the mixed Eulerian and Lagrangian concept is proposed: the computing area is partitioned by a moving mesh system which adjusts its location vertically to conform to the shape of the moving boundaries but keeps frozen in the horizontal direction. Accordingly, a modified dynamic condition is required to properly compute the boundary potentials. To demonstrate the effectiveness of the current method, two important problems for the interfacial wave dynamics, the generation and evolution processes, are investigated. Firstly, analytical solutions for the interfacial wave generations by the interaction between the barotropic tide and topography are derived and compared favorably with the numerical results. Furthermore simulations are performed for the nonlinear interfacial wave evolutions at various water depth ratios and satisfactory agreement is achieved with the existing asymptotical theories. (c) 2008 Elsevier Inc. All rights reserved.

Derivation Of A Nonlinear Reynolds Stress Model Using Renormalization Group Analysis And Two-Scale Expansion Technique

Adopting Yoshizawa's two-scale expansion technique, the fluctuating field is expanded around the isotropic field. The renormalization group method is applied for calculating the covariance of the fluctuating field at the lower order expansion. A nonlinear Reynolds stress model is derived and the turbulent constants inside are evaluated analytically. Compared with the two-scale direct interaction approximation analysis for turbulent shear flows proposed by Yoshizawa, the calculation is much more simple. The analytical model presented here is close to the Speziale model, which is widely applied in the numerical simulations for the complex turbulent flows.

Nonlinear analysis of pile foundation jacket platform

In this paper, the nonlinear collapse of the BOHAI-8 pile foundation jacket platform has been analyzed. The ultimate load and collapse process of two computational models of the structure are given. One model is of fixed support whose length is eight times the pile leg diameter and the other considers the nonlinearity of the soil-pile interaction.