Search NASA⌕ Search

SEARCH · Search NASA

Results for “collocation”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

On optimal control of hybrid dynamical systems using complementarity constraints

Optimal control for switch-based dynamical systems is a challenging problem in the process control literature. In this study, we model these systems as hybrid dynamical systems with finite number of unknown switching points and reformulate them using non-smooth and non-convex complementarity constraints as a mathematical program with complementarity constraints (MPCC). We utilize a moving finite element based strategy to discretize the differential equation system to accurately locate the unknown switching points at the finite element boundary and achieve high-order accuracy at intermediate non-collocation points. We propose a globalization approach to solve the discretized MPCC problem using a mixed NLP/MILP-based strategy to converge to a non-spurious first-order optimal solution. The method is tested on three dynamic optimization examples, including a gas–liquid tank model and an optimal control problem with a sliding mode solution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Performing a multi-unit level-3 PSA with MACCS

MACCS (MELCOR Accident Consequence Code System), WinMACCS, and MelMACCS now facilitate a multi-unit consequence analysis. MACCS evaluates the consequences of an atmospheric release of radioactive gases and aerosols into the atmosphere and is most commonly used to perform probabilistic safety assessments (PSAs) and related consequence analyses for nuclear power plants (NPPs). WinMACCS is a user-friendly preprocessor for MACCS. MelMACCS extracts source-term information from a MELCOR plot file. The current development can combine an arbitrary number of source terms, representing simultaneous releases from a multi-unit facility, into a single consequence analysis. The development supports different release signatures, fission product inventories, and accident initiation times for each unit. The treatment is completely general except that the model is currently limited to collocated units. A major practical consideration for performing a multi-unit PSA is that a comprehensive treatment for more than two units may involve an intractable number of combinations of source terms. This paper proposes and evaluates an approach for reducing the number of calculations to be tractable, even for sites with eight or ten units. The approximation error introduced by the approach is acceptable and is considerably less than other errors and uncertainties inherent in a Level 3 PSA.

42 ENGINEERING↗
DOE OSTI2024

Use of physics to improve solar forecast: Part III, impacts of different cloud types

Cloud-type impacts present a great challenge to solar forecasting due to diverse and complex cloud-radiation interactions. This third part of our paper sequence seeks to address this challenge by quantifying the forecast accuracies under eight cloud types: cumulus (Cu), stratified clouds (St), altocumulus (Ac), altostratus (As), cirrostratus/anvil (Cr), cirrus (Ci), congestus (Co), deep convective clouds (Dc) across four physics-informed persistence models reported in Part I. To generalize the cloud impacts, the eight cloud types are further grouped into three cloud categories based on their common features: weak convective clouds, stratiform clouds, and strong convective clouds. Here, the decade-long (2001 ~ 2014) collocated measurements of irradiances and cloud types at the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Program South Great Plain (SGP) Central Facility site are used for model evaluation. Results reveal a clear performance hierarchy for global horizontal irradiance (GHI) and direct normal irradiance (DNI): best for weak convective clouds and cirrus, intermediate for stratiform clouds, and worst for strong convective clouds. Performance for diffuse horizontal irradiance (DHI) is less influenced by cloud types. Cloud albedo dominates all three irradiances for Dc, while both cloud albedo and cloud fraction are influential for other cloud types. A 12 %~33 % improvement in accuracy at 6-hour lead time compared to the benchmark smart model confirms the effectiveness of incorporating physics into the models for various cloud types; further improvements are expected by directly integrating cloud type information into forecasting models by modifying the physical formulation of cloud-radiation interaction, and/or using more advanced machine learning models.

14 SOLAR ENERGY↗
DOE OSTI2024

Informative and non-informative decomposition of turbulent flow fields

Not all the information in a turbulent field is relevant for understanding particular regions or variables in the flow. Here, we present a method for decomposing a source field into its informative Φ I (x, t) and residual Φ R (x, t) components relative to another target field. The method is referred to as informative and non-informative decomposition (IND). All the necessary information for physical understanding, reduced-order modelling and control of the target variable is contained in Φ I (x, t), whereas Φ R (x, t) offers no substantial utility in these contexts. The decomposition is formulated as an optimisation problem that seeks to maximise the time-lagged mutual information of the informative component with the target variable while minimising the mutual information with the residual component. The method is applied to extract the informative and residual components of the velocity field in a turbulent channel flow, using the wall shear stress as the target variable. We demonstrate the utility of IND in three scenarios: (i) physical insight into the effect of the velocity fluctuations on the wall shear stress; (ii) prediction of the wall shear stress using velocities far from the wall; and (iii) development of control strategies for drag reduction in a turbulent channel flow using opposition control. In case (i), IND reveals that the informative velocity related to wall shear stress consists of wall-attached high- and low-velocity streaks, collocated with regions of vertical motions and weak spanwise velocity. This informative structure is embedded within a larger-scale streak–roll structure of residual velocity, which bears no information about the wall shear stress. In case (ii), the best-performing model for predicting wall shear stress is a convolutional neural network that uses the informative component of the velocity as input, while the residual velocity component provides no predictive capabilities. Finally, in case (iii), we demonstrate that the informative component of the wall-normal velocity is closely linked to the observability of the target variable and holds the essential information needed to develop successful control strategies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

3-D Geological Modeling for Numerical Flow Simulation Studies of Gas Hydrate Reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

58 GEOSCIENCES↗
DOE OSTI2021

Prediction of the Thermal Runaway Limit and Optimal Operation of Heat Transfer-Limited, Fixed-Bed Reactor Systems

Here, we derive a new prediction for thermal runaway starting from the alpha model for fixed-bed reactor systems. This method accounts for thermal resistance internal to the reactor tube and the radial temperature gradients that result. To showcase our method, we compare its predictions to other common criteria for thermal runaway using o-xylene oxidation as the example chemistry. Even in systems where internal heat transfer is negligible, the empirical practical design criterion for thermal runaway is inaccurate. For cases where internal heat transfer is relevant, our runaway limit is more stringent than limits derived from simpler 1-D models. To augment our work, we optimize the product yield with the thermal runaway constraint using orthogonal collocation. Using the alpha model, the results illustrate that the thermal runaway limit can be accurately determined using either numerical or analytical methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Ab Initio Computation of Rotationally-Averaged Pump–Probe X-ray and Electron Diffraction Signals

We develop a new algorithm for the computation of the rotationally averaged elastic molecular diffraction signal for the cases of perpendicular or parallel pump–probe geometries. The algorithm first collocates the charge density from an arbitrary ab initio wave function onto a Becke quadrature grid [A. Becke, J. Chem. Phys. 1988, 88, 2457], providing a high-fidelity multiresolution representation of the charge density. A double sum is then performed over the Becke grid points, and the interaction between points computed using the scattering kernels of Williamson and Zewail [J. C. Williamson and A. H. Zewail, J. Phys. Chem. 1994, 98, 2766]. These kernels analytically average over the molecular orientations with the cos 2 γ selection factor appropriate for one-photon dipole absorption in a perpendicular pump–probe geometry. We demonstrate that the method is converged with small grids containing <500 points/atom. We implement the algorithm on a GPU for increased efficiency and emonstrate the algorithm for molecules with up to a few dozen atoms. We explore the accuracy of the independent atom model (IAM) by comparison with our new and more accurate method. We also investigate the possibility of detecting signatures of electronic transitions in polyatomic pump–probe diffraction experiments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Contribution of Speciated Monoterpenes to Secondary Aerosol in the Eastern North Atlantic

Extension of laboratory results suggests that monoterpenes, a subset of marine biogenic volatile organic compounds, might play a competitive role in secondary aerosol formation in the remote marine atmosphere, where the oxidation of dimethyl sulfide has traditionally been thought to dominate. However, the current assessment of the role of monoterpenes in secondary aerosol formation in the remote marine atmosphere is limited by the scarcity of speciated monoterpene measurements and the complete absence of collocated measurements of particle chemical composition and monoterpene speciation. Here, in this paper, we present measurements of gas-phase volatile organic compounds, with particular focus on speciated monoterpenes, and commensurate measurements of aerosol chemical composition in the Eastern North Atlantic. The average monoterpene concentration in periods of marine air was 14 ± 10 ppt, and the dominant isomer was β-pinene (46 ± 10%), with other contributions from α-pinene, limonene, and β-ocimene. The total monoterpene concentration in marine air was greater than that of isoprene by a factor of 3.9 on average, in contrast to prior research and potentially due to the large β-pinene contribution. Monoterpenes were significantly smaller in concentration than dimethyl sulfide, which was also highlighted by the substantial contribution of sulfate to non-refractory submicron aerosol. The mean and interquartile range of the sulfate to organic aerosol mass ratio in clean marine air were 2.7 (0.7–3.4). Finally, a simple box model analysis showed that the BVOC precursors, dimethyl sulfide, methanethiol, monoterpenes, and isoprene, could sustain the measured sulfate to organic aerosol ratio in clean marine air.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Projecting Changes in Societally Impactful Northeastern U.S. Snowstorms

The northeastern United States is vulnerable to many impacts from snowfall-producing winter cyclones that are amplified by the proximity of population centers to storm tracks. Historically, climatic snowfall assessments have centered around seasonal means even though local impacts typically occur at scales of hours to days. To detect snowstorms at the event level, an objective algorithm is defined based on the Regional Snowfall Index. The metric collocates storm snowfall with population to produce statistics of snowstorms with societal impacts. When applied to the Community Earth System Model Large Ensemble, broad declines in snowstorm frequency are projected by the later 21 st century. These decreases are primarily due to a warmer atmosphere less conducive to snowfall as the predominant precipitation type. However, we note reductions are less significant for major events, since more hostile thermodynamic environments are partially offset by increased precipitation associated with cyclones that dynamically drive high-impact snowstorms.

58 GEOSCIENCES↗
DOE OSTI2019

Differences in Ice Cloud Optical Depth From CALIPSO and Ground-Based Raman Lidar at the ARM SGP and TWP Sites

Ice cloud column optical depths from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite and ground-based Raman lidars (RLs) at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) and Tropical Western Pacific (TWP) sites are compared for 8 and 4 years, respectively. The cloudy mean ice cloud column optical depths from CALIPSO Version 4 data product with a horizontal resolution of 5 km are 78% (63%) larger than those from the RLs with a resolution of 10 min/30 m at SGP (TWP) for collocated transparent profiles. The main difference at SGP is caused by the lidar ratio that for CALIPSO is related to its treatment of the multiple scattering factor. The main difference at TWP is caused by the averaging resolution. Large differences in the optical depth distribution between the CALIPSO and RL are found for small optical depths at both sites, which are due to optically very thin clouds only detectable by the RLs. The differences in ice cloud column optical depth distributions and their mean values between the CALIPSO and RL can largely be reconciled over both sites after accounting for the averaging resolutions, lidar ratios along with the CALIPSO multiple scattering factor treatment, sensitivity to optically very thin ice clouds, and definition of transparent profiles. This work also examines in detail the lidar ratios that are directly observed by the RLs, which does not support the temperature-dependent parameterizations of ice cloud lidar ratio and multiple scattering factor used in CALIPSO's Version 4 data product.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Investigation of Liquid Cloud Microphysical Properties of Deep Convective Systems: 2. Parameterization of Raindrop Size Distribution and its Application for Convective Rain Estimation

The liquid cloud microphysical properties for stratiform rain (SR) have been investigated in Part 1 of this series. Since convective rain (CR) has characteristics in raindrop size distribution (DSD) and precipitation that are distinct from SR, we investigate the CR properties in this study by using 20 hr of CR samples collected by 17 Automatic Parsivel Units disdrometers during the Midlatitude Continental Convective Clouds Experiment over the Atmospheric Radiation Measurement Southern Great Plains site. A full spectrum of DSD is constructed based on a total of 23 size channels (0.321 to 9.785 mm), and both Gamma and Exponential fitted functions are applied to extract the DSD shape parameters. Compared to SR properties, CR has distinct features including broader size range and narrower exponential slope parameter (λ E ). These results indicate that the assumption of constant exponential intercept parameter ( N 0E ) is inappropriate for CR rain rate estimates. Additionally, the subsetting scheme of the CR DSD spectra also has a strong impact on the Gamma/Exponential functions. Therefore, a new CR DSD parameterization scheme is developed by choosing the appropriate spectra subset with the constraint of rain rate and using the constant λ E instead of constant N 0E . With the input of the radar reflectivity at the lowest observed height, the newly calculated CR rain rates match well with the collocated surface rain gauge measurements (127 Mesonet stations and 17 Automatic Parsivel Units), while the rain rates calculated using traditional Z-R relationship are 3–4 times larger, indicating that constant λ E is a better assumption for CR DSD.

58 GEOSCIENCES↗
DOE OSTI2019

Analysis of Water Vapor Absorption in the Far-Infrared and Submillimeter Regions Using Surface Radiometric Measurements From Extremely Dry Locations

The second Radiative Heating in Underexplored Bands Campaign (RHUBC-II) was conducted in 2009 by the U.S. Department of Energy Atmospheric Radiation Measurement program to improve water vapor spectroscopy in the far-infrared spectral region. RHUBC-II was located in an extremely dry region of Chile to ensure very low opacities in this spectral region. Spectrally resolved measurements by a far-infrared spectrometer and a submillimeter interferometer from RHUBC-II are compared with line-by-line radiative transfer model calculations. Water vapor amounts and temperatures used in the calculations come from collocated radiosondes, with extensive adjustments to correct for issues due to the campaign's dry conditions and mountainous terrain. A reanalysis is also performed of far-infrared measurements taken at the Atmospheric Radiation Measurement North Slope of Alaska site before and during the first RHUBC campaign. These analyses determine that differences between the measurements and model calculations using existing spectroscopic parameters are significant in the far-infrared and submillimeter regions, leading to the derivation of improved water vapor continuum absorption coefficients and air-broadened widths of 74 water vapor lines. The foreign continuum is increased by more than 50% in part of the far-infrared and the widths of more than 20 lines are changed by more than 10%. The uncertainty in the foreign continuum coefficients is estimated as greater than 20% in some spectral regions, primarily a consequence of the uncertainty in the specification of water vapor. The improved far-infrared spectroscopic parameters have a notable impact on calculated spectral radiances and a modest impact on broadband radiative fluxes and heating rates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Fracture Caging to Limit Induced Seismicity

Geothermal resources offer a stable low-carbon energy source. However, geothermal sites can collocate with the hypocenters of large-magnitude seismic events. Large seismic events pose a risk to public safety and are therefore a liability for efforts to develop geothermal resources. In this work, we propose “fracture caging” to limit induced seismic event magnitudes and present evidence from numerical model predictions, laboratory experiments, and field observations. Fracture caging involves drilling tactical production wells around a geothermal injection zone to contain fluids in fracture-dominated flow systems. Prior to our work, the effect of small wells on the growth of large fractures and on flow through fractures was subject to debate. Our work shows that production wells can impede fracture growth and contain high-pressure fluids in fracture-dominated rocks. This containment offers a mechanism to limit induced seismicity.

58 GEOSCIENCES↗
DOE OSTI2020

A Satellite-Based Estimate of Convective Vertical Velocity and Convective Mass Flux: Global Survey and Comparison with Radar Wind Profiler Observations

Convective vertical velocity (w c ) and convective mass flux (M c ) lie at the heart of GCM cumulus parameterizations, but few observations of these critical parameters are available. In this paper, we develop and evaluate a novel, satellite-based method for estimating profiles of w c and M c . Here, comparisons with collocated ground-based radar wind profiler (RWP) observations show that satellite estimated median w c is slightly greater than the RWP estimates, but they show solid agreement when compared at the 95th percentiles (intense updrafts). RWP-derived and satellite estimated M c are broadly comparable in the lower and middle troposphere, with some differences in the upper troposphere due to differences in convective core sampling. A k-means cluster analysis of multiple years of w c data shows that convective characteristics are distinctly different among extratropical convection, tropical land convection, and tropical oceanic convection. Tropical land convection is significantly more intense and more variable than the oceanic counterpart.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Remote Sensing of Electric Fields Observed Within Winter Precipitation During the 2020 Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) Field Campaign

Aircraft electric fields from NASA's Lightning Instrument Package (LIP) were coupled with other airborne and ground-based, and in situ measurements to understand electrification in winter clouds that did not produce lightning. The measurements were made during seven research flights by a NASA ER-2 during the 2020 Investigation of Microphysics and Precipitation for Atlantic Coast-Threatening Snowstorms (IMPACTS) campaign. Observed total electric field magnitudes were as high as 80 V m -1 and variability in the electric field was observed along the flight path of the ER-2, indicating horizontal and/or vertical inhomogeneity in the cloud's electrical structure. X-band airborne radar data indicated 20-dBZ echo tops above 5 km in regions where electrification exceeded 10 V m -1 . In these regions, 85-GHz brightness temperatures (TB) from an airborne radiometer were lower than 265 K, with the lowest T B (~210 K) associated with ice scattering collocated with the strongest electric fields. In situ microphysical measurements from the NASA P-3 aircraft on February 7 indicated that regions near strong electric field contained supercooled water, rimed ice hydrometeors, ice water p- content as high as 1 g m -3 , liquid water content as high as 0.15 g m -3 , and supersaturation as high as 3.5%. Furthermore, these observations support the role of mixed phase microphysics in the generation of electric fields in clouds. In three case studies, ground based S-band polarimetric radar observed depolarization streaks in differential reflectivity near areas where the strongest electrification was observed. This observation reinforces the utility of depolarization streaks to identify areas of electrification prior to lightning occurrence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Performance of Cloud 3D Solvers in Ice Cloud Shortwave Radiation Closure Over the Equatorial Western Pacific Ocean

Abstract For retrieving cloud optical properties from satellite images or computing these properties from climate model output, computationally efficient treatments of cloud horizontal inhomogeneity include the Monte Carlo Independent Column Approximation (McICA) and the Tripleclouds method. Computationally efficient treatment of cloud horizontal radiation exchanges includes the SPeedy Algorithm for Radiative TrAnsfer through CloUd Sides (SPARTACUS). As a test to derive properties from satellite images, we collocate Moderate Resolution Imaging Spectroradiometer (MODIS) cloud retrievals with near‐nadir Cloud and the Earth's Radiant Energy System (CERES) footprints in July 2008 over an equatorial western Pacific Ocean region to compare the performance of the McICA, Tripleclouds, and SPARTACUS solvers to the conventional plane‐parallel homogeneous (PPH) treatment. PPH overestimates cloud albedo, and the three solvers effectively reduce overestimation with root mean square error of shortwave upwelling irradiance decreasing between 15.72 and 18.53 W m −2 , or about 22%–25%. Although cloud top variability does not get fed into the simulations, all three solvers also reduce the effect of cloud top variability on cloud albedo. Entrapment (energy reflected downward from clouds) and horizontal radiation transfer have opposite effects on the SPARTACUS cloud albedo simulation. The net effect depends on the cloud vertical extent, the unawareness of which limits the performance of the SPARTACUS solver.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The Dynamics of Megafire Smoke Plumes in Climate Models: Why a Converged Solution Matters for Physical Interpretations

As the climate system warms, megafires have become more frequent with devastating effects. A byproduct of these events is the creation of smoke plumes that can rise into the stratosphere and spread across the globe where they reside for many months. To gain a deeper understanding of the plume dynamics, global climate simulations of a megafire were performed at a wide range of grid spacings from 2.0° down to 7 km, including a 7 km nonhydrostatic experiment. The analysis focuses on how the resolved dynamics affects the specification of the plume characteristics such as injection height and black carbon (BC) mass. Prior studies initialize the smoke plume at one or a few grid points and this is shown here to produce severely dissipative dynamics. In order to validate such simulations with observations, enhancements of the plume characteristics to offset the dissipation is necessary. Using a numerically converged simulation, sensitivity tests show that to approximate the observed stratospheric lifetime, a reduction in BC fraction by 50% is necessary for external mixtures. The vorticity dynamics of the plume is also analyzed with a Lagrangian budget to understand the mechanisms responsible for the evolution of a collocated anticyclonic vortex. The results can be distilled down into a simple conceptual model. As the plume rises, the air diverges at the top of the updraft where the largest concentrations of smoke are found. This divergence induces a dilution of the background cyclonic absolute vorticity producing an anticyclonic vortex. Vortex decay occurs from opposite arguments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Advancing the Understanding of Snow Accumulation, Melting, and Associated Thermal Insulation Using Spatially Dense Snow Depth and Temperature Time Series

Snow thermal insulation is a critical factor influencing ground thermal dynamics and associated biogeochemical processes. We analyzed the spatiotemporal variability of snow accumulation, melting, and thermal insulation dynamics using spatially dense, collocated snow depth and ground interface temperature time series over two consecutive years. We demonstrated that considering late‐winter snow depth alone was insufficient to fully capture the complexity in snow and insulation dynamics. The influence of vegetation and topography on snow depth distribution varied over the season, across sites and years. We found that deep snow with a long melting period had a substantial impact on thawing n‐factors. To better predict snow insulation effects, we proposed a new weighted snow depth metric that integrates mean daily snow depth and air temperature throughout the cold season. Our results provide insights for developing space‐time remote sensing products and evaluating the representation of snow and permafrost processes in Earth system models.

54 ENVIRONMENTAL SCIENCES↗