The study on adaptive Cartesian grid methods for compressible flow and their applications

This research is mainly focused on the development of the adaptive Cartesian grid methods for compressibl  e flow. At first, the ghost cell method and its applications for inviscid compressible flow on adaptive tree Cartesian grid are developed. The proposed method is successfully used to evaluate various inviscid compressible flows around complex bodies. The mass conservation of the method is also studied by numerical analysis. The extension to three-dimensional flow is presented. Then, an h-adaptive Runge–Kutta discontinuous Galerkin (RKDG) method is presented in detail for the development of high accuracy numerical method under Cartesian grid. This method combined with the ghost cell immersed boundary method is also validated by well documented test problems involving both steady and unsteady compressible flows over complex bodies in a wide range of Mach numbers. In addition, in order to suppress the failure of preserving positivity of density or pressure, which may cause blow-ups of the high order numerical algorithms, a positivity-preserving limiter technique coupled with h-adaptive RKDG method is developed. Such a method has been successfully implemented to study flows with the large Mach number, strong shock/obstacle interactions and shock diffraction. The extension of the method to viscous flow under the adaptive Cartesian grid with hybrid overlapping bodyfitted grid is developed. The method is validated by benchmark problems and has been successfully implemented to study airfoil with ice accretion. Finally, based on an open source code, the detached eddy simulation (DES) is developed for massive separation flow, and it is used to perform the research on aerodynamic performance analysis over the wing with ice accretion

Kartezyen hesaplama ağları kullanılarak iki boyutlu sıkıştırılabilir akışlar için laminar navier-stokes çözücüsü geliştirilmesi

TÜBİTAK MAG01.05.201

Kartezyen hesaplama ağları için üç boyutlu Euler çözücüsü geliştirilmesi

TÜBİTAK MAG01.08.200

Topics in Energy Release and Particle Acceleration in the Heliosphere

This thesis investigates both the release of energy in solar flares, and the acceleration and transport of particles in various astrophysical situations. While numerical simulations are central to this thesis, these are always motivated by analytical arguments. A review of flare energy release is given in Chapter 2, with results presented in Chapters 3 and 4. The main goal of the flare work is to investigate the effect of viscosity on energy release rates. Scaling arguments and exact solutions of the magnetohydrodynamic equations are used to interpret the results of two-dimensional numerical simulations of magnetic reconnection. The results support viscous energy dissipation accounting for a significant fraction of flare energy release. Chapter 5 contains an introduction to astrophysical particle acceleration, using the Fokker-Planck formulation. The theory introduced in this chapter is used to study electron transport in solar flare loops (Section 5.5). A key aspect of the analysis is the expression of the Fokker-Planck equation as a system of stochastic differential equations. A generalisation to the flare loop hard X-ray emission prediction of Conway et al. (1998) is obtained, giving a stronger dependence on density for dispersed initial distributions. Chapter 6 uses the methods of the previous chapter to study the acceleration of cosmic-rays at the heliospheric termination shock. The applicability of the focused acceleration mechanism of Schlickeiser and Shalchi (2008) is examined using numerical simulations, which are interpreted using analytical arguments based on averaging the stochastic equations. The results show significant limitations in assuming a near-isotropic distribution, a requirement for the focused acceleration mechanism. In addition, momentum diffusion provides a significant effect that cannot be neglected. The theory is extended to include focused deceleration and pure momentum diffusion

Computational Simulation of Detonation Waves and Model Reduction for Reacting Flows

Ph.DDOCTOR OF PHILOSOPH

Fundamental studies of flame propagation in lean-burn natural gas engines

Lean-burn natural gas engines offer enhanced thermal efficiencies and reduced soot and NOx emissions. However, cycle-to-cycle variability in combustion that can result from unreliable ignition, variability in equivalence ratio and quenching is a challenge. Reliability of ignition can be improved by employing a dual-fuel ignition strategy in which a small quantity of diesel fuel is injected to initiate ignition. Computational studies of n-heptane/methane-air mixing layers are performed to provide insight into the fundamental physics of dual-fuel ignition. The results show that the characteristic time required for steady premixed flame propagation has three components: time for autoignition to occur, time for peak temperature to be achieved following autoignition, and time for steady flame propagation in the premixed fuel/air mixture to be achieved. The autoignition time correlates well with pressure and temperature of the unburned premixed charge. The time to achieve peak temperature is relatively short, but correlates with mixing layer thickness and premixed equivalence ratio. The time to achieve steady propagation correlates with mixing layer thickness and laminar flame speed and thickness. Subsequent work focuses on turbulent flame propagation in lean homogeneous mixtures by employing direct numerical simulations (DNS) under conditions that are relevant to lean-burn engines. Attention is specifically focused on the turbulent flame speed (ST) as a parameter of interest because of its importance in modeling combustion in engines. The studies are carried out in the thin reaction zone (TRZ) regime of turbulent premixed combustion. Normalized turbulence intensity (urms/SL) varies from 2 to 25 and the ratio of integral length scale to flame thickness (L o/δL) varies from 3.2 to 12.8. Initial studies show that the normalized turbulent flame speed (ST/SL) depends on more parameters than urms/SL suggested by some models. Although it is known that the turbulent flame speed varies with equivalence ratio, it is shown that the normalized turbulent flame speed does not change with equivalence ratio provided the Karlovitz (Ka) and Damköhler (Da) numbers are fixed. This suggests that Kaand/or Da are important parameters in characterizing the turbulent flame speed. Furthermore, ST/SL can be related to the flame area enhancement AT/AL and an efficiency factor Io which is close to unity. AT/AL is raised by increasing turbulent Reynolds number ReT and by reducing Ka. Increasing ReT leads to a broader spectrum of turbulent eddies that generate flame surface area. Increasing Karesults in fine wrinkling at the expense of larger scale wrinkling. This results in a net reduction in the effective surface area enhancement. Based on these insights, a correlation for ST that shows a dependence on Re T and Ka is proposed. Modeling of the Flame Surface Density (FSD) evolution is also considered. FSD is influenced by tangential strain rate and flame displacement speed. Surface averaged tangential strain rate is found to scale linearly with Ka. The effects of Ka on flame displacement speed are modeled using a Probability Density Function (PDF) based approach. The effects of premixed combustion on turbulence are investigated. For flames in the TRZ regime, the turbulence kinetic energy (TKE) decays monotonically across the flame brush. Scaling analyses of the terms in the transport equation of TKE reveal that viscous dissipation is the dominant contribution in the TKE equation. The relative importance of the other terms in the TKE equation decreases with increasing Ka

Development of an image processing method for automated, non-invasive and scale-independent monitoring of adherent cell cultures

Adherent cell culture is a key experimental method for biological investigations in diverse areas such as developmental biology, drug discovery and biotechnology. Light microscopy-based methods, for example phase contrast microscopy (PCM), are routinely used for visual inspection of adherent cells cultured in transparent polymeric vessels. However, the outcome of such inspections is qualitative and highly subjective. Analytical methods that produce quantitative results can be used but often at the expense of culture integrity or viability. In this work, an imaging-based strategy to adherent cell cultures monitoring was investigated. Automated image processing and analysis of PCM images enabled quantitative measurements of key cell culture characteristics. Two types of segmentation algorithms for the detection of cellular objects on PCM images were evaluated. The first one, based on contrast filters and dynamic programming was quick (<1s per 1280×960 image) and performed well for different cell lines, over a wide range of imaging conditions. The second approach, termed ‘trainable segmentation’, was based on machine learning using a variety of image features such as local structures and symmetries. It accommodated complex segmentation tasks while maintaining low processing times (<5s per 1280×960 image). Based on the output from these segmentation algorithms, imaging-based monitoring of a large palette of cell responses was demonstrated, including proliferation, growth arrest, differentiation, and cell death. This approach is non-invasive and applicable to any transparent culture vessel, including microfabricated culture devices where a lack of suitable analytical methods often limits their applicability. This work was a significant contribution towards the establishment of robust, standardised, and affordable monitoring methods for adherent cell cultures. Finally, automated image processing was combined with computer-controlled cultures in small-scale devices. This provided a first demonstration of how adaptive culture protocols could be established; i.e. culture protocols which are based on cellular response instead of arbitrary time points

Explorative approach to the dust evolution in binary star systems

This thesis discusses dust formation in binary systems, in particular for binary systems consisting of a Mira like star and a brown dwarf. A Mira-like star is an intermediate mass star in a late stage of their stellar evolution on the Asymptotic Giant Branch (AGB), and a brown dwarf, is a sub-stellar object with a mass below that necessary to maintain hydrogen-burning nuclear fusion reactions in their cores. In their radial pulsating elevated convective atmospheres, Mira-stars often develop strong stellar winds, which are driven by radiation pressure on the dust and lead to a substantial mass-loss of the star. Stellar winds are of central importance for the development of medium-heavy-AGB stars. Also, they are a reliable source for the production of dust particles and heavy elements for the interstellar medium and the chemical evolution of galaxies. In fact, most stars are in binary or multiple star systems.For a complete description of the processes of dust formation in binary star systems it is necessary to study the perturbative influence of a second star in the vicinity of a AGB-star with a strong stellar wind. This endeavour is embedded in a long standing tradition of scientific investigation of dust formation at the Zentrum für Astronomie und Astrophysik (ZAA) at the Technical University of Berlin...thesi

A novel dual-spin actuation mechanism for small calibre, spin stabilised, guided projectiles

© Cranfield University 2022. All rights reserved. No part of this publication may be reproduced without the written permission of the author and copyright holderSmall calibre projectiles are spin-stabilised to increase ballistic stability, often at high frequencies. Due to hardware limitations, conventional actuators and meth ods are unable to provide satisfactory control at such high frequencies. With the reduced volume for control hardware and increased financial cost, incorporating traditional guid ance methods into small-calibre projectiles is inherently difficult. This work presents a novel method of projectile control which addresses these issues and conducts a systems level analysis of the underlying actuation mechanism. The design is shown to be a viable alternative to traditional control methods, Firstly, a 7 Degree-of-Freedom (DoF) dynamic model is created for dual-spin pro jectiles, including aerodynamic coefficients. The stability of dual-spin projectiles, gov erned by the gyroscopic and dynamic stability factors is given, discussed and unified across available literature. The model is implemented in a Matlab/Simulink simulation environ ment, which is in turn validated against a range of academic literature and experimental test data. The novel design and fundamental operating principle are presented. The actuation mechanism (AM) is then mathematically formulated from both a velocity change (∆V ) and a lateral acceleration (a˜) perspective. A set of axioms are declared and verified using the 7-DoF model, showing that the inherently discrete system behaviour can be controlled continuously via these control variables, ∆V or a˜. Control state switching is simplified to be instantaneous, then expanded to be generically characterised by an arbitrarily complex mathematical function. A detailed investigation, parametric analysis and sensitivity study is undertaken to understand the system behaviour. A Monte Carlo procedure is described, which is used to compare the correction cap abilities of different guidance laws (GLs). A bespoke Zero-Effort-Miss (ZEM) based GLis synthesised from the mathematical formulation of the AM, with innately more know ledge of the system behaviour, which allows superior error correction. This bespoke GL is discussed in detail, a parametric study is undertaken, and both the GL parameters and PID controller gains are optimised using a genetic algorithm. Artificial Intelligence (AI) Reinforcement learning methods are used to emulate a GL, as well as controlling the AM and operating as a GL, simultaneously. The novel GLs are compared against a traditional proportional navigation GL in a nominal system and all GLs were able to control the AMs, reducing the miss distance to a satisfactory margin. The ZEM-based GL provided superior correction to the AI GL, which in turn provided superior correction over proportional navigation. Example CAD models are shown, and the stability analysis is conducted on the geometry. The CAD model is then used in CFD simulations to determine aerodynamic coefficients for use in the 7-DoF dynamic model. The novel control method was able to reduce the 95% dispersion diameter of a traditional ballistic 7.62mm projectile from 70mm to 33mm. Statistical data analysis showed there was no significant correlation or bias present in either the nominal or 7-DoF dispersion patterns. This project is co-sponsored by BAE Systems and ESPRC (ref. 1700064). The con tents of this thesis are covered by patent applications GB2011850.1, GB 2106035.5 and EP 20275128.5. Two papers are currently published (DOI: 10.1016/j.dt.2019.06.003, the second DOI is pending) and one is undergoing peer review..PH