On the Dirac equation with PT-symmetric potentials in the presence of position-dependent mass

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

The viscosity parameter's dependence on the profile of the disk scale height

It is shown that, for accretion disks, the height scale is a constant whenever hydrostatic equilibrium and the subsonic turbulence regime hold in the disk. In order to have a variable height scale, processes are needed that contribute an extra term to the continuity equation. This contribution makes the viscosity parameter much greater in the outer region and much smaller in the inner region. Under these circumstances, turbulence is the presumable source of viscosity in the disk.

The measurement of ice velocity, mass balance and thinning-rate on Johnsons Glacier, Livingston Island, South Shetland Islands, Antarctica

A network of twenty stakes was set up on Johnsons Glacier in order to determine its dynamics. During the austral summers from 1994-95 to 1997-98, we estimated surface velocities, mass balances and ice thickness variations. Horizontal velocity increased dow nstream from 1 m a- 1 near the ice divides to 40 m a- 1 near the ice terminus. The accumulation zone showed low accumulation rates (maximum of 0,6 m a- 1 (ice)), whereas in the lower part of the glacier, ablation rates were 4,3 m a- 1 (ice). Over the 3-year study period, both in the accumulation and ablation zones, we detected a reduction in the ice surface level ranging from 2 to 10 m from the annual ve rt ical velocities and ice-thinning data, the mass balance was obtained and compared with the mass balance field values, resulting in similar estimates. Flux values were calculated using cross-section data and horizontal velocities, and compared with the results obtained by means of mass balance and ice thinning data using the continuity equation. The two methods gave similar results.

Open solar flux estimates from near-Earth measurements of the interplanetary magnetic field: comparison of the first two perihelion passes of the Ulysses spacecraft

Results from all phases of the orbits of the Ulysses spacecraft have shown that the magnitude of the radial component of the heliospheric field is approximately independent of heliographic latitude. This result allows the use of near- Earth observations to compute the total open flux of the Sun. For example, using satellite observations of the interplanetary magnetic field, the average open solar flux was shown to have risen by 29% between 1963 and 1987 and using the aa geomagnetic index it was found to have doubled during the 20th century. It is therefore important to assess fully the accuracy of the result and to check that it applies to all phases of the solar cycle. The first perihelion pass of the Ulysses spacecraft was close to sunspot minimum, and recent data from the second perihelion pass show that the result also holds at solar maximum. The high level of correlation between the open flux derived from the various methods strongly supports the Ulysses discovery that the radial field component is independent of latitude. We show here that the errors introduced into open solar flux estimates by assuming that the heliospheric field’s radial component is independent of latitude are similar for the two passes and are of order 25% for daily values, falling to 5% for averaging timescales of 27 days or greater. We compare here the results of four methods for estimating the open solar flux with results from the first and second perehelion passes by Ulysses. We find that the errors are lowest (1–5% for averages over the entire perehelion passes lasting near 320 days), for near-Earth methods, based on either interplanetary magnetic field observations or the aa geomagnetic activity index. The corresponding errors for the Solanki et al. (2000) model are of the order of 9–15% and for the PFSS method, based on solar magnetograms, are of the order of 13–47%. The model of Solanki et al. is based on the continuity equation of open flux, and uses the sunspot number to quantify the rate of open flux emergence. It predicts that the average open solar flux has been decreasing since 1987, as Correspondence to: M. Lockwood (m.lockwood@rl.ac.uk) is observed in the variation of all the estimates of the open flux. This decline combines with the solar cycle variation to produce an open flux during the second (sunspot maximum) perihelion pass of Ulysses which is only slightly larger than that during the first (sunspot minimum) perihelion pass.

Genesis of quantum nonlocality

We revisit the problem of an otherwise classical particle immersed in the zero-point radiation field, with the purpose of tracing the origin of the nonlocality characteristic of Schrodinger`s equation. The Fokker-Planck-type equation in the particles phase-space leads to an infinite hierarchy of equations in configuration space. In the radiationless limit the first two equations decouple from the rest. The first is the continuity equation: the second one, for the particle flux, contains a nonlocal term due to the momentum fluctuations impressed by the field. These equations are shown to lead to Schrodinger`s equation. Nonlocality (obtained here for the one-particle system) appears thus as a property of the description, not of Nature. (C) 2011 Elsevier B.V. All rights reserved.

Understanding the tsunami with a simple model

In this paper, we use the approximation of shallow water waves (Margaritondo G 2005 Eur. J. Phys. 26 401) to understand the behaviour of a tsunami in a variable depth. We deduce the shallow water wave equation and the continuity equation that must be satisfied when a wave encounters a discontinuity in the sea depth. A short explanation about how the tsunami hit the west coast of India is given based on the refraction phenomenon. Our procedure also includes a simple numerical calculation suitable for undergraduate students in physics and engineering.

A marker-and-cell approach to free surface 2-D multiphase flows

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Electrofluid dynamic (EFD) energy conversion in gas flows with suspended nanoparticles

This work presents simulations of the Electrofluid Dynamic energy conversion process in slender channel devices having very small particles (in both micro and nano scales) as charge carriers. Solutions are discussed for a system composed by coupled differential equations, which includes the equation for the total current along the channel, the equations for total energy and momentum of the mixture (gas and solid particles), the continuity equation and the equations for energy and momentum of a single particle. Results for suspended particles of higher diameters have been previously published in the Literature, but the simulations here presented exhibit an appreciable increase in the values for output currents.

Electrofluid dynamic (EFD) energy conversion in gas flows with suspended nanoparticles

This work presents simulations of the Electrofluid Dynamic energy conversion process in slender channel devices having very small particles (in both micro and nano scales) as charge carriers. Solutions are discussed for a system composed by coupled differential equations, which includes the equation for the total current along the channel, the equations for total energy and momentum of the mixture (gas and solid particles), the continuity equation and the equations for energy and momentum of a single particle. Results for suspended particles of higher diameters have been previously published in the Literature, but the simulations here presented exhibit an appreciable increase in the values for output currents.

A marker-and-cell approach to free surface 2-D multiphase flows

This work describes a methodology to simulate free surface incompressible multiphase flows. This novel methodology allows the simulation of multiphase flows with an arbitrary number of phases, each of them having different densities and viscosities. Surface and interfacial tension effects are also included. The numerical technique is based on the GENSMAC front-tracking method. The velocity field is computed using a finite-difference discretization of a modification of the NavierStokes equations. These equations together with the continuity equation are solved for the two-dimensional multiphase flows, with different densities and viscosities in the different phases. The governing equations are solved on a regular Eulerian grid, and a Lagrangian mesh is employed to track free surfaces and interfaces. The method is validated by comparing numerical with analytic results for a number of simple problems; it was also employed to simulate complex problems for which no analytic solutions are available. The method presented in this paper has been shown to be robust and computationally efficient. Copyright (c) 2012 John Wiley & Sons, Ltd.

Studio Sperimentale e Modellazione della Separazione di Proteine con Membrane di Affinità

Chromatography represents one of the most important and widely used unit operation in the biotechnology industry. However this technique suffers from several limitations such as high pressure drop, slow mass transfer through the diffusive pores and strong dependence of the binding capacity on flow rate. In this work, affinity membranes with improved capacity have been considered as an alternative technology for the capturing step in antibody manufacturing. Several affinity membranes have been prepared starting from various membrane supports. Different affinity ligands have been utilized like Protein A, the natural ligand of choice for antibodies, as well as synthetic ligands that exhibit affinity for the Fc portion of antibodies. The membranes have been characterized in detail: binding and elution performance was evaluated in adsorption experiments using pure IgG solutions, while membrane selectivity was evaluated using complex solutions like a cell culture supernatant. The most promising affinity membranes were extensively tested in dynamic experiments. The effects of operating parameters like feed concentration and flow rate on separation performances like binding capacity, selectivity and process yield have been studied in detail in order to find the optimal conditions for binding and elution steps. The membranes have been used over several complete chromatographic cycles to evaluate the effects of ageing and of membrane regeneration on dynamic binding capacity. A novel mathematical model is proposed that can describe all the chromatographic steps involved in the membrane affinity chromatography process for protein purification. The mathematical description is based on the species continuity equation coupled with a proper binding kinetic equation, and suitable to describe adequately the dispersion phenomena occurring both in the micro-porous membranes as well as in the extra-column devices used in the system. The model considers specifically all the different chromatographic steps, namely adsorption, washing and elution. The few relevant fitting parameters of the model were derived from a calibration with the experimental affinity cycles performed with pure IgG solutions, then the model is used to describe experimental data obtained in chromatographic cycles carried out with complex feeds as the cell culture supernatant. Simulations reveal a good agreement with experimental data in all the chromatography steps, both in the case of pure IgG solutions and for the cell culture supernatant considered.

Numerical studies of colloidal systems

The interplay of hydrodynamic and electrostatic forces is of great importance for the understanding of colloidal dispersions. Theoretical descriptions are often based on the so called standard electrokinetic model. This Mean Field approach combines the Stokes equation for the hydrodynamic flow field, the Poisson equation for electrostatics and a continuity equation describing the evolution of the ion concentration fields. In the first part of this thesis a new lattice method is presented in order to efficiently solve the set of non-linear equations for a charge-stabilized colloidal dispersion in the presence of an external electric field. Within this framework, the research is mainly focused on the calculation of the electrophoretic mobility. Since this transport coefficient is independent of the electric field only for small driving, the algorithm is based upon a linearization of the governing equations. The zeroth order is the well known Poisson-Boltzmann theory and the first order is a coupled set of linear equations. Furthermore, this set of equations is divided into several subproblems. A specialized solver for each subproblem is developed, and various tests and applications are discussed for every particular method. Finally, all solvers are combined in an iterative procedure and applied to several interesting questions, for example, the effect of the screening mechanism on the electrophoretic mobility or the charge dependence of the field-induced dipole moment and ion clouds surrounding a weakly charged sphere. In the second part a quantitative data analysis method is developed for a new experimental approach, known as "Total Internal Reflection Fluorescence Cross-Correlation Spectroscopy" (TIR-FCCS). The TIR-FCCS setup is an optical method using fluorescent colloidal particles to analyze the flow field close to a solid-fluid interface. The interpretation of the experimental results requires a theoretical model, which is usually the solution of a convection-diffusion equation. Since an analytic solution is not available due to the form of the flow field and the boundary conditions, an alternative numerical approach is presented. It is based on stochastic methods, i. e. a combination of a Brownian Dynamics algorithm and Monte Carlo techniques. Finally, experimental measurements for a hydrophilic surface are analyzed using this new numerical approach.

Three dimensional circulation of water with variable wind in the Caspian Sea

This research is based on a numerical model for forecasting the three-dimensional behavior of (sea) water motion due to the effect of a variable wind velocity. The results obtained are then analyzed and compared with observation. This model is based on the equations that overcome the current and distribution of temperature by applying the method of finite difference with assuming Δx, Δy as constant and Δz, variable. The model is based on the momentum equation, continuity equation and thermodynamic energy equation and tension at the surface and middle layers and bottom stress. The horizontal and vertical eddy viscosity and thermal diffusivity coefficients we used in accordance with that of the Bennet on Outario Lake (1977). Considering the Caspian Sea dimension in numerical model the Coriolis parameter used with β effects and the approximation Boussines have been used. For the program controlling some simple experiment with boundary condition similar to that of the Caspian Sea have been done. For modeling the Caspian Sea the grid of the field was done as follows: At horizontal surface grid size is 10×10km extension and at vertical in 10 layers with varying thickness from surface to bed respectively as: 5, 10, 20, 3, 50, 100, 150, 200, 25, 500 and higher. The data of wind as velocity، direction and temperature of water related to 15th September 1995 at 6،12 and 18 o’clock were obtained from synoptic station at the Caspian Sea shore and the research marine of Haji Alief. The information concerning shore wind was measured and by the method of SPM (shore protection manual) was transferred to far shore winds through interpolation and by use of inverse square distance of position distribution of the wind velocity at the Caspian surface field was obtained. The model has been evaluated according to the reports and observations. Through studying the position of the current in different layers، the velocity in the cross section in the northern، southern and the middle layers، will be discussed. The results reveal the presence of the circulation cells in the three above mentioned areas. The circulation with depth is reduced too. The results obtained through the numerical solution of the temperature equation have been compared with the observation. The temperature change in different layers in cross section illustrates the relative accordance of the model mentioned.

Understanding the tsunami with a simple model

In this paper, we use the approximation of shallow water waves (Margaritondo G 2005 Eur. J. Phys. 26 401) to understand the behaviour of a tsunami in a variable depth. We deduce the shallow water wave equation and the continuity equation that must be satisfied when a wave encounters a discontinuity in the sea depth. A short explanation about how the tsunami hit the west coast of India is given based on the refraction phenomenon. Our procedure also includes a simple numerical calculation suitable for undergraduate students in physics and engineering.

Détection de mouvements dans des séquences d’images basée sur la dynamique de supraconductivité

Le but de ce travail est d’étudier la faisabilité de la détection de mouvements dans des séquences d’images en utilisant l’équation de continuité et la dynamique de supraconductivité. Notre approche peut être motivée par le fait que l’équation de continuité apparait dans plusieurs techniques qui estiment le flot optique. Un grand nombre de techniques qui utilisent les flots optiques utilisent une contrainte appelée contrainte de l’invariance lumineuse. La dynamique de supraconductivité nous permet de nous affranchir de la contrainte de l’invariance lumineuse. Les expériences se feront avec la base de données de séquences d’images CDNET 2014. Pour obtenir les résultats numériques en terme de score F1, une combinaison sera faite par la suite entre la dynamique de supraconductivité et un méchanisme d’attention qui est un résumé des vérites de terrain.