Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation model”

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 397 records · Page 22

A new method to construct a horizontal resolution-dependent wind speed adjustment factor for tropical cyclones in global climate model simulations

A new method to construct a horizontal resolution-dependent wind speed adjustment factor for evaluating tropical cyclones (TCs) in global climate models (GCMs) is presented here. In contrast to the previous studies that used idealized axisymmetric wind fields, this study analyzes 48 hours of 10-second surface wind fields from 1-km TC simulations. The adjustment factor is derived from filtering the simulated TC wind fields onto various horizontal grid spacings typical of those used in GCMs. The new adjustment factor leads to TCs with greater intensity than the existing adjustment factors for horizontal grid spacings smaller than 30 km. This difference is attributed to more realistic wind fields in the TC simulations that contain highly asymmetric, localized patches of higher wind speeds instead of axisymmetric wind fields. Applying the new adjustment factor to select GCM simulations suggests the common interpretation of low-intensity bias in GCM-simulated TCs might be slightly exaggerated.

54 ENVIRONMENTAL SCIENCES↗

Variable Lidar Ratios for Marine and Dusty Marine Using MODIS AOD Constrained Retrievals and GOCART Model Simulations and their Impact on Level 2 CALIPSO Data Products in Version 5

The Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) onboard CALIPSO provided global measurements of attenuated backscatter profiles of various tropospheric aerosol species from June 2006 through June 2023. Extinction profiles are retrieved from these backscatter measurements by assuming a value for the extinction-to-backscatter or lidar ratio (LR). Up until version V4.51 of the data products, only a constant global value for each of the species has been used, spatially and temporally. However, it has been well known from various ground-based measurements that LRs of each aerosol species can potentially vary spatially as well seasonally. For the forthcoming Version 5 of the CALIPSO products, the CALIPSO team will implement variable lidar ratios at 532 nm for the species classified as “marine” and “dusty marine”. Here we present the essential elements of this scheme, which involves estimating lidar ratios by constraining Fernald solutions for the CALIOP particulate backscatter profiles using aerosol optical depths (AOD) from collocated MODIS retrievals. To account for the lack of measurements in various regions and seasons (e.g., Arctic winters), we leverage collocated model sea-salt volume fractions (SSVF) simulated by the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model and use an empirical relationship between the SSVF and LR to build hybrid climatological maps for the entire globe, i.e., using the modelled values when reliable values of the constrained LRs are not available. Initial evaluation of the resulting extinction profiles and corresponding AODs, both globally and in specific regions are presented using a limited amount of data. The primary improvement in AOD occurs in the coastal areas, like the Arabian Sea and Bay of Bengal. The estimated marine LR values in these regions are significantly higher than the constant value of 23 sr used in previous data releases, likely due to mixing with offshore pollution which results in higher AODs. Initial validation results using the CALIPSO Ocean Derived Column Optical Depth (ODCOD) retrievals and ground-based AERONET AOD measurements are presented.

CALIPSO↗

Model simulations of rainout and washout from a warm stratiform cloud

Calculations utilizing a one-dimensional, time-dependent stratiform cloud model with parameterized microphysics to study the washout and rainout of soluble trace gases are presented. It is shown that the in-cloud concentrations of relatively insoluble species remain essentially unaffected by the presence of cloudwater and rainwater, while the concentrations of highly soluble species are rapidly depleted within the cloud. Due to this depletion it is shown that for highly soluble species the below-cloud scavenging via washout can dominate over in-cloud scavenging via rainout; in-cloud, aqueous-phase cloudwater and rainwater concentrations can differ considerably and, near the cloudbase, cloudwater concentrations can exceed rainwater concentrations by factors of two to ten.

Xing, Lin↗

Overview of Proposal on High Resolution Climate Model Simulations of Recent Hurricane and Typhoon Activity: The Impact of SSTs and the Madden Julian Oscillation

This talk gives an update on the progress and further plans for a coordinated project to carry out and analyze high-resolution simulations of tropical storm activity with a number of state-of-the-art global climate models. Issues addressed include, the mechanisms by which SSTs control tropical storm. activity on inter-annual and longer time scales, the modulation of that activity by the Madden Julian Oscillation on sub-seasonal time scales, as well as the sensitivity of the results to model formulation. The project also encourages companion coarser resolution runs to help assess resolution dependence, and. the ability of the models to capture the large-scale and long-terra changes in the parameters important for hurricane development. Addressing the above science questions is critical to understanding the nature of the variability of the Asian-Australian monsoon and its regional impacts, and thus CLIVAR RAMP fully endorses the proposed tropical storm simulation activity. The project is open to all interested organizations and investigators, and the results from the runs will be shared among the participants, as well as made available to the broader scientific community for analysis.

Schubert, Siegfried↗

Observations and wind-forced model simulations of the mean seasonal cycle in tropical Pacific sea surface topography

This paper quantitatively evaluates simulations of the mean seasonal cycle in the tropical Pacific Ocean, using a multiple vertical-mode version of the Busalacchi and O'Brien (1980) linear model. The period studied (1979-1981) includes several coincident oceanic and surface wind data sets, as well as a fairly regular seasonal cycle leading into the 1982-1983 El Nino. Simulated mean seasonal cycles in dynamic height and sea level are compared with observed variations in expendable bathythermograph and island tide gauge data averaged over the same period. It was found that, although all simulations show characteristic features of the mean meridional ridge-trough structure in surface topography, simulated north and south equatorial ridges at 20 deg N and 20 deg S are much higher than those observed. Only weak equatorial ridges were generated near 4 N, and none of the simulations exhibited a significant equatorial trough. These discrepancies are attributed to limitations in model physics and to the wind forcing.

Mcphaden, Michael J.↗

The Impacts of Microphysics and Planetary Boundary Layer Physics on Model Simulations of U.S. Deep South Summer Convection

Convection-allowing numerical weather simula- tions have often been shown to produce convective storms that have significant sensitivity to choices of model physical parameterizations. Among the most important of these sensitivities are those related to cloud microphysics, but planetary boundary layer parameterizations also have a significant impact on the evolution of the convection. Aspects of the simulated convection that display sensitivity to these physics schemes include updraft size and intensity, simulated radar reflectivity, timing and placement of storm initi- ation and decay, total storm rainfall, and other storm features derived from storm structure and hydrometeor fields, such as predicted lightning flash rates. In addition to the basic parameters listed above, the simulated storms may also exhibit sensitivity to im- posed initial conditions, such as the fields of soil temper- ature and moisture, vegetation cover and health, and sea and lake water surface temperatures. Some of these sensitivities may rival those of the basic physics sensi- tivities mentioned earlier. These sensitivities have the potential to disrupt the accuracy of short-term forecast simulations of convective storms, and thereby pose sig- nificant difficulties for weather forecasters. To make a systematic study of the quantitative impacts of each of these sensitivities, a matrix of simulations has been performed using all combinations of eight separate microphysics schemes, three boundary layer schemes, and two sets of initial conditions. The first version of initial conditions consists of the default data from large-scale operational model fields, while the second features specialized higher- resolution soil conditions, vegetation conditions and water surface temperatures derived from datasets created at NASA's Short-term Prediction and Operational Research Tran- sition (SPoRT) Center at the National Space Science and Technology Center (NSSTC) in Huntsville, AL. Simulations as outlined above, each 48 in number, were conducted for five midsummer weakly sheared coastal convective events each at two sites, Mobile, AL (MOB) and Houston, TX (HGX). Of special interest to operational forecasters at MOB and HGX were accuracy of timing and placement of convective storm initiation, reflectivity magnitudes and coverage, rainfall and inferred lightning threat.

McCaul, Eugene W., Jr.↗

Approximate simulation model for analysis and optimization in engineering system design

Computational support of the engineering design process routinely requires mathematical models of behavior to inform designers of the system response to external stimuli. However, designers also need to know the effect of the changes in design variable values on the system behavior. For large engineering systems, the conventional way of evaluating these effects by repetitive simulation of behavior for perturbed variables is impractical because of excessive cost and inadequate accuracy. An alternative is described based on recently developed system sensitivity analysis that is combined with extrapolation to form a model of design. This design model is complementary to the model of behavior and capable of direct simulation of the effects of design variable changes.

Sobieszczanski-Sobieski, Jaroslaw↗

Does Dynamical Downscaling Introduce Novel Information in Climate Model Simulations of Recipitation Change over a Complex Topography Region?

Current climate and future climate-warming runs with the RegCM Regional Climate Model (RCM) at 50 and 11 km-resolutions forced by the ECHAM GCM are used to examine whether the increased resolution of the RCM introduces novel information in the precipitation field when the models are run for the mountainous region of the Hellenic peninsula. The model results are inter-compared with the resolution of the RCM output degraded to match that of the GCM, and it is found that in both the present and future climate runs the regional models produce more precipitation than the forcing GCM. At the same time, the RCM runs produce increases in precipitation with climate warming even though they are forced with a GCM that shows no precipitation change in the region. The additional precipitation is mostly concentrated over the mountain ranges, where orographic precipitation formation is expected to be a dominant mechanism. It is found that, when examined at the same resolution, the elevation heights of the GCM are lower than those of the averaged RCM in the areas of the main mountain ranges. It is also found that the majority of the difference in precipitation between the RCM and the GCM can be explained by their difference in topographic height. The study results indicate that, in complex topography regions, GCM predictions of precipitation change with climate warming may be dry biased due to the GCM smoothing of the regional topography.

Tselioudis, George↗

The International Heat Stress Genotype Experiment for Modeling Wheat Response to Heat: Field Experiments and AgMIP-Wheat Multi-Model Simulations

The data set contains a portion of the International Heat Stress Genotype Experiment (IHSGE) data used in the AgMIP-Wheat project to analyze the uncertainty of 30 wheat crop models and quantify the impact of heat on global wheat yield productivity. It includes two spring wheat cultivars grown during two consecutive winter cropping cycles at hot, irrigated, and low latitude sites in Mexico (Ciudad Obregon and Tlaltizapan), Egypt (Aswan), India (Dharwar), the Sudan (Wad Medani), and Bangladesh (Dinajpur). Experiments in Mexico included normal (November-December) and late (January-March) sowing dates. Data include local daily weather data, soil characteristics and initial soil conditions, crop measurements (anthesis and maturity dates, anthesis and final total above ground biomass, final grain yields and yields components), and cultivar information. Simulations include both daily in-season and end-of-season results from 30 wheat models.

Wheat; Field Experimental Data; Heat Stress; Simul↗

How can biosphere models simulate enough vegetation biomass in the mountains of the western United States? Implications of meteorological forcing

Most carbon stocks and fluxes in the western United States are found in mountainous terrain, where observations and modeling are difficult. Terrestrial biosphere models generally underestimate above-ground biomass (AGB) over this region. In this report, we identify methods to reduce this underestimation by focusing upon 1) biases in meteorological datasets, 2) model representation of water stress, and 3) spatial resolution. We adopted the widely-used Community Land Model version 4.5 (CLM 4.5) with six different meteorological datasets and found a 6-fold variation in simulated AGB across Utah/Colorado. Simulations underestimated AGB because of warm and dry biases within the meteorological datasets that reduced water availability and restricted plant growth. To eliminate the AGB underestimation we adopted a meteorological dataset designed for complex terrain (gridMET), combined with a representation of plant hydraulic stress (CLM 5.0). Conversely, changes in spatial resolution (meteorological variables and land surface description) had negligible impact on simulated AGB.

54 ENVIRONMENTAL SCIENCES↗

Interactive Mesh Generation and Comparison Tools for Nuclear Energy Advanced Modeling & Simulation Phase I. Final Report

The promise of computer simulations for nuclear reactor designs is the opportunity to speed up the design process while also improving safety, providing more detail and higher fidelity, reducing cost, minimizing schedule risk, and avoiding many other potential problems; because simulations have the potential to be much less expensive and time-consuming than traditional development – where scale models and physical tests require construction and long lead times – it is possible to consider more alternatives and complete more comprehensive assessments before the design is finalized and constructed. As more design scenarios are studied by simulation, the speed that simulations promise is only available as long as bottlenecks are addressed for these scenarios. A frequent bottleneck is the amount of human labor required to describe the shape of the objects being designed and decompose the shape into simple elements, such as hexahedra, that are small enough to capture the physical phenomena of interest without being so small that even fast computers cannot perform the simulation in a reasonable amount of time. Recently, techniques for creating all-hexahedral decompositions have been developed, but they require some human input. Many simulation codes prefer all-hexahedral decompositions, so we proposed to evaluate the feasibility of these techniques on nuclear reactor geometries to see whether they required too much human input to be commercially viable. Our study concluded that they can be made viable with some additional software tools to reduce the amount of user input required.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Simulation Modeling and Performance Evaluation of Space Networks

In space exploration missions, the coordinated use of spacecraft as communication relays increases the efficiency of the endeavors. To conduct trade-off studies of the performance and resource usage of different communication protocols and network designs, JPL designed a comprehensive extendable tool, the Multi-mission Advanced Communications Hybrid Environment for Test and Evaluation (MACHETE). The design and development of MACHETE began in 2000 and is constantly evolving. Currently, MACHETE contains Consultative Committee for Space Data Systems (CCSDS) protocol standards such as Proximity-1, Advanced Orbiting Systems (AOS), Packet Telemetry/Telecommand, Space Communications Protocol Specification (SCPS), and the CCSDS File Delivery Protocol (CFDP). MACHETE uses the Aerospace Corporation s Satellite Orbital Analysis Program (SOAP) to generate the orbital geometry information and contact opportunities. Matlab scripts provide the link characteristics. At the core of MACHETE is a discrete event simulator, QualNet. Delay Tolerant Networking (DTN) is an end-to-end architecture providing communication in and/or through highly stressed networking environments. Stressed networking environments include those with intermittent connectivity, large and/or variable delays, and high bit error rates. To provide its services, the DTN protocols reside at the application layer of the constituent internets, forming a store-and-forward overlay network. The key capabilities of the bundling protocols include custody-based reliability, ability to cope with intermittent connectivity, ability to take advantage of scheduled and opportunistic connectivity, and late binding of names to addresses. In this presentation, we report on the addition of MACHETE models needed to support DTN, namely: the Bundle Protocol (BP) model. To illustrate the use of MACHETE with the additional DTN model, we provide an example simulation to benchmark its performance. We demonstrate the use of the DTN protocol and discuss statistics gathered concerning the total time needed to simulate numerous bundle transmissions

network protocols↗

Thermospheric dynamics during September 18-19, 1984. I - Model simulations

The September 18-19, 1984 Equinox Transient Study (ETS) interval was simulated using the National Center for Atmospheric Research thermospheric GCM with inputs guided by observations. High-latitude forcings for the model were defined with precision for this period. The importance of upward-propagating semidiurnal tides is emphasized, the first-realistic time-dependent simulation to incorporate such tides. A description of the simulated quiet time thermospheric structure for September 18 is contrasted with the storm time behavior of September 19. Several features of the storm response are described, including the generation of large-scale equatorward propagating disturbances. Longer-lived perturbations of the neutral temperature and composition are also discussed. A new feature discovered in the model predictions for altitudes around 200 km comprises a quasi-fixed, four-cell pattern of high-and low-density regions in the polar cap.

Crowley, G.↗

Enabling kilometer-scale E3SM land model simulation over North America: A new integrated framework solution

This study introduces a novel framework designed to enhance the performance, scalability, and portability of the kilometer-scale E3SM Land Model (km-ELM) within the E3SM modeling infrastructure. By seamlessly integrating cutting-edge data tools, we address existing challenges such as slow performance, limited scalability, and difficulties in software integration in current data-driven ELM simulation over large geographic areas. Our innovative approach leverages the KiloCraft data toolkit to generate unified inputs for simulations ranging from a single-cite case, to a 72,083-cell regional case to a continental configuration encompassing 21.6 million land grid cells at a 1 km × 1 km resolution. We conduct extensive strong- and weak-scaling experiments on three state-of-the-art supercomputers, utilizing up to 100,800 CPU cores across 2400 compute nodes to evaluate end-to-end metrics including wall-clock time, simulation-years-per-day (SYPD), initialization costs, and I/O throughput. Our results reveal the land (LND) component’s efficient scaling, demonstrating near-ideal weak scaling and strong-scaling parallel efficiencies reaching up to 87% at 50,400 cores. We confirm portability and reproducibility through bitwise-equivalent outputs across different machines using identical inputs over supported machines. Notably, at extreme scales, we identify I/O as a critical bottleneck and that leads to effective solution with the SCORPIO/ADIOS stack. Collectively, these findings validate the deployment of km-ELM at a continental scale with high parallel efficiency and provide essential guidance on configuration, decomposition, and I/O settings for optimized kilometer-scale land simulations in E3SM. This work emphasizes the innovative design and practical solutions that enhance the operational capabilities of km-ELM, focusing on software performance and scalability while leaving detailed scientific evaluations of simulated land processes for future investigations.

E3SM land model (ELM), km-ELM, scalability, perfor↗

Spatial Variability in Column CO2 Inferred from High Resolution GEOS-5 Global Model Simulations: Implications for Remote Sensing and Inversions

Column CO2 observations from current and future remote sensing missions represent a major advancement in our understanding of the carbon cycle and are expected to help constrain source and sink distributions. However, data assimilation and inversion methods are challenged by the difference in scale of models and observations. OCO-2 footprints represent an area of several square kilometers while NASA s future ASCENDS lidar mission is likely to have an even smaller footprint. In contrast, the resolution of models used in global inversions are typically hundreds of kilometers wide and often cover areas that include combinations of land, ocean and coastal areas and areas of significant topographic, land cover, and population density variations. To improve understanding of scales of atmospheric CO2 variability and representativeness of satellite observations, we will present results from a global, 10-km simulation of meteorology and atmospheric CO2 distributions performed using NASA s GEOS-5 general circulation model. This resolution, typical of mesoscale atmospheric models, represents an order of magnitude increase in resolution over typical global simulations of atmospheric composition allowing new insight into small scale CO2 variations across a wide range of surface flux and meteorological conditions. The simulation includes high resolution flux datasets provided by NASA s Carbon Monitoring System Flux Pilot Project at half degree resolution that have been down-scaled to 10-km using remote sensing datasets. Probability distribution functions are calculated over larger areas more typical of global models (100-400 km) to characterize subgrid-scale variability in these models. Particular emphasis is placed on coastal regions and regions containing megacities and fires to evaluate the ability of coarse resolution models to represent these small scale features. Additionally, model output are sampled using averaging kernels characteristic of OCO-2 and ASCENDS measurement concepts to create realistic pseudo-datasets. Pseudo-data are averaged over coarse model grid cell areas to better understand the ability of measurements to characterize CO2 distributions and spatial gradients on both short (daily to weekly) and long (monthly to seasonal) time scales

Ott, L.↗

Global Temperature Responses to Large Tropical Volcanic Eruptions in Paleo Data Assimilation Products and Climate Model Simulations Over the Last Millennium

Abstract Large volcanic eruptions are one of the dominant perturbations to global and regional atmospheric temperatures on timescales of years to decades. Discrepancies remain, however, in the estimated magnitude and persistence of the surface temperature cooling caused by volcanic eruptions, as characterized by paleoclimatic proxies and climate models. We investigate these discrepancies in the context of large tropical eruptions over the Last Millennium using two state‐of‐the‐art data assimilation products, the Paleo Hydrodynamics Data Assimilation product (PHYDA) and the Last Millennium Reanalysis (LMR), and simulations from the National Center for Atmospheric Research Community Earth System Model‐Last Millennium Ensemble (NCAR CESM‐LME). We find that PHYDA and LMR estimate mean global and hemispheric cooling that is similar in magnitude and persistence once effects from eruptions occurring in short succession are removed. The estimates also compare well to Northern‐Hemisphere reconstructions based solely or partially on tree‐ring density, which have been proposed as the most accurate proxy estimates of surface cooling due to volcanism. All proxy‐based estimates also agree well with the magnitude of the mean cooling simulated by the CESM‐LME. Differences remain, however, in the spatial patterns of the temperature responses in the PHYDA, LMR, and the CESM‐LME. The duration of cooling anomalies also persists for several years longer in the PHYDA and LMR relative to the CESM‐LME. Our results demonstrate progress in resolving discrepancies between proxy‐ and model‐based estimates of temperature responses to volcanism, but also indicate these estimates must be further reconciled to better characterize the risks of future volcanic eruptions.

Tejedor, E.↗

Diagnostic evaluation of numerical model simulations using the tendency equation

A procedure is proposed to expand the diagnostic capabilities of the pressure tendency equation of a primitive equation NWP model by computing the pressure tendency in physical coordinates. The advantages of isolating the density advection as a diagnostic tool to understand pressure changes is shown. By simple thermodynamic arguments it is demonstrated that in areas of synoptic-scale cyclonic development, the vertically integrated density advection is more than sufficient to explain the depletion of mass over a growing depression. Consequently, the joint contribution of the net divergence and vertical motion opposes the pressure fall. This is illustrated for a case of rapid cyclogenesis in southern South America.

Jusem, Juan C.↗