Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesoscale modeling”

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 199 records · Page 11

Hurricane Forecasts with a Global Mesoscale-resolving Model on the NASA Columbia Supercomputer Preliminary Simulations of Hurricane Katrina (2005)

It is known that the General Circulation Models (GCMs) have sufficient resolution to accurately simulate hurricane near-eye structure and intensity. To overcome this limitation, the mesoscale-resolving finite-element GCM (fvGCM) has been experimentally deployed on the NASA Columbia supercomputer, and its performance is evaluated choosing hurricane Katrina as an example in this study. On late August 2005 Katrina underwent two stages of rapid intensification and became the sixth most intense hurricane in the Atlantic. Six 5-day simulations of Katrina at both 0.25 deg and 0.125 deg show comparable track forecasts, but the 0,125 deg runs provide much better intensity forecasts, producing center pressure with errors of only +/- 12 hPa. The 0.125 deg simulates better near-eye wind distributions and a more realistic average intensification rate. A convection parameterization (CP) is one of the major limitations in a GCM, the 0.125 deg run with CP disabled produces very encouraging results.

Shen, B.-W.↗

Martian Atmospheric Circulation and Climate

This project was submitted with four work tasks. The reviews recommended focusing on the following three, with an appropriately reduced budget: 1. In the first, we proposed to use Thermal Emission Spectrometer (TES) data to examine the general circulation of the atmosphere and to better constrain and validate a Mars General Circulation Model (GCM). 2. In the second, we proposed to examine processes responsible for the injection of dust into the atmosphere. We will conduct a global, seasonal survey of dust lifting events, and specifically of dust devil activity, using MGS MOC, MGS TES, and Odyssey THEMIS data. The results of this survey will be compared with model predictions of dust devil and other dust lifting processes. It is expected that some of the observational results will motivate specific simulations with the GCM and mesoscale models. 3. In the third task, we proposed to extract seasonal and spatial near-surface wind information from aeolian streaks. Detailed analysis will focus on streaks in the polar regions and in regions of large topography (such as Tharsis). The streak information will be used along with the mesoscale model to examine the processes that control local and regional scale flow patterns.

Richardson, Mark I.↗

Analysis of the tropospheric water distribution during FIRE 2

The Penn State/NCAR mesoscale model, as adapted for use at ARC, was used as a testbed for the development and validation of cloud models for use in General Circulation Models (GCM's). This modeling approach also allows us to intercompare the predictions of the various cloud schemes within the same dynamical framework. The use of the PSU/NCAR mesoscale model also allows us to compare our results with FIRE-II (First International Satellite Cloud Climatology Project Regional Experiment) observations, instead of climate statistics. Though a promising approach, our work to date revealed several difficulties. First, the model by design is limited in spatial coverage and is only run for 12 to 48 hours at a time. Hence the quality of the simulation will depend heavily on the initial conditions. The poor quality of upper-tropospheric measurements of water vapor is well known and the situation is particularly bad for mid-latitude winter since the coupling with the surface is less direct than in summer so that relying on the model to spin-up a reasonable moisture field is not always successful. Though one of the most common atmospheric constituents, water vapor is relatively difficult to measure accurately, especially operationally over large areas. The standard NWS sondes have little sensitivity at the low temperatures where cirrus form and the data from the GOES 6.7 micron channel is difficult to quantify. For this reason, the goals of FIRE Cirrus II included characterizing the three-dimensional distribution of water vapor and clouds. In studying the data from FIRE Cirrus II, it was found that no single special observation technique provides accurate regional distributions of water vapor. The Raman lidar provides accurate measurements, but only at the Hub, for levels up to 10 km, and during nighttime hours. The CLASS sondes are more sensitive to moisture at low temperatures than are the NWS sondes, but the four stations only cover an area of two hundred kilometers on a side. The aircraft give the most accurate measurements of water vapor, but are limited in spatial and temporal coverage. This problem is partly alleviated by the use of the MAPS analyses, a four-dimensional data assimilation system that combines the previous 3-hour forecast with the available observations, but its upper-level moisture analyses are sometimes deficient because of the vapor measurement problem. An attempt was made to create a consistent four-dimensional description of the water vapor distribution during the second IFO by subjectively combining data from a variety of sources, including MAPS analyses, CLASS sondes, SPECTRE sondes, NWS sondes, GOES satellite analyses, radars, lidars, and microwave radiometers.

Westphal, Douglas L.↗

The Effect of Surface Heating Heterogeneity on Boundary Layer Height and Its Dependence on Background Wind Speed

The planetary boundary layer height (PBLH) is a fundamental variable of the planetary boundary layer (PBL). The effect of surface heterogeneity is challenging in the PBL parameterization for numerical weather and climate models. Here we use large eddy simulation data to investigate the impact of surface heating heterogeneity on the PBLH spatial variations and mean values over a domain of typical mesoscale model grid size. It is found that this impact depends on background wind condition. Variance decomposition and Fourier spectra both reveal that for the fully developed sheared convective PBL, surface heterogeneity contributes little to the spatial variation of PBLH, except for the case with heterogeneity scale of 14.4 km. This suggests that surface heterogeneity with length scales less than ~10 km, which is typical for mesoscale modeling, does not induce significant subgrid spatial variation of PBLH. The domain-averaged PBLH is found to generally increase as surface heterogeneity scale decreases, until reaching the scale on the order of PBLH. The findings can be used to inform the parameterization development of PBL processes over heterogeneous surface.

54 ENVIRONMENTAL SCIENCES↗

Routine High-Resolution Forecasts/Analyses for the Pacific Disaster Center: User Manual

Enclosed herein is our HWCMO user manual. This manual constitutes the final report for our NASA/PDC grant, NASA NAG5-8730, "Routine High Resolution Forecasts/Analysis for the Pacific Disaster Center". Since the beginning of the grant, we have routinely provided experimental high resolution forecasts from the RSM/MSM for the Hawaii Islands, while working to upgrade the system to include: (1) a more robust input of NCEP analyses directly from NCEP; (2) higher vertical resolution, with increased forecast accuracy; (3) faster delivery of forecast products and extension of initial 1-day forecasts to 2 days; (4) augmentation of our basic meteorological and simplified fireweather forecasts to firedanger and drought forecasts; (5) additional meteorological forecasts with an alternate mesoscale model (MM5); and (6) the feasibility of using our modeling system to work in higher-resolution domains and other regions. In this user manual, we provide a general overview of the operational system and the mesoscale models as well as more detailed descriptions of the models. A detailed description of daily operations and a cost analysis is also provided. Evaluations of the models are included although it should be noted that model evaluation is a continuing process and as potential problems are identified, these can be used as the basis for making model improvements. Finally, we include our previously submitted answers to particular PDC questions (Appendix V). All of our initially proposed objectives have basically been met. In fact, a number of useful applications (VOG, air pollution transport) are already utilizing our experimental output and we believe there are a number of other applications that could make use of our routine forecast/analysis products. Still, work still remains to be done to further develop this experimental weather, climate, fire danger and drought prediction system. In short, we would like to be a part of a future PDC team, if at all possible, to further develop and apply the system for the Hawaiian and other Pacific Islands as well as the entire Pacific Basin.

Roads, John↗

Numerical Simulations of Precipitation Processes, Microphysics, and Microwave Radiative Properties of flood Producing Storms in Mediterranean & Adriatic Basins

A comprehensive understanding of the meteorological and microphysical nature of Mediterranean storms requires a combination of in situ data analysis, radar data analysis, and satellite data analysis, effectively integrated with numerical modeling studies at various scales. An important aspect of understanding microphysical controls of severe storms, is first understanding the meteorological controls under which a storm has evolved, and then using that information to help characterize the dominant microphysical processes. For hazardous Mediterranean storms, highlighted by the October 5-6, 1998 Friuli flood event in northern Italy, a comprehensive microphysical interpretation requires an understanding of the multiple phases of storm evolution. This involves intense convective development, Sratiform decay, orographic lifting, and sloped frontal lifting processes, as well as the associated vertical motions and thermodynamical instabilities governing physical processes that effect details of the size distributions and fall rates of the various types of hydrometeors found within the storm environment. This talk overviews the microphysical elements of a severe Mediterranean storm in such a context, investigated with the aid of TRMM satellite and other remote sensing measurements, but guided by a nonhydrostatic mesoscale model simulation of the Friuli flood event. The data analysis for this paper was conducted by my research groups at the Global Hydrology and Climate Center in Huntsville, AL and Florida State University in Tallahassee, and in collaboration with Dr. Alberto Mugnai's research group at the Institute of Atmospheric Physics in Rome. The numerical modeling was conducted by Professor Oreg Tripoli and Ms. Giulia Panegrossi at the University of Wisconsin in Madison, using Professor Tripoli's nonhydrostatic modeling system (NMS). This is a scalable, fully nested mesoscale model capable of resolving nonhydrostatic circulations from regional scale down to cloud scale and below.

Smith, Eric A.↗

Bicrystallography-informed Frenkel–Kontorova model for interlayer dislocations in strained 2D heterostructures

In recent years, van der Waals (vdW) heterostructures and homostructures, which consist of stacks of two-dimensional (2D) materials, have risen to prominence due to their association with exotic quantum phenomena originating from correlated electronic states harbored by them. Atomistic scale relaxation effects play an extremely important role in the electronic scale quantum physics of these systems, providing means of manipulation of these materials and allowing them to be tailored for emergent technologies. We investigate such structural relaxation effects in this work using atomistic and mesoscale models, within the context of twisted bilayer graphene — a well-known heterostructure system that features moiré patterns arising from the lattices of the two graphene layers. For small twist angles, atomic relaxation effects in this system are associated with the natural emergence of interface dislocations or strain solitons, which result from the cyclic nature of the generalized stacking fault energy (GSFE), that measures the interface energy based on the relative movement of the two layers. Here, in this work, we first demonstrate using atomistic simulations that atomic reconstruction in bilayer graphene under a large twist also results from interface dislocations, although the Burgers vectors of such dislocations are considerably smaller than those observed in small-twist systems. To reveal the translational invariance of the heterointerface responsible for the formation of such dislocations, we derive the translational symmetry of the GSFE of a 2D heterostructure using the notions of coincident site lattices (CSLs) and displacement shift complete lattices (DSCLs). The workhorse for this exercise is a recently developed Smith normal form bicrystallography framework. Next, we construct a bicrystallography-informed and frame-invariant Frenkel–Kontorova model, which can predict the formation of strain solitons in arbitrary 2D heterostructures, and apply it to study a heterostrained, large-twist bilayer graphene system. Our mesoscale model is found to produce results consistent with atomistic simulations. We anticipate that the model will be invaluable in predicting structural relaxation and for providing insights into various heterostructure systems, especially in cases where the fundamental unit cell is large and therefore, atomistic simulations are computationally expensive.

2D heterostructures↗

Investigation of Aerodynamic and Radiometric Land Surface Temperatures

The surface temperature, T(sub s), of a land surface measured by a radiometer, T(sub s,r), and the temperature "felt" by the air, T(sub aero), often differ significantly and are difficult if not impossible to define rigorously. However, recent studies conducted by the principal investigators with several land surface models suggest that this problem can be largely resolved. The main goal of this project was to use model-based and empirical studies to improve understanding and reconcile the difference between T(sub s,r) and T(sub aero), while maintaining consistency within the models and with theory and data. The results from this effort have contributed progress towards the effective use of remotely sensed surface temperature measurements taken from an arbitrary view angle over a partial canopy cover for producing high quality sensible and latent heat flux estimates. In addition, we have developed parameterizations that are designed to improve the representation of the roughness length for heat in climate and mesoscale models. This not only provides improved representation of surface energy balance in such models, but should also facilitate the use of surface temperature measurements for validating or updating the surface temperature produced by SVATs (soil-vegetation-atmosphere schemes) in climate or mesoscale models.

Friedl, Mark↗

Utilization of VAS satellite data in the initialization of an oceanic-cyclogenesis simulation

A series of experiments was performed to test various method of incorporating Visible Infrared Spin Scan Radiometer Atmospheric Sounder (VAS)-sounding data into the initial conditions of the Penn State University/National Center for Atmospheric mesoscale model. The VAS data for this ocean-cyclogenesis case consist of 110 irregularly distributed temperature and humidity soundings located over the North Pacific Ocean and apply at approximately 1200 GMT 10 November 1981. Various methods of utilizing VAS data in the initial condition of a mesoscale model were evaluated.

Douglas, Sharon G.↗

Utilization of VAS satellite data in the initialization of an oceanic cyclogenesis simulation

A series of experiments was performed to test various methods of incorporating Visible Infrared Spin Scan Radiometer Atmospheric Sounder (VAS)-sounding data into the initial conditions of the Penn State University/National Center for Atmospheric mesoscale model. The VAS data for this ocean-cyclogenesis case consist of 110 irregularly distributed temperature and humidity soundings located over the North Pacific Ocean and apply at approximately 1200 GMT November 10, 1981. Various methods of utilizing VAS data in the initial condition of a mesoscale model were evaluated.

Douglas, Sharon G.↗

Atmospheric Constraints on Landing Site Selection

The Martian atmosphere is a significant part of the environment that the Mars Exploration Rovers (MER) will encounter. As such, it imposes important constraints on where the rovers can and cannot land. Unfortunately, as there are no meteorological instruments on the rovers, there is little atmospheric science that can be accomplished, and no scientific preference for landing sites. The atmosphere constrains landing site selection in two main areas, the entry descent and landing (EDL) process and the survivability of the rovers on the surface. EDL is influenced by the density profile and boundary layer winds (up to altitudes of 5 to 10 km). Surface survivability involves atmospheric dust, temperatures and winds. During EDL, the atmosphere is used to slow the lander down, both ballistically and on the parachute. This limits the maximum elevation of the landing site to -1.3 km below the MOLA reference aeroid. The landers need to encounter a sufficiently dense atmosphere to be able to stop, and the deeper the landing site, the more column integrated atmosphere the lander can pass through before reaching the surface. The current limit was determined both by a desire to be able to reach the hematite region and by a set of atmosphere models we developed for EDL simulations. These are based on Thermal Emission Spectrometer (TES) atmospheric profile measurements, Ames Mars General Circulation Model (MGCM) results, and the 1-D Ames GCM radiative/convective model by J. Murphy. The latter is used for the near surface diurnal cycle. The current version of our model encompasses representative latitude bands, but we intend to make specific models for the final candidate landing sites to insure that they fall within the general envelope. The second constraint imposed on potential landing sites through the EDL process is the near surface wind. The wind in the lower approximately 5 km determines the horizontal velocity that the landers have when they land. Due to the mechanics of the landing process, the total velocity (including both the horizontal and vertical components) determines whether or not the landers are successful. Unfortunately, the landing system has no easy way to nullify any horizontal velocity imparted by the wind, so the landing sites selected need to have as little wind as possible. In addition to the mean wind velocity, the landing system is sensitive to vertical wind shear in the lowest kilometer or so. Wind shear can deflect the retro rockets (RADs) from their nominal vertical orientation producing unwanted horizontal spacecraft velocities. Both mean velocity and wind shear are dominated by the the local topography and other surface properties (in particular albedo and thermal inertia which control the surface temperature). This is seen even in simplified 2-D mesoscale models. The effects in a fully 3-D model are expected to he even more topographically dependent. In particular there is potential for wind channeling in canyons and other terrain features. Boundary layer winds and wind shear are currently being modeled based on terrestrial data and boundary layer scaling laws modified for Martian conditions. We hope to supplement this with mesoscale model results (from several sources) once the number of landing sites is reduced to a manageable number.

Kass, David M.↗

Nonlocal Effect of Percolated Particle Networks on Viscoelasticity of Polymer–Filler Nanocomposites: A Mesoscale Simulation Study

With nanoparticles (NPs) as fillers, polymer nanocomposites (PNCs) usually exhibit enhanced mechanical properties. However, a direct connection between the microscopic structural relaxation and macroscopic mechanical properties of PNCs remains to be established. To investigate the micro-to-macro connection, we develop a mesoscale model, in which the NPbridging polymer chains are represented by a dynamic bonded interaction between NPs, and the bulk polymer matrix is implicitly modeled by overdamped Langevin dynamics. Extensive equilibrium simulations are performed to quantify the microscopic dynamics of model PNCs. Systematic analyses of modified Rouse dynamics, dynamic structure factor, and relaxation modulus uncover that the microscopic relaxation dynamics of PNCs are significantly decelerated across different length scales because of nonlocal effects of percolated particle networks a phenomenon that has not been adequately captured in prior simulation studies. We find that NPvolume- fraction and NP-bonding-energy barrier are the two critical variables that affect bulk viscoelasticity the most. The proposed mesoscale model is versatile and provides a powerful framework for studying structure−property relations of different PNCs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Use of Satellite Data Assimilation to Infer Land Surface Thermal Inertia

There are two important but observationally uncertain parameters in the grid averaged surface energy budgets of mesoscale models - surface moisture availability and thermal heat capacity. A technique has been successfully developed for assimilating Geostationary Operational Environmental Satellite (GOES) skin temperature tendencies during the mid-morning time frame to improve specification of surface moisture. In a new application of the technique, the use of satellite skin temperature tendencies in early evening is explored to improve specification of the surface thermal heat capacity. Together, these two satellite assimilation constraints have been shown to significantly improve the characterization of the surface energy budget of a mesoscale model on fine spatial scales. The GOES assimilation without the adjusted heat capacity was run operationally during the International H2O Project on a 12-km grid. This paper presents the results obtained when using both the moisture availability and heat capacity retrievals in concert. Preliminary results indicate that retrieved moisture availability alone improved the verification statistics of 2-meter temperature and dew point forecasts. Results from the 1.5 month long study period using the bulk heat capacity will be presented at the meeting.

Lapenta, William↗

A mesoscale 3D model of irradiated concrete informed via a 2.5 U-Net semantic segmentation

The concrete biological shield in light-water reactors is exposed to neutron and gamma irradiation, which deteriorates the concrete’s mechanical properties in the long term. To assess the irradiation-induced damage, predictive mechanical models are developed and used in parallel with the characterization of irradiated concrete samples. Realistic 3D simulation domains can drastically improve a model’s prediction. In this work, we utilized x-ray computed tomography (XCT) data of a concrete specimen to reconstruct its 3D microstructure. The XCT data shows low contrast between the concrete’s aggregates and cement paste, resulting in poor image segmentation when using traditional unsupervised techniques. To address this issue, we developed and trained a 2.5D U-Net model on only 24 pre-labeled XCT layers to segment 651 layers of the XCT data. The overall F1-score of the model is approximately 96%. Then, we created a 3D finite element (FE) mesh based on the stack of segmented images. The FE model contains radiation-induced expansion, damage, and creep. The constitutive equations are adapted to each phase (aggregates and cement paste). Here, we simulated the effects of neutron irradiation in the concrete specimen as well as the specimen’s mechanical response to uniaxial compression. Finally, model validation was performed using experimental data on similar concrete specimens in the literature.

2.5D U-Net↗

On the implementation of dislocation reactions in continuum dislocation dynamics modeling of mesoscale plasticity

The continuum dislocation dynamics framework for mesoscale plasticity is intended to capture the dislocation density evolution and the deformation of crystals when subjected to mechanical loading. It does so by solving a set of transport equations for dislocations concurrently with crystal mechanics equations, with the latter being cast in the form of an eigenstrain problem. Incorporating dislocation reactions in the dislocation transport equations is essential for making such continuum dislocation dynamics predictive. A formulation is proposed to incorporate dislocation reactions in the transport equations of the vector density-based continuum dislocation dynamics. This formulation aims to rigorously enforce dislocation line continuity using the concept of virtual dislocations that close all dislocation loops involved in cross slip, annihilation, and glissile and sessile junction reactions. The addition of virtual dislocations enables us to accurately enforce the divergence free condition upon the numerical solution of the dislocation transport equations for all slip systems individually. A set of tests were performed to illustrate the accuracy of the formulation and the solution of the transport equations within the vector density-based continuum dislocation dynamics. Comparing the results from these tests with an earlier approach in which the divergence free constraint was enforced on the total dislocation density tensor or the sum of two densities when only cross slip is considered shows that the new approach yields highly accurate results. Bulk simulations were performed for a face centered cubic crystal based on the new formulation and the results were compared with discrete dislocation dynamics predictions of the same. The microstructural features obtained from continuum dislocation dynamics were also analyzed with reference to relevant experimental observations.

36 MATERIALS SCIENCE↗

Coupled Ocean Wave-Atmosphere Models for Offshore Wind Energy Applications

Coupling ocean wave models to mesoscale atmospheric models is necessary to represent the effect of waves on wind turbine hub-height winds. In this report we provide a review of the most widely used ocean waves models and the phased-averaged spectral wave modeling paradigm that they are based on. Methodologies used to couple these wave models to mesoscale atmospheric models are described along with details of existing coupled modeling systems. We summarize impacts on offshore wind resource assessments that have been investigated with such coupled modeling systems to-date. Specifically, coupling in the North and Baltic Seas was shown to have a small negative effect on the offshore wind energy resource at weak-to-moderate wind speeds, with little impact at higher wind speeds. Finally, limitations of these existing coupled modeling systems and impact assessments are discussed, including examples of and the potential for improved parameterizations of the air-sea fluxes, and the use of fine-scaled simulations using phased-resolved wave models.

17 WIND ENERGY↗