Wave breaking on turbulent energy budget in the ocean surface mixed layer

As an important physical process at the air-sea interface, wave movement and breaking have a significant effect on the ocean surface mixed layer (OSML). When breaking waves occur at the ocean surface, turbulent kinetic energy (TKE) is input downwards, and a sublayer is formed near the surface and turbulence vertical mixing is intensively enhanced. A one-dimensional ocean model including the Mellor-Yamada level 2.5 turbulence closure equations was employed in our research on variations in turbulent energy budget within OSML. The influence of wave breaking could be introduced into the model by modifying an existing surface boundary condition of the TKE equation and specifying its input. The vertical diffusion and dissipation of TKE were effectively enhanced in the sublayer when wave breaking was considered. Turbulent energy dissipated in the sublayer was about 92.0% of the total depth-integrated dissipated TKE, which is twice higher than that of non-wave breaking. The shear production of TKE decreased by 3.5% because the mean flow fields tended to be uniform due to wave-enhanced turbulent mixing. As a result, a new local equilibrium between diffusion and dissipation of TKE was reached in the wave-enhanced layer. Below the sublayer, the local equilibrium between shear production and dissipation of TKE agreed with the conclusion drawn from the classical law-of-the-wall (Craig and Banner, 1994).

Traveling wave front for the Fisher equation on an infinite band region

Instead of discussing the existence of a one-dimensional traveling wave front solution which connects two constant steady states, the present work deals with the case connecting a constant and a nonhomogeneous steady state on an infinite band region. The corresponding model is the well-known Fisher equation with variational coefficient and Dirichlet boundary condition. (c) 2006 Elsevier Ltd. All rights reserved.

软土亚塑性模型与一维非线性固结理论研究

With the great development of Tianjing New Coastal District economy, people need more land to build and live. Land subsidence, which is caused by its special engineering geological conditions, has restricted the further development in the district. Soft soil consolidation is main factor of land subsidence ;thus , on the basis of consolidation theory, the paper make further study on soft soils one-dimension nonlinear consolidation which contains two parts:(1) the nonlinear consolidation of permeability coefficient and compressibility coefficient changing with time and depth, which means real one-dimension nonlinear consolidation;(2) the non-homogeneous consolidation of permeability coefficient and compressibility coefficient only changing with depth. Firstly, nonlinear characteristics of soft soils are elaborated. Hypoplastic theory is introduced to establish a modified soft soils nonlinear constitutive model; the nonlinear governing equation of compressibility coefficient is built, and the nonlinear characteristics of compressibility coefficient are analyzed. Secondly, Considering Load Fluctuation and soil thickness changing ,the consolidation characteristics of single layer is discussed in the paper; meanwhile, on the basis of the Davis and Raymond’s hypothesis and single layer nonlinear consolidation equation, the doubled-layer one-dimension nonlinear consolidation equation is also derived. The solution of the equation is obtained by analytical method, and the consolidation characteristics of doubled-layer soft soil nonlinear theory is also analyzed. Finally, based on assumption that permeability coefficient and compressibility coefficient is varying along depth, single layer soil one-dimension non-homogeneous consolidation differential equation is derived; and the approximate solution is obtained. Furthermore, the single layer non-homogeneous consolidation is extended to double layer non-homogeneous consolidation theory. By using parabolic differential scheme, the matrix equation is established; and the solution of the matrix equation is obtained by chase method. Consolidation characteristics of soil soft single (double) layer non-homogeneous consolidation theory and Terzaghi’s theory are also discussed.

用宽频带远震体波波形研究中国及邻近地区地壳构造及地质意义初探

The theory and approach of the broadband teleseismic body waveform inversion are expatiated in this paper, and the defining the crust structure's methods are developed. Based on the teleseismic P-wave data, the theoretic image of the P-wave radical component is calculated via the convolution of the teleseismic P-wave vertical component and the transform function, and thereby a P-wavefrom inversion method is built. The applied results show the approach effective, stable and its resolution high. The exact and reliable teleseismic P waveforms recorded by CDSN and IRIS and its geodynamics are utilized to obtain China and its vicinage lithospheric transfer functions, this region ithospheric structure is inverted through the inversion of reliable transfer functions, the new knowledge about the deep structure of China and its vicinage is obtained, and the reliable seismological evidence is provided to reveal the geodynamic evolution processes and set up the continental collisional theory. The major studies are as follows: Two important methods to study crustal and upper mantle structure -- body wave travel-time inversion and waveform modeling are reviewed systematically. Based on ray theory, travel-time inversion is characterized by simplicity, crustal and upper mantle velocity model can be obtained by using 1-D travel-time inversion preliminary, which introduces the reference model for studying focal location, focal mechanism, and fine structure of crustal and upper mantle. The large-scale lateral inhomogeneity of crustal and upper mantle can be obtained by three-dimensional t ravel-time seismic tomography. Based on elastic dynamics, through the fitting between theoretical seismogram and observed seismogram, waveform modeling can interpret the detail waveform and further uncover one-dimensional fine structure and lateral variation of crustal and upper mantle, especially the media characteristics of singular zones of ray. Whatever travel-time inversion and waveform modeling is supposed under certain approximate conditions, with respective advantages and disadvantages, and provide convincing structure information for elucidating physical and chemical features and geodynamic processes of crustal and upper mantle. Because the direct wave, surface wave, and refraction wave have lower resolution in investigating seismic velocity transitional zone, which is inadequate to study seismic discontinuities. On the contrary, both the converse and reflected wave, which sample the discontinuities directly, must be carefully picked up from seismogram to constrain the velocity transitional zones. Not only can the converse wave and reflected wave study the crustal structure, but also investigate the upper mantle discontinuities. There are a number of global and regional seismic discontinuities in the crustal and upper mantle, which plays a significant role in understanding physical and chemical properties and geodynamic processes of crustal and upper mantle. The broadband teleseismic P waveform inversion is studied particularly. The teleseismic P waveforms contain a lot of information related to source time function, near-source structure, propagation effect through the mantle, receiver structure, and instrument response, receiver function is isolated form teleseismic P waveform through the vector rotation of horizontal components into ray direction and the deconvolution of vertical component from the radial and tangential components of ground motion, the resulting time series is dominated by local receiver structure effect, and is hardly irrelevant to source and deep mantle effects. Receiver function is horizontal response, which eliminate multiple P wave reflection and retain direct wave and P-S converted waves, and is sensitive to the vertical variation of S wave velocity. Velocity structure beneath a seismic station has different response to radial and vertical component of an accident teleseismic P wave. To avoid the limits caused by a simplified assumption on the vertical response, the receiver function method is mended. In the frequency domain, the transfer function is showed by the ratio of radical response and vertical response of the media to P wave. In the time domain, the radial synthetic waveform can be obtained by the convolution of the transfer function with the vertical wave. In order to overcome the numerical instability, generalized reflection and transmission coefficient matrix method is applied to calculate the synthetic waveform so that all multi-reflection and phase conversion response can be included. A new inversion method, VFSA-LM method, is used in this study, which successfully combines very fast simulated annealing method (VFSA) with damped least square inversion method (LM). Synthetic waveform inversion test confirms its effectiveness and efficiency. Broadband teleseismic P waveform inversion is applied in lithospheric velocity study of China and its vicinage. According to the data of high quality CDSN and IRIS, we obtained an outline map showing the distribution of Asian continental crustal thickness. Based on these results gained, the features of distribution of the crustal thickness and outline of crustal structure under the Asian continent have been analyzed and studied. Finally, this paper advances the principal characteristics of the Asian continental crust. There exist four vast areas of relatively minor variations in the crustal thickness, namely, northern, eastern southern and central areas of Asian crust. As a byproduct, the earthquake location is discussed, Which is a basic issue in seismology. Because of the strong trade-off between the assumed initial time and focal depth and the nonlinear of the inversion problems, this issue is not settled at all. Aimed at the problem, a new earthquake location method named SAMS method is presented, In which, the objective function is the absolute value of the remnants of travel times together with the arrival times and use the Fast Simulated Annealing method is used to inverse. Applied in the Chi-Chi event relocation of Taiwan occurred on Sep 21, 2000, the results show that the SAMS method not only can reduce the effects of the trade-off between the initial time and focal depth, but can get better stability and resolving power. At the end of the paper, the inverse Q filtering method for compensating attenuation and frequency dispersion used in the seismic section of depth domain is discussed. According to the forward and inverse results of synthesized seismic records, our Q filtrating operator of the depth domain is consistent with the seismic laws in the absorbing media, which not only considers the effect of the media absorbing of the waves, but also fits the deformation laws, namely the frequency dispersion of the body wave. Two post stacked profiles about 60KM, a neritic area of China processed, the result shows that after the forward Q filtering of the depth domain, the wide of the wavelet of the middle and deep layers is compressed, the resolution and signal noise ratio are enhanced, and the primary sharp and energy distribution of the profile are retained.

地幔对流理论及其在解释地球物理观测的应用

With the improving of mantle convection theory, the developing of computing method and increasing of the measurement data, we can numerically simulate more clearly about the effects on some geophysical observed phenomenons such as the global heat flow and global lithospheric stress field in the Earth's surface caused by mantle convection, which is the primary mechanism for the transport of heat from the Earth's deep interior to its surface and the underlying force mechanism of dynamics in the Earth.Chapter 1 reviews the historical background and present research state of mantle convection theory.In Chapter 2, the basic conception of thermal convection and the basic theory about mantle flow.The effects on generation and distribution of global lithospheric stres s field induced by mantle flow are the subject of Chapter 3. Mantle convection causes normal stress and tangential stresses at the bottom of the lithosphere, and then the sublithospheric stress field induces the lithospheric deformation as sixrface force and results in the stress field within the lithosphere. The simulation shows that the agreement between predictions and observations is good in most regions. Most of subduction zones and continental collisions are under compressive. While ocean ridges, such as the east Pacific ridge, the Atlantic ridge and the east African rift valley, are under tensile. And most of the hotspots preferentially occur in regions where calculated stress is tensile. The calculated directions of the most compressive principal horizontal stress are largely in accord with that of the observation except for some regions such as the NW-Pacifie subduction zone and Qinghai-Tibet Plateau, in which the directions of the most compressive principal horizontal stress are different. It shows that the mantel flow plays an important role in causing or affecting the large-scale stress field within the lithosphere.The global heat flow simulation based on a kinematic model of mantle convection is given in Chapter 4. Mantle convection velocities are calculated based on the internal loading theory at first, the velocity field is used as the input to solve the thermal problem. Results show that calculated depth derivatives of the near surface temperature are closely correlated to the observed surface heat flow pattern. Higher heat flow values around midocean ridge systems can be reproduced very well. The predicted average temperature as a function of function of depth reveals that there are two thermal boundary layers, one is close to the surface and another is close to the core-mantle boundary, the rest of the mantle is nearly isothermal. Although, in most of the mantle, advection dominates the heat transfer, the conductive heat transfer is still locally important in the boundary layers and plays an important role for the surface heat flow pattern. The existence of surface plates is responsible for the long wavelength surface heat flow pattern.In Chapter 5, the effects on present-day crustal movement in the China Mainland resulted from the mantle convection are introduced. Using a dynamic method, we present a quantitative model for the present-day crustal movement in China. We consider not only the effect of the India-Eurasia collision, the gravitational potential energy difference of the Tibet Plateau, but also the contribution of the shear traction on the bottom of the lithosphere induced by the global mantle convection. The comparison between our results and the velocity field obtained from the GPS observation shows that our model satisfactorily reproduces the general picture of crustal deformation in China. Numerical modeling results reveal that the stress field on the base of the lithosphere induced by the mantle flow is probably a considerable factor that causes the movement and deformation of the lithosphere in continental China with its eflfcet focuing on the Eastern China A numerical research on the small-scale convection with variable viscosity in the upper mantle is introduced in Chapter 6. Based on a two-dimensional model, small-scale convection in the mantle-lithosphere system with variable viscosity is researched by using of finite element method. Variation of viscosity in exponential form with temperature is considered in this paper The results show that if viscosity is strongly temperature-dependent, the upper part of the system does not take a share in the convection and a stagnant lid, which is identified as lithosphere, is formed on the top of system because of low temperature and high viscosity. The calculated surface heat flow, topography and gravity anomaly are associated well with the convection pattern, namely, the regions with high heat flow and uplift correspond to the upwelling flow, and vice versa.In Chapter 7, we give a brief of future research subject: The inversion of lateral density heterogeneity in the mantle by minimizing the viscous dissipation.

Quantitative determination of compounds in tobacco essential oils by comprehensive two-dimensional gas chromatography coupled to time-of-flight mass spectrometry

A quantitative analysis of the individual compounds in tobacco essential oils is performed by comprehensive two-dimensional gas chromatography (GC x GC) combined with flame ionization detector (FID). A time-of-flight mass spectrometer (TOF/MS) was coupled to GC x GC for the identification of the resolved peaks. The response of a flame ionization detector to different compound classes was calibrated using multiple internal standards. In total, 172 compounds were identified with good match and 61 compounds with high probability value were reliably quantified. For comparative purposes, the essential oil sample was also quantified by one-dimensional gas chromatography-mass spectrometry (GC/MS) with multiple internal standards method. The results showed that there was close agreement between the two analysis methods when the peak purity and match quality in one-dimensional GC/MS are high enough. (c) 2005 Elsevier B.V. All rights reserved.

Multi-Area Visuotopic Map Complexes in Macaque Striate and Extra-striate Cortex

We propose that a simple, closed-form mathematical expression--the Wedge-Dipole mapping--provides a concise approximation to the full-field, two-dimensional topographic structure of macaque V1, V2, and V3. A single map function, which we term a map complex, acts as a simultaneous descriptor of all three areas. Quantitative estimation of the Wedge-Dipole parameters is provided via 2DG data of central-field V1 topography and a publicly available data set of full-field macaque V1 and V2 topography. Good quantitative agreement is obtained between the data and the model presented here. The increasing importance of fMRI-based brain imaging motivates the development of more sophisticated two-dimensional models of cortical visuotopy, in contrast to the one-dimensional approximations that have been in common use. One reason is that topography has traditionally supplied an important aspect of "ground truth", or validation, for brain imaging, suggesting that further development of high-resolution fMRI will be facilitated by this data analysis. In addition, several important insights into the nature of cortical topography follows from this work. The presence of anisotropy in cortical magnification factor is shown to follow mathematically from the shared boundary conditions at the V1-V2 and V2-V3 borders, and therefore may not causally follow from the existence of columnar systems in these areas, as is widely assumed. An application of the Wedge-Dipole model to localizing aspects of visual processing to specific cortical areas--extending previous work in correlating V1 cortical magnification factor to retinal anatomy or visual psychophysics data--is briefly discussed.

Normal and Amnesis Learning, Recognition, and Memory by a Neural Model of Cortico-Hippocampal Interactions

The processes by which humans and other primates learn to recognize objects have been the subject of many models. Processes such as learning, categorization, attention, memory search, expectation, and novelty detection work together at different stages to realize object recognition. In this article, Gail Carpenter and Stephen Grossberg describe one such model class (Adaptive Resonance Theory, ART) and discuss how its structure and function might relate to known neurological learning and memory processes, such as how inferotemporal cortex can recognize both specialized and abstract information, and how medial temporal amnesia may be caused by lesions in the hippocampal formation. The model also suggests how hippocampal and inferotemporal processing may be linked during recognition learning.

Two closely coupled lasers, dynamics & applications

The dynamics of two mutually coupled identical single-mode semi-conductor lasers are theoretically investigated. For small separation and large coupling between the lasers, symmetry-broken one-colour states are shown to be stable. In this case the light output of the lasers have significantly different intensities whilst at the same time the lasers are locked to a single common frequency. For intermediate coupling we observe stable two-colour states, where both single-mode lasers lase simultaneously at two optical frequencies which are separated by up to 150 GHz. For low coupling but possibly large separation, the frequency of the relaxation oscillations of the freerunning lasers defines the dynamics. Chaotic and quasi-periodic states are identified and shown to be stable. For weak coupling undamped relaxation oscillations dominate where each laser is locked to three or more odd number of colours spaced by the relaxation oscillation frequency. It is shown that the instabilities that lead to these states are directly connected to the two colour mechanism where the change in the number of optical colours due to a change in the plane of oscillation. At initial coupling, in-phase and anti-phase one colour states are shown to emerge from “on” uncoupled lasers using a perturbation method. Similarly symmetry-broken one-colour states come from considering one free-running laser initially “on” and the other laser initially “off”. The mechanism that leads to a bi-stability between in-phase and anti-phase one-colour states is understood. Due to an equivariant phase space symmetry of being able to exchange the identical lasers, a symmetric and symmetry-broken variant of all states mentioned above exists and is shown to be stable. Using a five dimensional model we identify the bifurcation structure which is responsible for the appearance of symmetric and symmetry-broken one-colour, symmetric and symmetry-broken two-colour, symmetric and symmetry-broken undamped relaxation oscillations, symmetric and symmetry-broken quasi-periodic, and symmetric and symmetry-broken chaotic states. As symmetry-broken states always exist in pairs, they naturally give rise to bi-stability. Several of these states show multistabilities between symmetric and symmetry-broken variants and among states. Three memory elements on the basis of bi-stabilities in one and two colour states for two coupled single-mode lasers are proposed. The switching performance of selected designs of optical memory elements is studied numerically.

Analytic expressions for the constitutive parameters of magnetoelectric metamaterials

Electromagnetic metamaterials are artificially structured media typically composed of arrays of resonant electromagnetic circuits, the dimension and spacing of which are considerably smaller than the free-space wavelengths of operation. The constitutive parameters for metamaterials, which can be obtained using full-wave simulations in conjunction with numerical retrieval algorithms, exhibit artifacts related to the finite size of the metamaterial cell relative to the wavelength. Liu showed that the complicated, frequency-dependent forms of the constitutive parameters can be described by a set of relatively simple analytical expressions. These expressions provide useful insight and can serve as the basis for more intelligent interpolation or optimization schemes. Here, we show that the same analytical expressions can be obtained using a transfer-matrix formalism applied to a one-dimensional periodic array of thin, resonant, dielectric, or magnetic sheets. The transfer-matrix formalism breaks down, however, when both electric and magnetic responses are present in the same unit cell, as it neglects the magnetoelectric coupling between unit cells. We show that an alternative analytical approach based on the same physical model must be applied for such structures. Furthermore, in addition to the intercell coupling, electric and magnetic resonators within a unit cell may also exhibit magnetoelectric coupling. For such cells, we find an analytical expression for the effective index, which displays markedly characteristic dispersion features that depend on the strength of the coupling coefficient. We illustrate the applicability of the derived expressions by comparing to full-wave simulations on magnetoelectric unit cells. We conclude that the design of metamaterials with tailored simultaneous electric and magnetic response-such as negative index materials-will generally be complicated by potentially unwanted magnetoelectric coupling. © 2010 The American Physical Society.

Model sensitivity and uncertainty analysis using roadside air quality measurements

Most of the air quality modelling work has been so far oriented towards deterministic simulations of ambient pollutant concentrations. This traditional approach, which is based on the use of one selected model and one data set of discrete input values, does not reflect the uncertainties due to errors in model formulation and input data. Given the complexities of urban environments and the inherent limitations of mathematical modelling, it is unlikely that a single model based on routinely available meteorological and emission data will give satisfactory short-term predictions. In this study, different methods involving the use of more than one dispersion model, in association with different emission simulation methodologies and meteorological data sets, were explored for predicting best CO and benzene estimates, and related confidence bounds. The different approaches were tested using experimental data obtained during intensive monitoring campaigns in busy street canyons in Paris, France. Three relative simple dispersion models (STREET, OSPM and AEOLIUS) that are likely to be used for regulatory purposes were selected for this application. A sensitivity analysis was conducted in order to identify internal model parameters that might significantly affect results. Finally, a probabilistic methodology for assessing urban air quality was proposed.

Study of Anisotropic Conductive Adhesive Joint Behaviour under 3-Point Bending

Flip chip interconnections using anisotropic conductive film (ACF) are now a very attractive technique for electronic packaging assembly. Although ACF is environmentally friendly, many factors may influence the reliability of the final ACF joint. External mechanical loading is one of these factors. Finite element analysis (FEA) was carried out to understand the effect of mechanical loading on the ACF joint. A 3-dimensional model of adhesively bonded flip chip assembly was built and simulations were performed for the 3-point bending test. The results show that the stress at its highest value at the corners, where the chip and ACF were connected together. The ACF thickness was increased at these corner regions. It was found that higher mechanical loading results in higher stress that causes a greater gap between the chip and the substrate at the corner position. Experimental work was also carried out to study the electrical reliability of the ACF joint with the applied bending load. As per the prediction from FEA, it was found that at first the corner joint failed. Successive open joints from the corner towards the middle were also noticed with the increase of the applied load.

Structure and denaturation of 4-chlorobenzoyl coenzyme A dehalogenase from Arthrobacter sp. strain TM-1

The secondary structure of the trimeric protein 4-chlorobenzoyl coenzyme A dehalogenase from Arthrobacter sp. strain TM-1, the second of three enzymes involved in the dechlorination of 4-chlorobenzoate to form 4-hydroxybenzoate, has been examined. E(mM) for the enzyme was 12.59. Analysis by circular dichroism spectrometry in the far uv indicated that 4-chlorobenzoyl coenzyme A dehalogenase was composed mostly of alpha-helix (56%) with lesser amounts of random coil (21%), beta-turn (13%) and beta-sheet (9%). These data are in close agreement with a computational prediction of secondary structure from the primary amino acid sequence, which indicated 55.8% alpha-helix, 33.7% random coil and 10.5% beta-sheet; the enzyme is, therefore, similar to the 4-chlorobenzoyl coenzyme A dehalogenase from Pseudomonas sp. CBS-3. The three-dimensional structure, including that of the presumed active site, predicted by computational analysis, is also closely similar to that of the Pseudomonas dehalogenase. Study of the stability and physicochemical properties revealed that at room temperature, the enzyme was stable for 24 h but was completely inactivated by heating to 60 degrees C for 5 min; thereafter by cooling at 1 degrees C min(-1) to 45 degrees C, 20.6% of the activity could be recovered. Mildly acidic (pH 5.2) or alkaline (pH 10.1) conditions caused complete inactivation, but activity was fully recovered on returning the enzyme to pH 7.4. Circular dichroism studies also indicated that secondary structure was little altered by heating to 60 degrees C, or by changing the pH from 7.4 to 6.0 or 9.2. Complete, irreversible destruction of, and maximal decrease in the fluorescence yield of the protein at 330-350 nm were brought about by 4.5 M urea or 1.1 M guanidinium chloride. Evidence was obtained to support the hypothetical three-dimensional model, that residues W140 and W167 are buried in a non-polar environment, whereas W182 appears at or close to the surface of the protein. At least one of the enzymes of the dehalogenase system (the combined 4-chlorobenzoate:CoA ligase, the dehalogenase and 4-hydroxybenzoyl coenzyme A thioesterase) appears to be capable of association with the cell membrane.

A Note On Tidal Asymmetry In The Severn Estuary

A one-dimensional, non-linear numerical model is used to investigate the tidally averaged frictional stress and set-up of water level due to tidal asymmetry in the Severn Estuary; these quantities depend on the overtides in the region. A linearized model of the overtides is applied to calculations of the M4 currents in order to delineate the mechanisms responsible for their generation. The relative importance of individual non-linear mechanisms to the tidally averaged stress and set-up is determined; these mechanisms are interactions between tidal flow and changes in depth or breadth over a cross-section, frictional interaction between the tidal flow and Stokes drift, interaction between the tidal fluctuations in water depth and frictional retardation and non-linear advection.

Drivetrain Effects on Small Engine Performance

Presented is a study that expands the body of knowledge on the effect of in-cycle speed fluctuations on performance of small engines. It uses the methods developed previously by Callahan, et al. (1) to examine a variety of two-stroke engines and one four-stroke engine. The two-stroke engines were: a high performance single-cylinder, a low performance single-cylinder, a high performance multi-cylinder, and a medium performance multi-cylinder. The four-stroke engine was a high performance single-cylinder unit. Each engine was modeled in Virtual Engines, which is a fully detailed one-dimensional thermodynamic engine simulator. Measured or predicted in-cycle speed data were input into the engine models. Predicted performance changes due to drivetrain effects are shown in each case, and conclusions are drawn from those results. The simulations for the high performance single-cylinder two-stroke engine predicted significant in-cycle crankshaft speed fluctuation amplitudes and significant changes in performance when the fluctuations were input into the engine model. This was validated experimentally on a firing test engine based on a Yamaha YZ250. The four-stroke engine showed significant changes in predicted performance compared to the prediction with zero speed fluctuation assumed in the model. Measured speed fluctuations from a firing Yamaha YZ400F engine were applied to the simulation in addition to data from a simple free mass model. Both methods predicted similar fluctuation profiles and changes in performance. It is shown that the gear reduction between the crankshaft and clutch allowed for this similar behavior. The multi-cylinder, high performance two-stroke engine also showed significant changes in performance, in this case depending on the firing configuration. The low output two-stroke engine simulation showed only a negligible change in performance in spite of high amplitude speed fluctuations. This was due to its flat torque versus speed characteristic. The medium performance multi-cylinder two-stroke engine also showed only a negligible change in performance, in this case due to a relatively high inertia rotating assembly and multiple cylinder firing events within the revolution. These smoothed the net torque pulsations and reduced the amplitude of the speed fluctuation itself.