Sviluppo di algoritmi di RCP per la stima della coppia indicata erogata da un motore a combustione interna.

Combustion control is one of the key factors to obtain better performances and lower pollutant emissions for diesel, spark ignition and HCCI engines. An algorithm that allows estimating, as an example, the mean indicated torque for each cylinder, could be easily used in control strategies, in order to carry out cylinders trade-off, control the cycle to cycle variation, or detect misfires. A tool that allows evaluating the 50% of Mass Fraction Burned (MFB50), or the net Cumulative Heat Release (CHRNET), or the ROHR peak value (Rate of Heat Release), could be used to optimize spark advance or to detect knock in gasoline engines and to optimize injection pattern in diesel engines. Modern management systems are based on the control of the mean indicated torque produced by the engine: they need a real or virtual sensor in order to compare the measured value with the target one. Many studies have been performed in order to obtain an accurate and reliable over time torque estimation. The aim of this PhD activity was to develop two different algorithms: the first one is based on the instantaneous engine speed fluctuations measurement. The speed signal is picked up directly from the sensor facing the toothed wheel mounted on the engine for other control purposes. The engine speed fluctuation amplitudes depend on the combustion and on the amount of torque delivered by each cylinder. The second algorithm processes in-cylinder pressure signals in the angular domain. In this case a crankshaft encoder is not necessary, because the angular reference can be obtained using a standard sensor wheel. The results obtained with these two methodologies are compared in order to evaluate which one is suitable for on board applications, depending on the accuracy required.

LES of flame acceleration and DDT in hydrogen-air mixture using artificially thickened flame approach and detailed chemical kinetics

A large eddy simulation is performed to study the deflagration to detonation transition phenomenon in an obstructed channel containing premixed stoichiometric hydrogen–air mixture. Two-dimensional filtered reactive Navier–Stokes equations are solved utilizing the artificially thickened flame approach (ATF) for modeling sub-grid scale combustion. To include the effect of induction time, a 27-step detailed mechanism is utilized along with an in situ adaptive tabulation (ISAT) method to reduce the computational cost due to the detailed chemistry. The results show that in the slow flame propagation regime, the flame–vortex interaction and the resulting flame folding and wrinkling are the main mechanisms for the increase of the flame surface and consequently acceleration of the flame. Furthermore, at high speed, the major mechanisms responsible for flame propagation are repeated reflected shock–flame interactions and the resulting baroclinic vorticity. These interactions intensify the rate of heat release and maintain the turbulence and flame speed at high level. During the flame acceleration, it is seen that the turbulent flame enters the ‘thickened reaction zones’ regime. Therefore, it is necessary to utilize the chemistry based combustion model with detailed chemical kinetics to properly capture the salient features of the fast deflagration propagation.

Vortices, sound and flames - a damaging combination

The interaction between vortices, sound and combustion can lead to self-excited oscillations of such large amplitudes that structural damage is done. These occur because any small unsteadiness in the rate of combustion is a source of sound, generating pressure and velocity fluctuations. However, the velocity fluctuations perturb the flame, thereby altering the instantaneous rate of heat release. Instability is then possible because while acoustic waves perturb the combustion, the unsteady combustion generates yet more sound! Combustion oscillations can occur in afterburners and at idle in conventional aeroengine combustors. Lean premixed, prevapourized technology has tremendous potential to reduce NOx emissions, but is proving highly susceptible to self-excited oscillations. An overview of the physics of the interaction between vortices, sound and flames is presented, and illustrated by examples of instability in generic premixed ducted flames and in aeroengine combustors. The potential for both passive and active control is discussed.

Instability in lean premixed combustors

A generalized acoustic equation is used to identify the mechanisms driving combustion instability. The relationship between the unsteady rate of heat release and the flow is found to influence significantly the frequency of oscillation. A kinematic flame model is reviewed and used to describe the unsteady combustion in a premixed ducted flame and in a typical lean premixed industrial gas turbine. Comparison is made between theory and experiment. | A generalized acoustic equation is used to identify the mechanisms driving combustion instability. The relationship between the unsteady rate of heat release and the flow is found to influence significantly the frequency of oscillation. A kinematic flame model is reviewed and used to describe the unsteady combustion in a premixed ducted flame and in a typical lean premixed industrial gas turbine. Comparison is made between theory and experiment.

Effects of compressibility and heat release on entrainment processes in mixing layers

Characteristics of the process of entrainment in plane mixing layers, and the changes with compressibility and heat release, were studied using temporal DNS with simultaneous fluid packet tracking. Convective Mach numbers of the simulations are 0.15, 0.7 and 1.1. The Reynolds number is quite high (between 11 000 and 37 000 based on layer width and velocity difference), and is above the mixing transition. The study agrees with recent findings in round jets: first, engulfed fluid volume and its growth rate are both very small compared with the volume of the turbulent region and its growth rate, respectively. Secondly, most often, the process occurs close to the turbulent-nonturbulent boundaries. A new finding is that both compressibility and heat release retard the entrainment process so that it takes an O(1) time for vorticity or scalar levels to grow even after growth has been initiated. This delay is manifested as the fall in mixing layer growth rates as compressibility and heat release levels increase.

The nonlinear heat release response of stratified lean-premixed flames to acoustic velocity oscillations

This paper analyzes the forced response of swirl-stabilized lean-premixed flames to high-amplitude acoustic forcing in a laboratory-scale stratified burner operated with CH4 and air at atmospheric pressure. The double-swirler, double-channel annular burner was specially designed to generate high-amplitude acoustic velocity oscillations and a radial equivalence ratio gradient at the inlet of the combustion chamber. Temporal oscillations of equivalence ratio along the axial direction are dissipated over a long distance, and therefore the effects of time-varying fuel/air ratio on the response are not considered in the present investigation. Simultaneous measurements of inlet velocity and heat release rate oscillations were made using a constant temperature anemometer and photomultiplier tubes with narrow-band OH*/CH* interference filters. Time-averaged and phase-synchronized CH* chemiluminescence intensities were measured using an intensified CCD camera. The measurements show that flame stabilization mechanisms vary depending on equivalence ratio gradients for a constant global equivalence ratio (φg=0.60). Under uniformly premixed conditions, an enveloped M-shaped flame is observed. In contrast, under stratified conditions, a dihedral V-flame and a toroidal detached flame develop in the outer stream and inner stream fuel enrichment cases, respectively. The modification of the stabilization mechanism has a significant impact on the nonlinear response of stratified flames to high-amplitude acoustic forcing (u'/U∼0.45 and f=60, 160Hz). Outer stream enrichment tends to improve the flame's stiffness with respect to incident acoustic/vortical disturbances, whereas inner stream stratification tends to enhance the nonlinear flame dynamics, as manifested by the complex interaction between the swirl flame and large-scale coherent vortices with different length scales and shedding points. It was found that the behavior of the measured flame describing functions (FDF), which depend on radial fuel stratification, are well correlated with previous measurements of the intensity of self-excited combustion instabilities in the stratified swirl burner. The results presented in this paper provide insight into the impact of nonuniform reactant stoichiometry on combustion instabilities, its effect on flame location and the interaction with unsteady flow structures. © 2011 The Combustion Institute.

Disruption of endogenous tidal rhythms of larval release linked to food supply and heat stress in an intertidal barnacle

The timing of larval release may greatly affect the survivorship and distribution of pelagic stages and reveal important aspects of life history tactics in marine invertebrates. Endogenous rhythms of breeding individuals and populations are valuable indicators of selected strategies because they are free of the neutral effect of stochastic environmental variation. The high-shore intertidal barnacle Chthamalus bisinuatus exhibits endogenous tidal and tidal amplitude rhythms in a way that larval release would more likely occur during fortnightly neap periods at high tide. Such timing would minimize larval loss due to stranding and promote larval retention close to shore. This fully explains temporal patterns in populations facing the open sea and inhabiting eutrophic areas. However, rhythmic activity breaks down to an irregular pattern in a population within the São Sebastião Channel subjected to large variation of food supply around a mesotrophic average. Peaks of chl a concentration precede release events by 6 d, suggesting resource limitation for egg production within the channel. Also, extreme daily temperatures imposing mortality risk correlate to release rate just 1 d ahead, suggesting a terminal reproductive strategy. Oceanographic conditions apparently dictate whether barnacles follow a rhythmic trend of larval release supported by endogenous timing or, alternatively, respond to the stochastic variation of key environmental factors, resulting in an erratic temporal pattern.

Technical report on intensification of turbulence by chemical heat release. Contract no. DA-04-495-ORD-3567(Z) Hypervelocity Range Research Program, a part of project "Defender."

"This research was supported by the Advanced Research Projects Agency, Department of defense, and was monitored by the U.S. Army Missile Command, Redstone Arsenal, Alabama."

Mathematical modelling of controlled drug release from polymer micro-spheres : incorporating the effects of swelling, diffusion and dissolution via moving boundary problems

Controlled drug delivery is a key topic in modern pharmacotherapy, where controlled drug delivery devices are required to prolong the period of release, maintain a constant release rate, or release the drug with a predetermined release profile. In the pharmaceutical industry, the development process of a controlled drug delivery device may be facilitated enormously by the mathematical modelling of drug release mechanisms, directly decreasing the number of necessary experiments. Such mathematical modelling is difficult because several mechanisms are involved during the drug release process. The main drug release mechanisms of a controlled release device are based on the device’s physiochemical properties, and include diffusion, swelling and erosion. In this thesis, four controlled drug delivery models are investigated. These four models selectively involve the solvent penetration into the polymeric device, the swelling of the polymer, the polymer erosion and the drug diffusion out of the device but all share two common key features. The first is that the solvent penetration into the polymer causes the transition of the polymer from a glassy state into a rubbery state. The interface between the two states of the polymer is modelled as a moving boundary and the speed of this interface is governed by a kinetic law. The second feature is that drug diffusion only happens in the rubbery region of the polymer, with a nonlinear diffusion coefficient which is dependent on the concentration of solvent. These models are analysed by using both formal asymptotics and numerical computation, where front-fixing methods and the method of lines with finite difference approximations are used to solve these models numerically. This numerical scheme is conservative, accurate and easily implemented to the moving boundary problems and is thoroughly explained in Section 3.2. From the small time asymptotic analysis in Sections 5.3.1, 6.3.1 and 7.2.1, these models exhibit the non-Fickian behaviour referred to as Case II diffusion, and an initial constant rate of drug release which is appealing to the pharmaceutical industry because this indicates zeroorder release. The numerical results of the models qualitatively confirms the experimental behaviour identified in the literature. The knowledge obtained from investigating these models can help to develop more complex multi-layered drug delivery devices in order to achieve sophisticated drug release profiles. A multi-layer matrix tablet, which consists of a number of polymer layers designed to provide sustainable and constant drug release or bimodal drug release, is also discussed in this research. The moving boundary problem describing the solvent penetration into the polymer also arises in melting and freezing problems which have been modelled as the classical onephase Stefan problem. The classical one-phase Stefan problem has unrealistic singularities existed in the problem at the complete melting time. Hence we investigate the effect of including the kinetic undercooling to the melting problem and this problem is called the one-phase Stefan problem with kinetic undercooling. Interestingly we discover the unrealistic singularities existed in the classical one-phase Stefan problem at the complete melting time are regularised and also find out the small time behaviour of the one-phase Stefan problem with kinetic undercooling is different to the classical one-phase Stefan problem from the small time asymptotic analysis in Section 3.3. In the case of melting very small particles, it is known that surface tension effects are important. The effect of including the surface tension to the melting problem for nanoparticles (no kinetic undercooling) has been investigated in the past, however the one-phase Stefan problem with surface tension exhibits finite-time blow-up. Therefore we investigate the effect of including both the surface tension and kinetic undercooling to the melting problem for nanoparticles and find out the the solution continues to exist until complete melting. The investigation of including kinetic undercooling and surface tension to the melting problems reveals more insight into the regularisations of unphysical singularities in the classical one-phase Stefan problem. This investigation gives a better understanding of melting a particle, and contributes to the current body of knowledge related to melting and freezing due to heat conduction.

Measurement of Heat Transfer Rate on Backward-Facing Steps at Hypersonic Mach Number

Two backward-facing models with step heights of 2 and 3 mm are used to measure the convective surface heat transfer rates by using platinum thin-film gauges, deposited on Macor inserts. Heat transfer rates have been theoretically calculated along the flat plate portion of a model using the Eckert reference temperature method. The experimentally determined surface heat transfer rate distributions are compared with theoretical and numerical estimations. Experimental heat flux distribution over a flat plate model showed good agreement with the reference temperature method at stagnation enthalpy range of 0.8-2 MJ/kg. Theoretical analysis has been used for downstream of a backward-facing step using Gai's nondimensional analysis. It has been found from the present study that approximately 10 and 8 step heights are required for the flow to reattach for 2 and 3 mm step height backward-facing step models, respectively, at a nominal Mach number of 7.6.

A new model for the energy release rate of fibre/matrix interfacial fracture

On the basis of the well-known shear-lag analysis of fibre/matrix interface stresses and the assumption of identical axial strains in the fibre and matrix, a new model for predicting the energy release rate of interfacial fracture of the fibre pull-out test model is attempted. The expressions for stresses in the fibre, matrix and interface are derived. The formula for interfacial debonding energy release rate is given. Numerical calculations are conducted and the results obtained are compared with those of the existing models.

Investigation on Depth of Heat Affected Zone of Discharge Spot by High Repetitive Rate YAG Laser-Induced Discharge Texturing

It is known that the press formability and the elongation of laser textured sheet are improved, and the service life of textured roll is longer than that of the un-textured roll due to hardening of the treated surface. One of the goals to develop high repetitive rate YAG laser-induced discharge texturing (LIDT) is to get deeper hardening zone. By observing and measuring cross-section of LIDT spots in different discharge conditions, it is found that the single-crater, which is formed by the discharge conditions of anode, which is covered by an oil film and with rectangular current waveform, has the most depth of heat affected zone (HAZ) comparing with other crater shapes when discharge energy is the same. The depth of HAZ is mainly depends on pulse duration when the discharge spot is single-crater. The results are analyzed.

Debonding and cracking energy release rate of the fiber/matrix interface

An investigation of fiber/matrix interfacial fracture energy is presented in this paper. Several existing theoretical expressions for the fracture energy of interfacial debonding are reviewed. For the single-fiber/matrix debonding and pull-out experimental model, a study is carried out on the effect of interfacial residual compressive stress and friction on interface cracking energy release rate.