General-domain compressible Navier-Stokes solvers exhibiting quasi-unconditional stability and high-order accuracy in space and time

This thesis presents a new class of solvers for the subsonic compressible Navier-Stokes equations in general two- and three-dimensional spatial domains. The proposed methodology incorporates: 1) A novel linear-cost implicit solver based on use of higher-order backward differentiation formulae (BDF) and the alternating direction implicit approach (ADI); 2) A fast explicit solver; 3) Dispersionless spectral spatial discretizations; and 4) A domain decomposition strategy that negotiates the interactions between the implicit and explicit domains. In particular, the implicit methodology is quasi-unconditionally stable (it does not suffer from CFL constraints for adequately resolved flows), and it can deliver orders of time accuracy between two and six in the presence of general boundary conditions. In fact this thesis presents, for the first time in the literature, high-order time-convergence curves for Navier-Stokes solvers based on the ADI strategy---previous ADI solvers for the Navier-Stokes equations have not demonstrated orders of temporal accuracy higher than one. An extended discussion is presented in this thesis which places on a solid theoretical basis the observed quasi-unconditional stability of the methods of orders two through six. The performance of the proposed solvers is favorable. For example, a two-dimensional rough-surface configuration including boundary layer effects at Reynolds number equal to one million and Mach number 0.85 (with a well-resolved boundary layer, run up to a sufficiently long time that single vortices travel the entire spatial extent of the domain, and with spatial mesh sizes near the wall of the order of one hundred-thousandth the length of the domain) was successfully tackled in a relatively short (approximately thirty-hour) single-core run; for such discretizations an explicit solver would require truly prohibitive computing times. As demonstrated via a variety of numerical experiments in two- and three-dimensions, further, the proposed multi-domain parallel implicit-explicit implementations exhibit high-order convergence in space and time, useful stability properties, limited dispersion, and high parallel efficiency.

The instability of fluids with time dependent heating

The stability of a fluid having a non-uniform temperature stratification is examined analytically for the response of infinitesimal disturbances. The growth rates of disturbances have been established for a semi-infinite fluid for Rayleigh numbers of 10³, 10⁴, and 10⁵ and for Prandtl numbers of 7.0 and 0.7.

The critical Rayleigh number for a semi-infinite fluid, based on the effective fluid depth, is found to be 32, while it is shown that for a finite fluid layer the critical Rayleigh number depends on the rate of heating. The minimum critical Rayleigh number, based on the depth of a fluid layer, is found to be 1340.

The stability of a finite fluid layer is examined for two special forms of heating. The first is constant flux heating, while in the second, the temperature of the lower surface is increased uniformly in time. In both cases, it is shown that for moderate rates of heating the critical Rayleigh number is reduced, over the value for very slow heating, while for very rapid heating the critical Rayleigh number is greatly increased. These results agree with published experimental observations.

The question of steady, non-cellular convection is given qualitative consideration. It is concluded that, although the motion may originate from infinitesimal disturbances during non-uniform heating, the final flow field is intrinsically non-linear.

Formation of jets by implusive acceleration of a curved free surface

The sudden axial acceleration of a column of liquid bounded at one end by a concave free surface has been found, experimentally, to produce a jet which issues from the free surface with a speed several times that imparted to the column.

Theoretical approximations to such flows, valid for small time, are formulated subject to the assumption that the fluid is inviscid and incompressible. In a special two-dimensional case, it is found that, for vanishingly small time, the velocity at the point on the free surface from which the jet emanates is π/2 times the velocity imparted to the column. The solutions to several problems in two and three dimensions assuming that the initial curvature of the free surface is small, lead to values for this ratio dependent upon the curvature—the initial velocity in the case of axial symmetry exceeding that of the analogous two-dimensional problem by approximately 25%.

Experiments conducted upon the phenomenon give values systematically in excess of those predicted by the theory, although theory and experiment are in qualitative agreement with respect to the displacement of the free surface. It is suggested that the discrepancy is attributable to effects of finite curvature having been imperfectly accounted for in the axially-symmetric analysis.

Photographic materials on pp. 115, 120, and 121 are essential and will not reproduce clearly on Xerox copies. Photographic copies should be ordered.

Hydrobiological service for observations and control of surface waters in the USSR

Technological progress, having reached in our time an unprecedented speed, is still increasing the rate of mineral extraction, industrial construction, and the mastering of new kinds of energy is growing. Correspondingly the anthropogenic load on the biosphere is increased and that requires the comprehensive development of monitoring the anthropogenic changes in the natural environment. Among problems resulting from the scientific-technological development, a noticeable place is given to the problem of pure water. Surface land waters proved to be a sensitive link in the natural environment. The hydrobiological service for observations and control of the surface waters is one of the subsystems of the State/Federal Service for Observations and Control of pollution levels in environmental objects, conducted by the USSR State Committee for Hydrometeor- ology and Control of the Natural Environment. This paper summarises the the main principles of the organisation and goals of the national system of monitoring of the state of the natural environment in the USSR.

Using Earth Deformation Caused by Surface Mass Loading to Constrain the Elastic Structure of the Crust and Mantle

Surface mass loads come in many different varieties, including the oceans, atmosphere, rivers, lakes, glaciers, ice caps, and snow fields. The loads migrate over Earth's surface on time scales that range from less than a day to many thousand years. The weights of the shifting loads exert normal forces on Earth's surface. Since the Earth is not perfectly rigid, the applied pressure deforms the shape of the solid Earth in a manner controlled by the material properties of Earth's interior. One of the most prominent types of surface mass loading, ocean tidal loading (OTL), comes from the periodic rise and fall in sea-surface height due to the gravitational influence of celestial objects, such as the moon and sun. Depending on geographic location, the surface displacements induced by OTL typically range from millimeters to several centimeters in amplitude, which may be inferred from Global Navigation and Satellite System (GNSS) measurements with sub-millimeter precision. Spatiotemporal characteristics of observed OTL-induced surface displacements may therefore be exploited to probe Earth structure. In this thesis, I present descriptions of contemporary observational and modeling techniques used to explore Earth's deformation response to OTL and other varieties of surface mass loading. With the aim to extract information about Earth's density and elastic structure from observations of the response to OTL, I investigate the sensitivity of OTL-induced surface displacements to perturbations in the material structure. As a case study, I compute and compare the observed and predicted OTL-induced surface displacements for a network of GNSS receivers across South America. The residuals in three distinct and dominant tidal bands are sub-millimeter in amplitude, indicating that modern ocean-tide and elastic-Earth models well predict the observed displacement response in that region. Nevertheless, the sub-millimeter residuals exhibit regional spatial coherency that cannot be explained entirely by random observational uncertainties and that suggests deficiencies in the forward-model assumptions. In particular, the discrepancies may reveal sensitivities to deviations from spherically symmetric, non-rotating, elastic, and isotropic (SNREI) Earth structure due to the presence of the South American craton.

A large-area photolithography technique based on surface plasmons leakage modes

An atomic force microscope (AFM) assisted surface plasmons leakage radiation photolithography technique has been numerically demonstrated by using two-dimensional finite-difference time-domain (2D-FDTD) method. With the aid of a metallic AFM tip, particular characteristic of the Kretstchmann configuration to excite surface plasmons (SPs) is utilized to achieve large-area patterns with high spatial resolution and contrast, the photoresist could be exposed with low power laser due to the remarkable local field enhancement at the metal/dielectric interface and the resonant localized SPs modes near the tip. Good tolerance on the film thickness and incident angle has been obtained, which provides a good practicability for experiments. This photolithography technique proposed here can realize large-area, high-resolution, high-contrast, nondestructive, arbitrary-structure fabrication of nanoscale devices. (c) 2007 Elsevier B.V. All rights reserved.

Quantum ricochets: surface capture, release and energy loss of fast ions hitting a polar surface at grazing incidence

A diffraction mechanism is proposed for the capture, multiple bouncing and final escape of a fast ion (keV) impinging on the surface of a polarizable material at grazing incidence. Capture and escape are effected by elastic quantum diffraction consisting of the exchange of a parallel surface wave vector G= 2p/ a between the ion parallel momentum and the surface periodic potential of period a. Diffraction- assisted capture becomes possible for glancing angles F smaller than a critical value given by Fc 2- 2./ a-| Vim|/ E, where E is the kinetic energy of the ion,. = h/ Mv its de Broglie wavelength and Vim its average electronic image potential at the distance from the surface where diffraction takes place. For F< Fc, the ion can fall into a selected capture state in the quasi- continuous spectrum of its image potential and execute one or several ricochets before being released by the time reversed diffraction process. The capture, ricochet and escape are accompanied by a large, periodic energy loss of several tens of eV in the forward motion caused by the coherent emission of a giant number of quanta h. of Fuchs- Kliewer surface phonons characteristic of the polar material. An analytical calculation of the energy loss spectrum, based on the proposed diffraction process and using a model ion-phonon coupling developed earlier (Lucas et al 2013 J. Phys.: Condens. Matter 25 355009), is presented, which fully explains the experimental spectrum of Villette et al (2000 Phys. Rev. Lett. 85 3137) for Ne+ ions ricocheting on a LiF(001) surface.

Study of the surface plasma transmission properties of a Fabry-Perot resonator by numerical simulation

It is demonstrated with powerful evidence that the extraordinary transmission of a metallic grating is undoubtedly due to the excitation of standing surface plasma waves in the Fabry-Perot like resonator. This is the first time that the strong standing waves set up in the groove of a sub-wavelength double-layer grating (SWDG) for the surface plasma waves have been reported. Moreover, about 90% transmission is gained with an SWDG, more easily fabricated than ordinary metallic gratings, in the first peak of transmission spectrum.

Surface plasmon resonance transmission filters at 1053 nm based on metallic grating with narrow slit

Metallic gratings with narrow slits can lead to special optical properties such as strongly enhancing the transmission and considerably strengthening the polarized effect. A narrow-band filter suitable for application in optical communication is designed by sandwiching a metallic grating between two identical dielectric films. The maximum transmission can reach 96% after optimizing the parameters of films and grating at a central wavelength of 1053 nm. It is the first time, to our knowledge, that such high transmission has been reported since the discovery of the extraordinarily high transmission through periodic holes or slits; moreover, the extremely polarized effect is also found in P mode of this symmetric grating.

Effect of surface plasmon resonance in TiO2/Au thin films on the fluorescence of self-assembled CdTe QDs structure

23rd Congress of the International Comission for Optics (ICO 23)

Tracking Pacific bluefin tuna (Thunnus thynnus orientalis) in the northeastern Pacific with an automated algorithm that estimates latitude by matching sea-surface-temperature data from satellites with temperature data from tags on fish

Data recovered from 11 popup satellite archival tags and 3 surgically implanted archival tags were used to analyze the movement patterns of juvenile northern bluefin tuna (Thunnus thynnus orientalis) in the eastern Pacific. The light sensors on archival and pop-up satellite transmitting archival tags (PSATs) provide data on the time of sunrise and sunset, allowing the calculation of an approximate geographic position of the animal. Light-based estimates of longitude are relatively robust but latitude estimates are prone to large degrees of error, particularly near the times of the equinoxes and when the tag is at low latitudes. Estimating latitude remains a problem for researchers using light-based geolocation algorithms and it has been suggested that sea surface temperature data from satellites may be a useful tool for refining latitude estimates. Tag data from bluefin tuna were subjected to a newly developed algorithm, called “PSAT Tracker,” which automatically matches sea surface temperature data from the tags with sea surface temperatures recorded by satellites. The results of this algorithm compared favorably to the estimates of latitude calculated with the lightbased algorithms and allowed for estimation of fish positions during times of the year when the lightbased algorithms failed. Three near one-year tracks produced by PSAT tracker showed that the fish range from the California−Oregon border to southern Baja California, Mexico, and that the majority of time is spent off the coast of central Baja Mexico. A seasonal movement pattern was evident; the fish spend winter and spring off central Baja California, and summer through fall is spent moving northward to Oregon and returning to Baja California.

A general framework for integrating environmental time series into stock assessment models: model description, simulation testing, and example

We present a method to integrate environmental time series into stock assessment models and to test the significance of correlations between population processes and the environmental time series. Parameters that relate the environmental time series to population processes are included in the stock assessment model, and likelihood ratio tests are used to determine if the parameters improve the fit to the data significantly. Two approaches are considered to integrate the environmental relationship. In the environmental model, the population dynamics process (e.g. recruitment) is proportional to the environmental variable, whereas in the environmental model with process error it is proportional to the environmental variable, but the model allows an additional temporal variation (process error) constrained by a log-normal distribution. The methods are tested by using simulation analysis and compared to the traditional method of correlating model estimates with environmental variables outside the estimation procedure. In the traditional method, the estimates of recruitment were provided by a model that allowed the recruitment only to have a temporal variation constrained by a log-normal distribution. We illustrate the methods by applying them to test the statistical significance of the correlation between sea-surface temperature (SST) and recruitment to the snapper (Pagrus auratus) stock in the Hauraki Gulf–Bay of Plenty, New Zealand. Simulation analyses indicated that the integrated approach with additional process error is superior to the traditional method of correlating model estimates with environmental variables outside the estimation procedure. The results suggest that, for the snapper stock, recruitment is positively correlated with SST at the time of spawning.

Snow accumulation time series and their climatic interpretation

EXTRACT (SEE PDF FOR FULL ABSTRACT): Several snow accumulation time series derived from ice cores and extending over 3 to 5 centuries are examined for spatial and temporal climatic information. ... A significant observation is the widespread depression of net snow accumulation during the latter part of the "Little Ice Age". This initially suggests sea surface temperatures were significantly depressed during the same period. However, prior to this, the available core records indicate generally higher than average precipitation rates. This also implies that influences such as shifted storm tracks or a dustier atmosphere may have been involved. Without additional spatial data coverage, these observations should properly be studied using a coupled (global) ocean/atmosphere GCM.

Variability in a radiolarian time series and its relationship to climate

EXTRACT (SEE PDF FOR FULL ABSTRACT): Recently, paleoceanographers have been challenged to produce reliable proxies of climate variables that can be incorporated into climate models. In developing proxies using time series of annual radiolarian species fluxes from Santa Barbara Basin, we identify groups of species associated with years of extreme sea surface temperatures and sea level heights.

ZnO thin film surface acoustic wave based lab-on-a-chip

Lab-on-a-chip (LOC) is one of the most important microsystem applications with promise for use in microanalysis, drug development, diagnosis of illness and diseases etc. LOC typically consists of two main components: microfluidics and sensors. Integration of microfluidics and sensors on a single chip can greatly enhance the efficiency of biochemical reactions and the sensitivity of detection, increase the reaction/detection speed, and reduce the potential cross-contamination, fabrication time and cost etc. However, the mechanisms generally used for microfluidics and sensors are different, making the integration of the two main components complicated and increases the cost of the systems. A lab-on-a-chip system based on a single surface acoustic wave (SAW) actuation mechanism is proposed. SAW devices were fabricated on nanocrystalline ZnO thin films deposited on Si substrates using sputtering. Coupling of acoustic waves into a liquid induces acoustic streaming and motion of droplets. A streaming velocity up to ∼ 5cm/s and droplet pumping speeds of ∼lcm/s were obtained. It was also found that a higher order mode wave, the Sezawa wave is more effective in streaming and transportation of microdroplets. The ZnO SAW sensor has been used for prostate antigen/antibody biorecognition systems, demonstrated the feasibility of using a single actuation mechanism for lab-on-a-chip applications. © 2010 Materials Research Society.