Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical core”

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 37 records · Page 2

An Intercomparison of the Dynamical Cores of Global Atmospheric Circulation Models for Mars

This is a Final Report for a Joint Research Interchange (JRI) between NASA Ames Research Cen- ter and San Jose State University, Department of Meteorology. The focus of this JRI has been to evaluate the dynamical "cores" of two global atmospheric circulation models for Mars that are in operation at the NASA Ames Research Center. ne two global circulation models in use are fundamentally different: one uses spherical harmonics in its horizontal representation of field variables; the other uses finite differences on a uniform longitude-latitude grid. Several simulations have been conducted to assess how the dynamical processors of each of these circulation models perform using identical "simple physics" parameterizations. A variety of climate statistics (e.g., time-mean flows and eddy fields) have been compared for realistic solstitial mean basic states. Results of this research have demonstrated that the two Mars circulation models with completely different spatial representations and discretizations produce rather similar circulation statistics for first-order meteorological fields, suggestive of a tendency for convergence of numerical solutions. Second and higher-order fields can, however, vary significantly between the two models.

Hollingsworth, Jeffery L.↗

Development of a River Dynamical Core for E3SM to simulate compound flooding on Exascale-class heterogeneous supercomputers

Flooding events pose significant risk to human life, property, and infrastructure. Physically-consistent quantification of altered flood risks in global models requires hyper-resolution (~1 km) or fine flood simulations using two-dimensional (2D) physics schemes, both of which are unavailable in the current generation Earth System Models. Here, in this work, we have developed the River Dynamical Core (RDycore), which is an open-source, 2D shallow water equation (SWE) library for the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM). RDycore uses PETSc and libCEED libraries that allows it to run efficiently on CPUs and GPUs, as well as select a time-integration algorithm at runtime without requiring any code modifications. RDycore achieves spatial error convergence rates for problems with analytical and manufactured solutions similar to those reported previously in the literature, or consistent with the implemented first-order spatial discretization scheme. RDycore's accuracy in predicting flooding for a well-studied dam break problem is comparable to existing SWE models. For a problem with 471 million grid cells, RDycore achieves a speedup of 6.6x and 7.6x on GPUs compared to CPUs when using 320 compute nodes on DOE's Perlmutter and Frontier supercomputers, respectively. The one-way coupling of the RDycore library within E3SM is demonstrated by performing multiple 5-day flooding simulations during Hurricane Harvey driven by five precipitation datasets. The E3SM--RDycore simulations at 30 m spatial resolution accurately simulate maximum water height during the hurricane when benchmarked against a previously published study and achieve a speedup of 15x (Perlmutter) and 21x (Frontier) on GPUs relative to CPUs. The work presented here is the foundational step in providing hardware and algorithmic portability framework for simulating kilometer-scale river dynamics within E3SM.

Flood Simulation↗

Progress Towards Integrating the Finite-Volume Cubed-Sphere (FV3) Dynamical Core Tangent Linear and Adjoint Models into JEDI

The Joint Effort for Data assimilation Integration (JEDI) -- led by the Joint Center for Satellite Data Assimilation (JCSDA) -- is an inter-organizational endeavor to develop a common framework for performing data assimilation. This extensive framework will ultimately provide solvers, observation operators, interpolation and model interfaces using object oriented modeling. Two partners involved in JEDI use or plan to use the Finite Volume Cubed-Sphere (FV3) dynamical core to produce weather forecasts; these are NASA's Global Modeling and Assimilation Office and NOAA's National Center for Environment Prediction. In this work we present an update on ongoing efforts to integrate the FV3 tangent linear and adjoint models into the prototype JEDI framework. We setup and run a simple cycled data assimilation experiment using 4DVAR on the cubed sphere grid and with the FV3 tangent linear and adjoint models. Development of the observation operators for JEDI is separately underway. Instead of using real observations a simplified set of simulated observations will be used. We discuss the steps required to bring the FV3 linearized model into the object oriented framework and consider what would be the computational requirements of running this configuration for an operational system. FV3 uses a small time-step to ensure that small scales are well resolved, however this presents design challenges when running 4DVAR with the adjoint. An approach to storing the FV3 model trajectory has been developed that maintains the flexibility of using automatic differentiation. We discuss how this approach is incorporated into the framework. Other important uses of adjoint models include computing observation impacts and singular vectors, we consider how these tools can be included in JEDI.

Holdaway, Daniel↗

Forced nutations of the earth: Influence of inner core dynamics. I - Theory. II - Numerical results and comparisons. III - Very long interferometry data analysis

A treatment is presented of the nutation problem for an oceanless, elastic, spheroidally stratified earth, with the dynamical role of the inner core explicitly included in the formulation. Solving the enlarged system of equations shows that a new almost diurnal eigenfrequency emerges. A rough estimate places it not far from the prograde annual tidal excitation frequency, so that possible resonance effects on nutation amplitudes need careful consideration. Tables are provided that exhibit the sensitivities of various relevant quantities, the eigenfrequencies and the coefficients which appear in the resonance expansion, as well as the nutation amplitudes at important tidal frequencies, to possible errors in the earth parameters which enter the theory set forth. Finally, the analysis of 798 VLBI experiments performed between July 1980 and February 1989 and the determination from this analysis of corrections to selected coefficients in the International Astronomical Union 1980 theory of the nutations of the earth are discussed.

Mathews, P. M.↗

Synchrotron radiography characterization of the liquid core dynamics in a canonical two-fluid coaxial atomizer

The liquid core of a canonical two-fluid coaxial atomizer has been characterized using synchrotron X-rays. The high energy photons allow for high-speed imaging of attenuation through the dense liquid-gas jet core, resolving the internal structures that include entrapped air bubbles and the formation of liquid ligaments and bags. When the gas-to-liquid momentum ratio increases, the liquid core transitions from an intact column, where primary break-up happens several liquid diameters downstream, to a hollow crown with a downstream span comparable to the liquid diameter that disintegrates by shedding ligaments from its rim. At high gas momentum ratios (limited by the sonic velocity at the gas nozzle exit), this crown suffers partial dewetting and, when angular momentum is added to the gas, it dewets on a large section of the liquid injection needle circumference. This partial crown exhibits azimuthal motions along the circumference, on timescales much longer than the relevant flow timescales. Here, the crown attachment to the liquid needle presents a bi-stable nature. The dramatic changes of the liquid core morphology, as the gas momentum and swirl ratios vary, have a strong impact on the gas-liquid boundary layers, which control the liquid break-up mechanisms and the resulting spray characteristics, such as droplet size distributions and the droplet volume fraction across the spray.

42 ENGINEERING↗

The source of the Earth's long wavelength geoid anomalies: Implications for mantle and core dynamics

The long wavelength components of the Earth's gravity field result mainly from density contrasts associated with convection in the mantle. Direct interpretation of the geoid for mantle convection is complicated by the fact that convective flow results in dynamically maintained deformation of the surface of the Earth, the core mantle boundary (CMB), and any interior chemical boundaries which might exist. These boundary deformations effect the geoid opposite in sign and are comparable in magnitude to those of the interior density contrasts driving the flow. The total difference of two relatively large quantities.

Hager, B. H.↗

MISR CMVs and Multiangular Views of Tropical Cyclone Inner-Core Dynamics

Multi-camera stereo imaging of cloud features from the MISR (Multiangle Imaging SpectroRadiometer) instrument on NASA's Terra satellite provides accurate and precise measurements of cloud top heights (CTH) and cloud motion vector (CMV) winds. MISR observes each cloudy scene from nine viewing angles (Nadir, +/-26(sup o), +/-46(sup o), +/-60(sup o), +/-70(sup o)) with approximatel 275-m pixel resolution. This paper provides an update on MISR CMV and CTH algorithm improvements, and explores a high-resolution retrieval of tangential winds inside the eyewall of tropical cyclones (TC). The MISR CMV and CTH retrievals from the updated algorithm are significantly improved in terms of spatial coverage and systematic errors. A new product, the 1.1-km cross-track wind, provides high accuracy and precision in measuring convective outflows. Preliminary results obtained from the 1.1-km tangential wind retrieval inside the TC eyewall show that the inner-core rotation is often faster near the eyewall, and this faster rotation appears to be related linearly to cyclone intensity.

hurricanes↗

Physically-Derived Dynamical Cores in Atmospheric General Circulation Models

The algorithm chosen to represent the advection in atmospheric models is often used as the primary attribute to classify the model. Meteorological models are generally classified as spectral or grid point, with the term grid point implying discretization using finite differences. These traditional approaches have a number of shortcomings that render them non-physical. That is, they provide approximate solutions to the conservation equations that do not obey the fundamental laws of physics. The most commonly discussed shortcomings are overshoots and undershoots which manifest themselves most overtly in the constituent continuity equation. For this reason many climate models have special algorithms to model water vapor advection. This talk focuses on the development of an atmospheric general circulation model which uses a consistent physically-based advection algorithm in all aspects of the model formulation. The shallow-water model is generalized to three dimensions and combined with the physics parameterizations of NCAR's Community Climate Model. The scientific motivation for the development is to increase the integrity of the underlying fluid dynamics so that the physics terms can be more effectively isolated, examined, and improved. The expected benefits of the new model are discussed and results from the initial integrations will be presented.

Rood, Richard B.↗

Physically-Derived Dynamical Cores in Atmospheric General Circulation Models

The algorithm chosen to represent the advection in atmospheric models is often used as the primary attribute to classify the model. Meteorological models are generally classified as spectral or grid point, with the term grid point implying discretization using finite differences. These traditional approaches have a number of shortcomings that render them non-physical. That is, they provide approximate solutions to the conservation equations that do not obey the fundamental laws of physics. The most commonly discussed shortcomings are overshoots and undershoots which manifest themselves most overtly in the constituent continuity equation. For this reason many climate models have special algorithms to model water vapor advection. This talk focuses on the development of an atmospheric general circulation model which uses a consistent physically-based advection algorithm in all aspects of the model formulation. The shallow-water model of Lin and Rood (QJRMS, 1997) is generalized to three dimensions and combined with the physics parameterizations of NCAR's Community Climate Model. The scientific motivation for the development is to increase the integrity of the underlying fluid dynamics so that the physics terms can be more effectively isolated, examined, and improved. The expected benefits of the new model are discussed and results from the initial integrations will be presented.

Rood, Richard B.↗

Dynamical Core Damping of Thermal Tides in Martian Atmosphere

Atmospheric oscillations with daily periodicity are observed in in-situ near-surface pressure, temperature, and winds observations, and also in remotely sensed temperature and pressure observations of the Martian atmosphere. Such oscillations are interpreted as thermal tides driven by the diurnal cycle of solar radiation and occur at various frequencies, with the most prominent being the diurnal, semidiurnal, terdiurnal and quadiurnal tides. Mars global circulation models reproduce these tides with varying levels of success. Until recently, both the MarsWRF and newly developed MarsMPAS models were able to produce realistic diurnal and semidiurnal tide amplitudes but predicted higher-order mode amplitudes that were significantly weaker than observed. We use linear wave analysis to show that the divergence damping applied within both MarsWRF and MarsMPAS is responsible for suppressing the amplitude of thermal tides with frequency greater than 2 per sol, despite being designed to suppress only acoustic wave modes. Decreasing the strength of the divergence damping in MarsWRF and MarsMPAS allows for excellent prediction of the higher order tidal modes. This finding demonstrates that care must be taken when applying numerical dampers and filters that may eliminate some desired dynamical features in planetary atmospheres.

Yuan Lian↗

Compact video content representation for video coding using low multi-linear tensor rank approximation with dynamic core tensor order

Tensor is a multi-way representation of an entity that is widely used in many applications of mathematics and physics. The wide range of tensor algebra operations attract the researchers of graphical and vision field to identify unique methods for various applications such as segmentation, compression, indexing and so on. In this paper, low multi-linear rank approximation of a tensor is used for decomposition so as to perform compression on a video that is represented as a 4D tensor. Thus, the encoder has blocks for (1) multilinear rank approximation method so as to aid decomposition, (2) sparsity removal of residual data to reduce the required memory storage to a small level, (3) quantization on the sparse less residual data and (4) LZ77 dictionary coding for the video. The decoder uses a tensor reconstruction block along with residual error correction block that reduces loss/error in the reconstructed video. The size of core tensor is an important factor in rank approximation. It is identified at run time depending on the video contents. This size is determined adaptively using Tikhonov’s regularization method. The best value of core tensor size is identified with the corner of L-curve which preserves the video contents from heavy loss/error. Experimental results of the proposed work on comparison with H.264/AVC compression method indicate that the proposed method achieves up to 10% improvement in compression ratio with the least amount of loss and a maximum signal to noise ratio of 0.14. The proposed system gives better compression ratio than H.265 for videos with more number of frames without compromising quality.

97 MATHEMATICS AND COMPUTING↗

Synergistic and Collaborative Development Strategies for FV3 Powered Next Generation Unified Global Modeling System

The GFDL (Geophysical Fluid Dynamics Laboratory, NOAA - National Oceanic and Atmospheric Administration) Finite-­Volume Cubed-Sphere Dynamical Core (FV3) is a scalable and flexible dynamical core capable of both hydrostatic and non-hydrostatic atmospheric simulations. FV3 has been chosen as the dynamical core for the Next Generation Global Prediction System project (NGGPS), designed to upgrade the current NCEP (National Centers for Environmental Prediction, NOAA) operational Global Forecast System (GFS) to run as a unified, fully-coupled system in NOAA's Environmental Modeling System infrastructure. FV3 dynamic core has a long history of serving as the main engine for global atmospheric models at various Government and Academic Research Laboratories including NOAA: GFDL Climate Modeling Suite (AM4 (Atmosphere Modeling 4), CM4 (Climate Model 4) ESM4 (Earth System Model 4), Hiram (HIgh Resolution Atmospheric Model)); NASA: GMAO (Global Modeling and Assimilation Office) Goddard Earth Observing System Model (GEOS); and NCAR (National Center for Atmospheric Research) Community Earth System Model (CESM). The three primary stakeholders in FV3 (GFDL, GMAO, EMC (Environmental Modeling Center - NOAA)) have embarked on synergistic and collaborative strategies with focus on the advancement of non-hydrostratic dynamic core, physics, chemistry, and data assimilation efforts, leveraging collective strengths of each of the partnering agencies. This talk presents the ongoing plans for model advancements at GFDL, GMAO and EMC, with special emphasis on contributions towards developing community based unified global modeling system for operational and research applications at respective organizations. Future plans include extending the collaborations to the development of earth system components including GMAO's aerosol chemistry models (GOCART/MAM (Goddard Chemistry Aerosol Radiation and Transport / Modal Aerosol Model)), Land Information System (LIS), and advanced data assimilation techniques; and GFDL's Modular Ocean Model (MOM) and Sea Ice Simulator (SIS) for transition to operations at NCEP.

Putman, William↗

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↗

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 kin 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↗