Numerical solutions of elliptic inverse problems via the equation error method

To estimate a parameter in an elliptic boundary value problem, the method of equation error chooses the value that minimizes the error in the PDE and boundary condition (the solution of the BVP having been replaced by a measurement). The estimated parameter converges to the exact value as the measured data converge to the exact value, provided Tikhonov regularization is used to control the instability inherent in the problem. The error in the estimated solution can be bounded in an appropriate quotient norm; estimates can be derived for both the underlying (infinite-dimensional) problem and a finite-element discretization that can be implemented in a practical algorithm. Numerical experiments demonstrate the efficacy and limitations of the method.

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.

Automation of shear-wave splitting parameter determination of local earthquakes at Yellowstone : application as indicator of crustal stress and temporal variation

Shear-wave splitting can be a useful technique for determining crustal stress fields in volcanic settings and temporal variations associated with activity. Splitting parameters were determined for a subset of local earthquakes recorded from 2000-2010 at Yellowstone. Analysis was automated using an unsupervised cluster analysis technique to determine optimum splitting parameters from 270 analysis windows for each event. Six stations clearly exhibit preferential fast polarization values sub-orthogonal to the direction of minimum horizontal compression. Yellowstone deformation results in a local crustal stress field differing from the regional field dominated by NE-SW extension, and fast directions reflect this difference rotating around the caldera maintaining perpendicularity to the rim. One station exhibits temporal variations concordant with identified periods of caldera subsidence and uplift. From splitting measurements, we calculated a crustal anisotropy of ~17-23% and crack density ~0.12-0.17 possibly resulting from stress-aligned fluid filled microcracks in the upper crust and an active hydrothermal system.

Quantum transport in a single molecular junction

The craze for faster and smaller electronic devices has never gone down and this has always kept researchers on their toes. Following Moore’s law, which states that the number of transistors in a single chip will double in every 18 months, today “30 million transistors can fit into the head of a 1.5 mm diameter pin”. But this miniaturization cannot continue indefinitely due to the ‘quantum leakage’ limit in the thickness of the insulating layer between the gate electrode and the current carrying channel. To bypass this limitation, scientists came up with the idea of using vastly available organic molecules as components in an electronic device. One of the primary challenges in this field was the ability to perform conductance measurements across single molecular junctions. Once that was achieved the focus shifted to a deeper understanding of the underlying physics behind the electron transport across these molecular scale devices. Our initial theoretical approach is based on the conventional Non-Equilibrium Green Function(NEGF) formulation, but the self-energy of the leads is modified to include a weighting factor that ensures negligible current in the absence of a molecular pathway as observed in a Mechanically Controlled Break Junction (MCBJ) experiment. The formulation is then made parameter free by a more careful estimation of the self-energy of the leads. The calculated conductance turns out to be atleast an order more than the experimental values which is probably due to a strong chemical bond at the metal-molecule junction unlike in the experiments. The focus is then shifted to a comparative study of charge transport in molecular wires of different lengths within the same formalism. The molecular wires, composed of a series of organic molecules, are sanwiched between two gold electrodes to make a two terminal device. The length of the wire is increased by sequentially increasing the number of molecules in the wire from 1 to 3. In the low bias regime all the molecular devices are found to exhibit Ohmic behavior. However, the magnitude of conductance decreases exponentially with increase in length of the wire. In the next study, the relative contribution of the ‘in-phase’ and the ‘out-of-phase’ components of the total electronic current under the influence of an external bias is estimated for the wires of three different lengths. In the low bias regime, the ‘out-of-phase’ contribution to the total current is minimal and the ‘in-phase’ elastic tunneling of the electrons is responsible for the net electronic current. This is true irrespective of the length of the molecular spacer. In this regime, the current-voltage characteristics follow Ohm’s law and the conductance of the wires is found to decrease exponentially with increase in length which is in agreement with experimental results. However, after a certain ‘off-set’ voltage, the current increases non-linearly with bias and the ‘out-of-phase’ tunneling of electrons reduces the net current substantially. Subsequently, the interaction of conduction electrons with the vibrational modes as a function of external bias in the three different oligomers is studied since they are one of the main sources of phase-breaking scattering. The number of vibrational modes that couple strongly with the frontier molecular orbitals are found to increase with length of the spacer and the external field. This is consistent with the existence of lowest ‘off-set’ voltage for the longest wire under study.

Parameter identification of a copper-zeolite SCR catalyst model using reactor data

The selective catalytic reduction system is a well established technology for NOx emissions control in diesel engines. A one dimensional, single channel selective catalytic reduction (SCR) model was previously developed using Oak Ridge National Laboratory (ORNL) generated reactor data for an iron-zeolite catalyst system. Calibration of this model to fit the experimental reactor data collected at ORNL for a copper-zeolite SCR catalyst is presented. Initially a test protocol was developed in order to investigate the different phenomena responsible for the SCR system response. A SCR model with two distinct types of storage sites was used. The calibration process was started with storage capacity calculations for the catalyst sample. Then the chemical kinetics occurring at each segment of the protocol was investigated. The reactions included in this model were adsorption, desorption, standard SCR, fast SCR, slow SCR, NH3 Oxidation, NO oxidation and N2O formation. The reaction rates were identified for each temperature using a time domain optimization approach. Assuming an Arrhenius form of the reaction rates, activation energies and pre-exponential parameters were fit to the reaction rates. The results indicate that the Arrhenius form is appropriate and the reaction scheme used allows the model to fit to the experimental data and also for use in real world engine studies.