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.
Dimensionless characteristics of turbines are discussed. On that basis, a universal method of mathematical approximation is formulated. The approximation of the turbine characteristics in a TKR-10ST turbocompressor is considered, as an example.
Angular two-point statistics of large-scale structure observables are important cosmological probes. To reach the high accuracy required by the statistical precision of future surveys, some of these statistics may need to be computed without the commonly employed Limber approximation; the exact computation however requires integration over Bessel functions, and a brute-force evaluation is slow to converge. Here, we present a new method based on our generalized FFTLog algorithm for the efficient computation of angular power spectra beyond the Limber approximation. The new method significantly simplifies the calculation and improves the numerical speed and stability. It is easily extended to handle integrals involving derivatives of Bessel functions, making it equally applicable to numerically more challenging cases such as contributions from redshift-space distortions and Doppler effects. We implement our method for galaxy clustering and galaxy-galaxy lensing power spectra. We find that using the Limber approximation for galaxy clustering in future analyses like LSST Year 1 and DES Year 6 may cause significant biases in cosmological parameters, indicating that going beyond the Limber approximation is necessary for these analyses.
Halide perovskites constitute a chemically diverse class of crystals with great promise as photovoltaic absorber materials, featuring band gaps between about 1 and 3.5 eV depending on composition. Their diversity calls for a general computational approach to predicting their band gaps. However, such an approach is still lacking. Here, we use density functional theory (DFT) and ab initio many-body perturbation theory within the GW approximation to compute the quasiparticle or fundamental band gap of a set of ten representative halide perovskites: CH 3 NH 3 PbI 3 ( MAPbI 3 ), MAPbBr 3 , CsSnBr 3 , ( MA ) 2 BiTlBr 6 , Cs 2 TlAgBr 6 , Cs 2 TlAgCl 6 , Cs 2 BiAgBr 6 , Cs 2 InAgCl 6 , Cs 2 SnBr 6 , and Cs 2 Au 2 I 6 . Comparing with recent measurements, we find that a standard generalized gradient exchange-correlation functional can significantly underestimate the experimental band gaps of these perovskites, particularly in cases with strong spin-orbit coupling (SOC) and highly dispersive band edges, to a degree that varies with composition. We show that these nonsystematic errors are inherited by one-shot G 0 W 0 and eigenvalue self-consistent GW 0 calculations, demonstrating that semilocal DFT starting points are insufficient for MAPbI 3 , MAPbBr 3 , CsSnBr 3 , ( MA ) 2 BiTlBr 6 , Cs 2 TlAgBr 6 , and Cs 2 TlAgCl 6 . On the other hand, we find that DFT with hybrid functionals leads to an improved starting point and G W 0 results in better agreement with experiment for these perovskites. Our results suggest that G W 0 with hybrid functional-based starting points are promising for predicting band gaps of systems with large SOC and dispersive bands in this technologically important class of semiconducting crystals.</
We derive a numerical approximation, namely L1–2 formula, to the Caputo–Fabrizio derivative by using a quadratic interpolation. Quadratic and cubic convergence rates are achieved for L1 and L1–2 formulas using Lagrange interpolation, respectively. We compute Caputo–Fabrizio derivatives of some known functions both theoretically and numerically. In addition, we solve non/linear sub-diffusion equations to test theoretical findings. Numerical results confirm the theoretically observed convergence rates.
It is shown that the local approximation for computing the electron distribution function depends both on the ratio between the energy relaxation length and a characteristic plasma length and on the ratio between heating and ambipolar electric fields. In particular, the local approximation is not valid at the discharge periphery even at high pressure due to the fact that the ambipolar electric field practically always is larger than the heating electric field.
In the present paper, concept of roughness for fuzzy filters with thresholds (u{sub 1},u{sub 2}) in ordered semigroups is introduced. Then, this concept is extended to fuzzy bi-filters with thresholds and fuzzy quasi-filters with thresholds. Further approximations of fuzzy ideals with thresholds, fuzzy bi-ideals with thresholds, and fuzzy interior ideals with thresholds are studied. Moreover, this concept is applied to study approximations of fuzzy quasi-ideals with thresholds and semiprime fuzzy quasi-ideals with thresholds.
Bayesian state estimation of a dynamical system from a stream of noisy measurements is important in many geophysical and engineering applications where high dimensionality of the state space, sparse observations, and model error pose key challenges. In this work, three computationally feasible, approximate Gaussian data assimilation/filtering algorithms are considered in various regimes of turbulent 2D Navier--Stokes dynamics in the presence of model error. The first source of error arises from the necessary use of reduced models for the forward dynamics of the filters, while a particular type of representation error arises from the finite resolution of observations which mix up information about resolved and unresolved dynamics. Two stochastically parameterized filtering algorithms, referred to as cSPEKF and GCF, are compared with 3DVAR---a prototypical time-sequential algorithm known to be accurate for filtering dissipative systems for a suitably inflated “background” covariance. We provide the first evidence that the stochastically parameterized algorithms, which do not rely on detailed knowledge of the underlying dynamics and do not require covariance inflation, can compete with or outperform an optimally tuned 3DVAR algorithm, and they can overcome competing sources of error in a range of dynamical scenarios.
This paper is devoted to the study of the solvability of random fuzzy fractional partial integro-differential equations under Caputo generalized Hukuhara differentiability. The notions of random fuzzy variables and fuzzy stochastic processes are developed for multivariable functions. The existence and uniqueness of two types of integral solutions generated from Darboux problem for nonlinear wave equations are proved using the successive approximations method and Gronwall’s inequality for stochastic processes. The continuous dependence on the data, the boundedness, and the stability with probability one of integral solutions are established to confirm the well posedness of our model. Some computational examples are presented to illustrate the theoretical results.
Hawking radiation is the evidence for the existence of black hole. What an observer can measure through Hawking radiation is the transmission probability. In the laboratory, miniature black holes can successfully be generated. The generated black holes are, most commonly, Myers-Perry black holes. In this paper, we will derive the rigorous bounds on the transmission probabilities for massless scalar fields of non-negative-angular-momentum modes emitted from a generated Myers-Perry black hole in six, seven, and eight dimensions. The results show that for low energy, the rigorous bounds increase with the increase in the energy of emitted particles. However, for high energy, the rigorous bounds decrease with the increase in the energy of emitted particles. When the black holes spin faster, the rigorous bounds decrease. For dimension dependence, the rigorous bounds also decrease with the increase in the number of extra dimensions. Furthermore, as comparison to the approximate transmission probability, the rigorous bound is proven to be useful.
The concept of local mass displacements (as a measure of nonlocality of the state of a physically small element of a body) is used to formulate a complete system of relations of the gradient-type mathematical model of viscous compressible fluid. In this model, the relationship between the thermomechanical processes and the process of local mass displacements, its irreversibility, and the inertia properties are taken into account. It is shown that, in the stationary approximation, the obtained system of relations enables us to justify the appearance of wedging pressure in thin interlayers of fluid.
We study uniform approximation of real-valued functions f, f(∞) = 0, on ℝ{sup +} and ℝ by real-valued simple partial fractions (the logarithmic derivatives of polynomials). We obtain a criterion for the best approximation on ℝ{sup +} and ℝ in terms of the Chebyshev alternance. This criterion is similar to the known criterion on finite segments. For the problem of approximating odd functions on ℝ we construct an alternance criterion with a weakened condition on the poles of fractions. We present a criterion for the best approximation by simple partial fractions on ℝ{sup +} and ℝ in terms of Kolmogorov. We prove analogs of the de la Vallee-Poussin alternation theorem.
Drawing on the adiabatic connection of density functional theory, exchange-correlation functionals of Kohn-Sham density functional theory are constructed which interpolate between the extreme limits of the electron-electron interaction strength. The first limit is the non-interacting one, where there is only exchange. The second limit is the strong correlated one, characterized as the minimum of the electron-electron repulsion energy. The exchange-correlation energy in the strong-correlation limit is approximated through a model for the exchange-correlation hole that is referred to as nonlocal-radius model [L. O. Wagner and P. Gori-Giorgi, Phys. Rev. A 90, 052512 (2014)]. Using the non-interacting and strong-correlated extremes, various interpolation schemes are presented that yield new approximations to the adiabatic connection and thus to the exchange-correlation energy. Some of them rely on empiricism while others do not. Several of the proposed approximations yield the exact exchange-correlation energy for one-electron systems where local and semi-local approximations often fail badly. Other proposed approximations generalize existing global hybrids by using a fraction of the exchange-correlation energy in the strong-correlation limit to replace an equal fraction of the semi-local approximation to the exchange-correlation energy in the strong-correlation limit. The performance of the proposed approximations is evaluated for molecular atomization energies, total atomic energies, and ionization potentials.
We prove an L{sup 2}-estimate for the homogenization of an elliptic operator A{sub ε} in a domain Ω with a Neumann boundary condition on the boundary ∂Ω. The coefficients of the operator A{sub ε} are rapidly oscillating over different groups of variables with periods of different orders of smallness as ε→0. We assume minimal regularity of the data, which makes it possible to impart to the result the meaning of an estimate in the operator (L{sup 2}(Ω)→L{sup 2}(Ω))-norm for the difference of the resolvents of the original and homogenized problems. We also find an approximation to the resolvent of the original problem in the operator (L{sup 2}(Ω)→H{sup 1}(Ω))-norm. Bibliography: 24 titles.
The Keller–Segel–Navier–Stokes system governs chemotaxis in liquid environments. This system is to be solved for the organism and chemoattractant densities and for the fluid velocity and pressure. It is known that if the total initial organism density mass is below 2π there exist globally defined generalised solutions, but what is less understood is whether there are blow-up solutions beyond such a threshold and its optimality. Motivated by this issue, a numerical blow-up scenario is investigated. Approximate solutions computed via a stabilised finite element method founded on a shock capturing technique are such that they satisfy a priori bounds as well as lower and L 1 (Ω) bounds for the organism and chemoattractant densities. In particular, these latter properties are essential in detecting numerical blow-up configurations, since the non-satisfaction of these two requirements might trigger numerical oscillations leading to non-realistic finite-time collapses into persistent Dirac-type measures. Our findings show that the existence threshold value 2π encountered for the organism density mass may not be optimal and hence it is conjectured that the critical threshold value 4π may be inherited from the fluid-free Keller–Segel equations. Additionally it is observed that the formation of singular points can be neglected if the fluid flow is intensified.
The primary computational challenge when simulating nonadiabatic ab initio molecular dynamics is the unfavourable compute costs of electronic structure calculations with molecular size. Simple electronic structure theories, like time-dependent density functional theory within the Tamm–Dancoff approximation (TDDFT/TDA), alleviate this cost for moderately sized molecular systems simulated on realistic time scales. Although TDDFT/TDA does have some limitations in accuracy, an appealing feature is that, in addition to including electron correlation through the use of a density functional, the cost of calculating analytic nuclear gradients and nonadiabatic coupling vectors is often computationally feasible even for moderately sized basis sets. Here in this work, some of the benefits and limitations of TDDFT/TDA are discussed and analysed with regard to its applicability as a ‘back-end’ electronic structure method for the symmetric quasi-classical Meyer–Miller model (SQC/MM). In order to investigate the benefits and limitations of TDDFT/TDA, SQC/MM is employed to predict and analyse a prototypical example of excited-state hydrogen transfer in gas-phase malonaldehyde. Then, the ring-opening dynamics of selenophene are simulated, which highlight some of the deficiencies of TDDFT/TDA. Additionally, some new algorithms are proposed that speed up the calculation of analytic nuclear gradients and nonadiabatic coupling vectors for a set of excited electronic states.
We study stochastic delay differential equations (SDDE) where the coefficients depend on the moving averages of the state process. As a first contribution, we provide sufficient conditions under which the solution of the SDDE and a linear path functional of it admit a finite-dimensional Markovian representation. As a second contribution, we show how approximate finite-dimensional Markovian representations may be constructed when these conditions are not satisfied, and provide an estimate of the error corresponding to these approximations. These results are applied to optimal control and optimal stopping problems for stochastic systems with delay.
We consider the incompressible fluid with a pressure-dependent viscosity flowing through a multiple pipe system. The viscosity–pressure relation is given by the Barus law commonly used in the engineering applications. Assuming that the ratio between pipes thickness and its length is small, we propose a rigorous asymptotic approach based on the concept of the transformed pressure. As a result, we obtain new macroscopic model describing the effective behavior of the fluid in the system. In particular, the generalized version of the Kirchhoff’s law is derived giving the explicit formula for the junction pressure. The error estimate for the asymptotic approximation is also provided. Mathematical analysis presented here can be applied to a general viscosity–pressure relation satisfied by other empiric laws.
In this paper, we interpret a fuzzy differential equation by using the strongly generalized differentiability concept. Utilizing the Generalized Characterization Theorem, we investigate the problem of finding a numerical approximation of solutions. The Runge–Kutta approximation method is implemented and its error analysis, which guarantees pointwise convergence, is given. The method applicability is illustrated by solving a linear first-order fuzzy differential equation.