DISTRIBUTED APPROACH FOR SOLVING OPTIMAL POWER FLOW PROBLEMS IN THREE-PHASE UNBALANCED DISTRIBUTION NETWORK

The objective of this report is to study distributed (decentralized) three phase optimal power flow (OPF) problem in unbalanced power distribution networks. A full three phase representation of the distribution networks is considered to account for the highly unbalance state of the distribution networks. All distribution network’s series/shunt components, and load types/combinations had been modeled on commercial version of General Algebraic Modeling System (GAMS), the high-level modeling system for mathematical programming and optimization. The OPF problem has been successfully implemented and solved in a centralized approach and distributed approach, where the objective is to minimize the active power losses in the entire system. The study was implemented on the IEEE-37 Node Test Feeder. A detailed discussion of all problem sides and aspects starting from the basics has been provided in this study. Full simulation results have been provided at the end of the report.

ACETYL RADICAL IN TOBACCO SMOKE: DETECTION, QUANTIFICATION AND SIMULATION

Free radicals are present in cigarette smoke and can have a negative effect on human health by attacking lipids, nucleic acids, proteins and other biologically important species. However, because of the complexity of the tobacco smoke system and the dynamic nature of radicals, little is known about the identity of the radicals, and debate continues on the mechanisms by which those radicals are produced. In this study, acetyl radicals were trapped from the gas phase using 3-amino-2, 2, 5, 5- tetramethyl-proxyl (3AP) on solid support to form stable 3AP adducts for later analysis by high performance liquid chromatography (HPLC), mass spectrometry/tandem mass spectrometry (MS-MS/MS) and liquid chromatography- mass spectrometry (LC-MS). Simulations of acetyl radical generation were performed using Matlab and the Master Chemical Mechanism (MCM) programs. A range of 10- 150 nmol/cigarette of acetyl radical was measured from gas phase tobacco smoke of both commerial and research cigarettes under several different smoking conditions. More radicals were detected from the puff smoking method compared to continuous flow sampling. Approximately twice as many acetyl radicals were trapped when a GF/F particle filter was placed before the trapping zone. Computational simulations show that NO/NO2 reacts with isoprene, initiating chain reactions to produce a hydroxyl radical, which abstracts hydrogen from acetaldehyde to generate acetyl radical. With initial concentrations of NO, acetaldehyde, and isoprene in a real-world cigarette smoke scenario, these mechanisms can account for the full amount of acetyl radical detected experimentally. This study contributes to the overall understanding of the free radical generation in gas phase cigarette smoke.

Assessments and Computational Simulations in Support of a Time-Varying Mass Flow Rate Measurement Technique for Pulsatile Gas Flow

This thesis covers the correction, and verification, development, and implementation of a computational fluid dynamics (CFD) model for an orifice plate meter. Past results were corrected and further expanded on with compressibility effects of acoustic waves being taken into account. One dynamic pressure difference transducer measures the time-varying differential pressure across the orifice meter. A dynamic absolute pressure measurement is also taken at the inlet of the orifice meter, along with a suitable temperature measurement of the mean flow gas. Together these three measurements allow for an incompressible CFD simulation (using a well-tested and robust model) for the cross-section independent time-varying mass flow rate through the orifice meter. The mean value of this incompressible mass flow rate is then corrected to match the mean of the measured flow rate( obtained from a Coriolis meter located up stream of the orifice meter). Even with the mean and compressibility corrections, significant differences in the measured mass flow rates at two orifice meters in a common flow stream were observed. This means that the compressibility effects associated with pulsatile gas flows is significant in the measurement of the time-varying mass flow rate. Future work (with the approach and initial runs covered here) will provide an indirect verification of the reported mass flow rate measurements.

MKS-Simulation eines Staplers mit aktivem Fahrwerk

In diesem Beitrag wird das Mehrkörpersystem-Modell (MKS-Modell) eines Gegengewichtsstap-lers mit aktivem Fahrwerk erläutert. Die wesentlichen Merkmale des MKS-Modells sind die neuentwi-ckelten Achsaufhängungen und die Hubgerüstlagerung des auf einen Gegengewichtsstapler mit 3 t - Tragkraft basierenden Fahrzeuges. Ein weiterer Aspekt ist die Kopplung der mechatronischen Modelltei-le Hydraulik und Regelung mit dem mechanischen Staplermodell. Abschließend zeigen die Ergebnisse der Co-Simulationen Verbesserungen der Fahrsicherheit und des Fahrkomforts.

Erforderlicher technischer Durchsatz von Materialfusssystemen zur Erreichung eines Solldurchsatzes - ermittelt durch Simulation mit Integrationsalgorithmus

Eine Analyse der Kundenauftragsprozesse in der Automobilindustrie zeigt, dass trotz großer An-strengungen nach wie vor Schwierigkeiten bei der Stabilisierung der komplexen Produktions- und Logistikprozesse mit negativen Auswirkungen auf die Zielgröße Termintreue auftreten. Innerhalb des mehrstufigen Wertschöpfungsnetzwerks wird dieser Planungs¬unsicherheit durch kurzfristige Bedarfs-schwankungen mit großen Sicherheits¬beständen und der Forderung nach hochflexiblen Logistik- und Produktionsprozessen bei den Zulieferern begegnet. In diesem Beitrag werden neue Methoden und Konzepte zur Planung und Umsetzung der Produktionssteuerung nach stabiler Auftragsfolge in der Automobilindustrie vorgestellt. Ziel ist es, den Produktionsprozess hinreichend stabilisieren zu können, sodass erhebliches Rationalisierungspotential sowohl bei den internen Prozessen als auch im gesamten Zuliefernetzwerk erreicht wird.

Setting Inventory Levels of CONWIP Flow Lines via Linear Programming

This paper treats the problem of setting the inventory level and optimizing the buffer allocation of closed-loop flow lines operating under the constant-work-in-process (CONWIP) protocol. We solve a very large but simple linear program that models an entire simulation run of a closed-loop flow line in discrete time to determine a production rate estimate of the system. This approach introduced in Helber, Schimmelpfeng, Stolletz, and Lagershausen (2011) for open flow lines with limited buffer capacities is extended to closed-loop CONWIP flow lines. Via this method, both the CONWIP level and the buffer allocation can be optimized simultaneously. The first part of a numerical study deals with the accuracy of the method. In the second part, we focus on the relationship between the CONWIP inventory level and the short-term profit. The accuracy of the method turns out to be best for such configurations that maximize production rate and/or short-term profit.

Hemodynamics and Flow Characteristics of a New Dialysis Port

Renal replacement therapy by hemodialysis requires a permanent vascular access. Implantable ports offer a potential alternative to standard vascular access strategies although their development is limited both in number and extent. We explored the fluid dynamics within two new percutaneous bone-anchored dialysis port prototypes, both by in vitro experiments and computer simulation. The new port is to be fixed to bone and allows the connection of a dialysis machine to a central venous catheter via a built-in valve. We found that the pressure drop induced by the two ports was between 20 and 50 mmHg at 500 ml/min, which is comparable with commercial catheter connectors (15–80 mmHg). We observed the formation of vortices in both geometries, and a shear rate in the physiological range (<10,000s-1), which is lower than maximal shear rates reported in commercial catheters (up to 13,000s-1). A difference in surface shear rate of 15% between the two ports was obtained.

Comparing three CPR feedback devices and standard BLS in a single rescuer scenario: a randomised simulation study

BACKGROUND Efficiently performed basic life support (BLS) after cardiac arrest is proven to be effective. However, cardiopulmonary resuscitation (CPR) is strenuous and rescuers' performance declines rapidly over time. Audio-visual feedback devices reporting CPR quality may prevent this decline. We aimed to investigate the effect of various CPR feedback devices on CPR quality. METHODS In this open, prospective, randomised, controlled trial we compared three CPR feedback devices (PocketCPR, CPRmeter, iPhone app PocketCPR) with standard BLS without feedback in a simulated scenario. 240 trained medical students performed single rescuer BLS on a manikin for 8min. Effective compression (compressions with correct depth, pressure point and sufficient decompression) as well as compression rate, flow time fraction and ventilation parameters were compared between the four groups. RESULTS Study participants using the PocketCPR performed 17±19% effective compressions compared to 32±28% with CPRmeter, 25±27% with the iPhone app PocketCPR, and 35±30% applying standard BLS (PocketCPR vs. CPRmeter p=0.007, PocketCPR vs. standard BLS p=0.001, others: ns). PocketCPR and CPRmeter prevented a decline in effective compression over time, but overall performance in the PocketCPR group was considerably inferior to standard BLS. Compression depth and rate were within the range recommended in the guidelines in all groups. CONCLUSION While we found differences between the investigated CPR feedback devices, overall BLS quality was suboptimal in all groups. Surprisingly, effective compression was not improved by any CPR feedback device compared to standard BLS. All feedback devices caused substantial delay in starting CPR, which may worsen outcome.

Real-time simulation of nonequilibrium transport of magnetization in large open quantum spin systems driven by dissipation

Using quantum Monte Carlo, we study the nonequilibrium transport of magnetization in large open strongly correlated quantum spin-12 systems driven by purely dissipative processes that conserve the uniform or staggered magnetization, disregarding unitary Hamiltonian dynamics. We prepare both a low-temperature Heisenberg ferromagnet and an antiferromagnet in two parts of the system that are initially isolated from each other. We then bring the two subsystems in contact and study their real-time dissipative dynamics for different geometries. The flow of the uniform or staggered magnetization from one part of the system to the other is described by a diffusion equation that can be derived analytically.

Computer simulation of electrophoretic aspects of enantiomer migration and separation in capillary electrochromatography with a neutral selector.

A computer simulation study describing the electrophoretic separation and migration of methadone enantiomers in presence of free and immobilized (2-hydroxypropyl)-β-CD is presented. The 1:1 interaction of methadone with the neutral CD was simulated by using experimentally determined mobilities and complexation constants for the complexes in a low-pH BGE comprising phosphoric acid and KOH. The use of complex mobilities represents free solution conditions with the chiral selector being a buffer additive, whereas complex mobilities set to zero provide data that mimic migration and separation with the chiral selector being immobilized, that is CEC conditions in absence of unspecific interaction between analytes and the chiral stationary phase. Simulation data reveal that separations are quicker, electrophoretic displacement rates are reduced, and sensitivity is enhanced in CEC with on-column detection in comparison to free solution conditions. Simulation is used to study electrophoretic analyte behavior at the interface between sample and the CEC column with the chiral selector (analyte stacking) and at the rear end when analytes leave the environment with complexation (analyte destacking). The latter aspect is relevant for off-column analyte detection in CEC and is described here for the first time via the dynamics of migrating analyte zones. Simulation provides insight into means to counteract analyte dilution at the column end via use of a BGE with higher conductivity. Furthermore, the impact of EOF on analyte migration, separation, and detection for configurations with the selector zone being displaced or remaining immobilized under buffer flow is simulated. In all cases, the data reveal that detection should occur within or immediately after the selector zone.

Impact of Taylor-Aris diffusivity on analyte and system zone dispersion in CZE assessed by computer simulation and experimental validation.

Application of pressure-driven laminar flow has an impact on zone and boundary dispersion in open tubular CE. The GENTRANS dynamic simulator for electrophoresis was extended with Taylor-Aris diffusivity which accounts for dispersion due to the parabolic flow profile associated with pressure-driven flow. Effective diffusivity of analyte and system zones as functions of the capillary diameter and the amount of flow in comparison to molecular diffusion alone were studied for configurations with concomitant action of imposed hydrodynamic flow and electroosmosis. For selected examples under realistic experimental conditions, simulation data are compared with those monitored experimentally using modular CE setups featuring both capacitively coupled contactless conductivity and UV absorbance detection along a 50 μm id fused-silica capillary of 90 cm total length. The data presented indicate that inclusion of flow profile based Taylor-Aris diffusivity provides realistic simulation data for analyte and system peaks, particularly those monitored in CE with conductivity detection.

Interstellar neutral helium in the heliosphere from Interstellar Boundary Explorer observations. III. Mach number of the flow, velocity vector, and temperature from the first six years of measurements

We analyzed observations of interstellar neutral helium (ISN He) obtained from the Interstellar Boundary Explorer (IBEX) satellite during its first six years of operation. We used a refined version of the ISN He simulation model, presented in the companion paper by Sokol et al. (2015b), along with a sophisticated data correlation and uncertainty system and parameter fitting method, described in the companion paper by Swaczyna et al. We analyzed the entire data set together and the yearly subsets, and found the temperature and velocity vector of ISN He in front of the heliosphere. As seen in the previous studies, the allowable parameters are highly correlated and form a four-dimensional tube in the parameter space. The inflow longitudes obtained from the yearly data subsets show a spread of similar to 6 degrees, with the other parameters varying accordingly along the parameter tube, and the minimum chi(2) value is larger than expected. We found, however, that the Mach number of the ISN He flow shows very little scatter and is thus very tightly constrained. It is in excellent agreement with the original analysis of ISN He observations from IBEX and recent reanalyses of observations from Ulysses. We identify a possible inaccuracy in the Warm Breeze parameters as the likely cause of the scatter in the ISN He parameters obtained from the yearly subsets, and we suppose that another component may exist in the signal or a process that is not accounted for in the current physical model of ISN He in front of the heliosphere. From our analysis, the inflow velocity vector, temperature, and Mach number of the flow are equal to lambda(ISNHe) = 255 degrees.8 +/- 0 degrees.5, beta(ISNHe) = 5 degrees.16 +/- 0 degrees.10, T-ISNHe = 7440 +/- 260 K, nu(SNHe) = 25.8 +/- 0.4 km s(-1), and M-ISNHe = 5.079 +/- 0.028, with uncertainties strongly correlated along the parameter tube.

Flow separation and roughness lengths over large bedforms in a tidal environment

This study characterises the shape of the flow separation zone (FSZ) and wake region over large asymmetric bedforms under tidal flow conditions. High resolution bathymetry, flow velocity and turbulence data were measured along two parallel transects in a tidal channel covered with bedforms. The field data are used to verify the applicability of a numerical model for a systematic study using the Delft3D modelling system and test the model sensitivity to roughness length. Three experiments are then conducted to investigate how the FSZ size and wake extent vary depending on tidally-varying flow conditions, water levels and bathymetry. During the ebb, a large FSZ occurs over the steep lee side of each bedform. During the flood, no flow separation develops over the bedforms having a flat crest; however, a small FSZ is observed over the steepest part of the crest of some bedforms, where the slope is locally up to 15°. Over a given bedform morphology and constant water levels, no FSZ occurs for velocity magnitudes smaller than 0.1 m s**-1; as the flow accelerates, the FSZ reaches a stable size for velocity magnitudes greater than 0.4 m s**-1. The shape of the FSZ is not influenced by changes in water levels. On the other hand, variations in bed morphology, as recorded from the high-resolution bathymetry collected during the tidal cycle, influence the size and position of the FSZ: a FSZ develops only when the maximum lee side slope over a horizontal distance of 5 m is greater than 10°. The height and length of the wake region are related to the length of the FSZ. The total roughness along the transect lines is an order of magnitude larger during the ebb than during the flood due to flow direction in relation to bedform asymmetry: during the ebb, roughness is created by the large bedforms because a FSZ and wake develops over the steep lee side. The results add to the understanding of hydrodynamics of natural bedforms in a tidal environment and may be used to better parameterise small-scale processes in large-scale studies.

3D Direct Simulation Monte Carlo Modelling of the Inner Gas Coma of Comet 67P/Churyumov-Gerasimenko: A Parameter Study

Direct Simulation Monte Carlo (DSMC) is a powerful numerical method to study rarefied gas flows such as cometary comae and has been used by several authors over the past decade to study cometary outflow. However, the investigation of the parameter space in simulations can be time consuming since 3D DSMC is computationally highly intensive. For the target of ESA's Rosetta mission, comet 67P/Churyumov-Gerasimenko, we have identified to what extent modification of several parameters influence the 3D flow and gas temperature fields and have attempted to establish the reliability of inferences about the initial conditions from in situ and remote sensing measurements. A large number of DSMC runs have been completed with varying input parameters. In this work, we present the simulation results and conclude on the sensitivity of solutions to certain inputs. It is found that among cases of water outgassing, the surface production rate distribution is the most influential variable to the flow field.

Numerical simulation of the upward propagation of a flame in a vertical tube filled with a very lean mixture.

Upwardpropagation of a premixed flame in averticaltubefilled with a very leanmixture is simulated numerically using a single irreversible Arrhenius reaction model with infinitely high activation energy. In the absence of heat losses and preferential diffusion effects, a curved flame with stationary shape and velocity close to those of an open bubble ascending in the same tube is found for values of the fuel mass fraction above a certain minimum that increases with the radius of the tube, while the numerical computations cease to converge to a stationary solution below this minimum mass fraction. The vortical flow of the gas behind the flame and in its transport region is described for tubes of different radii. It is argued that this flow may become unstable when the fuel mass fraction is decreased, and that this instability, together with the flame stretch due to the strong curvature of the flame tip in narrow tubes, may be responsible for the minimum fuel mass fraction. Radiation losses and a Lewis number of the fuel slightly above unity decrease the final combustion temperature at the flame tip and increase the minimum fuel mass fraction, while a Lewis number slightly below unity has the opposite effect.