Photon Management Structures for Absorption Enhancement in Intermediate Band Solar Cells and Crystalline Silicon Solar Cells

El objetivo de la tesis es investigar los beneficios que el atrapamiento de la luz mediante fenómenos difractivos puede suponer para las células solares de silicio cristalino y las de banda intermedia. Ambos tipos de células adolecen de una insuficiente absorción de fotones en alguna región del espectro solar. Las células solares de banda intermedia son teóricamente capaces de alcanzar eficiencias mucho mayores que los dispositivos convencionales (con una sola banda energética prohibida), pero los prototipos actuales se resienten de una absorción muy débil de los fotones con energías menores que la banda prohibida. Del mismo modo, las células solares de silicio cristalino absorben débilmente en el infrarrojo cercano debido al carácter indirecto de su banda prohibida. Se ha prestado mucha atención a este problema durante las últimas décadas, de modo que todas las células solares de silicio cristalino comerciales incorporan alguna forma de atrapamiento de luz. Por razones de economía, en la industria se persigue el uso de obleas cada vez más delgadas, con lo que el atrapamiento de la luz adquiere más importancia. Por tanto aumenta el interés en las estructuras difractivas, ya que podrían suponer una mejora sobre el estado del arte. Se comienza desarrollando un método de cálculo con el que simular células solares equipadas con redes de difracción. En este método, la red de difracción se analiza en el ámbito de la óptica física, mediante análisis riguroso con ondas acopladas (rigorous coupled wave analysis), y el sustrato de la célula solar, ópticamente grueso, se analiza en los términos de la óptica geométrica. El método se ha implementado en ordenador y se ha visto que es eficiente y da resultados en buen acuerdo con métodos diferentes descritos por otros autores. Utilizando el formalismo matricial así derivado, se calcula el límite teórico superior para el aumento de la absorción en células solares mediante el uso de redes de difracción. Este límite se compara con el llamado límite lambertiano del atrapamiento de la luz y con el límite absoluto en sustratos gruesos. Se encuentra que las redes biperiódicas (con geometría hexagonal o rectangular) pueden producir un atrapamiento mucho mejor que las redes uniperiódicas. El límite superior depende mucho del periodo de la red. Para periodos grandes, las redes son en teoría capaces de alcanzar el máximo atrapamiento, pero sólo si las eficiencias de difracción tienen una forma peculiar que parece inalcanzable con las herramientas actuales de diseño. Para periodos similares a la longitud de onda de la luz incidente, las redes de difracción pueden proporcionar atrapamiento por debajo del máximo teórico pero por encima del límite Lambertiano, sin imponer requisitos irrealizables a la forma de las eficiencias de difracción y en un margen de longitudes de onda razonablemente amplio. El método de cálculo desarrollado se usa también para diseñar y optimizar redes de difracción para el atrapamiento de la luz en células solares. La red propuesta consiste en un red hexagonal de pozos cilíndricos excavados en la cara posterior del sustrato absorbente de la célula solar. La red se encapsula en una capa dieléctrica y se cubre con un espejo posterior. Se simula esta estructura para una célula solar de silicio y para una de banda intermedia y puntos cuánticos. Numéricamente, se determinan los valores óptimos del periodo de la red y de la profundidad y las dimensiones laterales de los pozos para ambos tipos de células. Los valores se explican utilizando conceptos físicos sencillos, lo que nos permite extraer conclusiones generales que se pueden aplicar a células de otras tecnologías. Las texturas con redes de difracción se fabrican en sustratos de silicio cristalino mediante litografía por nanoimpresión y ataque con iones reactivos. De los cálculos precedentes, se conoce el periodo óptimo de la red que se toma como una constante de diseño. Los sustratos se procesan para obtener estructuras precursoras de células solares sobre las que se realizan medidas ópticas. Las medidas de reflexión en función de la longitud de onda confirman que las redes cuadradas biperiódicas consiguen mejor atrapamiento que las uniperiódicas. Las estructuras fabricadas se simulan con la herramienta de cálculo descrita en los párrafos precedentes y se obtiene un buen acuerdo entre la medida y los resultados de la simulación. Ésta revela que una fracción significativa de los fotones incidentes son absorbidos en el reflector posterior de aluminio, y por tanto desaprovechados, y que este efecto empeora por la rugosidad del espejo. Se desarrolla un método alternativo para crear la capa dieléctrica que consigue que el reflector se deposite sobre una superficie plana, encontrándose que en las muestras preparadas de esta manera la absorción parásita en el espejo es menor. La siguiente tarea descrita en la tesis es el estudio de la absorción de fotones en puntos cuánticos semiconductores. Con la aproximación de masa efectiva, se calculan los niveles de energía de los estados confinados en puntos cuánticos de InAs/GaAs. Se emplea un método de una y de cuatro bandas para el cálculo de la función de onda de electrones y huecos, respectivamente; en el último caso se utiliza un hamiltoniano empírico. La regla de oro de Fermi permite obtener la intensidad de las transiciones ópticas entre los estados confinados. Se investiga el efecto de las dimensiones del punto cuántico en los niveles de energía y la intensidad de las transiciones y se obtiene que, al disminuir la anchura del punto cuántico respecto a su valor en los prototipos actuales, se puede conseguir una transición más intensa entre el nivel intermedio fundamental y la banda de conducción. Tomando como datos de partida los niveles de energía y las intensidades de las transiciones calculados como se ha explicado, se desarrolla un modelo de equilibrio o balance detallado realista para células solares de puntos cuánticos. Con el modelo se calculan las diferentes corrientes debidas a transiciones ópticas entre los numerosos niveles intermedios y las bandas de conducción y de valencia bajo ciertas condiciones. Se distingue de modelos de equilibrio detallado previos, usados para calcular límites de eficiencia, en que se adoptan suposiciones realistas sobre la absorción de fotones para cada transición. Con este modelo se reproducen datos publicados de eficiencias cuánticas experimentales a diferentes temperaturas con un acuerdo muy bueno. Se muestra que el conocido fenómeno del escape térmico de los puntos cuánticos es de naturaleza fotónica; se debe a los fotones térmicos, que inducen transiciones entre los estados excitados que se encuentran escalonados en energía entre el estado intermedio fundamental y la banda de conducción. En el capítulo final, este modelo realista de equilibrio detallado se combina con el método de simulación de redes de difracción para predecir el efecto que tendría incorporar una red de difracción en una célula solar de banda intermedia y puntos cuánticos. Se ha de optimizar cuidadosamente el periodo de la red para equilibrar el aumento de las diferentes transiciones intermedias, que tienen lugar en serie. Debido a que la absorción en los puntos cuánticos es extremadamente débil, se deduce que el atrapamiento de la luz, por sí solo, no es suficiente para conseguir corrientes apreciables a partir de fotones con energía menor que la banda prohibida en las células con puntos cuánticos. Se requiere una combinación del atrapamiento de la luz con un incremento de la densidad de puntos cuánticos. En el límite radiativo y sin atrapamiento de la luz, se necesitaría que el número de puntos cuánticos de una célula solar se multiplicara por 1000 para superar la eficiencia de una célula de referencia con una sola banda prohibida. En cambio, una célula con red de difracción precisaría un incremento del número de puntos en un factor 10 a 100, dependiendo del nivel de la absorción parásita en el reflector posterior. Abstract The purpose of this thesis is to investigate the benefits that diffractive light trapping can offer to quantum dot intermediate band solar cells and crystalline silicon solar cells. Both solar cell technologies suffer from incomplete photon absorption in some part of the solar spectrum. Quantum dot intermediate band solar cells are theoretically capable of achieving much higher efficiencies than conventional single-gap devices. Present prototypes suffer from extremely weak absorption of subbandgap photons in the quantum dots. This problem has received little attention so far, yet it is a serious barrier to the technology approaching its theoretical efficiency limit. Crystalline silicon solar cells absorb weakly in the near infrared due to their indirect bandgap. This problem has received much attention over recent decades, and all commercial crystalline silicon solar cells employ some form of light trapping. With the industry moving toward thinner and thinner wafers, light trapping is becoming of greater importance and diffractive structures may offer an improvement over the state-of-the-art. We begin by constructing a computational method with which to simulate solar cells equipped with diffraction grating textures. The method employs a wave-optical treatment of the diffraction grating, via rigorous coupled wave analysis, with a geometric-optical treatment of the thick solar cell bulk. These are combined using a steady-state matrix formalism. The method has been implemented computationally, and is found to be efficient and to give results in good agreement with alternative methods from other authors. The theoretical upper limit to absorption enhancement in solar cells using diffractions gratings is calculated using the matrix formalism derived in the previous task. This limit is compared to the so-called Lambertian limit for light trapping with isotropic scatterers, and to the absolute upper limit to light trapping in bulk absorbers. It is found that bi-periodic gratings (square or hexagonal geometry) are capable of offering much better light trapping than uni-periodic line gratings. The upper limit depends strongly on the grating period. For large periods, diffraction gratings are theoretically able to offer light trapping at the absolute upper limit, but only if the scattering efficiencies have a particular form, which is deemed to be beyond present design capabilities. For periods similar to the incident wavelength, diffraction gratings can offer light trapping below the absolute limit but above the Lambertian limit without placing unrealistic demands on the exact form of the scattering efficiencies. This is possible for a reasonably broad wavelength range. The computational method is used to design and optimise diffraction gratings for light trapping in solar cells. The proposed diffraction grating consists of a hexagonal lattice of cylindrical wells etched into the rear of the bulk solar cell absorber. This is encapsulated in a dielectric buffer layer, and capped with a rear reflector. Simulations are made of this grating profile applied to a crystalline silicon solar cell and to a quantum dot intermediate band solar cell. The grating period, well depth, and lateral well dimensions are optimised numerically for both solar cell types. This yields the optimum parameters to be used in fabrication of grating equipped solar cells. The optimum parameters are explained using simple physical concepts, allowing us to make more general statements that can be applied to other solar cell technologies. Diffraction grating textures are fabricated on crystalline silicon substrates using nano-imprint lithography and reactive ion etching. The optimum grating period from the previous task has been used as a design parameter. The substrates have been processed into solar cell precursors for optical measurements. Reflection spectroscopy measurements confirm that bi-periodic square gratings offer better absorption enhancement than uni-periodic line gratings. The fabricated structures have been simulated with the previously developed computation tool, with good agreement between measurement and simulation results. The simulations reveal that a significant amount of the incident photons are absorbed parasitically in the rear reflector, and that this is exacerbated by the non-planarity of the rear reflector. An alternative method of depositing the dielectric buffer layer was developed, which leaves a planar surface onto which the reflector is deposited. It was found that samples prepared in this way suffered less from parasitic reflector absorption. The next task described in the thesis is the study of photon absorption in semiconductor quantum dots. The bound-state energy levels of in InAs/GaAs quantum dots is calculated using the effective mass approximation. A one- and four- band method is applied to the calculation of electron and hole wavefunctions respectively, with an empirical Hamiltonian being employed in the latter case. The strength of optical transitions between the bound states is calculated using the Fermi golden rule. The effect of the quantum dot dimensions on the energy levels and transition strengths is investigated. It is found that a strong direct transition between the ground intermediate state and the conduction band can be promoted by decreasing the quantum dot width from its value in present prototypes. This has the added benefit of reducing the ladder of excited states between the ground state and the conduction band, which may help to reduce thermal escape of electrons from quantum dots: an undesirable phenomenon from the point of view of the open circuit voltage of an intermediate band solar cell. A realistic detailed balance model is developed for quantum dot solar cells, which uses as input the energy levels and transition strengths calculated in the previous task. The model calculates the transition currents between the many intermediate levels and the valence and conduction bands under a given set of conditions. It is distinct from previous idealised detailed balance models, which are used to calculate limiting efficiencies, since it makes realistic assumptions about photon absorption by each transition. The model is used to reproduce published experimental quantum efficiency results at different temperatures, with quite good agreement. The much-studied phenomenon of thermal escape from quantum dots is found to be photonic; it is due to thermal photons, which induce transitions between the ladder of excited states between the ground intermediate state and the conduction band. In the final chapter, the realistic detailed balance model is combined with the diffraction grating simulation method to predict the effect of incorporating a diffraction grating into a quantum dot intermediate band solar cell. Careful optimisation of the grating period is made to balance the enhancement given to the different intermediate transitions, which occur in series. Due to the extremely weak absorption in the quantum dots, it is found that light trapping alone is not sufficient to achieve high subbandgap currents in quantum dot solar cells. Instead, a combination of light trapping and increased quantum dot density is required. Within the radiative limit, a quantum dot solar cell with no light trapping requires a 1000 fold increase in the number of quantum dots to supersede the efficiency of a single-gap reference cell. A quantum dot solar cell equipped with a diffraction grating requires between a 10 and 100 fold increase in the number of quantum dots, depending on the level of parasitic absorption in the rear reflector.

Realistic detailed balance study of the quantum efficiency of quantum dot solar cells

An attractive but challenging technology for high efficiency solar energy conversion is the intermediate band solar cell (IBSC), whose theoretical efficiency limit is 63%, yet which has so far failed to yield high efficiencies in practice. The most advanced IBSC technology is that based on quantum dots (QDs): the QD-IBSC. In this paper, k·p calculations of photon absorption in the QDs are combined with a multi-level detailed balance model. The model has been used to reproduce the measured quantum efficiency of a real QD-IBSC and its temperature dependence. This allows the analysis of individual sub-bandgap transition currents, which has as yet not been possible experimentally, yielding a deeper understanding of the failure of current QD-IBSCs. Based on the agreement with experimental data, the model is believed to be realistic enough to evaluate future QD-IBSC proposals.

Modelling of landslides propagation with SPH : effects of rheology and pore water pressure

El estudio desarrollado en este trabajo de tesis se centra en la modelización numérica de la fase de propagación de los deslizamientos rápidos de ladera a través del método sin malla Smoothed Particle Hydrodynamics (SPH). Este método tiene la gran ventaja de permitir el análisis de problemas de grandes deformaciones evitando operaciones costosas de remallado como en el caso de métodos numéricos con mallas tal como el método de los Elementos Finitos. En esta tesis, particular atención viene dada al rol que la reología y la presión de poros desempeñan durante estos eventos. El modelo matemático utilizado se basa en la formulación de Biot-Zienkiewicz v - pw, que representa el comportamiento, expresado en términos de velocidad del esqueleto sólido y presiones de poros, de la mezcla de partículas sólidas en un medio saturado. Las ecuaciones que gobiernan el problema son: • la ecuación de balance de masa de la fase del fluido intersticial, • la ecuación de balance de momento de la fase del fluido intersticial y de la mezcla, • la ecuación constitutiva y • una ecuación cinemática. Debido a sus propiedades geométricas, los deslizamientos de ladera se caracterizan por tener una profundidad muy pequeña frente a su longitud y a su anchura, y, consecuentemente, el modelo matemático mencionado anteriormente se puede simplificar integrando en profundidad las ecuaciones, pasando de un modelo 3D a 2D, el cual presenta una combinación excelente de precisión, sencillez y costes computacionales. El modelo propuesto en este trabajo se diferencia de los modelos integrados en profundidad existentes por incorporar un ulterior modelo capaz de proveer información sobre la presión del fluido intersticial a cada paso computacional de la propagación del deslizamiento. En una manera muy eficaz, la evolución de los perfiles de la presión de poros está numéricamente resuelta a través de un esquema explicito de Diferencias Finitas a cada nodo SPH. Este nuevo enfoque es capaz de tener en cuenta la variación de presión de poros debida a cambios de altura, de consolidación vertical o de cambios en las tensiones totales. Con respecto al comportamiento constitutivo, uno de los problemas principales al modelizar numéricamente deslizamientos rápidos de ladera está en la dificultad de simular con la misma ley constitutiva o reológica la transición de la fase de iniciación, donde el material se comporta como un sólido, a la fase de propagación donde el material se comporta como un fluido. En este trabajo de tesis, se propone un nuevo modelo reológico basado en el modelo viscoplástico de Perzyna, pensando a la viscoplasticidad como a la llave para poder simular tanto la fase de iniciación como la de propagación con el mismo modelo constitutivo. Con el fin de validar el modelo matemático y numérico se reproducen tanto ejemplos de referencia con solución analítica como experimentos de laboratorio. Finalmente, el modelo se aplica a casos reales, con especial atención al caso del deslizamiento de 1966 en Aberfan, mostrando como los resultados obtenidos simulan con éxito estos tipos de riesgos naturales. The study developed in this thesis focuses on the modelling of landslides propagation with the Smoothed Particle Hydrodynamics (SPH) meshless method which has the great advantage of allowing to deal with large deformation problems by avoiding expensive remeshing operations as happens for mesh methods such as, for example, the Finite Element Method. In this thesis, special attention is given to the role played by rheology and pore water pressure during these natural hazards. The mathematical framework used is based on the v - pw Biot-Zienkiewicz formulation, which represents the behaviour, formulated in terms of soil skeleton velocity and pore water pressure, of the mixture of solid particles and pore water in a saturated media. The governing equations are: • the mass balance equation for the pore water phase, • the momentum balance equation for the pore water phase and the mixture, • the constitutive equation and • a kinematic equation. Landslides, due to their shape and geometrical properties, have small depths in comparison with their length or width, therefore, the mathematical model aforementioned can then be simplified by depth integrating the equations, switching from a 3D to a 2D model, which presents an excellent combination of accuracy, computational costs and simplicity. The proposed model differs from previous depth integrated models by including a sub-model able to provide information on pore water pressure profiles at each computational step of the landslide's propagation. In an effective way, the evolution of the pore water pressure profiles is numerically solved through a set of 1D Finite Differences explicit scheme at each SPH node. This new approach is able to take into account the variation of the pore water pressure due to changes of height, vertical consolidation or changes of total stress. Concerning the constitutive behaviour, one of the main issues when modelling fast landslides is the difficulty to simulate with the same constitutive or rheological model the transition from the triggering phase, where the landslide behaves like a solid, to the propagation phase, where the landslide behaves in a fluid-like manner. In this work thesis, a new rheological model is proposed, based on the Perzyna viscoplastic model, thinking of viscoplasticity as the key to close the gap between the triggering and the propagation phase. In order to validate the mathematical model and the numerical approach, benchmarks and laboratory experiments are reproduced and compared to analytical solutions when possible. Finally, applications to real cases are studied, with particular attention paid to the Aberfan flowslide of 1966, showing how the mathematical model accurately and successfully simulate these kind of natural hazards.

A Hydrodynamic Method For Measuring Aqueous Nanoparticle Surface Interactions

The objectives of this research dissertation were to develop and present novel analytical methods for the quantification of surface binding interactions between aqueous nanoparticles and water-soluble organic solutes. Quantification of nanoparticle surface interactions are presented in this work as association constants where the solutes have interacted with the surface of the nanoparticles. By understanding these nanoparticle-solute interactions, in part through association constants, the scientific community will better understand how organic drugs and nanomaterials interact in the environment, as well as to understand their eventual environmental fate. The biological community, pharmaceutical, and consumer product industries also have vested interests in nanoparticle-drug interactions for nanoparticle toxicity research and in using nanomaterials as drug delivery vesicles. The presented novel analytical methods, applied to nanoparticle surface association chemistry, may prove to be useful in assisting the scientific community to understand the risks, benefits, and opportunities of nanoparticles. The development of the analytical methods presented uses a model nanoparticle, Laponite-RD (LRD). LRD was the proposed nanoparticle used to model the system and technique because of its size, 25 nm in diameter. The solutes selected to model for these studies were chosen because they are also environmentally important. Caffeine, oxytetracycline (OTC), and quinine were selected to use as models because of their environmental importance and chemical properties that can be exploited in the system. All of these chemicals are found in the environment; thus, how they interact with nanoparticles and are transported through the environment is important. The analytical methods developed utilize and a wide-bore hydrodynamic chromatography to induce a partial hydrodynamic separation between nanoparticles and dissolved solutes. Then, using deconvolution techniques, two separate elution profiles for the nanoparticle and organic solute can be obtained. Followed by a mass balance approach, association constants between LRD, our model nanoparticle, and organic solutes are calculated. These findings are the first of their kind for LRD and nanoclays in dilute dispersions.

Estudo cinético da lixiviação de metais de placas de circuito impresso obsoletas.

O processo tradicional de recuperação de metais de resíduos de equipamentos eletroeletrônicos (REEE) geralmente envolve processamento pirometalúrgico. Entretanto, o uso desta tecnologia para processar placas de circuito impresso (PCI) obsoletas pode levar à liberação de dioxinas e furanos, devido à decomposição térmica de retardantes de chama e resinas poliméricas presentes no substrato das placas. Portanto, este trabalho propõe uma rota hidrometalúrgica para recuperação de metais. O comportamento dos metais, com destaque para cobre, zinco e níquel, durante a lixiviação ácida, foi estudado em três temperaturas diferentes (35ºC, 65ºC e 75ºC), com e sem adição de um agente oxidante (peróxido de hidrogênio H2O2). A cinética de dissolução ácida desses metais foi estudada baseada na análise química por ICP-OES (Espectrometria de emissão ótica por plasma acoplado indutivamente) e EDX (Espectroscopia de fluorescência de raios-X por energia dispersiva). O balanço de massa e a análise química indicaram que a etapa de lixiviação sem adição de oxidante é pouco eficaz na extração dos metais, sendo responsável pela dissolução de menos do que 6% do total extraído. A 65ºC e H2SO4 1 mol/L, com adição de 5 mL de H2O2 (30%) a cada quinze minutos e densidade de polpa de 1 g / 10 mL, 98,1% do cobre, 99,9% do zinco e 99,0% do níquel foram extraídos após 4 horas. A cinética de dissolução desses metais é controlada pela etapa da reação química, seguindo, dependendo da temperatura, a equação 1 (1 XB)1/3 = k1.t ou a equação ln (1 XB) = k4.t.

A Construction Management Framework for Mass Customisation in Traditional Construction

A Mass Customisation model is discussed as a competitive positioning strategy in the marketplace adding value to the customer’s end-use. It includes the user as part of the construction process responding to the customer’s demands and wishes. To the present day, almost all proposals for Mass Customisation have been focused on the design phase and single family houses. The reality is that the processes carried out in the work execution are so inefficient that the costs of the Mass Customisation models are assumed by the customer and they do not offer solutions that support the change management. Furthermore, this inefficiency often makes Mass Customisation unfeasible in terms of deadlines and site management. Therefore, the present proposal focuses on achieving the paradigm of Mass Customisation in the traditional residential construction complementary to the existing proposals in the design phase. All this through the proposal of a framework for the integral management in the work execution, which will address change management introduced by the users offering an efficient and productive model that reduces costs in the process. This model will focus on the synergy between different strategies, techniques and technologies currently used in the construction management (such as Lean Construction or Six Sigma), together with, other strategies and technologies that have proven to be valid solutions in other fields (such as Business Process Management, Service Oriented Architecture, etc.).

Solar radiation over and under sea ice during the POLARSTERN cruise ARK-XXVI/3 (TransArc) in summer 2011

Arctic sea ice has declined and become thinner and younger (more seasonal) during the last decade. One consequence of this is that the surface energy budget of the Arctic Ocean is changing. While the role of surface albedo has been studied intensively, it is still widely unknown how much light penetrates through sea ice into the upper ocean, affecting sea-ice mass balance, ecosystems, and geochemical processes. Here we present the first large-scale under-ice light measurements, operating spectral radiometers on a remotely operated vehicle (ROV) under Arctic sea ice in summer. This data set is used to produce an Arctic-wide map of light distribution under summer sea ice. Our results show that transmittance through first-year ice (FYI, 0.11) was almost three times larger than through multi-year ice (MYI, 0.04), and that this is mostly caused by the larger melt-pond coverage of FYI (42 vs. 23%). Also energy absorption was 50% larger in FYI than in MYI. Thus, a continuation of the observed sea-ice changes will increase the amount of light penetrating into the Arctic Ocean, enhancing sea-ice melt and affecting sea-ice and upper-ocean ecosystems.

Greenland and Antarctic ice sheet topography, cavity geometry, and global bathymetry (RTopo-2), links to NetCDF files

The ocean plays an important role in modulating the mass balance of the polar ice sheets by interacting with the ice shelves in Antarctica and with the marine-terminating outlet glaciers in Greenland. Given that the flux of warm water onto the continental shelf and into the sub-ice cavities is steered by complex bathymetry, a detailed topography data set is an essential ingredient for models that address ice-ocean interaction. We followed the spirit of the global RTopo-1 data set and compiled consistent maps of global ocean bathymetry, upper and lower ice surface topographies and global surface height on a spherical grid with now 30-arc seconds resolution. We used the General Bathymetric Chart of the Oceans (GEBCO, 2014) as the backbone and added the International Bathymetric Chart of the Arctic Ocean version 3 (IBCAOv3) and the Interna- tional Bathymetric Chart of the Southern Ocean (IBCSO) version 1. While RTopo-1 primarily aimed at a good and consistent representation of the Antarctic ice sheet, ice shelves and sub-ice cavities, RTopo-2 now also contains ice topographies of the Greenland ice sheet and outlet glaciers. In particular, we aimed at a good representation of the fjord and shelf bathymetry sur- rounding the Greenland continent. We corrected data from earlier gridded products in the areas of Petermann Glacier, Hagen Bræ and Sermilik Fjord assuming that sub-ice and fjord bathymetries roughly follow plausible Last Glacial Maximum ice flow patterns. For the continental shelf off northeast Greenland and the floating ice tongue of Nioghalvfjerdsfjorden Glacier at about 79°N, we incorporated a high-resolution digital bathymetry model considering original multibeam survey data for the region. Radar data for surface topographies of the floating ice tongues of Nioghalvfjerdsfjorden Glacier and Zachariæ Isstrøm have been obtained from the data centers of Technical University of Denmark (DTU), Operation Icebridge (NASA/NSF) and Alfred Wegener Institute (AWI). For the Antarctic ice sheet/ice shelves, RTopo-2 largely relies on the Bedmap-2 product but applies corrections for the geometry of Getz, Abbot and Fimbul ice shelf cavities.

Volume changes of Jacobshavn Isbrae between 1985, 1997 and 2007

Following three decades of relative stability, Jakobshavn Isbrae, West Greenland, underwent dramatic thinning, retreat and speed-up starting in 1998. To assess the amount of ice loss, we analyzed 1985 aerial photos and derived a 40 m grid digital elevation model (DEM). We also obtained a 2007 40 m grid SPOT DEM covering the same region. Comparison of the two DEMs over an area of ~4000 km**2 revealed a total ice loss of 160 ± 4 km**3, with 107 ± 0.2 km**3 in grounded regions (0.27 mm eustatic sea-level rise) and 53 ± 4 km**3 from the disintegration of the floating tongue. Comparison of the DEMs with 1997 NASA Airborne Topographic Mapper data indicates that this ice loss essentially occurred after 1997, with +0.7 ± 5.6 km**3 between 1985 and 1997 and -160 ± 7 km**3 between 1997 and 2007. The latter is equivalent to an average specific mass balance of -3.7 ± 0.2 m/a over the study area. Previously reported thickening of the main glacier during the early 1990s was accompanied by similar-magnitude thinning outside the areas of fast flow, indicating that the land-based ice continued reacting to longer-term climate forcing.

Quantifying ice-shelf basal melt and ice-stream dynamics using high-resolution DEM and GPS time series

Thesis (Ph.D.)--University of Washington, 2016-06

Chemical kinetic investigation of a commercial batch reactor process

The aim of this investigation was to study the chemical reactions occurring during the batchwise production of a butylated melamine-formaldehyde resin, in order to optimise the efficiency and economics of the batch processes. The batch process models are largely empirical in nature as the reaction mechanism is unknown. The process chemistry and the commercial manufacturing method are described. A small scale system was established in glass and the ability to produce laboratory resins with the required quality was demonstrated, simulating the full scale plant. During further experiments the chemical reactions of methylolation, condensation and butylation were studied. The important process stages were identified and studied separately. The effects of variation of certain process parameters on the chemical reactions were also studied. A published model of methylolation was modified and used to simulate the methylolation stage. A major result of this project was the development of an indirect method for studying the condensation and butylation reactions occurring during the dehydration and acid reaction stages, as direct quantitative methods were not available. A mass balance method was devised for this purpose and used to collect experimental data. The reaction scheme was verified using this data. The reactions stages were simulated using an empirical model. This has revealed new information regarding the mechanism and kinetics of the reactions. Laboratory results were shown to be comparable with plant scale results. This work has improved the understanding of the batch process, which can be used to improve product consistency. Future work has been identified and recommended to produce an optimum process and plant design to reduce the batch time.

Osmotic dehydration in plant tissues

The primary aim of the thesis is to provide a comprehensive investigation of the osmotic dehydration processes in plant tissue. Effort has been concentrated on the modelling for simulating the processes. Two mathematical models for simulating the mass transfer during osmotic dehydration processes in plant tissues are developed and verified using existing experimental data. Both models are based on the mechanism of diffusion and convection of any mobile material that can transport in plant tissues. The mass balance equation for the transport of each constituent is established separately for intracellular and extra-cellular volumes with taking into account the mass transfer across the cell membrane the intracellular and extra-cellular volumes and the shrinkage of the whole tissue. The contribution from turgor pressure is considered in both models. Model two uses Darcy’s law to build the relation between shrinkage velocity and hydrostatic pressure in each volume because the plant tissue can be considered as the porous medium. Moreover, it has been extended to solve the multi-dimensional problems. A lot of efforts have been made to the parameter study and the sensitivity analyses. The parameters investigated including the concentration of the osmotic solution, diffusion coefficient, permeability of the cell membrane, elastic modulus of the cell wall, critical cell volume etc. The models allow us to quantitatively simulate the time evolution of intracellular and extra-cellular volumes as well as the time evolution of concentrations in each cross-section.

Concentration polarization model of spiral-wound membrane modules with application to batch-mode RO desalination of brackish water

Batch-mode reverse osmosis (batch-RO) operation is considered a promising desalination method due to its low energy requirement compared to other RO system arrangements. To improve and predict batch-RO performance, studies on concentration polarization (CP) are carried out. The Kimura-Sourirajan mass-transfer model is applied and validated by experimentation with two different spiral-wound RO elements. Explicit analytical Sherwood correlations are derived based on experimental results. For batch-RO operation, a new genetic algorithm method is developed to estimate the Sherwood correlation parameters, taking into account the effects of variation in operating parameters. Analytical procedures are presented, then the mass transfer coefficient models are developed for different operation processes, i.e., batch-RO and continuous RO. The CP related energy loss in batch-RO operation is quantified based on the resulting relationship between feed flow rates and mass transfer coefficients. It is found that CP increases energy consumption in batch-RO by about 25% compared to the ideal case in which CP is absent. For continuous RO process, the derived Sherwood correlation predicted CP accurately. In addition, we determined the optimum feed flow rate of our batch-RO system.

Modeling of the spectrum in a random distributed feedback fiber laser within the power balance modes

The simplest model for a description of the random distributed feedback (RDFB) Raman fiber laser is a power balance model describing the evolution of the intensities of the waves over the fiber length. The model predicts well the power performances of the RDFB fiber laser including the generation threshold, the output power and pump and generation wave intensity distributions along the fiber. In the present work, we extend the power balance model and modify equations in such a way that they describe now frequency dependent spectral power density instead of integral over the spectrum intensities. We calculate the generation spectrum by using the depleted pump wave longitudinal distribution derived from the conventional power balance model. We found the spectral balance model to be sufficient to account for the spectral narrowing in the RDFB laser above the threshold of the generation. © 2014 SPIE.

Biomass fast pyrolysis energy balance of a 1kg/h test rig

The present paper offers a methodological approach towards the estimation and definition of enthalpies constituting an energy balance around a fast pyrolysis experiment conducted in a laboratory scale fluid bed with a capacity of 1 kg/ h. Pure N2 was used as fluidization medium at atmospheric pressure and the operating temperature (∼500°C) was adjusted with electrical resistors. The biomass feedstock type that was used was beech wood. An effort was made to achieve a satisfying 92.5% retrieval of products (dry basis mass balance) with the differences mainly attributed to loss of some bio-oil constituents into the quenching medium, ISOPAR™. The chemical enthalpy recovery for bio-oil, char and permanent gases is calculated 64.6%, 14.5% and 7.1%, respectively. All the energy losses from the experimental unit into the environment, namely the pyrolyser, cooling unit etc. are discussed and compared to the heat of fast pyrolysis that was calculated at 1123.5 kJ per kg of beech wood. This only represents 2.4% of the biomass total enthalpy or 6.5% its HHV basis. For the estimation of some important thermo-physical properties such as heat capacity and density, it was found that using data based on the identified compounds from the GC/MS analysis is very close to the reference values despite the small fraction of the bio-oil components detected. The methodology and results can help as a starting point for the proper design of fast pyrolysis experiments, pilot and/or industrial scale plants.