Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulated current”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 973 records · Page 54

The potential role of methanesulfonic acid (MSA) in aerosol formation and growth and the associated radiative forcings

Atmospheric marine aerosol particles impact Earth's albedo and climate. These particles can be primary or secondary and come from a variety of sources, including sea salt, dissolved organic matter, volatile organic compounds, and sulfur-containing compounds. Dimethylsulfide (DMS) marine emissions contribute greatly to the global biogenic sulfur budget, and its oxidation products can contribute to aerosol mass, specifically as sulfuric acid and methanesulfonic acid (MSA). Further, sulfuric acid is a known nucleating compound, and MSA may be able to participate in nucleation when bases are available. As DMS emissions, and thus MSA and sulfuric acid from DMS oxidation, may have changed since pre-industrial times and may change in a warming climate, it is important to characterize and constrain the climate impacts of both species. Currently, global models that simulate aerosol size distributions include contributions of sulfate and sulfuric acid from DMS oxidation, but to our knowledge, global models typically neglect the impact of MSA on size distributions. In this study, we use the GEOS-Chem-TOMAS (GC-TOMAS) global aerosol microphysics model to determine the impact on aerosol size distributions and subsequent aerosol radiative effects from including MSA in the size-resolved portion of the model. The effective equilibrium vapor pressure of MSA is currently uncertain, and we use the Extended Aerosol Inorganics Model (E-AIM) to build a parameterization for GC-TOMAS of MSA's effective volatility as a function of temperature, relative humidity, and available gas-phase bases, allowing MSA to condense as an ideally nonvolatile or semivolatile species or too volatile to condense. We also present two limiting cases for MSA's volatility, assuming that MSA is always ideally nonvolatile (irreversible condensation) or that MSA is always ideally semivolatile (quasi-equilibrium condensation but still irreversible condensation). We further present simulations in which MSA participates in binary and ternary nucleation with the same efficacy as sulfuric acid whenever MSA is treated as ideally nonvolatile. When using the volatility parameterization described above (both with and without nucleation), including MSA in the model changes the global annual averages at 900 hPa of submicron aerosol mass by 1.2 %, N3 (number concentration of particles greater than 3 nm in diameter) by -3.9 % (non-nucleating) or 112.5 % (nucleating), N80 by 0.8 % (non-nucleating) or 2.1 % (nucleating), the cloud-albedo aerosol indirect effect (AIE) by -8.6 mW m -2 (non-nucleating) or -26 mW m -2 (nucleating), and the direct radiative effect (DRE) by -15 mW m -2 (non-nucleating) or -14  mW m -2 (nucleating). The sulfate and sulfuric acid from DMS oxidation produces 4–6 times more submicron mass than MSA does, leading to an ~10 times stronger cooling effect in the DRE. But the changes in N80 are comparable between the contributions from MSA and from DMS-derived sulfate/sulfuric acid, leading to comparable changes in the cloud-albedo AIE. Model–measurement comparisons with the Heintzenberg et al. (2000) dataset over the Southern Ocean indicate that the default model has a missing source or sources of ultrafine particles: the cases in which MSA participates in nucleation (thus increasing ultrafine number) most closely match the Heintzenberg distributions, but we cannot conclude nucleation from MSA is the correct reason for improvement. Model–measurement comparisons with particle-phase MSA observed with a customized Aerodyne high-resolution time-of-flight aerosol mass spectrometer (AMS) from the ATom campaign show that cases with the MSA volatility parameterizations (both with and without nucleation) tend to fit the measurements the best (as this is the first use of MSA measurements from ATom, we provide a detailed description of these measurements and their calibration). However, no one model sensitivity case shows the best model–measurement agreement for both Heintzenberg and the ATom campaigns. As there are uncertainties in both MSA's behavior (nucleation and condensation) and the DMS emissions inventory, further studies on both fronts are needed to better constrain MSA's past, current, and future impacts upon the global aerosol size distribution and radiative forcing.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Integrated modeling of temperature profiles in L-mode tokamak discharges

Simulations of doublet III-D, the joint European tokamak, and the tokamak fusion test reactor L-mode tokamak plasmas are carried out using the PTRANSP predictive integrated modeling code. The simulation and experimental temperature profiles are compared. The time evolved temperature profiles are computed utilizing the Multi-Mode anomalous transport model version 7.1 (MMM7.1) which includes transport associated with drift-resistive-inertial ballooning modes (the DRIBM model [T. Rafiq et al., Phys. Plasmas 17, 082511 (2010)]). The tokamak discharges considered involved a broad range of conditions including scans over gyroradius, ITER like current ramp-up, with and without neon impurity injection, collisionality, and low and high plasma current. The comparison of simulation and experimental temperature profiles for the discharges considered is shown for the radial range from the magnetic axis to the last closed flux surface. The regions where various modes in the Multi-Mode model contribute to transport are illustrated. In the simulations carried out using the MMM7.1 model it is found that: The drift-resistive-inertial ballooning modes contribute to the anomalous transport primarily near the edge of the plasma; transport associated with the ion temperature gradient and trapped electron modes contribute in the core region but decrease in the region of the plasma boundary; and neoclassical ion thermal transport contributes mainly near the center of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigating the impact of aerosol deposition on snowmelt over the Greenland Ice Sheet using a large-ensemble kernel

Accelerating surface melt on the Greenland Ice Sheet (GrIS) has led to a doubling of Greenland's contribution to global sea level rise during recent decades. Black carbon (BC), dust, and other light-absorbing impurities (LAIs) darken the surface and enhance snowmelt by boosting the absorption of solar energy. It is therefore important for coupled aerosol–climate and ice sheet models to include snow darkening effects from LAI, and yet most do not. In this study, we conduct several thousand simulations with the Community Land Model (CLM) component of the Community Earth System Model (CESM) to characterize changes in melt runoff due to variations in the amount, timing, and nature (wet or dry) of BC deposition on the GrIS. From this large matrix of simulations, we develop a kernel relating runoff to the location, month, year (from 2006 to 2015), and magnitudes of BC concentration within precipitation and dry deposition flux. BC deposition during June–August causes the largest increase in annually integrated runoff, but winter deposition events also exert large (roughly half as great) runoff perturbations due to reexposure of impurities at the snow surface during summer melt. Current BC deposition fluxes simulated with the atmosphere component of CESM induce a climatological-mean increase in GrIS-wide runoff of ∼ 8 Gt yr −1 , or +6.8 % relative to a paired simulation without BC deposition. We also provide linear equations that relate the increase in total runoff to GrIS-wide wet and dry BC deposition fluxes. It is our hope that the runoff kernel and simple equations provided here can be used to extend the utility of state-of-the-art aerosol models.

54 ENVIRONMENTAL SCIENCES↗

Global Sensitivity Analysis and Estimation of Model Error, Toward Uncertainty Quantification in Scramjet Computations

The development of scramjet engines is an important research area for advancing hypersonic and orbital flights. Progress toward optimal engine designs requires accurate flow simulations together with uncertainty quantification. However, performing uncertainty quantification for scramjet simulations is challenging due to the large number of uncertain parameters involved and the high computational cost of flow simulations. These difficulties are addressed in this paper by developing practical uncertainty quantification algorithms and computational methods, and deploying them in the current study to large-eddy simulations of a jet in crossflow inside a simplified HIFiRE Direct Connect Rig scramjet combustor. First, global sensitivity analysis is conducted to identify influential uncertain input parameters, which can help reduce the system’s stochastic dimension. Second, because models of different fidelity are used in the overall uncertainty quantification assessment, a framework for quantifying and propagating the uncertainty due to model error is presented. In conclusion, these methods are demonstrated on a nonreacting jet-in-crossflow test problem in a simplified scramjet geometry, with parameter space up to 24 dimensions, using static and dynamic treatments of the turbulence subgrid model, and with two-dimensional and three-dimensional geometries.

33 ADVANCED PROPULSION SYSTEMS↗

Nitrogen incorporated ultrananocrystalline diamond based field emitter array for a flat-panel x-ray source

A field emission based flat-panel transmission x-ray source is being developed as an alternative for medical and industrial imaging. A field emitter array (FEA) prototype based on nitrogen incorporated ultrananocrystalline diamond film has been fabricated to be used as the electron source of this flat panel x-ray source. The FEA prototype was developed using conventional microfabrication techniques. The field emission characteristics of the FEA prototype were evaluated. Results indicated that emission current densities of the order of 6 mA/cm{sup 2} could be obtained at electric fields as low as 10 V/μm to 20 V/μm. During the prototype microfabrication process, issues such as delamination of the extraction gate and poor etching of the SiO{sub 2} insulating layer located between the emitters and the extraction layer were encountered. Consequently, alternative FEA designs were investigated. Experimental and simulation data from the first FEA prototype were compared and the results were used to evaluate the performance of alternative single and double gate designs that would yield better field emission characteristics compared to the first FEA prototype. The best simulation results are obtained for the double gate FEA design, when the diameter of the collimator gate is around 2.6 times the diameter of the extraction gate.

07 ISOTOPE AND RADIATION SOURCES↗

Local measurements of the spatial magnetic field distribution in a z-pinch plasma during and near stagnation using polarization spectroscopy

We present here the detailed measurements of radial distribution of the magnetic field in a gas-puff z-pinch plasma at the final stages of the implosion phase and at stagnation. While the measurements are chordal, the radial distribution of different charge states was utilized to measure the magnetic field locally for certain radii, so that unlike chordal measurements in general, the magnetic field radial distribution was obtained with no need for the Abel inversion of the data. The distribution was measured using the Zeeman effect via a novel spectroscopic technique, at several axial locations, and demonstrates striking features such as the peak field remaining at a radius much larger than the stagnation radius at all times. Furthermore, while the distribution observed is sometimes monotonic with respect to the radius, it is often not, a behavior that can be linked to 2D features in the plasma column resulting from the Rayleigh–Taylor instability. The current flowing through the stagnating plasma was found to be a small fraction of the total current, resulting in clearly insufficient magnetic pressure to balance the plasma pressure at stagnation. The magnetic field data, taken over several axial positions, are used to obtain the true inductance in the imploding plasma for the first time; it is found that the data cannot explain the current turnover at stagnation. A simulation with the MACH2-Tabular Collisional-Radiative Equilibrium magnetohydrodynamics code in the r–z plane shows that the peak of the magnetic field pinches to a much smaller radius than is observed in the spectroscopic data. Furthermore, the depth of the computed current turnover at stagnation is smaller than the measured one. The two observed features of a radially extended magnetic field at stagnation together with a deep current turnover are a challenge to match in simulations. Various calculations and estimates of the inductive and resistive load voltages are examined to ascertain if they are responsible for the observed current notch. Lastly, the results demonstrate that the knowledge of the true inductance in the driven load requires such magnetic-field-distribution measurements and that imaging data or electrical measurements are insufficient.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multiscale modeling and characterization for performance and safety of lithium-ion batteries

Lithium-ion batteries are highly complex electrochemical systems whose performance and safety are governed by coupled nonlinear electrochemical-electrical-thermal-mechanical processes over a range of spatiotemporal scales. Gaining an understanding of the role of these processes as well as development of predictive capabilities for design of better performing batteries requires synergy between theory, modeling, and simulation, and fundamental experimental work to support the models. This paper presents the overview of the work performed by the authors aligned with both experimental and computational efforts. In this paper, we describe a new, open source computational environment for battery simulations with an initial focus on lithium-ion systems but designed to support a variety of model types and formulations. This system has been used to create a three-dimensional cell and battery pack models that explicitly simulate all the battery components (current collectors, electrodes, and separator). The models are used to predict battery performance under normal operations and to study thermal and mechanical safety aspects under adverse conditions. This paper also provides an overview of the experimental techniques to obtain crucial validation data to benchmark the simulations at various scales for performance as well as abuse. We detail some initial validation using characterization experiments such as infrared and neutron imaging and micro-Raman mapping. In addition, we identify opportunities for future integration of theory, modeling, and experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Calculation of Eddy Currents In the CTH Vacuum Vessel and Coil Frame

Knowledge of eddy currents in the vacuum vessel walls and nearby conducting support structures can significantly contribute to the accuracy of Magnetohydrodynamics (MHD) equilibrium reconstruction in toroidal plasmas. Moreover, the magnetic fields produced by the eddy currents could generate error fields that may give rise to islands at rational surfaces or cause field lines to become chaotic. In the Compact Toroidal Hybrid (CTH) device (R0 = 0.75 m, a = 0.29 m, B ≤ 0.7 T), the primary driver of the eddy currents during the plasma discharge is the changing flux of the ohmic heating transformer. Electromagnetic simulations are used to calculate eddy current paths and profile in the vacuum vessel and in the coil frame pieces with known time dependent currents in the ohmic heating coils. MAXWELL and SPARK codes were used for the Electromagnetic modeling and simulation. MAXWELL code was used for detailed 3D finite-element analysis of the eddy currents in the structures. SPARK code was used to calculate the eddy currents in the structures as modeled with shell/surface elements, with each element representing a current loop. In both cases current filaments representing the eddy currents were prepared for input into VMEC code for MHD equilibrium reconstruction of the plasma discharge. __________________________________________________

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Band alignment of epitaxial ZnS/Zn{sub 3}P{sub 2} heterojunctions

The energy-band alignment of epitaxial zb-ZnS(001)/{alpha}-Zn{sub 3}P{sub 2}(001) heterojunctions has been determined by measurement of shifts in the phosphorus 2p and sulfur 2p core-level binding energies for various thicknesses (0.6-2.2 nm) of ZnS grown by molecular beam epitaxy on Zn{sub 3}P{sub 2}. In addition, the position of the valence-band maximum for bulk ZnS and Zn{sub 3}P{sub 2} films was estimated using density functional theory calculations of the valence-band density-of-states. The heterojunction was observed to be type I, with a valence-band offset, {Delta}E{sub V,} of -1.19 {+-} 0.07 eV, which is significantly different from the type II alignment based on electron affinities that is predicted by Anderson theory. n{sup +}-ZnS/p-Zn{sub 3}P{sub 2} heterojunctions demonstrated open-circuit voltages of >750 mV, indicating passivation of the Zn{sub 3}P{sub 2} surface due to the introduction of the ZnS overlayer. Carrier transport across the heterojunction devices was inhibited by the large conduction-band offset, which resulted in short-circuit current densities of <0.1 mA cm{sup -2} under 1 Sun simulated illumination. Hence, constraints on the current density will likely limit the direct application of the ZnS/Zn{sub 3}P{sub 2} heterojunction to photovoltaics, whereas metal-insulator-semiconductor structures that utilize an intrinsic ZnS insulating layer appear promising.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Probing flow-induced nanostructure of complex fluids in arbitrary 2D flows using a fluidic four-roll mill (FFoRM)

Abstract Engineering flow processes to direct the microscopic structure of soft materials represents a growing area of materials research. In situ small-angle neutron scattering under flow (flow-SANS) is an attractive probe of fluid microstructure under simulated processing conditions, but current capabilities require many different sample environments to fully interrogate the deformations a fluid experiences in a realistic processing flow. Inspired by recent advances in microfluidics, we present a fluidic four-roll mill (FFoRM) capable of producing tunable 2D flow fields for in situ SANS measurements, that is intended to allow characterization of complex fluid nanostructure under arbitrary complex flows within a single sample environment. Computational fluid dynamics simulations are used to design a FFoRM that produces spatially homogeneous and sufficiently strong deformation fields. Particle tracking velocimetry experiments are then used to characterize the flows produced in the FFoRM for several classes of non-Newtonian fluids. Finally, a putative FFoRM-SANS workflow is demonstrated and validated through the characterization of flow-induced orientation in a semi-dilute cellulose nanocrystal dispersion under a range of 2D deformations. These novel experiments confirm that, for steady state straining flows at moderate strain rates, the nanocrystals orient along the principal strain-rate axis, in agreement with theories for rigid, rod-like Brownian particles in a homogeneous flow.

36 MATERIALS SCIENCE↗

A Novel Model of Thyristor in Reverse Recovery Process for Current Interruption Test of HVDC Circuit Breakers

The high voltage thyristor is widely used in high voltage DC (HVDC) grids due to its large current capability. Recently, the high voltage thyristor is implemented to produce a high current in the HVDC circuit breaker (CB) testing process. Unfortunately, the reverse recovery process (RRP) of a thyristor could induce significant turn-off overvoltage and power losses, jeopardizing the safety of the devices. The existing models of a thyristor cannot describe such hard-switching transient accurately. Therefore, this paper divides this process into four stages based on the physical mechanism and proposes a novel cosine exponential (CE), simulating the voltage and current of the thyristor with analytical equations, which are derived from the equivalent circuit. Moreover, the extraction procedure of the proposed model is illustrated in this paper. Finally, we develop an HVDC CB test platform to validate the proposed model. Experimental results show that the proposed CE model can accurately predict the overvoltage with a relative error of less than 8% at various current varying rates.

Dongye, Zhonghao↗

SIDSIM (Sub-Interval Distribution-function in Simulation) [SWR-17-22]

Time-series simulation is a technique used to estimate the energy delivery of a power generation system with output that changes under changing conditions. Such variable resources include photovoltaic energy systems that experience different levels of sunlight or wind farms with changing wind patterns. As distributed generation from these variable resources continues to flourish, having improved software tools to accurately measure, model, and estimate phenomena of interest to grid integration, in addition to energy output, will be crucial for safeguarding smooth integration and operation of these energy systems. Currently, time-series simulations assume conditions and power output are constant over the duration of a time-step (steady-state assumption). NREL researchers improve the accuracy with a new simulation software called SIDSIM (Sub-Interval Distribution Simulation) that replaces this steady value (solar or wind resource) with a distribution function that has a maximum value, a minimum value, and a shape to the curve within each time-step. The distribution function used for modeling varying values in the SIDSIM software is a novel equation and presents advantageous new benefits over the conventional methods. The first advantage comes in the ability to determine maximum instantaneous power within a time-step. This enables better protection of a system from overcurrent issues and it also allows for better forecasting and reservation of transmission capacity; such capabilities are essential for energy system integration. Second, a distribution function helps to determine the amount of energy loss due to “clipping” in cases where the photovoltaic direct current (DC) system capacity exceeds the inverter capacity, the AC interconnect capacity, or any other limiting component or threshold value. It improves estimates of battery throughput by allowing both charge and discharge of a battery within the duration of a time-step. Third, it provides a more refined estimate of utility-bill impacts by considering the delivered power may be higher than the site load or lower than the regular average power estimates that are made using a quasi-steady-state assumption. SEE ALSO OSTI ID 1,361,018

Walker, Andy↗

Hanford River Protection Project Life cycle Cost Modeling Tool to Enhance Mission Planning - 13396

The Life cycle Cost Model (LCM) Tool is an overall systems model that incorporates budget, and schedule impacts for the entire life cycle of the River Protection Project (RPP) mission, and is replacing the Hanford Tank Waste Operations Simulator (HTWOS) model as the foundation of the RPP system planning process. Currently, the DOE frequently requests HTWOS simulations of alternative technical and programmatic strategies for completing the RPP mission. Analysis of technical and programmatic changes can be performed with HTWOS; however, life cycle costs and schedules were previously generated by manual transfer of time-based data from HTWOS to Primavera P6. The LCM Tool automates the preparation of life cycle costs and schedules and is needed to provide timely turnaround capability for RPP mission alternative analyses. LCM is the simulation component of the LCM Tool. The simulation component is a replacement of the HTWOS model with new capability to support life cycle cost modeling. It is currently deployed in G22, but has been designed to work in any full object-oriented language with an extensive feature set focused on networking and cross-platform compatibility. The LCM retains existing HTWOS functionality needed to support system planning and alternatives studies going forward. In addition, it incorporates new functionality, coding improvements that streamline programming and model maintenance, and capability to input/export data to/from the LCM using the LCM Database (LCMDB). The LCM Cost/Schedule (LCMCS) contains cost and schedule data and logic. The LCMCS is used to generate life cycle costs and schedules for waste retrieval and processing scenarios. It uses time-based output data from the LCM to produce the logic ties in Primavera P6 necessary for shifting activities. The LCM Tool is evolving to address the needs of decision makers who want to understand the broad spectrum of risks facing complex organizations like DOE-RPP to understand how near-term programmatic decisions affect life cycle costs and commitments. (authors)

54 ENVIRONMENTAL SCIENCES↗

Corrosion Behavior of Pipeline Steel with Different Microstructures Under AC Interference in Acid Soil Simulation Solution

Corrosion behavior of X65 pipeline steels with different microstructures under alternating current (AC) interference was investigated in acid soil simulation solution by potentiodynamic polarization curve, potentiostatic polarization curve and immersion test. The results show that superimposed AC causes a sharp increase in corrosion current density of X65 steel. With the increase in i{sub AC}, the corrosion current densities of steels with various microstructures increase, especially at high i{sub AC}. Hot-rolled steel mainly experiences uniform corrosion, with very slight pit corrosion. Serious corrosion degrees with intensive corrosion pits can be observed on the surfaces of normalized and quenched microstructure steels. The annealed steel exhibits the feature of non-uniform corrosion with some pitting. The steels with various microstructures applied with AC have different corrosion resistance. The normalized steel shows the worst corrosion resistance, then the quenched microstructure, and the hot-rolled steel displays the optimum corrosion resistance. The difference in the microstructure can result in difference in corrosion degree and occurrence position of pitting corrosion of X65 steel. The normalized microstructure composed of polygonal ferrite and a large amount of pearlite and bainite is the most susceptible to AC corrosion.

36 MATERIALS SCIENCE↗

Role of substrate aspect ratio on the robustness of capillary alignment

Capillary forces associated with liquid bridges formed across solid substrates are routinely exploited to align and assemble micro- and nanoscale devices. The magnitude of these forces plays a critical role in minimizing substrate misalignment and therefore should be controlled for robust and reliable fabrication process. We explore through simulations and experiments the role of the substrate aspect ratio (L/W) on capillary restoring forces and torques. Simulations show that increasing the aspect ratio of the substrates increases the capillary torques and forces when the substrates are misaligned through either lateral or rotational perturbations. The effect of substrate area, perimeter, and liquid volume are also systematically explored to show that the increase in restoring torque is caused by an increase in aspect ratio. A simple theoretical model based on the geometry of the system shows excellent agreement with Surface Evolver simulations. Finally, parameters from experimental flip-chip devices [Josell, D. Wallace, W.E. Warren, J.A. Wheeler, D. Powell, A.C. J. Electron. Packag. 124, 227, (2002)] are used in our simulations to show how current capillary self-alignment schemes could benefit from using rectangular substrate shapes with aspect ratio greater than one.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Carbon Dioxide Removal Measurement, Reporting, and Verification Simulation toolkit (CDR MRVSim) v0.1

CDR MRVSim is a statistical software toolkit for techno-economic analysis of measurement, reporting, and verification of multiple carbon dioxide removal technologies. The current version applies Monte Carlo simulation to existing datasets estimate the cost and uncertainty of measuring the soil organic carbon (SOC) content of agricultural land, accounting for multiple sources of uncertainty, using data from field measurements. We are planning to incorporate enhanced rock weathering and other technologies into the software. The current model leverages existing data to characterize underlying variability in soil organic carbon and related parameters, including bulk density, in an agricultural field attempting to increase its SOC. We then simulate baseline and post-intervention "measurement campaigns" in which some number of SOC and bulk density measurements are conducted using a selected technology. This enables estimates of the cost and accuracy of measuring changes in SOC in the simulated field. The primary advantages over similar software are: 1) ability to quantify tradeoffs between cost and uncertainty in MRV across a range of possible MRV approaches 2) focus on guiding development of novel sensors by determining desirable sets of characteristics

Sherwin, Evan [Lawrence Berkeley National Laborato↗

Will Advanced Public Charging Infrastructure Speed Up Electrification of Future Transportation?

This study focuses on evaluating impacts of advanced, urban public charging infrastructure on the battery electric vehicle (BEV) adoption in the U.S market. Under various infrastructure scenarios (e.g., current and developed conditions), we investigate the infrastructure impacts on the near-term BEV adoption using a consumer-choice based market simulation approach. Our results suggest that current public charging infrastructure has continuous and significant impacts on the near-term BEV adoption. Additional infrastructure investment could further stimulate the public acceptance of BEVs. We also find that the actual infrastructure impact may vary depending on the assumption of how charging deployment could affect consumers' access of chargers and consumers' daily available charging time. However, this impact variability could be reduced when infrastructure is getting mature.

33 ADVANCED PROPULSION SYSTEMS↗

Numerical study on a canonized Hamiltonian system representing reduced magnetohydrodynamics and its comparison with two-dimensional Euler system

Introducing a Clebsch-like parameterization, we have formulated a canonical Hamiltonian system on a symplectic leaf of reduced magnetohydrodynamics. An interesting structure of the equations is in that the Lorentz-force, which is a quadratic nonlinear term in the conventional formulation, appears as a linear term −ΔQ, just representing the current density (Q is a Clebsch variable, and Δ is the two-dimensional Laplacian); omitting this term reduces the system into the two-dimensional Euler vorticity equation of a neutral fluid. A heuristic estimate shows that current sheets grow exponentially (even in a fully nonlinear regime) together with the action variable P that is conjugate to Q. By numerical simulation, the predicted behavior of the canonical variables, yielding exponential growth of current sheets, has been demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗