A detailed chemistry multi-cycle simulation of a gasoline fueled HCCI engine operated with NVO

A previously developed Stochastic Reactor Model (SRM) is used to simulate combustion in a four cylinder in-line four-stroke naturally aspirated direct injection Spark Ignition (SI) engine modified to run in Homogeneous Charge Compression Ignition (HCCI) mode with a Negative Valve Overlap (NVO). A portion of the fuel is injected during NVO to increase the cylinder temperature and enable HCCI combustion at a compression ratio of 12:1. The model is coupled with GT-Power, a one-dimensional engine simulation tool used for the open valve portion of the engine cycle. The SRM is used to model in-cylinder mixing, heat transfer and chemistry during the NVO and main combustion. Direct injection is simulated during NVO in order to predict heat release and internal Exhaust Gas Recycle (EGR) composition and mass. The NOx emissions and simulated pressure profiles match experimental data well, including the cyclic fluctuations. The model predicts combustion characteristics at different fuel split ratios and injection timings. The effect of fuel reforming on ignition timing is investigated along with the causes of cycle to cycle variations and unstable operation. A detailed flux analysis during NVO unearths interesting results regarding the effect of NOx on ignition timing compared with its effect during the main combustion. © 2009 SAE International.

Reconsideration On The Role Of The Specific Heat Ratio In Arrhenius Law Applications

Arrhenius law implicates that only those molecules which possess the internal energy greater than the activation energy E-a can react. However, the internal energy will not be proportional to the gas temperature if the specific heat ratio gamma and the gas constant R vary during chemical reaction processes. The varying gamma may affect significantly the chemical reaction rate calculated with the Arrhenius law under the constant gamma assumption, which has been widely accepted in detonation and combustion simulations for many years. In this paper, the roles of variable gamma and R in Arrhenius law applications are reconsidered, and their effects on the chemical reaction rate are demonstrated by simulating one-dimensional C-J and two-dimensional cellular detonations. A new overall one-step detonation model with variable gamma and R is proposed to improve the Arrhenius law. Numerical experiments demonstrate that this improved Arrhenius law works well in predicting detonation phenomena with the numerical results being in good agreement with experimental data.

Modeling, numerical methods, and simulation for particle-fluid two-phase flow problems

In this paper, we study the issues of modeling, numerical methods, and simulation with comparison to experimental data for the particle-fluid two-phase flow problem involving a solid-liquid mixed medium. The physical situation being considered is a pulsed liquid fluidized bed. The mathematical model is based on the assumption of one-dimensional flows, incompressible in both particle and fluid phases, equal particle diameters, and the wall friction force on both phases being ignored. The model consists of a set of coupled differential equations describing the conservation of mass and momentum in both phases with coupling and interaction between the two phases. We demonstrate conditions under which the system is either mathematically well posed or ill posed. We consider the general model with additional physical viscosities and/or additional virtual mass forces, both of which stabilize the system. Two numerical methods, one of them is first-order accurate and the other fifth-order accurate, are used to solve the models. A change of variable technique effectively handles the changing domain and boundary conditions. The numerical methods are demonstrated to be stable and convergent through careful numerical experiments. Simulation results for realistic pulsed liquid fluidized bed are provided and compared with experimental data. (C) 2004 Elsevier Ltd. All rights reserved.

Renormalization group approach on nanostructured systems

The electronic spectra of one-dimensional nanostructured systems are calculated within the pure hopping model on the tight-binding Hamiltonian. By means of the renormalization group Green's function method, the dependence of the density of states on the distributions of nanoscaled grains and the changes of values of hopping integrals in nanostructured systems are studied. It is found that the frequency shifts are dependent rather on the changes of the hopping integrals at nanoscaled grains than the distribution of nanoscaled grains.

A perturbational h4 exponential finite-difference scheme for the convective diffusion equation

A perturbational h4 compact exponential finite difference scheme with diagonally dominant coefficient matrix and upwind effect is developed for the convective diffusion equation. Perturbations of second order are exerted on the convective coefficients and source term of an h2 exponential finite difference scheme proposed in this paper based on a transformation to eliminate the upwind effect of the convective diffusion equation. Four numerical examples including one- to three-dimensional model equations of fluid flow and a problem of natural convective heat transfer are given to illustrate the excellent behavior of the present exponential schemes, the h4 accuracy of the perturbational scheme is verified using double precision arithmetic.

Statistical modeling of damage evolution in spallation

In order to understand the mechanism of the incipient spallation in rolled metals, a one dimensional statistical mode1 on evolution of microcracks in spallation was proposed. The crack length appears to be the fundamental variable in the statistical description. Two dynamic processes, crack nucleation and growth, were involved in the model of damage evolution. A simplified case was examined and preliminary correlation to experimental observations of spallation was made.

[G_Hondarribia_Muralla] Documentación geométrica de las murallas de Hondarribia

[ES] Al respecto del modelo virtual con texturas fotográficas generado en este proyecto, se puede consultar también el siguiente artículo:

Experimental Studies on Scramjet Tested in a Freejet Facility

Two different type scramjet models with side-wall compression and top-wall compression inlets have been tested in HPTF (Hypersonic Propulsion Test Facility) under the experimental conditions of Mach number 5.8, total temperature 1700K, total pressure 4.5MPa and mass flow rate 3.5kg/s. The liquid kerosene was used as main fuel for the scramjets. In order to get fast ignition in the combustor, a small amount of hydrogen was used as a pilot. A strut with alternative tail was employed for increasing the compression ratio and for mixing enhancement in the side-wall compression case. Recessed cavities were used as a flameholder for combustion stability. The combustion efficiency was estimated by one dimensional theory. The uniformity of the facility nozzle flow was verified by a scanning pitot rake. The experimental results showed that the kerosene fuel was successfully ignited and stable combustion was achieved for both scramjet models. However the thrusts were still less than the model drags due to the low combustion efficiencies.

Modeling the Chesapeake Bay and Tributaries: a synopsis

The last decade has seen the development and application of a spectrum of physical and numerical hydrographic models of the Chesapeake Bay and its tributaries. The success of the James River Hydraulic Model has initiated the construction of an estuarine hydraulic model of the entire Chesapeake System. Numerical analogues for hydrographic behavior and contaminant dispersion in one-, two-, and three dimensional model estuaries exist for various regions of the Bay. From an engineering viewpoint, one dimensional models are sufficiently advanced to be routinely employed in aiding management decisions. Bay investigators are playing leading roles in the development of two- and three-dimensional models of estuarine flows.

Simulation of Natural Gas Production in Hydrate Reservoirs

This paper simulates a one-dimensional physical model of natural gas production from hydrate dissociation in a reservoir by depressurization. According to the principles of solid hydrate decomposition in stratum and flow of natural gas in porous medium, the pressure governing equations for both gas zone and hydrate zone are set up based on the physical production model. Using the approximation reported by N. N. Verigin et al. (1980), the nonlinear governing equations are simplified and the self-similar solutions are obtained. Through calculation, for different reservoir parameters, the distribution characters of pressure are analyzed. The decline trend of natural gas production rate with time is also studied. The simulation results show that production of natural gas from a hydrate reservoir is very sensitive to several reservoir parameters, such as wellbore pressure and stratum porosity and permeability.

On the origin of shear banding instability in metallic glasses

To uncover the physical origin of shear-banding instability in metallic glass (MG), a theoretical description of thermo-mechanical deformation of MG undergoing one-dimensional simple shearing is presented. The coupled thermo-mechanical model takes into account the momentum balance, the energy balance and the dynamics of free volume. The interplay between free-volume production and temperature increase being two potential causes for shear-banding instability is examined on the basis of the homogeneous solution. It is found that the free-volume production facilitates the sudden increase in the temperature before instability and vice versa. A rigorous linear perturbation analysis is used to examine the inhomogeneous deformation, during which the onset criteria and the internal length and time scales for three types of instabilities, namely free-volume softening, thermal softening and coupling softening, are clearly revealed. The shear-banding instability originating from sole free-volume softening takes place easier and faster than that due to sole thermal softening, and dominates in the coupling softening. Furthermore, the coupled thermo-mechanical shear-band analysis does show that an initial slight distribution of local free volume can incur significant strain localization, producing a shear band. During such a localization process, the local free-volume creation occurs indeed prior to the increase in local temperature, indicating that the former is the cause of shear localization, whereas the latter is its consequence. Finally, extension of the above model to include the shear-induced dilatation shows that such dilatation facilitates the shear instability in metallic glasses.

High Order Multi-Moment Constrained Finite Volume Method. Part I: Basic Formulation

A new high-order finite volume method based on local reconstruction is presented in this paper. The method, so-called the multi-moment constrained finite volume (MCV) method, uses the point values defined within single cell at equally spaced points as the model variables (or unknowns). The time evolution equations used to update the unknowns are derived from a set of constraint conditions imposed on multi kinds of moments, i.e. the cell-averaged value and the point-wise value of the state variable and its derivatives. The finite volume constraint on the cell-average guarantees the numerical conservativeness of the method. Most constraint conditions are imposed on the cell boundaries, where the numerical flux and its derivatives are solved as general Riemann problems. A multi-moment constrained Lagrange interpolation reconstruction for the demanded order of accuracy is constructed over single cell and converts the evolution equations of the moments to those of the unknowns. The presented method provides a general framework to construct efficient schemes of high orders. The basic formulations for hyperbolic conservation laws in 1- and 2D structured grids are detailed with the numerical results of widely used benchmark tests. (C) 2009 Elsevier Inc. All rights reserved.

Analog VLSI-based, neuromorphic sensorimotor systems: modeling the primate oculomotor system

Using neuromorphic analog VLSI techniques for modeling large neural systems has several advantages over software techniques. By designing massively-parallel analog circuit arrays which are ubiquitous in neural systems, analog VLSI models are extremely fast, particularly when local interactions are important in the computation. While analog VLSI circuits are not as flexible as software methods, the constraints posed by this approach are often very similar to the constraints faced by biological systems. As a result, these constraints can offer many insights into the solutions found by evolution. This dissertation describes a hardware modeling effort to mimic the primate oculomotor system which requires both fast sensory processing and fast motor control. A one-dimensional hardware model of the primate eye has been built which simulates the physical dynamics of the biological system. It is driven by analog VLSI circuits mimicking brainstem and cortical circuits that control eye movements. In this framework, a visually-triggered saccadic system is demonstrated which generates averaging saccades. In addition, an auditory localization system, based on the neural circuits of the barn owl, is used to trigger saccades to acoustic targets in parallel with visual targets. Two different types of learning are also demonstrated on the saccadic system using floating-gate technology allowing the non-volatile storage of analog parameters directly on the chip. Finally, a model of visual attention is used to select and track moving targets against textured backgrounds, driving both saccadic and smooth pursuit eye movements to maintain the image of the target in the center of the field of view. This system represents one of the few efforts in this field to integrate both neuromorphic sensory processing and motor control in a closed-loop fashion.

I. Rigid body penetration into brittle materials. II. Phase change effect on shock wave propagation

Part I.

We have developed a technique for measuring the depth time history of rigid body penetration into brittle materials (hard rocks and concretes) under a deceleration of ~ 10⁵ g. The technique includes bar-coded projectile, sabot-projectile separation, detection and recording systems. Because the technique can give very dense data on penetration depth time history, penetration velocity can be deduced. Error analysis shows that the technique has a small intrinsic error of ~ 3-4 % in time during penetration, and 0.3 to 0.7 mm in penetration depth. A series of 4140 steel projectile penetration into G-mixture mortar targets have been conducted using the Caltech 40 mm gas/ powder gun in the velocity range of 100 to 500 m/s.

We report, for the first time, the whole depth-time history of rigid body penetration into brittle materials (the G-mixture mortar) under 10⁵ g deceleration. Based on the experimental results, including penetration depth time history, damage of recovered target and projectile materials and theoretical analysis, we find:

1. Target materials are damaged via compacting in the region in front of a projectile and via brittle radial and lateral crack propagation in the region surrounding the penetration path. The results suggest that expected cracks in front of penetrators may be stopped by a comminuted region that is induced by wave propagation. Aggregate erosion on the projectile lateral surface is < 20% of the final penetration depth. This result suggests that the effect of lateral friction on the penetration process can be ignored.

2. Final penetration depth, P_max, is linearly scaled with initial projectile energy per unit cross-section area, e_s , when targets are intact after impact. Based on the experimental data on the mortar targets, the relation is P_max(mm) 1.15e_s (J/mm² ) + 16.39.

3. Estimation of the energy needed to create an unit penetration volume suggests that the average pressure acting on the target material during penetration is ~ 10 to 20 times higher than the unconfined strength of target materials under quasi-static loading, and 3 to 4 times higher than the possible highest pressure due to friction and material strength and its rate dependence. In addition, the experimental data show that the interaction between cracks and the target free surface significantly affects the penetration process.

4. Based on the fact that the penetration duration, t_max, increases slowly with e_s and does not depend on projectile radius approximately, the dependence of t_max on projectile length is suggested to be described by t_max(μs) = 2.08e_s (J/mm² + 349.0 x m/(πR²), in which m is the projectile mass in grams and R is the projectile radius in mm. The prediction from this relation is in reasonable agreement with the experimental data for different projectile lengths.

5. Deduced penetration velocity time histories suggest that whole penetration history is divided into three stages: (1) An initial stage in which the projectile velocity change is small due to very small contact area between the projectile and target materials; (2) A steady penetration stage in which projectile velocity continues to decrease smoothly; (3) A penetration stop stage in which projectile deceleration jumps up when velocities are close to a critical value of ~ 35 m/s.

6. Deduced averaged deceleration, a, in the steady penetration stage for projectiles with same dimensions is found to be a(g) = 192.4v + 1.89 x 10⁴, where v is initial projectile velocity in m/s. The average pressure acting on target materials during penetration is estimated to be very comparable to shock wave pressure.

7. A similarity of penetration process is found to be described by a relation between normalized penetration depth, P/P_max, and normalized penetration time, t/t_max, as P/P_max = f(t/t_max, where f is a function of t/t_max. After f(t/t_max is determined using experimental data for projectiles with 150 mm length, the penetration depth time history for projectiles with 100 mm length predicted by this relation is in good agreement with experimental data. This similarity also predicts that average deceleration increases with decreasing projectile length, that is verified by the experimental data.

8. Based on the penetration process analysis and the present data, a first principle model for rigid body penetration is suggested. The model incorporates the models for contact area between projectile and target materials, friction coefficient, penetration stop criterion, and normal stress on the projectile surface. The most important assumptions used in the model are: (1) The penetration process can be treated as a series of impact events, therefore, pressure normal to projectile surface is estimated using the Hugoniot relation of target material; (2) The necessary condition for penetration is that the pressure acting on target materials is not lower than the Hugoniot elastic limit; (3) The friction force on projectile lateral surface can be ignored due to cavitation during penetration. All the parameters involved in the model are determined based on independent experimental data. The penetration depth time histories predicted from the model are in good agreement with the experimental data.

9. Based on planar impact and previous quasi-static experimental data, the strain rate dependence of the mortar compressive strength is described by σ_f/σ⁰_f = exp(0.0905(log(έ/έ_0) ^1.14, in the strain rate range of 10^-7/s to 10³/s (σ⁰_f and έ are reference compressive strength and strain rate, respectively). The non-dispersive Hugoniot elastic wave in the G-mixture has an amplitude of ~ 0.14 GPa and a velocity of ~ 4.3 km/s.

Part II.

Stress wave profiles in vitreous GeO₂ were measured using piezoresistance gauges in the pressure range of 5 to 18 GPa under planar plate and spherical projectile impact. Experimental data show that the response of vitreous GeO₂ to planar shock loading can be divided into three stages: (1) A ramp elastic precursor has peak amplitude of 4 GPa and peak particle velocity of 333 m/s. Wave velocity decreases from initial longitudinal elastic wave velocity of 3.5 km/s to 2.9 km/s at 4 GPa; (2) A ramp wave with amplitude of 2.11 GPa follows the precursor when peak loading pressure is 8.4 GPa. Wave velocity drops to the value below bulk wave velocity in this stage; (3) A shock wave achieving final shock state forms when peak pressure is > 6 GPa. The Hugoniot relation is D = 0.917 + 1.711u (km/s) using present data and the data of Jackson and Ahrens [1979] when shock wave pressure is between 6 and 40 GPa for ρ₀ = 3.655 gj cm³ . Based on the present data, the phase change from 4-fold to 6-fold coordination of Ge⁺⁴ with O^-2 in vitreous GeO₂ occurs in the pressure range of 4 to 15 ± 1 GPa under planar shock loading. Comparison of the shock loading data for fused SiO₂ to that on vitreous GeO₂ demonstrates that transformation to the rutile structure in both media are similar. The Hugoniots of vitreous GeO₂ and fused SiO₂ are found to coincide approximately if pressure in fused SiO₂ is scaled by the ratio of fused SiO₂to vitreous GeO₂ density. This result, as well as the same structure, provides the basis for considering vitreous Ge0₂ as an analogous material to fused SiO₂ under shock loading. Experimental results from the spherical projectile impact demonstrate: (1) The supported elastic shock in fused SiO₂ decays less rapidly than a linear elastic wave when elastic wave stress amplitude is higher than 4 GPa. The supported elastic shock in vitreous GeO₂ decays faster than a linear elastic wave; (2) In vitreous GeO₂ , unsupported shock waves decays with peak pressure in the phase transition range (4-15 GPa) with propagation distance, x, as α 1/x^-3.35 , close to the prediction of Chen et al. [1998]. Based on a simple analysis on spherical wave propagation, we find that the different decay rates of a spherical elastic wave in fused SiO₂ and vitreous GeO₂ is predictable on the base of the compressibility variation with stress under one-dimensional strain condition in the two materials.

Photodissociation and reaction dynamics of chlorine-containing species important in stratospheric ozone chemistry

Chlorine oxide species have received considerable attention in recent years due to their central role in the balance of stratospheric ozone. Many questions pertaining to the behavior of such species still remain unanswered and plague the ability of researchers to develop accurate chemical models of the stratosphere. Presented in this thesis are three experiments that study various properties of some specific chlorine oxide species.

In the first chapter, the reaction between ClONO_2 and protonated water clusters is investigated to elucidate a possible reaction mechanism for the heterogeneous reaction of chlorine nitrate on ice. The ionic products were various forms of protonated nitric acid, NO_2 +(H_20)_m, m = 0, 1, 2. These products are analogous to products previously reported in the literature for the neutral reaction occurring on ice surfaces. Our results support the hypothesis that the heterogeneous reaction is acid-catalyzed.

In the second chapter, the photochemistry of ClONO_2 was investigated at two wavelengths, 193 and 248 nm, using the technique of photofragmentation translational spectroscopy. At both wavelengths, the predominant dissociation pathways were Cl + NO_3 and ClO + NO_2. Channel assignments were confirmed by momentum matching the counterfragments from each channel. A one-dimensional stratospheric model using the new 248 nm branching ratio determined how our results would affect the predicted Cl_x and NO_x partitioning in the stratosphere.

Chapter three explores the photodissociation dynamics of Cl_2O at 193, 248 and 308 nm. At 193 nm, we found evidence for the concerted reaction channel, Cl_2 + O. The ClO + Cl channel was also accessed, however, the majority of the ClO fragments were formed with sufficient internal energies for spontaneous secondary dissociation to occur. At 248 and 308 nm, we only observed only the ClO + Cl channel. . Some of the ClO formed at 248 nm was formed internally hot and spontaneously dissociated. Bimodal translational energy distributions of the ClO and Cl products indicate two pathways leading to the same product exist.

Appendix A, B and C discuss the details of data analysis techniques used in Chapters 1 and 2. The development of a molecular beam source of ClO dimer is presented in Appendix D.