Velocity Gradient (velocity + gradient)

Distribution by Scientific Domains


Selected Abstracts


Evaluation of Left Ventricular Systolic and Diastolic Global Function: Peak Positive and Negative Myocardial Velocity Gradients in M-Mode Doppler Tissue Imaging

ECHOCARDIOGRAPHY, Issue 1 2002
Yoshiki Ueno M.D.
Objectives: To evaluate a new indicator of left ventricular global function: Myocardial velocity gradient (MVG) M-mode Doppler tissue imaging (DTI). Background: MVG is a new indicator of regional left ventricular function and global left ventricular diastolic function. However, it is unclear whether MVG also is an indicator of left ventricular global function in comparison with invasive indices. Methods: We performed conventional imaging and M-mode DTI in 85 subjects and calculated MVG at the posterior wall. We obtained satisfactory images in 65 subjects, who we divided into three groups: Noninvasive study group, invasive study group, and hemodialysis group. The noninvasive study group was divided into three subgroups (a younger normal subgroup, an older normal subgroup, and a cardiomyopathy subgroup), and MVG was compared with indices of conventional imaging. In the invasive study group, we compared MVG and indices of conventional imaging with hemodynamic data (peak positive and negative dp/dt, and the time constant T) using a high fidelity micromanometer-tipped catheter. In the hemodialysis group, we compared indices before hemodialysis with those after hemodialysis. Results: Peak positive MVG correlated well with peak positive dp/dt (r = 0.79), and this did not change with hemodialysis (P = 0.87). Peak negative MVG also correlated well with peak positive dp/dt and the time constant T (r = 0.88 and r = 0.80), and this did not change with hemodialysis (P = 0.97). Conclusions: Peak positive and negative MVG are sensitive and load-insensitive indicators of left ventricular function. [source]


Computational fluid dynamics modelling of boundary roughness in gravel-bed rivers: an investigation of the effects of random variability in bed elevation

EARTH SURFACE PROCESSES AND LANDFORMS, Issue 4 2001
A.P. Nicholas
Abstract Results from a series of numerical simulations of two-dimensional open-channel flow, conducted using the computational fluid dynamics (CFD) code FLUENT, are compared with data quantifying the mean and turbulent characteristics of open-channel flow over two contrasting gravel beds. Boundary roughness effects are represented using both the conventional wall function approach and a random elevation model that simulates the effects of supra-grid-scale roughness elements (e.g. particle clusters and small bedforms). Results obtained using the random elevation model are characterized by a peak in turbulent kinetic energy located well above the bed (typically at y/h,=,0·1,0·3). This is consistent with the field data and in contrast to the results obtained using the wall function approach for which maximum turbulent kinetic energy levels occur at the bed. Use of the random elevation model to represent supra-grid-scale roughness also allows a reduction in the height of the near-bed mesh cell and therefore offers some potential to overcome problems experienced by the wall function approach in flows characterized by high relative roughness. Despite these benefits, the results of simulations conducted using the random elevation model are sensitive to the horizontal and vertical mesh resolution. Increasing the horizontal mesh resolution results in an increase in the near-bed velocity gradient and turbulent kinetic energy, effectively roughening the bed. Varying the vertical resolution of the mesh has little effect on simulated mean velocity profiles, but results in substantial changes to the shape of the turbulent kinetic energy profile. These findings have significant implications for the application of CFD within natural gravel-bed channels, particularly with regard to issues of topographic data collection, roughness parameterization and the derivation of mesh-independent solutions. Copyright © 2001 John Wiley & Sons, Ltd. [source]


Evaluation of Left Ventricular Systolic and Diastolic Global Function: Peak Positive and Negative Myocardial Velocity Gradients in M-Mode Doppler Tissue Imaging

ECHOCARDIOGRAPHY, Issue 1 2002
Yoshiki Ueno M.D.
Objectives: To evaluate a new indicator of left ventricular global function: Myocardial velocity gradient (MVG) M-mode Doppler tissue imaging (DTI). Background: MVG is a new indicator of regional left ventricular function and global left ventricular diastolic function. However, it is unclear whether MVG also is an indicator of left ventricular global function in comparison with invasive indices. Methods: We performed conventional imaging and M-mode DTI in 85 subjects and calculated MVG at the posterior wall. We obtained satisfactory images in 65 subjects, who we divided into three groups: Noninvasive study group, invasive study group, and hemodialysis group. The noninvasive study group was divided into three subgroups (a younger normal subgroup, an older normal subgroup, and a cardiomyopathy subgroup), and MVG was compared with indices of conventional imaging. In the invasive study group, we compared MVG and indices of conventional imaging with hemodynamic data (peak positive and negative dp/dt, and the time constant T) using a high fidelity micromanometer-tipped catheter. In the hemodialysis group, we compared indices before hemodialysis with those after hemodialysis. Results: Peak positive MVG correlated well with peak positive dp/dt (r = 0.79), and this did not change with hemodialysis (P = 0.87). Peak negative MVG also correlated well with peak positive dp/dt and the time constant T (r = 0.88 and r = 0.80), and this did not change with hemodialysis (P = 0.97). Conclusions: Peak positive and negative MVG are sensitive and load-insensitive indicators of left ventricular function. [source]


A thermochemical boundary layer at the base of Earth's outer core and independent estimate of core heat flux

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2008
David Gubbins
SUMMARY Recent seismological observations suggest the existence of a ,150-km-thick density-stratified layer with a P -wave velocity gradient that differs slightly from PREM. Such a structure can only be caused by a compositional gradient, effects of a slurry or temperature being too small and probably the wrong sign. We propose a stably stratified, variable concentration layer on the liquidus. Heat is transported by conduction down the liquidus while the light and heavy components migrate through the layer by a process akin to zone refining, similar to the one originally proposed by Braginsky. The layer remains static in a frame of reference moving upwards with the expanding inner core boundary. We determine the gradient using estimates of co, the concentration in the main body of the outer core, and cb, the concentration of the liquid at the inner core boundary. We determine the depression of the melting point and concentrations using ideal solution theory and seismologically determined density jumps at the inner core boundary. We suppose that co determines ,,mod, the jump from normal mode eigenfrequencies that have long resolution lengths straddling the entire layer, and that cb determines ,,bod, the jump determined from body waves, which have fine resolution. A simple calculation then yields the seismic, temperature, and concentration profiles within the layer. Comparison with the distance to the C-cusp of PKP and normal mode eigenfrequencies constrain the model. We explore a wide range of possible input parameters; many fail to predict sensible seismic properties and heat fluxes. A model with ,,mod= 0.8 gm cc,1, ,,bod= 0.6 gm cc,1, and layer thickness 200 km is consistent with the seismic observations and can power the geodynamo with a reasonable inner core heat flux of ,2 TW and nominal inner core age of ,1 Ga. It is quite remarkable and encouraging that a model based on direct seismic observations and simple chemistry can predict heat fluxes that are comparable with those derived from recent core thermal history calculations. The model also provides plausible explanations of the observed seismic layer and accounts for the discrepancy between estimates of the inner core density jumps derived from body waves and normal modes. [source]


Geodetic observations of ice flow velocities over the southern part of subglacial Lake Vostok, Antarctica, and their glaciological implications

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2006
Jens Wendt
SUMMARY In the austral summer seasons 2001/02 and 2002/03, Global Positioning System (GPS) data were collected in the vicinity of Vostok Station to determine ice flow velocities over Lake Vostok. Ten GPS sites are located within a radius of 30 km around Vostok Station on floating ice as well as on grounded ice to the east and to the west of the lake. Additionally, a local deformation network around the ice core drilling site 5G-1 was installed. The derived ice flow velocity for Vostok Station is 2.00 m a,1± 0.01 m a,1. Along the flowline of Vostok Station an extension rate of about 10,5 a,1 (equivalent to 1 cm km,1 a,1) was determined. This significant velocity gradient results in a new estimate of 28 700 years for the transit time of an ice particle along the Vostok flowline from the bedrock ridge in the southwest of the lake to the eastern shoreline. With these lower velocities compared to earlier studies and, hence, larger transit times the basal accretion rate is estimated to be 4 mm a,1 along a portion of the Vostok flowline. An assessment of the local accretion rate at Vostok Station using the observed geodetic quantities yields an accretion rate in the same order of magnitude. Furthermore, the comparison of our geodetic observations with results inferred from ice-penetrating radar data indicates that the ice flow may not have changed significantly for several thousand years. [source]


Crustal structure of central and northern Iceland from analysis of teleseismic receiver functions

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 1 2000
Fiona A. Darbyshire
We present results from a teleseismic receiver function study of central and northern Iceland, carried out during the period 1995,1998. Data from eight broad-band seismometers installed in the SIL network operated by the Icelandic Meteorological Office were used for analysis. Receiver functions for each station were generated from events for a wide range of backazimuths and a combination of inversion and forward modelling was used to infer the crustal structure below each station. The models generated show a considerable variation in the nature and thickness of the crust across Iceland. The thinnest crust (20,21 km) is found in the northern half of the Northern Volcanic Zone approximately 120 km north of the centre of the Iceland mantle plume. Thicker crust (24,30 km) is found elsewhere in northern and central Iceland and the thickest crust (37 km) is found close to the plume centre. Velocity,depth profiles show a distinct division of the crust into two main sections, an upper high-velocity-gradient section of thickness 2,8 km and a lower crustal section with small or zero overall velocity gradient. The thickness of the upper crust correlates with the tectonic structure of Iceland; the upper crust is thickest on the flanks of the northern and central volcanic rift zones and thinnest close to active or extinct central volcanoes. Below the Krafla central volcano in northeastern Iceland the receiver function models show a prominent low-velocity zone at 10,15 km depth with minimum shear wave velocities of 2.0,2.5 km s,1. We suggest that this feature results from the presence of partially molten sills in the lower crust. Less prominent low-velocity zones found in other regions of Iceland may arise from locally high temperatures in the crust or from acidic intrusive bodies at depth. A combination of the receiver function results and seismic refraction results constrains the crustal thickness across a large part of Iceland. Melting by passive decompression of the hot mantle below the rift zone in northern Iceland forms a crust of thickness ,20 km. In contrast, the larger crustal thickness below central Iceland probably arises from enhanced melt production due to active upwelling in the plume core. [source]


The effects of near-surface conditions on anisotropy parameter estimations from 4C seismic data

GEOPHYSICAL PROSPECTING, Issue 1 2006
Bärbel Traub
ABSTRACT We present a study of anisotropic parameter estimation in the near-surface layers for P-wave and converted-wave (C-wave) data. Near-surface data is affected by apparent anisotropy due to a vertical velocity compaction gradient. We have carried out a modelling study, which showed that a velocity gradient introduces apparent anisotropy into an isotropic medium. Thus, parameter estimation will give anomalous values that affect the imaging of the target area. The parameter estimation technique is also influenced by phase reversals with diminishing amplitude, leading to erroneous parameters. In a modelling study using a near-surface model, we have observed phase reversals in near-surface PP reflections. The values of the P-wave anisotropy parameter , estimated from these events are about an order of magnitude larger than the model values. Next, we use C-wave data to estimate the effect of anisotropy (,) and compute , from these values. These calculated ,-values are closer to the model values, and NMO correction with both ,-values shows a better correction for the calculated value. Hence, we believe that calculating , from , gives a better representation of the anisotropy than picked , from the P-wave. Finally, we extract the anisotropy parameters , and , from real data from the Alba Field in the North Sea. Comparing the results with reference values from a model built according to well-log, VSP and surface data, we find that the parameters show differences of up to an order of magnitude. The ,-values calculated from the C-wave anisotropy parameter , fit the reference values much better and show values of the same order of magnitude. [source]


Coupled Navier,Stokes,Molecular dynamics simulations using a multi-physics flow simulation framework

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 10 2010
R. Steijl
Abstract Simulation of nano-scale channel flows using a coupled Navier,Stokes/Molecular Dynamics (MD) method is presented. The flow cases serve as examples of the application of a multi-physics computational framework put forward in this work. The framework employs a set of (partially) overlapping sub-domains in which different levels of physical modelling are used to describe the flow. This way, numerical simulations based on the Navier,Stokes equations can be extended to flows in which the continuum and/or Newtonian flow assumptions break down in regions of the domain, by locally increasing the level of detail in the model. Then, the use of multiple levels of physical modelling can reduce the overall computational cost for a given level of fidelity. The present work describes the structure of a parallel computational framework for such simulations, including details of a Navier,Stokes/MD coupling, the convergence behaviour of coupled simulations as well as the parallel implementation. For the cases considered here, micro-scale MD problems are constructed to provide viscous stresses for the Navier,Stokes equations. The first problem is the planar Poiseuille flow, for which the viscous fluxes on each cell face in the finite-volume discretization are evaluated using MD. The second example deals with fully developed three-dimensional channel flow, with molecular level modelling of the shear stresses in a group of cells in the domain corners. An important aspect in using shear stresses evaluated with MD in Navier,Stokes simulations is the scatter in the data due to the sampling of a finite ensemble over a limited interval. In the coupled simulations, this prevents the convergence of the system in terms of the reduction of the norm of the residual vector of the finite-volume discretization of the macro-domain. Solutions to this problem are discussed in the present work, along with an analysis of the effect of number of realizations and sample duration. The averaging of the apparent viscosity for each cell face, i.e. the ratio of the shear stress predicted from MD and the imposed velocity gradient, over a number of macro-scale time steps is shown to be a simple but effective method to reach a good level of convergence of the coupled system. Finally, the parallel efficiency of the developed method is demonstrated. Copyright © 2009 John Wiley & Sons, Ltd. [source]


Flocculation analysis and control system by laser diffractometry at industrial scale

AICHE JOURNAL, Issue 1 2008
Mohamed Bizi
Abstract The flocculant injection control system efficiency was evaluated on-site in an aggregate quarry, by means of laser diffraction and analysis of the size and texture of flocs. The configuration of the feed tank and the laser particle size measurement cell installed at the facility (280,000 t/year of aggregates) allowed characterization of flocs with particle size between 4.2 and 3473 ,m under hydrodynamic conditions that were highly favorable for the examination of large and very fragile flocs. Two days of analysis of the floc formation process along the path followed by the slurry showed that flocculation was optimal during standard operation of the facility when the flow rate of waste fines and concentrations of solids were close to those used to calibrate the flocculant injection control system. Conversely, when the concentration of solids in the flocculator feed slurry dropped by 57.3%, the flocculant dosing fluctuated during stabilization of the mechanism, the kinetics of flocculation slowed, the mean size of flocs arriving at the settling tank dropped by 69%, and the mode of smaller flocs shifted from 77.8 to 10.4 ,m in relation to normal operation. On-site analyses confirmed that the measurements made with laser diffraction (using a methodology developed in the laboratory) allow determination of the effects of conditioning on the characteristics of flocs in terms of particle size, porosity, density, and volume fraction in the slurry. Evolution of these characteristics according to the local parameters of conditioning (mean detention time, mean slurry velocity, and mean velocity gradient) provides a significant part of basic information necessary to a diagnosis of the operation of an industrial circuit of flocculation. © 2007 American Institute of Chemical Engineers AIChE J, 2008 [source]


Exponential attractor for a planar shear-thinning flow

MATHEMATICAL METHODS IN THE APPLIED SCIENCES, Issue 17 2007
Dalibor Pra
Abstract We study the dynamics of an incompressible, homogeneous fluid of a power-law type, with the stress tensor T = ,(1 + µ|Dv|)p,2Dv, where Dv is a symmetric velocity gradient. We consider the two-dimensional problem with periodic boundary conditions and p , (1, 2). Under these assumptions, we estimate the fractal dimension of the exponential attractor, using the so-called method of ,,-trajectories. Copyright © 2007 John Wiley & Sons, Ltd. [source]


Discovery of large-scale methanol and hydroxyl maser filaments in W3(OH)

MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY, Issue 4 2006
L. Harvey-Smith
ABSTRACT Images of the 6.7-GHz methanol maser emission from W3(OH) made at 50- and 100-mas angular resolution with the Multi-Element Radio-Linked Interferometer Network (MERLIN) are presented. The masers lie across the western face of the ultracompact H ii region in extended filaments which may trace large-scale shocks. There is a complex interrelation between the 6.7-GHz methanol masers and hydroxyl (OH) masers at 1.7 and 4.7 GHz. Together the two species trace an extended filamentary structure that stretches at least 3100 au across the face of the ultracompact H ii region. The dominant 6.7-GHz methanol emission coincides with the radio continuum peak and is populated by masers with broad spectral lines. The 6.7-GHz methanol emission is elongated at position angle 50° with a strong velocity gradient, and bears many similarities to the methanol maser disc structure reported in NGC 7538. It is surrounded by arcs of ground state OH masers at 1.7 GHz and highly excited OH masers at 13.44 GHz, some of which have the brightest methanol masers at their focus. We suggest that this region hosts the excitation centre for the ultracompact H ii region. [source]


OH megamasers, starburst and AGN activity in Markarian 231

MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY, Issue 2 2005
A. M. S. Richards
ABSTRACT We present Multi-Element Radio-Linked Interferometer Network (MERLIN) observations of OH maser and radio continuum emission within a few hundred pc of the core of the ultraluminous infrared galaxy (ULIRG) Markarian 231. This is the only known OH megamaser galaxy classed as a Seyfert 1. Maser emission is identified with the 1665- and 1667-MHz transitions over a velocity extent of 720 km s,1. Both lines show a similar position,velocity structure including a gradient of 1.7 km s,1 pc,1 from NW to SE along the 420-pc major axis. The (unresolved) inner few tens of pc possess a much steeper velocity gradient. The maser distribution is modelled as a torus rotating about an axis inclined at ,45° to the plane of the sky. We estimate the enclosed mass density to be 320 ± 90 M, pc,3 in a flattened distribution. This includes a central unresolved mass of ,8 × 106 M,. All the maser emission is projected against a region with a radio continuum brightness temperature ,105 K, giving a maser gain of ,2.2. The 1667:1665 MHz line ratio is close to 1.8 (the value predicted for thermal emission) consistent with radiatively pumped, unsaturated masers. This behaviour and the kinematics of the torus suggest that the size of individual masing regions is in the range 0.25,4 pc with a covering factor close to unity. There are no very bright compact masers, in contrast to galaxies such as the Seyfert 2 Markarian 273, where the masing torus is viewed nearer edge-on. The comparatively modest maser amplification seen from Markarian 231 is consistent with its classification in the unification scheme for Seyfert galaxies. Most of the radio continuum emission on 50,500 pc scales is probably of starburst origin but the compact peak is 0.4 per cent polarized by a magnetic field running north,south, similar to the jet direction on these scales. There is no close correlation between maser and continuum intensity, suggesting that much of the radio continuum must originate in the foreground and indeed the relative continuum brightness is slightly greater in the direction of the approaching jet. Comparisons with other data show that the jet changes direction close the nucleus and suggest that the sub-kpc disc hosting the masers and starburst activity is severely warped. [source]


A high-resolution radio study of neutral gas in the starburst galaxy NGC 520

MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY, Issue 2 2003
R. J. Beswick
ABSTRACT We present subarcsec angular resolution observations of the neutral gas in the nearby starburst galaxy NGC 520. The central kpc region of NGC 520 contains an area of significantly enhanced star formation. The radio continuum structure of this region resolves into ,10 continuum components. By comparing the flux densities of the brightest of these components at 1.4 GHz with published 15-GHz data we infer that these components detected at 1.4 and 1.6 GHz are related to the starburst and are most likely to be collections of several supernova remnants within the beam. None of these components is consistent with emission from an active galactic nuclei. Both neutral hydrogen (H i) and hydroxyl (OH) absorption lines are observed against the continuum emission, along with a weak OH maser feature probably related to the star formation activity in this galaxy. Strong H i absorption (NH, 1022 atoms cm,2) traces a velocity gradient of 0.5 km s,1 pc,1 across the central kpc of NGC 520. The H i absorption velocity structure is consistent with the velocity gradients observed in both the OH absorption and in CO emission observations. The neutral gas velocity structure observed within the central kpc of NGC 520 is attributed to a kpc-scale ring or disc. It is also noted that the velocity gradients observed for these neutral gas components appear to differ with the velocity gradients observed from optical ionized emission lines. This apparent disagreement is discussed and attributed to the extinction of the optical emission from the actual centre of this source hence implying that optical ionized emission lines are only detected from regions with significantly different radii to those sampled by the observations presented here. [source]


X-ray beaming caused by resonance scattering in the accretion column of magnetic cataclysmic variables

MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY, Issue 1 2001
Y. Terada
Extremely strong ionized Fe emission lines, with equivalent widths reaching ,4000 eV, were discovered by ASCA from a few Galactic compact objects, including AX J2315,0592, RX J1802.1+1804 and AX J1842.8,0423. These objects are thought to be binary systems containing magnetized white dwarfs (WDs). A possible interpretation of the strong Fe K line is the line-photon collimation in the WD accretion column, as a result of resonance scattering of line photons. The collimation occurs when the accretion column has a flat shape, and the effect is augmented by the vertical velocity gradient, which reduces the resonant trapping of resonant photons along the magnetic field lines. This effect was quantitatively confirmed with Monte Carlo simulations. Furthermore, with ASCA observations of the polar V834 Centauri, this collimation effect was clearly detected as a rotational modulation of the equivalent width of the Fe K emission line. The extremely strong emission lines mentioned above can be explained consistently by our interpretation. Combining this effect with other X-ray information, the geometry and plasma parameters in the accretion column were determined. [source]


Ultra-high-resolution observations of circumstellar K i and C2 around the post-AGB star HD 56126

MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY, Issue 2 2000
I. A. Crawford
We have used the Ultra-High-Resolution Facility (UHRF) at the AAT, operating at a resolution of 0.35 km s,1 (FWHM), to observe K i and C2 absorption lines arising in the circumstellar environment of the post-AGB star HD 56126. We find three narrow circumstellar absorption components in K i, two of which are also present in C2. We attribute this velocity structure to discrete shells resulting from multiple mass-loss events from the star. The very high spectral resolution has enabled us to resolve the intrinsic linewidths of these narrow lines for the first time, and we obtain velocity dispersions (b -values) of 0.2,0.3 km s,1 for the K i components, and 0.54±0.03 km s,1 for the strongest (and best defined) C2 component. These correspond to rigorous kinetic temperature upper limits of 211 K for K i and 420 K for C2, although the b -value ratio implies that these two species do not co-exist spatially. The observed degree of rotational excitation of C2 implies low kinetic temperatures (Tk,10 K) and high densities (n,106 to 107 cm,3) within the shell responsible for the main C2 component. Given this low temperature, the line profiles then imply either mildly supersonic turbulence or an unresolved velocity gradient through the shell. [source]


Regional GPS data confirm high strain accumulation prior to the 2000 June 4 Mw= 7.8 earthquake at southeast Sumatra

GEOPHYSICAL JOURNAL INTERNATIONAL, Issue 3 2001
G. W. Michel
Summary Site velocities derived from repeated measurements in a regional GPS network in Southeast Asia help to constrain the motion of tectonic blocks as well as slip rates along major faults in the area. Using 3-D forward dislocation modelling, the influence of seismic elastic loading and unloading on the measured site motions are approximated. Results suggest that the northwestern Sunda arc is fully coupled seismogenically, whereas its eastern part along Java shows localized deformation. Higher horizontal velocity gradients than expected from the modelling of a fully coupled plate interface west of Manila in the Philippines suggest that deformation may be localized there. Assuming that geodetically derived convergence represents long-term rates, accumulated geodetic moments are compared to those derived using seismic data from 1977 to 2000 (Harvard CMT catalogue). If areas displaying localized deformation are dominated by creep processes, the largest difference between accumulated and seismically released deformation is located where the 2000 June 4 Mw = 7.8 Sumatra earthquake occurred. [source]


Non-uniqueness with refraction inversion , the Mt Bulga shear zone

GEOPHYSICAL PROSPECTING, Issue 4 2010
Derecke Palmer
ABSTRACT The tau-p inversion algorithm is widely employed to generate starting models with many computer programs that implement refraction tomography. However, this algorithm can frequently fail to detect even major lateral variations in seismic velocities, such as a 50 m wide shear zone, which is the subject of this study. By contrast, the shear zone is successfully defined with the inversion algorithms of the generalized reciprocal method. The shear zone is confirmed with a 2D analysis of the head wave amplitudes, a spectral analysis of the refraction convolution section and with numerous closely spaced orthogonal seismic profiles recorded for a later 3D refraction investigation. Further improvements in resolution, which facilitate the recognition of additional zones with moderate reductions in seismic velocity, are achieved with a novel application of the Hilbert transform to the refractor velocity analysis algorithm. However, the improved resolution also requires the use of a lower average vertical seismic velocity, which accommodates a velocity reversal in the weathering. The lower seismic velocity is derived with the generalized reciprocal method, whereas most refraction tomography programs assume vertical velocity gradients as the default. Although all of the tomograms are consistent with the traveltime data, the resolution of each tomogram is comparable only with that of the starting model. Therefore, it is essential to employ inversion algorithms that can generate detailed starting models, where detailed lateral resolution is the objective. Non-uniqueness can often be readily resolved with head wave amplitudes, attribute processing of the refraction convolution section and additional seismic traverses, prior to the acquisition of any borehole data. It is concluded that, unless specific measures are taken to address non-uniqueness, the production of a single refraction tomogram that fits the traveltime data to sufficient accuracy does not necessarily demonstrate that the result is either correct, or even the most probable. [source]


A hybrid fast algorithm for first arrivals tomography

GEOPHYSICAL PROSPECTING, Issue 5 2009
Manuela Mendes
ABSTRACT A hybrid algorithm, combining Monte-Carlo optimization with simultaneous iterative reconstructive technique (SIRT) tomography, is used to invert first arrival traveltimes from seismic data for building a velocity model. Stochastic algorithms may localize a point around the global minimum of the misfit function but are not suitable for identifying the precise solution. On the other hand, a tomographic model reconstruction, based on a local linearization, will only be successful if an initial model already close to the best solution is available. To overcome these problems, in the method proposed here, a first model obtained using a classical Monte Carlo-based optimization is used as a good initial guess for starting the local search with the SIRT tomographic reconstruction. In the forward problem, the first-break times are calculated by solving the eikonal equation through a velocity model with a fast finite-difference method instead of the traditional slow ray-tracing technique. In addition, for the SIRT tomography the seismic energy from sources to receivers is propagated by applying a fast Fresnel volume approach which when combined with turning rays can handle models with both positive and negative velocity gradients. The performance of this two-step optimization scheme has been tested on synthetic and field data for building a geologically plausible velocity model. This is an efficient and fast search mechanism, which permits insertion of geophysical, geological and geodynamic a priori constraints into the grid model and ray path is completed avoided. Extension of the technique to 3D data and also to the solution of ,static correction' problems is easily feasible. [source]


Migration velocity analysis for tilted transversely isotropic media

GEOPHYSICAL PROSPECTING, Issue 1 2009
Laxmidhar Behera
ABSTRACT Tilted transversely isotropic formations cause serious imaging distortions in active tectonic areas (e.g., fold-and-thrust belts) and in subsalt exploration. Here, we introduce a methodology for P-wave prestack depth imaging in tilted transversely isotropic media that properly accounts for the tilt of the symmetry axis as well as for spatial velocity variations. For purposes of migration velocity analysis, the model is divided into blocks with constant values of the anisotropy parameters , and , and linearly varying symmetry-direction velocity VP0 controlled by the vertical (kz) and lateral (kx) gradients. Since determination of tilt from P-wave data is generally unstable, the symmetry axis is kept orthogonal to the reflectors in all trial velocity models. It is also assumed that the velocity VP0 is either known at the top of each block or remains continuous in the vertical direction. The velocity analysis algorithm estimates the velocity gradients kz and kx and the anisotropy parameters , and , in the layer-stripping mode using a generalized version of the method introduced by Sarkar and Tsvankin for factorized transverse isotropy with a vertical symmetry axis. Synthetic tests for several models typical in exploration (a syncline, uptilted shale layers near a salt dome and a bending shale layer) confirm that if the symmetry-axis direction is fixed and VP0 is known, the parameters kz, kx, , and , can be resolved from reflection data. It should be emphasized that estimation of , in tilted transversely isotropic media requires using nonhyperbolic moveout for long offsets reaching at least twice the reflector depth. We also demonstrate that application of processing algorithms designed for a vertical symmetry axis to data from tilted transversely isotropic media may lead to significant misfocusing of reflectors and errors in parameter estimation, even when the tilt is moderate (30°). The ability of our velocity analysis algorithm to separate the anisotropy parameters from the velocity gradients can be also used in lithology discrimination and geologic interpretation of seismic data in complex areas. [source]


Comparative study of lattice-Boltzmann and finite volume methods for the simulation of laminar flow through a 4:1 planar contraction

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 9 2004
Yarub Y. Al-Jahmany
Abstract In the present paper, a comparative study of numerical solutions for Newtonian fluids based on the lattice-Boltzmann method (LBM) and the classical finite volume method (FVM) is presented for the laminar flow through a 4:1 planar contraction at a Reynolds number of value one, Re=1. In this study, the stress field for LBM is directly obtained from the distribution function. The calculations of the stress based on the FVM-data use the evaluations of velocity gradients with finite differences. The stress field for both LBM and FVM is expressed in the present study in terms of the shear stress and the first normal stress difference. The lateral and axial profiles of the velocity, the shear stress and the first normal stress difference for both methods are investigated. It is shown that the LBM results for the velocity and the stresses are in excellent agreement with the FVM results. Copyright © 2004 John Wiley & Sons, Ltd. [source]


Godunov-type adaptive grid model of wave,current interaction at cuspate beaches

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 6 2004
Benedict D. Rogers
Abstract This paper presents a second-order accurate Godunov-type numerical scheme for depth- and period-averaged wave,current interaction. A flux Jacobian is derived for the wave conservation equations and its eigensystem determined, enabling Roe's approximate Riemann solver to be used to evaluate convective fluxes. Dynamically adaptive quadtree grids are used to focus on local hydrodynamic features, where sharp gradients occur in the flow variables. Adaptation criteria based on depth-averaged vorticity, wave-height gradient, wave steepness and the magnitude of velocity gradients are found to produce accurate solutions for nearshore circulation at a half-sinusoidal beach. However, the simultaneous combination of two or more separate criteria produces numerical instability and interference unless all criteria are satisfied for mesh depletion. Simulations of wave,current interaction at a multi-cusped beach match laboratory data from the United Kingdom Coastal Research Facility (UKCRF). A parameter study demonstrates the sensitivity of nearshore flow patterns to changes in relative cusp height, angle of wave incidence, bed roughness, offshore wave height and assumed turbulent eddy viscosity. Only a small deviation from normal wave incidence is required to initiate a meandering longshore current. Nearshore circulation patterns are highly dependent on the offshore wave height. Reduction of the assumed eddy viscosity parameter causes the primary circulation cells for normally incident waves to increase in strength whilst producing rip-like currents cutting diagonally across the surf zone. Copyright © 2004 John Wiley & Sons, Ltd. [source]


Transition from vortex to wall driven turbulence production in the Taylor,Couette system with a rotating inner cylinder

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 3 2002
W. M. J. Batten
Abstract Axisymmetrically stable turbulent Taylor vortices between two concentric cylinders are studied with respect to the transition from vortex to wall driven turbulent production. The outer cylinder is stationary and the inner cylinder rotates. A low Reynolds number turbulence model using the k - , formulation, facilitates an analysis of the velocity gradients in the Taylor,Couette flow. For a fixed inner radius, three radius ratios 0.734, 0.941 and 0.985 are employed to identify the Reynolds number range at which this transition occurs. At relatively low Reynolds numbers, turbulent production is shown to be dominated by the outflowing boundary of the Taylor vortex. As the Reynolds number increases, shear driven turbulence (due to the rotating cylinder) becomes the dominating factor. For relatively small gaps turbulent flow is shown to occur at Taylor numbers lower than previously reported. Copyright © 2002 John Wiley & Sons, Ltd. [source]


Hydrodynamic investigation of USP dissolution test apparatus II

JOURNAL OF PHARMACEUTICAL SCIENCES, Issue 9 2007
Ge Bai
Abstract The USP Apparatus II is the device commonly used to conduct dissolution testing in the pharmaceutical industry. Despite its widespread use, dissolution testing remains susceptible to significant error and test failures, and limited information is available on the hydrodynamics of this apparatus. In this work, laser-Doppler velocimetry (LDV) and computational fluid dynamics (CFD) were used, respectively, to experimentally map and computationally predict the velocity distribution inside a standard USP Apparatus II under the typical operating conditions mandated by the dissolution test procedure. The flow in the apparatus is strongly dominated by the tangential component of the velocity. Secondary flows consist of an upper and lower recirculation loop in the vertical plane, above and below the impeller, respectively. A low recirculation zone was observed in the lower part of the hemispherical vessel bottom where the tablet dissolution process takes place. The radial and axial velocities in the region just below the impeller were found to be very small. This is the most critical region of the apparatus since the dissolving tablet will likely be at this location during the dissolution test. The velocities in this region change significantly over short distances along the vessel bottom. This implies that small variations in the location of the tablet on the vessel bottom caused by the randomness of the tablet descent through the liquid are likely to result in significantly different velocities and velocity gradients near the tablet. This is likely to introduce variability in the test. © 2007 Wiley-Liss, Inc. and the American Pharmacists Association J Pharm Sci 96: 2327,2349, 2007 [source]


A high-resolution radio study of neutral gas in the starburst galaxy NGC 520

MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY, Issue 2 2003
R. J. Beswick
ABSTRACT We present subarcsec angular resolution observations of the neutral gas in the nearby starburst galaxy NGC 520. The central kpc region of NGC 520 contains an area of significantly enhanced star formation. The radio continuum structure of this region resolves into ,10 continuum components. By comparing the flux densities of the brightest of these components at 1.4 GHz with published 15-GHz data we infer that these components detected at 1.4 and 1.6 GHz are related to the starburst and are most likely to be collections of several supernova remnants within the beam. None of these components is consistent with emission from an active galactic nuclei. Both neutral hydrogen (H i) and hydroxyl (OH) absorption lines are observed against the continuum emission, along with a weak OH maser feature probably related to the star formation activity in this galaxy. Strong H i absorption (NH, 1022 atoms cm,2) traces a velocity gradient of 0.5 km s,1 pc,1 across the central kpc of NGC 520. The H i absorption velocity structure is consistent with the velocity gradients observed in both the OH absorption and in CO emission observations. The neutral gas velocity structure observed within the central kpc of NGC 520 is attributed to a kpc-scale ring or disc. It is also noted that the velocity gradients observed for these neutral gas components appear to differ with the velocity gradients observed from optical ionized emission lines. This apparent disagreement is discussed and attributed to the extinction of the optical emission from the actual centre of this source hence implying that optical ionized emission lines are only detected from regions with significantly different radii to those sampled by the observations presented here. [source]


The Effects of Power Law Fluid Rheology on Coil Heat Transfer for Agitated Vessel in the Transitional Flow Regime

THE CANADIAN JOURNAL OF CHEMICAL ENGINEERING, Issue 6 2004
Eric Ricci
Abstract A CFD model of heat transfer from power-law fluids to helical cooling coils in the transitional flow regime of a baffled tank mixed with a pitched blade turbine was developed with FluentTM. The model captured local temperature and velocity gradients. Simulations were run, varying Re, Pr, K and n. The results indicate that a Sieder-Tate type correlation, with the exponent on and the coefficient in front of the Reynolds number being a function of n, is recommended for estimating ho. Also, a new two coil bank design was found to be more efficient when 450 < Re < 650. On a développé à l'aide de FluentTM un modèle CFD de transfert de chaleur pour un serpentin hélicoïdal dans un réservoir à contrepales agité par une turbine à pales inclinées dans le régime de transition, appliqué au refroidissement de fluides de loi de puissance. Le modèle permet de déterminer la température locale et les gradients de vitesse. On a effectué des simulations en variant Re, Pr, K et n. Les résultats indiquent qu'une corrélation de type Sieder-Tate, où l'exposant et le coefficient prémultiplicateur du nombre de Reynolds sont fonction de n, est recommandée pour l'estimation de ho. On a également trouvé qu'une nouvelle conception de réservoir à deux serpentins était plus efficace lorsque 450 < Re < 650. [source]


Kinetic responses of Dunaliella in moving fluids

BIOTECHNOLOGY & BIOENGINEERING, Issue 1 2010
Ahammed Anwar Chengala
Abstract The objective of this work was to quantify the kinetic behavior of Dunaliella primolecta (D. primolecta) subjected to controlled fluid flow under laboratory conditions. In situ velocities of D. primolecta were quantified by micron-resolution particle image velocimetry and particle tracking velocimetry. Experiments were performed under a range of velocity gradients and corresponding energy dissipation levels at microscopic scales similar to the energy dissipation levels of natural aquatic ecosystems. An average swimming velocity of D. primolecta in a stagnant fluid was 41,µm/s without a preferential flow direction. In a moving fluid, the sample population velocities of D. primolecta follow a log-normal distribution. The variability of sample population velocities was maximal at the highest fluid flow velocity in the channel. Local fluid velocity gradients inhibited the accrual of D. primolecta by twofold 5 days after the initiation of the experiment in comparison to the non-moving fluid control experiment. Biotechnol. Bioeng. 2010;107: 65,75. © 2010 Wiley Periodicals, Inc. [source]