Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ASSIMILATION”

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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Using observational mean-flow data to drive large-eddy simulations of a diurnal cycle at the SWiFT site

Reproducing realistic date- and site-specific unsteady wind conditions in large-eddy simulations is becoming increasingly useful in wind energy. How to run a large-eddy simulation to match observed conditions, however, remains an open research question. One approach that has received considerable attention is mesoscale-to-microscale coupling, in which information about the mesoscale weather, most commonly acquired from a mesoscale numerical weather model, is passed on to a microscale model. In this paper, we demonstrate how the recently developed profile-assimilation technique, a form of mesoscale-to-microscale coupling, can be used to drive large-eddy simulations solely based on observed mean-flow profiles at a single location, bypassing the need for auxiliary mesoscale simulations. The new approach is evaluated for a diurnal cycle at the Scaled Wind Farm Technology site. Observed mean-flow profiles from the ground up to a height of 2 km are reconstructed by aggregating measurements from multiple instruments, and gaps in the data are infilled with natural neighbor interpolation. We perform nine simulations using various forcing approaches to deal with data limitations. The results show that it is indeed possible to drive microscale large-eddy simulation with observations using the profile-assimilation technique, notwithstanding large gaps in virtual potential temperature measurements. However, profile assimilation with vertical smoothing of the error between the desired and actual profiles is required. Without that smoothing, the microscale simulations develop unrealistically high turbulence levels under many situations. Finally, we show that simulated mesoscale data can account for missing observations, although care is needed as both data sources are not necessarily compatible.

17 WIND ENERGY↗
DOE OSTI2021

Development and application of marginal likelihood optimization for integral parameter adjustment

When adjusting nuclear data with integral experiments, care must be taken that spurious adjustments are not made by assimilating poorly characterized integral parameters. If there are unaccounted for biases or poorly estimated uncertainties in the calculated and experimental values for an integral parameter, the Bayesian data assimilation may adjust the nuclear data in a manner that does not reflect the physics of the integral parameter. To identify and lessen the impact of these inconsistent integral parameters, in this study we present a Marginal Likelihood Optimization algorithm. In a data-driven way, the marginalized likelihood is used to modulate hyperparameter terms that decrease the influence of inconsistent integral parameters on the adjustment. The advantage of this approach over other methods in the literature is that it incorporates correlation information and does not remove an integral parameter from the adjustment. Herein, we present and motivate the algorithm, and apply it to an integral data assimilation case study.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

NeuroSEM: A hybrid framework for simulating multiphysics problems by coupling PINNs and spectral elements

Multiphysics problems that are characterized by complex interactions among fluid dynamics, heat transfer, structural mechanics, and electromagnetics, are inherently challenging due to their coupled nature. While experimental data on certain state variables may be available, integrating these data with numerical solvers remains a significant challenge. Physics-informed neural networks (PINNs) have shown promising results in various engineering disciplines, particularly in handling noisy data and solving inverse problems in partial differential equations (PDEs). However, their effectiveness in forecasting nonlinear phenomena in multiphysics regimes, particularly involving turbulence, is yet to be fully established. Here, this study introduces NeuroSEM, a hybrid framework integrating PINNs with the highfidelity Spectral Element Method (SEM) solver, Nektar++. NeuroSEM leverages the strengths of both PINNs and SEM, providing robust solutions for multiphysics problems. PINNs are trained to assimilate data and model physical phenomena in specific subdomains, which are then integrated into the Nektar++ solver. We demonstrate the efficiency and accuracy of NeuroSEM for thermal convection in cavity flow and flow past a cylinder. The framework effectively handles data assimilation by addressing those subdomains and state variables where the data is available. We applied NeuroSEM to the Rayleigh-B´enard convection system, including cases with missing thermal boundary conditions and noisy datasets. Finally, we applied the proposed NeuroSEM framework to real particle image velocimetry (PIV) data to capture flow patterns characterized by horseshoe vortical structures. Our results indicate that NeuroSEM accurately models the physical phenomena and assimilates the data within the specified subdomains. The framework’s plug-and-play nature facilitates its extension to other multiphysics or multiscale problems. Furthermore, NeuroSEM is optimized for efficient execution on emerging integrated GPU-CPU architectures. This hybrid approach enhances the accuracy and efficiency of simulations, making it a powerful tool for tackling complex engineering challenges in various scientific domains.

42 ENGINEERING↗
DOE OSTI2018

Integrating population connectivity into pollution assessment: Overwintering mixing reveals flame retardant contamination in breeding areas in a migratory raptor

Highlights: • This multidisciplinary study integrates population connectivity into pollution assessment. • Flame retardants in inviable eggs of black kites were assessed in two populations. • A much higher occurrence and concentration of flame retardants were found in Madrid than Doñana. • Kites from Western Europe showed a strong winter intermixing in Africa. • Diffuse migratory connectivity reveals breeding areas as the main contamination grounds. Determining the exposure and magnitude at which various pollutants are differentially assimilated at the breeding and non-breeding grounds of migratory wildlife is challenging. Here, the possibility of applying the migratory connectivity framework to understanding contamination in birds is illustrated by considering flame retardants in inviable eggs of a migratory raptor, the black kite (Milvus migrans). The occurrence and concentration of legacy and emerging compounds in eggs from the southeastern peri-urban area of Madrid city, central Spain, were compared with those from Doñana National Park in southern Spain. A much higher occurrence and concentration of multiple polybrominated diphenyl ethers and Dechlorane 602 were found in Madrid than Doñana, but the opposite patterns were found for Dechlorane Plus. Individuals from these and other breeding areas in western Europe showed a strong intermixing pattern over widespread wintering areas in Africa, as assessed by ringing recoveries and movements tracked by satellite devices. This diffuse migratory connectivity reveals breeding areas as the main contamination grounds. High contamination burdens sequestered in eggs point to rapid assimilation of these compounds before laying, associated with important emission sources in Madrid, especially landfills of partially incinerated urban refuse, and other anthropogenic operations. Diet composition regarding aquatic vs. terrestrial prey, and bioaccumulation and biomagnification processes are suggested to explain differential assimilation of some compounds, especially Dechlorane Plus in Doñana, although a local emission source polluting this area cannot be ruled out. Insight from the migratory connectivity framework can help to disentangle large-scale patterns of contaminant uptake and refocus attention on key regions and potential causes of chemical hazards in declining migratory species and human populations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

A probabilistic drought forecasting framework: A combined dynamical and statistical approach

In order to improve drought forecasting skill, this study develops a probabilistic drought forecasting framework comprised of dynamical and statistical modeling components. The novelty of this study is to seek the use of data assimilation to quantify initial condition uncertainty with the Monte Carlo ensemble members, rather than relying entirely on the hydrologic model or land surface model to generate a single deterministic initial condition, as currently implemented in the operational drought forecasting systems. Next, the initial condition uncertainty is quantified through data assimilation and coupled with a newly developed probabilistic drought forecasting model using a copula function. The initial condition at each forecast start date are sampled from the data assimilation ensembles for forecast initialization. Finally, seasonal drought forecasting products are generated with the updated initial conditions. This study introduces the theory behind the proposed drought forecasting system, with an application in Columbia River Basin, Pacific Northwest, United States. Results from both synthetic and real case studies suggest that the proposed drought forecasting system significantly improves the seasonal drought forecasting skills and can facilitate the state drought preparation and declaration, at least three months before the official state drought declaration.

data assimilation↗
DOE OSTI2026

Hourly PM 2.5 Estimates across California from 2018 to 2023

This study presents a new data set of hourly PM 2.5 concentrations across California from 2018 to 2023 at a three-kilometer resolution. This data set was developed by assimilating observations from PurpleAir and the U.S. EPA Air Quality System monitors into wildfire smoke forecasts from the High-Resolution Rapid Refresh Smoke (HRRR-Smoke) model using the Gridpoint Statistical Interpolation (GSI) three-dimensional variational data assimilation framework. Archived forecasts of modeled wildfire smoke PM 2.5 from HRRR-Smoke create the background field for assimilation, which is then corrected using surface observations of total PM 2.5 . The resulting reanalysis from GSI provides an estimate of total PM 2.5 that minimizes error from both the observational and the model data. Validation results indicate strong performance, with monthly R 2 values ranging from 0.73 to 0.91 across the six-year data set, comparable to other PM 2.5 data sets. Case studies are presented for three major fire events, the 2018 Camp Fire, 2019 Kincade Fire, and 2020 Lightning Complex Fires to demonstrate the data set’s fidelity in resolving plume dynamics and local exposure patterns. Root-mean-squared error averaged over each month scales with average PM 2.5 concentrations, resulting in a low error under typical conditions but higher absolute errors during extreme smoke events. This is the first long-term, hourly PM 2.5 data set of its kind for California and enables the generation of subdaily exposure metrics, such as peak hourly concentrations, exceedance durations, and time-of-day exposure peaks. The novelty and strong validation of this data set make it a compelling resource for future studies on the impact and significance of subdaily PM 2.5 exposure.

PM2.5↗
DOE OSTI2023

Physics–Informed Neural Networks of the Saint–Venant Equations for Downscaling a Large–Scale River Model

Large-scale river models are being refined over coastal regions to improve the scientific understanding of coastal processes, hazards and responses to climate change. However, coarse mesh resolutions and approximations in physical representations of tidal rivers limit the performance of such models at resolving the complex flow dynamics near the river-ocean interface, resulting in inaccurate simulations of flood inundation. In this research, we propose a machine learning (ML) framework based on the state-of-the-art physics-informed neural network (PINN) to simulate the downscaled flow at the subgrid scale. First, we demonstrate that PINN is able to assimilate observations of various types and solve the one-dimensional (1-D) Saint-Venant equations (SVE) directly. We perform the flow simulations over a floodplain and along an open channel in several synthetic case studies. The PINN performance is evaluated against analytical solutions and numerical models. Our results indicate that the PINN solutions of water depth have satisfactory accuracy with limited observations assimilated. In the case of flood wave propagation induced by storm surge and tide, a new neural network architecture is proposed based on Fourier feature embeddings that seamlessly encodes the periodic tidal boundary condition in the PINN's formulation. Furthermore, we show that the PINN-based downscaling can produce more reasonable subgrid solutions of the along-channel water depth by assimilating observational data. The PINN solution outperforms the simple linear interpolation in resolving the topography and dynamic flow regimes at the subgrid scale. This study provides a promising path towards improving emulation capabilities in large-scale models to characterize fine-scale coastal processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Rapid CO2 changes cause oscillations in photosynthesis that implicate PSI acceptor-side limitations

Abstract Oscillations in CO2 assimilation rate and associated fluorescence parameters have been observed alongside the triose phosphate utilization (TPU) limitation of photosynthesis for nearly 50 years. However, the mechanics of these oscillations are poorly understood. Here we utilize the recently developed dynamic assimilation techniques (DATs) for measuring the rate of CO2 assimilation to increase our understanding of what physiological condition is required to cause oscillations. We found that TPU-limiting conditions alone were insufficient, and that plants must enter TPU limitation quickly to cause oscillations. We found that ramps of CO2 caused oscillations proportional in strength to the speed of the ramp, and that ramps induce oscillations with worse outcomes than oscillations induced by step change of CO2 concentration. An initial overshoot is caused by a temporary excess of available phosphate. During the overshoot, the plant outperforms steady-state TPU and ribulose 1,5-bisphosphate regeneration limitations of photosynthesis, but cannot exceed the rubisco limitation. We performed additional optical measurements which support the role of PSI reduction and oscillations in availability of NADP+ and ATP in supporting oscillations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

On the dual advantage of placing observations through forward sensitivity analysis

The four-dimensional variational data assimilation methodology for assimilating noisy observations into a deterministic model has been the workhorse in forecasting centres for over three decades. While this method provides a computationally efficient framework for dynamic data assimilation, it is largely silent on the important question concerning the minimum number and placement of observations. To answer this question, we demonstrate the dual advantage of placing the observations where the square of the sensitivity of the model solution with respect to the unknown control variables, called forward sensitivities, attains its maximum. By following this approach, we can force the observability Gramian to be of full rank, which in turn guarantees efficient recovery of the optimal values of the control variables, which is the first of the two advantages of this strategy. We further show that the proposed strategy of placing observations has another inherent optimality: the square of the sensitivity of the optimal estimates of the control with respect to the observations (used to obtain these estimates) attains its minimum value, a second advantage that is a direct consequence of the above strategy for placing observations. Our analytical framework and numerical experiments on linear and nonlinear systems confirm the effectiveness of our proposed strategy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Combining Observations and Models: A Review of the CARDAMOM Framework for Data‐Constrained Terrestrial Ecosystem Modeling

The rapid increase in the volume and variety of terrestrial biosphere observations (i.e., remote sensing data and in situ measurements) offers a unique opportunity to derive ecological insights, refine process‐based models, and improve forecasting for decision support. However, despite their potential, ecological observations have primarily been used to benchmark process‐based models, as many past and current models lack the capability to directly integrate observations and their associated uncertainties for parameterization. In contrast, data assimilation frameworks such as the CARbon DAta MOdel fraMework (CARDAMOM) and its suite of process‐based models, known as the Data Assimilation Linked Ecosystem Carbon Model (DALEC), are specifically designed for model‐data fusion. This review, motivated by a recent CARDAMOM community workshop, examines the development and applications of CARDAMOM, with an emphasis on its role in advancing ecosystem process understanding. CARDAMOM employs a Bayesian approach, using a Markov Chain Monte Carlo algorithm to enable data‐driven calibration of DALEC parameters and initial states (i.e., carbon pool sizes) through observation operators. CARDAMOM's unique ability to retrieve localized model process parameters from diverse datasets—ranging from in situ measurements to global satellite observations—makes it a highly flexible tool for analyzing spatially variable ecosystem responses to environmental change. However, assimilating these data also presents challenges, including data quality issues that propagate into model skill, as well as trade‐offs between model complexity, parameter equifinality, and predictive performance. We discuss potential solutions to these challenges, such as reducing parameter equifinality by incorporating new observations. This review also offers community recommendations for incorporating emerging datasets, integrating machine learning techniques, strengthening collaboration with remote sensing, field, and modeling communities, and expanding CARDAMOM's relevance for localized ecosystem monitoring and decision‐making. CARDAMOM enables a deep, mechanistic understanding of terrestrial ecosystem dynamics that cannot be achieved through empirical analyses of observational datasets or weakly constrained models alone.

Bayesian inference↗
DOE OSTI2022

Autodifferentiable Ensemble Kalman Filters

Data assimilation is concerned with sequentially estimating a temporally evolving state. This task, which arises in a wide range of scientific and engineering applications, is particularly challenging when the state is high-dimensional and the state-space dynamics are unknown. This paper introduces a machine learning framework for learning dynamical systems in data assimilation. Here, our auto-differentiable ensemble Kalman filters (AD-EnKFs) blend ensemble Kalman filters for state recovery with machine learning tools for learning the dynamics. In doing so, AD-EnKFs leverage the ability of ensemble Kalman filters to scale to high-dimensional states and the power of automatic differentiation to train high-dimensional surrogate models for the dynamics. Numerical results using the Lorenz-96 model show that AD-EnKFs outperform existing methods that use expectation-maximization or particle filters to merge data assimilation and machine learning. In addition, AD-EnKFs are easy to implement and require minimal tuning.

autodifferentiation↗
DOE OSTI2016

Final Technical Report for Department of Energy Award DE-SC0006625, “Predictability of the carbon-climate system on seasonal to decadal time scales.”

The project aims to investigate the feasibility of advancing our understanding of the carbon cycle, using a carbon-weather data assimilation system that updates the modeled carbon dioxide concentration and atmospheric circulation every six hours using CO 2 data (from the OCO 2 satellite) and weather data. At the core of the system is the DOE-NCAR-CAM5fv global circulation model coupled to the National Center for Atmospheric Research's Data Assimilation Testbed, running an ensemble of 30 models. This combination provides realistic vertical carbon dioxide gradients and conservation of dry air mass. A global four-dimensional distribution of atmospheric CO 2 concentration is produced. Our results show (1) that OCO 2 total precipitable water data are reliable and provide valuable uncertainty information for the OCO 2 data assimilation; and (2) that our approach is a promising method for monitoring national carbon dioxide emissions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2005

Diagnostic Comparison of Meteorological Analyses during the 2002 Antarctic Winter

Several meteorological datasets, including U.K. Met Office (MetO), European Centre for Medium-Range Weather Forecasts (ECMWF), National Centers for Environmental Prediction (NCEP), and NASA's Goddard Earth Observation System (GEOS-4) analyses, are being used in studies of the 2002 Southern Hemisphere (SH) stratospheric winter and Antarctic major warming. Diagnostics are compared to assess how these studies may be affected by the meteorological data used. While the overall structure and evolution of temperatures, winds, and wave diagnostics in the different analyses provide a consistent picture of the large-scale dynamics of the SH 2002 winter, several significant differences may affect detailed studies. The NCEP-NCAR reanalysis (REAN) and NCEP-Department of Energy (DOE) reanalysis-2 (REAN-2) datasets are not recommended for detailed studies, especially those related to polar processing, because of lower-stratospheric temperature biases that result in underestimates of polar processing potential, and because their winds and wave diagnostics show increasing differences from other analyses between similar to 30 and 10 hPa (their top level). Southern Hemisphere polar stratospheric temperatures in the ECMWF 40-Yr Re-analysis (ERA-40) show unrealistic vertical structure, so this long-term reanalysis is also unsuited for quantitative studies. The NCEP/Climate Prediction Center (CPC) objective analyses give an inferior representation of the upper-stratospheric vortex. Polar vortex transport barriers are similar in all analyses, but there is large variation in the amount, patterns, and timing of mixing, even among the operational assimilated datasets (ECMWF, MetO, and GEOS-4). The higher-resolution GEOS-4 and ECMWF assimilations provide significantly better representation of filamentation and small-scale structure than the other analyses, even when fields gridded at reduced resolution are studied. The choice of which analysis to use is most critical for detailed transport studies (including polar process modeling) and studies involving synoptic evolution in the upper stratosphere. The operational assimilated datasets are better suited for most applications than the NCEP/CPC objective analyses and the reanalysis datasets.

stratosphere↗
NASA NTRS2011

An Intercomparison of Tropospheric Ozone Retrievals Derived from Two Aura Instruments and Measurements in Western North America in 2006

Two recently developed methods for quantifying tropospheric ozone abundances based on Aura data, the Trajectory‐enhanced Tropospheric Ozone Residual (TTOR) and an assimilation of Aura data into Goddard Earth Observing System Version 4 (ASM), are compared to ozone measurements from ozonesonde data collected in April-May 2006 during the INTEX Ozonesonde Network Study 2006 (IONS‐06) campaign. Both techniques use Ozone Monitoring Instrument (OMI) and Microwave Limb Sounder (MLS) observations. Statistics on column ozone amounts for both products are presented. In general, the assimilation compares better to sonde integrated ozone to 200 hPa (28.6% difference for TTOR versus 2.7% difference for ASM), and both products are biased low. To better characterize the performance of ASM, ozone profiles based on the assimilation are compared to those from ozonesondes. We noted slight negative biases in the lower troposphere, and slight positive biases in the upper troposphere/lower stratosphere (UT/ LS), where we observed the greatest variability. Case studies were used to further understand ASM performance. We examine one case from 17 April 2006 at Bratt's Lake, Saskatchewan, where geopotential height gradients appear to be related to an underestimation in the ASM in the UT/LS region. A second case, from 21 April 2006 at Trinidad Head, California, is a situation where the overprediction of ozone in the UT/LS region does not appear to be due to current dynamic conditions but seems to be related to uncertainty in the flow pattern and large differences in MLS observations upstream.

UTLS↗
DOE OSTI2015

Development of fine-resolution analyses and expanded large-scale forcing properties. Part I: Methodology and evaluation

We produce fine-resolution, three-dimensional fields of meteorological and other variables for the U.S. Department of Energy’s Atmospheric Radiation Measurement (ARM) Southern Great Plains site. The Community Gridpoint Statistical Interpolation system is implemented in a multiscale data assimilation (MS-DA) framework that is used within the Weather Research and Forecasting model at a cloud-resolving resolution of 2 km. The MS-DA algorithm uses existing reanalysis products and constrains fine-scale atmospheric properties by assimilating high-resolution observations. A set of experiments show that the data assimilation analysis realistically reproduces the intensity, structure, and time evolution of clouds and precipitation associated with a mesoscale convective system. Evaluations also show that the large-scale forcing derived from the fine-resolution analysis has an overall accuracy comparable to the existing ARM operational product. For enhanced applications, the fine-resolution fields are used to characterize the contribution of subgrid variability to the large-scale forcing and to derive hydrometeor forcing, which are presented in companion papers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Formaldehyde as a carbon and electron shuttle between autotroph and heterotroph populations in acidic hydrothermal vents of Norris Geyser Basin, Yellowstone National Park

The Norris Geyser Basin in Yellowstone National Park contains a large number of hydrothermal systems, which host microbial populations supported by primary productivity associated with a suite of chemolithotrophic metabolisms. We demonstrate that Metallosphaera yellowstonesis MK1, a facultative autotrophic archaeon isolated from a hyperthermal acidic hydrous ferric oxide (HFO) spring in Norris Geyser Basin, excretes formaldehyde during autotrophic growth. To determine the fate of formaldehyde in this low organic carbon environment, we incubated native microbial mat (containing M. yellowstonensis) from a HFO spring with 13C-formaldehyde. Isotopic analysis of incubation-derived CO2 and biomass showed that formaldehyde was both oxidized and assimilated by members of the community. Autotrophy, formaldehyde oxidation, and formaldehyde assimilation displayed different sensitivities to chemical inhibitors, suggesting that distinct sub-populations in the mat selectively perform these functions. Our results demonstrate that electrons originally resulting from iron oxidation can energetically fuel autotrophic carbon fixation and associated formaldehyde excretion, and that formaldehyde is both oxidized and assimilated by different organisms within the native microbial community. Thus, formaldehyde can effectively act as a carbon and electron shuttle connecting the autotrophic, iron oxidizing members with associated heterotrophic members in the HFO community.

acidophilic↗
DOE OSTI2026

Improving coastal water level estimation by merging nadir-only satellite altimetry data into a hydrodynamic model

Providing robust real time flood warnings is of paramount importance to coastal communities. Although state-of-the-art hydrodynamic models are capable of robustly predicting Coastal Water Levels (CWL), unresolved drivers affecting level fluctuations are often not represented by the model governing equations. This work evaluates a novel method to improve the performance of the ADvanced CIRCulation (ADCIRC) hydrodynamic model by assimilating observations from four nadir-only satellite altimetry missions against a set of National Oceanic and Atmospheric Administration (NOAA) gauge stations located across the entire U.S. East Coast. Two different types of simulations were performed – Open Loop (OL) and Data Assimilation (DA). Five different simulations were performed where four different satellite altimetry observations were assimilated individually and combined with two different scenarios – with and without considering the data quality flags. Results indicate that, despite their limited spatial coverage, merging nadir-only observations into ADCIRC from the newly launched Surface Water and Ocean Topography (SWOT)’s nadir altimeter can improve the model performance at 76% of the gauge locations, whereas Sentinel-6 improves 73% of the total locations, Jason-3 74%, and SARAL 21%. Furthermore, combining observations from SWOT-nadir, Jason-3, and Sentinel-6 can improve the ADCIRC performance at more than 80% of the gauge locations for 107-day simulation. Nadir-only satellite altimetry observations can be useful for improving the model performance even if flagged as “poor quality” near the coast. When the flagged data are disregarded, SWOT can improve ADCIRC at 78%, Sentinel-6 at 73%, Jason-3 at 53%, and SARAL at 21% of the gauge locations. The ability to improve the model simulations largely depends on the availability of a satellite overpass nearby. Therefore, model performance can be further enhanced if satellite observations are available during a storm surge event, stressing the importance of frequent satellite overpasses.

Aafnan Bhuiyan, Soelem↗