MODELING AND SIMULATION OF HYDRODYNAMIC INTERACTION OF DNA IN A MICRO-FLUIDIC CHANNEL

The motion of DNA (in the bulk solution) and the non-Newtonian effective fluid behavior are considered separately and self-consistently with the fluid motion satisfying the no-slip boundary condition on the surface of the confining geometry in the presence of channel pressure gradients. A different approach has been developed to model DNA in the micro-channel. In this study the DNA is assumed as an elastic chain with its characteristic Young's modulus, Poisson's ratio and density. The force which results from the fluid dynamic pressure, viscous forces and electromotive forces is applied to the elastic chain in a coupled manner. The velocity fields in the micro-channel are influenced by the transport properties. Simulations are carried out for the DNAs attached to the micro-fluidic wall. Numerical solutions based on a coupled multiphysics finite element scheme are presented. The modeling scheme is derived based on mass conservation including biomolecular mass, momentum balance including stress due to Coulomb force field and DNA-fluid interaction, and charge transport associated to DNA and other ionic complexes in the fluid. Variation in the velocity field for the non-Newtonian flow and the deformation of the DNA strand which results from the fluid-structure interaction are first studied considering a single DNA strand. Motion of the effective center of mass is analyzed considering various straight and coil geometries. Effects of DNA statistical parameters (geometry and spatial distribution of DNAs along the channel) on the effective flow behavior are analyzed. In particular, the dynamics of different DNA physical properties such as radius of gyration, end-to-end length etc. which are obtained from various different models (Kratky-Porod, Gaussian bead-spring etc.) are correlated to the nature of interaction and physical properties under the same background fluid environment.

Partial delamination modeling in composite beams using a finite element method

A new method of modeling partial delamination in composite beams is proposed and implemented using the finite element method. Homogenized cross-sectional stiffness of the delaminated beam is obtained by the proposed analytical technique, including extension-bending, extension-twist and torsion-bending coupling terms, and hence can be used with an existing finite element method. A two noded C1 type Timoshenko beam element with 4 degrees of freedom per node for dynamic analysis of beams is implemented. The results for different delamination scenarios and beams subjected to different boundary conditions are validated with available experimental results in the literature and/or with the 3D finite element simulation using COMSOL. Results of the first torsional mode frequency for the partially delaminated beam are validated with the COMSOL results. The key point of the proposed model is that partial delamination in beams can be analyzed using a beam model, rather than using 3D or plate models. (c) 2013 Elsevier B.V. All rights reserved.

The problem of clear air turbulence: Changing perspectives in the understanding of the phenomenon

Due to rapid improvements in on-board instrumentation and atmospheric observation systems, in most cases, aircraft are able to steer clear of regions of adverse weather. However, they still encounter unexpected bumpy flight conditions in regions away from storms and clouds. This is the phenomenon of clear air turbulence (CAT), which has been a challenge to our understanding as well as efforts at prediction. While most of such cases result in mild discomfort, a few cases can be violent leading to serious injuries to passengers and damage to the aircraft. The underlying physical mechanisms have been sought to be explained in terms of fluid dynamic instabilities and waves in the atmosphere. The main mechanisms which have been proposed are: (i) Kelvin-Helmholtz instability of shear layers, (ii) waves generated from flow over mountains, (iii) inertia-gravity waves from clouds and other sources, (iv) spontaneous imbalance theory and (v) horizontal vortex tubes. This has also undergone a change over the years. We present an overview of the mechanisms proposed and their implications for prediction.

Thermal Efficiency of LPG and PNG-fired burners: Experimental and numerical studies

In the present work, the thermal efficiency of a conventional domestic burner is studied both experimentally and numerically for liquefied petroleum gas (LPG) and piped natural gas (PNG) fuels. Three-dimensional computational fluid dynamic (CFD) modeling of the steady-state flow, combustion and heat transfer to the vessel is reported for the first time in such burners. Based on the insights from the CFD model concerning the flow and heat transfer, design modifications in the form of a circular insert and a radiant sheet are proposed which are observed to increase thermal efficiency for LPG. For PNG, predictions showed that loading height was a much more important factor affecting efficiency than these design modifications and an optimal loading height could be identified. Experiments confirm these trends by showing an improvement in burner thermal efficiency of 2.5% for LPG with the modified design, and 10% for PNG with the optimal loading height, demonstrating that the CFD modeling approach developed in the present work is a useful tool to study domestic burners. (C) 2013 Elsevier Ltd. All rights reserved.

Effect of port gas injection on in-cylinder stratification in a CNG fuelled SI engine

The concept of barrel stratification of air-fuel mixture is evaluated for a port gas injection (PGI) single cylinder spark ignition (SI) internal combustion (IC) engine using a transient three-dimensional computational fluid dynamic (CFD) model. The gaseous fuel used in the study is compressed natural gas (CNG). It is observed that compared to the premixed gas carburettor case, a substantial amount of in-cylinder stratification can be achieved with port gas injection system. A detailed parametric study is reported to understand the effect of the various injection parameters such as injection location, injection orientation, start of injection (SOT) and its duration, and injection rate. Furthermore, the best injection timing is evaluated for various load and speed cases. It is observed that the best stratification pattern can be achieved at 50% engine load. The injection location is observed to have a profound effect on the in-cylinder stratification pattern, and injection towards the side of the spark plug is observed to give a rich fuel-air mixture near the spark plug. It is also shown that there exists an optimal injection pressure.

Flexibility in flapping foil suppresses meandering of induced jet in absence of free stream

Thrust-generating flapping foils are known to produce jets inclined to the free stream at high Strouhal numbers St = fA/U-infinity, where f is the frequency and A is the amplitude of flapping and U-infinity is the free-stream velocity. Our experiments, in the limiting case of St —> infinity (zero free-stream speed), show that a purely oscillatory pitching motion of a chordwise flexible foil produces a coherent jet composed of a reverse Benard-Karman vortex street along the centreline, albeit over a specific range of effective flap stiffnesses. We obtain flexibility by attaching a thin flap to the trailing edge of a rigid NACA0015 foil; length of flap is 0.79 c where c is rigid foil chord length. It is the time-varying deflections of the flexible flap that suppress the meandering found in the jets produced by a pitching rigid foil for zero free-stream condition. Recent experiments (Marais et al., J. Fluid Mech., vol. 710, 2012, p. 659) have also shown that the flexibility increases the St at which non-deflected jets are obtained. Analysing the near-wake vortex dynamics from flow visualization and particle image velocimetry (PIV) measurements, we identify the mechanisms by which flexibility suppresses jet deflection and meandering. A convenient characterization of flap deformation, caused by fluid-flap interaction, is through a non-dimensional effective stiffness', EI* = 8 EI/(rho V-TEmax(2) s(f) c(f)(3)/2), representing the inverse of the flap deflection due to the fluid-dynamic loading; here, EI is the bending stiffness of flap, rho is fluid density, V-TEmax is the maximum velocity of rigid foil trailing edge, s(f) is span and c(f) is chord length of the flexible flap. By varying the amplitude and frequency of pitching, we obtain a variation in EI* over nearly two orders of magnitude and show that only moderate EI*. (0.1 less than or similar to EI * less than or similar to 1 generates a sustained, coherent, orderly jet. Relatively `stiff' flaps (EI* greater than or similar to 1), including the extreme case of no flap, produce meandering jets, whereas highly `flexible' flaps (EI* less than or similar to 0.1) produce spread-out jets. Obtained from the measured mean velocity fields, we present values of thrust coefficients for the cases for which orderly jets are observed.

Multithreaded Test Synthesis for Deadlock Detection

Designing and implementing thread-safe multithreaded libraries can be a daunting task as developers of these libraries need to ensure that their implementations are free from concurrency bugs, including deadlocks. The usual practice involves employing software testing and/or dynamic analysis to detect. deadlocks. Their effectiveness is dependent on well-designed multithreaded test cases. Unsurprisingly, developing multithreaded tests is significantly harder than developing sequential tests for obvious reasons. In this paper, we address the problem of automatically synthesizing multithreaded tests that can induce deadlocks. The key insight to our approach is that a subset of the properties observed when a deadlock manifests in a concurrent execution can also be observed in a single threaded execution. We design a novel, automatic, scalable and directed approach that identifies these properties and synthesizes a deadlock revealing multithreaded test. The input to our approach is the library implementation under consideration and the output is a set of deadlock revealing multithreaded tests. We have implemented our approach as part of a tool, named OMEN1. OMEN is able to synthesize multithreaded tests on many multithreaded Java libraries. Applying a dynamic deadlock detector on the execution of the synthesized tests results in the detection of a number of deadlocks, including 35 real deadlocks in classes documented as thread-safe. Moreover, our experimental results show that dynamic analysis on multithreaded tests that are either synthesized randomly or developed by third-party programmers are ineffective in detecting the deadlocks.

Atmospheric pressure plasma chemical vapor deposition reactor for 100 mm wafers, optimized for minimum contamination at low gas flow rates

Gas discharge plasmas used for thinfilm deposition by plasma-enhanced chemical vapor deposition (PECVD) must be devoid of contaminants, like dust or active species which disturb the intended chemical reaction. In atmospheric pressure plasma systems employing an inert gas, the main source of such contamination is the residual air inside the system. To enable the construction of an atmospheric pressure plasma (APP) system with minimal contamination, we have carried out fluid dynamic simulation of the APP chamber into which an inert gas is injected at different mass flow rates. On the basis of the simulation results, we have designed and built a simple, scaled APP system, which is capable of holding a 100 mm substrate wafer, so that the presence of air (contamination) in the APP chamber is minimized with as low a flow rate of argon as possible. This is examined systematically by examining optical emission from the plasma as a function of inert gas flow rate. It is found that optical emission from the plasma shows the presence of atmospheric air, if the inlet argon flow rate is lowered below 300 sccm. That there is minimal contamination of the APP reactor built here, was verified by conducting an atmospheric pressure PECVD process under acetylene flow, combined with argon flow at 100 sccm and 500 sccm. The deposition of a polymer coating is confirmed by infrared spectroscopy. X-ray photoelectron spectroscopy shows that the polymer coating contains only 5% of oxygen, which is comparable to the oxygen content in polymer deposits obtained in low-pressure PECVD systems. (C) 2015 AIP Publishing LLC.

The hemodynamic effects of in-tandem carotid artery stenosis: implications for carotid endarterectomy.

OBJECTIVES: It remains controversial whether patients with severe disease of the internal carotid artery and a coexisting stenotic lesion downstream would benefit from a carotid endarterectomy (CEA) of the proximal lesion. The aim of this study was to simulate the hemodynamic and wall shear effects of in-tandem internal carotid artery stenosis using a computational fluid dynamic (CFD) idealized model to give insight into the possible consequences of CEA on these lesions. METHODS: A CFD model of steady viscous flow in a rigid tube with two asymmetric stenoses was introduced to simulate blood flow in arteries with multiple constrictions. The effect of varying the distance between the two stenoses, and the severity of the upstream stenosis on the pressure and wall shear stress (WSS) distributions on the second plaque, was investigated. The influence of the relative positions of the two stenoses was also assessed. RESULTS: The distance between the plaques was found to have minimal influence on the overall hemodynamic effect except for the presence of a zone of low WSS (range -20 to 30 dyne/cm2) adjacent to both lesions when the two stenoses were sufficiently close (<4 times the arterial diameter). The upstream stenosis was protective if it was larger than the downstream stenosis. The relative positions of the stenoses were found to influence the WSS but not the pressure distribution. CONCLUSIONS: The geometry and positions of the lesions need to be considered when considering the hemodynamic effects of an in-tandem stenosis. Low WSS is thought to cause endothelial dysfunction and initiate atheroma formation. The fact that there was a flow recirculation zone with low WSS in between the two stenoses may demonstrate how two closely positioned plaques may merge into one larger lesion. Decision making for CEA may need to take into account the hemodynamic situation when an in-tandem stenosis is found. CFD may aid in the risk stratification of patients with this problem.

Governing Equations and Numerical Solutions of Tension Leg Platform with Finite Amplitude Motion

It is demonstrated that when tension leg platform (TLP) moves with finite amplitude in waves, the inertia force, the drag force and the buoyancy acting on the platform are nonlinear functions of the response of TLP. The tensions of the tethers are also nonlinear functions of the displacement of TLP. Then the displacement, the velocity and the acceleration of TLP should be taken into account when loads are calculated. In addition, equations of motions should be set up on the instantaneous position. A theoretical model for analyzing the nonlinear behavior of a TLP with finite displacement is developed, in which multifold nonlinearities are taken into account, i.e., finite displacement, coupling of the six degrees of freedom, instantaneous position, instantaneous wet surface, free surface effects and viscous drag force. Based on the theoretical model, the comprehensive nonlinear differential equations are deduced. Then the nonlinear dynamic analysis of ISSC TLP in regular waves is performed in the time domain. The degenerative linear solution of the proposed nonlinear model is verified with existing published one. Furthermore, numerical results are presented, which illustrate that nonlinearities exert a significant influence on the dynamic responses of the TLP.

High Aspect Ratio (L/D) Riser Viv Prediction Using Wake Oscillator Model

A two-dimensional (2-D) vortex-induced vibration (VIV) prediction model for high aspect ratio (LID) riser subjected to uniform and sheared flow is studied in this paper. The nonlinear structure equations are considered. The near wake dynamics describing the fluctuating nature of vortex shedding is modeled using classical van der Pol equation. A new approach was applied to calibrate the empirical parameters in the wake oscillator model. Compared the predicted results with the experimental data and computational fluid dynamic (CFD) results. Good agreements are observed. It can be concluded that the present model can be used as simple computational tool in predicting some aspects of VIV of long flexible structures. (C) 2008 Elsevier Ltd. All rights reserved.

Self-consistent kinetic modeling of electron distribution function in a low-pressure inductively coupled plasma

Based upon the spatially inhomogeneous Boltzmann equation in two-term approximation coupled with electromagnetic and fluid model analysis for the recently developed inductively coupled plasma sources, a self-consistent electron kinetic model is developed. The electron distribution function, spatial distributions of the electron density and ionization rate are calculated and discussed.

基于变形修正的局部刚体化动力模型简化方法

提出一种基于变形修正的动力模型简化方法,以复杂工程结构的有限元模型为基础,根据运动同步性假设将结构分成若干同步性区域,采用适当的位移模式去逼近位移真值,再通过总体应变能最小来优化,从而得到整体位移与局部位移模式间的转换关系,算例结果表明了简化方法的有效性.

基于柔度修正的局部刚体化结构动力模型简化方法

提出一种基于柔度修正的局部刚体化结构动力模型简化方法.以复杂工程结构的有限元模型为基础,根据运动同步性假设将结构分成若干同步性区域,对每个区域引入准刚体模态,实现对结构的简化.同时对简化方法进行误差分析,利用柔度修正的思路实现对简化变化矩阵的修正,得到更精确的简化结果.以海洋平台结构作为算例应用该方法,得到了比较满意的简化结果.

三峡升船机系统动态响应分析

根据三峡升船机最新设计报告,采用有限元软件ANSYS构建了升船机系统有限元模型,包括承船厢、提升钢索、滑轮组、平衡重、厢内水体、驱动机构和刹车机构;通过数值模拟,对升船机系统进行了模态分析及动态响应研究,内容包括升船机系统的流固耦合模态及其在各种工况(起动、刹车、事故)下的动态响应和影响因素分析,如弹簧刚度、误载水深、阻尼、承船厢的位置等.计算结果表明:对于升船机这样一个复杂系统,该文所建立的有限元模型及数值算法是合理的,其结果可为升船机系统整体设计提供必要的参考.