The problem of the spreading of a liquid film along a solid surface: A new mathematical formulation

A new mathematical model is proposed for the spreading of a liquid film on a solid surface. The model is based on the standard lubrication approximation for gently sloping films (with the no-slip condition for the fluid at the solid surface) in the major part of the film where it is not too thin. In the remaining and relatively small regions near the contact lines it is assumed that the so-called autonomy principle holds—i.e., given the material components, the external conditions, and the velocity of the contact lines along the surface, the behavior of the fluid is identical for all films. The resulting mathematical model is formulated as a free boundary problem for the classical fourth-order equation for the film thickness. A class of self-similar solutions to this free boundary problem is considered.

Microsatellite variability and genetic distances.

We analyze the within- and between-population dynamics of the distribution of the number of repeats at multiple microsatellite DNA loci subject to stepwise mutation. Analytical expressions for moments up to the fourth order within a locus and the variance of between-locus variance at mutation-drift equilibrium have been obtained. These statistics may be used to test the appropriateness of the one-step mutation model and to detect between-locus variation in the mutation rate. Published data are compatible with the one-step mutation model, although they do not reject the two-step model. Using both multinomial sampling and diffusion approximations for the analysis of the genetic distance introduced by Goldstein et al. [Goldstein, D. B., Linares, A. R., Cavalli-Sforza, L. L. & Feldman, M. W. (1995) Proc. Natl. Acad. Sci. USA 92, 6723-6727], we show that this distance follows a chi 2 distribution with degrees of freedom equal to the number of loci when there is no variation in mutation rates among the loci. In the presence of such variation, the variance of the distance is obtained. We conclude that the number of microsatellite loci required for the construction of phylogenetic trees with reliable branch lengths may be several hundred. Also, mutations that change repeat scores by several units, even though extremely rare, may dramatically influence estimates of population parameters.

Equações de quarta ordem na modelagem de oscilações de pontes

Equações diferenciais de quarta ordem aparecem naturalmente na modelagem de oscilações de estruturas elásticas, como aquelas observadas em pontes pênseis. São considerados dois modelos que descrevem as oscilações no tabuleiro de uma ponte. No modelo unidimensional estudamos blow up em espaço finito de soluções de uma classe de equações diferenciais de quarta ordem. Os resultados apresentados solucionam uma conjectura apresentada em [F. Gazzola and R. Pavani. Wide oscillation finite time blow up for solutions to nonlinear fourth order differential equations. Arch. Ration. Mech. Anal., 207(2):717752, 2013] e implicam a não existência de ondas viajantes com baixa velocidade de propagação em uma viga. No modelo bidimensional analisamos uma equação não local para uma placa longa e fina, suportada nas extremidades menores, livre nas demais e sujeita a protensão. Provamos existência e unicidade de solução fraca e estudamos o seu comportamento assintótico sob amortecimento viscoso. Estudamos ainda a estabilidade de modos simples de oscilação, os quais são classificados como longitudinais ou torcionais.

Identification of nonlinear systems by correlation using pseudorandom signals

This thesis is concerned with the measurement of the characteristics of nonlinear systems by crosscorrelation, using pseudorandom input signals based on m sequences. The systems are characterised by Volterra series, and analytical expressions relating the rth order Volterra kernel to r-dimensional crosscorrelation measurements are derived. It is shown that the two-dimensional crosscorrelation measurements are related to the corresponding second order kernel values by a set of equations which may be structured into a number of independent subsets. The m sequence properties determine how the maximum order of the subsets for off-diagonal values is related to the upper bound of the arguments for nonzero kernel values. The upper bound of the arguments is used as a performance index, and the performance of antisymmetric pseudorandom binary, ternary and quinary signals is investigated. The performance indices obtained above are small in relation to the periods of the corresponding signals. To achieve higher performance with ternary signals, a method is proposed for combining the estimates of the second order kernel values so that the effects of some of the undesirable nonzero values in the fourth order autocorrelation function of the input signal are removed. The identification of the dynamics of two-input, single-output systems with multiplicative nonlinearity is investigated. It is shown that the characteristics of such a system may be determined by crosscorrelation experiments using phase-shifted versions of a common signal as inputs. The effects of nonlinearities on the estimates of system weighting functions obtained by crosscorrelation are also investigated. Results obtained by correlation testing of an industrial process are presented, and the differences between theoretical and experimental results discussed for this case;

Thin film bearings with phase change of lubricant

The possible evaporation of lubricant in fluid film bearings has been investigated theoretically and by experiment using a radial flow hydrostatic bearing supplied with liquid refrigerant R114. Good correlation between measured and theoretical values was obtained using a bespoke computational fluid dynamic model in which the flow was assumed to be laminar and adiabatic. The effects of viscous dissipation and vapour generation within the fluid film are fully accounted for by applying a fourth order Runge-Kutta routine to satisfy the radial and filmwise transverse constraints of momentum, energy and mass conservation. The results indicate that the radial velocity profile remains parabolic while the flow remains in the liquid phase and that the radial rate of enthalpy generation is then constant across the film at a given radius. The results also show that evaporation will commence at a radial location determined by geometry and flow conditions and in fluid layers adjacent to the solid boundaries. Evaporation is shown to progress in the radial direction and the load carrying capacity of such a bearing is reduced significantly. Expressions for the viscosity of the liquid/vapour mixture found in the literature survey have not been tested against experimental data. A new formulation is proposed in which the suitable choice of a characteristic constant yields close representation to any of these expressions. Operating constraints imposed by the design of the experimental apparatus limited the extent of the surface over which evaporation could be obtained, and prevented clear identification of the most suitable relationship for the viscosity of the liquid/vapour mixture. The theoretical model was extended to examine the development of two phase flow in a rotating shaft face seal of uniform thickness. Previous theoretical analyses have been based on the assumption that the radial velocity profile of the flow is always parabolic, and that the tangential component of velocity varies linearly from the value at the rotating surface, to zero at the stationary surface. The computational fluid dynamic analysis shows that viscous shear and dissipation in the fluid adjacent to the rotating surface leads to developing evaporation with a consequent reduction in tangential shear forces. The tangential velocity profile is predicted to decay rapidly through the film, exhibiting a profile entirely different to that assumed by previous investigators. Progressive evaporation takes place close to the moving wall and does not occur completely at a single radial location, as has been claimed in earlier work.

Line-by-line fiber Bragg grating made by femtosecond laser

In this letter, we report on the inscription of a fourth-order fiber Bragg grating made line-by-line in the optical fiber using a femtosecond laser. Strong Bragg resonance (~17 dB) and low insertion loss (~0.5 dB) were obtained with only 2000 periods. Measured refractive index change of these inscribed lines reaches up to 7 × 10-3. The grating was fully characterized and the low insertion loss together with low polarization-dependent loss were realized compared to gratings made by the point-by-point method. The high temperature annealing experiment shows the grating can survive up to at least 800°C.

New developments in the study of optical parabolic pulses in normally dispersive fibers

We report two recent studies dealing with the evolution of parabolic pulses in normally dispersive fibres. On the one hand, the nonlinear reshaping from a Gaussian intensity profile towards the asymptotic parabolic shape is experimentally investigated in a Raman amplifier. On the other hand, the significant impact of the fourth order dispersion on a passive propagation is theoretically discussed: we numerically demonstrate flat-top, coherent supercontinuum generation in an all-normal dispersion-flattened photonic crystal fiber. This shape is associated to a strong reshaping of the temporal profile what becomes triangular.

Quasi-separation of the biharmonic partial differential equation

In this paper, we consider analytical and numerical solutions to the Dirichlet boundary-value problem for the biharmonic partial differential equation on a disc of finite radius in the plane. The physical interpretation of these solutions is that of the harmonic oscillations of a thin, clamped plate. For the linear, fourth-order, biharmonic partial differential equation in the plane, it is well known that the solution method of separation in polar coordinates is not possible, in general. However, in this paper, for circular domains in the plane, it is shown that a method, here called quasi-separation of variables, does lead to solutions of the partial differential equation. These solutions are products of solutions of two ordinary linear differential equations: a fourth-order radial equation and a second-order angular differential equation. To be expected, without complete separation of the polar variables, there is some restriction on the range of these solutions in comparison with the corresponding separated solutions of the second-order harmonic differential equation in the plane. Notwithstanding these restrictions, the quasi-separation method leads to solutions of the Dirichlet boundary-value problem on a disc with centre at the origin, with boundary conditions determined by the solution and its inward drawn normal taking the value 0 on the edge of the disc. One significant feature for these biharmonic boundary-value problems, in general, follows from the form of the biharmonic differential expression when represented in polar coordinates. In this form, the differential expression has a singularity at the origin, in the radial variable. This singularity translates to a singularity at the origin of the fourth-order radial separated equation; this singularity necessitates the application of a third boundary condition in order to determine a self-adjoint solution to the Dirichlet boundary-value problem. The penultimate section of the paper reports on numerical solutions to the Dirichlet boundary-value problem; these results are also presented graphically. Two specific cases are studied in detail and numerical values of the eigenvalues are compared with the results obtained in earlier studies.

New developments in the study of optical parabolic pulses in normally dispersive fibers

We report two recent studies dealing with the evolution of parabolic pulses in normally dispersive fibres. On the one hand, the nonlinear reshaping from a Gaussian intensity profile towards the asymptotic parabolic shape is experimentally investigated in a Raman amplifier. On the other hand, the significant impact of the fourth order dispersion on a passive propagation is theoretically discussed: we numerically demonstrate flat-top, coherent supercontinuum generation in an all-normal dispersion-flattened photonic crystal fiber. This shape is associated to a strong reshaping of the temporal profile what becomes triangular.

Supersymmetry and Ghosts in Quantum Mechanics

2000 Mathematics Subject Classification: 81Q60, 35Q40.

Análise da biorremediação de compostos monoaromáticos em água através da pseudomonas aeruginosa

The monoaromatic compounds are toxic substances present in petroleum derivades and used broadly in the chemical and petrochemical industries. Those compounds are continuously released into the environment, contaminating the soil and water sources, leading to the possible unfeasibility of those hydrous resources due to their highly carcinogenic and mutagenic potentiality, since even in low concentrations, the BTEX may cause serious health issues. Therefore, it is extremely important to develop and search for new methodologies that assist and enable the treatment of BTEX-contaminated matrix. The bioremediation consists on the utilization of microbial groups capable of degrading hydrocarbons, promoting mineralization, or in other words, the permanent destruction of residues, eliminating the risks of future contaminations. This work investigated the biodegradation kinetics of water-soluble monoaromatic compounds (benzene, toluene and ethylbenzene), based on the evaluation of its consummation by the Pseudomonas aeruginosa bacteria, for concentrations varying from 40 to 200 mg/L. To do so, the performances of Monod kinetic model for microbial growth were evaluated and the material balance equations for a batch operation were discretized and numerically solved by the fourth order Runge-Kutta method. The kinetic parameters obtained using the method of least squares as statistical criteria were coherent when compared to those obtained from the literature. They also showed that, the microorganism has greater affinity for ethylbenzene. That way, it was possible to observe that Monod model can predict the experimental data for the individual biodegradation of the BTEX substrates and it can be applied to the optimization of the biodegradation processes of toxic compounds for different types of bioreactors and for different operational conditions.

Análise da biorremediação de compostos monoaromáticos em água através da pseudomonas aeruginosa

The monoaromatic compounds are toxic substances present in petroleum derivades and used broadly in the chemical and petrochemical industries. Those compounds are continuously released into the environment, contaminating the soil and water sources, leading to the possible unfeasibility of those hydrous resources due to their highly carcinogenic and mutagenic potentiality, since even in low concentrations, the BTEX may cause serious health issues. Therefore, it is extremely important to develop and search for new methodologies that assist and enable the treatment of BTEX-contaminated matrix. The bioremediation consists on the utilization of microbial groups capable of degrading hydrocarbons, promoting mineralization, or in other words, the permanent destruction of residues, eliminating the risks of future contaminations. This work investigated the biodegradation kinetics of water-soluble monoaromatic compounds (benzene, toluene and ethylbenzene), based on the evaluation of its consummation by the Pseudomonas aeruginosa bacteria, for concentrations varying from 40 to 200 mg/L. To do so, the performances of Monod kinetic model for microbial growth were evaluated and the material balance equations for a batch operation were discretized and numerically solved by the fourth order Runge-Kutta method. The kinetic parameters obtained using the method of least squares as statistical criteria were coherent when compared to those obtained from the literature. They also showed that, the microorganism has greater affinity for ethylbenzene. That way, it was possible to observe that Monod model can predict the experimental data for the individual biodegradation of the BTEX substrates and it can be applied to the optimization of the biodegradation processes of toxic compounds for different types of bioreactors and for different operational conditions.

A New Method for Modeling Free Surface Flows and Fluid-structure Interaction with Ocean Applications

The computational modeling of ocean waves and ocean-faring devices poses numerous challenges. Among these are the need to stably and accurately represent both the fluid-fluid interface between water and air as well as the fluid-structure interfaces arising between solid devices and one or more fluids. As techniques are developed to stably and accurately balance the interactions between fluid and structural solvers at these boundaries, a similarly pressing challenge is the development of algorithms that are massively scalable and capable of performing large-scale three-dimensional simulations on reasonable time scales. This dissertation introduces two separate methods for approaching this problem, with the first focusing on the development of sophisticated fluid-fluid interface representations and the second focusing primarily on scalability and extensibility to higher-order methods.

We begin by introducing the narrow-band gradient-augmented level set method (GALSM) for incompressible multiphase Navier-Stokes flow. This is the first use of the high-order GALSM for a fluid flow application, and its reliability and accuracy in modeling ocean environments is tested extensively. The method demonstrates numerous advantages over the traditional level set method, among these a heightened conservation of fluid volume and the representation of subgrid structures.

Next, we present a finite-volume algorithm for solving the incompressible Euler equations in two and three dimensions in the presence of a flow-driven free surface and a dynamic rigid body. In this development, the chief concerns are efficiency, scalability, and extensibility (to higher-order and truly conservative methods). These priorities informed a number of important choices: The air phase is substituted by a pressure boundary condition in order to greatly reduce the size of the computational domain, a cut-cell finite-volume approach is chosen in order to minimize fluid volume loss and open the door to higher-order methods, and adaptive mesh refinement (AMR) is employed to focus computational effort and make large-scale 3D simulations possible. This algorithm is shown to produce robust and accurate results that are well-suited for the study of ocean waves and the development of wave energy conversion (WEC) devices.

Ultra Compact Inline E-Plane Waveguide Bandpass Filters Using Cross Coupling

This paper presents novel ultra-compact waveguide bandpass filters that exhibit pseudo elliptic responses with ability to place transmission zeros on both sides of the passband to form sharp roll offs. The filters contain E plane extracted pole sections cascaded with cross-coupled filtering blocks. Compactness is achieved by the use of evanescent mode sections and closer arranged resonators modified to shrink in size. The filters containing non-resonating nodes are designed by means of the generalized coupling coefficients (GCC) extraction procedure for the cross-coupled filtering blocks and extracted pole sections. We illustrate the performance of the proposed structures through the design examples of a third and a fourth order filters with center frequencies of 9.2 GHz and 10 GHz respectively. The sizes of the proposed structures suitable for fabricating using the low cost E plane waveguide technology are 38% smaller than ones of the E plane extracted pole filter of the same order.

Design and implementation of generalized topologies of time-interleaved variable bandpass Σ−Δ modulators

In this thesis, novel analog-to-digital and digital-to-analog generalized time-interleaved variable bandpass sigma-delta modulators are designed, analysed, evaluated and implemented that are suitable for high performance data conversion for a broad-spectrum of applications. These generalized time-interleaved variable bandpass sigma-delta modulators can perform noise-shaping for any centre frequency from DC to Nyquist. The proposed topologies are well-suited for Butterworth, Chebyshev, inverse-Chebyshev and elliptical filters, where designers have the flexibility of specifying the centre frequency, bandwidth as well as the passband and stopband attenuation parameters. The application of the time-interleaving approach, in combination with these bandpass loop-filters, not only overcomes the limitations that are associated with conventional and mid-band resonator-based bandpass sigma-delta modulators, but also offers an elegant means to increase the conversion bandwidth, thereby relaxing the need to use faster or higher-order sigma-delta modulators. A step-by-step design technique has been developed for the design of time-interleaved variable bandpass sigma-delta modulators. Using this technique, an assortment of lower- and higher-order single- and multi-path generalized A/D variable bandpass sigma-delta modulators were designed, evaluated and compared in terms of their signal-to-noise ratios, hardware complexity, stability, tonality and sensitivity for ideal and non-ideal topologies. Extensive behavioural-level simulations verified that one of the proposed topologies not only used fewer coefficients but also exhibited greater robustness to non-idealties. Furthermore, second-, fourth- and sixth-order single- and multi-path digital variable bandpass digital sigma-delta modulators are designed using this technique. The mathematical modelling and evaluation of tones caused by the finite wordlengths of these digital multi-path sigmadelta modulators, when excited by sinusoidal input signals, are also derived from first principles and verified using simulation and experimental results. The fourth-order digital variable-band sigma-delta modulator topologies are implemented in VHDL and synthesized on Xilinx® SpartanTM-3 Development Kit using fixed-point arithmetic. Circuit outputs were taken via RS232 connection provided on the FPGA board and evaluated using MATLAB routines developed by the author. These routines included the decimation process as well. The experiments undertaken by the author further validated the design methodology presented in the work. In addition, a novel tunable and reconfigurable second-order variable bandpass sigma-delta modulator has been designed and evaluated at the behavioural-level. This topology offers a flexible set of choices for designers and can operate either in single- or dual-mode enabling multi-band implementations on a single digital variable bandpass sigma-delta modulator. This work is also supported by a novel user-friendly design and evaluation tool that has been developed in MATLAB/Simulink that can speed-up the design, evaluation and comparison of analog and digital single-stage and time-interleaved variable bandpass sigma-delta modulators. This tool enables the user to specify the conversion type, topology, loop-filter type, path number and oversampling ratio.