Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Adjoint Models”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2014

Investigating the Sensitivity of Surface-Level Nitrate Seasonality in Antarctica to Primary Sources Using a Global Model

Determining the sources of total nitrate (TNIT NO-3 + HNO3) reaching Antarctica is a long-standing challenge. Here we analyze the monthly sensitivity of surface-level TNIT in Antarctica to primary sources using a global 3-D chemical transport model, GEOS-Chem, and its adjoint. Modeled seasonal variation of TNIT concentrations shows good agreement with several measurement studies, given that the lack of post-depositional processing in the model leads to an expected underestimate of maximum values in November through January. Remote NOx sources have the greatest impact May-July, during when the model background concentrations are sensitive to NOx emissions from fossil fuel combustion, soil, and lightning originating from 25S to 65S. In this season, NOx is transported to Antarctica as TNIT, which is formed above continental source regions at an altitude of 5-11 km. In other seasons, more NOx is transported as a reservoir species (e.g., peroxyacetyl nitrate, PAN) through the free troposphere, transforming into TNIT within a cone of influence that extends to 35S and above 4 km altitude. Photolysis of PAN over Antarctica is the main driver of modeled NOx seasonality. Stratospheric production and loss of tracers are relatively unimportant in monthly sensitivities in GEOS-Chem, driving only a few percent of surface level variability of TNIT. A small peak concentration in August is captured by the model, although some measured values in August fall outside the range of simulated concentrations. Modifications to the model to represent sedimentation of polar stratospheric clouds (PSCs) lead to increased surface level August TNIT concentrations. However, this simple representation does not explicitly account for PSC particle deposition or disappearance of the tropopause in the middle of winter, and thus the influence of stratospheric nitrate sources estimated in this study is likely a lower bound.

stratosphere↗
NASA NTRS1987

Variational four-dimensional analysis using quasi-geostrophic constraints

A variational four-dimensional analysis technique using quasi-geostrophic models as constraints is examined using gridded fields as data. The analysis method uses a standard iterative nonlinear minimization technique to find the solution to the constraining forecast model which best fits the data as measured by a predefined functional. The minimization algorithm uses the derivative of the functional with respect to each of the initial condition values. This derivative vector is found by inserting the weighted differences between the model solution and the inserted data into a backwards integrating adjoint model. The four-dimensional analysis system was examined by applying it to fields created from a primitive equations model forecast and to fields created from satellite retrievals. The results show that the technique has several interesting characteristics not found in more traditional four-dimensional assimilation techniques. These features include a close fit of the model solution to the observations throughout the analysis interval and an insensitivity to the frequency of data insertion or the amount of data. The four-dimensional analysis technique is very versatile and can be extended to more complex problems with little theoretical difficulty.

Derber, John C.↗
NASA NTRS2009

The GMAO ADVAR System: Preliminary Results

The Goddard Earth Observing System (GEOS-S) Data Assimilation System is now configurable to run strong constraint 4DVAR. The system combines the 4DVAR-capable Grid-point Statistical Interpolation analysis with the GEOS-S general circulation model and an early version of the tangent linear and adjoint models of the finite-volume hydrodynamics with simplified physics. The GMAO 4DVAR uses a Lanczos-based conjugate gradient algorithm in a nested resolution inner-loop setting. A range of tests and experiments will have been carried out by the time of the Symposium. Results on balance-related issues, resolution configuration, and choice of assimilation time window will be shown. A comparison with the Global Modeling and Assimilation Office (GMAO) 3DVAR and a discussion of future directions of development will also be presented.

Todling, Ricardo↗
NASA NTRS2020

Robustness and Behavior of Adjoint Calculations of Observation Impacts in Numerical Weather Prediction

Adjoint models are powerful tools that can be used to estimate the impact of observations on a chosen norm for numerical weather prediction forecasts. In this study, the Global Modeling and Assimilation Office (NASA/GMAO) Observing System Simulation Experiment framework is employed to investigate the behavior of the adjoint tool in an environment where the 'true' state of the atmosphere is fully known. This allows for the calculation of adjoint estimates of observation impact for very short forecast times including the zero-hour analysis state. The adjoint calculations using self-analysis verification can also be compared to adjoint calculations using the 'truth' as verification in order to characterize the robustness of adjoint estimations in the operational setting. Results from a experiments exploring various aspects of performance of the adjoint tool will be presented.

Prive, Nikki C.↗
DOE OSTI2025

Technical Report on Subsurface Monitoring of the Brady Hot Spring Geothermal Site, Nevada, based upon Full Waveform Inversion

Abilities to accurately characterize the subsurface in a geothermal setting is key to assess and support production. An important element of geothermal reservoir monitoring is also the ability to investigate fluid transport within fracture network. This report focuses on improving subsurface imaging and monitoring in geothermal settings using full waveform inversion based on the adjoint method and time-lapse imaging. To assess our method, we rely on a dense seismic dataset collected in 2016 at the Brady Hot Springs geothermal site in Nevada for the DOE-funded project Poroelastic Tomography by Adjoint Inverse Modeling of Data from Seismology, Geodesy, and Hydrology. This dataset captures subsurface changes across four stages of geothermal power plant operations, which involve varying rates of fluid injection and extraction. Two velocity models were previously derived from this dataset using different methods: one based on travel times and another on sweep interferometry. Our first step is to refine these models using adjoint tomography, which has been applied successfully at global and regional-scales but is less common at the reservoir-scale. Two approaches are then explored for time-lapse analysis: directly comparing refined tomographic models from different stages or backpropagating waveform differences relative to a baseline tomographic model. The main take away is that both approaches highlight similar reservoir behaviors, but the latter approach is more computationally effective in capturing small-scale changes in subsurface properties. For this work, we leverage the use of Salvus (www.mondaic.com), an end-to-end seismic imaging solution, relying on the spectral element method to compute forward and adjoint simulations, and developed by Mondaic Ltd. It includes integrated workflow management that handles waveform and metadata, launches simulations, computes waveform misfits and adjoint sources, and iterates for model updates by nonlinear optimization.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Development and evaluation of a new 4DEnVar-based weakly coupled ocean data assimilation system in E3SMv2

The development, implementation, and evaluation of a new weakly coupled ocean data assimilation (WCODA) system for the fully coupled Energy Exascale Earth System Model version 2 (E3SMv2) utilizing the four-dimensional ensemble variational (4DEnVar) method are presented in this study. The 4DEnVar method, based on the dimension-reduced projection four-dimensional variational (DRP-4DVar) approach, replaces the adjoint model with the ensemble technique, thereby reducing computational demands. Monthly mean ocean temperature and salinity data from the EN4.2.1 reanalysis are integrated into the ocean component of E3SMv2 from 1950 to 2021 with the goal of providing realistic initial conditions for decadal predictions and predictability studies. The performance of the WCODA system is assessed using various metrics, including the reduction rate of the cost function, root mean square error (RMSE) differences, correlation differences, and model biases. Results indicate that the WCODA system effectively assimilates the reanalysis data into the climate model, consistently achieving negative reduction rates of the cost function and notable improvements in RMSE and correlation across various ocean layers and regions. Significant enhancements are observed in the upper ocean layers across the majority of global ocean regions, particularly in the north Atlantic, north Pacific, and Indian Ocean. Model biases in sea surface temperature and salinity are also substantially reduced. For sea surface temperature, cold biases in the north Pacific and north Atlantic are diminished by about 1–2 °C, and warm biases in the Southern Ocean are corrected by approximately 1.5–2.5 °C. In terms of salinity, improvements are observed with bias reductions of about 0.5–1 psu in the north Atlantic and north Pacific and up to 1.5 psu in parts of the Southern Ocean. The ultimate goal of the WCODA system is to advance the predictive capabilities of E3SM for subseasonal to decadal climate predictions, thereby supporting research on strategic energy-sector policies and planning.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2002

Mesoscale Assimilation of TRMM Data with 4DVAR: Preliminary Results

Surface rainfall data, derived from the TRMM Microwave Image (TMI), are assimilated into the PSU/NCAR MM5 model using a 4DVAR technique. Preliminary experiments are performed to incorporate TRMM rainfall data into a hurricane initialization. It is found that the rainfall data assimilation is sensitive to the error characteristics of the data and the physics in the adjoint model. In addition, assimilating the rainfall data alone produces a more realistic eye and rain bands in the hurricane but cannot ensure improvements of hurricane intensity forecasts. Numerical results indicate that it is necessary to incorporate TRMM rainfall data together with other types of data such as wind data into the model, in which case the inclusion of the rainfall data will further improve the intensity forecast of the hurricane. This fact might imply that some proper constraints will be needed for the rainfall assimilation. Relevant results and issues will be presented.

Pu, Zhoa-Xia↗
NASA NTRS2012

Spatially Refined Aerosol Direct Radiative Forcing Efficiencies

Global aerosol direct radiative forcing (DRF) is an important metric for assessing potential climate impacts of future emissions changes. However, the radiative consequences of emissions perturbations are not readily quantified nor well understood at the level of detail necessary to assess realistic policy options. To address this challenge, here we show how adjoint model sensitivities can be used to provide highly spatially resolved estimates of the DRF from emissions of black carbon (BC), primary organic carbon (OC), sulfur dioxide (SO2), and ammonia (NH3), using the example of emissions from each sector and country following multiple Representative Concentration Pathway (RCPs). The radiative forcing efficiencies of many individual emissions are found to differ considerably from regional or sectoral averages for NH3, SO2 from the power sector, and BC from domestic, industrial, transportation and biomass burning sources. Consequently, the amount of emissions controls required to attain a specific DRF varies at intracontinental scales by up to a factor of 4. These results thus demonstrate both a need and means for incorporating spatially refined aerosol DRF into analysis of future emissions scenario and design of air quality and climate change mitigation policies.

Climate↗
NASA NTRS2014

Carbon Monitoring System Flux Estimation and Attribution: Impact of ACOS-GOSAT X(CO2) Sampling on the Inference of Terrestrial Biospheric Sources and Sinks

Using an Observing System Simulation Experiment (OSSE), we investigate the impact of JAXA Greenhouse gases Observing SATellite 'IBUKI' (GOSAT) sampling on the estimation of terrestrial biospheric flux with the NASA Carbon Monitoring System Flux (CMS-Flux) estimation and attribution strategy. The simulated observations in the OSSE use the actual column carbon dioxide (X(CO2)) b2.9 retrieval sensitivity and quality control for the year 2010 processed through the Atmospheric CO2 Observations from Space algorithm. CMS-Flux is a variational inversion system that uses the GEOS-Chem forward and adjoint model forced by a suite of observationally constrained fluxes from ocean, land and anthropogenic models. We investigate the impact of GOSAT sampling on flux estimation in two aspects: 1) random error uncertainty reduction and 2) the global and regional bias in posterior flux resulted from the spatiotemporally biased GOSAT sampling. Based on Monte Carlo calculations, we find that global average flux uncertainty reduction ranges from 25% in September to 60% in July. When aggregated to the 11 land regions designated by the phase 3 of the Atmospheric Tracer Transport Model Intercomparison Project, the annual mean uncertainty reduction ranges from 10% over North American boreal to 38% over South American temperate, which is driven by observational coverage and the magnitude of prior flux uncertainty. The uncertainty reduction over the South American tropical region is 30%, even with sparse observation coverage. We show that this reduction results from the large prior flux uncertainty and the impact of non-local observations. Given the assumed prior error statistics, the degree of freedom for signal is approx.1132 for 1-yr of the 74 055 GOSAT X(CO2) observations, which indicates that GOSAT provides approx.1132 independent pieces of information about surface fluxes. We quantify the impact of GOSAT's spatiotemporally sampling on the posterior flux, and find that a 0.7 gigatons of carbon bias in the global annual posterior flux resulted from the seasonally and diurnally biased sampling when using a diagonal prior flux error covariance.

biased sampling↗
NASA NTRS1995

Sensitivity analysis of a model of CO2 exchange in tundra ecosystems by the adjoint method

A model of net primary production (NPP), decomposition, and nitrogen cycling in tundra ecosystems has been developed. The adjoint technique is used to study the sensitivity of the computed annual net CO2 flux to perturbation in initial conditions, climatic inputs, and model's main parameters describing current seasonal CO2 exchange in wet sedge tundra at Barrow, Alaska. The results show that net CO2 flux is most sensitive to parameters characterizing litter chemical composition and more sensitive to decomposition parameters than to NPP parameters. This underlines the fact that in nutrient-limited ecosystems, decomposition drives net CO2 exchange by controlling mineralization of main nutrients. The results also indicate that the short-term (1 year) response of wet sedge tundra to CO2-induced warming is a significant increase in CO2 emission, creating a positive feedback to atmosphreic CO2 accumulation. However, a cloudiness increase during the same year can severely alter this response and lead to either a slight decrease or a strong increase in emitted CO2, depending on its exact timing. These results demonstrate that the adjoint method is well suited to study systems encountering regime changes, as a single run of the adjoint model provides sensitivities of the net CO2 flux to perturbations in all parameters and variables at any time of the year. Moreover, it is shown that large errors due to the presence of thresholds can be avoided by first delimiting the range of applicability of the adjoint results.

Waelbroek, C.↗
DOE OSTI2024

Adjoint waveform tomography of East Asia for improved waveform prediction

We present a preliminary version of the East Asia Tomography (EAT) model, an adjoint waveform tomography model of East and Southeast Asia. We used SPiRaL (Simmons et al., 2021) as our starting model and source parameters for 238 earthquakes from the Global Centroid Moment Tensor catalogue (Ekström et al., 2012). After 50 iterations on Lawrence Livermore National Laboratory’s Lassen supercomputer, we converge on a model with a minimum period of 50 seconds. The preliminary model shows improved slab structure compared to SPiRaL and shows significantly reduced misfit. In later versions of the model, we aim to harness techniques proposed in other studies to improve waveform predictions that travel through the ocean (e.g., Wehner et al., 2022) and use relative amplitude-based misfit functions (e.g., Tao et al., 2018) to better constrain Earth structure. We hope to iterate the current extent of the model to 25 seconds minimum period before iterating to shorter periods for a smaller subregion of the full model.

58 GEOSCIENCES↗
DOE OSTI2026

Technical Report on Adjoint Waveform Tomography of East Asia for Improved Waveform Prediction

We present a preliminary version of the East Asia Tomography (EAT) model, an adjoint waveform tomography model of East and Southeast Asia. We used SPiRaL (Simmons et al., 2021) as our starting model and source parameters for 250 earthquakes from the Global Centroid Moment Tensor catalogue (Ekström et al., 2012). Over 198 iterations, we have iterated the EAT model down to a minimum period of 35 seconds. We plan on continuing our iteration technique down to 30 seconds period before updating our misfit function to use a normalized cross correlation-based misfit functions (e.g., Tao et al., 2018) to better constrain Earth structure. We hope to iterate the current extent of the model to 25 seconds minimum period before iterating to shorter periods for a smaller subregion of the full model.

58 GEOSCIENCES↗
NASA NTRSDate not supplied

Investigation of a Smooth Local Correlation-based Transition Model in a Discrete-Adjoint Aerodynamic Shape Optimization Algorithm

A smooth local correlation-based transition model is fully coupled to a RANS-based Newton-Krylov flow solver and discrete-adjoint gradient-based optimization algorithm. The free-transition optimization framework is evaluated using lift-constrained drag minimizations of airfoils at design conditions ranging from light to single-aisle aircraft and an infinite swept wing at design conditions representative of a transonic strut-braced wing aircraft. The impact of the streamwise grid resolution on the ability of the optimization algorithm to delay boundary-layer transition is investigated, with the results demonstrating that streamwise grid resolution requirements increase as the transition length decreases with increasing Reynolds number. The optimization problem at the light aircraft design conditions is demonstrated to be multi-modal, with the optimization algorithm producing two distinct designs: one with a thin, reflexed trailing edge and steep pressure recovery regions, the other with increased aft loading, with the latter design outperforming the former. A drag minimization of an airfoil at transonic design conditions demonstrates that the optimization algorithm successfully trades a decrease in viscous drag by delaying boundary-layer transition with an increase in wave drag, while the drag minimization of an infinite swept wing demonstrates the capability of the optimizational gorithm to delay both Tollmien-Schlichting and stationary crossflow instabilities.

AATT↗
NASA NTRS2003

Adjoint-Based Methods for Estimating CO2 Sources and Sinks from Atmospheric Concentration Data

Work to develop adjoint-based methods for estimating CO2 sources and sinks from atmospheric concentration data was initiated in preparation for last year's summer institute on Carbon Data Assimilation (CDAS) at the National Center for Atmospheric Research in Boulder, CO. The workshop exercises used the GSFC Parameterized Chemistry and Transport Model and its adjoint. Since the workshop, a number of simulations have been run to evaluate the performance of the model adjoint. Results from these simulations will be presented, along with an outline of challenges associated with incorporating a variety of disparate data sources, from sparse, but highly precise, surface in situ observations to less accurate, global future satellite observations.

Andrews, Arlyn E.↗
DOE OSTI2025

A direct-adjoint approach for material point model calibration with application to plasticity

Here, this paper proposes a new approach for the calibration of material parameters in local elastoplastic constitutive models. The calibration is posed as a constrained optimization problem, where the constitutive model evolution equations for a single material point serve as constraints. The objective function quantifies the mismatch between the stress predicted by the model and corresponding experimental measurements. To improve calibration efficiency, a novel direct-adjoint approach is presented to compute the Hessian of the objective function, which enables the use of second-order optimization algorithms. Automatic differentiation is used for gradient and Hessian computations. Two numerical examples are employed to validate the Hessian matrices and to demonstrate that the Newton–Raphson algorithm consistently outperforms gradient-based algorithms such as L-BFGS-B.

36 MATERIALS SCIENCE↗