Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model parameterization”

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 1,387 records · Page 77

Parameter Tuning and Calibration of RegCM3 with MIT-Emanuel Cumulus Parameterization Scheme over CORDEX East Asian Domain

In this study, we calibrated the performance of regional climate model RegCM3 with Massachusetts Institute of Technology (MIT)-Emanuel cumulus parameterization scheme over CORDEX East Asia domain by tuning the selected seven parameters through multiple very fast simulated annealing (MVFSA) sampling method. The seven parameters were selected based on previous studies, which customized the RegCM3 with MIT-Emanuel scheme through three different ways by using the sensitivity experiments. The responses of model results to the seven parameters were investigated. Since the monthly total rainfall is constrained, the simulated spatial pattern of rainfall and the probability density function (PDF) distribution of daily rainfall rates are significantly improved in the optimal simulation. Sensitivity analysis suggest that the parameter “relative humidity criteria” (RH), which has not been considered in the default simulation, has the largest effect on the model results. The responses of total rainfall over different regions to RH were examined. Positive responses of total rainfall to RH are found over northern equatorial western Pacific, which are contributed by the positive responses of explicit rainfall. Followed by an increase of RH, the increases of the low-level convergence and the associated increases in cloud water favor the increase of the explicit rainfall. The identified optimal parameters constrained by the total rainfall have positive effects on the low-level circulation and the surface air temperature. Furthermore, the optimized parameters based on the extreme case are suitable for a normal case and the model’s new version with mixed convection scheme.

emanual cumulus parameterization↗

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↗

Determinants of photochemical characteristics of the photosynthetic electron transport chain of maize

The photosynthetic electron transport chain (ETC) is the bridge that links energy harvesting during the photophysical reactions at one end and energy consumption during the biochemical reactions at the other. Its functioning is thus fundamental for the proper balance between energy supply and demand in photosynthesis. Currently, there is a lack of understanding regarding how the structural properties of the ETC are affected by nutrient availability and plant developmental stages, which is a major roadblock to comprehensive modeling of photosynthesis. Redox parameters reflect the structural controls of ETC on the photochemical reactions and electron transport. We conducted joint measurements of chlorophyll fluorescence (ChlF) and gas exchange under systematically varying environmental conditions and growth stages of maize and sampled foliar nutrient contents. We utilized the recently developed steady-state photochemical model to infer redox parameters of electron transport from these measurements. We found that the inferred values of these photochemical redox parameters varied with leaf macronutrient content. These variations may be caused either directly by these nutrients being components of protein complexes on the ETC or indirectly by their impacts on the structural integrity of the thylakoid and feedback from the biochemical reactions. Also, the redox parameters varied with plant morphology and developmental stage, reflecting seasonal changes in the structural properties of the ETC. Our findings will facilitate the parameterization and simulation of complete models of photosynthesis.

59 BASIC BIOLOGICAL SCIENCES↗

Analyzing Unfavored Fragmentation Functions Using NJL-Jet Model

The Nambu-Jona-Lasinio (NJL)-jet model provides a sound framework for calculating the fragmentation functions in an effective chiral quark theory, where the momentum and isospin sum rules are satisfied without the introduction of ad hoc parameters. The most recent version of the model includes the fragmentation of the light and strange quarks to pions, kaons, nucleons, and antinucleons; where the effects of the production of secondary pions and kaons from the vector mesons {rho}, K* and {phi} are also calculated. The results for the model fragmentation function exhibit a qualitative agreement with the empirical parameterizations. The results also allow to test, within the model assumptions, several assumption in parametrizations of the unfavored fragmentation functions used in empirical fits to the experimental data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Projecting Permafrost Thaw of Sub-Arctic Tundra With a Thermodynamic Model Calibrated to Site Measurements

Northern circumpolar permafrost thaw affects global carbon cycling, as large amounts of stored soil carbon becomes accessible to microbial breakdown under a warming climate. The magnitude of carbon release is linked to the extent of permafrost thaw, which is locally variable and controlled by soil thermodynamics. Soil thermodynamic properties, such as thermal diffusivity, govern the reactivity of the soil-atmosphere thermal gradient, and are controlled by soil composition and drainage. In order to project permafrost thaw for an Alaskan tundra experimental site, we used seven years of site data to calibrate a soil thermodynamic model using a data assimilation technique. The model reproduced seasonal and interannual temperature dynamics for shallow (5–40 cm) and deep soil layers (2–4 m), and simulations of seasonal thaw depth closely matched observed data. The model was then used to project permafrost thaw at the site to the year 2100 using climate forcing data for three future climate scenarios (RCP 4.5, 6.0, and 8.5). Minimal permafrost thawing occurred until mean annual air temperatures rose above the freezing point, after which we measured over a 1 m increase in thaw depth for every 1 °C rise in mean annual air temperature. Under no projected warming scenario was permafrost remaining in the upper 3 m of soil by 2100. We demonstrated an effective data assimilation method that optimizes parameterization of a soil thermodynamic model. The sensitivity of local permafrost to climate warming illustrates the vulnerability of sub-Arctic tundra ecosystems to significant and rapid soil thawing.

54 ENVIRONMENTAL SCIENCES↗

Soft particles and infinite-dimensional geometry

Abstract In the sigma model, soft insertions of moduli scalars enact parallel transport ofS-matrix elements about the finite-dimensional moduli space of vacua, and the antisymmetric double-soft theorem calculates the curvature of the vacuum manifold. We explore the analogs of these statements in gauge theory and gravity in asymptotically flat spacetimes, where the relevant moduli spaces are infinite-dimensional. These models have spaces of vacua parameterized by (trivial) flat connections on the celestial sphere, and soft insertions of photons, gluons, and gravitons parallel transportS-matrix elements about these infinite-dimensional manifolds. We argue that the antisymmetric double-soft gluon theorem ind + 2 bulk dimensions computes the curvature of a connection on the infinite-dimensional space Map ( S d , G ) / G , whereGis the global part of the gauge group. The analogous metrics in abelian gauge theory and gravity are flat, as indicated by the vanishing of the antisymmetric double-soft theorems in those models. In other words, Feynman diagram calculations not only know about the vacuum manifold of Yang–Mills theory, they can also be used to compute its curvature. The results have interesting implications for flat space holography.

Astronomy & Astrophysics↗

From material properties to multiscale modeling to improve lithium-ion energy storage safety

Energy storage using lithium-ion cells dominates consumer electronics and is rapidly becoming predominant in electric vehicles and grid-scale energy storage, but the high energy densities attained lead to the potential for release of this stored chemical energy. This article introduces some of the paths by which this energy might be unintentionally released, relating cell material properties to the physical processes associated with this potential release. The selected paths focus on the anode–electrolyte and cathode–electrolyte interactions that are of typical concern for current and near-future systems. Relevant material processes include bulk phase transformations, bulk diffusion, surface reactions, transport limitations across insulating passivation layers, and the potential for more complex material structures to enhance safety. We also discuss the development, parameterization, and application of predictive models for this energy release and give examples of the application of these models to gain further insight into the development of safer energy storage systems.

25 ENERGY STORAGE↗

Toward the Development of a Cold Regions Regional-Scale Hydrologic Model, Final Project Report

This project improves meso-scale hydrologic modeling in the boreal forest by: (1) demonstrating the importance of capturing the heterogeneity of the landscape using small scale datasets for parameterization for both small and large basins; (2) demonstrating that in drier parts of the landscape and as the boreal forest dries with climate change, modeling approaches must consider the sensitivity of simulations to soil hydraulic parameters - such as residual water content - that are usually held constant. Thus, variability / flexibility in residual water content must be considered for accurate simulation of hydrologic processes in the boreal forest; (3) demonstrating that assessing climate change impacts on boreal forest hydrology through multiple model integration must account for direct effects of climate change (temperature and precipitation), and indirect effects from climate impacts on landscape characteristics (permafrost and vegetation distribution). Simulations demonstrated that climate change will increase runoff, but will increase ET to a greater extent and result in a drying of the landscape; and (4) vegetation plays a significant role in boreal hydrologic processes in permafrost free areas that have deciduous trees. This landscape type results in a decoupling of ET and precipitation, a tight coupling of ET and temperature, low runoff, and overall soil drying.

54 ENVIRONMENTAL SCIENCES↗

SU-C-202-04: Adapting Biologically Optimized Dose Escalation Based On Mid-Treatment PET/CT for Non-Small-Cell Lung Cancer

Purpose: To develop a biological modeling strategy which incorporates the response observed on the mid-treatment PET/CT into a dose escalation design for adaptive radiotherapy of non-small-cell lung cancer. Method: FDG-PET/CT was acquired midway through standard fractionated treatment and registered to pre-treatment planning PET/CT to evaluate radiation response of lung cancer. Each mid-treatment PET voxel was assigned the median SUV inside a concentric 1cm-diameter sphere to account for registration and imaging uncertainties. For each voxel, the planned radiation dose, pre- and mid-treatment SUVs were used to parameterize the linear-quadratic model, which was then utilized to predict the SUV distribution after the full prescribed dose. Voxels with predicted post-treatment SUV≥2 were identified as the resistant target (response arm). An adaptive simultaneous integrated boost was designed to escalate dose to the resistant target as high as possible, while keeping prescription dose to the original target and lung toxicity intact. In contrast, an adaptive target volume was delineated based only on the intensity of mid-treatment PET/CT (intensity arm), and a similar adaptive boost plan was optimized. The dose escalation capability of the two approaches was compared. Result: Images of three patients were used in this planning study. For one patient, SUV prediction indicated complete response and no necessary dose escalation. For the other two, resistant targets defined in the response arm were multifocal, and on average accounted for 25% of the pre-treatment target, compared to 67% in the intensity arm. The smaller response arm targets led to a 6Gy higher mean target dose in the adaptive escalation design. Conclusion: This pilot study suggests that adaptive dose escalation to a biologically resistant target predicted from a pre- and mid-treatment PET/CT may be more effective than escalation based on the mid-treatment PET/CT alone. More plans and ultimately clinical protocols are needed to validate this approach. MSKCC has a research agreement with Varian Medical System.

60 APPLIED LIFE SCIENCES↗

Modeling the Winter Heat Conduction Through the Sea Ice System During MOSAiC

Abstract Models struggle to accurately simulate observed sea ice thickness changes, which could be partially due to inadequate representation of thermodynamic processes. We analyzed co‐located winter observations of the Arctic sea ice from the Multidisciplinary Drifting Observatory for the Study of the Arctic Climate for evaluating and improving thermodynamic processes in sea ice models, aiming to enable more accurate predictions of the warming climate system. We model the sea ice and snow heat conduction for observed transects forced by realistic boundary conditions to understand the impact of the non‐resolved meter‐scale snow and sea ice thickness heterogeneity on horizontal heat conduction. Neglecting horizontal processes causes underestimating the conductive heat flux of 10% or more. Furthermore, comparing model results to independent temperature observations reveals a ∼5 K surface temperature overestimation over ice thinner than 1 m, attributed to shortcomings in parameterizing surface turbulent and radiative fluxes rather than the conduction. Assessing the model deficiencies and parameterizing these unresolved processes is required for improved sea ice representation.

Geology↗

Atomistic measurement and modeling of intrinsic fracture toughness of two-dimensional materials

Quantifying the intrinsic mechanical properties of two-dimensional (2D) materials is essential to predict the long-term reliability of materials and systems in emerging applications ranging from energy to health to next-generation sensors and electronics. Currently, measurements of fracture toughness and identification of associated atomistic mechanisms remain challenging. Herein, we report an integrated experimental–computational framework in which in-situ high-resolution transmission electron microscopy (HRTEM) measurements of the intrinsic fracture energy of monolayer MoS 2 and MoSe 2 are in good agreement with atomistic model predictions based on an accurately parameterized interatomic potential. Changes in crystalline structures at the crack tip and crack edges, as observed in in-situ HRTEM crack extension tests, are properly predicted. Such a good agreement is the result of including large deformation pathways and phase transitions in the parameterization of the inter-atomic potential. The established framework emerges as a robust approach to determine the predictive capabilities of molecular dynamics models employed in the screening of 2D materials, in the spirit of the materials genome initiative. Moreover, it enables device-level predictions with superior accuracy (e.g., fatigue lifetime predictions of electro- and opto-electronic nanodevices).

Zhang, Xu↗

Distributary Channel Networks as Moving Boundaries: Causes and Morphodynamic Effects

In this study, we propose an exploratory model to describe the morphodynamics of distributary channel network growth on river deltas. The interface between deep channels and the shallow, unchannelized delta front deposits is modeled as a moving boundary. Steady flow over the unchannelized delta front is friction dominated and modeled by Laplace's equation. Shear stress along the network boundary produces nonlinear erosion rates at the interface, causing the boundary to move and network elements (channels and branches) to form. The model was run for boundary conditions resembling the Wax Lake Delta in coastal Louisiana, 20 parameterizations of sediment transport, and 3 parameterizations of discharge. In each case, the model produced a complex channel network with channel number, width, bifurcation angle, and channel shape depending on the sediment transport formula. For reasonable sediment transport parameters and gradually increasing water discharge, the model produced network characteristics and progradation rates similar to theWax Lake Delta. Lastly, this suggests that the model contains the processes responsible for network growth, despite its abstract formulation.

58 GEOSCIENCES↗

What Polarimetric Weather Radars Offer to Cloud Modelers: Forward Radar Operators and Microphysical/Thermodynamic Retrievals

The utilization of polarimetric weather radars for optimizing cloud models is a next frontier of research. It is widely understood that inadequacies in microphysical parameterization schemes in numerical weather prediction (NWP) models is a primary cause of forecast uncertainties. Due to its ability to distinguish between hydrometeors with different microphysical habits and to identify “polarimetric fingerprints” of various microphysical processes, polarimetric radar emerges as a primary source of needed information. There are two approaches to leverage this information for NWP models: (1) radar microphysical and thermodynamic retrievals and (2) forward radar operators for converting the model outputs into the fields of polarimetric radar variables. In this paper, we will provide an overview of both. Polarimetric measurements can be combined with cloud models of varying complexity, including ones with bulk and spectral bin microphysics, as well as simplified Lagrangian models focused on a particular microphysical process. Combining polarimetric measurements with cloud modeling can reveal the impact of important microphysical agents such as aerosols or supercooled cloud water invisible to the radar on cloud and precipitation formation. Some pertinent results obtained from models with spectral bin microphysics, including the Hebrew University cloud model (HUCM) and 1D models of melting hail and snow coupled with the NSSL forward radar operator, are illustrated in the paper.

54 ENVIRONMENTAL SCIENCES↗

Calibration of cloud and aerosol related parameters for solar irradiance forecasts in WRF-solar

Model parameters are a major source of uncertainty in numerical weather prediction. Recently, the Weather Research and Forecasting model with Solar extensions (WRF-Solar) has been upgraded by enhancing the treatment of sub-grid scale cloud and aerosols with augmentations of a sub-grid scale cloud scheme (CLD3) and an upgraded aerosol-aware Thompson-Eidhammer scheme (TE14). However, the value of model parameters associated with these parameterizations are assigned based on limited measurements or theoretical calculations. Calibrating the most sensitive parameters has the potential to improve solar irradiance predictions. Here, we adopted a multiobjective surrogate-based optimization (SBO) framework to calibrate nine parameters used in CLD3 and TE14 that lead to the largest sensitivity in simulated irradiance. The normalized mean-absolute-error (NMAE) of global horizontal irradiance (GHI) and direct normal irradiance (DNI) are minimized by calibrating WRF-Solar over two regions including the Southern Great Plains (SGP) and Central California, in order to focus on parameter calibration under cloudy conditions with different aerosol loading. The results show that generalized linear model (GLM)-based surrogate models approximate physical models well, particularly when the third order and three-way interaction terms are considered. The SBO framework efficiently searches the parameter space for optimal solutions with less computational costs than directly calibrating the physical model. We first calibrate CLD3 parameters over the less-polluted SGP region. Optimized CLD3 parameters alone result in NMAE reduction by 14% for the site-mean and up to 33% for individual cases over the SGP region. With further calibration of TE14 parameters over the Central California during active fire periods, the optimized parameters lead to over 20% reductions of NMAE. Our investigation reveals, however, that optimizing TE14 has a limited impact on irradiance simulations under less-polluted conditions in the SGP.

14 SOLAR ENERGY↗

Model Structure and Climate Data Uncertainty in Historical Simulations of the Terrestrial Carbon Cycle (1850–2014)

Abstract The divergence among Earth system models in the terrestrial carbon cycle has prompted interest in how to reduce uncertainty. Previous studies have identified model structural uncertainty arising from process parameterizations and parameter values. The current study highlights the importance of climate forcing in generating carbon cycle uncertainty. We use simulations in which three models (Community Land Model version 4 (CLM4), CLM4.5, CLM5) with substantially different carbon cycles are forced with two climate reconstructions (CRUNCEPv7, Global Soil Wetness Project 3 version 1 (GSWP3v1)) to examine the contributions of model structure and climate to uncertainty in the carbon cycle over the period 1850–2014. Climate uncertainty for global annual net biome production exceeds one third of total uncertainty (defined as the sum of climate and model structure uncertainty) in the first half of the twentieth century, but declines after the 1950s. Global annual gross primary productivity, net primary productivity, heterotrophic respiration, and vegetation and soil carbon stocks have substantial climate uncertainty (relative to total uncertainty) throughout the simulation period. Climate forcing contributes more than one half of total uncertainty for these carbon cycle fluxes and stocks throughout boreal North America and Eurasia, some midlatitude regions, and in eastern Amazonia and western equatorial Africa during the decade 2000–2009. Comparison with observationally based data sets of the carbon cycle using model benchmarking methods provides insight into strengths and deficiencies among models and climate forcings, but we caution against overreliance on benchmarking to discriminate among models. The conceptualization of uncertainty arising from this study implies embracing multiple feasible model simulations rather than focusing on which model or simulation is best.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The ellipticity parameterization for an NFW profile: An overlooked angular structure in strong lens modeling

Galaxy-scale gravitational lenses are often modeled with two-component mass profiles where one component represents the stellar mass and the second is a Navarro Frenk White (NFW) profile representing the dark matter. Outside of the spherical case, the NFW profile is costly to implement, and so it is approximated via two different methods; ellipticity can be introduced via the lensing potential (NFWp) or via the mass by approximating the NFW profile as a sum of analytical profiles (NFWm). While the NFWp method has been the default for lensing applications, it gives a different prescription of the azimuthal structure, which we show introduces ubiquitous gradients in ellipticity and boxiness in the mass distribution rather than having a constant elliptical shape. Because an unmodeled azimuthal structure has been shown to be able to bias lens model results, we explored the degree to which this azimuthal structure that was introduced can affect the model accuracy. We constructed input profiles using composite models using both the NFWp and NFWm methods and fit these mocks with a power-law elliptical mass distribution (PEMD) model with external shear. As a measure of the accuracy of the recovered lensing potential, we calculated the value of the Hubble parameter H 0 one would determine from the lensing fit. We found that the fits to the NFWp input return H 0 values that are systematically biased by about 3% lower than the NFWm counterparts. We explored whether such an effect is attributable to the mass sheet transformation (MST) by using an MST-independent quantity, ξ 2 . We show that, as expected, the NFWm mocks are degenerate with PEMD through an MST. For the NFWp, an additional bias was found beyond the MST due to the azimuthal structure exterior to the Einstein radius. We recommend modelers use an NFWm prescription in the future, such that the azimuthal structure can be introduced explicitly rather than implicitly.

79 ASTRONOMY AND ASTROPHYSICS↗