用光学干涉方法对液相扩散传质过程的研究

液相中扩散传质过程研究无论对于基础理论还是生产实践都具有重要的意义.但在通常的重力环境中,传质过程不是单因素地由浓度梯度来决定,对流和沉降会对实验研究产生重要的干扰.通过光学干涉技术的应用,实现对于扩散过程中的传质系数进行测量,实验中采用Math-Zehnder干涉仪对于整个传质过程进行监控并记录相关图像信息.这些图像信息是贯穿在整个实验过程中的连续录像,不同于其它实验记录的静态图像.通过计算这些随时间改变的干涉条纹的变化,就可以推导出传质系数的结果.实验液体选用水/葡萄糖溶液,之后还将进一步将该实验装置搭载TF-1火箭,进行微重力实验,以排除重力产生的影响.

Investigation on mechanism of critical cavitating flow in liquid jet pumps under operating limits

The critical cavitating flow in liquid jet pumps under operating limits is investigated in this paper. Measurements on the axial pressure distribution along the wall of jet pumps indicate that two-phase critical flow occurs in the throat pipe under operating limits. The entrained flow rate and the distribution of the wall pressure upstream lowest pressure section does not change when the outlet pressure is lower than a critical value. A liquid-vapor mixing shockwave is also observed under operating limits. The wave front moves back and forth in low frequency around the position of the lowest pressure. With the measured axial wall pressures, the Mach number of the two-phase cavitating flow is calculated. It's found that the maximum Mach number is very close to I under operating limits. Further analysis infers a cross-section where Mach number approaches to I near the wave front. Thus, the liquid-vapor mixture velocity should reach the local sound velocity and resulting in the occurrence of operating limits.

氨热法生长氮化镓晶体中传热传质的研究

本文利用基于非正交网格的二阶精度有限体积法,对氨热法生长过程中温度场和流场进行了模拟,其中隔板开孔率分别为10%(中心开孔5%,侧壁与隔板边缘开孔5%)和20%(中心开孔10%,侧壁与隔板边缘开孔10%).通过对流场和温度场的分析,了解了高压釜内部营养素的输运及溶液结晶的本质.结果显示在釜底的多孔介质层流动较弱,在流体层流动较强.在多孔介质层热量的输运主要通过热传导;在流体层中,流体与原料的分界处以及流体和高压釜的侧壁出现了大的温度梯度.

I. The effects of variations in certain parameters upon the efficiency of in vivo hemodialysis with a Kiil dialyzer. II. Interfacial structure of concurrent air-water flow in a two-inch diameter horizontal tube.

Part I

Regression analyses are performed on in vivo hemodialysis data for the transfer of creatinine, urea, uric acid and inorganic phosphate to determine the effects of variations in certain parameters on the efficiency of dialysis with a Kiil dialyzer. In calculating the mass transfer rates across the membrane, the effects of cell-plasma mass transfer kinetics are considered. The concept of the effective permeability coefficient for the red cell membrane is introduced to account for these effects. A discussion of the consequences of neglecting cell-plasma kinetics, as has been done to date in the literature, is presented.

A physical model for the Kiil dialyzer is presented in order to calculate the available membrane area for mass transfer, the linear blood and dialysate velocities, and other variables. The equations used to determine the independent variables of the regression analyses are presented. The potential dependent variables in the analyses are discussed.

Regression analyses were carried out considering overall mass-transfer coefficients, dialysances, relative dialysances, and relative permeabilities for each substance as the dependent variables. The independent variables were linear blood velocity, linear dialysate velocity, the pressure difference across the membrane, the elapsed time of dialysis, the blood hematocrit, and the arterial plasma concentrations of each substance transferred. The resulting correlations are tabulated, presented graphically, and discussed. The implications of these correlations are discussed from the viewpoint of a research investigator and from the viewpoint of patient treatment.

Recommendations for further experimental work are presented.

Part II

The interfacial structure of concurrent air-water flow in a two-inch diameter horizontal tube in the wavy flow regime has been measured using resistance wave gages. The median water depth, r.m.s. wave height, wave frequency, extrema frequency, and wave velocity have been measured as functions of air and water flow rates. Reynolds numbers, Froude numbers, Weber numbers, and bulk velocities for each phase may be calculated from these measurements. No theory for wave formation and propagation available in the literature was sufficient to describe these results.

The water surface level distribution generally is not adequately represented as a stationary Gaussian process. Five types of deviation from the Gaussian process function were noted in this work. The presence of the tube walls and the relatively large interfacial shear stresses precludes the use of simple statistical analyses to describe the interfacial structure. A detailed study of the behavior of individual fluid elements near the interface may be necessary to describe adequately wavy two-phase flow in systems similar to the one used in this work.

I. The attenuation and dispersion of sound in a condensing medium. II. The flow of a gas-particle mixture downstream of a normal shock in a nozzle

I. The attenuation of sound due to particles suspended in a gas was first calculated by Sewell and later by Epstein in their classical works on the propagation of sound in a two-phase medium. In their work, and in more recent works which include calculations of sound dispersion, the calculations were made for systems in which there was no mass transfer between the two phases. In the present work, mass transfer between phases is included in the calculations.

The attenuation and dispersion of sound in a two-phase condensing medium are calculated as functions of frequency. The medium in which the sound propagates consists of a gaseous phase, a mixture of inert gas and condensable vapor, which contains condensable liquid droplets. The droplets, which interact with the gaseous phase through the interchange of momentum, energy, and mass (through evaporation and condensation), are treated from the continuum viewpoint. Limiting cases, for flow either frozen or in equilibrium with respect to the various exchange processes, help demonstrate the effects of mass transfer between phases. Included in the calculation is the effect of thermal relaxation within droplets. Pressure relaxation between the two phases is examined, but is not included as a contributing factor because it is of interest only at much higher frequencies than the other relaxation processes. The results for a system typical of sodium droplets in sodium vapor are compared to calculations in which there is no mass exchange between phases. It is found that the maximum attenuation is about 25 per cent greater and occurs at about one-half the frequency for the case which includes mass transfer, and that the dispersion at low frequencies is about 35 per cent greater. Results for different values of latent heat are compared.

II. In the flow of a gas-particle mixture through a nozzle, a normal shock may exist in the diverging section of the nozzle. In Marble’s calculation for a shock in a constant area duct, the shock was described as a usual gas-dynamic shock followed by a relaxation zone in which the gas and particles return to equilibrium. The thickness of this zone, which is the total shock thickness in the gas-particle mixture, is of the order of the relaxation distance for a particle in the gas. In a nozzle, the area may change significantly over this relaxation zone so that the solution for a constant area duct is no longer adequate to describe the flow. In the present work, an asymptotic solution, which accounts for the area change, is obtained for the flow of a gas-particle mixture downstream of the shock in a nozzle, under the assumption of small slip between the particles and gas. This amounts to the assumption that the shock thickness is small compared with the length of the nozzle. The shock solution, valid in the region near the shock, is matched to the well known small-slip solution, which is valid in the flow downstream of the shock, to obtain a composite solution valid for the entire flow region. The solution is applied to a conical nozzle. A discussion of methods of finding the location of a shock in a nozzle is included.

Modelling algal blooms

This is the final report on the research project to develop predictive models to quantify algal blooms in relation to environmental variables. The project's objectives were to develop models simulating the impact of vertical structure and mass transfer upon the dynamics of planktonic algae, including cyanobacteria, in lakes and reservoirs, to assess the potential of sedimentary phosphorus to sustain algal growth following reduction in external loading and to expand and enhance formulations to predict behaviour of blue-green algal populations and to incorporate these into a model software package. As part of the project a strategy for the production of a user-friendly packaging for the modelling software PROTEC-2 adaptable to particular sites was developed.

Advanced Reaction Systems for Hydrogen Production

[EN] This PhD work started in March 2010 with the support of the University of the Basque Country (UPV/EHU) under the program named “Formación de Personal Investigador” at the Chemical and Environmental Engineering Department in the Faculty of Engineering of Bilbao. The major part of the Thesis work was carried out in the mentioned department, as a member of the Sustainable Process Engineering (SuPrEn) research group. In addition, this PhD Thesis includes the research work developed during a period of 6 months at the Institut für Mikrotechnik Mainz GmbH, IMM, in Germany. During the four years of the Thesis, conventional and microreactor systems were tested for several feedstocks renewable and non-renewable, gases and liquids through several reforming processes in order to produce hydrogen. For this purpose, new catalytic formulations which showed high activity, selectivity and stability were design. As a consequence, the PhD work performed allowed the publication of seven scientific articles in peer-reviewed journals. This PhD Thesis is divided into the following six chapters described below. The opportunity of this work is established on the basis of the transition period needed for moving from a petroleum based energy system to a renewable based new one. Consequently, the present global energy scenario was detailed in Chapter 1, and the role of hydrogen as a real alternative in the future energy system was justified based on several outlooks. Therefore, renewable and non-renewable hydrogen production routes were presented, explaining the corresponding benefits and drawbacks. Then, the raw materials used in this Thesis work were described and the most important issues regarding the processes and the characteristics of the catalytic formulations were explained. The introduction chapter finishes by introducing the concepts of decentralized production and process intensification with the use of microreactors. In addition, a small description of these innovative reaction systems and the benefits that entailed their use were also mentioned. In Chapter 2 the main objectives of this Thesis work are summarized. The development of advanced reaction systems for hydrogen rich mixtures production is the main objective. In addition, the use and comparison between two different reaction systems, (fixed bed reactor (FBR) and microreactor), the processing of renewable raw materials, the development of new, active, selective and stable catalytic formulations, and the optimization of the operating conditions were also established as additional partial objectives. Methane and natural gas (NG) steam reforming experimental results obtained when operated with microreactor and FBR systems are presented in Chapter 3. For these experiments nickel-based (Ni/Al2O3 and Ni/MgO) and noble metal-based (Pd/Al2O3 and Pt/Al2O3) catalysts were prepared by wet impregnation and their catalytic activity was measured at several temperatures, from 973 to 1073 K, different S/C ratios, from 1.0 to 2.0, and atmospheric pressure. The Weight Hourly Space Velocity (WHSV) was maintained constant in order to compare the catalytic activity in both reaction systems. The results obtained showed a better performance of the catalysts operating in microreactors. The Ni/MgO catalyst reached the highest hydrogen production yield at 1073 K and steam-to-carbon ratio (S/C) of 1.5 under Steam methane Reforming (SMR) conditions. In addition, this catalyst also showed good activity and stability under NG reforming at S/C=1.0 and 2.0. The Ni/Al2O3 catalyst also showed high activity and good stability and it was the catalyst reaching the highest methane conversion (72.9 %) and H2out/CH4in ratio (2.4) under SMR conditions at 1073 K and S/C=1.0. However, this catalyst suffered from deactivation when it was tested under NG reforming conditions. Regarding the activity measurements carried out with the noble metal-based catalysts in the microreactor systems, they suffered a very quick deactivation, probably because of the effects attributed to carbon deposition, which was detected by Scanning Electron Microscope (SEM). When the FBR was used no catalytic activity was measured with the catalysts under investigation, probably because they were operated at the same WHSV than the microreactors and these WHSVs were too high for FBR system. In Chapter 4 biogas reforming processes were studied. This chapter starts with an introduction explaining the properties of the biogas and the main production routes. Then, the experimental procedure carried out is detailed giving concrete information about the experimental set-up, defining the parameters measured, specifying the characteristics of the reactors used and describing the characterization techniques utilized. Each following section describes the results obtained from activity testing with the different catalysts prepared, which is subsequently summarized: Section 4.3: Biogas reforming processes using γ-Al2O3 based catalysts The activity results obtained by several Ni-based catalysts and a bimetallic Rh-Ni catalyst supported on magnesia or alumina modified with oxides like CeO2 and ZrO2 are presented in this section. In addition, an alumina-based commercial catalyst was tested in order to compare the activity results measured. Four different biogas reforming processes were studied using a FBR: dry reforming (DR), biogas steam reforming (BSR), biogas oxidative reforming (BOR) and tri-reforming (TR). For the BSR process different steam to carbon ratios (S/C) from 1.0 to 3.0, were tested. In the case of BOR process the oxygen-to-methane (O2/CH4) ratio was varied from 0.125 to 0.50. Finally, for TR processes different S/C ratios from 1.0 to 3.0, and O2/CH4 ratios of 0.25 and 0.50 were studied. Then, the catalysts which achieved high activity and stability were impregnated in a microreactor to explore the viability of process intensification. The operation with microreactors was carried out under the best experimental conditions measured in the FBR. In addition, the physicochemical characterization of the fresh and spent catalysts was carried out by Inductively Coupled Plasma Atomic Emission Spectroscopy (ICP-AES), N2 physisorption, H2 chemisorption, Temperature Programmed Reduction (TPR), SEM, X-ray Photoelectron Spectroscopy (XPS) and X-ray powder Diffraction (XRD). Operating with the FBR, conversions close to the ones predicted by thermodynamic calculations were obtained by most of the catalysts tested. The Rh-Ni/Ce-Al2O3 catalyst obtained the highest hydrogen production yield in DR. In BSR process, the Ni/Ce-Al2O3 catalyst achieved the best activity results operating at S/C=1.0. In the case of BOR process, the Ni/Ce-Zr-Al2O3 catalyst showed the highest reactants conversion values operating at O2/CH4=0.25. Finally, in the TR process the Rh-Ni/Ce-Al2O3 catalyst obtained the best results operating at S/C=1.0 and O2/CH4=0.25. Therefore, these three catalysts were selected to be coated onto microchannels in order to test its performance under BOR and TR processes conditions. Although the operation using microreactors was carried out under considerably higher WHSV, similar conversions and yields as the ones measured in FBR were measured. Furthermore, attending to other measurements like Turnover Frequency (TOF) and Hydrogen Productivity (PROD), the values calculated for the catalysts tested in microreactors were one order of magnitude higher. Thus, due to the low dispersion degree measured by H2-chemisorption, the Ni/Ce-Al2O3 catalyst reached the highest TOF and PROD values. Section 4.4: Biogas reforming processes using Zeolites L based catalysts In this section three type of L zeolites, with different morphology and size, were synthesized and used as catalyst support. Then, for each type of L zeolite three nickel monometallic and their homologous Rh-Ni bimetallic catalysts were prepared by the wetness impregnation method. These catalysts were tested using the FBR under DR process and different conditions of BSR (S/C ratio of 1.0 and 2.0), BOR (O2/CH4 ratio of 0.25 and 0.50) and TR processes (at S/C=1.0 and O2/CH4=0.25). The characterization of these catalysts was also carried out by using the same techniques mentioned in the previous section. Very high methane and carbon dioxide conversion values were measured for almost all the catalysts under investigation. The experimental results evidenced the better catalytic behavior of the bimetallic catalysts as compared to the monometallic ones. Comparing the catalysts behavior with regards to their morphology, for the BSR process the Disc catalysts were the most active ones at the lowest S/C ratio tested. On the contrary, the Cylindrical (30–60 nm) catalysts were more active under BOR conditions at O2/CH4=0.25 and TR processes. By the contrary, the Cylindrical (1–3 µm) catalysts showed the worst activity results for both processes. Section 4.5: Biogas reforming processes using Na+ and Cs+ doped Zeolites LTL based catalysts A method for the synthesis of Linde Type L (LTL) zeolite under microwave-assisted hydrothermal conditions and its behavior as a support for heterogeneously catalyzed hydrogen production is described in this section. Then, rhodium and nickel-based bimetallic catalysts were prepared in order to be tested by DR process and BOR process at O2/CH4=0.25. Moreover, the characterization of the catalysts under investigation was also carried out. Higher activities were achieved by the catalysts prepared from the non-doped zeolites, Rh-Ni/D and Rh-Ni/N, as compared to the ones supported on Na+ and Cs+ exchanged supports. However, the differences between them were not very significant. In addition, the Na+ and Cs+ incorporation affected mainly to the Disc catalysts. Comparing the results obtained by these catalysts with the ones studied in the section 4.4, in general worst results were achieved under DR conditions and almost the same results when operated under BOR conditions. In Chapter 5 the ethylene glycol (EG) as feed for syngas production by steam reforming (SR) and oxidative steam reforming (OSR) was studied by using microchannel reactors. The product composition was determined at a S/C of 4.0, reaction temperatures between 625°C and 725°C, atmospheric pressure and Volume Hourly Space Velocities (VHSV) between 100 and 300 NL/(gcath). This work was divided in two sections. The first one corresponds to the introduction of the main and most promising EG production routes. Then, the new experimental procedure is detailed and the information about the experimental set-up and the measured parameters is described. The characterization was carried out using the same techniques as for the previous chapter. Then, the next sections correspond to the catalytic activity and catalysts characterization results. Section 5.3: xRh-cm and xRh-np catalysts for ethylene glycol reforming Initially, catalysts with different rhodium loading, from 1.0 to 5.0 wt. %, and supported on α-Al2O3 were prepared by two different preparation methods (conventional impregnation and separate nanoparticle synthesis). Then, the catalysts were compared regarding their measured activity and selectivity, as well as the characterization results obtained before and after the activity tests carried out. The samples prepared by a conventional impregnation method showed generally higher activity compared to catalysts prepared from Rh nanoparticles. By-product formation of species such as acetaldehyde, ethane and ethylene was detected, regardless if oxygen was added to the feed or not. Among the catalysts tested, the 2.5Rh-cm catalyst was considered the best one. Section 5.4: 2.5Rh-cm catalyst support modification with CeO2 and La2O3 In this part of the Chapter 5, the catalyst showing the best performance in the previous section, the 2.5Rh-Al2O3 catalyst, was selected in order to be improved. Therefore, new Rh based catalysts were designed using α-Al2O3 and being modified this support with different contents of CeO2 or La2O3 oxides. All the catalysts containing additives showed complete conversion and selectivities close to the equilibrium in both SR and OSR processes. In addition, for these catalysts the concentrations measured for the C2H4, CH4, CH3CHO and C2H6 by-products were very low. Finally, the 2.5Rh-20Ce catalyst was selected according to its catalytic activity and characterization results in order to run a stability test, which lasted more than 115 hours under stable operation. The last chapter, Chapter 6, summarizes the main conclusions achieved throughout this Thesis work. Although very high reactant conversions and rich hydrogen mixtures were obtained using a fixed bed reaction system, the use of microreactors improves the key issues, heat and mass transfer limitations, through which the reforming reactions are intensified. Therefore, they seem to be a very interesting and promising alternative for process intensification and decentralized production for remote application.

Desenvolvimento e avaliação de novas abordagens de modelagem de processos de separação em leito móvel simulado

O Leito Móvel Simulado (LMS) é um processo de separação de compostos por adsorção muito eficiente, por trabalhar em um regime contínuo e também possuir fluxo contracorrente da fase sólida. Dentre as diversas aplicações, este processo tem se destacado na resolução de petroquímicos e principalmente na atualidade na separação de misturas racêmicas que são separações de um grau elevado de dificuldade. Neste trabalho foram propostas duas novas abordagens na modelagem do LMS, a abordagem Stepwise e a abordagem Front Velocity. Na modelagem Stepwise as colunas cromatográficas do LMS foram modeladas com uma abordagem discreta, onde cada uma delas teve seu domínio dividido em N células de mistura interligadas em série, e as concentrações dos compostos nas fases líquida e sólida foram simuladas usando duas cinéticas de transferência de massa distintas. Essa abordagem pressupõe que as interações decorrentes da transferência de massa entre as moléculas do composto nas suas fases líquida e sólida ocorram somente na superfície, de forma que com essa suposição pode-se admitir que o volume ocupado por cada molécula nas fases sólida e líquida é o mesmo, o que implica que o fator de residência pode ser considerado igual a constante de equilíbrio. Para descrever a transferência de massa que ocorre no processo cromatográfico a abordagem Front Velocity estabelece que a convecção é a fase dominante no transporte de soluto ao longo da coluna cromatográfica. O Front Velocity é um modelo discreto (etapas) em que a vazão determina o avanço da fase líquida ao longo da coluna. As etapas são: avanço da fase líquida e posterior transporte de massa entre as fases líquida e sólida, este último no mesmo intervalo de tempo. Desta forma, o fluxo volumétrico experimental é utilizado para a discretização dos volumes de controle que se deslocam ao longo da coluna porosa com a mesma velocidade da fase líquida. A transferência de massa foi representada por dois mecanismos cinéticos distintos, sem (tipo linear) e com capacidade máxima de adsorção (tipo Langmuir). Ambas as abordagens propostas foram estudadas e avaliadas mediante a comparação com dados experimentais de separação em LMS do anestésico cetamina e, posteriormente, com o fármaco Verapamil. Também foram comparados com as simulações do modelo de equilíbrio dispersivo para o caso da Cetamina, usado por Santos (2004), e para o caso do Verapamil (Perna 2013). Na etapa de caracterização da coluna cromatográfica as novas abordagens foram associadas à ferramenta inversa R2W de forma a determinar os parâmetros globais de transferência de massa apenas usando os tempos experimentais de residência de cada enantiômero na coluna de cromatografia líquida de alta eficiência (CLAE). Na segunda etapa os modelos cinéticos desenvolvidos nas abordagens foram aplicados nas colunas do LMS com os valores determinados na caracterização da coluna cromatográfica, para a simulação do processo de separação contínua. Os resultados das simulações mostram boa concordância entre as duas abordagens propostas e os experimentos de pulso para a caracterização da coluna na separação enantiomérica da cetamina ao longo do tempo. As simulações da separação em LMS, tanto do Verapamil quando da Cetamina apresentam uma discrepância com os dados experimentais nos primeiros ciclos, entretanto após esses ciclos iniciais a correlação entre os dados experimentais e as simulações. Para o caso da separação da cetamina (Santos, 2004), a qual a concentração da alimentação era relativamente baixa, os modelos foram capazes de predizer o processo de separação com as cinéticas Linear e Langmuir. No caso da separação do Verapamil (Perna, 2013), onde a concentração da alimentação é relativamente alta, somente a cinética de Langmuir representou o processo, devido a cinética Linear não representar a saturação das colunas cromatográficas. De acordo como o estudo conduzido ambas as abordagens propostas mostraram-se ferramentas com potencial na predição do comportamento cromatográfico de uma amostra em um experimento de pulso, assim como na simulação da separação de um composto no LMS, apesar das pequenas discrepâncias apresentadas nos primeiros ciclos de trabalho do LMS. Além disso, podem ser facilmente implementadas e aplicadas na análise do processo, pois requer um baixo número de parâmetros e são constituídas de equações diferenciais ordinárias.