Search NASA⌕ Search

SEARCH · Search NASA

Results for “Uncertainty Preservation”

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 127 records · Page 7

Enabling Autonomous Rover Science through Dynamic Planning and Scheduling

This paper describes how dynamic planning and scheduling techniques can be used onboard a rover to autonomously adjust rover activities in support of science goals. These goals could be identified by scientists on the ground or could be identified by onboard data-analysis software. Several different types of dynamic decisions are described, including the handling of opportunistic science goals identified during rover traverses, preserving high priority science targets when resources, such as power, are unexpectedly over-subscribed, and dynamically adding additional, ground-specified science targets when rover actions are executed more quickly than expected. After describing our specific system approach, we discuss some of the particular challenges we have examined to support autonomous rover decision-making. These include interaction with rover navigation and path-planning software and handling large amounts of uncertainty in state and resource estimations.

planning↗

Event generators for high-energy physics experiments

We provide an overview of the status of Monte-Carlo event generators for high-energy particle physics. Guided by the experimental needs and requirements, we highlight areas of active development, and opportunities for future improvements. Particular emphasis is given to physics models and algorithms that are employed across a variety of experiments. These common themes in event generator development lead to a more comprehensive understanding of physics at the highest energies and intensities, and allow models to be tested against a wealth of data that have been accumulated over the past decades. A cohesive approach to event generator development will allow these models to be further improved and systematic uncertainties to be reduced, directly contributing to future experimental success. Event generators are part of a much larger ecosystem of computational tools. They typically involve a number of unknown model parameters that must be tuned to experimental data, while maintaining the integrity of the underlying physics models. Making both these data, and the analyses with which they have been obtained accessible to future users is an essential aspect of open science and data preservation. It ensures the consistency of physics models across a variety of experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Spent Nuclear Fuel Exploratory Roadmaps - 20054

The U.S. Department of Energy manages nearly 2,500 metric tons of heavy metal (MTHM) of spent nuclear fuel (SNF) resulting from several decades of research, testing, and production reactors [DOE 1995]. This SNF is managed at the DOE Hanford Site in Washington State (Hanford), Idaho National Laboratory (INL), and Savannah River Site in South Carolina (SRS). From 1995 to 2004, the Department of Energy (DOE) made several key programmatic decisions, supported by appropriate documentation in accordance with the National Environmental Policy Act (NEPA) for the management of SNF [DOE 1995]. These decisions have provided an overarching framework for SNF management, as well as site-specific and SNF-related management decisions for the past two decades. In the years since these decisions were made, with the notable exception of the successful drying and packaging of the production reactor fuel at Hanford, a majority of the decisions have been largely unimplemented. Also, since these decisions were made, a number of changes that bear on considerations relative to the path forward have occurred. A SNF Exploratory Road map activity, identified reasonable alternate pathways for DoE's inventory of SNF. Three fundamental pathways were identified in this activity for the long-term management and disposition of DOE SNF: Direct Disposal Pathway, Existing Processing Pathway and Alternate Processing Pathway. Further actions would be required (e.g., evaluation of NEPA analysis, technology maturity evaluations, refinement of cost and risks, benefits and advantages relative to other alternate paths) to define and inform many key decisions that will result in the selection of the disposition pathways that are most beneficial to the US Government in dispositioning DoE's SNF inventory. DOE is currently managing all SNF safely; however, the age of DOE-owned SNF and facilities for storing and processing SNF, coupled with the uncertainty of the storage duration, necessitates decisions and actions to ensure that the infrastructure will be in place to ensure continued safe and effective long-term management and eventual disposal. The identification of reasonable alternate pathways took these considerations into account in an effort to proactively manage conditions that could challenge the safety of storing and managing DOE SNF over the time periods now contemplated and also to ensure that flexibility is preserved to ensure that the DOE SNF remains compatible with final disposition pathways when they become available. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Oxygen, Magnesium, and Aluminum Isotopes in the Ivuna CAI: Re-Examining High-Temperature Fractionations in CI Chondrites

CI chondrites are thought to approximate the bulk solar system composition since they closely match the composition of the solar photosphere. Thus, chemical differences between a planetary object and the CI composition are interpreted to result from fractionations of a CI starting composition. This interpretation is often made despite the secondary mineralogy of CI chondrites, which resulted from low-T aqueous alteration on the parent asteroid(s). Prevalent alteration and the relatively large uncertainties in the photospheric abundances (approx. +/-5-10%) permit chemical fractionation of CI chondrites from the bulk solar system, if primary chondrules and/or CAIs have been altered beyond recognition. Isolated olivine and pyroxene grains that range from approx. 5 microns to several hundred microns have been reported in CI chondrites, and acid residues of Orgueil were found to contain refractory oxides with oxygen isotopic compositions matching CAIs. However, the only CAI found to be unambiguously preserved in a CI chondrite was identified in Ivuna. The Ivuna CAI's primary mineralogy, small size (approx.170 microns), and fine-grained igneous texture classify it as a compact type A. Aqueous alteration infiltrated large portions of the CAI, but other regions remain pristine. The major primary phases are melilite (Ak 14-36 ), grossmanite (up to 20.8 wt.% TiO 2 ), and spinel. Both melilite and grossmanite have igneous textures and zoning patterns. An accretionary rim consists primarily of olivine (Fa 2-17 ) and low-Ca pyroxene (Fs 2-10 ), which could be either surviving CI2 material or a third lithology.

Frank, D. R.↗

Harmonized and high-quality datasets of aerosol optical depth at a US continental site, 1997-2018

Aerosol optical depth (AOD) characterizes the aerosol burden in the atmosphere, while its wavelength dependence is a sign of particle size. Long-term records of wavelength-resolved AOD with high quality and suitable continuity are required for climate change assessment. Typically, climate-related studies use AOD products provided by several, and perhaps different, ground-based instruments. The measurements from these instruments often have different accuracy and temporal resolution. To preserve the advantages of these products (high quality) and to reduce their disadvantages (patchy records), we generate a merged dataset obtained from several instruments in which a nearly-continuous AOD record is found at two wavelengths (500 and 870 nm) with high quality and high temporal resolution (1-min) for a 21-yr period (1997-2018). The combined dataset addresses: (1) varying data quality and resolution mismatch of the individual AOD records, and (2) the uncertainty of the merged AOD and its relevance for user-specified needs. The generated dataset will be beneficial for a wide range of applications including aerosol-radiation interactions and their impact in forcing climate change.

54 ENVIRONMENTAL SCIENCES↗

Sulfur Pellets Responses to a Bare and Steel Reflected Pulse of the Oak Ridge National Laboratory Health Physics Research Reactor

The experiments analyzed in this report were conducted at the Health Physics Research Reactor (HPRR), also known as the $\textit{Fast Burst Reactor}$. The reactor was designed and built at Oak Ridge National Laboratory (ORNL) in 1961. The HPRR was an unmoderated, unshielded fast reactor that used highly enriched uranium and molybdenum alloy as fuel. The reactor was initially sent to the Nevada Test Site in 1962, where it was used to evaluate radiation doses received as a result of the Hiroshima and Nagasaki bombings during World War II. A few years later, the reactor was sent back to ORNL to be part of the Dosimetry Application Research (DOSAR) facility shown in Figure 1, which included a reactor building shown on the left (west) of the picture and a control and laboratory building in the upper right corner (northeast). The critical assembly was used for numerous technical studies, including systems calibration, dosimetry, radiobiology of plants and animals, testing of radiation alarms, as well as teaching and training in radiation dosimetry and nuclear engineering. Between 1963 and 1987, the HPRR was operated for thousands of hours, achieving criticality close to 10,000 times and motivating many publications. The HPRR was decommissioned in 1987. The goal of this effort was to use historical data from operation of the HPRR to create a criticality accident alarm system (CAAS) benchmark to be included in the $\textit{International Handbook of Evaluated Criticality Safety Benchmark Experiments}$ (ICSBEP Handbook). A thorough inspection was performed of all available documentation and information available. The most promising experiments that were selected for evaluation were those described in the 1987 ORNL report entitled $\textit{Health Physics Research Reactor Reference Dosimetry}$, ORNL-6240. The report includes reference dosimetry results of the shielded and unshielded configurations of the HPRR after burst operations. Because of changes to the reactor positioning and storage systems that were made in 1985, the previous dosimetry reports became obsolete, and the newly designed experiments were needed to create the HPRR’s adjusted dosimetry data. The various results reported in ORNL-6240 include reference doses and dose equivalents from different conventions at different distances and elevations as determined using the detected neutron fluence and conversion factors. The HPRR neutron fluence was obtained through different methods, including sulfur pellet analysis and threshold detector unit data. Information about the HPRR spectrum was also obtained through Bonner sphere measurements. This benchmark is focused on a part of the measured sulfur fluences reported in Appendix H of ORNL-6240. Standard commercial sulfur pellets were placed at different distances from the HPRR centerline during burst operation and were activated due to the 32 S(n,p) 32 P reaction. The resulting 32 P activity was then measured and the information about the corresponding sulfur fluence and/or neutron dose could be extracted. Many of those measurements have 7 been performed with the HPRR in its bare configuration or with different shields (combinations of Lucite, concrete, steel). All the necessary, precise information about material and/or dimensions of the different shields was not found, so it was decided to focus only on the unshielded and steel-shielded configurations to minimize the benchmark uncertainty. A total of 31 cases (24 unshielded and 7 shielded cases at different positions) of sulfur fluence were selected before evaluation to develop the benchmark.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Evaluation of Oak Ridge National Laboratory Health Physics Research Reactor Operation Data for Critical Benchmark Creation [Abstract]

The Oak Ridge National Laboratory (ORNL) Health Physics Research Reactor (HPRR) was a research reactor designed and built at ORNL in 1961. The critical assembly was using a highly enriched uranium and molybdenum alloy as the fuel, and it could be operated in steady-state or burst modes. The reactor was used for about 25 years to produce a lot of publications related to dosimetry, radiobiology and radiation detectors testing before its decommissioning in 1987. In recent years, the idea of using legacy operation data from the to create a valuable critical accident alarm system shielding benchmark arose. Such a benchmark has been submitted to the International Criticality Safety Benchmark Experiment Project (ICSBEP) Technical Review Group for a potential inclusion in the 2022 version of the handbook. Another way to use the valuable data from the operation of the HPRR is to evaluate the feasibility of the creation of a subcritical or prompt supercritical benchmark for inclusion in the ICSBEP or the International Reactor Physics Experiments Evaluation Project (IRPhEP) handbooks. To initiate a burst, the HPRR had to be operated in a slightly subcritical state for a few minutes. Then, the insertion of the burst control rod would greatly increase the reactivity of the system and start the burst. No critical configuration of the HPRR critical assembly could be located. The only information available concerns stable subcritical and prompt supercritical states, found in a burst experiments’ logbook. In the recovered logbook pages, information about 8 different bursts is available. The information includes the rods positions before and during a burst, the recorded subcritical reactor period and reactivity, and the burst fission yield derived from the temperature elevation sulfur pellet irradiation analysis. By using the HPRR logbook information and the as-built drawings of the critical assembly, a highly detailed model of the HPRR was created with SCALE 6.2.4/KENO-VI. Eight KENO-VI models were created to replicate the sub-critical assembly configurations described in the eight bursts from the recovered logbook pages. KENO-VI calculates k eff and it can be linked to a reactivity value in cents by using the delayed neutron fraction B eff , also calculated by KENO-VI. KENO-VI can also be used to model the prompt super-critical configurations of the HPRR and to assess the similarity with the burst measurements by comparing the calculated k eff and the measured fission yields between each burst. Unfortunately, high uncertainty exist and the obtained discrepancies between experiments and calculation results are high, compromising the creation of a valuable critical benchmark from HPRR operation data. The reasons of the discrepancies and potential ways to solve them are explored.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

How May the Choice of Downscaling Techniques and Meteorological Reference Observations Affect Future Hydroclimate Projections?

Abstract We present an intercomparison of a suite of high‐resolution downscaled climate projections based on a six‐member General Circulation Model (GCM) ensemble from Coupled Models Intercomparison Project (CMIP6). The CMIP6 GCMs have been downscaled using dynamical and statistical downscaling techniques based on two meteorological reference observations over the conterminous United States. We use the regional climate model, RegCM4, for dynamical downscaling, double bias correction constructed analogs method for statistical downscaling, and Daymet and Livneh datasets as the reference observations for statistical training and bias‐correction. We evaluate the performances of downscaled data in both historical and future periods under the SSP585 scenario. While dynamical downscaling improves the simulation of some performance evaluation indices, it adds an extra bias in others, highlighting the need for statistical correction before its use in impact assessments. Downscaled datasets after bias‐correction compare exceptionally well with observations. However, the choice of downscaling techniques and the underlying reference observations influence the hydroclimate characteristics of downscaled data. For instance, the statistical downscaling generally preserves the GCMs climate change signal but overestimates the frequency of hot extremes. Similarly, simulated future changes are sensitive to the choice of reference observations, particularly for precipitation extremes that exhibit a higher projected increase in the ensembles trained and/or corrected by Daymet than Livneh. Overall, these results demonstrate that multiple factors, including downscaling techniques and reference observations, can substantially influence the outcome of downscaled climate projections and stress the need for a comprehensive understanding of such method‐based uncertainties.

54 ENVIRONMENTAL SCIENCES↗

Drill/borescope System for the Mars Polar Pathfinder

The primary goals of the Mars Polar Pathfinder (MPP) Discovery Mission are to characterize the composition and structure of Mars' north polar ice cap, and to determine whether a climate record may be preserved in layers of ice and dust. The MPP would land as close as possible to the geographic north pole of Mars and use a set of instruments similar to those used by glaciologists to study polar ice caps on Earth: a radar sounder, a drill/borescope system, and a thermal probe. The drill/borescope system will drill approximately 50 cm into the surface and image the sides of the hole at 10 micron resolution for compositional and stratigraphic analysis. Several uncertainties have guided the development of this instrument, and they are discussed.

Paige, D. A.↗

In-Situ Oxygen Isotopic Composition of Tagish Lake: An Ungrouped Type 2 Carbonaceous Chondrite

We have measured the oxygen isotopic composition of several components of Tagish Lake by ion microprobe. This meteorite constitutes the best preserved sample of C2 matter presently available for study. It presents two different lithologies (carbonate-poor and -rich) which have fairly comparable oxygen isotopic composition, with regard to both the primary or secondary minerals. For the olivine and pyroxene grains, their delta O-18 values range from - 10.5% to + 7.4% in the carbonate-poor lithology, with a mean Delta O-17 value of - 3.7 2.4%. In the carbonate-rich lithology, delta O-18 varies from - 7.9% to + 3.3%, and the mean Delta O-17 value is - 4.7 +/- 1.4%. Olivine inclusions (Fo(sub >99)) with extreme O-16-enrichment were found in both lithologies: delta O-18 = - 46.1 %, delta O-187= - 48.3% and delta O-18 = - 40.6%, delta O-17 = - 41.2% in the carbonate-rich lithology; delta O-18 = - 41.5%, delta O-17 = -43.4%0 in the carbonate-poor lithology. Anhydrous minerals in the carbonate-poor lithology are slightly more O-16-rich than in the carbonate-rich one. Four low-iron manganese-rich (LIME) olivine grains do not have an oxygen isotopic composition distinct from the other "normal" olivines. The phyllosilicate matrix presents the same range of oxygen isotopic compositions in both lithologies: delta O-18 from approximately 11 % to approximately 6%, with an average Delta. O-17 approximately 0%. Because the bulk Tagish Lake oxygen isotopic composition given by Brown et al. is on the high end of our matrix analyses, we assume that this "bulk Tagish Lake" composition probably only represents that of the carbonate-rich lithology. Calcium carbonates have delta O-18 values up to 35%, with Delta O-17 approximately 0.5%0. Magnetite grains present very high Delta O-17 values approximately + 3.4%0 +/- 1.2%. Given our analytical uncertainties and our limited carbonate data, the matrix and the carbonate seem to have formed in isotopic equilibrium. In that case, their large isotopic fractionation would argue for a low temperature (CM-like, T approximately 0 deg) formation. Magnetite probably formed during a separate event. Tagish Lake magnetite data is surprisingly compatible with that of R-chondrites and unequilibrated ordinary (LL3) chondrites. Our oxygen isotope data strongly supports the hypothesis of a single precursor for both lithologies. Drastic mineralogical changes between the two lithologies not being accompanied with isotopic fractionation seem compatible with the alteration model presented by Young et aI. Tagish Lake probably represents the first well preserved large sample of the C2 matter that dominates interplanetary matter since the formation of the solar system.

Zolensky, Michael E.↗

SCALE Code System

The SCALE Code System is a widely used modeling and simulation suite for nuclear safety analysis and design that is developed, maintained, tested, and managed by the Reactor and Nuclear Systems Division (RNSD) of Oak Ridge National Laboratory (ORNL). SCALE provides a comprehensive, verified and validated, user-friendly tool set for criticality safety, reactor and lattice physics, radiation shielding, spent fuel and radioactive source term characterization, and sensitivity and uncertainty analysis. Since 1980, regulators, licensees, and research institutions around the world have used SCALE for safety analysis and design. SCALE provides an integrated framework with dozens of computational modules, including three deterministic and three Monte Carlo radiation transport solvers that are selected based on the desired solution strategy. SCALE includes current nuclear data libraries and problem-dependent processing tools for continuous-energy (CE) and multigroup (MG) neutronics and coupled neutron-gamma calculations, as well as activation, depletion, and decay calculations. SCALE includes unique capabilities for automated variance reduction for shielding calculations, as well as sensitivity and uncertainty analysis. SCALE’s graphical user interfaces assist with accurate system modeling, visualization of nuclear data, and convenient access to desired results.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Compositional inhomogeneities in a single Icelandic tholeiite flow

New trace element analyses of 25 basalt samples from a vertical traverse in an 11 m thick flow of Icelandic tholeiite show compositional differences which exceed analytical uncertainties and are random with respect to height. These variations suggest a short-range segregation model in which the amounts of phenocrysts, groundmass materials, and residual liquid vary randomly in different samples of the flow; a least-squares method is applied to determine if these variations reflect different proportions of crystallizing phases and residual liquid. Monte Carlo computations indicate that short-range segregation results in inhomogeneity of a magnitude of interflow compositional differences; this causes incoherence in elements with different geochemical behaviors while preserving coherence in elements of similar behavior.

Lindstrom, M. M.↗

Selecting Appropriate Model Complexity: An Example of Tracer Inversion for Thermal Prediction in Enhanced Geothermal Systems

Abstract A major challenge in the inversion of subsurface parameters is the ill‐posedness issue caused by the inherent subsurface complexities and the generally spatially sparse data. Appropriate simplifications of inversion models are thus necessary to make the inversion process tractable and meanwhile preserve the predictive ability of the inversion results. In this study, we investigate the effect of model complexity on fracture aperture inversion and thermal performance prediction in a field‐scale EGS model. Principal component analysis was used to map the aperture field to a low‐dimensional latent space. The complexity of the inversion model was quantitatively represented by the percentage of total variance in the original aperture fields preserved by the latent space. Tracer, pressure and flow rate data were used to invert for fracture aperture through an ensemble‐based inversion method, and the inferred aperture field was used to predict thermal performance. With an over‐simplified aperture model, ensemble collapse occurred. The inverted aperture models failed to resolve necessary flow and transport features, leading to a biased thermal performance prediction. A complex aperture model involved excessive features and was prone to overinterpreting the inversion data. Both the tracer/pressure/flow rate data reproduction and thermal prediction showed significant uncertainties, making it difficult to properly estimate long‐term thermal performance. Fortunately, our results indicate that there exists an appropriate model complexity which can simultaneously match inversion data and predict thermal performance with an acceptable uncertainty. The quality of the fit of tracer data appears to be a useful indicator of such an appropriate model complexity.

15 GEOTHERMAL ENERGY↗

Desmearing Bonse–Hart USANS data using Bayesian Gaussian process regression

Ultra-small-angle neutron scattering (USANS) enables access to micrometer-scale structures but is intrinsically affected by strong, anisotropic resolution smearing arising from slit-geometry optics. As a result, recovery of the intrinsic scattering intensity constitutes an ill-posed inverse problem, and commonly used iterative desmearing methods lack rigorous uncertainty quantification. We present a Bayesian desmearing framework for slit-geometry USANS based on Gaussian process regression. In this approach, the scattering intensity is modeled as a smooth random function, and the instrumental point spread function is incorporated explicitly as a forward operator. The resulting formulation yields a closed-form maximum a posteriori solution with well-defined credibility intervals. Computational benchmarks and experimental validation using combined USANS and small-angle neutron scattering (SANS) measurements demonstrate that the framework enables stable desmearing, suppresses experimental noise, and preserves physically meaningful structural features under realistic conditions.

Tung, Chi-Huan [Oak Ridge National Laboratory (ORN↗

Extending Conformal Prediction to Hidden Markov Models with Exact Validity via de Finetti’s Theorem for Markov Chains

Conformal prediction is a widely used method to quantify uncertainty in settings where the data is independent and identically distributed (IID), or more generally, exchangeable. Conformal prediction takes in a pre-trained classifier and a calibration dataset as inputs, and returns a function which maps feature vectors to subsets of classes. The output of the returned function for a new feature vector is guaranteed to contain the true class with a pre-specified confidence. Despite its success and usefulness in IID settings, extending conformal prediction to non-exchangeable (e.g., Markovian) data in a manner that provably preserves all desirable theoretical properties has largely remained an open problem. As a solution, we extend conformal prediction to the setting of a Hidden Markov Model (HMM) with unknown parameters. The key idea behind the proposed method is to partition the non-exchangeable Markovian data from the HMM into exchangeable blocks by exploiting the de Finetti’s Theorem for Markov Chains discovered by Diaconis and Freedman (1980). The permutations of the exchangeable blocks are then viewed as randomizations of the observed Markovian data from the HMM. The proposed method provably retains all desirable theoretical guarantees offered by the classical conformal prediction framework and is general enough to be useful in many sequential prediction problems.

Nettasinghe, Don Buddhika Wijayantha↗

MrHyDE v.1.0

SAND2024-01324O MrHyDE, which stands for Multi-resolution Hybridized Differential Equations, is a general-purpose C++ package for the solution of coupled multiphysics and multiscale systems on massively parallel computing systems. MrHyDE is designed to enable moving beyond forward simulation for multiscale applications which includes optimization, control, uncertainty quantification, and stochastic inversion. The framework provides interfaces to several packages within the Trilinos framework and leverages automatic differentiation to enable adjoint capabilities for large-scale, gradient-based optimization. MrHyDE provides automated multiscale capabilities through a subgrid model interface and multiscale Dirichlet-to-Neumann maps. For extreme-scale applications, MrHyDE provides in situ data-compression algorithms to reduce memory requirements while maintaining performance. MrHyDE is a general-purpose, computational framework for the solution of multiscale and multiphysics applications. It uses a combination of structure-preserving, physics-compatible discretizations, fully implicit methods, multi-resolution schemes, or fully explicit methods. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗

Development of SCALE Tools for Sensitivity and Uncertainty Analysis Methodology Implementation (TSUNAMI) from SCALE 5 through SCALE 6.2

This review paper provides an abbreviated historical review of the development of the TSUNAMI capabilities in SCALE, beginning with the first initiatives, to provide new insights into validation, uncertainty quantification, and safety assessment. The tools of TSUNAMI have evolved substantially over two decades to meet the growing application analysis and experiment design needs. Even more advancements are emerging as the next generation of innovators continue to use and develop these tools.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

High-Efficiency High-Resolution Global Model Developments at the NASA Goddard Data Assimilation Office

The Data Assimilation Office (DAO) has been developing a new generation of ultra-high resolution General Circulation Model (GCM) that is suitable for 4-D data assimilation, numerical weather predictions, and climate simulations. These three applications have conflicting requirements. For 4-D data assimilation and weather predictions, it is highly desirable to run the model at the highest possible spatial resolution (e.g., 55 km or finer) so as to be able to resolve and predict socially and economically important weather phenomena such as tropical cyclones, hurricanes, and severe winter storms. For climate change applications, the model simulations need to be carried out for decades, if not centuries. To reduce uncertainty in climate change assessments, the next generation model would also need to be run at a fine enough spatial resolution that can at least marginally simulate the effects of intense tropical cyclones. Scientific problems (e.g., parameterization of subgrid scale moist processes) aside, all three areas of application require the model's computational performance to be dramatically improved as compared to the previous generation. In this talk, I will present the current and future developments of the "finite-volume dynamical core" at the Data Assimilation Office. This dynamical core applies modem monotonicity preserving algorithms and is genuinely conservative by construction, not by an ad hoc fixer. The "discretization" of the conservation laws is purely local, which is clearly advantageous for resolving sharp gradient flow features. In addition, the local nature of the finite-volume discretization also has a significant advantage on distributed memory parallel computers. Together with a unique vertically Lagrangian control volume discretization that essentially reduces the dimension of the computational problem from three to two, the finite-volume dynamical core is very efficient, particularly at high resolutions. I will also present the computational design of the dynamical core using a hybrid distributed-shared memory programming paradigm that is portable to virtually any of today's high-end parallel super-computing clusters.

Lin, Shian-Jiann↗