Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed state estimation”

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 721 records · Page 40

Examining the Impacts of Great Lakes Temperature Perturbations on Simulated Precipitation in the Northeastern United States

Most inland water bodies are not resolved by general circulation models, requiring that lake surface temperatures be estimated. Given the large spatial and temporal variability of the surface temperatures of the North American Great Lakes, such estimations can introduce errors when used as lower boundary conditions for dynamical downscaling. Lake surface temperatures (LSTs) influence moisture and heat fluxes, thus impacting precipitation within the immediate region and potentially in regions downwind of the lakes. For this study, the Advanced Research version of the Weather Research and Forecasting Model (WRF-ARW) was used to simulate precipitation over the six New England states during a 5-yr historical period. The model simulation was repeated with perturbed LSTs, ranging from 10°C below to 10°C above baseline values obtained from reanalysis data, to determine whether the inclusion of erroneous LST values has an impact on simulated precipitation and synoptic-scale features. Results show that simulated precipitation in New England is statistically correlated with LST perturbations, but this region falls on a wet–dry line of a larger bimodal distribution. Wetter conditions occur to the north and drier conditions occur to the south with increasing LSTs, particularly during the warm season. Additionally, the precipitation differences coincide with large-scale anomalous temperature, pressure, and moisture patterns. Care must therefore be taken to ensure reasonably accurate Great Lakes surface temperatures when simulating precipitation, especially in southeastern Canada, Maine, and the mid-Atlantic region.

54 ENVIRONMENTAL SCIENCES↗

Calculating Methane Emissions from Offshore Facilities Using Bottom-Up Methods

With changing demands in regulation, understanding methane emissions from offshore oil and gas production infrastructure has become increasingly important. Reported emissions from facilities in the Gulf of Mexico range from zero to thousands of tons of methane per hour, but these is currently no clear understanding of how this range compares to expected emissions from normally operating facilities. To generate realistic emission estimates, we create two bottom-up models that simulate emissions from facilities operating in the Gulf of Mexico. We estimate type 1 prototypical facilities (typically unmanned, older, lower-producing platforms in shallow water with little processing equipment, compressors, or storage tanks) to emit an average of 13 kg CH 4 h −1 , which corresponds to a loss of 2.7% of the average facility production. Type 2 prototypical facilities (continuously manned, higher production and operate in deeper water with processing equipment, oil storage tanks, compressors and power generation) emit an average of 88 kg CH 4 h −1 , which corresponds to a loss of 2.5% of production. The average measured emission from type 1 facilities was 18 kg CH 4 h −1 with a median production loss estimated at 8%. The average measured emission from type 2 facilities was 36 kg CH 4 h −1 with a median production loss estimated at 2.4%. Using emission factors that consider the long-tail emission distribution partly reconciles the difference between modelled and measured emission estimates, but we suggest the current the fugitive emission estimate may be an underestimate and more data on the number and size of fugitive emissions could explain differences between the modelled and measured emission estimate. We suggest the bottom-up approach described here that uses production data coupled with facility equipment could be used to identify facilities that have abnormally large measured emissions, caused by methodological failure or larger than expected fugitive emissions, which should be targeted for further evaluation resulting in remeasurement or identification of source type so that a more accurate estimates can be made on the absolute emission.

bottom-up↗

The Economic Accessibility of CO2 Sequestration through Bioenergy with Carbon Capture and Storage (BECCS) in the US

Bioenergy with carbon capture and storage (BECCS) is one strategy to remove CO2 from the atmosphere. To assess the potential scale and cost of CO2 sequestration from BECCS in the US, this analysis models carbon sequestration net of supply chain emissions and costs of biomass production, delivery, power generation, and CO2 capture and sequestration in saline formations. The analysis includes two biomass supply scenarios (near-term and long-term), two biomass logistics scenarios (conventional and pelletized), and two generation technologies (pulverized combustion and integrated gasification combined cycle). Results show marginal cost per tonne CO2 (accounting for costs of electricity and CO2 emissions of reference power generation scenarios) as a function of CO2 sequestered (simulating capture of up to 90% of total CO2 sequestration potential) and associated spatial distribution of resources and generation locations for the array of scenario options. Under a near-term scenario using up to 206 million tonnes per year of biomass, up to 181 million tonnes CO2 can be sequestered annually at scenario-average costs ranging from $62 to $137 per tonne CO2; under a long-term scenario using up to 740 million tonnes per year of biomass, up to 737 million tonnes CO2 can be sequestered annually at scenario-average costs ranging from $42 to $92 per tonne CO2. These estimates of CO2 sequestration potential may be reduced if future competing demand reduces resource availability or may be increased if displaced emissions from conventional power sources are included. Results suggest there are large-scale opportunities to implement BECCS at moderate cost in the US, particularly in the Midwest, Plains States, and Texas.

09 BIOMASS FUELS↗

Inverse rendering of fusion plasmas: Inferring plasma composition from imaging systems

In this work, we develop a differentiable rendering pipeline for visualising plasma emission within tokamaks, and estimating the gradients of the emission and estimating other physical quantities. Unlike prior work, we are able to leverage arbitrary representations of plasma quantities and easily incorporate them into a non-linear optimisation framework. The efficiency of our method enables not only estimation of a physically plausible image of plasma, but also recovery of the neutral Deuterium distribution from imaging and midplane measurements alone. We demonstrate our method with three different levels of complexity showing first that a poloidal neutrals density distribution can be recovered from imaging alone, second that the distributions of neutral Deuterium, electron density and electron temperature can be recovered jointly, and finally, that this can be done in the presence of realistic imaging systems that incorporate sensor cropping and quantisation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Matrix Approach to Land Carbon Cycle Modeling

Land ecosystems contribute to climate change mitigation by taking up approximately 30% of anthropogenically emitted carbon. However, estimates of the amount and distribution of carbon uptake across the world's ecosystems or biomes display great uncertainty. The latter hinders a full understanding of the mechanisms and drivers of land carbon uptake, and predictions of the future fate of the land carbon sink. The latter is needed as evidence to inform climate mitigation strategies such as afforestation schemes. To advance land carbon cycle modeling, we have developed a matrix approach. Land carbon cycle models use carbon balance equations to represent carbon exchanges among pools. Our approach organizes this set of equations into a single matrix equation without altering any processes of the original model. The matrix equation enables the development of a theoretical framework for understanding the general, transient behavior of the land carbon cycle. While carbon input and residence time are used to quantify carbon storage capacity at steady state, a third quantity, carbon storage potential, integrates fluxes with time to define dynamic disequilibrium of the carbon cycle under global change. The matrix approach can help address critical contemporary issues in modeling, including pinpointing sources of model uncertainty and accelerating spin-up of land carbon cycle models by tens of times. The accelerated spin-up liberates models from the computational burden that hinders comprehensive parameter sensitivity analysis and assimilation of observational data to improve model accuracy. Such computational efficiency offered by the matrix approach enables substantial improvement of model predictions using ever-increasing data availability. Overall, the matrix approach offers a step change forward for understanding and modeling the land carbon cycle.

54 ENVIRONMENTAL SCIENCES↗

Astrometric microlensing probes of the isolated neutron star population with Roman

Notoriously hard to detect and study, isolated neutron stars (NSs) might provide valuable answers to fundamental questions about stellar evolution and explosion physics. With the upcoming Roman Space Telescope, scheduled for launch in 2026, a new and powerful channel for their detection will become available: astrometric microlensing. We set out to create a realistic sample of simulated gravitational microlensing events as observed by Roman with the Galactic Bulge Time Domain Survey. We focus in particular on the population of NS lenses, which has until now been largely understudied. We used dedicated Galactic models tailored for application to microlensing by compact objects. In addition to populations of stars, white dwarfs, and black holes, we simulated four different NS populations with Maxwellian natal kick distributions: v = (150, 250, 350, 450) km/s. For each simulation, we applied projected Roman precision, cadence, and detectability criteria. We found that the parameter space log10 t E –log 10 θ E , which will be accessible to Roman observations, is efficient for the classification of stellar remnants. We found a feature in this space that is characteristic of NSs; using this feature, optimal samples of NS candidates can be constructed from Roman -like datasets. We describe the dependence of the observable parameter distributions on the assumed mean kick velocities. As the effects of natal kicks are very complex and mutually counteracting, we suggest that more detailed studies focused on the dynamics of NSs are needed in anticipation of Roman and future surveys. We estimate that Roman will observe approximately 11 000 microlensing events, including ~100 with NS lenses, whose photometric and astrometric signals are detectable; the event yield decreases by 38% when gap-filling low-cadence observations are not included. We make all simulated microlensing event datasets publicly available in preparation for Roman data.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigation of Anchor Link Magnetohydrodynamic Effects in the Toroidally Symmetric Lead-Lithium (TSLL) Blanket Concept

This paper presents a comprehensive study on the magnetohydrodynamic (MHD) flow in a slotted channel with a cylindrical anchor link, which is an imperative component of the toroidally symmetric lead-lithium (TSLL) blanket concept. Following the validation of the MHD solver implemented in COMSOL and a proposed subtraction approach to compute the anchor link pressure drop, the effects of computational domain size on the pressure drop and velocity distribution are examined. The results show that the pressure drop associated with the anchor link and the maximum flow velocity follow an asymptotic trend as the domain width increases, with wall-induced pressure drop being more dominant than that of the anchor link. The velocity distribution analysis revealed the formation of an internal boundary layer extending along the magnetic field direction, which is a unique feature of the investigated MHD flow. An estimation of the total MHD pressure drop associated with the array of anchor links under the TSLL blanket conditions suggests ~0.3 MPa, which is significantly lower than the recommended maximum allowable blanket pressure drop of 2 MPa. In conclusion, this work offers valuable insights into the anchor link–associated MHD phenomena and serves as a foundation for further development of the TSLL blanket concept for future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The role of high-order electron correlation effects in a model system for non-valence correlation-bound anions

The diffusion Monte Carlo (DMC), auxiliary field quantum Monte Carlo (AFQMC), and equation-of-motion coupled cluster (EOM-CC) methods are used to calculate the electron binding energy (EBE) of the non-valence anion state of a model (H2O)4 cluster. Two geometries are considered, one at which the anion is unbound and the other at which it is bound in the Hartree–Fock (HF) approximation. It is demonstrated that DMC calculations can recover from the use of a HF trial wave function that has collapsed onto a discretized continuum solution, although larger EBEs are obtained when using a trial wave function for the anion that provides a more realistic description of the charge distribution and, hence, of the nodal surface. For the geometry at which the cluster has a non-valence correlation-bound anion, both the inclusion of triples in the EOM-CC method and the inclusion of supplemental diffuse d functions in the basis set are important. DMC calculations with suitable trial wave functions give EBE values in good agreement with our best estimate EOM-CC result. AFQMC using a trial wave function for the anion with a realistic electron density gives a value of the EBE nearly identical to the EOM-CC result when using the same basis set. For the geometry at which the anion is bound in the HF approximation, the inclusion of triple excitations in the EOM-CC calculations is much less important. The best estimate EOM-CC EBE value is in good agreement with the results of DMC calculations with appropriate trial wave functions.

Upadhyay, Shiv (ORCID:0000000285010501)↗

Assessing the Viability of Geothermal Microgrid Deployment: A Geospatial Analysis Across the United States

Geothermal microgrids hold a potential of supplying clean and dependable power to communities throughout the United States (US), all while sidestepping the expenses associated with connecting to strained or isolated power grids. Nonetheless, their implementation is still in its early stages in the country. The objective of this analysis is to leverage available data to pinpoint regions across the US that exhibit favorable conditions for the development of geothermal microgrids. Drawing from a variety of sources, including estimates of geothermal resources, the costs associated with geothermal energy generation and electricity transmission, existing microgrid locations, and subsidy programs, we aim to identify promising areas for further exploration. By mapping out the contiguous US, Alaska, and Hawaii, we delineate regions with high relative favorability for geothermal microgrid deployment. Our findings reveal the presence of highly favorable regions across the Western states of the contiguous US, as well as isolated areas in Alaska and Hawaii. Furthermore, we delve into a discussion on state policies and incentive programs, considering their role in fostering favorable conditions or posing barriers to geothermal microgrid development.

Alaska↗

Scale effects on core design, fuel costs, and spent fuel volume of pressurized water reactors

The desire to improve the economic competitiveness and deployment pace of nuclear energy through modularization, manufacturing, and series production had led to the development of smaller size reactors. As the standard 17x17 fuel technology is mainly maintained in the pressurized water reactors (PWRs) category, this translates into a lower number of fuel assemblies in the core and sometimes a reduced fuel height. To assess the impact of such scale change in core design on fuel cycle cost and spent fuel volume, a scoping analysis tool is developed based on infinite lattice calculations, leakage, fuel management reduced models, and levelized unit cost of electricity (LCOE) estimate. As such, cost dynamics driven by fuel specific power, burnup, core leakage, feed, cycle length, fuel assembly height as well as uranium market data are captured with consistent set of assumptions and analysis methods. A selection of 5 reactor designs representative of leading PWR developers is assessed and compared. Pursuing higher specific powers and optimal burnups are highlighted as the main fuel cost reduction drivers, nevertheless, practical limitations and opportunities must be evaluated to establish the feasibility of such enhanced fuel operation. In consequence, a detailed core design is performed using SIMULATE3 code for 5 PWR variations including natural and forced coolant circulation modes, two reactor scales, power densities of 73, 112, and 123 kW/l and higher discharge burnups. Design and optimization are performed at the lattice level, for the reflector, and at the core loading level. Satisfactory steady-state operation including power distribution, coolant operating limits, and reactivity requirements are analyzed and reported in this paper. The fuel economics of the detailed core designs confirm the scoping analysis findings. Despite the unlocked power uprates in small PWRs, the achievable burnup for a given fuel specific power requires more enrichment and shorter fuel height results in higher fabrication costs per mass of fuel, which makes scaling down core size a more expensive endeavor on the fuel cycle front. Spent fuel volumes are reported for the PWRs designed in this paper. Furthermore, these volumes are driven by the core average discharge burnup regardless of the scale in consideration. Additional cost and core performance aspects related to heavy reflector gains, fuel-reflector substitution, and disposal cost policy in the U.S. are examined.

42 ENGINEERING↗

Transverse single-spin asymmetries in $γ$ SIDIS as a direct probe of quark-gluon-quark longitudinal momentum structure

Transverse single-spin asymmetries in the semi-inclusive deep-inelastic production of isolated photons ($\gamma$SIDIS), $A_{UT}^{\gamma {\rm SIDIS}}$, provide an unprecedented opportunity to extract the quark-gluon-quark correlators $F_{FT}(x,x')$ and $G_{FT}(x,x')$ point-by-point in their full support $x,x'$. We utilize realistic models for these functions, based on input from the Sivers transverse momentum dependent parton distribution function and imposing constraints from the $d_2$ matrix element calculated in lattice QCD, in order to provide numerical estimates for $A_{UT}^{\gamma {\rm SIDIS}}$ at the Electron-Ion Collider (EIC). We thoroughly explore the EIC phase space in order to isolate in which regions the asymmetry can be sizable, finding it can be as much as $10\%$ or larger for certain kinematics. Given that $F_{FT}(x,x')$ and $G_{FT}(x,x')$ are basically unknown, $A_{UT}^{\gamma {\rm SIDIS}}$ will be an important future measurement to learn about multi-parton correlations in the nucleon.

Electron-Ion Collider↗

Quantification of Shallow Defect Suppression in Aminosilane-Passivated Lead Halide Perovskites via Optical and Electrochemical Methods

We combine fluence-dependent photoluminescence (PL) and electrochemical characterization to estimate defect-dependent trapping rates, energetic distributions, and densities of defects, in lead halide perovskite films Cs 0.1 FA 0.9 PbI 3 and Cs 0.17 FA 0.83 Pb­(I 0.75 Br 0.25 ) 3 , before and after trap passivation with vapor-deposited (3-aminopropyl)­trimethoxysilane (APTMS). For both compositions, the PL studies show that energetically shallow subgap near-valence defect densities lead to bimolecular nonradiative recombination which is reduced by 92% after APTMS treatment. Electrochemical characterization of these same active layers shows for the first time that APTMS passivation impacts mainly on near-valence (mobile anionic) defects that decrease in density by 4 orders of magnitude, from 10 18 cm –3 to 10 14 cm –3 , as a result of silane modification. Here, the combination of PL and electrochemical characterization promises a unique approach to speciation and quantification of a broader distribution of trap states in perovskites and provides straightforward assessments of the efficacy of defect mitigation strategies.

Defects↗

High-temperature 2D ferromagnetism in conjugated microporous porphyrin-type polymers

The need for magnetic 2D materials that are stable to the enviroment and have high Curie temperatures is very important for various electronic and spintronic applications. We have found that two-dimensional porphyrin-type aza-conjugated microporous polymer crystals are such a material (Fe-aza-CMPs). Fe-aza-CMPs are stable to CO, CO 2 , and O 2 atmospheres and show unusual adsorption, electronic, and magnetic properties. Indeed, they are semiconductors with small energy band gaps ranging from 0.27 eV to 0.626 eV. CO, CO 2 , and O 2 molecules can be attached in three different ways where single, double, or triple molecules are bound to iron atoms in Fe-aza-CMPs. For different attachment configurations we find that for CO and CO 2 a uniform distribution of the molecules is most energetically favorable while for O 2 molecules aggregation is most energetically preferable. The magnetic moments decrease from 4 to 2 to 0 for singly, doubly, triply occupied configurations for all gasses respectively. The most interesting magnetic properties are found for O 2 molecules attached to the Fe-aza-CMP. For a single attachment configuration we find that an antiferromagnetic state is favorable. When two O 2 molecules are attached, the calculations show the highest exchange integral with a value of J = 1071 μeV. This value has been verified by two independent methods where in the first method J is calculated by the energy difference between ferromagnetic and anitferromagnetic configurations. The second method is based on the frozen magnon approach where the magnon dispersion curve has been fitted by the Ising model. Finally, for the second method J has been estimated at J = 1100 μeV in excellent agreement with the first method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rising infrastructure inequalities accompany urbanization and economic development

Impending global urban population growth is expected to occur with considerable infrastructure expansion. However, our understanding of attendant infrastructure inequalities is limited, highlighting a critical knowledge gap in the sustainable development implications of urbanization. Using satellite data from 2000 to 2019, we examine country-level population-adjusted biases in infrastructure distribution within and between regions of varying urbanization levels and derive four key findings. First, we find long-run positive associations between infrastructure inequalities and both urbanization and economic development. Second, our estimates highlight increasing infrastructure inequalities across most of the countries examined. Third, we find greater future infrastructure inequality increases in the global south, where inequalities will rise more in countries with substantial urban primacy. Fourth, we find that infrastructure inequality may evolve differently than economic inequalities. Overall, advancing sustainable development vis-à-vis urbanization and economic development will require intentional infrastructure planning for spatial equity.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Single-shot in-line x-ray phase-contrast imaging of void-shockwave interactions in fusion energy materials

Recent breakthroughs in nuclear fusion, specifically the report of reactions exceeding scientific breakeven at the National Ignition Facility (NIF), highlight the potential of inertial fusion energy (IFE) as a sustainable and virtually limitless energy source. However, further progress in IFE requires characterization of defects in ablator materials and how they affect fuel capsule compression. Voids within the ablator can degrade energy yield, but their impact on the density distribution has primarily been studied through simulations, with limited high-resolution experimental validation. To address this, we used the x-ray free-electron laser (XFEL) at the matter in extreme conditions (MECs) instrument at the Linac coherent light source (LCLS) to capture 2D x-ray phase-contrast (XPC) images of a void-bearing sample with a composition similar to inertial confinement fusion (ICF) ablators. By driving a compressive shockwave through the sample using MEC's long-pulse laser system, we analyzed how voids influence shockwave propagation and density distribution during compression. To quantify this impact, we extracted phase information using two phase retrieval algorithms. First, we applied the contrast transfer function (CTF) method, paired with Tikhonov regularization and a fast optimization approach to generate an initial phase estimate. We then refined the result using a projected gradient descent (PGD) method that works directly with the sample's refractive index. Comparing these results with radiation adaptive grid Eulerian (xRAGE) radiation hydrodynamic simulations enables identification of model validation needs or improvements. By calculating phase maps in situ, it becomes possible to reconstruct areal density maps, improving understanding of laser-capsule interactions and advancing IFE research.

Hodge, D. S. [Colorado State Univ., Fort Collins, ↗

A Statistical Interpolation Code for Ocean Analysis and Forecasting

Abstract We present a data assimilation package for use with ocean circulation models in analysis, forecasting, and system evaluation applications. The basic functionality of the package is centered on a multivariate linear statistical estimation for a given predicted/background ocean state, observations, and error statistics. Novel features of the package include support for multiple covariance models, and the solution of the least squares normal equations either using the covariance matrix or its inverse—the information matrix. The main focus of this paper, however, is on the solution of the analysis equations using the information matrix, which offers several advantages for solving large problems efficiently. Details of the parameterization of the inverse covariance using Markov random fields are provided and its relationship to finite-difference discretizations of diffusion equations are pointed out. The package can assimilate a variety of observation types from both remote sensing and in situ platforms. The performance of the data assimilation methodology implemented in the package is demonstrated with a yearlong global ocean hindcast with a 1/4° ocean model. The code is implemented in modern Fortran, supports distributed memory, shared memory, multicore architectures, and uses climate and forecasts compliant Network Common Data Form for input/output. The package is freely available with an open source license from www.tendral.com/tsis/ .

Srinivasan, Ashwanth↗

Studies of laser stimulated photodetachment from nanoparticles for particle charge measurements

Determining nanoparticle charge is more challenging than that for microparticles due to change in the particle size during the synthesis substantial plasma property variations, and difficulties in visualizing individual particles, rendering conventional microparticle charge diagnostics ineffective in dusty plasma. In this work, we utilized laser-stimulated photodetachment (LSPD) to deduce the mean charge of nanoparticles. Nanoparticles were grown in an Ar/C 2 H 2 mixture using a capacitively coupled RF discharge and the LSPD induced changes in the electron current monitored by a cylindrical Langmuir probe. LSPD signals were obtained and analyzed across different dust growth phases. The prolonged decay of electron current pulses was attributed to the presence of residual negative ions, caused by the effective electrostatic trapping of these ions and the potential post—LSPD re-formation of new ones. The particle charge was estimated by combining the laser-stimulated photodetachment signal from the probe with the dust density obtained from laser-light extinction using the measured nanoparticle size distribution. For a nanoparticles size range of approximately 100–250 nm and mean diameter of $d$ p ∼154.37 nm, the effective mean charge was estimated to be $\langle$$Q$$\rangle$ d $≈$ 37 elementary charge units. The measured charge values are lower than those predicted by orbital motion limited theory, which may be attributed to significant electron depletion in the nanodusty plasma. LSPD results in Ar/C 2 H 2 nano-dusty plasma confirm the applicability of this method for estimating individual nanoparticle charges. However, it has also been demonstrated that electron detachment from residual background negative ions can influence the detachment current decay and must be carefully considered.

dust particle charge↗