Development of a model-based transient calibration process for diesel engine electronic control module tables – Part 1: data requirements, processing, and analysis

This is the first part of a study investigating a model-based transient calibration process for diesel engines. The motivation is to populate hundreds of parameters (which can be calibrated) in a methodical and optimum manner by using model-based optimization in conjunction with the manual process so that, relative to the manual process used by itself, a significant improvement in transient emissions and fuel consumption and a sizable reduction in calibration time and test cell requirements is achieved. Empirical transient modelling and optimization has been addressed in the second part of this work, while the required data for model training and generalization are the focus of the current work. Transient and steady-state data from a turbocharged multicylinder diesel engine have been examined from a model training perspective. A single-cylinder engine with external air-handling has been used to expand the steady-state data to encompass transient parameter space. Based on comparative model performance and differences in the non-parametric space, primarily driven by a high engine difference between exhaust and intake manifold pressures (ΔP) during transients, it has been recommended that transient emission models should be trained with transient training data. It has been shown that electronic control module (ECM) estimates of transient charge flow and the exhaust gas recirculation (EGR) fraction cannot be accurate at the high engine ΔP frequently encountered during transient operation, and that such estimates do not account for cylinder-to-cylinder variation. The effects of high engine ΔP must therefore be incorporated empirically by using transient data generated from a spectrum of transient calibrations. Specific recommendations on how to choose such calibrations, how many data to acquire, and how to specify transient segments for data acquisition have been made. Methods to process transient data to account for transport delays and sensor lags have been developed. The processed data have then been visualized using statistical means to understand transient emission formation. Two modes of transient opacity formation have been observed and described. The first mode is driven by high engine ΔP and low fresh air flowrates, while the second mode is driven by high engine ΔP and high EGR flowrates. The EGR fraction is inaccurately estimated at both modes, while EGR distribution has been shown to be present but unaccounted for by the ECM. The two modes and associated phenomena are essential to understanding why transient emission models are calibration dependent and furthermore how to choose training data that will result in good model generalization.

Development of Measurement Methods for Testing of Hydrokinetic Devices to Evaluate the Environmental Effect on Local Substrate

Modifications and upgrades to the hydraulic flume facility in the Environmental Fluid Mechanics and Hydraulics Laboratory (EFM&H) at Bucknell University are described. These changes enable small-scale testing of model marine hydrokinetic(MHK) devices. The design of the experimental platform provides a controlled environment for testing of model MHK devices to determine their effect on localsubstrate. Specifically, the effects being studied are scour and erosion around a cylindrical support structure and deposition of sediment downstream from the device.

Absence of the calcium-binding protein calretinin, not of calbindin D-28k, causes a permanent impairment of murine adult hippocampal neurogenesis

Calretinin (CR) and calbindin D-28k (CB) are cytosolic EF-hand Ca(2+)-binding proteins and function as Ca(2+) buffers affecting the spatiotemporal aspects of Ca(2+) transients and possibly also as Ca(2+) sensors modulating signaling cascades. In the adult hippocampal circuitry, CR and CB are expressed in specific principal neurons and subsets of interneurons. In addition, CR is transiently expressed within the neurogenic dentate gyrus (DG) niche. CR and CB expression during adult neurogenesis mark critical transition stages, onset of differentiation for CR, and the switch to adult-like connectivity for CB. Absence of either protein during these stages in null-mutant mice may have functional consequences and contribute to some aspects of the identified phenotypes. We report the impact of CR- and CB-deficiency on the proliferation and differentiation of progenitor cells within the subgranular zone (SGZ) neurogenic niche of the DG. Effects were evaluated (1) two and four weeks postnatally, during the transition period of the proliferative matrix to the adult state, and (2) in adult animals (3 months) to trace possible permanent changes in adult neurogenesis. The absence of CB from differentiated DG granule cells has no retrograde effect on the proliferative activity of progenitor cells, nor affects survival or migration/differentiation of newborn neurons in the adult DG including the SGZ. On the contrary, lack of CR from immature early postmitotic granule cells causes an early loss in proliferative capacity of the SGZ that is maintained into adult age, when it has a further impact on the migration/survival of newborn granule cells. The transient CR expression at the onset of adult neurogenesis differentiation may thus have two functions: (1) to serve as a self-maintenance signal for the pool of cells at the same stage of neurogenesis contributing to their survival/differentiation, and (2) it may contribute to retrograde signaling required for maintenance of the progenitor pool.

Agonist specific L-type Ca(2+)-current stimulation in ventricular myocytes by a novel steroid-like compound

In cardiac muscle the amplitude of Ca(2+) transients can be increased by enhancing Ca(2+) influx. Among the processes leading to increased Ca(2+) influx, agonists of the L-type Ca(2+)-channel can play an important role. Known pharmacological Ca(2+)-channel agonists act on different binding sites on the channel protein, which may lead not only to enhanced peak currents, but also to distinct changes in other biophysical characteristics of the current. In this study, membrane currents were recorded with the patch-clamp technique in the whole-cell configuration in guinea pig isolated ventricular myocytes in combination with confocal fluorescence Ca(2+) imaging techniques and a variety of pharmacological tools. Testing a new positive inotropic steroid-like compound, we found that it increased the L-type Ca(2+)-current by 2.5-fold by shifting the voltage-dependence of activation by 20.2 mV towards negative potentials. The dose-response relationship revealed two vastly different affinities (EC(50(high-affinity))=4.5+/-1.7 nM, EC(50(low-affinity))=8.0+/-1.1 microM) exhibiting differential pharmacological interactions with three classes of Ca(2+)-current antagonists, suggesting more than one binding site on the channel protein. Therefore, we identified and characterized a novel positive inotropic compound (F90927) as a member of a new class of Ca(2+)-channel agonists exhibiting unique features, which set it apart from other presently known L-type Ca(2+)-channel agonists.

Spine Ca2+ signaling in spike-timing-dependent plasticity

Calcium is a second messenger, which can trigger the modification of synaptic efficacy. We investigated the question of whether a differential rise in postsynaptic Ca2+ ([Ca2+]i) alone is sufficient to account for the induction of long-term potentiation (LTP) and long-term depression (LTD) of EPSPs in the basal dendrites of layer 2/3 pyramidal neurons of the somatosensory cortex. Volume-averaged [Ca2+]i transients were measured in spines of the basal dendritic arbor for spike-timing-dependent plasticity induction protocols. The rise in [Ca2+]i was uncorrelated to the direction of the change in synaptic efficacy, because several pairing protocols evoked similar spine [Ca2+]i transients but resulted in either LTP or LTD. The sequence dependence of near-coincident presynaptic and postsynaptic activity on the direction of changes in synaptic strength suggested that LTP and LTD were induced by two processes, which were controlled separately by postsynaptic [Ca2+]i levels. Activation of voltage-dependent Ca2+ channels before metabotropic glutamate receptors (mGluRs) resulted in the phospholipase C-dependent (PLC-dependent) synthesis of endocannabinoids, which acted as a retrograde messenger to induce LTD. LTP required a large [Ca2+]i transient evoked by NMDA receptor activation. Blocking mGluRs abolished the induction of LTD and uncovered the Ca2+-dependent induction of LTP. We conclude that the volume-averaged peak elevation of [Ca2+]i in spines of layer 2/3 pyramids determines the magnitude of long-term changes in synaptic efficacy. The direction of the change is controlled, however, via a mGluR-coupled signaling cascade. mGluRs act in conjunction with PLC as sequence-sensitive coincidence detectors when postsynaptic precede presynaptic action potentials to induce LTD. Thus presumably two different Ca2+ sensors in spines control the induction of spike-timing-dependent synaptic plasticity.

Tyrosine kinase modulation of protein kinase C activity regulates G protein-linked Ca2+ signaling in leukemic hematopoietic cells

We have used a recombinant mouse pre-B cell line (TonB210.1, expressing Bcr/Abl under the control of an inducible promoter) and several human leukemia cell lines to study the effect of high tyrosine kinase activity on G protein-coupled receptor (GPCR) agonist-stimulated cellular Ca(2+) release and store-operated Ca(2+) entry (SOCE). After induction of Bcr/Abl expression, GPCR-linked SOCE increased. The effect was reverted in the presence of the specific Abl inhibitor imatinib (1microM) and the Src inhibitor PP2 (10microM). In leukemic cell lines constitutively expressing high tyrosine kinase activity, Ca(2+) transients were reduced by imatinib and/or PP2. Ca(2+) transients were enhanced by specific inhibitors of PKC subtypes and this effect was amplified by tyrosine kinase inhibition in Bcr/Abl expressing TonB210.1 and K562 cells. Under all conditions Ca(2+) transients were essentially blocked by the PKC activator PMA. In Bcr/Abl expressing (but not in native) TonB210.1 cells, tyrosine kinase inhibitors enhanced PKCalpha catalytic activity and PKCalpha co-immunoprecipitated with Bcr/Abl. Unlike native TonB210.1 cells, Bcr/Abl expressing cells showed a high rate of cell death if Ca(2+) influx was reduced by complexing extracellular Ca(2+) with BAPTA. Our data suggest that tonic inhibition of PKC represents a mechanism by which high tyrosine kinase activity can enhance cellular Ca(2+) transients and thus exert profound effects on the proliferation, apoptosis and chemotaxis of leukemic cells.

Deficiency in parvalbumin, but not in calbindin D-28k upregulates mitochondrial volume and decreases smooth endoplasmic reticulum surface selectively in a peripheral, subplasmalemmal region in the soma of Purkinje cells

The Ca(2+)-binding proteins parvalbumin (PV) and calbindin D-28k (CB) are key players in the intracellular Ca(2+)-buffering in specific cells including neurons and have profound effects on spatiotemporal aspects of Ca(2+) transients. The previously observed increase in mitochondrial volume density in fast-twitch muscle of PV-/- mice is viewed as a specific compensation mechanism to maintain Ca(2+) homeostasis. Since cerebellar Purkinje cells (PC) are characterized by high expression levels of the Ca(2+) buffers PV and CB, the question was raised, whether homeostatic mechanisms are induced in PC lacking these buffers. Mitochondrial volume density, i.e. relative mitochondrial mass was increased by 40% in the soma of PV-/- PC. Upregulation of mitochondrial volume density was not homogenous throughout the soma, but was selectively restricted to a peripheral region of 1.5 microm width underneath the plasma membrane. Accompanied was a decreased surface of subplasmalemmal smooth endoplasmic reticulum (sPL-sER) in a shell of 0.5 microm thickness underneath the plasma membrane. These alterations were specific for the absence of the "slow-onset" buffer PV, since in CB-/- mice neither changes in peripheral mitochondria nor in sPL-sER were observed. This implicates that the morphological alterations are aimed to specifically substitute the function of the slow buffer PV. We propose a novel concept that homeostatic mechanisms of components involved in Ca(2+) homeostasis do not always occur at the level of similar or closely related molecules. Rather the cell attempts to restore spatiotemporal aspects of Ca(2+) signals prevailing in the undisturbed (wildtype) situation by subtly fine tuning existing components involved in the regulation of Ca(2+) fluxes.

Calcium indicator loading of neurons using single-cell electroporation

Studies of subcellular Ca(2+) signaling rely on methods for labeling cells with fluorescent Ca(2+) indicator dyes. In this study, we demonstrate the use of single-cell electroporation for Ca(2+) indicator loading of individual neurons and small neuronal networks in rat neocortex in vitro and in vivo. Brief voltage pulses were delivered through glass pipettes positioned close to target cells. This approach resulted in reliable and rapid (within seconds) loading of somata and subsequent complete labeling of dendritic and axonal arborizations. By using simultaneous whole-cell recordings in brain slices, we directly addressed the effect of electroporation on neurons. Cell viability was high (about 85%) with recovery from the membrane permeabilization occurring within a minute. Electrical properties of recovered cells were indistinguishable before and after electroporation. In addition, Ca(2+) transients with normal appearance could be evoked in dendrites, spines, and axonal boutons of electroporated cells. Using negative-stains of somata, targeted single-cell electroporation was equally applicable in vivo. We conclude that electroporation is a simple approach that permits Ca(2+) indicator loading of multiple cells with low background staining within a short amount of time, which makes it especially well suited for functional imaging of subcellular Ca(2+) dynamics in small neuronal networks.

Dendritic spikes in apical dendrites of neocortical layer 2/3 pyramidal neurons

Layer 2/3 (L2/3) pyramidal neurons are the most abundant cells of the neocortex. Despite their key position in the cortical microcircuit, synaptic integration in dendrites of L2/3 neurons is far less understood than in L5 pyramidal cell dendrites, mainly because of the difficulties in obtaining electrical recordings from thin dendrites. Here we directly measured passive and active properties of the apical dendrites of L2/3 neurons in rat brain slices using dual dendritic-somatic patch-clamp recordings and calcium imaging. Unlike L5 cells, L2/3 dendrites displayed little sag in response to long current pulses, which suggests a low density of I(h) in the dendrites and soma. This was also consistent with a slight increase in input resistance with distance from the soma. Brief current injections into the apical dendrite evoked relatively short (half-width 2-4 ms) dendritic spikes that were isolated from the soma for near-threshold currents at sites beyond the middle of the apical dendrite. Regenerative dendritic potentials and large concomitant calcium transients were also elicited by trains of somatic action potentials (APs) above a critical frequency (130 Hz), which was slightly higher than in L5 neurons. Initiation of dendritic spikes was facilitated by backpropagating somatic APs and could cause an additional AP at the soma. As in L5 neurons, we found that distal dendritic calcium transients are sensitive to a long-lasting block by GABAergic inhibition. We conclude that L2/3 pyramidal neurons can generate dendritic spikes, sharing with L5 pyramidal neurons fundamental properties of dendritic excitability and control by inhibition.

Visual and quantitative electroencephalographic analysis of healthy young and adult cats under medetomidine sedation

A study was designed to investigate the effect of medetomidine sedation on quantitative electroencephalography (q-EEG) in healthy young and adult cats to determine objective guidelines for diagnostic EEG recordings and interpretation. Preliminary visual examination of EEG recordings revealed high-voltage low-frequency background activity. Spindles, k-complexes and vertex sharp transients characteristic of sleep or sedation were superimposed on a low background activity. Neither paroxysmal activity nor EEG burst-suppression were observed. The spectral analysis of q-EEG included four parameters, namely, relative power (%), and mean, median and peak frequency (Hz) of all four frequency bands (delta, theta, alpha and beta). The findings showed a prevalence of slow delta and theta rhythms as opposed to fast alpha and beta rhythms in both young (group A) and adult (group B) cats. A posterior gradient was reported for the theta band and an anterior gradient for the alpha and beta bands in both groups, respectively. The relative power value in group B compared to group A was significantly higher for theta, alpha and beta bands, and lower for the delta band. The mean and median frequency values in group B was significantly higher for delta, theta and beta bands and lower for the alpha band. The study has shown that a medetomidine sedation protocol for feline EEG may offer a method for investigating bio-electrical cortical activity. The use of q-EEG analysis showed a decrease in high frequency bands and increased activity of the low frequency band in healthy cats under medetomidine sedation.

Apparent survival of the salamander

Background Understanding the demographic processes underlying population dynamics is a central theme in ecology. Populations decline if losses from the population (i.e., mortality and emigration) exceed gains (i.e., recruitment and immigration). Amphibians are thought to exhibit little movement even though local populations often fluctuate dramatically and are likely to go exinct if there is no rescue effect through immigration from nearby populations. Terrestrial salamanders are generally portrayed as amphibians with low migratory activity. Our study uses demographic analysis as a key to unravel whether emigration or mortality is the main cause of "losses" from the population. In particular, we use the analysis to challenge the common belief that terrestrial salamanders show low migratory activity. Results The mark-recapture analysis of adult salamanders showed that monthly survival was high (> 90%) without a seasonal pattern. These estimates, however, translate into rather low rates of local annual survival of only ~40% and suggest that emigration was important. The estimated probability of emigration was 49%. Conclusion Our analysis shows that terrestrial salamanders exhibit more migratory activity than commonly thought. This may be due either because the spatial extent of salamander populations is underestimated or because there is a substantial exchange of individuals between populations. Our current results are in line with several other studies that suggest high migratory activity in amphibians. In particular, many amphibian populations may be characterized by high proportions of transients and/or floaters.

Transport model of the human Na+-coupled L-ascorbic acid (vitamin C) transporter SVCT1

Vitamin C (L-ascorbic acid) is an essential micronutrient that serves as an antioxidant and as a cofactor in many enzymatic reactions. Intestinal absorption and renal reabsorption of the vitamin is mediated by the epithelial apical L-ascorbic acid cotransporter SVCT1 (SLC23A1). We explored the molecular mechanisms of SVCT1-mediated L-ascorbic acid transport using radiotracer and voltage-clamp techniques in RNA-injected Xenopus oocytes. L-ascorbic acid transport was saturable (K(0.5) approximately 70 microM), temperature dependent (Q(10) approximately 5), and energized by the Na(+) electrochemical potential gradient. We obtained a Na(+)-L-ascorbic acid coupling ratio of 2:1 from simultaneous measurement of currents and fluxes. L-ascorbic acid and Na(+) saturation kinetics as a function of cosubstrate concentrations revealed a simultaneous transport mechanism in which binding is ordered Na(+), L-ascorbic acid, Na(+). In the absence of L-ascorbic acid, SVCT1 mediated pre-steady-state currents that decayed with time constants 3-15 ms. Transients were described by single Boltzmann distributions. At 100 mM Na(+), maximal charge translocation (Q(max)) was approximately 25 nC, around a midpoint (V(0.5)) at -9 mV, and with apparent valence approximately -1. Q(max) was conserved upon progressive removal of Na(+), whereas V(0.5) shifted to more hyperpolarized potentials. Model simulation predicted that the pre-steady-state current predominantly results from an ion-well effect on binding of the first Na(+) partway within the membrane electric field. We present a transport model for SVCT1 that will provide a framework for investigating the impact of specific mutations and polymorphisms in SLC23A1 and help us better understand the contribution of SVCT1 to vitamin C metabolism in health and disease.

ATP modeling of internal transformer faults for relay performance testing

Transformers are very important elements of any power system. Unfortunately, they are subjected to through-faults and abnormal operating conditions which can affect not only the transformer itself but also other equipment connected to the transformer. Thus, it is essential to provide sufficient protection for transformers as well as the best possible selectivity and sensitivity of the protection. Nowadays microprocessor-based relays are widely used to protect power equipment. Current differential and voltage protection strategies are used in transformer protection applications and provide fast and sensitive multi-level protection and monitoring. The elements responsible for detecting turn-to-turn and turn-to-ground faults are the negative-sequence percentage differential element and restricted earth-fault (REF) element, respectively. During severe internal faults current transformers can saturate and slow down the speed of relay operation which affects the degree of equipment damage. The scope of this work is to develop a modeling methodology to perform simulations and laboratory tests for internal faults such as turn-to-turn and turn-to-ground for two step-down power transformers with capacity ratings of 11.2 MVA and 290 MVA. The simulated current waveforms are injected to a microprocessor relay to check its sensitivity for these internal faults. Saturation of current transformers is also studied in this work. All simulations are performed with the Alternative Transients Program (ATP) utilizing the internal fault model for three-phase two-winding transformers. The tested microprocessor relay is the SEL-487E current differential and voltage protection relay. The results showed that the ATP internal fault model can be used for testing microprocessor relays for any percentage of turns involved in an internal fault. An interesting observation from the experiments was that the SEL-487E relay is more sensitive to turn-to-turn faults than advertized for the transformers studied. The sensitivity of the restricted earth-fault element was confirmed. CT saturation cases showed that low accuracy CTs can be saturated with a high percentage of turn-to-turn faults, where the CT burden will affect the extent of saturation. Recommendations for future work include more accurate simulation of internal faults, transformer energization inrush, and other scenarios involving core saturation, using the newest version of the internal fault model. The SEL-487E relay or other microprocessor relays should again be tested for performance. Also, application of a grounding bank to the delta-connected side of a transformer will increase the zone of protection and relay performance can be tested for internal ground faults on both sides of a transformer.

Time-domain models for power system stability and unbalance

It is an important and difficult challenge to protect modern interconnected power system from blackouts. Applying advanced power system protection techniques and increasing power system stability are ways to improve the reliability and security of power systems. Phasor-domain software packages such as Power System Simulator for Engineers (PSS/E) can be used to study large power systems but cannot be used for transient analysis. In order to observe both power system stability and transient behavior of the system during disturbances, modeling has to be done in the time-domain. This work focuses on modeling of power systems and various control systems in the Alternative Transients Program (ATP). ATP is a time-domain power system modeling software in which all the power system components can be modeled in detail. Models are implemented with attention to component representation and parameters. The synchronous machine model includes the saturation characteristics and control interface. Transient Analysis Control System is used to model the excitation control system, power system stabilizer and the turbine governor system of the synchronous machine. Several base cases of a single machine system are modeled and benchmarked against PSS/E. A two area system is modeled and inter-area and intra-area oscillations are observed. The two area system is reduced to a two machine system using reduced dynamic equivalencing. The original and the reduced systems are benchmarked against PSS/E. This work also includes the simulation of single-pole tripping using one of the base case models. Advantages of single-pole tripping and comparison of system behavior against three-pole tripping are studied. Results indicate that the built-in control system models in PSS/E can be effectively reproduced in ATP. The benchmarked models correctly simulate the power system dynamics. The successful implementation of a dynamically reduced system in ATP shows promise for studying a small sub-system of a large system without losing the dynamic behaviors. Other aspects such as relaying can be investigated using the benchmarked models. It is expected that this work will provide guidance in modeling different control systems for the synchronous machine and in representing dynamic equivalents of large power systems.

Parameter estimation for transformer modeling

Large Power transformers, an aging and vulnerable part of our energy infrastructure, are at choke points in the grid and are key to reliability and security. Damage or destruction due to vandalism, misoperation, or other unexpected events is of great concern, given replacement costs upward of $2M and lead time of 12 months. Transient overvoltages can cause great damage and there is much interest in improving computer simulation models to correctly predict and avoid the consequences. EMTP (the Electromagnetic Transients Program) has been developed for computer simulation of power system transients. Component models for most equipment have been developed and benchmarked. Power transformers would appear to be simple. However, due to their nonlinear and frequency-dependent behaviors, they can be one of the most complex system components to model. It is imperative that the applied models be appropriate for the range of frequencies and excitation levels that the system experiences. Thus, transformer modeling is not a mature field and newer improved models must be made available. In this work, improved topologically-correct duality-based models are developed for three-phase autotransformers having five-legged, three-legged, and shell-form cores. The main problem in the implementation of detailed models is the lack of complete and reliable data, as no international standard suggests how to measure and calculate parameters. Therefore, parameter estimation methods are developed here to determine the parameters of a given model in cases where available information is incomplete. The transformer nameplate data is required and relative physical dimensions of the core are estimated. The models include a separate representation of each segment of the core, including hysteresis of the core, λ-i saturation characteristic, capacitive effects, and frequency dependency of winding resistance and core loss. Steady-state excitation, and de-energization and re-energization transients are simulated and compared with an earlier-developed BCTRAN-based model. Black start energization cases are also simulated as a means of model evaluation and compared with actual event records. The simulated results using the model developed here are reasonable and more correct than those of the BCTRAN-based model. Simulation accuracy is dependent on the accuracy of the equipment model and its parameters. This work is significant in that it advances existing parameter estimation methods in cases where the available data and measurements are incomplete. The accuracy of EMTP simulation for power systems including three-phase autotransformers is thus enhanced. Theoretical results obtained from this work provide a sound foundation for development of transformer parameter estimation methods using engineering optimization. In addition, it should be possible to refine which information and measurement data are necessary for complete duality-based transformer models. To further refine and develop the models and transformer parameter estimation methods developed here, iterative full-scale laboratory tests using high-voltage and high-power three-phase transformer would be helpful.