Gradient-Driven Vortex Motion in Nonneutral Plasmas and Ideal 2D Fluids
NASA Astrophysics Data System (ADS)
Schecter, David A.
2000-10-01
Two-dimensional (2D) turbulent flows can relax to metastable patterns without dissipation of kinetic energy. This ``rapid'' relaxation has been observed in computer simulations of ideal 2D fluids, and more recently in experiments with pure electron plasmas, which can obey similar dynamics. The late stage of relaxation often involves small vortices moving in a larger ``background'' shear-flow.(X.P. Huang et al., Phys. Rev. Lett. 74), 4424 (1995). In time, positive vortices (rotating counter-clockwise) move to peaks in background vorticity, whereas negative vortices (rotating clockwise) move to minima.(C.G. Rossby, J. Mar. Res. 7), 175 (1948); C.H. Liu and L. Ting, Comp. & Fluids 15, 77 (1987). In general, the rate of this migration increases with the magnitude of the background vorticity gradient, whereas it decreases as the background shear intensifies.\\vspace12pt Positive and negative vortices can also be classified as either prograde or retrograde, depending on whether they rotate with or against the local background shear. Surprisingly, a retrograde vortex moves up or down a background vorticity gradient orders of magnitude faster than a prograde vortex of equal strength.(D.A. Schecter and D.H.E. Dubin, Phys. Rev. Lett. 83), 2191 (1999). An accurate expression for the velocity of a weak retrograde vortex is obtained from an analytic calculation, in which the response of the background flow to the vortex is linearized. However, this linear theory fails for prograde vortices of any strength. Interestingly, the velocity of a prograde vortex can be obtained from a simple estimate, which accounts for the nonlinear ``trapping'' of background fluid around the vortex. The analytic expressions for the velocities of both prograde and retrograde vortices are in good quantitative agreement with vortex-in-cell simulations, and with electron plasma experiments, when the background shear is below a critical level. When the ratio of background shear to background vorticity
Rheological Properties of Quasi-2D Fluids in Microgravity
NASA Technical Reports Server (NTRS)
Stannarius, Ralf; Trittel, Torsten; Eremin, Alexey; Harth, Kirsten; Clark, Noel; Maclennan, Joseph; Glaser, Matthew; Park, Cheol; Hall, Nancy; Tin, Padetha
2015-01-01
In recent years, research on complex fluids and fluids in restricted geometries has attracted much attention in the scientific community. This can be attributed not only to the development of novel materials based on complex fluids but also to a variety of important physical phenomena which have barely been explored. One example is the behavior of membranes and thin fluid films, which can be described by two-dimensional (2D) rheology behavior that is quite different from 3D fluids. In this study, we have investigated the rheological properties of freely suspended films of a thermotropic liquid crystal in microgravity experiments. This model system mimics isotropic and anisotropic quasi 2D fluids [46]. We use inkjet printing technology to dispense small droplets (inclusions) onto the film surface. The motion of these inclusions provides information on the rheological properties of the films and allows the study of a variety of flow instabilities. Flat films have been investigated on a sub-orbital rocket flight and curved films (bubbles) have been studied in the ISS project OASIS. Microgravity is essential when the films are curved in order to avoid sedimentation. The experiments yield the mobility of the droplets in the films as well as the mutual mobility of pairs of particles. Experimental results will be presented for 2D-isotropic (smectic-A) and 2D-nematic (smectic-C) phases.
ENERGY LANDSCAPE OF 2D FLUID FORMS
Y. JIANG; ET AL
2000-04-01
The equilibrium states of 2D non-coarsening fluid foams, which consist of bubbles with fixed areas, correspond to local minima of the total perimeter. (1) The authors find an approximate value of the global minimum, and determine directly from an image how far a foam is from its ground state. (2) For (small) area disorder, small bubbles tend to sort inwards and large bubbles outwards. (3) Topological charges of the same sign repel while charges of opposite sign attract. (4) They discuss boundary conditions and the uniqueness of the pattern for fixed topology.
Collective excitations in 2D hard-disc fluid.
Huerta, Adrian; Bryk, Taras; Trokhymchuk, Andrij
2015-07-01
Collective dynamics of a two-dimensional (2D) hard-disc fluid was studied by molecular dynamics simulations in the range of packing fractions that covers states up to the freezing. Some striking features concerning collective excitations in this system were observed. In particular, the short-wavelength shear waves while being absent at low packing fractions were observed in the range of high packing fractions, just before the freezing transition in a 2D hard-disc fluid. In contrast, the so-called "positive sound dispersion" typically observed in dense Lennard-Jones-like fluids, was not detected for the 2D hard-disc fluid. The ratio of specific heats in the 2D hard-disc fluid shows a monotonic increase with density approaching the freezing, resembling in this way the similar behavior in the vicinity of the Widom line in the case of supercritical fluids. PMID:25595625
A new inversion method for (T2, D) 2D NMR logging and fluid typing
NASA Astrophysics Data System (ADS)
Tan, Maojin; Zou, Youlong; Zhou, Cancan
2013-02-01
One-dimensional nuclear magnetic resonance (1D NMR) logging technology has some significant limitations in fluid typing. However, not only can two-dimensional nuclear magnetic resonance (2D NMR) provide some accurate porosity parameters, but it can also identify fluids more accurately than 1D NMR. In this paper, based on the relaxation mechanism of (T2, D) 2D NMR in a gradient magnetic field, a hybrid inversion method that combines least-squares-based QR decomposition (LSQR) and truncated singular value decomposition (TSVD) is examined in the 2D NMR inversion of various fluid models. The forward modeling and inversion tests are performed in detail with different acquisition parameters, such as magnetic field gradients (G) and echo spacing (TE) groups. The simulated results are discussed and described in detail, the influence of the above-mentioned observation parameters on the inversion accuracy is investigated and analyzed, and the observation parameters in multi-TE activation are optimized. Furthermore, the hybrid inversion can be applied to quantitatively determine the fluid saturation. To study the effects of noise level on the hybrid method and inversion results, the numerical simulation experiments are performed using different signal-to-noise-ratios (SNRs), and the effect of different SNRs on fluid typing using three fluid models are discussed and analyzed in detail.
A discrete simulation of 2-D fluid flow on TERASYS
Mullins, P.G.; Krolak, P.D.
1995-12-01
A discrete simulation of two-dimensional (2-D) fluid flow, on a recently designed novel architecture called TERASYS is presented. The simulation uses a cellular automaton approach, implemented in a new language called data-parallel bit C (dbC). A performance comparison between our implementation on TERASYS and an implementation on the Connection Machine is discussed. We comment briefly on the suitability of the TERASYS system for modeling fluid flow using cellular automata.
2-D linear motion system. Innovative technology summary report
1998-11-01
The US Department of Energy's (DOE's) nuclear facility decontamination and decommissioning (D and D) program requires buildings to be decontaminated, decommissioned, and surveyed for radiological contamination in an expeditious and cost-effective manner. Simultaneously, the health and safety of personnel involved in the D and D activities is of primary concern. D and D workers must perform duties high off the ground, requiring the use of manlifts or scaffolding, often, in radiologically or chemically contaminated areas or in areas with limited access. Survey and decontamination instruments that are used are sometimes heavy or awkward to use, particularly when the worker is operating from a manlift or scaffolding. Finding alternative methods of performing such work on manlifts or scaffolding is important. The 2-D Linear Motion System (2-D LMS), also known as the Wall Walker{trademark}, is designed to remotely position tools and instruments on walls for use in such activities as radiation surveys, decontamination, and painting. Traditional (baseline) methods for operating equipment for these tasks require workers to perform duties on elevated platforms, sometimes several meters above the ground surface and near potential sources of contamination. The Wall Walker 2-D LMS significantly improves health and safety conditions by facilitating remote operation of equipment. The Wall Walker 2-D LMS performed well in a demonstration of its precision, accuracy, maneuverability, payload capacity, and ease of use. Thus, this innovative technology is demonstrated to be a viable alternative to standard methods of performing work on large, high walls, especially those that have potential contamination concerns. The Wall Walker was used to perform a final release radiological survey on over 167 m{sup 2} of walls. In this application, surveying using a traditional (baseline) method that employs an aerial lift for manual access was 64% of the total cost of the improved technology
A 2D electrohydrodynamic model for electrorotation of fluid drops.
Feng, James Q
2002-02-01
A theoretical analysis of spontaneous electrorotation of deformable fluid drops in a DC electric field is presented with a 2D electrohydrodynamic model. The fluids in the system are assumed to be leaky dielectric and Newtonian. If the rotating flow is dominant over the cellular convection type of electrohydrodynamic flow, closed-form solutions for drops of small deformations can be obtained. Because the governing equations are in general nonlinear even when drop deformations are ignored, the general solution for even undeformed drop takes a form of infinite series and can only be evaluated by numerical means. Both closed-form solutions for special cases and numerical solutions for more general cases are obtained here to describe steady-state field variables and first-order drop deformations. In a DC electric field of strength beyond the threshold value, spontaneous electrorotation of a drop is shown to occur when charge relaxation in the surrounding fluid is faster than the fluid inside the drop. With increasing the strength of the applied electric field from the threshold for onset of electrorotation, the axis of drop contraction deviates from from that of the applied electric field in the direction of the rotating flow with an angle increasing with the field strength. PMID:16290391
Validation of a 2-D semi-coupled numerical model for fluid-structure-seabed interaction
NASA Astrophysics Data System (ADS)
Ye, Jianhong; Jeng, Dongsheng; Wang, Ren; Zhu, Changqi
2013-10-01
A 2-D semi-coupled model PORO-WSSI 2D (also be referred as FSSI-CAS 2D) for the Fluid-Structure-Seabed Interaction (FSSI) has been developed by employing RANS equations for wave motion in fluid domain, VARANS equations for porous flow in porous structures; and taking the dynamic Biot's equations (known as "u - p" approximation) for soil as the governing equations. The finite difference two-step projection method and the forward time difference method are adopted to solve the RANS, VARANS equations; and the finite element method is adopted to solve the "u - p" approximation. A data exchange port is developed to couple the RANS, VARANS equations and the dynamic Biot's equations together. The analytical solution proposed by Hsu and Jeng (1994) and some experiments conducted in wave flume or geotechnical centrifuge in which various waves involved are used to validate the developed semi-coupled numerical model. The sandy bed involved in these experiments is poro-elastic or poro-elastoplastic. The inclusion of the interaction between fluid, marine structures and poro-elastoplastic seabed foundation is a special point and highlight in this paper, which is essentially different with other previous coupled models The excellent agreement between the numerical results and the experiment data indicates that the developed coupled model is highly reliablefor the FSSI problem.
NASA Astrophysics Data System (ADS)
Wang, Zenghui; Feng, Philip X.-L.
2016-07-01
Atomic layer crystals are emerging building blocks for enabling new two-dimensional (2D) nanomechanical systems, whose motions can be coupled to other attractive physical properties in such 2D systems. Optical interferometry has been very effective in reading out the infinitesimal motions of these 2D structures and spatially resolving different modes. To quantitatively understand the detection efficiency and its dependence on the device parameters and interferometric conditions, here we present a systematic study of the intrinsic motion responsivity in 2D nanomechanical systems using a Fresnel-law-based model. We find that in monolayer to 14-layer structures, MoS2 offers the highest responsivity among graphene, h-BN, and MoS2 devices and for the three commonly used visible laser wavelengths (633, 532, and 405 nm). We also find that the vacuum gap resulting from the widely used 300 nm-oxide substrate in making 2D devices, fortunately, leads to close-to-optimal responsivity for a wide range of 2D flakes. Our results elucidate and graphically visualize the dependence of motion transduction responsivity upon 2D material type and number of layers, vacuum gap, oxide thickness, and detecting wavelength, thus providing design guidelines for constructing 2D nanomechanical systems with optimal optical motion readout.
Wang, Zenghui; Feng, Philip X-L
2016-01-01
Atomic layer crystals are emerging building blocks for enabling new two-dimensional (2D) nanomechanical systems, whose motions can be coupled to other attractive physical properties in such 2D systems. Optical interferometry has been very effective in reading out the infinitesimal motions of these 2D structures and spatially resolving different modes. To quantitatively understand the detection efficiency and its dependence on the device parameters and interferometric conditions, here we present a systematic study of the intrinsic motion responsivity in 2D nanomechanical systems using a Fresnel-law-based model. We find that in monolayer to 14-layer structures, MoS2 offers the highest responsivity among graphene, h-BN, and MoS2 devices and for the three commonly used visible laser wavelengths (633, 532, and 405 nm). We also find that the vacuum gap resulting from the widely used 300 nm-oxide substrate in making 2D devices, fortunately, leads to close-to-optimal responsivity for a wide range of 2D flakes. Our results elucidate and graphically visualize the dependence of motion transduction responsivity upon 2D material type and number of layers, vacuum gap, oxide thickness, and detecting wavelength, thus providing design guidelines for constructing 2D nanomechanical systems with optimal optical motion readout. PMID:27464908
Wang, Zenghui; Feng, Philip X.-L.
2016-01-01
Atomic layer crystals are emerging building blocks for enabling new two-dimensional (2D) nanomechanical systems, whose motions can be coupled to other attractive physical properties in such 2D systems. Optical interferometry has been very effective in reading out the infinitesimal motions of these 2D structures and spatially resolving different modes. To quantitatively understand the detection efficiency and its dependence on the device parameters and interferometric conditions, here we present a systematic study of the intrinsic motion responsivity in 2D nanomechanical systems using a Fresnel-law-based model. We find that in monolayer to 14-layer structures, MoS2 offers the highest responsivity among graphene, h-BN, and MoS2 devices and for the three commonly used visible laser wavelengths (633, 532, and 405 nm). We also find that the vacuum gap resulting from the widely used 300 nm-oxide substrate in making 2D devices, fortunately, leads to close-to-optimal responsivity for a wide range of 2D flakes. Our results elucidate and graphically visualize the dependence of motion transduction responsivity upon 2D material type and number of layers, vacuum gap, oxide thickness, and detecting wavelength, thus providing design guidelines for constructing 2D nanomechanical systems with optimal optical motion readout. PMID:27464908
Collective motion of microswimmers in viscoelastic fluids
NASA Astrophysics Data System (ADS)
Li, Gaojin; Ardekani, Arezoo
2015-11-01
The dynamics of suspension of self-propelled microorganisms show fascinating hydrodynamic phenomena, such as, large scale swarming motion, locally correlated motion, enhanced particle diffusion, and enhanced fluid mixing. Even though many studies have been conducted in a Newtonian fluid, the collective motion of microorganisms in non-Newtonian fluids is less understood. The non-Newtonian fluid rheological properties, such as viscoelasticity and shear-dependent viscosity in saliva, mucus and biofilm, significantly affect the swimming properties and hydrodynamic interaction of microorganisms. In this work, we use direct numerical simulation to investigate the collective motion of rod-like swimmers in viscoelastic fluids. Two swimming types, pusher and puller, are investigated. The background viscoelastic fluid is modeled using an Oldroyd-B constitutive equation. This work is supported by NSF CBET-1445955 and Indiana CTSI TR001108.
2-D traveling-wave patterns in binary fluid convection
Surko, C.M.; Porta, A.L.
1996-12-31
An overview is presented of recent experiments designed to study two-dimensional traveling-wave convection in binary fluid convection in a large aspect ratio container. Disordered patterns are observed when convection is initiated. As time proceeds, they evolve to more ordered patterns, consisting of several domains of traveling-waves separated by well-defined domain boundaries. The detailed character of the patterns depends sensitively on the Rayleigh number. Numerical techniques are described which were developed to provide a quantitative characterization of the traveling-wave patterns. Applications of complex demodulation techniques are also described, which make a detailed study of the structure and dynamics of the domain boundaries possible.
NASA Astrophysics Data System (ADS)
Olmez, O.; Ozbulut, M.; Yildiz, M.; Goren, O.
2016-06-01
The present study investigates the vortical and nonlinear effects in the roll motion of a 2-D body with square cross-sections by using Smoothed Particle Hydrodynamics (SPH). A 2-D rigid body with square cross-section is taken into account for the benchmark study and subjected to the oscillatory roll motion with a given angular frequency. The governing equations are continuity equation and Euler's equation with artificial viscosity term. Weakly Compressible SPH (WCSPH) scheme is employed for the discretization of the governing equations. Velocities of the fluid particles are updated by means of XSPH+Artificial Particle Displacement (VXSPH+APD) algorithm. In this method only the free surface fluid particles are subjected to VXSPH algorithm while the APD algorithm is employed for the fully populated flow regions. The hybrid usage of numerical treatment keeps free surface particles together by creating an artificial surface tension on the free surface. VXSPH+APD is a proven numerical treatment to provide the most accurate results for this type of free surface flows (Ozbulut et al. 2014). The results of the present study are compared with those of the experimental studies as well as with those of the numerical methods obtained from the current literature.
Edge preserving motion estimation with occlusions correction for assisted 2D to 3D conversion
NASA Astrophysics Data System (ADS)
Pohl, Petr; Sirotenko, Michael; Tolstaya, Ekaterina; Bucha, Victor
2014-02-01
In this article we propose high quality motion estimation based on variational optical flow formulation with non-local regularization term. To improve motion in occlusion areas we introduce occlusion motion inpainting based on 3-frame motion clustering. Variational formulation of optical flow proved itself to be very successful, however a global optimization of cost function can be time consuming. To achieve acceptable computation times we adapted the algorithm that optimizes convex function in coarse-to-fine pyramid strategy and is suitable for modern GPU hardware implementation. We also introduced two simplifications of cost function that significantly decrease computation time with acceptable decrease of quality. For motion clustering based motion inpaitning in occlusion areas we introduce effective method of occlusion aware joint 3-frame motion clustering using RANSAC algorithm. Occlusion areas are inpainted by motion model taken from cluster that shows consistency in opposite direction. We tested our algorithm on Middlebury optical flow benchmark, where we scored around 20th position, but being one of the fastest method near the top. We also successfully used this algorithm in semi-automatic 2D to 3D conversion tool for spatio-temporal background inpainting, automatic adaptive key frame detection and key points tracking.
In situ fluid typing and quantification with 1D and 2D NMR logging.
Sun, Boqin
2007-05-01
In situ nuclear magnetic resonance (NMR) fluid typing has recently gained momentum due to data acquisition and inversion algorithm enhancement of NMR logging tools. T(2) distributions derived from NMR logging contain information on bulk fluids and pore size distributions. However, the accuracy of fluid typing is greatly overshadowed by the overlap between T(2) peaks arising from different fluids with similar apparent T(2) relaxation times. Nevertheless, the shapes of T(2) distributions from different fluid components are often different and can be predetermined. Inversion with predetermined T(2) distributions allows us to perform fluid component decomposition to yield individual fluid volume ratios. Another effective method for in situ fluid typing is two-dimensional (2D) NMR logging, which results in proton population distribution as a function of T(2) relaxation time and fluid diffusion coefficient (or T(1) relaxation time). Since diffusion coefficients (or T(1) relaxation time) for different fluid components can be very different, it is relatively easy to separate oil (especially heavy oil) from water signal in a 2D NMR map and to perform accurate fluid typing. Combining NMR logging with resistivity and/or neutron/density logs provides a third method for in situ fluid typing. We shall describe these techniques with field examples. PMID:17466778
Rupture dynamics and ground motions from earthquakes in 2-D heterogeneous media
NASA Astrophysics Data System (ADS)
Bydlon, Samuel A.; Dunham, Eric M.
2015-03-01
We perform 2-D simulations of earthquakes on rough faults in media with random heterogeneities (with von Karman distribution) to study the effects of geometric and material heterogeneity on the rupture process and resulting high-frequency ground motions in the near-fault region (out to ˜20 km). Variations in slip and rupture velocity can arise from material heterogeneity alone but are dominantly controlled by fault roughness. Scattering effects become appreciable beyond ˜3 km from the fault. Near-fault scattering extends the duration of incoherent, high-frequency ground motions and, at least in our 2-D simulations, elevates root-mean-square accelerations (i.e., Arias intensity) with negligible reduction in peak velocities. We also demonstrate that near-fault scattering typically occurs in the power law tail of the power spectral density function, quantified by the Hurst exponent and another parameter combining standard deviation and correlation length.
Numerical simulation of ( T 2, T 1) 2D NMR and fluid responses
NASA Astrophysics Data System (ADS)
Tan, Mao-Jin; Zou, You-Long; Zhang, Jin-Yan; Zhao, Xin
2012-12-01
One-dimensional nuclear magnetic resonance (1D NMR) logging technology is limited for fluid typing, while two-dimensional nuclear magnetic resonance (2D NMR) logging can provide more parameters including longitudinal relaxation time ( T 1) and transverse relaxation time ( T 2) relative to fluid types in porous media. Based on the 2D NMR relaxation mechanism in a gradient magnetic field, echo train simulation and 2D NMR inversion are discussed in detail. For 2D NMR inversion, a hybrid inversion method is proposed based on the damping least squares method (LSQR) and an improved truncated singular value decomposition (TSVD) algorithm. A series of spin echoes are first simulated with multiple waiting times ( T W s) in a gradient magnetic field for given fluid models and these synthesized echo trains are inverted by the hybrid method. The inversion results are consistent with given models. Moreover, the numerical simulation of various fluid models such as the gas-water, light oil-water, and vicious oil-water models were carried out with different echo spacings ( T E s) and T W s by this hybrid method. Finally, the influences of different signal-to-noise ratios (SNRs) on inversion results in various fluid models are studied. The numerical simulations show that the hybrid method and optimized observation parameters are applicable to fluid typing of gas-water and oil-water models.
Nonrigid 2D registration of fluoroscopic coronary artery image sequence with layered motion
NASA Astrophysics Data System (ADS)
Park, Taewoo; Jung, Hoyup; Yun, Il Dong
2016-03-01
We present a new method for nonrigid registration of coronary artery models with layered motion information. 2D nonrigid registration method is proposed that brings layered motion information into correspondence with fluoroscopic angiograms. The registered model is overlaid on top of interventional angiograms to provide surgical assistance during image-guided chronic total occlusion procedures. The proposed methodology is divided into two parts: layered structures alignments and local nonrigid registration. In the first part, inpainting method is used to estimate a layered rigid transformation that aligns layered motion information. In the second part, a nonrigid registration method is implemented and used to compensate for any local shape discrepancy. Experimental evaluation conducted on a set of 7 fluoroscopic angiograms results in a reduced target registration error, which showed the effectiveness of the proposed method over single layered approach.
Principles of the motion of fluids
NASA Astrophysics Data System (ADS)
Euler, Leonhard
2008-08-01
The elements of the theory of the motion of fluids in general are treated here, the whole matter being reduced to this: given a mass of fluid, either free or confined in vessels, upon which an arbitrary motion is impressed, and which in turn is acted upon by arbitrary forces, to determine the motion carrying forward each particle, and at the same time to ascertain the pressure exerted by each part, acting on it as well as on the sides of the vessel. At first in this memoir, before undertaking the investigation of these effects of the forces, the Most Famous Author carefully evaluates all the possible motions which can actually take place in the fluid. Indeed, even if the individual particles of the fluid are free from each other, motions in which the particles interpenetrate are nevertheless excluded, since we are dealing with fluids that do not permit any compression into a narrower volume. Thus it is clear that an arbitrary small portion of fluid cannot receive a motion other than the one which constantly conserves the same volume; even though meanwhile the shape is changed in any way. It would hold indeed, as long as no elementary portion would be compressed at any time into a smaller volume; furthermore if the portion expanded into a larger volume, the continuity of the particles was violated, these were dispersed and no longer clung together, such a motion would no longer pertain to the science of the motion of fluids; but individual droplets would separately perform their motion. Therefore, this case being excluded, motion of the fluids must be restricted by this rule that each small portion must retain for ever the same volume; and this principle restricts the general expressions of motion for elements of the fluid. Plainly, considering an arbitrary small portion of the fluid, its individual points have to be carried by such a motion that, when at a moment of time they arrive at the next location, until then they occupy a volume equal to the previous one
The use of 2D ultrasound elastography for measuring tendon motion and strain.
Chernak Slane, Laura; Thelen, Darryl G
2014-02-01
The goal of the current study was to investigate the fidelity of a 2D ultrasound elastography method for the measurement of tendon motion and strain. Ultrasound phantoms and ex vivo porcine flexor tendons were cyclically stretched to 4% strain while cine ultrasound radiofrequency (RF) data and video data were simultaneously collected. 2D ultrasound elastography was used to estimate tissue motion and strain from RF data, and surface tissue motion and strain were separately estimated using digital image correlation (DIC). There were strong correlations (R(2)>0.97) between DIC and RF measurements of phantom displacement and strain, and good agreement in estimates of peak phantom strain (DIC: 3.5±0.2%; RF: 3.7±0.1%). For tendon, elastographic estimates of displacement profiles also correlated well with DIC measurements (R(2)>0.92), and exhibited similar estimated peak tendon strain (DIC: 2.6±1.4%; RF: 2.2±1.3%). Elastographic tracking with B-Mode images tended to under-predict peak strain for both the phantom and tendon. This study demonstrates the capacity to use quantitative elastographic techniques to measure tendon displacement and strain within an ultrasound image window. The approach may be extendible to in vivo use on humans, which would allow for the non-invasive analysis of tendon deformation in both normal and pathological states. PMID:24388164
The Use of 2D Ultrasound Elastography for Measuring Tendon Motion and Strain
Slane, Laura Chernak; Thelen, Darryl G.
2014-01-01
The goal of the current study was to investigate the fidelity of a 2D ultrasound elastography method for the measurement of tendon motion and strain. Ultrasound phantoms and ex vivo porcine flexor tendons were cyclically stretched to 4% strain while cine ultrasound radiofrequency (RF) data and video data were simultaneously collected. 2D ultrasound elastography was used to estimate tissue motion and strain from RF data, and surface tissue motion and strain were separately estimated using digital image correlation (DIC). There were strong correlations (R2 > 0.97) between DIC and RF measurements of phantom displacement and strain, and good agreement in estimates of peak phantom strain (DIC: 3.5 ± 0.2%; RF: 3.7 ± 0.1%). For tendon, elastographic estimates of displacement profiles also correlated well with DIC measurements (R2 > 0.92), and exhibited similar estimated peak tendon strain (DIC: 2.6 ± 1.4%; RF: 2.2 ± 1.3%). Elastographic tracking with B-Mode images tended to under-predict peak strain for both the phantom and tendon. This study demonstrates the capacity to use quantitative elastographic techniques to measure tendon displacement and strain within an ultrasound image window. The approach may be extendible to in vivo use on humans, which would allow for the non-invasive analysis of tendon deformation in both normal and pathological states. PMID:24388164
Coronary arteries motion modeling on 2D x-ray images
NASA Astrophysics Data System (ADS)
Gao, Yang; Sundar, Hari
2012-02-01
During interventional procedures, 3D imaging modalities like CT and MRI are not commonly used due to interference with the surgery and radiation exposure concerns. Therefore, real-time information is usually limited and building models of cardiac motion are difficult. In such case, vessel motion modeling based on 2-D angiography images become indispensable. Due to issues with existing vessel segmentation algorithms and the lack of contrast in occluded vessels, manual segmentation of certain branches is usually necessary. In addition, such occluded branches are the most important vessels during coronary interventions and obtaining motion models for these can greatly help in reducing the procedure time and radiation exposure. Segmenting different cardiac phases independently does not guarantee temporal consistency and is not efficient for occluded branches required manual segmentation. In this paper, we propose a coronary motion modeling system which extracts the coronary tree for every cardiac phase, maintaining the segmentation by tracking the coronary tree during the cardiac cycle. It is able to map every frame to the specific cardiac phase, thereby inferring the shape information of the coronary arteries using the model corresponding to its phase. Our experiments show that our motion modeling system can achieve promising results with real-time performance.
Brownian motion of particles in nematic fluids
NASA Astrophysics Data System (ADS)
Yao, Xuxia; Nayani, Karthik; Park, Jung; Srinivasarao, Mohan
2011-03-01
We studied the brownian motion of both charged and neutral polystyrene particles in two nematic fluids, a thermotropic liquid crystal, E7, and a lyotropic chromonic liquid crystal, Sunset Yellow FCF (SSY). Homogeneous planar alignment of E7 was easliy achieved by using rubbed polyimide film coated on the glass. For SSY planar mondomain, we used the capillary method recently developed in our lab. By tracking a single particle, the direction dependent diffussion coefficients and Stokes drag were measured in the nematic phase and isotropic phase for both systems.
Chen, Chia-Hsiung; Azari, David; Hu, Yu Hen; Lindstrom, Mary J.; Thelen, Darryl; Yen, Thomas Y.; Radwin, Robert G.
2015-01-01
Objective Marker-less 2D video tracking was studied as a practical means to measure upper limb kinematics for ergonomics evaluations. Background Hand activity level (HAL) can be estimated from speed and duty cycle. Accuracy was measured using a cross correlation template-matching algorithm for tracking a region of interest on the upper extremities. Methods Ten participants performed a paced load transfer task while varying HAL (2, 4, and 5) and load (2.2 N, 8.9 N and 17.8 N). Speed and acceleration measured from 2D video were compared against ground truth measurements using 3D infrared motion capture. Results The median absolute difference between 2D video and 3D motion capture was 86.5 mm/s for speed, and 591 mm/s2 for acceleration, and less than 93 mm/s for speed and 656 mm/s2 for acceleration when camera pan and tilt were within ±30 degrees. Conclusion Single-camera 2D video had sufficient accuracy (< 100 mm/s) for evaluating HAL. Practitioner Summary This study demonstrated that 2D video tracking had sufficient accuracy to measure HAL for ascertaining the American Conference of Government Industrial Hygienists Threshold Limit Value® for repetitive motion when the camera is located within ±30 degrees off the plane of motion when compared against 3D motion capture for a simulated repetitive motion task. PMID:25978764
Chen, Chia-Hsiung; Azari, David P; Hu, Yu Hen; Lindstrom, Mary J; Thelen, Darryl; Yen, Thomas Y; Radwin, Robert G
2015-01-01
Marker-less 2D video tracking was studied as a practical means to measure upper limb kinematics for ergonomics evaluations. Hand activity level (HAL) can be estimated from speed and duty cycle. Accuracy was measured using a cross-correlation template-matching algorithm for tracking a region of interest on the upper extremities. Ten participants performed a paced load transfer task while varying HAL (2, 4, and 5) and load (2.2 N, 8.9 N and 17.8 N). Speed and acceleration measured from 2D video were compared against ground truth measurements using 3D infrared motion capture. The median absolute difference between 2D video and 3D motion capture was 86.5 mm/s for speed, and 591 mm/s(2) for acceleration, and less than 93 mm/s for speed and 656 mm/s(2) for acceleration when camera pan and tilt were within ± 30 degrees. Single-camera 2D video had sufficient accuracy (< 100 mm/s) for evaluating HAL. Practitioner Summary: This study demonstrated that 2D video tracking had sufficient accuracy to measure HAL for ascertaining the American Conference of Government Industrial Hygienists Threshold Limit Value(®) for repetitive motion when the camera is located within ± 30 degrees off the plane of motion when compared against 3D motion capture for a simulated repetitive motion task. PMID:25978764
Chang, F.H.; Santee, G.E. Jr.; Mortensen, G.A.; Brockett, G.F.; Gross, M.B.; Silling, S.A.; Belytschko, T.
1981-03-01
This report, the second in a series of reports for RP-1065, describes the second step in the stepwise approach for developing the three-dimensional, nonlinear, fluid/structure interaction methodology to assess the hydroloads on a large PWR during the subcooled portions of a hypothetical LOCA. The second step in the methodology considers enhancements and special modifications to the 2D STEALTH-HYDRO computer program and the 2D WHAMSE computer program. The 2D STEALTH-HYDRO enhancements consist of a fluid-fluid coupling control-volume model and an orifice control-volume model. The enhancements to 2D WHAMSE include elimination of the implicit integration routines, material models, and structural elements not required for the hydroloads application. In addition the logic for coupling the 2D STEALTH-HYDRO computer program to the 2D WHAMSE computer program is discussed.
An Integrative Model of Excitation Driven Fluid Flow in a 2D Uterine Channel
NASA Astrophysics Data System (ADS)
Maggio, Charles; Fauci, Lisa; Chrispell, John
2009-11-01
We present a model of intra-uterine fluid flow in a sagittal cross-section of the uterus by inducing peristalsis in a 2D channel. This is an integrative multiscale computational model that takes as input fluid viscosity, passive tissue properties of the uterine channel and a prescribed wave of membrane depolarization. This voltage pulse is coupled to a model of calcium dynamics inside a uterine smooth muscle cell, which in turn drives a kinetic model of myosin phosphorylation governing contractile muscle forces. Using the immersed boundary method, these muscle forces are communicated to a fluid domain to simulate the contractions which occur in a human uterus. An analysis of the effects of model parameters on the flow properties and emergent geometry of the peristaltic channel will be presented.
The oscillatory motion of a surfactant-laden liquid plug in a 2D-channel
NASA Astrophysics Data System (ADS)
Fujioka, Hideki; Grotberg, James B.
2004-11-01
Liquid plugs can form in the lung's small airways near the end of expiration. This happens more frequently when the amount of pulmonary surfactant is reduced. In medical treatments such as surfactant replacement therapy, partial liquid ventilation, and drug delivery, the formation of plugs in an airway is important to deliver the instilled liquid uniformly throughout the lung. In this study, we investigate numerically the oscillatory motion of a surfactant-laden liquid plug within a two-dimensional channel lined by a thin liquid film. The viscosity of both the left and right air phases is assumed to be negligible, so that the only fluid dynamics of the liquid phase is considered. The plug motion is regulated by the flow rate in the left air phase, which is prescribed as a sinusoidal function of time. The pressure drop between the left and right air phases varies for time with a different phase of the flow rate. The plug length and the film thickness oscillate with an average value during a cycle. These behaviors changes by system parameters, Reynolds number, Womersley number, Capillary number, and surfactant properties. The significance of this study on mechanical stresses acting on airway epithelial cells caused by the motion of a liquid plug during normal breath, conventional or high-frequency ventilation is discussed. Supported by NIH grant HL41126, NASA grant NAG3-2740.
Roton Excitations and the Fluid-Solid Phase Transition in Superfluid 2D Yukawa Bosons
NASA Astrophysics Data System (ADS)
Molinelli, S.; Galli, D. E.; Reatto, L.; Motta, M.
2016-05-01
We compute several ground-state properties and the dynamical structure factor of a zero-temperature system of Bosons interacting with the 2D screened Coulomb (2D-SC) potential. We resort to the exact shadow path integral ground state (SPIGS) quantum Monte Carlo method to compute the imaginary-time correlation function of the model, and to the genetic algorithm via falsification of theories (GIFT) to retrieve the dynamical structure factor. We provide a detailed comparison of ground-state properties and collective excitations of 2D-SC and ^4 He atoms. The roton energy of the 2D-SC system is an increasing function of density, and not a decreasing one as in ^4 He. This result is in contrast with the view that the roton is the soft mode of the fluid-solid transition. We uncover a remarkable quasi-universality of backflow and of other properties when expressed in terms of the amount of short-range order as quantified by the height of the first peak of the static structure factor.
NASA Astrophysics Data System (ADS)
Shaari, M. F.; Abu Bakar, H.; Nordin, N.; Saw, S. K.; Samad, Z.
2013-12-01
Contractile body is an alternative mechanism instead of rotating blade propeller to generate water jet for locomotion. The oscillating motion of the actuator at different frequencies varies the pressure and volume of the pressure chamber in time to draw in and jet out the water at a certain mass flow rate. The aim of this research was to analyze the influence of the actuating frequency of the fluid flow in the pressure chamber of the thruster during this inflation-deflation process. A 70mm × 70mm × 18mm (L × W × T) 2D water jet thruster was fabricated for this purpose. The contractile function was driven using two lateral pneumatic actuators where the fluid flow analysis was focused on the X-Y plane vector. Observation was carried out using a video camera and Matlab image measurement technique to determine the volume of the flowing mass. The result demonstrated that the greater actuating frequency decreases the fluid flow rate and the Reynolds number. This observation shows that the higher frequency would give a higher mass flow rate during water jet generation.
Computational Fluid Dynamics Demonstration of Rigid Bodies in Motion
NASA Technical Reports Server (NTRS)
Camarena, Ernesto; Vu, Bruce T.
2011-01-01
The Design Analysis Branch (NE-Ml) at the Kennedy Space Center has not had the ability to accurately couple Rigid Body Dynamics (RBD) and Computational Fluid Dynamics (CFD). OVERFLOW-D is a flow solver that has been developed by NASA to have the capability to analyze and simulate dynamic motions with up to six Degrees of Freedom (6-DOF). Two simulations were prepared over the course of the internship to demonstrate 6DOF motion of rigid bodies under aerodynamic loading. The geometries in the simulations were based on a conceptual Space Launch System (SLS). The first simulation that was prepared and computed was the motion of a Solid Rocket Booster (SRB) as it separates from its core stage. To reduce computational time during the development of the simulation, only half of the physical domain with respect to the symmetry plane was simulated. Then a full solution was prepared and computed. The second simulation was a model of the SLS as it departs from a launch pad under a 20 knot crosswind. This simulation was reduced to Two Dimensions (2D) to reduce both preparation and computation time. By allowing 2-DOF for translations and 1-DOF for rotation, the simulation predicted unrealistic rotation. The simulation was then constrained to only allow translations.
NASA Astrophysics Data System (ADS)
Choquard, Ph.; Vuffray, M.
2014-10-01
The coupling between dilatation and vorticity, two coexisting and fundamental processes in fluid dynamics (Wu et al., 2006, pp. 3, 6) is investigated here, in the simplest cases of inviscid 2D isotropic Burgers and pressureless Euler-Coriolis fluids respectively modeled by single vortices confined in compressible, local, inertial and global, rotating, environments. The field equations are established, inductively, starting from the equations of the characteristics solved with an initial Helmholtz decomposition of the velocity fields namely a vorticity free and a divergence free part (Wu et al., 2006, Sects. 2.3.2, 2.3.3) and, deductively, by means of a canonical Hamiltonian Clebsch like formalism (Clebsch, 1857, 1859), implying two pairs of conjugate variables. Two vector valued fields are constants of the motion: the velocity field in the Burgers case and the momentum field per unit mass in the Euler-Coriolis one. Taking advantage of this property, a class of solutions for the mass densities of the fluids is given by the Jacobian of their sum with respect to the actual coordinates. Implementation of the isotropy hypothesis entails a radial dependence of the velocity potentials and of the stream functions associated to the compressible and to the rotational part of the fluids and results in the cancellation of the dilatation-rotational cross terms in the Jacobian. A simple expression is obtained for all the radially symmetric Jacobians occurring in the theory. Representative examples of regular and singular solutions are shown and the competition between dilatation and vorticity is illustrated. Inspired by thermodynamical, mean field theoretical analogies, a genuine variational formula is proposed which yields unique measure solutions for the radially symmetric fluid densities investigated. We stress that this variational formula, unlike the Hopf-Lax formula, enables us to treat systems which are both compressible and rotational. Moreover in the one
Turbulent solutions of equations of fluid motion
NASA Technical Reports Server (NTRS)
Deissler, R. G.
1985-01-01
Some turbulent solutions of the unaveraged Navier-Stokes equations (equations of fluid motion) are reviewed. Those equations are solved numerically in order to study the nonlinear physics of incompressible turbulent flow. The three components of the mean-square velocity fluctuations are initially equal for the conditions chosen. The resulting solutions show characteristics of turbulence, such as the linear and nonlinear excitation of small-scale fluctuations. For the stronger fluctuations the initially nonrandom flow develops into an apparently random turbulence. The cases considered include turbulence that is statistically homogeneous or inhomogeneous and isotropic or anisotropic. A statistically steady-state turbulence is obtained by using a spatially periodic body force. Various turbulence processes, including the transfer of energy between eddy sizes and between directional components and the production, dissipation, and spatial diffusion of turbulence, are considered. It is concluded that the physical processes occurring in turbulence can be profitably studied numerically.
An experimental method for eliminating effect of rigid out-of-plane motion on 2D-DIC
NASA Astrophysics Data System (ADS)
Zhiqiang, Wang; Fengzhou, Fang; Bing, Liu; Zhiyong, Wang
2015-10-01
The out-of-plane motion is one of the most important factors that affect the precision of two-dimensional digital image correlation (2D-DIC). In this paper, a novel solution is presented to improve conventional 2D-DIC by eliminating the effect of out-of-plane motion, including translation and rotation. Firstly, an experimental technique using two projected laser strips is proposed to measure the out-of-plane motion of a planar specimen. A theoretical model is then established to predict the pseudostrains caused by out-of-plane motion based on the pin-hole imaging model. Using the measured out-of-plane displacement, the captured deformed images used in 2D-DIC are amended to eliminate the effect of out-of-plane motion by the theoretical model. Finally, two experiments were conducted to validate the effectiveness of the proposed method. Results indicate that application of the proposed method can effectively eliminate the errors caused by out-of-plane motion.
2D numerical modelling of fluid percolation in the subduction zone
NASA Astrophysics Data System (ADS)
Dymkova, D.; Gerya, T.; Podladchikov, Y.
2012-04-01
Subducting slab dehydration and resulting aqueous fluid percolation triggers partial melting in the mantle wedge and is accompanied with the further melt percolation through the porous space to the region above the slab. This problem is a complex coupled chemical, thermal and mechanical process responsible for the magmatic arcs formation and change of the mantle wedge properties. We have created a two-dimensional model of a two-phase flow in a porous media solving a coupled Darcy-Stokes system of equations for two incompressible media for the case of nonlinear visco-plastic rheology of solid matrix. Our system of equation is expanded for the high-porosity limits and stabilized for the case of high porosity contrasts. We use a finite-difference method with fully staggered grid in a combination with marker-in-cell technique for advection of fluid and solid phase. We performed a comparison with a benchmark of a thermal convection in a porous media in a bottom-heated box to verify the interdependency of Rayleigh and Nusselt numbers with earlier obtained ones (Cherkaoui & Wilcock, 1999). We have demonstrated the stability and robustness of the algorithm in case of strongly non-linear visco-plastic rheology of solid including cases with localization of both deformation and porous flow along spontaneously forming shear bands. We have checked our model for the forming of localized porous channels under a simple shear stress (Katz et al, 2006). We have developed a setup of a self-initiating due to gravitational instability subduction. With our coupled fluid-solid flow we have achieved a self-consistent water downward suction by a slab bending predicted by the other models with a simplified fluid kinematical motion implementation (Faccenda et al, 2009). With this setup we have obtained a self-consistent upper crust weakening by a porous fluid pressure which was theoretically assumed in the previously existing subduction models (Gerya & Meilick, 2011; Faccenda et al, 2009
Chen, Mingqing; Bai, Junjie; Siochi, R Alfredo C
2013-02-01
To present a new method of estimating 3D positions of the ipsi-lateral hemi-diaphragm apex (IHDA) from 2D projection images of mega-voltage cone beam CT (MVCBCT). The detection framework reconstructs a 3D volume from all the 2D projection images. An initial estimated 3D IHDA position is determined in this volume based on an imaging processing pipeline, including Otsu thresholding, connected component labeling and template matching. This initial position is then projected onto each 2D projection image to create a region of interest (ROI). To accurately detect the IHDA position in 2D projection space, two methods, dynamic Hough transform (DHT) and a tracking approach based on a joint probability density function (PDF) are developed. Both methods utilize a double-parabola model to fit the 2D diaphragm boundary. The 3D IHDA motion in the superior-inferior (SI) direction is estimated from the initial static 3D position and the detected 2D positions in projection space. The two Hough-based detection methods are tested on 35 MVCBCT scans from 15 patients. The detection is compared to manually identified IHDA positions in 2D projection space by three clinicians. An average and standard deviation of 4.252 ± 3.354 and 2.485 ± 1.750 mm was achieved for DHT and tracking-based approaches respectively, compared with the inter-expert variance among three experts of 1.822 ± 1.106 mm. Based on the results of the scans, the PDF tracking-based approach appears more robust than the DHT. The combination of the automatic ROI localization and the tracking-based approach is a quicker and more accurate method of extracting 3D IHDA motion from 2D projection images. PMID:23321998
Dynamics and cortical distribution of neural responses to 2D and 3D motion in human
McKee, Suzanne P.; Norcia, Anthony M.
2013-01-01
The perception of motion-in-depth is important for avoiding collisions and for the control of vergence eye-movements and other motor actions. Previous psychophysical studies have suggested that sensitivity to motion-in-depth has a lower temporal processing limit than the perception of lateral motion. The present study used functional MRI-informed EEG source-imaging to study the spatiotemporal properties of the responses to lateral motion and motion-in-depth in human visual cortex. Lateral motion and motion-in-depth displays comprised stimuli whose only difference was interocular phase: monocular oscillatory motion was either in-phase in the two eyes (lateral motion) or in antiphase (motion-in-depth). Spectral analysis was used to break the steady-state visually evoked potentials responses down into even and odd harmonic components within five functionally defined regions of interest: V1, V4, lateral occipital complex, V3A, and hMT+. We also characterized the responses within two anatomically defined regions: the inferior and superior parietal cortex. Even harmonic components dominated the evoked responses and were a factor of approximately two larger for lateral motion than motion-in-depth. These responses were slower for motion-in-depth and were largely independent of absolute disparity. In each of our regions of interest, responses at odd-harmonics were relatively small, but were larger for motion-in-depth than lateral motion, especially in parietal cortex, and depended on absolute disparity. Taken together, our results suggest a plausible neural basis for reduced psychophysical sensitivity to rapid motion-in-depth. PMID:24198326
A case study of fluid flow in fractured rock mass based on 2-D DFN modeling
NASA Astrophysics Data System (ADS)
Han, Jisu; Noh, Young-Hwan; Um, Jeong-Gi; Choi, Yosoon
2014-05-01
A two dimensional steady-state fluid flow through fractured rock mass of an abandoned copper mine in Korea is addressed based on discrete fracture network modeling. An injection well and three observation wells were installed at the field site to monitor the variations of total heads induced by injection of fresh water. A series of packer tests were performed to estimate the rock mass permeability. First, the two dimensional stochastic fracture network model was built and validated for a granitic rock mass using the geometrical and statistical data obtained from surface exposures and borehole logs. This validated fracture network model was combined with the fracture data observed on boreholes to generate a stochastic-deterministic fracture network system. Estimated apertures for each of the fracture sets using permeability data obtained from borehole packer tests were discussed next. Finally, a systematic procedure for fluid flow modeling in fractured rock mass in two dimensional domain was presented to estimate the conductance, flow quantity and nodal head in 2-D conceptual linear pipe channel network. The results obtained in this study clearly show that fracture geometry parameters (orientation, density and size) play an important role in the hydraulic behavior of fractured rock masses.
Dynamic Linkages Between the Transition Zone & Surface Plate Motions in 2D Models of Subduction
NASA Astrophysics Data System (ADS)
Arredondo, K.; Billen, M. I.
2012-12-01
Descending subducted slabs affect both plate tectonics at the surface and overall mantle flow (e.g. Conrad and Lithgow-Bertelloni, 2002). For time-dependent numerical models, the potential evolution of these slabs, ranging from immediate penetration into the lower mantle to prior buckling and stagnation, are affected by parameters such as the plate age, the viscosity jump into the lower mantle, the presence of phase transitions, trench motion and the chosen governing equation approximation (e.g. Billen and Hirth, 2007). Similarly, the overall deviatoric stress within the slab, especially where modified by the phase transitions, may explain the uneven distribution of deep earthquakes with depth (e.g. Bina, 1997). Better understanding of these processes may arise from a more realistic 2-D model that is fully-dynamic, with an overriding plate, freely-moving trench, compositionally-layered slab and seven major phase transitions, in addition to using the compressible (TALA) form of the governing equations. Though the thermodynamic parameters of certain phase transitions may be uncertain, this study aims to test the latest data and encourage further mineralogical research. We will present fully-dynamic models, which explore the importance of the phase transitions, especially those that have been previously excluded such as the wadsleyite to ringwoodite and the pyroxene and garnet phase transitions. These phase transitions, coupled with the modeled compositionally distinct crust, harzburgite, and pyrolite lithosphere layers, may produce new large-scale dynamic behavior not seen in past numerical models, as well as stress variations within the slab related to deep slab seismicity. Feedback from the compositionally complex slab to the dynamic trench may provide further insight on the mechanics of slab stagnation and behavior in the upper and lower mantle. Billen, M. I., and G. Hirth, Rheologic controls on slab dynamics, Geochemistry, Geophysics and Geosystems, 8 (Q08012
Numerical study of fluid motion in bioreactor with two mixers
NASA Astrophysics Data System (ADS)
Zheleva, I.; Lecheva, A.
2015-10-01
Numerical study of hydrodynamic laminar behavior of a viscous fluid in bioreactor with multiple mixers is provided in the present paper. The reactor is equipped with two disk impellers. The fluid motion is studied in stream function-vorticity formulation. The calculations are made by a computer program, written in MATLAB. The fluid structure is described and numerical results are graphically presented and commented.
Numerical study of fluid motion in bioreactor with two mixers
Zheleva, I.; Lecheva, A.
2015-10-28
Numerical study of hydrodynamic laminar behavior of a viscous fluid in bioreactor with multiple mixers is provided in the present paper. The reactor is equipped with two disk impellers. The fluid motion is studied in stream function-vorticity formulation. The calculations are made by a computer program, written in MATLAB. The fluid structure is described and numerical results are graphically presented and commented.
2-D Three Fluid Simulation of Upstreaming Ions Above Auroral Precipitation
NASA Astrophysics Data System (ADS)
Danielides, M. A.; Lummerzheim, D.; Otto, A.; Stevens, R. J.
2006-12-01
The ionosphere is a rich reservoir of charged particles from which a variable fraction is transported to the magnetosphere. An important transport phenomena is the formation of upward ion flow above auroral structure. A primary region of the outflow is not known, but contributions come from polar cap, dayside cusp/cleft region, auroral oval, or even from mid-latitudes. In the past global magnetospheric models and fluid codes were used to simulate large scale ion outflow above, e.g., the polar-cap aurora. However, satellites orbiting at low- altitudes have repeatingly detected localized ion outflow above the auroral oval. Ionosphere-magnetosphere coupling simulations gave first insides into the small-scale dynamics of aurora. The aim of this study is the investigation of coupled plasma and neutral dynamics in smaller scale aurora to explain the generation, structure, and dynamics of vertical ion upstream. We consider auroral electron precipitation at ionospheric heights in a 2-D three fluid ionospheric-magnetospheric coupling code (Otto and Zhu, 2003). Specially we examine the effects of the electron precipitation, heat conduction and heating in field- aligned current through coulomb collisions or turbulence causing: i) electron heating, ii) electron pressure gradients, and iii) upstreaming of ions through a resulting ambipolar electric field. Our first case studies are performed for different boundary conditions and for different auroral electron precipitation parameters (variation in characteristic auroral energy, auroral energy flux and horizontal scale). The results shall clarify how auroral precipitation can drive ions upwards. Finally we discuss the effect of ion drag and the interaction of the upstreaming ions with a stable neutral constituent. Otto, O. and H. Zhu, Fluid plasma simulation of coupled systems: Ionosphere and magnetosphere, Space Plasma Simulation. Edited by J. Buechner, C. Dum, and M. Scholer., Lecture Notes in Physics, vol. 615, p.193
Kim, Young-Keun; Kim, Kyung-Soo
2014-10-15
Maritime transportation demands an accurate measurement system to track the motion of oscillating container boxes in real time. However, it is a challenge to design a sensor system that can provide both reliable and non-contact methods of 6-DOF motion measurements of a remote object for outdoor applications. In the paper, a sensor system based on two 2D laser scanners is proposed for detecting the relative 6-DOF motion of a crane load in real time. Even without implementing a camera, the proposed system can detect the motion of a remote object using four laser beam points. Because it is a laser-based sensor, the system is expected to be highly robust to sea weather conditions.
NASA Astrophysics Data System (ADS)
Kim, Young-Keun; Kim, Kyung-Soo
2014-10-01
Maritime transportation demands an accurate measurement system to track the motion of oscillating container boxes in real time. However, it is a challenge to design a sensor system that can provide both reliable and non-contact methods of 6-DOF motion measurements of a remote object for outdoor applications. In the paper, a sensor system based on two 2D laser scanners is proposed for detecting the relative 6-DOF motion of a crane load in real time. Even without implementing a camera, the proposed system can detect the motion of a remote object using four laser beam points. Because it is a laser-based sensor, the system is expected to be highly robust to sea weather conditions.
Petasecca, M. Newall, M. K.; Aldosari, A. H.; Fuduli, I.; Espinoza, A. A.; Porumb, C. S.; Guatelli, S.; Metcalfe, P.; Lerch, M. L. F.; Rosenfeld, A. B.; Booth, J. T.; Colvill, E.; Duncan, M.; Cammarano, D.; Carolan, M.; Oborn, B.; Perevertaylo, V.; Keall, P. J.
2015-06-15
Purpose: Spatial and temporal resolutions are two of the most important features for quality assurance instrumentation of motion adaptive radiotherapy modalities. The goal of this work is to characterize the performance of the 2D high spatial resolution monolithic silicon diode array named “MagicPlate-512” for quality assurance of stereotactic body radiation therapy (SBRT) and stereotactic radiosurgery (SRS) combined with a dynamic multileaf collimator (MLC) tracking technique for motion compensation. Methods: MagicPlate-512 is used in combination with the movable platform HexaMotion and a research version of radiofrequency tracking system Calypso driving MLC tracking software. The authors reconstruct 2D dose distributions of small field square beams in three modalities: in static conditions, mimicking the temporal movement pattern of a lung tumor and tracking the moving target while the MLC compensates almost instantaneously for the tumor displacement. Use of Calypso in combination with MagicPlate-512 requires a proper radiofrequency interference shielding. Impact of the shielding on dosimetry has been simulated by GEANT4 and verified experimentally. Temporal and spatial resolutions of the dosimetry system allow also for accurate verification of segments of complex stereotactic radiotherapy plans with identification of the instant and location where a certain dose is delivered. This feature allows for retrospective temporal reconstruction of the delivery process and easy identification of error in the tracking or the multileaf collimator driving systems. A sliding MLC wedge combined with the lung motion pattern has been measured. The ability of the MagicPlate-512 (MP512) in 2D dose mapping in all three modes of operation was benchmarked by EBT3 film. Results: Full width at half maximum and penumbra of the moving and stationary dose profiles measured by EBT3 film and MagicPlate-512 confirm that motion has a significant impact on the dose distribution. Motion
Theory of Brownian motion in a Jeffreys fluid
Raikher, Yu. L.; Rusakov, V. V.
2010-11-15
We have constructed a kinetic theory of Brownian motion in a rheologically complex medium-a Jeffreys fluid that is characterized by a combination of two viscosity mechanisms: ordinary and delayed. This model is shown to be much better suited for the interpretation of experiments on the microrheology of viscoelastic media than the standard Maxwell model. In particular, no oscillations of the mean-square particle displacement arise in a Jeffreys fluid, which is a nonremovable artifact of the theory of Brownian motion in a Maxwell fluid. The developed approach can to be used also consider the diffusion of particles in other complex fluids whose rheology is described by phenomenological schemes.
2D fluid simulations of acoustic waves in pulsed ICP discharges: Comparison with experiments
NASA Astrophysics Data System (ADS)
Despiau-Pujo, Emilie; Cunge, Gilles; Sadeghi, Nader; Braithwaite, N. St. J.
2012-10-01
Neutral depletion, which is mostly caused by gas heating under typical material processing conditions, is an important phenomenon in high-density plasmas. In low pressure pulsed discharges, experiments show that additional depletion due to electron pressure (Pe) may have a non-negligible influence on radical transport [1]. To evaluate this effect, comparisons between 2D fluid simulations and measurements of gas convection in Ar/Cl2 pulsed ICP plasmas are reported. In the afterglow, Pe drops rapidly by electron cooling which generates a neutral pressure gradient between the plasma bulk and the reactor walls. This in turn forces the cold surrounding gas to move rapidly towards the center, thus launching an acoustic wave in the reactor. Time-resolved measurements of atoms drift velocity and gas temperature by LIF and LAS in the early afterglow are consistent with gas drifting at acoustic wave velocity followed by rapid gas cooling. Similar results are predicted by the model. The ion flux at the reactor walls is also shown to oscillate in phase with the acoustic wave due to ion-neutral friction forces. Finally, during plasma ignition, experiments show opposite phenomena when Pe rises.[4pt] [1] Cunge et al, APL 96, 131501 (2010)
The Accuracy of Webcams in 2D Motion Analysis: Sources of Error and Their Control
ERIC Educational Resources Information Center
Page, A.; Moreno, R.; Candelas, P.; Belmar, F.
2008-01-01
In this paper, we show the potential of webcams as precision measuring instruments in a physics laboratory. Various sources of error appearing in 2D coordinate measurements using low-cost commercial webcams are discussed, quantifying their impact on accuracy and precision, and simple procedures to control these sources of error are presented.…
On the coupling between fluid flow and mesh motion in the modelling of fluid structure interaction
NASA Astrophysics Data System (ADS)
Dettmer, Wulf G.; Perić, Djordje
2008-12-01
Partitioned Newton type solution strategies for the strongly coupled system of equations arising in the computational modelling of fluid solid interaction require the evaluation of various coupling terms. An essential part of all ALE type solution strategies is the fluid mesh motion. In this paper, we investigate the effect of the terms which couple the fluid flow with the fluid mesh motion on the convergence behaviour of the overall solution procedure. We show that the computational efficiency of the simulation of many fluid solid interaction processes, including fluid flow through flexible pipes, can be increased significantly if some of these coupling terms are calculated exactly.
Vehicular motion in 2D city traffic network with signals controlled by phase shift
NASA Astrophysics Data System (ADS)
Komada, Kazuhito; Kojima, Kengo; Nagatani, Takashi
2011-03-01
We study the dynamic behavior of vehicular traffic through the series of traffic lights controlled by phase shift in two-dimensional (2D) city traffic network. The nonlinear-map model is presented for the vehicular traffic. The city traffic network is made of one-way perpendicular streets arranged in a square lattice with traffic signals where vertical streets are oriented upwards and horizontal streets are oriented rightwards. There are two traffic lights for the movement to north or that to east at each crossing. The traffic lights are controlled by the cycle time, split, and phase shift. The vehicle moves through the series of signals on a path selected by the driver. The city traffic with a heterogeneous density distribution is also studied. The dependence of the arrival time on cycle time, split, phase shift, selected path, and density is clarified for 2D city traffic. It is shown that the vehicular traffic is efficiently controlled by the phase shift.
Multi-level model for 2D human motion analysis and description
NASA Astrophysics Data System (ADS)
Foures, Thomas; Joly, Philippe
2003-01-01
This paper deals with the proposition of a model for human motion analysis in a video. Its main caracteristic is to adapt itself automatically to the current resolution, the actual quality of the picture, or the level of precision required by a given application, due to its possible decomposition into several hierarchical levels. The model is region-based to address some analysis processing needs. The top level of the model is only defined with 5 ribbons, which can be cut into sub-ribbons regarding to a given (or an expected) level of details. Matching process between model and current picture consists in the comparison of extracted subject shape with a graphical rendering of the model built on the base of some computed parameters. The comparison is processed by using a chamfer matching algorithm. In our developments, we intend to realize a platform of interaction between a dancer and tools synthetizing abstract motion pictures and music in the conditions of a real-time dialogue between a human and a computer. In consequence, we use this model in a perspective of motion description instead of motion recognition: no a priori gestures are supposed to be recognized as far as no a priori application is specially targeted. The resulting description will be made following a Description Scheme compliant with the movement notation called "Labanotation".
1 kHz 2D Visual Motion Sensor Using 20 × 20 Silicon Retina Optical Sensor and DSP Microcontroller.
Liu, Shih-Chii; Yang, MinHao; Steiner, Andreas; Moeckel, Rico; Delbruck, Tobi
2015-04-01
Optical flow sensors have been a long running theme in neuromorphic vision sensors which include circuits that implement the local background intensity adaptation mechanism seen in biological retinas. This paper reports a bio-inspired optical motion sensor aimed towards miniature robotic and aerial platforms. It combines a 20 × 20 continuous-time CMOS silicon retina vision sensor with a DSP microcontroller. The retina sensor has pixels that have local gain control and adapt to background lighting. The system allows the user to validate various motion algorithms without building dedicated custom solutions. Measurements are presented to show that the system can compute global 2D translational motion from complex natural scenes using one particular algorithm: the image interpolation algorithm (I2A). With this algorithm, the system can compute global translational motion vectors at a sample rate of 1 kHz, for speeds up to ±1000 pixels/s, using less than 5 k instruction cycles (12 instructions per pixel) per frame. At 1 kHz sample rate the DSP is 12% occupied with motion computation. The sensor is implemented as a 6 g PCB consuming 170 mW of power. PMID:25879969
Development of models for the two-dimensional, two-fluid code for sodium boiling NATOF-2D. [LMFBR
Zielinski, R.G.; Kazimi, M.S.
1981-09-01
Several features were incorporated into NATOF-2D, a two-dimensional, two fluid code developed at MIT for the purpose of analysis of sodium boiling transients under LMFBR conditions. They include improved interfacial mass, momentum and energy exchange rate models, and a cell-to-cell radial heat conduction mechanism which was calibrated by simulation of Westinghouse Blanket Heat Transfer Test Program Runs 544 and 545. Finally, a direct method of pressure field solution was implemented into a direct method of pressure field solution was implemented into NATOF-2D, replacing the iterative technique previously available, and resulted in substantially reduced computational costs.
Binocular Perception of 2D Lateral Motion and Guidance of Coordinated Motor Behavior.
Fath, Aaron J; Snapp-Childs, Winona; Kountouriotis, Georgios K; Bingham, Geoffrey P
2016-04-01
Zannoli, Cass, Alais, and Mamassian (2012) found greater audiovisual lag between a tone and disparity-defined stimuli moving laterally (90-170 ms) than for disparity-defined stimuli moving in depth or luminance-defined stimuli moving laterally or in depth (50-60 ms). We tested if this increased lag presents an impediment to visually guided coordination with laterally moving objects. Participants used a joystick to move a virtual object in several constant relative phases with a laterally oscillating stimulus. Both the participant-controlled object and the target object were presented using a disparity-defined display that yielded information through changes in disparity over time (CDOT) or using a luminance-defined display that additionally provided information through monocular motion and interocular velocity differences (IOVD). Performance was comparable for both disparity-defined and luminance-defined displays in all relative phases. This suggests that, despite lag, perception of lateral motion through CDOT is generally sufficient to guide coordinated motor behavior. PMID:26614099
Energy Exchange during Plunge/Surge Motions of a 2D Wing
NASA Astrophysics Data System (ADS)
Kerstens, Wesley; Choi, Jeesoon; Colonius, Tim; Williams, David
2011-11-01
The rate of energy transfer between an NACA-0006 wing and an unsteady flow is examined at pre-stall and post-stall conditions using numerical simulations and wind tunnel experiments. The plunge and surge motions simulate the fluctuating vertical (wz) and longitudinal (wx) velocity components of a wind gust. In a steady flow the wing loses energy to the flow through the drag power term, but in an unsteady flow the wing may gain energy from the fluctuating lift power and fluctuating drag power terms. The net energy transfer averaged over the period of oscillation depends on the phase angle between the plunge and surge motions. The largest increase of energy occurs when wx and wz are in-phase. When the fluctuations are large enough, then it is possible for the net energy gain to be positive. The numerical simulations conducted at Reynolds numbers near the critical value for vortex shedding show qualitative agreement with the experiments. The simulations highlight the role of vortex shedding in determining the optimal frequency and phase for energy extraction from the gust. Support of the AFOSR through grant FA9550-09-1-0189 managed by Dr. Douglas Smith is gratefully acknnowledged.
Nonlinear state-space modeling of human motion using 2-D marker observations.
Vartiainen, Paavo; Bragge, Timo; Arokoski, Jari P; Karjalainen, Pasi A
2014-07-01
A novel method for the estimation of human kinematics, based on state-space modeling, is proposed. The state consists of the positions, orientations, velocities, and accelerations of an articulated model. Estimation is performed using the unscented Kalman filter (UKF) algorithm with a fixed-interval smoother. Impulsive acceleration at floor contact of the foot is estimated by implementing a contact constraint in the UKF evolution model. The constraint inserts an acceleration impulse into the model state. The estimation method was applied to marker-based motion analysis in a motion laboratory. Validation measurements were performed with a rigid test device and with human gait. A triaxial accelerometer was used to evaluate acceleration estimates. Comparison between the proposed method and the extended Kalman smoother showed a clear difference in the quality of estimates during impulsive accelerations. The proposed approach enables estimation of human kinematics during both continuous and transient accelerations. The approach provides a novel way of estimating acceleration at foot initial contact, and thus enables more accurate evaluation of loading from the beginning of the floor contact. PMID:24760898
Robust 2D/3D registration for fast-flexion motion of the knee joint using hybrid optimization.
Ohnishi, Takashi; Suzuki, Masahiko; Kobayashi, Tatsuya; Naomoto, Shinji; Sukegawa, Tomoyuki; Nawata, Atsushi; Haneishi, Hideaki
2013-01-01
Previously, we proposed a 2D/3D registration method that uses Powell's algorithm to obtain 3D motion of a knee joint by 3D computed-tomography and bi-plane fluoroscopic images. The 2D/3D registration is performed consecutively and automatically for each frame of the fluoroscopic images. This method starts from the optimum parameters of the previous frame for each frame except for the first one, and it searches for the next set of optimum parameters using Powell's algorithm. However, if the flexion motion of the knee joint is fast, it is likely that Powell's algorithm will provide a mismatch because the initial parameters are far from the correct ones. In this study, we applied a hybrid optimization algorithm (HPS) combining Powell's algorithm with the Nelder-Mead simplex (NM-simplex) algorithm to overcome this problem. The performance of the HPS was compared with the separate performances of Powell's algorithm and the NM-simplex algorithm, the Quasi-Newton algorithm and hybrid optimization algorithm with the Quasi-Newton and NM-simplex algorithms with five patient data sets in terms of the root-mean-square error (RMSE), target registration error (TRE), success rate, and processing time. The RMSE, TRE, and the success rate of the HPS were better than those of the other optimization algorithms, and the processing time was similar to that of Powell's algorithm alone. PMID:23138929
Towards real-time 2D/3D registration for organ motion monitoring in image-guided radiation therapy
NASA Astrophysics Data System (ADS)
Gendrin, C.; Spoerk, J.; Bloch, C.; Pawiro, S. A.; Weber, C.; Figl, M.; Markelj, P.; Pernus, F.; Georg, D.; Bergmann, H.; Birkfellner, W.
2010-02-01
Nowadays, radiation therapy systems incorporate kV imaging units which allow for the real-time acquisition of intra-fractional X-ray images of the patient with high details and contrast. An application of this technology is tumor motion monitoring during irradiation. For tumor tracking, implanted markers or position sensors are used which requires an intervention. 2D/3D intensity based registration is an alternative, non-invasive method but the procedure must be accelerate to the update rate of the device, which lies in the range of 5 Hz. In this paper we investigate fast CT to a single kV X-ray 2D/3D image registration using a new porcine reference phantom with seven implanted fiducial markers. Several parameters influencing the speed and accuracy of the registrations are investigated. First, four intensity based merit functions, namely Cross-Correlation, Rank Correlation, Mutual Information and Correlation Ratio, are compared. Secondly, wobbled splatting and ray casting rendering techniques are implemented on the GPU and the influence of each algorithm on the performance of 2D/3D registration is evaluated. Rendering times for a single DRR of 20 ms were achieved. Different thresholds of the CT volume were also examined for rendering to find the setting that achieves the best possible correspondence with the X-ray images. Fast registrations below 4 s became possible with an inplane accuracy down to 0.8 mm.
2D-3D rigid registration to compensate for prostate motion during 3D TRUS-guided biopsy
NASA Astrophysics Data System (ADS)
De Silva, Tharindu; Fenster, Aaron; Bax, Jeffrey; Gardi, Lori; Romagnoli, Cesare; Samarabandu, Jagath; Ward, Aaron D.
2012-02-01
Prostate biopsy is the clinical standard for prostate cancer diagnosis. To improve the accuracy of targeting suspicious locations, systems have been developed that can plan and record biopsy locations in a 3D TRUS image acquired at the beginning of the procedure. Some systems are designed for maximum compatibility with existing ultrasound equipment and are thus designed around the use of a conventional 2D TRUS probe, using controlled axial rotation of this probe to acquire a 3D TRUS reference image at the start of the biopsy procedure. Prostate motion during the biopsy procedure causes misalignments between the prostate in the live 2D TRUS images and the pre-acquired 3D TRUS image. We present an image-based rigid registration technique that aligns live 2D TRUS images, acquired immediately prior to biopsy needle insertion, with the pre-acquired 3D TRUS image to compensate for this motion. Our method was validated using 33 manually identified intrinsic fiducials in eight subjects and the target registration error was found to be 1.89 mm. We analysed the suitability of two image similarity metrics (normalized cross correlation and mutual information) for this task by plotting these metrics as a function of varying parameters in the six degree-of-freedom transformation space, with the ground truth plane obtained from registration as the starting point for the parameter exploration. We observed a generally convex behaviour of the similarity metrics. This encourages their use for this registration problem, and could assist in the design of a tool for the detection of misalignment, which could trigger the execution of a non-real-time registration, when needed during the procedure.
NASA Astrophysics Data System (ADS)
Sharapov, V. N.; Cherepanov, A. N.; Popov, V. N.; Bykova, V. G.
2012-11-01
A model describing two-dimensional (2D) dynamics of heat transfer in the fluid systems with a localized sink of a magmatic fluid into local fractured zones above the roof of crystallizing crustal intrusions is suggested. Numerical modeling of the migration of the phase boundaries in 2D intrusive chambers under retrograde boiling of magma with relatively high initial water content in the melt shows that, depending on the character of heat dissipation from a magmatic fluid into the host rock, two types of fluid magmatic systems can arise. (1) At high heat losses, the zoning of fluidogenic ore formation is determined by the changes in temperature of the rocks within the contact aureole of the intrusive bodies. These temperature variations are controlled by the migration of the phase boundaries in the cooling melt towards the center of the magmatic bodies from their contacts. (2) In the case of a localized sink of the magmatic fluid in different parts of the top of the intrusive chambers, a specific characteristic scenario of cooling of the magmatic bodies is probably implemented. In 2D systems with a heat transfer coefficient α k < 5 × 104 W/m2 K, an area with quasi-stationary phase boundaries develops close to the region of fluid drainage through the fractured zone in the intrusion. Therefore, as the phase boundaries contract to the sink zone of a fluid, specific thermal tubes arise, whose characteristics depend on the width of the fluid-conductive zone and the heat losses into the side rocks. (3) The time required for the intrusion to solidify varies depending on the particular position of the fluid conductor above the top of the magmatic body.
SAGE 2D and 3D Simulations of the Explosive Venting of Supercritical Fluids Through Porous Media
NASA Astrophysics Data System (ADS)
Weaver, R.; Gisler, G.; Svensen, H.; Mazzini, A.
2008-12-01
Magmatic intrusive events in large igneous provinces heat sedimentary country rock leading to the eventual release of volatiles. This has been proposed as a contributor to climate change and other environmental impacts. By means of numerical simulations, we examine ways in which these volatiles can be released explosively from depth. Gases and fluids cooked out of country rock by metamorphic heating may be confined for a time by impermeable clays or other barriers, developing high pressures and supercritical fluids. If confinement is suddenly breached (by an earthquake for example) in such a way that the fluid has access to porous sediments, a violent eruption of a non-magmatic mixture of fluid and sediment may result. Surface manifestations of these events could be hydrothermal vent complexes, kimberlite pipes, pockmarks, or mud volcanoes. These are widespread on Earth, especially in large igneous provinces, as in the Karoo Basin of South Africa, the North Sea off the Norwegian margin, and the Siberian Traps. We have performed 2D and 3D simulations with the Sage hydrocode (from Los Alamos and Science Applications International) of supercritical venting in a variety of geometries and configurations. The simulations show several different patterns of propagation and fracturing in porous or otherwise weakened overburden, dependent on depth, source conditions (fluid availability, temperature, and pressure), and manner of confinement breach. Results will be given for a variety of 2D and 3D simulations of these events exploring the release of volatiles into the atmosphere.
Physiological flow of Carreau fluid due to ciliary motion
NASA Astrophysics Data System (ADS)
Nadeem, S.; Munim, A.; Shaheen, A.; Hussain, S.
2016-03-01
In this paper, we have explained the ciliary motion of a Carreau fluid inside a symmetrical channel with ciliated walls. Carreau fluid equations for the two dimensional channel are interpreted by employing the low Reynolds number and long wave-length approximations. The analytical series solutions are found by employing the Homotopy perturbation method. Solutions are presented through graphs. The results for the distinct values of the constants are described by generating graphs of velocity, pressure gradient, pressure rise and stream lines.
Modeling Selective Local Interactions with Memory: Motion on a 2D Lattice.
Weinberg, Daniel; Levy, Doron
2014-06-15
We consider a system of particles that simultaneously move on a two-dimensional periodic lattice at discrete times steps. Particles remember their last direction of movement and may either choose to continue moving in this direction, remain stationary, or move toward one of their neighbors. The form of motion is chosen based on predetermined stationary probabilities. Simulations of this model reveal a connection between these probabilities and the emerging patterns and size of aggregates. In addition, we develop a reaction diffusion master equation from which we derive a system of ODEs describing the dynamics of the particles on the lattice. Simulations demonstrate that solutions of the ODEs may replicate the aggregation patterns produced by the stochastic particle model. We investigate conditions on the parameters that influence the locations at which particles prefer to aggregate. This work is a two-dimensional generalization of [Galante & Levy, Physica D, http://dx.doi.org/10.1016/j.physd.2012.10.010], in which the corresponding one-dimensional problem was studied. PMID:25045193
Dynamic Linkages Between the Transition Zone & Surface Plate Motions in 2D Models of Subduction
NASA Astrophysics Data System (ADS)
Arredondo, K.; Billen, M. I.
2013-12-01
While slab pull is considered the dominant force controlling plate motion and speed, its magnitude is controlled by slab behavior in the mantle, where tomographic studies show a wide range of possibilities from direct penetration to folding, or stagnation directly above the lower mantle (e.g. Fukao et al., 2009). Geodynamic studies have investigated various parameters, such as plate age and two phase transitions, to recreate observed behavior (e.g. Běhounková and Cízková, 2008). However, past geodynamic models have left out known slab characteristics that may have a large impact on slab behavior and our understanding of subduction processes. Mineral experiments and seismic observations have indicated the existence of additional phase transitions in the mantle transition zone that may produce buoyancy forces large enough to affect the descent of a subducting slab (e.g. Ricard et al., 2005). The current study systematically tests different common assumptions used in geodynamic models: kinematic versus free-slip boundary conditions, the effects of adiabatic heating, viscous dissipation and latent heat, compositional layering and a more complete suite of phase transitions. Final models have a complete energy equation, with eclogite, harzburgite and pyrolite lithosphere compositional layers, and seven composition-dependent phase transitions within the olivine, pyroxene and garnet polymorph minerals. Results show important feedback loops between different assumptions and new behavior from the most complete models. Kinematic models show slab weakening or breaking above the 660 km boundary and between compositional layers. The behavior in dynamic models with a free-moving trench and overriding plate is compared to the more commonly found kinematic models. The new behavior may have important implications for the depth distribution of deep earthquakes within the slab. Though the thermodynamic parameters of certain phase transitions may be uncertain, their presence and
Beginning Continuous Fluid Motion in the Music Classroom.
ERIC Educational Resources Information Center
Westervelt, Todd G.
2002-01-01
Focuses on how to include movement education into the music classroom. Discusses the techniques developed by three movement specialists: (1) Edwin Gordon's Continuous Fluid Motion (CFM); (2) Rudolf von Laban's effort/flow elements; and (3) Phyllis Weikart's taxonomy of movement. Includes a bibliography of resources. (CMK)
Large-amplitude inviscid fluid motion in an accelerating container
NASA Technical Reports Server (NTRS)
Perko, L. M.
1968-01-01
Study of dynamic behavior of the liquid-vapor interface of an inviscid fluid in an accelerating cylindrical container includes an analytical-numerical method for determining large amplitude motion. The method is based on the expansion of the velocity potential in a series of harmonic functions with time dependent coefficients.
Microgravity and its effects on residual motions in fluids
NASA Technical Reports Server (NTRS)
Alexander, J. Iwan D.; Lundquist, Charles A.
1990-01-01
The primary reason for conducting many materials science experiments in space is to minimize or eliminate undesirable effects that might result owing to convective motions in fluids that are driven by buoyancy effects. Of particular concern are the low frequency accelerations caused by the Earth's gravity gradient field, spacecraft attitude motions, and atmospheric drag. In order to gain a limited understanding of the effects of these accelerations, researchers calculated the Stokes' motion of a spherical particle in a fluid for various types of spacecraft attitudes. Researchers assessed the effect of slowly rotating the experimental system relative to the spacecraft in order to reduce the rate at which the particles accumulate against the container wall.
NASA Astrophysics Data System (ADS)
Pan, Li-Hua; Hou, Peng-Fei; Chen, Jia-Yun
2016-08-01
The 2D steady-state solutions regarding the expressions of stress and strain for fluid-saturated, orthotropic, poroelastic plane are derived in this paper. For this object, the general solutions of the corresponding governing equation are first obtained and expressed in harmonic functions. Based on these compact general solutions, the suitable harmonic functions with undetermined constants for line fluid source in the interior of infinite poroelastic body and a line fluid source on the surface of semi-infinite poroelastic body are presented, respectively. The fundamental solutions can be obtained by substituting these functions into the general solution, and the undetermined constants can be obtained by the continuous conditions, equilibrium conditions and boundary conditions.
Vlasov Fluid stability of a 2-D plasma with a linear magnetic field null
Kim, J.S.
1984-01-01
Vlasov Fluid stability of a 2-dimensional plasma near an O type magnetic null is investigated. Specifically, an elongated Z-pinch is considered, and applied to Field Reversed Configurations at Los Alamos National Laboratory by making a cylindrical approximation of the compact torus. The orbits near an elliptical O type null are found to be very complicated; the orbits are large and some are stochastic. The kinetic corrections to magnetohydrodynamics (MHD) are investigated by evaluating the expectation values of the growth rates of a Vlasov Fluid dispersion functional by using a set of trial functions based on ideal MHD. The dispersion functional involves fluid parts and orbit dependent parts. The latter involves phase integral of two time correlations. The phase integral is replaced by the time integral both for the regular and for the stochastic orbits. Two trial functions are used; one has a large displacement near the null and the other away from the null.
Comparing a 2D fluid model of the DC planar magnetron cathode to experiments
Garcia, M.
1996-05-01
Planar magnetron cathodes have arching magnetic field lines which concentrate plasma density near the electrode surface. This enhances the ion bombardment of the surface and the yield of sputtered atoms. Magnetron cathodes are used in the Plasma Electrode Pockels Cell (PEPC) devices of the Laser Program because they provide for significantly higher conduction than do glow discharges. An essential feature of magnetron cathodes is that the vector product of the perpendicular electric field, E[sub y], with the parallel component of the magnetic field, B[sub x], forms a closed track with a circulating current along the cathode surface. An analytical, 2D, two component, quasi-neutral, continuum model yields formulas for the plasma density, the total and component current densities, the electric field, and the positive electrical potential, between the cathode surface and a distant, uniform plasma. For a specific gas, the free parameters are electron temperature, gas number density, and total current. The model is applied to the interpretation of experimental data from the PEPC device, as well as a small vacuum facility for testing magnetron cathodes. Finally, the model has been applied to generate cross sectional views of a PEPC magnetron cathode track.
Anomalous diffusion of an ellipsoid in quasi-2D active fluids
NASA Astrophysics Data System (ADS)
Peng, Yi; Yang, Ou; Tang, Chao; Cheng, Xiang
Enhanced diffusion of a tracer particle is a unique feature in active fluids. Here, we studied the diffusion of an ellipsoid in a free-standing film of E. coli. Particle diffusion is linearly enhanced at low bacterial concentrations, whereas a non-linear enhancement is observed at high bacterial concentrations due to the giant fluctuation. More importantly, we uncover an anomalous coupling between the translational and rotational degrees of freedom that is strictly prohibited in the classical Brownian diffusion. Combining experiments with theoretical modeling, we show that such an anomaly arises from the stretching flow induced by the force dipole of swimming bacteria. Our work illustrates a novel universal feature of active matter and transforms the understanding of fundamental transport processes in microbiological systems. ACS Petroleum Research Fund #54168-DNI9, NSF Faculty Early Career Development Program, DMR-1452180.
Regional Multi-Fluid-Based Geophysical Excitation of Polar Motion
NASA Technical Reports Server (NTRS)
Nastula, Jolanta; Salstein, David A.; Gross, Richard
2011-01-01
By analyzing geophysical fluids geographic distribution, we can isolate the regional provenance for some of the important signals in polar motion. An understanding of such will enable us to determine whether certain climate signals can have an impact on polar motion. Here we have compared regional patterns of three surficial fluids: the atmosphere, ocean and land-based hydrosphere. The oceanic excitation function of polar motion was estimated with the ECCO/JPL data - assimilating model, and the atmospheric excitation function was determined from NCEP/NCAR reanalyses. The excitation function due to land hydrology was estimated from the Gravity Recovery and Climate Experiment (GRACE) data by an indirect approach that determines water thickness. Our attention focuses on the regional distribution of atmospheric and oceanic excitation of the annual and Chandler wobbles during 1993-2010, and on hydrologic excitation of these wobbles during 2002.9-2011.5. It is found that the regions of maximum fractional covariance (those exceeding a value of 3 .10 -3) for the annual band are over south Asia, southeast Asia and south central Indian ocean, for hydrology, atmosphere and ocean respectively; and for the Chandler period, areas over North America, Asia, and South America; and scattered across the southern oceans for the atmosphere and oceans respectively
The motion of ellipsoids in a second order fluid
NASA Astrophysics Data System (ADS)
Kim, S.
1985-09-01
The rigid body motion of an ellipsoid in a second order fluid (SOF) under the action of specified (time independent) external forces and torques have been obtained to first order in the Weissenberg number by inverting the resistance relations for the force an torque under specified rigid body motions. The reciprocal theorem of Lorentz was used to bypass the calculation of the O(W) velocity field. The results agree with known analytic solutions for SOF with the secondary to primary normal stress ratio of -1/2. The solution procedure was also tested by computing the torque on a translating prolate spheroid with aspect ratios ranging from slender bodies to near-spheres. One result is that for a SOF with zero secondary normal stress (Weissenberg fluid), previous asymptotic results for near-spheres were found to be accurate even at fairly large aspect ratios. New results of nondegenerate ellipsoids suggest that the orientation (as monitored by Euler angles) and trajectory of sedimenting, nonaxisymmetric particles such as ellipsoids provide useful information on the rheology of the suspending fluid.
Bi-planar 2D-to-3D registration in Fourier domain for stereoscopic x-ray motion tracking
NASA Astrophysics Data System (ADS)
Zosso, Dominique; Le Callennec, Benoît; Bach Cuadra, Meritxell; Aminian, Kamiar; Jolles, Brigitte M.; Thiran, Jean-Philippe
2008-03-01
In this paper we present a new method to track bone movements in stereoscopic X-ray image series of the knee joint. The method is based on two different X-ray image sets: a rotational series of acquisitions of the still subject knee that allows the tomographic reconstruction of the three-dimensional volume (model), and a stereoscopic image series of orthogonal projections as the subject performs movements. Tracking the movements of bones throughout the stereoscopic image series means to determine, for each frame, the best pose of every moving element (bone) previously identified in the 3D reconstructed model. The quality of a pose is reflected in the similarity between its theoretical projections and the actual radiographs. We use direct Fourier reconstruction to approximate the three-dimensional volume of the knee joint. Then, to avoid the expensive computation of digitally rendered radiographs (DRR) for pose recovery, we develop a corollary to the 3-dimensional central-slice theorem and reformulate the tracking problem in the Fourier domain. Under the hypothesis of parallel X-ray beams, the heavy 2D-to-3D registration of projections in the signal domain is replaced by efficient slice-to-volume registration in the Fourier domain. Focusing on rotational movements, the translation-relevant phase information can be discarded and we only consider scalar Fourier amplitudes. The core of our motion tracking algorithm can be implemented as a classical frame-wise slice-to-volume registration task. Results on both synthetic and real images confirm the validity of our approach.
On steady motion of viscoelastic fluid of Oldroyd type
Baranovskii, E. S.
2014-06-01
We consider a mathematical model describing the steady motion of a viscoelastic medium of Oldroyd type under the Navier slip condition at the boundary. In the rheological relation, we use the objective regularized Jaumann derivative. We prove the solubility of the corresponding boundary-value problem in the weak setting. We obtain an estimate for the norm of a solution in terms of the data of the problem. We show that the solution set is sequentially weakly closed. Furthermore, we give an analytic solution of the boundary-value problem describing the flow of a viscoelastic fluid in a flat channel under a slip condition at the walls. Bibliography: 13 titles. (paper)
Motion of a hot particle in viscous fluids
NASA Astrophysics Data System (ADS)
Oppenheimer, Naomi; Navardi, Shahin; Stone, Howard A.
2016-05-01
We study the motion of a hot particle in a viscous liquid at low Reynolds numbers, which is inspired by recent experiments with Brownian particles heated by a laser. The difference in temperature between a particle and the ambient fluid causes a spatial variation of the viscosity in the vicinity of the solid body. We derive a general analytical expression determining the force and the torque on a particle for low Péclet numbers by exploiting the Lorentz reciprocal theorem. For small temperature and viscosity variations, a perturbation analysis is implemented to evaluate the leading-order correction to the hydrodynamic force and torque on the particle. The results are applied to describe dynamics of a uniformly hot spherical particle and to spherical particles with a nonuniform surface temperature described by dipole and quadrupole moments. Among other results, we find for dipolar thermal fields that there is coupling of the translational and rotational motions when there are local viscosity variations; such coupling is absent in an isothermal fluid.
Flow in left atrium using MR fluid motion estimation
NASA Astrophysics Data System (ADS)
Wong, Kelvin K. L.; Kelso, Richard M.; Worthley, Steve M.; Sanders, Prash; Mazumdar, Jagannath; Abbott, Derek
2007-12-01
A recent development based on optical flow applied onto Fast Imaging in Steady State Free Precession (TrueFISP) magnetic resonance imaging is able to deliver good estimation of the flow profile in the human heart chamber. The examination of cardiac flow based on tracking of MR signals emitted by moving blood is able to give medical doctors insight into the flow patterns within the human heart using standard MRI procedure without specifically subjecting the patient to longer scan times using more dedicated scan protocols such as phase contrast MRI. Although MR fluid motion estimation has its limitations in terms of accurate flow mapping, the use of a comparatively quick scan procedure and computational post-processing gives satisfactory flow quantification and can assist in management of cardiac patients. In this study, we present flow in the left atria of five human subjects using MR fluid motion tracking. The measured flow shows that vortices exist within the atrium of heart. Although the scan is two-dimensional, we have produced multiple slices of flow maps in a spatial direction to show that the vortex exist in a three-dimensional space.
Ladstein, Jarle; Evensmoen, Hallvard R.; Håberg, Asta K.; Kristoffersen, Anders; Goa, Pål E.
2016-01-01
Purpose: To compare 2D and 3D echo-planar imaging (EPI) in a higher cognitive level fMRI paradigm. In particular, to study the link between the presence of task-correlated physiological fluctuations and motion and the fMRI contrast estimates from either 2D EPI or 3D EPI datasets, with and without adding nuisance regressors to the model. A signal model in the presence of partly task-correlated fluctuations is derived, and predictions for contrast estimates with and without nuisance regressors are made. Materials and Methods: Thirty-one healthy volunteers were scanned using 2D EPI and 3D EPI during a virtual environmental learning paradigm. In a subgroup of 7 subjects, heart rate and respiration were logged, and the correlation with the paradigm was evaluated. FMRI analysis was performed using models with and without nuisance regressors. Differences in the mean contrast estimates were investigated by analysis-of-variance using Subject, Sequence, Day, and Run as factors. The distributions of group level contrast estimates were compared. Results: Partially task-correlated fluctuations in respiration, heart rate and motion were observed. Statistically significant differences were found in the mean contrast estimates between the 2D EPI and 3D EPI when using a model without nuisance regressors. The inclusion of nuisance regressors for cardiorespiratory effects and motion reduced the difference to a statistically non-significant level. Furthermore, the contrast estimate values shifted more when including nuisance regressors for 3D EPI compared to 2D EPI. Conclusion: The results are consistent with 3D EPI having a higher sensitivity to fluctuations compared to 2D EPI. In the presence partially task-correlated physiological fluctuations or motion, proper correction is necessary to get expectation correct contrast estimates when using 3D EPI. As such task-correlated physiological fluctuations or motion is difficult to avoid in paradigms exploring higher cognitive functions, 2
ANFIS modeling for prediction of particle motions in fluid flows
NASA Astrophysics Data System (ADS)
Safdari, Arman; Kim, Kyung Chun
2015-11-01
Accurate dynamic analysis of parcel of solid particles driven in fluid flow system is of interest for many natural and industrial applications such as sedimentation process, study of cloud particles in atmosphere, etc. In this paper, numerical modeling of solid particles in incompressible flow using Eulerian-Lagrangian approach is carried out to investigate the dynamic behavior of particles in different flow conditions; channel and cavity flow. Although modern computers have been well developed, the high computational time and costs for this kind of problems are still demanded. The Lattice Boltzmann Method (LBM) is used to simulate fluid flows and combined with the Lagrangian approach to predict the motion of particles in the range of masses. Some particles are selected, and subjected to Adaptive-network-based fuzzy inference system (ANFIS) to predict the trajectory of moving solid particles. Using a hybrid learning procedure from computational particle movement, the ANFIS can construct an input-output mapping based on fuzzy if-then rules and stipulated computational fluid dynamics prediction pairs. The obtained results from ANFIS algorithm is validated and compared with the set of benchmark data provided based on point-like approach coupled with the LBM method.
Goksel, Orcun; Zahiri-Azar, Reza; Salcudean, Septimiu E
2007-01-01
Motion estimation in sequences of ultrasound echo signals is essential for a wide range of applications. In time domain cross correlation, which is a common motion estimation technique, the displacements are typically not integral multiples of the sampling period. Therefore, to estimate the motion with sub-sample accuracy, 1D and 2D interpolation methods such as parabolic, cosine, and ellipsoid fitting have been introduced in the literature. In this paper, a simulation framework is presented in order to compare the performance of currently available techniques. First, the tissue deformation is modeled using the finite element method (FEM) and then the corresponding pre-/post-deformation radio-frequency (RF) signals are generated using Field II ultrasound simulation software. Using these simulated RF data of deformation, both axial and lateral tissue motion are estimated with sub-sample accuracy. The estimated displacements are then evaluated by comparing them to the known displacements computed by the FEM. This simulation approach was used to evaluate three different lateral motion estimation techniques employing (i) two separate 1D sub-sampling, (ii) two consecutive 1D sub-sampling, and (iii) 2D joint sub-sampling estimators. The estimation errors during two different tissue compression tests are presented with and without spatial filtering. Results show that RF signal processing methods involving tissue deformation can be evaluated using the proposed simulation technique, which employs accurate models. PMID:18002416
A new energy harvester for fluids in motion
NASA Astrophysics Data System (ADS)
Boragno, Corrado; Boccalero, Gregorio
2015-04-01
A new energy harvester, based on the fluttering phenomenon, is presented. The device is done with a wing connected to a support via two elastomers. When a fluid in motion impinges on this elastic structure, an amount of kinetic energy is transferred to the system, inducing large amplitude oscillations if few mechanical parameters are correctly set. In order to transform the mechanical energy in electrical energy, an electromagnetic coupling is adopted. In this way, it is possible to produce several mW in a wind of 4 m/s with a centimeter-sized device. The device is conceived as an autonomous power source for distributed sensors to be used in Internet of Things.
NASA Astrophysics Data System (ADS)
Kawamura, E.; Lichtenberg, A. J.; Lieberman, M. A.; Marakhtanov, A. M.
2016-06-01
A fast 2D axisymmetric fluid-analytical multifrequency capacitively coupled plasma (CCP) reactor code is used to study center high nonuniformity in a low pressure electronegative chlorine discharge. In the code, a time-independent Helmholtz wave equation is used to solve for the capacitive fields in the linearized frequency domain. This eliminates the time dependence from the electromagnetic (EM) solve, greatly speeding up the simulations at the cost of neglecting higher harmonics. However, since the code allows up to three driving frequencies, we can add the two most important harmonics to the CCP simulations as the second and third input frequencies. The amplitude and phase of these harmonics are estimated by using a recently developed 1D radial nonlinear transmission line (TL) model of a highly asymmetric cylindrical discharge (Lieberman et al 2015 Plasma Sources Sci. Technol. 24 055011). We find that at higher applied frequencies, the higher harmonics contribute significantly to the center high nonuniformity due to their shorter plasma wavelengths.
Bonanno, Gabriele; Puy, Gilles; Wiaux, Yves; van Heeswijk, Ruud B.; Piccini, Davide; Stuber, Matthias
2014-01-01
Purpose Respiratory motion correction remains a challenge in coronary magnetic resonance imaging (MRI) and current techniques, such as navigator gating, suffer from sub-optimal scan efficiency and ease-of-use. To overcome these limitations, an image-based self-navigation technique is proposed that uses “sub-images” and compressed sensing (CS) to obtain translational motion correction in 2D. The method was preliminarily implemented as a 2D technique and tested for feasibility for targeted coronary imaging. Methods During a 2D segmented radial k-space data acquisition, heavily undersampled sub-images were reconstructed from the readouts collected during each cardiac cycle. These sub-images may then be used for respiratory self-navigation. Alternatively, a CS reconstruction may be used to create these sub-images, so as to partially compensate for the heavy undersampling. Both approaches were quantitatively assessed using simulations and in vivo studies, and the resulting self-navigation strategies were then compared to conventional navigator gating. Results Sub-images reconstructed using CS showed a lower artifact level than sub-images reconstructed without CS. As a result, the final image quality was significantly better when using CS-assisted self-navigation as opposed to the non-CS approach. Moreover, while both self-navigation techniques led to a 69% scan time reduction (as compared to navigator gating), there was no significant difference in image quality between the CS-assisted self-navigation technique and conventional navigator gating, despite the significant decrease in scan time. Conclusions CS-assisted self-navigation using 2D translational motion correction demonstrated feasibility of producing coronary MRA data with image quality comparable to that obtained with conventional navigator gating, and does so without the use of additional acquisitions or motion modeling, while still allowing for 100% scan efficiency and an improved ease-of-use. In
NASA Astrophysics Data System (ADS)
Chen, K.; You, Y.; Noblesse, F.
2016-07-01
Experiments are conducted in a linear stratified fluid with a momentum source modeled via a nozzle jet moving horizontally. The generation mechanism of the quasi-two-dimensional dipolar vortex streets is investigated and their evolution characteristics are analyzed. Observation shows that the formation of a dipolar vortex street requires a nonzero motion of the nozzle in addition to conditions of the Reynolds and Froude number (Re, Fr). The (Re, Fr) condition that the dipolar vortex streets can be generated is determined via experimental measurements. The explanation for the absence of such a vortex street can be the low energy of the jet and the strong body-effect disturbance of the solid nozzle. The dependence of the vortex street dimensionless formation time τ and the Strouhal number St on the Froude number Fr or the Reynolds number Re is analyzed. This analysis shows that τ and St appear to be independent of Re and approximately have power-law relations with Fr via data fitting. The exponents of Fr in the two power-law functions are -0.27 for τ and -0.21 for St, while the constant coefficients are 65 and 0.21.
Spoerk, Jakob; Gendrin, Christelle; Weber, Christoph; Figl, Michael; Pawiro, Supriyanto Ardjo; Furtado, Hugo; Fabri, Daniella; Bloch, Christoph; Bergmann, Helmar; Gröller, Eduard; Birkfellner, Wolfgang
2012-01-01
A common problem in image-guided radiation therapy (IGRT) of lung cancer as well as other malignant diseases is the compensation of periodic and aperiodic motion during dose delivery. Modern systems for image-guided radiation oncology allow for the acquisition of cone-beam computed tomography data in the treatment room as well as the acquisition of planar radiographs during the treatment. A mid-term research goal is the compensation of tumor target volume motion by 2D/3D registration. In 2D/3D registration, spatial information on organ location is derived by an iterative comparison of perspective volume renderings, so-called digitally rendered radiographs (DRR) from computed tomography volume data, and planar reference x-rays. Currently, this rendering process is very time consuming, and real-time registration, which should at least provide data on organ position in less than a second, has not come into existence. We present two GPU-based rendering algorithms which generate a DRR of 512 × 512 pixels size from a CT dataset of 53 MB size at a pace of almost 100 Hz. This rendering rate is feasible by applying a number of algorithmic simplifications which range from alternative volume-driven rendering approaches – namely so-called wobbled splatting – to sub-sampling of the DRR-image by means of specialized raycasting techniques. Furthermore, general purpose graphics processing unit (GPGPU) programming paradigms were consequently utilized. Rendering quality and performance as well as the influence on the quality and performance of the overall registration process were measured and analyzed in detail. The results show that both methods are competitive and pave the way for fast motion compensation by rigid and possibly even non-rigid 2D/3D registration and, beyond that, adaptive filtering of motion models in IGRT. PMID:21782399
Study of the Motion of a Vertically Falling Sphere in a Viscous Fluid
ERIC Educational Resources Information Center
Soares, A. A.; Caramelo, L.; Andrade, M. A. P. M.
2012-01-01
This paper aims at contributing to a better understanding of the motion of spherical particles in viscous fluids. The classical problem of spheres falling through viscous fluids for small Reynolds numbers was solved taking into account the effects of added mass. The analytical solution for the motion of a falling sphere, from the beginning to the…
Shape matters: Near-field fluid mechanics dominate the collective motions of ellipsoidal squirmers.
Kyoya, K; Matsunaga, D; Imai, Y; Omori, T; Ishikawa, T
2015-12-01
Microswimmers show a variety of collective motions. Despite extensive study, questions remain regarding the role of near-field fluid mechanics in collective motion. In this paper, we describe precisely the Stokes flow around hydrodynamically interacting ellipsoidal squirmers in a monolayer suspension. The results showed that various collective motions, such as ordering, aggregation, and whirls, are dominated by the swimming mode and the aspect ratio. The collective motions are mainly induced by near-field fluid mechanics, despite Stokes flow propagation over a long range. These results emphasize the importance of particle shape in collective motion. PMID:26764823
NASA Astrophysics Data System (ADS)
Gao, Yan; Sui, Fusheng; Muggleton, Jennifer M.; Yang, Jun
2016-08-01
The dispersion characteristics of axisymmetric (n=0) waves offer a way to gain physical insight into the low-frequency vibrational behaviour of underground pipe systems. Whilst these can be found in the literature, they are generally calculated numerically. Coupled equations of motion for the n=0 waves that propagate in a buried fluid-filled pipe are presented in this paper and, from this, an analytical solution is developed for the fluid-dominated (s=1) wavenumber. The effect of the frictional stress at the pipe-soil interface on the dispersion behaviour of the s=1 wave is characterised by adopting a soil loading matrix. Overall, the fluid loading has a greater effect on the propagation wavespeed compared with the soil loading: for metal pipes, the effect of soil loading is negligible; for plastic pipes, however, simply neglecting the effect of soil loading can lead to a considerable underestimation in the calculation of the wavespeed. The wave attenuation increases significantly at higher frequencies regardless of pipe material resulting from the added damping due to radiation into the soil. Theoretical predictions of the s=1 wavenumber are compared with experimental data measured on an MDPE water pipe. The degree of agreement between prediction and experiment makes clear that, although the wavespeed is only slightly affected by the presence of the frictional stress, the frictional stress at the pipe-soil interface needs to be appropriately taken into account for attenuation predictions.
NASA Astrophysics Data System (ADS)
Haris, L.; Khotimah, S. N.; Haryanto, F.; Viridi, S.
2014-02-01
Molecular dynamics has been widely used to numerically solve equation of motion of classical many-particle system. It can be used to simulate many systems including biophysics, whose complexity level is determined by the involved elements. Based on this method, a numerical model had been constructed to mimic the behaviour of malaria-infected red blood cells within capillary vessel. The model was governed by three forces namely Coulomb force, normal force, and Stokes force. By utilizing two dimensional four-cells scheme, theoretical observation was carried out to test its capability. Although the parameters were chosen deliberately, all of the quantities were given arbitrary value. Despite this fact, the results were quite satisfactory. Combined with the previous results, it can be said that the proposed model were sufficient enough to mimic the malaria-infected red blood cells motion within obstructed capillary vessel.
Real-time ultrasound-tagging to track the 2D motion of the common carotid artery wall in vivo
Zahnd, Guillaume; Salles, Sébastien; Liebgott, Hervé; Vray, Didier; Sérusclat, André; Moulin, Philippe
2015-02-15
Purpose: Tracking the motion of biological tissues represents an important issue in the field of medical ultrasound imaging. However, the longitudinal component of the motion (i.e., perpendicular to the beam axis) remains more challenging to extract due to the rather coarse resolution cell of ultrasound scanners along this direction. The aim of this study is to introduce a real-time beamforming strategy dedicated to acquire tagged images featuring a distinct pattern in the objective to ease the tracking. Methods: Under the conditions of the Fraunhofer approximation, a specific apodization function was applied to the received raw channel data, in real-time during image acquisition, in order to introduce a periodic oscillations pattern along the longitudinal direction of the radio frequency signal. Analytic signals were then extracted from the tagged images, and subpixel motion tracking of the intima–media complex was subsequently performed offline, by means of a previously introduced bidimensional analytic phase-based estimator. Results: The authors’ framework was applied in vivo on the common carotid artery from 20 young healthy volunteers and 6 elderly patients with high atherosclerosis risk. Cine-loops of tagged images were acquired during three cardiac cycles. Evaluated against reference trajectories manually generated by three experienced analysts, the mean absolute tracking error was 98 ± 84 μm and 55 ± 44 μm in the longitudinal and axial directions, respectively. These errors corresponded to 28% ± 23% and 13% ± 9% of the longitudinal and axial amplitude of the assessed motion, respectively. Conclusions: The proposed framework enables tagged ultrasound images of in vivo tissues to be acquired in real-time. Such unconventional beamforming strategy contributes to improve tracking accuracy and could potentially benefit to the interpretation and diagnosis of biomedical images.
NASA Astrophysics Data System (ADS)
Velioǧlu, Deniz; Cevdet Yalçıner, Ahmet; Zaytsev, Andrey
2016-04-01
Tsunamis are huge waves with long wave periods and wave lengths that can cause great devastation and loss of life when they strike a coast. The interest in experimental and numerical modeling of tsunami propagation and inundation increased considerably after the 2011 Great East Japan earthquake. In this study, two numerical codes, FLOW 3D and NAMI DANCE, that analyze tsunami propagation and inundation patterns are considered. Flow 3D simulates linear and nonlinear propagating surface waves as well as long waves by solving three-dimensional Navier-Stokes (3D-NS) equations. NAMI DANCE uses finite difference computational method to solve 2D depth-averaged linear and nonlinear forms of shallow water equations (NSWE) in long wave problems, specifically tsunamis. In order to validate these two codes and analyze the differences between 3D-NS and 2D depth-averaged NSWE equations, two benchmark problems are applied. One benchmark problem investigates the runup of long waves over a complex 3D beach. The experimental setup is a 1:400 scale model of Monai Valley located on the west coast of Okushiri Island, Japan. Other benchmark problem is discussed in 2015 National Tsunami Hazard Mitigation Program (NTHMP) Annual meeting in Portland, USA. It is a field dataset, recording the Japan 2011 tsunami in Hilo Harbor, Hawaii. The computed water surface elevation and velocity data are compared with the measured data. The comparisons showed that both codes are in fairly good agreement with each other and benchmark data. The differences between 3D-NS and 2D depth-averaged NSWE equations are highlighted. All results are presented with discussions and comparisons. Acknowledgements: Partial support by Japan-Turkey Joint Research Project by JICA on earthquakes and tsunamis in Marmara Region (JICA SATREPS - MarDiM Project), 603839 ASTARTE Project of EU, UDAP-C-12-14 project of AFAD Turkey, 108Y227, 113M556 and 213M534 projects of TUBITAK Turkey, RAPSODI (CONCERT_Dis-021) of CONCERT
NASA Astrophysics Data System (ADS)
Rank, Christopher M.; Heußer, Thorsten; Flach, Barbara; Brehm, Marcus; Kachelrieß, Marc
2015-03-01
We propose a new method for PET/MR respiratory motion compensation, which is based on a 3D-2D registration of strongly undersampled MR data and a) runs in parallel with the PET acquisition, b) can be interlaced with clinical MR sequences, and c) requires less than one minute of the total MR acquisition time per bed position. In our simulation study, we applied a 3D encoded radial stack-of-stars sampling scheme with 160 radial spokes per slice and an acquisition time of 38 s. Gated 4D MR images were reconstructed using a 4D iterative reconstruction algorithm. Based on these images, motion vector fields were estimated using our newly-developed 3D-2D registration framework. A 4D PET volume of a patient with eight hot lesions in the lungs and upper abdomen was simulated and MoCo 4D PET images were reconstructed based on the motion vector fields derived from MR. For evaluation, average SUVmean values of the artificial lesions were determined for a 3D, a gated 4D, a MoCo 4D and a reference (with ten-fold measurement time) gated 4D reconstruction. Compared to the reference, 3D reconstructions yielded an underestimation of SUVmean values due to motion blurring. In contrast, gated 4D reconstructions showed the highest variation of SUVmean due to low statistics. MoCo 4D reconstructions were only slightly affected by these two sources of uncertainty resulting in a significant visual and quantitative improvement in terms of SUVmean values. Whereas temporal resolution was comparable to the gated 4D images, signal-to-noise ratio and contrast-to-noise ratio were close to the 3D reconstructions.
CFD simulation of the vertical motion characteristics of the moonpool fluid for the truss spar
NASA Astrophysics Data System (ADS)
Wang, Bin; Liu, Liqin; Tang, Yougang
2014-03-01
The research purpose of this paper is to estimate the impacts of the parameters of the guide plate on the vertical motion characteristics of the moonpool fluid. With the volume of fluid (VOF) method, three-dimensional models of the moonpool fluid motions of the truss spar platform are established. Simulation results are then presented for the moonpool forced oscillation by employing the dynamic mesh method and user-defined functions in FLUENT. The motions of the moonpool fluid and the loads on the guide plates are obtained for both cases of square-ring and crisscross. The results show that the shape and area of the guide plate at the bottom of the moonpool have a significant impact on the physical parameters of the moonpool, including the load on the moonpool guide plate, motion form of the moonpool fluid and the mass flow rate.
NASA Astrophysics Data System (ADS)
Furtado, H.; Steiner, E.; Stock, M.; Georg, D.; Birkfellner, W.
2014-03-01
Intra-fractional respiratorymotion during radiotherapy is one of themain sources of uncertainty in dose application creating the need to extend themargins of the planning target volume (PTV). Real-time tumormotion tracking by 2D/3D registration using on-board kilo-voltage (kV) imaging can lead to a reduction of the PTV. One limitation of this technique when using one projection image, is the inability to resolve motion along the imaging beam axis. We present a retrospective patient study to investigate the impact of paired portal mega-voltage (MV) and kV images, on registration accuracy. We used data from eighteen patients suffering from non small cell lung cancer undergoing regular treatment at our center. For each patient we acquired a planning CT and sequences of kV and MV images during treatment. Our evaluation consisted of comparing the accuracy of motion tracking in 6 degrees-of-freedom(DOF) using the anterior-posterior (AP) kV sequence or the sequence of kV-MV image pairs. We use graphics processing unit rendering for real-time performance. Motion along cranial-caudal direction could accurately be extracted when using only the kV sequence but in AP direction we obtained large errors. When using kV-MV pairs, the average error was reduced from 3.3 mm to 1.8 mm and the motion along AP was successfully extracted. The mean registration time was of 190+/-35ms. Our evaluation shows that using kVMV image pairs leads to improved motion extraction in 6 DOF. Therefore, this approach is suitable for accurate, real-time tumor motion tracking with a conventional LINAC.
Feng, Yongqiang; Max, Ludo
2014-01-01
Purpose Studying normal or disordered motor control requires accurate motion tracking of the effectors (e.g., orofacial structures). The cost of electromagnetic, optoelectronic, and ultrasound systems is prohibitive for many laboratories, and limits clinical applications. For external movements (lips, jaw), video-based systems may be a viable alternative, provided that they offer high temporal resolution and sub-millimeter accuracy. Method We examined the accuracy and precision of 2D and 3D data recorded with a system that combines consumer-grade digital cameras capturing 60, 120, or 240 frames per second (fps), retro-reflective markers, commercially-available computer software (APAS, Ariel Dynamics), and a custom calibration device. Results Overall mean error (RMSE) across tests was 0.15 mm for static tracking and 0.26 mm for dynamic tracking, with corresponding precision (SD) values of 0.11 and 0.19 mm, respectively. The effect of frame rate varied across conditions, but, generally, accuracy was reduced at 240 fps. The effect of marker size (3 vs. 6 mm diameter) was negligible at all frame rates for both 2D and 3D data. Conclusion Motion tracking with consumer-grade digital cameras and the APAS software can achieve sub-millimeter accuracy at frame rates that are appropriate for kinematic analyses of lip/jaw movements for both research and clinical purposes. PMID:24686484
Feng, Yongqiang; Max, Ludo
2014-04-01
PURPOSE Studying normal or disordered motor control requires accurate motion tracking of the effectors (e.g., orofacial structures). The cost of electromagnetic, optoelectronic, and ultrasound systems is prohibitive for many laboratories and limits clinical applications. For external movements (lips, jaw), video-based systems may be a viable alternative, provided that they offer high temporal resolution and submillimeter accuracy. METHOD The authors examined the accuracy and precision of 2-D and 3-D data recorded with a system that combines consumer-grade digital cameras capturing 60, 120, or 240 frames per second (fps), retro-reflective markers, commercially available computer software (APAS, Ariel Dynamics), and a custom calibration device. RESULTS Overall root-mean-square error (RMSE) across tests was 0.15 mm for static tracking and 0.26 mm for dynamic tracking, with corresponding precision (SD) values of 0.11 and 0.19 mm, respectively. The effect of frame rate varied across conditions, but, generally, accuracy was reduced at 240 fps. The effect of marker size (3- vs. 6-mm diameter) was negligible at all frame rates for both 2-D and 3-D data. CONCLUSION Motion tracking with consumer-grade digital cameras and the APAS software can achieve submillimeter accuracy at frame rates that are appropriate for kinematic analyses of lip/jaw movements for both research and clinical purposes. PMID:24686484
Theory and Validation of Magnetic Resonance Fluid Motion Estimation Using Intensity Flow Data
Wong, Kelvin Kian Loong; Kelso, Richard Malcolm; Worthley, Stephen Grant; Sanders, Prashanthan; Mazumdar, Jagannath; Abbott, Derek
2009-01-01
Background Motion tracking based on spatial-temporal radio-frequency signals from the pixel representation of magnetic resonance (MR) imaging of a non-stationary fluid is able to provide two dimensional vector field maps. This supports the underlying fundamentals of magnetic resonance fluid motion estimation and generates a new methodology for flow measurement that is based on registration of nuclear signals from moving hydrogen nuclei in fluid. However, there is a need to validate the computational aspect of the approach by using velocity flow field data that we will assume as the true reference information or ground truth. Methodology/Principal Findings In this study, we create flow vectors based on an ideal analytical vortex, and generate artificial signal-motion image data to verify our computational approach. The analytical and computed flow fields are compared to provide an error estimate of our methodology. The comparison shows that the fluid motion estimation approach using simulated MR data is accurate and robust enough for flow field mapping. To verify our methodology, we have tested the computational configuration on magnetic resonance images of cardiac blood and proved that the theory of magnetic resonance fluid motion estimation can be applicable practically. Conclusions/Significance The results of this work will allow us to progress further in the investigation of fluid motion prediction based on imaging modalities that do not require velocity encoding. This article describes a novel theory of motion estimation based on magnetic resonating blood, which may be directly applied to cardiac flow imaging. PMID:19270756
This study is a part of an ongoing research project that aims at assessing the environmental benefits of DNAPL removal. The laboratory part of the research project is to examine the functional relationship between DNAPL architecture, mass removal and contaminant mass flux in 2-D ...
Motion-dependent fluid forces acting on tube arrays in crossflow
Chen, S.S.; Zhu, S.; Jendrzejczyk, J.A.
1993-06-01
Motion-dependent fluids forces acting on a tube array were measured as a function of excitation frequency, excitation amplitude, and flow velocity. Fluid-damping and fluid-stiffness coefficients were obtained from measured motion-dependent fluid forces as a function of reduced flow velocity and excitation amplitude. The water channel and test setup provide a sound facility for obtaining key coefficients for fluidelastic instability of tube arrays in crossflow. Once the guideline, based on the unsteady flow theory, can be developed for fluidelastic instability of tube arrays crossflow.
Wang, Zhirui; Xu, Jia; Huang, Zuzhen; Zhang, Xudong; Xia, Xiang-Gen; Long, Teng; Bao, Qian
2016-01-01
To detect and estimate ground slowly moving targets in airborne single-channel synthetic aperture radar (SAR), a road-aided ground moving target indication (GMTI) algorithm is proposed in this paper. First, the road area is extracted from a focused SAR image based on radar vision. Second, after stationary clutter suppression in the range-Doppler domain, a moving target is detected and located in the image domain via the watershed method. The target’s position on the road as well as its radial velocity can be determined according to the target’s offset distance and traffic rules. Furthermore, the target’s azimuth velocity is estimated based on the road slope obtained via polynomial fitting. Compared with the traditional algorithms, the proposed method can effectively cope with slowly moving targets partly submerged in a stationary clutter spectrum. In addition, the proposed method can be easily extended to a multi-channel system to further improve the performance of clutter suppression and motion estimation. Finally, the results of numerical experiments are provided to demonstrate the effectiveness of the proposed algorithm. PMID:26999140
Wang, Zhirui; Xu, Jia; Huang, Zuzhen; Zhang, Xudong; Xia, Xiang-Gen; Long, Teng; Bao, Qian
2016-01-01
To detect and estimate ground slowly moving targets in airborne single-channel synthetic aperture radar (SAR), a road-aided ground moving target indication (GMTI) algorithm is proposed in this paper. First, the road area is extracted from a focused SAR image based on radar vision. Second, after stationary clutter suppression in the range-Doppler domain, a moving target is detected and located in the image domain via the watershed method. The target's position on the road as well as its radial velocity can be determined according to the target's offset distance and traffic rules. Furthermore, the target's azimuth velocity is estimated based on the road slope obtained via polynomial fitting. Compared with the traditional algorithms, the proposed method can effectively cope with slowly moving targets partly submerged in a stationary clutter spectrum. In addition, the proposed method can be easily extended to a multi-channel system to further improve the performance of clutter suppression and motion estimation. Finally, the results of numerical experiments are provided to demonstrate the effectiveness of the proposed algorithm. PMID:26999140
Exploiting the color of Brownian motion for high-frequency microrheology of Newtonian fluids
NASA Astrophysics Data System (ADS)
Domínguez-García, Pablo; Mor, Flavio M.; Forró, László; Jeney, Sylvia
2013-09-01
Einstein's stochastic description of the random movement of small objects in a fluid, i.e. Brownian motion, reveals to be quite different, when observed on short timescales. The limitations of Einstein's theory with respect to particle inertia and hydrodynamic memory yield to the apparition of a colored frequency-dependent component in the spectrum of the thermal forces, which is called "the color of Brownian motion". The knowledge of the characteristic timescales of the motion of a trapped microsphere motion in a Newtonian fluid allowed to develop a high-resolution calibration method for optical interferometry. Well-calibrated correlation quantities, such as the mean square displacement or the velocity autocorrelation function, permit to study the mechanical properties of fluids at high frequencies. These properties are estimated by microrheological calculations based on the theoretical relations between the complex mobility of the beads and the rheological properties of a complex fluid.
Motion of a viscous fluid and a wall in the presence of a stationary wall
NASA Astrophysics Data System (ADS)
Sennitskii, V. L.
2016-03-01
The problem of damping motion of a hydromechanical system consisting of a viscous fluid and its bounding rigid walls is solved. A condition under which there is an abrupt deceleration of the hydromechanical system is determined.
Investigating the fluid mechanics behind red blood cell-induced lateral platelet motion
NASA Astrophysics Data System (ADS)
Crowl Erickson, Lindsay; Fogelson, Aaron
2009-11-01
Platelets play an essential role in blood clotting; they adhere to damaged tissue and release chemicals that activate other platelets. Yet in order to adhere, platelets must first come into contact with the injured vessel wall. Under arterial flow conditions, platelets have an enhanced concentration near blood vessel walls. This non-uniform cell distribution depends on the fluid dynamics of blood as a heterogeneous medium. We use a parallelized lattice Boltzmann-immersed boundary method to solve the flow dynamics of red cells and platelets in a periodic 2D vessel with no-slip boundary conditions. Red cells are treated as biconcave immersed boundary objects with isotropic Skalak membrane tension and an internal viscosity five times that of the surrounding plasma. Using this method we analyze the influence of shear rate, hematocrit, and red cell membrane properties on lateral platelet motion. We find that the effective diffusion of platelets is significantly lower near the vessel wall compared to the center of the vessel. Insight gained from this work could lead to significant improvements to current models for platelet adhesion where the presence of red blood cells is neglected due to computational intensity.
Incompressible wave motion of inhomogeneous, compressible fluids in a gravity field
NASA Astrophysics Data System (ADS)
Godin, O. A.
2012-04-01
We consider a particular class of linear and non-linear wave motions in fluids, in which pressure remains constant in each moving fluid parcel. The fluid is assumed to be inviscid, and wave motion is considered as an adiabatic thermodynamic process. An exact, analytic solution of linearized hydrodynamics equations is obtained that describes the wave motion in inhomogeneous, compressible, rotating fluids with piece-wise continuous parameters in a uniform gravity field. The solution is valid under surprisingly general assumptions about the environment and reduces to some classical wave types in appropriate limiting cases. Free waves in bounded and unbounded domains as well as excitation of wave fields by a point source are considered. Edge waves propagating along vertical and inclined rigid boundaries are found in rotating and non-rotating fluids. Allowance for three-dimensional variation of the sound speed and for arbitrary density stratification, including density discontinuities, makes the exact solution an attractive model of acoustic-gravity waves in a coupled ocean-atmosphere system. The new wave type complements classical exact solutions of linearized equations of fluid mechanics known as the Rossby, Lamb, Kelvin, and Poincaré waves, which provide much of the conceptual foundation of geophysical fluid dynamics. In addition to a wide class of exact solutions for linear waves, an exact solution of full non-linear hydrodynamics equations is found that describes a propagating wave in inhomogeneous, compressible fluids with piece-wise continuous parameters in a uniform gravity field. The fluid may have a free surface and a rigid boundary. Depending on the geometry of the problem, the solution has the meaning of either surface or edge wave. The exact solution describes a finite-amplitude wave in an otherwise quiescent fluid. Extensions to finite-amplitude waves in fluids with background currents are considered. Relation of the new exact solution for the non
NASA Astrophysics Data System (ADS)
Lopes Filho, Milton C.; Nussenzveig Lopes, Helena J.; Titi, Edriss S.; Zang, Aibin
2015-06-01
The second-grade fluid equations are a model for viscoelastic fluids, with two parameters: α > 0, corresponding to the elastic response, and , corresponding to viscosity. Formally setting these parameters to 0 reduces the equations to the incompressible Euler equations of ideal fluid flow. In this article we study the limits of solutions of the second-grade fluid system, in a smooth, bounded, two-dimensional domain with no-slip boundary conditions. This class of problems interpolates between the Euler- α model (), for which the authors recently proved convergence to the solution of the incompressible Euler equations, and the Navier-Stokes case ( α = 0), for which the vanishing viscosity limit is an important open problem. We prove three results. First, we establish convergence of the solutions of the second-grade model to those of the Euler equations provided , as α → 0, extending the main result in (Lopes Filho et al., Physica D 292(293):51-61, 2015). Second, we prove equivalence between convergence (of the second-grade fluid equations to the Euler equations) and vanishing of the energy dissipation in a suitably thin region near the boundary, in the asymptotic regime , as α → 0. This amounts to a convergence criterion similar to the well-known Kato criterion for the vanishing viscosity limit of the Navier-Stokes equations to the Euler equations. Finally, we obtain an extension of Kato's classical criterion to the second-grade fluid model, valid if , as . The proof of all these results relies on energy estimates and boundary correctors, following the original idea by Kato.
Hou, Gary Y.; Provost, Jean; Grondin, Julien; Wang, Shutao; Marquet, Fabrice; Bunting, Ethan; Konofagou, Elisa E.
2015-01-01
Harmonic Motion Imaging for Focused Ultrasound (HMIFU) is a recently developed High-Intensity Focused Ultrasound (HIFU) treatment monitoring method. HMIFU utilizes an Amplitude-Modulated (fAM = 25 Hz) HIFU beam to induce a localized focal oscillatory motion, which is simultaneously estimated and imaged by confocally-aligned imaging transducer. HMIFU feasibilities have been previously shown in silico, in vitro, and in vivo in 1-D or 2-D monitoring of HIFU treatment. The objective of this study is to develop and show the feasibility of a novel fast beamforming algorithm for image reconstruction using GPU-based sparse-matrix operation with real-time feedback. In this study, the algorithm was implemented onto a fully integrated, clinically relevant HMIFU system composed of a 93-element HIFU transducer (fcenter = 4.5MHz) and coaxially-aligned 64-element phased array (fcenter = 2.5MHz) for displacement excitation and motion estimation, respectively. A single transmit beam with divergent beam transmit was used while fast beamforming was implemented using a GPU-based delay-and-sum method and a sparse-matrix operation. Axial HMI displacements were then estimated from the RF signals using a 1-D normalized cross-correlation method and streamed to a graphic user interface. The present work developed and implemented a sparse matrix beamforming onto a fully-integrated, clinically relevant system, which can stream displacement images up to 15 Hz using a GPU-based processing, an increase of 100 fold in rate of streaming displacement images compared to conventional CPU-based conventional beamforming and reconstruction processing. The achieved feedback rate is also currently the fastest and only approach that does not require interrupting the HIFU treatment amongst the acoustic radiation force based HIFU imaging techniques. Results in phantom experiments showed reproducible displacement imaging, and monitoring of twenty two in vitro HIFU treatments using the new 2D system showed a
An IPOT meshless method using DC PSE approximation for fluid flow equations in 2D and 3D geometries
NASA Astrophysics Data System (ADS)
Bourantas, G. C.; Loukopoulos, V. C.; Skouras, E. D.; Burganos, V. N.; Nikiforidis, G. C.
2016-06-01
Navier-Stokes (N-S) equations, in their primitive variable (u-v-p) formulation, are numerically solved using the Implicit Potential (IPOT) numerical scheme in the context of strong form Meshless Point Collocation (MPC) method. The unknown field functions are computed using the Discretization Correction Particle Strength Exchange (DC PSE) approximation method. The latter makes use of discrete moment conditions to derive the operator kernels, which leads to low condition number for the moment matrix compared to other meshless interpolation methods and increased stability for the numerical solution. The proposed meshless scheme is applied on 2D and 3D spatial domains, using uniform or irregular set of nodes to represent the domain. The numerical results obtained are compared against those obtained using well-established methods.
Motion of spheres along a fluid-gas interface.
Cichocki, Bogdan; Ekiel-Jezewska, Maria L; Nägele, Gerhard; Wajnryb, Eligiusz
2004-08-01
A system of many spherical particles, suspended in a quiescent fluid and touching a planar free fluid-gas interface, is considered. Stick fluid boundary conditions at the sphere surfaces are assumed. The free surface boundary conditions are taken into account with the use of the method of images. For such a quasi-two-dimensional system, the one-sphere resistance operator is calculated numerically. Moreover, the corresponding friction and mobility tensors are constructed from irreducible multipole expansion. Finally, the long-distance terms of the two-sphere mobility tensor are evaluated explicitly up to the order of 1/r3, where r is the interparticle distance. Experiments which have motivated this work are outlined. PMID:15260785
Forward-in-time differencing for fluids: Nonhydrostatic modeling of fluid motions on a sphere
Smolarkiewicz, P.K.; Grubisic, V.; Margolin, L.G.; Wyszogrodzki, A.A.
1998-12-31
Traditionally, numerical models for simulating planetary scale weather and climate employ the hydrostatic primitive equations--an abbreviated form of Navier-Stokes` equations that neglect vertical accelerations and use simplified Coriolis forces. Although there is no evidence so far that including nonhydrostatic effects in global models has any physical significance for large scale solutions, there is an emerging trend in the community toward restoring Navier-Stokes` equations (or at least their less constrained forms) in global models of atmospheres and oceans. The primary motivation is that state-of-the-art computers already admit resolutions where local nonhydrostatic effects become noticeable. much of this present research aims to improve the design of a high-performance numerical model for simulating the flows of moist (and precipitating), rotating, stratified fluids past a specified time-dependent irregular lower boundary. This model is representative of a class of nonhydrostatic atmospheric codes that employs the anelastic equations of motion in a terrain-following curvilinear framework, and contains parallel implementations of semi-Lagrangian and Eulerian approximations selectable by the user. The model has been employed in a variety of application; the quality of results suggest that modern nonoscillatory forward-in-time (NFT) methods are superior to the more traditional centered-in-time-and-space schemes, in terms of accuracy, computational efficiency, flexibility and robustness. The authors have extended the Cartesian NFT model to a mountainous sphere and, consequently, have dispensed with the traditional geophysical simplifications of hydrostaticity, gentle terrain slopes, and weak rotation. In this paper, they discuss the algorithmic design, relative efficiency and accuracy of several different variants (hydrostatic, nonhydrostatic, implicit, explicit, etc.) of the NFT global model. They substantiate their theoretical discussions with the results of
Motion of a bubble ring in a viscous fluid
NASA Astrophysics Data System (ADS)
Cheng, M.; Lou, J.; Lim, T. T.
2013-06-01
In this paper, lattice Boltzmann method was undertaken to study the dynamics of a vortex ring bubble (or bubble ring) in a viscous incompressible fluid. The study is motivated partly by our desire to assess whether a bubble ring keeps increasing its radius and decreasing its rise velocity as it rises through fluid as was predicted by Turner ["Buoyant vortex rings," Proc. R. Soc. London, Ser. A 239, 61 (1957)], 10.1098/rspa.1957.0022 and Pedley ["The toroidal bubble," J. Fluid Mech. 32, 97 (1968)], 10.1017/S0022112068000601, or does the ring like a rising bubble, eventually reaches a steady state where its radius and velocity remain constant as was predicted by Joseph et al. [Potential Flows of Viscous and Viscoelastic Fluids (Cambridge University Press, 2008)]. The parameters investigated included ring circulation, Reynolds number, density ratio and Bond number. Our numerical results show that a rising bubble ring increases its radius and decreases its velocity, but the process is interrupted by ring instability that eventually causes it to break up into smaller bubbles. This finding is consistent with the stability analysis by Pedley, who predicted that a bubble ring has a finite lifespan and is ultimately destroyed by surface tension instability. Furthermore, it is found that increasing initial circulation has a stabilizing effect on a bubble ring while increasing Reynolds number or Bond number hastens ring instability, resulting in an earlier break up into smaller bubbles; the number of bubbles depends on the wavenumber of the perturbation.
Fluid forces on rotating centrifugal impeller with whirling motion
NASA Technical Reports Server (NTRS)
Shoji, H.; Ohashi, H.
1980-01-01
Fluid forces on a centrifugal impeller, whose rotating axis whirls with a constant speed, were calculated by using unsteady potential theory. Calculations were performed for various values of whirl speed, number of impeller blades and angle of blades. Specific examples as well as significant results are given.
Near field fluid coupling between internal motion of the organ of Corti and the basilar membrane
Elliott, Stephen J.; Ni, Guangjian
2015-12-31
The pressure distribution in each of the fluid chambers of the cochlea can be decomposed into a 1D, or plane wave, component and a near field component, which decays rapidly away from the excitation point. The transverse motion of the basilar membrane, BM, for example, generates both a 1D pressure field, which couples into the slow wave, and a local near field pressure, proportional to the BM acceleration, that generates an added mass on the BM due to the fluid motion. When the organ of Corti, OC, undergoes internal motion, due for example to outer hair cell activity, this motion will not itself generate any 1D pressure if the OC is incompressible and the BM is constrained not to move volumetrically, and so will not directly couple into the slow wave. This motion will, however, generate a near field pressure, proportional to the OC acceleration, which will act on the OC and thus increases its effective mass. The near field pressure due to this OC motion will also act on the BM, generating a force on the BM proportional to the acceleration of the OC, and thus create a “coupling mass” effect. By reciprocity, this coupling mass is the same as that acting on the OC due to the motion of the BM. This near field fluid coupling is initially observed in a finite element model of a slice of the cochlea. These simulations suggest a simple analytical formulation for the fluid coupling, using higher order beam modes across the width of the cochlear partition. It is well known that the added mass due to the near field pressure dominates the overall mass of the BM, and thus significantly affects the micromechanical dynamics. This work not only quantifies the added mass of the OC due its own motion in the fluid, and shows that this is important, but also demonstrates that the coupling mass effect between the BM and OC significantly affects the dynamics of simple micromechanical models.
Near field fluid coupling between internal motion of the organ of Corti and the basilar membrane
NASA Astrophysics Data System (ADS)
Elliott, Stephen J.; Ni, Guangjian
2015-12-01
The pressure distribution in each of the fluid chambers of the cochlea can be decomposed into a 1D, or plane wave, component and a near field component, which decays rapidly away from the excitation point. The transverse motion of the basilar membrane, BM, for example, generates both a 1D pressure field, which couples into the slow wave, and a local near field pressure, proportional to the BM acceleration, that generates an added mass on the BM due to the fluid motion. When the organ of Corti, OC, undergoes internal motion, due for example to outer hair cell activity, this motion will not itself generate any 1D pressure if the OC is incompressible and the BM is constrained not to move volumetrically, and so will not directly couple into the slow wave. This motion will, however, generate a near field pressure, proportional to the OC acceleration, which will act on the OC and thus increases its effective mass. The near field pressure due to this OC motion will also act on the BM, generating a force on the BM proportional to the acceleration of the OC, and thus create a "coupling mass" effect. By reciprocity, this coupling mass is the same as that acting on the OC due to the motion of the BM. This near field fluid coupling is initially observed in a finite element model of a slice of the cochlea. These simulations suggest a simple analytical formulation for the fluid coupling, using higher order beam modes across the width of the cochlear partition. It is well known that the added mass due to the near field pressure dominates the overall mass of the BM, and thus significantly affects the micromechanical dynamics. This work not only quantifies the added mass of the OC due its own motion in the fluid, and shows that this is important, but also demonstrates that the coupling mass effect between the BM and OC significantly affects the dynamics of simple micromechanical models.
Teleparallelism, Brownian Motion, Quantum Mechanics and Fluid-Dynamics I
NASA Astrophysics Data System (ADS)
Rapoport, Diego
2002-12-01
Extending the rules of teleparallelism for the introduction of a metric and a connection with torsion on a smooth manifold, M, we define generalized Brownian motions on M starting with a standard Wiener process. The laplacian operator generating this diffusion is the square of the teleparallelism connection on M, yet it is found to depend on the trace-torsion, and thus we restrict to Riemann-Cartan-Weyl connections. We extend these constructions to the generalized Brownian motions of differential forms. We apply this to give random covariant implicit solutions of the Navier-Stokes equations. We give the constitutive equations for the trace-torsion Q, and obtain a non-linear wave equation with quantum potential term for a scalar ψ appearing in the term d lnψ of Q. We relate the diffusion with drift ∇lnψ, to the heat kernel of quantum gravity for a scalar field. In Q appear two electromagnetic potentials which are proved to produce the time-evolution irreversibility of the Brownian motions. They appear related to the rotational degrees of freedom of a massive non-linear Dirac-Hestenes spinor field which defines a global spinor structure on M and a solution of the Clifford-Maxwell equation.
Laminar Motion of the Incompressible Fluids in Self-Acting Thrust Bearings with Spiral Grooves
Velescu, Cornel; Popa, Nicolae Calin
2014-01-01
We analyze the laminar motion of incompressible fluids in self-acting thrust bearings with spiral grooves with inner or external pumping. The purpose of the study is to find some mathematical relations useful to approach the theoretical functionality of these bearings having magnetic controllable fluids as incompressible fluids, in the presence of a controllable magnetic field. This theoretical study approaches the permanent motion regime. To validate the theoretical results, we compare them to some experimental results presented in previous papers. The laminar motion of incompressible fluids in bearings is described by the fundamental equations of fluid dynamics. We developed and particularized these equations by taking into consideration the geometrical and functional characteristics of these hydrodynamic bearings. Through the integration of the differential equation, we determined the pressure and speed distributions in bearings with length in the “pumping” direction. These pressure and speed distributions offer important information, both quantitative (concerning the bearing performances) and qualitative (evidence of the viscous-inertial effects, the fluid compressibility, etc.), for the laminar and permanent motion regime. PMID:24526896
NASA Astrophysics Data System (ADS)
Kim, Ho Jun; Lee, Hae June
2016-06-01
The wide applicability of capacitively coupled plasma (CCP) deposition has increased the interest in developing comprehensive numerical models, but CCP imposes a tremendous computational cost when conducting a transient analysis in a three-dimensional (3D) model which reflects the real geometry of reactors. In particular, the detailed flow features of reactive gases induced by 3D geometric effects need to be considered for the precise calculation of radical distribution of reactive species. Thus, an alternative inclusive method for the numerical simulation of CCP deposition is proposed to simulate a two-dimensional (2D) CCP model based on the 3D gas flow results by simulating flow, temperature, and species fields in a 3D space at first without calculating the plasma chemistry. A numerical study of a cylindrical showerhead-electrode CCP reactor was conducted for particular cases of SiH4/NH3/N2/He gas mixture to deposit a hydrogenated silicon nitride (SiN x H y ) film. The proposed methodology produces numerical results for a 300 mm wafer deposition reactor which agree very well with the deposition rate profile measured experimentally along the wafer radius.
Determining effects of turbine blades on fluid motion
Linn, Rodman Ray; Koo, Eunmo
2012-05-01
Disclosed is a technique for simulating wind interaction with wind turbines. A turbine blade is divided into radial sections. The effect that each of these radial sections has on the velocities in Eulerian computational cells they overlap is determined. The effect is determined using Lagrangian techniques such that the calculations need not include wind components in the radial direction. A force on each radial section of turbine blade is determined. This force depends on the axial and azimuthal components of the fluid flow in the computational cell and the geometric properties of the turbine blade. The force on the turbine blade is fed back to effect the fluid flow in the computational cell for the next time step.
Determining effects of turbine blades on fluid motion
Linn, Rodman Ray; Koo, Eunmo
2011-05-31
Disclosed is a technique for simulating wind interaction with wind turbines. A turbine blade is divided into radial sections. The effect that each of these radial sections has on the velocities in Eulerian computational cells they overlap is determined. The effect is determined using Lagrangian techniques such that the calculations need not include wind components in the radial direction. A force on each radial section of turbine blade is determined. This force depends on the axial and azimuthal components of the fluid flow in the computational cell and the geometric properties of the turbine blade. The force on the turbine blade is fed back to effect the fluid flow in the computational cell for the next time step.
NASA Astrophysics Data System (ADS)
Yang, Xiaochen; Zhang, Qinghe; Hao, Linnan
2015-03-01
A water-fluid mud coupling model is developed based on the unstructured grid finite volume coastal ocean model (FVCOM) to investigate the fluid mud motion. The hydrodynamics and sediment transport of the overlying water column are solved using the original three-dimensional ocean model. A horizontal two-dimensional fluid mud model is integrated into the FVCOM model to simulate the underlying fluid mud flow. The fluid mud interacts with the water column through the sediment flux, current, and shear stress. The friction factor between the fluid mud and the bed, which is traditionally determined empirically, is derived with the assumption that the vertical distribution of shear stress below the yield surface of fluid mud is identical to that of uniform laminar flow of Newtonian fluid in the open channel. The model is validated by experimental data and reasonable agreement is found. Compared with numerical cases with fixed friction factors, the results simulated with the derived friction factor exhibit the best agreement with the experiment, which demonstrates the necessity of the derivation of the friction factor.
Two-Phase Acto-Cytosolic Fluid Flow in a Moving Keratocyte: A 2D Continuum Model.
Nikmaneshi, M R; Firoozabadi, B; Saidi, M S
2015-09-01
The F-actin network and cytosol in the lamellipodia of crawling cells flow in a centripetal pattern and spout-like form, respectively. We have numerically studied this two-phase flow in the realistic geometry of a moving keratocyte. Cytosol has been treated as a low viscosity Newtonian fluid flowing through the high viscosity porous medium of F-actin network. Other involved phenomena including myosin activity, adhesion friction, and interphase interaction are also discussed to provide an overall view of this problem. Adopting a two-phase coupled model by myosin concentration, we have found new accurate perspectives of acto-cytosolic flow and pressure fields, myosin distribution, as well as the distribution of effective forces across the lamellipodia of a keratocyte with stationary shape. The order of magnitude method is also used to determine the contribution of forces in the internal dynamics of lamellipodia. PMID:26403420
Frictional heating, fluid pressure, and the resistance to fault motion
NASA Astrophysics Data System (ADS)
Lachenbruch, Arthur H.
1980-11-01
Expansion of pore fluid caused by frictional heating might have an important effect on the factional resistance and temperature during an earthquake and a controlling influence on the physics of the earthquake process. When confined water is heated, the pressure increases rapidly (≳10 bars/°C). As Sibson (1973) has pointed out, this could cause a sharp reduction of effective normal stress and dynamic friction on the fault surface. Whether or not this transient stress reduction occurs depends upon the tandem operation of several processes, any of which can break the chain that links frictional heat to frictional stress: the friction must cause an appreciable temperature rise (imposing conditions on the width of the shear zone and rate of conductive transport); the temperature rise must cause an appreciable fluid pressure rise (imposing conditions on the rate of pore dilatation or hydrofracturing, and the rate of Darcian transport); the fluid pressure rise must cause an appreciable reduction of friction (requiring the presence of a continuous fluid phase). Each process depends upon event duration, particle velocity, and the initial value of dynamic friction. With the present uncertainty in the controlling parameters (principally permeability, width of the shear zone, initial stress, and factors controlling transient hydrofracture and pore dilatation) a wide variety of fault behavior is possible. Limits to fault behavior for various ranges of the controlling parameters can be estimated from the governing equations, however, and results can be summarized graphically. If the effective stress law applies and pore dilatation is unimportant, dynamic friction would drop from an initial value of 1 kbar to ˜100 bars when shear strain reached 10 for most earthquakes if the permeability were less than 0.1 μdarcy; the maximum temperature rise would be only ˜150°C irrespective of final strain. If the permeability were ≳100 mdarcies, however, friction would be unaffected
Simultaneous Measurement of Fluid and Particle Motion in Shear Induced Erosion
NASA Astrophysics Data System (ADS)
Krueger, Paul S.; An, Zhongfeng
2015-11-01
Fluid particle interaction is fundamental to shear induced particle erosion, but experimental measurements of this interaction are challenging due to differing optical characteristics of the fluid and particles and because of the high particle volume fraction in the particle bed. To address these challenges, monodisperse glass beads were used with a refractive-index matched aqueous solution of NaI flowing horizontally over the particle bed. Two cameras separately imaged the fluid and particle phase motion using optical filters to isolate the emission bands of the fluorescent fluid tracer particles and dye added to the fluid for the fluid and particle phase cameras, respectively. Then digital particle image velocimetry and particle tracking were used to obtain the full-field, time-varying evolution of the fluid and particle motion simultaneously. The results showed rapid, fluctuating particle transport near flow initiation for sufficiently high fluid flow rates. Increased slip in mean particle velocities was observed above from the particle bed surface and an approximately linear relationship was observed between particle and fluid velocity fluctuations. This material is based on work supported by the National Science Foundation under grant no. 1000908.
On the motion of a heavy rigid body in an ideal fluid with circulation.
Borisov, Alexey V; Mamaev, Ivan S
2006-03-01
We consider Chaplygin's equations [Izd. Akad. Nauk SSSR 3, 3 (1933)] describing the planar motion of a rigid body in an unbounded volume of an ideal fluid while circulation around the body is not zero. Hamiltonian structures and new integrable cases are revealed; certain remarkable partial solutions are found and their stability is examined. The nonintegrability of the system describing the motion of a body in the field of gravity is proved and the chaotic behavior of the system is illustrated. PMID:16599749
Armstrong, William D.; Naughton, Jonathan; Lindberg, William R.
2008-09-02
A shear stress sensor for measuring fluid wall shear stress on a test surface is provided. The wall shear stress sensor is comprised of an active sensing surface and a sensor body. An elastic mechanism mounted between the active sensing surface and the sensor body allows movement between the active sensing surface and the sensor body. A driving mechanism forces the shear stress sensor to oscillate. A measuring mechanism measures displacement of the active sensing surface relative to the sensor body. The sensor may be operated under periodic excitation where changes in the nature of the fluid properties or the fluid flow over the sensor measurably changes the amplitude or phase of the motion of the active sensing surface, or changes the force and power required from a control system in order to maintain constant motion. The device may be operated under non-periodic excitation where changes in the nature of the fluid properties or the fluid flow over the sensor change the transient motion of the active sensor surface or change the force and power required from a control system to maintain a specified transient motion of the active sensor surface.
An analysis of peristaltic motion of compressible convected Maxwell fluid
NASA Astrophysics Data System (ADS)
Abbasi, A.; Ahmad, I.; Ali, N.; Hayat, T.
2016-01-01
This paper presents a theoretical study for peristaltic flow of a non-Newtonian compressible Maxwell fluid through a tube of small radius. Constitutive equation of upper convected Maxwell model is used for the non-Newtonian rheology. The governing equations are modeled for axisymmetric flow. A regular perturbation method is used for the radial and axial velocity components up to second order in dimensionless amplitude. Exact expressions for the first-order radial and axial velocity components are readily obtained while second-order mean axial velocity component is obtained numerically due to presence of complicated non-homogenous term in the corresponding equation. Based on the mean axial velocity component, the net flow rate is calculated through numerical integration. Effects of various emerging parameters on the net flow rate are discussed through graphical illustrations. It is observed that the net flow rate is positive for larger values of dimensionless relaxation time λ1. This result is contrary to that of reported by [D. Tsiklauri and I. Beresnev, "Non-Newtonian effects in the peristaltic flow of a Maxwell fluid," Phys. Rev. E. 64 (2001) 036303]." i.e. in the extreme non-Newtonian regime, there is a possibility of reverse flow.
NASA Astrophysics Data System (ADS)
Arredondo, K.; Billen, M. I.
2015-12-01
Observations of seismicity and seismic tomography provide constraints on the geometry of slabs within mantle, while compression/tension axis derived from moment tensor solutions provide constraints on the internal deformation of slabs. However, since these observations provide only a somewhat blurred or incomplete snapshot of the slab in time, it is difficult to directly relate these observations to the evolution of the slab geometry and the forces acting on and within the slab. In contrast, plate tectonic reconstructions provide time-dependent constraints on the surface motion of plates and the trench at subduction zones, which are related to the dynamical evolution of the slab. We use 2D geodynamical simulations of subduction to explore the relationship between dynamical process within the deforming slab and the observations of surface plate motion and the state-of-stress in slabs. Specifically we utilize models that include the extended Boussinesq approximation (shear heating and latent heat terms in the energy equation), a layered lithosphere with pyrolite, harzburgite and basalt/eclogite, compositionally-dependent phase transitions, and a composite rheology with yielding. The models employ a weak crustal layer that decouples the overriding and subducting plates and allows for dynamically determined trench motion. Here we show that, 1) multiple phase transitions increase slab folding, 2) ridge push significantly increases trench retreat, and 3) strength of the weak crustal layer influences slab detachment. Compared to past studies a more realistic treatment of the phase transitions makes trench retreat more difficult to generate: a weaker plate may encourage slab retreat but detaches once the slab tip crosses into the transition zone due to the rapid increase in slab density. As suggested by previous studies, slab folding within the transition zone changes the direction of forces on the slab and causes periodic changes from trench retreat to trench advance. We
NASA Astrophysics Data System (ADS)
Nissen-Meyer, Tarje; Fournier, Alexandre; Dahlen, F. A.
2008-09-01
We portray a dedicated spectral-element method to solve the elastodynamic wave equation upon spherically symmetric earth models at the expense of a 2-D domain. Using this method, 3-D wavefields of arbitrary resolution may be computed to obtain Fréchet sensitivity kernels, especially for diffracted arrivals. The meshing process is presented for varying frequencies in terms of its efficiency as measured by the total number of elements, their spacing variations and stability criteria. We assess the mesh quantitatively by defining these numerical parameters in a general non-dimensionalized form such that comparisons to other grid-based methods are straightforward. Efficient-mesh generation for the PREM example and a minimum-messaging domain decomposition and parallelization strategy lay foundations for waveforms up to frequencies of 1 Hz on moderate PC clusters. The discretization of fluid, solid and respective boundary regions is similar to previous spectral-element implementations, save for a fluid potential formulation that incorporates the density, thereby yielding identical boundary terms on fluid and solid sides. We compare the second-order Newmark time extrapolation scheme with a newly implemented fourth-order symplectic scheme and argue in favour of the latter in cases of propagation over many wavelengths due to drastic accuracy improvements. Various validation examples such as full moment-tensor seismograms, wavefield snapshots, and energy conservation illustrate the favourable behaviour and potential of the method.
Oscillatory motion of a viscous fluid in a porous medium
Siraev, R. R.
2015-08-15
An oscillatory flow of an incompressible fluid in a saturated porous medium in the presence of a solid inclusion has been theoretically studied. Unsteady filtration has been described by the Brinkman–Forchheimer equation, where inertial effects and terms with acceleration characteristic of high filtration rates and the presence of pulsations are taken into account. The convective part of the acceleration is responsible for nonlinear effects near macroinhomogeneities. These effects can play a noticeable role in unsteady flows in the porous medium, as is shown for the problem of a solid ball streamed by an oscillatory flow having a given velocity at infinity. The results indicate that a secondary averaged flow appears in the case of high frequencies and cannot be described by Darcy’s or Forchheimer’s filtration laws.
Thermal diffusion by Brownian-motion-induced fluid stress
NASA Astrophysics Data System (ADS)
Kreft, Jennifer; Chen, Yeng-Long
2007-08-01
The Ludwig-Soret effect, the migration of a species due to a temperature gradient, has been extensively studied without a complete picture of its cause emerging. Here we investigate the dynamics of DNA and spherical particles subjected to a thermal gradient using a combination of Brownian dynamics and the lattice Boltzmann method. We observe that the DNA molecules will migrate to colder regions of the channel, an observation also made in experiments. In fact, the thermal diffusion coefficient found agrees quantitatively with the experimentally measured value. We also observe that the thermal diffusion coefficient decreases as the radius of the studied spherical particles increases. Furthermore, we observe that the thermal-fluctuation-fluid-momentum-flux coupling induces a gradient in the stress which leads to thermal migration in both systems.
Turbulent fluid motion 2: Scalars, vectors, and tensors
NASA Technical Reports Server (NTRS)
Deissler, Robert G.
1991-01-01
The author shows that the sum or difference of two vectors is a vector. Similarly the sum of any two tensors of the same order is a tensor of that order. No meaning is attached to the sum of tensors of different orders, say u(sub i) + u(sub ij); that is not a tensor. In general, an equation containing tensors has meaning only if all the terms in the equation are tensors of the same order, and if the same unrepeated subscripts appear in all the terms. These facts will be used in obtaining appropriate equations for fluid turbulence. With the foregoing background, the derivation of appropriate continuum equations for turbulence should be straightforward.
Cerebrospinal fluid constituents of cat vary with susceptibility to motion sickness
NASA Technical Reports Server (NTRS)
Lucot, James B.; Crampton, George H.; Matson, Wayne R.; Gamache, Paul H.
1989-01-01
The cerebrospinal fluid drawn from the fourth ventricles of the brains of cats during and after the development of motion sickness was studied to determine what neurotransmitters may be involved in the development of the sickness. The analytical procedure, which uses HPLC coupled with n-electrode coulometric electrochemical detection to measure many compounds with picogram sensitivity, is described. Baseline levels of DOPAC, MHPGSO4, uric acid, DA, 5-HIAA, and HVA were lower on motion and control days in cats which became motion sick when compared with cats which did not. None of the total of 36 identified compounds identified in the samples varied as a function of either exposure to motion or provocation of emesis. It is concluded that susceptibility to motion sickness is a manifestation of individual differences related to fundamental neurochemical composition.
Scaglione, S; Wendt, D; Miggino, S; Papadimitropoulos, A; Fato, M; Quarto, R; Martin, I
2008-08-01
In this study, we investigated the effect of the long-term (10 days) application of a defined and uniform level of fluid flow (uniform shear stress of 1.2 x 10(-3) N/m(2)) on human bone marrow stromal cells (BMSC) cultured on different substrates (i.e., uncoated glass or calcium phosphate coated glass, Osteologictrade mark) in a 2D parallel plate model. Both exposure to flow and culture on Osteologic significantly reduced the number of cell doublings. BMSC cultured under flow were more intensely stained for collagen type I and by von Kossa for mineralized matrix. BMSC exposed to flow displayed an increased osteogenic commitment (i.e., higher mRNA expression of cbfa-1 and osterix), although phenotype changes in response to flow (i.e., mRNA expression of osteopontin, osteocalcin and bone sialoprotein) were dependent on the substrate used. These findings highlight the importance of the combination of physical forces and culture substrate to determine the functional state of differentiating osteoblastic cells. The results obtained using a simple and controlled 2D model system may help to interpret the long-term effects of BMSC culture under perfusion within 3D porous scaffolds, where multiple experimental variables cannot be easily studied independently, and shear stresses cannot be precisely computed. PMID:17969030
Measurements of Coupled Fluid and Sediment Motion Over Mobile Sand Dunes in a Laboratory Flume
Technology Transfer Automated Retrieval System (TEKTRAN)
The relation between turbulent fluid motions and sediment particles over mobile sand dunes may be better understood by examining the time scales over which the quantities fluctuate. In laboratory experiments performed at the USDA-ARS-National Sedimentation Laboratory, profiles of acoustic backscatt...
NASA Astrophysics Data System (ADS)
Li, Zhe; Leduc, Julien; Nunez-Ramirez, Jorge; Combescure, Alain; Marongiu, Jean-Christophe
2015-04-01
We propose a non-intrusive numerical coupling method for transient fluid-structure interaction (FSI) problems simulated by means of different discretization methods: smoothed particle hydrodynamics (SPH) and finite element (FE) methods for the fluid and the solid sub-domains, respectively. As a partitioned coupling method, the present algorithm can ensure a zero interface energy during the whole period of numerical simulation, even in the presence of large interface motion. In other words, the time integrations of the two sub-domains (second order Runge-Kutta scheme for fluid and Newmark integrator for solid) are synchronized. Thanks to this energy-conserving feature, one can preserve the minimal order of accuracy in time and the numerical stability of the FSI simulations, which are validated with a 1D and a 2D trivial numerical test cases. Additionally, some other 2D FSI simulations involving large interface motion have also been carried out with the proposed SPH-FE coupling method. Finally, an example of aquaplaning problem is given in order to show the feasibility of such coupling method in multi-dimensional applications with complicated structural geometries.
Salt-Finger Convection in a Stratified Fluid Layer Induced by Thermal and Solutal Capillary Motion
NASA Technical Reports Server (NTRS)
Chen, Chuan F.; Chan, Cho Lik
1996-01-01
Salt-finger convection in a double-diffusive system is a motion driven by the release of gravitational potential due to different diffusion rates. Normally, when the gravitational field is reduced, salt-finger convection together with other convective motions driven by buoyancy forces will be rapidly suppressed. However, because the destabilizing effect of the concentration gradient is amplified by the Lewis number, with values varying from 10(exp 2) for aqueous salt solutions to 10 (exp 4) for liquid metals, salt-finger convection may be generated at much reduced gravity levels. In the microgravity environment, the surface tension gradient assumes a dominant role in causing fluid motion. In this paper, we report on some experimental results showing the generation of salt-finger convection due to capillary motio on the surface of a stratified fluid layer. A numerical simulation is presented to show the cause of salt-finger convection.
Starting solutions for oscillating motions of an Oldroyd-B fluid over a plane wall
NASA Astrophysics Data System (ADS)
Anjum, Asia; Ayub, Muhammad; Khan, Masood
2012-01-01
In this paper, we establish the starting solutions for oscillating motions of an Oldroyd-B fluid between two side walls perpendicular to a plane wall. The expressions for the velocity field and the associated tangential stress at the bottom wall are obtained, presented under integral and series form. These satisfy all imposed initial and boundary conditions. The obtained solutions are graphically analyzed for the variations of interesting flow parameters. In the absence of side walls, all solutions that have been obtained reduce to those corresponding to the motion over an infinite plate. Moreover, the obtained solutions can be specialized to give similar solutions for Maxwell, second grade and Newtonian fluids performing the same motions.
An incompressible smoothed particle hydrodynamics method for the motion of rigid bodies in fluids
NASA Astrophysics Data System (ADS)
Tofighi, N.; Ozbulut, M.; Rahmat, A.; Feng, J. J.; Yildiz, M.
2015-09-01
A two-dimensional incompressible smoothed particle hydrodynamics scheme is presented for simulation of rigid bodies moving through Newtonian fluids. The scheme relies on combined usage of the rigidity constraints and the viscous penalty method to simulate rigid body motion. Different viscosity ratios and interpolation schemes are tested by simulating a rigid disc descending in quiescent medium. A viscosity ratio of 100 coupled with weighted harmonic averaging scheme has been found to provide satisfactory results. The performance of the resulting scheme is systematically tested for cases with linear motion, rotational motion and their combination. The test cases include sedimentation of a single and a pair of circular discs, sedimentation of an elliptic disc and migration and rotation of a circular disc in linear shear flow. Comparison with previous results at various Reynolds numbers indicates that the proposed method captures the motion of rigid bodies driven by flow or external body forces accurately.
Bacterial migration and motion in a fluid phase and near a solid surface
Frymier, P.D. Jr.
1995-01-01
An understanding of the migration and motion of bacteria in a fluid phase and near solid surfaces is necessary to characterize processes such as the bioremediation of hazardous waste, the pathogenesis of infection, industrial biofouling and wastewater treatment, among others. This study addresses three questions concerning the prediction of the distribution of a population of bacteria in a fluid phase and the motion of bacteria near a solid surface: Under what conditions does a one-dimensional phenomenological model for the density of a population of chemotactic bacteria yield an adequate representation of the migration of bacteria subject to a one-dimensional attractant gradient? How are the values of transport coefficients obtained from experimental data affected by the use of the one-dimensional phenomenological model and also by the use of different descriptions of bacterial swimming behavior in a mathematically rigorous balance equation? How is the characteristic motion of bacteria swimming in a fluid affected by the presence of a solid phase? A computer simulation that rigorously models the movement of a large population of individual chemotactic bacteria in three dimensions is developed to test the validity of a one-dimensional phenomenological model for bacterial migration in a fluid.
NASA Astrophysics Data System (ADS)
Bhatti, M. M.; Zeeshan, A.
2016-06-01
In this paper, effects of variable viscosity with heat transfer on solid particle motion of dusty Jeffrey fluid model through a planar channel has been examined. The governing flow problem for fluid phase and dusty phase is formulated with the help of momentum and energy equation. The resulting coupled ordinary differential equations have been solved analytically and closed form solutions are presented. The influence of all the physical parameters are sketched for velocity profile, pressure rise and temperature profile. Numerical computation is used to evaluate the expression for pressure rise. The present analysis is also presented for Newtonian fluid by taking λ1 → 0 as a special case of our study. It is found that due to the influence of variable viscosity, the fluid velocity changes in the center of the channel and shows opposite behavior near the walls. It is also found that temperature profile increases for larger values of Prandtl number (Pr) and Eckert number (Ec).
Droplet motion in one-component fluids on solid substrates with wettability gradients.
Xu, Xinpeng; Qian, Tiezheng
2012-05-01
Droplet motion on solid substrates has been widely studied not only because of its importance in fundamental research but also because of its promising potentials in droplet-based devices developed for various applications in chemistry, biology, and industry. In this paper, we investigate the motion of an evaporating droplet in one-component fluids on a solid substrate with a wettability gradient. As is well known, there are two major difficulties in the continuum description of fluid flows and heat fluxes near the contact line of droplets on solid substrates, namely, the hydrodynamic (stress) singularity and thermal singularity. To model the droplet motion, we use the dynamic van der Waals theory [Phys. Rev. E 75, 036304 (2007)] for the hydrodynamic equations in the bulk region, supplemented with the boundary conditions at the fluid-solid interface. In this continuum hydrodynamic model, various physical processes involved in the droplet motion can be taken into account simultaneously, e.g., phase transitions (evaporation or condensation), capillary flows, fluid velocity slip, and substrate cooling or heating. Due to the use of the phase field method (diffuse interface method), the hydrodynamic and thermal singularities are resolved automatically. Furthermore, in the dynamic van der Waals theory, the evaporation or condensation rate at the liquid-gas interface is an outcome of the calculation rather than a prerequisite as in most of the other models proposed for evaporating droplets. Numerical results show that the droplet migrates in the direction of increasing wettability on the solid substrates. The migration velocity of the droplet is found to be proportional to the wettability gradients as predicted by Brochard [Langmuir 5, 432 (1989)]. The proportionality coefficient is found to be linearly dependent on the ratio of slip length to initial droplet radius. These results indicate that the steady migration of the droplets results from the balance between the
Korecka, Magdalena; Waligorska, Teresa; Figurski, Michal; Toledo, Jon B; Arnold, Steven E; Grossman, Murray; Trojanowski, John Q; Shaw, Leslie M
2014-01-01
The primary aims of this work were to: 1) establish a calibrator surrogate matrix for quantification of amyloid-β (Aβ)42 in human cerebrospinal fluid (CSF) and preparation of quality control samples for LC-MS-MS methodology, 2) validate analytical performance of the assay, and 3) evaluate its diagnostic utility and compare it with the AlzBio3 immunoassay. The analytical methodology was based on a 2D-UPLC-MS-MS platform. Sample pretreatment used 5 M guanidine hydrochloride and extraction on μElution SPE columns as previously described. A column cleaning procedure involved gradual removal of aqueous solvents by acetonitrile assured consistent long-term chromatography performance. Receiver-operator characteristic (ROC) curve and correlation analyses evaluated the diagnostic utility of UPLC-MS-MS compared to AlzBio3 immunoassay for detection of Alzheimer's disease (AD). The surrogate matrix, artificial CSF containing 4 mg/mL of BSA, provides linear and reproducible calibration comparable to human pooled CSF as calibration matrix. Appropriate cleaning of the trapping and analytical columns provided every-day, trouble-free runs. Analyses of CSF Aβ42 showed that UPLC-MS-MS distinguished neuropathologically-diagnosed AD subjects from healthy controls with at least equivalent diagnostic utility to AlzBio3. Comparison of ROC curves for these two assays showed no statistically significant difference (p = 0.2229). Linear regression analysis of Aβ42 concentrations measured by this mass spectrometry-based method compared to the AlzBio3 immunoassay showed significantly higher but highly correlated results. In conclusion, the newly established surrogate matrix for 2D-UPLC-MS-MS measurement of Aβ42 provides selective, reproducible, and accurate results. The documented analytical performance and diagnostic performance for AD versus controls supports consideration as a candidate reference method. PMID:24625802
Mitri, F G
2015-09-01
The classical Resonance Scattering Theory (RST) for plane waves in acoustics is generalized for the case of a 2D arbitrarily-shaped beam incident upon an elastic cylinder with arbitrary location that is immersed in a nonviscous fluid. The formulation is valid for an elastic (or viscoelastic) cylinder (or a cylindrical shell, a layered cylinder/shell, or a multilayered cylindrical shell, etc.) of any size and material. Partial-wave series expansions (PWSEs) for the incident, internal and scattered fields are derived, and numerical examples illustrate the theory. The wave-fields are expressed using a generalized PWSE involving the beam-shape coefficients (BSCs) and the scattering coefficients of the cylinder. When the beam is shifted off the center of the cylinder, the off-axial BSCs are evaluated by performing standard numerical integration. Acoustic resonance scattering directivity diagrams are calculated by subtracting an appropriate background from the expression of the scattered pressure field. The properties related to the arbitrary scattering of a zeroth-order quasi-Gaussian cylindrical beam (chosen as an example) by an elastic brass cylinder centered on the axis of wave propagation of the beam, and shifted off-axially are analyzed and discussed. Moreover, the total and resonance backscattering form function moduli are numerically computed, and the results discussed with emphasis on the contribution of the surface waves circumnavigating the cylinder circular surface to the resonance backscattering. Furthermore, the analysis is extended to derive general expressions for the axial and transverse acoustic radiation force functions for the cylinder in any 2D beam of arbitrary shape. Examples are provided for a zeroth-order quasi Gaussian cylindrical beam with different waist. Potential applications are in underwater and physical acoustics, however, ongoing research in biomedical ultrasound, non-destructive evaluation, imaging, manufacturing, instrumentation, and
Brownian motion in a rotating fluid: Diffusivity is a function of the rotation rate
NASA Astrophysics Data System (ADS)
Ryskin, Gregory
1988-09-01
The phenomenological relations between thermodynamic fluxes and forces are normally assumed to be invariant with respect to arbitrary motion of the frame of reference. We describe a breakdown of this invariance strong enough to be observable. It is shown that the diffusivity in a rotating fluid is anisotropic and also smaller in magnitude than in a fluid at rest in an inertial frame, giving rise to a diffusion analog of the Hall effect. For large Brownian particles (e.g., biological macromolecules) the diffusivity may decrease by 50% at the rotation speeds achievable in ultracentrifuges.
Simulation of single microorganism motion in fluid based on granular model
NASA Astrophysics Data System (ADS)
Viridi, S.; Nuraini, N.
2016-04-01
Microorganism model for simulating its motion is proposed in this work. It consists of granular particles which can interact to each other through linear spring mimicking microorganism muscles, which is simpler than other model. As a part of the organism organ is moving, while the other remains at its position, it will push the surrounding fluid through Stoke's force and as reaction the fluid pushes back the microorganism. Contracting force is used to change the distance between two points in the organ. Gravity influence is simply neglected in this work. All the considered forces are used to get motion parameters of organism through molecular dynamics method. It is observed that the use of contracting (push-pull) organ constructs slightly more effective model than shrink- and swell-organs as previously investigated, if weighted effectiveness formula is used as function of number of considered forces and involved particles.
On motion of fluid in boundary layer near line of intersection of two planes
NASA Technical Reports Server (NTRS)
Loitsianskii, L G; Bolshakov, V P
1951-01-01
In the paper "The Mutual Interference of Boundary Layers," the authors investigated the problem of the interference of two planes intersecting at right angles on the boundary layers formed by the motion of fluid along the line of intersection of these planes. In the present paper, the results of the preceding one are generalized to the case of planes intersecting at any angle. The motion of a fluid in an angle less than 180 degrees is discussed and the enlargement of the boundary layers near the line of intersection of the planes, the limits of the interference effects of the boundary layers, and the corrections on the drag are determined. All computations are conducted by the Karman-Pohlhausen method for laminar and turbulent boundary layers. The results are reduced to tabulated form.
Xu, Xinpeng; Qian, Tiezheng
2012-06-01
Using a continuum model capable of describing the one-component liquid-gas hydrodynamics down to the contact line scale, we carry out numerical simulation and physical analysis for the droplet motion driven by thermal singularity. For liquid droplets in one-component fluids on heated or cooled substrates, the liquid-gas interface is nearly isothermal. Consequently, a thermal singularity occurs at the contact line and the Marangoni effect due to temperature gradient is suppressed. Through evaporation or condensation in the vicinity of the contact line, the thermal singularity makes the contact angle increase with the increasing substrate temperature. This effect on the contact angle can be used to move the droplets on substrates with thermal gradients. Our numerical results for this kind of droplet motion are explained by a simple fluid dynamical model at the droplet length scale. Since the mechanism for droplet motion is based on the change of contact angle, a separation of length scales is exhibited through a comparison between the droplet motion induced by a wettability gradient and that by a thermal gradient. It is shown that the flow field at the droplet length scale is independent of the statics or dynamics at the contact line scale. PMID:23005105
NASA Astrophysics Data System (ADS)
Baudon, Catherine; Gillet, Hervé; Cremer, Michel
2013-04-01
High-quality bathymetric, 2D seismic and Chirp data located in the southern parts of the Bay of Biscay, France, collected by the University of Bordeaux 1 (Cruises ITSAS 2, 2001; PROSECAN 3, 2006 and SARGASS, 2010) have recently been compiled. The survey area widely covers the Capbreton Canyon, which lies on the boundary between two major structural zones: the Aquitanian passive margin to the North, and the Basque-Cantabrian margin to the South which corresponds to the offshore Pyrenean front. The dataset revealed a large number of key seafloor features potentially associated with focused fluid-flow processes and subsurface sediment-remobilization. Focused fluid migration through sub-seabed sediments is a common phenomenon on continental margins worldwide and has widespread implications from both industrial and fundamental perspectives, from seafloor marine environmental issues to petroleum exploration and hazard assessments. Our study analyses the relationships between seafloor features, deeper structures and fluid migration through the Plio-Quaternary sedimentary pile. The geometrical characteristics, mechanisms of formation and kinematics of four main groups of seabed features have been investigated. (i) A 150km2 field of pockmarks can be observed on the Basque margin. These features are cone-shaped circular or elliptical depressions that are either randomly distributed as small pockmarks (diameter < 20m) or aligned in trains of large pockmarks (ranging from 200 to 600m in diameter) along shallow troughs leading downstream to the Capbreton Canyon. Seismic data show that most pockmarks reach the seabed through vertically staked V-shaped features but some are buried and show evidence of lateral migration through time. (ii) A second field of widely-spaced groups of pockmarks pierce the upper slope of the Aquitanian margin. These depressions are typically a few hundred meters in diameter and seem to be preferentially located in the troughs or on the stoss sides of
Onset and cessation of grain motion in fluid-sheared beds
NASA Astrophysics Data System (ADS)
Clark, Abe; Salevan, Julia; Shattuck, Mark; Ouellette, Nick; O'Hern, Corey
2015-11-01
We performed molecular dynamics simulations of granular beds driven by a model hydrodynamic shear flow to elucidate general grain-scale mechanisms that determine the onset and cessation of sediment transport. By varying the Shields number (the nondimensional shear stress at the top of the bed) and particle Reynolds number (the ratio of particle inertia to viscous damping), we explore how variations of the fluid flow rate, particle inertia, and fluid viscosity affect the onset and cessation of bed motion. For low to moderate particle Reynolds numbers, a critical boundary separates mobile and static states. Transition times between these states diverge as this boundary is approached both from above and below. At high particle Reynolds number, inertial effects become dominant, and particle motion can be sustained well below flow rates at which mobilization of a static bed occurs. We also find that the onset of bed motion (for both low and high particle Reynolds numbers) is described by Weibullian weakest-link statistics, and thus is crucially dependent on the packing structure of the granular bed, even deep beneath the surface. This work was supported by the US Army Research Office under Grant No. W911NF-14-1-0005.
Visualization study of the normal-fluid motion in superfluid helium-4
NASA Astrophysics Data System (ADS)
Gao, Jian; Marokov, Alex; Guo, Wei; van Sciver, Steven W.; Ihas, Gary G.; McKinsey, Daniel N.; Vinen, William F.
2014-03-01
Flow visualization in superfluid 4He is challenging, yet crucial for attaining a detailed understandingof quantum turbulence. Two problems have impeded progress: finding and introducing suitable tracersthat are small yet visible; and unambiguous interpretation of the tracer motion. Metastable He2 triplet molecules form angstrom-sized bubbles in helium and can be imaged using a laser-induced-fluorescence technique. At temperatures above 1 K, helium molecules solely follow the motion of the normal-fluid component without being affected by quantized vortices. In our recent experiments on thermal counterflow, by tracing a thin molecular line created via femtosecond-laser field-ionization technique, we are able to measure the instantaneous normal-fluid velocity field. We show that the obtained velocity probability density function (PDF) obeys a Gaussian distribution. We also discuss the calculated structure function of the novel normal-fluid turbulence in thermal counterflow. The work is supported by the start-up grant of W.G. provided by the National High Magnetic Field Laboratory and Florida State University.
Small Particle Response to Fluid Motion Using Tethered Particles to Simulate Microgravity
NASA Technical Reports Server (NTRS)
Trolinger, James; L'Esperance, Drew; Rangel, Roger; Coimbra, Carlos; Witherow, William K.; Rogers, Jan; Lal, Ravindra
2003-01-01
This paper reports on ground based work conducted to support the Spaceflight Definition project SHIVA (Spaceflight Holography Investigation in a Virtual Apparatus). SHIVA will advance our understanding of the movement of a particle in a fluid. Gravity usually dominates the equations of motion, but in microgravity as well as on earth other terms can become important. Through an innovative application of fractional differential equations, two members of our team produced the first analytical solution of a fundamental equation of motion, which had only been solved numerically or by approximation before. The general solution predicts that the usually neglected history term becomes important in particle response to a sinusoidal fluid movement when the characteristic viscous time is in the same order as the fluid oscillation period and peaks when the two times are equal. In this case three force terms, the Stokes drag, the added mass, and the history drag must all be included in predicting particle movement. We have developed diagnostic recording methods using holography to save all of the particle field data, allowing the experiment to essentially be transferred from space back to earth in what we call the virtual apparatus for on-earth microgravity experimentation. We can quantify precisely the three-dimensional motion of sets of particles, allowing us to test and apply the new analytical solutions. We are examining the response of particles up to 2 mm radius to fluid oscillation at frequencies up to 80 Hz with amplitudes up to 200 microns. Ground studies to support the flight development program have employed various schemes to simulate microgravity. One of the most reliable and meaningful methods uses spheres tethered to a fine hair suspended in the fluid. We have also investigated particles with nearly neutral buoyancy. Recordings are made at the peak amplitudes of vibration of the cell providing a measure of the ratio of fluid to particle amplitude. The experiment
NASA Technical Reports Server (NTRS)
Cho, Y. I.; Crawford, D. W.; Back, L. H.; Back, M. R.
1987-01-01
A flow visualization study using selective dye injection and frame by frame analysis of a movie provided qualitative and quantitative data on the motion of marked fluid particles in a 60 degree artery branch model for simulation of physiological femoral artery flow. Physical flow features observed included jetting of the branch flow into the main lumen during the brief reverse flow period, flow separation along the main lumen wall during the near zero flow phase of diastole when the core flow was in the downstream direction, and inference of flow separation conditions along the wall opposite the branch later in systole at higher branch flow ratios. There were many similarities between dye particle motions in pulsatile flow and the comparative steady flow observations.
On the stability of a convective motion generated by a chemically reacting fluid in a pipe
NASA Astrophysics Data System (ADS)
Koliskina, V.; Kolyshkin, A.; Volodko, I.; Kalis, H.
2016-06-01
Linear stability analysis of a chemically reacting fluid motion in a pipe is performed in the present paper. The reaction rate has an Arrhenius form. The base flow and temperature distribution is obtained from the nonlinear heat equation coupled with the equations of motion. The stability of the flow with respect to asymmetric (spiral) perturbations is investigated numerically. The critical Grasshof number of the flow depends on two dimensionless parameters: the Prandtl number and the Frank-Kamenetsky parameter. The increase of both parameters has a destabilizing influence on the flow. It is shown that the second branch of a marginal stability curve corresponding to smaller critical Grasshof numbers appears as the Prandtl number increases.
Mixed initial-boundary value problem for equations of motion of Kelvin-Voigt fluids
NASA Astrophysics Data System (ADS)
Baranovskii, E. S.
2016-07-01
The initial-boundary value problem for equations of motion of Kelvin-Voigt fluids with mixed boundary conditions is studied. The no-slip condition is used on some portion of the boundary, while the impermeability condition and the tangential component of the surface force field are specified on the rest of the boundary. The global-in-time existence of a weak solution is proved. It is shown that the solution is unique and depends continuously on the field of external forces, the field of surface forces, and initial data.
On the relative rotational motion between rigid fibers and fluid in turbulent channel flow
NASA Astrophysics Data System (ADS)
Marchioli, C.; Zhao, L.; Andersson, H. I.
2016-01-01
In this study, the rotation of small rigid fibers relative to the surrounding fluid in wall-bounded turbulence is examined by means of direct numerical simulations coupled with Lagrangian tracking. Statistics of the relative (fiber-to-fluid) angular velocity, referred to as slip spin in the present study, are evaluated by modelling fibers as prolate spheroidal particles with Stokes number, St, ranging from 1 to 100 and aspect ratio, λ, ranging from 3 to 50. Results are compared one-to-one with those obtained for spherical particles (λ = 1) to highlight effects due to fiber length. The statistical moments of the slip spin show that differences in the rotation rate of fibers and fluid are influenced by inertia, but depend strongly also on fiber length: Departures from the spherical shape, even when small, are associated with an increase of rotational inertia and prevent fibers from passively following the surrounding fluid. An increase of fiber length, in addition, decouples the rotational dynamics of a fiber from its translational dynamics suggesting that the two motions can be modelled independently only for long enough fibers (e.g., for aspect ratios of order ten or higher in the present simulations).
An Unstructured Finite Volume Approach for Structural Dynamics in Response to Fluid Motions
Xia, Guohua; Lin, Ching-Long
2008-01-01
A new cell-vortex unstructured finite volume method for structural dynamics is assessed for simulations of structural dynamics in response to fluid motions. A robust implicit dual-time stepping method is employed to obtain time accurate solutions. The resulting system of algebraic equations is matrix-free and allows solid elements to include structure thickness, inertia, and structural stresses for accurate predictions of structural responses and stress distributions. The method is coupled with a fluid dynamics solver for fluid-structure interaction, providing a viable alternative to the finite element method for structural dynamics calculations. A mesh sensitivity test indicates that the finite volume method is at least of second-order accuracy. The method is validated by the problem of vortex-induced vibration of an elastic plate with different initial conditions and material properties. The results are in good agreement with existing numerical data and analytical solutions. The method is then applied to simulate a channel flow with an elastic wall. The effects of wall inertia and structural stresses on the fluid flow are investigated. PMID:18496602
Massoudi, Mehrdad; Tran, P.X.
2008-09-15
Unsteady problems involving the second grade fluids have received considerable attention in recent years. The present study is an attempt to look at the motion of an oscillating rod in a second grade fluid. Specifically, we solve numerically for the flow of a second grade fluid surrounding a solid cylindrical rod that is suddenly set into longitudinal and torsional motion. The equations are made dimensionless. The results are presented for the shear stresses at the wall, related to the drag force; these are physical quantities of interest, especially in oil-drilling applications.
ERIC Educational Resources Information Center
Feng, Yongqiang; Max, Ludo
2014-01-01
Purpose: Studying normal or disordered motor control requires accurate motion tracking of the effectors (e.g., orofacial structures). The cost of electromagnetic, optoelectronic, and ultrasound systems is prohibitive for many laboratories and limits clinical applications. For external movements (lips, jaw), video-based systems may be a viable…
Motion of Two Compressible Fluids With Interface in a Porous Reservoir
NASA Astrophysics Data System (ADS)
Wooding, R. A.; Weir, G. J.
1984-07-01
This paper describes the development of a physical model of a reservoir involving two compressible fluids, typically gas over water. The reservoir is assumed to be of constant thickness, but its height varies with position, leading to an undulating topography. Withdrawal of the upper fluid begins and proceeds at a constant rate. Topics of interest are the distribution of fluxes throughout the reservoir, the pressure field as a function of time and space, and the movement of the interface between the fluids. The problem is formulated using a version of the two-fluid layer concept. One of the space variables is eliminated by vertical integration, assuming that the fluid velocity vector is almost parallel to the upper and lower confining boundaries (Boussinesq reservoir). The resulting equation of motion for the horizontal components of the flow resembles that for a horizontal reservoir, but terms involving the gradient of the undulating reservoir modify the vertical pressure gradient. The problem is solved in two space dimensions and time. The difficulty presented by intersections of the fluid-fluid interface with one or other of the confining boundaries is resolved by introducing a mathematical continuation of the fluids and their interface through and past the boundary. This is particularly appropriate for "leaky" aquifers; in the case of impermeable boundaries, a concern of the present paper, approximate solutions may be obtained by assuming a virtual porosity and permeability which are extremely small. The method resembles that of Wilson and Sa da Costa (1982). Some useful results have been obtained for Boussinesq reservoirs of fairly arbitrary geometry. Approximate numerical solutions are presented for two-dimensional flow in a rectangular reservoir, the outer part of which is completely submerged, and for axisymmetric and three-dimensional flows in a circular quadrant. A constriction upon the gas flow due to a "dip" in the reservoir is found to produce enhanced
Small Particle Response to Fluid Motion using Tethered Particles to Simulate Microgravity
NASA Technical Reports Server (NTRS)
Trolinger, J. D.; Rangel, R.; Coimbra, C.; Witherow, W.; Rogers, J. R.; Lal, R. B.; Curreri, Peter A. (Technical Monitor)
2002-01-01
This paper reports on ground based work conducted to support the Spaceflight Definition project SHIVA (Spaceflight Holography Investigation in a Virtual Apparatus). SHIVA will advance our understanding of the movement of a particle in a fluid. Gravity usually dominates the equations of motion, but in microgravity as well as on earth other terms can become important. Before two members of our team found an analytical solution of the equations, numerical methods and/or neglecting terms were required. The general solution predicts that the usually neglected history term becomes important when the characteristic viscous time is in the same order as the vibration period and peaks when the two times are equal. In this case three force terms, the Stokes drag, the added mass, and the history drag must all be included in predicting particle movement. We also developed diagnostic recording methods using holography to save all of the particle field data, allowing the experiment to essentially be transferred from space back to earth in what we call the "virtual apparatus". Using state-of-the-art methods in holography we will quantify the three-dimensional motion of sets of particles, allowing us to test and apply the new analytical solutions. The motion of particles up to 4 mm in diameter in a fluid that oscillates at frequencies up to 100 Hz with amplitudes up to 200 microns is being examined. Ground studies to support the flight development program have employed various schemes to simulate microgravity. One of the most reliable and meaningful methods uses spheres tethered to a fine hair suspended in the fluid. We have also investigated particles with nearly neutral buoyancy. Recordings are made at the peak amplitudes of vibration of the cell providing a measure of the ratio of fluid to particle amplitude. The experiment requires precise location of the particle at the time of recording. The hologram of the particle provides microscopic images of the particle that are used
Motion of a Probe Ball in the Fluid under Centrifugal Acceleration
NASA Astrophysics Data System (ADS)
Nyrkova, I. A.; Semenov, A. N.; Khokhlov, A. R.; Linliu, K.; Chu, B.
1997-11-01
The viscosity of a fluid can be measured by observing the motion of a probe sphere (or ball) in a centrifuge tube filled with this fluid. The hydrodynamic behavior of the probe ball moving in the centrifuge tube has been solved theoretically. We have got the universal relationship (for balls of a given material andsize in a given tube) between the terminal ball velocity, the fluid viscosity and the centrifuge acceleration using the only adjustable parameter — the rotational friction coefficient between the ball and the tube. The rotation of the centrifuge tube in the horizontal plane induces an inertia force which is counterbalanced by the friction force acting on the ball. As a result, the ball moves along the tube with some characteristic speed, which is a measure of the viscosity of the fluid. This speed was calculated in the lubrication approximation. The gravitational acceleration causes the ball to move very close to the bottom of the centrifuge tube. In this situation, the gravity is balanced by a “levitation” force introduced and calculated in the present paper. The origin of this force is the formation of the “bubble” behind and below the moving ball. The theoretical development on the terminal velocity for the ball moving very near the bottom of the horizontal centrifuge tube is tested by using a specially designed centrifuge for two types of balls and a wide set of viscosity standards. Excellent agreement between theory and experiment suggests that we have developed a new approach to measure high viscosities of fluids at low shear rates which might be especially useful for the investigation of polymer melts.
Development of a suspension type sliding planar motion table using magnetic fluid lubrication
NASA Astrophysics Data System (ADS)
Li, Xinghui; Shinshi, Tadahiko; Hijikata, Wataru; Morimoto, Yoshihiro
2016-06-01
A sliding planar motion table system that can be used for the lens driving actuator of a laser cutting machine was developed. The system uses magnetic fluid as the lubricant to avoid the leakage of lubricating oil under the table and reduce environmental pollution. The motion table is suspended from the guide surface by an attractive force generated by electromagnets to reduce the contact and frictional forces between the table and the guide surface. The table is capable of movement in one rotational and two translational directions over the guide surface using six electromagnets and three non-contact displacement sensors. Experimental results showed that the magnetic suspension of the table reduced the friction by 82.1% compared to the friction that would otherwise be generated by the dead weight of the table. Circular motion within a diameter of 2 mm was achieved with resolutions of 5 μm and 20 μrad in the translational and rotational directions, respectively. A bandwidth of higher than 100 Hz was also achieved in the three movement directions.
Development of a suspension type sliding planar motion table using magnetic fluid lubrication.
Li, Xinghui; Shinshi, Tadahiko; Hijikata, Wataru; Morimoto, Yoshihiro
2016-06-01
A sliding planar motion table system that can be used for the lens driving actuator of a laser cutting machine was developed. The system uses magnetic fluid as the lubricant to avoid the leakage of lubricating oil under the table and reduce environmental pollution. The motion table is suspended from the guide surface by an attractive force generated by electromagnets to reduce the contact and frictional forces between the table and the guide surface. The table is capable of movement in one rotational and two translational directions over the guide surface using six electromagnets and three non-contact displacement sensors. Experimental results showed that the magnetic suspension of the table reduced the friction by 82.1% compared to the friction that would otherwise be generated by the dead weight of the table. Circular motion within a diameter of 2 mm was achieved with resolutions of 5 μm and 20 μrad in the translational and rotational directions, respectively. A bandwidth of higher than 100 Hz was also achieved in the three movement directions. PMID:27370485
Felderhof, B U
2006-10-14
The motion of a particle immersed in a fluid near a fluid-fluid interface is studied on the basis of the linearized Navier-Stokes equations. The motion is influenced by surface tension, dilatational surface elasticity modulus, and surface shear modulus, as well as by gravity. The backflow at the location of the particle after a sudden impulse has some universal features that are the same as for a rigid wall with stick boundary conditions. At short times the flow depends only on the mass densities of the two fluids. The nature of the short-time flow is calculated from potential flow theory. At a somewhat later time the particle shows a pronounced rebound. The maximum value of the rebound and the time at which the maximum occurs depend on the elastic properties of the interface. PMID:17042642
NASA Technical Reports Server (NTRS)
Martin, J. J.; Holt, J. B.
2000-01-01
This report details the results of a series of fluid motion experiments to investigate the use of magnets to orient fluids in a low-gravity environment. The fluid of interest for this project was liquid oxygen (LO2) since it exhibits a paramagnetic behavior (is attracted to magnetic fields). However, due to safety and handling concerns, a water-based ferromagnetic mixture (produced by Ferrofluidics Corporation) was selected to simplify procedures. Three ferromagnetic fluid mixture strengths and a nonmagnetic water baseline were tested using three different initial fluid positions with respect to the magnet. Experiment accelerometer data were used with a modified computational fluid dynamics code termed CFX-4 (by AEA Technologies) to predict fluid motion. These predictions compared favorably with experiment video data, verifying the code's ability to predict fluid motion with and without magnetic influences. Additional predictions were generated for LO2 with the same test conditions and geometries used in the testing. Test hardware consisted of a cylindrical Plexiglas tank (6-in. bore with 10-in. length), a 6,000-G rare Earth magnet (10-in. ring), three-axis accelerometer package, and a video recorder system. All tests were conducted aboard the NASA Reduced-Gravity Workshop, a KC-135A aircraft.
NASA Astrophysics Data System (ADS)
Polukhin, V. A.; Kurbanova, E. D.
2016-02-01
Molecular dynamics simulation is used to study the thermal stability of the interfacial states of metallic Al, Ag, Sn, Pb, and Hg films (i.e., the structural elements of superconductor composites and conducting electrodes) reinforced by 2D graphene and silicene crystals upon heating up to disordering and to analyze the formation of nonautonomous fluid pseudophases in interfaces. The effect of perforation defects in reinforcing 2D-C and 2D-Si planes with passivated edge covalent bonds on the atomic dynamics is investigated. As compared to Al and Ag, the diffusion coefficients in Pd and Hg films increase monotonically with temperature during thermally activated disordering processes, the interatomic distances decrease, the sizes decrease, drops form, and their density profile grows along the normal. The coagulation of Pb and Hg drops is accompanied by a decrease in the contact angle, the reduction of the interface contact with graphene, and the enhancement of its corrugation (waviness).
Granular spirals on erodible sand bed submitted to a circular fluid motion.
Caps, H; Vandewalle, N
2003-09-01
An experimental study of a granular surface submitted to a circular fluid motion is presented. The appearance of an instability along the sand-water interface is observed beyond a critical radius r(c). This creates ripples with a spiral shape on the granular surface. A phase diagram of such patterns is constructed and discussed as a function of the rotation speed omega of the flow and as a function of the height of water h above the surface. The study of r(c) as a function of h, omega, and r parameters is reported. Thereafter, r(c) is shown to depend on the rotation speed according to a power law. The ripple wavelength is found to decrease when the rotation speed increases and is proportional to the radial distance r. The azimuthal angle epsilon of the spiral arms is studied. It is found that epsilon scales with homegar. This lead to the conclusion that epsilon depends on the fluid momentum. Comparison with experiments performed with fluids allows us to state that the spiral patterns are not the signature of an instability of the boundary layer. PMID:14524759
NASA Astrophysics Data System (ADS)
Takada, Naoki; Misawa, Masaki; Tomiyama, Akio; Fujiwara, Shinya
2000-07-01
This study describes the numerical simulations of two-phase fluid motions under gravity by the lattice Boltzmann method (LBM), in which the fluid motions result from collision and translation of mesoscopic particles and the interface interaction in multiphase fluids can be reproduced in a self-organizing way. Our aims are to examine the applicability of LBM to the numerical analysis of bubble motions in comparison with the two-dimensional results by the Volume Of Fluid (VOF) method based on the Navier-Stokes and the liquid-volume convective equations, and to develop the three-dimensional binary fluids model, consisting of two sets of distribution functions to represent the total fluid density and the density difference, which introduces the repulsive interaction consistent with a free energy function between fluid particles. We included the buoyancy terms due to the density difference between two phases in the lattice Boltzmann equations, and simulated the motions of single bubble and two bubbles rising in a duct, calculating the surface tension from the Laplace's law represented by the non-dimensional numbers, Eotvos and Morton numbers. In the two-dimensional simulations, the results by LBM agree with those by the VOF method. The three-dimensional simulation of two bubble interaction shows that the upper bubble takes a shape of skirt as the lower bubble approaches due to the wake formation, and they coalesce into a single bubble eventually. These results prove the validity of the buoyancy model proposed here and the applicability of LBM to the quantitative numerical analysis of two-phase fluid motions.
Fast algorithms for visualizing fluid motion in steady flow on unstructured grids
NASA Technical Reports Server (NTRS)
Ueng, S. K.; Sikorski, K.; Ma, Kwan-Liu
1995-01-01
The plotting of streamlines is an effective way of visualizing fluid motion in steady flows. Additional information about the flowfield, such as local rotation and expansion, can be shown by drawing in the form of a ribbon or tube. In this paper, we present efficient algorithms for the construction of streamlines, streamribbons and streamtubes on unstructured grids. A specialized version of the Runge-Kutta method has been developed to speed up the integration of particle paths. We have also derived closed-form solutions for calculating angular rotation rate and radius to construct streamribbons and streamtubes, respectively. According to our analysis and test results, these formulations are two to four times better in performance than previous numerical methods. As a large number of traces are calculated, the improved performance could be significant.
Control of self-motion in dynamic fluids: fish do it differently from bees.
Scholtyssek, Christine; Dacke, Marie; Kröger, Ronald; Baird, Emily
2014-05-01
To detect and avoid collisions, animals need to perceive and control the distance and the speed with which they are moving relative to obstacles. This is especially challenging for swimming and flying animals that must control movement in a dynamic fluid without reference from physical contact to the ground. Flying animals primarily rely on optic flow to control flight speed and distance to obstacles. Here, we investigate whether swimming animals use similar strategies for self-motion control to flying animals by directly comparing the trajectories of zebrafish (Danio rerio) and bumblebees (Bombus terrestris) moving through the same experimental tunnel. While moving through the tunnel, black and white patterns produced (i) strong horizontal optic flow cues on both walls, (ii) weak horizontal optic flow cues on both walls and (iii) strong optic flow cues on one wall and weak optic flow cues on the other. We find that the mean speed of zebrafish does not depend on the amount of optic flow perceived from the walls. We further show that zebrafish, unlike bumblebees, move closer to the wall that provides the strongest visual feedback. This unexpected preference for strong optic flow cues may reflect an adaptation for self-motion control in water or in environments where visibility is limited. PMID:24872463
Visiometrics: From Solitons to Vortex Projectiles- Art and Science of Fluid Motions
NASA Astrophysics Data System (ADS)
Zabusky, Norman J.
2003-11-01
This talk is derived from a personal and interdisciplinary half-century journey in the science and art of fluids in motion with emphasis on waves and vortices. The talk is mostly interconnected with images and animations. They are from real astrophysical, geophysical and laboratory images and scientific computer simulations. These images are linked to artistic modes of expression, e.g. painting, photography, sculpture and digital animations, some with musical accompaniment. They are also connected to historical, scientific, mathematical and literary sources. Abstract scientific computing images are created by projecting data to lower dimensions - for visiometric use. The visiometric approach is explained in the context of the discovery of solitons and recent work in vortex dynamics of inhomogeneous media, e.g. the Richtmyer-Meshkov instability [1,2]. [1] Vortex paradigm for accelerated inhomogeneous flows: Visiometrics for the Rayleigh-Taylor and Richtmyer-Meshkov environments. N.J. Zabusky. Ann. Review of Fluid Mechanics, 1999. 31, 495-535. [2] Scientific computing visualization- a new venue in the arts. N.J. Zabusky. Proc. Science and Art Symposium 2000. Eds A. Gyr, P.D. Koumoutsakos & U. Burr. pp 1-11. Kluwer Academic Pub., 2000.
Dynamic simulation of a peristaltic micropump considering coupled fluid flow and structural motion
NASA Astrophysics Data System (ADS)
Lin, Qiao; Yang, Bozhi; Xie, Jun; Tai, Yu-Chong
2007-02-01
This paper presents lumped-parameter simulation of dynamic characteristics of peristaltic micropumps. The pump consists of three pumping cells connected in series, each of which is equipped with a compliant diaphragm that is electrostatically actuated in a peristaltic sequence to mobilize the fluid. Diaphragm motion in each pumping cell is first represented by an effective spring subjected to hydrodynamic and electrostatic forces. These cell representations are then used to construct a system-level model for the entire pump, which accounts for both cell- and pump-level interactions of fluid flow and diaphragm vibration. As the model is based on first principles, it can be evaluated directly from the device's geometry, material properties and operating parameters without using any experimentally identified parameters. Applied to an existing pump, the model correctly predicts trends observed in experiments. The model is then used to perform a systematic analysis of the impact of geometry, materials and pump loading on device performance, demonstrating its utility as an efficient tool for peristaltic micropump design.
NASA Astrophysics Data System (ADS)
Turner, J. Stewart; Gustafson, Lewis B.
1981-12-01
The results of a continuing series of laboratory experiments, designed to model the fluid motions which accompany crystallization, are both described and related in a preliminary way to prototype flows in magma chambers. Previous experiments have demonstrated the importance of compositional inhomogeneity, produced by crystallization and melting in a thermal gradient and coupled with double-diffusive effects, in driving convective flows which result in thermal and compositional stratification in an originally homogeneous fluid. The present experiments examine effects produced in tanks cooled at the side, by the upward flow of a less dense boundary layer depleted in the crystallizing component as crystals grow on the side wall. These processes are examined in simple two and three component aqueous systems (H 2O-Na 2CO 3, H 2O-Na 2CO 3-K 2CO 3, H 2O-CuSO 4-Na 2SO 4) with one and two crystallizing phases. In each of these systems, an initially downward flow of a cooled boundary layer against the side wall is reversed as crystallization commences and depletes the boundary layer in the crystallizing component. Accumulation of this cooler but lighter depleted fluid at the top of the chamber produces thermal and compositional layering by a "filling box" mechanism, partly modified by interchange between the boundary layer and the convecting layers outside. When more than one component is present in the solution, the crystallization process produces a differentiated fluid column, i.e. one with compositional gradients which are different for each of the components. The compositional and thermal distributions within the fluid change with time, but finally appear to reach a steady state. These distributions are the integrated result of compositional changes produced by crystallization from a thin boundary layer, a small proportion of the bulk fluid which evolves in composition and temperature independently of the bulk fluid, in a manner controlled by the dynamics of the system
Uma, B; Radhakrishnan, R; Eckmann, D M; Ayyaswamy, P S
2013-01-01
A hybrid scheme based on Markovian fluctuating hydrodynamics of the fluid and a non-Markovian Langevin dynamics with the Ornstein-Uhlenbeck noise perturbing the translational and rotational equations of motion of a nanoparticle is employed to study the thermal motion of a nearly neutrally buoyant nanoparticle in an incompressible Newtonian fluid medium. A direct numerical simulation adopting an arbitrary Lagrangian-Eulerian based finite element method is employed in simulating the thermal motion of the particle suspended in the fluid contained in a cylindrical vessel. The instantaneous flow around the particle and the particle motion are fully resolved. The numerical results show that (a) the calculated temperature of the nearly neutrally buoyant Brownian particle in a quiescent fluid satisfies the equipartition theorem; (b) the translational and rotational decay of the velocity autocorrelation functions result in algebraic tails, over long time; (c) the translational and rotational mean square displacements of the particle obeys Stokes-Einstein and Stokes-Einstein-Debye relations, respectively; and (d) the parallel and perpendicular diffusivities of the particle closer to the wall are consistent with the analytical results, where available. The study has important implications for designing nanocarriers for targeted drug delivery. PMID:23814315
Liu, Zhongzheng; Kim, Yong-Joe; Wang, Han; Han, Arum
2016-01-01
A numerical modeling method for accurately predicting the acoustophoretic motion of compressible microparticles in microfluidic devices is presented to consider the effects of fluid medium flow and spatial temperature variation that can significantly influence the acoustophoretic motion. In the proposed method, zeroth-order fluid medium flow and temperature, and first- and second-order acoustic fields in the microfluidic devices are first calculated by applying quadratic mapping functions and a second-order finite difference method (FDM) to perturbed mass, momentum, and energy conservation equations and state equation. Then, the acoustic radiation force is obtained based on the Gorkov's acoustic radiation force equation and applied to the Newton's Equation of Motion to calculate the microparticle motion. The proposed method was validated by comparing its results to a commercial software package, COMSOL Multiphysics results, one-dimensional, analytical modeling results, and experimental results. It is shown that the fluid medium flow affects the acoustic radiation force and streaming significantly, resulting in the acoustic radiation force and streaming prediction errors of 10.9% and 67.4%, respectively, when the fluid medium flow speed is increased from 0 to 1 m/s. A local temperature elevation from 20 °C to 22 °C also results in the prediction errors of 88.4% and 73.4%. PMID:26827029
NASA Astrophysics Data System (ADS)
Shin, Dong-Youn; Grassia, Paul; Derby, Brian
2003-09-01
A simple oscillatory, slightly compressible, fluid flow model in a thick-walled piezoelectric tube used in a drop-on-demand inkjet print head is developed from the point of view of fluid-structure interaction to take account of pressure wave propagation and pressure loading opposing wall motion. A frequency sweep is performed computationally using the model revealing the first acoustic fluid-structure resonance frequency and the influence of fluid viscosity. The validity of the model, with given information on the speed of sound in a fluid, is evaluated by comparing the theoretically predicted resonance frequency to the experimentally measured resonance frequency. In addition, the intrinsic speed of sound can be easily computed using the measured acoustic resonance frequency and this computed speed of sound agrees closely with speeds of sound reported in the literature.
NASA Astrophysics Data System (ADS)
Glass, Olivier; Lacave, Christophe; Sueur, Franck
2016-02-01
In this paper we consider the motion of a rigid body immersed in a two dimensional unbounded incompressible perfect fluid with vorticity. We prove that when the body shrinks to a massless pointwise particle with fixed circulation, the "fluid+rigid body" system converges to the vortex-wave system introduced by Marchioro and Pulvirenti (Mathematical theory of incompressible nonviscous fluids. Applied Mathematical Sciences 96, Springer-Verlag, 1994). This extends both the paper (Glass et al. Bull Soc Math France 142(3):489-536, 2014) where the case of a solid tending to a massive pointwise particle was tackled and the paper (Glass et al. Dynamics of a point vortex as limits of a shrinking solid in an irrotational fluid, 2014) where the massless case was considered but in a bounded cavity filled with an irrotational fluid.
Bubble motion through a generalized power-law fluid flowing in a vertical tube.
Mukundakrishnan, Karthik; Eckmann, David M; Ayyaswamy, P S
2009-04-01
Intravascular gas embolism may occur with decompression in space flight, as well as during cardiac and vascular surgery. Intravascular bubbles may be deposited into any end organ, such as the heart or the brain. Surface interactions between the bubble and the endothelial cells lining the vasculature result in serious impairment of blood flow and can lead to heart attack, stroke, or even death. To develop effective therapeutic strategies, there is a need for understanding the dynamics of bubble motion through blood and its interaction with the vessel wall through which it moves. Toward this goal, we numerically investigate the axisymmetric motion of a bubble moving through a vertical circular tube in a shear-thinning generalized power-law fluid, using a front-tracking method. The formulation is characterized by the inlet Reynolds number, capillary number, Weber number, and Froude number. The flow dynamics and the associated wall shear stresses are documented for a combination of two different inlet flow conditions (inlet Reynolds numbers) and three different effective bubble radii (ratio of the undeformed bubble radii to the tube radii). The results of the non-Newtonian model are then compared with that of the model assuming a Newtonian blood viscosity. Specifically, for an almost occluding bubble (effective bubble radius = 0.9), the wall shear stress and the bubble residence time are compared for both Newtonian and non-Newtonian cases. Results show that at low shear rates, for a given pressure gradient the residence time for a non-Newtonian flow is higher than that for a Newtonian flow. PMID:19426324
On the influence of the Basset history force on the motion of a particle through a fluid
NASA Astrophysics Data System (ADS)
Thomas, P. J.
1992-09-01
A parameter study on the development of the Basset integral appearing in the equation of motion of a particle moving through a fluid is presented. The particle motion is investigated numerically for the flow across an aerodynamic shock. It is found that, for this type of flow, the Basset history force acting upon the particle described by the Basset integral can be many times larger than the viscous drag in the immediate shock region. It is demonstrated how this force affects the particle motion across the shock for different particles. Nevertheless, the results obtained show that it is justified to neglect the Basset integral for the theoretical description of the motion of the types of particles commonly used in flow measurement tracer techniques for the type of flow considered here.
Motion of beads in an oscillatory rotating fluid: micro-bead-beating
NASA Astrophysics Data System (ADS)
Nadim, Ali; Sterling, James; Doebler, Robert
2008-11-01
One method for mechanical lysis of biological cells and spores is to mix them with a suspension of beads and vigorously ``shake'' the mixture. The precise mechanisms of lysis are not understood but lysis is thought to result from collisions between the beads and the cells and the associated stresses exerted on the cells. For instance, in the micro-bead-beater^TM instrument from Claremont BioSolutions LLC (Upland, CA), the ``shaking'' occurs when a small cartridge filled with a mixture of cells/spores and 100-micron beads is driven at high frequencies in a small arc trajectory. In this presentation, we describe our initial modeling effort aimed at understanding this system via analysis of the trajectories of beads within such an instrument. The equations governing the motion of non-neutrally-buoyant spherical beads in an oscillatory rotating flow are derived and analyzed numerically. The resulting trajectories are found to be quite complex and very different from those in a steadily rotating fluid. A catalog of possible trajectories at various values of the governing dimensionless parameters is presented.
Effect of fluid motion on the impact erosion by a micro-particle on quartz crystals
NASA Astrophysics Data System (ADS)
Cai, D. H.; Qi, H.; Wen, D. H.; Zhang, L.; Yuan, Q. L.; Chen, Z. Z.
2016-08-01
Abrasive slurry jet (ASJ) is a promising technology to process a variety of materials with advantages of high flexibility, no heat affected zone and high cutting efficiency. In this paper, the impressions generated on a quartz crystal specimen by the impacts of micro-particles laden in a water flow and the associated impact erosion mechanisms are presented and discussed in order to effectively and efficiently control the machining quality. Both brittle and ductile mode erosions coexist in the machining process due to the influence of the fluid motion on the trajectories of particles near the target surface. Large-scale craters produced by brittle conchoidal fractures associated with crashed zone, radial and lateral cracks, dominate the erosion process at large jet impact angles while small-scale craters involving micro-ploughing and micro-cutting are produced by the ductile mode erosion at small jet impact angles. The relation between the process parameters and the overall average volume of craters has also been quantitatively analyzed. A combination of small jet impact angle and abrasive particles and low water pressure is preferred for improving the surface quality after the ASJ machining process caused by the more formation of ductile mode induced craters on the target material, but it is at the sacrifice of the material removal rate as well.
Effect of mass ratio on fluid induced motions of a circular cylinder with strips
NASA Astrophysics Data System (ADS)
Vinod, Ashwin; Banerjee, Arindam
2014-11-01
The objective of the current experimental work is to investigate the effects of mass ratio on Fluid Induced Motions, such as vortex induced vibration (VIV) and galloping, of elastically mounted circular cylinders attached with strips to their outer surface. Although the effect of mass ratio on VIV of a smooth circular cylinder is well documented in literature, however, their effects on circular cylinders with strips, capable of inciting galloping oscillations haven't been investigated and could have potential applications in the domain of vibration based energy harvesters. In the current work, three different mass ratios were tested, out of which, one falls below the critical mass in vortex induced vibration of a circular cylinder. The strips used for the experiments included sandpaper strips of prescribed roughness and smooth strips with no roughness, both of which served as surface protrusion based mechanisms of altering the flow around the cylinder. Interesting variations were observed in the amplitude, frequency response and the power spectrum, depending on the mass ratio of the oscillating system tested. The authors acknowledge support of the Office of Naval Research (Grant # ONR-000141210495 - Dr. Ron Joslin).
Híjar, Humberto
2015-02-01
We study the Brownian motion of a particle bound by a harmonic potential and immersed in a fluid with a uniform shear flow. We describe this problem first in terms of a linear Fokker-Planck equation which is solved to obtain the probability distribution function for finding the particle in a volume element of its associated phase space. We find the explicit form of this distribution in the stationary limit and use this result to show that both the equipartition law and the equation of state of the trapped particle are modified from their equilibrium form by terms increasing as the square of the imposed shear rate. Subsequently, we propose an alternative description of this problem in terms of a generalized Langevin equation that takes into account the effects of hydrodynamic correlations and sound propagation on the dynamics of the trapped particle. We show that these effects produce significant changes, manifested as long-time tails and resonant peaks, in the equilibrium and nonequilibrium correlation functions for the velocity of the Brownian particle. We implement numerical simulations based on molecular dynamics and multiparticle collision dynamics, and observe a very good quantitative agreement between the predictions of the model and the numerical results, thus suggesting that this kind of numerical simulations could be used as complement of current experimental techniques. PMID:25768490
ERIC Educational Resources Information Center
Gerhart, James B.; Nussbaum, Rudi H.
This monograph was written for the Conference on the New Instructional Materials in Physics held at the University of Washington in summer, 1965. It is intended for use in an introductory course in college physics. It consists of an extensive qualitative discussion of motion followed by a detailed development of the quantitative methods needed to…
ERIC Educational Resources Information Center
Brand, Judith, Ed.
2002-01-01
This issue of Exploratorium Magazine focuses on the topic of motion. Contents include: (1) "First Word" (Zach Tobias); (2) "Cosmic Collisions" (Robert Irion); (3) "The Mobile Cell" (Karen E. Kalumuck); (4) "The Paths of Paths" (Steven Vogel); (5) "Fragments" (Pearl Tesler); (6) "Moving Pictures" (Amy Snyder); (7) "Plants on the Go" (Katharine…
Fluid Motion and the Toroidal Magnetic Field Near the Top of Earth's Liquid Outer Core.
NASA Astrophysics Data System (ADS)
Celaya, Michael Augustine
This work considers two unresolved problems central to the study of Earth's deep interior: (1) What is the surface flow of the complete three dimensional motion sustaining the geomagnetic field in the fluid outer core? (2) How strong is the toroidal component of that field just beneath the mantle inside the core?. A solution of these problems is necessary to achieve even a basic understanding of magnetic field generation, and core-mantle interactions. Progress in solving (1) is made by extending previous attempts to resolve the core surface flow, and identifying obstacles which lead to distorted solutions. The extension relaxes the steady motions constraint. This permits more realistic solutions which should resemble more closely the real Earth flow. A difficulty with the assumption of steady flow is that if the real motion is unsteady, as it is likely to be, then steady models will suffer from aliasing. Aliased solutions can be highly corrupted. The effects of aliasing incurred through model underparametrization are explored. It is found that flow spectral energy must fall rapidly with increasing degree to escape aliasing's distortion. Damping does not appear to remedy the problem, but in fact obscures it by forcing the solution to converge upon a single, but possibly still aliased estimate. Inversions of a magnetic field model for unsteady motions, indicate steady flows are indeed aliased in time. By comparison, unsteady flows appear free of aliasing and show significant temporal variation, changing by about 30% of their magnitude over 20 years. However, it appears that noise in the high degree secular variation (SV) data used to determine the flow acts as a further impediment to solving (1). Damping is shown to be effective in removing noise, but only once aliasing is no longer a factor and noise is restricted to that part of the SV which makes only a small contribution to the solution. To solve (2) the radial component of Ohm's law is inverted for the toroidal
NASA Astrophysics Data System (ADS)
Mahmood, A.; Fetecau, C.; Khan, N. A.; Jamil, M.
2010-08-01
The velocity field and the associated shear stress corresponding to the longitudinal oscillatory flow of a generalized second grade fluid, between two infinite coaxial circular cylinders, are determined by means of the Laplace and Hankel transforms. Initially, the fluid and cylinders are at rest and at t = 0+ both cylinders suddenly begin to oscillate along their common axis with simple harmonic motions having angular frequencies Ω1 and Ω2. The solutions that have been obtained are presented under integral and series forms in terms of the generalized G and R functions and satisfy the governing differential equation and all imposed initial and boundary conditions. The respective solutions for the motion between the cylinders, when one of them is at rest, can be obtained from our general solutions. Furthermore, the corresponding solutions for the similar flow of ordinary second grade fluid and Newtonian fluid are also obtained as limiting cases of our general solutions. At the end, the effect of different parameters on the flow of ordinary second grade and generalized second grade fluid are investigated graphically by plotting velocity profiles.
2-D Animation's Not Just for Mickey Mouse.
ERIC Educational Resources Information Center
Weinman, Lynda
1995-01-01
Discusses characteristics of two-dimensional (2-D) animation; highlights include character animation, painting issues, and motion graphics. Sidebars present Silicon Graphics animations tools and 2-D animation programs for the desktop computer. (DGM)
Modeling on Fluid Flow and Inclusion Motion in Centrifugal Continuous Casting Strands
NASA Astrophysics Data System (ADS)
Wang, Qiangqiang; Zhang, Lifeng; Sridhar, Seetharaman
2016-05-01
During the centrifugal continuous casting process, unreasonable casting parameters can cause violent level fluctuation, serious gas entrainment, and formation of frozen shell pieces at the meniscus. Thus, in the current study, a three-dimensional multiphase turbulent model was established to study the transport phenomena during centrifugal continuous casting process. The effects of nozzle position, casting and rotational speed on the flow pattern, centrifugal force acting on the molten steel, level fluctuation, gas entrainment, shear stress on mold wall, and motion of inclusions during centrifugal continuous casting process were investigated. Volume of Fluid model was used to simulate the molten steel-air two-phase. The level fluctuation and the gas entrainment during casting were calculated by user-developed subroutines. The trajectory of inclusions in the rotating system was calculated using the Lagrangian approach. The results show that during centrifugal continuous casting, a large amount of gas was entrained into the molten steel, and broken into bubbles of various sizes. The greater the distance to the mold wall, the smaller the centrifugal force. Rotation speed had the most important influence on the centrifugal force distribution at the side region. Angular moving angle of the nozzle with 8° and keeping the rotation speed with 60 revolutions per minute can somehow stabilize the level fluctuation. The increase of angular angle of nozzle from 8 to 18 deg and rotation speed from 40 to 80 revolutions per minute favored to decrease the total volume of entrained bubbles, while the increase of distance of nozzle moving left and casting speed had reverse effects. The trajectories of inclusions in the mold were irregular, and then rotated along the strand length. After penetrating a certain distance, the inclusions gradually moved to the center of billet and gathered there. More work, such as the heat transfer, the solidification, and the inclusions entrapment
Modeling on Fluid Flow and Inclusion Motion in Centrifugal Continuous Casting Strands
NASA Astrophysics Data System (ADS)
Wang, Qiangqiang; Zhang, Lifeng; Sridhar, Seetharaman
2016-08-01
During the centrifugal continuous casting process, unreasonable casting parameters can cause violent level fluctuation, serious gas entrainment, and formation of frozen shell pieces at the meniscus. Thus, in the current study, a three-dimensional multiphase turbulent model was established to study the transport phenomena during centrifugal continuous casting process. The effects of nozzle position, casting and rotational speed on the flow pattern, centrifugal force acting on the molten steel, level fluctuation, gas entrainment, shear stress on mold wall, and motion of inclusions during centrifugal continuous casting process were investigated. Volume of Fluid model was used to simulate the molten steel-air two-phase. The level fluctuation and the gas entrainment during casting were calculated by user-developed subroutines. The trajectory of inclusions in the rotating system was calculated using the Lagrangian approach. The results show that during centrifugal continuous casting, a large amount of gas was entrained into the molten steel, and broken into bubbles of various sizes. The greater the distance to the mold wall, the smaller the centrifugal force. Rotation speed had the most important influence on the centrifugal force distribution at the side region. Angular moving angle of the nozzle with 8° and keeping the rotation speed with 60 revolutions per minute can somehow stabilize the level fluctuation. The increase of angular angle of nozzle from 8 to 18 deg and rotation speed from 40 to 80 revolutions per minute favored to decrease the total volume of entrained bubbles, while the increase of distance of nozzle moving left and casting speed had reverse effects. The trajectories of inclusions in the mold were irregular, and then rotated along the strand length. After penetrating a certain distance, the inclusions gradually moved to the center of billet and gathered there. More work, such as the heat transfer, the solidification, and the inclusions entrapment
NASA Astrophysics Data System (ADS)
Leonov, G. A.; Kuznetsov, N. V.; Mokaev, T. N.
2015-11-01
In this paper a Lorenz-like system, describing convective fluid motion in rotating cavity, is considered. It is shown numerically that this system, like the classical Lorenz system, possesses a homoclinic trajectory and a chaotic self-excited attractor. However, for the considered system, unlike the classical Lorenz system, along with self-excited attractor a hidden attractor can be localized. Analytical-numerical localization of hidden attractor is demonstrated.
NASA Astrophysics Data System (ADS)
Bartzke, Gerhard; Rogers, Benedict D.; Fourtakas, Georgios; Mokos, Athanasios; Huhn, Katrin
2016-04-01
The processes that cause the creation of a variety of sediment morphological features, e.g. laminated beds, ripples, or dunes, are based on the initial motion of individual sediment grains. However, with experimental techniques it is difficult to measure the flow characteristics, i.e., the velocity of the pore water flow in sediments, at a sufficient resolution and in a non-intrusive way. As a result, the role of fluid infiltration at the surface and in the interior affecting the initiation of motion of a sediment bed is not yet fully understood. Consequently, there is a strong need for numerical models, since these are capable of quantifying fluid driven sediment transport processes of complex sediment beds composed of irregular shapes. The numerical method Smoothed Particle Hydrodynamics (SPH) satisfies this need. As a meshless and Lagrangian technique, SPH is ideally suited to simulating flows in sediment beds composed of various grain shapes, but also flow around single grains at a high temporal and spatial resolution. The solver chosen is DualSPHysics (www.dual.sphysics.org) since this is validated for a range of flow conditions. For the present investigation a 3-D numerical flume model was generated using SPH with a length of 4.0 cm, a width of 0.05 cm and a height of 0.2 cm where mobile sediment particles were deposited in a recess. An experimental setup was designed to test sediment configurations composed of irregular grain shapes (grain diameter, D50=1000 μm). Each bed consisted of 3500 mobile objects. After the bed generation process, the entire domain was flooded with 18 million fluid particles. To drive the flow, an oscillating motion perpendicular to the bed was applied to the fluid, reaching a peak value of 0.3 cm/s, simulating 4 seconds of real time. The model results showed that flow speeds decreased logarithmically from the top of the domain towards the surface of the beds, indicating a fully developed boundary layer. Analysis of the fluid
Wada, Yuji; Kundu, Tribikram; Nakamura, Kentaro
2014-08-01
The distributed point source method (DPSM) is extended to model wave propagation in viscous fluids. Appropriate estimation on attenuation and boundary layer formation due to fluid viscosity is necessary for the ultrasonic devices used for acoustic streaming or ultrasonic levitation. The equations for DPSM modeling in viscous fluids are derived in this paper by decomposing the linearized viscous fluid equations into two components-dilatational and rotational components. By considering complex P- and S-wave numbers, the acoustic fields in viscous fluids can be calculated following similar calculation steps that are used for wave propagation modeling in solids. From the calculations reported the precision of DPSM is found comparable to that of the finite element method (FEM) for a fundamental ultrasonic field problem. The particle velocity parallel to the two bounding surfaces of the viscous fluid layer between two rigid plates (one in motion and one stationary) is calculated. The finite element results agree well with the DPSM results that were generated faster than the transient FEM results. PMID:25096081
Propagator-resolved 2D exchange in porous media in the inhomogeneous magnetic field.
Burcaw, Lauren M; Hunter, Mark W; Callaghan, Paul T
2010-08-01
We present a propagator-resolved 2D exchange spectroscopy technique for observing fluid motion in a porous medium. The susceptibility difference between the matrix and the fluid is exploited to produce an inhomogeneous internal magnetic field, causing the Larmor frequency to change as molecules migrate. We test our method using a randomly packed monodisperse 100 microm diameter glass bead matrix saturated with distilled water. Building upon previous 2D exchange spectroscopy work we add a displacement dimension which allows us to obtain 2D exchange spectra that are defined by both mixing time and spatial displacement rather than by mixing time alone. We also simulate our system using a Monte Carlo process in a random nonpenetrating monodisperse bead pack, finding good agreement with experiment. A simple analytic model is used to interpret the NMR data in terms of a characteristic length scale over which molecules must diffuse to sample the inhomogeneous field distribution. PMID:20554230
Inertial solvation in femtosecond 2D spectra
NASA Astrophysics Data System (ADS)
Hybl, John; Albrecht Ferro, Allison; Farrow, Darcie; Jonas, David
2001-03-01
We have used 2D Fourier transform spectroscopy to investigate polar solvation. 2D spectroscopy can reveal molecular lineshapes beneath ensemble averaged spectra and freeze molecular motions to give an undistorted picture of the microscopic dynamics of polar solvation. The transition from "inhomogeneous" to "homogeneous" 2D spectra is governed by both vibrational relaxation and solvent motion. Therefore, the time dependence of the 2D spectrum directly reflects the total response of the solvent-solute system. IR144, a cyanine dye with a dipole moment change upon electronic excitation, was used to probe inertial solvation in methanol and propylene carbonate. Since the static Stokes' shift of IR144 in each of these solvents is similar, differences in the 2D spectra result from solvation dynamics. Initial results indicate that the larger propylene carbonate responds more slowly than methanol, but appear to be inconsistent with rotational estimates of the inertial response. To disentangle intra-molecular vibrations from solvent motion, the 2D spectra of IR144 will be compared to the time-dependent 2D spectra of the structurally related nonpolar cyanine dye HDITCP.
Enhancement of biomixing by swimming cells in 2D films
NASA Astrophysics Data System (ADS)
Gollub, Jerry; Kurtuldu, Huseyin; Guasto, Jeffrey; Johnson, Karl
2011-11-01
Fluid mixing in active suspensions of microorganisms is important to ecological phenomena and shows surprising statistical behavior. We investigate the mixing produced by swimming unicellular algal cells (Chlamydomonas) in quasi-2D films by tracking the motions of cells and of microscopic passive tracer particles advected by the fluid. The reduced spatial dimension of the system leads to long-range flows and a surprisingly strong dependence of tracer transport on the swimmer concentration. The mean square displacements are well described by a stochastic Langevin model, with an effective diffusion coefficient D growing as the 3/2 power of the swimmer concentration, due to the interaction of tracer particles with multiple swimmers. We also discuss the anomalous probability distributions of tracer displacements, which become Gaussian at high concentration, but show strong power-law tails at low concentration. Supported by NSF Grant DMR-0803153.
Kent, J C; Eaton, A R
1982-03-01
A new technique has been developed for studies of fluid motion within the cylinder of a reciprocating piston engine during the air induction process. Helium-filled bubbles, serving as neutrally buoyant flow tracer particles, enter the cylinder along with the inducted air charge. The bubble motion is recorded by stereo cine photography through the transparent cylinder of a specially designed research engine. Quantitative data on the 3-D velocity field generated during induction is obtained from frame-to-frame analysis of the stereo images, taking into account refraction of the rays due to the transparent cylinder. Other applications for which this technique appears suitable include measurements of velocity fields within intake ports and flow-field dynamics within intake manifolds of multicylinder engines. PMID:20372559
Friedel, Michael J.
2001-01-01
This report describes a model for simulating transient, Variably Saturated, coupled water-heatsolute Transport in heterogeneous, anisotropic, 2-Dimensional, ground-water systems with variable fluid density (VST2D). VST2D was developed to help understand the effects of natural and anthropogenic factors on quantity and quality of variably saturated ground-water systems. The model solves simultaneously for one or more dependent variables (pressure, temperature, and concentration) at nodes in a horizontal or vertical mesh using a quasi-linearized general minimum residual method. This approach enhances computational speed beyond the speed of a sequential approach. Heterogeneous and anisotropic conditions are implemented locally using individual element property descriptions. This implementation allows local principal directions to differ among elements and from the global solution domain coordinates. Boundary conditions can include time-varying pressure head (or moisture content), heat, and/or concentration; fluxes distributed along domain boundaries and/or at internal node points; and/or convective moisture, heat, and solute fluxes along the domain boundaries; and/or unit hydraulic gradient along domain boundaries. Other model features include temperature and concentration dependent density (liquid and vapor) and viscosity, sorption and/or decay of a solute, and capability to determine moisture content beyond residual to zero. These features are described in the documentation together with development of the governing equations, application of the finite-element formulation (using the Galerkin approach), solution procedure, mass and energy balance considerations, input requirements, and output options. The VST2D model was verified, and results included solutions for problems of water transport under isohaline and isothermal conditions, heat transport under isobaric and isohaline conditions, solute transport under isobaric and isothermal conditions, and coupled water
Energy Science and Technology Software Center (ESTSC)
2005-07-01
Aniso2d is a two-dimensional seismic forward modeling code. The earth is parameterized by an X-Z plane in which the seismic properties Can have monoclinic with x-z plane symmetry. The program uses a user define time-domain wavelet to produce synthetic seismograms anrwhere within the two-dimensional media.
Preece, D.S. Perkins, E.D.
1999-02-10
Techniques for modeling oil well sand production have been developed using the formulations for superquadric discrete elements and Darcy fluid flow. Discrete element models are generated using the new technique of particle cloning. Discrete element sources and sinks allow simulation of sand production from the initial state through the transition to an equilibrium state where particles are created and removed at the same rate.
ERIC Educational Resources Information Center
Vollmer, Michael; Mollmann, Klaus-Peter
2012-01-01
We present fascinating simple demonstration experiments recorded with high-speed cameras in the field of fluid dynamics. Examples include oscillations of falling droplets, effects happening upon impact of a liquid droplet into a liquid, the disintegration of extremely large droplets in free fall and the consequences of incompressibility. (Contains…
NASA Technical Reports Server (NTRS)
Nastula, J.; Ponte, R. M.; Salstein, D. A.
2007-01-01
Three sets of degree-2, order-1 harmonics of the gravity field, derived from the Gravity Recovery and Climate Experiment (GRACE) data processed at the Center for Space Research (CSR), Jet Propulsion Laboratory (JPL) and GeoforschungsZentrum (GFZ), are used to compute polar motion excitation functions X1 and X2. The GFZ and JPL excitations and the CSR X2, excitation compare generally well with geodetically observed excitation after removal of effects of oceanic currents and atmospheric winds. The agreement considerably exceeds that from previous GRACE data releases. For the JPL series, levels of correlation with the geodetic observations and the variance explained are comparable to, but still lower than, those obtained independently from available models and analyses of the atmosphere, ocean, and land hydrology. Improvements in data quality of gravity missions are still needed to deliver even tighter constraints on mass-related excitation of polar motion.
2015-01-01
We report a simple yet highly efficient chemical motor that can be controlled with visible light. The motor made from a noble metal and doped silicon acts as a pump, which is driven through a light-activated catalytic reaction process. We show that the actuation is based on electro-osmosis with the electric field generated by chemical reactions at the metal and silicon surfaces, whereas the contribution of diffusio-osmosis to the actuation is negligible. Surprisingly, the pump can be operated using water as fuel. This is possible because of the large ζ-potential of silicon, which makes the electro-osmotic fluid motion sizable even though the electric field generated by the reaction is weak. The electro-hydrodynamic process is greatly amplified with the addition of reactive species, such as hydrogen peroxide, which generates higher electric fields. Another remarkable finding is the tunability of silicon-based pumps. That is, it is possible to control the speed of the fluid with light. We take advantage of this property to manipulate the spatial distribution of colloidal microparticles in the liquid and to pattern colloidal microparticle structures at specific locations on a wafer surface. Silicon-based pumps hold great promise for controlled mass transport in fluids. PMID:26349036
Esplandiu, Maria J; Farniya, Ali Afshar; Bachtold, Adrian
2015-11-24
We report a simple yet highly efficient chemical motor that can be controlled with visible light. The motor made from a noble metal and doped silicon acts as a pump, which is driven through a light-activated catalytic reaction process. We show that the actuation is based on electro-osmosis with the electric field generated by chemical reactions at the metal and silicon surfaces, whereas the contribution of diffusio-osmosis to the actuation is negligible. Surprisingly, the pump can be operated using water as fuel. This is possible because of the large ζ-potential of silicon, which makes the electro-osmotic fluid motion sizable even though the electric field generated by the reaction is weak. The electro-hydrodynamic process is greatly amplified with the addition of reactive species, such as hydrogen peroxide, which generates higher electric fields. Another remarkable finding is the tunability of silicon-based pumps. That is, it is possible to control the speed of the fluid with light. We take advantage of this property to manipulate the spatial distribution of colloidal microparticles in the liquid and to pattern colloidal microparticle structures at specific locations on a wafer surface. Silicon-based pumps hold great promise for controlled mass transport in fluids. PMID:26349036
NASA Technical Reports Server (NTRS)
Deissler, Robert G.
1996-01-01
Background material on Fourier analysis and on the spectral form of the continuum equations, both averaged and unaveraged, are given. The equations are applied to a number of cases of homogeneous turbulence with and without mean gradients. Spectral transfer of turbulent activity between scales of motion is studied in some detail. The effects of mean shear, heat transfer, normal strain, and buoyancy are included in the analyses.
Comparisons between measurement and analysis of fluid motion in internal combustion engines
Witze, P.O.
1981-10-01
The Engine Combustion Technology Project was created for the purpose of promoting the development of advanced piston engine concepts by the development of techniques to measure, analyze, and understand the combustion process. The technologies emphasized in the project include laser-based measurement techniques and large-scale computer simulations. Considerable progress has already been achieved by project participants in modeling engine air motion, fuel sprays, and engine combustion phenomena. This milestone report covers one part of that progress, summarizing the current capabilities of multi-dimensional computer codes being developed by the project to predict the behavior of turbulent air motion in an engine environment. Computed results are compared directly with experimental data in six different areas of importance to internal combustion engines: (1) Induction-generated ring-vortex structures; (2) Piston-induced vortex roll-up; (3) Behavior of turbulence during compression; (4) Decay of swirling flow during compression; (5) Decay of swirling flow in a constant volume engine simulator; (6) Exhaust-pipe flow. The computational procedures used include vortex dynamics, rapid distortion theory, and finite difference models employing two-equation and subgrid-scale turbulence models. Although the capability does not yet exist to predict the air motion in an engine from its geometric configuration alone, the results presented show that many flowfield sub-processes can be predicted given well-specified initial and boundary conditions.
Greg Flach, Frank Smith
2011-12-31
Mesh2d is a Fortran90 program designed to generate two-dimensional structured grids of the form [x(i),y(i,j)] where [x,y] are grid coordinates identified by indices (i,j). The x(i) coordinates alone can be used to specify a one-dimensional grid. Because the x-coordinates vary only with the i index, a two-dimensional grid is composed in part of straight vertical lines. However, the nominally horizontal y(i,j0) coordinates along index i are permitted to undulate or otherwise vary. Mesh2d also assigns an integer material type to each grid cell, mtyp(i,j), in a user-specified manner. The complete grid is specified through three separate input files defining the x(i), y(i,j), and mtyp(i,j) variations.
Energy Science and Technology Software Center (ESTSC)
2011-12-31
Mesh2d is a Fortran90 program designed to generate two-dimensional structured grids of the form [x(i),y(i,j)] where [x,y] are grid coordinates identified by indices (i,j). The x(i) coordinates alone can be used to specify a one-dimensional grid. Because the x-coordinates vary only with the i index, a two-dimensional grid is composed in part of straight vertical lines. However, the nominally horizontal y(i,j0) coordinates along index i are permitted to undulate or otherwise vary. Mesh2d also assignsmore » an integer material type to each grid cell, mtyp(i,j), in a user-specified manner. The complete grid is specified through three separate input files defining the x(i), y(i,j), and mtyp(i,j) variations.« less
The dispersion of the jet fluid due to the large scale motion in bluffbody flows
NASA Astrophysics Data System (ADS)
Ghoniem, Ahmed F.; Martins, Luis-Filipe; Kelly, John; Rotman, Douglas
1991-06-01
The unsteady flow established behind the bluff body, simulated numerically using the vortex method, is found to depend strongly on the velocity ratio between the two streams. At low jet velocity, the recirculation zone is long and quiet and the jet flow is stagnated close to the bluffbody. As the jet velocity approaches the annular flow velocity, the recirculation zone becomes unsteady with strong oscillations associated with the shedding of large eddies from both sides of the bluffbody, and the jet stagnation point approaches that of the annular flow. At even higher jet velocities, the jet penetrates through the recirculation zone and the unsteadiness is weakened. Simulations of the dispersion of the inner jet fluid into the recirculation zone show that these large scale phenomena play an important role in the mixing between the two streams. At low jet velocities, the jet fluid is dispersed almost uniformly within the recirculation zone of the bluffbody, while at high jet velocities, the jet fluid remains confined within a narrow zone around the centerline of the bluffbody. For intermediate values of the jet velocity, mixed zones appear intermittently in the form of large toroidal eddies which are shed from the downstream end of the recirculation zone.
Electroosmotic fluid motion and late-time solute transport at non-negligible zeta potentials
S. K. Griffiths; R. H. Nilson
1999-12-01
Analytical and numerical methods are employed to determine the electric potential, fluid velocity and late-time solute distribution for electroosmotic flow in a tube and channel when the zeta potential is not small. The electric potential and fluid velocity are in general obtained by numerical means. In addition, new analytical solutions are presented for the velocity in a tube and channel in the extremes of large and small Debye layer thickness. The electroosmotic fluid velocity is used to analyze late-time transport of a neutral non-reacting solute. Zeroth and first-order solutions describing axial variation of the solute concentration are determined analytically. The resulting expressions contain eigenvalues representing the dispersion and skewness of the axial concentration profiles. These eigenvalues and the functions describing transverse variation of the concentration field are determined numerically using a shooting technique. Results are presented for both tube and channel geometries over a wide range of the normalized Debye layer thickness and zeta potential. Simple analytical approximations to the eigenvalues are also provided for the limiting cases of large and small values of the Debye layer thickness. The methodology developed here for electroosmotic flow is also applied to the Taylor problem of late-time transport and dispersion in pressure-driven flows.
NASA Astrophysics Data System (ADS)
Lotsch, Bettina V.
2015-07-01
Graphene's legacy has become an integral part of today's condensed matter science and has equipped a whole generation of scientists with an armory of concepts and techniques that open up new perspectives for the postgraphene area. In particular, the judicious combination of 2D building blocks into vertical heterostructures has recently been identified as a promising route to rationally engineer complex multilayer systems and artificial solids with intriguing properties. The present review highlights recent developments in the rapidly emerging field of 2D nanoarchitectonics from a materials chemistry perspective, with a focus on the types of heterostructures available, their assembly strategies, and their emerging properties. This overview is intended to bridge the gap between two major—yet largely disjunct—developments in 2D heterostructures, which are firmly rooted in solid-state chemistry or physics. Although the underlying types of heterostructures differ with respect to their dimensions, layer alignment, and interfacial quality, there is common ground, and future synergies between the various assembly strategies are to be expected.
Fluid motion within the cylinder of internal combustion engines - The 1986 Freeman Scholar Lecture
NASA Astrophysics Data System (ADS)
Heywood, John B.
1987-03-01
Aspects of gas motion into, within, and out of the engine cylinder which govern the combustion characteristics and capabilities of spark-ignition engines and compression-ignition or diesel engines are considered. Flow characteristics through inlet and exhaust valves in four-stroke cycle engines, and through ports in the cylinder liner in two-stroke cycle engines, are discussed. Features and turbulence characteristics of common in-cylinder flows including the large scale rotating flows precipitated by the conical intake jet and two-stroke scavenger flows are reviewed. The flow phenomenon near walls are then discussed, with application to heat transfer and hydrocarbon emissions phenomena.
Fluid motion within the cylinder of internal combustion engines - The 1986 Freeman Scholar Lecture
Heywood, J.B.
1987-03-01
The flow field within the cylinder of internal combustion engines is the most important factor controlling the combustion process. Thus it has a major impact on engine operation. This paper reviews those aspects of gas motion into, within, and out of the engine cylinder that govern the combustion characteristics and breathing capabilities of spark-ignition engines and compression-ignition or diesel engines. Necessary background information and reciprocating engine operating cycles, the primary effect of piston motion and the spark-ignition and diesel engine combustion processes is first summarized. Then the characteristics of flow through inlet and exhaust valves in four-stroke cycle engines, and through ports in the cylinder liner in two-stroke cycle engines are reviewed. The essential features of common in-cylinder flows - the large scale rotating flows set up by the conical intake jet, the creation and development of swirl about the cylinder axis, the flows produced during compression due to combustion chamber shape called squish, flow during the combustion process, and two-stroke scavenging flows - are then described. The turbulence characteristics of these flows are then defined and discussed. Finally, flow phenomena which occur near the walls, which are important to heat transfer and hydrocarbon emissions phenomena, are reviewed.
Uma, B.; Radhakrishnan, R.; Eckmann, D.M.
2014-01-01
A direct numerical simulation adopting an arbitrary Lagrangian-Eulerian based finite element method is employed to simulate the motion of a nanocarrier in a quiescent fluid contained in a cylindrical tube. The nanocarrier is treated as a solid sphere. Thermal fluctuations are implemented using two different approaches: (1) fluctuating hydrodynamics; (2) generalized Langevin dynamics (Mittag-Leffler noise). At thermal equilibrium, the numerical predictions for temperature of the nanoparticle, velocity distribution of the particle, decay of the velocity autocorrelation function, diffusivity of the particle and particle-wall interactions are evaluated and compared with analytical results, where available. For a neutrally buoyant nanoparticle of 200 nm radius, the comparisons between the results obtained from the fluctuating hydrodynamics and the generalized Langevin dynamics approaches are provided. Results for particle diffusivity predicted by the fluctuating hydrodynamics approach compare very well with analytical predictions. Ease of computation of the thermostat is obtained with the Langevin approach although the dynamics gets altered. PMID:25621317
Jenny, Patrick Torrilhon, Manuel; Heinz, Stefan
2010-02-20
In this paper, a stochastic model is presented to simulate the flow of gases, which are not in thermodynamic equilibrium, like in rarefied or micro situations. For the interaction of a particle with others, statistical moments of the local ensemble have to be evaluated, but unlike in molecular dynamics simulations or DSMC, no collisions between computational particles are considered. In addition, a novel integration technique allows for time steps independent of the stochastic time scale. The stochastic model represents a Fokker-Planck equation in the kinetic description, which can be viewed as an approximation to the Boltzmann equation. This allows for a rigorous investigation of the relation between the new model and classical fluid and kinetic equations. The fluid dynamic equations of Navier-Stokes and Fourier are fully recovered for small relaxation times, while for larger values the new model extents into the kinetic regime. Numerical studies demonstrate that the stochastic model is consistent with Navier-Stokes in that limit, but also that the results become significantly different, if the conditions for equilibrium are invalid. The application to the Knudsen paradox demonstrates the correctness and relevance of this development, and comparisons with existing kinetic equations and standard solution algorithms reveal its advantages. Moreover, results of a test case with geometrically complex boundaries are presented.
On the collision of impulsive gravitational waves when coupled with fluid motions
NASA Astrophysics Data System (ADS)
Chandrasekhar, S.; Xanthopoulos, B. C.
1985-11-01
An exact solution of Einstein's equations, with a source derived from a perfect fluid in which the energy density, epsilon, is equal to the pressure, p, is obtained. The solution describes the space-time following the collision of plane impulsive gravitational waves and is the natural generalization of the Nutku-Hali solution of the vacuum equations, in the region of interaction under similar basic conditions. A consistent extension of the solution, prior to the instant of collision, requires that the fluid in the region of interaction is the direct result of a transformation of incident null-dust (i.e., of massless particles describing null trajectories). The ultimate result of the collision is the development of a space-time singularity, the nature of which is strongly dependent on the amplitude and the character of the sound waves that are present. The distribution of epsilon that follows the collision has many intriguing features. The solution obtained in this paper provides the first example of an induced transformation of a massless into a massive particle.
An experimental and numerical study of wave motion and upstream influence in a stratified fluid
NASA Technical Reports Server (NTRS)
Hurdis, D. A.
1974-01-01
A system consisting of two superimposed layers of liquid of different densities, with a thin transition layer at the interface, provides a good laboratory model of an ocean thermocline or of an atmospheric inversion layer. This research was to gain knowledge about the propagation of disturbances within these two geophysical systems. The technique used was to observe the propagation of internal waves and of upstream influence within the density-gradient region between the two layers of liquid. The disturbances created by the motion of a vertical flat plate, which was moved longitudinally through this region, were examined both experimentally and numerically. An upstream influence, which resulted from a balance of inertial and gravitational forces, was observed, and it was possible to predict the behavior of this influence with the numerical model. The prediction included a description of the propagation of the upstream influence to steadily increasing distances from the flat plate and the shapes and magnitudes of the velocity profiles.
Turbulent fluid motion IV-averages, Reynolds decomposition, and the closure problem
NASA Technical Reports Server (NTRS)
Deissler, Robert G.
1992-01-01
Ensemble, time, and space averages as applied to turbulent quantities are discussed, and pertinent properties of the averages are obtained. Those properties, together with Reynolds decomposition, are used to derive the averaged equations of motion and the one- and two-point moment or correlation equations. The terms in the various equations are interpreted. The closure problem of the averaged equations is discussed, and possible closure schemes are considered. Those schemes usually require an input of supplemental information unless the averaged equations are closed by calculating their terms by a numerical solution of the original unaveraged equations. The law of the wall for velocities and temperatures, the velocity- and temperature-defect laws, and the logarithmic laws for velocities and temperatures are derived. Various notions of randomness and their relation to turbulence are considered in light of ergodic theory.
von Kármán-Howarth and Corrsin equations closure based on Lagrangian description of the fluid motion
NASA Astrophysics Data System (ADS)
de Divitiis, Nicola
2016-05-01
A new approach to obtain the closure formulas for the von Kármán-Howarth and Corrsin equations is presented, which is based on the Lagrangian representation of the fluid motion, and on the Liouville theorem associated to the kinematics of a pair of fluid particles. This kinematics is characterized by the finite scale separation vector which is assumed to be statistically independent from the velocity field. Such assumption is justified by the hypothesis of fully developed turbulence and by the property that this vector varies much more rapidly than the velocity field. This formulation leads to the closure formulas of von Kármán-Howarth and Corrsin equations in terms of longitudinal velocity and temperature correlations following a demonstration completely different with respect to the previous works. Some of the properties and the limitations of the closed equations are discussed. In particular, we show that the times of evolution of the developed kinetic energy and temperature spectra are finite quantities which depend on the initial conditions.
Peristaltic Motion of Johnson-Segalman Fluid in a Curved Channel with Slip Conditions
Hina, Sadia; Mustafa, Meraj; Hayat, Tasawar
2014-01-01
Slip effects on the peristaltic transport of Johnson-Segalman fluid through a curved channel have been addressed. The influence of wall properties is also analyzed. Long wavelength and low Reynolds number assumptions have been utilized in the mathematical formulation of the problem. The equations so formed have been solved numerically by shooting method through computational software Mathematica 8. In addition the analytic solution for small Weissenberg number (elastic parameter) is computed through a regular perturbation method. An excellent agreement is noticed between the two solutions. The results indicate an increase in the magnitude of velocity with an intensification in the slip effect. Moreover the size and circulation of the trapped boluses increase with an increase in the slip parameter. Unlike the planar channel, the profiles of axial velocity are not symmetric about the central line of the channel. PMID:25474212
Peristaltic motion of Johnson-Segalman fluid in a curved channel with slip conditions.
Hina, Sadia; Mustafa, Meraj; Hayat, Tasawar
2014-01-01
Slip effects on the peristaltic transport of Johnson-Segalman fluid through a curved channel have been addressed. The influence of wall properties is also analyzed. Long wavelength and low Reynolds number assumptions have been utilized in the mathematical formulation of the problem. The equations so formed have been solved numerically by shooting method through computational software Mathematica 8. In addition the analytic solution for small Weissenberg number (elastic parameter) is computed through a regular perturbation method. An excellent agreement is noticed between the two solutions. The results indicate an increase in the magnitude of velocity with an intensification in the slip effect. Moreover the size and circulation of the trapped boluses increase with an increase in the slip parameter. Unlike the planar channel, the profiles of axial velocity are not symmetric about the central line of the channel. PMID:25474212
Uma, B.; Eckmann, D.M.; Ayyaswamy, P.S.; Radhakrishnan, R.
2012-01-01
A novel hybrid scheme based on Markovian fluctuating hydrodynamics of the fluid and a non-Markovian Langevin dynamics with the Ornstein-Uhlenbeck noise perturbing the translational and rotational equations of motion of the nanoparticle is employed to study the thermal motion of a nanoparticle in an incompressible Newtonian fluid medium. A direct numerical simulation adopting an arbitrary Lagrangian-Eulerian (ALE) based finite element method (FEM) is employed in simulating the thermal motion of a particle suspended in the fluid confined in a cylindrical vessel. The results for thermal equilibrium between the particle and the fluid are validated by comparing the numerically predicted temperature of the nanoparticle with that obtained from the equipartition theorem. The nature of the hydrodynamic interactions is verified by comparing the velocity autocorrelation function (VACF) and mean squared displacement (MSD) with well-known analytical results. For nanoparticle motion in an incompressible fluid, the fluctuating hydrodynamics approach resolves the hydrodynamics correctly but does not impose the correct equipartition of energy based on the nanoparticle mass because of the added mass of the displaced fluid. In contrast, the Langevin approach with an appropriate memory is able to show the correct equipartition of energy, but not the correct short- and long-time hydrodynamic correlations. Using our hybrid approach presented here, we show for the first time, that we can simultaneously satisfy the equipartition theorem and the (short- and long-time) hydrodynamic correlations. In effect, this results in a thermostat that also simultaneously preserves the true hydrodynamic correlations. The significance of this result is that our new algorithm provides a robust computational approach to explore nanoparticle motion in arbitrary geometries and flow fields, while simultaneously enabling us to study carrier adhesion mediated by biological reactions (receptor
Fluids in micropores. V. Effects of thermal motion in the walls of a slit-micropore
Diestler, D.J.; Schoen, M.
1996-05-01
Previous articles in this series have concerned the prototypal slit-pore with {ital rigid} walls, in which a Lennard-Jones (12,6) monatomic film is constrained between two plane-parallel walls comprising like atoms fixed in the face-centered-cubic (fcc) (100) configuration. The behavior of molecularly thin films in the rigid-wall prototype is governed by the template effect, whereby solid films can form epitaxially when the walls are properly aligned in the lateral directions. In this article the influence of thermal motion of the wall atoms on the template effect is investigated. The walls are treated as Einstein solids, the atoms moving independently in harmonic potentials centered on rigidly fixed equilibrium positions in the fcc (100) configuration. The force constant {ital f}{sub {ital c}} is a measure of the stiffness of the walls, the rigid-wall limit being {ital f}{sub {ital c}}={infinity}. Formal thermodynamic and statistical mechanical analyses of the system are carried out. The results of grand canonical ensemble Monte Carlo simulations indicate that for values of {ital f}{sub {ital c}} characteristic of a soft (e.g., noble-gas) crystal dynamic coupling between wall and film has a substantial influence on such equilibrium properties as normal stress (load) and interfacial tensions. In general, the softer the walls (i.e., the smaller the value of {ital f}{sub {ital c}}), the weaker the template effect and hence the softer and more disordered the confined film. {copyright} {ital 1996 American Institute of Physics.}
NASA Astrophysics Data System (ADS)
Leary, K. P.; Schmeeckle, M. W.
2014-12-01
Despite numerous experimental and numerical studies investigating transport over ripples and dunes in rivers, the spatiotemporal details of the pattern of transport over bedforms remain largely unknown. Here we report turbulence-resolving, simultaneous measurements of bedload motion and near-bed fluid velocity downstream of a backward facing step in a laboratory flume. Details are compared to a coupled large eddy simulation and distinct element simulation (LES-DEM) of the same geometry. Two synchronized high-speed video cameras simultaneously observed bed load motion and the motion of neutrally buoyant particles in a laser light sheet 6 mm above the bed at 250 frames/s downstream of a 3.8 cm backward-facing step. Particle imaging velocimetry algorithms were applied to the laser sheet images to obtain two-dimensional field of two-dimensional vectors while manual particle tracking techniques were applied to the video images of the bed. As expected, there is a strong positive correlation between sediment flux and near-bed fluid velocity. Sediment flux was determined by manually tracking grains that passed over a 6 cm long line in the middle of the field of view on the bedload images. Sediment flux increased monotonically downstream of flow reattachment. Localized, intermittent, high-magnitude transport events were more apparent near flow reattachment than further downstream. Often, these high-magnitude events were seen to have significant cross-stream particle velocities. These events are consistent with permeable "splat events" visualized in the LES-DEM numerical simulations, wherein a volume of fluid moves toward and impinges on the bed. Fluid impingement and penetration of the bed results in outward flow and sediment motion from the center of the splat. Work is ongoing to quantify spatial and temporal autocorrelations and covariances of the fluid velocity and sediment motions.
NASA Astrophysics Data System (ADS)
Elkins, Christopher J.; Alley, Marcus T.
2007-12-01
Magnetic resonance velocimetry (MRV) is a non-invasive technique capable of measuring the three-component mean velocity field in complex three-dimensional geometries with either steady or periodic boundary conditions. The technique is based on the phenomenon of nuclear magnetic resonance (NMR) and works in conventional magnetic resonance imaging (MRI) magnets used for clinical imaging. Velocities can be measured along single lines, in planes, or in full 3D volumes with sub-millimeter resolution. No optical access or flow markers are required so measurements can be obtained in clear or opaque MR compatible flow models and fluids. Because of its versatility and the widespread availability of MRI scanners, MRV is seeing increasing application in both biological and engineering flows. MRV measurements typically image the hydrogen protons in liquid flows due to the relatively high intrinsic signal-to-noise ratio (SNR). Nonetheless, lower SNR applications such as fluorine gas flows are beginning to appear in the literature. MRV can be used in laminar and turbulent flows, single and multiphase flows, and even non-isothermal flows. In addition to measuring mean velocity, MRI techniques can measure turbulent velocities, diffusion coefficients and tensors, and temperature. This review surveys recent developments in MRI measurement techniques primarily in turbulent liquid and gas flows. A general description of MRV provides background for a discussion of its accuracy and limitations. Techniques for decreasing scan time such as parallel imaging and partial k-space sampling are discussed. MRV applications are reviewed in the areas of physiology, biology, and engineering. Included are measurements of arterial blood flow and gas flow in human lungs. Featured engineering applications include the scanning of turbulent flows in complex geometries for CFD validation, the rapid iterative design of complex internal flow passages, velocity and phase composition measurements in
Effect of particle geometry on triple line motion of nano-fluid drops and deposit nano-structuring.
Askounis, Alexandros; Sefiane, Khellil; Koutsos, Vasileios; Shanahan, Martin E R
2015-08-01
We illustrate the importance of particle geometry on droplet contact line pinning, 'coffee-stain' formation and nano-structuring within the resulting rings. We present the fundamentals of pure liquid droplet evaporation and then discuss the effect of particles on the evaporation process. The resulting coffee-stain patterns and particle structuring within them are presented and discussed. In the second part, we turn our attention to the effect of particle geometry on the evaporation process. A wide range of particle shapes, categorised according to aspect ratio, from the simple shape of a sphere to the highly irregular shapes of platelets and tubes is discussed. Particle geometry effect on evaporation behaviour was quantified in terms of change in contact angle and contact radius for the stick-slip cases. Consequently the hysteretic energy barrier pinning the droplets was estimated, showing an increasing trend with particle aspect ratio. The three-phase contact line (TL) motion kinetics are complemented with analysis of the nano-structuring behaviour of each shape, leading to the identification of the two main parameters affecting nanoparticle self-assembly behaviour at the wedge. Flow velocity and wedge constraints were found to have antagonist effects on particle deposition, although these varied with particle shape. This description should help in understanding the drying behaviour of more complex fluids. Furthermore, knowing the fundamentals of this simple and inexpensive surface patterning technique should permit its tailoring to the needs of many potential applications. PMID:24927853
NASA Astrophysics Data System (ADS)
Mettier, Ralph; Pfiffner, O. Adrian
2010-05-01
Surface motion is, apart from the obvious topography, the most easily accessible and best quantifiable characteristic of a typical alpine-style orogen. While it is understood that several different processes, such as i.e. isostatic unloading and thermodynamic effects contribute to the overall motion, it is mostly unclear how large the individual contributions are, and how much of the observed motion is a consequence of ongoing tectonic shortening. A number of methods, such as enhanced GPS measurements, Interferometric Synthetic Aperture Radar (InSAR) and fission track (FT) dating, as well as precise leveling can now provide us with a good description of the vertical motion at present as well as in the fairly recent history of the orogen. This in turn, provides us with reliable, and often much needed, criteria for calibrating conceptual and numerical models of orogenesis and the involved processes. We present a series of finite element models, that attempt to reproduce the observed vertical surface motion on a roughly north-south cross section of the Swiss Alps in the 'ABAQUS' commercial FEM package. Unlike most comparable modeling approaches, we apply a fairly simple formulation of rheology, and focus on a highly complex geometrical representation of the cross section, constructed of individual tectonomorphic units such as the Aar- and Gotthard massifs, the Helvetic and Penninic nappe structures as well as the underlying subduction of the European crust. The models simulate a short timespan, with a fixed rate of shortening prescribed by the boundary conditions and the various interactions between the tectonomorphic units being the dominant adjustable parameters. The resulting motion at the surface of the model, as well as the internal deformation of the individual tectonomorphic units is then examined, interpreted and compared to their real-world counterparts. The models incorporate variations in the chosen physical descriptions of the materials, deforming in
NASA Astrophysics Data System (ADS)
Shi, Juan; Qiu, Bing; Tan, Hui-Li
2009-06-01
A lattice Boltzmann model is presented to simulate the deformation and motions of a red blood cell (RBC) in a shear flow. The curvatures of the membrane of a static RBC with different chemical potential drops calculated by our model agree with those computed by a shooting method very well. Our simulation results show that in a shear flow, a biconcave RBC becomes highly flattened and undergoes tank-treading motion. With intrinsically parallel dynamics, this lattice Boltzmann method is expected to find wide applications to both single and multi-vesicles suspension as well as complex open membranes in various fluid flows for a wide range of Reynolds numbers.
NASA Astrophysics Data System (ADS)
Najafi, Amin
2014-05-01
Using the Monte Carlo simulations, we have calculated mean-square fluctuations in statistical mechanics, such as those for colloids energy configuration are set on square 2D periodic substrates interacting via a long range screened Coulomb potential on any specific and fixed substrate. Random fluctuations with small deviations from the state of thermodynamic equilibrium arise from the granular structure of them and appear as thermal diffusion with Gaussian distribution structure as well. The variations are showing linear form of the Fluctuation-Dissipation Theorem on the energy of particles constitutive a canonical ensemble with continuous diffusion process of colloidal particle systems. The noise-like variation of the energy per particle and the order parameter versus the Brownian displacement of sum of large number of random steps of particles at low temperatures phase are presenting a markovian process on colloidal particles configuration, too.
Kim, Sungheon; Decarlo, Lindsey; Cho, Gene Y.; Jensen, Jens H.; Sodickson, Daniel K.; Moy, Linda; Formenti, Silvia; Schneider, Robert J.; Goldberg, Judith D.; Sigmund, Eric E.
2013-01-01
Effective delivery of therapeutic drug to the core of a tumor is often impeded by physiological barriers, such as interstitial fluid pressure (IFP). There are a number of therapies to lower IFP and induce tumor vascular normalization. However, lack of a non-invasive means to measure IFP hinders utilization of such a window of opportunity for maximizing the treatment response. Thus, the purpose of this study was to investigate the feasibility of using intravoxel incoherent motion (IVIM) diffusion parameters as noninvasive imaging biomarkers for IFP. Mice bearing the 4T1 mammary carcinoma model were studied with diffusion weighted magnetic resonance imaging (DWI) immediately followed by wick-in-needle IFP measurement. Voxelwise analysis was conducted with a conventional monoexponential diffusion model as well as a biexponential model taking IVIM into account. There was no significant correlation of IFP with either median apparent diffusion coefficient from the monoexponential model (r = 0.11, p = 0.78) or median tissue diffusivity from the biexponential model (r = 0.30, p = 0.44). On the other hand, IFP was correlated with the median pseudo-diffusivity (Dp) of apparent vascular voxels (r = 0.76, p = 0.02) and with the median product of perfusion-fraction and pseudo-diffusivity (fp·Dp) of apparent vascular voxels (r = 0.77, p = 0.02). Although the effect of IVIM in tumors has been reported previously, to our knowledge, this study represents the first direct comparison of IVIM metrics with IFP, with the results supporting the feasibility of using IVIM-DWI metrics as noninvasive biomarkers for tumor IFP. PMID:22072561
NASA Astrophysics Data System (ADS)
Wang, Jin; Ma, Jianyong; Zhou, Changhe
2014-11-01
A 3×3 high divergent 2D-grating with period of 3.842μm at wavelength of 850nm under normal incidence is designed and fabricated in this paper. This high divergent 2D-grating is designed by the vector theory. The Rigorous Coupled Wave Analysis (RCWA) in association with the simulated annealing (SA) is adopted to calculate and optimize this 2D-grating.The properties of this grating are also investigated by the RCWA. The diffraction angles are more than 10 degrees in the whole wavelength band, which are bigger than the traditional 2D-grating. In addition, the small period of grating increases the difficulties of fabrication. So we fabricate the 2D-gratings by direct laser writing (DLW) instead of traditional manufacturing method. Then the method of ICP etching is used to obtain the high divergent 2D-grating.
Koch, Jon; Borg, John; Mattson, Abby; Olsen, Kris; Bahcall, James
2012-01-01
Objective. This in vitro study compared the flow pattern and shear stress of an irrigant induced by ultrasonic and polymer rotary finishing file activation in an acrylic root canal model. Flow visualization analysis was performed using an acrylic canal filled with a mixture of distilled water and rheoscopic fluid. The ultrasonic and polymer rotary finishing file were separately tested in the canal and activated in a static position and in a cyclical axial motion (up and down). Particle movement in the fluid was captured using a high-speed digital camera and DaVis 7.1 software. The fluid shear stress analysis was performed using hot film anemometry. A hot-wire was placed in an acrylic root canal and the canal was filled with distilled water. The ultrasonic and polymer rotary finishing files were separately tested in a static position and in a cyclical axial motion. Positive needle irrigation was also tested separately for fluid shear stress. The induced wall shear stress was measured using LabVIEW 8.0 software. PMID:22461994
NASA Astrophysics Data System (ADS)
Celigueta, Miguel Angel; Deshpande, Kedar M.; Latorre, Salvador; Oñate, Eugenio
2016-04-01
We present a procedure for coupling the finite element method (FEM) and the discrete element method (DEM) for analysis of the motion of particles in non-Newtonian fluids. Particles are assumed to be spherical and immersed in the fluid mesh. A new method for computing the drag force on the particles in a non-Newtonian fluid is presented. A drag force correction for non-spherical particles is proposed. The FEM-DEM coupling procedure is explained for Eulerian and Lagrangian flows, and the basic expressions of the discretized solution algorithm are given. The usefulness of the FEM-DEM technique is demonstrated in its application to the transport of drill cuttings in wellbores.
Rise characteristics of gas bubbles in a 2D rectangular column: VOF simulations vs experiments
Krishna, R.; Baten, J.M. van
1999-10-01
About five centuries ago, Leonardo da Vinci described the sinuous motion of gas bubbles rising in water. The authors have attempted to simulate the rise trajectories of bubbles of 4, 5, 7, 8, 9, 12, and 20 mm in diameter rising in a 2D rectangular column filled with water. The simulations were carried out using the volume-of-fluid (VOF) technique developed by Hirt and Nichols (J. Computational Physics, 39, 201--225 (1981)). To solve the Navier-Stokes equations of motion the authors used a commercial solver, CFX 4.1c of AEA Technology, UK. They developed their own bubble-tracking algorithm to capture sinuous bubble motions. The 4 and 5 mm bubbles show large lateral motions observed by Da Vinci. The 7, 8 and 9 mm bubble behave like jellyfish. The 12 mm bubble flaps its wings like a bird. The extent of lateral motion of the bubbles decreases with increasing bubble size. Bubbles larger than 20 mm in size assume a spherical cap form and simulations of the rise characteristics match experiments exactly. VOF simulations are powerful tools for a priori determination of the morphology and rise characteristics of bubbles rising in a liquid. Bubble-bubble interactions are also properly modeled by the VOF technique.
Gross, M.B.
1984-10-01
STEALTH is a family of computer codes that can be used to calculate a variety of physical processes in which the dynamic behavior of a continuum is involved. The version of STEALTH described in this volume is designed for calculations of fluid-structure interaction. This version of the program consists of a hydrodynamic version of STEALTH which has been coupled to a finite-element code, WHAMSE. STEALTH computes the transient response of the fluid continuum, while WHAMSE computes the transient response of shell and beam structures under external fluid loadings. The coupling between STEALTH and WHAMSE is performed during each cycle or step of a calculation. Separate calculations of fluid response and structural response are avoided, thereby giving a more accurate model of the dynamic coupling between fluid and structure. This volume provides the theoretical background, the finite-difference equations, the finite-element equations, a discussion of several sample problems, a listing of the input decks for the sample problems, a programmer's manual and a description of the input records for the STEALTH/WHAMSE computer program.
Efficient framework for deformable 2D-3D registration
NASA Astrophysics Data System (ADS)
Fluck, Oliver; Aharon, Shmuel; Khamene, Ali
2008-03-01
Using 2D-3D registration it is possible to extract the body transformation between the coordinate systems of X-ray and volumetric CT images. Our initial motivation is the improvement of accuracy of external beam radiation therapy, an effective method for treating cancer, where CT data play a central role in radiation treatment planning. Rigid body transformation is used to compute the correct patient setup. The drawback of such approaches is that the rigidity assumption on the imaged object is not valid for most of the patient cases, mainly due to respiratory motion. In the present work, we address this limitation by proposing a flexible framework for deformable 2D-3D registration consisting of a learning phase incorporating 4D CT data sets and hardware accelerated free form DRR generation, 2D motion computation, and 2D-3D back projection.
NASA Technical Reports Server (NTRS)
Hide, Raymond
1995-01-01
General expressions (with potential applications in several areas of geophysical fluid dynamics) are derived for all three components of the contribution made by the geostrophic part of the pressure field associated with flow in a rotating gravitating fluid to the topographic torque exerted by the fluid on a rigid impermeable bounding surface of any shape. When applied to the Earth's liquid metallic core, which is bounded by nearly spherical surfaces and can be divided into two main regions, the "torosphere" and "polosphere," the expressions reduce to formulae given previously by the author, thereby providing further support for his work and that of others on the role of topographic coupling at the core-mantle boundary in the excitation by core motions of Earth rotation fluctuations on decadal time scales. They also show that recent criticisms of that work are vitiated by mathematical and physical errors. Contrary to these criticisms, the author's scheme for exploiting Earth rotation and other geophysical data (either real or simulated in computer models) in quantitative studies of the topography of the core-mantle boundary (CMB) by intercomparing various models of (a) motions in the core based on geomagnetic secular variation data and (b) CMB topography based on seismological and gravity data has a sound theoretical basis. The practical scope of the scheme is of course limited by the accuracy of real data, but this is a matter for investigation, not a priori assessment.
Energy Science and Technology Software Center (ESTSC)
2004-08-01
AnisWave2D is a 2D finite-difference code for a simulating seismic wave propagation in fully anisotropic materials. The code is implemented to run in parallel over multiple processors and is fully portable. A mesh refinement algorithm has been utilized to allow the grid-spacing to be tailored to the velocity model, avoiding the over-sampling of high-velocity materials that usually occurs in fixed-grid schemes.
Characterization of Porous Medium Properties Using 2D NMR
NASA Astrophysics Data System (ADS)
Sun, Boqin; Dunn, Keh-Jim
2003-03-01
We have successfully applied the concept of 2D NMR to the characterization of properties of fluid-saturated porous medium. Using a two-windowed modified CPMG pulse sequence, we were able to explore the magnetic internal filed gradient distribution within the pore space of a fluid-saturated porous medium due to magnetic susceptibility contrast between the solid matrix and pore fluid. Similar scheme is used to identify and quantify different types of pore fluids, such as oil, water, and gas, based on the contrast in their diffusion coefficients. The magic angle spinning technique (MAS) can also be applied in the 2D NMR framework for delineating the chemical shift spectra of the pore fluids in a porous medium at different T1 or T2 relaxation times. The results can be displayed in a two-dimensional plot, with one axis being the T1 or T2 relaxation times, the other axis being the internal field gradient, diffusion coefficient, or chemical shift, and the third axis being the proton population. Our preliminary laboratory work indicates that the 2D NMR approach can be a powerful tool for the characterization of properties of fluid-saturated porous medium, such as fluid typing, oil viscosity determination, surface wettability, etc.
NASA Astrophysics Data System (ADS)
Mayor, Louise
2016-05-01
Graphene might be the most famous example, but there are other 2D materials and compounds too. Louise Mayor explains how these atomically thin sheets can be layered together to create flexible “van der Waals heterostructures”, which could lead to a range of novel applications.
NASA Astrophysics Data System (ADS)
Rionero, Salvatore
2013-05-01
Either for its great geophysical relevance or the frequent occurrence of porous materials in real life, research on convective-diffusive fluid motions in porous horizontal layers has a notable relevance, which is increasing with the number of salts dissolved in the fluid. In the present paper, porous horizontal layers heated from below and salted by m salts partly from above and partly from below are studied forall min {N}. In the Darcy-Boussinesq scheme it is shown that: (i) the L2 solutions are bounded, uniquely determined, and asymptotically converging toward an absorbing set; (ii) for each Fourier component of the perturbations to the thermal conduction solution, there exists an own nonlinear admissible evolution system; (iii) subcritical instabilities do not exist and the conditions of linear stability also guarantee the global nonlinear stability; (iv) global nonlinear stability is guaranteed by the general condition (1.2) holding forall min {N}; (v) condition (1.2) is hidden in the Darcy-Boussinesq equations, it can be found by substituting the salt concentration fields via new suitable unknown fields and looking for symmetries and skew-symmetries in the new system of equations. The present paper - originating from Rionero ["Absence of subcritical instabilities and global nonlinear stability for porous ternary diffusive-convective fluid mixtures," Phys. Fluids 24, 104101 (2012)], 10.1063/1.4757858 - generalizes the properties (ii)-(iv) (obtained for m = 2) to any min {N} and furnishes the newly obtained properties (i) and (v). We stress the relevant physical meaning of (1.2). In fact (1.2) - in simple algebraic closed form - guarantees that the onset of convection cannot occur and appears to be useful not only for theoreticians but also for experimentalists in the research field of physics of fluids. Analogously, conditions guaranteeing the onset of convection - in simple algebraic closed form (cf. d6.18 d6.19">(6.18) and (6.19) reversed) - are furnished.
ERIC Educational Resources Information Center
Lavenda, Bernard H.
1985-01-01
Explains the phenomenon of Brownian motion, which serves as a mathematical model for random processes. Topics addressed include kinetic theory, Einstein's theory, particle displacement, and others. Points out that observations of the random course of a particle suspended in fluid led to the first accurate measurement of atomic mass. (DH)
Roy, C.; Ohana, M.; Host, Ph.; Alemann, G.; Labani, A.; Wattiez, A.; Lang, H.
2014-01-01
Objective The goal of this prospective study was to compare the efficiency of two types of MRU after diuretic administration to identify the non-dilated ureter. Methods MR pelvic examinations were performed in 126 patients after receiving furosemide. Each patient underwent in addition to their protocol for context, two types of MRU: 2D T2-weighted FSE (T2w-MRU) and 3D Gd T1-weighted GE (CE-MRU). Four segments were checked for each ureter. For the first part of the analysis, readers evaluated the whole image quality using a four points subjective scale and for the second part, they were asked to score separately each ureteral segment as present or absent. Results 1008 ureteral segments were checked. For the image quality, readers did not find any significant difference (3.8 ± 0.5 vs 3.6 ± 0.7, p value: 0.13) between MRU methods. The interobserver agreement was excellent with a κ correlation coefficient as high as 0.89 for T2w-MRU and 0.92 for CE-MRU, respectively. For the detection of the segments and considering the 9 rotations for the T2W MRU, there were no statistically significant differences between the two groups. Conclusion T2-weighted MRU with multiple orientations and diuretic is sufficient to identify the non-dilated ureter. It offers information on ureteral peristaltism. It can be suggested that this sequence is able to detect an initial obstruction before hydronephrosis occurs. PMID:26937423
Martin, Bryn A; Yiallourou, Theresia I; Pahlavian, Soroush Heidari; Thyagaraj, Suraj; Bunck, Alexander C; Loth, Francis; Sheffer, Daniel B; Kröger, Jan Robert; Stergiopulos, Nikolaos
2016-05-01
For the first time, inter-operator dependence of MRI based computational fluid dynamics (CFD) modeling of cerebrospinal fluid (CSF) in the cervical spinal subarachnoid space (SSS) is evaluated. In vivo MRI flow measurements and anatomy MRI images were obtained at the cervico-medullary junction of a healthy subject and a Chiari I malformation patient. 3D anatomies of the SSS were reconstructed by manual segmentation by four independent operators for both cases. CFD results were compared at nine axial locations along the SSS in terms of hydrodynamic and geometric parameters. Intraclass correlation (ICC) assessed the inter-operator agreement for each parameter over the axial locations and coefficient of variance (CV) compared the percentage of variance for each parameter between the operators. Greater operator dependence was found for the patient (0.19 < ICC < 0.99) near the craniovertebral junction compared to the healthy subject (ICC > 0.78). For the healthy subject, hydraulic diameter and Womersley number had the least variance (CV = ~2%). For the patient, peak diastolic velocity and Reynolds number had the smallest variance (CV = ~3%). These results show a high degree of inter-operator reliability for MRI-based CFD simulations of CSF flow in the cervical spine for healthy subjects and a lower degree of reliability for patients with Type I Chiari malformation. PMID:26446009
Silt motion simulation using finite volume particle method
NASA Astrophysics Data System (ADS)
Jahanbakhsh, E.; Vessaz, C.; Avellan, F.
2014-03-01
In this paper, we present a 3-D FVPM which features rectangular top-hat kernels. With this method, interaction vectors are computed exactly and efficiently. We introduce a new method to enforce the no-slip boundary condition. With this boundary enforcement, the interaction forces between fluid and wall are computed accurately. We employ the boundary force to predict the motion of rigid spherical silt particles inside the fluid. To validate the model, we simulate the 2-D sedimentation of a single particle in viscous fluid tank and compare results with benchmark data. The particle resolution is verified by convergence study. We also simulate the sedimentation of two particles exhibiting drafting, kissing and tumbling phenomena in 2-D and 3-D. We compare the results with other numerical solutions.
Energy Science and Technology Software Center (ESTSC)
2001-01-31
This software reduces the data from two-dimensional kSA MOS program, k-Space Associates, Ann Arbor, MI. Initial MOS data is recorded without headers in 38 columns, with one row of data per acquisition per lase beam tracked. The final MOSS 2d data file is reduced, graphed, and saved in a tab-delimited column format with headers that can be plotted in any graphing software.
NASA Astrophysics Data System (ADS)
Hide, Raymond
1998-01-01
Previous theoretical and laboratory studies of mechanically driven fluids in general rotation relative to an inertial frame have shown that there is a special class of flows for which the (Eulerian) flow field u( r, t) relative to the rotating frame of reference is unaffected by gyroscopic (Coriolis) forces, and therefore remains the same for all values of the rotation vector Ω. (Here t denotes time and r the position of a general point R in a reference frame attached to the rotating apparatus.) Such flows occur when (a) Ω is independent of time t; (b) u( r, t) is independent of the coordinate z (say) parallel to Ω, (c) the fluid has constant density and is therefore 'barotropic' (i.e. no density variations on horizontal surfaces) and (d) the topology of the cross-section of the (cylindrical) container, in planes z = constant, is such that the bounding surfaces can support the concomitant field of (kinematic) pressure P1 satisfying ▿ P 1 + 2 Ω × u = 0 Condition (d) is equivalent to the requirement that any fluid sources or siks within the system be multipole in character, but not monopole. In the present study the 'baroclinic' case is treated, where buoyancy forces due to the action of gravity (and centripetal forces) on horizontal density variations have to be taken into account. These include investigations of flows due entirely to buoyancy forces, such as thermal convection in fluids in rotating cylindrical containers of various shapes and topological characteristics subject to horizontal temperature gradients. The implications for the impressed temperature field of the mathematical requirements that the fields of kinematic pressure P1 and density overlineϱϑ (where overlineϱ denotes the mean density) be everywhere single-valued are guiding such investigations and facilitating the interpretation of their findings. The investigations include laboratory studies, reported elsewhere, of convection in a rotating fluid annulus with a circular cross
Nanoimprint lithography: 2D or not 2D? A review
NASA Astrophysics Data System (ADS)
Schift, Helmut
2015-11-01
Nanoimprint lithography (NIL) is more than a planar high-end technology for the patterning of wafer-like substrates. It is essentially a 3D process, because it replicates various stamp topographies by 3D displacement of material and takes advantage of the bending of stamps while the mold cavities are filled. But at the same time, it keeps all assets of a 2D technique being able to pattern thin masking layers like in photon- and electron-based traditional lithography. This review reports about 20 years of development of replication techniques at Paul Scherrer Institut, with a focus on 3D aspects of molding, which enable NIL to stay 2D, but at the same time enable 3D applications which are "more than Moore." As an example, the manufacturing of a demonstrator for backlighting applications based on thermally activated selective topography equilibration will be presented. This technique allows generating almost arbitrary sloped, convex and concave profiles in the same polymer film with dimensions in micro- and nanometer scale.