Search NASA⌕ Search

SEARCH · Search NASA

Results for “inverter availability”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Utilities Perspective on Protection Challenges with High IBR Penetration

Utilities have seen rapid increase in Solar, Wind, and Battery Energy Storage Resources that interconnect to their electric system through Inverters. Inverter-Based Resources (IBRs) have fault current characteristics that are unlike the fault current response of traditional rotating-machine-based generators, which is well known and repeatable. IBR’s non-traditional fault current behavior is due to the IBR control scheme, which is configured to provide a clean AC output but also protect the inverter’s sensitive power electronics devices from damage, one source of which is overcurrent. This results in low fault current magnitude, low or no negative sequence current injection, the variability of sequence component currents, the variability of voltage with respect to current angles, and the lack of inertia. The control scheme also results in a fault current response that can vary between manufactures and between models of the same manufacturer. High penetration of IBRs can adversely affect the protection schemes applied in areas with high penetration of IBRs. With the proliferation of IBRs, utilities are finding out that conventional protection schemes are not adequately equipped to protect the electric systems. This is mainly because the existing protection elements and practices have been designed based on the fault current response of conventional rotating machines. In several cases, the available literature does not provide any clear solution for the issues when the protection scheme does not operate properly near IBRs. This presentation identifies various protection challenges due to IBRs that industry is facing, from the utility perspective. Instead of facing on one issue, we are looking broadly on all the challenges that system protection has experienced with high penetration of IBRs. Based on the IBR response from various utilities during real fault events and gathering perspective from different utility SMEs via questionnaire, the presentation summarizes on gathered data and internal experiences.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Utility Perspective on Protection Challenges with High IBR Penetration [Slides]

Utilities have seen rapid increase in Solar, Wind, and Battery Energy Storage Resources that interconnect to their electric system through Inverters. Inverter-Based Resources (IBRs) have fault current characteristics that are unlike the fault current response of traditional rotating-machine-based generators, which is well known and repeatable. IBR’s non-traditional fault current behavior is due to the IBR control scheme, which is configured to provide a clean AC output but also protect the inverter’s sensitive power electronics devices from damage, one source of which is overcurrent. This results in low fault current magnitude, low or no negative sequence current injection, the variability of sequence component currents, the variability of voltage with respect to current angles, and the lack of inertia. The control scheme also results in a fault current response that can vary between manufactures and between models of the same manufacturer. High penetration of IBRs can adversely affect the protection schemes applied in areas with high penetration of IBRs. With the proliferation of IBRs, utilities are finding out that conventional protection schemes are not adequately equipped to protect the electric systems. This is mainly because the existing protection elements and practices have been designed based on the fault current response of conventional rotating machines. In several cases, the available literature does not provide any clear solution for the issues when the protection scheme does not operate properly near IBRs. This report identifies various protection challenges due to IBRs that industry is facing, from the utility perspective. Instead of facing on one issue, this report looks broadly on all the challenges that system protection has experienced with high penetration of IBRs. Based on the IBR response from various utilities during real fault events and gathering perspective from different utility SMEs via questionnaire, the report summarizes on gathered data and internal experiences.

14 SOLAR ENERGY↗
DOE OSTI2022

Field Experience Detecting PV Underperformance in Real Time Using Existing Instrumentation

Maintenance at large-scale photovoltaic plants employs a mix of preventative and corrective maintenance practices. Large outages, such as an inverter tripping offline, are often easy to detect. More subtle sub-inverter faults and failures can accumulate and go unnoticed for months or years. A software-based fault detection method has been developed to analyze commonly measured data from large-scale PV plants for more timely detection of subtle underperformance. The method has been demonstrated on eight datasets from large-scale plants with high accuracy of detection. Results are validated using aerial infrared scanning. String outages are detected with a true positive rate of 73 percent and tracker issues are detected with a true positive rate of 88 percent. The developed method can be uniformly applied to photovoltaic plants across a range of scales and configurations to assess performance, quickly detect underperformance, and determine the source and location of failures. The results inform and improve operations and maintenance at PV plants, ultimately aiding in improved affordability, reliability, availability, and resiliency of solar electricity.

14 SOLAR ENERGY↗
DOE OSTI2024

Assessing Inverter-Based Resources Modeling Gaps in Commonly Used Short-Circuit Programs

Future bulk power systems are expected to operate with a high penetration of inverter-based resources (IBR). This is anticipated to fundamentally change the system’s short-circuit behaviors and demand a change in existing transmission line protection settings or schemes. Phasor-domain short circuit programs (e.g., ASPEN, CAPE, and CYME) remain the mainstream tools for developing protection settings. However, the system protection community has shared concerns about the accuracy of IBR models in short-circuit programs. This report compares the performance of existing generic IBR models in short-circuit programs with detailed electro-magnetic transient (EMT) models provided by several different IBR manufacturers in PSCAD, to quantify the short-circuit program’s IBR model accuracy, assess their ability to replicate individual IBR controls such as negative sequence current injections, and identify specific gaps in the existing generic phasor-domain IBR models that are currently available.

42 ENGINEERING↗
NASA NTRS2000

Airborne Aerosol Closure Studies During PRIDE

The Puerto Rico Dust Experiment (PRIDE) was conducted during June/July of 2000 to study the properties of Saharan dust aerosols transported across the Atlantic Ocean to the Caribbean Islands. During PRIDE, the NASA Ames Research Center six-channel (380 - 1020 nm) airborne autotracking sunphotometer (AATS-6) was operated aboard a Piper Navajo airplane alongside a suite of in situ aerosol instruments. The in situ aerosol instrumentation relevant to this paper included a Forward Scattering Spectrometer Probe (FSSP-100) and a Passive Cavity Aerosol Spectrometer Probe (PCASP), covering the radius range of approx. 0.05 to 10 microns. The simultaneous and collocated measurement of multi-spectral aerosol optical depth and in situ particle size distribution data permits a variety of closure studies. For example, vertical profiles of aerosol optical depth obtained during local aircraft ascents and descents can be differentiated with respect to altitude and compared to extinction profiles calculated using the in situ particle size distribution data (and reasonable estimates of the aerosol index of refraction). Additionally, aerosol extinction (optical depth) spectra can be inverted to retrieve estimates of the particle size distributions, which can be compared directly to the in situ size distributions. In this paper we will report on such closure studies using data from a select number of vertical profiles at Cabras Island, Puerto Rico, including measurements in distinct Saharan Dust Layers. Preliminary results show good agreement to within 30% between mid-visible aerosol extinction derived from the AATS-6 optical depth profiles and extinction profiles forward calculated using 60s-average in situ particle size distributions and standard Saharan dust aerosol refractive indices published in the literature. In agreement with tendencies observed in previous studies, our initial results show an underestimate of aerosol extinction calculated based on the in situ size distributions relative to the extinction obtained from the sunphotometer measurements. However, a more extensive analysis of all available AATS-6 and in situ size distribution data is necessary to ascertain whether the preliminary results regarding the degree of extinction closure is representative of the entire range of dust conditions encountered in PRIDE. Finally, we will compare the spectral extinction measurements obtained in PRIDE to similar data obtained in Saharan dust layers encountered above the Canary Islands during ACE-2 (Aerosol Characterization Experiment) in July 1997. Thus, the evolution of Saharan dust spectral properties during its transport across the Atlantic can be investigated, provided the dust origin and microphysical properties are found to be comparable.

Redemann, Jens↗
DOE OSTI2022

Making Invisible Visible: Data-Driven Seismic Inversion With Spatio-Temporally Constrained Data Augmentation

Deep learning and data-driven approaches have shown great potential in scientific domains. The promise of data-driven techniques relies on the availability of a large volume of high-quality training datasets. Due to the high cost of obtaining data through expensive physical experiments, instruments, and simulations, data augmentation techniques for scientific applications have emerged as a new direction for obtaining scientific data recently. However, existing data augmentation techniques originating from computer vision yield physically unacceptable data samples that are not helpful for the domain problems that we are interested in. In this article, we develop new data augmentation techniques based on convolutional neural networks. Specifically, our generative models leverage different physics knowledge (such as governing equations, observable perception, and physics phenomena) to improve the quality of the synthetic data. To validate the effectiveness of our data augmentation techniques, we apply them to solve a subsurface seismic full-waveform inversion using simulated CO 2 leakage data. Our interest is to invert for subsurface velocity models associated with very small CO 2 leakage. We validate the performance of our methods using comprehensive numerical tests. Here via comparison and analysis, we show that data-driven seismic imaging can be significantly enhanced by using our data augmentation techniques. Particularly, the imaging quality has been improved by 15% in test scenarios of general-sized leakage and 17% in small-sized leakage when using an augmented training set obtained with our techniques.

58 GEOSCIENCES↗
NASA NTRS1990

Recursive dynamics for flexible multibody systems using spatial operators

Due to their structural flexibility, spacecraft and space manipulators are multibody systems with complex dynamics and possess a large number of degrees of freedom. Here the spatial operator algebra methodology is used to develop a new dynamics formulation and spatially recursive algorithms for such flexible multibody systems. A key feature of the formulation is that the operator description of the flexible system dynamics is identical in form to the corresponding operator description of the dynamics of rigid multibody systems. A significant advantage of this unifying approach is that it allows ideas and techniques for rigid multibody systems to be easily applied to flexible multibody systems. The algorithms use standard finite-element and assumed modes models for the individual body deformation. A Newton-Euler Operator Factorization of the mass matrix of the multibody system is first developed. It forms the basis for recursive algorithms such as for the inverse dynamics, the computation of the mass matrix, and the composite body forward dynamics for the system. Subsequently, an alternative Innovations Operator Factorization of the mass matrix, each of whose factors is invertible, is developed. It leads to an operator expression for the inverse of the mass matrix, and forms the basis for the recursive articulated body forward dynamics algorithm for the flexible multibody system. For simplicity, most of the development here focuses on serial chain multibody systems. However, extensions of the algorithms to general topology flexible multibody systems are described. While the computational cost of the algorithms depends on factors such as the topology and the amount of flexibility in the multibody system, in general, it appears that in contrast to the rigid multibody case, the articulated body forward dynamics algorithm is the more efficient algorithm for flexible multibody systems containing even a small number of flexible bodies. The variety of algorithms described here permits a user to choose the algorithm which is optimal for the multibody system at hand. The availability of a number of algorithms is even more important for real-time applications, where implementation on parallel processors or custom computing hardware is often necessary to maximize speed.

Jain, A.↗
DOE OSTI2021

Data-driven Optimal Control Strategy for Virtual Synchronous Generator via Deep Reinforcement Learning Approach

This paper aims at developing a data-driven optimal control strategy for virtual synchronous generator (VSG) in the scenario where no expert knowledge or requirement for system model is available. Firstly, the optimal and adaptive control problem for VSG is transformed into a reinforcement learning task. Specifically, the control variables, i.e., virtual inertia and damping factor, are defined as the actions. Meanwhile, the active power output, angular frequency and its derivative are considered as the observations. Moreover, the reward mechanism is designed based on three preset characteristic functions to quantify the control targets: (1) maintaining the deviation of angular frequency within special limits; (2) preserving well-damped oscillations for both the angular frequency and active power output; (3) obtaining slow frequency drop in the transient process. Next, to maximize the cumulative rewards, a decentralized deep policy gradient algorithm, which features model-free and faster convergence, is developed and employed to find the optimal control policy. With this effort, a data-driven adaptive VSG controller can be obtained. By using the proposed controller, the inverter-based distributed generator can adaptively adjust its control variables based on current observations to fulfill the expected targets in model-free fashion. Finally, simulation results validate the feasibility and effectiveness of the proposed approach.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Advanced Distributed Wind Turbine Controls Series: Part 1-Flatirons Campus Model Overview – Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL)

Wind turbines are typically deployed to provide energy, reduce diesel-fuel consumption, reduce carbon emissions, and reduce costs for energy and fuel transportation. However, in addition to solely providing energy to the power system, wind turbines contain rotating masses and inverter-based controls that can enable various reliability and resilience services through advance controls. As part of the Microgrids, Infrastructure Resilience, and Advanced Controls Launchpad (MIRACL), it is demonstrated that advanced wind turbine controls can be employed to support higher contributions of wind, and to demonstrate ways that wind can play a role in supporting grid stability in islanded or grid-connected configurations. This paper documents models of various subsystem comprising a portion of NREL's Flatirons campus that will be used in three subsequent reports to demonstrate capabilities of advanced wind turbine controls. The series of reports will detail advanced capabilities of distributed wind turbines to provide support to isolated grids, distribution grids, and microgrids. We developed models to simulate a wind turbine (600 kW), solar PV (430 kW), battery energy storage system (1 MW/1MWh), a diesel generator (2 MW) and various types of loads (critical, dynamic). The model of the subsystems in MATLAB/Simulink are validated with available data from real-world components on NREL's Flatirons Campus. These validated models can be configured for various studies including four MIRACL use cases: 1) isolated grids, 2) microgrids, and 3) behind-the-meter, and 4) front-of-the-meter wind turbine deployments.

17 WIND ENERGY↗
DOE OSTI2026

Grid and Wind Transient Ride Through Assessment of Type-5 Synchronous Wind Turbine

The integration of the inverter-based renewable energy resources is increasing and replacing conventional synchronous generators in the evolving bulk power system (BPS), impacting the system inertia along the way. Also, the limited over-current capability of power electronic converters continues to deteriorate the grid strength of the evolving power systems. Type-5 wind turbines have the potential to mitigate these problems because this configuration interfaces with the grid through a synchronous generator; hence, its operation and consequent grid impacts might mimic the conventional synchronous generator. The drivetrain of this type of turbines, however, is different from the more common Type-3 and Type-4 turbines. A hydraulic torque converter and gearbox with torque limiting feature are integral parts of a Type-5 wind turbine drivetrain, to isolate the wind speed variations at the wind turbine rotor from the grid connected generator side. This supports the synchronous operation of the generator. However, high fidelity models of Type-5 wind turbine drivetrain including these core components are not openly and widely available for grid integration and transient stability studies. Moreover, Type-5 turbines are not commonly deployed and hence its power generation and frequency stability characteristics are not widely publicized. This hinders appropriate assessment of Type-5 wind power plant's contribution to bulk grid resilience. This work provides transient stability assessment of a standalone Type-5 wind turbine. Normal operation and grid transient studies are performed to see the step load change and wind speed disturbance impacts on the Type-5 drivetrain and generation dispatch. Preliminary results show that Type-5 wind turbines can support load changes of between 0.6 p.u to 1.15 p.u. Beyond this point, the synchronization of the generator starts to deviate from its set point.

17 - WIND ENERGY↗
DOE OSTI2024

Next-Generation, High-temperature, High-frequency, High-efficiency, High-power-density Traction System

To meet performance and reliability requirements necessary for broader adoption of electric drive vehicles, the Electrical and Electronics Technical Team of the U.S. Drive partnership has established aggressive design goals for next-generation electric vehicle drivetrains. Specifically, the 2025 roadmap stipulates a 100 kW/L power density target and a $\$$2.7/kW cost target for power electronics, in addition to high-voltage operation (i.e., greater than 800 VDC). The additional targets for traction motor and the overall system performance impose further challenges on the power electronics design. For example, many high specific power machines have reduced iron content, and therefore reduced intrinsic filtering, thus requiring the inverter to supply a low-distortion drive current. These machines also typically have a high pole count, thus requiring drive current at a higher electrical frequency. Other motors, such as brush-less dc and switch reluctance machines, require a carefully-shaped, non-sinusoidal drive current (Yang, Shang, Brown, & Krishnamurthy, 2015), (Zhang, Bowman, O'Connel, & Haran, 2018), (Anderson, et al., 2018). Two- and three-level inverter topologies are the conventional framework for the power electronics design of the drivetrain, and some demonstrations have shown recent progress towards addressing cost, power density and efficiency goals (Gurpinar & Ozpineci, 2018), (Zhu, Kim, Chen, Erickson, & Maksimović, 2018), (Deshpande, Chen, Narayanasamy, Sathyanarayanan, & Luo, 2018), (Alizadeh, et al., 2019). However, an unconventional approach may be necessary to take the dramatic leap in power density necessitated by the roadmap—while simultaneously addressing the other system needs. Therefore, this project leverages the flying capacitor multilevel (FCML) topology, together with a scalable, modular approach, to address these needs. This type of hybrid converter has several advantages: lower voltage (i.e., less than 300 V) transistors can be used, energy-dense capacitors process most of the power, and the output current waveform is multilevel and exhibits a frequency multiplying effect—in other words, the output has reduced dv/dt and filtering requirements for the same high voltage dc bus. For example, in an electric vehicle with an 800 V bus, a 10-level FCML could leverage 100 V, commercially available GaN devices switching at 115 kHz to produce a ~1 MHz switching waveform (modulated according to the motor drive requirements) with one ninth of the dv/dt of a two-level converter. Prior work has already demonstrated promising performance and gravimetric power density figures for more electric aircraft applications (Pallo, Foulkes, Modeer, Coday, & Pilawa-Podgurski, 2018). This project leverages lessons learned to achieve the volumetric power density of 100 kW/L by employing advanced liquid cooling, address the 300,000 mile reliability challenge with redundant design, topology failure studies and online health monitoring, and reduce costs to $\$$2.7/kW through the use of low-cost GaN devices, modular converter assemblies, and modest modifications to traditional manufacturing methods. The project involved several hardware designs, each achieving increasing performance. At the conclusion of the project, a volumetric power density of 380 kW/L was achieved, in a 800V dc-ac converter, greatly surpassing even the aggressive target goal.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

The Dependence of the Type Ia Supernova Host Bias on Observation or Fitting Technique

More luminous Type Ia supernovae prefer less massive hosts and regions of higher star formation. This correlation is inverted during width–color–luminosity light-curve standardization resulting in step-like biases of distance measurements with respect to host properties. Using the PMAS/PPak Integral-field Supernovahosts COmpilation (PISCO) supernova host sample and Sloan Digital Sky Survey, Galaxy Evolution Explorer, and Two Micron All Sky Survey photometry, we compare host stellar mass and specific star-formation rate (sSFR) from different observation methods, including local versus global, and fitting techniques to measure their impact on the host step biases. Mass-step measurements for all our mass samples are consistent within a 1σ significance from –0.03 ± 0.02 mag to –0.04 ± 0.02 mag. Including or excluding UV information had no effect on measured mass-step size or location. sSFR step sizes are more significant than mass-step measurements and varied from 0.05 ± 0.03 mag (Hα) and 0.06 ± 0.02 mag (UV) for a 51 host sample. The sSFR step location is influenced by the mass sample used to normalize star formation and by sSFR tracer choice. The step size is reduced to 0.04 ± 0.03 mag when using all available 73 hosts with Hα measurements. This 73 PISCO host subsample overall lacked a clear step signal, but here we are searching for whether different choices of mass or sSFR estimation can create a step signal. We find no evidence that different observation or fitting techniques choices can create a distance measurement step in either mass or sSFR.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1999

A Study of Nucleate Boiling with Forced Convection in Microgravity

The ultimate objective of basic studies of flow boiling in microgravity is to improve the understanding of the processes involved, as manifested by the ability to predict its behavior. This is not yet the case for boiling heat transfer even in earth gravity, despite the considerable research activity over the past 30 years. The elements that constitute the nucleate boiling process - nucleation, growth, motion, and collapse of the vapor bubbles (if the bulk liquid is subcooled) - are common to both pool and flow boiling. It is well known that the imposition of bulk liquid motion affects the vapor bubble behavior relative to pool boiling, but does not appear to significantly influence the heat transfer. Indeed, it has been recommended in the past that empirical correlations or experimental data of pool boiling be used for design purposes with forced convection nucleate boiling. It is anticipated that such will most certainly not be possible for boiling in microgravity, based on observations made with pool boiling in microgravity. In earth gravity buoyancy will act to remove the vapor bubbles from the vicinity of the heater surface regardless of how much the imposed bulk velocity is reduced, depending, of course, on the geometry of the system. Vapor bubbles have been observed to dramatically increase in size in pool boiling in microgravity, and the heat flux at which dryout took place was reduced considerably below what is generally termed the critical heat flux (CHF) in earth gravity, depending on the bulk liquid subcooling. However, at heat flux levels below dryout, the nucleate pool boiling process was enhanced considerably over that in earth gravity, in spite of the large vapor bubbles formed in microgravity and perhaps as a consequence. These large vapor bubbles tended to remain in the vicinity of the heater surface, and the enhanced heat transfer appeared to be associated with the presence of what variously has been referred to as a liquid microlayer between the bubble and the heater surface. The enhancement of the boiling process with low velocities in earth gravity for those orientations producing the formation of a liquid macrolayer described above, accompanied by "sliding" vapor bubbles, has been demonstrated. The enhancement was presented as a function of orientation, subcooling, and heated length, while a criterion for the heat transfer for mixed natural/forced convection nucleate boiling was given previously. A major unknown in the prediction and application of flow boiling heat transfer in microgravity is the upper limit of the heat flux for the onset of dryout (or critical heat flux - CHF), for given conditions of fluid-heater surfaces, including geometry, system pressure and bulk liquid subcooling. It is clearly understood that the behavior in microgravity will be no different than on earth with sufficiently high flow velocities, and would require no space experimentation. However, the boundary at which this takes place is still an unknown. Previous results of CHF measurements were presented for low velocity flow boiling at various orientations in earth gravity as a function of flow velocity and bulk liquid subcooling, along with preliminary measurements of bubble residence times on a flat heater surface. This showed promise as a parameter to be used in modeling the CHF, both in earth gravity and in microgravity. The objective of the work here is to draw attention to and show results of current modeling efforts for the CHF, with low velocities in earth gravity at different orientations and subcoolings. Many geometrical possibilities for a heater surface exist in flowing boiling, with boiling on the inner and outer surfaces of tubes perhaps being the most common. If the vapor bubble residence time on and departure size from the heater surface bear a relationship to the CHF, as results to be given indicate, it is important that visualization of and access to vapor bubble growth be conveniently available for research purposes. In addition, it is desirable to reduce the number of variables as much as possible in a fundamental study. These considerations dictated the use of a flat heater surface, which is rectangular in shape, 1.91 cm by 3.81 cm (0.75 x 1.5 inches), consisting either of a 400 Angstrom thick semi-transparent gold film sputtered on a quartz substrate which serves simultaneously as a heater and a resistance thermometer, or a copper substrate of the same size. The heater substrate is a disc which can be rotated so that the heated length in the flow direction can be changed from 1.91 to 3.81 cm (0.75 to 1.5 inches). The fluid is R-113, and the velocities can be varied between 0.5 cm/s and 60 cm/s. For a sufficiently low velocity the CHF can be modeled reasonably well at various orientations by the correlation for pool boiling corrected for the influence of bulk liquid subcooling, multiplied by the square root of q, the angle relative to horizontal. This arises from equating buoyancy and drag forces in the inverted positions where the vapor bubbles are held against the heater surface as they slide. A distortion of the measurements relative to pool boiling occurs as the flow velocity increases. In modeling this effect at different levels of subcooling it appeared appropriate to estimate the volumetric rate of vapor generation, using measurements of bubble frequency (or residence time), void fraction and average bubble boundary layer thickness. These were determined with the use of a platinum hot wire probe 0.025 mm in diameter by 1.3 mm long, applying a constant current to distinguish between contact with liquid or vapor. Two-dimensional spatial variations are obtained with a special mechanism to resolve displacements in increments of 0.025 mm. From such measurements it was determined that the fraction of the surface heat transfer resulting in evaporation varies inversely with the subcooling correction factor for the CHF. The measured inverse bubble residence time is normalized relative to that predicted for an infinite horizontal flat plate at the CHF, and is correlated well with the CHF normalized relative to that for pool boiling, for various orientation angles and subcooling levels. This correspondence is then combined with a normalizing factor for the energy flux leaving the heater surface at the CHF and the computed bubble radius at departure, determined from the balance between the outward velocity of the interface due to evaporation and the buoyance induced velocity of the center of mass of the bubble. The product of the CHF and the corresponding residence time was determined to be a constant for all orientations at a given bulk flow velocity and liquid subcooling, and must be determined empirically for each velocity and subcooling at present. It then becomes possible to predict the CHF for the different orientations, velocities, and subcoolings. These are compared with normalized measurements of the CHF for velocities ranging from 4 cm/s to 55 cm/s, subcoolings from 2.8 to 22.2 K, over orientations angles of 360 degrees.

Merte, Herman, Jr.↗
DOE OSTI2022

Module-Level Thermal Interface Material Degradation in HALT

Abstract In this study, thermal interface material (TIM) degradation is driven through highly accelerated life test (HALT) using temperature cycling with a prescribed vibrational acceleration for two commercially available materials having thermal conductivities of 6.0 and 8.5 W/m K. HALT specimens were prepared by applying TIM through a 4-mil stencil over AlSiC baseplates in the shape of those used in Wolfspeed CAS325M12HM2 power electronics modules. Baseplates were mounted onto aluminum carrier blocks with embedded thermocouples to characterize the thermal resistance across the baseplate and TIM layer. Thermal dissipation into the top of the baseplates was provided by a custom heating block, which mimics the size and placement of the die junctions in CAS325 modules, applying power loads of 200, 300, and 400 W. After initial characterization, samples were transferred to the HALT chamber with one set of samples exposed to temperature cycling only (TCO) and the other temperature cycling and vibration (TCV). Both sample sets were cycled between temperature extremes of −40 °C and 180 °C with vibrations applied at a peak acceleration of 3.21 Grms. After hundreds of cycles, samples were reevaluated to assess changes in thermal resistance to provide an accelerated measure of TIM degradation. This allows for reliability prediction of useful lifetime (illustrated in a solar inverter case study herein), as well as to provide a basis for developing an accelerated testing method to related temperature cycling to faster methods of degradation. Such techniques provide a means to develop maintenance schedules for power modules for ensuring sufficient thermal performance over the operating lifetime.

Engineering↗
DOE OSTI2021

Geothermal Play-Fairway Analysis of Washington State Prospects: Final Report

The Washington State Geothermal Play-Fairway Analysis overcomes the exploration challenges posed by dense vegetation, glacial deposits, and extreme precipitation. The geothermal play-fairways we target are locations where heat, permeability, and saturated porosity are present in sufficient volume to provide adequate heat exchange at depths accessible by modern drilling technology. The three study areas lie along the Cascade Range magmatic arc and are near Mount Baker, Mount St. Helens, and the Wind River Valley. The seven-year project is divided into three phases. In Phase 1 we build on a previous statewide assessment of geothermal resources and develop an initial modeling approach. The results are a series of favorability, uncertainty, and risk maps for three targeted study areas. Based on these initial results, we collect new geologic and geophysical data to further refine our modeling and reduce exploration uncertainty in Phase 2. We improve the modeling method to handle the new data and update the favorability, uncertainty, and risk maps. We also update the conceptual geothermal resource models. In Phase 3 we validate our modeling approach by drilling two temperature-gradient holes and collecting and analyzing core, image logs, and new geochemistry. Our modeling approach improves on an earlier statewide method through a more-rigorous and detailed assessment of heat and permeability. Permeability potential is assessed through geomechanical modeling of the deformation that can generate and maintain reservoir porosity and permeability. Metrics to inform heat potential include temperature-gradient wells, which are sparse in Washington; proximity of Quaternary volcanic vents and young intrusive rock; spring temperature; and reservoir temperature inferred from geothermometry. We weight the individual components using an expert-guided approach known as the Analytical Hierarchy Process. During Phase 2 we also develop a fluid-filled fracture model, and an infrastructure model that helps to delineate areas which are more favorable for geothermal development based on proximity to transmission lines, elevation, land ownership and use restrictions, and availability of process water. New geologic and geophysical data is collected during Phase 2 in each of our three main study areas. At Mount Baker and north of Mount St. Helens we conduct 1:24,000-scale geologic mapping and lidar analysis to better constrain the location and character of surface faults; detailed mapping in the Wind River Valley was completed just prior to the start of this project. Ages of intrusive rocks are determined with 40 Ar/ 39 Ar geochronology, though all of our samples are Miocene or older. We collect ground based gravity observations (a total of 1,580 new stations) in all of our study areas and ground-based magnetic lines (a total of 93 km) at Mount Baker. These data are combined with existing gravity and aeromagnetic data and used to constrain fault locations and geometry. Two to three cross sections are constructed at each study area using the mapped surface geology and forward-modeling of the gravity and magnetic data; these cross sections form the basis for our updated conceptual models. We collect magnetotelluric surveys at Mount Baker and Mount St. Helens and these data are inverted to form a resistivity model from the surface to about 10 km depth; each model shows conductive zones that can be interpreted as upwelling geothermal fluids. At Mount St. Helens we deploy a passive seismic array and use the newly detected events to refine the location of the Saint Helens seismic zone. We also employ ambient-noise tomography to develop a detailed seismic-velocity model for the study area and use this model to help constrain our cross sections and conceptual model. Based on the new data collected during Phase 2—and our updated models—we develop a campaign of temperature-gradient holes and core analysis to validate our modeling in Phase 3. Drill hole MB76-31 is located near Little Park Creek, 11 km west-southwest of the summit of Mount Baker, and is 1,471 ft deep. About 410 ft of core from the lower portion of the hole—and image logs from ~175 ft below ground surface to the bottom—are collected and analyzed. Water samples are collected and processed for geothermometry. Drill hole MSH17-24 is located along upper Schultz Creek, 16 km north-northeast of Mount St. Helens and has core from 470 ft to the bottom at 1,053 ft. We did not collect image logs due to borehole stability concerns, but water samples are collected and analyzed for geothermometry. Repeat temperature-gradient measurements are made at both sites and thermal conductivity is measured from core samples. At MB76-31, the equilibrated temperature gradient of 64°C/km and calculated heat flow of 141–159 mW/m 2 is more than twice the regional average. Detailed mapping and analysis of the core, coupled with correlation to the image logs, indicates a history of permeability generation consistent with our predictions of high permeability. Because the site has high favorability in the Phase 2 model, we consider the results a positive validation of the modeling. At site MSH17-24, the equilibrated temperature gradient of ~15°C/km and calculated heat flow of 41–43 mW/m 2 are similar to regional. Geochemical analysis of the water samples indicates a meteoric source without any geothermal component. Detailed outcrop-based mapping of fault exposures near the drill site and analysis of image logs from nearby boreholes indicates a history of permeability generation consistent with our predictions. Because the site has low favorability in the Phase 2 model, we consider the results a positive validation of the modeling. Together, the two sites provide a reasonably positive validation of the Phase 2 modeling and should encourage future use of this modeling approach.

15 GEOTHERMAL ENERGY↗
NASA NTRS2000

Estimates of Ground Temperature and Atmospheric Moisture from CERES Observations

A method is developed to retrieve surface ground temperature (Tg) and atmospheric moisture using clear sky fluxes (CSF) from CERES-TRMM observations. In general, the clear sky outgoing long-wave radiation (CLR) is sensitive to upper level moisture (q(sub h)) over wet regions and Tg over dry regions The clear sky window flux from 800 to 1200 /cm (RadWn) is sensitive to low level moisture (q(sub j)) and Tg. Combining these two measurements (CLR and RadWn), Tg and q(sub h) can be estimated over land, while q(sub h) and q(sub t) can be estimated over the oceans. The approach capitalizes on the availability of satellite estimates of CLR and RadWn and other auxiliary satellite data. The basic methodology employs off-line forward radiative transfer calculations to generate synthetic CSF data from two different global 4-dimensional data assimilation products. Simple linear regression is used to relate discrepancies in CSF to discrepancies in Tg, q(sub h) and q(sub t). The slopes of the regression lines define sensitivity parameters that can be exploited to help interpret mismatches between satellite observations and model-based estimates of CSF. For illustration, we analyze the discrepancies in the CSF between an early implementation of the Goddard Earth Observing System Data Assimilation System (GEOS-DAS) and a recent operational version of the European Center for Medium-Range Weather Prediction data assimilation system. In particular, our analysis of synthetic total and window region SCF differences (computed from two different assimilated data sets) shows that simple linear regression employing (Delta)Tg and broad layer (Delta)q(sub l) from 500 hPa to surface and (Delta)q(sub h) from 200 to 500 hPa provides a good approximation to the full radiative transfer calculations, typically explaining more than 90% of the 6-hourly variance in the flux differences. These simple regression relations can be inverted to "retrieve" the errors in the geophysical parameters. Uncertainties (normalized by standard deviation) in the monthly mean retrieved parameters range from 7% for (Delta)T to about 20% for (Delta)q(sub t). Our initial application of the methodology employed an early CERES-TRMM data set (CLR and Radwn) to assess the quality of the GEOS2 data. The results showed that over the tropical and subtropical oceans GEOS2 is, in general, too wet in the upper troposphere (mean bias of 0.99 mm) and too dry in the lower troposphere (mean bias of -4.7 mm). We note that these errors, as well as a cold bias in the Tg, have largely been corrected in the current version of GEOS-2 with the introduction of a land surface model, a moist turbulence scheme and the assimilation of SSTM/I total precipitable water.

Wu, Man Li C.↗
NASA NTRS2010

Launch Pad Flame Trench Refractory Materials

The launch complexes at NASA's John F. Kennedy Space Center (KSC) are critical support facilities for the successful launch of space-based vehicles. These facilities include a flame trench that bisects the pad at ground level. This trench includes a flame deflector system that consists of an inverted, V-shaped steel structure covered with a high temperature concrete material five inches thick that extends across the center of the flame trench. One side of the "V11 receives and deflects the flames from the orbiter main engines; the opposite side deflects the flames from the solid rocket boosters. There are also two movable deflectors at the top of the trench to provide additional protection to shuttle hardware from the solid rocket booster flames. These facilities are over 40 years old and are experiencing constant deterioration from launch heat/blast effects and environmental exposure. The refractory material currently used in launch pad flame deflectors has become susceptible to failure, resulting in large sections of the material breaking away from the steel base structure and creating high-speed projectiles during launch. These projectiles jeopardize the safety of the launch complex, crew, and vehicle. Post launch inspections have revealed that the number and frequency of repairs, as well as the area and size of the damage, is increasing with the number of launches. The Space Shuttle Program has accepted the extensive ground processing costs for post launch repair of damaged areas and investigations of future launch related failures for the remainder of the program. There currently are no long term solutions available for Constellation Program ground operations to address the poor performance and subsequent failures of the refractory materials. Over the last three years, significant liberation of refractory material in the flame trench and fire bricks along the adjacent trench walls following Space Shuttle launches have resulted in extensive investigations of failure mechanisms, load response, ejected material impact evaluation, and repair design analysis (environmental and structural assessment, induced environment from solid rocket booster plume, loads summary, and repair integrity), assessment of risk posture for flame trench debris, and justification of flight readiness rationale. Although the configuration of the launch pad, water and exhaust direction, and location of the Mobile Launcher Platform between the flame trench and the flight hardware should protect the Space Vehicle from debris exposure, loss of material could cause damage to a major element of the ground facility (resulting in temporary usage loss); and damage to other facility elements is possible. These are all significant risks that will impact ground operations for Constellation and development of new refractory material systems is necessary to reduce the likelihood of the foreign object debris hazard during launch. KSC is developing an alternate refractory material for the launch pad flame trench protection system, including flame deflector and flame trench walls, that will withstand launch conditions without the need for repair after every launch, as is currently the case. This paper will present a summary of the results from industry surveys, trade studies, life cycle cost analysis, and preliminary testing that have been performed to support and validate the development, testing, and qualification of new refractory materials.

Calle, Luz M.↗
NASA NTRS2000

Estimates of Ground Temperature and Atmospheric Moisture from CERES Observations

A method is developed to retrieve surface ground temperature (T(sub g)) and atmospheric moisture using clear sky fluxes (CSF) from CERES-TRMM observations. In general, the clear sky outgoing longwave radiation (CLR) is sensitive to upper level moisture (q(sub l)) over wet regions and (T(sub g)) over dry regions The clear sky window flux from 800 to 1200/cm (RadWn) is sensitive to low level moisture (q(sub t)) and T(sub g). Combining these two measurements (CLR and RadWn), Tg and q(sub h) can be estimated over land, while q(sub h) and q(sub l) can be estimated over the oceans. The approach capitalizes on the availability of satellite estimates of CLR and RadWn and other auxiliary satellite data. The basic methodology employs off-line forward radiative transfer calculations to generate synthetic CSF data from two different global 4-dimensional data assimilation products. Simple linear regression is used to relate discrepancies in CSF to discrepancies in T(sub g), q(sub h) and q(sub l). The slopes of the regression lines define sensitivity parameters that can be exploited to help interpret mismatches between satellite observations and model-based estimates of CSF. For illustration, we analyze the discrepancies in the CSF between an early implementation of the Goddard Earth Observing System Data Assimilation System (GEOS-DAS) and a recent operational version of the European Center for Medium-Range Weather Prediction data assimilation system. In particular, our analysis of synthetic total and window region SCF differences (computed from two different assimilated data sets) shows that simple linear regression employing Delta(T(sub g)) and broad layer Delta(q(sub l) from .500 hPa to surface and Delta(q(sub h)) from 200 to .300 hPa provides a good approximation to the full radiative transfer calculations. typically explaining more than 90% of the 6-hourly variance in the flux differences. These simple regression relations can be inverted to "retrieve" the errors in the geophysical parameters. Uncertainties (normalized by standard deviation) in the monthly mean retrieved parameters range from 7% for Delta(T(sub g)) to about 20% for Delta(q(sub l)). Our initial application of the methodology employed an early CERES-TRMM data set (CLR and Radwn) to assess the quality of the GEOS2 data. The results showed that over the tropical and subtropical oceans GEOS2 is, in general, too wet in the upper troposphere (mean bias of 0.99 mm) and too dry in the lower troposphere (mean bias of -4.7 min). We note that these errors, as well as a cold bias in the T(sub g). have largely been corrected in the current version of GEOS-2 with the introduction of a land surface model, a moist turbulence scheme and the assimilation of SSM/I total precipitable water.

Wu, Man Li C.↗