Numerical Investigation On Evolution Of Cylindrical Cellular Detonation

Cylindrical cellular detonation is numerically investigated by solving two-dimensional reactive Euler equations with a finite volume method on a two-dimensional self-adaptive unstructured mesh. The one-step reversible chemical reaction model is applied to simplify the control parameters of chemical reaction. Numerical results demonstrate the evolution of cellular cell splitting of cylindrical cellular detonation explored in experimentas. Split of cellular structures shows different features in the near-field and far-field from the initiation zone. Variation of the local curvature is a key factor in the behavior of cell split of cylindrical cellular detonation in propagation. Numerical results show that split of cellular structures comes from the self-organization of transverse waves corresponding to the development of small disturbances along the detonation front related to detonation instability.

First-passage failure of quasi-integrable Hamiltonian systems

The first-passage failure of quasi-integrable Hamiltonian si-stems (multidegree-of-freedom integrable Hamiltonian systems subject to light dampings and weakly random excitations) is investigated. The motion equations of such a system are first reduced to a set of averaged Ito stochastic differential equations by using the stochastic averaging method for quasi-integrable Hamiltonian systems. Then, a backward Kolmogorov equation governing the conditional reliability function and a set of generalized Pontryagin equations governing the conditional moments of first-passage time are established. Finally, the conditional reliability function, and the conditional probability density and moments of first-passage time are obtained by solving these equations with suitable initial and boundary conditions. Two examples are given to illustrate the proposed procedure and the results from digital simulation are obtained to verify the effectiveness of the procedure.

First-passage time of Duffing oscillator under combined harmonic and white-noise excitations

The first-passage time of Duffing oscillator under combined harmonic and white-noise excitations is studied. The equation of motion of the system is first reduced to a set of averaged Ito stochastic differential equations by using the stochastic averaging method. Then, a backward Kolmogorov equation governing the conditional reliability function and a set of generalized Pontryagin equations governing the conditional moments of first-passage time are established. Finally, the conditional reliability function, and the conditional probability density and moments of first-passage time are obtained by solving the backward Kolmogorov equation and generalized Pontryagin equations with suitable initial and boundary conditions. Numerical results for two resonant cases with several sets of parameter values are obtained and the analytical results are verified by using those from digital simulation.

Magnetohydrodynamics equilibrium of a self-confined elliptic plasma ball

A variational principle is applied to the problem of magnetohydrodynamics (MHD) equilibrium of a self-contained elliptical plasma ball, such as elliptical ball lightning. The principle is appropriate for an approximate solution of partial differential equations with arbitrary boundary shape. The method reduces the partial differential equation to a series of ordinary differential equations and is especially valuable for treating boundaries with nonlinear deformations. The calculations conclude that the pressure distribution and the poloidal current are more uniform in an oblate self-confined plasma ball than that of an elongated plasma ball. The ellipticity of the plasma ball is obviously restricted by its internal pressure, magnetic field, and ambient pressure. Qualitative evidence is presented for the absence of sighting of elongated ball lightning.

The thermal and mechanical structure of a two-dimensional plume in the earths mantle

On the condition that the distribution of velocity and temperature at the mid-plane of a mantle plume has been obtained (pages 213–218, this issue), the problem of determining the lateral structure of the plume at a given depth is reduced to solving an eigenvalue problem of a set of ordinary differential equations with five unknown functions, with an eigenvalue being related to the thermal thickness of the plume at this depth. The lateral profiles of upward velocity, temperature and viscosity in the plume and the thickness of the plume at various depths are calculated for two sets of Newtonian rheological parameters. The calculations show that the precondition for the existence of the plume, δT/L 1 (L = the height of the plume, δT = lateral distance from the mid-plane), can be satisfied, except for the starting region of the plume or near the base of the lithosphere. At the lateral distance, δT, the upward velocity decreases to 0.1 – 50% of its maximum value at different depths. It is believed that this model may provide an approach for a quantitative description of the detailed structure of a mantle plume.

Monte Carlo Modeling of Micro Scale Gas Flows

An information preservation (IP) method has been used to simulate many micro scale gas flows. It may efficiently reduce the statistical scatter inherent in conventional particle approaches such as the direct simulation Monte Carlo (DSMC) method. This paper reviews applications of IP to some benchmark problems. Comparison of the IP results with those given by experiment, DSMC, and the linearized Boltzmann equation, as well as the Navier-Stokes equations with a slip boundary condition, and the lattice Boltzmann equation, shows that the IP method is applicable to micro scale gas flows over the entire flow regime from continuum to free molecular.

An observer-based vaccination control law for an SEIR epidemic model based on feedback linearization techniques for nonlinear systems

This paper presents a vaccination strategy for fighting against the propagation of epidemic diseases. The disease propagation is described by an SEIR (susceptible plus infected plus infectious plus removed populations) epidemic model. The model takes into account the total population amounts as a refrain for the illness transmission since its increase makes the contacts among susceptible and infected more difficult. The vaccination strategy is based on a continuous-time nonlinear control law synthesised via an exact feedback input-output linearization approach. An observer is incorporated into the control scheme to provide online estimates for the susceptible and infected populations in the case when their values are not available from online measurement but they are necessary to implement the control law. The vaccination control is generated based on the information provided by the observer. The control objective is to asymptotically eradicate the infection from the population so that the removed-by-immunity population asymptotically tracks the whole one without precise knowledge of the partial populations. The model positivity, the eradication of the infection under feedback vaccination laws and the stability properties as well as the asymptotic convergence of the estimation errors to zero as time tends to infinity are investigated.

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.

A Hybrid Scheme Based on Finite Element/Volume Methods for Two Immiscible Fluid Flows

We have successfully extended our implicit hybrid finite element/volume (FE/FV) solver to flows involving two immiscible fluids. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centered FV method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. This updating strategy can be rigorously proven to be able to eliminate the unphysical pressure boundary layer and is crucial for the correct temporal convergence rate. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centers and the auxiliary variable at vertices. The fluid interface is captured by solving an advection equation for the volume fraction of one of the fluids. The same matrix-free FV method, as the one used for momentum equations, is used to solve the advection equation. We will focus on the interface sharpening strategy to minimize the smearing of the interface over time. We have developed and implemented a global mass conservation algorithm that enforces the conservation of the mass for each fluid.

Nonlinear oscillations in discrete and continuous systems

The first thesis topic is a perturbation method for resonantly coupled nonlinear oscillators. By successive near-identity transformations of the original equations, one obtains new equations with simple structure that describe the long time evolution of the motion. This technique is related to two-timing in that secular terms are suppressed in the transformation equations. The method has some important advantages. Appropriate time scalings are generated naturally by the method, and don't need to be guessed as in two-timing. Furthermore, by continuing the procedure to higher order, one extends (formally) the time scale of valid approximation. Examples illustrate these claims. Using this method, we investigate resonance in conservative, non-conservative and time dependent problems. Each example is chosen to highlight a certain aspect of the method.

The second thesis topic concerns the coupling of nonlinear chemical oscillators. The first problem is the propagation of chemical waves of an oscillating reaction in a diffusive medium. Using two-timing, we derive a nonlinear equation that determines how spatial variations in the phase of the oscillations evolves in time. This result is the key to understanding the propagation of chemical waves. In particular, we use it to account for certain experimental observations on the Belusov-Zhabotinskii reaction.

Next, we analyse the interaction between a pair of coupled chemical oscillators. This time, we derive an equation for the phase shift, which measures how much the oscillators are out of phase. This result is the key to understanding M. Marek's and I. Stuchl's results on coupled reactor systems. In particular, our model accounts for synchronization and its bifurcation into rhythm splitting.

Finally, we analyse large systems of coupled chemical oscillators. Using a continuum approximation, we demonstrate mechanisms that cause auto-synchronization in such systems.

Phase control of higher spectral components in the presence of a static electric field

We investigate the higher spectral component generations driven by a few-cycle laser pulse in a dense medium when a static electric field is present. Our results show that, when assisted by a static electric field, the dependence of the transmitted laser spectrum on the carrier-envelope phase (CEP) is significantly increased. Continuum and distinct peaks can be achieved by controlling the CEP of the few-cycle ultrashort laser pulse. Such a strong variation is due to the fact that the presence of the static electric field modifies the waveform of the combined electric field, which further affects the spectral distribution of the generated higher spectral components.

Surface effects on spinwave resonance in thin magnetic films

Over the past few decades, ferromagnetic spinwave resonance in magnetic thin films has been used as a tool for studying the properties of magnetic materials. A full understanding of the boundary conditions at the surface of the magnetic material is extremely important. Such an understanding has been the general objective of this thesis. The approach has been to investigate various hypotheses of the surface condition and to compare the results of these models with experimental data. The conclusion is that the boundary conditions are largely due to thin surface regions with magnetic properties different from the bulk. In the calculations these regions were usually approximated by uniform surface layers; the spins were otherwise unconstrained except by the same mechanisms that exist in the bulk (i.e., no special "pinning" at the surface atomic layer is assumed). The variation of the ferromagnetic spinwave resonance spectra in YIG films with frequency, temperature, annealing, and orientation of applied field provided an excellent experimental basis for the study.

This thesis can be divided into two parts. The first part is ferromagnetic resonance theory; the second part is the comparison of calculated with experimental data in YIG films. Both are essential in understanding the conclusion that surface regions with properties different from the bulk are responsible for the resonance phenomena associated with boundary conditions.

The theoretical calculations have been made by finding the wave vectors characteristic of the magnetic fields inside the magnetic medium, and then combining the fields associated with these wave vectors in superposition to match the specified boundary conditions. In addition to magnetic boundary conditions required for the surface layer model, two phenomenological magnetic boundary conditions are discussed in detail. The wave vectors are easily found by combining the Landau-Lifshitz equations with Maxwell's equations. Mode positions are most easily predicted from the magnetic wave vectors obtained by neglecting damping, conductivity, and the displacement current. For an insulator where the driving field is nearly uniform throughout the sample, these approximations permit a simple yet accurate calculation of the mode intensities. For metal films this calculation may be inaccurate but the mode positions are still accurately described. The techniques necessary for calculating the power absorbed by the film under a specific excitation including the effects of conductivity, displacement current and damping are also presented.

In the second part of the thesis the properties of magnetic garnet materials are summarized and the properties believed associated with the two surface regions of a YIG film are presented. Finally, the experimental data and calculated data for the surface layer model and other proposed models are compared. The conclusion of this study is that the remarkable variety of spinwave spectra that arises from various preparation techniques and subsequent treatments can be explained by surface regions with magnetic properties different from the bulk.

Upper hybrid resonance absorption and scattering from a plasma column

The electromagnetic scattering and absorption properties of small (kr~1/2) inhomogeneous magnetoplasma columns are calculated via the full set of Maxwell's equations with tensor dielectric constitutive relation. The cold plasma model with collisional damping is used to describe the column. The equations are solved numerically, subject to boundary conditions appropriate to an infinite parallel strip line and to an incident plane wave. The results are similar for several density profiles and exhibit semiquantitative agreement with measurements in waveguide. The absorption is spatially limited, especially for small collision frequency, to a narrow hybrid resonant layer and is essentially zero when there is no hybrid layer in the column. The reflection is also enhanced when the hybrid layer is present, but the value of the reflection coefficient is strongly modified by the presence of the glass tube. The nature of the solutions and an extensive discussion of the conditions under which the cold collisional model should yield valid results is presented.

Effect of microcosmic optical parameters of dopants on the signal-to-noise ratio in doubly doped LiNbO3 crystals

We propose a united theory that describes the two-center recording system by taking scattering noise into account. The temporal evolution of the signal-to-noise ratio in doubly doped photorefractive crystals is described based on jointly solving material equations and coupled-wave equations with the fourth-order Runge-Kutta method. Roles of microcosmic optical parameters of dopants on the signal-to-noise ratio are discussed in detail. The theoretical results can confirm and predict experimental results. (c) 2005 Elsevier GmbH. All rights reserved.

Improvement in refractive-index change in LiNbO3 : Ce : Cu by applying an external electric field

By jointly solving two-centre material equations with a nonzero external electric field and coupled-wave equations, we have numerically studied the dependence of the non-volatile holographic recording in LiNbO3:Ce:Cu crystals on the external electric field. The dominative photovoltaic effect of the non-volatile holographic recording in doubly doped LiNbO3 crystals is directly verified. And an external electric field that is applied in the positive direction along the c-axis (or a large one in the negative direction of the c-axis) in the recording phase and another one that is applied in the negative direction of the c-axis in the fixing phase are both proved to benefit strong photorefractive performances. Experimental verifications are given with a small electric field applied externally.