PROPAGATION OF SONIC BOOMS THROUGH A REAL, STRATIFIED ATMOSPHERE

Sonic boom propagation in a quiet) stratified) lossy atmosphere is the subject of this dissertation. Two questions are considered in detail: (1) Does waveform freezing occur? (2) Are sonic booms shocks in steady state? Both assumptions have been invoked in the past to predict sonic boom waveforms at the ground. A very general form of the Burgers equation is derived and used as the model for the problem. The derivation begins with the basic conservation equations. The effects of nonlinearity) attenuation and dispersion due to multiple relaxations) viscosity) and heat conduction) geometrical spreading) and stratification of the medium are included. When the absorption and dispersion terms are neglected) an analytical solution is available. The analytical solution is used to answer the first question. Geometrical spreading and stratification of the medium are found to slow down the nonlinear distortion of finite-amplitude waves. In certain cases the distortion reaches an absolute limit) a phenomenon called waveform freezing. Judging by the maturity of the distortion mechanism, sonic booms generated by aircraft at 18 km altitude are not frozen when they reach the ground. On the other hand, judging by the approach of the waveform to its asymptotic shape, N waves generated by aircraft at 18 km altitude are frozen when they reach the ground. To answer the second question we solve the full Burgers equation and for this purpose develop a new computer code, THOR. The code is based on an algorithm by Lee and Hamilton (J. Acoust. Soc. Am. 97, 906-917, 1995) and has the novel feature that all its calculations are done in the time domain, including absorption and dispersion. Results from the code compare very well with analytical solutions. In a NASA exercise to compare sonic boom computer programs, THOR gave results that agree well with those of other participants and ran faster. We show that sonic booms are not steady state waves because they travel through a varying medium, suffer spreading, and fail to approximate step shocks closely enough. Although developed to predict sonic boom propagation, THOR can solve other problems for which the extended Burgers equation is a good propagation model.

Multiphysics modelling of the metals casting process

Metals casting is a process governed by the interaction of a range of physical phenomena. Most computational models of this process address only what are conventionally regarded as the primary phenomena-heat conduction and solidification. However, to predict the formation of porosity (a factor of key importance in cast quality) requires the modelling of the interaction of the fluid flow, heat transfer, solidification and the development of stress-deformation in the solidified part of a component. In this paper, a model of the casting process is described which addresses all the main continuum phenomena involved in a coupled manner. The model is solved numerically using novel finite volume unstructured mesh techniques, and then applied to both the prediction of shape deformation (plus the subsequent formation of a gap at the metal-mould interface and its impact on the heat transfer behaviour) and porosity formation in solidifying metal components. Although the porosity prediction model is phenomenologically simplistic it is based on the interaction of the continuum phenomena and yields good comparisons with available experimental results. This work represents the first of the next generation of casting simulation tools to predict aspects of the structure of cast components.

Explicit finite volume modeling of aircraft anti-icing and de-icing

In this paper, a runback water and ice prediction model is extended to anti-icing and thermal de-icing situations. The resulting coupled equations that govern thin-film flow, ice accretion, and heat conduction in the multilayered system substrate-ice-water are solved using an explicit finite volume approach. The procedure is implemented in the three-dimensional icing code ICECREMO2, and both structured and unstructured grids can be considered. Numerical results are presented to compare the present code simulations to some data provided by other ice prediction codes and to show the capabilities of the present numerical tool.

Thermal transport in out-of-equilibrium quantum harmonic chains

We address the problem of heat transport in a chain of coupled quantum harmonic oscillators, exposed to the influences of local environments of various nature, stressing the effects that the specific nature of the environment has on the phenomenology of the transport process. We study in detail the behavior of thermodynamically relevant quantities such as heat currents and mean energies of the oscillators, establishing rigorous analytical conditions for the existence of a steady state, whose features we analyze carefully. In particular, we assess the conditions that should be faced to recover trends reminiscent of the classical Fourier law of heat conduction and highlight how such a possibility depends on the environment linked to our system.

Thermal transport across incommensurate phases in potassium selenate: photo-pyroelectric and calorimetric measurements

The thermal transport properties—thermal diffusivity, thermal conductivity and specific heat capacity—of potassium selenate crystal have been measured through the successive phase transitions, following the photo-pyroelectric thermal wave technique. The variation of thermal conductivity with temperature through the incommensurate (IC) phase of this crystal is measured. The enhancement in thermal conductivity in the IC phase is explained in terms of heat conduction by phase modes, and the maxima in thermal conductivity during transitions is due to enhancement in the phonon mean free path and the corresponding reduction in phonon scattering. The anisotropy in thermal conductivity and its variation with temperature are reported. The variation of the specific heat with temperature through the high temperature structural transition at 745 K is measured, following the differential scanning calorimetric method. By combining the results of photo-pyroelectric thermal wave methods and differential scanning calorimetry, the variation of the specific heat capacity with temperature through all the four phases of K2SeO4 is reported. The results are discussed in terms of phonon mode softening during transitions and phonon scattering by phase modes in the IC phase.

Thermoaktive Bauteilsysteme - Ein neuer simulationstechnischer Berechnungsansatz

Thermoaktive Bauteilsysteme sind Bauteile, die als Teil der Raumumschließungsflächen über ein integriertes Rohrsystem mit einem Heiz- oder Kühlmedium beaufschlagt werden können und so die Beheizung oder Kühlung des Raumes ermöglichen. Die Konstruktionenvielfalt reicht nach diesem Verständnis von Heiz, bzw. Kühldecken über Geschoßtrenndecken mit kern-integrierten Rohren bis hin zu den Fußbodenheizungen. Die darin enthaltenen extrem trägen Systeme werden bewußt eingesetzt, um Energieangebot und Raumenergiebedarf unter dem Aspekt der rationellen Energieanwendung zeitlich zu entkoppeln, z. B. aktive Bauteilkühlung in der Nacht, passive Raumkühlung über das kühle Bauteil am Tage. Gebäude- und Anlagenkonzepte, die träge reagierende thermoaktive Bauteilsysteme vorsehen, setzen im kompetenten und verantwortungsvollen Planungsprozeß den Einsatz moderner Gebäudesimulationswerkzeuge voraus, um fundierte Aussagen über Behaglichkeit und Energiebedarf treffen zu können. Die thermoaktiven Bauteilsysteme werden innerhalb dieser Werkzeuge durch Berechnungskomponenten repräsentiert, die auf mathematisch-physikalischen Modellen basieren und zur Lösung des bauteilimmanenten mehrdimensionalen instationären Wärmeleitungsproblems dienen. Bisher standen hierfür zwei unterschiedliche prinzipielle Vorgehensweisen zur Lösung zur Verfügung, die der physikalischen Modellbildung entstammen und Grenzen bzgl. abbildbarer Geometrie oder Rechengeschwindigkeit setzen. Die vorliegende Arbeit dokumentiert eine neue Herangehensweise, die als experimentelle Modellbildung bezeichnet wird. Über den Weg der Systemidentifikation können aus experimentell ermittelten Datenreihen die Parameter für ein kompaktes Black-Box-Modell bestimmt werden, das das Eingangs-Ausgangsverhalten des zugehörigen beliebig aufgebauten thermoaktiven Bauteils mit hinreichender Genauigkeit widergibt. Die Meßdatenreihen lassen sich über hochgenaue Berechnungen generieren, die auf Grund ihrer Detailtreue für den unmittelbaren Einsatz in der Gebäudesimulation ungeeignet wären. Die Anwendung der Systemidentifikation auf das zweidimensionale Wärmeleitungsproblem und der Nachweis ihrer Eignung wird an Hand von sechs sehr unterschiedlichen Aufbauten thermoaktiver Bauteilsysteme durchgeführt und bestätigt sehr geringe Temperatur- und Energiebilanzfehler. Vergleiche zwischen via Systemidentifikation ermittelten Black-Box-Modellen und physikalischen Modellen für zwei Fußbodenkonstruktionen zeigen, daß erstgenannte auch als Referenz für Genauigkeitsabschätzungen herangezogen werden können. Die Praktikabilität des neuen Modellierungsansatzes wird an Fallstudien demonstriert, die Ganzjahressimulationen unter Bauteil- und Betriebsvariationen an einem exemplarischen Büroraum betreffen. Dazu erfolgt die Integration des Black-Box-Modells in das kommerzielle Gebäude- und Anlagensimulationsprogramm CARNOT. Die akzeptablen Rechenzeiten für ein Einzonen-Gebäudemodell in Verbindung mit den hohen Genauigkeiten bescheinigen die Eignung der neuen Modellierungsweise.

Ein Beitrag zur Modellierung und Simulation von elektromechanischen Systemen

Die in dieser Arbeit durchgeführten Untersuchungen zeigen, daß es möglich ist, komplexe thermische Systeme, unter Verwendung der thermisch-, elektrischen Analogien, mit PSpice zu simulieren. Im Mittelpunkt der Untersuchungen standen hierbei Strangkühlkörper zur Kühlung von elektronischen Bauelementen. Es konnte gezeigt werden,daß alle Wärmeübertragungsarten, (Wärmeleitung, Konvektion und Wärmestrahlung) in der Simulation berücksichtigt werden können. Für die Berechnung der Konvektion wurden verschiedene Methoden hergeleitet. Diese gelten zum einen für verschiedene Kühlkörpergeometrien, wie z.B. ebene Flächen und Kühlrippenzwischenräume, andererseits unterscheiden sie sich, je nachdem, ob freie oder erzwungene Konvektion betrachtet wird. Für die Wärmestrahlung zwischen den Kühlrippen wurden verschiedenen Berechnungsmethoden entwickelt. Für die Simulation mit PSpice wurde die Berechnung der Wärmestrahlung zwischen den Kühlrippen vereinfacht. Es konnte gezeigt werden, daß die Fehler, die durch die Vereinfachung entstehen, vernachlässigbar klein sind. Für das thermische Verhalten einer zu kühlenden Wärmequelle wurde ein allgemeines Modell entworfen. Zur Bestimmung der Modellparameter wurden verschiedene Meßverfahren entwickelt. Für eine im Fachgebiet Elektromechanik entwickelte Wärmequelle zum Test von Kühlvorrichtungen wurde mit Hilfe dieser Meßverfahren eine Parameterbestimmung durchgeführt. Die Erstellung des thermischen Modells eines Kühlkörpers für die Simulation in PSpice erfordert die Analyse der Kühlkörpergeometrie. Damit diese Analyse weitestgehend automatisiert werden kann, wurden verschiedene Algorithmen unter Matlab entwickelt. Es wurde ein Algorithmus entwickelt, der es ermöglicht, den Kühlkörper in Elementarzellen zu zerlegen, die für die Erstellung des Simulationsmodells benötigt werden. Desweiteren ist es für die Simulation notwendig zu wissen, welche der Elementarzellen am Rand des Kühlkörpers liegen, welche der Elementarzellen an einem Kühlrippenzwischenraum liegen und welche Kühlkörperkanten schräg verlaufen. Auch zur Lösung dieser Aufgaben wurden verschiedene Algorithmen entwickelt. Diese Algorithmen wurden zu einem Programm zusammengefaßt, das es gestattet, unterschiedliche Strangkühlkörper zu simulieren und die Simulationsergebnisse in Form der Temperaturverteilung auf der Montagefläche des Kühlkörpers grafisch darzustellen. Es können stationäre und transiente Simulationen durchgeführt werden. Desweiteren kann der thermische Widerstand des Kühlkörpers RthK als Funktion der Verlustleistung der Wärmequelle dargestellt werden. Zur Verifikation der Simulationsergebnisse wurden Temperaturmessungen an Kühlkörpern durchgeführt und mit den Simulationsergebnissen verglichen. Diese Vergleiche zeigen, daß die Abweichungen im Bereich der Streuung der Temperaturmessung liegen. Das hier entwickelte Verfahren zur thermischen Simulation von Strangkühlkörpern kann somit als gut bewertet werden.

Atomistic-continuum modeling of ultrafast laser-induced melting of silicon targets

In this work, we present an atomistic-continuum model for simulations of ultrafast laser-induced melting processes in semiconductors on the example of silicon. The kinetics of transient non-equilibrium phase transition mechanisms is addressed with MD method on the atomic level, whereas the laser light absorption, strong generated electron-phonon nonequilibrium, fast heat conduction, and photo-excited free carrier diffusion are accounted for with a continuum TTM-like model (called nTTM). First, we independently consider the applications of nTTM and MD for the description of silicon, and then construct the combined MD-nTTM model. Its development and thorough testing is followed by a comprehensive computational study of fast nonequilibrium processes induced in silicon by an ultrashort laser irradiation. The new model allowed to investigate the effect of laser-induced pressure and temperature of the lattice on the melting kinetics. Two competing melting mechanisms, heterogeneous and homogeneous, were identified in our big-scale simulations. Apart from the classical heterogeneous melting mechanism, the nucleation of the liquid phase homogeneously inside the material significantly contributes to the melting process. The simulations showed, that due to the open diamond structure of the crystal, the laser-generated internal compressive stresses reduce the crystal stability against the homogeneous melting. Consequently, the latter can take a massive character within several picoseconds upon the laser heating. Due to the large negative volume of melting of silicon, the material contracts upon the phase transition, relaxes the compressive stresses, and the subsequent melting proceeds heterogeneously until the excess of thermal energy is consumed. A series of simulations for a range of absorbed fluences allowed us to find the threshold fluence value at which homogeneous liquid nucleation starts contributing to the classical heterogeneous propagation of the solid-liquid interface. A series of simulations for a range of the material thicknesses showed that the sample width we chosen in our simulations (800 nm) corresponds to a thick sample. Additionally, in order to support the main conclusions, the results were verified for a different interatomic potential. Possible improvements of the model to account for nonthermal effects are discussed and certain restrictions on the suitable interatomic potentials are found. As a first step towards the inclusion of these effects into MD-nTTM, we performed nanometer-scale MD simulations with a new interatomic potential, designed to reproduce ab initio calculations at the laser-induced electronic temperature of 18946 K. The simulations demonstrated that, similarly to thermal melting, nonthermal phase transition occurs through nucleation. A series of simulations showed that higher (lower) initial pressure reinforces (hinders) the creation and the growth of nonthermal liquid nuclei. For the example of Si, the laser melting kinetics of semiconductors was found to be noticeably different from that of metals with a face-centered cubic crystal structure. The results of this study, therefore, have important implications for interpretation of experimental data on the kinetics of melting process of semiconductors.

Crumpet structures: experimental and modelling studies

Crumpets are made by heating fermented batter on a hot plate at around 230°C. The characteristic structure dominated by vertical pores develops rapidly: structure has developed throughout around 75% of the product height within 30s, which is far faster than might be expected from transient heat conduction through the batter. Cooking is complete within around 3 min. Image analysis based on results from X-ray tomography shows that the voidage fraction is approximately constant and that there is continual coalescence between the larger pores throughout the product although there is also a steady level of small bubbles trapped within the solidified batter. We report here experimental studies which shed light on some of the mechanisms responsible for this structure, together with some models of key phenomena.Three aspects are discussed here: the role of gas (carbon dioxide and nitrogen) nuclei in initiating structure development; convective heat transfer inside the developing pores; and the kinetics of setting the batter into an elastic solid structure. It is shown conclusively that the small bubbles of carbon dioxide resulting from the fermentation stage play a crucial role as nuclei for pore development: without these nuclei, the result is not a porous structure, but rather a solid, elastic, inedible, gelatinized product. These nuclei are also responsible for the tiny bubbles which are set in the final product. The nuclei form the source of the dominant pore structure which is largely driven by the, initially explosive, release of water vapour from the batter together with the desorption of dissolved carbon dioxide. It is argued that the rapid evaporation, transport and condensation of steam within the growing pores provides an important mechanism, as in a heat pipe, for rapid heat transfer, and models for this process are developed and tested. The setting of the continuous batter phase is essential for final product quality: studies using differential scanning calorimetry and on the kinetics of change in the visco-elastic properties of the batter suggest that this process is driven by the kinetics of gelatinization. Unlike many thermally driven food processes the rates of heating are such that gelatinization kinetics cannot be neglected. The implications of these results for modelling and for the development of novel structures are discussed.

An urban parameterization for a global climate model. Part I: Formulation and evaluation for two cities

Urbanization, the expansion of built-up areas, is an important yet less-studied aspect of land use/land cover change in climate science. To date, most global climate models used to evaluate effects of land use/land cover change on climate do not include an urban parameterization. Here, the authors describe the formulation and evaluation of a parameterization of urban areas that is incorporated into the Community Land Model, the land surface component of the Community Climate System Model. The model is designed to be simple enough to be compatible with structural and computational constraints of a land surface model coupled to a global climate model yet complex enough to explore physically based processes known to be important in determining urban climatology. The city representation is based upon the “urban canyon” concept, which consists of roofs, sunlit and shaded walls, and canyon floor. The canyon floor is divided into pervious (e.g., residential lawns, parks) and impervious (e.g., roads, parking lots, sidewalks) fractions. Trapping of longwave radiation by canyon surfaces and solar radiation absorption and reflection is determined by accounting for multiple reflections. Separate energy balances and surface temperatures are determined for each canyon facet. A one-dimensional heat conduction equation is solved numerically for a 10-layer column to determine conduction fluxes into and out of canyon surfaces. Model performance is evaluated against measured fluxes and temperatures from two urban sites. Results indicate the model does a reasonable job of simulating the energy balance of cities.

A thermoelastic system of memory type in noncylindrical domains

In this paper, we consider a hyperbolic thermoelastic system of memory type in domains with moving boundary. The problem models vibrations of an elastic bar under thermal effects according to the heat conduction law of Gurtin and Pipkin. Global existence is proved by using the penalty method of Lions. (c) 2007 Elsevier Inc. All rights reserved.

LOSS of EXPONENTIAL STABILITY FOR A THERMOELASTIC SYSTEM WITH MEMORY

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Experimental and theoretical study of workpiece temperature when end milling hardened steels using (TiAl)N-coated and PcBN-tipped tools

The present work shows an experimental and theoretical study on heat flow when end milling, at high-speed, hardened steels applied to moulds and dies. AISI H13 and AISI D2 steels were machined with two types of ball nose end mills: coated with (TiAl)N and tipped with PcBN. The workpiece geometry was designed to simulate tool-workpiece interaction in real situations found in mould industries, in which complex surfaces and thin walls are commonly machined. The compressed and cold air cooling systems were compared to dry machining Results indicated a relatively small temperature variation, with higher range when machining AISI D2 with PcBN-tipped end mill. All cooling systems used demonstrated good capacity to remove heat from the machined surface, especially the cold air. Compressed air was the most indicated to keep workpiece at relatively stable temperature. A theoretical model was also proposed to estimate the energy transferred to the workpiece (Q) and the average convection coefficient ((h) over bar) for the cooling systems used. The model used a FEM simulation and a steepest decent method to find the best values for both variables. (c) 2007 Elsevier B.V. All rights reserved.

Heating time required for an MgO-graphite refractory during a hot modulus of rupture test

In the work described in the present paper, an analytical solution of the general heat conduction equation was employed to assay the temperature profile inside a solid slab which is initially at room temperature and is suddenly plunged into a fluid maintained at a high temperature. The results were then extrapolated to a simulation of a hot modulus of rupture test of typical MgO-graphite refractory samples containing different amounts of graphite in order to evaluate how fast the temperature equilibrates inside the test specimens. Calculations indicated that, depending on the graphite content, the time to full temperature homogenization was in the range of 80 to 200 s. These findings are relevant to the high temperature testing of such refractories in oxidizing conditions in view of the graphite oxidation risks in the proper evaluation of the hot mechanical properties.

Magnetothermal instabilities in magnetized anisotropic plasmas

Using the 16-moment transport equations for an ideal anisotropic collisionless plasma we analyze the influence of pressure anisotropy on the magnetothermal (MTI) and heat-flux-driven buoyancy (HBI) instabilities. We calculate the dispersion relation and the growth rates for these instabilities in the presence of a background heat flux and for configurations with static pressure anisotropy, finding that when the frequency at which heat conduction acts is much larger than any other frequency in the system (i.e. weak magnetic field) the pressure anisotropy has no effect on the MTI/HBI, provided the degree of anisotropy is small. In contrast, when this ordering of timescales does not apply the instability criteria depend on pressure anisotropy.