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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

First-principles Modeling and Design of Solid-State Interfaces for the Protection and Use of Lithium Metal Anodes

Li-ion batteries are one of the most advanced energy storage technologies in use today. Li-ion batteries are used in a multitude of applications ranging from consumer electronics, medical devices, sensors and grid storage. However, improving the capacity and energy density delivered by current Li-ion technology requires advanced materials research into novel chemical systems. In this project we have focused particularly in the use of solid-state electrolytes with lithium metal electrodes. Research into all solid-state batteries (ASSB) with Li metal electrodes has significantly expanded in recent years, however most studies reported experimental findings, which left substantial room for theoretical and modeling work as a tool to understand and determine design principles allowing reliable and safe use of ASSBs with Li metal. Among the remaining obstacles preventing reliable use of ASSBs with a Li metal electrode, the stability of the interface between the solid electrolyte and Li metal, and the propagation/dendrite formation of Li metal and resulting mechanical degradation of the electrolyte are key phenomenon that are yet to be fully understood. In the current project we have addressed these two coupled phenomena using first principles calculations and mesoscale continuum modeling. We have obtained chemical and electrochemical stability windows for several solid electrolyte materials. Additionally, from mathematical and numerical modeling of Li protrusion and dendrite initiation during plating and stripping we have determined design criteria in terms of chemical, electrochemical, and mechanical properties and operating conditions for which stable deposition can occur. We also considered the effects of mixed electronic-ionic conduction in solid electrolytes, which has more recently been suggested as another important mechanism involved in ASSB failure. Throughout our work we have successfully addressed important questions necessary for the use of ASSB’s. We have determined guiding principles for materials properties and operating conditions necessary to operate ASSB’s. And have proposed novel solid electrolyte materials with predicted chemical stability an ionic conductivity. Although this represents significant progress in our understanding, open questions remain in order to fully develop reliable and safely operate ASSBs with Li metal. Future work, building on this project will require further experimental, theoretical and simulation efforts to address remaining questions.

25 ENERGY STORAGE↗
NASA NTRS2005

Remote Sensing Characterization of the Urban Landscape for Improvement of Air Quality Modeling

The urban landscape is inherently complex and this complexity is not adequately captured in air quality models, particularly the Community Multiscale Air Quality (CMAQ) model that is used to assess whether urban areas are in attainment of EPA air quality standards, primarily for ground level ozone. This inadequacy of the CMAQ model to sufficiently respond to the heterogeneous nature of the urban landscape can impact how well the model predicts ozone pollutant levels over metropolitan areas and ultimately, whether cities exceed EPA ozone air quality standards. We are exploring the utility of high-resolution remote sensing data and urban growth projections as improved inputs to the meteorology component of the CMAQ model focusing on the Atlanta, Georgia metropolitan area as a case study. These growth projections include "business as usual" and "smart growth" scenarios out to 2030. The growth projections illustrate the effects of employing urban heat island mitigation strategies, such as increasing tree canopy and albedo across the Atlanta metro area, in moderating ground-level ozone and air temperature, compared to "business as usual" simulations in which heat island mitigation strategies are not applied. The National Land Cover Dataset at 30m resolution is being used as the land use/land cover input and aggregated to the 4km scale for the MM5 mesoscale meteorological model and the (CMAQ) modeling schemes. Use of these data has been found to better characterize low densityhburban development as compared with USGS 1 km land use/land cover data that have traditionally been used in modeling. Air quality prediction for fiture scenarios to 2030 is being facilitated by land use projections using a spatial growth model. Land use projections were developed using the 2030 Regional Transportation Plan developed by the Atlanta Regional Commission, the regional planning agency for the area. This allows the state Environmental Protection agency to evaluate how these transportation plans will affect fbture air quality.

Quattrochi, Dale A.↗
DOE OSTI2022

Scientific challenges to characterizing the wind resource in the marine atmospheric boundary layer

Abstract. With the increasing level of offshore wind energy investment, it is correspondingly important to be able to accurately characterize the wind resource in terms of energy potential as well as operating conditions affecting wind plant performance, maintenance, and lifespan. Accurate resource assessment at a particular site supports investment decisions. Following construction, accurate wind forecasts are needed to support efficient power markets and integration of wind power with the electrical grid. To optimize the design of wind turbines, it is necessary to accurately describe the environmental characteristics, such as precipitation and waves, that erode turbine surfaces and generate structural loads as a complicated response to the combined impact of shear, atmospheric turbulence, and wave stresses. Despite recent considerable progress both in improvements to numerical weather prediction models and in coupling these models to turbulent flows within wind plants, major challenges remain, especially in the offshore environment. Accurately simulating the interactions among winds, waves, wakes, and their structural interactions with offshore wind turbines requires accounting for spatial (and associated temporal) scales from O(1 m) to O(100 km). Computing capabilities for the foreseeable future will not be able to resolve all of these scales simultaneously, necessitating continuing improvement in subgrid-scale parameterizations within highly nonlinear models. In addition, observations to constrain and validate these models, especially in the rotor-swept area of turbines over the ocean, remains largely absent. Thus, gaining sufficient understanding of the physics of atmospheric flow within and around wind plants remains one of the grand challenges of wind energy, particularly in the offshore environment. This paper provides a review of prominent scientific challenges to characterizing the offshore wind resource using as examples phenomena that occur in the rapidly developing wind energy areas off the United States. Such phenomena include horizontal temperature gradients that lead to strong vertical stratification; consequent features such as low-level jets and internal boundary layers; highly nonstationary conditions, which occur with both extratropical storms (e.g., nor'easters) and tropical storms; air–sea interaction, including deformation of conventional wind profiles by the wave boundary layer; and precipitation with its contributions to leading-edge erosion of wind turbine blades. The paper also describes the current state of modeling and observations in the marine atmospheric boundary layer and provides specific recommendations for filling key current knowledge gaps.

17 WIND ENERGY↗
NASA NTRS1994

Data assimilation and model evaluation experiment datasets

The Institute for Naval Oceanography, in cooperation with Naval Research Laboratories and universities, executed the Data Assimilation and Model Evaluation Experiment (DAMEE) for the Gulf Stream region during fiscal years 1991-1993. Enormous effort has gone into the preparation of several high-quality and consistent datasets for model initialization and verification. This paper describes the preparation process, the temporal and spatial scopes, the contents, the structure, etc., of these datasets. The goal of DAMEE and the need of data for the four phases of experiment are briefly stated. The preparation of DAMEE datasets consisted of a series of processes: (1) collection of observational data; (2) analysis and interpretation; (3) interpolation using the Optimum Thermal Interpolation System package; (4) quality control and re-analysis; and (5) data archiving and software documentation. The data products from these processes included a time series of 3D fields of temperature and salinity, 2D fields of surface dynamic height and mixed-layer depth, analysis of the Gulf Stream and rings system, and bathythermograph profiles. To date, these are the most detailed and high-quality data for mesoscale ocean modeling, data assimilation, and forecasting research. Feedback from ocean modeling groups who tested this data was incorporated into its refinement. Suggestions for DAMEE data usages include (1) ocean modeling and data assimilation studies, (2) diagnosis and theoretical studies, and (3) comparisons with locally detailed observations.

Lai, Chung-Cheng A.↗
DOE OSTI2021

From classical thermodynamics to phase-field method

Phase-field method is a density-based computational method at the mesoscale for modeling and predicting the temporal microstructure and property evolution during materials processes. The focus of this article is on connecting the most common phase-field equations to the very basic first and second laws of classical thermodynamics through rudimentary irreversible thermodynamics. It briefly discusses the relations of the continuum phase-field equations to their counter parts at the microscopic and atomic levels. It attempts to clarify the contributions of long-range elastic, electrostatic, and magnetic interactions to domain structure evolution during structural, ferroelectric, and ferromagnetic phase transformations by separating order parameter changes due to the presence of quasi-static fields and those arising from phase transformations. A few examples are presented to demonstrate the possibility of employing the phase-field method to provide guidance to designing materials for optimum properties or discovering novel mesoscale phenomena or new materials functionalities. Here, the article ends with a brief perspective on a number of potential future directions on the development and applications of phase-field method beyond its traditional applications to structural alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

TRITIUM TRANSPORT MODELING AT THE PORE SCALE IN CERAMIC BEEDER MATERIALS USING TMAP8

Fusion reactors depend on the blanket material to breed and release tritium at the same rate or faster than it is consumed by the fusion reaction. Cellular ceramic breeders are dense materials that can maintain a high tritium breeding ratio while promoting tritium release thanks to highly connected pores. Assessing the tritium breeding capabilities of these materials requires a combination of extensive experimental and modeling efforts. In this work, we develop and calibrate a multiphysics model of tritium transport. This novel model accounts for ceramic and pore diffusion, trapping and detrapping, and several surface reactions at the pore surface. We perform a sensitivity analysis and calibrate the model by comparing its predictions against experimental measurements of deuterium absorption. The calibrated model is then used to model tritium absorption in samples with different pore microstructures to investigate the effect of pore interconnectivity on tritium absorption. The model is part of the development of the multiscale, multiphysics framework for tritium transport (i.e., the Tritium Migration Analysis Program [TMAP8]), which is itself built on top of the finite element multiphysics framework MOOSE. This study demonstrates some of TMAP8’s capabilities, and is the first step toward assessing the tritium breeding capabilities of ceramic breeder material designs.

08 HYDROGEN↗
NASA NTRS2014

Introducing Multisensor Satellite Radiance-Based Evaluation for Regional Earth System Modeling

Earth System modeling has become more complex, and its evaluation using satellite data has also become more difficult due to model and data diversity. Therefore, the fundamental methodology of using satellite direct measurements with instrumental simulators should be addressed especially for modeling community members lacking a solid background of radiative transfer and scattering theory. This manuscript introduces principles of multisatellite, multisensor radiance-based evaluation methods for a fully coupled regional Earth System model: NASA-Unified Weather Research and Forecasting (NU-WRF) model. We use a NU-WRF case study simulation over West Africa as an example of evaluating aerosol-cloud-precipitation-land processes with various satellite observations. NU-WRF-simulated geophysical parameters are converted to the satellite-observable raw radiance and backscatter under nearly consistent physics assumptions via the multisensor satellite simulator, the Goddard Satellite Data Simulator Unit. We present varied examples of simple yet robust methods that characterize forecast errors and model physics biases through the spatial and statistical interpretation of various satellite raw signals: infrared brightness temperature (Tb) for surface skin temperature and cloud top temperature, microwave Tb for precipitation ice and surface flooding, and radar and lidar backscatter for aerosol-cloud profiling simultaneously. Because raw satellite signals integrate many sources of geophysical information, we demonstrate user-defined thresholds and a simple statistical process to facilitate evaluations, including the infrared-microwave-based cloud types and lidar/radar-based profile classifications.

Planetary Boundary Layer↗
NASA NTRS2018

Comparisons of Bin and Bulk Microphysics Schemes in Simulations of Topographic Winter Precipitation with Radar and Radiometer Measurements

The Hebrew University Cloud Model (HUCM) bin scheme and the Thompson bulk scheme in the Weather Research and Forecasting (WRF) model are compared to assess biases often found in simulated brightness temperature and radar reflectivity. Compared to our preceding study that evaluated several bulk schemes in the WRF model, the current study obtains a reduction of the bias from excessive microwave scattering by precipitation ice for both HUCM bin and the Thompson bulk microphysics schemes for a topographic winter precipitation event associated with an atmospheric river. The Thompson particle size distributions (PSDs) and snow particle density assumption are implemented into the Goddard Satellite Data Simulator Unit (G‐SDSU) and have produced improvements. Despite the greater sophistication of the bin scheme in representing cloud and precipitation processes, the simulation with the Thompson bulk scheme is generally in better agreement with observations for this winter event. The explicitly resolved hydrometeor PSDs in HUCM enable analysis of mass spectra variations in response to changes in microphysics assumptions. Two HUCM sensitivity runs tested the enhancement of snow particle breakup and the influence of ice nuclei (IN) concentration. Higher IN concentration resulted in increased snow mass and broadened the spectrum toward small‐size particles. Modified snow mass spectra and resultant changes in graupel contributed to modifications in scattering and reflectivity simulations. The article demonstrates the bin scheme's capability to provide a new means to improve our understanding of uncertainties in mesoscale weather models and radiative transfer models.

microwave remote sensing↗
DOE OSTI2026

Assessing the Impact of a Forest Canopy on Near-Surface Wind Statistics

Representing the forest canopy in atmospheric numerical models should improve simulated winds within and above the canopy up to a few hundred meters above the ground. Here, in this study, we implement a forest canopy parameterization into the Weather Research and Forecasting (WRF) Model in a large-eddy simulation (LES) mode by applying drag forces across multiple layers within the canopy height. We use unique observations from the Lidar Experiments for Assessing Flow over Forests (LEAFF) field campaign at the Wind River Experimental Forest (WREF) in the U.S. Pacific Northwest to evaluate model performance. In a 2-day case study, the canopy parameterization improved wind predictions both within and above the canopy, particularly during the daytime and at finer grid resolution. Without it, winds were frequently overpredicted above the canopy. Similarly, derived quantities such as the wind shear index also yielded estimates closer to observations with the canopy parameterization implemented. These findings suggest that representing the canopy using drag forces alone can improve simulated mean winds up to 200 m above the surface. Furthermore, second-order statistical moments of wind were more sensitive to canopy density than first-order moments, especially during the daytime. This increased sensitivity and the improved daytime performance in wind speed—evidenced by the lowest bias from observations (3% compared to 20% over diurnal cycle)—imply that winds above the canopy layer are strongly influenced by how well turbulence above the canopy is modeled. The results of this study can serve as a foundation for parameterizing forest canopy effects in coarser weather forecast models.

Energy - Wind↗
DOE OSTI2021

Evolution of Idealized Vortices in Monsoon-Like Shears: Application to Monsoon Depressions

Abstract This study examines processes fundamental to the development of South Asian monsoon depressions using an array of integrations of an idealized convection-permitting numerical model. In each integration, a wave of initially small amplitude is subjected to a different amount of vertical and meridional wind shear, with temperature and moisture fields constructed according to realistic constraints. Based on the evolution of this disturbance into monsoon depression–like vortices, two features of the background environment emerge as important: the low-level gradient of moist static energy (MSE) and the low-level meridional shear. As the low-level MSE gradient steepens, the disturbance becomes stronger and produces more rain. This strengthening results from the interaction of the vortex with latent heat release by convection that is in turn organized by positive MSE advection in the northerly flow west of the vortex. In this region of advection, moister air from the north ascends along upward-sloping isentropes, driving moist convection. The disturbance also becomes stronger with increasing meridional shear, which makes the environment more barotropically unstable. The absence of either of these two features of the background environment prevents substantial growth of the disturbance. Our results suggest that monsoon depression growth in South Asia is fostered by the coexistence of a strong low-level MSE gradient with strong meridional wind shear associated with the monsoon trough.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Lake Huron Enhances Snowfall Downwind of Lake Erie: A Modeling Study of the 2010 New Year’s Lake Effect Snowfall Event

In the Northeast and Great Lakes regions of the United States, the influence of multiple lakes on the overlying air can greatly affect lake-effect snowfall over downwind communities. To assess the impact of Lake Huron on snowfall downwind of Lake Erie, we simulated a lake-effect snow event which occurred from 1 to 6 January 2010 using the Regional Atmospheric Modeling System (RAMS). We found that the presence of Lake Huron enhances snowfall downwind of Lake Erie by almost 20% and leads to much heavier local snowfall totals than when Lake Huron is not present. This increase in snowfall is due to a lake-to-lake (L2L) convective band, as secondary circulations associated with lake-effect convection form over Lake Huron and persist overland between the lakes before reintensifying over Lake Erie. As these secondary circulations move over Lake Erie, lifting within the upward branch of the circulations accelerates the development of convection within the L2L band. Surface fluxes and convection over Lake Huron deepen the boundary layer, promoting deeper vertical development of this L2L band over Lake Erie. However, although this boundary layer modification strengthens the L2L band, we found that it actually reduced lake-effect snowfall produced by wind-parallel bands (WPB) in other parts of the lake. This indicates that boundary layer modification from upstream lakes may impact L2L bands differently than WPB which are not connected to upstream secondary circulations.

Lake effects↗
DOE OSTI2022

Assimilation of a Coordinated Fleet of Uncrewed Aircraft System Observations in Complex Terrain: Observing System Experiments

Uncrewed aircraft system (UAS) observations from the Lower Atmospheric Profiling Studies at Elevation–A Remotely-Piloted Aircraft Team Experiment (LAPSE-RATE) field campaign were assimilated into a high-resolution configuration of the Weather Research and Forecasting (WRF) Model. The impact of assimilating targeted UAS observations in addition to surface observations was compared to that obtained when assimilating surface observations alone using observing system experiments (OSEs) for a terrain-driven flow case and a convection initiation (CI) case observed within Colorado’s San Luis Valley (SLV). The assimilation of UAS observations in addition to surface observations results in a clear increase in skill for both flow regimes over that obtained when assimilating surface observations alone. For the terrain-driven flow case, the UAS observations improved the representation of thermal stratification across the northern SLV, which produced stronger upvalley flow over the eastern half of the SLV that better matched the observations. For the CI case, the UAS observations improved the representation of the pre-convective environment by reducing dry biases across the SLV and over the surrounding terrain. This led to earlier CI and more organized convection over the foothills that spilled outflows into the SLV, ultimately helping to increase low-level convergence and CI there. In addition, the importance of UAS capturing an outflow that originated over the Sangre de Cristo Mountains and triggered CI is discussed. These outflows and subsequent CI were not well captured in the simulation that assimilated surface observations alone. We report that observations obtained with a fleet of UAS are shown to notably improve high-resolution analyses and short-term predictions of two very different mesogamma-scale weather events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Influence of the distribution of expansive sites in aggregates on microscopic damage caused by alkali-silica reaction: Insights into the mechanical origin of expansion

The origin of damage in concrete due to the alkali-silica reaction (ASR) is attributed to the expansion site in the aggregate. To investigate the cracking process of the aggregate during ASR and its consequences on concrete damage, the effect of the distribution of the expansive sites in the aggregate on ASR expansion and the crack patterns must be evaluated. Thus, in this study, a mesoscale discrete model was applied to ASR modeling to represent the propagation of cracking during ASR accurately. The distribution of the expansive sites in the aggregate was based on the gel pocket and reaction rim models, which are two ASR mechanisms reported in the literature. These two expansion models highlight the different crack patterns obtained based on the aggregate characteristics. Further, the expansion cracking processes determined based on the gel pocket and reaction rim models are consistent with the evolution of cracking with the expansion level.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Sensitivity of the strength and toughness of concrete to the properties of the interfacial transition zone

We report Civil infrastructure is made primarily of concrete structures or components and therefore understanding durability and fracture behavior of concrete is of utmost importance. Concrete contains an interfacial transition zone (ITZ), a porous region surrounding the aggregates, that is often considered to be the weakest region in the concrete. The ITZ is poorly characterized and property estimates for the ITZ differ considerably. In this simulation study, representative concrete mesostructures are produced by packing coarse aggregates with realistic geometries into a mortar matrix. A meshless numerical method, peridynamics, is utilized to simulate the mechanical response including fracture under uniaxial compression and tension. The sensitivity of the stiffness and fracture toughness of the samples to the ITZ properties is computed, showing strong relationships between the ITZ properties and the effective modulus and effective yield strength of the concrete. These results provides insight into the influence of the poorly characterized ITZ on the stiffness and strength of concrete. This work showcases the applicability of peridynamics to concrete systems, matching experimental strength and modulus values. Additionally, relationships between the ITZ’s mechanical properties and the overall concrete strength and stiffness are presented to enable future design decisions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

On the computational solution of vector-density based continuum dislocation dynamics models: A comparison of two plastic distortion and stress update algorithms

Continuum dislocation dynamics models of mesoscale plasticity consist of dislocation transport-reaction equations coupled with crystal mechanics equations. The coupling between these two sets of equations is such that dislocation transport gives rise to the evolution of plastic distortion (strain), while the evolution of the latter fixes the stress from which the dislocation velocity field is found via a mobility law. Earlier solutions of these equations employed a staggered solution scheme for the two sets of equations in which the plastic distortion was updated via time integration of its rate, as found from Orowan's law. In this work, we show that such a direct time integration scheme can suffer from accumulation of numerical errors. We introduce an alternative scheme based on field dislocation mechanics that ensures consistency between the plastic distortion and the dislocation content in the crystal. The new scheme is based on calculating the compatible and incompatible parts of the plastic distortion separately, and the incompatible part is calculated from the current dislocation density field. Stress field and dislocation transport calculations were implemented within a finite element based discretization of the governing equations, with the crystal mechanics part solved by a conventional Galerkin method and the dislocation transport equations by the least squares method. A simple test was first performed to show the accuracy of the two schemes for updating the plastic distortion, which shows that the solution method based on field dislocation mechanics is more accurate. This method then was used to simulate an austenitic steel crystal under uniaxial loading and multiple slip conditions. By considering dislocation interactions caused by junctions, a hardening rate similar to discrete dislocation dynamics simulation results was obtained. Finally, the simulations show that dislocations exhibit some self-organized structures as the strain is increased.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Micro-beam bending of FCC bicrystals: A comparison between defect dynamics simulations and experiments

To understand the role of the grain boundary (GB) in plasticity at small scale, a concurrently coupled mesoscale plasticity model was developed to simulate micro-bending of bicrystalline micron-sized beams. By coupling dislocation dynamics (DD) with a finite element model (FEM), a novel defect dynamics model provides the means to investigate intricate interactions between dislocations and GBs under various loading conditions. Our simulations of micro-bending agree well with corresponding micro-bending experiments, and they show that mechanical response of bicrystals could have not only hardening but also softening depending on the characters of the GB. In addition, changing the location of the GB in the microbeams results in different mechanical responses; GBs located at the neutral plane show softening compared to single crystals, while inclined GBs located halfway along the length of the beam show little effect. In conclusion, simulation results could provide a clear picture on detailed dislocation-GB interactions, and quantitative resolved shear stress analysis supplemented by dislocation density distribution is used to analyze the mechanical response of bicrystalline samples.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Impacts of the morphology of new neighborhoods on microclimate and building energy

In anticipation of emerging global urbanization and consequent increases in energy use and carbon dioxide emissions, better understanding and quantification of climate effects on energy use in cities are needed, requiring coordinated research into large-scale, regional, and microclimate impacts to and from the city structure. The methodology described here addresses this need by (1) demonstrating a process for creating and testing example morphologies for new neighborhoods for their impact on local and regional meteorology within a two-way-coupled four-domain nested mesoscale weather model (6 km horizontal resolution outer domain, 90 m horizontal innermost domain) and (2) allocating resulting building-level meteorological profiles to each building in a neighborhood for parallel computation of building-by-building energy use. Our Chicago Loop test case shows that the morphology of even a small new added development to a neighborhood affects not only its own microclimate, but also the microclimate of the original neighborhood to which the development was added, and that these changes in microclimate affect both neighborhoods’ building energy use. Furthermore, this method represents an important step toward quantifying and analyzing the relationships among climatic conditions, urban morphology, and energy use and using these relationships to inform energy-efficient urban development and planning.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗