Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “core model proposal”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Extending the Licklider Transmission Protocol to Multi-Band Links

Most deep space missions return data to Earth using links operating at a single frequency band. Indeed, their data requirements are low enough that bandwidth regulations do not constrain the system. In contrast, spacecraft such as Kepler or Europa Clipper are transitioning to a new operational paradigm where engineering and science data are transmitted through simultaneous links operating at different frequency bands (henceforth termed multi-band links). This ensures, for instance, that critical data is correctly received using a well characterized X-band link, while science data at a much larger data rate can be returned efficiently (both in terms of bandwidth and energy) through a Ka-band link.Having a spacecraft establish two simultaneous links with a ground station opens a large span of potential improvements for space communications and mission operations. In this paper, we consider the problem of running a Licklider Transmission Protocol (LTP) session over a multi-band link. LTP is an implementation of a selective Automatic Repeat reQuest (ARQ) protocol, i.e. it ensures correct delivery of data over an error prone link with potentially long propagation delays. To maximize its efficiency in deep space environments, LTP operates in deferred-ACK mode and is typically included as one of the core protocols in the Delay Tolerant Networking (DTN) suite.The contributions of this paper are as follows: First, we propose an extension to LTP for multi-band links (denoted MBLTP) and sketch how it can be implemented without modifying the current definition of LTP data units. Next, we develop bounds on the performance of MBLTP when transmitting a single data file over a multi-band link. Three metrics are considered, file expected delivery time, total energy spent and bundle jitter. The results of the analytic model are first benchmarked against simulations to ensure validity, and then compared against the performance of both traditional LTP and Parallel LTP (PLTP).We demonstrate that MBLTP can significantly reduce the latency and jitter with which data products are delivered to destination over a deep space link compared to LTP at moderate energy cost. Similarly, we also demonstrate that MBLTP outperforms PLTP in all considered metrics.

Sanchez Net, Marc↗
DOE OSTI2021

Computational design of microarchitected porous electrodes for redox flow batteries

Porous electrodes are used as the core reactive component across electrochemical technologies. In flowing systems, controlling the fluid distribution, species transport, and reactive environment is critical to attaining high performance. However, conventional electrode materials like felts and papers provide few opportunities for precise engineering of the electrode and its microstructure. To address these limitations, architected electrodes composed of unit cells with spatially varying geometry determined via computational optimization are proposed. Resolved simulation is employed to develop a homogenized description of the constituent unit cells. These effective properties serve as inputs to a continuum model for the electrode when used in the negative half-cell of a vanadium redox flow battery. Porosity distributions minimizing power loss are then determined via computational design optimization to generate architected porosity electrodes. The architected electrodes are compared to bulk, uniform porosity electrodes and found to lead to increased power efficiency across operating flow rates and currents. The design methodology is further used to generate a scaled-up electrode with comparable power efficiency to the bench-scale systems. Finally, the variable porosity architecture and computational design methodology presented here thus offers a novel pathway for automatically generating spatially engineered electrode structures with improved power performance.

25 ENERGY STORAGE↗
DOE OSTI2021

The role of edge plasma parameters in H-mode density limit on the JET-ILW

A study of a dataset of JET H-mode plasma with the Be/W ITER-like wall (JET-ILW) shows that reaching the edge MHD ballooning limit leads to confinement degradation. However, unlike JET plasmas with a carbon wall (JET-C), the JET-ILW plasmas stay in a marginal dithering phase for a relatively long period, associated with a higher (≈20%) H-mode density limit (HDL) than JET-C equivalents. This suggests that ITER could be operated in H-mode with higher density than the scaling based on carbon wall devices, but likely with a dithering phase plasma with lower confinement. A new, reliable estimator for JET E r, min has been derived by combining HRTS measurements of pedestal gradient and edge-SOL decay lengths. JET radial E r ETB wells are observed in the range of -15 to -60 kV m -1 in high performance H-modes, consistent with previous CXRS results in ASDEX Upgrade. The results imply that a higher positive E × B shear in the near SOL plays a role in sustaining a marginal phase in JET-ILW which leads to a higher HDL than that in JET-C. The results of the JET-ILW dataset show agreement with the Goldston finite collisionality HD model for SOL broadening at high collisionality. A hypothesis for the dithering H-mode phase is proposed: as n e,SOL increases, ν *,SOL increases, SOL broadens, E r shear decreases, triggers L-mode; n e drops, ν *,SOL decreases, SOL becomes narrower, and E r shear increases, triggering H-mode, resulting in a cycle of H–L–H- oscillations. For burning plasma devices, such as ITER, operating just below the MHD limit for the dithering phase could be a promising regime for maximising core density, and fusion performance while minimising plasma-material interaction. The oscillatory signal during the dithering phase could be used as a precursor of undesirable plasma performance for control purposes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Rare earth element enrichment in coal and coal-adjacent strata of the Uinta Region, Utah and Colorado

This study aims to quantify rare earth element enrichment within coal and coal-adjacent strata in the Uinta Region of central Utah and western Colorado. Rare earth elements are a subset of critical minerals as defined by the U.S. Geological Survey. These elements are used for a wide variety of applications, including renewable energy technology in the transition toward carbon-neutral energy. While rare earth element enrichment has been associated with Appalachian coals, there has been a more limited evaluation of western U.S. coals. Here, samples from six active mines, four idle/historical mines, four mine waste piles, and seven stratigraphically complete cores within the Uinta Region were geochemically evaluated using portable X-ray fluorescence ( n = 3,113) and inductively coupled plasma-mass spectrometry ( n = 145) elemental analytical methods. Results suggest that 24%–45% of stratigraphically coal-adjacent carbonaceous shale and siltstone units show rare earth element enrichment (>200 ppm), as do 100% of sampled igneous material. A small subset (5%–8%) of coal samples display rare earth element enrichment, specifically in cases containing volcanic ash. This study proposes two multi-step depositional and diagenetic models to explain the enrichment process, requiring the emplacement and mobilization of rare earth element source material due to hydrothermal and other external influences. Historical geochemical evaluations of Uinta Region coal and coal-adjacent data are sparse, emphasizing the statistical significance of this research. These results support the utilization of active mines and coal processing waste piles for the future of domestic rare earth element extraction, offering economic and environmental solutions to pressing global demands.

Coe, Haley H.↗
NASA NTRS1986

Lunar nitrogen: Secular variation or mixing?

The two current models to explain the nearly 40% variation of the lunar nitrogen isotopic composition are: (1) secular variation of solar wind nitrogen; and (2) a two component mixing model having a constant, heavy solar wind admixed with varying amounts of indigenous light lunar N (LLN). Both models are needed to explain the step pyrolysis extraction profile. The secular variation model proposes that the low temperature release is modern day solar wind implanted into grain surfaces, the 900 C to 1100 C release is from grain surfaces which were once exposed to the ancient solar wind but which are now trapped inside agglutinates, and the >1100 C release as spallogenic N produced by cosmic rays. The mixing model ascribes the components to solar wind, indigenous lunar N and spallogenic N respectively. An extension of either interpretation is that the light N seen in lunar breccias or deep drill cores represent conditions when more N-14 was available to the lunar surface.

Norris, S. J.↗
NASA NTRS1988

Arcs from gravitational lensing

The proposal made by Paczynski (1987) that the arcs of blue light found recently in two cluster cores are gravitationally lensed elongated images of background galaxies is investigated. It is shown that lenses that are circularly symmetric in projection produce pairs of arcs, in conflict with the observations. However, more realistic asymmetric lenses produce single arcs, which can become as elongated as the observed ones whenever the background galaxy is located on or close to a cusp caustic. Detailed computer simulations of lensing by clusters using a reasonable model of the mass distribution are presented. Elongated and curved lensed images longer than 10 arcsec occur in 12 percent of the simulated clusters. It is concluded that the lensing hypothesis must be taken seriously.

Grossman, Scott A.↗
NASA NTRSDate not supplied

Jarosite Formation in Deep Antarctic Ice, an Earth Analogue to Unveil the Origin of Jarosite on Mars?

On Earth, jarosite is a weathering product forming in acidic-oxidative environments from the alteration of iron-bearing minerals in presence of liquid water. Typical settings where this iron-potassium hydrated sulphate is found, are weathering zones of pyrite-rich deposits, evaporative basins and fumaroles. Jarosite is not only known on Earth, it also occurs on Mars, where it has been identified by the Opportunity rover. The mineral was in fact recognized in the finely layered formations outcropping at Meridiani Planum and that were accurately investigated by the rover (Klingelhöfer et al. 2004). Since jarosite requires liquid water to form, its occurrence on Mars has been regarded as an evidence for the occurrence of liquid water in the geologic past of Mars (Elwood-Madden et al., 2004). Since then, many models have been proposed to describe the environments where the precipitation of Martian jarosite took place. The most accepted ones deal with evaporative basins similar to Earth’s playas, others concern volcanic activity and hydrothermal processes. An alternative proposal predicted that jarosite may have formed as a consequence of weathering of mineral dust trapped in massive ice deposits, i.e. the ice-weathering model (Niles & Michalsky, 2009). The hypothesis that jarosite formed on Mars because of low-temperature, acidic and water limited weathering, is not new (Burns, 1987), but until now no direct evidences were available to support it.A potential Earth analogue to investigate such processes is deep Antarctic ice. We present a first investigation of deep ice samples from the Talos Dome ice core (East Antarctica) aimed at the identification of englacial jarosite, so as to support the ice-weathering model. Evidences gathered through independent techniques showed that jarosite is actually present in deep Antarctic ice and results from the weathering of dust trapped into ice. The process is controlled by the re-crystallization of ice grains and the concurrent re-location of impurities at grain-junctions, which both depend on ice depth. This study demonstrates that the deep englacial environment is suitable for jarosite precipitation. Our findings support the hypothesis that, as originally predicted by the ice-weathering model, paleo ice-related processes have been important in thegeologic and geochemical history of Mars.

Lindsay P Keller↗
DOE OSTI2017

Insufficient Energy From MgO Exsolution to Power Early Geodynamo

The origin of Earth's ancient magnetic field is an outstanding problem. It has recently been proposed that exsolution of MgO from the core may provide sufficient energy to drive an early geodynamo. Here, in this work, we present new experiments on Mg partitioning between iron-rich liquids and silicate/oxide melts. Our results indicate that Mg partitioning depends strongly on the oxygen content in the iron-rich liquid, in contrast to previous findings that it depends only on temperature. Consequently, MgO exsolution during core cooling is drastically reduced and insufficient to drive an early geodynamo alone. Using the new experimental data, our thermal model predicts inner core nucleation at ~850 Ma and a nearly constant paleointensity.

58 GEOSCIENCES↗
DOE OSTI2016

Adsorption of ethylene on Sn and In terminated Si(001) surface studied by photoelectron spectroscopy and scanning tunneling microscopy

Interaction of ethylene (C{sub 2}H{sub 4}) with Si(001)-Sn-2 × 2 and Si(001)-In-2 × 2 at room temperature has been studied using core level (C 1s) X-ray photoelectron spectroscopy with synchrotron radiation and scanning tunneling microscopy. Sn and In form similar dimer chains on Si(001)2 × 1, but exhibit different interaction with ethylene. While ethylene adsorbs on top of Sn dimers of the Si(001)-Sn-2 × 2 surface, the Si(001)-In-2 × 2 surface turned out to be inert. Furthermore, the reactivity of the Sn terminated surface is found to be considerably decreased in comparison with Si(001)2 × 1. According to the proposed adsorption model ethylene bonds to Sn dimers via [2 + 2] cycloaddition by interacting with their π dimer bonds. In contrast, indium dimers do not contain π bonds, which renders the In terminated Si(001) surface inert for ethylene adsorption.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Why is there a dearth of close-in planets around fast-rotating stars?

We propose that the reported dearth of Kepler objects of interest (KOIs) with orbital periods P {sub orb} ≲ 2-3 days around stars with rotation periods P {sub rot} ≲ 5-10 days can be attributed to tidal ingestion of close-in planets by their host stars. We show that the planet distribution in this region of the log P {sub orb}-log P {sub rot} plane is qualitatively reproduced with a model that incorporates tidal interaction and magnetic braking as well as the dependence on the stellar core-envelope coupling timescale. We demonstrate the consistency of this scenario with the inferred break in the P {sub orb} distribution of close-in KOIs and point out a potentially testable prediction of this interpretation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Consistent SPH Simulations of Protostellar Collapse and Fragmentation

We study the consistency and convergence of smoothed particle hydrodynamics (SPH) as a function of the interpolation parameters, namely the number of particles N, the number of neighbors n, and the smoothing length h, using simulations of the collapse and fragmentation of protostellar rotating cores. The calculations are made using a modified version of the GADGET-2 code that employs an improved scheme for the artificial viscosity and power-law dependences of n and h on N, as was recently proposed by Zhu et al., which comply with the combined limit N→∞, h→0, and n→∞ with n/N→0 for full SPH consistency as the domain resolution is increased. We apply this realization to the “standard isothermal test case” in the variant calculated by Burkert and Bodenheimer and the Gaussian cloud model of Boss to investigate the response of the method to adaptive smoothing lengths in the presence of large density and pressure gradients. The degree of consistency is measured by tracking how well the estimates of the consistency integral relations reproduce their continuous counterparts. In particular, C {sup 0} and C {sup 1} particle consistency is demonstrated, meaning that the calculations are close to second-order accuracy. As long as n is increased with N, mass resolution also improves as the minimum resolvable mass M{sub min}∼n{sup −1}. This aspect allows proper calculation of small-scale structures in the flow associated with the formation and instability of protostellar disks around the growing fragments, which are seen to develop a spiral structure and fragment into close binary/multiple systems as supported by recent observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Effect of Reactor Radiation on the Thermal Conductivity of Graphite-based Dispersion Fuel

The Fukushima Daiichi nuclear accident in 2011 stimulated materials development efforts to improve the current fuel design for Light Water Reactors (LWRs). From fuel element to fuel cladding, new concepts in advanced nuclear fuel development have been proposed, studied, and tested. The new fuel systems are expected to have better resistance to damage, melting, and chemical reaction during accident scenarios. To test the performance and accident tolerance of the advanced fuel, it is highly preferable to simulate off-normal transients in addition to standard in-core fuel performance tests (available at High Flux Isotope Reactor at Oak Ridge National Laboratory and Advanced Test Reactor at Idaho National Laboratory). The Transient Reactor Test Facility (TREAT) is now being prepared to restart its operation, performing transient testing to assist advanced fuel development. TREAT was operated from 1959 to 1994, and was extensively used to generate thousands of reactor transients for different programs. The reactor has not operated since 1994. The original fuel used at TREAT was graphite-based dispersion fuel with highly-enriched uranium (HEU, 93.1% enriched uranium). For nonproliferation purposes, a fuel conversion program, supported by U.S. Department of Energy, National Nuclear Safety Administration (NNSA), Office of Material Management and Minimization (NA-23) Reactor Conversion Program, is underway to develop a replacement core in which the UO{sub 2} contains only low-enriched uranium (LEU). The preliminary LEU core design has shown a much lower carbon/uranium atomic ratio and enrichment than in the HEU core design. To achieve performance similar to the HEU core, the prospective LEU core will contain a much higher density of UO{sub 2} particles. As a result, a higher level of radiation damage, particularly from fission fragments emerging from the UO{sub 2} particles, can be expected dependent on the UO{sub 2} particle size within the fuel. To support LEU fuel design, the present study investigates the effect of reactor radiations on the thermal conductivity of graphite-based dispersion fuel. Both neutron- and fission-radiation induced degradations of the fuel are modeled, and validated based on the available experimental data. In this study, the reactor-radiation induced thermal conductivity degradation of graphite-based dispersion fuel was analyzed. A semi-empirical method to model the relative thermal conductivity with reactor radiations was proposed and validated. The proposed method can be further applied to predict the changes in the thermal conductivity of the prospective LEU TREAT fuel during long-term operation. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Integrated Dynamic Facade Control with an Agent-based Architecture for Commercial Buildings

Dynamic façades have significant technical potential to minimize heating, cooling, and lighting energy use and peak electric demand in the perimeter zone of commercial buildings, but the performance of these systems is reliant on being able to balance complex trade-offs between solar control, daylight admission, comfort, and view over the life of the installation. As the context for controllable energy-efficiency technologies grows more complex with the increased use of intermittent renewable energy resources on the grid, it has become increasingly important to look ahead towards more advanced approaches to integrated systems control in order to achieve optimum life-cycle performance at a lower cost. This study examines the feasibility of a model predictive control system for low-cost autonomous dynamic façades. A system architecture designed around lightweight, simple agents is proposed. The architecture accommodates whole building and grid level demands through its modular, hierarchical approach. Automatically-generated models for computing window heat gains, daylight illuminance, and discomfort glare are described. The open source Modelica and JModelica software tools were used to determine the optimum state of control given inputs of window heat gains and lighting loads for a 24-hour optimization horizon. Penalty functions for glare and view/ daylight quality were implemented as constraints. The control system was tested on a low-power controller (1.4 GHz single core with 2 GB of RAM) to evaluate feasibility. The target platform is a low-cost ($35/unit) embedded controller with 1.2 GHz dual-core cpu and 1 GB of RAM. Configuration and commissioning of the curtainwall unit was designed to be largely plug and play with minimal inputs required by the manufacturer through a web-based user interface. An example application was used to demonstrate optimal control of a three-zone electrochromic window for a south-facing zone. The overall approach was deemed to be promising. Further engineering is required to enable scalable, turnkey solutions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1979

How tidal heating in Io drives the Galilean orbital resonance locks

The mechanisms by which orbital resonance locks are maintained among Io, Europa and Ganymede are analyzed, and the effects of the dissipative tides in Jupiter and Io on their establishment and evolution are considered. According to the proposed model, initially all three satellites are in orbits far from the 2:1 commensurabilities or the three body lock. The tide raised on Io damps down the free eccentricity; only modest tidal heating occurs. Subsequently the dissipative tide raised on Jupiter by Io causes Io's orbit to spiral outwards; Io approaches the 2:1 commensurability with Europa. Io's forced eccentricity increases rapidly to a critical value, and thereafter the resonant interaction causes Europa's orbit to expand at half that of Io's orbit. A fluid core is probably formed as the result of tidal heating. Finally Europa approaches the 2:1 commensurability, angular momentum is transferred from Europa's orbit to Ganymede's, and a steady state is attained. On the basis of the observed three-body resonance amplitude (0.066 deg), it is concluded that three-body resonance may have formed less than 500 million years ago; a measurement of the rate of tidal heating in Io may determine the Q value of both Io and Jupiter. Upper and lower limits for the Jovian Q value of 2 million and 20,000 are suggested.

Yoder, C. F.↗
NASA NTRS1997

Pre-Launch Tasks Proposed in our Contract of December 1991

We propose, during the pre-EOS phase to: (1) develop, with other MODIS Team Members, a means of discriminating different major biome types with NDVI and other AVHRR-based data. (2) develop a simple ecosystem process model for each of these biomes, BIOME-BGC (3) relate the seasonal trend of weekly composite NDVI to vegetation phenology and temperature limits to develop a satellite defined growing season for vegetation; and (4) define physiologically based energy to mass conversion factors for carbon and water for each biome. Our final core at-launch product will be simplified, completely satellite driven biome specific models for net primary production. We will build these biome specific satellite driven algorithms using a family of simple ecosystem process models as calibration models, collectively called BIOME-BGC, and establish coordination with an existing network of ecological study sites in order to test and validate these products. Field datasets will then be available for both BIOME-BGC development and testing, use for algorithm developments of other MODIS Team Members, and ultimately be our first test point for MODIS land vegetation products upon launch. We will use field sites from the National Science Foundation Long-Term Ecological Research network, and develop Glacier National Park as a major site for intensive validation.

Running, Steven W.↗
NASA NTRS1998

Pre-Launch Tasks Proposed in our Contract of December 1991

We propose, during the pre-EOS phase to: (1) develop, with other MODIS Team Members, a means of discriminating different major biome types with NDVI and other AVHRR-based data; (2) develop a simple ecosystem process model for each of these biomes, BIOME-BGC; (3) relate the seasonal trend of weekly composite NDVI to vegetation phenology and temperature limits to develop a satellite defined growing season for vegetation; and (4) define physiologically based energy to mass conversion factors for carbon and water for each biome. Our final core at-launch product will be simplified, completely satellite driven biome specific models for net primary production. We will build these biome specific satellite driven algorithms using a family of simple ecosystem process models as calibration models, collectively called BIOME-BGC, and establish coordination with an existing network of ecological study sites in order to test and validate these products. Field datasets will then be available for both BIOME-BGC development and testing, use for algorithm developments of other MODIS Team Members, and ultimately be our first test point for MODIS land vegetation products upon launch. We will use field sites from the National Science Foundation Long-Term Ecological Research network, and develop Glacier National Park as a major site for intensive validation.

Source record↗
DOE OSTI2016

Collaborative Project: Understanding Climate Model Biases in Tropical Atlantic and Their Impact on Simulations of Extreme Climate Events

Recent studies have revealed that among all the tropical oceans, the tropical Atlantic has experienced the most pronounced warming trend over the 20th century. Many extreme climate events affecting the U.S., such as hurricanes, severe precipitation and drought events, are influenced by conditions in the Gulf of Mexico and the Atlantic Ocean. It is therefore imperative to have accurate simulations of the climatic mean and variability in the Atlantic region to be able to make credible projections of future climate change affecting the U.S. and other countries adjoining the Atlantic Ocean. Unfortunately, almost all global climate models exhibit large biases in their simulations of tropical Atlantic climate. The atmospheric convection simulation errors in the Amazon region and the associated errors in the trade wind simulations are hypothesized to be a leading cause of the tropical Atlantic biases in climate models. As global climate models have resolutions that are too coarse to resolve some of the atmospheric and oceanic processes responsible for the model biases, we propose to use a high-resolution coupled regional climate model (CRCM) framework to address the tropical bias issue. We propose to combine the expertise in tropical coupled atmosphere-ocean modeling at Texas A&M University (TAMU) and the coupled land-atmosphere modeling expertise at Pacific Northwest National Laboratory (PNNL) to develop a comprehensive CRCM for the Atlantic sector within a general and flexible modeling framework. The atmospheric component of the CRCM will be the NCAR WRF model and the oceanic component will be the Rutgers/UCLA ROMS. For the land component, we will use CLM modified at PNNL to include more detailed representations of vegetation and soil hydrology processes. The combined TAMU-PNNL CRCM model will be used to simulate the Atlantic climate, and the associated land-atmosphere-ocean interactions at a horizontal resolution of 9 km or finer. A particular focus of the model development effort will be to optimize the performance of WRF and ROMS over several thousand of cores by focusing on both the parallel communication libraries and the I/O interfaces, in order to achieve the sustained throughput needed to perform simulations on such fine resolution grids. The CRCM model will be developed within the framework of the Coupler (CPL7) software that is part of the NCAR Community Earth System Model (CESM). Through efforts at PNNL and within the community, WRF and CLM have already been coupled via CPL7. Using the flux coupler approach for the whole CRCM model will allow us to flexibly couple WRF, ROMS, and CLM with each model running on its own grid at different resolutions. In addition, this framework will allow us to easily port parameterizations between CESM and the CRCM, and potentially allow partial coupling between the CESM and the CRCM. TAMU and PNNL will contribute cooperatively to this research endeavor. The TAMU team led by Chang and Saravanan has considerable experience in studying atmosphere-ocean interactions within tropical Atlantic sector and will focus on modeling issues that relate to coupling WRF and ROMS. The PNNL team led by Leung has extensive expertise in atmosphere-land interaction and will be responsible for improving the land surface parameterization. Both teams will jointly work on integrating WRF-ROMS and WRF-CLM to couple WRF, ROMS, and CLM through CPL7. Montuoro of the TAMU Supercomputing Center will be responsible for improving the MPI and Parallel IO interfaces of the CRCM. Both teams will contribute to the design and execution of the proposed numerical experiments and jointly perform analysis of the numerical experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Performance and Accuracy Implications of Parallel Split Physics-Dynamics Coupling in the Energy Exascale Earth System Atmosphere Model

Simultaneous calculation of atmospheric processes is faster than calculating processes one at a time. This type of parallelism is beneficial or perhaps even necessary to provide good performance on modern supercomputers, which achieve faster performance through increased processor count rather than improved clock speed. The scalability of the Energy Exascale Earth System Model (E3SM) Atmosphere Model (EAM) is limited by the fluid dynamics which scales up to the number of mesh cells in the global mesh. In contrast, the suite of physics parameterizations in EAM is scalable up to the total number of physics columns, which is an order of magnitude greater than the number of mesh cells. A proposed solution to unlocking the greater potential performance from the physics suite is to solve the physics and dynamics in parallel. This work represents a first attempt at parallel splitting of the grid-scale fluid dynamics model and the subgrid-scale physics parameterizations in a global atmosphere model. We will demonstrate that switching to parallel physics-dynamics coupling extends the scalability of the EAM to up to 3 times the previous peak scalability limit and is up to 20% faster than the sequentially split coupling at the highest core counts and the same time step. Decadal simulations of both coupling approaches show very little impact to the model climate. This improved performance does not come without drawbacks, however. Parallel splitting requires a shorter time step and other modifications which largely offset performance gains. A mass fixer is required for conservation. Techniques for mitigating these issues are also discussed.

97 MATHEMATICS AND COMPUTING↗