Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complexity theory”

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 523 records · Page 29

Effective permittivity of saline ice under thermal variation

A model for calculating the effective permittivity of saline ice under thermal variation is presented. The model includes multiphase inhomogeneities with multiple species characterized by orientation, size and shape distributions. The model is used to derive the effective permittivity as a function of temperature under the strong fluctuation theory which is extended to account for the complexity. The results calculated from the model are compared with experimental data at 4.8 GHz for saline ice grown at the US Army Cold Regions Research and Engineering Laboratory (CRREL). The comparison between measured and calculated complex permittivities is good for the imaginary part, and the difference is within 10 percent for the real part.

Nghiem, S. V.↗

Molecular Spectroscopy by Ab Initio Methods

Due to recent advances in methods and computers, the accuracy of ab calculations has reached a point where these methods can be used to provide accurate spectroscopic constants for small molecules; this will be illustrated with several examples. We will show how ab initio calculations where used to identify the Hermann infrared system in N2 and two band systems in CO. The identification of all three of these band systems relied on very accurate calculations of quintet states. The analysis of the infrared spectra of cool stars requires knowledge of the intensity of vibrational transitions in SiO for high nu and J levels. While experiment can supply very accurate dipole moments for nu = 0 to 3, this is insufficient to construct a global dipole moment function. We show how theory, combined by the experiment, can be used to generate the line intensities up to nu = 40 and J = 250. The spectroscopy of transition metal containing systems is very difficult for both theory and experiment. We will discuss the identification of the ground state of Ti2 and the spectroscopy of AlCu as examples of how theory can contribute to the understanding of these complex systems.

Bauschlicher, Charles W., Jr.↗

An Efficient Scheme for Updating Sparse Cholesky Factors

Raghavan had earlier developed the software package DCSPACK which can be used for solving sparse linear systems where the coefficient matrix is symmetric and positive definite (this project was not funded by NASA but by agencies such as NSF). DSCPACK-S is the serial code and DSCPACK-P is a parallel implementation suitable for multiprocessors or networks-of-workstations with message passing using MCI. The main algorithm used is the Cholesky factorization of a sparse symmetric positive positive definite matrix A = LL(T). The code can also compute the factorization A = LDL(T). The complexity of the software arises from several factors relating to the sparsity of the matrix A. A sparse N x N matrix A has typically less that cN nonzeroes where c is a small constant. If the matrix were dense, it would have O(N2) nonzeroes. The most complicated part of such sparse Cholesky factorization relates to fill-in, i.e., zeroes in the original matrix that become nonzeroes in the factor L. An efficient implementation depends to a large extent on complex data structures and on techniques from graph theory to reduce, identify, and manage fill. DSCPACK is based on an efficient multifrontal implementation with fill-managing algorithms and implementation arising from earlier research by Raghavan and others. Sparse Cholesky factorization is typically a four step process: (1) ordering to compute a fill-reducing numbering, (2) symbolic factorization to determine the nonzero structure of L, (3) numeric factorization to compute L, and, (4) triangular solution to solve L(T)x = y and Ly = b. The first two steps are symbolic and are performed using the graph of the matrix. The numeric factorization step is of dominant cost and there are several schemes for improving performance by exploiting the nested and dense structure of groups of columns in the factor. The latter are aimed at better utilization of the cache-memory hierarchy on modem processors to prevent cache-misses and provide execution rates (operations/second) that are close to the peak rates for dense matrix computations. Currently, EPISCOPACY is being used in an application at NASA directed by J. Newman and M. James. We propose the implementation of efficient schemes for updating the LL(T) or LDL(T) factors computed in DSCPACK-S to meet the computational requirements of their project. A brief description is provided in the next section.

Raghavan, Padma↗

Pythium invasion of plant-based life support systems: biological control and sources

Invasion of plant-based life support systems by plant pathogens could cause plant disease and disruption of life support capability. Root rot caused by the fungus, Pythium, was observed during tests of prototype plant growth systems containing wheat at the Kennedy Space Center (KSC). We conducted experiments to determine if the presence of complex microbial communities in the plant root zone (rhizosphere) resisted invasion by the Pythium species isolated from the wheat root. Rhizosphere inocula of different complexity (as assayed by community-level physiological profile: CLPP) were developed using a dilution/extinction approach, followed by growth in hydroponic rhizosphere. Pythium growth on wheat roots and concomitant decreases in plant growth were inversely related to the complexity of the inocula during 20-day experiments in static hydroponic systems. Pythium was found on the seeds of several different wheat cultivars used in controlled environmental studies, but it is unclear if the seed-borne fungal strain(s) were identical to the pathogenic strain recovered from the KSC studies. Attempts to control pathogens and their effects in hydroponic life support systems should include early inoculation with complex microbial communities, which is consistent with ecological theory.

NASA Discipline Life Support Systems↗

Nonlinear dynamics in cardiac conduction

Electrical conduction in the heart shows many phenomena familiar from nonlinear dynamics. Among these phenomena are multiple basins of attraction, phase locking, and perhaps period-doubling bifurcations and chaos. We describe a simple cellular-automation model of electrical conduction which simulates normal conduction patterns in the heart as well as a wide range of disturbances of heart rhythm. In addition, we review the application of percolation theory to the analysis of the development of complex, self-sustaining conduction patterns.

NASA Discipline Cardiopulmonary↗

Efficient Web Services Policy Combination

Large-scale Web security systems usually involve cooperation between domains with non-identical policies. The network management and Web communication software used by the different organizations presents a stumbling block. Many of the tools used by the various divisions do not have the ability to communicate network management data with each other. At best, this means that manual human intervention into the communication protocols used at various network routers and endpoints is required. Developing practical, sound, and automated ways to compose policies to bridge these differences is a long-standing problem. One of the key subtleties is the need to deal with inconsistencies and defaults where one organization proposes a rule on a particular feature, and another has a different rule or expresses no rule. A general approach is to assign priorities to rules and observe the rules with the highest priorities when there are conflicts. The present methods have inherent inefficiency, which heavily restrict their practical applications. A new, efficient algorithm combines policies utilized for Web services. The method is based on an algorithm that allows an automatic and scalable composition of security policies between multiple organizations. It is based on defeasible policy composition, a promising approach for finding conflicts and resolving priorities between rules. In the general case, policy negotiation is an intractable problem. A promising method, suggested in the literature, is when policies are represented in defeasible logic, and composition is based on rules for non-monotonic inference. In this system, policy writers construct metapolicies describing both the policy that they wish to enforce and annotations describing their composition preferences. These annotations can indicate whether certain policy assertions are required by the policy writer or, if not, under what circumstances the policy writer is willing to compromise and allow other assertions to take precedence. Meta-policies are specified in defeasible logic, a computationally efficient non-monotonic logic developed to model human reasoning. One drawback of this method is that at one point the algorithm starts an exhaustive search of all subsets of the set of conclusions of a defeasible theory. Although the propositional defeasible logic has linear complexity, the set of conclusions here may be large, especially in real-life practical cases. This phenomenon leads to an inefficient exponential explosion of complexity. The current process of getting a Web security policy from combination of two meta-policies consists of two steps. The first is generating a new meta-policy that is a composition of the input meta-policies, and the second is mapping the meta-policy onto a security policy. The new algorithm avoids the exhaustive search in the current algorithm, and provides a security policy that matches all requirements of the involved metapolicies.

Vatan, Farrokh↗

Initiating Molecular Growth in the Interstellar Medium via Dimeric Complexes of Observed Ions and Molecules

A feasible initiation step for particle growth in the interstellar medium (ISM) is simulated by means of ab quantum chemistry methods. The systems studied are dimer ions formed by pairing nitrogen containing small molecules known to exist in the ISM with ions of unsaturated hydrocarbons or vice versa. Complexation energies, structures of ensuing complexes and electronic excitation spectra of the encounter complexes are estimated using various quantum chemistry methods. Moller-Plesset perturbation theory (MP2, Z-averaged perturbation theory (ZAP2), coupled cluster singles and doubles with perturbative triples corrections (CCSD(T)), and density functional theory (DFT) methods (B3LYP) were employed along with the correlation consistent cc-pVTZ and aug-cc-pVTZ basis sets. Two types of complexes are predicted. One type of complex has electrostatic binding with moderate (7-20 kcal per mol) binding energies, that are nonetheless significantly stronger than typical van der Waals interactions between molecules of this size. The other type of complex develops strong covalent bonds between the fragments. Cyclic isomers of the nitrogen containing complexes are produced very easily by ion-molecule reactions. Some of these complexes show intense ultraviolet visible spectra for electronic transitions with large oscillator strengths at the B3LYP, omegaB97, and equations of motion coupled cluster (EOM-CCSD) levels. The open shell nitrogen containing carbonaceous complexes especially exhibit a large oscillator strength electronic transition in the visible region of the electromagnetic spectrum.

Bera, Partha P.↗

Towards Accurate and Efficient Predictions of Martensitic Transition Temperatures for Shape Memory Alloys from First Principles

Shape memory alloys (SMAs) can remember and recover their original shapes upon heating due to the existence of a reversible martensitic transition (MT) between the high-temperature austenite (A) and low-temperature martensite (M) phases. The martensitic transition temperature (MTT) is a crucial characteristic of an SMA. SMAs have a wide range of potential applications in aerospace, civil engineering, bioengineering, etc., but their operating temperatures are limited by the available SMAs. MTT can be tuned by alloying a binary with other metals, and the multicomponent NiTi-based SMAs have attracted tremendous research efforts recently. It is not efficient to employ the trial-and-error method alone due to the dramatically increased complexity and possibilities in compositions, and thus reliable theory and accurate computations play an indispensable role in creating SMAs with desirable properties.

Zhigang Wu↗

Approach for Inferring Full-Scope Human Reliability Data Based on Simplified Simulator Data

This paper proposes a method for inferring full-scope human reliability data based on the Simplified Human Error Experimental Program (SHEEP) data. It mainly focuses on the human errors observed when using simulators with different complexity levels. In the proposed method, the manner in which human error probabilities (HEPs) change as a result of increasing simulator complexity and how simulator complexity levels are quantified represent key information for inferring full-scope data. In the present study, SHEEP error data pertaining to actual professional operators using Rancor Microworld (Rancor) (i.e., a more simplified simulator) and Compact Nuclear Simulator (CNS) (i.e., a less simplified simulator) were compared with the HuREX error data. An approach to quantifying simulator complexity levels was then proposed based on information theory and acquired eye-tracker data.

99 - GENERAL AND MISCELLANEOUS↗

Fidelity Analysis of Sampled Imaging Systems

Many modeling, simulation and performance analysis studies of sampled imaging systems are inherently incomplete because they are conditioned on a discrete-input, discrete-output model that only accounts for blurring during image acquisition and additive noise. For those sampled imaging systems where the effects of digital image acquisition, digital filtering and reconstruction are significant, the modeling, simulation and performance analysis should be based on a more comprehensive continuous-input, discrete-processing, continuous-output end-to-end model. This more comprehensive model should properly account for the low-pass filtering effects of image acquisition prior to sampling, the potentially important noiselike effects of the aliasing caused by sampling, additive noise due to device electronics and quantization, the generally high-boost filtering effects of digital processing, and the low-pass filtering effects of image reconstruction. This model should not, however, be so complex as to preclude significant mathematical analysis, particularly the mean-square (fidelity) type of analysis so common in linear system theory. We demonstrate that, although the mathematics of such a model is more complex, the increase in complexity is not so great as to prevent a complete fidelity-metric analysis at both the component level and at the end-to-end system level: that is, computable mean-square-based fidelity metrics are developed by which both component-level and system-level performance can be quantified. In addition, we demonstrate that system performance can be assessed qualitatively by visualizing the output image as the sum of three component images, each of which relates to a corresponding fidelity metric. The cascaded, or filtered, component accounts for the end-to-end system filtering of image acquisition, digital processing, and image reconstruction; the random noise component accounts for additive random noise, modulated by digital processing and image reconstruction filtering; and the aliased noise component accounts for the frequency folding effect of sampling, modulated by digital processing and image reconstruction filtering.

Park, Stephen K.↗

Overdestabilization vs Overstabilization in the Theoretical Analysis of f-Orbital Covalency

The complex nature of the f-orbital electronic structures and their interaction with the chemical environment pose significant computational challenges. Advanced computational techniques that variationally include scalar relativities and spin–orbit coupling directly at the molecular orbital level have been developed to address this complexity. Among these, variational relativistic multiconfigurational multireference methods stand out for their high accuracy and systematic improvement in studies of f-block complexes. Additionally, these advanced methods offer the potential for calibrating low-scaling electronic structure methods such as density functional theory. However, studies on the Cl K-edge X-ray absorption spectra of the [Ce(III)Cl 6 ] 3– and [Ce(IV)Cl 6 ] 2– complexes show that time-dependent density functional theory with approximate exchange–correlation kernels can lead to inaccuracies, resulting in an overstabilization of 4f orbitals and incorrect assessments of covalency. In contrast, approaches utilizing small active space wave function methods may understate the stability of these orbitals. The results herein demonstrate the need for large active space, multireference, and variational relativistic methods in studying f-block complexes.

chemical calculations↗

Three-dimensional potential flows from functions of a 3D complex variable

Potential, or ideal, flow velocities can be found from the gradient of an harmonic function. An ordinary complex valued analytic function can be written as the sum of two real valued functions, both of which are harmonic. Thus, 2D complex valued functions serve as a source of functions that describe two-dimensional potential flows. However, this use of complex variables has been limited to two-dimensions. Recently, a new system of three-dimensional complex variables has been developed at the NASA Ames Research Center. As a step toward application of this theory to the analysis of 3D potential flow, several functions of a three-dimensional complex variable have been investigated. The results for two such functions, the 3D exponential and 3D logarithm, are presented in this paper. Potential flows found from these functions are investigated. Important characteristics of these flows fields are noted.

Kelly, Patrick↗

Application of Quantum Computing to Biochemical Systems: A Look to the Future

Chemistry is considered as one of the more promising applications to science of near-term quantum computing. Recent work in transitioning classical algorithms to a quantum computer has led to great strides in improving quantum algorithms and illustrating their quantum advantage. Because of the limitations of near-term quantum computers, the most effective strategies split the work over classical and quantum computers. There is a proven set of methods in computational chemistry and materials physics that has used this same idea of splitting a complex physical system into parts that are treated at different levels of theory to obtain solutions for the complete physical system for which a brute force solution with a single method is not feasible. These methods are variously known as embedding, multi-scale, and fragment techniques and methods. We review these methods and then propose the embedding approach as a method for describing complex biochemical systems, with the parts not only treated with different levels of theory, but computed with hybrid classical and quantum algorithms. Such strategies are critical if one wants to expand the focus to biochemical molecules that contain active regions that cannot be properly explained with traditional algorithms on classical computers. While we do not solve this problem here, we provide an overview of where the field is going to enable such problems to be tackled in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A comparison of Gaussian processes and neural networks for computer model emulation and calibration

The Department of Energy relies on complex physics simulations for prediction in domains like cosmology, nuclear theory, and materials science. These simulations are often extremely computationally intensive, with some requiring days or weeks for a single simulation. In order to assure their accuracy, these models are calibrated against observational data in order to estimate inputs and systematic biases. Because of their great computational complexity, this process typically requires the construction of an emulator, a fast approximation to the simulation. In this paper, two emulator approaches are compared: Gaussian process regression and neural networks. Their emulation accuracy and calibration performance on three real problems of Department of Energy interest is considered. On these problems, the Gaussian process emulator tends to be more accurate with narrower, but still well-calibrated uncertainty estimates. The neural network emulator is accurate, but tends to have large uncertainty on its predictions. Finally, as a result, calibration with the Gaussian process emulator produces more constrained posteriors that still perform well in prediction.

97 MATHEMATICS AND COMPUTING↗

The lifetime of ice particles in the solar system

The behavior of spherical water-ice particles in interplanetary space is investigated theoretically using an equation which gives the sublimation rate as a function of the solar energy distribution at 1 AU, solar distance, the temperature of the ice-particle surface, the temperature of the surrounding background, particle absorption efficiency, particle radius, and the complex refractive index of ice. The variation of sublimation rate with particle size is examined, and it is shown that the sublimation rate depends critically on the absorptivity of ice in the near-IR spectral region. Two different spectral-absorption curves are plotted on the basis of two different choices for the complex refractive index of ice, and an attempt is made to decide experimentally which energy distribution is more appropriate. Results of direct laboratory measurements of the sublimation rate of smooth spherical ice particles are compared with the theoretical results obtained from the cited sublimation-rate equation. Agreement between theory and experiment is found to be excellent when the complex refractive indices given by Bertie et al., (1969) are employed. Sublimation rates and particle lifetimes are then computed as a function of particle radius for heliocentric distances of 0.5, 0.75, and 1.0 AU.

Patashnick, H.↗

Influence of Ga doping on magnetic properties, magnetocaloric effect, and electronic structure of pseudobinary GdZn 1 - x Ga x ( x = 0 - 0.1 )

Here, we explore the impact of introducing IIIA-group element Ga in place of IIB-group element Zn in binary intermetallic GdZn on its magnetic and magnetocaloric properties, as well as explicate the modified electronic band structure of the compound. The magnetic transition temperature of the compound decreases with the increase of Ga concentration in GdZn 1-x ⁢Ga x ⁢(x = 0–0.1) while the crystal structure (CsCl prototype) and lattice parameters remain unchanged. Our detailed analysis of magnetization and magnetocaloric data conclusively proves that long-ranged magnetic ordering exists in the sample, despite the magnetic interaction considerably weakening with the increase of Ga. The experimental data are rationalized using both a theoretical machine-learning model and first-principle density functional theory. The electronic band structure of GdZn is manifested with some unusual complex features which gradually diminish with Ga doping and the conventional sinusoidal feature of Ruderman-Kittel-Kasuya-Yosida-type interactions also disappears. A mean-field theory model is developed and can successfully describe the overall magnetocaloric behavior of the GdZn 1-x⁢ Ga x series of samples.

36 MATERIALS SCIENCE↗

Fifth-degree elastic energy for predictive continuum stress–strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

36 MATERIALS SCIENCE↗

Fifth-degree elastic energy for predictive continuum stress-strain relations and elastic instabilities under large strain and complex loading in silicon

Materials under complex loading develop large strains and often phase transformation via an elastic instability, as observed in both simple and complex systems. Here, we represent a material (exemplified for Si I) under large Lagrangian strains within a continuum description by a 5th-order elastic energy found by minimizing error relative to density functional theory (DFT) results. The Cauchy stress—Lagrangian strain curves for arbitrary complex loadings are in excellent correspondence with DFT results, including the elastic instability driving the Si I → II phase transformation (PT) and the shear instabilities. PT conditions for Si I → II under action of cubic axial stresses are linear in Cauchy stresses in agreement with DFT predictions. Such continuum elastic energy permits study of elastic instabilities and orientational dependence leading to different PTs, slip, twinning, or fracture, providing a fundamental basis for continuum physics simulations of crystal behavior under extreme loading.

Levitas, Valery↗