Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “global refinement”

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 577 records · Page 32

Preliminary result of teleseismic double-difference relocation of earthquakes in the Molucca collision zone with a 3D velocity model

We have relocated hypocenters of earthquakes occurring in the Molucca collision zone and surrounding region taken from the BMKG catalog using teleseismic double-difference relocation algorithm (teletomoDD). We used P-wave arrival times of local, regional, and teleseismic events recorded at 304 recording stations. Over 7,000 earthquakes were recorded by the BMKG seismographicnetworkin the study region from April, 2009 toJune, 2014. We used a 3D regional-global nested velocity modelresulting fromprevious global tomographystudy. In this study, the3D seismic velocity model was appliedto theIndonesian region, whilethe1D seismicvelocity model (ak135)wasused for regions outside of Indonesia. Our relocation results show a better improvement in travel-time RMS residuals comparedto those of the BMKG catalog.Ourresultsalso show that relocation shifts were dominated intheeast-west direction, whichmaybeinfluenced by theexistingvelocity anomaly related to the reversed V-shaped slabbeneaththestudy region. Our eventrelocation results refine the geometry of slabs beneath the Halmahera and Sangihe arcs.

58 GEOSCIENCES↗

Photodynamics of oxybenzone sunscreen: Nonadiabatic dynamics simulations

Herein we have used combined static electronic structure calculations and “on-the-fly” global-switching trajectory surface-hopping dynamics simulations to explore the photochemical mechanism of oxybenzone sunscreen. We have first employed the multi-configurational CASSCF method to optimize minima, conical intersections, and minimum-energy reaction paths related to excited-state intramolecular proton transfer (ESIPT) and excited-state decays in the {sup 1}ππ{sup ∗}, {sup 1}nπ{sup ∗}, and S{sub 0} states (energies are refined at the higher MS-CASPT2 level). According to the mapped potential energy profiles, we have identified two ultrafast excited-state deactivation pathways for the initially populated {sup 1}ππ{sup ∗} system. The first is the diabatic ESIPT process along the {sup 1}ππ{sup ∗} potential energy profile. The generated {sup 1}ππ{sup ∗} keto species then decays to the S{sub 0} state via the keto {sup 1}ππ{sup ∗}/gs conical intersection. The second is internal conversion to the dark {sup 1}nπ{sup ∗} state near the {sup 1}ππ{sup ∗} /{sup 1}nπ{sup ∗} crossing point in the course of the diabatic {sup 1}ππ{sup ∗} ESIPT process. Our following dynamics simulations have shown that the ESIPT and {sup 1}ππ{sup ∗} → S{sub 0} internal conversion times are 104 and 286 fs, respectively. Finally, our present work demonstrates that in addition to the ESIPT process and the {sup 1}ππ{sup ∗} → S{sub 0} internal conversion in the keto region, the {sup 1}ππ{sup ∗} → {sup 1}nπ{sup ∗} internal conversion in the enol region plays as well an important role for the excited-state relaxation dynamics of oxybenzone.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Enhancing Least-Squares Finite Element Methods Through a Quantity-of-Interest

Here, we introduce an approach that augments least-squares finite element formulations with user-specified quantities-of-interest. The method incorporates the quantity-of-interest into the least-squares functional and inherits the global approximation properties of the standard formulation as well as increased resolution of the quantity-of-interest. We establish theoretical properties such as optimality and enhanced convergence under a set of general assumptions. Central to the approach is that it offers an element-level estimate of the error in the quantity-of-interest. As a result, we introduce an adaptive approach that yields efficient, adaptively refined approximations. Several numerical experiments for a range of situations are presented to support the theory and highlight the effectiveness of our methodology. Notably, the results show that the new approach is effective at improving the accuracy per total computational cost.

97 MATHEMATICS AND COMPUTING↗

Technology Strategy Assessment: Findings from Storage Innovations 2030 Bidirectional Hydrogen Storage

Hydrogen is the most common element in the universe, comprising nearly 75% of all normal matter, and it has been used by scientists for centuries, but it was not fully recognized as an element until 1766, when it was isolated by Henry Cavendish. Early work focused on the generation of hydrogen through the oxidation of metals in water, which released hydrogen gas. Hydrogen’s lighter-than-air and flammable properties were immediately used in engines, zeppelins, and as feedstock for a wide variety of chemical reactions. Several approaches were developed for the production of hydrogen with the most common being associated with the production and conversion of hydrocarbon-based fuels. Coal gasification, steam methane reforming, and other reformation processes provide the majority of current hydrogen production due to the relatively low cost of hydrogen produced through these processes. More than 95% of hydrogen production is used for industrial processes rather than energy storage. To facilitate affordable decarbonization of these industrial processes and to advance the use of hydrogen as a fuel in transportation, DOE launched the Hydrogen Shot as part of the Energy Earthshots Initiative. The goal of the Hydrogen Shot is to reduce the cost of clean hydrogen by 80% to $1/kg of clean hydrogen production within one decade (known as the “1 1 1” goal). This is distinct from the Long-Duration Storage Shot, which is the primary focus of this report; however, it is intrinsically linked to bidirectional hydrogen storage. Several important chemical synthesis processes are dependent upon hydrogen, and the production and use of hydrogen is generally driven by its connection to one of these markets. For example, ammonia is one of the most highly produced chemicals in the world and it depends chiefly on hydrogen. Ammonia is primarily used for agricultural fertilizer and is considered to be largely responsible for a doubling of agricultural production per unit of land over the last century. Another one of hydrogen’s primary uses is as a catalyst in petroleum refining during the desulfurization process. Beyond chemical production, hydrogen is used as a reductant in the production of steel and has been demonstrated as a substitute for metallurgical coal in the production of raw iron. It is even used in the hydrogenation reaction for food products to create more shelf-stable semi-solid fats. However, while hydrogen is produced on the order of 100 million metric tons/year globally to feed these industries, more than 95% of hydrogen is produced from hydrocarbons that emit CO2 during the process. Conversely, electrolysis is a process by which electricity is used to separate hydrogen and oxygen in water molecules, usually across a membrane. Hydrogen production via electrolysis lowers the carbon intensity of produced hydrogen when coupled with low-carbon electricity. Currently, global electrolysis capacity is on the order of 1 GW, which equates to about 500 metric tons/day of hydrogen production. To support large-scale industrial decarbonization, capacity will likely need to increase by two to three orders of magnitude. Electrolysis technology is broadly separated into groups that are defined by the electrolyte used, with further subdivision based on the operating characteristics. The majority of commercial electrolyzer systems are based around three main technology groups: liquid alkaline, proton exchange membrane, and solid oxide. Liquid Alkaline (LA) electrolysis is the oldest, most mature, least expensive, and most common commercial technology, with 400 plants in operation by 1902. Its hydrogen output is low relative to the size of the system due to a low current density. LA electrolysis utilizes a liquid potassium hydroxide solution as the electrolyte. Proton exchange membrane (PEM) electrolysis (also known as polymer electrolyte membrane electrolysis), described in 1960, relies on an acid-impregnated polymer membrane as the electrolyte and typically offers three to six times higher hydrogen production per unit cell area than LA electrolysis. Solid oxide electrolysis, or high-temperature electrolysis, utilizes a ceramic cell as the electrolyte and operates on steam rather than liquid water, enabling electrical efficiencies of more than 90%, which is up from 60% with PEM. Two pre-commercial electrolyzer technologies to note are alkaline exchange membrane (AEM) and proton-conducting solid oxide electrolysis cell (SOEC). AEM potentially has the advantages of both LA and PEM technologies in that it is able to use low-cost materials like LA but with the ability to operate at higher output pressures with a smaller footprint like PEM. Proton-conducting SOEC is similar to commercial SOEC, which uses an oxide-conducting ceramic; however, it uses a proton-conducting ceramic that has the potential to operate at lower temperatures and has lower capital costs. Each of these technologies is experiencing a rapid improvement in performance and a reduction in installed cost, and each appears to be well suited to specific applications. Besides differences in the type of electrolyzer used, the main difference in the architecture of bidirectional hydrogen systems is how the hydrogen is stored. Currently, the most cost-effective way to store large amounts of hydrogen gas is underground, such as in large salt caverns that have been hollowed out. These salt caverns are geographically concentrated in small portions of the United States and are not generally near large metropolitan areas; however, other subsurface architectures are being investigated to expand this reach. A more widely deployable option is aboveground pressurized tanks. These systems are about 10 times as expensive because of the materials and safety margins required to hold hydrogen at high pressures. A third option is using materials-based storage, such as liquid organic hydrogen carriers. By reversibly attaching the produced hydrogen to other molecules, it can be stored at near atmospheric pressure and room temperature. This has the potential to reduce the material cost of storage but may result in a reduction in the efficiency of the process because there are both hydrogen uptake and release processes. While materials-based storage has not been used extensively for large-scale hydrogen storage in the past, there is currently significant activity regarding developing materials and processes for use in large-scale hydrogen storage applications. Electrolysis-produced hydrogen offers an unusual opportunity for energy storage applications. Unlike more conventional energy storage approaches, such as batteries, which operate entirely within electrical markets, hydrogen is a valuable product beyond the electric market and can be directed to the most lucrative use. Hydrogen also can be directly converted back to electricity using either a fuel cell or turbine, or it can be sold to other markets, such as chemical synthesis, steel production, or even export. In this way, excess electricity can be upgraded to the most valuable product. Finally, its use can be actively managed between multiple off-takers; for example, local hydrogen storage can provide a specific amount of stored electricity and any excess can be exported to ammonia production. This flexibility is amplified by the fact that hydrogen storage has fully decoupled power and energy components, which allows for affordable scaling options. Together, this allows a substantial amount of creativity to enable the economic utilization of variable power resources while supporting decarbonization of the industry.

08 HYDROGEN↗

A simplified hydrocarbon reaction mechanism for combustion applications

A simplified chemical kinetic reaction mechanism for the combustion of a hydrocarbon fuel is presented and discussed. The observed kinetic behavior of propane combustion as determined in shock tube studies and the calculated kinetic behavior of propane/air mixtures as derived from a detailed propane combustion mechanism were used as a data base for constructing and refining the mechanism. The mechanism contains 13 chemical species and 26 chemical reactions. Numerical studies revealed that the simplified mechanism can reproduce the observed kinetic behavior of propane/air mixtures, including laminar flame speeds, over a wide range of temperatures and pressures. It is anticipated that the mechanism can be used in various combustion applications and provide a more realistic description of the combustion process than can be achieved using global reaction mechanisms.

Jachimowski, C. J.↗

Visualization of High Latitude Ion Upflow in Support of the Image Mission

The study of the magnetosphere is a 400 year old science that began with the publication by Gilbert, in 1600, of his hypotheses that the Earth was a giant magnet. Since then we have learned many things about the magnetosphere, particularly in the last 40 years of the space age, but we still have many unanswered questions. In spite of the many thousands of observations of this system we still lack a global understanding of how it works. This is due to its large size and tenuous nature that mean that any measurement made of the fields or particles involved only give one a knowledge of the local conditions at a given time. To gain a global perspective through such observations would require the simultaneous operation of thousands of satellites spread throughout the magnetospheric system in addition to observations made on the ground. Such a program would be impractical at least from financial considerations. What is needed for the advancement of magnetospheric physics is to develop the same capabilities that astrophysicists, solar physicists and meteorologists have been using for years --- the ability to stand back from the object under study and see it in its entirety. The challenge for doing this for the magnetosphere is that the particle densities are very low and the material is, for the most part, not luminous. In the last 25 years several ideas have been proposed that would allow at least the imaging of certain portions of the magnetosphere. These include imaging of the plasmasphere through the resonant scattering of solar 304 A from He+ ions, imaging of various hot plasma populations (i.e. the ring current, plasmasheet, upflowing ionospheric ions, etc.) from the neutral atoms that result when ions of these populations charge exchange with the hydrogen geocorona, and imaging the aurora at various wavelengths in the far ultraviolet. In addition, a novel technique for probing various boundaries in the magnetosphere by bouncing low frequency radio waves off of them has been extensively studied. Such a technique is analogous to the way the under water world can be probed with sonar. About five years ago NASA convened a science working group to study the possibility of flying a magnetospheric imaging mission. This resulted in a number of proposals for such a mission, one of which was selected to be the first MIDEX mission, to be launched in early 2000. The mission is called IMAGE (Imager for Magnetopause to Aurora Global Exploration) and its P.I. is J. Burch at SwRI. The IMAGE spacecraft will carry imagers to view the plasmasphere, aurora, ring current, inner plasmasheet, and upflowing ionospheric ions as well as a radio sounder to probe the location, shape and dynamics of the magnetopause, plasmapause, etc. Between its selection last April and the non advocacy mission review, which takes place next spring, the IMAGE teams needs to further refine the design of the mission and its instruments. The theory and modeling (T&M) subgroup of this team has the task of demonstrating what kind of images the instruments on IMAGE will see as well as showing that useful scientific information can be extracted from such images. As a central element to the efforts of the T&M subgroup we have decided to simulate and create synthetic images for the magnetic cloud event of October, 1995. In this event a large cloud, with high plasma densities and strong magnetic fields, ejected from the sun collided with the earth's magnetosphere triggering a three day period of intense magnetic storms and substorms. This event was observed from a number of different spacecraft and on the ground so we have a good data set to work with. In our work we will place the IMAGE spacecraft in the magnetosphere on its proposed orbit, with its proposed instruments, to see what it would see had it been there. Existing models of the plasmasphere, ring current and magnetopause will be run for this event to give the structures for the imaging instruments. There are several models which are lacking and which need to be developed. These include a model for the cusp, the inner plasmasheet and the upflowing ions. My task this summer was to develop the upflowing ion model and use it to create synthetic images.

Wilson, Gordon R.↗

On Europa's Magnetospheric Interaction: A MHD Simulation of the E4 Flyby

The global three-dimensional interaction of Europa with the Jovian magnetosphere is modeled by using a complete set of ideal magnetohydrodynamic (MHD) equations. The model accounts for exospheric mass loading, ion-neutral charge exchange, recombination, and a possible intrinsic dipole magnetic field of Europa. The single-fluid MHD equations are solved by using a modem, finite volume, higher-order, Godunov-type method on an adoptively refined unstructured grid, which allows detailed modeling of the region near Europa while still resolving both the upstream region and the satellite's wake. The magnetic field and plasma density measured during Galileo's EGA flyby of December 19, 1996, are reproduced reasonably well in the simulation. We find the agreement between the data and our model particularly convincing if we assume that the plasma velocity during the EGA flyby deviated from the nominal corotation direction by approximately 20 deg. Evidence from the Galileo energetic particle detector also supports this assumption. In this case, we can fit the data using a dipole with orientation close to that of an induced dipole arising from the interaction of a hypothetical conducting subsurface layer on Europa with the periodically changing magnetic field of Jupiter. However, the magnitude of the dipole in our model is somewhat smaller (70%) than that suggested by Khurana et al. The total mass loading and ion-neutral charge exchange rates are consistent with the estimates of Europa's atmosphere and ionosphere.

Kabin, K.↗

Aerocapture Design Study for a Titan Polar Orbiter

In 2014 a team at NASA Goddard Space Flight Center (GSFC) studied the feasibility of using active aerocapture to reduce the chemical Delta V requirements for inserting a small scientific satellite into Titan polar orbit. The scientific goals of the mission would be multi-spectral imaging and active radar mapping of Titan's surface and subsurface. The study objectives were to: (i) identify and select from launch window opportunities and refine the trajectory to Titan; (ii) study the aerocapture flight path and refine the entry corridor; (iii) design a carrier spacecraft and systems architecture; (iv) develop a scientific and engineering plan for the orbital portion of the mission. Study results include: (i) a launch in October 2021 on an Atlas V vehicle, using gravity assists from Earth and Venus to arrive at Titan in January 2031; (ii) initial aerocapture via an 8-km wide entry corridor to reach an initial 350X6000 km orbit, followed by aerobraking to reach a 350X1500 km orbit, and a periapse raise maneuver to reach a final 1500 km circular orbit; (iii) a three-part spacecraft system consisting of a cruise stage, radiator module, and orbiter inside a heat shield; (iv) a 22-month mission including station keeping to prevent orbital decay due to Saturn perturbations, with 240 Gb of compressed data returned. High-level issues identified include: (i) downlink capability - realistic downlink rates preclude the desired multi-spectral, global coverage of Titan's surface; (ii) power - demise of the NASA ASRG (Advanced Stirling Radioisotope Generator) program, and limited availability at present of MMRTGs (Multi-Mission Radioisotope Generators) needed for competed outer planet missions; (iii) thermal - external radiators must be carried to remove 4 kW of waste heat from MMRTGs inside the aeroshell, requiring heat pipes that pass through the aeroshell lid, compromising shielding ability; (iv) optical navigation to reach the entry corridor; (v) the NASA requirement of continuous critical event coverage for the orbiter, especially during the peak heating of the aerocapture when the radio link will be broken. In conclusion, although Titan aerocapture allows for considerable savings in propellant mass, this comes at a cost of increased mission complexity. Further architecture study and refinement is required to reduce high-level mission risks and to elucidate the optimum architecture.

Nixon, Conor A.↗

A novel multistage estimation of signal parameters

A multistage estimation scheme is presented for estimating the parameters of a received carrier signal possibly phase-modulated by unknown data and experiencing very high Doppler, Doppler rate, etc. Such a situation arises, for example, in the case of the Global Positioning Systems (GPS). In the proposed scheme, the first-stage estimator operates as a coarse estimator of the frequency and its derivatives, resulting in higher rms estimation errors but with a relatively small probability of the frequency estimation error exceeding one-half of the sampling frequency (an event termed cycle slip). The second stage of the estimator operates on the error signal available from the first stage, refining the overall estimates, and in the process also reduces the number of cycle slips. The first-stage algorithm is a modified least-squares algorithm operating on the differential signal model and referred to as differential least squares (DLS). The second-stage algorithm is an extended Kalman filter, which yields the estimate of the phase as well as refining the frequency estimate. A major advantage of the is a reduction in the threshold for the received carrier power-to-noise power spectral density ratio (CNR) as compared with the threshold achievable by either of the algorithms alone.

Kumar, Rajendra↗

On the hydrodynamic interaction of shock waves with interstellar clouds. 1: Nonradiative shocks in small clouds

The interstellar medium (ISM) is inhomogeneous, with clouds of various temperatures and densities embedded in a tenuous intercloud medium. Shocks propagating through the ISM can ablate or destroy the clouds, at the same time significantly altering the properties of the intercloud medium. This paper presents a comprehensive numerical study of the simplest case of the interaction between a shock wave and a spherical cloud, in which the shock far from the cloud is steady and planar, and in which radiative losses, thermal conduction, magnetic fields, and gravitational forces are all neglected. As a result, the problem is completely specified by two numbers: the Mach number of the shock, M, and the ratio of the density of the cloud to that of the intercloud medium, Chi. For strong shocks we show that the dependence on M scales out, so the primary independent parameter is Chi. Variations from this simple case are also considered: the potential effect of radiative losses is assessed by calculations in which the ratio of specific heats in the cloud is 1.1 instead of 5/3; the effect of the initial shape of the cloud is studied by using a cylindrical cloud instead of a spherical one; and the role of the initial shock is determined by considering the case of a cloud embedded in a wind. Local adaptive mesh refinement techniques with a second-order, two-fluid, two-dimensional Godunov hydrodynamic scheme are used to address these problems, allowing heretofore unobtainable numerical resolution. Convergence studies to be described in a subsequent paper demonstrate that about 100 zones per cloud radius are needed for accurate results; previous calculations have generally used about a third of this number. The results of the calculations are analyzed in terms of global quantities which provide an overall description of te shocked cloud: the size and shape of the cloud, the mean density, the mean pressure, the mean velocity, the velocity dispersion, and the total circulation.

Klein, Richard I.↗

Formation and Reconnection of Three-Dimensional Current Sheets in the Solar Corona

Current-sheet formation and magnetic reconnection are believed to be the basic physical processes responsible for much of the activity observed in astrophysical plasmas, such as the Sun s corona. We investigate these processes for a magnetic configuration consisting of a uniform background field and an embedded line dipole, a topology that is expected to be ubiquitous in the corona. This magnetic system is driven by a uniform horizontal flow applied at the line-tied photosphere. Although both the initial field and the driver are translationally symmetric, the resulting evolution is calculated using a fully three-dimensional magnetohydrodynamic (3D MHD) simulation with adaptive mesh refinement that resolves the current sheet and reconnection dynamics in detail. The advantage of our approach is that it allows us to apply directly the vast body of knowledge gained from the many studies of 2D reconnection to the fully 3D case. We find that a current sheet forms in close analogy to the classic Syrovatskii 2D mechanism, but the resulting evolution is different than expected. The current sheet is globally stable, showing no evidence for a disruption or a secondary instability even for aspect ratios as high as 80:1. The global evolution generally follows the standard Sweet- Parker 2D reconnection model except for an accelerated reconnection rate at a very thin current sheet, due to the tearing instability and the formation of magnetic islands. An interesting conclusion is that despite the formation of fully 3D structures at small scales, the system remains close to 2D at global scales. We discuss the implications of our results for observations of the solar corona. Subject Headings: Sun: corona Sun: magnetic fields Sun: reconnection

Edmondson, J. K.↗

Refinement of the Transition-edge Sensor Design for ATHENA X-IFU

The X-ray Integral Field Unit (X-IFU) instrument on the Advanced Telescope for High ENergy Astrophysics (ATHENA) is baselined to have 2376 transition-edge sensor (TES) microcalorimeter pixels in a single array. The required performance for X-IFU has been demonstrated on a kilo-pixel array of square TES’s with 50 m side length, and this is considered the baseline pixel design. However, over the last few years we have explored small modifications to this design in search of a globally optimized instrument performance. We have previously reported on investigations of extending the length of the TES’s and variations in the number of X-ray absorber support stems. Here we report on investigations of TES designs with a length of 50 m but narrower width. We will discuss the consequence of this on magnetic field sensitivity, resistive transition parameters, spectral performance, and ease of multiplexing. We will also discuss how the number and position of the absorber attachments influences the vibrational modes of the pixels, and the impact this may have on performance. Finally, we will present measurements of more substantial change options of our TES design, including more extreme TES geometries, the use of metal islands or etched holes in the silicon nitride membranes, and the position of wiring to minimize current-induced magnetic field effects. The exploration of all these changes to the design are not only useful for optimizing performance on X-IFU, but also guide our fundamental understanding of the key physics of TES microcalorimeters.

Nick Wakeham↗

Structural investigation of the temperature-stable relaxor dielectric Ba 0.8 Ca 0.2 TiO 3 -Bi(Mg 0.5 Ti 0.5 )O 3

Aberration corrected scanning transmission electron microscopy (STEM) and electron diffraction have been used to disclose local structure and nano-chemistry in a Ca modified BaTiO 3 -Bi(Mg 0.5 Ti 0.5 )O 3 relax or dielectric ceramic which exhibits high and near-invariant relative permittivity over a wide temperature range. High resolution, synchrotron X-ray diffraction indicated a globally cubic structure ($Pm$$\overline{3}$$m$), but direct atomic-scale imaging by STEM revealed local tetragonal distortions. Nanopolar clusters were identified from B-site atomic displacement vectors measured relative to oxygen ion positions along < 100 > and < 110 > zone axes of integrated differential phase contrast (iDPC) STEM images, highlighting cluster sizes of 2–5 nm. Chemical analysis by STEM-energy dispersive X-ray spectroscopy and full pattern refinements of X-ray powder diffraction data each implied high levels of Bi vacancies within the matrix. The possibility that these A-site vacancies modulate the nanopolar structure and promote flattening of the permittivity-temperature response in this class of dielectric is discussed.

36 MATERIALS SCIENCE↗

EMMA: an adaptive mesh refinement cosmological simulation code with radiative transfer

EMMA is a cosmological simulation code aimed at investigating the reionization epoch. It handles simultaneously collisionless and gas dynamics, as well as radiative transfer physics using a moment-based description with the M1 approximation. Field quantities are stored and computed on an adaptive three-dimensional mesh and the spatial resolution can be dynamically modified based on physically motivated criteria. Physical processes can be coupled at all spatial and temporal scales. We also introduce a new and optional approximation to handle radiation: the light is transported at the resolution of the non-refined grid and only once the dynamics has been fully updated, whereas thermo-chemical processes are still tracked on the refined elements. Such an approximation reduces the overheads induced by the treatment of radiation physics. A suite of standard tests are presented and passed by EMMA, providing a validation for its future use in studies of the reionization epoch. The code is parallel and is able to use graphics processing units (GPUs) to accelerate hydrodynamics and radiative transfer calculations. Depending on the optimizations and the compilers used to generate the CPU reference, global GPU acceleration factors between ×3.9 and ×16.9 can be obtained. As a result, vectorization and transfer operations currently prevent better GPU performance and we expect that future optimizations and hardware evolution will lead to greater accelerations.

79 ASTRONOMY AND ASTROPHYSICS↗

Separating Spatial and Temporal Variations of the Aurora Using Two Nearly Colocated Satellites

This final report describes the efforts accomplished during the grant's period of performance, covering the period of 1 May 1997 to 30 April 2001, of a NASA Supporting Research and Technology Program grant under the Ionospheric, Thermospheric, and Mesospheric Physics component of the Sun-Earth Connections program. We have met and exceeded the goals set forth in the proposed research objectives. Referred publications have appeared in the scientific literature and several others are in the review process. In addition, numerous invited and contributed presentations of these studies were presented at national and international meetings during the performance period. One graduate student completed his PhD and won two AGU Best Student Paper awards based on research funded by this grant. These studies are summarized below. The science goal delineated in the initial proposal was "to systematically explore the temporal and spatial characteristics of the aurora in a way heretofore impossible, using data from two coplanar DMSP spacecraft." We accomplished this goal through a series of related studies. One study used these unique data to establish the role of Ps6 waves in coupling between the magnetosphere and the auroral ionosphere (omega bands) during the recovery phase of a magnetic storm; the published paper demonstrated the causal relationships between geospace processes occurring in different regions and established a simple conceptual model based on the fortuitous constellation of observations. In the second string of papers, we used these data to explore velocity-dispersed ions (VDIS) in and near the cusp, to test region identification models, and to look at space/time structure of auroral precipitation. On the first topic, the unique DMSP data revealed a remarkable double VDIS with a latitudinal overlap. This could only be explained in terms of a unified reconnection geometry that builds on several earlier unrelated models. The paper outlining this discovery has drawn considerable attention from the community and is currently in press - it adds significantly to the debate over whether reconnection is study state versus bursty and patchy versus global. The second paper develops the model further by incorporating the electron signature - these ionospheric particle precipitation signatures reveal the presence of magnetospheric "fossilized" FTEs, demonstrating the power of ionospheric measurements as a remote diagnostic of magnetospheric processes. Finally, the general nature of aurora] stability and coherence and region identification by particle characteristics were fully explored in a final paper. We identify candidate mechanisms controlling coherence time scales and length scales and refine boundary region identification criteria. We also use the dual-DMSP observations to identify the open and closed LLBL region and related its significance to the generalized bursty, multiple x-line model developed in the first paper. All of these topics are chapters of Dr. Boudouridis' recently completed PhD thesis.

Spence, Harlan E.↗

Sensitivity of Mountain Hydroclimate Simulations in Variable-Resolution CESM to Microphysics and Horizontal Resolution

Mountains are natural dams that impede atmospheric moisture transport and water towers that cool, condense, and store precipitation. They are essential in the western United States where precipitation is seasonal, and snowpack is needed to meet water demand. With anthropogenic climate change increasingly threatening mountain snowpack, there is a pressing need to better understand the driving climatological processes. However, the coarse resolution typical of modern global climate models renders them largely insufficient for this task, and signals a need for an advanced strategy. This paper continues the assessment of variable-resolution in the Community Earth System Model (VR-CESM) in modeling mountain hydroclimatology to understand the role of grid-spacing at 55, 28, 14, and 7 km and microphysics, specifically the Morrison and Gettelman (2008, MG1, https://doi.org/10.1175/2008JCLI2105.1) scheme versus the Gettelman and Morrison (2015, MG2, https://doi.org/10.1175/JCLI-D-14-00102.1) scheme. Eight VR-CESM simulations were performed from 1999 to 2015 with the F_AMIP_CAM5 component set, which couples the atmosphere-land models and prescribes ocean data. Refining horizontal grid-spacing from 28 to 7 km with the MG1 scheme did not improve the simulated mountain hydroclimatology. Substantial improvements occurred with the use of MG2 at grid-spacings ≤28 km compared to MG1 as shown with subsequent statistics. Average SWE bias diminished by 9.4X, 4.9X, and 3.5X from 55 to 7 km. The range in minimum (maximum) DJF spatial correlations increased by 0.1–0.2 in both precipitation and SWE. Mountain windward/leeward distributions and elevation profiles improved across hydroclimate variables, however not always with model resolution alone. Disconcertingly, all VR-CESM simulations exhibited a systemic mountain cold bias that worsened with elevation and will require further examination.

58 GEOSCIENCES↗

Quantifying the Impact of Vertical Resolution on the Representation of Marine Boundary Layer Physics for Global-Scale Models

Abstract While GCM horizontal resolution has received the majority of scale improvements in recent years, ample evidence suggests that a model’s vertical resolution exerts a strong control on its ability to accurately simulate the physics of the marine boundary layer. Here we show that, regardless of parameter tuning, the ability of a single-column model (SCM) to simulate the subtropical marine boundary layer improves when its vertical resolution is improved. We introduce a novel objective tuning technique to optimize the parameters of an SCM against profiles of temperature and moisture and their turbulent fluxes, horizontal winds, cloud water, and rainwater from large-eddy simulations (LES). We use this method to identify optimal parameters for simulating marine stratocumulus and shallow cumulus. The novel tuning method utilizes an objective performance metric that accounts for the uncertainty in the LES output, including the covariability between model variables. Optimization is performed independently for different vertical grid spacings and value of time step, ranging from coarse scales often used in current global models (120 m, 180 s) to fine scales often used in parameterization development and large-eddy simulations (10 m, 15 s). Uncertainty-weighted disagreement between the SCM and LES decreases by a factor of ∼5 when vertical grid spacing is improved from 120 to 10 m, with time step reductions being of secondary importance. Model performance is shown to converge at a vertical grid spacing of 20 m, with further refinements to 10 m leading to little further improvement. Significance Statement In successive generations of computer models that simulate Earth’s atmosphere, improvements have been mainly accomplished by reducing the horizontal sizes of discretized grid boxes, while the vertical grid spacing has seen comparatively lesser refinements. Here we advocate for additional attention to be paid to the number of vertical layers in these models, especially in the model layers closest to Earth’s surface where climatologically important marine stratocumulus and shallow cumulus clouds reside. Our experiments show that the ability of a one-dimensional model to represent physical processes important to these clouds is strongly dependent on the model’s vertical grid spacing.

Meteorology & Atmospheric Sciences↗

Quantitative analysis of Indonesia’s reserves and energy security as an evaluation by the nation in facing global competition

Energy security has become a serious concern for all countries in the world and each country has its own definiton for measuring its energy security. The objective of this study was to measure energy security of Indonesia quantitatively by comparing it with other countries and provide some recommendations for enhancing the energy security. In this study, the database was developed from various sources and was cross-checked to confirm validity of the data. Then the parameters of energy security were defined, where all of data will be processed towards the selected parameters. These parameters (e.g. Primary Energy mix, TPES/capita, FEC/capita, Self Sufficiency, Refining capacity, Overseas Energy Resources, Resources diversification) are the standards used to produce an analysis or evaluation of national energy management. Energy balances for Indonesia and 10 selected countries (USA, Germany, Russia, England, Japan, China, South Korea, Singapore, Thailand and India) were presented from 2009 to 2013. With a base index of 1.0 for Indonesia, calculated energy security index capable of representing Indonesia energy security compared relatively to other countries were also presented and discussed in detail. In 2012, Indonesia security index is ranked 11 from 11 countries, while USA and South Korea are the highest with security index of 3.36 and 2.89, respectively. According to prediction for 2025, Indonesia energy security is ranked 10 from 11 countries with only Thailand has lower security index (0.98). This result shows that Indonesia energy security was vulnerable to crisis and must be improved. Therefore this study proposed some recommendations to improve Indonesia energy security. Indonesia need to increase oil production by constructing new refinery plants, developing infrastructure for energy distribution to reduce the potential of energy shortage and accelerating the utilization of renewable energy to reduce the excessive use of primary energy. From energy policy proposed in this study, Indonesia energy security for 2025 could be improved to ranked 8 of 11 countries, better than Malaysia, Thailand and Singapore.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗