Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ENFORCEMENT”

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 235 records · Page 13

Security Policies for Mitigating the Risk of Load Altering Attacks on Smart Grid Systems

While demand response programs implement energy efficiency and power quality objectives, they bring potential security threats to the Smart Grid. The ability to influence load in a system enables attackers to cause system failures and impacts the quality and integrity of power delivered to customers. This paper presents a security mechanism to monitor and control load according to a set of security policies during normal system operation. The mechanism monitors, detects, and responds to load altering attacks. We examined the security requirements of Smart Grid stakeholders and constructed a set of load control policies enforced by the mechanism. We implemented a proof of concept prototype and tested it using the simulation environment. By enforcing the proposed policies in this prototype, the system is maintained in a safe state in the presence of load drop attacks.

Ryutov, Tatyana↗

Distribution-Agnostic Stochastic Optimal Power Flow for Distribution Grids

This paper outlines a data-driven, distributionally robust approach to solve chance-constrained AC optimal power flow problems in distribution networks. Uncertain forecasts for loads and power generated by photovoltaic (PV) systems are considered, with the goal of minimizing PV curtailment while meeting power flow and voltage regulation constraints. A data- driven approach is utilized to develop a distributionally robust conservative convex approximation of the chance-constraints; particularly, the mean and covariance matrix of the forecast errors are updated online, and leveraged to enforce voltage regulation with predetermined probability via Chebyshev-based bounds. By combining an accurate linear approximation of the AC power flow equations with the distributionally robust chance constraint reformulation, the resulting optimization problem becomes convex and computationally tractable.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A Blockchain Based Co-Simulation Framework for Integrating DERs to Wholesale Electricity Markets

As the number of distributed energy resources (DERs) continue to increase across energy-delivery systems, there remains a need for integrating their capabilities into traditional grid operations. In this paper, a blockchain-based solution is proposed to facilitate FERC's Order No. 2222 implementations. The presented use-case enables small-scale DERs to participate in wholesale market operations through DER aggregators, while also enabling local distribution system operators to enforce distribution system constraints in a secure and traceable manner. The presented use case is built around the Unified Testing Platform (UTP) being developed as a part of the Blockchain for Optimized Security and Energy Management (BLOSEM) project. This is a multi-lab effort intended to simplify the deployment of blockchain-powered grid solutions by enabling the integration of simulation tools, and blockchain technologies through the use of system-agnostic interfaces that provide a modular, interoperable, and reusable connectivity layer.

blockchain interoperability↗

A Hierarchical Local Electricity Market for a DER-rich Grid Edge

We report with increasing penetration of distributed energy resources (DER) in the distribution system, it is critical to design market structures that enable smooth integration of DERs. A hierarchical local electricity market (LEM) structure is proposed in this paper with a secondary market (SM) at the lower level representing secondary feeders and a primary market (PM) at the upper level, representing primary feeders, in order to effectively use DERs to increase grid efficiency and resilience. The lower level SM enforces budget, power balance, and flexibility constraints and accounts for costs related to consumers, such as their disutility, flexibility limits, and commitment reliability, while the upper level PM enforces grid physics constraints such as power balance and capacity limits, and also minimizes line losses. The hierarchical LEM is extensively evaluated using a modified IEEE-123 bus with high DER penetration, with each primary feeder consisting of at least three secondary feeders. Data from a GridLAB-D model is used to emulate realistic power injections and load profiles over the course of 24 hours. The performance of the LEM is illustrated by delineating the family of power-injection profiles across the primary and secondary feeders as well as corresponding local electricity tariffs that vary across the distribution grid. Through numerical simulations, the hierarchical LEM is shown to improve the efficiency of the market in terms of lowering overall costs, including both the distribution-level locational marginal prices (d-LMP) as well as retail tariffs paid by customers. Together, it represents an overall framework for a Distribution System Operator (DSO) who can provide the oversight for the entire LEM.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Radiation dose reduction in medical x-ray CT via Fourier-based iterative reconstruction

Purpose: A Fourier-based iterative reconstruction technique, termed Equally Sloped Tomography (EST), is developed in conjunction with advanced mathematical regularization to investigate radiation dose reduction in x-ray CT. The method is experimentally implemented on fan-beam CT and evaluated as a function of imaging dose on a series of image quality phantoms and anonymous pediatric patient data sets. Numerical simulation experiments are also performed to explore the extension of EST to helical cone-beam geometry. Methods: EST is a Fourier based iterative algorithm, which iterates back and forth between real and Fourier space utilizing the algebraically exact pseudopolar fast Fourier transform (PPFFT). In each iteration, physical constraints and mathematical regularization are applied in real space, while the measured data are enforced in Fourier space. The algorithm is automatically terminated when a proposed termination criterion is met. Experimentally, fan-beam projections were acquired by the Siemens z-flying focal spot technology, and subsequently interleaved and rebinned to a pseudopolar grid. Image quality phantoms were scanned at systematically varied mAs settings, reconstructed by EST and conventional reconstruction methods such as filtered back projection (FBP), and quantified using metrics including resolution, signal-to-noise ratios (SNRs), and contrast-to-noise ratios (CNRs). Pediatric data sets were reconstructed at their original acquisition settings and additionally simulated to lower dose settings for comparison and evaluation of the potential for radiation dose reduction. Numerical experiments were conducted to quantify EST and other iterative methods in terms of image quality and computation time. The extension of EST to helical cone-beam CT was implemented by using the advanced single-slice rebinning (ASSR) method. Results: Based on the phantom and pediatric patient fan-beam CT data, it is demonstrated that EST reconstructions with the lowest scanner flux setting of 39 mAs produce comparable image quality, resolution, and contrast relative to FBP with the 140 mAs flux setting. Compared to the algebraic reconstruction technique and the expectation maximization statistical reconstruction algorithm, a significant reduction in computation time is achieved with EST. Finally, numerical experiments on helical cone-beam CT data suggest that the combination of EST and ASSR produces reconstructions with higher image quality and lower noise than the Feldkamp Davis and Kress (FDK) method and the conventional ASSR approach. Conclusions: A Fourier-based iterative method has been applied to the reconstruction of fan-bean CT data with reduced x-ray fluence. This method incorporates advantageous features in both real and Fourier space iterative schemes: using a fast and algebraically exact method to calculate forward projection, enforcing the measured data in Fourier space, and applying physical constraints and flexible regularization in real space. Our results suggest that EST can be utilized for radiation dose reduction in x-ray CT via the readily implementable technique of lowering mAs settings. Numerical experiments further indicate that EST requires less computation time than several other iterative algorithms and can, in principle, be extended to helical cone-beam geometry in combination with the ASSR method.

60 APPLIED LIFE SCIENCES↗

Cardiac C-arm computed tomography using a 3D + time ROI reconstruction method with spatial and temporal regularization

Purpose: Reconstruction of the beating heart in 3D + time in the catheter laboratory using only the available C-arm system would improve diagnosis, guidance, device sizing, and outcome control for intracardiac interventions, e.g., electrophysiology, valvular disease treatment, structural or congenital heart disease. To obtain such a reconstruction, the patient's electrocardiogram (ECG) must be recorded during the acquisition and used in the reconstruction. In this paper, the authors present a 4D reconstruction method aiming to reconstruct the heart from a single sweep 10 s acquisition. Methods: The authors introduce the 4D RecOnstructiOn using Spatial and TEmporal Regularization (short 4D ROOSTER) method, which reconstructs all cardiac phases at once, as a 3D + time volume. The algorithm alternates between a reconstruction step based on conjugate gradient and four regularization steps: enforcing positivity, averaging along time outside a motion mask that contains the heart and vessels, 3D spatial total variation minimization, and 1D temporal total variation minimization. Results: 4D ROOSTER recovers the different temporal representations of a moving Shepp and Logan phantom, and outperforms both ECG-gated simultaneous algebraic reconstruction technique and prior image constrained compressed sensing on a clinical case. It generates 3D + time reconstructions with sharp edges which can be used, for example, to estimate the patient's left ventricular ejection fraction. Conclusions: 4D ROOSTER can be applied for human cardiac C-arm CT, and potentially in other dynamic tomography areas. It can easily be adapted to other problems as regularization is decoupled from projection and back projection.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Prior data assisted compressed sensing: A novel MR imaging strategy for real time tracking of lung tumors

Purpose: Hybrid radiotherapy-MRI devices promise real time tracking of moving tumors to focus the radiation portals to the tumor during irradiation. This approach will benefit from the increased temporal resolution of MRI's data acquisition and reconstruction. In this work, the authors propose a novel spatial-temporal compressed sensing (CS) imaging strategy for the real time MRI–-prior data assisted compressed sensing (PDACS), which aims to improve the image quality of the conventional CS without significantly increasing reconstruction times. Methods: Conventional 2D CS requires a random sampling of partial k-space data, as well as an iterative reconstruction that simultaneously enforces the image's sparsity in a transform domain as well as maintains the fidelity to the acquired k-space. PDACS method requires the additional acquisition of the prior data, and for reconstruction, it additionally enforces fidelity to the prior k-space domain similar to viewsharing. In this work, the authors evaluated the proposed PDACS method by comparing its results to those obtained from the 2D CS and viewsharing methods when performed individually. All three methods are used to reconstruct images from lung cancer patients whose tumors move and who are likely to benefit from lung tumor tracking. The patients are scanned, using a 3T MRI, under free breathing using the fully sampled k-space with 2D dynamic bSSFP sequence in a sagittal plane containing lung tumor. These images form a reference set for the evaluation of the partial k-space methods. To create partial k-space, the fully sampled k-space is retrospectively undersampled to obtain a range of acquisition acceleration factors, and reconstructed with 2D-CS, PDACS, and viewshare methods. For evaluation, metrics assessing global image artifacts as well as tumor contour shape fidelity are determined from the reconstructed images. These analyses are performed both for the original 3T images and those at a simulated 0.5T equivalent noise level. Results: In the 3.0T images, the PDACS strategy is shown to give superior results compared to viewshare and conventional 2D CS using all metrics. The 2D-CS tends to perform better than viewshare at the low acceleration factors, while the opposite is true at the high acceleration factors. At simulated 0.5T images, PDACS method performs only marginally better than the viewsharing method, both of which are superior compared to 2D CS. The PDACS image reconstruction time (0.3 s/image) is similar to that of the conventional 2D CS. Conclusions: The PDACS method can potentially improve the real time tracking of moving tumors by significantly increasing MRI's data acquisition speeds. In 3T images, the PDACS method does provide a benefit over the other two methods in terms of both the overall image quality and the ability to accurately and automatically contour the tumor shape. MRI's data acquisition may be accelerated using the simpler viewsharing strategy at the lower, 0.5T magnetic field, as the marginal benefit of the PDACS method may not justify its additional reconstruction times.

60 APPLIED LIFE SCIENCES↗

SU-E-I-01: Iterative CBCT Reconstruction with a Feature-Preserving Penalty

Purpose: Low-dose CBCT is desired in various clinical applications. Iterative image reconstruction algorithms have shown advantages in suppressing noise in low-dose CBCT. However, due to the smoothness constraint enforced during the reconstruction process, edges may be blurred and image features may lose in the reconstructed image. In this work, we proposed a new penalty design to preserve image features in the image reconstructed by iterative algorithms. Methods: Low-dose CBCT is reconstructed by minimizing the penalized weighted least-squares (PWLS) objective function. Binary Robust Independent Elementary Features (BRIEF) of the image were integrated into the penalty of PWLS. BRIEF is a general purpose point descriptor that can be used to identify important features of an image. In this work, BRIEF distance of two neighboring pixels was used to weigh the smoothing parameter in PWLS. For pixels of large BRIEF distance, weaker smooth constraint will be enforced. Image features will be better preserved through such a design. The performance of the PWLS algorithm with BRIEF penalty was evaluated by a CatPhan 600 phantom. Results: The image quality reconstructed by the proposed PWLS-BRIEF algorithm is superior to that by the conventional PWLS method and the standard FDK method. At matched noise level, edges in PWLS-BRIEF reconstructed image are better preserved. Conclusion: This study demonstrated that the proposed PWLS-BRIEF algorithm has great potential on preserving image features in low-dose CBCT.

60 APPLIED LIFE SCIENCES↗

Conservative Multimoment Transport along Characteristics for Discontinuous Galerkin Methods

A set of algorithms based on characteristic discontinuous Galerkin methods is presented for tracer transport on the sphere. The algorithms are designed to reduce message passing interface communication volume per unit of simulated time relative to current methods generally, and to the spectral element scheme employed by the U.S. Department of Energy's Exascale Earth System Model (E3SM) specifically. Two methods are developed to enforce discrete mass conservation when the transport schemes are coupled to a separate dynamics solver; constrained transport and Jacobian-combined transport. A communication-efficient method is introduced to enforce tracer consistency between the transport scheme and dynamics solver; this method also provides the transport scheme's shape preservation capability. A subset of the algorithms derived here is implemented in E3SM and shown to improve transport performance by a factor of 2.2 for the model's standard configuration with 40 tracers at the strong scaling limit of one element per core. Read More: https://epubs.siam.org/doi/10.1137/18M1165943

97 MATHEMATICS AND COMPUTING↗

Analysis of Electric Vehicle Charging Behavior Patterns with Function Principal Component Analysis Approach

This manuscript focused on analyzing electric vehicles’ (EV) charging behavior patterns with a functional data analysis (FDA) approach, with the goal of providing theoretical support to the EV infrastructure planning and regulation, as well as the power grid load management. 5-year real-world charging log data from a total of 455 charging stations in Kansas City, Missouri, was used. The focuses were placed on analyzing the daily usage occupancy variability, daily energy consumption variability, and station-level usage variability. Compared with the traditional discrete-based analysis models, the proposed FDA modeling approach had unique advantages in preserving the smooth function behavior of the data, bringing more flexibility in the modeling process with little required assumptions or background knowledge on independent variables, as well as the capability of handling time series data with different lengths or sizes. In addition to the patterns revealed in the EV charging station’s occupancy and energy consumption, the differences between EV driver’s charging time and parking time were analyzed and called for the needs for parking regulation and enforcement. The different usage patterns observed at charging stations located on different land-use types were also analyzed.

Engineering↗

Charliecloud

Charliecloud is a set of scripts to let users run a virtual cluster of virtual machines (VMs) on a desktop or supercomputer. Key functions include: 1. Creating (typically by installing an operating system from vendor media) and updating VM images; 2. Running a single VM; 3. Running multiple VMs in a virtual cluster. The virtual machines can talk to one another over the network and (in some cases) the outside world. This is accomplished by calling external programs such as QEMU and the Virtual Distributed Ethernet (VDE) suite. The goal is to let users have a virtual cluster containing nodes where they have privileged access, while isolating that privilege within the virtual cluster so it cannot affect the physical compute resources. Host configuration enforces security; this is not included in Charliecloud, though security guidelines are included in its documentation and Charliecloud is designed to facilitate such configuration. Charliecloud manages passing information from host computers into and out of the virtual machines, such as parameters of the virtual cluster, input data specified by the user, output data from virtual compute jobs, VM console display, and network connections (e.g., SSH or X11). Parameters for the virtual cluster (number of VMs, RAM and disk per VM, etc.) are specified by the user or gathered from the environment (e.g., SLURM environment variables). Example job scripts are included. These include computation examples (such as a "hello world" MPI job) as well as performance tests. They also include a security test script to verify that the virtual cluster is appropriately sandboxed. Tests include: 1. Pinging hosts inside and outside the virtual cluster to explore connectivity; 2. Port scans (again inside and outside) to see what services are available; 3. Sniffing tests to see what traffic is visible to running VMs; 4. IP address spoofing to test network functionality in this case; 5. File access tests to make sure host access permissions are enforced. This test script is not a comprehensive scanner and does not test for specific vulnerabilities. Importantly, no information about physical hosts or network topology is included in this script (or any of Charliecloud); while part of a sensible test, such information is specified by the user when the test is run. That is, one cannot learn anything about the LANL network or computing infrastructure by examining Charliecloud code.

Priedhorsky, Reid↗

A directional microcrack damage theory for brittle solids based on continuous hyperplasticity

We present a novel microcrack-damage theory for brittle solids under compression. Instead of using internal variables like zeroth, second or fourth rank damage tensors, the state of material damage is represented by an internal function that encapsulates the information of direction, density and size of microcracks. Just like other internal variables, the evolution of this state function must obey the second law of thermodynamics for arbitrary loading paths. This is done by casting the model in the framework of continuous hyperplasticity and enforcing a non-negative dissipation rate functional. The proposed framework offers predictions on the continuous evolution of microcrack density and the induced material anisotropy along with the macroscopic stress-strain curves. The use of continuous damage function grants the model significantly enhanced resolution in characterizing the direction-dependent response of cracked solids compared to classical models that are based on damage tensors. Two scenarios are considered in developing the theory, one assumes frictionless cracks and the other incorporates friction between crack surfaces. The results highlight that inelasticity, pressure dependence, and loading-unloading hysteresis exhibited by brittle solids are natural consequences of frictional microcracks. The proposed theory offers a generic and versatile framework to upscale micromechanical processes operating at individual crack scale to explain the macroscopic behavior of cracked solids.

Materials Science↗

Analysis of the Chinese Market for Building Energy Efficiency

China will account for about half of the new construction globally in the coming decade. Its floorspace doubled from 1996 to 2011, and Chinese rural buildings alone have as much floorspace as all of U.S. residential buildings. Building energy consumption has also grown, increasing by over 40% since 1990. To curb building energy demand, the Chinese government has launched a series of policies and programs. Combined, this growth in buildings and renovations, along with the policies to promote green buildings, are creating a large market for energy efficiency products and services. This report assesses the impact of China’s policies on building energy efficiency and on the market for energy efficiency in the future. The first chapter of this report introduces the trends in China, drawing on both historical analysis, and detailed modeling of the drivers behind changes in floorspace and building energy demand such as economic and population growth, urbanization, policy. The analysis describes the trends by region, building type and energy service. The second chapter discusses China’s policies to promote green buildings. China began developing building energy codes in the 1980s. Over time, the central government has increased the stringency of the code requirements and the extent of enforcement. The codes are mandatory in all new buildings and major renovations in China’s cities, and they have been a driving force behind the expansion of China’s markets for insulation, efficient windows, and other green building materials. China also has several other important policies to encourage efficient buildings, including the Three-Star Rating System (somewhat akin to LEED), financial incentives tied to efficiency, appliance standards, a phasing out of incandescent bulbs and promotion of efficient lighting, and several policies to encourage retrofits in existing buildings. In the third chapter, we take “deep dives” into the trends affecting key building components. This chapter examines insulation in walls and roofs; efficient windows and doors; heating, air conditioning and controls; and lighting. These markets have seen significant growth because of the strength of the construction sector but also the specific policies that require and promote efficient building components. At the same time, as requirements have become more stringent, there has been fierce competition, and quality has at time suffered, which in turn has created additional challenges. Next we examine existing buildings in chapter four. China has many Soviet-style, inefficient buildings built before stringent requirements for efficiency were more widely enforced. As a result, there are several specific market opportunities related to retrofits. These fall into two or three categories. First, China now has a code for retrofitting residential buildings in the north. Local governments have targets of the number of buildings they must retrofit each year, and they help finance the changes. The requirements focus on insulation, windows, and heat distribution. Second, the Chinese government recently decided to increase the scale of its retrofits of government and state-owned buildings. It hopes to achieve large scale changes through energy service contracts, which creates an opportunity for energy service companies. Third, there is also a small but growing trend to apply energy service contracts to large commercial and residential buildings. This report assesses the impacts of China’s policies on building energy efficiency. By examining the existing literature and interviewing stakeholders from the public, academic, and private sectors, the report seeks to offer an in-depth insights of the opportunities and barriers for major market segments related to building energy efficiency. The report also discusses trends in building energy use, policies promoting building energy efficiency, and energy performance contracting for public building retrofits.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Performance Analysis: Issues Tracking System Data through June 2013

This report presents the results from an analysis of issues and assessments in LLNL’s Issues Tracking System (ITS). The analysis is conducted to identify issues that may require additional management attention and noncompliances that may not have been previously identified that meet the threshold for reporting to the DOE Noncompliance Tracking System (NTS) or to the DOE Safeguards and Security Information Management System (SSIMS). This report includes ITS data within the DOE Office of Enforcement regulated subjects through June 2013. The analysis in this report primarily focuses on deficiencies identified in the last twelve months. The analysis of issues concluded that data for 17 of the 25 Office of Enforcement regulated safety/security subjects were within expected variation. The data for eight met a common test and were discussed further. Three of the eight safety subjects met an action limit and were analyzed to resolution. Analysis of data from one of the safety subjects where a common test was met identified a potential safety related significant or programmatic (systemic) noncompliance with the LLNL work planning and control process. Additional analysis is needed to make a final determination. Note: Since the analysis of WPC deficiencies (through June 2013) was completed, a noncompliance evaluation was also completed in November 2013 and included more recent WPC deficiencies. The analysis of ITS data through June 2013, and the results of the noncompliance evaluation completed in November 2013 led to LLNL filing a programmatic noncompliance report to the NTS in January 2014 titled, Programmatic noncompliance with topical areas within LLNL’s Work Planning and Control process.

99 GENERAL AND MISCELLANEOUS↗

AVTA Federal Fleet PEV Readiness Data Logging and Characterization Study for the National Park Service: Golden Gate National Recreation Area

Battelle Energy Alliance, LLC, managing and operating contractor for the U.S. Department of Energy's Idaho National Laboratory, is the lead laboratory for U.S. Department of Energy Advanced Vehicle Testing. Battelle Energy Alliance, LLC contracted with Intertek Testing Services, North America (ITSNA) to collect data on federal fleet operations as part of the Advanced Vehicle Testing Activity's Federal Fleet Vehicle Data Logging and Characterization study. The Advanced Vehicle Testing Activity study seeks to collect data to validate the utilization of advanced electric drive vehicle transportation. This report focuses on the Golden Gate National Recreation Area (GGNRA) fleet to identify daily operational characteristics of select vehicles and report findings on vehicle and mission characterizations to support the successful introduction of plug-in electric vehicles (PEVs) into the agencies' fleets. Individual observations of these selected vehicles provide the basis for recommendations related to electric vehicle adoption and whether a battery electric vehicle or plug-in hybrid electric vehicle (PHEV) (collectively PEVs) can fulfill the mission requirements. GGNRA identified 182 vehicles in its fleet, which are under the management of the U.S. General Services Administration. Fleet vehicle mission categories are defined in Section 4, and while the GGNRA vehicles conduct many different missions, only two (i.e., support and law enforcement missions) were selected by agency management to be part of this fleet evaluation. The selected vehicles included sedans, trucks, and sport-utility vehicles. This report will show that battery electric vehicles and/or PHEVs are capable of performing the required missions and providing an alternative vehicle for support vehicles and PHEVs provide the same for law enforcement, because each has a sufficient range for individual trips and time is available each day for charging to accommodate multiple trips per day. These charging events could occur at the vehicle home base, high-use work areas, or intermediately along routes that the vehicles frequently travel. Replacement of vehicles in the current fleet would result in significant reductions in the emission of greenhouse gases and petroleum use, while also reducing fuel costs. The San Francisco Bay Area is a leader in the adoption of PEVs in the United States. PEV charging stations, or more appropriately identified as electric vehicle supply equipment, located on the GGNRA facility would be a benefit for both GGNRA fleets and general public use. Fleet drivers and park visitors operating privately owned PEVs benefit by using the charging infrastructure. ITSNA recommends location analysis of the GGNRA site to identify the optimal placement of the electric vehicle supply equipment station. ITSNA recognizes the support of Idaho National Laboratory and ICF International for their efforts to initiate communication with the National Parks Service and GGNRA for participation in the study. ITSNA is pleased to provide this report and is encouraged by the high interest and support from the National Park Service and GGNRA personnel.

99 GENERAL AND MISCELLANEOUS↗

Theoretical Plasma Physics

Lattice Boltzmann algorithms are a mesoscopic method to solve problems in nonlinear physics which are highly parallelized – unlike the direction solution of the original problem. These methods are applied to both fluid and magnetohydrodynamic turbulence. By introducing entropic constraints one can enforce the positive definiteness of the distribution functions and so be able to simulate fluids at high Reynolds numbers without numerical instabilities. By introducing a vector distribution function for the magnetic field one can enforce the divergence free condition on the magnetic field automatically, without the need of divergence cleaning as needed in most direct numerical solutions of the resistive magnetohydrodynamic equations. The principal reason for the high parallelization of lattice Boltzmann codes is that they consist of a kinetic collisional relaxation step (which is purely local) followed by a simple shift of the relaxed data to neighboring lattice sites. In large eddy simulations, the closure schemes are highly nonlocal – the most famous of these schemes is that due to Smagorinsky. Under a lattice Boltzmann representation the Smagorinsky closure is purely local – being simply a particular moment on the perturbed distribution fucntions. After nonlocal fluid moment models were discovered to represent Landau damping, it was found possible to model these fluid models using an appropriate lattice Boltzmann algorithm. The close to ideal parallelization of the lattice Boltzmann codes permitted us to be Gordon Bell finalists on using the Earth Simulation in Japan. We have also been involved in the radio frequency propagation of waves into a tokamak and into a spherical overdense tokamak plasma. Initially we investigated the use of a quasi-optical grill for the launching of lower hybrid waves into a tokamak. It was found that the conducting walls do not prevent the rods from being properly irradiated, the overloading of the quasi-optical grill is not severe with the electric field only being about three times higher than in the ideal case. Moreover, the quasi-optical grill was significantly fewer structural elements that the multijunction grill. Nevertheless there has not been much interest from experimental fusion groups to implementing these structures. Hence we have returned to optimizing the multijunction grill so that the large number of coupling matrix elements can be efficiently evaluated using symmetry arguments. In overdense plasmas, the standard electromagnetic waves cannot propagate into the plasma center, but are reflected at the plasma edge. By optimizing mode conversion processes (in particular, the O-X-B wave propagation of Ordinary Mode converting to an Extraordinary mode which then converts into an electrostatic Bernstein wave) one can excite within the plasma an electrostatic Bernstein wave that does not suffer density cutoffs and is absorbed on the electron cyclotron harmonics. Finally we have started looking at other mesoscopic lattice algorithms that involve unitary collision and streaming steps. Because these algorithms are unitary they can be run on quantum computers when they become available – unlike their computational cousin of lattice Boltzmann which is a purely classical code. These quantum lattice gas algorithms have been tested successfully on exact analytic soliton collision solution. These calculations are hoped to be able to study Bose Einstein condensed atomic gases and their ground states in an optical lattice.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Federal Oversight of Hydrogen Systems

The application of hydrogen as an energy carrier has been expanding into industrial and transportation sectors enabling sustainable energy resources and providing a zero-emission energy infrastructure. The hydrogen supply infrastructure includes processes from production and storage, to transportation and distribution, to end use. Each portion of the hydrogen supply infrastructure is regulated by international, federal, state, and local entities. Regulations are enforced by entities which provide guidance and updates as necessary. While energy sources such as natural gas are currently regulated via the Code of Federal Regulations and United States Code, there might be some ambiguity as to which regulations are applicable to hydrogen and where regulatory gaps may exist. This report contains an overview of the regulations that apply to hydrogen, and those that may indirectly cover hydrogen as an energy carrier participating in a sustainable zero emission global energy system. As part of this effort, the infrastructure of hydrogen systems and regulation enforcement entities are defined, and a visual map and reference table are developed. This regulatory map and table can be used to identify the boundaries of federal oversight for each component of the hydrogen supply value chain which includes production, storage, distribution, and use.

08 HYDROGEN↗

Efficient Integration of Algebraic Constraint for Exponential Time Integration

This work is a continuation of a previous project where the high efficiency of exponential integration methods for DRMHD systems was demonstrated. Often algebraic constraints must also be enforced for DRMHD systems of interest. Straightforward application of exponential methods to DRMHD equations with constraints leads to prohibitively computationally expensive methods. In this work, we propose new exponential schemes that allow the constraints to be removed from evaluation of exponential matrix functions which drastically reduces computational complexity by eliminating the need to perform computations enforcing constraints in exponential calculations while still preserving a high order of accuracy and allowing for a large time step even when the problem is stiff. This idea is similar to the W-methods and is achieved by carefully designing a method with the desired order of accuracy, even with an incomplete Jacobian matrix used as an argument of exponential-like functions. The constraints are accounted for by including them in the evaluation of the right-hand-side forcing function of the spatially discretized system. We study performance of the new methods on test problems and outline future research directions that this work opens.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗