Numerical CFD modelling of non-neutral atmospheric boundary layers for offshore wind resource assessment based on Monin-Obukhov theory

The presented works aim at proposing a methodology for the simulation of offshore wind conditions using CFD. The main objective is the development of a numerical model for the characterization of atmospheric boundary layers of different stability levels, as the most important issue in offshore wind resource assessment. Based on Monin-Obukhov theory, the steady k-ε Standard turbulence model is modified to take into account thermal stratification in the surface layer. The validity of Monin-Obukhov theory in offshore conditions is discussed with an analysis of a three day episode at FINO-1 platform.

Modelling wake effects using two CFD techniques

Comparativa de dos técnica de simulación de rotor: BEM (Blade Element Momentum) y MRF (Moving Reference Frame) para el perfil S809 de NREL

CFD modelling of wind farms in complex terrain

Modelling of entire wind farms in flat and complex terrain using a full 3D Navier–Stokes solver for incompressible flow is presented in this paper. Numerical integration of the governing equations is performed using an implicit pressure correction scheme, where the wind turbines (W/Ts) are modelled as momentum absorbers through their thrust coefficient. The k–ω turbulence model, suitably modified for atmospheric flows, is employed for closure. A correction is introduced to account for the underestimation of the near wake deficit, in which the turbulence time scale is bounded using a general “realizability” constraint for the fluctuating velocities. The second modelling issue that is discussed in this paper is related to the determination of the reference wind speed for the thrust calculation of the machines. Dealing with large wind farms and wind farms in complex terrain, determining the reference wind speed is not obvious when a W/T operates in the wake of another WT and/or in complex terrain. Two alternatives are compared: using the wind speed value at hub height one diameter upstream of the W/T and adopting an induction factor-based concept to overcome the utilization of a wind speed at a certain distance upwind of the rotor. Application is made in two wind farms, a five-machine one located in flat terrain and a 43-machine one located in complex terrain.

CFD modelling of the interaction between the Surface Boundary Layer and rotor wake. Comparison of results obtained with different turbulence models and mesh strategies

A simplified CFD wake model based on the actuator disk concept is used to simulate the wind turbine, represented by a disk upon which a distribution of forces, defined as axial momentum sources, are applied on the incoming non-uniform flow. The rotor is supposed to be uniformly loaded, with the exerted forces function of the incident wind speed, the thrust coefficient and the rotor diameter. The model is tested under different parameterizations of turbulence models and validated through experimental measurements downwind of a wind turbine in terms of wind speed deficit and turbulence intensity.

Analysis of cogeneration in the present energy framework

In this paper, a general vision of cogeneration penetration in the European Union is shown; after this, a case study is included, evaluating as a function of two factors (electricity and emission allowance prices) the suitability of installing, for an industry with a determined thermal demand, two different options. The first one is a gas turbine cogeneration plant generating steam through a heat recovery steam generator (HRSG). The second one consists of installing a natural gas boiler for steam production covering the electricity demand from the grid. The CO2 emissions from both options are compared regarding different kinds of generation mixes from the electricity grid in the case of using the industrial boiler; taking into account the advantages of using biomass in relation to emissions, a last comparison has been carried out considering a biomass boiler instead of the natural gas boiler.

Application of gas-cooled Accelerator Driven System (ADS) transmutation devices to sustainable nuclear energy development

The conceptual design of a pebble bed gas-cooled transmutation device is shown with the aim to evaluate its potential for its deployment in the context of the sustainable nuclear energy development, which considers high temperature reactors for their operation in cogeneration mode, producing electricity, heat and Hydrogen. As differential characteristics our device operates in subcritical mode, driven by a neutron source activated by an accelerator that adds clear safety advantages and fuel flexibility opening the possibility to reduce the nuclear stockpile producing energy from actual LWR irradiated fuel with an efficiency of 45?46%, either in the form of Hydrogen, electricity, or both.

Numerical Error Prediction and its applications in CFD using tau-estimation

Nowadays, Computational Fluid Dynamics (CFD) solvers are widely used within the industry to model fluid flow phenomenons. Several fluid flow model equations have been employed in the last decades to simulate and predict forces acting, for example, on different aircraft configurations. Computational time and accuracy are strongly dependent on the fluid flow model equation and the spatial dimension of the problem considered. While simple models based on perfect flows, like panel methods or potential flow models can be very fast to solve, they usually suffer from a poor accuracy in order to simulate real flows (transonic, viscous). On the other hand, more complex models such as the full Navier- Stokes equations provide high fidelity predictions but at a much higher computational cost. Thus, a good compromise between accuracy and computational time has to be fixed for engineering applications. A discretisation technique widely used within the industry is the so-called Finite Volume approach on unstructured meshes. This technique spatially discretises the flow motion equations onto a set of elements which form a mesh, a discrete representation of the continuous domain. Using this approach, for a given flow model equation, the accuracy and computational time mainly depend on the distribution of nodes forming the mesh. Therefore, a good compromise between accuracy and computational time might be obtained by carefully defining the mesh. However, defining an optimal mesh for complex flows and geometries requires a very high level expertize in fluid mechanics and numerical analysis, and in most cases a simple guess of regions of the computational domain which might affect the most the accuracy is impossible. Thus, it is desirable to have an automatized remeshing tool, which is more flexible with unstructured meshes than its structured counterpart. However, adaptive methods currently in use still have an opened question: how to efficiently drive the adaptation ? Pioneering sensors based on flow features generally suffer from a lack of reliability, so in the last decade more effort has been made in developing numerical error-based sensors, like for instance the adjoint-based adaptation sensors. While very efficient at adapting meshes for a given functional output, the latter method is very expensive as it requires to solve a dual set of equations and computes the sensor on an embedded mesh. Therefore, it would be desirable to develop a more affordable numerical error estimation method. The current work aims at estimating the truncation error, which arises when discretising a partial differential equation. These are the higher order terms neglected in the construction of the numerical scheme. The truncation error provides very useful information as it is strongly related to the flow model equation and its discretisation. On one hand, it is a very reliable measure of the quality of the mesh, therefore very useful in order to drive a mesh adaptation procedure. On the other hand, it is strongly linked to the flow model equation, so that a careful estimation actually gives information on how well a given equation is solved, which may be useful in the context of _ -extrapolation or zonal modelling. The following work is organized as follows: Chap. 1 contains a short review of mesh adaptation techniques as well as numerical error prediction. In the first section, Sec. 1.1, the basic refinement strategies are reviewed and the main contribution to structured and unstructured mesh adaptation are presented. Sec. 1.2 introduces the definitions of errors encountered when solving Computational Fluid Dynamics problems and reviews the most common approaches to predict them. Chap. 2 is devoted to the mathematical formulation of truncation error estimation in the context of finite volume methodology, as well as a complete verification procedure. Several features are studied, such as the influence of grid non-uniformities, non-linearity, boundary conditions and non-converged numerical solutions. This verification part has been submitted and accepted for publication in the Journal of Computational Physics. Chap. 3 presents a mesh adaptation algorithm based on truncation error estimates and compares the results to a feature-based and an adjoint-based sensor (in collaboration with Jorge Ponsín, INTA). Two- and three-dimensional cases relevant for validation in the aeronautical industry are considered. This part has been submitted and accepted in the AIAA Journal. An extension to Reynolds Averaged Navier- Stokes equations is also included, where _ -estimation-based mesh adaptation and _ -extrapolation are applied to viscous wing profiles. The latter has been submitted in the Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering. Keywords: mesh adaptation, numerical error prediction, finite volume Hoy en día, la Dinámica de Fluidos Computacional (CFD) es ampliamente utilizada dentro de la industria para obtener información sobre fenómenos fluidos. La Dinámica de Fluidos Computacional considera distintas modelizaciones de las ecuaciones fluidas (Potencial, Euler, Navier-Stokes, etc) para simular y predecir las fuerzas que actúan, por ejemplo, sobre una configuración de aeronave. El tiempo de cálculo y la precisión en la solución depende en gran medida de los modelos utilizados, así como de la dimensión espacial del problema considerado. Mientras que modelos simples basados en flujos perfectos, como modelos de flujos potenciales, se pueden resolver rápidamente, por lo general aducen de una baja precisión a la hora de simular flujos reales (viscosos, transónicos, etc). Por otro lado, modelos más complejos tales como el conjunto de ecuaciones de Navier-Stokes proporcionan predicciones de alta fidelidad, a expensas de un coste computacional mucho más elevado. Por lo tanto, en términos de aplicaciones de ingeniería se debe fijar un buen compromiso entre precisión y tiempo de cálculo. Una técnica de discretización ampliamente utilizada en la industria es el método de los Volúmenes Finitos en mallas no estructuradas. Esta técnica discretiza espacialmente las ecuaciones del movimiento del flujo sobre un conjunto de elementos que forman una malla, una representación discreta del dominio continuo. Utilizando este enfoque, para una ecuación de flujo dado, la precisión y el tiempo computacional dependen principalmente de la distribución de los nodos que forman la malla. Por consiguiente, un buen compromiso entre precisión y tiempo de cálculo se podría obtener definiendo cuidadosamente la malla, concentrando sus elementos en aquellas zonas donde sea estrictamente necesario. Sin embargo, la definición de una malla óptima para corrientes y geometrías complejas requiere un nivel muy alto de experiencia en la mecánica de fluidos y el análisis numérico, así como un conocimiento previo de la solución. Aspecto que en la mayoría de los casos no está disponible. Por tanto, es deseable tener una herramienta que permita adaptar los elementos de malla de forma automática, acorde a la solución fluida (remallado). Esta herramienta es generalmente más flexible en mallas no estructuradas que con su homóloga estructurada. No obstante, los métodos de adaptación actualmente en uso todavía dejan una pregunta abierta: cómo conducir de manera eficiente la adaptación. Sensores pioneros basados en las características del flujo en general, adolecen de una falta de fiabilidad, por lo que en la última década se han realizado grandes esfuerzos en el desarrollo numérico de sensores basados en el error, como por ejemplo los sensores basados en el adjunto. A pesar de ser muy eficientes en la adaptación de mallas para un determinado funcional, este último método resulta muy costoso, pues requiere resolver un doble conjunto de ecuaciones: la solución y su adjunta. Por tanto, es deseable desarrollar un método numérico de estimación de error más asequible. El presente trabajo tiene como objetivo estimar el error local de truncación, que aparece cuando se discretiza una ecuación en derivadas parciales. Estos son los términos de orden superior olvidados en la construcción del esquema numérico. El error de truncación proporciona una información muy útil sobre la solución: es una medida muy fiable de la calidad de la malla, obteniendo información que permite llevar a cabo un procedimiento de adaptación de malla. Está fuertemente relacionado al modelo matemático fluido, de modo que una estimación precisa garantiza la idoneidad de dicho modelo en un campo fluido, lo que puede ser útil en el contexto de modelado zonal. Por último, permite mejorar la precisión de la solución resolviendo un nuevo sistema donde el error local actúa como término fuente (_ -extrapolación). El presenta trabajo se organiza de la siguiente manera: Cap. 1 contiene una breve reseña de las técnicas de adaptación de malla, así como de los métodos de predicción de los errores numéricos. En la primera sección, Sec. 1.1, se examinan las estrategias básicas de refinamiento y se presenta la principal contribución a la adaptación de malla estructurada y no estructurada. Sec 1.2 introduce las definiciones de los errores encontrados en la resolución de problemas de Dinámica Computacional de Fluidos y se examinan los enfoques más comunes para predecirlos. Cap. 2 está dedicado a la formulación matemática de la estimación del error de truncación en el contexto de la metodología de Volúmenes Finitos, así como a un procedimiento de verificación completo. Se estudian varias características que influyen en su estimación: la influencia de la falta de uniformidad de la malla, el efecto de las no linealidades del modelo matemático, diferentes condiciones de contorno y soluciones numéricas no convergidas. Esta parte de verificación ha sido presentada y aceptada para su publicación en el Journal of Computational Physics. Cap. 3 presenta un algoritmo de adaptación de malla basado en la estimación del error de truncación y compara los resultados con sensores de featured-based y adjointbased (en colaboración con Jorge Ponsín del INTA). Se consideran casos en dos y tres dimensiones, relevantes para la validación en la industria aeronáutica. Este trabajo ha sido presentado y aceptado en el AIAA Journal. También se incluye una extensión de estos métodos a las ecuaciones RANS (Reynolds Average Navier- Stokes), en donde adaptación de malla basada en _ y _ -extrapolación son aplicados a perfiles con viscosidad de alas. Este último trabajo se ha presentado en los Actas de la Institución de Ingenieros Mecánicos, Parte G: Journal of Aerospace Engineering. Palabras clave: adaptación de malla, predicción del error numérico, volúmenes finitos

Modelling wake effects in large wind farms in complex terrain: the problem, the methods and the issues

Computational fluid dynamic (CFD) methods are used in this paper to predict the power production from entire wind farms in complex terrain and to shed some light into the wake flow patterns. Two full three-dimensional Navier–Stokes solvers for incompressible fluid flow, employing k − ϵ and k − ω turbulence closures, are used. The wind turbines are modeled as momentum absorbers by means of their thrust coefficient through the actuator disk approach. Alternative methods for estimating the reference wind speed in the calculation of the thrust are tested. The work presented in this paper is part of the work being undertaken within the UpWind Integrated Project that aims to develop the design tools for next generation of large wind turbines. In this part of UpWind, the performance of wind farm and wake models is being examined in complex terrain environment where there are few pre-existing relevant measurements. The focus of the work being carried out is to evaluate the performance of CFD models in large wind farm applications in complex terrain and to examine the development of the wakes in a complex terrain environment.

Modelling and forecasting fossil fuels, CO2 and electricity prices and their volatilities

In the current uncertain context that affects both the world economy and the energy sector, with the rapid increase in the prices of oil and gas and the very unstable political situation that affects some of the largest raw materials’ producers, there is a need for developing efficient and powerful quantitative tools that allow to model and forecast fossil fuel prices, CO2 emission allowances prices as well as electricity prices. This will improve decision making for all the agents involved in energy issues. Although there are papers focused on modelling fossil fuel prices, CO2 prices and electricity prices, the literature is scarce on attempts to consider all of them together. This paper focuses on both building a multivariate model for the aforementioned prices and comparing its results with those of univariate ones, in terms of prediction accuracy (univariate and multivariate models are compared for a large span of days, all in the first 4 months in 2011) as well as extracting common features in the volatilities of the prices of all these relevant magnitudes. The common features in volatility are extracted by means of a conditionally heteroskedastic dynamic factor model which allows to solve the curse of dimensionality problem that commonly arises when estimating multivariate GARCH models. Additionally, the common volatility factors obtained are useful for improving the forecasting intervals and have a nice economical interpretation. Besides, the results obtained and methodology proposed can be useful as a starting point for risk management or portfolio optimization under uncertainty in the current context of energy markets.

Distribution of the coal flow in the mill-duct system of the As Pontes Power Plant using CFD modeling

The efficiency of a Power Plant is affected by the distribution of the pulverized coal within the furnace. The coal, which is pulverized in the mills, is transported and distributed by the primary gas through the mill-ducts to the interior of the furnace. This is done with a double function: dry and enter the coal by different levels for optimizing the combustion in the sense that a complete combustion occurs with homogeneous heat fluxes to the walls. The mill-duct systems of a real Power Plant are very complex and they are not yet well understood. In particular, experimental data concerning the mass flows of coal to the different levels are very difficult to measure. CFD modeling can help to determine them. An Eulerian/Lagrangian approach is used due to the low solid–gas volume ratio.

Estudio de la refrigeración de una piscina de combustible gastado con un código CFD

Tras el accidente nuclear de Fukushima se demostró que las piscinas de combustible gastado en las centrales nucleares ven comprometida su refrigeración a largo plazo en caso de una pérdida total de energía eléctrica (SBO), ya que si experimentan un SBO de larga duración no existen a priori sistemas para mantener la refrigeración de los elementos combustibles que no dependan de los diésel de emergencia o de la red externa. En este trabajo se ha estudiado la refrigeración de una piscina de combustible gastado con el programa CFD STAR-CCM+, tanto en condiciones normales como en caso de pérdida del sistema de refrigeración. Posteriormente se ha evaluado la misma mediante el empleo de sistemas pasivos que permiten refrigerar los elementos combustibles durante cierto tiempo tras la pérdida del sistema de refrigeración y de una manera pasiva. De esta manera se consigue cierto margen antes de la entrada en ebullición del agua de la piscina, mejorándose por tanto la refrigeración de la misma. ABSTRACT. After the Fukushima nuclear accident, it was proved that the cooling of the current spent fuel pools are not sure for long term in case of a Station Blackout (SBO) Accident. If a long lasting blackout SBO occurs there are no systems available to keep cooling the spent fuel assemblies that do not rely on diesel generators or the external grid. During this thesis, the author has studied the spent fuel pool cooling, in ordinary conditions and if the spent fuel pool loses its cooling system, using the CFD program STAR-CCM+. Afterwards, the spent fuel pool cooling has been studied through the use of passive systems. Those two systems are able to cool the spent fuel assemblies in a passive way during a certain period of time after losing the cooling system. As a consequence, the pool´s water would boil later and the spent fuel pools safety would be enhanced.

Long Term Performance Study of a Direct Methanol Fuel Cell Fed with Alcohol Blends

The use of alcohol blends in direct alcohol fuel cells may be a more environmentally friendly and less toxic alternative to the use of methanol alone in direct methanol fuel cells. This paper assesses the behaviour of a direct methanol fuel cell fed with aqueous methanol, aqueous ethanol and aqueous methanol/ethanol blends in a long term experimental study followed by modelling of polarization curves. Fuel cell performance is seen to decrease as the ethanol content rises, and subsequent operation with aqueous methanol only partly reverts this loss of performance. It seems that the difference in the oxidation rate of these alcohols may not be the only factor affecting fuel cell performance.

Conceptual design of offshore platform supply vessel based on hybrid diesel generator-fuel cell power plant

Nowadays increasing fuel prices and upcoming pollutant emission regulations are becoming a growing concern for the shipping industry worldwide. While fuel prices will keep rising in future years, the new International Convention for the Prevention of Pollution from Ships (MARPOL) and Sulphur Emissions Control Areas (SECA) regulations will forbid ships to use heavy fuel oils at certain situations. To fulfil with these regulations, the next step in the marine shipping business will comprise the use of cleaner fuels on board as well as developing new propulsion concept. In this work a new conceptual marine propulsion system is developed, based on the integration of diesel generators with fuel cells in a 2850 metric tonne of deadweight platform supply vessel. The efficiency of the two 250 kW methanol-fed Solid Oxide Fuel Cell (SOFC) system installed on board combined with the hydro dynamically optimized design of the hull of the ship will allow the ship to successfully operate at certain modes of operation while notably reduce the pollutant emissions to the atmosphere. Besides the cogeneration heat obtained from the fuel cell system will be used to answer different heating needs on board the vessel

Simulación CFD de un accidente de caída de un elemento combustible en piscina

It is intended to provide a methodology of analysis more realistic this accident referred to in calculations of the license that requires fuel catastrophic break regardless of the height of the fall, with the consequent release of inventory analysers. Accidents that occurred in the past indicate that this hypothesis could be too conservative.

Análisis de viabilidad e implantación de sistemas de cogeneración en el sector residencial

En la actualidad la generación y utilización eficientes de la energía es el vector principal que permite el desarrollo sostenible en el marco ambiental, económico, seguro y rentable. Todo ello genera una necesidad en el ser humano de guiar a los avances tecnológicos hacia una manera cada vez más eficiente de generar nuestras necesidades básicas, como es el caso de la energía. La cogeneración ha sido uno de los resultados positivos en la búsqueda de la eficiencia energética, debido a tratarse de un sistema de producción simultánea de calor y electricidad partiendo inicialmente de un combustible como energía primaria. Es por ello, que en el presente proyecto se estudia, analiza y propone la posibilidad de implantar sistemas de cogeneración en el sector residencial, un sector que podría beneficiarse enormemente de los beneficios que ofrecen dichos sistemas. En una primera parte se analiza la tecnología de cogeneración y sus variantes, como son, la microcogeneración y la trigeneración. También se muestra la evolución legislativa que han sufrido estos sistemas. En una segunda parte se ha tomado un caso modelo, un edificio de 72 viviendas con sistema de calderas centralizado convencional, y se ha estudiado la posibilidad de implantar un sistema de cogeneración. Para ello se han calculado previamente las demandas energéticas del edificio y se han ido proponiendo diferentes modos de operación para cubrir dichas demandas por medio de sistemas de microcogeneración o cogeneración. Finalmente, una vez valoradas las opciones se muestra la elegida y se efectúa un análisis económico ABSTRACT Nowadays the efficient generation of energy is the main vector that allows sustainable development in environmental, economic, safety and cost effectiveness. All this generates a need in humans to lead to new technological advances towards an even more efficient way to generate our basic needs, such as energy. Cogeneration has been one of the positive results in the search for energy efficiency, due to the fact that it is a system of simultaneous production of heat and electricity initially starting from a primary energy fuel. It is for this reason that this project studies, analyzes and proposes the possibility of introducing cogeneration systems in the residential sector, a sector that could benefit greatly from the benefits offered by these systems. In the first part, cogeneration technology and its variants are analyzed, like, micro-cogeneration and trigeneration. The legislative evolutions that have suffered these systems are also displayed. In a second part, a model case has been taken; a building of 72 flats with conventional centralized boiler system, the possibility of introducing a cogeneration system has been studied. Previously the energy demands of the building have been calculated proposing different operating modes to meet those demands through micro-CHP or cogeneration systems. Finally, once the options are valued the chosen one is shown and an economic analysis is performed.