Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by $\delt T_{Lp}$ and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of $\delt T_{Lp}$ on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading.

A Multiwavelength Study of the Intracluster Medium and the Characterization of the Multiwavelength Sub/millimeter Inductance Camera

The first part of this thesis combines Bolocam observations of the thermal Sunyaev-Zel’dovich (SZ) effect at 140 GHz with X-ray observations from Chandra, strong lensing data from the Hubble Space Telescope (HST), and weak lensing data from HST and Subaru to constrain parametric models for the distribution of dark and baryonic matter in a sample of six massive, dynamically relaxed galaxy clusters. For five of the six clusters, the full multiwavelength dataset is well described by a relatively simple model that assumes spherical symmetry, hydrostatic equilibrium, and entirely thermal pressure support. The multiwavelength analysis yields considerably better constraints on the total mass and concentration compared to analysis of any one dataset individually. The subsample of five galaxy clusters is used to place an upper limit on the fraction of pressure support in the intracluster medium (ICM) due to nonthermal processes, such as turbulent and bulk flow of the gas. We constrain the nonthermal pressure fraction at r500c to be less than 0.11 at 95% confidence, where r500c refers to radius at which the average enclosed density is 500 times the critical density of the Universe. This is in tension with state-of-the-art hydrodynamical simulations, which predict a nonthermal pressure fraction of approximately 0.25 at r500c for the clusters in this sample.

The second part of this thesis focuses on the characterization of the Multiwavelength Sub/millimeter Inductance Camera (MUSIC), a photometric imaging camera that was commissioned at the Caltech Submillimeter Observatory (CSO) in 2012. MUSIC is designed to have a 14 arcminute, diffraction-limited field of view populated with 576 spatial pixels that are simultaneously sensitive to four bands at 150, 220, 290, and 350 GHz. It is well-suited for studies of dusty star forming galaxies, galaxy clusters via the SZ Effect, and galactic star formation. MUSIC employs a number of novel detector technologies: broadband phased-arrays of slot dipole antennas for beam formation, on-chip lumped element filters for band definition, and Microwave Kinetic Inductance Detectors (MKIDs) for transduction of incoming light to electric signal. MKIDs are superconducting micro-resonators coupled to a feedline. Incoming light breaks apart Cooper pairs in the superconductor, causing a change in the quality factor and frequency of the resonator. This is read out as amplitude and phase modulation of a microwave probe signal centered on the resonant frequency. By tuning each resonator to a slightly different frequency and sending out a superposition of probe signals, hundreds of detectors can be read out on a single feedline. This natural capability for large scale, frequency domain multiplexing combined with relatively simple fabrication makes MKIDs a promising low temperature detector for future kilopixel sub/millimeter instruments. There is also considerable interest in using MKIDs for optical through near-infrared spectrophotometry due to their fast microsecond response time and modest energy resolution. In order to optimize the MKID design to obtain suitable performance for any particular application, it is critical to have a well-understood physical model for the detectors and the sources of noise to which they are susceptible. MUSIC has collected many hours of on-sky data with over 1000 MKIDs. This work studies the performance of the detectors in the context of one such physical model. Chapter 2 describes the theoretical model for the responsivity and noise of MKIDs. Chapter 3 outlines the set of measurements used to calibrate this model for the MUSIC detectors. Chapter 4 presents the resulting estimates of the spectral response, optical efficiency, and on-sky loading. The measured detector response to Uranus is compared to the calibrated model prediction in order to determine how well the model describes the propagation of signal through the full instrument. Chapter 5 examines the noise present in the detector timestreams during recent science observations. Noise due to fluctuations in atmospheric emission dominate at long timescales (less than 0.5 Hz). Fluctuations in the amplitude and phase of the microwave probe signal due to the readout electronics contribute significant 1/f and drift-type noise at shorter timescales. The atmospheric noise is removed by creating a template for the fluctuations in atmospheric emission from weighted averages of the detector timestreams. The electronics noise is removed by using probe signals centered off-resonance to construct templates for the amplitude and phase fluctuations. The algorithms that perform the atmospheric and electronic noise removal are described. After removal, we find good agreement between the observed residual noise and our expectation for intrinsic detector noise over a significant fraction of the signal bandwidth.

Large eddy simulation of turbine internal cooling ducts

Abstract Large-Eddy Simulation (LES) and hybrid Reynolds-averaged Navier–Stokes–LES (RANS–LES) methods are applied to a turbine blade ribbed internal duct with a 180° bend containing 24 pairs of ribs. Flow and heat transfer predictions are compared with experimental data and found to be in agreement. The choice of LES model is found to be of minor importance as the flow is dominated by large geometric scale structures. This is in contrast to several linear and nonlinear RANS models, which display turbulence model sensitivity. For LES, the influence of inlet turbulence is also tested and has a minor impact due to the strong turbulence generated by the ribs. Large scale turbulent motions destroy any classical boundary layer reducing near wall grid requirements. The wake-type flow structure makes this and similar flows nearly Reynolds number independent, allowing a range of flows to be studied at similar cost. Hence LES is a relatively cheap method for obtaining accurate heat transfer predictions in these types of flows.

Columnar structures and stress relaxation in thick GaN films grown on sapphire by HVPE

Thick GaN films with high quality are directly grown on sapphire in a home-built vertical hydride vapour phase epitaxy (HVPE) reactor. The optical and structural properties of large scale columnar do-mains near the interface are studied using cathodoluminescence and micro-Raman scattering. These columnar do-mains Show a strong emission intensity due to extremely high free carrier concentration up to 2 x 10(19) cm(-3), which are related with impurities trapped in structural defects. The compressive stress in GaN Elm clearly decreases with increasing distance from interface. The quasi-continuous columnar domains play an important role in the stress relaxation for the upper high quality layer.

Temperature dependence of the formation of nano-scale indium clusters in InAlGaN alloys on Si(111) substrates

The temperature dependence of the formation of nano-scale indium clusters in InAlGaN quaternary alloys, which are grown by metalorganic chemical vapour deposition on GaN/Si(111) epilayers, is investigated. Firm evidence is provided to support the existence of phase separation, or nano-scale In-rich clusters, by the combined results of high-resolution transmission electron microscopy (HRTEM), high-resolution x-ray diffraction (HRXRD) and micro-Raman spectra. The results of HRXRD and Raman spectra indicate that the degree of phase separation is strong and the number of In clusters in the InAlGaN layers on silicon substrate is higher at lower growth temperatures than that at higher growth temperatures, which limits the In and Al incorporated into the InAlGaN quaternary alloys. The detailed mechanism of luminescence in this system is studied by low temperature photoluminescence (LT-PL). We conclude that the ultraviolet (UV) emission observed in the quaternary InAlGaN alloys arises from the matrix of a random alloy, and the second emission peak in the blue-green region results from the nano-scale indium clusters.

Constrain intergalactic medium from the SZ effect map

In this paper, we try to detect the SZ effect in the 2MASS DWT clusters and less bound objects in order to constrain the warm-hot intergalactic medium distribution on large scales by cross-correlation analysis. The results of both observed WMAP and mock SZ effect map indicate that the hot gas distributes from inside as well as outside of the high density regions of galaxy clusters, which is consistent with the results of both observation and hydro simulation. Therefore, the DWT measurement of the cross-correlation would be a powerful tool to probe the missing of baryons in the universe.

Cross-correlation between WMAP and 2MASS: non-Gaussianity induced by the SZ effect

We study the non-Gaussianity induced by the Sunyaev-Zel'dovich (SZ) effect in cosmic microwave background (CMB) fluctuation maps. If a CMB map is contaminated by the SZ effect of galaxies or galaxy clusters, the CMB maps should have similar non-Gaussian features to the galaxy and cluster fields. Using the WMAP data and 2MASS galaxy catalogue, we show that the non-Gaussianity of the 2MASS galaxies is imprinted on WMAP maps. The signature of non-Gaussianity can be seen with the fourth-order cross-correlation between the wavelet variables of the WMAP maps and 2MASS clusters. The intensity of the fourth-order non-Gaussian features is found to be consistent with the contamination of the SZ effect of 2MASS galaxies. We also show that this non-Gaussianity can not be seen by the high-order autocorrelation of the WMAP. This is because the SZ signals in the autocorrelations of the WMAP data generally are weaker than the WMAP-2MASS cross-correlations by a factor f(2), which is the ratio between the powers of the SZ-effect map and the CMB fluctuations on the scale considered. Therefore, the ratio of high-order autocorrelations of CMB maps to cross-correlations of the CMB maps and galaxy field would be effective to constrain the powers of the SZ effect on various scales.

Atomic Structures and Radiative Properties of He-like Ions in Debye Plasmas

A detailed investigation of plasma screening effects on atomic structure and transition properties are presented for He-like ions embedded in dense plasma environment. Multi-configuration Dirac-Fock calculations were carried out for these ions by considering a Debye-Huckel potential. A large-scale relativistic configuration-interaction method is adopted to calculate transition energies and transition probabilities and to allow for a systematic improvement of the calculations. Comparison of the presently calculated results with others, when available, is made.

电离层行进式扰动的GPS台网监测研究

In this dissertation, we investigated two types of traveling ionospheric disturbances (TIDs)/gravity waves (GWs) triggered separately by auroral energy input during super geomagnetic storms and solar terminator (ST) under quiet geomagnetic conditions (kp<3+) using TEC measurements from the global network of GPS receivers. Research into the generation and propagation of TIDs/GWs during storms greatly enhance our understandings on the evolution processes of energy transportation from the high-latitude’s magnetosphere to the low-latitude ionosphere and the conjugated effect of TIDs propagation between the northern and southern hemispheres. Our results revealed that the conjugacy of propagation direction between the northern and southern hemispheres was subject to the influence of Coriolis force. We also figure out the evolution processes of ionospheric disturbances at the global scale. These are important topics that had not been well addressed previously. In addition, we also obtained thee wave structures of medium scale TIDs excited by the solar terminator (ST) moving over the northern America and physical mechanisms involved. Our observations confirm that the ST is a stable and repetitive source of ionospheric wave disturbances and the evidence of solar terminator generated disturbances has been demonstrated experimentally via the GPS TEC measurement. The main researches and results of this dissertation are as follows. First, the global traveling ionospheric disturbances (TIDs) during the drastic magnetic storms of October 29–31, 2003 were analyzed using the Global Position System (GPS) total electron content (TEC) data observed in the Asian-Australian, European and North American sectors. We collected the most comprehensive set of the TEC data from more than 900 GPS stations on the International GNSS Services (IGS) website and introduce here a strategy that combines polynomial fitting and multi-channel maximum entropy spectral analysis to obtain TID parameters. Moreover, in collaboration with my thesis advisor, I have developed an imaging technique of 2-dimensional map of TIDs structures to obtain spatial and temporal maps of large scale traveling ionospheric disturbances (LSTIDs). The clear structures of TEC perturbations map during the passage of TIDs were displayed. The results of our study are summarized as follows: (1) Large-scale TIDs (LSTIDs) and medium-scale TIDs (MSTIDs) were detected in all three sectors after the sudden commencement (SC) of the magnetic storm, and their features showed longitudinal and latitudinal dependences. The duration of TIDs was longer at higher latitudes than at middle latitudes, with a maximum of about 16 h. The TEC variation amplitude of LSTIDs was larger in the North American sector than in the two other sectors. At the lower latitudes, the ionospheric perturbations were more complicated, and their duration and amplitude were relatively longer and larger. (2) The periods and phase speeds of TIDs were different in these three sectors. In Europe, the TIDs propagated southward; in North America and Asia, the TIDs propagated southwestward; in the near-equator region, the disturbances propagated with the azimuth (the angle of the propagation direction of the LSTIDs measured clockwise from due north with 0°) of 210° showing the influence of Coriolis force; in the Southern Hemisphere, the LSTIDs propagated conjugatedly northwestward. Both the southwestward and northeastward propagating LSTIDs are found in the equatorial region. These results mean that the Coriolis effect cannot be ignored for the wave propagation of LSTIDs and that the propagation direction is correlated with the polar magnetic activity. (3) The day (day of year: 301) before the SC (sudden commencement) of magnetic storm, we observed a sudden TEC skip disturbances (±10 TECU). It should be a response for the high flux of proton during the solar flare event, but not the magnetic storms. Next, the most comprehensive and dense GPS network’s data from North-America region were used in this paper to analyze the medium scale traveling ionospheric disturbances (MSTIDs) which were generated by the moving solar terminator during the quiet days in 2005. We applied the multi-channel maximum entropy spectral analysis to calculated TID parameters, and found that the occurrence of ST-MSTIDs depends on the seasonal variations. The results of our study are summarized as follows: (1) MSTIDs stimulated by the moving ST (ST-MSTIDs) are detected at mid-latitudes after the passage of the solar terminator with the life time of 2～3 hours and the variation amplitude of 0.2～0.8 TECU. Spectral analysis indicated that the horizontal wavelength, average period, horizontal phase velocity of the MSTIDs are around 300±150 km，150±80 m/s and 25±15 min, respectively. In addition, ST-MSTIDs have wave fronts elongating the moving ST direction and almost parallel to ST. (2) The statistical results demonstrate that the dusk MSTIDs stimulated by ST is more obvious than the dawn MSTIDs in summer. On the contrary, the more-pronounced dawn MSTIDs occurs in winter. (3) Further analysis indicates that the seasonal variations of ST-MSTIDs occurrence frequency are most probably related to the seasonal differences of the variations of EUV flux in the ionosphere region and recombination process during sunrise and sunset period at mid-latitudes. Statistical study of occurrence characteristics of TIDs using the GPS network in North-American and European during solar maximum, In conclusion, statistical studies of the propagation characteristics of TIDs, which excited by the two common origins including geomagnetic storms and moving solar terminator, were involved with global GPS TEC databasein this thesis. We employed the multichannel maximum entropy spectral analysis method to diagnose the characteristics of propagation and evolvement of ionospheric disturbances, also, the characteristics of their regional distribution and climatological variations were revealed by the statistic analysis. The results of these studies can improve our knowledge about the energy transfer in the solar-terrestrial system and the coupling process between upper and lower atmosphere (thermosphere-ionosphere-mesosphere). On the other hand, our results of the investigation on TIDs generated by particular linear origin such as ST are important for developing ionospheric irregularity physics and modeling the transionosphere radio wave propagation. Besides, the GPS TEC representation of the ST-generated ionospheric structure suggests a better possibility for investigating this phenomenon. Subsequently, there are scientific meaning of the result of this dissertation to deeply discuss the energy transfer and coupling in the ionosphere, as well as realistic value to space weather forecast in the ionosphere region.

Columnar-to-disk structural transition in nanoscale (SiO2)(N) clusters

Bromley, S. T., Flikkema, E. (2005) 'Columnar-to-disk structural transition in nanoscale (SiO2)(N) clusters'. Physical Review Letters 95 (18), article no: 185505 RAE2008

The dayside high-latitude trough under quiet geomagnetic conditions: Radio tomography and the CTIP model

Pryse, Sian; Dewis, K.L.; Middleton, H.R.; Balthazor, R.L., (2005) 'The dayside high-latitude trough under quiet geomagnetic conditions: Radio tomography and the CTIP model', Annales Geophysicae 23(4) pp.1199-1206 RAE2008

Large molecules in the envelope surrounding IRC+10216

A new chemical model of the circumstellar envelope surrounding the carbon-rich star IRC+10216 is developed that includes carbon-containing molecules with up to 23 carbon atoms. The model consists of 3851 reactions involving 407 gas-phase species. Sizeable abundances of a variety of large molecules - including carbon clusters, unsaturated hydrocarbons and cyanopolyynes - have been calculated. Negative molecular ions of chemical formulae C-n(-) and CnH- (7 less than or equal to n less than or equal to 23) exist in considerable abundance, with peak concentrations at distances from the central star somewhat greater than their neutral counterparts. The negative ions might be detected in radio emission, or even in the optical absorption of background field stars. The calculated radial distributions of the carbon-chain CnH radicals are looked at carefully and compared with interferometric observations.

Self-organized electromagnetic field structures in laser-produced counter-streaming plasmas

Self-organization(1,2) occurs in plasmas when energy progressively transfers from smaller to larger scales in an inverse cascade(3). Global structures that emerge from turbulent plasmas can be found in the laboratory(4) and in astrophysical settings; for example, the cosmic magnetic field(5,6,) collisionless shocks in supernova remnants(7) and the internal structures of newly formed stars known as Herbig-Haro objects(8). Here we show that large, stable electromagnetic field structures can also arise within counter-streaming supersonic plasmas in the laboratory. These surprising structures, formed by a yet unexplained mechanism, are predominantly oriented transverse to the primary flow direction, extend for much larger distances than the intrinsic plasma spatial scales and persist for much longer than the plasma kinetic timescales. Our results challenge existing models of counter-streaming plasmas and can be used to better understand large-scale and long-time plasma self-organization.

Behavioural context of multi-scale species distribution models assessed by radio-tracking

Incorporating ecological processes and animal behaviour into Species Distribution Models (SDMs) is difficult. In species with a central resting or breeding place, there can be conflict between the environmental requirements of the 'central place' and foraging habitat. We apply a multi-scale SDM to examine habitat trade-offs between the central place, roost sites, and foraging habitat in . Myotis nattereri. We validate these derived associations using habitat selection from behavioural observations of radio-tracked bats. A Generalised Linear Model (GLM) of roost occurrence using land cover variables with mixed spatial scales indicated roost occurrence was positively associated with woodland on a fine scale and pasture on a broad scale. Habitat selection of radio-tracked bats mirrored the SDM with bats selecting for woodland in the immediate vicinity of individual roosts but avoiding this habitat in foraging areas, whilst pasture was significantly positively selected for in foraging areas. Using habitat selection derived from radio-tracking enables a multi-scale SDM to be interpreted in a behavioural context. We suggest that the multi-scale SDM of . M. nattereri describes a trade-off between the central place and foraging habitat. Multi-scale methods provide a greater understanding of the ecological processes which determine where species occur and allow integration of behavioural processes into SDMs. The findings have implications when assessing the resource use of a species at a single point in time. Doing so could lead to misinterpretation of habitat requirements as these can change within a short time period depending on specific behaviour, particularly if detectability changes depending on behaviour. © 2011 Gesellschaft für ökologie.

Developed turbulence and nonlinear amplification of magnetic fields in laboratory and astrophysical plasmas

The visible matter in the universe is turbulent and magnetized. Turbulence in galaxy clusters is produced by mergers and by jets of the central galaxies and believed responsible for the amplification of magnetic fields. We report on experiments looking at the collision of two laser-produced plasma clouds, mimicking, in the laboratory, a cluster merger event. By measuring the spectrum of the density fluctuations, we infer developed, Kolmogorov-like turbulence. From spectral line broadening, we estimate a level of turbulence consistent with turbulent heating balancing radiative cooling, as it likely does in galaxy clusters. We show that the magnetic field is amplified by turbulent motions, reaching a nonlinear regime that is a precursor to turbulent dynamo. Thus, our experiment provides a promising platform for understanding the structure of turbulence and the amplification of magnetic fields in the universe.