Search NASA⌕ Search

SEARCH · Search NASA

Results for “determining the coordination number”

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 253 records · Page 14

Kepler Mission's Focal Plane Characterization Models Implementation

The Kepler Mission photometer is an unusually complex array of CCDs. A large number of time-varying instrumental and systemic effects must be modeled and removed from the Kepler pixel data to produce light curves of sufficiently high quality for the mission to be successful in its planet-finding objective. After the launch of the spacecraft, many of these effects are difficult to remeasure frequently, and various interpolations over a small number of sample measurements must be used to determine the correct value of a given effect at different points in time. A library of software modules, called Focal Plane Characterization (FC) Models, is the element of the Kepler Science Data Pipeline (hereafter "pipeline") that handles this. FC, or products generated by FC, are used by nearly every element of the SOC processing chain. FC includes Java components: database persistence classes, operations classes, model classes, and data importers; and MATLAB code: model classes, interpolation methods, and wrapper functions. These classes, their interactions, and the database tables they represent, are discussed. This paper describes how these data and the FC software work together to provide the pipeline with the correct values to remove non-photometric effects caused by the photometer and its electronics from the Kepler light curves. The interpolation mathematics is reviewed, as well as the special case of the sky-to-pixel,pixel-to-sky coordinate transformation code, which incorporates a compound model that is unique in the SOC software.

mission↗

Metagenome-assembled genomes measured at 3 depths during snowmelt period in East River, CO (March, May, and June, September 2017)

Snowmelt is a critical biogeochemical period that accounts for large nitrogen (N) export events from high-elevation watersheds. Soil microbial populations bloom and immobilize N during snowmelt, yet the population size crashes in spring, which releases a pulse of soil N. We sought to discover the N sources fueling this microbial bloom and determine the fate of N following microbial die-off. Here, focusing on the snowmelt period within a headwater catchment of the Upper Colorado River Basin (East River, CO), we deployed strain-resolved metagenomics to identify the metabolic pathways and processes that mobilize soil N during and after snowmelt. Soil metagenome samples were taken from 6 snowpits from 3 depths (0-5cm, 5-15cm, >15cm) at 4 time points during snowmelt period (March 2017, May 2017, and June 2017, September 2017) generating 48 metagenomes. We reconstructed 474 metagenome-assembled genomes (MAGs) across all metagenomes.All 48 metagenomes were sequenced at JGI and raw data can be found under JGI (Joint Genome Institute) GOLD Study Gs0135149. Metagenome assemblies from IMG under the same study were used for genome binning. This dataset (1) a zip file of 474 MAGs (as fasta files, Gs0135149_bins_tar.gz), (2) sample metadata file with sample IGSNs (International Generic Sample Numbers) (samples.csv), (3) bounding box coordinates for the sampled locations (Gs0135149.kml), (4) metagenome metadata file listing IMG/M (Integrated Microbial Genomes/Metagenomes) metagenome accessions linking samples to metagenomes (metagenomes.csv), (5) location metadata file (locations.csv), (6) file-level metadata file (flmd.csv) and (7) data dictionary (dd.csv) file.This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

54 ENVIRONMENTAL SCIENCES↗

Numerical Modeling of a Vortex Stabilized Arcjet

Arcjet thrusters are being actively considered for use in Earth orbit maneuvering applications. Experimental studies are currently the chief means of determining an optimal thruster configuration. Earlier numerical studies have failed to include all of the effects found in typical arcjets including complex geometries, viscosity, and swirling flow. Arcjet geometries are large area ratio converging nozzles with centerbodies in the subsonic portion of the nozzle. The nozzle walls serve as the anode while the centerbody functions as the cathode. Viscous effects are important because the Reynolds number, based on the throat radius, is typically less than 1,000. Experimental studies have shown that a swirl or circumferential velocity component stabilizes a constricted arc. This dissertation describes the equations governing flow through a constricted arcjet thruster. An assumption that the flowfield is in local thermodynamic equilibrium leads to a single fluid plasma temperature model. An order of magnitude analysis reveals the governing fluid mechanics equations are uncoupled from the electromagnetic field equations. A numerical method is developed to solve the governing fluid mechanics equations, the Thin Layer Navier-Stokes equations. A coordinate transformation is employed in deriving the governing equations to simplify the application of boundary conditions in complex geometries. An axisymmetric formulation is employed to include the swirl velocity component as well as the axial and radial velocity components. The numerical method is an implicit finite-volume technique and allows for large time steps to reach a converged steady-state solution. The inviscid fluxes are flux-split, and Gauss-Seidel line relaxation is used to accelerate convergence. Converging-diverging nozzles with exit-to-throat area ratios up to 100:1 and annular nozzles were examined. Quantities examined included Mach number and static wall pressure distributions, and oblique shock structures. As the level of swirl and viscosity in the flowfield increased the mass flow rate and thrust decreased. The technique was used to predict the flow through a typical arcjet thruster geometry. Results indicate swirl and viscosity play an important role in the complex geometry of an arcjet.

Pawlas, Gary E.↗

Unravelling the mechanism of pressure induced polyamorphic transition in an inorganic molecular glass

The atomic structure of a germanium doped phosphorous selenide glass of composition Ge 2.8 P 57.7 Se 39.5 is determined as a function of pressure from ambient to 24 GPa using Monte-Carlo simulations constrained by high energy x-ray scattering data. The ambient pressure structure consists primarily of P 4 Se 3 molecules and planar edge shared phosphorus rings, reminiscent of those found in red phosphorous as well as a small fraction of locally clustered corner-sharing GeSe 4 tetrahedra. This low-density amorphous phase transforms into a high-density amorphous phase at ~6.3 GPa. The high-pressure phase is characterized by an extended network structure. The polyamorphic transformation between these two phases involves opening of the P 3 ring at the base of the P 4 Se 3 molecules and subsequent reaction with red phosphorus type moieties to produce a cross linked structure. The compression mechanism of the low-density phase involves increased molecular packing, whereas that of the high pressure phase involves an increase in the nearest-neighbor coordination number while the bond angle distributions broaden and shift to smaller angles. The entropy and volume changes associated with this polyamorphic transformation are positive and negative, respectively, and consequently the corresponding Clapeyron slope for this transition would be negative. This result has far reaching implications in our current understanding of the thermodynamics of polyamorphic transitions in glasses and glass-forming liquids.

42 ENGINEERING↗

In Situ Determination of Speciation and Local Structure of NaCl–SrCl 2 and LiF–ZrF 4 Molten Salts

Understanding the local environment of the metal atoms in salt melts is important for modeling the properties of melts and predicting their behavior and thus helping enable the development of technologies such as molten salt reactors and solar-thermal power systems and new approaches to recycling rare-earth metals. Toward that end, we have developed an in situ approach for measuring the coordination of metals in molten salt coupling X-ray absorption spectroscopy (XAS) and Raman spectroscopy. Our approach was demonstrated for two salt mixtures (1.9 and 5 mol % SrCl 2 in NaCl, 0.8 and 5 mol % ZrF 4 in LiF) at up to 1100 °C. Near-edge (X-ray absorption near-edge structure, XANES) and extended X-ray absorption fine structure (EXAFS) spectra were measured. The EXAFS response was modeled using ab initio FEFF calculations. Strontium’s first shell is observed to be coordinated with chlorine (Sr 2+ –Cl – ) and zirconium’s first shell is coordinated by fluorine (Zr 4+ –F – ), both having coordination numbers that decrease with increasing temperature. Multiple zirconium complexes are believed to be present in the melt, which may interfere and distort the EXAFS spectra and result in an anomalously low zirconium first shell coordination number. Finally, the use of boron nitride (BN) powder as a salt diluent for XAFS measurements was found to not interfere with measurements and thus can be used for investigations of such systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impingement of Droplets in 60 Deg Elbows with Potential Flow

Trajectories were determined for water droplets or other aerosol particles in air flowing through 600 elbows especially designed for two-dimensional potential motion. The elbows were established by selecting as walls of each elbow two streamlines of a flow field produced by a complex potential function that establishes a two-dimensional flow around. a 600 bend. An unlimited number of elbows with slightly different shapes can be established by selecting different pairs of streamlines as walls. Some of these have a pocket on the outside wall. The elbows produced by the complex potential function are suitable for use in aircraft air-inlet ducts and have the following characteristics: (1) The resultant velocity at any point inside the elbow is always greater than zero but never exceeds the velocity at the entrance. (2) The air flow field at the entrance and exit is almost uniform and rectilinear. (3) The elbows are symmetrical with respect to the bisector of the angle of bend. These elbows should have lower pressure losses than bends of constant cross-sectional area. The droplet impingement data derived from the trajectories are presented along with equations so that collection efficiency, area, rate, and distribution of droplet impingement can be determined for any elbow defined by any pair of streamlines within a portion of the flow field established by the complex potential function. Coordinates for some typical streamlines of the flow field and velocity components for several points along these streamlines are presented in tabular form. A comparison of the 600 elbow with previous calculations for a comparable 90 elbow indicated that the impingement characteristics of the two elbows were very similar.

Hacker, Paul T.↗

Reining in Radium for Nuclear Medicine: Extra-Large Chelator Development for an Extra-Large Ion

Targeted α therapy (TAT) of soft-tissue cancers using the α particle-emitting radionuclide 223 Ra holds great potential because of its favorable nuclear properties, adequate availability, and established clinical use for treating metastatic prostate cancer of the bone. Despite these advantages, the use of 223 Ra has been largely overshadowed by other α emitters due to its challenging chelation chemistry. A key criterion that needs to be met for a radionuclide to be used in TAT is its stable attachment to a targeting vector via a bifunctional chelator. The low charge density of Ra 2+ arising from its large ionic radius weakens its electrostatic binding interactions with chelators, leading to insufficient complex stability in vivo. In this study, we synthesized and evaluated macropa-XL as a novel chelator for 223 Ra. It bears a large 21-crown-7 macrocyclic core and two picolinate pendent groups, which we hypothesized would effectively saturate the large coordination sphere of the Ra 2+ ion. The structural chemistry of macropa-XL was first established with the nonradioactive Ba 2+ ion using X-ray diffraction and X-ray absorption spectroscopy, which revealed the formation of an 11-coordinate complex in a rare anti pendent-arm configuration. Subsequently, the stability constant of the [Ra(macropa-XL)] complex was determined via competitive cation exchange with 223 Ra and 224 Ra radiotracers and compared with that of macropa, the current state-of-the-art chelator for Ra 2+ . A moderate log K ML value of 8.12 was measured for [Ra(macropa-XL)], which is approximately 1.5 log K units lower than the stability constant of [Ra(macropa)]. This relative decrease in Ra 2+ complex stability for macropa-XL versus macropa was further probed using density functional theory calculations. Additionally, macropa-XL was radiolabeled with 223 Ra, and the kinetic stability of the resulting complex was evaluated in human serum. Furthermore, although macropa-XL could effectively bind 223 Ra under mild conditions, the complex appeared to be unstable to transchelation. Collectively, this study sheds additional light on the chelation chemistry of the exotic Ra 2+ ion and contributes to the small, but growing, number of chelator development efforts for 223 Ra-based TAT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Human Research Program Human Health Countermeasures Element Cardiovascular Risks Standing Review Panel (SRP)

The Cardiovascular Risk Standing Review Panel (SRP) evaluated several cardiovascular risks associated with space flight along with the ongoing and emerging plans to study these issues and potentially propose and/or develop countermeasures. The areas of focus included: 1) The risk of cardiac rhythm problems during prolonged space flight, and 2) Issues related to the risk of orthostatic intolerance during re-exposure to gravity. An emerging area of concern is radiation associated vascular injury. The risk of cardiac rhythm disturbances has emerged based on case reports only. No systematic study of this risk has been published. However, concerns about this risk are heightened by the age range of astronauts, the structural changes in the heart that occur during space flight, and the potential shifts in fluids and electrolytes. The current plan is to use prolonged Holter monitor EKG records made as part of the "Integrated Cardiovascular SMO" in space to determine more about the frequency and magnitude of this problem and to link this data to complementary data from the nutrition group on electrolytes. The SRP was supportive of this approach. The SRP also felt that any data related to cardiovascular risk in space should be better coordinated with the medical screening data that all astronauts undergo at regular intervals. Additionally, while there are potential privacy issues related to this suggestion, many of the current barriers to better coordination of experimental and clinical data appear to reflect longstanding cultural traditions at NASA that need rethinking. The risk of orthostatic intolerance during re-exposure to gravity was seen by the SRP as an area supported by a wealth of published physiological evidence. The SRP also felt that moving forward with the planned approach to countermeasures was reasonable and that extensive additional hypothesis testing on the physiology of orthostatic intolerance was not needed at this time. There was support for developing ground based models of limited (e.g. 1/6 th) G environments on Earth that generated a number of ideas for consideration by NASA investigators.

Joyner, Michael↗

Differential Processing of Low and High LET Radiation Induced DNA Damage: Investigation of Switch from ATM to ATR Signaling

The members of the phosphatidylinositol kinase-like kinase family of proteins namely ataxia-telangiectasia mutated (ATM) and ATM- and Rad3-related (ATR) are directly responsible for the maintenance of genomic integrity by mounting DDR through signaling and facilitating the recruitment of repair factors at the sites of DNA damage along with coordinating the deployment of cell cycle checkpoints to permit repair by phosphorylating Checkpoint kinase Chk1, Chk2 and p53. High LET radiation from GCR (Galactic Cosmic Rays) consisting mainly of protons and high energy and charged (HZE) particles from SPE (Solar Particle Event) pose a major health risk for astronauts on their space flight missions. The determination of these risks and the design of potential safeguards require sound knowledge of the biological consequences of lesion induction and the capability of the cells to counter them. We here strive to determine the coordination of ATM and ATR kinases at the break sites directly affecting checkpoint signaling and DNA repair and whether differential processing of breaks induced by low and high LET radiation leads to possible augmentation of swap of these damage sensors at the sites of DNA damage. Exposure of cells to IR triggers rapid autophosphorylation of serine-1981 that causes dimer dissociation and initiates monomer formation of ATM. ATM kinase activity depends on the disruption of the dimer, which allows access and phosphorylation of downstream ATM substrates like Chk2. Evidence suggests that ATM is activated by the alterations in higher-order chromatin structure although direct binding of ATM to DSB ends may be a crucial step in its activation. On the other hand, in case of ATR, RPA (replication protein A)-coated ssDNA (single-stranded DNA) generated as a result of stalled DNA replication or during processing of chromosomal lesions is crucial for the localization of ATR to sites of DNA damage in association with ATR-interacting protein (ATRIP). Although the majority of RPA-coated ssDNA is generally present only during DNA replication, ATR activation in G1 and G2-phase might still require formation of RPA-coated ssDNA, probably initiated by the MRN-CtIP complex and then extended by the Exo1- or BLM-dependent mechanisms at the sites of DSBs. Evidence accumulates that activation of ATM and ATR are oppositely regulated by the length of single stranded overhangs generated at the break sites by processes mentioned above and these stretches of single stranded overhangs hold the clue for ATM to ATR switch at broken DNA ends. We irradiated 82-6hTERT human fibroblast cells with low LET gamma-rays and high LET Fe and Si particles. Preliminary results with cells exposed to 1Gy gamma-rays show that the kinetics of pChk2-pT68 foci formation is comparable to that of gamma-H2AX although they appear to recede quicker. The number and intensity of observed foci reaches a maximum at 30 min and 60 min post IR for Chk2-pT68 and gamma-H2AX foci respectively and all Chk2-pT68 foci colocalize with gamma-H2AX foci. The kinetics of Chk1-pS345 and ATRIP are being determined. Results of Chk2-pT68 foci kinetics was also corroborated by western blot experiments, although phosphorylation was detected as early as 10 min and started receding 30 min post IR with 2Gy of gamma-rays. On the other hand, level of ATR-pS428 reached its maximum between 60 and 120 min and was maintained until the last measured time point of 4 hours post IR as determined by western blotting. Experiments performed with high LET Fe and Si particles will be reported.

Saha, Janapriya↗

Reflector Surface Error Compensation in Dual-Reflector Antennas

By probing the field on a small subreflector at a minimal number of points, the main reflector surface errors can be obtained and subsequently used to design a phase-correction subreflector that can compensate for main reflector errors. The compensating phase-error profile across the subreflector can be achieved either by a surface deformation or by the use of an array of elements such as patch antennas that can cause a phase shift between the incoming and outgoing fields. The second option is of primary interest here, but the methodology can be applied to either case. The patch array is most easily implemented on a planar surface. Therefore, the example of a flat subreflector and a parabolic main reflector (a Newtonian dual reflector system) is considered in this work. The subreflector is assumed to be a reflector array covered with patch elements. The phase variation on a subreflector can be detected by a small number of receiving patch elements (probes). By probing the phase change at these few selected positions on the subreflector, the phase error over the entire surface can be recovered and used to change the phase of all the patch elements covering the subreflector plane to compensate for main reflector errors. This is accomplished by using a version of sampling theorem on the circular aperture. The sampling is performed on the phase-error function on the circular aperture of the main reflector by a method developed using Zernike polynomials. This method is based upon and extended from a theory previously proposed and applied to reflector aperture integration. This sampling method provides for an exact retrieval of the coefficients of up to certain orders in the expansion of the phase function, from values on a specifically calculated set of points in radial and azimuthal directions in the polar coordinate system, on the circular reflector aperture. The corresponding points on the subreflector are then obtained and, by probing the fields at these points, a set of phase values is determined that is then transferred back to the main reflector aperture for recovering the phase function. Once this function is recovered, the corresponding phase function on the subreflector is calculated and used to compensate for main reflector surface errors. In going back and forth between sub and main reflectors, geometrical (ray) optics is employed, which even though it ignores edge diffraction and other effects, is shown to be accurate for phase recovery.

Jamnejad, Vahraz↗

Stability and Control Characteristics at Subsonic Speeds of a Flat-Top Arrowhead Wing-Body Combination

A wind-tunnel investigation was made to determine the longitudinal- and lateral-stability derivatives of a flat-top wing-body configuration at Mach numbers from 0.22 to 0.90 and Reynolds numbers of 3.5 and 17 million. The wing had a leading-edge sweepback of 78.9 deg and a cathedral of 45 deg on the outer panels. The tests included the determination of the effectiveness of elevon and rudder controls and also an investigation of ground effects. The model was tested at angles of attack up to 28 deg and angles of sideslip up to 18 deg. The dynamic response of this configuration has been determined from the wind-tunnel data for a simulated airplane having a wing loading of 17.7 pounds per square foot. The longitudinal data show a forward shift in aerodynamic center of 10 percent of the mean aerodynamic chord as the lift coefficient is increased above 0.1. Although flown in the lift range of decreasing stability, the simulated airplane did not encounter pitch-up in maneuvers initiated from steady level flight with zero static margin unless a load factor of 2.2 was exceeded. This maneuver margin was provided by a large value of pitching moment due to pitching velocity. The number of cycles to damp the Dutch roll mode to half amplitude, the time constants of the roll subsidence and spiral divergence modes, and control effectiveness in roll are computed. The lateral stability is shown to be positive but is marginal in meeting the military specifications for today's aircraft. An analog computer study has been made in five degrees of freedom (constant velocity) which illustrates that the handling characteristics are satisfactory. Several programed rolling maneuvers and coordinated turns also illustrate the handling qualities of the airplane.

Buell, Donald A.↗

GT Flex: A Coordinated Multi-Building Pilot Study

Buildings are a significant and untapped resource for providing utility electric grid services. Recent studies have estimated that buildings could reduce the peak demand on the electric grid in the U.S. by almost 25% through effective combinations of energy efficiency (EE) measures and load flexibility strategies (Langevin et al. 2021). The U.S. Department of Energy (DOE) has established a goal to triple energy efficiency and demand flexibility in both residential and commercial sections by 2030 compared to 2020 levels (Satchwell et al. 2021). Such findings place buildings alongside electric vehicles, photovoltaics, electric batteries, and other distributed energy resources (DERs) as primary technologies needed for supporting high renewable energy generation grids. Coordinating and optimizing multiple buildings and other DERs is more beneficial and valuable when compared with individual buildings and DERs operating as siloed resources, uncoordinated with others (Olgyay et al. 2020). A pilot study at the Georgia Institute of Technology (GIT) was conducted to evaluate value propositions of a multi-building scale project seeking carbon reduction, energy efficiency and grid-interactive capabilities, by demonstrating the means by which stakeholders can determine the technical and financial merits of grid-interactivity and energy efficiency technologies coordinated across multiple assets. The study focused on analyzing technical feasibility of deploying thermal load flexibility strategies at the multi-building scale, coordinated to not exceed existing infrastructure constraints at the pilot site. Results show that campus can provide 3-3.5 MW of potential load shed over a 4-hour event window through coordinated dispatch of thermal cooling load flexibility without exceeding existing infrastructure capacities. Under future high renewable scenarios, this thermal flexibility resource is also valuable when coordinated to reduce curtailment of intermittent renewables. Economic analyses were performed to effectively communicate various value propositions of grid-interactive and efficient building (GEB) thermal flexibility strategies. Load flexibility presents a financial value proposition to the campus today. By conducting rationalized, coordinated dispatch in response to real time price (RTP) fluctuations, the campus can benefit materially from daily price arbitrage. The RTP signal acts as an aggregating mechanism between the utility and customer to call on demand flexibility resources, with a large portion of the benefit deriving from a relatively small number of days. Realizing and maximizing this benefit with thermal load flexibility requires careful attention to the timing of pricing signals and parameterization of dispatch to overcome efficiency penalties. Grid value and signals are expected to evolve over time, and thermal load flexibility shows potential to adapt dispatch logic to support intermittent renewable generation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

The Photogrammetry Cube

We can determine distances between objects and points of interest in 3-D space to a useful degree of accuracy from a set of camera images by using multiple camera views and reference targets in the camera s field of view (FOV). The core of the software processing is based on the previously developed foreign-object debris vision trajectory software (see KSC Research and Technology 2004 Annual Report, pp. 2 5). The current version of this photogrammetry software includes the ability to calculate distances between any specified point pairs, the ability to process any number of reference targets and any number of camera images, user-friendly editing features, including zoom in/out, translate, and load/unload, routines to help mark reference points with a Find function, while comparing them with the reference point database file, and a comprehensive output report in HTML format. In this system, scene reference targets are replaced by a photogrammetry cube whose exterior surface contains multiple predetermined precision 2-D targets. Precise measurement of the cube s 2-D targets during the fabrication phase eliminates the need for measuring 3-D coordinates of reference target positions in the camera's FOV, using for example a survey theodolite or a Faroarm. Placing the 2-D targets on the cube s surface required the development of precise machining methods. In response, 2-D targets were embedded into the surface of the cube and then painted black for high contrast. A 12-inch collapsible cube was developed for room-size scenes. A 3-inch, solid, stainless-steel photogrammetry cube was also fabricated for photogrammetry analysis of small objects.

Source record↗

Hawai‘i Play Fairway (Final Report)

Most of Hawai'i's geothermal resources are blind—their manifestations, such as hot springs and steam vents, do not appear on the ground surface because the heated water flows far below. With the exception of $K\bar{i}lauea East Rift Zone$, in most areas of Hawai'i, high lateral permeability in the first kilometer below ground surface prevents surface thermal features from developing. As a methodology for discovering these blind resources, Play Fairway Analysis (PFA) involves finding potential locations of blind hydrothermal systems and describing potential geothermal sources in rift-zone settings. Using the PFA to find Hawai'i's geothermal resources, the University of Hawai'i (UH) conducted the Hawai‘i Play Fairway Project, Hawai'i's first statewide geothermal resource assessment since 1985. Sponsored by the U.S. Department of Energy, the Hawai'i Play Fairway Project provided an updated resource assessment, a roadmap for additional exploration activities, and the identification of areas for further exploration. Benefitting from UH's core competency in earth sciences and experienced geothermal researchers, the project comprised three phases. During the first phase, the team identified, compiled, and ranked existing geologic, groundwater, and geophysical datasets relevant to subsurface heat, fluid and permeability. Using a Bayesian statistical approach, the team developed a statistical methodology to integrate these data into a resource probability map. The team evaluated the confidence in the probability value and considered development viability of areas with geothermal resources. With these analyses, the team identified 10 locations in the Hawaiian Islands for exploration activities. For the second phase, the team collected new groundwater data in 10 locations across the state and new geophysical data on $L\bar{a}na‘i, Maui$, and central Hawai'i Island and modeled topographically induced stress to better characterize subsurface permeability. Analyzing the subsurface stresses, the team evaluated the potential for fracture-induced permeability. The team inverted the MT and gravity data to produce 3D models of resistivity and density, respectively, on $L\bar{a}na‘i$, across $Haleakal\bar{a}'s$ SW rift (Maui), and surrounding Mauna Kea (Hawai‘i Island). The team developed and applied a new method for incorporating depth information about resistivity, density, and potential for fracture-induced permeability into the statistical method for computing resource probability in these three focus areas. The team incorporated the new groundwater results with the new geophysical results and the calculations of potential for fracture-induced permeability to produce updated maps of resource probability and confidence. Through combining data from the first and second phases, the team determined locations for further exploration during the third phase. For MT and gravity surveys, the team recommended $Kaua'i's$ $L\bar{i}hu'e$ $Basin$, the east rift of $Maui's$ $Haleakal\bar{a}$ volcano, and the southwest rift of Hawai'i Island's Mauna Loa volcano. The MT and gravity surveys aimed to enable improved confidence in the resource potential in these locations. For drilling deep groundwater well(s), the team recommended Southeast Mauna Kea and $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$. During the third phase, further exploration involved drilling a groundwater well in $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$ and performing more geophysical surveys. We deepened an existing water well proximal to our target area on $L\bar{a}na'i$ due to funding constraints that precluded us from spudding a new well that would exceed 1km depth. Drilling was preceded by a number of substantial elements including: writing an Environmental Assessment and the subsequent legal process, performance of deviation logging, lowering a camera down the well, coordinating site preparation with $P\bar{u}lama$ $L\bar{a}na'i$, shipping the UH-owned rig interisland, procuring supplies, and leading 3 community meetings on $L\bar{a}na'i$. Drilling occurred 24/7 the entire month of June 2019 over which time $L\bar{a}na'i$ $Well$ 10 was deepened from 427 m to 1057 m, with continuous core collected. We measured a roughly linear temperature gradient averaging 42°C/km and a maximum bottom hole temperature of 66°C. This gradient is more than twice the background for Hawai'i and within a range of gradients measured in this depth range for some exploration wells within KERZ. We consider these results encouraging for $L\bar{a}na'i's$ resource potential and recommend following with a slim hole within $L\bar{a}na'i's$ caldera (our target zone) to ~ 2 km. Further, the positive implications such results have for the island of O‘ahu are substantial - the shield stage of O'ahu's volcanoes ended 1-2 My earlier. However, O'ahu uses more electricity than the rest of the islands combined, and the utility recently called for 500-700MW of firm, dispatchable renewable electricity on O'ahu by 2033. In Phase 3, we also collected limited new encouraging groundwater data, and updated our thoughts on the probabilities of fluid and permeability at resource depths (Pr F = 1; Pr P = mostly unconstrained). Ultimately, we advocate for using our final probability of heat, and confidence in this probability, to drive the next phase of exploration. We contend further development of geothermal in Hawai‘i will enable the state to achieve its 100% renewable policy objective and Hawai'i to transition off of fossil fuels through geothermal discovery and development. The project not only produced a large amount of data and expanded the existing knowledge of Hawai'i's geothermal resources, but also produced publications, theses, presentations, core photos, datasets, media reports, television interviews, community events, and a blog. Students and new professionals benefitted from the project's hands-on research experiences and educational opportunities and earned awards and recognition.

15 GEOTHERMAL ENERGY↗

Inventory of forest and rangeland and detection of forest stress

The author has identified the following significant results. Three small scales of CIR photography were interpreted to determine the number of bark beetle-killed trees detected in each of six spot size categories. A procedure was developed to predict the probability of detecting spots in each spot size category and in turn to estimate the number of infestations and dead trees even on the smallest scale. Statistical tests of the data indicated that the linear model did not fit the data and that other models should be tested. As a result of daily monitoring of Black Hills radiometric instruments it was possible to show the spectral energy relationships in the ponderosa pine ecosystems over time. These data have been helpful for comparison with radiance signatures extracted from ERTS-1 bulk 70mm using precision microdensitometry. Effects of atmospheric interference were shown by a 30 percent increase in scene radiance on channel 4 of the satellite imagery. A calibration and scaling technique was developed and tested to enable interpretation of ERTS-1 bulk and precision data for the Atlanta test site. The technique includes calibration of a photographic copy system for the I2S image combiner and the production of scaled overlays of grid coordinate systems, study area locations, and outline maps of county boundaries.

Heller, R. C.↗

Interplanetary Fast Shocks and Associated Drivers Observed through the Twenty-Third Solar Minimum by WIND Over its First 2.5 Years

A list of the interplanetary (IP) shocks observed by WIND from its launch (in November 1994) to May 1997 is presented. Forty two shocks were identified. The magnetohydrodynamic nature of the shocks is investigated, and the associated shock parameters and their uncertainties are accurately computed using a practical scheme which combines two techniques. These techniques are a combination of the "pre-averaged" magnetic-coplanarity, velocity-coplanarity, and the Abraham-Schrauner-mixed methods, on the one hand, and the Vinas and Scudder [1986] technique for solving the non-linear least-squares Rankine-Hugoniot shock equations, on the other. Within acceptable limits these two techniques generally gave the same results, with some exceptions. The reasons for the exceptions are discussed. It is found that the mean strength and rate of occurrence of the shocks appears to correlated with the solar cycle. Both showed a decrease in 1996 coincident with the time of the lowest ultraviolet solar radiance, indicative of solar minimum and start of solar cycle 23, which began around June 1996. Eighteen shocks appeared to be associated with corotating interaction regions (CIRs). The distribution of their shock normals showed a mean direction peaking in the ecliptic plane and with a longitude (phi(sub n)) in that plane between perpendicular to the Parker spiral and radial from the Sun. When grouped according to the sense of the direction of propagation of the shocks the mean azimuthal (longitude) angle in GSE coordinates was approximately 194 deg for the fast-forward and approximately 20 deg for the fast-reverse shocks. Another 16 shocks were determined to be driven by solar transients, including magnetic clouds. These shocks had a broader distribution of normal directions than those of the CIR cases with a mean direction close to the Sun-Earth line. Eight shocks of unknown origin had normal orientation well off the ecliptic plane. No shock propagated with longitude phi(sub n) >= 220 +/- 10 deg, this would suggest strong hindrance to the propagation of shocks contra a rather tightly winding Parker spiral. Examination of the obliquity angle theta(sub Bn) (that between the shock normal and the upstream interplanetary magnetic field) for the full set of shocks revealed that about 58% was quasi-perpendicular, and some were very nearly perpendicular. About 32% of the shocks were oblique, and the rest (only 10%) were quasi-parallel, with one on Dec. 9, 1996 that showed field pulsations. Small uncertainty in the estimated angle theta(sub Bn) was obtained for about 10 shocks with magnetosonic Mach numbers between 1 and 2, hopefully significantly contributing to studies researching particle acceleration mechanisms at IP shocks, and to investigations where accurate values of theta(sub Bn) are crucial.

Mariani, F.↗

Selection Effects in Identifying Magnetic Clouds and the Importance of the Closest Approach Parameter

This study is motivated by the unusually low number of magnetic clouds (MCs) that are strictly identified within interplanetary coronal mass ejections (ICMEs), as observed at 1 AU; this is usually estimated to be around 30% or lower. But a looser definition of MCs may significantly increase this percentage. Another motivation is the unexpected shape of the occurrence distribution of the observers' "closest approach distances" (measured from a MC's axis, and called CA) which drops off somewhat rapidly as |CA| (in % of MC radius) approaches 100%, based on earlier studies. We suggest, for various geometrical and physical reasons, that the |CA|-distribution should be somewhere between a uniform one and the one actually observed, and therefore the 30% estimate should be higher. So we ask, When there is a failure to identify a MC within an ICME, is it occasionally due to a large |CA| passage, making MC identification more difficult, i.e., is it due to an event selection effect? In attempting to answer this question we examine WIND data to obtain an accurate distribution of the number of MCs vs. |CA| distance, whether the event is ICME-related or not, where initially a large number of cases (N=98) are considered. This gives a frequence distribution that is far from uniform, confirming earlier studies. This along with the fact that there are many ICME identification-parameters that do not depend on |CA| suggest that, indeed an MC event selection effect may explain at least part of the low ratio of (No. MCs)/(No. ICMEs). We also show that there is an acceptable geometrical and physical consistency in the relationships for both average "normalized" magnetic field intensity change and field direction change vs. |CA| within a MC, suggesting that our estimates of |CA|, B(sub 0) (magnetic field intensity on the axis), and choice of a proper "cloud coordinate" system (all needed in the analysis) are acceptably accurate. Therefore the MC fitting model (Lepping et al., 1990) is adequate, on average, for our analysis. However, this selection effect is not likely to completely answer our original question, on the unexpected ratio of MCs to ICMEs, so we must look for other factors, such as peculiarities of CME birth conditions. As a by-product of this analysis, we determine that the first order structural effects within a MC due to its interaction with the solar wind, plus the MC's usual expansion at 1 AU (i.e., the non-force free components of the MC's field) are, on average, weakly dependent on radial distance from the MC's axis; that is, in the outer reaches of a typical MC the non-force free effects show up, but even there they are rather weak. Finally, we show that it is not likely that a MC's size distribution statistically controls the occurrence distribution of the estimated |CA|s.

Lepping, R. P.↗

An Explicit Upwind Algorithm for Solving the Parabolized Navier-Stokes Equations

An explicit, upwind algorithm was developed for the direct (noniterative) integration of the 3-D Parabolized Navier-Stokes (PNS) equations in a generalized coordinate system. The new algorithm uses upwind approximations of the numerical fluxes for the pressure and convection terms obtained by combining flux difference splittings (FDS) formed from the solution of an approximate Riemann (RP). The approximate RP is solved using an extension of the method developed by Roe for steady supersonic flow of an ideal gas. Roe's method is extended for use with the 3-D PNS equations expressed in generalized coordinates and to include Vigneron's technique of splitting the streamwise pressure gradient. The difficulty associated with applying Roe's scheme in the subsonic region is overcome. The second-order upwind differencing of the flux derivatives are obtained by adding FDS to either an original forward or backward differencing of the flux derivative. This approach is used to modify an explicit MacCormack differencing scheme into an upwind differencing scheme. The second order upwind flux approximations, applied with flux limiters, provide a method for numerically capturing shocks without the need for additional artificial damping terms which require adjustment by the user. In addition, a cubic equation is derived for determining Vegneron's pressure splitting coefficient using the updated streamwise flux vector. Decoding the streamwise flux vector with the updated value of Vigneron's pressure splitting improves the stability of the scheme. The new algorithm is applied to 2-D and 3-D supersonic and hypersonic laminar flow test cases. Results are presented for the experimental studies of Holden and of Tracy. In addition, a flow field solution is presented for a generic hypersonic aircraft at a Mach number of 24.5 and angle of attack of 1 degree. The computed results compare well to both experimental data and numerical results from other algorithms. Computational times required for the upwind PNS code are approximately equal to an explicit PNS MacCormack's code and existing implicit PNS solvers.

Korte, John J.↗