Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fixed grid”

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 379 records · Page 21

Crack Offset Measurement With the Projected Laser Target Device

The device and associated analysis methodology summarized in this report were developed for the purpose of estimating the size of discontinuities in the surface of the foam that covers the Space Shuttle External Tank. These surface offsets are thought to be due to subsurface cracks in the foam insulation. The mathematical analysis and procedure described here provide a method to quantity the dimensions of the crack offset in a direction perpendicular to the surface, making use of the projected laser target device (PLTD) tool and a laser line projector. The keys to the construction and use of the PLTD are the following geometrical design requirements: Laser dots are on a square grid: length on a side. Laser beams are perpendicular to projected surface. Beams are parallel out to the distance being projected. The PLTD can be used to (1) calibrate fixed cameras of unknown magnification and orientation (far-field solution); (2) provide equivalent calibration to multiple cameras, previously achieved only by the use of known target points, for example, in 3.D foreign-object debris tracking on a fixed launch platform; (3) compute scaling for conventional 2.D images, and depth of field for 3.D images (near-field solution); and (4) in conjunction with a laser line projector, achieve accurate measurements of surface discontinuity (cracks) in a direction perpendicular to the surface.

Source record↗

Methods and Results of a Search for Gravitational Waves Associated with Gamma-Ray Bursts Using the GEO 600, LIGO, and Virgo Detectors

In this paper we report on a search for short-duration gravitational wave bursts in the frequency range 64 Hz-1792 Hz associated with gamma-ray bursts (GRBs), using data from GEO600 and one of the LIGO or Virgo detectors. We introduce the method of a linear search grid to analyze GRB events with large sky localization uncertainties such as the localizations provided by the Fermi Gamma-ray Burst Monitor (GBM). Coherent searches for gravitational waves (GWs) can be computationally intensive when the GRB sky position is not well-localized, due to the corrections required for the difference in arrival time between detectors. Using a linear search grid we are able to reduce the computational cost of the analysis by a factor of O(10) for GBM events. Furthermore, we demonstrate that our analysis pipeline can improve upon the sky localization of GRBs detected by the GBM, if a high-frequency GW signal is observed in coincidence. We use the linear search grid method in a search for GWs associated with 129 GRBs observed satellite-based gamma-ray experiments between 2006 and 2011. The GRBs in our sample had not been previously analyzed for GW counterparts. A fraction of our GRB events are analyzed using data from GEO600 while the detector was using squeezed-light states to improve its sensitivity; this is the first search for GWs using data from a squeezed-light interferometric observatory. We find no evidence for GW signals, either with any individual GRB in this sample or with the population as a whole. For each GRB we place lower bounds on the distance to the progenitor, assuming a fixed GW emission energy of 10(exp -2)Stellar Mass sq c, with a median exclusion distance of 0.8 Mpc for emission at 500 Hz and 0.3 Mpc at 1 kHz. The reduced computational cost associated with a linear search grid will enable rapid searches for GWs associated with Fermi GBM events in the Advanced detector era.

methods↗

Power System Recovery Coordinated with (Non-)Black-Start Generators

Power restoration is an urgent task after a black-out, and recovery efficiency is critical when quantifying system resilience. Multiple elements should be considered to restore the power system quickly and safely. This paper proposes a recovery model to solve a direct-current optimal power flow (DCOPF) based on mixed-integer linear programming (MILP). Since most of the generators cannot start independently, the interaction between black-start (BS) and non-black-start (NBS) generators must be modeled appropriately. The energization status of the NBS is coordinated with the recovery status of transmission lines, and both of them are modeled as binary variables. Also, only after an NBS unit receives the cranking power through connected transmission lines, will it be allowed to participate in the following system dispatch. The amount of cranking power is estimated as a fixed proportion to the maximum generation capacity. The proposed model is validated on several test systems, as well as a 1393-bus representation system of the Puerto Rican electric power grid. Test results demonstrate how the recovery of NBS units and damaged transmission lines can be optimized, resulting in an efficient and well-coordinated recovery procedure.

Zhao, Meng↗

Numerical simulation of the effect of spatial disturbances on vortex asymmetry

The steady asymmetric vortex pattern observed on slender bodies of revolution at large angle of attack was investigated using fine-grid Navier-Stokes computations. The computed results demonstrate the marked asymmetry which has been observed in experiments. To obtain asymmetry, it was found essential to introduce a space-fixed time-invariant perturbation into the computation. The computational results suggest that vortex asymmetry is forced by amplification of small disturbances, such as those due to surface roughness, occurring within the body viscous boundary layer.

Degani, David↗

Calculation of unsteady linearized Euler flows in cascades using harmonically deforming grids

A method for calculating unsteady, inviscid, compressible flows in cascades is presented. Using the linearized Euler technique, the flow is decomposed into a steady or mean flow plus a harmonically varying small disturbance flow. The equations that describe the small disturbance flow are linear variable coefficient equations, and are solved using a pseudo-time time marching Lax-Wendroff technique. Unlike previous linearized methods, however, the solution is computed on a harmonically deforming computational grid that conforms to the motion of the vibrating airfoils. The mean flow and perturbation flow solutions are defined in the deforming coordinate system rather than in a coordinate system fixed in space. Hence, no extrapolation terms are required to implement the upwash boundary conditions at the airfoil surfaces, significantly improving the accuracy of the method. For transonic flow calculations, unsteady shock motions are modelled using shock capturing. The unsteady loads due to the shock motion are then seen as pressure impulses. Representative computational results are presented for transonic channel flows and subsonic and transonic cascade flows.

Hall, Kenneth C.↗

Supersonic vortex breakdown over a delta wing in transonic flow

The effects of freestream Mach number and angle of attack on the leading-edge vortex breakdown due to the terminating shock on a 65-degree, sharp-edged, cropped delta wing are investigated computationally, using the time-accurate solution of the laminar unsteady compressible full Navier-Stokes equations with the implicit upwind flux-difference splitting, finite-volume scheme. A fine O-H grid consisting of 125 x 85 x 84 points in the wrap-around, normal, and axial directions, respectively, is used for all the flow cases. Keeping the Reynolds number fixed at 3.23 x 10 exp 6, the Mach number is varied from 0.85 to 0.9 and the angle of attack is varied from 20 to 24 deg. The results show that, at 20-deg angle of attack, the increase of the Mach number from 0.85 to 0.9 results in moving the location of the terminating shock downstream. The results also show that, at 0.85 Mach number, the increase of the angle of attack from 20 to 24 deg results in moving the location of the terminating shock upstream. The results are in good agreement with the experimental data.

Kandil, Hamdy A.↗

Performance Tests of Snow-Related Variables Over the Tibetan Plateau and Himalayas Using a New Version of NASA GEOS-5 Land Surface Model that Includes the Snow Darkening Effect

The snow darkening effect, i.e. the reduction of snow albedo, is caused by absorption of solar radiation by absorbing aerosols (dust, black carbon, and organic carbon) deposited on the snow surface. This process is probably important over Himalayan and Tibetan glaciers due to the transport of highly polluted Atmospheric Brown Cloud (ABC) from the Indo-Gangetic Plain (IGP). This effect has been incorporated into the NASA Goddard Earth Observing System model, version 5 (GEOS-5) atmospheric transport model. The Catchment land surface model (LSM) used in GEOS-5 considers 3 snow layers. Code was developed to track the mass concentration of aerosols in the three layers, taking into account such processes as the flushing of the compounds as liquid water percolates through the snowpack. In GEOS-5, aerosol emissions, transports, and depositions are well simulated in the Goddard Chemistry Aerosol Radiation and Transport (GO CART) module; we recently made the connection between GOCART and the GEOS-5 system fitted with the revised LSM. Preliminary simulations were performed with this new system in "replay" mode (i.e., with atmospheric dynamics guided by reanalysis) at 2x2.5 degree horizontal resolution, covering the period 1 November 2005 - 31 December 2009; we consider the final three years of simulation here. The three simulations used the following variants of the LSM: (1) the original Catchment LSM with a fixed fresh snowfall density of 150 kg m-3 ; (2) the LSM fitted with the new snow albedo code, used here without aerosol deposition but with changes in density formulation and melting water effect on snow specific surface area, (3) the LSM fitted with the new snow albedo code as same as (2) but with fixed aerosol deposition rates (computed from GOCART values averaged over the Tibetan Plateau domain [Ion.: 60-120E; lat.: 20-50N] during March-May 2008) applied to all grid points at every time step. For (2) and (3), the same setting on the fresh snowfall density as in (1) was used.

Yasunari, Tppei J.↗

Enhancing the Flexibility of Generation of Run-Of-the-River Hydro Power Plants

Run-of-the-River Hydro Power Plants (ROR HPPs) output can be characterized as fixed and predictable output with low cost of both installation and operation. ROR HPP generation is proportional to the water flow and cannot be adjusted to provide inertia to the grid, which is critical to power systems stability. Hence ROR HPP can only participate in low profit energy markets and not in the high profit ancillary service market. With the advances in energy storage, they can help build a synthetic energy reservoir. The combined operation of ROR HPPs and energy storage are capable of both ramping and variable power output, allowing participation in the ancillary service markets. However, energy storage technologies such as supercapacitors and flywheels, which are power dense can quickly charge/discharge within a short duration. Energy dense storage technologies such as batteries can provide long duration energy injection while the maximum ramping rate is much lower than power dense energy storage. In order to provide long duration fast ramping capability, Hybrid Energy Storage Systems (HESSs) including both power and energy dense energy storage are proposed to be integrated to enable the flexible operation of ROR HPP. Power dense energy storage with quick burst of energy can assist in regulation, whereas energy dense storage can provide energy balance over larger time. In this paper, the properties of different types of energy storage including the ramping rate, affordable energy storage capacity, thermal management, etc. are analyzed. The size of each type of energy storage is determined to ensure optimal combination of HESS and ROR HPP. A generic Front End Controller (FEC) is developed and customized for ROR HPP to provide real time control and cohesive operation of multiple units. FEC serves as the interface between utility management systems and ROR HPPs to ensure their integration with the centralized controls. Using the Automatic Generation Control (AGC) signal received from grid operator, FEC optimizes the ROR HPP turbine and HESS so that the aggregated generation can satisfy the market needs. Based on the study, system operator can acquire extra inertia support for system stability, and the ROR HPP operators can enhance their revenues by participating in wider market avenues.

13 HYDRO ENERGY↗

Regulatory Considerations Associated with the Expanded Adoption of Distributed Solar

Increased adoption of distributed PV, and other forms of distributed generation, have the potential to affect utility-customer interactions, system costs recovery, and utility revenue streams. If a greater number of electricity customers choose to self-generate, demand for system power will decrease and utility fixed costs will have to be recovered over fewer kilowatt hours of sales. As such, regulators will need to determine the value and cost of additional distributed PV and determine the appropriate allocation of the costs and benefits among consumers. The potential for new business models to emerge also has implications for regulation and rate structures that ensure equitable solutions for all electricity grid users. This report examines regulatory tools and rate designs for addressing emerging issues with the expanded adoption of distributed PV and evaluates the potential effectiveness and viability of these options going forward. It offers the groundwork needed in order for regulators to explore mechanisms and ensure that utilities can collect sufficient revenues to provide reliable electric service, cover fixed costs, and balance cost equity among ratepayers -- while creating a value proposition for customers to adopt distributed PV.

14 SOLAR ENERGY↗

ECAR-6585 MARVEL Core Support Structure Analysis: Project #33526

Idaho National Laboratory (INL) is developing a small microreactor to produce electrical power utilizing a small nuclear core and stirling engines for power output under the Microreactor Application Research Validation and Evaluation (MARVEL) program. The MARVEL Reactor (Figure 1) will be installed in the Transient Reactor Test (TREAT) Facility. Walsh Engineering Services, PC (WES) has contracted with the INL to detail the final design. There are six main systems in the MARVEL, one of which is the Fuel and Core System (FCS). The FCS contains a subsystem–the Core Support Structures Subsystem (CSS)–that supports and orients each of the 36 nuclear fuel pins in the MARVEL Reactor System. The FCS (apart from the fixed Beryllium Oxide radial reflectors) is shown in light green and is depicted in more detail in Figure 1 and Figure 2. The main subject of this analysis is the lower grid plate, which is the principal item in the CSS which supports the fuel, see Figure 3. Additionally, the upper grid plate thermal expansion of the upper grid plate is presented in Appendix A.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Wind turbine generator interaction with diesel generators on an isolated power system

The results of a dynamic interaction investigation to characterize any disturbances caused by interfacing the Mod 0A wind turbine (150 kW configuration) with the Block Island utility diesel generator grid are reported. The tests were run when only two diesel generators were on line, and attention was given to power, frequency, and voltage time profiles. The interconnected system was examined in the start-up and synchronization phase, normal shutdown and cut-out of the wind turbine, during fixed pitch generation, and during variable pitch operation. Governors were installed on the diesel generators to accommodate the presence of wind-derived electricity. The blade pitch control was set to maintain power at 150 kW or below. Power and voltage transients were insignificant during start-up and shutdown, and frequency aberrations were within the range caused by load fluctuations. It is concluded that wind turbine generation can be successfully implemented by an isolated utility, even with a significant penetration to the total grid output.

Scott, G. W.↗

DSD2D-FLS 2010: Bdzil's 2010 DSD Code Base; Computing tb and Dn with Edits to Reduce the Noise in the Dn Field Near HE Boundaries

The full level-set function code, DSD3D, is fully described in LA-14336 (2007) [1]. This ASCI-supported, DSD code project was the last such LANL DSD code project that I was involved with before my retirement in 2007. My part in the project was to design and build the core DSD3D solver, which was to include a robust DSD boundary condition treatment. A robust boundary condition treatment was required, since for an important local “customer,” the only description of the explosives’ boundary was through volume fraction data. Given this requirement, the accuracy issues I had encountered with our “fast-tube,” narrowband, DSD2D solver, and the difficulty we had building an efficient MPI-parallel version of the narrowband DSD2D, I decided DSD3D should be built as a full level-set function code, using a totally local DSD boundary condition algorithm for the level-­set function, phi, which did not rely on the gradient of the level-­set function being one, |grad(phi)| = 1. The narrowband DSD2D solver was built on the assumption that |grad(phi)| could be driven to one, and near the boundaries of the explosive this condition was not being satisfied. Since the narrowband is typically no more than10*dx wide, narrowband methods are discrete methods with a fixed, non-­resolvable error, where the error is related to the thickness of the band: the narrower the band the larger the errors. Such a solution represents a discrete approximation to the true solution and does not limit to the solution of the underlying PDEs under grid resolution.The full level-­set function code, DSD3D, is fully described in LA-14336 (2007) [1]. This ASCI-­supported, DSD code project was the last such LANL DSD code project that I was involved with before my retirement in 2007. My part in the project was to design and build the core DSD3D solver, which was to include a robust DSD boundary condition treatment. A robust boundary condition treatment was required, since for an important local “customer,” the only description of the explosives’ boundary was through volume fraction data. Given this requirement, the accuracy issues I had encountered with our “fast-­tube,” narrowband, DSD2D solver, and the difficulty we had building an efficient MPI-parallel version of the narrowband DSD2D, I decided DSD3D should be built as a full level-­set function code, using a totally local DSD boundary condition algorithm for the level-­set function, phi, which did not rely on the gradient of the level-­set function being one, |grad(phi)| = 1. The narrowband DSD2D solver was built on the assumption that |grad(phi)| could be driven to one, and near the boundaries of the explosive this condition was not being satisfied. Since the narrowband is typically no more than10*dx wide, narrowband methods are discrete methods with a fixed, non-resolvable error, where the error is related to the thickness of the band: the narrower the band the larger the errors. Such a solution represents a discrete approximation to the true solution and does not limit to the solution of the underlying PDEs under grid resolution.

97 MATHEMATICS AND COMPUTING↗

Field evaluation of remote wind sensing technologies: Shore-based and buoy mounted LIDAR systems

In developing a national energy strategy, the United States has a number of objectives, including increasing economic growth, improving environmental quality, and enhancing national energy security. Wind power contributes to these objectives through the deployment of clean, affordable and reliable domestic energy. To achieve U.S. wind generation objectives, the Wind and Water Power Program within the Department of Energy’s (DOE) Office of Energy Efficiency and Renewable Energy (EERE) instituted the U.S. Offshore Wind: Removing Market Barriers Program in FY 2011. Accurate and comprehensive information on offshore wind resource characteristics across a range of spatial and temporal scales is one market barrier that needs to be addressed through advanced research in remote sensing technologies. There is a pressing need for reliable offshore wind-speed measurements to assess the availability of the potential wind energy resource in terms of power production and to identify any frequently occurring spatial variability in the offshore wind resource that may impact the operational reliability and lifetime of wind turbines and their components and to provide a verification program to validate the “bankability” of the output of these alternative technologies for use by finance institutions for the financing of offshore wind farm construction. The application of emerging remote sensing technologies is viewed as a means to cost-effectively meet the data needs of the offshore wind industry. In particular, scanning and buoy mounted LIDAR have been proposed as a means to obtain accurate offshore wind data at multiple locations without the high cost and regulatory hurdles associated with the construction of offshore meteorological towers. However; before these remote sensing technologies can be accepted the validity of the measured data must be evaluated to ensure their accuracy. The proposed research will establish a unique coastal ocean test-bed in the Mid-Atlantic for the evaluation of LIDAR-based wind measurement systems to validate the accuracy of remotely measured wind data in marine applications. Specifically, the test-bed will be utilized to systematically evaluate the capability of emerging scanning LIDAR and buoy mounted vertically profiling LIDAR by: (1) Evaluating a fixed scanning LIDAR against land-based 50 and 60 meter high meteorological masts fitted with research quality cup-vane and/or sonic anemometers; (2) Evaluating a buoy mounted vertically profiling LIDAR fixed on land and floating in a sheltered bay against a co-located 60 meter high meteorological mast fitted with a research quality cup-vane and/or sonic anemometers and the fixed scanning LIDAR; and (3) Offshore field evaluation of both LIDAR platforms through a comparison of the fixed scanning LIDAR data and data obtained by the buoy mounted LIDAR located 10 miles offshore. The proposed research will systematically validate Light Detection and Ranging (LIDAR) based wind measurement systems and assess the temporal and spatial variability of the offshore wind resource in the Mid-Atlantic east of New Jersey. The goal of the proposed project is to address the technical and commercial challenges of the offshore wind energy industry by validating and assessing cost-effective, over ocean wind resource characterization technologies. The objective is to systematically evaluate the capability of both scanning and vertically profiling LIDARs to accurately measure 3D wind fields through comparison with fixed met masts and intercomparison among LIDAR platforms. Once validated, data collected by both buoy mounted vertically profiling LIDARs and shore-based, pulsed horizontally scanning LIDARs can be used to accurately assess offshore wind resources and to quantify the spatial and temporal variability in the offshore wind fields. One of the fundamental research questions to be addressed in phase 1 is the assessment of various measurement and data processing schemes to retrieve accurate wind vectors in the marine environment over large sampling ranges (10 to 12 km) and varying atmospheric aerosol levels. Atmospheric conditions and aerosol content within the coastal ocean region of the Mid-Atlantic seaboard of the US can vary significantly over short time periods in response to frontal passages and extratropical and tropical low pressure system passage offshore of the coast. Since aerosols provide the scattering medium for the determination of LIDAR Doppler shifts in the atmosphere the accuracy and range of LIDAR derived velocity measurements as a function of variation in aerosol content in the marine environment is a key research question to be addressed. In phase 1, it is desired to capture as much variation in atmospheric conditions and aerosol content as possible. To this end, collocated measurements of LIDAR and standard anemometer wind fields will be captured by the project PIs over all four seasons and during specific events (e.g., coastal low pressure system passage) in year 1. Additionally, since the meteorological masts are permanent structures, additional events can be captured over the three year duration of the field research project. All research instruments are owned by Fishermen’s Energy and made available to the PIs though a lease agreement as part of the DOE grant. Energy Fishermen’s Energy will be responsible for the operation and maintenance of the scanning LIDAR and met mast anemometers. On a daily basis, environmental data and systems performance indicators will be transmitted from each measurement station to the Fishermen’s project team consisting of both in-house personnel and equipment manufacturer engineers. Data sets include compiled LIDAR files as well as data sets from ancillary sensors. Diagnostic parameters to be monitored include standard deviations of measured values, battery levels and charging systems output, and the operational status. Once data have been confirmed as complete and reliable, files will be transferred to the Garrad Hassan (a subcontractor to Fishermen’s Energy) for incorporation in to the validation database, which is accessible to other scientific team members. Data collection times and durations will be determined by the PI and Co-PIs in consultation with instrument engineers to ensure the capture of data representative of the expected range of mid-Atlantic atmospheric conditions (e.g., temperature, moisture, coastal low pressure systems, tropical systems, rain, snow, fog). The collection and processing of the data is a function of site specific measurement requirements (Kelley et.al. 2007; Hannon et.al. 2008). To determine the optimal profiles of wind speed and direction from the LIDAR radial velocities as a function of azimuth angle, rigorous estimates of the bias and random error of each radial velocity estimate are required. Lockheed Martin Coherent Technologies, Inc., under contract with Fishermen’s Energy, will provide analyses of raw and processed data using various scan patterns to determine optimal performance settings for the pulsed scanning LIDAR. Once optimized, appropriate processing and analyses techniques will be evaluated by Garrad Hassan for use in validating the accuracy of the LIDAR wind field measurements against the standard anemometer measurements from the meteorological masts. The most attractive capability of the scanning LIDAR is the ability to provide high spatial resolution observations in a three-dimensional volume which provides superior statistical accuracy due to the large number of samples obtained. Each radial scan provides measurements in 100 range gates over a distance of 10 to 12 km at an update rate of 5 to 10 Hz and rotation of 2.5° per second. Each rotation at a fixed azimuth requires 2.4 minutes. Depending on the number of azimuths desired a complete scan can take up to 10 minutes or longer to complete. Once collected the radial velocities are processed to produce vector wind velocity estimates based on a set of data distributed in angle and range around points of interest, typically a standard grid within the radial wind map. To calculate wind vectors over a limited spatial area of interest for the comparison of data with other measurement platforms a localized least-squares approach has been applied by Hannon et.al. (2008) and a Variation Assimilation (VAR) processing technique has been applied by Chan and al Assimilation (VAR) processing technique has been applied by Chan and Shao (2006). Additionally, Kelley et.al. (2007) applied a “stare” technique that fixed a scanning LIDAR in both azimuth and elevation angles to measure over collocated sampling volumes of the LIDAR and a 3D sonic anemometer mounted to a fixed mast. Although the stare technique is limited to wind directions aligned with the sampling radial of the LIDAR, it does provide a direct comparison of sampling volumes. Each of the processing techniques described above (and possibly others) will be evaluated to determine the validity of the LIDAR derived wind fields in the marine environment. Numerical methods such as linear regression and comparison of probability density functions of wind fields measured by each instrument platform will be used to assess the processing techniques. Linear regression has the advantage of directly evaluating corresponding pairs of wind data measured by each instrument and can lend insight into deviations and bias between instruments as a function of wind speed. Assessment of the coherency between the probability density function of the wind measured by each instrument provides insight into processes that may not be accurately resolved by each instrument at specific frequencies. Once processed and assessed the most appropriate technique will be utilized to provide valid wind measurements from the pulsed scanning LIDAR. A detailed analysis of the measurement data from the LIDAR and the three meteorological towers and comparison of the coastal wind characteristics from the different systems will be performed by the CO_PI at NREL. This analysis will evaluate how the WT LIDAR performance and measurement of the wind characteristics vary with distance from the LIDAR and by atmospheric conditions; using the tower measurements at different distances (4.8, 9.6, and 19 km) from the LIDAR as a reference. The comparative analysis will include, to the extent possible, evaluation of parameters such as wind speed and direction distributions, wind shear, turbulence intensity and their variations by atmospheric conditions, month or season, and time of day.

17 WIND ENERGY↗

s-Step Krylov Subspace Methods as Bottom Solvers for Geometric Multigrid

Geometric multigrid solvers within adaptive mesh refinement (AMR) applications often reach a point where further coarsening of the grid becomes impractical as individual sub domain sizes approach unity. At this point the most common solution is to use a bottom solver, such as BiCGStab, to reduce the residual by a fixed factor at the coarsest level. Each iteration of BiCGStab requires multiple global reductions (MPI collectives). As the number of BiCGStab iterations required for convergence grows with problem size, and the time for each collective operation increases with machine scale, bottom solves in large-scale applications can constitute a significant fraction of the overall multigrid solve time. In this paper, we implement, evaluate, and optimize a communication-avoiding s-step formulation of BiCGStab (CABiCGStab for short) as a high-performance, distributed-memory bottom solver for geometric multigrid solvers. This is the first time s-step Krylov subspace methods have been leveraged to improve multigrid bottom solver performance. We use a synthetic benchmark for detailed analysis and integrate the best implementation into BoxLib in order to evaluate the benefit of a s-step Krylov subspace method on the multigrid solves found in the applications LMC and Nyx on up to 32,768 cores on the Cray XE6 at NERSC. Overall, we see bottom solver improvements of up to 4.2x on synthetic problems and up to 2.7x in real applications. This results in as much as a 1.5x improvement in solver performance in real applications.

97 MATHEMATICS AND COMPUTING↗

SU-E-T-631: Commissioning and Comprehensive Evaluation of the ArcCHECK Cylindrical Diode Array for VMAT QA

Purpose: Commissioning and comprehensive evaluation of ArcCHECK phantom for dosimetry of VMAT QA, using 6MV photon beam with and without the flattening filter. Methods: ArcCHECK was evaluated for response dependency on linac dose rate, instantaneous dose rate, radiation field size, beam angle and couch insertion. Scatter dose characterization, consistency and symmetry of response, dosimetric accuracy of fixed aperture arcs and clinical VMAT plans were investigated. Measurements were done using TrueBeam™ STx accelerator (Console version 1.6) with a 6 MV beam with and without flattening filter. Reference dose-grids were calculated using Eclipse TPS Analytical Anisotropic Algorithm (AAA version 10.0.39). Planned doses were calculated using symmetric 2mm 3D dose grids with 4 degree angular resolution defaulted to each control point. Gamma evaluations were performed in absolute dose mode, with default normalization to maximum dose in the curved plane and a low dose threshold of 10% to restrict the analysis to clinically relevant areas. Global and local gamma indices at 3mm/3% and 2mm/2% level were computed using SNC software (version 6.0). Results: Results of gamma analysis demonstrated an overall agreement between ArcCHECK measured and TPS calculated reference doses. Field size dependency was within 0.5% of the reference. Dose-rate based dependency was well within 1% of the TPS reference and the angular dependency was ±3% of the reference, as tested for BEV angles. At the level of 3%/3mm, narrow and wide open arcs as well as clinical VMAT cases demonstrated high level of dosimetry accuracy in global gamma passing rates for both 6X and 6F beams. At the level of 2%/2mm two VMAT cases involving the narrow heavily modulated arcs demonstrated lower passing rates. Conclusion: ArcCHECK phantom with latest software and hardware upgrades is suitable for VMAT QA. For higher sensitivity of 2%/2mm gamma analysis, we intend to use it as one of the VMAT QA evaluation metrics.

07 ISOTOPE AND RADIATION SOURCES↗

Renewable Energy Potential Model: Hawaii Geothermal Supply Curves

This dataset extends the development of the Renewable Energy Potential (reV) model to include geothermal energy, with a specific focus on Hawaii. Provided here are the results of two scenarios that were modeled for geothermal energy in Hawaii: binary enhanced geothermal systems (EGS) at a depth of 2.5 km and hydrothermal binary systems at a depth of 1.5 km. The resource data for both scenarios were derived from Lautze and Haskins (2024) using an exponential method. The PFA probability of heat map was used as a look up table for which temperature gradient to use (Lautze and Haskins, 2024). The dataset provides geospatial and techno-economic details for evaluating geothermal energy potential. It includes spatial coordinates, estimated capacity factors, developable area, resource potential, and annual energy production metrics. Economic details such as levelized cost of electricity (LCOE), site development costs, transmission costs, and fixed-charge rates are also included. The reV model, originally developed for wind and solar energy, incorporates these variables to evaluate deployment constraints related to land use, environmental and cultural factors, and grid integration.

15 GEOTHERMAL ENERGY↗

Multiobjective Optimization of Regional and Teleseismic Data to Constrain the Source of the 12 September 2016 Mw 5.4 Earthquake in South Korea

In this paper, we propose a new approach for source inversion with automatic criteria for the station selection process and using a multiobjective optimization scheme combining teleseismic body-wave and regional surface-wave observations. To illustrate the approach, we determine the source parameters of the 12 September 2016 M w 5.4 South Korea earthquake, as well as its foreshock (same day) and its main aftershock on 19 September 2016. Available teleseismic (short period and broadband) and regional (broadband) seismograms from the Incorporated Research Institutions for Seismology International–Federation of Digital Seismograph Networks (IRIS-FDSN) webservice are selected (see Data and Resources). First, all teleseismic data with signal-to-noise ratio (SNR) above 2 in the 0.3–1.5 Hz frequency band are selected and a cepstrum function is blindly applied to help detect reflected wave arrivals inside a P coda. This method allows the identification of records with weak P-coda and clear depth phase (pP, sP) arrivals and can be used to select an optimal set of teleseismic seismograms to be used for the inversion. Second, regional seismograms are selected according to the fit between their observed and theoretical dispersion curves extracted from local 1D velocity model. Depending on the epicentral distance, two frequency bands are proposed for regional data: 15–50 s for epicentral distances smaller than 2° and 30–80 s otherwise. Finally, we explore the space parameters from exhaustive grid searches for the two objective functions Tf and Rf for teleseismic and regional data, respectively. We show that the seismic moment parameters can be exclusively derived from Rf and can be fixed for Tf. However, the depth parameter is not sensitive to Rf and can be simultaneously inverted from both. The selection of the optimal solution satisfying both cost functions uses the optimum Pareto front that allows a decrease in epistemic uncertainty. Lastly, we discuss the advantage of this approach for improving the characterization of moderate earthquake sources in areas with limited regional instrumentation.

58 GEOSCIENCES↗

Gridded uncertainty in fossil fuel carbon dioxide emission maps, a CDIAC example

Due to a current lack of physical measurements at appropriate spatial and temporal scales, all current global maps and distributions of fossil fuel carbon dioxide (FFCO2) emissions use one or more proxies to distribute those emissions. These proxies and distribution schemes introduce additional uncertainty into these maps. This paper examines the uncertainty associated with the magnitude of gridded FFCO2 emissions. This uncertainty is gridded at the same spatial and temporal scales as the mass magnitude maps. This gridded uncertainty includes uncertainty contributions from the spatial, temporal, proxy, and magnitude components used to create the magnitude map of FFCO2 emissions. Throughout this process, when assumptions had to be made or expert judgment employed, the general tendency in most cases was toward overestimating or increasing the magnitude of uncertainty. The results of the uncertainty analysis reveal a range of 4–190 %, with an average of 120 % (2 σ ) for populated and FFCO2-emitting grid spaces over annual timescales. This paper also describes a methodological change specific to the creation of the Carbon Dioxide Information Analysis Center (CDIAC) FFCO2 emission maps: the change from a temporally fixed population proxy to a temporally varying population proxy.

54 ENVIRONMENTAL SCIENCES↗