SEARCH · Search NASA
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.
We study various classes of nonlinear equations containing operators of potential type (Riesz potential). By the method of monotone operators in the Lebesgue spaces of real-valued functions L{sub p}(a, b) we prove global theorems on the existence, uniqueness, estimates, and methods of construction of their solutions. We present applications that illustrate the results obtained.
Hydropower is one of the most mature forms of renewable energy generation. The United States (US) has almost 103 GW of installed, with 80 GW of conventional generation and 23 GW of pumped hydropower [1]. Moreover, the potential for future development on Non-Powered Dams is up to 10 GW. With the US setting its goals to become carbon neutral [2], more renewable energy in the form of hydropower needs to be integrated with the grid. Currently, there are no publicly available tool that can estimate the hydropower potential for existing hydropower dams or other non-powered dams. The HydroGenerate is an open-source python library that has the capability of estimating hydropower generation based on flow rate either provided by the user or received from United States Geological Survey (USGS) water data services. The tool calculates the efficiency as a function of flow based on the turbine type either selected by the user or estimated based on the “head” provided by the user.
Here, we present a novel minimum entropy control algorithm for a class of stochastic nonlinear systems subjected to non-Gaussian noises. The entropy control can be considered as an optimization problem for the system randomness attenuation, but the mean value has to be considered separately. To overcome this disadvantage, a new representation of the system stochastic properties was given using the cumulant-generating function based on the moment-generating function, in which the mean value and the entropy was reflected by the shape of the cumulant-generating function. Based on the samples of the system output and control input, a time-variant linear model was identified, and the minimum entropy optimization was transformed to system stabilization. Then, an optimal control strategy was developed to achieve the randomness attenuation, and the boundedness of the controlled system output was analyzed. The effectiveness of the presented control algorithm was demonstrated by a numerical example. In this paper, a data-driven minimum entropy design is presented without pre-knowledge of the system model; entropy optimization is achieved by the system stabilization approach in which the stochastic distribution control and minimum entropy are unified using the same identified structure; and a potential framework is obtained since all the existing system stabilization methods can be adopted to achieve the minimum entropy objective.
Naturally Fractured Reservoirs (NFR’s) have received little attention as potential CO 2 storage sites. Two main facts deter from storage projects in fractured reservoirs: (1) CO 2 tends to be nonwetting in target formations and capillary forces will keep CO 2 in the fractures, which typically have low pore volume; and (2) the high conductivity of the fractures may lead to increased spatial spreading of the CO 2 plume. Numerical simulations are a powerful tool to understand the physics behind brine-CO 2 flow in NFR’s. Dual-porosity models are typically used to simulate multiphase flow in fractured formations. However, existing dual-porosity models are based on crude approximations of the matrix-fracture fluid transfer processes and often fail to capture the dynamics of fluid exchange accurately. Therefore, more accurate transfer functions are needed in order to evaluate the CO 2 transfer to the matrix. This work presents an assessment of CO 2 storage potential in NFR’s using dual-porosity models. We investigate the impact of a system of fractures on storage in a saline aquifer, by analyzing the time scales of brine drainage by CO 2 in the matrix blocks and the maximum CO 2 that can be stored in the rock matrix. A new model to estimate drainage time scales is developed and used in a transfer function for dual-porosity simulations. Furthermore, we then analyze how injection rates should be limited in order to avoid early spill of CO 2 (lost control of the plume) on a conceptual anticline model. Numerical simulations on the anticline show that naturally fractured reservoirs may be used to store CO 2 .
Correlations in the Lyman-α (Lyα) forest, both as a function of line of sight separation (1D) and 3D separation, provide a unique window to the distribution of matter at redshifts not accessible by current galaxy surveys. While optimal quadratic estimators have been used to measure 1D correlations, they are computationally expensive and difficult to extend to 3D analyses. On the other hand, estimators based on the Fast Fourier Transform (FFT) are significantly faster, but are affected by missing data in the spectra (masked pixels) and so far have not used pixel weights to reduce the uncertainties in the measurement. In this publication we describe how to compute the window matrix that enables forward-modelling the impact of masked pixels and weights on the FFT-based estimators. Here, we use Gaussian and hydrodynamical simulations with artificially masked pixels to validate the method on the measurement of 1D correlations. Finally, we show that the formalism can be extended to model the impact on 3D correlations, in particular on the cross-spectrum, the correlation of 1D Fourier modes as a function of transverse separation. This work will enable more precise clustering measurements with the Lyα forest dataset recently collected by the Dark Energy Spectroscopic Instrument (DESI).
Doubly-fed induction generators (DFIGs) are prone to subsynchronous resonance (SSR) with series-compensated transmission lines and experiencing SSR-generated harmonics of wide-ranging magnitudes. This paper presents a large-signal impedance-based approach for the prediction of the magnitude of SSR-generated harmonics in DFIGs. Note that the large-signal impedance of a network shows variation in its impedance response with the magnitude of perturbation at its terminals. It is discovered in this paper that the impedance of a series-compen- sated transmission line starts changing with the magnitude of SSR because of the metal-oxide varistor-based protection of series capacitors in the line. The changing line impedance halts the growth of the SSR beyond a certain point, where the DFIG forms a limit-cycle mode and emits harmonics at the SSR frequency. The discovered role of the protection of series capacitors in SSR not only enables the prediction of SSR-generated harmonics but also can be used for the timely detection of an SSR event to avoid system disruption and sudden loss of bulk generation. The large- signal impedance-based prediction of SSR-generated harmonics in DFIG-based wind and hydro generators is demonstrated using numerical simulations of a 2.5-MVA, 690-V DFIG supplying to a series-compensated transmission line.
Here, we present ComCTQMC, a GPU accelerated quantum impurity solver. It uses the continuous-time quantum Monte Carlo (CTQMC) algorithm wherein the partition function is expanded in terms of the hybridisation function (CT-HYB). ComCTQMC supports both partition and worm-space measurements, and it uses improved estimators and the reduced density matrix to improve observable measurements whenever possible. ComCTQMC efficiently measures all one and two-particle Green's functions, all static observables which commute with the local Hamiltonian, and the occupation of each impurity orbital. ComCTQMC can solve complex-valued impurities with crystal fields that are hybridized to both fermionic and bosonic baths. Most importantly, ComCTQMC utilizes graphical processing units (GPUs), if available, to dramatically accelerate the CTQMC algorithm when the Hilbert space is sufficiently large. We demonstrate acceleration by a factor of over 600 (100) in a simulation of δ-Pu at 600 K with (without) crystal fields. In easier problems, the GPU offers less impressive acceleration or even decelerates the CTQMC. Here we describe the theory, algorithms, and structure used by ComCTQMC in order to achieve this set of features and level of acceleration.
The use of pulsars as astrophysical clocks for gravitational wave (GW) experiments demands the highest possible timing precision. Pulse times of arrival (TOAs) are limited by stochastic processes that occur in the pulsar itself, along the line of sight through the interstellar medium, and in the measurement process. On timescales of seconds to hours, the TOA variance exceeds that from template-fitting errors due to additive noise. We assess contributions to the total variance from two additional effects: amplitude and phase jitter intrinsic to single pulses and changes in the interstellar impulse response from scattering. The three effects have different dependencies on time, frequency, and pulse signal-to-noise ratio. We use data on 37 pulsars from the North American Nanohertz Observatory for GWs to assess the individual contributions to the overall intraday noise budget for each pulsar. We detect jitter in 22 pulsars and estimate the average value of rms jitter in our pulsars to be ∼1% of pulse phase. We examine how jitter evolves as a function of frequency and find evidence for evolution. Finally, we compare our measurements with previous noise parameter estimates and discuss methods to improve GW detection pipelines.
Buildings rarely perform as designed/simulated and there are numerous tangible benefits if this gap is reconciled. A new scientifically rigorous yet pragmatic methodology for calibrating building energy simulations - called Enhanced Parameter Estimation (EPE) - is proposed that allows physically relevant parameter estimation rather than a blind force-fit to energy use data. Starting with a rapidly created simulation model, calibration is performed in two stages: (a) building shell calibration with the HVAC system replaced by an ideal system that meets the loads (b) HVAC system calibration with the building shell and all internal loads replaced by a box with only process loads. In the first stage, EPE identifies a small number of high-level heat flows in the energy balance, calculates them with specifically tailored individual driving functions, introduces physically significant parameters to best accomplish energy balance, and, estimates the parameters and their uncertainty bounds. Calibration is thus done with corrective heat flows without any arbitrary tuning of input parameters. Calibration accuracy is enhanced by machine learning of the residual errors. The EPE methodology is demonstrated by means of: a synthetic building and an actual 75,000 sq. ft. building in Pennsylvania. Future work needed for widespread application is discussed.
ABSTRACT We present the MARD-Y3 catalogue of between 1086 and 2171 galaxy clusters (52 per cent and 65 per cent new) produced using multicomponent matched filter (MCMF) follow-up in 5000 deg2 of DES-Y3 optical data of the ∼20 000 overlapping ROSAT All-Sky Survey source catalogue (2RXS) X-ray sources. Optical counterparts are identified as peaks in galaxy richness as a function of redshift along the line of sight towards each 2RXS source within a search region informed by an X-ray prior. All peaks are assigned a probability fcont of being a random superposition. The clusters lie at 0.02 < z < 1.1 with more than 100 clusters at z > 0.5. Residual contamination is 2.6 per cent and 9.6 per cent for the cuts adopted here. For each cluster we present the optical centre, redshift, rest frame X-ray luminosity, M500 mass, coincidence with NWAY infrared sources, and estimators of dynamical state. About 2 per cent of MARD-Y3 clusters have multiple possible counterparts, the photo-z’s are high quality with σΔz/(1 + z) = 0.0046, and ∼1 per cent of clusters exhibit evidence of X-ray luminosity boosting from emission by cluster active galactic nuclei. Comparison with other catalogues (MCXC, RM, SPT-SZ, Planck) is performed to test consistency of richness, luminosity, and mass estimates. We measure the MARD-Y3 X-ray luminosity function and compare it to the expectation from a fiducial cosmology and externally calibrated luminosity- and richness–mass relations. Agreement is good, providing evidence that MARD-Y3 has low contamination and can be understood as a simple two step selection – X-ray and then optical – of an underlying cluster population described by the halo mass function.
The largest source of uncertainty in any source inversion is the velocity model used to construct the transfer function employed in the forward model that relates observed ground motion to the seismic moment tensor. However, standard inverse procedures often does not quantify uncertainty in the seismic moment tensor due to error in the Green’s functions from uncertain event location and Earth structure. We attempt to incorporate this uncertainty into an estimation of the seismic moment tensor using a distribution of velocity models calculated in a prior effort based on different and complementary data sets. The posterior distribution of velocity models is then used to construct Green’s functions for use in Bayesian inference of an unknown seismic moment tensor using regional waveform data. The combined likelihood is estimated using data-specific error models and the posterior of the seismic moment tensor is estimated and can be interpreted in terms of most-probable source-type.
Best approximation and approximation by families of linear polynomial operators (FLPO) in the spaces L{sub p}, 0<p⩽+∞, are investigated for periodic functions of an arbitrary number of variables in terms of the generalized modulus of smoothness generated by a periodic generator which, near the origin, is assumed to be close in a certain sense to some homogeneous function of positive order. Direct and inverse theorems (Jackson- and Bernstein-type estimates) are proved; conditions on the generators are obtained under which the approximation error by an FLPO is equivalent to an appropriate modulus of smoothness. These problems are solved by going over from the modulus to an equivalent K-functional. The general results obtained are applied to classical objects in the theory of approximation and smoothness. In particular, they are applied to the methods of approximation generated by Fejér, Riesz and Bochner-Riesz kernels, and also to the moduli of smoothness and K-functionals corresponding to the conventional, Weyl and Riesz derivatives and to the Laplace operator. Bibliography: 24 titles. (paper)
The value of time (VOT) is a fundamental component used in transportation modeling, policy analysis, and economic appraisal. Decades of research and practice have empirically estimated the VOT across many factors (e.g., mode, purpose, time, comfort, etc.), yet little is known about its underlying form. Although it is well established that VOT can vary, it is still unclear whether patterns exist in this variation. The objective of this paper is not to merely estimate the VOT, but to model the VOT across multiple continuous and interacting variables. The purpose is to reveal its functional form with respect to mode, age, gender, purpose, income, and time of day to provide a generalizable understanding for future research and practice. Such an understanding can help develop simpler models and reduce the need for bespoke estimations for every conceivable variable perturbation. This research utilized a household travel survey containing 14,159 reported trips with imputed travel time and costs for the alternative mode choices. The average overall estimated VOT is 40.32 $/h, with results showing VOT varying log-linearly with income and trip distance, but following a Gaussian function (normal curve) with age. Overall, the results show that travel distance dominates VOT variation, which increases exponentially at a rate that is 3.61 times higher per mile of distance than per $10,000 of income, and that VOT by age peaks at age 54. This basic understanding of how the VOT varies sets the foundation for answering the subsequent question for why it might vary.
Purpose: The purpose of this work is to estimate effective and lung doses from a low-dose lung cancer screening CT protocol using Tube Current Modulation (TCM) across patient models of different sizes. Methods: Monte Carlo simulation methods were used to estimate effective and lung doses from a low-dose lung cancer screening protocol for a 64-slice CT (Sensation 64, Siemens Healthcare) that used TCM. Scanning parameters were from the AAPM protocols. Ten GSF voxelized patient models were used and had all radiosensitive organs identified to facilitate estimating both organ and effective doses. Predicted TCM schemes for each patient model were generated using a validated method wherein tissue attenuation characteristics and scanner limitations were used to determine the TCM output as a function of table position and source angle. The water equivalent diameter (WED) was determined by estimating the attenuation at the center of the scan volume for each patient model. Monte Carlo simulations were performed using the unique TCM scheme for each patient model. Lung doses were tallied and effective doses were estimated using ICRP 103 tissue weighting factors. Effective and lung dose values were normalized by scanspecific 32 cm CTDIvol values based upon the average tube current across the entire simulated scan. Absolute and normalized doses were reported as a function of WED for each patient. Results: For all ten patients modeled, the effective dose using TCM protocols was below 1.5 mSv. Smaller sized patient models experienced lower absolute doses compared to larger sized patients. Normalized effective and lung doses showed some dependence on patient size (R2 = 0.77 and 0.78, respectively). Conclusion: Effective doses for a low-dose lung screening protocol using TCM were below 1.5 mSv for all patient models used in this study. Institutional research agreement, Siemens Healthcare; Past recipient, research grant support, Siemens Healthcare; Consultant, Toshiba America Medical Systems; Consultant, Samsung Electronics.
Abstract We introduce the Weak-form Estimation of Nonlinear Dynamics (WENDy) method for estimating model parameters for non-linear systems of ODEs. Without relying on any numerical differential equation solvers, WENDy computes accurate estimates and is robust to large (biologically relevant) levels of measurement noise. For low dimensional systems with modest amounts of data, WENDy is competitive with conventional forward solver-based nonlinear least squares methods in terms of speed and accuracy. For both higher dimensional systems and stiff systems, WENDy is typically both faster (often by orders of magnitude) and more accurate than forward solver-based approaches. The core mathematical idea involves an efficient conversion of the strong form representation of a model to its weak form, and then solving a regression problem to perform parameter inference. The core statistical idea rests on the Errors-In-Variables framework, which necessitates the use of the iteratively reweighted least squares algorithm. Further improvements are obtained by using orthonormal test functions, created from a set of $$C^{\infty }$$ C ∞ bump functions of varying support sizes.We demonstrate the high robustness and computational efficiency by applying WENDy to estimate parameters in some common models from population biology, neuroscience, and biochemistry, including logistic growth, Lotka-Volterra, FitzHugh-Nagumo, Hindmarsh-Rose, and a Protein Transduction Benchmark model. Software and code for reproducing the examples is available at https://github.com/MathBioCU/WENDy .
In a typical finite temperature quantum Monte Carlo (QMC) simulation, estimators for simple static observables such as specific heat and magnetization are known. With a great deal of system-specific manual labor, one can sometimes also derive more complicated non-local or even dynamic observable estimators. In contrast, we show that arbitrary static observables can be estimated within the permutation matrix representation (PMR) flavor for any Hamiltonian. We then generalize these results to general imaginary-time correlation functions and non-trivial integrated susceptibilities thereof. Finally, we demonstrate the practical versatility of our method by estimating various non-local, random observables for the transverse-field Ising model on a square lattice and a toy random model.
Exploring biodiversity linkages between above-ground and below-ground biota is a core topic in ecology, and can have implications on our understanding of ecosystem process stability. Yet, this topic still remains underexplored. Here, we explored diversity linkages, in terms of both alpha- and beta- diversity, between plant and top soil microbial communities in a semi-arid grassland ecosystem. Soil microbial community structure was assessed based on both 16S rRNA and functional genes, and plant community composition was evaluated by traditional “species composition” and a newly-defined “biomass composition”, which includes the information on the biomass of each species. The bacterial alpha-diversity, expressed as the richness and Shannon diversity of 16S rRNA genes, was significantly correlated with plant species richness and Shannon diversity, whereas the alpha-diversity of microbial functional genes showed marginal association with total plant biomass. Microbial beta-diversity, evaluated by 16S rRNA genes, showed close relationship with plant beta-diversity estimated by both “species composition” and “biomass composition”, while the microbial beta-diversity based on functional genes was only associated with the compositional variation in above-ground plant biomass. These results showed that the differences in metabolic potential of soil microbial communities, which is closely related with ecosystem functions, can be better predicted by the variation of plant-derived resources returned to soil, than merely by the species composition of the macro-organism communities.
The β-decay strength function of nuclides produced in fission is important as it dictates the distribution of decay energy between electrons, neutrinos, and γ rays and so is critical for calculating decay heat in reactors and for estimating the reactor antineutrino spectrum. Several experimental techniques are available to determine this strength function, including electron spectroscopy, γ-ray calorimetry (TAGS spectroscopy), and detailed, high-resolution spectroscopy with modern large high-purity germanium arrays. This work investigates the decay of the well-known and strongly produced fission fragment 141 Ba. A beam of 141 Cs was implanted at the target position of the Gammasphere and the subsequent decay of the daughter 141 Ba was studied. Extensive decay spectroscopy was possible up to the decay Q value of 3.197(7) MeV, including a significant extension of the level scheme and detailed angular correlation measurements for all levels with greater than 0.25% β feeding. The distribution of the β-decay strength was then inferred and compared to previous calorimetric studies. In conclusion, the agreement was excellent and provides a benchmark for comparing strength function methods and data for a more detailed understanding of the structure of 141 La.