Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2025

Hemispheric Asymmetry of Phase Partition in Mixed‐Phase Clouds Based on Near Global‐Scale Airborne Observations

Mixed-phase clouds contribute to substantial uncertainties in global climate models due to their complex microphysical properties. Former model evaluations almost exclusively rely on satellite observations to assess cloud phase distributions globally. This study investigated mixed-phase cloud properties using near global-scale in situ observation data sets from 14 flight campaigns in combination with collocated output from a global climate model. The Southern Hemisphere (SH) shows significantly higher occurrence frequencies and higher mass fractions of supercooled liquid water than Northern Hemisphere (NH) based on observations at 0.2 and 100 km horizontal scales. Such hemispheric asymmetry is not captured by the model. The model also consistently overestimates liquid water content (LWC) in all cloud phases but shows ice water content (IWC) biases that vary with phase. Key processes contributing to model biases in phase partition can be identified through the combination of evaluation of phase frequency, liquid mass fraction, LWC and IWC.

Yang, Ching An [San Jose State University, CA (Uni↗
DOE OSTI2026

Precipitation‐Buoyancy Relationships in the Life Cycle of Tropical Mesoscale Convective Systems

This study aims to establish process-level benchmarks linking Mesoscale Convective Systems (MCSs) at various stages of their life cycle to their thermodynamic environment. The relationship between MCS precipitation and an empirical buoyancy measure (B 𝐿 ) is examined using collocated satellite-observed MCS tracks and reanalysis data. A positive relationship is identified between the frequency of tropical MCSs and that of high B 𝐿 conditions. The buoyancy measure, integrating instability and entrainment, helps elucidate thermodynamic characteristics throughout the MCS life cycle. Environments with high instability and moderate subsaturation are frequently linked to the initial stage, while environments with low instability and near saturation are frequently linked to the mature stage. Stable and highly subsaturated environments are more likely associated with the termination of the life cycle. These associations are qualitatively similar for oceanic and land MCSs. Overall, the MCS-environment relationships can serve as observational benchmarks with which to diagnose MCS-resolving models.

Tsai, Wei‐Ming [University of California, Los Ange↗
DOE OSTI2021

Significance of the organic aerosol driven climate feedback in the boreal area

Aerosol particles cool the climate by scattering solar radiation and by acting as cloud condensation nuclei. Higher temperatures resulting from increased greenhouse gas levels have been suggested to lead to increased biogenic secondary organic aerosol and cloud condensation nuclei concentrations creating a negative climate feedback mechanism. Here, we present direct observations on this feedback mechanism utilizing collocated long term aerosol chemical composition measurements and remote sensing observations on aerosol and cloud properties. Summer time organic aerosol loadings showed a clear increase with temperature, with simultaneous increase in cloud condensation nuclei concentration in a boreal forest environment. Remote sensing observations revealed a change in cloud properties with an increase in cloud reflectivity in concert with increasing organic aerosol loadings in the area. The results provide direct observational evidence on the significance of this negative climate feedback mechanism.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Highly active, ultra-low loading single-atom iron catalysts for catalytic transfer hydrogenation

Highly effective and selective noble metal-free catalysts attract significant attention. Here, a single-atom iron catalyst is fabricated by saturated adsorption of trace iron onto zeolitic imidazolate framework-8 (ZIF-8) followed by pyrolysis. Its performance toward catalytic transfer hydrogenation of furfural is comparable to state-of-the-art catalysts and up to four orders higher than other Fe catalysts. Isotopic labeling experiments demonstrate an intermolecular hydride transfer mechanism. First principles simulations, spectroscopic calculations and experiments, and kinetic correlations reveal that the synthesis creates pyrrolic Fe(II)-plN 3 as the active center whose flexibility manifested by being pulled out of the plane, enabled by defects, is crucial for collocating the reagents and allowing the chemistry to proceed. The catalyst catalyzes chemoselectively several substrates and possesses a unique trait whereby the chemistry is hindered for more acidic substrates than the hydrogen donors. This work paves the way toward noble-metal free single-atom catalysts for important chemical reactions.

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

Offshore wind and wave energy can reduce total installed capacity required in zero-emissions grids

Abstract As the world races to decarbonize power systems to mitigate climate change, the body of research analyzing paths to zero emissions electricity grids has substantially grown. Although studies typically include commercially available technologies, few of them consider offshore wind and wave energy as contenders in future zero-emissions grids. Here, we model with high geographic resolution both offshore wind and wave energy as independent technologies with the possibility of collocation in a power system capacity expansion model of the Western Interconnection with zero emissions by 2050. In this work, we identify cost targets for offshore wind and wave energy to become cost effective, calculate a 17% reduction in total installed capacity by 2050 when offshore wind and wave energy are fully deployed, and show how curtailment, generation, and transmission change as offshore wind and wave energy deployment increase.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Quantifying local and global mass balance errors in physics-informed neural networks

Physics-informed neural networks (PINN) have recently become attractive for solving partial differential equations (PDEs) that describe physics laws. By including PDE-based loss functions, physics laws such as mass balance are enforced softly in PINN. This paper investigates how mass balance constraints are satisfied when PINN is used to solve the resulting PDEs. We investigate PINN’s ability to solve the 1D saturated groundwater flow equations (diffusion equations) for homogeneous and heterogeneous media and evaluate the local and global mass balance errors. We compare the obtained PINN’s solution and associated mass balance errors against a two-point finite volume numerical method and the corresponding analytical solution. We also evaluate the accuracy of PINN in solving the 1D saturated groundwater flow equation with and without incorporating hydraulic heads as training data. We demonstrate that PINN’s local and global mass balance errors are significant compared to the finite volume approach. Tuning the PINN’s hyperparameters, such as the number of collocation points, training data, hidden layers, nodes, epochs, and learning rate, did not improve the solution accuracy or the mass balance errors compared to the finite volume solution. Mass balance errors could considerably challenge the utility of PINN in applications where ensuring compliance with physical and mathematical properties is crucial.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Multimodal quantification of degradation pathways during extreme fast charging of lithium-ion batteries

Enabling fast charging of Li-ion batteries will be a key step towards realizing the technology's full potential in electric vehicles. Currently, fast charging is limited by a variety of processes that reduce cell capacity upon extended cycling. Using a multimodal approach combining incremental capacity analysis (dQ/dV), high energy X-ray diffraction (HEXRD), and mass spectrometry titration (MST), we identify specific degradation mechanisms—including Li plating, dead Li x C 6 formation, Li 2 C 2 formation, solid carbonate solid-electrolyte interphase (SEI) deposition, and loss of positive electrode active material (LAM PE )—that occur during extended fast charge cycling. We find that Li plating is the major source of capacity loss in cells cycled at 6C, while non-carbonate SEI species deposition on the graphite anode is the main source of capacity loss when cycled at 4C. We also study local degradative phenomena by examining specific ~1–5 cm 2 regions of the cells using HEXRD and MST. Here, we find that plated Li is often collocated with dead Li x C 6 , Li 2 C 2 , and solid carbonate SEI species, and these additional species cumulatively account for ~20% of the capacity lost during 6C cycling. Finally, in a cell with an anomalously high amount of LAM PE (quantified via dQ/dV), we find that regions of cathode degradation were accompanied by non-carbonate SEI products on the adjacent region of the anode. We postulate that this phenomenon arises due to crosstalk between the electrodes, wherein soluble electrolyte oxidation products formed at the delithiated cathode migrate to the graphite anode and are ultimately deposited on the graphite surface. This work demonstrates the utility of combining multiple characterization techniques to reveal a more holistic understanding of degradative phenomena that occur across multiple length scales during fast charge.

25 ENERGY STORAGE↗
DOE OSTI2021

Well-balanced treatment of gravity in astrophysical fluid dynamics simulations at low Mach numbers

We report accurate simulations of flows in stellar interiors are crucial to improving our understanding of stellar structure and evolution. Because the typically slow flows are merely tiny perturbations on top of a close balance between gravity and the pressure gradient, such simulations place heavy demands on numerical hydrodynamics schemes. We demonstrate how discretization errors on grids of reasonable size can lead to spurious flows orders of magnitude faster than the physical flow. Well-balanced numerical schemes can deal with this problem. Three such schemes were applied in the implicit, finite-volume SEVEN-LEAGUE HYDRO code in combination with a low-Mach-number numerical flux function. We compare how the schemes perform in four numerical experiments addressing some of the challenges imposed by typical problems in stellar hydrodynamics. We find that the α-β and deviation well-balancing methods can accurately maintain hydrostatic solutions provided that gravitational potential energy is included in the total energy balance. They accurately conserve minuscule entropy fluctuations advected in an isentropic stratification, which enables the methods to reproduce the expected scaling of convective flow speed with the heating rate. The deviation method also substantially increases accuracy of maintaining stationary orbital motions in a Keplerian disk on long timescales. The Cargo–LeRoux method fares substantially worse in our tests, although its simplicity may still offer some merits in certain situations. Overall, we find the well-balanced treatment of gravity in combination with low Mach number flux functions essential to reproducing correct physical solutions to challenging stellar slow-flow problems on affordable collocated grids.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Stability of the steady three-dimensional lid-driven flow in a cube and the supercritical flow dynamics

The stability of the steady flow in a lid-driven cube is investigated by a collocation method making use of asymptotic solutions for the singular edges of the cavity up- and downstream of the moving wall. Owing to the rapid convergence of the method high-accuracy critical data are obtained. To determine the critical point subcritical growth rates of small perturbations are extrapolated to zero. We find the bifurcation to be of Hopf-type and slightly subcritical. Above the critical point, the oscillatory flow is symmetric with respect to the symmetric midplane of the cavity and characterized by nearly streamwise vortices in the boundary layer on the wall upstream of the moving wall. The oscillation amplitude grows slowly and seems to saturate. On a long time scale, however, the constant-amplitude oscillations are unstable. The periodic oscillations are interrupted by short bursts during which the oscillation amplitude grows substantially and the spatial structure of the oscillating streamwise vortices changes. Towards the end of each burst the mirror symmetry of the oscillatory flow is lost, the flow returns to the vicinity of the unstable steady state and the growth of symmetric oscillations starts again leading to an intermittent chaotic flow.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Transitional regimes of natural convection in a differentially heated cubical cavity under the effects of wall and molecular gas radiation

The transition to unsteadiness and the dynamics of weakly turbulent natural convection, coupled to wall or gas radiation in a differentially heated cubical cavity with adiabatic lateral walls, are studied numerically. The working fluid is air with small contents of water vapor and carbon dioxide whose infrared spectral radiative properties are modelled by the absorption distribution function model. A pseudo spectral Chebyshev collocation method is used to solve the flow field equations and is coupled to a direct ray tracing method for radiation transport. Flow structures are identified by means of either the proper orthogonal decomposition or the dynamic mode decomposition methods. We first retrieve the classical mechanism of transition to unsteadiness without radiation, characterized by counter-rotating streamwise-oriented vortices generated at the exit of the vertical boundary layers. Wall radiation through a transparent medium leads to a homogenization of lateral wall temperatures and the resulting transition mechanism is similar to that obtained with perfectly conducting lateral walls. The transition is due to an unstable stratification upstream the vertical boundary layers and is characterized by periodically oscillating transverse rolls of axis perpendicular to the main flow. When molecular gas radiation is accounted for, no periodic solution is found and the transition to unsteadiness displays complex structures with chimneys-like rolls whose axes are again parallel to the main flow. The origin of this instability is probably due to centrifugal forces, as suggested previously for the case without radiation. Above the transition to unsteadiness, at Ra = 3 × 10{sup 8}, it is shown that both wall and gas radiation significantly intensify turbulent fluctuations, decrease the thermal stratification in the core of the cavity, and increase the global circulation.

42 ENGINEERING↗
DOE OSTI2014

Phase-field modeling of switchable diode-like current-voltage characteristics in ferroelectric BaTiO{sub 3}

A self-consistent model has been proposed to study the switchable current-voltage (I-V) characteristics in Cu/BaTiO{sub 3}/Cu sandwiched structure combining the phase-field model of ferroelectric domains and diffusion equations for ionic/electronic transport. The electrochemical transport equations and Ginzburg-Landau equations are solved using the Chebyshev collocation algorithm. We considered a single parallel plate capacitor configuration which consists of a single layer BaTiO{sub 3} containing a single tetragonal domain orientated normal to the plate electrodes (Cu) and is subject to a sweep of ac bias from −1.0 to 1.0 V at 25 °C. Our simulation clearly shows rectifying I-V response with rectification ratios amount to 10{sup 2}. The diode characteristics are switchable with an even larger rectification ratio after the polarization direction is flipped. The effects of interfacial polarization charge, dopant concentration, and dielectric constant on current responses were investigated. The switchable I-V behavior is attributed to the polarization bound charges that modulate the bulk conduction.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Radiation effects on bifurcation and dual solutions in transient natural convection in a horizontal annulus

Transitions and bifurcations of transient natural convection in a horizontal annulus with radiatively participating medium are numerically investigated using the coupled lattice Boltzmann and direct collocation meshless (LB-DCM) method. As a hybrid approach based on a common multi-scale Boltzmann-type model, the LB-DCM scheme is easy to implement and has an excellent flexibility in dealing with the irregular geometries. Separate particle distribution functions in the LBM are used to calculate the density field, the velocity field and the thermal field. In the radiatively participating medium, the contribution of thermal radiation to natural convection must be taken into account, and it is considered as a radiative term in the energy equation that is solved by the meshless method with moving least-squares (MLS) approximation. The occurrence of various instabilities and bifurcative phenomena is analyzed for different Rayleigh number Ra and Prandtl number Pr with and without radiation. Then, bifurcation diagrams and dual solutions are presented for relevant radiative parameters, such as convection-radiation parameter Rc and optical thickness τ. Numerical results show that the presence of volumetric radiation changes the static temperature gradient of the fluid, and generally results in an increase in the flow critical value. Besides, the existence and development of dual solutions of transient convection in the presence of radiation are greatly affected by radiative parameters. Finally, the advantage of LB-DCM combination is discussed, and the potential benefits of applying the LB-DCM method to multi-field coupling problems are demonstrated.

42 ENGINEERING↗
DOE OSTI2015

Numerical simulation of scattering of acoustic waves by inelastic bodies using hypersingular boundary integral equation

A numerical method for solving a problem of diffraction of acoustic waves by system of solid and thin objects based on the reduction the problem to a boundary integral equation in which the integral is understood in the sense of finite Hadamard value is proposed. To solve this equation we applied piecewise constant approximations and collocation methods numerical scheme. The difference between the constructed scheme and earlier known is in obtaining approximate analytical expressions to appearing system of linear equations coefficients by separating the main part of the kernel integral operator. The proposed numerical scheme is tested on the solution of the model problem of diffraction of an acoustic wave by inelastic sphere.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗