Search NASA⌕ Search

SEARCH · Search NASA

Results for “calculation method”

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 739 records · Page 41

A Mixed Finite Volume Element Method for Flow Calculations in Porous Media

A key ingredient in the simulation of flow in porous media is the accurate determination of the velocities that drive the flow. The large scale irregularities of the geology, such as faults, fractures, and layers suggest the use of irregular grids in the simulation. Work has been done in applying the finite volume element (FVE) methodology as developed by McCormick in conjunction with mixed methods which were developed by Raviart and Thomas. The resulting mixed finite volume element discretization scheme has the potential to generate more accurate solutions than standard approaches. The focus of this paper is on a multilevel algorithm for solving the discrete mixed FVE equations. The algorithm uses a standard cell centered finite difference scheme as the 'coarse' level and the more accurate mixed FVE scheme as the 'fine' level. The algorithm appears to have potential as a fast solver for large size simulations of flow in porous media.

Jim E Jones↗

The photoionization cross section of magnesium near threshold.

Description of a method for calculating the photoionization cross section of magnesium, using the technique of configuration reaction in the continuum. The two series of levels which give rise to autoionization near the threshold are explicitly included in the calculation. The described method is a logical extension of the basic technique developed by Atlick and Moore (1966) for similar calculations in application to helium.

Bates, G. N.↗

Calculation of the Vertical Strata Load of Utility Tunnel Crossing Ground Fissure Zone

A ground fissure is a geological disaster in which the vertical dislocation of strata causes surface rupture. Ground fissures can cause extreme harm to the surface and underground buildings. Ground fissure activity can result in different settlement on the two sides of the strata, which will generate additional stress (pressure) that differs from the stress of the general stratum on underground structures across the ground fissure zone. It is essential to assess the effective stress of strata in the design of underground engineering structures across a ground fissure zone. The Xi’an ground fissure through a utility tunnel was focus of the research, and a physical model and data for oblique crossing of the 45° ground fissure were analyzed. A model of the utility tunnel structure was established, including the surrounding soil load as an active ground fissure environment. This model was used to calculate the vertical formation pressure of the overlying soil on the utility tunnel. A method to calculate the overlying load on the utility tunnel caused by ground fissure activity was proposed and compared with the calculation based on the A. Marston principle. The results showed that the ground fissure load calculation method based on the strata-holding effect can effectively calculate the earth pressure of the surrounding soil layer of the utility tunnel in the cross-ground fissure section. The results of this work provide guidance and reference value for the design of a utility tunnel in an area with the potential for a ground fissure.

Zhang, Dan↗

Pressure algorithm for elliptic flow calculations with the PDF method

An algorithm to determine the mean pressure field for elliptic flow calculations with the probability density function (PDF) method is developed and applied. The PDF method is a most promising approach for the computation of turbulent reacting flows. Previous computations of elliptic flows with the method were in conjunction with conventional finite volume based calculations that provided the mean pressure field. The algorithm developed and described here permits the mean pressure field to be determined within the PDF calculations. The PDF method incorporating the pressure algorithm is applied to the flow past a backward-facing step. The results are in good agreement with data for the reattachment length, mean velocities, and turbulence quantities including triple correlations.

Anand, M. S.↗

A computationally-efficient method for flamelet calculations

A new open-source code for the simulation of the diffusion flamelet equations is proposed. Emphasis is placed on using an approximate Jacobian to reduce the computational cost of the matrix operations. Performance of the proposed solvers is tested by performing flamelet calculations with kinetic mechanisms of varying sizes. For the unity Lewis number equations, the present iterative Newton solver using an approximate Jacobian greatly outperforms direct Newton solvers using exact Jacobians. The computation cost scales linearly with the number of species, leading to a reduction in solution times by two orders of magnitude for mechanisms containing thousands of species. The applicability of the Jacobian approximations to the solution of the non-unity Lewis number flamelet equations is assessed. The approximations are generally inadequate to solve the full non-unity Lewis number equations but can be used in some applications depending on the balance of terms in the flamelet equations. As an example, the flamelet solver is applied to the study of sooting tendencies in laminar co-flow diffusion flames where modified non-unity Lewis number flamelet equations, previously shown to accurately reproduce experimentally-measured Yield Sooting Indices (YSI), are solved. Here, the accelerated flamelet solver is well suited for sensitivity analysis and uncertainty quantification with large detailed kinetic mechanisms, tasks for which the computational cost was previously prohibitive.

42 ENGINEERING↗

Benchmarking Density Functional Theory Methods for Efficient Calculations of a Strongly Correlated Li 1– x Ni 1– y O 2−δ System

Transition metal oxides (TMOs), such as LiNiO 2 , are promising candidates for energy storage and electronic devices due to their unique electronic properties, exceptional physical and chemical characteristics, and ability to adopt multiple oxidation states. However, accurately predicting their properties using mean-field density functional theory (DFT) is challenging due to the presence of strongly correlated d-electrons and the complex interplay between their structural, electronic, and magnetic responses. These challenges are further exacerbated by the need to model defects, surfaces, and interfaces, which require computationally efficient, large-scale simulations. To address these issues, we carry out a benchmark study on the Li 1–x NiO 2 system, evaluating the performance of several popular functionals. Our findings demonstrate that combining SCAN functional relaxation with single-step HSE calculations provides a practical and scalable computational strategy. This approach balances accuracy and efficiency, enabling high-throughput simulations of strongly correlated TMOs and improved predictive modeling capability of TMOs for practical applications.

25 ENERGY STORAGE↗

On the uncertainty of long-period return values of extreme daily precipitation

Methods for calculating return values of extreme precipitation and their uncertainty are compared using daily precipitation rates over the Western U.S. and Southwestern Canada from a large ensemble of climate model simulations. The roles of return-value estimation procedures and sample size in uncertainty are evaluated for various return periods. We compare two different generalized extreme value (GEV) parameter estimation techniques, namely L-moments and maximum likelihood (MLE), as well as empirical techniques. Even for very large datasets, confidence intervals calculated using GEV techniques are narrower than those calculated using empirical methods. Furthermore, the more efficient L-moments parameter estimation techniques result in narrower confidence intervals than MLE parameter estimation techniques at small sample sizes, but similar best estimates. It should be noted that we do not claim that either parameter fitting technique is better calibrated than the other to estimate long period return values. While a non-stationary MLE methodology is readily available to estimate GEV parameters, it is not for the L-moments method. Comparison of uncertainty quantification methods are found to yield significantly different estimates for small sample sizes but converge to similar results as sample size increases. Finally, practical recommendations about the length and size of climate model ensemble simulations and the choice of statistical methods to robustly estimate long period return values of extreme daily precipitation statistics and quantify their uncertainty.

54 ENVIRONMENTAL SCIENCES↗

Numerical methods for the calculation of three-dimensional nozzle exhaust flow fields

Numerical codes developed for the calculation of three-dimensional nozzle exhaust flow fields associated with hypersonic airbreathing aircraft are described. Both codes employ reference plane grid networks with respect to three coordinate systems. Program CHAR3D is a characteristic code utilizing a new wave preserving network within the reference planes, while program BIGMAC is a finite difference code utilizing conservation variables and a one-sided difference algorithm. Secondary waves are numerically captured by both codes, while the underexpansion shock and plume boundary are treated discretely. The exhaust gas properties consist of hydrogen-air combustion product mixtures in local chemical equilibrium. Nozzle contours are treated by a newly developed geometry package based on dual cubic splines. Results are presented for simple configurations demonstrating two- and three-dimensional multiple wave interactions.

Dash, S. M.↗

A simplified method for the calculation of diurnally averaged photodissociation rates

Expressions are given for the calculation of the diurnally averaged transmission of solar flux from a series summation involving relatively few terms and readily available mathematical functions. It is shown that these expressions can be scaled to remove much of the dependence on latitude and solar declination so that at large optical depth, they become asymptotic to a function which depends only on the optical depth at noon and displays an inverse-square-root variation. From this variation, a simple analytic form containing one adjustable parameter is proposed. The functional dependence of this parameter on latitude and solar declination is determined by numerical fitting procedures. The accuracy of the analytic approximation is better than 4.5% for all altitudes, latitudes and seasons of the year.

Boughner, R. E.↗