Manufacturing of uniformly-sized silicon particles for solar cell from molten metal jet by electromagnetic pinch force

Spherical silicon solar cells are expected to serve as a technology to reduce silicon usage of photovoltaic (PV) power systems[1, 2, 3]. In order to establish the spherical silicon solar cell, a manufacturing method of uniformly sized silicon particles of 1mm in diameter is required. However, it is difficult to mass-produce the mono-sized silicon particles at low cost by existent processes now. We proposed a new method to generate liquid metal droplets uniformly by applying electromagnetic pinch force to a liquid metal jet[4]. The electromagnetic force was intermittently applied to the liquid metal jet issued from a nozzle in order to fluctuate the surface of the jet. As the fluctuation grew, the liquid jet was broken up into small droplets according to a frequency of the intermittent electromagnetic force. Firstly, a preliminary experiment was carried out. A single pulse current was applied instantaneously to a single turn coil around a molten gallium jet. It was confirmed that the jet could be split up by pinch force generated by the current. And then, electromagnetic pinch force was applied intermittently to the jet. It was found that the jet was broken up into mono-sized droplets in the case of a force frequency was equal to a critical frequency[5], which corresponds to a natural disturbance wave length of the jet. Numerical simulations of the droplet generation from the liquid jet were then carried out, which consisted of an electromagnetic analysis and a fluid flow calculation with a free surface of the jet. The simulation results were compared with the experiments and the agreement between the two was quite good.

Fundamentals of steady state, non-Newtonian rimming flow

Rimming flow on the inner surface of a horizontal rotating cylinder is investigated. Using a scale analysis, a theoretical description is obtained for steady-state non-Newtonian flow. Simple lubrication theory is applied since the Reynolds number is small and the liquid film is thin. Since the Deborah number is very small the flow is viscometric. The shear-thinning number, which characterizes the shear-thinning effect, may be small or large. A general constitutive law for this kind of flow requires only a single function relating shear stress and shear rate that corresponds to a generalized Newtonian liquid. For this case the run-off condition for rimming flow is derived. Provided the run-off condition is satisfied, the existence of a continuous steady-state solution is proved. The rheological models, which show Newtonian behavior at low shear rates with transition to power-law shear thinning at moderate shear rates, are considered. Numerical results are carried out for the Carreau and Ellis models, which exhibit Newtonian behavior near the free surface and power-law behavior near the wall of the rotating cylinder.

Biaxial Characterisation of Materials for Thermoforming and Blow Moulding

During free surface moulding processes such as thermoforming and blow moulding heated polymer materials are subjected to rapid biaxial deformation as they are drawn into the shape of a mould. In the development of process simulations it is therefore essential to be able to accurately measure and model this behaviour. Conventional uniaxial test methods are generally inadequate for this purpose and this has led to the development of specialised biaxial test rigs. In this paper the results of several programmes of biaxial tests conducted at Queen’s University are presented and discussed. These have included tests on high impact polystyrene (HIPS), polypropylene (PP) and aPET, and the work has involved a wide variety of experimental conditions. In all cases the results clearly demonstrate the unique characteristics of materials when subjected to biaxial deformation. PP draws the highest stresses and it is the most temperature sensitive of the materials. aPET is initially easier to form but exhibits strain hardening at higher strains. This behaviour is increased with increasing strain rate but at very high strain rates these effects are increasingly mollified by adiabatic heating. Both aPET and PP (to a lesser degree) draw much higher stresses in sequential stretching showing that this behaviour must be considered in process simulations. HIPS showed none of these effects and it is the easiest material to deform.

The effects of lichen cover upon the rate of solutional weathering of limestone

The contribution of lichens to the biomodification of limestone surfaces is an area of conflict within bioweathering studies, with some researchers suggesting a protective effect induced by lichen coverage and others a deteriorative effect induced by the same organisms.Data are reported demonstrating the potential role of endolithic lichen, in particular of Bagliettoa baldensis, in the active protection of Carboniferous limestone surfaces from rainfall-induced solutional weathering. During a 12-month microcatchment exposure period in the west of Northern Ireland, average dissolutional losses of calciumare greater from a lichen-free limestone surface compared with a predominantly endolithic lichen-covered surface by just under 1.25 times. During colderwintermonths, the lichen free surface experiences calcium loss almost 1.5 times greater than the lichen-covered surface. Using extrapolation to upscale from the micro-catchment sample scale, for the year of sample exposure, the rate of calcium loss is 1.001 g m−2 a−1 from lichen-covered limestone surfaces and 1.228 gm−2 a−1 from lichen-free bare limestone surfaces. This research has implications for our understanding of karst environments, the contribution of lichens to karren development and the conservation of lichen-colonised dimension stone within a cultural setting.

Predicting axial velocity profiles within a diffusing marine propeller jet

A full understanding of the hydrodynamic processes within the jet produced by a manoeuvring ship’s propeller is essential in the development and maintenance of ports, docks and harbours. In this study the predominant axial velocity component of a diffusing propeller jet was studied. The flow fields formed by four propellers, each operating at four power levels (speeds of rotation), were investigated under bollard pull conditions within a large free surface tank using Laser Doppler Anemometry. Comparison were made to existing methodologies by which a prediction of the magnitudes of the axial velocity can be made, and where deficient modifications to the methodologies have been developed. The jets were found to produce a maximum axial velocity along the initial efflux plane at a location near the blade mid-span. The position and magnitude of the axial velocity was seen to decrease as the jet entrained more flow and transitioned from the zone of flow establishment into the zone of established flow.

Development of numerical methodologies for parameter identification and shape optimization in metal forming simulations

Por parte da indústria de estampagem tem-se verificado um interesse crescente em simulações numéricas de processos de conformação de chapa, incluindo também métodos de engenharia inversa. Este facto ocorre principalmente porque as técnicas de tentativa-erro, muito usadas no passado, não são mais competitivas a nível económico. O uso de códigos de simulação é, atualmente, uma prática corrente em ambiente industrial, pois os resultados tipicamente obtidos através de códigos com base no Método dos Elementos Finitos (MEF) são bem aceites pelas comunidades industriais e científicas Na tentativa de obter campos de tensão e de deformação precisos, uma análise eficiente com o MEF necessita de dados de entrada corretos, como geometrias, malhas, leis de comportamento não-lineares, carregamentos, leis de atrito, etc.. Com o objetivo de ultrapassar estas dificuldades podem ser considerados os problemas inversos. No trabalho apresentado, os seguintes problemas inversos, em Mecânica computacional, são apresentados e analisados: (i) problemas de identificação de parâmetros, que se referem à determinação de parâmetros de entrada que serão posteriormente usados em modelos constitutivos nas simulações numéricas e (ii) problemas de definição geométrica inicial de chapas e ferramentas, nos quais o objetivo é determinar a forma inicial de uma chapa ou de uma ferramenta tendo em vista a obtenção de uma determinada geometria após um processo de conformação. São introduzidas e implementadas novas estratégias de otimização, as quais conduzem a parâmetros de modelos constitutivos mais precisos. O objetivo destas estratégias é tirar vantagem das potencialidades de cada algoritmo e melhorar a eficiência geral dos métodos clássicos de otimização, os quais são baseados em processos de apenas um estágio. Algoritmos determinísticos, algoritmos inspirados em processos evolucionários ou mesmo a combinação destes dois são usados nas estratégias propostas. Estratégias de cascata, paralelas e híbridas são apresentadas em detalhe, sendo que as estratégias híbridas consistem na combinação de estratégias em cascata e paralelas. São apresentados e analisados dois métodos distintos para a avaliação da função objetivo em processos de identificação de parâmetros. Os métodos considerados são uma análise com um ponto único ou uma análise com elementos finitos. A avaliação com base num único ponto caracteriza uma quantidade infinitesimal de material sujeito a uma determinada história de deformação. Por outro lado, na análise através de elementos finitos, o modelo constitutivo é implementado e considerado para cada ponto de integração. Problemas inversos são apresentados e descritos, como por exemplo, a definição geométrica de chapas e ferramentas. Considerando o caso da otimização da forma inicial de uma chapa metálica a definição da forma inicial de uma chapa para a conformação de um elemento de cárter é considerado como problema em estudo. Ainda neste âmbito, um estudo sobre a influência da definição geométrica inicial da chapa no processo de otimização é efetuado. Este estudo é realizado considerando a formulação de NURBS na definição da face superior da chapa metálica, face cuja geometria será alterada durante o processo de conformação plástica. No caso dos processos de otimização de ferramentas, um processo de forjamento a dois estágios é apresentado. Com o objetivo de obter um cilindro perfeito após o forjamento, dois métodos distintos são considerados. No primeiro, a forma inicial do cilindro é otimizada e no outro a forma da ferramenta do primeiro estágio de conformação é otimizada. Para parametrizar a superfície livre do cilindro são utilizados diferentes métodos. Para a definição da ferramenta são também utilizados diferentes parametrizações. As estratégias de otimização propostas neste trabalho resolvem eficientemente problemas de otimização para a indústria de conformação metálica.

Estimating depth and explosive charge weight for an extremely shallow underwater explosion of the ROKS Cheonan sinking in the Yellow Sea

We estimated the detonation depth and net explosive weight for a very shallow underwater explosion using cutoff frequencies and spectral analysis. With detonation depth and a bubble pulse the net explosive weight for a shallow underwater explosion could simply be determined. The ray trace modeling confirms the detonation depth as a source of the hydroacoustic wave propagation in a shallow channel. We found cutoff frequencies of the reflection off the ocean bottom to be 8.5 Hz, 25 Hz, and 43 Hz while the cutoff frequency of the reflection off the free surface to be 45 Hz including 1.01 Hz for the bubble pulse, and also found the cutoff frequency of surface reflection to well fit the ray-trace modeling. We also attempted to corroborate our findings using a 3D bubble shape modeling and boundary element method. Our findings led us to the net explosive weight of the underwater explosion offshore of Baengnyeong-do for the ROKS Cheonan sinking to be approximately 136 kg TNT at a depth of about 8 m within an ocean depth of around 44 m. © 2015 Elsevier B.V.

Development of new technologies for the selective removal of volatile compounds with applications in forensic science and natural aromas

Tese de doutoramento, Farmácia (Química Farmacêutica e Terapêutica), Universidade de Lisboa, Faculdade de Farmácia, 2015

Propriétés des valeurs propres de ballotement pour contenants symétriques

Le problème d’oscillation de fluides dans un conteneur est un problème classique d’hydrodynamique qui est etudié par des mathématiciens et ingénieurs depuis plus de 150 ans. Le présent travail est lié à l’étude de l’alternance des fonctions propres paires et impaires du problème de Steklov-Neumann pour les domaines à deux dimensions ayant une forme symétrique. On obtient des résultats sur la parité de deuxième et troisième fonctions propres d’un tel problème pour les trois premiers modes, dans le cas de domaines symétriques arbitraires. On étudie aussi la simplicité de deux premières valeurs propres non nulles d’un tel problème. Il existe nombre d’hypothèses voulant que pour le cas des domaines symétriques, toutes les valeurs propres sont simples. Il y a des résultats de Kozlov, Kuznetsov et Motygin [1] sur la simplicité de la première valeur propre non nulle obtenue pour les domaines satisfaisants la condition de John. Dans ce travail, il est montré que pour les domaines symétriques, la deuxième valeur propre non-nulle du problème de Steklov-Neumann est aussi simple.

Modélisation de l'évolution du réseau magnétique au cours du cycle solaire

Le réseau magnétique consiste en un ensemble de petites concentrations de flux magnétique sur la photosphère solaire. Vu sa petite échelle de taille et de flux, à la limite de détection, son comportement n'est connu que depuis récemment. Les interactions du réseau sont pourtant cruciales afin de comprendre la dynamo et l'irradiance solaires, car beaucoup de caractérisques du réseau dépendent de ces interactions. De plus, le réseau est la principale contribution magnétique surfacique à l'irradiance solaire. Les modèles existants du réseau ne tenaient jusqu'à maintenant pas compte des interactions du réseau. Nous avons tenté de combler cette lacune avec notre modèle. Nos simulations impliquent une marche aléatoire en 2D de tubes de flux magnétiques sur la photosphère solaire. Les tubes de flux sont injectés puis soumis à des règles de déplacement et d'interaction. L'injection se fait à deux échelles, respectivement la plus petite et la plus grande observables: les tubes de flux élémentaires et les taches solaires. Des processus de surface imitant ceux observés sont inclus, et consistent en l'émergence, la coalescence, l'annulation et la submergence de flux. La fragmentation des concentrations n'est présente que pour les taches, sous forme de désintégration libérant des tubes de flux. Le modèle est appliqué au cycle solaire 21 (1976-1986, le mieux documenté en termes de caractéristiques de taches solaires. Il en résulte des réponses à deux questions importantes en physique solaire. La première est: l'injection de flux magnétique à deux échelles très distinctes peut-elle conduire à une distribution de flux en loi de puissance comme on l'observe, si l'on inclut des processus de surface qui retraitent le flux? Cette question est étroitement liée à l'origine de la dynamo solaire, qui pourrait produire ladite distribution. Nous trouvons qu'on peut effectivement produire une telle distribution avec ce type d'injection et ce type de processus de surface. Cela implique que la distribution de flux observée ne peut servir à déterminer quel type de dynamo opère dans le Soleil. La deuxième question à laquelle nous avons apporté un élément de réponse est celle à savoir combien de temps il faut au réseau pour retrouver son état d'activité de base. Cet état a été observé lors du minimum de Maunder en 1645-1715 et touche de près la question de l'influence de l'activité solaire sur le climat terrestre. Le récent minimum d'activité est considéré par certains comme ayant atteint cet état. Nous trouvons plutôt que ça n'a pas été le cas. En effet, le temps de relaxation du réseau que nous avons calculé est supérieur au temps écoulé entre la fin du dernier cycle solaire d'activité et celui de l'amorce du présent cycle.

Modeling the interplay between land use, bioenergy production and climate change

Land use has become a force of global importance, considering that 34% of the Earth’s ice-free surface was covered by croplands or pastures in 2000. The expected increase in global human population together with eminent climate change and associated search for energy sources other than fossil fuels can, through land-use and land-cover changes (LUCC), increase the pressure on nature’s resources, further degrade ecosystem services, and disrupt other planetary systems of key importance to humanity. This thesis presents four modeling studies on the interplay between LUCC, increased production of biofuels and climate change in four selected world regions. In the first study case two new crop types (sugarcane and jatropha) are parameterized in the LPJ for managed Lands dynamic global vegetation model for calculation of their potential productivity. Country-wide spatial variation in the yields of sugarcane and jatropha incurs into substantially different land requirements to meet the biofuel production targets for 2015 in Brazil and India, depending on the location of plantations. Particularly the average land requirements for jatropha in India are considerably higher than previously estimated. These findings indicate that crop zoning is important to avoid excessive LUCC. In the second study case the LandSHIFT model of land-use and land-cover changes is combined with life cycle assessments to investigate the occurrence and extent of biofuel-driven indirect land-use changes (ILUC) in Brazil by 2020. The results show that Brazilian biofuels can indeed cause considerable ILUC, especially by pushing the rangeland frontier into the Amazonian forests. The carbon debt caused by such ILUC would result in no carbon savings (from using plant-based ethanol and biodiesel instead of fossil fuels) before 44 years for sugarcane ethanol and 246 years for soybean biodiesel. The intensification of livestock grazing could avoid such ILUC. We argue that such an intensification of livestock should be supported by the Brazilian biofuel sector, based on the sector’s own interest in minimizing carbon emissions. In the third study there is the development of a new method for crop allocation in LandSHIFT, as influenced by the occurrence and capacity of specific infrastructure units. The method is exemplarily applied in a first assessment of the potential availability of land for biogas production in Germany. The results indicate that Germany has enough land to fulfill virtually all (90 to 98%) its current biogas plant capacity with only cultivated feedstocks. Biogas plants located in South and Southwestern (North and Northeastern) Germany might face more (less) difficulties to fulfill their capacities with cultivated feedstocks, considering that feedstock transport distance to plants is a crucial issue for biogas production. In the fourth study an adapted version of LandSHIFT is used to assess the impacts of contrasting scenarios of climate change and conservation targets on land use in the Brazilian Amazon. Model results show that severe climate change in some regions by 2050 can shift the deforestation frontier to areas that would experience low levels of human intervention under mild climate change (such as the western Amazon forests or parts of the Cerrado savannas). Halting deforestation of the Amazon and of the Brazilian Cerrado would require either a reduction in the production of meat or an intensification of livestock grazing in the region. Such findings point out the need for an integrated/multicisciplinary plan for adaptation to climate change in the Amazon. The overall conclusions of this thesis are that (i) biofuels must be analyzed and planned carefully in order to effectively reduce carbon emissions; (ii) climate change can have considerable impacts on the location and extent of LUCC; and (iii) intensification of grazing livestock represents a promising venue for minimizing the impacts of future land-use and land-cover changes in Brazil.

Water wave scattering by finite arrays of circular structures

The scattering of small amplitude water waves by a finite array of locally axisymmetric structures is considered. Regions of varying quiescent depth are included and their axisymmetric nature, together with a mild-slope approximation, permits an adaptation of well-known interaction theory which ultimately reduces the problem to a simple numerical calculation. Numerical results are given and effects due to regions of varying depth on wave loading and free-surface elevation are presented.

Benchmarking 2D hydraulic models for urban flooding

This paper describes benchmark testing of six two-dimensional (2D) hydraulic models (DIVAST, DIVASTTVD, TUFLOW, JFLOW, TRENT and LISFLOOD-FP) in terms of their ability to simulate surface flows in a densely urbanised area. The models are applied to a 1·0 km × 0·4 km urban catchment within the city of Glasgow, Scotland, UK, and are used to simulate a flood event that occurred at this site on 30 July 2002. An identical numerical grid describing the underlying topography is constructed for each model, using a combination of airborne laser altimetry (LiDAR) fused with digital map data, and used to run a benchmark simulation. Two numerical experiments were then conducted to test the response of each model to topographic error and uncertainty over friction parameterisation. While all the models tested produce plausible results, subtle differences between particular groups of codes give considerable insight into both the practice and science of urban hydraulic modelling. In particular, the results show that the terrain data available from modern LiDAR systems are sufficiently accurate and resolved for simulating urban flows, but such data need to be fused with digital map data of building topology and land use to gain maximum benefit from the information contained therein. When such terrain data are available, uncertainty in friction parameters becomes a more dominant factor than topographic error for typical problems. The simulations also show that flows in urban environments are characterised by numerous transitions to supercritical flow and numerical shocks. However, the effects of these are localised and they do not appear to affect overall wave propagation. In contrast, inertia terms are shown to be important in this particular case, but the specific characteristics of the test site may mean that this does not hold more generally.

Variation in dynamic elastic shear modulus of sandstone upon fluid saturation and substitution

Experimental acoustic measurements on sandstone rocks at both sonic and ultrasonic frequencies show that fluid saturation can cause a noticeable change in both the dynamic bulk and shear elastic moduli of sandstones. We observed that the change in dynamic shear modulus upon fluid saturation is highly dependent on the type of saturant, its viscosity, rock microstructure, and applied pressures. Frequency dispersion has some influence on dynamic elastic moduli too, but its effect is limited to the ultrasonic frequency ranges and above. We propose that viscous coupling, reduction in free surface energy, and, to a limited extent, frequency dispersion due to both local and global flow are the main mechanisms responsible for the change in dynamic shear elastic modulus upon fluid saturation and substitution, and we quantify influences.

Approximations to the scattering of water waves by steep topography

A new method is developed for approximating the scattering of linear surface gravity waves on water of varying quiescent depth in two dimensions. A conformal mapping of the fluid domain onto a uniform rectangular strip transforms steep and discontinuous bed profiles into relatively slowly varying, smooth functions in the transformed free-surface condition. By analogy with the mild-slope approach used extensively in unmapped domains, an approximate solution of the transformed problem is sought in the form of a modulated propagating wave which is determined by solving a second-order ordinary differential equation. This can be achieved numerically, but an analytic solution in the form of a rapidly convergent infinite series is also derived and provides simple explicit formulae for the scattered wave amplitudes. Small-amplitude and slow variations in the bedform that are excluded from the mapping procedure are incorporated in the approximation by a straightforward extension of the theory. The error incurred in using the method is established by means of a rigorous numerical investigation and it is found that remarkably accurate estimates of the scattered wave amplitudes are given for a wide range of bedforms and frequencies.