This paper compares four different formulations of model predictive control that attempt to maximise electrical power generated by a wave energy converter (WEC). Control laws include (1) pure maximisation of mechanical power, (2) maximisation of mechanical power with a control penalty factor, (3) maximisation of electrical power using power conversion efficiency, and (4) maximisation of electrical power using the full electro-mechanical model of a system. For this study, a wave-to-wire model is developed for a floating spherical buoy connected to a permanent magnet synchronous generator. The performance of the controllers, including the mechanical and electrical power outputs, is compared in irregular wave conditions for the unconstrained and force-constrained scenarios. The results demonstrate that the controller designed to maximise mechanical power is not suitable for practical applications and may lead to negative electrical power output due to the non-ideal power take-off efficiency. Moreover, the replacement of the power take-off dynamics by the efficiency coefficient does not guarantee the maximum electrical power production.
The progression of wind turbine technology has led to wind turbines being incredibly optimized machines often approaching their theoretical maximum production capabilities. When placed together in arrays to make wind farms, however, they are subject to wake interference that greatly reduces downstream turbines' power production, increases structural loading and maintenance, reduces their lifetimes, and ultimately increases the levelized cost of energy. Development of techniques to manage wakes and operate larger and larger arrays of turbines more efficiently is now a crucial field of research. Herein, four wake management techniques in various states of development are reviewed. These include axial induction control, wake steering, the latter two combined, and active wake control. Each of these is reviewed in terms of its control strategies and use for power maximization, load reduction, and ancillary services. By evaluating existing research, several directions for future research are suggested.
Magnetized Liner Inertial Fusion (MagLIF) [Slutz et al., Phys. Plasmas 17, 056303 (2010)] experiments driven by the Z machine produce >1013 deuterium-deuterium fusion reactions [Gomez et al., Phys. Rev. Lett. 125, 155002 (2020)]. Simulations indicate high yields and gains (1000) with increased current and deuterium-tritium layers for burn propagation [Slutz et al., Phys. Plasmas 23, 022702 (2016)]. Such a coating also isolates the metal liner from the gaseous fuel, which should reduce mixing of liner material into the fuel. However, the vapor density at the triple point is only 0.3 kg/m3, which is not high enough for MagLIF operation. We present two solutions to this problem. First, a fuel wetted low-density plastic foam can be used to form a layer on the inside of the liner. The desired vapor density can be obtained by controlling the temperature. This does however introduce carbon into the layer which will enhance radiation losses. Simulations indicate that this wetted foam layer can significantly contribute to the fusion yield when the foam density is less than 35 kg/m3. Second, we show that a pure frozen fuel layer can first be formed on the inside of the liner and then low temperature gaseous fuel can be introduced just before the implosion without melting a significant amount of the ice layer. This approach is the most promising for MagLIF to produce high yield and gain.
Fast detection and isolation of faults in a DC microgrid is of particular importance. Fast tripping protection (i) increases the lifetime of power electronics (PE) switches by avoiding high fault current magnitudes and (ii) enhances the controllability of PE converters. This paper proposes a traveling wave (TW) based scheme for fast tripping protection of DC microgrids. The proposed scheme utilizes a discrete wavelet transform (DWT) to calculate the high-frequency components of DC fault currents. Multiresolution analysis (MRA) using DWT is utilized to detect TW components for different frequency ranges. The Parseval energy of the MRA coefficients are then calculated to demonstrate a quantitative relationship between the fault current signal energy and coefficients’ energy. The calculated Parseval energy values are used to train a Support Vector Machine classifier to identify the fault type and a Gaussian Process regression engine to estimate the fault location on the DC cables. The proposed approach is verified by simulating two microgrid test systems in PSCAD/EMTDC.
With a wide variety of wave energy device archetypes currently under consideration, it is a major challenge to ensure that research findings and methods are broadly applicable. In particular, the design and testing of wave energy control systems, a process which includes experimental design, empirical modeling, control design, and performance evaluation, is of interest. This goal motivated the redesign and testing of a floating dual flap wave energy converter. As summarized in this paper, the steps taken in the design, testing, and analysis of the device mirrored those previously demonstrated on a three-degree of freedom point absorber device. The method proposed does not require locking WEC degrees of freedom to develop an excitation model, and presents a more attainable system identification procedure for at-sea deployments. The results show that the methods employed work well for this dual flap device, lending additional support for the broad applicability of the design and testing methods applied here. The aim of this paper is to demonstrate that these models are particularly useful for deducing areas of device design or controller implementation that can be reasonably improved to increase device power capture.
We explore the character angle dependence of dislocation-solute interactions in a face-centered cubic random Fe0.70Ni0.11Cr0.19 alloy through molecular dynamics (MD) simulations of dislocation mobility. Using the MD mobility data, we determine the phonon and thermally activated solute drag parameters which govern mobility for each dislocation character angle. The resulting parameter set indicates that, surprisingly, the solute energy barrier does not depend on character angle. Instead, only the zero-temperature flow stress—which is dictated by the activation area for thermal activation—is dependent on character angle. By analyzing the line roughness from MD simulations and the geometry of a bowing dislocation line undergoing thermal activation, we conclude that the character angle dependence of the activation area in this alloy is governed by the dislocation line tension, rather than the dislocation-solute interaction itself. Our findings motivate further investigation into the line geometry of dislocations in solid solutions.
Kiani, Mehrdad T.; Gan, Lucia T.; Traylor, Rachel; Yang, Rui; Barr, Christopher M.; Hattar, Khalid; Fan, Jonathan A.; Wendy Gu, X.
Grain boundaries have complex structural features that influence strength, ductility and fracture in metals and alloys. Grain boundary misorientation angle has been identified as a key parameter that controls their mechanical behavior, but the effect of misorientation angle has been challenging to isolate in polycrystalline materials. Here, we describe the use of bicrystal Au thin films made using a rapid melt growth process to study deformation at a single grain boundary. Tensile testing is performed on bicrystals with different misorientation angles using in situ TEM, as well as on a single crystalline sample. Plastic deformation is initiated through dislocation nucleation from free surfaces. Grain boundary sliding is not observed, and failure occurs away from the grain boundary through plastic collapse in all cases. The failure behavior in these nanoscale bicrystals does not appear to depend on the misorientation angle or grain boundary energy but instead has a more complex dependence on sample surface structure and dislocation activity.
Dual-fuel (DF) engines, in which premixed natural gas and air in an open-type combustion chamber is ignited by diesel-fuel pilot sprays, have been more popular for marine use than pre-chamber spark ignition (PCSI) engines because of their superior durability. However, control of ignition and combustion in DF engines is more difficult than in PCSI engines. In this context, this study focuses on the ignition stability of n-heptane pilot-fuel jets injected into a compressed premixed charge of natural gas and air at low-load conditions. To aid understanding of the experimental data, chemical-kinetics simulations were carried out in a simplified engine-environment that provided insight into the chemical effects of methane (CH4) on pilot-fuel ignition. The simulations reveal that CH4 has an effect on both stages of n-heptane autoignition: the small, first-stage, cool-flame-type, low-temperature ignition (LTI) and the larger, second-stage, high-temperature ignition (HTI). As the ratio of pilot-fuel to CH4 entrained into the spray decreases, the initial oxidization of CH4 consumes the OH radicals produced by pilot-fuel decomposition during LTI, thereby inhibiting its progression to HTI. Using imaging diagnostics, the spatial and temporal progression of LTI and HTI in DF combustion are measured in a heavy-duty optical engine, and the imaging data are analyzed to understand the cause of severe fluctuations in ignition timing and combustion completeness at low-load conditions. Images of cool-flame and hydroxyl radical (OH*) chemiluminescence serve as indicators of LTI and HTI, respectively. The cycle-to-cycle and spatial variation in ignition extracted from the imaging data are used as key metrics of comparison. The imaging data indicate that the local concentration of the pilot-fuel and the richness of the surrounding natural-gas air mixture are important for LTI and HTI, but in different ways. In particular, higher injection pressures and shorter injection durations increase the mixing rate, leading to lower concentrations of pilot-fuel more quickly, which can inhibit HTI even as LTI remains relatively robust. Decreasing the injection pressure from 80 MPa to 40 MPa and increasing the injection duration from 500 µs to 760 µs maintained constant pilot-fuel mass, while promoting robust transition from LTI to HTI by effectively slowing the mixing rate. This allows enough residence time for the OH radicals, produced by the two-stage ignition chemistry of the pilot-fuel, to accelerate the transition from LTI to HTI before being consumed by CH4 oxidation. Thus from a practical perspective, for a premixed natural gas fuel–air equivalence-ratio, it is possible to improve the “stability” of the combustion process by solely manipulating the pilot-fuel injection parameters while maintaining constant mass of injected pilot-fuel. This allows for tailoring mixing trajectories to offset changes in fuel ignition chemistry, so as to promote a robust transition from LTI to HTI by changing the balance between the local concentration of the pilot-fuel and richness of the premixed natural gas and air. This could prove to be a valuable tool for combustion design to improve fuel efficiency or reduce noise or perhaps even reduce heat-transfer losses by locating early combustion away from in-cylinder walls.
Material produced by current metal additive manufacturing processes is susceptible to variable performance due to imprecise control of internal porosity, surface roughness, and conformity to designed geometry. Using a double U-notched specimen, we investigate the interplay of nominal geometry and porosity in determining ductile failure characteristics during monotonic tensile loading. We simulate the effects of distributed porosity on plasticity and damage using a statistical model based on populations of pores visible in computed tomography scans and additional sub-threshold voids required to match experimental observations of deformation and failure. We interpret the simulation results from a physical viewpoint and provide a statistical model of the probability of failure near stress concentrations. We provide guidance for designs where material defects could cause unexpected failures depending on the relative importance of these defects with respect to features of the nominal geometry.
Electromagnetic (EM) methods are among the original techniques for subsurface characterization in exploration geophysics because of their particular sensitivity to the earth electrical conductivity, a physical property of rocks distinct yet complementary to density, magnetization, and strength. However, this unique ability also makes them sensitive to metallic artifacts - infrastructure such as pipes, cables, and other forms of cultural clutter - the EM footprint of which often far exceeds their diminutive stature when compared to that of bulk rock itself. In the hunt for buried treasure or unexploded ordnance, this is an advantage; in the long-term monitoring of mature oil fields after decades of production, it is quite troublesome indeed. Here we consider the latter through the lens of an evolving energy industry landscape in which the traditional methods of EM characterization for the exploration geophysicist are applied toward emergent problems in well-casing integrity, carbon capture and storage, and overall situational awareness in the oil field. We introduce case studies from these exemplars, showing how signals from metallic artifacts can dominate those from the target itself and impose significant burdens on the requisite simulation complexity. We also show how recent advances in numerical methods mitigate the computational explosivity of infrastructure modeling, providing feasible and real-time analysis tools for the desktop geophysicist. Lastly, we demonstrate through comparison of field data and simulation results that incorporation of infrastructure into the analysis of such geophysical data is, in a growing number of cases, a requisite but now manageable step.
Poondla, Yasvanth; Goldstein, David; Varghese, Philip; Clarke, Peter; Moore, Christopher H.
The goal of this work is to build up the capability of quasi-particle simulation (QuiPS), a novel flow solver, such that it can adequately model the rarefied portion of an atmospheric reentry trajectory. Direct simulation Monte Carlo (DSMC) is the conventional solver for such conditions, but struggles to resolve transient flows, trace species, and high-level internal energy states due to stochastic noise. Quasi-particle simulation (QuiPS) is a novel Boltzmann solver that describes a system with a discretized, truncated velocity distribution function. The resulting fixed-velocity, variable weight quasi-particles enable smooth variation of macroscopic properties. The distribution function description enables the use of a variance-reduced collision model, greatly minimizing expense near equilibrium. This work presents the addition of a neutral air chemistry model to QuiPS and some demonstrative 0D simulations. The explicit representation of internal distributions in QuiPS reveals some of the flaws in existing physics models. Variance reduction, a key feature of QuiPS, can greatly reduce expense of multi-dimensional calculations, but is only cheaper when the gas composition is near chemical equilibrium.
Cookoff experiments of powdered and pressed TATB-based plastic bonded explosives (PBXs) have been modeled using a pressure-dependent universal cookoff model (UCM) in combination with a micromechanics pressurization (MMP) model described in a companion paper. The MMP model is based on the accumulation of decomposition gases at nucleation sites that load the surrounding TATB crystals and binder. This is the first cookoff model to use an analytical mechanics solution for compressibility and thermal expansion to describe internal pressurization caused by both temperature and decomposition occurring within closed-pore explosives. This approach produces more accurate predictions of ignition time and pressurization within high-density explosives than simple equation-of-state models. The current paper gives details of the reaction chemistry, model parameters, predicted uncertainty, and validation using experiments from multiple laboratories with errors less than 6 %. The UCM/MMP model framework gives more accurate thermal ignition predictions for high density explosives that are initially impermeable to decomposition gases.
Neuromorphic computing, which aims to replicate the computational structure and architecture of the brain in synthetic hardware, has typically focused on artificial intelligence applications. What is less explored is whether such brain-inspired hardware can provide value beyond cognitive tasks. Here we show that the high degree of parallelism and configurability of spiking neuromorphic architectures makes them well suited to implement random walks via discrete-time Markov chains. These random walks are useful in Monte Carlo methods, which represent a fundamental computational tool for solving a wide range of numerical computing tasks. Using IBM’s TrueNorth and Intel’s Loihi neuromorphic computing platforms, we show that our neuromorphic computing algorithm for generating random walk approximations of diffusion offers advantages in energy-efficient computation compared with conventional approaches. We also show that our neuromorphic computing algorithm can be extended to more sophisticated jump-diffusion processes that are useful in a range of applications, including financial economics, particle physics and machine learning.
Triangle counting is a fundamental building block in graph algorithms. In this article, we propose a block-based triangle counting algorithm to reduce data movement during both sequential and parallel execution. Our block-based formulation makes the algorithm naturally suitable for heterogeneous architectures. The problem of partitioning the adjacency matrix of a graph is well-studied. Our task decomposition goes one step further: it partitions the set of triangles in the graph. By streaming these small tasks to compute resources, we can solve problems that do not fit on a device. We demonstrate the effectiveness of our approach by providing an implementation on a compute node with multiple sockets, cores and GPUs. The current state-of-the-art in triangle enumeration processes the Friendster graph in 2.1 seconds, not including data copy time between CPU and GPU. Using that metric, our approach is 20 percent faster. When copy times are included, our algorithm takes 3.2 seconds. This is 5.6 times faster than the fastest published CPU-only time.
Hsu, Nai W.; Hou, Wei C.; Chen, Yen Y.; Wu, Yu J.; Kao, Hsiang S.; Harris, Charles T.; Lu, T.M.; Li, Jiun Y.
Capacitance-voltage ( {C} - {V} ) characteristics and carrier transport properties of 2-D electron gases (2DEGs) in an undoped Si/SiGe heterostructure at {T}= {4} - {35} K are presented. Two capacitance plateaus due to density saturation of the 2DEG in the buried Si quantum well (QW) are observed and explained by a model of surface tunneling. The peak mobility at 4 K is 4.1 \times 10^{{5}} cm2/ \text{V}\cdot \text{s} and enhanced by a factor of 1.97 at an even lower carrier density compared to the saturated carrier density, which is attributed to the effect of remote carrier screening. At {T}\,\,=35 K, the mobility enhancement with a factor of 1.35 is still observed, which suggests the surface tunneling is still dominant.
Agencies that monitor for underground nuclear tests are interested in techniques that automatically characterize mining blasts to reduce the human analyst effort required to produce high-quality event bulletins. Waveform correlation is effective in finding similar waveforms from repeating seismic events, including mining blasts. We report the results of an experiment to detect and identify mining blasts for two regions, Wyoming (U.S.A.) and Scandinavia, using waveform templates recorded by multiple International Monitoring System stations of the Preparatory Commission for the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO PrepCom) for up to 10 yr prior to the time of interest. We discuss approaches for template selection, threshold setting, and event detection that are specialized for characterizing mining blasts using a sparse, global network. We apply the approaches to one week of data for each of the two regions to evaluate the potential for establishing a set of standards for waveform correlation processing of mining blasts that can be generally applied to operational monitoring systems with a sparse network. We compare candidate events detected with our processing methods to the Reviewed Event Bulletin of the International Data Centre to assess potential reduction in analyst workload.
The self-magnetic pinch (SMP) diode is a type of radiographic diode used to generate an intense electron beam for radiographic applications. At Sandia National Laboratories, SMP was the diode load for the six-cavity radiographic integrated test stand inductive voltage adder (IVA) driver operated in a magnetically insulated transmission line (MITL). The MITL contributes a flow current in addition to the current generated within the diode itself. Extensive experiments with a MITL of 40 ω load impedance [T. J. Renk et al., Phys. Plasmas 29, 023105 (2022)] indicate that the additional flow current leads to results similar to what might be expected from a conventional high-voltage interface driver, where flow current is not present. However, when the MITL flow impedance was increased to 80 ω, qualitatively different diode behavior was observed. This includes large retrapping waves suggestive of an initial coupling to low impedance as well as diode current decreasing with time even as the total current does not. A key observation is that the driver generates total current (flow + diode) consistent with the flow impedance of the MITL used. The case is made in this paper that the 80 ω MITL experiments detailed here can only be understood when the IVA-MITL-SMP diode is considered as a total system. The constraint of fixed total current plus the relatively high flow impedance limits the ability of the diode (whether SMP or other type) to act as an independent load. An unexpected new result is that in tracking the behavior of the electron strike angle on the converter as a function of time, we observed that the conventional cIVx "Radiographic"radiation scaling (where x ∼2.2) begins to break down for voltages above 8 MV, and cubic scaling is required to recover accurate angle tracking.
The Multi-Fidelity Toolkit (MFTK) is a simulation tool being developed at Sandia National Laboratories for aerodynamic predictions of compressible flows over a range of physics fidelities and computational speeds. These models include the Reynolds-Averaged-Navier-Stokes (RANS) equations, the Euler equations, and modified Newtonian aerodynamics (MNA) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal code- and solution-verification activities as well as model validation with uncertainty considerations. Code verification is performed on the MNA model by comparing with an analytical solution for flat-plate and inclined-plate geometries. Solution-verification activities include grid-refinement studies of HIFiRE-1 wind tunnel measurements, which are used for validation, for all model fidelities. A thorough treatment of the validation comparison with prediction error and validation uncertainty is also presented.
