Coupling of Chemical Kinetics with Computational Fluid Dynamics in a Three-Dimensional Engine Model

The role of computer modeling has grown recently to integrate itself as an inseparable tool to experimental studies for the optimization of automotive engines and the development of future fuels. Traditionally, computer models rely on simplified global reaction steps to simulate the combustion and pollutant formation inside the internal combustion engine. With the current interest in advanced combustion modes and injection strategies, this approach depends on arbitrary adjustment of model parameters that could reduce credibility of the predictions. The purpose of this study is to enhance the combustion model of KIVA, a computational fluid dynamics code, by coupling its fluid mechanics solution with detailed kinetic reactions solved by the chemistry solver, CHEMKIN. As a result, an engine-friendly reaction mechanism for n-heptane was selected to simulate diesel oxidation. Each cell in the computational domain is considered as a perfectly-stirred reactor which undergoes adiabatic constant- volume combustion. The model was applied to an ideally-prepared homogeneous- charge compression-ignition combustion (HCCI) and direct injection (DI) diesel combustion. Ignition and combustion results show that the code successfully simulates the premixed HCCI scenario when compared to traditional combustion models. Direct injection cases, on the other hand, do not offer a reliable prediction mainly due to the lack of turbulent-mixing model, inherent in the perfectly-stirred reactor formulation. In addition, the model is sensitive to intake conditions and experimental uncertainties which require implementation of enhanced predictive tools. It is recommended that future improvements consider turbulent-mixing effects as well as optimization techniques to accurately simulate actual in-cylinder process with reduced computational cost. Furthermore, the model requires the extension of existing fuel oxidation mechanisms to include pollutant formation kinetics for emission control studies.

Unexpected controllable pair-structure in ferroelectric nanodomains.

The imminent inability of silicon-based memory devices to satisfy Moore's Law is approaching rapidly. Controllable nanodomains of ferroic systems are anticipated to enable future high-density nonvolatile memory and novel electronic devices. We find via piezoresponse force microscopy (PFM) studies on lead zirconate titanate (PZT) films an unexpected nanostructuring of ferroelectric-ferroelastic domains. These consist of c-nanodomains within a-nanodomains in proximity to a-nanodomains within c-domains. These structures are created and annihilated as pairs, controllably. We treat these as a new kind of vertex-antivertex pair and consider them in terms of the Srolovitz-Scott 4-state Potts model, which results in pairwise domain vertex instabilities that resemble the vortex-antivortex mechanism in ferromagnetism, as well as dislocation pairs (or disclination pairs) that are well-known in nematic liquid crystals. Finally, we show that these nanopairs can be scaled up to form arrays that are engineered at will, paving the way toward facilitating them to real technologies.

Freshwater influences on hydrology and seagrass dynamics of intermittent estuaries

Many temperate estuaries have intermittently open and closed mouths, a feature that is often related to intermittent freshwater input. These systems, often overlooked due to their small size, can have large hydrological variability over medium-term time scales.

This variability presents potential difficulties for estuarine species particularly where anthropogenic alterations to freshwater flows can cause large deviations from natural patterns of tidal influence and inundation of habitat.

Influences of natural and hydrological variability on seagrasses were examined in two central Victorian estuaries with anthropogenically-modified but naturally-intermittent freshwater flows and mouth openings. Comparisons were focused on differences between an estuary with artificially-augmented freshwater inflow and an adjacent system, in which the volume and timing of inflows were altered by a reservoir. Eight additional estuaries in the region were also used to provide a context for these two main sites.

Hydrological changes during the three-year field component were affected by the ending of a drought and then a major flood a year later as well as by ongoing anthropogenic flow reduction and augmentation. These influences on hydrology were associated with an initially high seagrass coverage that was substantially reduced and showed signs of recovery only in the system that was affected by lower inflows. Such influences and responses also changed seasonally but to a much lesser extent than the responses to stochastic climatic events.

Natural flows were intermittent and varied substantially between years. Flooding flows represented up to 89% of the long-term annual average flow. Water quality was broadly typical of the region, with the exception of low pH in some tributaries, especially those of Anglesea estuary. Anthropogenic changes to flow were most evident at times of low natural flows and resulted in longer and more frequent periods of zero inflow to Painkalac estuary and a continual base flow to Anglesea. This base flow, from ponds containing coal ash, neutralised waters flowing from upstream and increased conductivity, except at times of high natural flow.

A three-state conceptual model of the magnitude and variability of water levels, based largely on the degree of tidal influence was identified and quantitatively assessed for the two estuaries that were the main focus of the study. These states in turn had a large influence on the area and inundation of benthic habitat. Floods tended to open the mouths of estuaries, which then remained tidal given sufficient flow to overcome sedimentary processes at the mouths. Low and zero inflow was a precondition for closure of the mouths of the estuaries. When closed, differences in inflow resulted in different endpoints in salinity patterns. From an initial pattern similar to a classic ‘salt wedge’, Painkalac estuary, with reduced inflow, quickly destratified and gradually became more saline, at times hypersaline. Anglesea estuary, with augmented flow, tended to remain stratified for longer until becoming completely fresh, given a long enough period of closure.

Episodic changes in the water quality of the estuaries were associated with different components of the freshwater flow regimes. At high flows, fresh waters of low pH with a high metal load entered Anglesea estuary. Except during the largest flood, when the estuary was completely flushed, this water was neutralised at the halocline and resulting in precipitation of metals. High flows into Painkalac were associated with elevated concentrations of clay-sourced suspended solids. During a closed period, with zero flow, a release of sediment-bound nutrients triggered by anoxia was observed in Painkalac, followed by an algal bloom.

The large decline in seagrass extent that was observed in both estuaries was closely related to floods and the subsequent reductions in potential habitat associated with the tidal states that followed. Analysis of historical patterns of extent against rainfall records suggested that periods of drought and extended mouth closures were related to establishment and expansion of beds. This model was similar to that described for South African estuaries and contrasted with more-seasonal patterns reported for local marine embayments.

Rates of in situ decomposition of seagrass detritus showed a mix of seasonal and disturbance-driven patterns of change, depending on estuary. Variability of these rates on a scale of 100s of metres was typically not significant, but there were a few episodes that were highly significant. A negative correlation between decomposition rate and seagrass extent was also observed. A novel technique for assessing cellulose decomposition potential in sediment, adapted from soil science, proved to be a useful tool for estuarine research. Results from this component of the study highlighted both small-scale variability that was inconsistent through time, and also stable differences in decomposition potential between depths and estuaries that were consistent with differences in hydrological state and salinity.

Given the relative lack of knowledge about processes in intermittent estuaries, particularly those relating to changes in freshwater inflow, results from this study will be of value both locally and for similar systems elsewhere. Locally, it is likely that flow regimes to both Anglesea and Painkalac estuaries will be reduced, following closure of the mine power station at Anglesea and due to increased demand from the reservoir above Painkalac. There is potential to manage flows from each of these sources to minimise downstream effects. Regionally, and globally, there are many intermittent estuaries in areas with Mediterranean-type climates. It has been predicted that the climates of these regions will become drier but with an increase in intensity of storm events, both of which have ramifications for flow regimes to estuaries. It is hoped that results of this study will contribute to more informed management of intermittent estuaries in the context of these likely changes.

Propagação de danos em sistemas cooperativos: propriedades termodinâmicas

High-precision calculations of the correlation functions and order parameters were performed in order to investigate the critical properties of several two-dimensional ferro- magnetic systems: (i) the q-state Potts model; (ii) the Ashkin-Teller isotropic model; (iii) the spin-1 Ising model. We deduced exact relations connecting specific damages (the difference between two microscopic configurations of a model) and the above mentioned thermodynamic quanti- ties which permit its numerical calculation, by computer simulation and using any ergodic dynamics. The results obtained (critical temperature and exponents) reproduced all the known values, with an agreement up to several significant figures; of particular relevance were the estimates along the Baxter critical line (Ashkin-Teller model) where the exponents have a continuous variation. We also showed that this approach is less sensitive to the finite-size effects than the standard Monte-Carlo method. This analysis shows that the present approach produces equal or more accurate results, as compared to the usual Monte Carlo simulation, and can be useful to investigate these models in circumstances for which their behavior is not yet fully understood

Generalized partition function zeros of 1D spin models and their critical behavior at edge singularities

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Estimation of electronic coupling in π -stacked donor-bridge-acceptor systems: correction of the two-state model

Comparison of donor-acceptor electronic couplings calculated within two-state and three-state models suggests that the two-state treatment can provide unreliable estimates of Vda because of neglecting the multistate effects. We show that in most cases accurate values of the electronic coupling in a π stack, where donor and acceptor are separated by a bridging unit, can be obtained as Ṽ da = (E2 - E1) μ12 Rda + (2 E3 - E1 - E2) 2 μ13 μ23 Rda2, where E1, E2, and E3 are adiabatic energies of the ground, charge-transfer, and bridge states, respectively, μij is the transition dipole moments between the states i and j, and Rda is the distance between the planes of donor and acceptor. In this expression based on the generalized Mulliken-Hush approach, the first term corresponds to the coupling derived within a two-state model, whereas the second term is the superexchange correction accounting for the bridge effect. The formula is extended to bridges consisting of several subunits. The influence of the donor-acceptor energy mismatch on the excess charge distribution, adiabatic dipole and transition moments, and electronic couplings is examined. A diagnostic is developed to determine whether the two-state approach can be applied. Based on numerical results, we showed that the superexchange correction considerably improves estimates of the donor-acceptor coupling derived within a two-state approach. In most cases when the two-state scheme fails, the formula gives reliable results which are in good agreement (within 5%) with the data of the three-state generalized Mulliken-Hush model

Positronium-hydrogen-atom scattering in a five-state model

The scattering of orthopositronium (Ps) by hydrogen atoms has been investigated in a five-state coupled-channel model allowing for Ps(1s)H(2s,2p) and Ps(2s,2p)H(1s) excitations using a recently proposed electron-exchange model potential. The higher (n greater than or equal to 3) excitations and ionization of the Ps atom are calculated using the first Born approximation. Calculations are reported of scattering lengths, phase shifts. elastic, Ps and H excitation, and total cross sections. Remarkable correlations are observed between the S-wave Ps-H binding energy and the singlet scattering length, effective range, and resonance energy obtained in various model calculations. These correlations suggest that if a Ps-H dynamical model yields the correct result for one of these four observables, it is expected to lead to the correct result for the other three. The present model, which is constructed so as to reproduce the Ps-H resonance at 4.01 eV, automatically yields a Ps-H bound state at - 1.05 eV that compares well with the accurate value of - 1.067 eV. The model leads to a singlet scattering length of 3.72a(0) and effective range of 1.67a(0), whereas the correlations suggest the precise values of 3.50a(0) and 1.65a(0) for these observables, respectively. [S1050-2947(99)07703-3].

Critical behavior of three-dimensional disordered Potts models with many states

We study the 3D Disordered Potts Model with p = 5 and p = 6. Our numerical simulations (that severely slow down for increasing p) detect a very clear spin glass phase transition. We evaluate the critical exponents and the critical value of the temperature, and we use known results at lower p values to discuss how they evolve for increasing p. We do not find any sign of the presence of a transition to a ferromagnetic regime.

A three species model to simulate application of hyperbaric oxygen therapy to chronic wounds

Chronic wounds are a significant socioeconomic problem for governments worldwide. Approximately 15% of people who suffer from diabetes will experience a lower-limb ulcer at some stage of their lives, and 24% of these wounds will ultimately result in amputation of the lower limb. Hyperbaric Oxygen Therapy (HBOT) has been shown to aid the healing of chronic wounds; however, the causal reasons for the improved healing remain unclear and hence current HBOT protocols remain empirical. Here we develop a three-species mathematical model of wound healing that is used to simulate the application of hyperbaric oxygen therapy in the treatment of wounds. Based on our modelling, we predict that intermittent HBOT will assist chronic wound healing while normobaric oxygen is ineffective in treating such wounds. Furthermore, treatment should continue until healing is complete, and HBOT will not stimulate healing under all circumstances, leading us to conclude that finding the right protocol for an individual patient is crucial if HBOT is to be effective. We provide constraints that depend on the model parameters for the range of HBOT protocols that will stimulate healing. More specifically, we predict that patients with a poor arterial supply of oxygen, high consumption of oxygen by the wound tissue, chronically hypoxic wounds, and/or a dysfunctional endothelial cell response to oxygen are at risk of nonresponsiveness to HBOT. The work of this paper can, in some way, highlight which patients are most likely to respond well to HBOT (for example, those with a good arterial supply), and thus has the potential to assist in improving both the success rate and hence the costeffectiveness of this therapy.

A semantics-based approach to querying process model repositories

As business process management technology matures, organisations acquire more and more business process models. The resulting collections can consist of hundreds, even thousands of models and their management poses real challenges. One of these challenges concerns model retrieval where support should be provided for the formulation and efficient execution of business process model queries. As queries based on only structural information cannot deal with all querying requirements in practice, there should be support for queries that require knowledge of process model semantics. In this paper we formally define a process model query language that is based on semantic relationships between tasks. This query language is independent of the particular process modelling notation used, but we will demonstrate how it can be used in the context of Petri nets by showing how the semantic relationships can be determined for these nets in such a way that state space explosion is avoided as much as possible. An experiment with three large process model repositories shows that queries expressed in our language can be evaluated efficiently.

Numerical simulation of idealised three-dimensional downburst wind fields

Convective downburst wind storms generate the peak annual gust wind speed for many parts of the non-cyclonic world at return periods of importance for ultimate limit state design. Despite this there is little clear understanding of how to appropriately design for these wind events given their significant dissimilarities to boundary layer winds upon which most design is based. To enhance the understanding of wind fields associated with these storms a three-dimensional numerical model was developed to simulate a multitude of idealised downburst scenarios and to investigate their near-ground wind characteristics. Stationary and translating downdraft wind events in still and sheared environments were simulated with baseline results showing good agreement with previous numerical work and full-scale observational data. Significant differences are shown in the normalised peak wind speed velocity profiles depending on the environmental wind conditions in the vicinity of the simulated event. When integrated over the height of mid- to high rise structures, all simulated profiles are shown to produce wind loads smaller than an equivalent 10 m height matched open terrain boundary layer profile. This suggests that for these structures the current design approach is conservative from an ultimate loading standpoint. Investigating the influence of topography on the structure of the simulated near-ground downburst wind fields, it is shown that these features amplify wind speeds in a manner similar to that expected for boundary layer winds, but the extent of amplification is reduced. The level of reduction is shown to be dependent on the depth of the simulated downburst outflow.

State space formulation of ammonia reactor dynamics

The specific objective of this paper is to develop a state space model of a tubular ammonia reactor which is the heart of an ammonia plant in a fertiliser complex. A ninth order model with three control inputs and two disturbance inputs is generated from the nonlinear distributed model using linearization and lumping approximations. The lumped model is chosen such that the steady state temperature at the exit of the catalyst bed computed from the simplified state space model is close enough to the one computed from the nonlinear steady state model. The model developed in this paper is very useful for the design of continuous/discrete versions of single variable/multivariable control algorithms.

Flap‐Lag Damping in Hover and Forward Flight with a Three‐Dimensional Wake

Prediction of lag damping is difficult owing to the delicate balance of drag, induced drag and Coriolis forces in the in‐plane direction. Moreover, induced drag” is sensitive to dynamic wake, bath shed and trailing components, and thus its prediction requires adequate unsteady‐wake representation. Accordingly, rigid‐blade flap‐lag equations are coupled with a three‐dimensional finite‐state wake model; three isolatcd rotor canfigurations with three, four and five blades are treated over a range of thrust levels, tack numbers, lag frequencies and advance ratios. The investigation includes convergence characteristics of damping with respect to the number of radial shape functions and harmonics of the wake model for multiblade modes of low frequency (< 1/ rev.) to high frequency (> 1/rev.). Predicted flap and lag damping levels are then compared with similar predictions with 1) rigid wake (no unsteady induced now), 2) Loewy lift deficiency and 3) dynamic inflow. The coverage also includes correlations with the measured lag regressive‐mode damping in hover and forward flight and comparisons with similar correlations with dynamic inflow. Lag‐damping predictions with the dynamic wake model are consistently higher than the predictions with the dynamic inflow model; even for the low frequency lag regressive mode, the number of wake harmonics should at least be equal to twice the number of blades.

125 GeV Higgs state in the context of four generations with two Higgs doublets

We interpret the recent discovery of a 125 GeV Higgs-like state in the context of a two-Higgs-doublet model with a heavy fourth sequential generation of fermions, in which one Higgs doublet couples only to the fourth-generation fermions, while the second doublet couples to the lighter fermions of the first three families. This model is designed to accommodate the apparent heaviness of the fourth-generation fermions and to effectively address the low-energy phenomenology of a dynamical electroweak-symmetry-breaking scenario. The physical Higgs states of the model are, therefore, viewed as composites primarily of the fourth-generation fermions. We find that the lightest Higgs, h, is a good candidate for the recently discovered 125 GeV spin-zero particle, when tan beta similar to O(1), for typical fourth-generation fermion masses of M-4G = 400-600 GeV, and with a large t-t' mixing in the right-handed quark sector. This, in turn, leads to BR(t' -> th) similar to O(1), which drastically changes the t' decay pattern. We also find that, based on the current Higgs data, this two-Higgs-doublet model generically predicts an enhanced production rate (compared to the Standard Model) in the pp -> h -> tau tau channel, and reduced rates in the VV -> h -> gamma gamma and p (p) over bar /pp -> V -> hV -> Vbb channels. Finally, the heavier CP-even Higgs is excluded by the current data up to m(H) similar to 500 GeV, while the pseudoscalar state, A, can be as light as 130 GeV. These heavier Higgs states and the expected deviations from the Standard Model din some of the Higgs production channels can be further excluded or discovered with more data.

A multi-sheeted three-dimensional potential-energy surface for the H-atom photodissociation of phenol

Extending the previous work of Lan et al. J. Chem. Phys., 122, 224315 (2005)], a multi-state potential model for the H atom photodissociation is presented. All three ``disappearing coordinates'' of the departing H atom have been considered. Ab initio CASSCF computations have been carried out for the linear COH geometry of C-2v symmetry, and for several COH angles with the OH group in the ring plane and also perpendicular to the ring plane. By keeping the C6H5O fragment frozen in a C-2v-constrained geometry throughout, we have been able to apply symmetry-based simplifications in the constructions of a diabatic model. This model is able to capture the overall trends of twelve adiabats at both torsional limits for a wide range of COH bend angles.