Two-scale computational modelling of water flow in unsaturated soils containing irregular-shaped inclusions

The focus of this paper is two-dimensional computational modelling of water flow in unsaturated soils consisting of weakly conductive disconnected inclusions embedded in a highly conductive connected matrix. When the inclusions are small, a two-scale Richards’ equation-based model has been proposed in the literature taking the form of an equation with effective parameters governing the macroscopic flow coupled with a microscopic equation, defined at each point in the macroscopic domain, governing the flow in the inclusions. This paper is devoted to a number of advances in the numerical implementation of this model. Namely, by treating the micro-scale as a two-dimensional problem, our solution approach based on a control volume finite element method can be applied to irregular inclusion geometries, and, if necessary, modified to account for additional phenomena (e.g. imposing the macroscopic gradient on the micro-scale via a linear approximation of the macroscopic variable along the microscopic boundary). This is achieved with the help of an exponential integrator for advancing the solution in time. This time integration method completely avoids generation of the Jacobian matrix of the system and hence eases the computation when solving the two-scale model in a completely coupled manner. Numerical simulations are presented for a two-dimensional infiltration problem.

Intrinsic equations of steady rotating, incompressible viscous fluid flow

The equations governing the flow of a steady rotating incompressible viscous fluid are expressed in intrinsic form along the vortex lines and their normals. Using these equations the effects of rotation on the geometric properties of viscous fluid flows are studied. A particular flow in which the vortex lines are right circular helices is discussed.

Velocity ratio-cum-transfer matrix method for the evaluation of a muffler with mean flow

The paper deals with a method for the evaluation of exhaust muffers with mean flow. A new set of variables, convective pressure and convective mass velocity, have been defined to replace the acoustic variables. An expression for attenuation (insertion loss) of a muffler has been proposed in terms of convective terminal impedances and a velocity ratio, on the lines of the one existing for acoustic filters. In order to evaluate the velocity ratio in terms of convective variables, transfer matrices for various muffler elements have been derived from the basic relations of energy, mass and momentum. Finally, the velocity ratiocum-transfer matrix method is illustrated for a typical straight-through muffler.

Non-linear couette flow with heat transfer using the B-G-K model

In recent years a large number of investigators have devoted their efforts to the study of flow and heat transfer in rarefied gases, using the BGK [1] model or the Boltzmann kinetic equation. The velocity moment method which is based on an expansion of the distribution function as a series of orthogonal polynomials in velocity space, has been applied to the linearized problem of shear flow and heat transfer by Mott-Smith [2] and Wang Chang and Uhlenbeck [3]. Gross, Jackson and Ziering [4] have improved greatly upon this technique by expressing the distribution function in terms of half-range functions and it is this feature which leads to the rapid convergence of the method. The full-range moments method [4] has been modified by Bhatnagar [5] and then applied to plane Couette flow using the B-G-K model. Bhatnagar and Srivastava [6] have also studied the heat transfer in plane Couette flow using the linearized B-G-K equation. On the other hand, the half-range moments method has been applied by Gross and Ziering [7] to heat transfer between parallel plates using Boltzmann equation for hard sphere molecules and by Ziering [83 to shear and heat flow using Maxwell molecular model. Along different lines, a moment method has been applied by Lees and Liu [9] to heat transfer in Couette flow using Maxwell's transfer equation rather than the Boltzmann equation for distribution function. An iteration method has been developed by Willis [10] to apply it to non-linear heat transfer problems using the B-G-K model, with the zeroth iteration being taken as the solution of the collisionless kinetic equation. Krook [11] has also used the moment method to formulate the equivalent continuum equations and has pointed out that if the effects of molecular collisions are described by the B-G-K model, exact numerical solutions of many rarefied gas-dynamic problems can be obtained. Recently, these numerical solutions have been obtained by Anderson [12] for the non-linear heat transfer in Couette flow,

Very large array detection of radio recombination lines from the radio nucleus of NGC 253: Ionization by a weak active galactic nucleus, an obscured super star cluster, or a compact supernova remnant?

We have imaged the H92alpha and H75alpha radio recombination line (RRL) emissions from the starburst galaxy NGC 253 with a resolution of similar to4 pc. The peak of the RRL emission at both frequencies coincides with the unresolved radio nucleus. Both lines observed toward the nucleus are extremely wide, with FWHMs of similar to200 km s(-1). Modeling the RRL and radio continuum data for the radio nucleus shows that the lines arise in gas whose density is similar to10(4) cm(-3) and mass is a few thousand M., which requires an ionizing flux of (6-20) x 10(51) photons s(-1). We consider a supernova remnant (SNR) expanding in a dense medium, a star cluster, and also an active galactic nucleus (AGN) as potential ionizing sources. Based on dynamical arguments, we rule out an SNR as a viable ionizing source. A star cluster model is considered, and the dynamics of the ionized gas in a stellar-wind driven structure are investigated. Such a model is only consistent with the properties of the ionized gas for a cluster younger than similar to10(5) yr. The existence of such a young cluster at the nucleus seems improbable. The third model assumes the ionizing source to be an AGN at the nucleus. In this model, it is shown that the observed X-ray flux is too weak to account for the required ionizing photon flux. However, the ionization requirement can be explained if the accretion disk is assumed to have a big blue bump in its spectrum. Hence, we favor an AGN at the nucleus as the source responsible for ionizing the observed RRLs. A hybrid model consisting of an inner advection-dominated accretion flow disk and an outer thin disk is suggested, which could explain the radio, UV, and X-ray luminosities of the nucleus.

Parallel Flow-Sensitive Pointer Analysis by Graph-Rewriting

Precise pointer analysis is a problem of interest to both the compiler and the program verification community. Flow-sensitivity is an important dimension of pointer analysis that affects the precision of the final result computed. Scaling flow-sensitive pointer analysis to millions of lines of code is a major challenge. Recently, staged flow-sensitive pointer analysis has been proposed, which exploits a sparse representation of program code created by staged analysis. In this paper we formulate the staged flow-sensitive pointer analysis as a graph-rewriting problem. Graph-rewriting has already been used for flow-insensitive analysis. However, formulating flow-sensitive pointer analysis as a graph-rewriting problem adds additional challenges due to the nature of flow-sensitivity. We implement our parallel algorithm using Intel Threading Building Blocks and demonstrate considerable scaling (upto 2.6x) for 8 threads on a set of 10 benchmarks. Compared to the sequential implementation of staged flow-sensitive analysis, a single threaded execution of our implementation performs better in 8 of the benchmarks.

Framework for morphometric classification of cells in imaging flow cytometry

Imaging flow cytometry is an emerging technology that combines the statistical power of flow cytometry with spatial and quantitative morphology of digital microscopy. It allows high-throughput imaging of cells with good spatial resolution, while they are in flow. This paper proposes a general framework for the processing/classification of cells imaged using imaging flow cytometer. Each cell is localized by finding an accurate cell contour. Then, features reflecting cell size, circularity and complexity are extracted for the classification using SVM. Unlike the conventional iterative, semi-automatic segmentation algorithms such as active contour, we propose a noniterative, fully automatic graph-based cell localization. In order to evaluate the performance of the proposed framework, we have successfully classified unstained label-free leukaemia cell-lines MOLT, K562 and HL60 from video streams captured using custom fabricated cost-effective microfluidics-based imaging flow cytometer. The proposed system is a significant development in the direction of building a cost-effective cell analysis platform that would facilitate affordable mass screening camps looking cellular morphology for disease diagnosis. Lay description In this article, we propose a novel framework for processing the raw data generated using microfluidics based imaging flow cytometers. Microfluidics microscopy or microfluidics based imaging flow cytometry (mIFC) is a recent microscopy paradigm, that combines the statistical power of flow cytometry with spatial and quantitative morphology of digital microscopy, which allows us imaging cells while they are in flow. In comparison to the conventional slide-based imaging systems, mIFC is a nascent technology enabling high throughput imaging of cells and is yet to take the form of a clinical diagnostic tool. The proposed framework process the raw data generated by the mIFC systems. The framework incorporates several steps: beginning from pre-processing of the raw video frames to enhance the contents of the cell, localising the cell by a novel, fully automatic, non-iterative graph based algorithm, extraction of different quantitative morphological parameters and subsequent classification of cells. In order to evaluate the performance of the proposed framework, we have successfully classified unstained label-free leukaemia cell-lines MOLT, K562 and HL60 from video streams captured using cost-effective microfluidics based imaging flow cytometer. The cell lines of HL60, K562 and MOLT were obtained from ATCC (American Type Culture Collection) and are separately cultured in the lab. Thus, each culture contains cells from its own category alone and thereby provides the ground truth. Each cell is localised by finding a closed cell contour by defining a directed, weighted graph from the Canny edge images of the cell such that the closed contour lies along the shortest weighted path surrounding the centroid of the cell from a starting point on a good curve segment to an immediate endpoint. Once the cell is localised, morphological features reflecting size, shape and complexity of the cells are extracted and used to develop a support vector machine based classification system. We could classify the cell-lines with good accuracy and the results were quite consistent across different cross validation experiments. We hope that imaging flow cytometers equipped with the proposed framework for image processing would enable cost-effective, automated and reliable disease screening in over-loaded facilities, which cannot afford to hire skilled personnel in large numbers. Such platforms would potentially facilitate screening camps in low income group countries; thereby transforming the current health care paradigms by enabling rapid, automated diagnosis for diseases like cancer.

Structure of Vortex Dislocations in a Wake-type Flow

Structure and dynamical processes of vortex dislocations in a kind of wake-type flow are described clearly by vortex lines, which are directly constructed from data of three-dimensional direct numerical simulations of the flow evolution.

Experimental Study the Flow Pattern in Thermal Capillary Convection

In the present research work, the thermal capillary convection has been investigated and measured by particle image velocimetry (PIV) technique. There is one liquid layer in a rectangular cavity with different temperature’s sidewalls. The cavity is 52mm,42mm,20mm, 4mm in height of the silicon oil liquid layer. A sidewall of the cavity is heated by electro-thermal film, another sidewall is cooled by the semiconductor cooling sheet. The velocity field and the stream lines in cross section in liquid layer have been obtained at different temperature difference. The present experiment demonstrates that the pattern of the convection mainly relates with temperature difference.

Characteristics of flow fields induced by interfacial waves in two-layer fluid

When designing deep ocean structures, it is necessary to estimate the effects of internal waves on the platform and auxiliary parts such as tension leg, riser and mooring lines. Up to now, only a few studies are concerned with the internal wave velocity fields. By using the most representative two-layer model, we have analyzed the behavior of velocity field induced by interfacial wave in the present paper. We find that there may exist velocity shear of fluid particles in the upper and lower layers so that any structures in the ocean are subjected to shear force nearby the interface. In the meantime, the magnitude of velocity for long internal wave appears spatially uniform in the respective layer although they still decay exponentially. Finally, the temporal variation for Stokes and solitary waves are shown to be of periodical and pulse type.

Vortex dislocations in wake-type flow induced by spanwise disturbances

Vortex dislocations in wake-type flow induced by three types of spanwise disturbances superimposed on an upstream velocity profile are investigated by direct numerical simulations. Three distinct modes of vortex dislocations and flow transitions have been found. A local spanwise exponential decay disturbance leads to the appearance of a twisted chainlike mode of vortex dislocation. A stepped spanwise disturbance causes a streamwise periodic spotlike mode of vortex dislocation. A spanwise sinusoidal wavy disturbance with a moderate waviness causes a strong unsteadiness of wake behavior. This unsteadiness starts with a systematic periodic mode of vortex dislocation in the spanwise direction followed by the spanwise vortex shedding suppressed completely with increased time and the near wake becoming a steady shear flow. Characteristics of these modes of vortex dislocation and complex vortex linkages over the dislocation, as well as the corresponding dynamic processes related to the appearance of dislocations, are described by examining the variations of vortex lines and vorticity distribution. The nature of the vortex dislocation is demonstrated by the substantial vorticity modification of the spanwise vortex from the original spanwise direction to streamwise and vertical directions, accompanied by the appearance of noticeable vortex branching and complex vortex linking, all of which are produced at the locations with the biggest phase difference or with a frequency discontinuity between shedding cells. The effect of vortex dislocation on flow transition, either to an unsteady irregular vortex flow or suppression of the Kaacutermaacuten vortex shedding making the wake flow steady state, is analyzed. Distinct similarities are found in the mechanism and main flow phenomena between the present numerical results obtained in wake-type flows and the experimental-numerical results of cylinder wakes reported in previous studies.

Effects of exhaust inlet size on the fluid flow patterns created by an Aaberg exhaust hood

The aim of this paper is to investigate the effects of the size of the exhaust inlet on the fluid flow patterns created by an Aaberg exhaust hood. Streamlines and lines of constant speed modelling those created by an Aaberg hood are used to examine the effect of the inlet size on the size and profile of the hood's effective suction area. © 1992.

Wave dispersion and power flow analysis of an infinite aerofoil-shaped beam

Classic flutter analysis models an aerofoil as a two degree-of-freedom rigid body supported by linear and torsional springs, which represent the bending and torsional stiffness of the aerofoil section. In this classic flutter model, no energy transfer or dissipation can occur in the span-wise direction of the aerofoil section. However, as the aspect ratio of an aerofoil section increases, this span-wise energy transfer - in the form of travelling waves - becomes important to the overall system dynamics. This paper extends the classic flutter model to include travelling waves in the span-wise direction. Namely, wave dispersion and power flow analysis of an infinite, aerofoil-shaped beam, subject to bending, torsion, tension and a constant wind excitation, is used to investigate the overall system stability. Examples of potential applications for these high aspect ratio aerofoil sections include high-altitude balloon tethers, towed cables, offshore risers and mooring lines.

Effect of 80.55 MeV/u(12)C(6+) ions on radiosensitivity and cell cycle of human hepatoma cell lines

In this paper, the relationship between radiosensitivity, cell cycle alteration and the change of apoptosis in different human hepatoma cell lines irradiated by heavy ions were studied with the aim of building up the base data for clinical therapy. Exponentially growing hepatoma cell lines were irradiated by 80.55 MeV/u(12)C(6+) ions at a dose of 0 Gy, 0.5 Gy, 1 Gy, 2 Gy, 4 Gy and 8 Gy. The radiosensitivity was assessed by means of the colony-forming assay. The DNA content, the percentage of each cell-cycle phase and the apoptosis rate were obtained with flow cytometry methods. After the irradiation, the SF2 (survival fraction at 2 gray) of SMMC-7721 cells were evidently lower than that of HepG2 cells. The S phase arrest, G2/M phase arrest delay and the apoptosis in the two hepatoma cell lines varied with the increase of the dose and repair time. The heavy ions could obviously kill the human hepatoma cell lines. Compared to HepG2 cells, SMMC-7721 cells were more radiosensitive to C-12(6+) ions.

Molecular characterization of chordoma xenografts generated from a novel primary chordoma cell source and two chordoma cell lines.

OBJECT: Chordoma cells can generate solid-like tumors in xenograft models that express some molecular characteristics of the parent tumor, including positivity for brachyury and cytokeratins. However, there is a dearth of molecular markers that relate to chordoma tumor growth, as well as the cell lines needed to advance treatment. The objective in this study was to isolate a novel primary chordoma cell source and analyze the characteristics of tumor growth in a mouse xenograft model for comparison with the established U-CH1 and U-CH2b cell lines. METHODS: Primary cells from a sacral chordoma, called "DVC-4," were cultured alongside U-CH1 and U-CH2b cells for more than 20 passages and characterized for expression of CD24 and brachyury. While brachyury is believed essential for driving tumor formation, CD24 is associated with healthy nucleus pulposus cells. Each cell type was subcutaneously implanted in NOD/SCID/IL2Rγ(null) mice. The percentage of solid tumors formed, time to maximum tumor size, and immunostaining scores for CD24 and brachyury (intensity scores of 0-3, heterogeneity scores of 0-1) were reported and evaluated to test differences across groups. RESULTS: The DVC-4 cells retained chordoma-like morphology in culture and exhibited CD24 and brachyury expression profiles in vitro that were similar to those for U-CH1 and U-CH2b. Both U-CH1 and DVC-4 cells grew tumors at rates that were faster than those for U-CH2b cells. Gross tumor developed at nearly every site (95%) injected with U-CH1 and at most sites (75%) injected with DVC-4. In contrast, U-CH2b cells produced grossly visible tumors in less than 50% of injected sites. Brachyury staining was similar among tumors derived from all 3 cell types and was intensely positive (scores of 2-3) in a majority of tissue sections. In contrast, differences in the pattern and intensity of staining for CD24 were noted among the 3 types of cell-derived tumors (p < 0.05, chi-square test), with evidence of intense and uniform staining in a majority of U-CH1 tumor sections (score of 3) and more than half of the DVC-4 tumor sections (scores of 2-3). In contrast, a majority of sections from U-CH2b cells stained modestly for CD24 (scores of 1-2) with a predominantly heterogeneous staining pattern. CONCLUSIONS: This is the first report on xenografts generated from U-CH2b cells in which a low tumorigenicity was discovered despite evidence of chordoma-like characteristics in vitro. For tumors derived from a primary chordoma cell and U-CH1 cell line, similarly intense staining for CD24 was observed, which may correspond to their similar potential to grow tumors. In contrast, U-CH2b tumors stained less intensely for CD24. These results emphasize that many markers, including CD24, may be useful in distinguishing among chordoma cell types and their tumorigenicity in vivo.