Library

Notes: Aims: To describe the clinicopathological and immunophenotypic features of 25 cases of Kikuchi–Fujimoto disease (K–F), which remains a poorly recognized entity and is still frequently confused with malignant lymphoma, and to discuss the main diagnostic problems experienced by the referring pathologist.〈section xml:id="abs1-1"〉〈title type="main"〉Methods and resultsHaematoxylin and eosin sections of 27 lymph node biopsies were re-examined. Immunostains for B-lymphocytes, T-lymphocytes and macrophages were performed. Clinical and follow-up data were obtained through a questionnaire to the referring pathologist or from the patients' notes where available. The suggested initial diagnoses are discussed. The lymph nodes showed a necrotizing process characterized by patchy or confluent areas of necrosis associated with karyorrhexis and absence or paucity of granulocytes. This was associated with a proliferation of large blastic cells consisting of a mixture of T-lymphocytes and histiocytes. Fragmentation of the biopsy was a frequent feature. The diagnosis of K–F was suggested by the referring pathologist in three cases only. The most common suggested diagnosis was that of a non-Hodgkin's lymphoma.〈section xml:id="abs1-2"〉〈title type="main"〉ConclusionThis series documents continuing difficulties in the diagnosis of Kikuchi–Fujimoto disease in the UK and emphasizes that cases are still being mistakenly diagnosed as malignant lymphomas. The diagnosis of Kikuchi–Fujimoto disease merits active consideration in any nodal biopsy showing fragmentation, necrosis and karyorrhexis, especially in young women presenting with cervical lymphadenopathy.

Notes: Particle-laden turbulent flows, at average volume fraction less than 4×10−4, in open channels are numerically simulated by using a pseudospectral method. The motion of particles, that are large compared with the dissipative length scale, is coupled to the fluid motion by a method that generates a "virtual" no-slip boundary on the particle surface by imposition of an external force field on the grid-points enclosed by the particle. Cases for both moving and stationary particles, lying on the wall, are simulated. The investigations focus on particle-turbulence interaction. It is found that particles increase turbulence intensities and Reynolds stress. By examining higher order turbulence statistics and doing a quadrant analysis of the Reynolds stress, it is found that the ejection-sweep cycle is affected—primarily through suppression of sweeps by the smaller particles and enhancement of sweep activity by the larger particles. An assessment of the impact of these findings on scalar transfer is made, as enhancement of wall heat/mass transfer rates is a motivation of the overall work on this subject. In the cases considered, comparison of the calculations with an existing experiment was possible, and shows good agreement. At present, due to limitations in available computational resources, this method cannot be used when the particle diameter is smaller than the smallest turbulence scale (e.g. the Kolmogorov length scale) and the volume fraction is of the same order as studied in this paper, i.e. between 10−3 and 10−4. © 1997 American Institute of Physics.

Notes: This paper presents the results of a numerical investigation of the effects of near-neutral density solid particles on turbulent liquid flow in a channel. Interactions of particles, in a size range about the dissipative length scale, with wall turbulence have been simulated at low volume fractions (average volume fraction less than 4×10−4). Fluid motion is calculated by directly solving the Navier-Stokes equations by a pseudo-spectral method and resolving all scales of motion. Particles are moved in a Lagrangian frame through the action of forces imposed by the fluid and gravity. Particle effects on fluid motion are fed back at each time step by calculating the velocity disturbance caused by the particles assuming the flow around them is locally Stokesian. Particle-particle interactions are not considered. The slightly heavier-than-fluid particles of the size range considered are found to preferentially accumulate in the low-speed streaks, as reported in several other investigations. It is also found that particles smaller than the dissipative length scale reduce turbulence intensities and Reynolds stress, whereas particles that are somewhat larger increase intensities and stress. By examining higher order turbulence statistics and doing a quadrant analysis of the Reynolds stress, it is found that the ejection-sweep cycle is affected—primarily through suppression of sweeps by the smaller particles and enhancement of sweep activity by the larger particles. A preliminary assessment of the impact of these findings on scalar transfer is made, as enhancement of transfer rate is a motivation of the overall work on this subject. For the case investigated, comparison of the calculations with an existing experiment was possible, and shows good agreement. © 1996 American Institute of Physics.

Notes: Direct numerical simulations of open-channel flow indicate that turbulence at the free surface contains large-scale persistent structures. They are "upwellings'' caused by impingement of bursts emanating from the bottom boundary; "downdrafts'' in regions where adjacent upwellings interact, and whirlpool-like "attached vortices'' which form at the edge of upwellings. The attached vortices are particularly long-lived in the sense that once formed, unless destroyed by other upwellings, they tend to interact with each other and dissipate only slowly. If turbulence generation at the bottom wall is turned off by changing the boundary condition to free slip, then the upwellings (related to bursts) and downdrafts no longer form. The dominant structures at the free surface become the attached vortices which pair, merge, and slowly dissipate. In the central regions, as expected, the structure remains three dimensional throughout the decay process. Near the free surface, the structure appears to be quasi- two dimensional, as indicated by quantitative measures such as energy spectra, interwave number energy transfer, invariants of the anisotropy tensor, and length scales. In the decaying case, the quasi-two-dimensional region increases in thickness, with decay time, though the structure in the central regions of the flow remains three dimensional. © 1995 American Institute of Physics.

Notes: The velocity, flux, and concentration distribution of solid particles in a turbulent boundary layer of a horizontal water flume were investigated experimentally by means of LDA and visualization techniques. The particles were of polystyrene (specific density ∼1.05). Results show that coherent wall structures are responsible for most of the characteristics of particle behavior throughout the boundary layer. Particles are often concentrated in regions of low velocity, associated with wall structures, and as a result the average particle velocity is lower than the fluid's. This was also noted previously by Rashidi et al., but not explained. The actual relative velocity between the particles and the surrounding fluid is often small. In addition, the data suggest that as the shear rate increases, the particle flux profiles asymptotically approach a shape where a strong gradient of particle flux exists in the lower part of the boundary layer (y+≤250), while it is relatively constant at higher elevations. This phenomenon may also be attributable to interactions with the wall structures. © 1995 American Institute of Physics.

Notes: The structural features of turbulence at the free surface of a channel flow have been experimentally investigated. The experiments were conducted in a horizontal channel of large aspect ratio in the (depth based) Reynolds number range of 2800–8800. The results indicate that the persistent structures on the free surface can be classified as upwellings, downdrafts, and spiral eddies. Upwellings are shown to be related to the bursts originating in the sheared region at the channel bottom and the eddies are seen to be generated at the edges of the upwellings. The eddies often merge if rotating in the same direction, and form "pairs" if rotating in opposite directions—though there are occasional mergers of such counter-rotating ones. The spiral eddies decay slowly and are sometimes annihilated by fresh upwellings. The population densities and the persistence times of the various structures were measured for different flow conditions. The resulting data show that the physical parameters characterizing the structures at the interface, scale with a mix of inner (wall shear stress and viscosity) and outer variables. Measurement of the streamwise and spanwise velocities at the free-surface were made by particle imaging velocimetry (PIV) and the surface normal velocity near the free-surface estimated by continuity. The results indicate that the upwellings and spiral eddy regions would be expected to dominate scalar transport rates at high Prandtl/Schmidt numbers. The one-dimensional energy spectra of the flow field at the free-surface compare well with direct numerical simulations and show a region with −5/3 slope at low wave numbers. This experimentally confirms a previous result regarding the two-dimensionality of turbulence near the free surface, based on numerical simulations by Pan and Banerjee [Phys. Fluids 7, 1649 (1995)]. © 1998 American Institute of Physics.

Notes: A method for particle image velocimetry (PIV) is presented which improves upon the accuracy, computational efficiency and dynamic range (i.e., the difference between the largest and smallest resolvable particle displacement vectors) of conventional PIV techniques. The technique is applied to free-surface turbulence to resolve energy spectra for motions with a wide dynamic range. The methodology—based on multi-grid image processing algorithms for rigid body motion analysis, estimates the displacement vectors at discrete particle locations. The essence of this technique is to estimate large scale motions from image intensity patterns of low spatial frequencies and small scale motions from intensity patterns of high spatial frequencies. Cross-correlation between a pair of time separated particle images is implemented by the hierarchical computational scheme of Burt ["Fast filter transforms for image processing," Int. J. Comput. Vision 16, 20 (1981)]. Each image is convolved with a series of band-pass filters and subsampled to obtain a set of images progressively decreasing in resolution and size. A coarse estimate of the displacement field obtained from pairs of lower resolution images are used to obtain more accurate estimates at the next (finer) level. Processing starts at the level of lowest resolution and stops at the highest resolution level, which contains the original image pair. Due to subsampling of low resolution images, the match template size can be kept constant for all stages of computation, thus eliminating the dependence of the largest resolvable displacement on the size of match template. In the present work, the search area at each level is kept constant at 3×3 pixels and the match template size at 5×5 pixels for all levels of computation. The algorithm has been implemented using simple thresholding based on the confidence level of an estimated displacement vector, as suggested by Anandan ["A computational framework and an algorithm for measurement of visual motion," Int. J. Comput. Vision, 2, 283, (1987)]. However, the confidence-level-based smoothing technique for rigid body motions (continuous velocity fields) could not be applied to displacement estimates obtained at discrete points i.e., the particle locations. Instead, smoothing was performed over the area covered by each particle. The algorithm has been tested against direct numerical simulations of turbulent flows when the flow field is known and particle images have been generated from these with the addition of noise. Both the accuracy of motion estimation and the computation time are seen to improve as compared to conventional PIV methods. Finally, video images taken of particle motion on the free-surface of a channel flow have been used to determine the capabilities of the technique in an experimental study. The resulting spectra show a quasi-two-dimensional character of the free-surface turbulent flow field, which corresponds well with the direct numerical simulations. © 1998 American Institute of Physics.

Notes: This work focuses on subgrid-scale (SGS) modeling for finite-difference large-eddy simulations, employing filters in physical space. When a filter in physical space is used, an overlap is allowed between the unresolved and the resolved scales. For such a filter, all the three terms in the classical decomposition of the SGS stress tensor are present: the Leonard and cross-terms, due to the overlap between scales, and the true SGS Reynolds tensor, expressing the pure effect of the small scales. A dynamic subgrid-scale stress model is proposed, for finite-difference large-eddy simulation of incompressible and compressible flows in which the Leonard and cross-parts of the SGS stress tensor are assumed to be proportional to the resolved part (the "modified Leonard term''), which is computed explicity. The SGS Reynolds stress is modeled by the eddy-viscosity Smagorinsky model. The two unknown parameters in this model are computed dynamically, as in Germano et al. [Phys. Fluids A 3, 1790 (1991)], but using a least squares technique. The model is tested using direct numerical simulation data for fully developed turbulent incompressible flows in presence of solid boundaries and free surfaces, and for compressible homogeneous turbulence. A "box filter'' in physical space is used. Other SGS models are also tested, viz. the dynamic model of Germano et al. (DSM), and its compressible extension by Moin et al. [Phys. Fluids A 3, 2746 (1991)], and the dynamic mixed model in Zang et al. [Phys. Fluids A 5, 3186 (1993)] (DMM) and its compressible version developed here. Results on the behavior of the different models with regard to energy exchanges and correlation with the exact SGS stresses are presented for different filter widths. In particular high correlation is found between the modified Leonard and cross-terms thus justifying the basic assumption made in the model. © 1995 American Institute of Physics.

Notes: The motion of solid particles near the wall in a turbulent boundary layer was investigated experimentally in a water flume by flow visualization techniques and by LDA. The particles were of polystyrene (specific density ∼1.05) with diameters ranging from 100 to 900 μm. Results show that particle motion, as well as entrainment and deposition processes, are controlled by the action of coherent wall structures, which appear to be funnel vortices. The behavior of the particles is consistent with the motion and effects of such vortices. The vortices appear to cause the formation of particle streaks near the wall, to create suitable conditions for particle entrainment, and to assist in particle deposition by conveying them from the outer flow to the wall region. © 1995 American Institute of Physics.

Notes: The Navier–Stokes equations have been solved, by a pseudospectral method, for pressure-driven flows between a no-slip wavy wall and a slip flat wall. Periodic boundary conditions were used in the streamwise and spanwise directions. The physical domain is mapped into a computational domain that is a rectangular parallelepiped using a nonorthogonal transformation. The pseudospectral solution procedure employed in previous studies, for example, Lam and Banerjee [Phys. Fluids A 4, 306 (1992)], eliminated the pressure and solved for the wall–normal velocity and vorticity. The other velocity components were calculated using the definition of vorticity, and the continuity equation. This procedure leads to oscillations in the pressure field when solutions were attempted in the mapped computational domain. To overcome the problem, the procedure had to be modified and the pressure solved for directly using a fractional time step technique. For the cases examined here, these modifications resulted in spectral accuracy being maintained. Flow over sinusoidal wave trains has been simulated and the results compare well with available experiments. The simulations show significant effects of the wavy boundary on the mean flow and the turbulence statistics. The mean velocity profile differs substantially from the profile for the flat-wall case, particularly in the buffer region where the fluid is under the influence of both the wavy wall and the slip boundary. The velocity fluctuations in the streamwise direction decrease in the buffer region. This effect becomes more pronounced when the wave amplitude increases. Most of the redistribution of energy, from the streamwise direction to the spanwise and wall–normal directions, occurs in a thin layer close to the boundary, downstream of the wave troughs. The energy primarily redistributes into spanwise fluctuations. High shear stress regions form downstream of the wave troughs, and streaky structures and quasi-streamwise vortices are also seen to initiate in these regions. The length of the streaks, and the extent of the quasi-streamwise vortices, scale with wave length for the two cases investigated. © 1997 American Institute of Physics.