Theoretical investigation of turbulent boundary layer over a wavy surface

The important features of the two-dimensional incompressible turbulent flow over a wavy surface of wavelength comparable with the boundary layer thickness are analyzed.

A turbulent field method using model equation for turbulent shear stress similar to the scheme of Bradshaw, Ferriss and Atwell (1967) is employed with suitable modification to cover the viscous sublayer. The governing differential equations are linearized based on the small but finite amplitude to wavelength ratio. An orthogonal wavy coordinate system, accurate to the second order in the amplitude ratio, is adopted to avoid the severe restriction to the validity of linearization due to the large mean velocity gradient near the wall. Analytic solution up to the second order is obtained by using the method of matched-asymptotic-expansion based on the large Reynolds number and hence the small skin friction coefficient.

In the outer part of the layer, the perturbed flow is practically "inviscid." Solutions for the velocity, Reynolds stress and also the wall pressure distributions agree well with the experimental measurement. In the wall region where the perturbed Reynolds stress plays an important role in the process of momentum transport, only a qualitative agreement is obtained. The results also show that the nonlinear second-order effect is negligible for amplitude ratio of 0.03. The discrepancies in the detailed structure of the velocity, shear stress, and skin friction distributions near the wall suggest modifications to the model are required to describe the present problem.

Energy inequalities and error estimates for axisymmetric torsion of thin elastic shells of revolution

The problem motivating this investigation is that of pure axisymmetric torsion of an elastic shell of revolution. The analysis is carried out within the framework of the three-dimensional linear theory of elastic equilibrium for homogeneous, isotropic solids. The objective is the rigorous estimation of errors involved in the use of approximations based on thin shell theory.

The underlying boundary value problem is one of Neumann type for a second order elliptic operator. A systematic procedure for constructing pointwise estimates for the solution and its first derivatives is given for a general class of second-order elliptic boundary-value problems which includes the torsion problem as a special case.

The method used here rests on the construction of “energy inequalities” and on the subsequent deduction of pointwise estimates from the energy inequalities. This method removes certain drawbacks characteristic of pointwise estimates derived in some investigations of related areas.

Special interest is directed towards thin shells of constant thickness. The method enables us to estimate the error involved in a stress analysis in which the exact solution is replaced by an approximate one, and thus provides us with a means of assessing the quality of approximate solutions for axisymmetric torsion of thin shells.

Finally, the results of the present study are applied to the stress analysis of a circular cylindrical shell, and the quality of stress estimates derived here and those from a previous related publication are discussed.

Biossorção de tório com emprego de Sargassum filipendula

Com o aumento do rigor para descarte de efluentes contaminado com metais pesados, as pesquisas têm se intensificado na busca de métodos de remoção, que tragam bons resultados de captação dos metais, aliado a um baixo custo. O uso de biomassas como bactérias, fungos e algas marinhas como material adsorvente, tem se apresentado como uma alternativa, principalmente para soluções com baixo teor de metais. Neste trabalho, a alga marinha Sargassum filipendula foi avaliada na sua capacidade de remoção do metal tório de uma solução sintética e do efluente dos laboratórios de análises ambientais do IRD. A cinética desta biossorção foi estudada em regime de batelada e o equilíbrio foi alcançado com 180 min de reação. Dois modelos cinéticos foram utilizados nesta avaliação, um de primeira ordem e um modelo de pseudo segunda ordem, tendo o modelo de segunda ordem apresentado um melhor ajuste dos dados. Na avaliação da capacidade máxima de captação do tório pela biomassa marinha em regime de batelada, foi construída a isoterma que apresentou um perfil crescente na captação alcançando um máximo de 2,59 mol/g. Os modelos de Langmuir e Freundlich foram utilizados para ajustar os dados da isoterma, tendo apresentado maior correlação com os dados o modelo de Langmuir, resultando num valor de captação máxima calculado pelo modelo de 2,92 mol/g. A capacidade de remoção do metal tório da alga Sargassum filipendula também foi avaliada em regime contínuo. Um estudo de altura crítica de leito foi realizado preenchendo-se uma coluna com diferentes massas de alga correspondendo a diferentes alturas de leito. A concentração de tório residual foi quantificada na solução de saída e a menor concentração na saída foi alcançada com 40 cm de leito ou 96 g de biomassa. Após este estudo um sistema contínuo com duas colunas, com 96 g de biomassa cada, e uma bomba peristáltica foi utilizado para o tratamento do efluente real do IRD, contendo não apenas o metal tório, mas outros metais como urânio, cálcio, cromo, ferro, chumbo, etc. Cento e cinco litros de efluente foram tratados numa concentração de 6 mg/L, a concentração do efluente de saída foi de 3,75 mg/L. A caracterização do efluente quanto aos metais presentes foi realizada em ICP-MS, os resultados demonstraram que não houve competição pelos sítios ligantes da biomassa entre o tório e os demais metais. Além disso, alguns metais como cálcio, ferro e magnésio, tiveram um aumento na concentração de saída indicando a presença do mecanismo de troca iônica na biossorção do tório por Sargassum filipendula

Analysis of dispersion compensation by chirped fiber grating in directly and externally modulated CATV systems

The dispersion compensation effect of the chirped fiber grating (CFG) is analyzed theoretically, and analytic expressions are derived for composite second-order (CSO) distortion in analog modulated sub-carrier multiplexed (AM-SCM) cable television (CATV) systems with externally and directly modulated transmitters. Simulations are given for the two kinds of modulations and for standard single mode fiber and non-zero dispersion shift fiber (NZDSF) systems. The results show that CFG could be used as a dispersion compensator in directly modulated systems, but its dispersion coefficient should be adjusted much more precisely than the externally modulated system. The requirements for the NZDSF system could be loosened much. It is proposed that directly modulated source may be used as a transmitter in CATV systems combined with tunable CFG dispersion compensator being adjusted precisely, which may be more cost-effective than externally modulation technology. (c) 2006 Elsevier GmbH. All rights reserved.

Synthesis and characterization of photorefractive polymer based on chemically hybridized CdS-PVK nanocomposite with a new azo chromophore

We report an organic/inorganic polymer composite based on the chemically hybridized photoconductor CdS-PVK nanocomposite doped with a new second-order optically nonlinear chromophore 1-n-butoxy-2-methyl-(4-p-nitrophenylazo)benzene (BMNPAB) and plasticizer 9-ethyl-carbazole (ECZ) to manifest a photorefractive (PR) effect. A detailed description of the synthesis and characterization of BMNPAB is presented. The poled film including PVK-10-CdS nanocomposite and BMNPAB exhibits a high second harmonic generation (SHG) coefficient of 31 pm/V The photoconductivity of PVK-CdS nanocomposite also was studied here. Two-beam coupling experiment clearly indicated an asymmetric optical energy exchange between two beams on the polymer composite at zero electrical field, and the two-beam coupling gain of 50.0 cm(-1) and diffraction efficiencv of 4.2% were obtained at 647.1 nm wavelength. (C) 2007 Elsevier Ltd. All rights reserved.

A kinetic theory description for external spherical flows with arbitrary Knudsen number by a moment method

The Maxwell integral equations of transfer are applied to a series of problems involving flows of arbitrary density gases about spheres. As suggested by Lees a two sided Maxwellian-like weighting function containing a number of free parameters is utilized and a sufficient number of partial differential moment equations is used to determine these parameters. Maxwell's inverse fifth-power force law is used to simplify the evaluation of the collision integrals appearing in the moment equations. All flow quantities are then determined by integration of the weighting function which results from the solution of the differential moment system. Three problems are treated: the heat-flux from a slightly heated sphere at rest in an infinite gas; the velocity field and drag of a slowly moving sphere in an unbounded space; the velocity field and drag torque on a slowly rotating sphere. Solutions to the third problem are found to both first and second-order in surface Mach number with the secondary centrifugal fan motion being of particular interest. Singular aspects of the moment method are encountered in the last two problems and an asymptotic study of these difficulties leads to a formal criterion for a "well posed" moment system. The previously unanswered question of just how many moments must be used in a specific problem is now clarified to a great extent.

A multiple scattering problem

The present work deals with the problem of the interaction of the electromagnetic radiation with a statistical distribution of nonmagnetic dielectric particles immersed in an infinite homogeneous isotropic, non-magnetic medium. The wavelength of the incident radiation can be less, equal or greater than the linear dimension of a particle. The distance between any two particles is several wavelengths. A single particle in the absence of the others is assumed to scatter like a Rayleigh-Gans particle, i.e. interaction between the volume elements (self-interaction) is neglected. The interaction of the particles is taken into account (multiple scattering) and conditions are set up for the case of a lossless medium which guarantee that the multiple scattering contribution is more important than the self-interaction one. These conditions relate the wavelength λ and the linear dimensions of a particle a and of the region occupied by the particles D. It is found that for constant λ/a, D is proportional to λ and that |Δχ|, where Δχ is the difference in the dielectric susceptibilities between particle and medium, has to lie within a certain range.

The total scattering field is obtained as a series the several terms of which represent the corresponding multiple scattering orders. The first term is a single scattering term. The ensemble average of the total scattering intensity is then obtained as a series which does not involve terms due to products between terms of different orders. Thus the waves corresponding to different orders are independent and their Stokes parameters add.

The second and third order intensity terms are explicitly computed. The method used suggests a general approach for computing any order. It is found that in general the first order scattering intensity pattern (or phase function) peaks in the forward direction Θ = 0. The second order tends to smooth out the pattern giving a maximum in the Θ = π/2 direction and minima in the Θ = 0 , Θ = π directions. This ceases to be true if ka (where k = 2π/λ) becomes large (> 20). For large ka the forward direction is further enhanced. Similar features are expected from the higher orders even though the critical value of ka may increase with the order.

The first order polarization of the scattered wave is determined. The ensemble average of the Stokes parameters of the scattered wave is explicitly computed for the second order. A similar method can be applied for any order. It is found that the polarization of the scattered wave depends on the polarization of the incident wave. If the latter is elliptically polarized then the first order scattered wave is elliptically polarized, but in the Θ = π/2 direction is linearly polarized. If the incident wave is circularly polarized the first order scattered wave is elliptically polarized except for the directions Θ = π/2 (linearly polarized) and Θ = 0, π (circularly polarized). The handedness of the Θ = 0 wave is the same as that of the incident whereas the handedness of the Θ = π wave is opposite. If the incident wave is linearly polarized the first order scattered wave is also linearly polarized. The second order makes the total scattered wave to be elliptically polarized for any Θ no matter what the incident wave is. However, the handedness of the total scattered wave is not altered by the second order. Higher orders have similar effects as the second order.

If the medium is lossy the general approach employed for the lossless case is still valid. Only the algebra increases in complexity. It is found that the results of the lossless case are insensitive in the first order of k_imD where k_im = imaginary part of the wave vector k and D a linear characteristic dimension of the region occupied by the particles. Thus moderately extended regions and small losses make (k_imD)² ≪ 1 and the lossy character of the medium does not alter the results of the lossless case. In general the presence of the losses tends to reduce the forward scattering.

Distributed Control and Optimization for Communication and Power Systems

We are at the cusp of a historic transformation of both communication system and electricity system. This creates challenges as well as opportunities for the study of networked systems. Problems of these systems typically involve a huge number of end points that require intelligent coordination in a distributed manner. In this thesis, we develop models, theories, and scalable distributed optimization and control algorithms to overcome these challenges.

This thesis focuses on two specific areas: multi-path TCP (Transmission Control Protocol) and electricity distribution system operation and control. Multi-path TCP (MP-TCP) is a TCP extension that allows a single data stream to be split across multiple paths. MP-TCP has the potential to greatly improve reliability as well as efficiency of communication devices. We propose a fluid model for a large class of MP-TCP algorithms and identify design criteria that guarantee the existence, uniqueness, and stability of system equilibrium. We clarify how algorithm parameters impact TCP-friendliness, responsiveness, and window oscillation and demonstrate an inevitable tradeoff among these properties. We discuss the implications of these properties on the behavior of existing algorithms and motivate a new algorithm Balia (balanced linked adaptation) which generalizes existing algorithms and strikes a good balance among TCP-friendliness, responsiveness, and window oscillation. We have implemented Balia in the Linux kernel. We use our prototype to compare the new proposed algorithm Balia with existing MP-TCP algorithms.

Our second focus is on designing computationally efficient algorithms for electricity distribution system operation and control. First, we develop efficient algorithms for feeder reconfiguration in distribution networks. The feeder reconfiguration problem chooses the on/off status of the switches in a distribution network in order to minimize a certain cost such as power loss. It is a mixed integer nonlinear program and hence hard to solve. We propose a heuristic algorithm that is based on the recently developed convex relaxation of the optimal power flow problem. The algorithm is efficient and can successfully computes an optimal configuration on all networks that we have tested. Moreover we prove that the algorithm solves the feeder reconfiguration problem optimally under certain conditions. We also propose a more efficient algorithm and it incurs a loss in optimality of less than 3% on the test networks.

Second, we develop efficient distributed algorithms that solve the optimal power flow (OPF) problem on distribution networks. The OPF problem determines a network operating point that minimizes a certain objective such as generation cost or power loss. Traditionally OPF is solved in a centralized manner. With increasing penetration of volatile renewable energy resources in distribution systems, we need faster and distributed solutions for real-time feedback control. This is difficult because power flow equations are nonlinear and kirchhoff's law is global. We propose solutions for both balanced and unbalanced radial distribution networks. They exploit recent results that suggest solving for a globally optimal solution of OPF over a radial network through a second-order cone program (SOCP) or semi-definite program (SDP) relaxation. Our distributed algorithms are based on the alternating direction method of multiplier (ADMM), but unlike standard ADMM-based distributed OPF algorithms that require solving optimization subproblems using iterative methods, the proposed solutions exploit the problem structure that greatly reduce the computation time. Specifically, for balanced networks, our decomposition allows us to derive closed form solutions for these subproblems and it speeds up the convergence by 1000x times in simulations. For unbalanced networks, the subproblems reduce to either closed form solutions or eigenvalue problems whose size remains constant as the network scales up and computation time is reduced by 100x compared with iterative methods.

The rate and mechanism of the partial oxidation of acetaldehyde by PPM concentrations of nitrogen dioxide

The thermal reaction between nitrogen dioxide and acetaldehyde in the gas phase was investigated at room temperature and atmospheric pressure. The initial rate of disappearance of nitrogen dioxide was 1.00 ± 0.03 order with respect to nitrogen dioxide and 1.00 ± 0.07 order with respect to acetaldehyde. An initial second order rate constant of (8.596 ± 0.189) x 10^-3 1.mole^-1 sec^-1 was obtained at 22.0 ± 0.1 °C and a total pressure of one atmosphere. The activation energy of the reaction was 12,900 cal/mole in the temperature range between 22°C and 122°C.

The products of the reaction were nitric oxide, carbon dioxide, methyl nitrite, nitromethane and a trace amount of trans-dimeric nitrosomethane. The addition of nitric oxide increased the rate of formation of nitromethane and decreased the rate of formation of methyl nitrite. There were no measurable surface effects due to the addition of glass wool or glass beads to the reactor.

Reactants and products were analyzed by gas chromatography. A mechanism was proposed incorporating the principal features of the reaction.

Photorefractive performance of a novel multifunctional inorganic-organic hybridized nanocomposite sensitized by CdS nanoparticles

A novel multifunctional inorganic-organic photorefractive (PR) poly(N-vinyl)-3-[p-nitrophenylazolcarbazolyl-CdS nanocomposites with different molar ratios of US to poly(N-vinyl)-3-[p-nitrophenylazo]carbazolyl (PVNPAK) were synthesized via a postazo-coupling reaction and chemically hybridized approach, respectively. The nanocomposites are highly soluble and could be obtained as film-forming materials with appreciably high molecular weights and low glass transition temperature (T,) due to the flexible spacers. The PVNPAK matrix possesses a highest-occupied molecular orbital value of about -5.36 eV determined from cyclic voltammetry. Second harmonic generation (SHG) could be observed in PVNPAK film without any poling procedure and 4.7 pm/V of effective second-order nonlinear optical susceptibility is obtained. The US particles as photosensitizers had a nanoscale size in PVNPAK adopting transmission electron microscopy. The improvement of interface quality between US and polymer matrix is responsible for efficient photoinduced charge generation efficiency in the nanocomposites. An asymmetric optical energy exchange between two beams on the polymer composites PVNPAK-CdS/ECZ has been found even without an external field in two-beam coupling (TBC) experiment, and the TBC gain and diffraction efficiency of 14.26 cm(-1) and 3.4% for PVNPAK-5-CdS/ECZ, 16.43 cm(-1) and 4.4% for PVNPAK-15-CdS/ECZ were measured at a 647.1 nm wavelength, respectively.

Measurements of mantle velocities of P waves with a large array

A large array has been used to investigate the P-wave velocity structure of the lower mantle. Linear array processing methods are reviewed and a method of nonlinear processing is presented. Phase velocities, travel times, and relative amplitudes of P waves have been measured with the large array at the Tonto Forest Seismological Observatory in Arizona for 125 earthquakes in the distance range of 30 to 100 degrees. Various models are assumed for the upper 771 km of the mantle and the Wiechert-Herglotz method applied to the phase velocity data to obtain a velocity depth structure for the lower mantle. The phase velocity data indicates the presence of a second-order discontinuity at a depth of 840 km, another at 1150 km, and less pronounced discontinuities at 1320, 1700 and 1950 km. Phase velocities beyond 85 degrees are interpreted in terms of a triplication of the phase velocity curve, and this results in a zone of almost constant velocity between depths of 2670 and 2800 km. Because of the uncertainty in the upper mantle assumptions, a final model cannot be proposed, but it appears that the lower mantle is more complicated than the standard models and there is good evidence for second-order discontinuities below a depth of 1000 km. A tentative lower bound of 2881 km can be placed on the depth to the core. The importance of checking the calculated velocity structure against independently measured travel times is pointed out. Comparisons are also made with observed PcP times and the agreement is good. The method of using measured values of the rate of change of amplitude with distances shows promising results.

On the uniqueness of singular solutions to boundary-initial value problems in linear elastodynamics

This investigation deals with certain generalizations of the classical uniqueness theorem for the second boundary-initial value problem in the linearized dynamical theory of not necessarily homogeneous nor isotropic elastic solids. First, the regularity assumptions underlying the foregoing theorem are relaxed by admitting stress fields with suitably restricted finite jump discontinuities. Such singularities are familiar from known solutions to dynamical elasticity problems involving discontinuous surface tractions or non-matching boundary and initial conditions. The proof of the appropriate uniqueness theorem given here rests on a generalization of the usual energy identity to the class of singular elastodynamic fields under consideration.

Following this extension of the conventional uniqueness theorem, we turn to a further relaxation of the customary smoothness hypotheses and allow the displacement field to be differentiable merely in a generalized sense, thereby admitting stress fields with square-integrable unbounded local singularities, such as those encountered in the presence of focusing of elastic waves. A statement of the traction problem applicable in these pathological circumstances necessitates the introduction of "weak solutions'' to the field equations that are accompanied by correspondingly weakened boundary and initial conditions. A uniqueness theorem pertaining to this weak formulation is then proved through an adaptation of an argument used by O. Ladyzhenskaya in connection with the first boundary-initial value problem for a second-order hyperbolic equation in a single dependent variable. Moreover, the second uniqueness theorem thus obtained contains, as a special case, a slight modification of the previously established uniqueness theorem covering solutions that exhibit only finite stress-discontinuities.

Corrections to the paraxial approximation solutions in transversely nonuniform refractive-index media

The coupled differential recurrence equations for the corrections to the paraxial approximation solutions in transversely nonuniform refractive-index media are established in terms of the perturbation method. All the corrections (including the longitudinal field corrections) to the paraxial approximation solutions are presented in the weak-guidance approximation. As a concrete application, the first-order longitudinal field correction and the second-order transverse field correction to the paraxial approximation of a Gaussian beam propagating in a transversely quadratic refractive index medium are analytically investigated. (C) 1999 Optical Society of America [S0740-3232(99)00310-5].

Efeito da temperatura de incubação na cinética e no equilíbrio da biossorção do cádmio por Sargassum filipendula

O uso de biomassas para biossorção de metais pesados é bem documentado na literatura e vários tipos de espécies de microrganismos e algas já foram testados. A maior parte destes trabalhos foi realizada com biomassa seca para prevenir qualquer resposta metabólica indesejável. Vários estudos na literatura sugerem o uso de biomassa seca sobre condições moderadas, tais como secagem ao sol; por outro lado, vários trabalhos recomendam a faixa de 313K a 353K para garantir completa inativação da biomassa. O efeito da biomassa seca ao sol a 303K e seca a 333K em estufa na remoção de Cd2+ é aqui reportado. A avaliação dos resultados foi baseada na cinética e capacidade de remoção do metal pela alga Sargassum filipendula. Os resultados indicam que a adsorção máxima de metal não foi notadamente reduzida quando a biomassa seca em estufa foi usada, para concentrações de cádmio na faixa de 10,0 a 500,0 mg L-1. O estudo cinético realizado indicou que o modelo de pseudo segunda ordem ajustou melhor os dados experimentais, tanto para uma solução diluída (10 mg L-1) quanto para a concentrada (100 mg L-1). Em ambos os casos, os efeitos da secagem em estufa, a 60C refletiu-se suavemente na remoção do metal. Os dados experimentais foram melhor ajustados pelo modelo de Langmuir em comparação com o modelo de Freundlich. Análises termogravimétricas mostraram que não havia dano estrutural no biossorvente devido à secagem em estufa. O espectro de infravermelho não indicou diferença entre a biomassa in natura e seca. O efeito da temperatura na biossorção do metal significativo na faixa de 303K a 328K, refletindo-se na capacidade de remoção do cádmio

Aplicações de métodos de classificação e calibração multivariada acoplados com técnicas espectroscópicas em amostras ambientais e alimentos

Este trabalho de pesquisa descreve dois estudos de caso de métodos quimiométricos empregados para a quantificação de hidrocarbonetos policíclicos aromáticos HPAs (naftaleno, fluoreno, fenantreno e fluoranteno) em água potável usando espectroscopia de fluorescência molecular e a classificação e caracterização de sucos de uva e seus parâmetros de qualidade através de espectroscopia de infravermelho próximo. O objetivo do primeiro estudo é a aplicação combinada de métodos quimiométricos de segunda ordem (N-PLS, U-PLS, U-PLS/RBL e PARAFAC) e espectrofluorimetria para determinação direta de HPAs em água potável, visando contribuir para o conhecimento do potencial destas metodologias como alternativa viável para a determinação tradicional por cromatografia univariada. O segundo estudo de caso destinado à classificação e determinação de parâmetros de qualidade de sucos de uva, densidade relativa e teor de sólidos solúveis totais, foi medida por espectroscopia de infravermelho próximo e métodos quimiométricos. Diversos métodos quimiométricos, tais como HCA, PLS-DA, SVM-DA e SIMCA foram investigados para a classificação amostras de sucos de uva ao mesmo tempo que métodos de calibração multivariada de primeira ordem, tais como PLS, iPLS e SVM-LS foram usadas para a predição dos parâmetros de qualidade. O princípio orientador para o desenvolvimento dos estudos aqui descritos foi a necessidade de metodologias analíticas com custo, tempo de execução e facilidade de operação melhores e menor produção de resíduos do que os métodos atualmente utilizados para a quantificação de HPAs, em água de torneira, e classificação e caracterização das amostras de suco de uva e seus parâmetros de qualidade