SEARCH · Engineering Papers

Results for “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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,567 records · Page 87

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Biocrust contribution to ecosystem carbon fluxes varies along an elevational gradient

Abstract Understanding how each component of an ecosystem contributes to carbon fluxes across a range of abiotic conditions enables accurate forecasts for future emission scenarios. In drylands, biological soil crust (biocrust) contribution to ecosystem carbon fluxes may vary at a regional scale but is rarely quantified due to the difficulty of parameterizing process‐based models or parsing biocrust (on the soil surface) from whole‐soil flux measurements. We measured excised biocrust carbon fluxes across current and future predicted warmer summertime temperatures from dryland sites (grassland, shrubland, savanna, and woodland) and compared those to the ecosystem and soil fluxes from eddy flux towers. Overall, biocrust fluxes showed net carbon loss during the warm growing season temperatures, suggesting that cool‐season photosynthesis is likely critical for maintaining positive biocrust carbon balance in these sites. Biocrust flux temperature responses differed by site: Grassland and shrubland biocrust gross photosynthesis was relatively invariant, while respiration increased with temperature; in the woodland and savanna, biocrust gross photosynthesis and respiration increased with temperature. Biocrust fluxes contributed <3% to the observed ecosystem gross primary productivity at 26°C and ecosystem respiration at 19°C in grassland, shrubland, and woodland, but contributed >10% of observed soil respiration in grassland and shrubland at 19°C, reinforcing the need to separate biocrust activity from root and subsurface heterotroph activity to understand drivers of fluxes at different sites. Regional resolution of biocrust type and cover will improve predictions of biocrust contribution to global carbon flux with changing temperatures.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Establishing connectivity patterns of eastern oysters ( Crassostrea virginica ) on regional oceanographic scales

Localized depletion of marine and estuarine populations often results from large-scale natural and anthropogenic disturbances (e.g., hurricanes, oil spills) as well as overharvest of fisheries resources. Understanding how such localized depletions may affect populations at larger regional scales requires knowledge of connectivity among local populations within the larger regional landscape or metapopulation. Efforts to restore populations following such decreases require similar knowledge. During the 2010 Deepwater Horizon oil spill, dramatic declines of oysters occurred throughout estuaries west and east of the Mississippi River. We examined trajectories of particles, which were parameterized to mimic oyster larvae, using the ADvanced CIRCulation (ADCIRC; https://adcirc.org ) model to evaluate potential connectivity within and among embayments from Western Louisiana to Alabama. Patterns of larval settlement, which we defined as the intersection of a larval particle with known or expected oyster habitat at any point 13–21 days post-release, reflected much greater local contributions, with 50%–90% of settled particles originating within the same subdivision of embayments. Exchange among subdivision was much less (0%–40%) and settlement originating from outside the embayment of release was trivial under most scenarios (0%–14%). Connectivity between adjacent basins was greatest for larvae released in the southern portions of the embayment, whereas connectivity among nonadjacent basins was not predicted under the scenarios modeled. Because most local populations are relatively isolated on ecological timescales, areas suffering from extensive local depletion are likely to require extensive time to recover due to the lack of larval subsidy from the overall regional population. Restoration would require building stepping-stone populations or reefs within and among basin to restore a high degree of connectivity.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Simulation-based characterization of the variability of earthquake risk to buildings in the near-field

Recent advancements in high performance computing platforms and computational workflow for regional-scale simulations are enabling unprecedented modeling of fault-to-structure earthquake processes. Regional simulations resolving ground motions at frequencies relevant to engineered systems are becoming computationally viable and provide a new capability to improve understanding of the geographical distribution and intensity of risk to buildings and critical infrastructure. As computational capabilities advance, it is essential to move beyond illustrative single rupture realizations for scenario earthquake events towards the development of a full suite of rupture realizations that appropriately characterize the range of risk to building systems. The work described in this article investigates the application of a suite of fault rupture realizations with the objective of assessing near-fault, site-specific seismic demand variability for building structures. A representative high-performance regional-scale computational model is utilized to execute ground motion and building response simulations based on 18 kinematic rupture realizations of an M7 strike-slip scenario earthquake. The fault rupture models for the scenario earthquake are created by systematically perturbing the hypocenter location and stochastically generating rupture parameters (slip, rise time, rake angle) to represent a breadth of ground motion intensities resulting from the spatial and temporal variabilities of an earthquake rupture process. The resulting seismic demand variability for three-story (short period) and forty-story (long period) steel moment-resisting frame buildings is characterized in terms of the median and distribution of peak inter-story drift ratio for a range of near-fault sites. The full suite of 18 fault rupture realizations and approximately 280,000 nonlinear dynamic building simulations indicate that the three-story building undergoes higher median seismic demand and significantly greater variability of demand at a given site than the forty-story building, which has important implications for the level of certainty in predicting building performance during an earthquake. The simulations performed provide deeper insight into the relationship between fault rupture parameterization and building response, which is essential information for developing a representative suite of rupture realizations for specific earthquake scenarios.

58 GEOSCIENCES↗

DOE OSTI2023

Connecting suborganismal data to bioenergetic processes: killifish embryos exposed to a dioxin-like compound

A core challenge for ecological risk assessment is to integrate molecular responses into a chain of causality to organismal or population level outcomes. Bioenergetic theory may be a useful approach for integrating suborganismal responses to predict organismal level responses that influence population dynamics. In this work, we describe a novel application of Dynamic Energy Budget (DEB), theory in the context of a toxicity framework (Adverse Outcome Pathways, AOP) to make quantitative predictions of chemical exposures to individuals, starting from suborganismal data. We use early life stage exposure of Fundulus heteroclitus to dioxin-like chemicals (DLCs) and connect AOP Key Events (KEs) to DEB processes through “damage” that is produced at a rate proportional to the internal toxicant concentration. We use transcriptomic data of fish embryos exposed to DLCs to translate molecular indicators of damage into changes in DEB parameters (damage increases somatic maintenance costs) and use DEB models to predict sublethal and lethal effects of young fish. By changing a small subset of model parameters, we predict the evolved tolerance to DLCs in some wild F. heteroclitus populations, a data set not used in model parameterization. The differences in model parameters points to reduced sensitivity and altered damage repair dynamics as contributing to this evolved resistance. Our methodology has potential extrapolation to untested chemicals of ecological concern.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

DOE OSTI2025

A Data Library of Liquid Clouds Modelled With a Large Eddy Simulation Framework

We describe a library of atmospheric large eddy simulations (LES) of liquid-phase boundary layer clouds constructed to enable aerosol–cloud–turbulence interaction studies, support parameterization evaluation and development, and provide training data for machine learning applications. The simulations use a modern LES framework designed for high numerical accuracy, coupled to a detailed spectral bin microphysical scheme. Case studies are configured to represent observed conditions in four key global cloud regions—the Northeastern Atlantic, Northeastern Pacific, Continental United States and Southern Ocean—following a semi-idealised approach. The library also includes aerosol concentration halving and doubling experiments to expose the sensitivities of the case studies to aerosol perturbations. Simulation results are compared to observations on a case-by-case basis, then the library's coverage is evaluated in terms of spreads in meteorological factors and atmospheric boundary layer attributes.

aerosol↗

DOE OSTI2017

Evaluation of energy balance closure adjustment methods by independent evapotranspiration estimates from lysimeters and hydrological simulations

Here, we report non-closure of the surface energy balance is a frequently observed phenomenon of hydrometeorological field measurements, when using the eddy-covariance method, which can be ascribed to an underestimation of the turbulent fluxes. Several approaches have been proposed in order to adjust the measured fluxes for this apparent systematic error. However, there are uncertainties about partitioning of the energy balance residual between the sensible and latent heat flux and whether such a correction should be applied on 30-minute data or longer time scales. The data for this study originate from two grassland sites in southern Germany, where measurements from weighable lysimeters are available as reference. The adjusted evapotranspiration rates are also compared with joint energy and water balance simulations using a physically-based distributed hydrological model. We evaluate two adjustment methods: the first one preserves the Bowen ratio and the correction factor is determined on a daily basis. The second one attributes a smaller portion of the residual energy to the latent heat flux than to the sensible heat flux for closing the energy balance for every 30-minute flux integration interval. Both methods lead to an improved agreement of the eddy-covariance based fluxes with the independent lysimeter estimates and the physically-based model simulations. The first method results in a better comparability of evapotranspiration rates, and the second method leads to a smaller overall bias. These results are similar between both sites despite considerable differences in terrain complexity and grassland management. Moreover, we found that a daily adjustment factor leads to less scatter than a complete partitioning of the residual for every half-hour time interval. Lastly, the vertical temperature gradient in the surface layer and friction velocity were identified as important predictors for a potential future parameterization of the energy balance residual.

58 GEOSCIENCES↗

DOE OSTI2023

The importance of explicitly representing the streambed in watershed models

Abstract The streambed is the critical interface between the aquatic and terrestrial systems and hosts important biogeochemical hot spots within river corridors. Although the streambed characteristics are significantly different from those of its surrounding soil, the streambed itself has not been explicitly represented in watershed models. Here, we explicitly incorporated a streambed layer into an integrated hydrologic model through model parameterization and discretization. We examined the hydrological effects of streambed characteristics, including hydraulic conductivity ( K ), layer thickness, and resolution, on the exchange fluxes across the streambed as well as the streamflow at the watershed outlet. The numerical experiments were performed in the American River Watershed, a headwater, mountainous watershed within the Yakima River Basin in central Washington. Despite having a negligible effect on the watershed streamflow, an explicit representation of the streambed with distinctive properties dramatically changed the magnitude and variability of the exchange flux. In general, a larger streambed K along with a thicker streambed layer induced larger exchange fluxes. The exchange flux was most sensitive to the streambed resolution. A finer streambed resolution increased exchange fluxes per unit area while reducing the overall exchange volumes across the entire streambed. The amount of baseflow decreased by 6% as the streambed resolution increased from 250 to 50 m. This finding is important because these hydrological changes may, in turn, affect the exchange of nutrients and contaminants between surface water and groundwater and the associated biogeochemical processes. Our work demonstrated the importance of representing streambeds in fully distributed, process‐based watershed models to better capture the exchange flow dynamics in river corridors.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Modelling the Effects of Wetland Restoration on Coastal Hydrology: A Case Study of Elkhorn Slough Watershed, California

ABSTRACT Coastal wetlands, some of the most productive ecosystems on Earth, provide critical ecosystem services, including support of biodiversity, carbon sequestration and flood protection. In recent decades, these ecosystems have experienced extensive coastal wetland loss. Coastal wetland restoration provides a beacon of hope, offering a chance to reclaim these important habitats. However, even with billions of dollars invested worldwide in restoring coastal wetlands, we still lack comprehensive knowledge about the effectiveness of these restoration efforts in recovering wetland ecosystem functions and how future climate change may affect these efforts. The ability to evaluate how these ecosystems will function in the future is vital for examining current investments and developing future protection and management plans. We selected Elkhorn Slough, a tidal estuary, in California, to investigate the impact of wetland restoration and sea level rise (SLR) on coastal hydrology using the process‐based coastal hydrologic model, Advanced Terrestrial Simulator (ATS), informed by site‐specific data. We designed a novel modelling workflow for incorporating wetland restoration features into land cover and soil properties for the model parameterization. The validation results demonstrate a strong agreement between modelled and observed data. We studied the characteristics of coastal watershed hydrology, then focused on the surface water dynamics at two wetland sites within Elkhorn Slough, a reference site and a restored site. Our simulation results indicate that the restored site successfully maintains surface elevation, resulting in reduced surface inundation. We also examined the impact of wetland restoration under expected SLR over the next few decades. The low‐lying Yampah Marsh, the reference site, is likely to be inundated due to future SLR when highest tides arrive, while a higher percentage of Hester Marsh, the restored site, would retain marsh vegetation in coming decades, regardless of tidal conditions. Our study provides important information for examining the outcome of restoration practices that include surface elevation in tidal wetlands under climate changes.

advanced terrestrial simulator↗

DOE OSTI2022

Estimating dermal contact soil exposure for amphibians

Abstract Chemical exposure estimation through the dermal route is an underemphasized area of ecological risk assessment for terrestrial animals. Currently, there are efforts to create exposure models to estimate doses from this pathway for use in ecological risk assessment. One significant limitation has been insufficient published data to characterize exposure and to support the selection and parameterization of appropriate models, particularly for amphibians in terrestrial habitats. Recent publications measuring pesticide doses to terrestrial‐phase amphibians have begun to rectify this situation. We collated and summarized available measurements of terrestrial amphibian dermal exposure to pesticides from 11 studies in which researchers measured tissue concentrations associated with known pesticide experimental application rates. This data set included tissue concentrations in 11 amphibian species and 14 different pesticides. We then compared the results of two screening exposure models that differed based on surface area scaling approaches as a function of body weight (one based on birds as surrogates for amphibians and another amphibian‐specific) to the measured tissue residue concentrations. We define a false‐negative rate for each screening model as the proportion of amphibians for which the predicted concentration is less than the observed concentration (i.e., underestimate), contrary to the intent of screening models, which are intended to have a bias for higher exposure concentrations. The screening model that uses birds as surrogates did not have any instances where estimated expected avian doses were less than measured amphibian body burdens. When using the amphibian‐specific exposure model that corrected for differences between avian and amphibian surface area, measured concentrations were greater than model estimates for 11.3% of the 1158 comparisons. The database of measured pesticide concentrations in terrestrial amphibians is provided for use in calculating bioconcentration factors and for future amphibian dermal exposure model development. Integr Environ Assess Manag 2023;19:9–16. © 2022 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Environmental Sciences & Ecology↗

DOE OSTI2023

Number of sites-based solver for determining coverages from steady-state mean-field micro-kinetic models

Kinetic models parameterized by ab-initio calculations have led to significant improvements in understanding chemical reactions in heterogeneous catalysis. These studies have been facilitated by implementations which determine steady-state coverages and rates of mean-field micro-kinetic models. As implemented in the open-source kinetic modeling program, CatMAP, the conventional solution strategy is to use a root-finding algorithm to determine the coverage of all intermediates through the steady-state expressions, constraining all coverages to be non-negative and to properly sum to unity. Though intuitive, this root-finding strategy causes issues with convergence to solution due to these imposed constraints. In this work, we avoid explicitly imposing these constraints, solving the mean-field steady-state micro-kinetic model in the space of number of sites instead of solving it in the space of coverages. Here, we transform the constrained root-finding problem to an unconstrained least-squares minimization problem, leading to significantly improved convergence in solving micro-kinetic models and thus enabling the efficient study of more complex catalytic reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2013

The Mechanism of First Raindrops Formation in Deep Convective Clouds

The formation of first raindrops in deep convective clouds is investigated. A combination of observational data analysis and 2-D and 3-D numerical bin microphysical simulations of deep convective clouds suggests that the first raindrops form at the top of undiluted or slightly diluted cores. It is shown that droplet size distributions in these regions are wider and contain more large droplets than in diluted volumes. The results of the study indicate that the initial raindrop formation is determined by the basic microphysical processes within ascending adiabatic volumes. It allows one to predict the height of the formation of first raindrops considering the processes of nucleation, diffusion growth and collisions. The results obtained in the study explain observational results reported by Freud and Rosenfeld (2012) according to which the height of first raindrop formation depends linearly on the droplet number concentration at cloud base. The results also explain why a simple adiabatic parcel model can reproduce this dependence. The present study provides a physical basis for retrieval algorithms of cloud microphysical properties and aerosol properties using satellites proposed by Rosenfeld et al. ( 2012). The study indicates that the role of mixing and entrainment in the formation of the first raindrops is not of crucial importance. It is also shown that low variability of effective and mean volume radii along horizontal traverses, as regularly observed by in situ measurements, can be simulated by high-resolution cloud models, in which mixing is parameterized by a traditional 1.5 order turbulence closure scheme.

clouds↗

DOE OSTI2013

Spatial and temporal variations of new particle formation in East Asia using an NPF-explicit WRF-chem model: North-south contrast in new particle formation frequency

The new particle formation (NPF)-explicit version of the WRF-chem model, which we developed recently, can calculate the growth and sink of nucleated clusters explicitly with 20 aerosol size bins from 1 nm to 10 μm. In this study, the model is used to understand spatial and temporal variations of the frequency of NPF events and the concentrations of aerosols (condensation nuclei, CN) and cloud condensation nuclei (CCN) within the boundary layer in East Asia in spring 2009. Model simulations show distinct north-south contrast in the frequency and mechanism of NPF in East Asia. NPF mostly occurred over limited periods and regions between 30° and 45°N, such as northeast China, Korea, and Japan, including regions around active volcanoes (Miyakejima and Sakurajima). At these latitudes, NPF was considerably suppressed by high concentrations of preexisting particles under stagnant air conditions associated with high-pressure systems, while nucleation occurred more extensively on most days during the simulation period. Conversely, neither nucleation nor NPF occurred frequently south of 30°N because of lower SO2 emissions and H2SO4 concentrations. The period-averaged NPF frequency was 3 times higher at latitudes of 30° - 45°N than at latitudes of 20° - 30°N. The north-south contrast of NPF frequency is validated by surface measurements in outflow regions in East Asia. The period- and domain-averaged contribution of secondary particles is estimated to be 44% for CN (> 10 nm) and 26% for CCN at a supersaturation of 1.0% in our simulation, though the contribution is highly sensitive to the magnitudes and size distributions of primary aerosol emissions and the coefficients in the nucleation parameterizations.

Matsui, H.↗

DOE OSTI2017

Projected precipitation changes within the Great Lakes and Western Lake Erie Basin: a multi‐model analysis of intensity and seasonality

ABSTRACT The Great Lakes region encompasses the largest freshwater lake network in the world and supports a diverse network of agriculture, transportation, and tourism. Recently, Lake Erie has experienced increased hypoxia events, which have been attributed to agricultural practices and changes in run‐off. Here we examine the projected changes in extreme precipitation events to address concerns regarding regional agriculture, surface run‐off, and subsequent water quality. Precipitation projections within the overall Great Lakes Basin and the Western Lake Erie Basin subregion are examined using climate model simulations of varying spatial resolutions to understand historical precipitation and projected future precipitation. We develop three model ensembles for the historical period (1980–1999) and the mid‐century (2041–2060) that cover a range of spatial resolutions and future emissions scenarios, including: (1) 12 global model members from the fifth Climate Model Intercomparison Project ( CMIP5 ) using Representative Concentration Pathway ( RCP ) 8.5, (2) ten regional climate model ( RCM ) members from the North American Regional Climate Change Assessment Program driven by CMIP3 global models using the A2 emissions scenario, and (3) two high resolution RCM simulations ( RCM4 ) driven by CMIP5 global models using the RCP 8.5 scenario. For the historical period, all model ensembles overestimate winter and spring precipitation, and many of the models simulate a summer drying that is not observed. At mid‐century, most of the models predict a 10–20% increase in precipitation depending on the time of year. Daily probability distribution functions from three model ensembles reveal spring seasonal increases in high precipitation event probabilities when compared to the historical period, suggesting an increase in the frequency of high intensity precipitation at mid‐century. Overall, the presence of lakes or higher spatial resolution does not ensure improved representation of historical processes, and more complex interactions between large‐scale dynamics, local feedbacks, and physical parameterizations drive the model spread.

Basile, Samantha J.↗

DOE OSTI2020

Plant‐mediated methane transport in emergent and floating‐leaved species of a temperate freshwater mineral‐soil wetland

Abstract Methane flux from freshwater mineral‐soil (FWMS) wetlands and its variability among sites is largely modulated by plant‐mediated transport. However, plant‐mediated transport processes are rarely resolved in land surface models and are poorly parametrized for plants commonly found in FWMS wetlands. Here, relationships between methane flux and CO 2 uptake, as well as plant conductance of methane were evaluated for three plant species and two characteristic functional types: emergent (narrow‐leaved cattail) and floating‐leaved (American lotus and water lily). We found significant but contrasting correlations between methane flux and CO 2 uptake in cattails ( r 2 = 0.51, slope = −0.16, during morning) and water lily ( r 2 = 0.32, slope = 0.064, after midday). This relationship was not significant in American lotus, showing that stomata regulation of methane fluxes is species‐specific and not generalizable across the floating‐leaved plant functional type. Conductance of methane per leaf area showed distinct seasonal dynamics across species. Conductance was similar among the floating‐leaved species (6.2 × 10 −3 m d −1 in lotus and 7.2 × 10 −3 m d −1 in water lily) and higher than conductance in the emergent species (2.7 × 10 −3 m d −1 ). Our results provide direct observations of plant conductance rates and identify the vegetation parameters (leaf area, stomatal conductance) that modify them. Our results further suggest that models of methane emissions from FWMS should parameterize plant‐mediated transport in different plant functional types, scaled by leaf area and with variable seasonal phenological dynamics, and consider possible species‐specific mechanisms that control methane transport through plants.

Villa, Jorge A.↗

DOE OSTI2025

Demonstrating Hierarchical System Development With the Common Community Physics Package Single‐Column Model: A Case Study Over the Southern Great Plains

This study demonstrates a specific application of the hierarchical system development (HSD) approach to investigate, analyze, and attribute model issues within the Unified Forecast System (UFS), with a focus on process isolation. By evaluating a non‐precipitating, shallow cumulus case at the Atmospheric Radiation Measurement Southern Great Plains site in the UFS global forecast against the observation, the investigation identifies a warmer and deeper daytime convective planetary boundary layer (PBL) and misrepresented nocturnal PBL transition. Hypothesis testing, which employs the Common Community Physics Package (CCPP) single‐column model (SCM) and uses the same physics as the UFS global model, confirms that these issues are attributed to the model physics and initialization. Specifically, misrepresented PBL processes are linked to problematic surface condition and a lack of cloud formation, which may stem from deficiencies in PBL and cloud microphysics parameterizations and their interactions. The UFS initial condition contributes to an earlier, excessively collapsed daytime convective boundary layer and a lack of decoupling between the stable boundary layer and residual layer late in the afternoon. This work introduces an avenue for the community to engage with the application of HSD, along with the CCPP and CCPP SCM, to understand the interplay of model physics, disentangle the roles of model components, as well as facilitate model and forecast improvement.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Methyl α‐D‐galactopyranosyl‐(1→3)‐β‐D‐galactopyranoside and methyl β‐D‐galactopyranosyl‐(1→3)‐β‐D‐galactopyranoside: Glycosidic linkage conformation determined from MA'AT analysis

Abstract MA'ATanalysis has been applied to two biologically‐importantO‐glycosidic linkages in two disaccharides, α‐D‐Galp‐(1→3)‐β‐D‐GalpOMe (3) and β‐D‐Galp‐(1→3)‐β‐D‐GalpOMe (4). Using density functional theory (DFT) to obtain parameterized equations relating a group of trans‐O‐glycosidic NMR spin‐couplings to eitherphi(ϕ') orpsi(ψ'), and experimental 3 J COCH , 2 J COC , and 3 J COCC spin‐couplings measured in aqueous solution in 13 C‐labeled isotopomers, probability distributions ofϕ'andψ'in each linkage were determined and compared to those determined by aqueous 1‐μs molecular dynamics (MD) simulation. Good agreement was found between theMA'ATand single‐state MD conformational models of these linkages for the most part, with modest (approximately <15°) differences in the mean values ofϕ'andψ', although the envelope of allowed angles (encoded in circular standard deviations or CSDs) is consistently larger forϕ'determined fromMA'ATanalysis than from MD for both linkages. TheMA'ATmodel of the α‐Galp‐(1→3)‐β‐Galplinkage agrees well with those determined previously using conventional NMR methods ( 3 J COCH values and/or 1 H‐ 1 H NOEs), but some discrepancy was observed for the β‐Galp‐(1→3)‐β‐Galplinkage, which may arise from errors in the conventions used to describe the linkage torsion angles. Statistical analyses of X‐ray crystal structures show ranges ofϕ'andψ'for both linkages that include the mean angles determined fromMA'ATanalyses, although both angles adopt a wide range of values in the crystalline state, withϕ'in β‐Galp‐(1→3)‐β‐Galplinkages showing greater‐than‐expected conformational variability.

Chemistry↗

DOE OSTI2014

Adaptive h -refinement for reduced-order models: ADAPTIVE h -refinement for reduced-order models

Our work presents a method to adaptively refine reduced-order models a posteriori without requiring additional full-order-model solves. The technique is analogous to mesh-adaptive h-refinement: it enriches the reduced-basis space online by ‘splitting’ a given basis vector into several vectors with disjoint support. The splitting scheme is defined by a tree structure constructed offline via recursive k-means clustering of the state variables using snapshot data. This method identifies the vectors to split online using a dual-weighted-residual approach that aims to reduce error in an output quantity of interest. The resulting method generates a hierarchy of subspaces online without requiring large-scale operations or full-order-model solves. Furthermore, it enables the reduced-order model to satisfy any prescribed error tolerance regardless of its original fidelity, as a completely refined reduced-order model is mathematically equivalent to the original full-order model. Experiments on a parameterized inviscid Burgers equation highlight the ability of the method to capture phenomena (e.g., moving shocks) not contained in the span of the original reduced basis.

42 ENGINEERING↗

DOE OSTI2016

Exploration of efficient reduced-order modeling and a posteriori error estimation

Efficient algorithms are considered for the computation of a reduced-order model based on the proper orthogonal decomposition methodology for the solution of parameterized elliptic partial differential equations. The method relies on partitioning the parameter space into subdomains based on the properties of the solution space and then forming a reduced basis for each of the subdomains. This yields more efficient offline and online stages for the proper orthogonal decomposition method. We extend these ideas for inexpensive adjoint based a posteriori error estimation of both the expensive finite element method solutions and the reduced-order model solutions, for a single and multiple quantities of interest. Finally, various numerical results indicate the efficacy of the approach.

42 ENGINEERING↗