Search NASA⌕ Search

SEARCH · Search NASA

Results for “collocated”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Tensor Network Space-Time Spectral Collocation Method for Time-Dependent Convection-Diffusion-Reaction Equations

Emerging tensor network techniques for solutions of partial differential equations (PDEs), known for their ability to break the curse of dimensionality, deliver new mathematical methods for ultra-fast numerical solutions of high-dimensional problems. Here, we introduce a Tensor Train (TT) Chebyshev spectral collocation method, in both space and time, for the solution of the time-dependent convection-diffusion-reaction (CDR) equation with inhomogeneous boundary conditions, in Cartesian geometry. Previous methods for numerical solution of time-dependent PDEs often used finite difference for time, and a spectral scheme for the spatial dimensions, which led to a slow linear convergence. Spectral collocation space-time methods show exponential convergence; however, for realistic problems they need to solve large four-dimensional systems. We overcome this difficulty by using a TT approach, as its complexity only grows linearly with the number of dimensions. We show that our TT space-time Chebyshev spectral collocation method converges exponentially, when the solution of the CDR is smooth, and demonstrate that it leads to a very high compression of linear operators from terabytes to kilobytes in TT-format, and a speedup of tens of thousands of times when compared to a full-grid space-time spectral method. These advantages allow us to obtain the solutions at much higher resolutions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Investigating Inherent Numerical Stabilization for the Moist, Compressible, Non‐Hydrostatic Euler Equations on Collocated Grids

This study investigates inherent numerical dissipation due to upwind fluxes and reconstruction strategies for collocated Finite-Volume integration of the Euler equations. Idealized supercell simulations are used without any explicit dissipation. Flux terms are split into: mass flux, pressure, and advected quantities. They are computed with the following upwind strategies: central, advectively upwind, and acoustically upwind. This is performed for third and ninth-order-accurate reconstructions with and without Weighted Essentially Non-Oscillatory limiting. Acoustic-only upwinding for pressure and mass flux terms and advective-only upwinding for advected quantities is the most flexible simplification found. It reduces data movement and computations. Assuming a constant speed of sound in acoustic upwinding gives similar results to using the true speed of sound. Dissipation from upwind adapts automatically to grid spacing, time step, reconstruction accuracy, and flow smoothness. While stability is maintained even at 21st-order spatial accuracy, there is a limit to the spatial order of accuracy for which upwinding alone can create a realizable solution in the conditions of this study. Convex combinations of upwind and central solutions for flux terms also reduced dissipation, but as the central proportion grows, solutions become physically unrealizable. The range of length scales of the kinetic energy spectra can be extended along k −5/3 to smaller spatial scales by reducing dissipation either with higher-order reconstructions or using convex combinations of upwind and central fluxes. However, not all extensions of the length scale range along k −5/3 exhibit physically realizable solutions, even though the spectra appear to be physical.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Application of Ocean Thermal Energy Conversion (OTEC) Systems for Powering Safety Monitoring Systems of Offshore Oil and Gas Operations, OESI 2.0 M-1 T-1-P1.1 – Objective 3 Report: Collocation of OTEC with Offshore Oil and Gas Activities

This report is a component of a comprehensive research initiative aimed at evaluating the technoeconomic feasibility of deploying Ocean Thermal Energy Conversion (OTEC) systems to power safety monitoring systems for offshore oil and gas operations in the Gulf of America (GOA). The overarching study has three primary objectives: summarizing existing oil and gas equipment monitoring systems, modeling OTEC as a renewable marine energy source, and exploring the collocation of OTEC power sources with offshore oil and gas activities to offset power demand. This report specifically addresses Objective 3, focusing on the integration of OTEC systems with floating oil and gas facilities.

02 PETROLEUM↗
DOE OSTI2024

Collocation methods for nonlinear differential equations on low-rank manifolds

We introduce new methods for integrating nonlinear differential equations on low-rank manifolds. These methods rely on interpolatory projections onto the tangent space, enabling low-rank time integration of vector fields that can be evaluated entry-wise. A key advantage of our approach is that it does not require the vector field to exhibit low-rank structure, thereby overcoming significant limitations of traditional dynamical low-rank methods based on orthogonal projection. To construct the interpolatory projectors, we develop a sparse tensor sampling algorithm based on the discrete empirical interpolation method (DEIM) that parameterizes tensor train manifolds and their tangent spaces with cross interpolation. Using these projectors, we propose two time integration schemes on low-rank tensor train manifolds. The first scheme integrates the solution at selected interpolation indices and constructs the solution with cross interpolation. The second scheme generalizes the well-known orthogonal projector-splitting integrator to interpolatory projectors. We demonstrate the proposed methods with applications to several tensor differential equations arising from the discretization of partial differential equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Additive Polynomial Time Integrators, Part I: Framework and Fully Implicit-Explicit Collocation Methods

In this paper we generalize the polynomial time integration framework to additively partitioned initial value problems. The framework we present is general and enables the construction of many new families of additive integrators with arbitrary order-of-accuracy and varying degree of implicitness. In this first work, we focus on a new class of implicit-explicit polynomial block methods that are based on fully implicit Runge–Kutta methods with Radau nodes and that possess high stage order. Here, we show that the new fully implicit-explicit (FIMEX) integrators have improved stability compared to existing IMEX Runge–Kutta methods, while also being more computationally efficient due to recent developments in preconditioning techniques for solving the associated systems of nonlinear equations. For PDEs on periodic domains where the implicit component is trivial to invert, we will show how parallelization of the right-hand side evaluations can be exploited to obtain significant speedup compared to existing serial IMEX Runge–Kutta methods. For parallel (in space) finite element discretizations, the new methods can achieve orders of magnitude better accuracy than existing IMEX Runge–Kutta methods and/or achieve a given accuracy several times times faster in terms of computational runtime.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Multi-link spherical joint with collocated centers of rotation

A joint includes a shell that is hollow and at least partially spherical. The shell includes a plurality of shell rail sections including a first shell rail section. The first shell rail section includes a portion of an outer surface of the shell, a portion of an inner surface of the shell, and an opening edge section that defines a shell opening. The joint also includes a plurality of rotatable members including a first rotatable member. The first rotatable member includes an outer cup positioned at least partially outside of the shell and having an outer track surface, and an inner cup positioned at least partially inside of the shell and having an inner track surface. The first shell rail section is positioned at least partially between the outer and inner track surfaces.

Hunsaker, Dakota↗
DOE OSTI2023

Cantilever with a collocated piezoelectric actuator-sensor pair

Illustrative embodiments provide an apparatus comprising a substrate comprising a cantilever, a bottom electrode on the substrate, a bottom piezoelectric transducer on the bottom electrode such that the bottom electrode is between the substrate and the bottom piezoelectric transducer, a middle electrode on the bottom piezoelectric transducer such that the bottom piezoelectric transducer is between the bottom electrode and the middle electrode, a top piezoelectric transducer on the middle electrode such that the middle electrode is between the bottom piezoelectric transducer and the top piezoelectric transducer, and a top electrode on the top piezoelectric transducer, such that the top piezoelectric transducer is between the middle electrode and the top electrode. Illustrative embodiments also provide a method of making the apparatus and a method of using the apparatus for atomic force microscopy.

Moheimani, Seyed Omid Reza↗
DOE OSTI2024

HybridSystemsSimulations.jl - Solving the Merchant Collocated Facilities with JuMP

The development of new clean-generation technologies also leads to new plant-level architectures that combine several generation and storage assets behind the point of connection. These co-located generation resources (Hybrid Systems) primarily operate as merchant assets that employ automated market bidding models and internal Energy Management Systems (EMS) to comply with the operator's signals. Formulating an optimal bidding model requires embedding the EMS control model into the bidding algorithm, resulting in a bi-level optimization problem. In this presentation, we first showcase using JuMP to formulate and solve this problem effectively for multiple merchant systems and the bidding outcomes considering different model formulations. Second, the bidding outcomes are later integrated into a PowerSimulations.jl (also built with JuMP) simulation to study the system-level effects of the various merchant bidding and the interactions between market-clearing models and the embedded EMS model. We will showcase simulations conducted in the RTS system considering different levels of merchant hybrid systems participation. The presentation provides the following specific insights on JuMP usage: 1) the Formulation of specialized bi-level problems with custom cuts to solve the merchant hybrid system bidding problem; 2) the integration of a modular model within a complex simulation workflow supported by JuMP in PowerSimulations.jl; 3) Accelerating the solution of power systems operations simulation that employ agent optimization problems using JuMP.

energy markets↗
DOE OSTI2025

Relationships Between Mesoscale Convective System Properties and Midlevel Dynamic Perturbations

Abstract Past studies implicate dynamic anomalies operating on subsynoptic scales as a possible initiation source of summertime (July–August) mesoscale convective systems (MCSs) in the central United States during northwesterly flow regimes. To improve our understanding of warm season MCSs occurring over a variety of flow regimes, we track midlevel (600 hPa) vorticity perturbations (“MPs”) as 2D objects comprising wavelengths of 500–2,500 km over the central US from May–August of 2004–2021. We perform statistical analysis of relationships between metrics of MP objects (e.g., duration, size, intensity, and origin) and high‐resolution MCS precipitation characteristics (e.g., duration, total rainfall, rain coverage area, and motion) that occur while collocated with or in the absence of MPs to discern predictive capability of background dynamic features on storm precipitation potential. Although the majority of MPs collocated with MCS initiation occur during July–August, a significant number (40%) occur between May and June. Northwesterly flow MPs comprise a relative minority of our events, suggesting that MPs can affect MCSs across a variety of warm season flow regimes. MPs affecting MCSs initiated primarily over the high plains near the central Rockies. Only approximately 20% of tracked MCS initiation events were collocated with MPs, but these storms produced ∼25% greater lifetime rainfall and coverage area, and ∼29% more stratiform rain than non‐MP‐induced MCSs. In general, larger and more vigorous MPs resulted in more hydrologically impactful MCSs. The most directly attributable benefit to MCS initiation was from MP‐enhanced background vertical motion and thermodynamic instability (e.g., increased CAPE).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaluation of the hyperspectral radiometer (HSR1) at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site

Abstract. The Peak Design Ltd hyperspectral radiometer (HSR1) was tested at the Atmospheric Radiation Measurement (ARM) user facility Southern Great Plains (SGP) site in Lamont, Oklahoma, for 2 months from May to July 2022. The HSR1 is a prototype instrument that measures total (Ftotal) and diffuse (Fdiffuse) spectral irradiance from 360 to 1100 nm with a spectral resolution of 3 nm. The HSR1 spectral irradiance measurements are compared to nearby collocated spectral radiometers, including two multifilter rotating shadowband radiometers (MFRSRs) and the Shortwave Array Spectroradiometer–Hemispheric (SASHe) radiometer. The Ftotal at 500 nm for the HSR1 compared to the MFRSRs has a mean (relative) difference of 0.01 W m−2 nm−1 (1 %–2 %). The HSR1 mean Fdiffuse at 500 nm is smaller than the MFRSRs' by 0.03–0.04 (10 %) W m−2 nm−1. The HSR1 clear-sky aerosol optical depth (AOD) is also retrieved by considering Langley regressions and compared to collocated instruments such as the Cimel sunphotometer (CSPHOT), MFRSRs, and SASHe. The mean HSR1 AOD at 500 nm is larger than the CSPHOT's by 0.010 (8 %) and larger than the MFRSRs' by 0.007–0.017 (6 %–18 %). In general, good agreement between the HSR1 and other instruments is found in terms of the Ftotal, Fdiffuse, and AODs at 500 nm. The HSR1 quantities are also compared at other wavelengths to the collocated instruments. The comparisons are within ∼ 10 % for the Ftotal and Fdiffuse, except for 940 nm, where there is relatively larger disagreement. The AOD comparisons are within ∼ 10 % at 415 and 440 nm; however, a relatively larger disagreement in the AOD comparison is found for higher wavelengths.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Wet-radome attenuation in ARM cloud radars and its utilization in radar calibration using disdrometer measurements

Abstract. A relative calibration technique has been developed for the US Department of Energy's (DOE's) Atmospheric Radiation Measurement (ARM) user facility Ka-band ARM Zenith Radars (KAZRs). This method uses the signal attenuation caused by water on the radome to estimate reflectivity factor (Ze) offsets. The wet-radome attenuation (WRA) is assumed to follow a log-linear relationship with rainfall rate during light and moderate rain, as measured by a collocated surface disdrometer. The technique has an uncertainty of approximately 3 dB, due to factors such as disdrometer measurement error, rain variability between radar and disdrometer sample volumes, and the fitting function's uncertainty for the WRA behavior. A practical advantage of this WRA-based approach to shorter-wavelength radar monitoring is that, while it requires a reference disdrometer, it proves feasible for a wider range of collocated disdrometer measurements compared to traditional direct disdrometer comparison at the onset of light rain. This technique thus offers a cost-effective monitoring tool for remote or long-term radar deployments. This calibration technique was applied during the ARM Tracking Aerosol Convection Interactions Experiment (TRACER) from October 2021 through September 2022. The estimated Ze offsets were compared against traditional radar calibration and monitoring methods using available datasets from this campaign. Results show that the WRA-based offsets align closely with mean offsets found between cloud radars and from direct disdrometer comparison near the onset of rain, while also reflecting similar offset and campaign-long trends when compared to collocated, independently calibrated radar wind profilers. Nevertheless, overall, the KAZR Ze offsets estimated during TRACER remained stable at approximately 2 dB lower than the disdrometer estimates from the campaign start until the end of June 2022; afterward, the offsets increased to around 7 dB by the campaign's end. This increase is linked to a drop of about 1 dB in transmitter power toward the end of the project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Marine and continental stratocumulus cloud microphysical properties obtained from routine ARM Cimel sunphotometer observations

This study investigates marine and continental stratocumulus (Sc) cloud properties obtained from an automated implementation of a multispectral photometer retrieval. Photometer methods simultaneously retrieve cloud optical depth (τ) and cloud droplet effective radius (r e ), with estimates for liquid water path (LWP) calculated on the availability of those quantities. These applied methods evaluate retrieved cloud properties for Sc identified during a recent 6 year period over the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) program sites in Oklahoma, USA (SGP) and in the Azores, Portugal (ENA). Modest agreement in key quantity retrievals is found between the routine photometer products and multisensor collocated profiling references. Cumulative breakdowns contingent on cloud thickness indicate increases in all retrieved quantities in thicker clouds, with larger discrepancies in the relative performance between the retrievals collected in the presence of drizzle. Under continental cloud conditions, the clouds of a similar thickness and r e to those sampled under marine conditions report a factor of 1.5 larger τ and LWP. An r 2 ≅0.65 is found between photometer τ retrievals and shadowband radiometer measurements, with photometer retrievals reporting a high (relative) bias. The τ intercomparisons indicate that variability between retrievals is a factor of three larger than errors reported from individual retrieval input perturbation tests. Photometer r e retrievals suggest a low r 2 (< 0.1) having a standard deviation ≅ 3 µm when compared to ARM baseline multi-sensor radar/radiometer references (accounting for offsets in the cloud droplet number concentration assumptions of the latter). However, photometer LWP calculations remain relatively unbiased in non-drizzling conditions, with errors O (50 g m −2 ) and r 2 ≅0.5 to collocated radiometer and interferometer references. Additional sensitivity tests for island influences on marine Sc properties suggest that while island-influenced winds may promote larger cloud LWP or thickness, the influence could be within retrieval method uncertainty and/or collocated instrument variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Constrained or unconstrained? Neural-network-based equation discovery from data

Throughout many fields, practitioners often rely on differential equations to model systems. Yet, for many applications, the theoretical derivation of such equations and/or the accurate resolution of their solutions may be intractable. Instead, recently developed methods, including those based on parameter estimation, operator subset selection, and neural networks, allow for the data-driven discovery of both ordinary and partial differential equations (PDEs), on a spectrum of interpretability. The success of these strategies is often contingent upon the correct identification of representative equations from noisy observations of state variables and, as importantly and intertwined with that, the mathematical strategies utilized to enforce those equations. Specifically, the latter has been commonly addressed via unconstrained optimization strategies. Representing the PDE as a neural network, we propose to discover the PDE (or the associated operator) by solving a constrained optimization problem and using an intermediate state representation similar to a physics-informed neural network (PINN). The objective function of this constrained optimization problem promotes matching the data, while the constraints require that the discovered PDE is satisfied at a number of spatial collocation points. We present a penalty method and a widely used trust-region barrier method to solve this constrained optimization problem, and we compare these methods on numerical examples. Our results on several example problems demonstrate that the latter constrained method outperforms the penalty method, particularly for higher noise levels or fewer collocation points. This work motivates further exploration into using sophisticated constrained optimization methods in scientific machine learning, as opposed to their commonly used, penalty-method or unconstrained counterparts. For both of these methods, we solve these discovered neural network PDEs with classical methods, such as finite difference methods, as opposed to PINNs-type methods relying on automatic differentiation. Here, we briefly highlight how simultaneously fitting the data while discovering the PDE improves the robustness to noise and other small, yet crucial, implementation details.

Data-driven discovery↗
DOE OSTI2025

Gradient flow based phase-field modeling using separable neural networks

Allen–Cahn equation is a reaction–diffusion equation and is widely used for modeling phase separation. Machine learning methods for solving the Allen–Cahn equation in its strong form suffer from inaccuracies in collocation techniques, errors in computing higher-order spatial derivatives, and the large system size required by the space–time approach. To overcome these challenges, we propose solving the gradient flow of the Ginzburg–Landau free energy functional, which is equivalent to the Allen–Cahn equation, thereby avoiding the second-order spatial derivatives associated with the Allen–Cahn equation. A minimizing movement scheme is employed to solve the gradient flow problem, eliminating the complexities of a space–time approach. We utilize a separable neural network that efficiently represents the phase field through low-rank tensor decomposition. As we use the minimizing movement scheme to numerically solve the gradient flow problem, we thus, refer to the proposed method as the Separable Deep Minimizing Movement (SDMM) method. The evaluation of the functional in the minimizing movement scheme using the Gauss quadrature technique bypasses the inaccuracies associated with collocation techniques traditionally used to solve partial differential equations. A hyperbolic tangent transformation is introduced on the phase field prior to the evaluation of the functional to ensure that it remains strictly bounded within the values of the two phases. For this transformation, theoretical guarantee for energy stability of the minimizing movement scheme is established. Our results suggest that this transformation helps to improve the accuracy and efficiency significantly. The proposed method resolves the challenges faced by state-of-the-art machine learning techniques, outperforming them in both accuracy and efficiency. It is also the first machine learning method to achieve an order of magnitude speed improvement over the finite element method. In addition to its formulation and computational implementation, several case studies illustrate the applicability of the proposed method.

42 ENGINEERING↗