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 1,495 records · Page 83

Reliable electrochemical phase diagrams of magnetic transition metals and related compounds from high-throughput ab initio calculations

Abstract Magnetic transition metals (mTM = Cr, Mn, Fe, Co, and Ni) and their complex compounds (oxides, hydroxides, and oxyhydroxides) are highly important material platforms for diverse technologies, where electrochemical phase diagrams with respect to electrode potential and solution pH can be used to effectively understand their corrosion and oxidation behaviors in relevant aqueous environments. Many previous decades-old mTM–Pourbaix diagrams are inconsistent with various direct electrochemical observations, because experimental complexities associated with extracting reliable free energies of formation (Δ f G ) lead to inaccuracies in the data used for modeling. Here, we develop a high-throughput simulation approach based on density-functional theory (DFT), which quickly screens structures and compounds using efficient DFT methods and calculates accurate Δ f G values, using high-level exchange-correlation functions to obtain ab initio Pourbaix diagrams in comprehensive and close agreement with various important electrochemical, geological, and biomagnetic observations reported over the last few decades. We also analyze the microscopic mechanisms governing the chemical trends among the Δ f G values and Pourbaix diagrams to further understand the electrochemical behaviors of mTM-based materials. Last, we provide probability profiles at variable electrode potential and solution pH to show quantitatively the likely coexistence of multiple-phase areas and diffuse phase boundaries.

Huang, Liang-Feng (ORCID:0000000329375183)↗

Analytical Identification Method of Generalized Short‐Circuit Ratio Using Phasor Measurement Units

This paper introduces a novel analytical approach for the identification of the admittance matrix and the generalized short-circuit ratio (gSCR) in power systems integrated with renewable energy sources. The proposed method leverages voltage and current measurements from phasor measurement units (PMUs) to construct a least squares objective function, which is then solved using matrix calculus and partial derivatives. Unlike conventional optimization algorithms, this approach provides an analytical solution that substantially reduces data requirements, enabling the efficient and accurate identification of the gSCR with smaller datasets. Additionally, its fixed computational complexity allows for real-time updates as new data are collected, ensuring continuous refinement of the system of equations and enabling rapid, precise gSCR calculations. The method also exhibits strong robustness against measurement noise, making it well-suited for practical applications in dynamic power systems. The combination of reduced data requirements, real-time adaptability, noise robustness and fixed computational load establishes this method as a highly efficient and reliable tool for real-time power system stability analysis. Case studies on an EPRI 36-bus system demonstrate the method's effectiveness, highlighting its accuracy in closely matching true gSCR values, even under diverse disturbances and noisy conditions.

Han, Zelei [Hohai University, Nanjing (China)] (OR↗

Static penetration resistance of soils

Model test results were used to define the failure mechanism associated with the static penetration resistance of cohesionless and low-cohesion soils. Knowledge of this mechanism has permitted the development of a new analytical method for calculating the ultimate penetration resistance which explicitly accounts for penetrometer base apex angle and roughness, soil friction angle, and the ratio of penetration depth to base width. Curves relating the bearing capacity factors to the soil friction angle are presented for failure in general shear. Strength parameters and penetrometer interaction properties of a fine sand were determined and used as the basis for prediction of the penetration resistance encountered by wedge, cone, and flat-ended penetrometers of different surface roughness using the proposed analytical method. Because of the close agreement between predicted values and values measured in laboratory tests, it appears possible to deduce in-situ soil strength parameters and their variation with depth from the results of static penetration tests.

Durgunoglu, H. T.↗

Atomistic determination of Peierls barriers of dislocation glide in nickel

The Peierls barrier measures the lattice resistance to dislocation glide in crystalline solids. We use the nudged elastic band (NEB) method to calculate the Peierls barriers for screw and edge dislocation glide in a face-centered cubic (FCC) metal of Ni. The minimum energy paths (MEPs) across single or sequential Peierls barriers are determined under shear loading. The NEB results show the decreasing Peierls barrier with increasing shear stress, giving the Peierls stress at which the Peierls barrier vanishes. The effects of boundary condition and system size on Peierls barriers are studied by comparing strain- and stress-controlled NEB results. Furthermore, the free-end NEB methods are applied to determine MEPs with improved computational efficiency. The NEB results are also used to evaluate the energetic driving force of dislocation glide, which is consistent with that determined from the Peach-Koehler force. The accuracy of the present NEB results based on an empirical interatomic potential is assessed by comparison with a machine-learning potential. This work demonstrates the robust and efficient quantification of Peierls barriers to dislocation glide in an FCC metal, and it lays a solid foundation for the atomistic determination of Peierls barriers in compositionally complex alloys with the FCC structure in future studies.

42 ENGINEERING↗

Thermodynamics Modeling for Actinide Monocarbides and Mononitrides from First Principles

The high-temperature thermodynamical properties for the actinide monocarbides and mononitrides ThC, ThN, UC, UN, PuC, and PuN are calculated from first-principles electronic-structure theory. The electronic structure is modeled with density-functional theory (DFT) and is fully relativistic, including the spin-orbit interaction. Furthermore, the DFT is extended to account for orbital–orbital interactions, by means of a parameter-free orbital-polarization (OP) technique, that has proven to be essential for the 5f electrons in plutonium. Strong anharmonicity and the temperature dependence of the lattice vibrations are captured with the self-consistent ab initio lattice dynamics (SCAILD) method. The calculated free energies and heat capacities are compared to published results from quasi-harmonic (QH) theory, and experiments, where available. For the uranium and plutonium compounds, we make use of CALPHAD assessments to help evaluate the theory. Generally, our anharmonic relativistic approach compares well with both CALPHAD and experiments. For the thorium compounds, our theory is in good accord with QH modeling of the free energy at lower temperatures but for the heat capacity the comparison is less favorable.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Computationally efficient control allocation

A computationally efficient method for calculating near-optimal solutions to the three-objective, linear control allocation problem is disclosed. The control allocation problem is that of distributing the effort of redundant control effectors to achieve some desired set of objectives. The problem is deemed linear if control effectiveness is affine with respect to the individual control effectors. The optimal solution is that which exploits the collective maximum capability of the effectors within their individual physical limits. Computational efficiency is measured by the number of floating-point operations required for solution. The method presented returned optimal solutions in more than 90% of the cases examined; non-optimal solutions returned by the method were typically much less than 1% different from optimal and the errors tended to become smaller than 0.01% as the number of controls was increased. The magnitude of the errors returned by the present method was much smaller than those that resulted from either pseudo inverse or cascaded generalized inverse solutions. The computational complexity of the method presented varied linearly with increasing numbers of controls; the number of required floating point operations increased from 5.5 i, to seven times faster than did the minimum-norm solution (the pseudoinverse), and at about the same rate as did the cascaded generalized inverse solution. The computational requirements of the method presented were much better than that of previously described facet-searching methods which increase in proportion to the square of the number of controls.

Durham, Wayne↗

A Lagrangian dual method for two-stage robust optimization with binary uncertainties

This report presents a new exact method to calculate worst-case parameter realizations in two-stage robust optimization problems with categorical or binary-valued uncertain data. Traditional exact algorithms for these problems, notably Benders decomposition and column-and-constraint generation, compute worst-case parameter realizations by solving mixed-integer bilinear optimization subproblems. However, their numerical solution can be computationally expensive not only due to their resulting large size after reformulating the bilinear terms, but also because decision-independent bounds on their variables are typically unknown. We propose an alternative Lagrangian dual method that circumvents these difficulties and is readily integrated in either algorithm. We specialize the method to problems where the binary parameters switch on or off constraints as these are commonly encountered in applications, and discuss extensions to problems that lack relatively complete recourse and to those with integer recourse. Numerical experiments provide evidence of significant computational improvements over existing methods.

42 ENGINEERING↗

Skyshine Calculations for a Large Spent Nuclear Fuel Storage Facility with SCALE 6.2.3

The SCALE code system developed at Oak Ridge National Laboratory includes state-of-the-art capabilities for radiation source term and radiation transport simulations that can be used in numerous applications, including dose rate analyses of complex consolidated interim storage facilities (CISFs). A licensed CISF could be used to store tens of thousands of tonnes of spent nuclear fuel discharged from commercial power reactors using various cask and storage pad designs. A CISF design must comply with the regulatory requirements provided in 10 CFR Part 72, including requirements related to annual dose limits applicable to real individuals located beyond the area controlled by the licensee. Therefore, calculating a dose to the public is a necessary part of the licensing process for the construction of a CISF. These calculations are very challenging because of the complexity of the CISF design and the low magnitude of dose rate at large distances from the facility. This paper describes detailed far-field dose rate calculations performed for a proposed CISF using MAVRIC, the Monte Carlo radiation shielding sequence in SCALE 6.2.3, with automated variance reduction based on discrete ordinates calculations. The method presented in this paper uses a detailed Monte Carlo radiation transport simulation in one step from source to dose rate. A series of independent simulations was made using the complete site geometry (all casks present), but with only one cask containing radiation sources to obtain the dose rate maps produced by each storage cask. The CISF dose rate map was obtained by adding the dose rate maps produced by the independent individual cask simulations. Ample volumes of air and soil extending beyond the location of interest for dose rate calculation were included in the calculation model to properly simulate important radiation attenuation and scattering events that affect far-field dose rates. Furthermore, a comprehensive sensitivity study is included in this paper to illustrate the importance of selecting appropriate air volume, mass density, and composition for CISF skyshine dose rate calculations. Dry soil and soil containing water were analyzed to determine their effects on groundshine radiation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Costs of measuring leaf area index of corn

The magnitude of plant-to-plant variability of leaf area of corn plants selected from uniform plots was examined and four representative methods for measuring leaf area index (LAI) were evaluated. The number of plants required and the relative costs for each sampling method were calculated to detect 10, 20, and 50% differences in LAI using 0.05 and 0.01 tests of significance and a 90% probability of success (beta = 0.1). The natural variability of leaf area per corn plant was nearly 10%. Additional variability or experimental error may be introduced by the measurement technique employed and by nonuniformity within the plot. Direct measurement of leaf area with an electronic area meter had the lowest CV, required that the fewest plants be sampled, but required approximately the same amount of time as the leaf area/weight ratio method to detect comparable differences. Indirect methods based on measurements of length and width of leaves required more plants but less total time than the direct method. Unless the coefficients for converting length and width to area are verified frequently, the indirect methods may be biased. When true differences in LAI among treatments exceed 50% of mean, all four methods are equal. The method of choice depends on the resources available, the differences to be detected, and what additional information, such as leaf weight or stalk weight, is also desired.

Daughtry, C. S. T.↗

An in-line burn-up algorithm: proposal and preliminary analysis

A new burn-up algorithm is proposed for coupling neutronics and isotopic depletion. The algorithm is analogous to those proposed for enforcing xenon equilibrium and converging neutronics and thermal-hydraulics. During the corrector step of standard predictor-corrector methods, the corrector nuclide density is updated as neutronics converges. This is expected to improve stability while reducing computational expense compared to other implicit algorithms. The stability of the algorithm is investigated by Fourier analysis for a simple case whose results have been successfully compared with numerical experiments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

PV O&M Cost Model [SWR-21-14]

This tool presents a method for calculating costs associated with the operation and maintenance (O&M) of photovoltaic (PV) systems. The tool compiles details regarding the cost and frequency of multiple O&M services to estimate annual O&M costs ($/year) for each year of an analysis period, the net present value ($) of life cycle costs accumulated over the analysis period, and the reserve account amount ($) that might be required to fund unexpected repairs. The reserve account includes parts inventory and is important for providing a source of funds to make repairs quickly and avoid lost production. This method allows a detailed selection of services to perform based on system size, market served (e.g., residential, commercial, or utility), type and configuration of system components (e.g., micro-, string, or central inverter), and site and environmental conditions (e.g., pollen, bird populations) which is an improvement over simple per unit valuations of O&M costs ($/kW/year). This model also distinguishes costs that vary from year to year and increase at different rates over time as modeled by heuristic failure distributions (e.g., Weibull or Lognormal distribution) based on actuarial data for many of the services.

Walker, Andy↗

A user's guide for V174, a program using a finite difference method to analyze transonic flow over oscillating wings

The design and usage of a pilot program using a finite difference method for calculating the pressure distributions over harmonically oscillating wings in transonic flow are discussed. The procedure used is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equation for small disturbances. The steady velocity potential which must be obtained from some other program, is required for input. The unsteady differential equation is linear, complex in form with spatially varying coefficients. Because sinusoidal motion is assumed, time is not a variable. The numerical solution is obtained through a finite difference formulation and a line relaxation solution method.

Butler, T. D.↗

Estimation of effective hydrologic properties of soils from observations of vegetation density

An existing one-dimensional model of the annual water balance is reviewed. Slight improvements are made in the method of calculating the bare soil component of evaporation, and in the way surface retention is handled. A natural selection hypothesis, which specifies the equilibrium vegetation density for a given, water limited, climate-soil system, is verified through comparisons with observed data and is employed in the annual water balance of watersheds in Clinton, Ma., and Santa Paula, Ca., to estimate effective areal average soil properties. Comparison of CDF's of annual basin yield derived using these soil properties with observed CDF's provides excellent verification of the soil-selection procedure. This method of parameterization of the land surface should be useful with present global circulation models, enabling them to account for both the non-linearity in the relationship between soil moisture flux and soil moisture concentration, and the variability of soil properties from place to place over the Earth's surface.

Tellers, T. E.↗

Computing the projected area of a spacecraft solar array with matrices

Matrix methods for computing the projected area of a solar array as a function of the rotational position of a spacecraft and of array position on a spacecraft are presented. Formulas are derived which provide the optimum solar array pitch, cant and tracking angles for a given spacecraft configuration and orbit. These formulas are general and applicable to many spacecraft. Formulas are also provided for determining the energy output from an array for a given orbit and the resultant energy available to spacecraft loads. Results are simply obtained and realizable with a hand calculator. The methods above can be extended to the case of computer analysis of solar array shadowing. How this can be done is outlined and results are presented from a spacecraft study program.

E. M. Gaddy↗

THERMAL SAFETY MARGIN CALCULATION OF THE MP-2 EXPERIMENT IN THE ADVANCED TEST REACTOR

The Mini-Plate 2 (MP-2) irradiation test is a fueled experiment designed for irradiation in multiple test locations in the Advanced Test Reactor (ATR) at the Idaho National Laboratory (INL). The experiment is a drop-in test where small aluminum-clad fuel plate samples (mini plates) are cooled directly by the ATR Primary Coolant System (PCS) water. The MP-2 fuel plate experiment will be irradiated in several different irradiation locations of the ATR. This fueled experiment contains aluminum-clad fuel mini plates consisting of monolithic U-Mo. Four different types of fuel plates with fuel meat thickness and cladding are part of the MP-2 test. A thermal analysis has been performed on the MP-2 experiment. A method for calculating Departure from Nucleate Boiling Ratio (DNBR) and Flow Instability Ratio (FIR) during a reactivity transient using the commercial finite element and heat transfer code ABAQUS is discussed. At the start of an ATR cycle the heat generation rate of the fueled experiment is high and the heat rate multiplier from the outer shim control cylinders is low, while the reverse is true at the end of the ATR cycle. Thermal analyses at 10 day increments during the cycle calculate the DNBR and FIR during a reactivity transient. This technique calculates DNBR for the fuel plate surfaces and FIR for all water components for each finite element surface and node at various times during the ATR cycle. Heat rates vary with time during the transient calculation that are provided by a detailed physics analysis. Oxide growth on the fuel plates is also incorporated Results from the transient calculations are displayed with the ABAQUS post processor. By calculating these parameters at each location in the finite element model, conservatism is replaced with accuracy. This allows for a greater margin for the thermal hydraulic safety parameters.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Weak gravitational lensing by Kerr-MOG black hole and Gauss–Bonnet theorem

Highlights: • The deflection angle of the Kerr-MOG black hole is significantly greater than the original Kerr black hole. • In the future astrophysical observation might shed light on the effects of modified gravity theory. • MOG effect could be taken into account in the gravitational lensing experiments. • Our results put forward that many features of weak gravitational lensing arise from the topological properties of an underlying optical manifold. • For weak gravitational lensing of (at least) spherically symmetric distributions, the Gauss–Bonnet approach holds several advantages over common methods. - Abstract: The deflection angle of Kerr-MOG black holes is studied for different values of the parameter in modified gravity (MOG). To this end, we employ the Gauss–Bonnet theorem, which was first studied by Gibbons and Werner and then extended by Ono, Ishihara and Asada, who use a generalized optical metric where the deflection of light for an observer and source at finite distance. By using this method, we study the weak gravitational lensing by Kerr-MOG black hole. Our computations show that with an increase in the MOG parameter (α), the deflection angle becomes significantly larger than that of Kerr black hole. The results obtained show that MOG effect could be taken into account in the gravitational lensing experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

Improved Method for Evaluating and Specifying the Chromaticity of Light Sources

This article recommends a method for calculating and specifying light source chromaticity using the CIE 2015 10° color matching functions (CMFs), which, according to analysis of existing psychophysical experiment data, can reduce visual mismatch compared to specifications based on the traditional CIE 1931 2° CMFs. Specifically, this work evaluates, documents, and recommends for adoption a supporting system of measures to be used with the CIE 2015 10° CMFs: a new uniform chromaticity scale (UCS) diagram with coordinates (s, t), a measure of correlated color temperature (CCTst), and a measure of distance from the Planckian locus (Dst). It also presents options for updating nominal classification quadrangles. A complete method of this nature has not yet been standardized, which may be contributing to the slow uptake of the CIE 2015 CMFs. The proposed tools are analogous to u, v, CCT, Duv, and the ANSI C78.377 nominal classification quadrangles that are currently specified in the CIE 1960 UCS using the CIE 1931 2° CMFs. While conceptually equivalent, the differences between the current standard method and the new st system are important for reducing unintended mismatch in the chromaticity of light. Furthermore, the implications of changing chromaticity specification methods are identified by a comparison over a diverse set of real light source SPDs. Finally, as a step toward improving the design and specification of light sources for illumination purposes, we propose a process for implementing the recommendations.

42 ENGINEERING↗