Transparent boundary conditions for wave propagation on unbounded domains

The numerical solution of the time dependent wave equation in an unbounded domain generally leads to a truncation of this domain, which requires the introduction of an artificial boundary with associated boundary conditions. Such nonreflecting conditions ensure the equivalence between the solution of the original problem in the unbounded region and the solution inside the artificial boundary. We consider the acoustic wave equation and derive exact transparent boundary conditions that are local in time and can be directly used in explicit methods. These conditions annihilate wave harmonics up to a given order on a spherical artificial boundary, and we show how to combine the derived boundary condition with a finite difference method. The analysis is complemented by a numerical example in two spatial dimensions that illustrates the usefulness and accuracy of transparent boundary conditions.

Fluctuations and single molecules

The fluorescence of single molecules coupled to a thermal bath is studied both experimentally and theoretically. The effect of different fluctuations on the coherence properties of resonance fluorescence is considered first. Coherence is measured in an interference experiment where a single molecule is used as a light source. A standard approach based on the optical Bloch equations apparently provides quite an accurate description of the interference experiment. Systems with long correlation times (where spectra are time dependent on any timescale) are considered next. It is shown that intensity-time-frequency correlation spectroscopy, which provides both high signal-to-noise ratio and high time resolution, is very suitable for such a case. The Bloch equations are further tested in an experiment where the shape of an excitation spectral line of a single molecule is accurately measured over six orders of magnitude of the exciting laser power. Significant deviations from the predictions of the Bloch equations are found. The role of critical parameters-the correlation time of the bath, the Rabi oscillation period, and the coupling constant between the bath and the molecule-is discussed. The paper also includes a short general introduction to the methodology of single-molecule studies.

A time-harmonic target-field method for designing unshielded RF coils in MRI

Time-harmonic methods are required in the accurate design of RF coils as operating frequency increases. This paper presents such a method to find a current density solution on the coil that will induce some desired magnetic field upon an asymmetrically located target region within. This inverse method appropriately considers the geometry of the coil via a Fourier series expansion, and incorporates some new regularization penalty functions in the solution process. A new technique is introduced by which the complex, time-dependent current density solution is approximated by a static coil winding pattern. Several winding pattern solutions are given, with more complex winding patterns corresponding to more desirable induced magnetic fields.

DOSY NMR studies of chemical exchange behavior of poly(2-hydroxyethyl methacrylate)

PFG-NMR was used to study the chemical exchange of linear PHEMA having a range of molecular weights with water in DMSO containing varying quantities of water. The aim was to investigate the use of PFG-NMR to study chemical exchange between a polymer with exchangeable protons and a small fast diffusing molecule to provide insight into the conformation adopted by a polymer in solution. The experimental data were simulated closely for the two-site exchange case using the Bloch equations modified for chemical exchange and diffusion. The exchange rate could be used to detect changes in polymer conformation resulting from changes in the solvent. PHEMA of MW 10 000 showed significant time-dependent changes in exchange rate, resulting from preferential solvation of the OH sites by water, and subsequent conformational changes which altered accessibility of the OH sites to water. This behavior was not observed for larger MW PHEMA, which adopted a stable conformation immediately. Large changes in the exchange rate were not reflected in changes to the hydrodynamic radius, suggesting that a minimal overall change in the chain dimensions occurred. DMSO was found to be a poor solvent for PHEMA, which adopts a compact conformation in DMSO. This work has demonstrated that PFG-NMR is a sensitive method for detecting subtle changes in polymer conformation in polymers with exchangeable protons.

A basis function approach to Bayesian inference in diffusion processes

In this paper, we present a framework for Bayesian inference in continuous-time diffusion processes. The new method is directly related to the recently proposed variational Gaussian Process approximation (VGPA) approach to Bayesian smoothing of partially observed diffusions. By adopting a basis function expansion (BF-VGPA), both the time-dependent control parameters of the approximate GP process and its moment equations are projected onto a lower-dimensional subspace. This allows us both to reduce the computational complexity and to eliminate the time discretisation used in the previous algorithm. The new algorithm is tested on an Ornstein-Uhlenbeck process. Our preliminary results show that BF-VGPA algorithm provides a reasonably accurate state estimation using a small number of basis functions.

Theoretical studies of the mid-latitude ionosphere

Following a brief description of the atmosphere and ionosphere in Chapter I we describe how the equations of continuity and momentum for 0+, H+, He+, 0++ are derived from the formulations of St-Maurice and Schunk(1977) and Quegan et al.(1981) in Chapter II. In Chapter III we investigate the nature of the downward flow of protons in a collapsing post-sunset ionosphere. We derive an analytical form for the limiting temperature, we also note the importance of the polarization field term and concluded that the flow will remain subsonic for realistic conditions. The time-dependent behaviour of He+ under sunspot minimum conditions is investigated in Chapter IV. This is achieved by numerical solution of the 0+, H+ and,He+ continuity and momentum equations, treating He+ as a minor ion with 0+ , H+ as major ions. We found that He+ flows upwards during the day-time and downwards during the nighttime. He+ flux tube content reached a maximum on the 8th day of the integration period and started to decreasing. This is due to the large amount of H+ present at the late stages of the integration period which makes He+ unable to diffuse through the H+ layer away from the loss region. In Chapter V we investigate the behaviour of 0++ using sunspot maximum parameters. Although our results support the findings of Geis and Young (1981) that the large amounts of 0++ at the equator are caused mainly by thermal diffusion, the model used by Geis and Young overemphesizes the effect of thermal diffusion. The importance of 0++ - 0+ collision frequency is also noted. In Chapter VI we extend the work of Chapter IV, presenting a comparative study of H and He at sunspot minimum and sunspot maximum.In this last Chapter all three ions, O+ ,H+ and He+ , are treated theoretically as major ions and we concentrate mainly on light ion contents and fluxes. The results of this Chapter indicate that by assuming He+ as a minor ion we under-estimate He+ and over-estimate. H+. Some interesting features concerning the day to day behaviour of the light ion fluxes arise. In particular the day-time H+ fluxes decrease from day to day in contrast to the work of Murphy et al.(1976). In appendix.A we derive some analytical forms for the optical depth so that the models can include a realistic description of photoionization.

Time-resolved Q-factor measurement and its application in performance analysis of 42.6 Gbit/s packets generated by SGDBR lasers

We demonstrate a novel time-resolved Q-factor measurement technique and demonstrate its application in the analysis of optical packet switching systems with high information spectral density. For the first time, we report the time-resolved Q-factor measurement of 42.6 Gbit/s AM-PSK and DQPSK modulated packets, which were generated by a SGDBR laser under wavelength switching. The time dependent degradation of Q-factor performance during the switching transient was analyzed and was found to be correlated with different laser switching characteristics in each case.

Benefits of retailer-supplier partnership initiatives under time-varying demand:a comparative analytical study

This paper aims to help supply chain managers to determine the value of retailer-supplier partnership initiatives beyond information sharing (IS) according to their specific business environment under time-varying demand conditions. For this purpose, we use integer linear programming models to quantify the benefits that can be accrued by a retailer, a supplier and system as a whole from shift in inventory ownership and shift in decision-making power with that of IS. The results of a detailed numerical study pertaining to static time horizon reveal that the shift in inventory ownership provides system-wide cost benefits in specific settings. Particularly, when it induces the retailer to order larger quantities and the supplier also prefers such orders due to significantly high setup and shipment costs. We observe that the relative benefits of shift in decision-making power are always higher than the shift in inventory ownership under all the conditions. The value of the shift in decision-making power is greater than IS particularly when the variability of underlying demand is low and time-dependent variation in production cost is high. However, when the shipment cost is negligible and order issuing efficiency of the supplier is low, the cost benefits of shift in decision-making power beyond IS are not significant. © 2012 Taylor & Francis.

Time-spatial drift of decelerating electromagnetic pulses

A time dependent electromagnetic pulse generated by a current running laterally to the direction of the pulse propagation is considered in paraxial approximation. It is shown that the pulse envelope moves in the time-spatial coordinates on the surface of a parabolic cylinder for the Airy pulse and a hyperbolic cylinder for the Gaussian. These pulses propagate in time with deceleration along the dominant propagation direction and drift uniformly in the lateral direction. The Airy pulse stops at infinity while the asymptotic velocity of the Gaussian is nonzero. © 2013 Optical Society of America.

Local Energy Decay in Even Dimensions for the Wave Equation with a Time-Periodic Non-Trapping Metric and Applications to Strichartz Estimates

2000 Mathematics Subject Classification: 35B40, 35L15.

Activated escape from polarization pulling in fibre Raman amplifiers

Fluctuation-induced escape (FIE) from a metastable state with probability controlled by external force is a process inherent in many physical phenomena such as diffusion in crystals, protein folding, activated chemical reactions etc. [1-3]. In this work we present a novel example of FIE problem, considering a very practical nonlinear system recently emerged in the area of fibre telecommunications. Unlike the standard FIE problems where noise is time-dependent, in fibre Raman amplifier (FRA) the role of noise is played by frozen fluctuations of parameters (random birefringence) along the fibre span which result from the breaking of cylindrical symmetry during the fibre drawing [4-6]. The role of periodic forcing in this problem is played by the periodic fibre spinning, leading to key model that is formally similar to the time-domain equations for periodically forced escape [1-3]. © 2011 IEEE.

Strongly localized moving discrete dissipative breather-solitons in Kerr nonlinear media supported by intrinsic gain

We investigate the mobility of nonlinear localized modes in a generalized discrete Ginzburg-Landau-type model, describing a one-dimensional waveguide array in an active Kerr medium with intrinsic, saturable gain and damping. It is shown that exponentially localized, traveling discrete dissipative breather-solitons may exist as stable attractors supported only by intrinsic properties of the medium, i.e., in the absence of any external field or symmetry-breaking perturbations. Through an interplay by the gain and damping effects, the moving soliton may overcome the Peierls-Nabarro barrier, present in the corresponding conservative system, by self-induced time-periodic oscillations of its power (norm) and energy (Hamiltonian), yielding exponential decays to zero with different rates in the forward and backward directions. In certain parameter windows, bistability appears between fast modes with small oscillations and slower, large-oscillation modes. The velocities and the oscillation periods are typically related by lattice commensurability and exhibit period-doubling bifurcations to chaotically "walking" modes under parameter variations. If the model is augmented by intersite Kerr nonlinearity, thereby reducing the Peierls-Nabarro barrier of the conservative system, the existence regime for moving solitons increases considerably, and a richer scenario appears including Hopf bifurcations to incommensurately moving solutions and phase-locking intervals. Stable moving breathers also survive in the presence of weak disorder. © 2014 American Physical Society.

Dissipative parametric modulation instability and pattern formation in nonlinear optical systems

We present the essential features of the dissipative parametric instability, in the universal complex Ginzburg- Landau equation. Dissipative parametric instability is excited through a parametric modulation of frequency dependent losses in a zig-zag fashion in the spectral domain. Such damping is introduced respectively for spectral components in the +ΔF and in the -ΔF region in alternating fashion, where F can represent wavenumber or temporal frequency depending on the applications. Such a spectral modulation can destabilize the homogeneous stationary solution of the system leading to growth of spectral sidebands and to the consequent pattern formation: both stable and unstable patterns in one- and in two-dimensional systems can be excited. The dissipative parametric instability provides an useful and interesting tool for the control of pattern formation in nonlinear optical systems with potentially interesting applications in technological applications, like the design of mode- locked lasers emitting pulse trains with tunable repetition rate; but it could also find realizations in nanophotonics circuits or in dissipative polaritonic Bose-Einstein condensates.

Mechanistic Models of Anti-HIV Microbicide Drug Delivery

A new modality for preventing HIV transmission is emerging in the form of topical microbicides. Some clinical trials have shown some promising results of these methods of protection while other trials have failed to show efficacy. Due to the relatively novel nature of microbicide drug transport, a rigorous, deterministic analysis of that transport can help improve the design of microbicide vehicles and understand results from clinical trials. This type of analysis can aid microbicide product design by helping understand and organize the determinants of drug transport and the potential efficacies of candidate microbicide products.

Microbicide drug transport is modeled as a diffusion process with convection and reaction effects in appropriate compartments. This is applied here to vaginal gels and rings and a rectal enema, all delivering the microbicide drug Tenofovir. Although the focus here is on Tenofovir, the methods established in this dissertation can readily be adapted to other drugs, given knowledge of their physical and chemical properties, such as the diffusion coefficient, partition coefficient, and reaction kinetics. Other dosage forms such as tablets and fiber meshes can also be modeled using the perspective and methods developed here.

The analyses here include convective details of intravaginal flows by both ambient fluid and spreading gels with different rheological properties and applied volumes. These are input to the overall conservation equations for drug mass transport in different compartments. The results are Tenofovir concentration distributions in time and space for a variety of microbicide products and conditions. The Tenofovir concentrations in the vaginal and rectal mucosal stroma are converted, via a coupled reaction equation, to concentrations of Tenofovir diphosphate, which is the active form of the drug that functions as a reverse transcriptase inhibitor against HIV. Key model outputs are related to concentrations measured in experimental pharmacokinetic (PK) studies, e.g. concentrations in biopsies and blood. A new measure of microbicide prophylactic functionality, the Percent Protected, is calculated. This is the time dependent volume of the entire stroma (and thus fraction of host cells therein) in which Tenofovir diphosphate concentrations equal or exceed a target prophylactic value, e.g. an EC50.

Results show the prophylactic potentials of the studied microbicide vehicles against HIV infections. Key design parameters for each are addressed in application of the models. For a vaginal gel, fast spreading at small volume is more effective than slower spreading at high volume. Vaginal rings are shown to be most effective if inserted and retained as close to the fornix as possible. Because of the long half-life of Tenofovir diphosphate, temporary removal of the vaginal ring (after achieving steady state) for up to 24h does not appreciably diminish Percent Protected. However, full steady state (for the entire stromal volume) is not achieved until several days after ring insertion. Delivery of Tenofovir to the rectal mucosa by an enema is dominated by surface area of coated mucosa and whether the interiors of rectal crypts are filled with the enema fluid. For the enema 100% Percent Protected is achieved much more rapidly than for vaginal products, primarily because of the much thinner epithelial layer of the mucosa. For example, 100% Percent Protected can be achieved with a one minute enema application, and 15 minute wait time.

Results of these models have good agreement with experimental pharmacokinetic data, in animals and clinical trials. They also improve upon traditional, empirical PK modeling, and this is illustrated here. Our deterministic approach can inform design of sampling in clinical trials by indicating time periods during which significant changes in drug concentrations occur in different compartments. More fundamentally, the work here helps delineate the determinants of microbicide drug delivery. This information can be the key to improved, rational design of microbicide products and their dosage regimens.

Time-dependency of Rocks and Implications Associated with Tunnelling

The formulation of a geotechnical model and the associated prediction of the mechanical behaviour is a challenge engineers need to overcome in order to optimize tunnel design and meet project requirements. Special challenges arise in cases where rocks and rockmasses are susceptible to time-effects and time-dependent processes govern. Progressive rockmass deformation and instability, time-dependent overloading of support and delayed failures are commonly the result of time-dependent phenomena. The research work presented in this thesis serves as an attempt to provide more insight into the time-dependent behaviour of rocks. Emphasis is given on investigating and analyzing creep deformation and time-dependent stress relaxation phenomenon at the laboratory scale and in-depth analyses are presented. This thesis further develops the understanding of these phenomena and practical yet scientific tools for estimating and predicting the long-term strength and the maximum stress relaxation of rock materials are proposed. The identification of the existence of three distinct behavioural stages during stress relaxation is presented and discussed. The main observations associated with time-dependent behaviour are employed in numerical analyses and applied at the tunnel scale. A new approach for simulating and capturing the time-dependent behaviour coupled with the tunnel advancement effect is also developed and analyzed. Guidance is provided to increase the understanding of the support-rockmass interaction and the main implications and significance of time-dependent behaviour associated with rock tunnelling are discussed. The work presented in this thesis advances the scientific understanding of time-dependent rock and rockmass behaviour, increases the awareness of how such phenomena are captured numerically, and lays out a framework for dealing with such deformations when predicting tunnel deformations. Practical aspects of this thesis are also presented, which will increase their usage in the associated industries and close the gap between the scientific and industry communities.