Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computing Methodologies”

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 325 records · Page 18

An overview of methods to identify and manage uncertainty for modelling problems in the water-environment-agriculture cross-sector

Uncertainty pervades the representation of systems in the water–environment–agriculture cross-sector. Successful methods to address uncertainties have largely focused on standard mathematical formulations of biophysical processes in a single sector, such as partial or ordinary differential equations. More attention to integrated models of such systems is warranted. Model components representing the different sectors of an integrated model can have less standard, and different, formulations to one another, as well as different levels of epistemic knowledge and data informativeness. Thus, uncertainty is not only pervasive but also crosses boundaries and propagates between system components. Uncertainty assessment (UA) cries out for more eclectic treatment in these circumstances, some of it being more qualitative and empirical. Here in this paper, we discuss the various sources of uncertainty in such a cross-sectoral setting and ways to assess and manage them. We have outlined a fast-growing set of methodologies, particularly in the computational mathematics literature on uncertainty quantification (UQ), that seem highly pertinent for uncertainty assessment. There appears to be considerable scope for advancing UA by integrating relevant UQ techniques into cross-sectoral problem applications. Of course this will entail considerable collaboration between domain specialists who often take first ownership of the problem and computational methods experts.

97 MATHEMATICS AND COMPUTING↗

Convergence of micro-geochemistry and micro-geomechanics towards understanding proppant shale rock interaction: A Caney shale case study in southern Oklahoma, USA

As a direct outcome of economic development coupled with an increase in population, global energy demand will continue to rise in the coming decades. Although renewable energy sources are increasingly investigated for optimal production, the immediate needs require focus on energy sources that are currently available and reliable, with a minimal environmental impact; the efficient exploration and production of unconventional hydrocarbon resources is bridging the energy needs and energy aspirations, during the current energy transition period. The main challenges are related to the accurate quantification of the critical rock properties that influence production, their heterogeneity and the multiscale driven physico-chemical nature of rock–fluid interactions. A key feature of shale reservoirs is their low permeability due to dominating nanoporosity of the clay-rich matrix. As a means of producing these reservoirs in a cost-effective manner, a prerequisite is creation of hydraulic fracture networks capable of the highest level of continued conductivity. Fracturing fluid chemical design, formation brine geochemical composition, and rock mineralogy all contribute to swelling-induced conductivity damage. The Caney Shale is an organic-rich, often calcareous mudrock. Many studies have examined the impact that clay has on different kinds of shale productivity but there is currently no data reported on the Caney Shale in relation to horizontal drilling; all reported data on the Caney Shale is on vertical wells which are shallow, compared to an emerging play that is at double the depth. Here, in this work we develop geochemical–geomechanical integration of rock properties at micro-and nanoscales that can provide insights into the potential proppant embedment and its mitigation. The novel methodology amalgamates the following: computed X-ray tomography, scanning electron microscopy, energy dispersive spectroscopy, micro-indentation, and Raman spectroscopy techniques. Here, our results show that due to the multiscale heterogeneity in the Caney Shale, these geochemical and structural properties translate into a variation in mechanical properties that will impact interaction between the proppant and the host shale rock.

03 NATURAL GAS↗

Automated high-throughput tensile testing reveals stochastic process parameter sensitivity

The mechanical properties of additively manufactured metals tend to show high variability, due largely to the stochastic nature of defect formation during the printing process. This study seeks to understand how automated high throughput testing can be utilized to understand the variable nature of additively manufactured metals at different print conditions, and to allow for statistically meaningful analysis. This is demonstrated by analyzing how different processing parameters, including laser power, scan velocity, and scan pattern, influence the tensile behavior of additively manufactured stainless steel 316L utilizing a newly developed automated test methodology. Microstructural characterization through computed tomography and electron backscatter diffraction is used to understand some of the observed trends in mechanical behavior. Specifically, grain size and morphology are shown to depend on processing parameters and influence the observed mechanical behavior. In the current study, laser-powder bed fusion, also known as selective laser melting or direct metal laser sintering, is shown to produce 316L over a wide processing range without substantial detrimental effect on the tensile properties. Ultimate tensile strengths above 600 MPa, which are greater than that for typical wrought annealed 316L with similar grain sizes, and elongations to failure greater than 40% were observed. As a result, it is demonstrated that this process has little sensitivity to minor intentional or unintentional variations in laser velocity and power.

36 MATERIALS SCIENCE↗

Reliable p K a Prediction through Efficient Incorporation of Anharmonicity within the Nuclear–Electronic Orbital Framework

Accurate pK a prediction is critical for understanding chemical reactivity and molecular properties across a wide range of applications. Computational approaches usually invoke a harmonic treatment of the vibrational modes for zero-point energies, as well as thermal and entropic contributions. Herein, we present a general protocol for relative pK a prediction that incorporates the significant anharmonic effects using nuclear–electronic orbital (NEO) theory. This protocol is validated against experimental data for a range of molecules in acetonitrile, including protonated nitrogen bases, nitrophenols, anilines, and diamines, as well as cobalt electrocatalysts. For simple acids, the NEO approach offers only a slight improvement over conventional density functional theory with the standard harmonic vibrational treatment, whereas for hydrogen-bonded acids, the NEO approach offers more significantly improved performance at a comparable computational cost. This accessible methodology provides a practical route for accurate pKa prediction in challenging systems and is extendable to related thermodynamic properties such as hydricities and proton-coupled redox potentials.

Density functional theory↗

Data-driven quantum approximate optimization algorithm for power systems

Abstract Quantum technology provides a ground-breaking methodology to tackle challenging computational issues in power systems. It is especially promising for Distributed Energy Resources (DERs) dominant systems that have been widely developed to promote energy sustainability. In those systems, knowing the maximum sections of power and data delivery is essential for monitoring, operation, and control. However, high computational effort is required. By leveraging quantum resources, Quantum Approximate Optimization Algorithm (QAOA) provides a means to search for these sections efficiently. However, QAOA performance relies heavily on critical parameters, especially for weighted graphs. Here we present a data-driven QAOA, which transfers quasi-optimal parameters between weighted graphs based on the normalized graph density. We verify the strategy with 39,774 expectation value calculations. Without parameter optimization, our data-driven QAOA is comparable with the Goemans-Williamson algorithm. This work advances QAOA and pilots its practical application to power systems in noisy intermediate-scale quantum devices.

Jing, Hang (ORCID:0000000256343489)↗

Capacity Design Method for Integrated Convective/Radiant Terminals to Guarantee Overall and Local Environment

To achieve energy efficiency in the operation of heating devices and ensure thermal comfort in indoor environments, the integrated convective/radiant terminals have become an important development direction for heating systems. This study selected a 14 m 2 bedroom hot summer and cold winter region in China as a case study. Computational Fluid Dynamics (CFD) methodology was used to investigate the design of integrated convective/radiant terminals to ensure thermal comfort for both the entire room and partial space. The results indicate that focusing on partial space during the steady-state stage can achieve energy savings of 31.1% compared to guaranteeing the entire room. Additionally, during the start-up stage, there is a significant reduction of convective unit capacity by 17.6% when the start-up time is 15 minutes. These findings provide data support for the design and engineering applications of products of integrated convective/radiant terminals. Besides, the excess convection ratios for the entire room and partial space are specified as 0.69 to 0.80 and 0.73 to 0.83, respectively.

Chi, Junjie↗

Generalized correlation integral vectors: A distance concept for chaotic dynamical systems

Several concepts of fractal dimension have been developed to characterise properties of attractors of chaotic dynamical systems. Numerical approximations of them must be calculated by finite samples of simulated trajectories. In principle, the quantities should not depend on the choice of the trajectory, as long as it provides properly distributed samples of the underlying attractor. In practice, however, the trajectories are sensitive with respect to varying initial values, small changes of the model parameters, to the choice of a solver, numeric tolerances, etc. The purpose of this paper is to present a statistically sound approach to quantify this variability. We modify the concept of correlation integral to produce a vector that summarises the variability at all selected scales. The distribution of this stochastic vector can be estimated, and it provides a statistical distance concept between trajectories. Here, we demonstrate the use of the distance for the purpose of estimating model parameters of a chaotic dynamic model. The methodology is illustrated using computational examples for the Lorenz 63 and Lorenz 95 systems, together with a framework for Markov chain Monte Carlo sampling to produce posterior distributions of model parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Improving accuracy of electrochemical capacitance and solvation energetics in first-principles calculations

Here, reliable first-principles calculations of electrochemical processes require accurate prediction of the interfacial capacitance, a challenge for current computationally efficient continuum solvation methodologies. We develop a model for the double layer of a metallic electrode that reproduces the features of the experimental capacitance of Ag(100) in a non-adsorbing, aqueous electrolyte, including a broad hump in the capacitance near the potential of zero charge and a dip in the capacitance under conditions of low ionic strength. Using this model, we identify the necessary characteristics of a solvation model suitable for first-principles electrochemistry of metal surfaces in non-adsorbing, aqueous electrolytes: dielectric and ionic nonlinearity, and a dielectric-only region at the interface. The dielectric nonlinearity, caused by the saturation of dipole rotational response in water, creates the capacitance hump, while ionic nonlinearity, caused by the compactness of the diffuse layer, generates the capacitance dip seen at low ionic strength. We show that none of the previously developed solvation models simultaneously meet all these criteria. We design the nonlinear electrochemical soft-sphere solvation model which both captures the capacitance features observed experimentally and serves as a general-purpose continuum solvation model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

THE NASA AMES PAH IR SPECTROSCOPIC DATABASE VERSION 2.00: UPDATED CONTENT, WEB SITE, AND ON(OFF)LINE TOOLS

A significantly updated version of the NASA Ames PAH IR Spectroscopic Database, the first major revision since its release in 2010, is presented. The current version, version 2.00, contains 700 computational and 75 experimental spectra compared, respectively, with 583 and 60 in the initial release. The spectra span the 2.5-4000 μm (4000-2.5 cm{sup -1}) range. New tools are available on the site that allow one to analyze spectra in the database and compare them with imported astronomical spectra as well as a suite of IDL object classes (a collection of programs utilizing IDL's object-oriented programming capabilities) that permit offline analysis called the AmesPAHdbIDLSuite. Most noteworthy among the additions are the extension of the computational spectroscopic database to include a number of significantly larger polycyclic aromatic hydrocarbons (PAHs), the ability to visualize the molecular atomic motions corresponding to each vibrational mode, and a new tool that allows one to perform a non-negative least-squares fit of an imported astronomical spectrum with PAH spectra in the computational database. Finally, a methodology is described in the Appendix, and implemented using the AmesPAHdbIDLSuite, that allows the user to enforce charge balance during the fitting procedure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theory of non-Gaussianity in warm inflation

The theory and methodology is developed to compute the bispectrum in warm inflation, leading to results for the non-linearity parameter and the shape of the bispectrum. Particular attention is paid to the study of the bispectrum in the regime of weak dissipation and how stochastic fluctuations affect the bispectrum. It is shown that, in contrast to the strong dissipative regime, the amplitude of non-Gaussianity is strongly dependent on the parameters governing the microscopic physics in the intermediate and weak dissipation warm inflation regimes. The most important results concern the shape of the bispectrum, which has two different, but distinct, forms in the weak and strong dissipative regimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Predicting nonequilibrium Green’s function dynamics and photoemission spectra via nonlinear integral operator learning

Understanding the dynamics of nonequilibrium quantum many-body systems is an important research topic in a wide range of fields across condensed matter physics, quantum optics, and high-energy physics. However, numerical studies of large-scale nonequilibrium phenomena in realistic materials face serious challenges due to intrinsic high-dimensionality of quantum many-body problems and the absence of time-invariance. The nonequilibrium properties of many-body systems can be described by the dynamics of the correlator, or the Green's function of the system, whose time evolution is given by a high-dimensional system of integro-differential equations, known as the Kadanoff–Baym equations (KBEs). The time-convolution term in KBEs, which needs to be recalculated at each time step, makes it difficult to perform long-time numerical simulation. In this paper, we develop an operator-learning framework based on recurrent neural networks (RNNs) to address this challenge. We utilize RNNs to learn the nonlinear mapping between Green's functions and convolution integrals in KBEs. By using the learned operators as a surrogate model in the KBE solver, we obtain a general machine-learning scheme for predicting the dynamics of nonequilibrium Green's functions. Besides significant savings per each time step, the new methodology reduces the temporal computational complexity from $O(N_t^3)$ to $O(N_t)$ where N t is the number of steps taken in a simulation, thereby making it possible to study large many-body problems which are currently infeasible with conventional KBE solvers. Through various numerical examples, we demonstrate the effectiveness of the operator-learning based approach in providing accurate predictions of physical observables such as the reduced density matrix and time-resolved photoemission spectra. Moreover, our framework exhibits clear numerical convergence and can be easily parallelized, thereby facilitating many possible further developments and applications.

97 MATHEMATICS AND COMPUTING↗

Algorithm sensitivity analysis and parameter tuning for tissue image segmentation pipelines

Motivation: Sensitivity analysis and parameter tuning are important processes in large-scale image analysis. They are very costly because the image analysis workflows are required to be executed several times to systematically correlate output variations with parameter changes or to tune parameters. An integrated solution with minimum user interaction that uses effective methodologies and high performance computing is required to scale these studies to large imaging datasets and expensive analysis workflows. Results: The experiments with two segmentation workflows show that the proposed approach can (i) quickly identify and prune parameters that are non-influential; (ii) search a small fraction (about 100 points) of the parameter search space with billions to trillions of points and improve the quality of segmentation results (Dice and Jaccard metrics) by as much as 1.42 x compared to the results from the default parameters; (iii) attain good scalability on a high performance cluster with several effective optimizations. Conclusions: Our work demonstrates the feasibility of performing sensitivity analyses, parameter studies and auto-tuning with large datasets. The proposed framework can enable the quantification of error estimations and output variations in image segmentation pipelines.

Biochemistry & Molecular Biology↗

New Approach Methodology for Assessing Inhalation Risks of a Contact Respiratory Cytotoxicant: Computational Fluid Dynamics-Based Aerosol Dosimetry Modeling for Cross-Species and In Vitro Comparisons

Regulatory agencies are considering alternative approaches to assessing inhalation toxicity that utilizes in vitro studies with human cells and in silico modeling in lieu of additional animal studies. In support of this goal, computational fluid-particle dynamics models were developed to estimate site-specific deposition of inhaled aerosols containing the fungicide, chlorothalonil, in the rat and human for comparisons to prior rat inhalation studies and new human in vitro studies. Under bioassay conditions, the deposition was predicted to be greatest at the front of the rat nose followed by the anterior transitional epithelium and larynx corresponding to regions most sensitive to local contact irritation and cytotoxicity. For humans, simulations of aerosol deposition covering potential occupational or residential exposures (1–50 µm diameter) were conducted using nasal and oral breathing. Aerosols in the 1–5 µm range readily penetrated the deep region of the human lung following both oral and nasal breathing. Under actual use conditions (aerosol formulations >10 µm), the majority of deposited doses were in the upper conducting airways. Beyond the nose or mouth, the greatest deposition in the pharynx, larynx, trachea, and bronchi was predicted for aerosols in the 10–20 µm size range. Only small amounts of aerosols >20 µm penetrated past the pharyngeal region. Using the ICRP clearance model, local retained tissue dose metrics including maximal concentrations and areas under the curve were calculated for each airway region following repeated occupational exposures. These results are directly comparable with benchmark doses from in vitro toxicity studies in human cells leading to estimated human equivalent concentrations that reduce the reliance on animals for risk assessments.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Temperature dependence of magnetic anisotropy and magnetoelasticity from classical spin-lattice calculations

Here we present a classical molecular-spin dynamics (MSD) methodology that enables accurate computations of the temperature dependence of the magnetocrystalline anisotropy as well as magnetoelastic properties of magnetic materials. The nonmagnetic interactions are accounted for by a spectral neighbor analysis potential (SNAP) machine-learned interatomic potential, whereas the magnetoelastic contributions are accounted for using a combination of an extended Heisenberg Hamiltonian and a Néel pair interaction model, representing both the exchange interaction and spin-orbit-coupling effects, respectively. All magnetoelastic potential components are parameterized using a combination of first-principles and experimental data. Our framework is applied to the α phase of iron. Initial testing of our MSD model is done using a 0 K parametrization of the Néel interaction model. After this, we examine how individual Néel parameters impact the $B$ 1 and $B$ 2 magnetostrictive coefficients using a moment-independent δ sensitivity analysis. The results from this study are then used to initialize a genetic algorithm optimization which explores the Néel parameter phase space and tries to minimize the error in the B 1 and B 2 magnetostrictive coefficients in the range of 0–1200 K. Our results show that while both the 0 K and genetic algorithm optimized parametrization provide good experimental agreement for $B$ 1 and $B$ 2 , only the genetic algorithm optimized results can capture the second peak in the $B$ 1 magnetostrictive coefficient which occurs near approximately 800 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Predictions of first passage times in sparse discrete fracture networks using graph-based reductions

Here, we present a graph-based methodology to reduce the computational cost of obtaining first passage times through sparse fracture networks. We also derive graph representations of generic three-dimensional discrete fracture networks (DFNs) using the DFN topology and flow boundary conditions. Subgraphs corresponding to the union of the k shortest paths between the inflow and outflow boundaries are identified and transport on their equivalent subnetworks is compared to transport through the full network. The number of paths included in the subgraphs is based on the scaling behavior of the number of edges in the graph with the number of shortest paths. First passage times through the subnetworks are in good agreement with those obtained in the full network, both for individual realizations and in distribution. We obtain accurate estimates of first passage times with an order of magnitude reduction of CPU time and mesh size using the proposed method.

54 ENVIRONMENTAL SCIENCES↗