Heating the Early Universe : Numerical Methods and Their Analysis

During the epoch when the first collapsed structures formed (6<z<50) our Universe went through an extended period of changes. Some of the radiation from the first stars and accreting black holes in those structures escaped and changed the state of the Intergalactic Medium (IGM). The era of this global phase change in which the state of the IGM was transformed from cold and neutral to warm and ionized, is called the Epoch of Reionization.In this thesis we focus on numerical methods to calculate the effects of this escaping radiation. We start by considering the performance of the cosmological radiative transfer code C2-Ray. We find that although this code efficiently and accurately solves for the changes in the ionized fractions, it can yield inaccurate results for the temperature changes. We introduce two new elements to improve the code. The first element, an adaptive time step algorithm, quickly determines an optimal time step by only considering the computational cells relevant for this determination. The second element, asynchronous evolution, allows different cells to evolve with different time steps. An important constituent of methods to calculate the effects of ionizing radiation is the transport of photons through the computational domain or ``ray-tracing''. We devise a novel ray tracing method called PYRAMID which uses a new geometry - the pyramidal geometry. This geometry shares properties with both the standard Cartesian and spherical geometries. This makes it on the one hand easy to use in conjunction with a Cartesian grid and on the other hand ideally suited to trace radiation from a radially emitting source. A time-dependent photoionization calculation not only requires tracing the path of photons but also solving the coupled set of photoionization and thermal equations. Several different solvers for these equations are in use in cosmological radiative transfer codes. We conduct a detailed and quantitative comparison of four different standard solvers in which we evaluate how their accuracy depends on the choice of the time step. This comparison shows that their performance can be characterized by two simple parameters and that the C2-Ray generally performs best.

Domain decomposition methods for some aerodynamic noise problems

Abstract not available

Resistive Embedded Heating for Homogeneous Curing of Composite Damage Repair

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

Coupled CFD Shape Optimization for aerodynamic profiles

The present document deals with the optimization of shape of aerodynamic profiles -- The objective is to reduce the drag coefficient on a given profile without penalising the lift coefficient -- A set of control points defining the geometry are passed and parameterized as a B-Spline curve -- These points are modified automatically by means of CFD analysis -- A given shape is defined by an user and a valid volumetric CFD domain is constructed from this planar data and a set of user-defined parameters -- The construction process involves the usage of 2D and 3D meshing algorithms that were coupled into own- code -- The volume of air surrounding the airfoil and mesh quality are also parametrically defined -- Some standard NACA profiles were used by obtaining first its control points in order to test the algorithm -- Navier-Stokes equations were solved for turbulent, steady-state ow of compressible uids using the k-epsilon model and SIMPLE algorithm -- In order to obtain data for the optimization process an utility to extract drag and lift data from the CFD simulation was added -- After a simulation is run drag and lift data are passed to the optimization process -- A gradient-based method using the steepest descent was implemented in order to define the magnitude and direction of the displacement of each control point -- The control points and other parameters defined as the design variables are iteratively modified in order to achieve an optimum -- Preliminary results on conceptual examples show a decrease in drag and a change in geometry that obeys to aerodynamic behavior principles

3D Ionospheric Effects on HF Propagation and Heating

The thesis uses a three-dimensional, first-principles model of the ionosphere in combination with High Frequency (HF) raytracing model to address key topics related to the physics of HF propagation and artificial ionospheric heating. In particular: 1. Explores the effect of the ubiquitous electron density gradients caused by Medium Scale Traveling Ionospheric Disturbances (MSTIDs) on high-angle of incidence HF radio wave propagation. Previous studies neglected the all-important presence of horizontal gradients in both the cross- and down-range directions, which refract the HF waves, significantly changing their path through the ionosphere. The physics-based ionosphere model SAMI3/ESF is used to generate a self-consistently evolving MSTID that allows for the examination of the spatio-temporal progression of the HF radio waves in the ionosphere. 2. Tests the potential and determines engineering requirements for ground- based high power HF heaters to trigger and control the evolution of Equatorial Spread F (ESF). Interference from ESF on radio wave propagation through the ionosphere remains a critical issue on HF systems reliability. Artificial HF heating has been shown to create plasma density cavities in the ionosphere similar to those that may trigger ESF bubbles. The work explores whether HF heating may trigger or control ESF bubbles. 3. Uses the combined ionosphere and HF raytracing models to create the first self-consistent HF Heating model. This model is utilized to simulate results from an Arecibo experiment and to provide understanding of the physical mechanism behind observed phenomena. The insights gained provide engineering guidance for new artificial heaters that are being built for use in low to middle latitude regions. In accomplishing the above topics: (i) I generated a model MSTID using the SAMI3/ESF code, and used a raytrace model to examine the effects of the MSTID gradients on radio wave propagation observables; (ii) I implemented a three- dimensional HF heating model in SAMI3/ESF and used the model to determine whether HF heating could artificially generate an ESF bubble; (iii) I created the first self-consistent model for artificial HF heating using the SAMI3/ESF ionosphere model and the MoJo raytrace model and ran a series of simulations that successfully modeled the results of early artificial heating experiments at Arecibo.

The use of solar colletor for heating and electricity for a single family house

The presented work is related to the use of solar energy for the needs of heating and electricity for a single house located in Poland. Electricity will provided by energy conversion in the turbine by means of Organic Rankine Cycle (ORC), in which the operating medium (water heated in solar collector) is heating refrigerator in the heating exchanger. The solar installation is integrated with heat accumulator and wood boiler, which is used in the situation that collector is not enough to fill requirements of thermal comfort. There are chosen also all the necessary components of the system. In the work is also performed the economic assessment, by F chart method, to evaluate the profitability of the project, taking into total costs and savings.

Effect of Geometric Uncertainties on the Aerodynamic Characteristic of Offshore Wind Turbine Blades

Offshore wind turbines operate in a complex unsteady flow environment which causes unsteady aerodynamic loads. The unsteady flow environment is characterized by a high degree of uncertainty. In addition, geometry variations and material imperfections also cause uncertainties in the design process. Probabilistic design methods consider these uncertainties in order to reach acceptable reliability and safety levels for offshore wind turbines. Variations of the rotor blade geometry influence the aerodynamic loads which also affect the reliability of other wind turbine components. Therefore, the present paper is dealing with geometric uncertainties of the rotor blades. These can arise from manufacturing tolerances and operational wear of the blades. First, the effect of geometry variations of wind turbine airfoils on the lift and drag coefficients are investigated using a Latin hypercube sampling. Then, the resulting effects on the performance and the blade loads of an offshore wind turbine are analyzed. The variations of the airfoil geometry lead to a significant scatter of the lift and drag coefficients which also affects the damage-equivalent flapwise bending moments. In contrast to that, the effects on the power and the annual energy production are almost negligible with regard to the assumptions made.

Aerodynamic behavior of an airfoil with morphing trailing edge for wind turbine applications

The length of wind turbine rotor blades has been increased during the last decades. Higher stresses arise especially at the blade root because of the longer lever arm. One way to reduce unsteady blade-root stresses caused by turbulence, gusts, or wind shear is to actively control the lift in the blade tip region. One promising method involves airfoils with morphing trailing edges to control the lift and consequently the loads acting on the blade. In the present study, the steady and unsteady behavior of an airfoil with a morphing trailing edge is investigated. Two-dimensional Reynolds-Averaged Navier-Stokes (RANS) simulations are performed for a typical thin wind turbine airfoil with a morphing trailing edge. Steady-state simulations are used to design optimal geometry, size, and deflection angles of the morphing trailing edge. The resulting steady aerodynamic coefficients are then analyzed at different angles of attack in order to determine the effectiveness of the morphing trailing edge. In order to investigate the unsteady aerodynamic behavior of the optimal morphing trailing edge, time-resolved RANS-simulations are performed using a deformable grid. In order to analyze the phase shift between the variable trailing edge deflection and the dynamic lift coefficient, the trailing edge is deflected at four different reduced frequencies for each different angle of attack. As expected, a phase shift between the deflection and the lift occurs. While deflecting the trailing edge at angles of attack near stall, additionally an overshoot above and beyond the steady lift coefficient is observed and evaluated.

Geothermal district heating networks: modelling novel operational strategies incorporating heat storage

The value of integrating a heat storage into a geothermal district heating system has been investigated. The behaviour of the system under a novel operational strategy has been simulated focusing on the energetic, economic and environmental effects of the new strategy of incorporation of the heat storage within the system. A typical geothermal district heating system consists of several production wells, a system of pipelines for the transportation of the hot water to end-users, one or more re-injection wells and peak-up devices (usually fossil-fuel boilers). Traditionally in these systems, the production wells change their production rate throughout the day according to heat demand, and if their maximum capacity is exceeded the peak-up devices are used to meet the balance of the heat demand. In this study, it is proposed to maintain a constant geothermal production and add heat storage into the network. Subsequently, hot water will be stored when heat demand is lower than the production and the stored hot water will be released into the system to cover the peak demands (or part of these). It is not intended to totally phase-out the peak-up devices, but to decrease their use, as these will often be installed anyway for back-up purposes. Both the integration of a heat storage in such a system as well as the novel operational strategy are the main novelties of this thesis. A robust algorithm for the sizing of these systems has been developed. The main inputs are the geothermal production data, the heat demand data throughout one year or more and the topology of the installation. The outputs are the sizing of the whole system, including the necessary number of production wells, the size of the heat storage and the dimensions of the pipelines amongst others. The results provide several useful insights into the initial design considerations for these systems, emphasizing particularly the importance of heat losses. Simulations are carried out for three different cases of sizing of the installation (small, medium and large) to examine the influence of system scale. In the second phase of work, two algorithms are developed which study in detail the operation of the installation throughout a random day and a whole year, respectively. The first algorithm can be a potentially powerful tool for the operators of the installation, who can know a priori how to operate the installation on a random day given the heat demand. The second algorithm is used to obtain the amount of electricity used by the pumps as well as the amount of fuel used by the peak-up boilers over a whole year. These comprise the main operational costs of the installation and are among the main inputs of the third part of the study. In the third part of the study, an integrated energetic, economic and environmental analysis of the studied installation is carried out together with a comparison with the traditional case. The results show that by implementing heat storage under the novel operational strategy, heat is generated more cheaply as all the financial indices improve, more geothermal energy is utilised and less fuel is used in the peak-up boilers, with subsequent environmental benefits, when compared to the traditional case. Furthermore, it is shown that the most attractive case of sizing is the large one, although the addition of the heat storage most greatly impacts the medium case of sizing. In other words, the geothermal component of the installation should be sized as large as possible. This analysis indicates that the proposed solution is beneficial from energetic, economic, and environmental perspectives. Therefore, it can be stated that the aim of this study is achieved in its full potential. Furthermore, the new models for the sizing, operation and economic/energetic/environmental analyses of these kind of systems can be used with few adaptations for real cases, making the practical applicability of this study evident. Having this study as a starting point, further work could include the integration of these systems with end-user demands, further analysis of component parts of the installation (such as the heat exchangers) and the integration of a heat pump to maximise utilisation of geothermal energy.

Error estimation and adaptive mesh refinement for aerodynamic flows

This lecture course covers the theory of so-called duality-based a posteriori error estimation of DG finite element methods. In particular, we formulate consistent and adjoint consistent DG methods for the numerical approximation of both the compressible Euler and Navier-Stokes equations; in the latter case, the viscous terms are discretized based on employing an interior penalty method. By exploiting a duality argument, adjoint-based a posteriori error indicators will be established. Moreover, application of these computable bounds within automatic adaptive finite element algorithms will be developed. Here, a variety of isotropic and anisotropic adaptive strategies, as well as $hp$-mesh refinement will be investigated.

Analysis and improvement of the envelope, heating and ventilation of a low energy building with district heating

In this thesis project, a building in Vegagatan 12, Gävle has been analysed in order to see why it does consume more energy than it was expected. This building is a low energy building certified by Miljöbyggnad and it should use less than 55kWh/m2 year and nowadays it is using 62.23 kWh/m2. To get the needed data, some information about the building has been gathered, some measurements have been done in the building and some calculations have been done with those measurements. Finally, some possible solutions have been offered to reduce the energy use of the building. Insulating the floor, the pipes and the walls, reducing the indoor temperature in winter... All of these changes need the help of environmentally friendly attitudes, which is a very important fact in low energy buildings.

Aerodynamic parametric analysis of an unconventional joined-wing aircraft configuration

In this paper, the experimental results of an unconventional joined-wing aircraft configuration are presented. The test model uses two different wings, forward and rear, both joined in tandem and forming diamond shapes both in plant and front views. The wings are joined in such a way that it is possible to change the rear wing dihedral angle values and the rear wing sweep angle values in 25 different positions that modify the relative distance and the relative height between the wings. To measure the system aerodynamic coefficients itis necessary to perform wind tunnel tests. The datapresented corresponds to the lift, drag and induced drag aerodynamic coefficients, as well as the aerodynamic efficiency and the parameter for minimum required power, from the calculated values of the lift and drag time series measured by a 6-axis force and torque sensor. The results show the influence on the aerodynamic coefficients of the rear wing sweep and dihedral angles parameters. As a main result, it can be concluded that, in general terms, the lift and induced drag aerodynamic coefficients values decrease as both the distance and height between the wings increase, on the other hand, the total drag aerodynamic coefficient decreases if the distance between the wings increases, but nevertheless shows a slight tendency to increase if the height of the rear wing increases, whereas the aerodynamic efficiency and the parameter for minimum required power increase if the distance between the wings increases

Investigation of capacitive discharge heating of metallic glasses

In recent years, the discovery of bulk metallic glasses with exceptional properties has generated much interest. One of their most intriguing features is their capacity for viscous flow above the glass transition temperature. This characteristic allows metallic glasses to be formed like plastics at modest temperatures. However, crystallization of supercooled metallic liquids in the best bulk metallic glass-formers is much more rapid than in most polymers and silicate glass-forming liquids. The short times to crystallization impairs experimentation on and processing of supercooled glass-forming metallic liquids. A technique to rapidly and uniformly heat metallic glasses at rates of 10⁵ to 10⁶ kelvin per second is presented. A capacitive discharge is used to ohmically heat metallic glasses to temperatures in the super cooled liquid region in millisecond time-scales. By heating samples rapidly, the most time-consuming step in experiments on supercooled metallic liquids is reduced orders of magnitude in length. This allows for experimentation on and processing of metallic liquids in temperature ranges that were previously inaccessible because of crystallization.

A variety of forming techniques, including injection molding and forging, were coupled with capacitive discharge heating to produce near net-shaped metallic glass parts. In addition, a new forming technique, which combines a magnetic field with the heating current to produce a forming force, was developed. Viscosities were measured in previously inaccessible temperature ranges using parallel plate rheometry combined with capacitive discharge heating. Lastly, a rapid pulse calorimeter was developed with this technique to investigate the thermophysical behavior of metallic glasses at these rapid heating rates.