Small scale turbulence in high Karlovitz number premixed flames

The purpose of this thesis is to characterize the behavior of the smallest turbulent scales in high Karlovitz number (Ka) premixed flames. These scales are particularly important in the two-way coupling between turbulence and chemistry and better understanding of these scales will support future modeling efforts using large eddy simulations (LES). The smallest turbulent scales are studied by considering the vorticity vector, ω, and its transport equation.

Due to the complexity of turbulent combustion introduced by the wide range of length and time scales, the two-dimensional vortex-flame interaction is first studied as a simplified test case. Numerical and analytical techniques are used to discern the dominate transport terms and their effects on vorticity based on the initial size and strength of the vortex. This description of the effects of the flame on a vortex provides a foundation for investigating vorticity in turbulent combustion.

Subsequently, enstrophy, ω² = ω • ω, and its transport equation are investigated in premixed turbulent combustion. For this purpose, a series of direct numerical simulations (DNS) of premixed n-heptane/air flames are performed, the conditions of which span a wide range of unburnt Karlovitz numbers and turbulent Reynolds numbers. Theoretical scaling analysis along with the DNS results support that, at high Karlovitz number, enstrophy transport is controlled by the viscous dissipation and vortex stretching/production terms. As a result, vorticity scales throughout the flame with the inverse of the Kolmogorov time scale, τ_η, just as in homogeneous isotropic turbulence. As τ_η is only a function of the viscosity and dissipation rate, this supports the validity of Kolmogorov’s first similarity hypothesis for sufficiently high Ka numbers (Ka ≳ 100). These conclusions are in contrast to low Karlovitz number behavior, where dilatation and baroclinic torque have a significant impact on vorticity within the flame. Results are unaffected by the transport model, chemical model, turbulent Reynolds number, and lastly the physical configuration.

Next, the isotropy of vorticity is assessed. It is found that given a sufficiently large value of the Karlovitz number (Ka ≳ 100) the vorticity is isotropic. At lower Karlovitz numbers, anisotropy develops due to the effects of the flame on the vortex stretching/production term. In this case, the local dynamics of vorticity in the strain-rate tensor, S, eigenframe are altered by the flame. At sufficiently high Karlovitz numbers, the dynamics of vorticity in this eigenframe resemble that of homogeneous isotropic turbulence.

Combined, the results of this thesis support that both the magnitude and orientation of vorticity resemble the behavior of homogeneous isotropic turbulence, given a sufficiently high Karlovitz number (Ka ≳ 100). This supports the validity of Kolmogorov’s first similarity hypothesis and the hypothesis of local isotropy under these condition. However, dramatically different behavior is found at lower Karlovitz numbers. These conclusions provides/suggests directions for modeling high Karlovitz number premixed flames using LES. With more accurate models, the design of aircraft combustors and other combustion based devices may better mitigate the detrimental effects of combustion, from reducing CO₂ and soot production to increasing engine efficiency.

Offshore coho salmon populations near the Pacific Northwest and large-scale atmospheric events

EXTRACT (SEE PDF FOR FULL ABSTRACT): Catch of coho salmon off the coast of Washington and Oregon since 1925 appears to be related to large-scale events in the atmosphere, which in turn affect ocean currents and coastal upwelling intensities in the northeastern Pacific. At least two time scales of variations can be identified. The first is that of the El Nino/Southern Oscillation phenomenon giving rise to an irregular cycle of between 3 to 7 years. ... The second time scale of variation seems to have a periodicity of about 20 years, although this is based on a limited dataset. ... This paper endeavors to describe how, if real, these atmospheric/oceanic effects are integrated and might affect the salmon catch. The possibility must also be considered that the atmospheric events are symbiotically related to the oceanic events and, further, that both may be enmeshed in even longer-term variability of climate.

TIME-RESOLVED RESONANCE RAMAN-SPECTROSCOPY AND RAMAN SPECTROELECTROCHEMISTRY OF (CO)(5)W[4,4'-BPY]W(CO)(5), PROBED IN THE VISIBLE AND NEAR-INFRARED

Time-resolved resonance Raman spectroscopy of the lowest energy excited state of the 4,4'-bipyridyl ligand-bridged complex, [(CO)(5)W(L)W(CO5] (1), and Raman spectroscopy of electrochemically reduced 1, both give bands characteristic of the the L(.-) species. This confirms that the ligand L is negatively charged in the lowest energy exicited state which is therefore metal-ligand charge transfer (MLCT) in character. Raman spectra of the radical anion of 1 excited in the far red (800 nm) exhibited a band near 2050 cm(-1) due to a vco symmetric CO stretching mode, compared to the corresponding band at 2070 cm(-1) in the spectrum of the parent, uncharged complex. The lower vco in the reduced complex supports the recent finding by time-resolved IR spectroscopy of a similar frequency decrease for nu(CO) in the longest lived (MLCT) excited state of 1 which was attributed to electron/hole localisation in this state on the IR time scale.

Calculus of variations on time scales and applications to economics

We consider some problems of the calculus of variations on time scales. On the beginning our attention is paid on two inverse extremal problems on arbitrary time scales. Firstly, using the Euler-Lagrange equation and the strengthened Legendre condition, we derive a general form for a variation functional that attains a local minimum at a given point of the vector space. Furthermore, we prove a necessary condition for a dynamic integro-differential equation to be an Euler-Lagrange equation. New and interesting results for the discrete and quantum calculus are obtained as particular cases. Afterwards, we prove Euler-Lagrange type equations and transversality conditions for generalized infinite horizon problems. Next we investigate the composition of a certain scalar function with delta and nabla integrals of a vector valued field. Euler-Lagrange equations in integral form, transversality conditions, and necessary optimality conditions for isoperimetric problems, on an arbitrary time scale, are proved. In the end, two main issues of application of time scales in economic, with interesting results, are presented. In the former case we consider a firm that wants to program its production and investment policies to reach a given production rate and to maximize its future market competitiveness. The model which describes firm activities is studied in two different ways: using classical discretizations; and applying discrete versions of our result on time scales. In the end we compare the cost functional values obtained from those two approaches. The latter problem is more complex and relates to rate of inflation, p, and rate of unemployment, u, which inflict a social loss. Using known relations between p, u, and the expected rate of inflation π, we rewrite the social loss function as a function of π. We present this model in the time scale framework and find an optimal path π that minimizes the total social loss over a given time interval.

Diatom-based characterization of iberian coastal environments at different time scales

Tese de doutoramento, Geologia (Geologia Económica e do Ambiente), Universidade de Lisboa, Faculdade de Ciências, 2014

Significance of Time Scales in Nonlinear Dynamical Analysis of Electroencephalogram Signals

We propose to show in this paper, that the time series obtained from biological systems such as human brain are invariably nonstationary because of different time scales involved in the dynamical process. This makes the invariant parameters time dependent. We made a global analysis of the EEG data obtained from the eight locations on the skull space and studied simultaneously the dynamical characteristics from various parts of the brain. We have proved that the dynamical parameters are sensitive to the time scales and hence in the study of brain one must identify all relevant time scales involved in the process to get an insight in the working of brain.

Femtosecond real-time probing of reactions. 12. Vectorial dynamics of transition states

Femtosecond time-resolved techniques with KETOF (kinetic energy time-of-flight) detection in a molecular beam are developed for studies of the vectorial dynamics of transition states. Application to the dissociation reaction of IHgI is presented. For this system, the complex [I---Hg---I](++)* is unstable and, through the symmetric and asymmetric stretch motions, yields different product fragments: [I---Hg---I](++)* -> HgI(X^2/sigma^+) + I(^2P_3/2) [or I*(^2P_l/2)] (1a); [I---Hg---I](++)* -> Hg(^1S_0) + I(^2P_3/2) + I(^2P_3/2) [or I* (^2P_1/2)] (1 b). These two channels, (1a) and (1b), lead to different kinetic energy distributions in the products. It is shown that the motion of the wave packet in the transition-state region can be observed by MPI mass detection; the transient time ranges from 120 to 300 fs depending on the available energy. With polarized pulses, the vectorial properties (transition moments alignment relative to recoil direction) are studied for fragment separations on the femtosecond time scale. The results indicate the nature of the structure (symmetry properties) and the correlation to final products. For 311-nm excitation, no evidence of crossing between the I and I* potentials is found at the internuclear separations studied. (Results for 287-nm excitation are also presented.) Molecular dynamics simulations and studies by laser-induced fluorescence support these findings.

Femtosecond time-resolved wave packet motion in molecular multiphoton ionization and fragmentation

The dynamics of molecular multiphoton ionization and fragmentation of a diatomic molecule (Na_2) have been studied in molecular beam experiments. Femtosecond laser pulses from an amplified colliding-pulse mode-locked (CPM) ring dye laser are employed to induce and probe the molecular transitions. The final continuum states are analyzed by photoelectron spectroscopy, by ion mass spectrometry and by measuring the kinetic energy of the formed ionic fragments. Pump-probe spectra employing 70-fs laser pulses have been measured to study the time dependence of molecular multiphoton ionization and fragmentation. The oscillatory structure of the transient spectra showing the dynamics on the femtosecond time scale can best be understood in terms of the motion of wave packets in bound molecular potentials. The transient Na_2^+ ionization and the transient Na^+ fragmentation spectra show that contributions from direct photoionization of a singly excited electronic state and from excitation and autoionization of a bound doubly excited molecular state determine the time evolution of molecular multiphoton ionization.

Femtosecond real-time probing of reactions

Femtosecond reaction dynamics of OClO in a supersonic molecular beam are reported. The system is excited to the A^2A_2 state with a femtosecond pulse, covering a range of excitation in the symmetric stretch between v_1 = 17 to v_1 = 11 (308-352 nm). A time-delayed femtosecond probe pulse ionizes the OClO, and OClO^+ is detected. This ion has not been observed in previous experiments because of its ultrafast fragmentation. Transients are reported for the mass of the parent OClO as well as the mass of the ClO. Apparent biexponential decays are observed and related to the fragmentation dynamics: OClO+hv \rightarrow (OClO)^{(++)*} \rightarrow ClO+O \rightarrow Cl+O_2. Clusters of OClO with water (OClO)_n (H_2 0)_m with n from 1 to 3 and m from 0 to 3 are also observed. The dynamics of the fragmentation reveal the nuclear motions and the electronic coupling between surfaces. The time scale for bond breakage is in the range of 300-500 fs, depending on v_1; surface crossing to form new intermediates is a pathway for the two channels of fragmentation: ClO+O (primary) and Cl+O_2 (minor). Comparisons with results of ab initio calculations are made.

Horizontal potential vorticity dipoles on the convective storm scale

The structure and dynamics of potential vorticity (PV) anomalies generated by convective storms is investigated both theoretically and in a numerical model case study. Linear theory suggests that if the storm-induced heating is on a sufficiently small scale (relative to the Rossby radius of deformation), and the environment contains moderate vertical wind shear (of order 1 m s(-1) km(-1)), then the dominant mode of a diabatically generated PV anomaly is a horizontally oriented dipole. The horizontal dipoles are typically of O(10 PVU), compared with the O(1 PVU) vertical dipoles that have been studied extensively throughout the literature. Furthermore, the horizontal PV dipoles are realized almost entirely as relative vorticity anomalies (on a time-scale of the order of tens of minutes after the heating has been turned on). The analysis of horizontal PV dipoles offers a new perspective on the vorticity dynamics of individual convective cells, implying that moist processes play a role in the maintenance of vertical vorticity in the convective storm environment.

Assessing levels of uncertainty in recent temperature time series

We examine to what degree we can expect to obtain accurate temperature trends for the last two decades near the surface and in the lower troposphere. We compare temperatures obtained from surface observations and radiosondes as well as satellite-based measurements from the Microwave Soundings Units (MSU), which have been adjusted for orbital decay and non-linear instrument-body effects, and reanalyses from the European Centre for Medium-Range Weather Forecasts (ERA) and the National Centre for Environmental Prediction (NCEP). In regions with abundant conventional data coverage, where the MSU has no major influence on the reanalysis, temperature anomalies obtained from microwave sounders, radiosondes and from both reanalyses agree reasonably. Where coverage is insufficient, in particular over the tropical oceans, large differences are found between the MSU and either reanalysis. These differences apparently relate to changes in the satellite data availability and to differing satellite retrieval methodologies, to which both reanalyses are quite sensitive over the oceans. For NCEP, this results from the use of raw radiances directly incorporated into the analysis, which make the reanalysis sensitive to changes in the underlying algorithms, e.g. those introduced in August 1992. For ERA, the bias-correction of the one-dimensional variational analysis may introduce an error when the satellite relative to which the correction is calculated is biased itself or when radiances change on a time scale longer than a couple of months, e.g. due to orbit decay. ERA inhomogeneities are apparent in April 1985, October/November 1986 and April 1989. These dates can be identified with the replacements of satellites. It is possible that a negative bias in the sea surface temperatures (SSTs) used in the reanalyses may have been introduced over the period of the satellite record. This could have resulted from a decrease in the number of ship measurements, a concomitant increase in the importance of satellite-derived SSTs, and a likely cold bias in the latter. Alternately, a warm bias in SSTs could have been caused by an increase in the percentage of buoy measurements (relative to deeper ship intake measurements) in the tropical Pacific. No indications for uncorrected inhomogeneities of land surface temperatures could be found. Near-surface temperatures have biases in the boundary layer in both reanalyses, presumably due to the incorrect treatment of snow cover. The increase of near-surface compared to lower tropospheric temperatures in the last two decades may be due to a combination of several factors, including high-latitude near-surface winter warming due to an enhanced NAO and upper-tropospheric cooling due to stratospheric ozone decrease.

Combined use of satellite estimates and rain gauge observations to generate high-quality historical rainfall time series over Ethiopia

Climate data are used in a number of applications including climate risk management and adaptation to climate change. However, the availability of climate data, particularly throughout rural Africa, is very limited. Available weather stations are unevenly distributed and mainly located along main roads in cities and towns. This imposes severe limitations to the availability of climate information and services for the rural community where, arguably, these services are needed most. Weather station data also suffer from gaps in the time series. Satellite proxies, particularly satellite rainfall estimate, have been used as alternatives because of their availability even over remote parts of the world. However, satellite rainfall estimates also suffer from a number of critical shortcomings that include heterogeneous time series, short time period of observation, and poor accuracy particularly at higher temporal and spatial resolutions. An attempt is made here to alleviate these problems by combining station measurements with the complete spatial coverage of satellite rainfall estimates. Rain gauge observations are merged with a locally calibrated version of the TAMSAT satellite rainfall estimates to produce over 30-years (1983-todate) of rainfall estimates over Ethiopia at a spatial resolution of 10 km and a ten-daily time scale. This involves quality control of rain gauge data, generating locally calibrated version of the TAMSAT rainfall estimates, and combining these with rain gauge observations from national station network. The infrared-only satellite rainfall estimates produced using a relatively simple TAMSAT algorithm performed as good as or even better than other satellite rainfall products that use passive microwave inputs and more sophisticated algorithms. There is no substantial difference between the gridded-gauge and combined gauge-satellite products over the test area in Ethiopia having a dense station network; however, the combined product exhibits better quality over parts of the country where stations are sparsely distributed.

The population dynamics of disease on short and long time-scales

Observational evidence is scarce concerning the distribution of plant pathogen population sizes or densities as a function of time-scale or spatial scale. For wild pathosystems we can only get indirect evidence from evolutionary patterns and the consequences of biological invasions.We have little or no evidence bearing on extermination of hosts by pathogens, or successful escape of a host from a pathogen. Evidence over the last couple of centuries from crops suggest that the abundance of particular pathogens in the spectrum affecting a given host can vary hugely on decadal timescales. However, this may be an artefact of domestication and intensive cultivation. Host-pathogen dynamics can be formulated mathematically fairly easily–for example as SIR-type differential equation or difference equation models, and this has been the (successful) focus of recent work in crops. “Long-term” is then discussed in terms of the time taken to relax from a perturbation to the asymptotic state. However, both host and pathogen dynamics are driven by environmental factors as well as their mutual interactions, and both host and pathogen co-evolve, and evolve in response to external factors. We have virtually no information about the importance and natural role of higher trophic levels (hyperpathogens) and competitors, but they could also induce long-scale fluctuations in the abundance of pathogens on particular hosts. In wild pathosystems the host distribution cannot be modelled as either a uniform density or even a uniform distribution of fields (which could then be treated as individuals). Patterns of short term density-dependence and the detail of host distribution are therefore critical to long-term dynamics. Host density distributions are not usually scale-free, but are rarely uniform or clearly structured on a single scale. In a (multiply structured) metapopulation with coevolution and external disturbances it could well be the case that the time required to attain equilibrium (if it exists) based on conditions stable over a specified time-scale is longer than that time-scale. Alternatively, local equilibria may be reached fairly rapidly following perturbations but the meta-population equilibrium be attained very slowly. In either case, meta-stability on various time-scales is a more relevant than equilibrium concepts in explaining observed patterns.

Response time of the high-latitude dayside ionosphere to sudden changes in the north-south component of the IMF

The time scale of the response of the high-latitude dayside ionospheric flow to changes in the North-South component of the interplanetary magnetic field (IMF) has been investigated by examining the time delays between corresponding sudden changes. Approximately 40 h of simultaneous IMF and ionospheric flow data have been examined, obtained by the AMPTE-UKS and -IRM spacecraft and the EISCAT “Polar” experiment, respectively, in which 20 corresponding sudden changes have been identified. Ten of these changes were associated with southward turnings of the IMF, and 10 with northward turnings. It has been found that the corresponding flow changes occurred simultaneously over the whole of the “Polar” field-of-view, extending more than 2° in invariant latitude, and that the ionospheric response delay following northward turnings is the same as that following southward turnings, though the form of the response is different in the two cases. The shortest response time, 5.5 ± 3.2 min, is found in the early- to mid-afternoon sector, increasing to 9.5 ± 3.0 min in the mid-morning sector, and to 9.5 ± 3.1 min near to dusk. These times represent the delays in the appearance of perturbed flows in the “Polar” field-of-view following the arrival of IMF changes at the subsolar magnetopause. Overall, the results agree very well with those derived by Etemadi et al. (1988, Planet. Space Sci.36, 471) from a general cross-correlation analysis of the IMF Bz and “Polar” beam-swinging vector flow data.

Cauchy-Stieltjes Integral on Time Scales in Banach Spaces and Hysteresis Operators

The play operator has a fundamental importance in the theory of hysteresis. It was studied in various settings as shown by P. Krejci and Ph. Laurencot in 2002. In that work it was considered the Young integral in the frame of Hilbert spaces. Here we study the play in the frame of the regulated functions (that is: the ones having only discontinuities of the first kind) on a general time scale T (that is: with T being a nonempty closed set of real numbers) with values in a Banach space. We will be showing that the dual space in this case will be defined as the space of operators of bounded semivariation if we consider as the bilinearity pairing the Cauchy-Stieltjes integral on time scales.