Sharply curved turn around duct flow predictions using spectral partitioning of the turbulent kinetic energy and a pressure modified wall law

Computational predictions of turbulent flow in sharply curved 180 degree turn around ducts are presented. The CNS2D computer code is used to solve the equations of motion for two-dimensional incompressible flows transformed to a nonorthogonal body-fitted coordinate system. This procedure incorporates the pressure velocity correction algorithm SIMPLE-C to iteratively solve a discretized form of the transformed equations. A multiple scale turbulence model based on simplified spectral partitioning is employed to obtain closure. Flow field predictions utilizing the multiple scale model are compared to features predicted by the traditional single scale k-epsilon model. Tuning parameter sensitivities of the multiple scale model applied to turn around duct flows are also determined. In addition, a wall function approach based on a wall law suitable for incompressible turbulent boundary layers under strong adverse pressure gradients is tested. Turn around duct flow characteristics utilizing this modified wall law are presented and compared to results based on a standard wall treatment

Turbulent-laminar patterns in shear flows without walls

Turbulent-laminar intermittency, typically in the form of bands and spots, is a ubiquitous feature of the route to turbulence in wall-bounded shear flows. Here we study the idealised shear between stress-free boundaries driven by a sinusoidal body force and demonstrate quantitative agreement between turbulence in this flow and that found in the interior of plane Couette flow -- the region excluding the boundary layers. Exploiting the absence of boundary layers, we construct a model flow that uses only four Fourier modes in the shear direction and yet robustly captures the range of spatiotemporal phenomena observed in transition, from spot growth to turbulent bands and uniform turbulence. The model substantially reduces the cost of simulating intermittent turbulent structures while maintaining the essential physics and a direct connection to the Navier-Stokes equations. We demonstrate the generic nature of this process by introducing stress-free equivalent flows for plane Poiseuille and pipe flows which again capture the turbulent-laminar structures seen in transition.Comment: 13 pages, 9 figure

Simplified Two-Time Step Method for Calculating Combustion and Emission Rates of Jet-A and Methane Fuel With and Without Water Injection

A simplified kinetic scheme for Jet-A, and methane fuels with water injection was developed to be used in numerical combustion codes, such as the National Combustor Code (NCC) or even simple FORTRAN codes. The two time step method is either an initial time averaged value (step one) or an instantaneous value (step two). The switch is based on the water concentration in moles/cc of 1x10(exp -20). The results presented here results in a correlation that gives the chemical kinetic time as two separate functions. This two time step method is used as opposed to a one step time averaged method previously developed to determine the chemical kinetic time with increased accuracy. The first time averaged step is used at the initial times for smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, initial water to fuel mass ratio, temperature, and pressure. The second instantaneous step, to be used with higher water concentrations, gives the chemical kinetic time as a function of instantaneous fuel and water mole concentration, pressure and temperature (T4). The simple correlations would then be compared to the turbulent mixing times to determine the limiting rates of the reaction. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates are used to calculate the necessary chemical kinetic times. Chemical kinetic time equations for fuel, carbon monoxide and NOx are obtained for Jet-A fuel and methane with and without water injection to water mass loadings of 2/1 water to fuel. A similar correlation was also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium concentrations of carbon monoxide and nitrogen oxide as functions of overall equivalence ratio, water to fuel mass ratio, pressure and temperature (T3). The temperature of the gas entering the turbine (T4) was also correlated as a function of the initial combustor temperature (T3), equivalence ratio, water to fuel mass ratio, and pressure

Limitations of rapid myelin water quantification using 3D bSSFP

Object: Imaging of the myelin water fraction (MWF) is conventionally performed using a multi-echo spin-echo sequence. This technique requires long acquisition times and therefore often suffers from a lack of volume coverage. In this work, the application of 3D balanced steady-state free precession (bSSFP) sequences to extract high-resolution myelin water maps is discussed. Materials and Methods: Based on a two-pool water exchange model, an approximate bSSFP signal equation is derived and applied to fit the flip angle dependence of the invivo bSSFP signal. Thereby, the MWF and signal amplitude are fitted, while a priori assumptions are made for the other parameters of the two-pool system. Results: The effects of magnetization transfer, finite RF pulses, B 0 and B 1 inhomogeneities, as well as variation of the constant fit parameters, are investigated. Acquisition and calculation of quantitative, high-resolution MWF maps from white matter of healthy volunteers based on bSSFP is feasible and averaged MWF fit results agree with literature. However, results from numerical simulations indicate a severe dependence of the derived MWF values on the constant two-pool parameters. Conclusion: The demonstrated dependence of the MWF on the two-pool parameters considerably impairs the applicability of the proposed metho

Comparative Analysis of Small-Scale Organic Rankine Cycle Systems for Solar Energy Utilisation

Small-scale organic Rankine cycle (ORC) systems driven by solar energy are compared in this paper, which aims to explore the potential of power generation for domestic utilisation. A solar thermal collector was used as the heat source for a hot water storage tank. Thermal performance was then evaluated in terms of both the conventional ORC and an ORC using thermal driven pump (TDP). It is established that the solar ORC using TDP has a superior performance to the conventional ORC under most working conditions. Results demonstrate that power output of the ORC using TDP ranges from 72 W to 82 W with the increase of evaporating temperature, which shows an improvement of up to 3.3% at a 100 °C evaporating temperature when compared with the power output of the conventional ORC. Energy and exergy efficiencies of the ORC using TDP increase from 11.3% to 12.6% and from 45.8% to 51.3% when the evaporating temperature increases from 75 °C to 100 °C. The efficiency of the ORC using TDP is improved by up to 3.27%. Additionally, the exergy destruction using TDP can be reduced in the evaporator and condenser. The highest exergy efficiency in the evaporator is 96.9%, an improvement of 62% in comparison with that of the conventional ORC, i.e., 59.9%. Thus, the small-scale solar ORC system using TDP is more promising for household application

A new framework for large strain electromechanics based on convex multi-variable strain energies: Finite Element discretisation and computational implementation

In Gil and Ortigosa (2016), Gil and Ortigosa introduced a new convex multi-variable framework for the numerical simulation of Electro Active Polymers (EAPs) in the presence of extreme deformations and electric fields. This extends the concept of polyconvexity to strain energies which depend on non-strain based variables. The consideration of the new concept of multi-variable convexity guarantees the well posedness of generalised Gibbs’ energy density functionals and, hence, opens up the possibility of a new family of mixed variational principles. The aim of this paper is to present, as an example, the Finite Element implementation of two of these mixed variational principles. These types of enhanced methodologies are known to be necessary in scenarios in which the simpler displacement-potential based formulation yields non-physical results, such as volumetric locking, bending and shear locking, pressure oscillations and electro-mechanical locking, to name but a few. Crucially, the use of interpolation spaces in which some of the unknown fields are described as piecewise discontinuous across elements can be used in order to efficiently condense these fields out. This results in mixed formulations with a computational cost comparable to that of the displacement-potential based approach, yet far more accurate. Finally, a series of very challenging numerical examples are presented in order to demonstrate the accuracy, robustness and efficiency of the proposed methodology

On the track of the brain's microstructure : myelin water imaging using quantitative MRI

Conventional quantitative magnetic resonance imaging (MRI), for example monoexponential determination of the relaxation times T1 and T2, is sensitive to the various pathologies of myelinated tissue in the brain. However, it gives relatively unspecific information about the underlying nature of the disease. A parameter that directly correlates with the integrity of the myelin sheath is the so-called myelin water fraction (MWF). Based on multi-component analysis of non-invasive quantitative MRI measurements, mapping of the MWF becomes feasible and proved to be useful for studying demyelination and remyelination processes in the course of multiple sclerosis (MS) and other myelin related pathologies. Common myelin water imaging techniques often suffer from a lack of volume coverage due to their 2D acquisition schemes. This thesis focuses on the development of new myelin water mapping procedures, especially on fast 3D MRI measurements that provide whole brain coverage. In chapter 2, an MWF mapping technique based on balanced steady-state free precession (bSSFP) sequences is introduced. An extended bSSFP signal equation, which is based on a two-pool water model describing brain tissue, is derived to determine typical multi-compartment parameters, including the MWF, of healthy subjects. Possible influences of magnetization transfer effects, infinite radiofrequency pulses and B0/B1 inhomogeneities are discussed extensively. Chapter 3 introduces a 3D acquisition scheme based on multi-gradient-echo (mGRE) pulse sequences that is applied for sampling multi-component T2* decays in the human brain of healthy volunteers and MS patients. Quantitative myelin water maps are generated based on analysis of T2* spectra. Chapter 4 discusses possible adaptations and modifications of the proposed procedure from chapter 3 when moving to higher main magnetic field strengths. The effects of B0 inhomogeneities on the data sets and possible correction methods are additionally covered in this part of the thesis. Finally, the crucial role of accurate B1 and B0 imaging and the influences on myelin water imaging are revisited in chapter 5. A solution to simultaneous mapping of B1 and B0 is presented that might help to overcome systematic error sources in MWF mapping in the future

Multiphase Neutral Interstellar Medium: Analyzing Simulation with H I 21cm Observational Data Analysis Techniques

Several different methods are regularly used to infer the properties of the neutral interstellar medium (ISM) using atomic hydrogen (H I) 21cm absorption and emission spectra. In this work, we study various techniques used for inferring ISM gas phase properties, namely the correlation between brightness temperature and optical depth $(T_B(v)$, $\tau(v))$ at each channel velocity $(v)$, and decomposition into Gaussian components, by creating mock spectra from a 3D magnetohydrodynamic simulation of a two-phase, turbulent ISM. We propose a physically motivated model to explain the $T_B(v)-\tau(v)$ distribution and relate the model parameters to properties like warm gas spin temperature and cold cloud length scales. Two methods based on Gaussian decomposition -- using only absorption spectra and both absorption and emission spectra -- are used to infer the column density distribution as a function of temperature. In observations, such analysis reveals the puzzle of large amounts (significantly higher than in simulations) of gas with temperature in the thermally unstable range of $\sim200\mathrm{\ K}$ to $\sim2000\mathrm{\ K}$ and a lack of the expected bimodal (two-phase) temperature distribution. We show that, in simulation, both methods are able to recover the true gas distribution till temperatures $\lesssim2500\mathrm{\ K}$ (and the two-phase distribution in general) reasonably well. We find our results to be robust to a range of effects such as noise, varying emission beam size, and simulation resolution. This shows that the observational inferences are unlikely to be artifacts, thus highlighting a tension between observations and simulations. We discuss possible reasons for this tension and ways to resolve it.Comment: 21 pages (including appendixes), 15 figures, 3 tables, Submitted to MNRAS, Comments are welcom