Le présent travail de thèse met l'accent sur divers aspects de la distillation membranaire dans l'objectif de concevoir des procédés de dessalement proches du " zéro effluent liquide ". De manière générale, deux sujets sont discutés en détail: (i) la corrélation entre les caractéristiques de la membrane et les performances du procédé de distillation membranaire (ii) la compréhension et le contrôle de la polarisation thermique en DM. L'analyse de l'état de l'art en distillation membranaire porte notamment sur les progrès dans le développement des membranes, dans la compréhension des phénomènes de transport, les récents développements dans la conception des modules et le colmatage. Des phénomènes annexes et les applications innovantes sont également discutés dans la partie introductive de la thèse. L'effet des conditions de fabrication et de la composition des collodions sur les caractéristiques des membranes et la corrélation entre ces dernières et leurs performances a été discuté dans la section suivante. Il est établi que la morphologie de la membrane joue un rôle crucial dans ses performances pour des applications sur des fluides réels. En outre, on met en évidence que l'impact de la morphologie de la membrane est différente selon que la procédé fonctionne avec une phase liquide froide du côté distillat (Direct Contact Membrane Distillation - DCMD)) ou avec un courant d'air sec ou le vide (Air Gap ou Vacuum Membrane Distillation). Dans une deuxième partie, les aspects théoriques et expérimentaux de la polarisation thermique en distillation membranaire (DCMD) ont également été étudiés. Les phénomènes de polarisation thermique sur une membrane plane ont été étudiés en utilisant une cellule spécialement conçue. L'effet des conditions de fonctionnement et de la concentration de la solution sur la polarisation thermique a été étudié expérimentalement. Nous avons observé que l'augmentation de concentration de la solution favorise la polarisation thermique à cause d'une détérioration de l'hydrodynamique résultante à la surface de la membrane. Certaines techniques actives et passives pour réduire la polarisation thermique et le colmatage en distillation membranaire ont également été examinées dans l'étude cette étude. Nous avons montré que la polarisation thermique peut être considérablement réduite en générant des écoulements secondaires dans le fluide circulant à l'intérieur du canal d'alimentation, donc à l'intérieur de la fibre creuse si c'est cette configuration qui est retenue. Dans la présente étude, l'induction d'un écoulement secondaire a été réalisée en utilisant les fibres torsadées en hélice et une configuration ondulée. En raison de l'amélioration du niveau de polarisation thermique du côté de l'alimentation et du distillat, les géométries de fibres ondulées fournissent des flux et des taux de rendement supérieurs à ceux des autres configurations. La mise en œuvre d'un écoulement pulsé et intermittent pour contrôler la polarisation en distillation membranaire a également été examinée. Notre étude permet de conclure que ces types d'écoulements ont un impact positif sur les taux de rendement et le facteur d'amélioration volumique (gain en flux ramené par rapport à l'augmentation du volume de l'équipement) sans compromis sur le taux de remplissage des carters de fibres creuses. / The current PhD work emphasizes on various aspects of membrane distillation for approaching zero liquid discharge in seawater desalination. In broader sense, two themes have been discussed in detail: (i) correlation between membrane features and their performance in MD (ii) understanding and control of thermal polarization in MD. Introduction and state-of-the-art studies of MD including progress in membrane development, understanding the transport phenomenon, recent developments in module fabrication, fouling and related phenomenon and innovative applications have been discussed in introductory part of the thesis. The effect of operating conditions and dope compositions on membrane characteristics and correlation between membrane features and their performance has been discussed in subsequent section. It has been established that membrane morphology plays a crucial role in performance of the membrane for real applications. Furthermore, it has been demonstrated that the effect of membrane morphology is different for direct contact and vacuum configurations. Theoretical and experimental aspects of thermal polarization in direct contact membrane distillation have also been investigated. Thermal polarization phenomenon in a flat sheet membrane has been studied by using a specifically designed cell. The effect of operating conditions and solution concentration on thermal polarization has been explored experimentally. It has been observed that increased solution concentration favors the thermal polarization due to resulting poor hydrodynamic at the membrane surface and increase in diffusion resistance to the water vapors migrating from bulk feed phase to the membrane surface. Some active and passive techniques to decrease thermal polarization and possible fouling in membrane distillation have also been discussed in the current study. Thermal polarization can be greatly reduced by inducing secondary flows in the fluid flowing inside the fiber. The induction of secondary flows in the current study has been realized by using the fibers twisted in helical and wavy configurations. Due to improvement of thermal polarization coefficient on up and downstream, the undulating fiber geometries provide high flux and superior performance ratio. Application of intermittent and pulsatile flow to control thermal polarization in MD has also been discussed. It has been inferred that these flows have positive impact on performance ratio and volume based enhancement factors without compromising on packing density of the system. The application of MD for treatment of produced water has also been studied.

Orientador: Rosângela Barros Zanoni Lopes Moreno / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-20T08:21:26Z (GMT). No. of bitstreams: 1 Silveira_BrunoMarcodeOliveira_M.pdf: 1762377 bytes, checksum: cbdcb4c20dc00b869c37775090896cf4 (MD5) Previous issue date: 2012 / Resumo: As recentes descobertas do pré-sal no Brasil trouxeram novos desafios para as atividades de explotação de petróleo. Perfuração em formações carbonáticas apresentam grandes diferenças em relação a perfurações em rochas areníticas. Enquanto arenitos são mais homogeneos, rochas carbonáticas podem apresentar significativas variações em relação ao tamanho e à distribuição de poros, refletindo na produtividade de formação, como também na operação e manutenção de produção. A produtividade dos poços é avaliada logo após a perfuração dos mesmos, uma vez que danos à formação produtora podem alterar suas características originais. Diante disso, fluidos especificamente formulados para perfurar o reservatório, conhecidos como fluidos drill-in, vem sendo estudados com vistas à otimização da perfuração e à minimização de dano. Estes fluidos são desenvolvidos de acordo com as características do reservatório a ser perfurado e testes em laboratórios são necessários para que se possa adequar a formulação do fluido de perfuração com o menor custo, além de avaliar a invasão na rocha e verificar o possível dano residual para os fluidos candidatos. Este trabalho visa a comparação da invasão de fluidos drill-in e fluxo reverso de óleo em amostras de carbonatos e de arenitos de forma a avaliar a influência do tipo de rocha no dano à formação e sua interação com o fluido de perfuração testado. Os testes foram realizados usando dois aparatos experimentais diferentes, um para amostras com cerca de 7 cm de comprimento e outro para amostras medindo cerca de 30 cm, ambas com diâmetro de aproximadamente 3,7 cm. Os ensaios foram executados em amostras de carbonato e de arenito com permeabilidades absolutas ao gás semelhantes e submetendo-as ao respectivo protocolo de teste, cujas etapas incluíam a preparação das amostras, a caracterização de escoamento multifásico, a invasão com fluido de perfuração e o fluxo reverso de óleo. A invasão foi induzida por meio de um diferencial de pressão positivo através das amostras saturadas com óleo na condição de água conata, simulando o reservatório. O fluxo reverso ou deslocamento de óleo no sentido oposto ao de invasão representou a produção de óleo bem como a limpeza natural do reservatório. Durante a caracterização foi possível observar que as amostras de carbonatos utilizadas são menos molháveis à água em comparação aos arenitos. Com os resultados obtidos com o aparato utilizado para as amostras mais longas, foi possível identificar o deslocamento da frente de fluido invasor por meio dos pontos de medida de pressão distribuídos ao longo da amostra. A dinâmica de escoamento apontou que a invasão do fluido de perfuração foi mais rápida para as amostras de carbonatos em comparação com as amostras de arenito. Na etapa de fluxo reverso, pode-se observar que as amostras de carbonato apresentaram um retorno de produtividade maior e mais rápido relativamente aos arenitos / Abstract: The recent pre-salt discoveries in Brazil have brought new challenges to oil exploitation activities. Drilling in carbonate formations differs widely of drilling in sandstone rocks. While sandstones are more homogeneous, carbonate rocks are characterized by a large range of size and distribution of pores, which reflects in the reservoir productivity, as well as in the operation activities and production maintenance. The productivity and changes of original reservoir's characteristics due to formation damage are evaluated after drilling. Then fore, drilling fluids specifically designed to cross pay-zone, known as drill-in fluids, have been studied aiming at drilling optimization and damage minimizing. These fluids are developed according to the drilled reservoir characteristics and laboratory tests are performed to obtain proper drilling fluids formulations as fluid invasion and clean up are evaluated. The objective of this work is to compare drill-in fluids invasion (water-based polymer fluid) and backflow of oil in carbonates and sandstones samples to assess the rock type influence in damage formation and their interaction with the tested drilling fluid. Tests were run using two different apparatus, one for samples ~7,0 cm long and another for samples ~30 cm long, both with ~3,7 cm diameter. Lab tests were performed using carbonate and sandstone samples, with similar absolute permeability. These samples were submitted to a proper test protocol, that includes samples preparation, multiphase flow characterization, invasion and backflow processes. Damage formation was induced by injecting drill-in fluid with a positive differential pressure through samples saturated with oil and connate water, simulating the reservoir. Oil displaced in the opposite flow direction was performed to simulated natural cleanup or oil production. During the samples characterization, it was possible to observe that tested carbonate samples were less water wet than sandstones ones. From the results obtained with the apparatus used with longer samples, it was possible to identify the drill-in fluid displacement through pressure observation points distributed throughout the sample. The flow dynamics pointed out that the drilling fluid invasion was faster for carbonate samples comparing to sandstone samples. However, in backflow stage, a higher and faster productivity restoration was observed for carbonate samples than for sandstone ones / Mestrado / Reservatórios e Gestão / Mestre em Ciências e Engenharia de Petróleo

As possíveis soluções da equação integro-diferencial de Boltzmann constituem uma importante ferramenta para o estudo de gases e plasmas. No entanto, suas soluções analíticas são difíceis de serem encontradas. Uma abordagem bastante utilizada na literatura para obter soluções aproximadas da equação de Boltzmann é através de hipóteses que simplificam a forma da integral de colisão. Nesta dissertação, discutimos dois modelos colisionais alternativos que generalizam o método originalmente proposto por Bhatnagar, Gross e Krook, usualmente referido na literatura como aproximação BGK. O primeiro deles é um modelo de relaxação de segunda ordem, no qual introduzimos um segundo tempo de relaxação, 2 , relacionado com efeitos não lineares. O segundo modelo é baseado em outra generalização do modelo BGK obtida através de uma lei de potência parametrizada por um índice . No limite 1 o modelo BGK padrão é recuperado. As duas aproximações são fisicamente interpretadas. Além disso, para ilustrar nossos resultados com algumas aplicações mais quantitativas, obtemos as expressões analíticas para diversos coeficientes de transporte, a saber: a condutividade térmica (), a viscosidade de cisalhamento () e a condutividade elétrica (). Em particular, no modelo de relaxação de segunda ordem, as correções dependem da razão 1/2, onde 1 é a escala de tempo característica do modelo BGK padrão e 2 a nova escala de tempo associada aos efeitos não lineares. Finalmente, como um resultado geral, mostramos também que todas as correções nos coeficientes de transporte dependem numa certa potência do chamado número de Knudsen. / The possible solutions of the integro-diferential Boltzmann equation constitute an important tool for studying gases and plasmas. |However, its analytical solutions are hardly derived. An approach often adopted in the literature for obtaining approximate solutions of the Boltzmann equation is to consider some simplifying hypothesis on the collisional term. In this Dissertation, we discuss two diferent alternative collisional models which generalize the method originally proposed by Bhatnagar, Gross e Krook, and usually referred to as BGK approximation. The first one is a second order relaxation model in which a second relaxation time, 2, related with the nonlinear efects, is introduced. The second one is based on a diferent generalization of the BGK model which is obtained through a power law parameterized by a index . In the limit 1, the BGK model is recovered. Both approximations are physically interpreted. Further, in order to illustrate our results with some more quantitative applications, we derive the analytical expressions for several transport coefficients, among them the thermal conductivity (), the shear viscosity () and the electric conductivity (). In particular, for the second order relaxation model, we find that the corrections depend on the ratio 1/2 where 1 is the characteristic time scale of the BGK model and 2 describe the nonlinear efects. Finally, as a general result, it is also shown that all the corrections on the transport coefficients depend on a given power of the so-called Knudsen number.

As possíveis soluções da equação integro-diferencial de Boltzmann constituem uma importante ferramenta para o estudo de gases e plasmas. No entanto, suas soluções analíticas são difíceis de serem encontradas. Uma abordagem bastante utilizada na literatura para obter soluções aproximadas da equação de Boltzmann é através de hipóteses que simplificam a forma da integral de colisão. Nesta dissertação, discutimos dois modelos colisionais alternativos que generalizam o método originalmente proposto por Bhatnagar, Gross e Krook, usualmente referido na literatura como aproximação BGK. O primeiro deles é um modelo de relaxação de segunda ordem, no qual introduzimos um segundo tempo de relaxação, 2 , relacionado com efeitos não lineares. O segundo modelo é baseado em outra generalização do modelo BGK obtida através de uma lei de potência parametrizada por um índice . No limite 1 o modelo BGK padrão é recuperado. As duas aproximações são fisicamente interpretadas. Além disso, para ilustrar nossos resultados com algumas aplicações mais quantitativas, obtemos as expressões analíticas para diversos coeficientes de transporte, a saber: a condutividade térmica (), a viscosidade de cisalhamento () e a condutividade elétrica (). Em particular, no modelo de relaxação de segunda ordem, as correções dependem da razão 1/2, onde 1 é a escala de tempo característica do modelo BGK padrão e 2 a nova escala de tempo associada aos efeitos não lineares. Finalmente, como um resultado geral, mostramos também que todas as correções nos coeficientes de transporte dependem numa certa potência do chamado número de Knudsen. / The possible solutions of the integro-diferential Boltzmann equation constitute an important tool for studying gases and plasmas. |However, its analytical solutions are hardly derived. An approach often adopted in the literature for obtaining approximate solutions of the Boltzmann equation is to consider some simplifying hypothesis on the collisional term. In this Dissertation, we discuss two diferent alternative collisional models which generalize the method originally proposed by Bhatnagar, Gross e Krook, and usually referred to as BGK approximation. The first one is a second order relaxation model in which a second relaxation time, 2, related with the nonlinear efects, is introduced. The second one is based on a diferent generalization of the BGK model which is obtained through a power law parameterized by a index . In the limit 1, the BGK model is recovered. Both approximations are physically interpreted. Further, in order to illustrate our results with some more quantitative applications, we derive the analytical expressions for several transport coefficients, among them the thermal conductivity (), the shear viscosity () and the electric conductivity (). In particular, for the second order relaxation model, we find that the corrections depend on the ratio 1/2 where 1 is the characteristic time scale of the BGK model and 2 describe the nonlinear efects. Finally, as a general result, it is also shown that all the corrections on the transport coefficients depend on a given power of the so-called Knudsen number.

Mapas simpléticos têm sido amplamente utilizados para modelar o transporte caótico em plasmas e fluidos. Neste trabalho, propomos três tipos de mapas simpléticos que descrevem o movimento de deriva elétrica em plasmas magnetizados. Efeitos de raio de Larmor finito são incluídos em cada um dos mapas. No limite do raio de Larmor tendendo a zero, o mapa com frequência monotônica se reduz ao mapa de Chirikov-Taylor, e, nos casos com frequência não-monotônica, os mapas se reduzem ao mapa padrão não-twist. Mostramos como o raio de Larmor finito pode levar à supressão de caos, modificar a topologia do espaço de fases e a robustez de barreiras de transporte. Um método baseado na contagem dos tempos de recorrência é proposto para analisar a influência do raio de Larmor sobre os parâmetros críticos que definem a quebra de barreiras de transporte. Também estudamos um modelo para um sistema de partículas onde a deriva elétrica é descrita pelo mapa de frequência monotônica, e o raio de Larmor é uma variável aleatória que assume valores específicos para cada partícula do sistema. A função densidade de probabilidade para o raio de Larmor é obtida a partir da distribuição de Maxwell-Boltzmann, que caracteriza plasmas na condição de equilíbrio térmico. Um importante parâmetro neste modelo é a variável aleatória gama, definida pelo valor da função de Bessel de ordem zero avaliada no raio de Larmor da partícula. Resultados analíticos e numéricos descrevendo as principais propriedades estatísticas do parâmetro gama são apresentados. Tais resultados são então aplicados no estudo de duas medidas de transporte: a taxa de escape e a taxa de aprisionamento por ilhas de período um. / Area-preserving maps have been extensively used to model chaotic transport in plasmas and fluids. In this work we propose three types of maps describing electric drift motion in magnetized plasmas. Finite Larmor radius effects are included in all maps. In the limit of zero Larmor radius, the monotonic frequency map reduces to the Chirikov-Taylor map, and, in cases with non-monotonic frequencies, the maps reduce to the standard nontwist map. We show how the finite Larmor radius can lead to chaos suppression, modify the phase space topology and the robustness of transport barriers. A method based on counting the number of recurrence times is used to quantify the dependence on the Larmor radius of the threshold for the breakup of transport barriers. We also study a model for a system of particles where the electric drift is described by the monotonic frequency map, and the Larmor radius is a random variable that takes a specific value for each particle of the system. The Larmor radius\' probability density function is obtained from the Maxwell-Boltzmann distribution, which characterizes plasmas in thermal equilibrium. An important parameter in this model is the random variable gamma, defined by the zero-order Bessel function evaluated at the Larmor radius\'particle. We show analytical and numerical computations related to the statistics of gamma. The set of analytical results obtained here is then applied to the study of two numerical transport measures: the escape rate and the rate of trapping by period-one islands.

Mapas simpléticos têm sido amplamente utilizados para modelar o transporte caótico em plasmas e fluidos. Neste trabalho, propomos três tipos de mapas simpléticos que descrevem o movimento de deriva elétrica em plasmas magnetizados. Efeitos de raio de Larmor finito são incluídos em cada um dos mapas. No limite do raio de Larmor tendendo a zero, o mapa com frequência monotônica se reduz ao mapa de Chirikov-Taylor, e, nos casos com frequência não-monotônica, os mapas se reduzem ao mapa padrão não-twist. Mostramos como o raio de Larmor finito pode levar à supressão de caos, modificar a topologia do espaço de fases e a robustez de barreiras de transporte. Um método baseado na contagem dos tempos de recorrência é proposto para analisar a influência do raio de Larmor sobre os parâmetros críticos que definem a quebra de barreiras de transporte. Também estudamos um modelo para um sistema de partículas onde a deriva elétrica é descrita pelo mapa de frequência monotônica, e o raio de Larmor é uma variável aleatória que assume valores específicos para cada partícula do sistema. A função densidade de probabilidade para o raio de Larmor é obtida a partir da distribuição de Maxwell-Boltzmann, que caracteriza plasmas na condição de equilíbrio térmico. Um importante parâmetro neste modelo é a variável aleatória gama, definida pelo valor da função de Bessel de ordem zero avaliada no raio de Larmor da partícula. Resultados analíticos e numéricos descrevendo as principais propriedades estatísticas do parâmetro gama são apresentados. Tais resultados são então aplicados no estudo de duas medidas de transporte: a taxa de escape e a taxa de aprisionamento por ilhas de período um. / Area-preserving maps have been extensively used to model chaotic transport in plasmas and fluids. In this work we propose three types of maps describing electric drift motion in magnetized plasmas. Finite Larmor radius effects are included in all maps. In the limit of zero Larmor radius, the monotonic frequency map reduces to the Chirikov-Taylor map, and, in cases with non-monotonic frequencies, the maps reduce to the standard nontwist map. We show how the finite Larmor radius can lead to chaos suppression, modify the phase space topology and the robustness of transport barriers. A method based on counting the number of recurrence times is used to quantify the dependence on the Larmor radius of the threshold for the breakup of transport barriers. We also study a model for a system of particles where the electric drift is described by the monotonic frequency map, and the Larmor radius is a random variable that takes a specific value for each particle of the system. The Larmor radius\' probability density function is obtained from the Maxwell-Boltzmann distribution, which characterizes plasmas in thermal equilibrium. An important parameter in this model is the random variable gamma, defined by the zero-order Bessel function evaluated at the Larmor radius\'particle. We show analytical and numerical computations related to the statistics of gamma. The set of analytical results obtained here is then applied to the study of two numerical transport measures: the escape rate and the rate of trapping by period-one islands.

Cette thèse apporte des éléments sur la compréhension de la gestion de l'eau et de ses effets sur les performances électriques d'une PEMFC au moyen de modélisations multi-échelle des transferts. Une analyse du transport couplé de charges et de matière dans les pores de la membrane est proposée. La présence d'eau liquide est prise en compte dans les GDL (écoulements diphasiques) et les couches actives (noyage). Le couplage de ces modèles à une description des transferts de matière le long des canaux d’alimentation permet de mettre en évidence une répartition non-uniforme des concentrations en eau, des flux et donc de la densité de courant. Les résultats numériques sont comparés à des données expérimentales (coefficient de partage de l'eau et performance électrique locale) obtenues au laboratoire sur deux piles. Ceci permet de valider les modèles de fonctionnement du cœur de pile et d'alimenter la réflexion sur la connaissance et la modélisation des transferts d'eau dans le cœur de pile / This works contributes to the understanding of water management of polymer electrolyte membrane fuel cell and of its links with the electrical performances. More specifically, the manuscript deals with the multi-scale modelling of transport phenomena. An analysis of coupled mass and charge transfer in the pores of a polymer membrane is presented. The presence of liquid water is considered in the GDL (two-phase flow) and in the active layers (flooding). The description of these phenomena is associated with that of gas depletion along the bipolar plate channels. This allows to emphasize the non-uniformity of water concentration, of the fluxes and as a consequence, of current density. The numerical results are compared with experimental data (water transport coefficient, local electrical performances) measured on two different fuel cells. This comparison validates at least partially the numerical models and provides further information for the analysis of water management within PEMFC

本研究旨在應用書目計量學方法分析1900年至2007年間，輸送現象文獻之特性，藉以探討輸送現象之發展情形。研究樣本取自SCIE資料庫，計104854篇相關書目資料。期望本研究結果能提供輸送現象學界與圖書館或資料單位參考，協助其做有效之館藏規劃及期刊資訊服務。 本研究結果歸納如下：(1)輸送現象文獻成長狀況目前呈穩定成長狀態，1965年後之文獻累積成長曲線符合指數成長。(2)輸送現象文獻之主要出版國為美國、英國、荷蘭、瑞士、德國、日本等國。台灣文獻出版量之排名則為第22名。(3)英文為輸送現象文獻最重要之寫作語文，德文、法文、俄文、日文依序於後。(4)輸送現象文獻之出版類型以期刊文章為主，佔總文獻之91.37%。(5)本研究不符合布萊德福定律，但符合布萊德福-齊夫定律曲線，求得之核心期刊有八種，分別為美國及英國所出版。(6)以高被引期刊、影響係數以及被資料庫收錄情況等方式，驗證布萊德福分區法所求之核心期刊品質，發現核心期刊品質值得信任。(7)以洛卡定律最小平方法求得n值為-2.5141(平等法)、-2.5389(第一作者法)，c值為0.7488(平等法)、0.7546(第一作者法)，進一步利用K-S檢定法加以檢定，證明洛卡定律並不適用於本研究。(8)本研究之作者生產力分佈不符合普萊斯平方根定律與80/20定律。(9)印度理工學院為輸送現象文獻重要的生產機構；台灣則為國立成功大學。總生產力台灣排名世界第22。(10)美國、法國、德國、日本、俄羅斯為相關研究機構最多的國家；英國、美國、荷蘭、瑞士則為最具吸收他國研究成果的學術傳播環境。 本研究最後依據研究發現之結果，提出以下建議：(1)本研究之各項研究成果，可作為圖書資訊界館藏發展與管理之參考。(2)輸送現象產業界應加強對於研究的重視。(3)輸送現象相關研究人員，應具備英文、德文、法文、俄文等外語能力。(4)資料庫公司應提升書目資料精確性，以利研究。(5)定期對輸送現象文獻進行分析，有助於輸送現象相關研究人員瞭解該領域之研究發展。 / This study analyzed the characteristics of transport phenomenon literature during 1900-2007 by a bibliometric approach. A total of 104854 bibliographic data were retrieved from SCIE database. The results attempt to have a better understanding of the developments and trends on transport phenomenon reaserch. The results of the present study reveal that: (1) Transport Phenomenon literature increased steadily. The growth curve during 1965-2007 fits exponential growth. (2)USA, UK, Netherlands, Switzerland, Germany and Japan are the major countries of Transport Phenomenon publications. Taiwan ranked 22. (3)English is the most often used language. German ranked number two and the third most often used language is French. (4)Journal article is the major of publication and accounted for 91.37% of total publication. (5)The journal article distribution doesn't fit Bradford Law. However, it fits the typical Bradford-Zipf’s S-shaped curve. (6)Core journals are found to be associated with highly cited times and high impact factor. (7)For only counting first author, Lotka’s law is found to be inapplicable to author productivity distribution with n=-2.5389, c=0.7488. The K-S test also explored the invalidity of observed distribution at 0.01 level of significance. (8) Author productivity distribution also does’t fit Price Square Root Law and 80/20 Rule. (9) Indian Institute of Technology is the most important institute of Transport Phenomenon publications. In Taiwan, National Cheng Kung University published the most Transport Phenomenon literature and it ranked 22. (10)USA, French, Germany, Japan and Russia have the most relevant research institutions than other countries.