Elerctron Transport Coefficients in N2 in Crossed Fields at Low E/p

The numerical solutions of Boltzmann transpott equation for the energy distribution of electrons moving in crossed fields in nitrogen have been obtained for 100 ÿ E/p ÿ 1000 V M-1 Torr-1 and for 0ÿ B/p ÿ 0.02 Tesla Torr-1 using the concept of energy dependent effective field intensity. From the derived distribution functions the electron mean energy, the tranaverse and perpendicular drift velocities and the averaged effective field intensity (Eavef) which signifies the average field intensity experienced by electron swarms in E àB field have been derived. The maximum difference between the electron mean energy for a given E ÃÂB field and that corresponding to Eavef/p (p is the gas pressure) is found to be within ñ3.5%.

Study of the vertical transport in p-doped superlattices based on group III-V semiconductors

The electrical conductivity σ has been calculated for p-doped GaAs/Al0.3Ga0.7As and cubic GaN/Al0.3Ga0.7N thin superlattices (SLs). The calculations are done within a self-consistent approach to the k → ⋅ p → theory by means of a full six-band Luttinger-Kohn Hamiltonian, together with the Poisson equation in a plane wave representation, including exchange correlation effects within the local density approximation. It was also assumed that transport in the SL occurs through extended minibands states for each carrier, and the conductivity is calculated at zero temperature and in low-field ohmic limits by the quasi-chemical Boltzmann kinetic equation. It was shown that the particular minibands structure of the p-doped SLs leads to a plateau-like behavior in the conductivity as a function of the donor concentration and/or the Fermi level energy. In addition, it is shown that the Coulomb and exchange-correlation effects play an important role in these systems, since they determine the bending potential.

Role of different scattering mechanisms on the temperature dependence of transport in graphene

Detailed experimental and theoretical studies of the temperature dependence of the effect of different scattering mechanisms on electrical transport properties of graphene devices are presented. We find that for high mobility devices the transport properties are mainly governed by completely screened short range impurity scattering. On the other hand, for the low mobility devices transport properties are determined by both types of scattering potentials - long range due to ionized impurities and short range due to completely screened charged impurities. The results could be explained in the framework of Boltzmann transport equations involving the two independent scattering mechanisms.

Um método SN híbrido direto para cálculos de sistemas combustível-moderador em geometria unidimensional

Desenvolvemos nesta dissertação um método híbrido direto para o cálculo do fator de desvantagem e descrição da distribuição do fluxo de nêutrons em sistemas combustível-moderador. Na modelagem matemática, utilizamos a equação de transporte de Boltzmann independente do tempo, considerando espalhamento linearmente anisotrópico no modelo monoenergético e espalhamento isotrópico no modelo multigrupo, na formulação de ordenadas discretas (SN), em geometria unidimensional. Desenvolvemos nesta dissertação um método híbrido direto para o cálculo do fator de desvantagem e descrição da distribuição do fluxo de nêutrons em sistemas combustível-moderador. Na modelagem matemática, utilizamos a equação de transporte de Boltzmann independente do tempo, considerando espalhamento linearmente anisotrópico no modelo monoenergético e espalhamento isotrópico no modelo multigrupo, na formulação de ordenadas discretas (SN), em geometria unidimensional. Descrevemos uma análise espectral das equações de ordenadas discretas (SN)a um grupo e a dois grupos de energia, onde seguimos uma analogia com o método de Case. Utilizamos, neste método, quadraturas angulares diferentes no combustível (NC) e no moderador (NM), onde em geral assumimos que NC > NM . Condições de continuidade especiais que acoplam os fluxos angulares que emergem do combustível (moderador) e incidem no moderador (combustível), foram utilizadas com base na equivalência entre as equações SN e PN-1, o que caracteriza a propriedade híbrida do modelo proposto. Sendo um método híbrido direto, utilizamos as NC + NM equações lineares e algébricas constituídas pelas (NC + NM)/2 condições de contorno reflexivas e (NC + NM)/2 condições de continuidade para determinarmos as NC + NM constantes. Com essas constantes podemos calcular os valores dos fluxos angulares e dos fluxos escalares em qualquer ponto do domínio. Apresentamos resultados numéricos para ilustrar a eficiência e a precisão do método proposto.

Solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico e aproximação sintética de difusão

Nesta dissertação, são apresentados os seguintes modelos matemáticos de transporte de nêutrons: a equação linearizada de Boltzmann e a equação da difusão de nêutrons monoenergéticos em meios não-multiplicativos. Com o objetivo de determinar o período fluxo escalar de nêutrons, é descrito um método espectronodal que gera soluções numéricas para o problema de difusão em geometria planar de fonte fixa, que são livres de erros de truncamento espacial, e que conjugado com uma técnica de reconstrução espacial intranodal gera o perfil detalhado da solução. A fim de obter o valor aproximado do fluxo angular de nêutrons em um determinado ponto do domínio e em uma determinada direção de migração, descreve-se também um método de reconstrução angular baseado na solução analítica da equação unidimensional de transporte de nêutrons monoenergéticos com espalhamento linearmente anisotrópico com aproximação sintética de difusão nos termos de fonte por espalhamento. O código computacional desenvolvido nesta dissertação foi implementado na plataforma livre Scilab, e para ilustrar a eficiência do código criado,resultados numéricos obtidos para três problemas-modelos são apresentados

Intermittency transport modeling of separated flow transition

An intermittency transport model is proposed for modeling separated-flow transition. The model is based on earlier work on prediction of attached flow bypass transition and is applied for the first time to model transition in a separation bubble at various degrees of free-stream turbulence. The model has been developed so that it takes into account the entrainment of the surrounding fluid. Experimental investigations suggest that it is this phenomena which ultimately determines the extent of the separation bubble. Transition onset is determined via a boundary layer correlation based on momentum thickness at the point of separation. The intermittent flow characteristic of the transition process is modeled via an intermittency transport equation. This accounts for both normal and streamwise variation of intermittency and hence models the entrainment of surrounding flow in a more accurate manner than alternative prescribed intermittency models. The model has been validated against the well established T3L semicircular leading edge flat plate test case for three different degrees of free-stream turbulence characteristic of turbomachinery blade applications.

Statistics and modelling of turbulent kinetic energy transport in different regimes of premixed combustion

The statistical behaviour of turbulent kinetic energy transport in turbulent premixed flames is analysed using data from three-dimensional Direct Numerical Simulation (DNS) of freely propagating turbulent premixed flames under decaying turbulence. For flames within the corrugated flamelets regime, it is observed that turbulent kinetic energy is generated within the flame brush. By contrast, for flames within the thin reaction zones regime it has been found that the turbulent kinetic energy decays monotonically through the flame brush. Similar trends are observed also for the dissipation rate of turbulent kinetic energy. Within the corrugated flamelets regime, it is demonstrated that the effects of the mean pressure gradient and pressure dilatation within the flame are sufficient to overcome the effects of viscous dissipation and are responsible for the observed augmentation of turbulent kinetic energy in the flame brush. In the thin reaction zones regime, the effects of the mean pressure gradient and pressure dilatation terms are relatively much weaker than those of viscous dissipation, resulting in a monotonic decay of turbulent kinetic energy across the flame brush. The modelling of the various unclosed terms of the turbulent kinetic energy transport equation has been analysed in detail. The predictions of existing models are compared with corresponding quantities extracted from DNS data. Based on this a-priori DNS assessment, either appropriate models are identified or new models are proposed where necessary. It is shown that the turbulent flux of turbulent kinetic energy exhibits counter-gradient (gradient) transport wherever the turbulent scalar flux is counter-gradient (gradient) in nature. A new model has been proposed for the turbulent flux of turbulent kinetic energy, and is found to capture the qualitative and quantitative behaviour obtained from DNS data for both the corrugated flamelets and thin reaction zones regimes without the need to adjust any of the model constants. © 2010 Springer Science+Business Media B.V.

Local strain rate and curvature dependences of scalar dissipation rate transport in turbulent premixed flames: A direct numerical simulation analysis

The statistical behaviours of the instantaneous scalar dissipation rate Nc of reaction progress variable c in turbulent premixed flames have been analysed based on three-dimensional direct numerical simulation data of freely propagating statistically planar flame and V-flame configurations with different turbulent Reynolds number Ret. The statistical behaviours of N c and different terms of its transport equation for planar and V-flames are found to be qualitatively similar. The mean contribution of the density-variation term T1 is positive, whereas the molecular dissipation term (-D2) acts as a leading order sink. The mean contribution of the strain rate term T2 is predominantly negative for the cases considered here. The mean reaction rate contribution T3 is positive (negative) towards the unburned (burned) gas side of the flame, whereas the mean contribution of the diffusivity gradient term (D) assumes negative (positive) values towards the unburned (burned) gas side. The local statistical behaviours of Nc, T1, T2, T 3, (-D2), and f(D) have been analysed in terms of their marginal probability density functions (pdfs) and their joint pdfs with local tangential strain rate aT and curvature km. Detailed physical explanations have been provided for the observed behaviour. © 2014 Y. Gao et al.

Nonlinear Dynamics for Local Fractional Burgers Equation Arising in Fractal Flow

The local fractional Burgers’ equation (LFBE) is investigated from the point of view of local fractional conservation laws envisaging a nonlinear local fractional transport equation with a linear non-differentiable diffusion term. The local fractional derivative transformations and the LFBE conversion to a linear local fractional diffusion equation are analyzed.

Teoria cinética não extensiva: efeitos físicos em gases e plasmas

The standard kinetic theory for a nonrelativistic diluted gas is generalized in the spirit of the nonextensive statistic distribution introduced by Tsallis. The new formalism depends on an arbitrary q parameter measuring the degree of nonextensivity. In the limit q = 1, the extensive Maxwell-Boltzmann theory is recovered. Starting from a purely kinetic deduction of the velocity q-distribution function, the Boltzmann H-teorem is generalized for including the possibility of nonextensive out of equilibrium effects. Based on this investigation, it is proved that Tsallis' distribution is the necessary and sufficient condition defining a thermodynamic equilibrium state in the nonextensive context. This result follows naturally from the generalized transport equation and also from the extended H-theorem. Two physical applications of the nonextensive effects have been considered. Closed analytic expressions were obtained for the Doppler broadening of spectral lines from an excited gas, as well as, for the dispersion relations describing the eletrostatic oscillations in a diluted electronic plasma. In the later case, a comparison with the experimental results strongly suggests a Tsallis distribution with the q parameter smaller than unity. A complementary study is related to the thermodynamic behavior of a relativistic imperfect simple fluid. Using nonequilibrium thermodynamics, we show how the basic primary variables, namely: the energy momentum tensor, the particle and entropy fluxes depend on the several dissipative processes present in the fluid. The temperature variation law for this moving imperfect fluid is also obtained, and the Eckart and Landau-Lifshitz formulations are recovered as particular cases

Dirac fermions in strong electric field and quantum transport in graphene

Our previous results on the nonperturbative calculations of the mean current and of the energy-momentum tensor in QED with the T-constant electric field are generalized to arbitrary dimensions. The renormalized mean values are found, and the vacuum polarization contributions and particle creation contributions to these mean values are isolated in the large T limit; we also relate the vacuum polarization contributions to the one-loop effective Euler-Heisenberg Lagrangian. Peculiarities in odd dimensions are considered in detail. We adapt general results obtained in 2 + 1 dimensions to the conditions which are realized in the Dirac model for graphene. We study the quantum electronic and energy transport in the graphene at low carrier density and low temperatures when quantum interference effects are important. Our description of the quantum transport in the graphene is based on the so-called generalized Furry picture in QED where the strong external field is taken into account nonperturbatively; this approach is not restricted to a semiclassical approximation for carriers and does not use any statistical assumptions inherent in the Boltzmann transport theory. In addition, we consider the evolution of the mean electromagnetic field in the graphene, taking into account the backreaction of the matter field to the applied external field. We find solutions of the corresponding Dirac-Maxwell set of equations and with their help we calculate the effective mean electromagnetic field and effective mean values of the current and the energy-momentum tensor. The nonlinear and linear I-V characteristics experimentally observed in both low-and high-mobility graphene samples are quite well explained in the framework of the proposed approach, their peculiarities being essentially due to the carrier creation from the vacuum by the applied electric field. DOI: 10.1103/PhysRevD.86.125022

Fast and accurate numerical solutions in some problems of particle and radiation transport: synthetic acceleration for the method of short characteristics, Doppler-broadened scattering kernel, remote sensing of the cryosphere

The aim of this work is to present various aspects of numerical simulation of particle and radiation transport for industrial and environmental protection applications, to enable the analysis of complex physical processes in a fast, reliable, and efficient way. In the first part we deal with speed-up of numerical simulation of neutron transport for nuclear reactor core analysis. The convergence properties of the source iteration scheme of the Method of Characteristics applied to be heterogeneous structured geometries has been enhanced by means of Boundary Projection Acceleration, enabling the study of 2D and 3D geometries with transport theory without spatial homogenization. The computational performances have been verified with the C5G7 2D and 3D benchmarks, showing a sensible reduction of iterations and CPU time. The second part is devoted to the study of temperature-dependent elastic scattering of neutrons for heavy isotopes near to the thermal zone. A numerical computation of the Doppler convolution of the elastic scattering kernel based on the gas model is presented, for a general energy dependent cross section and scattering law in the center of mass system. The range of integration has been optimized employing a numerical cutoff, allowing a faster numerical evaluation of the convolution integral. Legendre moments of the transfer kernel are subsequently obtained by direct quadrature and a numerical analysis of the convergence is presented. In the third part we focus our attention to remote sensing applications of radiative transfer employed to investigate the Earth's cryosphere. The photon transport equation is applied to simulate reflectivity of glaciers varying the age of the layer of snow or ice, its thickness, the presence or not other underlying layers, the degree of dust included in the snow, creating a framework able to decipher spectral signals collected by orbiting detectors.

Analysis of nonlinear gain-induced effects on short-pulse amplification in doped fibers by use of an extended power equation

Optical pulse amplification in doped fibers is studied using an extended power transport equation for the coupled pulse spectral components. This equation includes the effects of gain saturation, gain dispersion, fiber dispersion, fiber nonlinearity, and amplified spontaneous emission. The new model is employed to study nonlinear gain-induced effects on the spectrotemporal characteristics of amplified subpicosecond pulses, in both the anomalous and the normal dispersion regimes.

A numerical model for inert gas transport in the lung based on fractal airway morphology

Patients suffering from cystic fibrosis (CF) show thick secretions, mucus plugging and bronchiectasis in bronchial and alveolar ducts. This results in substantial structural changes of the airway morphology and heterogeneous ventilation. Disease progression and treatment effects are monitored by so-called gas washout tests, where the change in concentration of an inert gas is measured over a single or multiple breaths. The result of the tests based on the profile of the measured concentration is a marker for the severity of the ventilation inhomogeneity strongly affected by the airway morphology. However, it is hard to localize underlying obstructions to specific parts of the airways, especially if occurring in the lung periphery. In order to support the analysis of lung function tests (e.g. multi-breath washout), we developed a numerical model of the entire airway tree, coupling a lumped parameter model for the lung ventilation with a 4th-order accurate finite difference model of a 1D advection-diffusion equation for the transport of an inert gas. The boundary conditions for the flow problem comprise the pressure and flow profile at the mouth, which is typically known from clinical washout tests. The natural asymmetry of the lung morphology is approximated by a generic, fractal, asymmetric branching scheme which we applied for the conducting airways. A conducting airway ends when its dimension falls below a predefined limit. A model acinus is then connected to each terminal airway. The morphology of an acinus unit comprises a network of expandable cells. A regional, linear constitutive law describes the pressure-volume relation between the pleural gap and the acinus. The cyclic expansion (breathing) of each acinus unit depends on the resistance of the feeding airway and on the flow resistance and stiffness of the cells themselves. Special care was taken in the development of a conservative numerical scheme for the gas transport across bifurcations, handling spatially and temporally varying advective and diffusive fluxes over a wide range of scales. Implicit time integration was applied to account for the numerical stiffness resulting from the discretized transport equation. Local or regional modification of the airway dimension, resistance or tissue stiffness are introduced to mimic pathological airway restrictions typical for CF. This leads to a more heterogeneous ventilation of the model lung. As a result the concentration in some distal parts of the lung model remains increased for a longer duration. The inert gas concentration at the mouth towards the end of the expirations is composed of gas from regions with very different washout efficiency. This results in a steeper slope of the corresponding part of the washout profile.

Modelling of chloride transport in non-saturated concrete : from microscale to macroscale

Corrosion of steel bars embedded in concrete has a great influence on structural performance and durability of reinforced concrete. Chloride penetration is considered to be a primary cause of concrete deterioration in a vast majority of structures. Therefore, modelling of chloride penetration into concrete has become an area of great interest. The present work focuses on modelling of chloride transport in concrete. The differential macroscopic equations which govern the problem were derived from the equations at the microscopic scale by comparing the porous network with a single equivalent pore whose properties are the same as the average properties of the real porous network. The resulting transport model, which accounts for diffusion, migration, advection, chloride binding and chloride precipitation, consists of three coupled differential equations. The first equation models the transport of chloride ions, while the other two model the flow of the pore water and the heat transfer. In order to calibrate the model, the material parameters to determine experimentally were identified. The differential equations were solved by means of the finite element method. The classical Galerkin method was employed for the pore solution flow and the heat transfer equations, while the streamline upwind Petrov Galerkin method was adopted for the transport equation in order to avoid spatial instabilities for advection dominated problems. The finite element codes are implemented in Matlab® . To retrieve a good understanding of the influence of each variable and parameter, a detailed sensitivity analysis of the model was carried out. In order to determine the diffusive and hygroscopic properties of the studied concretes, as well as their chloride binding capacity, an experimental analysis was performed. The model was successfully compared with experimental data obtained from an offshore oil platform located in Brazil. Moreover, apart from the main objectives, numerous results were obtained throughout this work. For instance, several diffusion coefficients and the relation between them are discussed. It is shown how the electric field set up between the ionic species depends on the gradient of the species’ concentrations. Furthermore, the capillary hysteresis effects are illustrated by a proposed model, which leads to the determination of several microstructure properties, such as the pore size distribution and the tortuosity-connectivity of the porous network. El fenómeno de corrosión del acero de refuerzo embebido en el hormigón ha tenido gran influencia en estructuras de hormigón armado, tanto en su funcionalidad estructural como en aspectos de durabilidad. La penetración de cloruros en el interior del hormigón esta considerada como el factor principal en el deterioro de la gran mayoría de estructuras. Por lo tanto, la modelización numérica de dicho fenómeno ha generado gran interés. El presente trabajo de investigación se centra en la modelización del transporte de cloruros en el interior del hormigón. Las ecuaciones diferenciales que gobiernan los fenómenos a nivel macroscópico se deducen de ecuaciones planteadas a nivel microscópico. Esto se obtiene comparando la red porosa con un poro equivalente, el cual mantiene las mismas propiedades de la red porosa real. El modelo está constituido por tres ecuaciones diferenciales acopladas que consideran el transporte de cloruros, el flujo de la solución de poro y la transferencia de calor. Con estas ecuaciones se tienen en cuenta los fenómenos de difusión, migración, advección, combinación y precipitación de cloruros. El análisis llevado a cabo en este trabajo ha definido los parámetros necesarios para calibrar el modelo. De acuerdo con ellas, se seleccionaron los ensayos experimentales a realizar. Las ecuaciones diferenciales se resolvieron mediante el método de elementos finitos. El método clásico de Galerkin se empleó para solucionar las ecuaciones de flujo de la solución de poro y de la transferencia de calor, mientras que el método streamline upwind Petrov-Galerkin se utilizó para resolver la ecuación de transporte de cloruros con la finalidad de evitar inestabilidades espaciales en problemas con advección dominante. El código de elementos finitos está implementado en Matlab® . Con el objetivo de facilitar la comprensión del grado de influencia de cada variable y parámetro, se realizó un análisis de sensibilidad detallado del modelo. Se llevó a cabo una campaña experimental sobre los hormigones estudiados, con el objeto de obtener sus propiedades difusivas, químicas e higroscópicas. El modelo se contrastó con datos experimentales obtenidos en una plataforma petrolera localizada en Brasil. Las simulaciones numéricas corroboraron los datos experimentales. Además, durante el desarrollo de la investigación se obtuvieron resultados paralelos a los planteados inicialmente. Por ejemplo, el análisis de diferentes coeficientes de difusión y la relación entre ellos. Así como también se observó que el campo eléctrico establecido entre las especies iónicas disueltas en la solución de poro depende del gradiente de concentración de las mismas. Los efectos de histéresis capilar son expresados por el modelo propuesto, el cual conduce a la determinación de una serie de propiedades microscópicas, tales como la distribución del tamaño de poro, además de la tortuosidad y conectividad de la red porosa.