Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time distribution”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Plasma thermal transport with a generalized 8-moment distribution function

Moment equations that model plasma transport require an ansatz distribution function to close the system of equations. The resulting transport is sensitive to the specific closure used, and several options have been proposed in the literature. Two different 8-moment distribution functions can be generalized to form a single-parameter family of distribution functions. The transport coefficients resulting from this generalized distribution function can be expressed in terms of this free parameter. This provides the flexibility of matching the 8-moment model to some validating result at a given magnetization value, such as Braginskii’s transport, or the more recent results of Davies et al. [Physics of Plasma, 28, 012305 (2021)]. Here, this process can be thought of as solving for the 8-moment distribution function that matches the value of a transport coefficient given by a Chapman-Enskog expansion, while retaining the improved physical properties, such as finite propagation speeds and time dependence which belong to the hyperbolic moment models. Since the presented generalized distribution function only has a single free parameter, only a single transport coefficient can be matched at a time. However, this generalization process may be extended to provide multiple free parameters. The focus of this Brief Communication is on the dramatically improved thermal conductivity of the proposed model compared to the two base moment models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The Kimberlina synthetic multiphysics dataset for CO 2 monitoring investigations

Abstract We present a synthetic multi‐scale, multi‐physics dataset constructed from the Kimberlina 1.2 CO 2 reservoir model based on a potential CO 2 storage site in the Southern San Joaquin Basin of California. Among 300 models, one selected reservoir‐simulation scenario produces hydrologic‐state models at the onset and after 20 years of CO 2 injection. Subsequently, these models were transformed into geophysical properties, including P‐ and S‐wave seismic velocities, saturated density where the saturating fluid can be a combination of brine and supercritical CO 2 , and electrical resistivity using established empirical petrophysical relationships. From these 3D distributions of geophysical properties, we have generated synthetic time‐lapse seismic, gravity and electromagnetic responses with acquisition geometries that mimic realistic monitoring surveys and are achievable in actual field situations. We have also created a series of synthetic well logs of CO 2 saturation, acoustic velocity, density and induction resistivity in the injection well and three monitoring wells. These were constructed by combining the low‐frequency trend of the geophysical models with the high‐frequency variations of actual well logs collected at the potential storage site. In addition, to better calibrate our datasets, measurements of permeability and pore connectivity have been made on cores of Vedder Sandstone, which forms the primary reservoir unit. These measurements provide the range of scales in the otherwise synthetic dataset to be as close to a real‐world situation as possible. This dataset consisting of the reservoir models, geophysical models, simulated time‐lapse geophysical responses and well logs forms a multi‐scale, multi‐physics testbed for designing and testing geophysical CO 2 monitoring systems as well as for imaging and characterization algorithms. The suite of numerical models and data have been made publicly available for downloading on the National Energy Technology Laboratory's (NETL) Energy Data Exchange (EDX) website.

58 GEOSCIENCES↗
DOE OSTI2020

Ring-opening and hydrodenitrogenation of indole under hydrothermal conditions over Ni, Pt, Ru, and Ni-Ru bimetallic catalysts

Here, the activity and selectivity of activated-carbon-supported Ni, Pt, Ru, and Ni-Ru bimetallic catalysts was examined for hydrothermal denitrogenation of indole. The molar yield of pyrrole ring-opening compounds, without an added hydrogen source, are in the order: Ni < Pt < Ni 90 Ru 10 < Ni 75 Ru 25 < Ni 50 Ru 50 ~ Ni 25 Ru 75 ~ Ru. Ru-containing catalysts facilitated production of hydrocarbons (hydrodenitrogenation (HDN) products) when used with added formic acid (hydrogen source). We elucidated catalytic hydrothermal HDN pathways for indole based on the product distributions and the variation of their yields with time. Hydrogenation of indole to indoline is the primary pathway and ring-opening of indoline to form alkyl anilines is faster than forming HDN products (alkyl benzenes). DFT calculations confirmed experimental activity trends, showing Ru is more active than Ni for indole ring opening and for o-toluidine deamination. If no hydrogen source is present, directly breaking the N-C bond in the pyrrole ring is more favorable than breaking the C-N bond with an aromatic carbon. If a H source is provided, the pyrrole ring hydrogenates first, forming indoline, followed by cleavage of the C-N bond.

42 ENGINEERING↗
DOE OSTI2022

Qualitative trend analysis based on a mixed-integer representation

Shape constrained spline fitting is a useful method to impose prior knowledge onto flexible semi-parametric models during parameter estimation. Most typically, the function shape is imposed through order restrictions on the regression coefficients. The intended shape is considered known or selected based on heuristic rules. In this study, we present a method to estimate the optimal set of order restrictions to segment a univariate data series into episodes with distinct shapes. This is also known as the qualitative trend analysis (QTA) problem. The obtained solution uses a trade-off between lack-of-fit and model complexity. Further, our practical implementation takes inspiration from the generalized order restricted information criterion (GORIC) for inequality-constrained model selection. From this, one learns (a) that QTA can be formulated as a mixed-integer quadratic program (MIQP) and (b) that the newly proposed mixed order restricted information criterion (MORIC) enables optimal segmentation. This is illustrated through didactic case studies.

42 ENGINEERING↗
DOE OSTI2022

A nuclear-hydrogen hybrid energy system with large-scale storage: A study in optimal dispatch and economic performance in a real-world market

Here a study in optimal dispatch and economic analysis of a novel nuclear hybrid energy system (NHES) with large-scale hydrogen storage and a novel control scheme is conducted. The control scheme makes charge and discharge decisions by using the real-time electricity price relative to the historical price distribution. The pricing data is taken from the ERCOT Houston and West load zones with different price oscillations. Six case control studies with different levels of aggressiveness were chosen for each load zone for a total of 12 case studies. It is found that in a high price volatility load zone, a more aggressive control strategy maximizes revenue, while in a less oscillatory load zone, a moderate control strategy performs better. The NHES can store energy when the price is low and sell additional electricity when the price is high, as well as participate in the ancillary services market. Due to these advantages, the NHES can generate 10–40 million USD more revenue annually than a stand-alone nuclear plant. Due to higher costs, the NHES has a higher LCOE (81.36 USD/MWh) than a stand-alone plant (68.66 USD/MWh); however, the LCOE of the NHES is still comparable to other new forms of electricity generation. It is found that both the NHES and the stand-alone plant are not economical in the existing electricity market. Due to the extra revenue generated, the hybrid system has a positive annual cashflow whereas the stand-alone plant does not. It is found that the hybrid plant needs less cost subsidy than a stand-alone plant to become economical.

08 HYDROGEN↗
DOE OSTI2022

In situ characterization of phase transformations in petroleum pitch by high temperature X-ray diffraction

The phase transformation from isotropic to mesophase of a petroleum-derived pitch was monitored by in situ X-ray diffraction at 410 °C. The kinetics of the transformation were characterized by monitoring the growth of a broad peak at 2θ~24° and could be described by Avrami’s equation. The volatilization of low molecular weight components originally present in the pitch, or produced during the condensation reaction, was assessed by thermal gravimetry analysis, following the same thermal history. Ex situ characterizations of the pitch precursor before and after thermal treatment provided insight into the evolution of the molecular structure associated with mesophase generation (optical polarized microscopy, molecular weight distribution by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry, X-ray photoelectron spectroscopy, Auger spectroscopy, and thermal stability). Here, this investigation demonstrates the value of using complementary in situ and ex situ experimental techniques to characterize complex mechanisms during the processing of organic materials.

02 PETROLEUM↗
DOE OSTI2023

A conservative Galerkin solver for the quasilinear diffusion model in magnetized plasmas

We propose a conservative Galerkin scheme for the quasilinear model in three-dimensional momentum space and three-dimensional spectral space, with cylindrical symmetry. We construct an unconditionally conservative weak form and use a discretization that preserves conservation properties independent of the wave emission probability. The discrete operators, combined with a consistent quadrature rule, preserve all the conservation laws rigorously. The proposed scheme is quite general: it works for both relativistic and non-relativistic systems, for both magnetized and unmagnetized plasmas, and even for problems with time-dependent dispersion relations. We represent the particle distribution by continuous basis functions and use discontinuous basis functions for the wave spectral energy density, which enables the application of a positivity-preserving technique. We adopt the marching simplex algorithm, designed initially for computer graphics, for numerical integration on the resonance manifold. Furthermore, the numerical examples with a bump-on-tail initial configuration show how the unstable waves produce strong momentum space diffusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Optimization of La 2 NiO 4+δ Electrolysis Cell Oxygen Electrode through Surfactant-Enabled LaCoO 3±δ Nanocatalyst Deposition

Lanthanum nickelate (LNO) has shown promise as a Cr-resistant air electrode material for SOECs but has suboptimal surface oxygen exchange properties. Nanocoating of the LNO surface with lanthanum cobaltite (LCO) was chosen to improve cell performance as a surface oxygen conductor. The work focused on the implementation of a two-step nano-LCO film deposition utilizing catechol molecules in a porous LNO electrode. The subgoals of the work were to maintain nanosized LCO particles/ grains to increase active surface area and to control the regularity/ homogeneity of the coating across the microstructure. To achieve these goals, a novel surfactant-enhanced liquid infiltration method was utilized, where nucleation sites were spread across the electrode structure to control the location and size of LCO particles. Various catechol surfactant compositions were evaluated for their ability to control the kinetics of nanoparticle deposition and the homogeneity of the coating. Chelated LCO was characterized by X-ray diffraction (XRD), which found a substantial improvement in LCO formation with surfactant addition and determined polymerized norepinephrine to be the best-performing surfactant, with 88.4% pure LCO formed at low temperature. X-ray photoelectron spectroscopy (XPS) confirmed LCO nanostructures formed by the two-step infiltration process, showing no impurities and a stable perovskite structure. Deposition kinetics were analyzed using atomic force microscopy (AFM), correlating infiltration times and solution molarity to nanoparticle size and distribution, the results of which were confirmed in symmetrical cell samples by scanning electron microscopy (SEM). Electrochemical impedance spectroscopy (EIS) testing demonstrated substantial improvements in polarization resistance, where the nanocoating reduced the resistance by ∼55% to 0.152 Ω·cm 2 at 700 °C and 0.039 Ω·cm 2 at 800 °C. Electrical conductivity relaxation (ECR) at this temperature confirmed an improved surface oxygen exchange coefficient of the LCO + LNO heterostructure predicted by the Bode data from EIS, alongside a reduction in activation energy by about 30%.

Deposition↗
DOE OSTI2020

Simultaneous Observations of Electromagnetic Ion Cyclotron (EMIC) Waves and Pitch Angle Scattering During a Van Allen Probes Conjunction

On 22 December 2015, the two Van Allen Probes observed two sets of electromagnetic ion cyclotron (EMIC) wave bursts during a close conjunction when both Probe A and Probe B were separated by 0.57 to 0.68 R E . Here, the EMIC waves occurred during an active period in the recovery phase of a coronal mass ejection-driven geomagnetic storm. Both spacecraft observed EMIC wave bursts that had similar spatial structure within a 1–2 min time delay. The EMIC waves occurred outside the plasmasphere, within ΔL ≈ 1–2 of the plasmapause and within a few degrees in magnetic latitude of the equatorial plane. The spatial structure of the EMIC wave bursts may have been related to the proton drift paths outside the plasmasphere and influenced by total magnetic field strength variations associated with solar wind pressure enhancements. The EMIC waves were observed in a narrow L shell region from L ≈ 4.55–5.32 between 10 and 11 magnetic local time (MLT) on the outbound halves of the spacecraft orbits and from L ≈ 4.82–5.51 between 13 and 14 MLT on the inbound halves of the spacecraft orbits. However, Pc1 pulsations were observed on the ground over a broad range of local times. The anisotropy of the proton pitch angle distributions was enhanced when the EMIC waves were observed. Although the overall radiation belt response during this storm was dominated by acceleration and transport processes, the EMIC waves produced local pitch angle scattering of 13–15 keV protons and 2.1–2.6 MeV electrons, consistent with calculations of the expected resonant energies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Energy Budget Constraints on the Time History of Aerosol Forcing and Climate Sensitivity

Abstract An observationally constrained time series of historical aerosol effective radiative forcing (ERF) from 1750 to 2019 is developed in this study. We find that the time history of aerosol ERFs diagnosed in CMIP6 models exhibits considerable variation and explore how the time history of aerosol forcing influences the probability distributions of present‐day aerosol forcing and emergent metrics such as climate sensitivity. Using a simple energy balance model, trained on CMIP6 climate models and constrained by observed near‐surface warming and ocean heat uptake, we derive estimates for the historical aerosol forcing. We find 2005–2014 mean aerosol ERF to be −1.1 (−1.8 to −0.5) W m −2 relative to 1750. Assuming recently published historical emissions from fossil fuel and industrial sectors and biomass burning emissions from SSP2‐4.5, aerosol ERF in 2019 is −0.9 (−1.5 to −0.4) W m −2 . There is a modest recovery in aerosol forcing (+0.025 W m −2 decade −1 ) between 1980 and 2014. This analysis also gives a 5%–95% range of equilibrium climate sensitivity of 1.8°C –5.1°C (best estimate 3.1°C) with a transient climate response of 1.2°C –2.6°C (best estimate 1.8°C).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Using TRMM Latent Heat as a Source to Estimate Convection Induced Gravity Wave Momentum Flux in the Lower Stratosphere

Due to a lack of observations, it is a big challenge to quantify the momentum flux transported from the troposphere into the lower stratosphere by gravity waves generated by convection. This limits our understanding of the dynamics of general circulation in the stratosphere and the realism of simulations of the Quasi Biennial Oscillation. In the past decade, some general circulation models have linked the momentum flux to the latent heating from subgrid scale convective precipitation with some success. However, there is still a large uncertainty in the sources of subgrid scale convection in these models. This study applies the links between the momentum flux and convective latent heating directly to the convective precipitation derived from 16 years of TRMM precipitation radar observations in the tropics and subtropics. The total and directional momentum flux at 100 hPa are derived at individual convective pixels by using the TRMM Spectrum Latent Heating product and large-scale wind profiles from the ERA-Interim dataset. For the first time, we are able to estimate the geographical distribution of momentum flux at 100 hPa in tropics and subtropics from the observed convective sources. The diurnal, seasonal, and interannual variations of the derived momentum flux are presented and discussed. These results could provide a reference to validate the dynamic coupling between the tropical troposphere and the stratosphere in the general circulation models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Impact of Coastal Marsh Eco‐Geomorphologic Change on Saltwater Intrusion Under Future Sea Level Rise

Abstract Coastal saltwater intrusion (SWI) is one key factor that affects the hydrology, ecology, and biogeochemistry of coastal ecosystems. Future climate change, especially intensified sea level rise (SLR), is expected to trigger SWI to encroach on coastal freshwater aquifers more intensively. Numerous studies have investigated decadal/century scale SWI under SLR by assuming a static coastal landscape topography. However, coastal landscapes are highly dynamic in response to SLR, and the impact of coastal landscape evolution on SWI has received very little attention. Thus, this study used a coastal marsh landscape as an example and investigated how coastal marsh evolution affects future SWI with a physically‐based coastal hydro‐eco‐geomorphologic model, Advanced Terrestrial Simulator. Our numerical experiments showed that it is very likely that the marsh elevation increases with future SLR due to sediment deposition, and a depression zone is formed due to different marsh accretion rates between the ocean boundary and the inland. We found that marsh accretion may significantly reduce the surface saltwater inflow at the ocean boundary, and the evolved topographic depression zone may prolong the residence time of surface ponded saltwater, affecting subsurface salinity distribution differently. We also predicted that marshlands may become more sensitive to upland freshwater supply under future SLR, compared with previous predictions without marsh evolution. This study demonstrates the importance of coastal evolution to coastal freshwater‐saltwater interaction. The eco‐geomorphologic effect may not be ignored when evaluating coastal SWI under SLR at decadal or century scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Applying Magnetic Curvature to MMS Data to Identify Thin Current Sheets Relative to Tail Reconnection

Thin current sheets (TCSs) have been postulated to be a necessary precondition for reconnection onset. Magnetic reconnection X-lines in the magnetotail have been observed to be more common duskward of midnight. We take advantage of the MMS tetrahedral formation during the 2017–2020 MMS tail seasons to calculate the thickness of the cross-tail neutral sheet relative to ion gyroradius. While a similar technique was applied to Cluster data, current sheet thickness over a broader range of radial distances has not been robustly explored before this study. We compare our analysis to recent theories regarding mechanisms of tail current sheet thinning and to recent simulations. We find MMS spent more than twice as long in ion-scale TCSs in the pre-midnight sector than post-midnight, despite nearly even plasma sheet dwell time. The dawn-dusk asymmetry in the distribution of Ion Diffusion Regions, as previously reported in relation to regions of TCSs, is also analyzed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A Prefire Approach for Probabilistic Assessments of Postfire Debris‐Flow Inundation

Increases in wildfire activity and rainfall intensification are driving more postfire debris flows (PFDF) in many regions around the world. PFDFs are most common in the first postfire year and may even occur before a fire is fully controlled. This underscores the importance of assessing postfire hazards before a fire starts. Evaluation of PFDF hazards prior to fire can help strategize interventions lessening the negative effects of future fires. However, debris-flow runout and inundation analyses are not routine in PFDF hazard assessments, partially due to time constraints and substantial uncertainties in boundary conditions. Here, we propose a prefire PFDF inundation assessment framework using a debris-flow runout model based on the Herschel-Bulkley (HB) rheology (HEC-RAS v6.1). We constrain model inputs and parameters using Bayesian posterior analysis, rainfall-runoff simulations, and a debris-flow volume model. We use observations from recent PFDF incidents in northern Arizona, USA, to calibrate model components and then apply our prefire inundation assessment framework in a nearby unburned area. Specifically, we (a) identify yield stress as the most influential factor on inundation extent and arrival time in a HB model, (b) establish posterior distributions for model parameters suitable for forward modeling by leveraging uncertainties in field observations, and (c) implement a predictive forward analysis in an area that has not burned recently to evaluate PFDF inundation under several future fire scenarios. This study improves our ability to assess postfire debris-flow hazards before a fire begins and provides guidance for future applications of single-phase rheological models when assessing PFDF hazards.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Gyrokinetic particle simulations of interactions between energetic particles and magnetic islands induced by neoclassical tearing modes

Interactions between energetic particles (EPs) and neoclassical tearing mode (NTM) islands in the DIII-D tokamak are studied using the global gyrokinetic toroidal code (GTC). GTC simulations find that the EP radial profile is partially flattened within the magnetic island regions and that there are stochastic regions in the particle phase space. Radial particle flux is induced mainly around the magnetic island regions and decreases with time to almost zero when the initial EP distribution achieves a new steady-state in the absence of EP sources. Stochastic regions of magnetic field lines induced by the superposition of multiple islands have weak effects on the particle flux when the width of stochastic regions is smaller than the EP drift orbit width. In conclusion, the perturbed parallel EP current induced by the magnetic islands has weak stabilizing effects on the linear growth rate of the NTM instability in this DIII-D experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Whistler waves generated by nongyrotropic and gyrotropic electron beams during asymmetric guide field reconnection

Using a two-dimensional particle-in-cell simulation of asymmetric reconnection with a guide field whose strength is 0.3 times the reconnecting magnetic field, we study electron distribution functions and wave intensities in the diffusion region, focusing on the electron diffusion region (EDR). Wave activities with frequencies below the electron cyclotron frequency are observed, and these are whistler waves propagating almost anti-parallel to the magnetic field. The waves are concentrated near the magnetospheric separatrix away from the X line, but the wave activity also spreads through the EDR near the X line. The reconnection outflows are asymmetric in the outflow direction in the magnetospheric side, and the wave intensity is stronger in the side of the faster electron outflow. We study the whistler waves using the fast Fourier transform, analyses of electron velocity distribution functions, and the dispersion solver calculation. Along the magnetospheric separatrix in the stronger outflow side, highly anisotropic electron beams exist with super-Alfvénic drift speeds. The dispersion analysis shows that there are two modes: a temperature anisotropy mode and a beam mode. Outside the EDR, the whistler wave intensity is highest near the separatrix, but the wave intensity decreases if we move away from the separatrix toward the magnetic neutral line because of the increase in the electron population near zero parallel velocity. In the EDR, in the velocity plane perpendicular to the magnetic field, ring/crescent electron distribution functions are observed. Near the X-line, the wave power is enhanced where nongyrotropic electrons contribute to increase the perpendicular temperature anisotropy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Kinetics of Carbon Condensation in Detonation of High Explosives: First-Order Phase Transition Theory Perspective

We report the kinetics of carbon condensation, or carbon clustering, in detonation of carbon-rich high explosives is modeled by solving a system of rate equations for concentrations of carbon particles. Unlike previous efforts, the rate equations account not only for the aggregation of particles but also for their fragmentation in a thermodynamically consistent manner. Numerical simulations are performed, yielding the distribution of particle concentrations as a function of time. In addition to that, analytical expressions are obtained for all the distinct steps and regimes of the condensation kinetics, which facilitates the analysis of the numerical results and allows one to study the sensitivity of the kinetic behavior to the variation of system parameters. The latter is important because the numerical values of many parameters are not reliably known at present. The theory of the kinetics of first-order phase transitions is found adequate to describe the general kinetic trends of carbon condensation, as described by the rate equations. Such physical phenomena and processes as the coagulation, nucleation, growth, and Ostwald ripening are observed, and their dependence on various system parameters is studied and reported. It is believed that the present work will become useful when analyzing the present and future results for the kinetics of carbon condensation, obtained from experiments or atomistic simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗