Scheme Construction with Numerical Flux Residual Correction (NFRC) and Group Velocity Control (GVC)

For simulating multi-scale complex flow fields like turbulent flows, the high order accurate schemes are preferred. In this paper, a scheme construction with numerical flux residual correction (NFRC) is presented. Any order accurate difference approximation can be obtained with the NFRC. To improve the resolution of the shock, the constructed schemes are modified with group velocity control (GVC) and weighted group velocity control (WGVC). The method of scheme construction is simple, and it is used to solve practical problems.

A Constrained Particle Dynamics for Continuum-Particle Hybrid Method in Micro- and Nano-Fluidics

A hybrid method of continuum and particle dynamics is developed for micro- and nano-fluidics, where fluids are described by a molecular dynamics (MD) in one domain and by the Navier-Stokes (NS) equations in another domain. In order to ensure the continuity of momentum flux, the continuum and molecular dynamics in the overlap domain are coupled through a constrained particle dynamics. The constrained particle dynamics is constructed with a virtual damping force and a virtual added mass force. The sudden-start Couette flows with either non-slip or slip boundary condition are used to test the hybrid method. It is shown that the results obtained are quantitatively in agreement with the analytical solutions under the non-slip boundary conditions and the full MD simulations under the slip boundary conditions.

Rapidly Infrared-Assisted Cooperatively Self-Assembled Highly Ordered Multiscale Porous Materials

In this paper, cooperative self-assembly (CSA) of colloidal spheres with different sizes was studied. It was found that a complicated jamming effect makes it difficult to achieve an optimal self-assembling condition for construction of a well-ordered stacking of colloidal spheres in a relatively short growth time by CSA. Through the use of a characteristic infrared (IR) technique to significantly accelerate local evaporation on the growing interface without changing the bulk growing environment, a concise three-parameter (temperature, pressure, and IR intensity) CSA method to effectively overcome the jamming effect has been developed. Mono- and multiscale inverse opals in a large range of lattice scales can be prepared within a growth time (15-30 min) that is remarkably shorter than the growth times of several hours for previous methods. Scanning electron microscopy images and transmittance spectra demonstrated the superior crystalline and optical qualities of the resulting materials. More importantly, the new method enables optimal conditions for CSA without limitations on sizes and materials of multiple colloids. This strategy not only makes a meaningful advance in the applicability and universality of colloidal crystals and ordered porous materials but also can be an inspiration to the self-assembly systems widely used in many other fields, such as nanotechnology and molecular bioengineering.

A Cip/Multi-Moment Finite Volume Method For Shallow Water Equations With Source Terms

A novel finite volume method has been presented to solve the shallow water equations. In addition to the volume-integrated average (VIA) for each mesh cell, the surface-integrated average (SIA) is also treated as the model variable and is independently predicted. The numerical reconstruction is conducted based on both the VIA and the SIA. Different approaches are used to update VIA and SIA separately. The SIA is updated by a semi-Lagrangian scheme in terms of the Riemann invariants of the shallow water equations, while the VIA is computed by a flux-based finite volume formulation and is thus exactly conserved. Numerical oscillation can be effectively avoided through the use of a non-oscillatory interpolation function. The numerical formulations for both SIA and VIA moments maintain exactly the balance between the fluxes and the source terms. 1D and 2D numerical formulations are validated with numerical experiments. Copyright (c) 2007 John Wiley & Sons, Ltd.

Equilibrium of isolated magnetostatic flux structure

Two-dimensional magnetostatic models of flux structure confined in stratified atmosphere are discussed in the present paper. The magnetic field in the flux structure is assumed to be force-free at the first step. Numerical solutions for this nonlinear free boundary problem are obtained by finite element method. Results show clearly the relation between the inside fields and outside pressure, especially the influence of atmospheric pressure distribution on the flux structure.

A model of an isolated magnetic-flux tube in a stratified atmosphere

A two-dimensional model of a magnetic flux tube confined in a gravitational stratified atmosphere is discussed. The magnetic field in the flux tube is assumed to be force-free. By using the approximation of large scale height, the problem of a free boundary with nonlinear conditions may be reduced to one involving a fixed boundary. The two-dimensional features are obtained by applying the perturbation method and adopting the Luest-Schlueter model as the basic state. The results show that the configuration of a flux tube confined in a gravitational stratified atmosphere is divergent, and the more twisted the magnetic field, the more divergent is the flux tube.

High Order Multi-Moment Constrained Finite Volume Method. Part I: Basic Formulation

A new high-order finite volume method based on local reconstruction is presented in this paper. The method, so-called the multi-moment constrained finite volume (MCV) method, uses the point values defined within single cell at equally spaced points as the model variables (or unknowns). The time evolution equations used to update the unknowns are derived from a set of constraint conditions imposed on multi kinds of moments, i.e. the cell-averaged value and the point-wise value of the state variable and its derivatives. The finite volume constraint on the cell-average guarantees the numerical conservativeness of the method. Most constraint conditions are imposed on the cell boundaries, where the numerical flux and its derivatives are solved as general Riemann problems. A multi-moment constrained Lagrange interpolation reconstruction for the demanded order of accuracy is constructed over single cell and converts the evolution equations of the moments to those of the unknowns. The presented method provides a general framework to construct efficient schemes of high orders. The basic formulations for hyperbolic conservation laws in 1- and 2D structured grids are detailed with the numerical results of widely used benchmark tests. (C) 2009 Elsevier Inc. All rights reserved.

Angular distributions of multiphoton detachment of H- in various infrared laser fields

Using a nonperturbative quantum scattering theory, the photoelectron angular distributions (PADs) from the multiphoton detachment of H- ions in strong, linearly polarized infrared laser fields are obtained to interpret recent experimental observations. In our theoretical treatment, the PADs in n-photon detachment are determined by the nth-order generalized phased Bessel (GPB) functions X-n(Z(f),eta). The advantage of using the GPB scenario to calculate PADs is its simplicity: a single special function (GPB) without any mixing coefficient can express PADs observed by recent experiments. Thus, the GPB scenario can be called a parameterless scenario.

Ultraviolet-infrared femtosecond laser-induced damage in fused silica and CaF2 crystals

The damage in fused silica and CaF2 crystals induced by wavelength tunable femtosecond lasers is studied. The threshold fluence is observed to increase rapidly with laser wavelength lambda in the region of 250-800 nm, while it is nearly a constant for 800 method. The reflectivity increases rapidly in the latter half of pump pulse, which supports that impact ionization plays an important role in the generation of conduction band electrons (CBEs). We study the CBEs absorption via subconduction-band (sub-CB) transition, and develop a coupled avalanche model. Our results indicate that the CBEs absorption via sub-CB transition plays an important role in the damage in dielectrics irradiated by the visible and near ultraviolet femtosecond lasers. Our theory explains well the experiments.

Fluorescent Ag nanoclusters in glass induced by an infrared femtosecond laser

We report a method for the selective introduction of fluorescent Ag nanoclusters in glass. Extinction and photoluminescence spectra show that a fraction of the Ag atoms are generated through femtosecond laser induced multiphoton reduction and then aggregate to form Ag nanoclusters after heat treatment. Red luminescence from the irradiated region is observed under blue or green laser excitation. The fluorescence can be attributed to interband transitions within Ag nanoclusters. This method provides a novel route to fabricate fluorescent nanomaterials in 3D transparent materials. (c) 2007 Elsevier B.V. All rights reserved.

The heat capacity of superfluid ^4He in the presence of a constant heat flux near Tλ

We present the first experimental evidence that the heat capacity of superfluid ⁴He, at temperatures very close to the lambda transition temperature, T_λ,is enhanced by a constant heat flux, Q. The heat capacity at constant Q, C_Q,is predicted to diverge at a temperature T_c(Q) < T_λ at which superflow becomes unstable. In agreement with previous measurements, we find that dissipation enters our cell at a temperature, T_DAS(Q),below the theoretical value, T_c(Q). Our measurements of C_Q were taken using the discrete pulse method at fourteen different heat flux values in the range 1µW/cm² ≤ Q≤ 4µW /cm². The excess heat capacity ∆C_Q we measure has the predicted scaling behavior as a function of T and Q:∆C_Q • t^α ∝ (Q/Q_c)², where Q_cT) ~ t^2ν is the critical heat current that results from the inversion of the equation for T_c(Q). We find that if the theoretical value of T_c( Q) is correct, then ∆C_Q is considerably larger than anticipated. On the other hand,if T_c(Q)≈ T_DAS(Q),then ∆C_Q is the same magnitude as the theoretically predicted enhancement.

Did galaxies reionize the universe?

The epoch of reionization remains one of the last uncharted eras of cosmic history, yet this time is of crucial importance, encompassing the formation of both the first galaxies and the first metals in the universe. In this thesis, I present four related projects that both characterize the abundance and properties of these first galaxies and uses follow-up observations of these galaxies to achieve one of the first observations of the neutral fraction of the intergalactic medium during the heart of the reionization era.

First, we present the results of a spectroscopic survey using the Keck telescopes targeting 6.3 < z < 8.8 star-forming galaxies. We secured observations of 19 candidates, initially selected by applying the Lyman break technique to infrared imaging data from the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST). This survey builds upon earlier work from Stark et al. (2010, 2011), which showed that star-forming galaxies at 3 < z < 6, when the universe was highly ionized, displayed a significant increase in strong Lyman alpha emission with redshift. Our work uses the LRIS and NIRSPEC instruments to search for Lyman alpha emission in candidates at a greater redshift in the observed near-infrared, in order to discern if this evolution continues, or is quenched by an increase in the neutral fraction of the intergalactic medium. Our spectroscopic observations typically reach a 5-sigma limiting sensitivity of < 50 AA. Despite expecting to detect Lyman alpha at 5-sigma in 7-8 galaxies based on our Monte Carlo simulations, we only achieve secure detections in two of 19 sources. Combining these results with a similar sample of 7 galaxies from Fontana et al. (2010), we determine that these few detections would only occur in < 1% of simulations if the intrinsic distribution was the same as that at z ~ 6. We consider other explanations for this decline, but find the most convincing explanation to be an increase in the neutral fraction of the intergalactic medium. Using theoretical models, we infer a neutral fraction of X_HI ~ 0.44 at z = 7.

Second, we characterize the abundance of star-forming galaxies at z > 6.5 again using WFC3 onboard the HST. This project conducted a detailed search for candidates both in the Hubble Ultra Deep Field as well as a number of additional wider Hubble Space Telescope surveys to construct luminosity functions at both z ~ 7 and 8, reaching 0.65 and 0.25 mag fainter than any previous surveys, respectively. With this increased depth, we achieve some of the most robust constraints on the Schechter function faint end slopes at these redshifts, finding very steep values of alpha_{z~7} = -1.87 +/- 0.18 and alpha_{z~8} = -1.94 +/- 0.23. We discuss these results in the context of cosmic reionization, and show that given reasonable assumptions about the ionizing spectra and escape fraction of ionizing photons, only half the photons needed to maintain reionization are provided by currently observable galaxies at z ~ 7-8. We show that an extension of the luminosity function down to M_{UV} = -13.0, coupled with a low level of star-formation out to higher redshift, can fit all available constraints on the ionization history of the universe.

Third, we investigate the strength of nebular emission in 3 < z < 5 star-forming galaxies. We begin by using the Infrared Array Camera (IRAC) onboard the Spitzer Space Telescope to investigate the strength of H alpha emission in a sample of 3.8 < z < 5.0 spectroscopically confirmed galaxies. We then conduct near-infrared observations of star-forming galaxies at 3 < z < 3.8 to investigate the strength of the [OIII] 4959/5007 and H beta emission lines from the ground using MOSFIRE. In both cases, we uncover near-ubiquitous strong nebular emission, and find excellent agreement between the fluxes derived using the separate methods. For a subset of 9 objects in our MOSFIRE sample that have secure Spitzer IRAC detections, we compare the emission line flux derived from the excess in the K_s band photometry to that derived from direct spectroscopy and find 7 to agree within a factor of 1.6, with only one catastrophic outlier. Finally, for a different subset for which we also have DEIMOS rest-UV spectroscopy, we compare the relative velocities of Lyman alpha and the rest-optical nebular lines which should trace the cites of star-formation. We find a median velocity offset of only v_{Ly alpha} = 149 km/s, significantly less than the 400 km/s observed for star-forming galaxies with weaker Lyman alpha emission at z = 2-3 (Steidel et al. 2010), and show that this decrease can be explained by a decrease in the neutral hydrogen column density covering the galaxy. We discuss how this will imply a lower neutral fraction for a given observed extinction of Lyman alpha when its visibility is used to probe the ionization state of the intergalactic medium.

Finally, we utilize the recent CANDELS wide-field, infra-red photometry over the GOODS-N and S fields to re-analyze the use of Lyman alpha emission to evaluate the neutrality of the intergalactic medium. With this new data, we derive accurate ultraviolet spectral slopes for a sample of 468 3 < z < 6 star-forming galaxies, already observed in the rest-UV with the Keck spectroscopic survey (Stark et al. 2010). We use a Bayesian fitting method which accurately accounts for contamination and obscuration by skylines to derive a relationship between the UV-slope of a galaxy and its intrinsic Lyman alpha equivalent width probability distribution. We then apply this data to spectroscopic surveys during the reionization era, including our own, to accurately interpret the drop in observed Lyman alpha emission. From our most recent such MOSFIRE survey, we also present evidence for the most distant galaxy confirmed through emission line spectroscopy at z = 7.62, as well as a first detection of the CIII]1907/1909 doublet at z > 7.

We conclude the thesis by exploring future prospects and summarizing the results of Robertson et al. (2013). This work synthesizes many of the measurements in this thesis, along with external constraints, to create a model of reionization that fits nearly all available constraints.

I. The effect of droplet solidification upon two-phase flow in a rocket nozzle. II. A kinetic theory investigation of some condensation-evaporation phenomena by a moment method

Part I:

The perturbation technique developed by Rannie and Marble is used to study the effect of droplet solidification upon two-phase flow in a rocket nozzle. It is shown that under certain conditions an equilibrium flow exists, where the gas and particle phases have the same velocity and temperature at each section of the nozzle. The flow is divided into three regions: the first region, where the particles are all in the form of liquid droplets; a second region, over which the droplets solidify at constant freezing temperature; and a third region, where the particles are all solid. By a perturbation about the equilibrium flow, a solution is obtained for small particle slip velocities using the Stokes drag law and the corresponding approximation for heat transfer between the particle and gas phases. Singular perturbation procedure is required to handle the problem at points where solidification first starts and where it is complete. The effects of solidification are noticeable.

Part II:

When a liquid surface, in contact with only its pure vapor, is not in the thermodynamic equilibrium with it, a net condensation or evaporation of fluid occurs. This phenomenon is studied from a kinetic theory viewpoint by means of moment method developed by Lees. The evaporation-condensation rate is calculated for a spherical droplet and for a liquid sheet, when the temperatures and pressures are not too far removed from their equilibrium values. The solutions are valid for the whole range of Knudsen numbers from the free molecule to the continuum limit. In the continuum limit, the mass flux rate is proportional to the pressure difference alone.

A kinetic theory description for external spherical flows with arbitrary Knudsen number by a moment method

The Maxwell integral equations of transfer are applied to a series of problems involving flows of arbitrary density gases about spheres. As suggested by Lees a two sided Maxwellian-like weighting function containing a number of free parameters is utilized and a sufficient number of partial differential moment equations is used to determine these parameters. Maxwell's inverse fifth-power force law is used to simplify the evaluation of the collision integrals appearing in the moment equations. All flow quantities are then determined by integration of the weighting function which results from the solution of the differential moment system. Three problems are treated: the heat-flux from a slightly heated sphere at rest in an infinite gas; the velocity field and drag of a slowly moving sphere in an unbounded space; the velocity field and drag torque on a slowly rotating sphere. Solutions to the third problem are found to both first and second-order in surface Mach number with the secondary centrifugal fan motion being of particular interest. Singular aspects of the moment method are encountered in the last two problems and an asymptotic study of these difficulties leads to a formal criterion for a "well posed" moment system. The previously unanswered question of just how many moments must be used in a specific problem is now clarified to a great extent.

Transparent Ni2+-doped MgO-Al2O3-SiO2 glass ceramics with broadband infrared luminescence

We report on transparent Ni2+-doped MgO-Al2O3-SiO2 glass ceramics with broadband infrared luminescence. Ni2+-doped MgO-Al2O3-SiO2 glass is Prepared by using the conventional method. After heat treatment at high temperature, MgAl2O4 crystallites are precipitated, and their average size is about 4.3nm. No luminescence is detected in the as-prepared glass sample, while broadband infrared luminescence centred at around 1315nm with full width at half maximum (FWHM) of about 300nm is observed from the glass ceramics. The observed infrared emission could be attributed to the T-3(2g)(F-3) -> (3)A(2g)(F-3) transition of octahedral Ni2+ ions in the MgAl2O4 crystallites of the transparent glass ceramics. The product of the fluorescence lifetime and the stimulated emission cross section is about 1.6 X 10(-24) s cm(2).