Analytic Signals and Zero-locus in Multiplexing Systems
Analytic functions theory used to study conventional multiplexing systems involving double and triple modulation techniques
ANALYTIC FUNCTION↗
SEARCH · Search NASA
Results for “functional 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.
Electronic Transport Through Carbon Nanotubes: Effects of Structural Deformation and the Tube Chirality
Atomistic simulations using a combination of classical force field and Density-Functional-Theory (DFT) show that carbon atoms remain essentially sp2 coordinated in either bent tubes or tubes pushed by an atomically sharp AFM tip. Subsequent Green's-function-based transport calculations reveal that for armchair tubes there is no significant drop in conductance, while for zigzag tubes the conductance can drop by several orders of magnitude in AFM-pushed tubes. The effect can be attributed to simple stretching of the tube under tip deformation, which opens up an energy gap at the Fermi surface.
Maiti, Amitesh↗
Heats of Segregation of BCC Metals Using Ab Initio and Quantum Approximate Methods
Many multicomponent alloys exhibit surface segregation, in which the composition at or near a surface may be substantially different from that of the bulk. A number of phenomenological explanations for this tendency have been suggested, involving, among other things, differences among the components' surface energies, molar volumes, and heats of solution. From a theoretical standpoint, the complexity of the problem has precluded a simple, unified explanation, thus preventing the development of computational tools that would enable the identification of the driving mechanisms for segregation. In that context, we investigate the problem of surface segregation in a variety of bcc metal alloys by computing dilute-limit heats of segregation using both the quantum-approximate energy method of Bozzolo, Ferrante and Smith (BFS), and all-electron density functional theory. In addition, the composition dependence of the heats of segregation is investigated using a BFS-based Monte Carlo procedure, and, for selected cases of interest, density functional calculations. Results are discussed in the context of a simple picture that describes segregation behavior as the result of a competition between size mismatch and alloying effects
Good, Brian↗
Experimental and Theoretical Study of Thermodynamics of the Reaction of Titania and Water at High Temperatures
The transpiration method was used to determine the volatility of titanium dioxide (TiO2) in water vapor-containing environments at temperatures between 1473 and 1673 K. Water contents ranged from 0 to 76 mole % in oxygen or argon carrier gases for 20 to 250 hr exposure times. Results indicate that oxygen is not a key contributor to volatilization and the primary reaction for volatilization in this temperature range is: TiO2(s) + H2O(g) = TiO(OH)2(g). Data were analyzed with both the second and third law methods to extract an enthalpy and entropy of formation. The geometry and vibrational frequencies of TiO(OH)2(g) were computed using B3LYP density functional theory, and the enthalpy of formation was computed using the coupled-cluster singles and doubles method with a perturbative correction for connected triple substitutions [CCSD(T)]. Thermal functions are calculated using both a structure with bent and linear hydroxyl groups. Calculated second and third heats show closer agreement with the linear hydroxyl group, suggesting more experimental and computational spectroscopic and structural work is needed on this system.
equilibrium↗
A Fast Method of Deriving the Kirchhoff Formula for Moving Surfaces
The Kirchhoff formula for a moving surface is very useful in many wave propagation problems, particularly in the prediction of noise from rotating machinery. Several publications in the last two decades have presented derivations of the Kirchhoff formula for moving surfaces in both time and frequency domains. Here we present a method originally developed by Farassat and Myers in time domain that is both simple and direct. It is based on generalized function theory and the useful concept of imbedding the problem in the unbounded three-dimensional space. We derive an inhomogeneous wave equation with the source terms that involve Dirac delta functions with their supports on the moving data surface. This wave equation is then solved using the simple free space Green's function of the wave equation resulting in the Kirchhoff formula. The algebraic manipulations are minimal and simple. We do not need the Green's theorem in four dimensions and there is no ambiguity in the interpretation of any terms in the final formulas. Furthermore, this method also gives the simplest derivation of the classical Kirchhoff formula which has a fairly lengthy derivation in physics and applied mathematics books. The Farassat-Myers method can be used easily in frequency domain.
Farassat, F.↗
The 3-D inelastic analysis methods for hot section components: Brief description
Advanced 3-D inelastic structural/stress analysis methods and solution strategies for more accurate yet more cost-effective analysis of components subjected to severe thermal gradients and loads in the presence of mechanical loads, with steep stress and strain gradients are being developed. Anisotropy, time and temperature dependent plasticity and creep effects are also addressed. The approach is to develop four different theories, one linear and three higher order theories (polynomial function, special function, general function). The theories are progressively more complex from linear to general function in order to provide streamlined analysis capability with increasing accuracy for each hot section component and for different parts of the same component according to the severity of the local stress, strain and temperature gradients associated with hot spots, cooling holes and surface coating cracks. To further enhance the computational effectiveness, the higher order theories will have embedded singularities (cooling passages, for example) in the generic modeling region. Each of the four theories consists of three formulation models derivable from independent theoretical formulations. These formulation models are based on: (1) mechanics of materials; (2) special finite elements; and (3) an advanced formulation to be recommended by the contractor.
Chamis, C. C.↗
Modeling Ring-Opening Polymerization with Machine-Learned Committors
Many important industrial and biological processes depend on ring-opening polymerization (ROP), including large-scale production polymers, such as nylon 6, and a variety of biodegradable polymers and high-performance resins. Large-scale simulations (>200K atoms) of ROP were performed using a deep neural network-based method that captured the effect of the local chemical environment of the reactive sites using training data from quantum chemical methods. The utility of this method was demonstrated by modeling the cationic ROP of bis-benzoxazine, a promising resin for use in advanced high-temperature application composites (>200°C). Density functional theory calculations were used to map out the various ring opening and polymerization reaction paths that occur during the thermal processing of bis-benzoxazine. The probabilities of these reactions, as a function of local atomic configurations, were used to construct committor functions using a neural network approach. The learned committor functions were then used to determine when a particular reaction occurred within the REACTER framework for reactive molecular dynamics simulations. The experimentally observed exotherm that occurs during bis-benzoxazine ROP was captured in an effective manner in the simulations by imposing a local heating of the system that corresponds to the enthalpy of reaction for each reaction. The dependence of the final polymer morphology on the ROP reaction temperature was analyzed at the molecular scale. The results presented in this work indicate that machine-learned committor functions are a promising approach for incorporating high-fidelity reactivity criteria that capture the effect of local chemical environments into large-scale polymerization simulations.
machine learning, polymer simulations, molecular d↗
High-Performance Computing Optimization for Aladyn – Adaptive Neural Network Molecular Dynamics Mini-Application
This report provides a description and performance evaluation of the optimization techniques for high performance computing (HPC) implementation of the open source Computational Materials mini-application Aladyn (https://github.com/nasa/aladyn). Aladyn is a basic molecular dynamics code written in FORTRAN 2003, which is designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory (DFT) method. While achieving orders of magnitude faster computational performance than DFT, the ANN-based approach was still very computationally demanding compared to the conventional approach of using empirically fitted energy functions. After its initial development, Aladyn was evaluated and optimized by experts at the NASA Advanced Supercomputing (NAS) division to exploit modern supercomputer architectures. The code has been optimized for execution on multicore central processing units (CPUs), including Intel® Skylake microarchitecture, and on graphic accelerators, such as Nvidia® V100 graphic processing units (GPUs), using Open Multi-Processing (OpenMP) and Open Accelerators (OpenACC) programming interfaces. The optimization achieved a speedup of 4.7 times the baseline version on CPU performance and an additional 2.4 times on CPU+GPU performance. Atomistic computer simulations are a fundamental tool in materials research to model material properties form physics-based first principles. Atomic interaction, governed by Quantum Mechanics (QM) require sophisticated and highly computationally demanding mathematical models to calculate [1]. Classical methods use approximate functional forms, empirically fitted through a set of variable parameters to emulate atomic energies as direct functions of atomic coordinates [2]. While empirical potentials are computationally much simpler, allowing simulations of large-scale systems of up to a trillion (1012) atoms [3], they are substantially less accurate compared to quantum calculations and applicable only to very specific atomic configurations or predefined crystallographic phases. A recently suggested approach is to use heuristic machine learning methods [4], such as those based on Adaptive Neural Networks (ANNs) to predict atomic energies, after being trained on a sufficiently large database of QM-calculated structures [5,6]. This approach reduces significantly the computational complexity, allowing for simulations of orders of magnitude larger systems compared to QM-based methods without compromising accuracy. Still, compared to classical methods using empirical energy functions, ANN methods remain two- to three orders of magnitude more computationally demanding. Hence, the computational cost of simulations, together with the need for extensive training of ANNs, still makes the practical implementation of ANN-based methods quite challenging. The purpose of the Aladyn mini-application software [7], available as open source at https://github.com/nasa/aladyn, is to be a testbed for exploring possible optimization strategies to develop highly scalable parallel algorithms for ANN-based atomistic simulations. Aladyn is aimed at utilizing the architecture of the high-end modern highperformance computing (HPC) hardware based on multicore central processing units (CPUs) equipped with graphic processing unit (GPU) accelerators. Specifically, the goal is to optimize the performance on a single HPC compute node, before implementing scaling to multi-node parallelization using message passing interface (MPI). At the same time, the open source code of Aladyn can serve as a training model for students and professors in academia.
Yamakov, Vesselin I.↗
Impact of Functionally Graded Cylinders: Theory
This final report summarizes the work funded under the Grant NAG3-2411 during the 04/05/2000-04/04/2001 period. The objective of this one-year project was to generalize the theoretical framework of the two-dimensional higher-order theory for the analysis of cylindrical functionally graded materials/structural components employed in advanced aircraft engines developed under past NASA Glenn funding. The completed generalization significantly broadens the theory's range of applicability through the incorporation of dynamic impact loading capability into its framework. Thus, it makes possible the assessment of the effect of damage due to fuel impurities, or the presence of submicron-level debris, on the life of functionally graded structural components. Applications involving advanced turbine blades and structural components for the reusable-launch vehicle (RLV) currently under development will benefit from the completed work. The theory's predictive capability is demonstrated through a numerical simulation of a one-dimensional wave propagation set up by an impulse load in a layered half-plane. Full benefit of the completed generalization of the higher-order theory described in this report will be realized upon the development of a related computer code.
Aboudi, Jacob↗
The Kirchhoff Formulas for Moving Surfaces in Aeroacoustics - The Subsonic and Supersonic Cases
One of the active areas of computational aeroacoustics is the application of the Kirchhoff formulas to the problems of the rotating machinery noise predictions. The original Kirchhoff formula was derived for a stationary surface. In 1988, Farassat and Myers derived a Kirchhoff Formula obtained originally by Morgans using modem mathematics. These authors gave a formula particularly useful for applications in aeroacoustics. This formula is for a surface moving at subsonic speed. Later in 1995 these authors derived the Kirchhoff formula for a super-sonically moving surface. This technical memorandum presents the viewgraphs of a day long workshop by the author on the derivation of the Kirchhoff formulas. All necessary background mathematics such as differential geometry and multidimensional generalized function theory are discussed in these viewgraphs. Abstraction is kept at minimum level here. These viewgraphs are also suitable for understanding the derivation and obtaining the solutions of the Ffowcs Williams-Hawkings equation. In the first part of this memorandum, some introductory remarks are made on generalized functions, the derivation of the Kirchhoff formulas and the development and validation of Kirchhoff codes. Separate lists of references by Lyrintzis, Long, Strawn and their co-workers are given in this memorandum. This publication is aimed at graduate students, physicists and engineers who are in need of the understanding and applications of the Kirchhoff formulas in acoustics and electromagnetics.
Farassat, F.↗
On the application of Mathieu functions in the theory of subsonic compressible flow past oscillating airfoils
Explore the source record for details and available documents.
THEORIES - FLOW↗
Machine Learning the COSMO Model for Predicting Thermodynamics of Electrolyte Mixtures
Bottom-up design of electrolyte mixtures for battery systems requires predicting macro thermodynamic properties from molecular constituents. For instance, molten salt electrolyte batteries require conditions far above room temperature to operate. Therefore, discovering mixtures with increasingly lower eutectic melting points is desirable. A model that can approximate chemical activity is a valuable tool to search through the vast compositional design space. Machine learning can predict properties of materials such as vibrational free energies, electronic energy gaps, and thermal conductivities. Moreover, they can learn physical models such as interatomic potentials. The COSMO-SAC model uses theory and empirical parameterization to predict liquid-vapor and liquid-solid properties using first-principles calculations. However, obtaining activity coefficients required for parameterizing the COSMO-SAC model is costly and limited to a select chemical space. In this work, we explored if machine learning methods could improve the COSMO-SAC model and bridge density functional theory calculations to liquid phase thermodynamic properties. Our data-driven approach uses existing databases for sigma-profiles of organic solvents and reconciles their methodological differences via ensemble averaging. First, an optimal machine learning model is constructed for each dataset. Our machine learning algorithms use the sigma-profile as an input feature to predict binary mixtures' activity coefficients using multi-output regression. Each dataset uses different choices of functionals, methods, and basis sets. Therefore, our ensemble model attempts to predict corrected activity coefficients given the combination of all the model outputs. The activity coefficients used for training are generated using the COSMO-SAC model. This approach enables the extraction of meaningful information from the existing datasets to improve the COSMO-SAC model for obtaining thermodynamic properties of electrolyte mixtures. With the liquid phase activities, we can identify electrolyte mixtures that meet desired phase equilibria conditions.
Thermodynamics↗
AladynPi – Adaptive Neural Network Molecular Dynamics Simulation Code with Physically Informed Potential: Computational Materials Mini-Application
This report provides an overview and description of commands used in the Computational Materials mini-application, AladynPi. AladynPi is an extension of a previously released mini-application, Aladyn (https://github.com/nasa/aladyn; Yamakov, V.I., and Glaessgen, E.H., NASA/TM-2018-220104). Aladyn and AladynPi are basic molecular dynamics codes written in FORTRAN 2003, which are designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory method. An input for the ANN is a set of structure coefficients, characterizing the local atomic environment of each atom, for which the atomic energy is obtained in the ANN inference process. In Aladyn, the ANN gives directly the energy of interatomic interactions. In AladynPi, the ANN gives optimized parameters for a predefined empirical function, known as bond-order-potential (BOP). The parameterized BOP function is then used to calculate the energy. AladynPi code is being released to serve as a training testbed for students and professors in academia to explore possible optimization algorithms for parallel computing on multicore central processing unit (CPU) computers or computers utilizing manycore architectures based on graphic processing units (GPUs). The effort is supported by the High Performance Computing incubator (HPCi) project at NASA Langley Research Center.
Yamakov, Vesselin I.↗
Effective Forces Between Colloidal Particles
Colloidal suspensions have proven to be excellent model systems for the study of condensed matter and its phase behavior. Many of the properties of colloidal suspensions can be investigated with a systematic variation of the characteristics of the systems and, in addition, the energy, length and time scales associated with them allow for experimental probing of otherwise inaccessible regimes. The latter property also makes colloidal systems vulnerable to external influences such as gravity. Experiments performed in micro-ravity by Chaikin and Russell have been invaluable in extracting the true behavior of the systems without an external field. Weitz and Pusey intend to use mixtures of colloidal particles with additives such as polymers to induce aggregation and form weak, tenuous, highly disordered fractal structures that would be stable in the absence of gravitational forces. When dispersed in a polarizable medium, colloidal particles can ionize, emitting counterions into the solution. The standard interaction potential in these charged colloidal suspensions was first obtained by Derjaguin, Landau, Verwey and Overbeek. The DLVO potential is obtained in the mean-field linearized Poisson-Boltzmann approximation and thus has limited applicability. For more precise calculations, we have used ab initio density functional theory. In our model, colloidal particles are charged hard spheres, the counterions are described by a continuum density field and the solvent is treated as a homogeneous medium with a specified dielectric constant. We calculate the effective forces between charged colloidal particles by integrating over the solvent and counterion degrees of freedom, taking into account the direct interactions between the particles as well as particle-counterion, counterion-counterion Coulomb, counterion entropic and correlation contributions. We obtain the effective interaction potential between charged colloidal particles in different configurations. We evaluate two- and three-body forces in the bulk as well as study the influence of soft walls. We qualitatively explain the effects of the walls on the forces and demonstrate that many-body effects are negligible in our system. With adjustments in the parameters, the DLVO pair-potential can describe the results quantitatively. Besides electrostatic interactions, entropic depletion effects that arise from (hard-core) exclusion play an important role in determining the behavior of multi-component colloidal suspensions. A standard theory for depletion forces is due to Asakura and Oosawa and is based on the ideal gas approximation. To go beyond this approximation, we have studied entropic forces in molecular dynamics simulations of systems of hard spheres (the effects of the solvent have been ignored). The effective depletion forces for these systems can be found either from equilibrium distribution functions or from direct momentum transfer calculations. Our results obtained by either method show qualitative differences from the Asakura-Oosawa forces, indicating a longer range, higher value at contact and most importantly a more complicated structure, comprising of several maxima and minima. Our calculations include the determination of effective forces between two spheres, a hard sphere and a wall, and the behavior of a hard sphere near a step-edge and a corner. We also demonstrate that such entropic forces do not necessarily satisfy pairwise additivity.
Tehver, Riina↗
Aerodynamics via acoustics - Application of acoustic formulas for aerodynamic calculations
Prediction of aerodynamic loads on bodies in arbitrary motion is considered from an acoustic point of view, i.e., in a frame of reference fixed in the undisturbed medium. An inhomogeneous wave equation which governs the disturbance pressure is constructed and solved formally using generalized function theory. When the observer is located on the moving body surface there results a singular linear integral equation for surface pressure. Two different methods for obtaining such equations are discussed. Both steady and unsteady aerodynamic calculations are considered. Two examples are presented, the more important being an application to propeller aerodynamics. Of particular interest for numerical applications is the analytical behavior of the kernel functions in the various integral equations.
Farassat, F.↗
Aerodynamics Via Acoustics: Application of Acoustic Formulas for Aerodynamic Calculations
Prediction of aerodynamic loads on bodies in arbitrary motion is considered from an acoustic point of view, i.e., in a frame of reference fixed in the undisturbed medium. An inhomogeneous wave equation which governs the disturbance pressure is constructed and solved formally using generalized function theory. When the observer is located on the moving body surface there results a singular linear integral equation for surface pressure. Two different methods for obtaining such equations are discussed. Both steady and unsteady aerodynamic calculations are considered. Two examples are presented, the more important being an application to propeller aerodynamics. Of particular interest for numerical applications is the analytical behavior of the kernel functions in the various integral equations.
Farassat, F.↗
Formation of Carbon Nanotube Based Gears: Quantum Chemistry and Molecular Mechanics Study of the Electrophilic Addition of o-Benzyne to Fullerenes, Graphene, and Nanotubes
Considerable progress has been made in recent years in chemical functionalization of fullerene molecules. In some cases, the predominant reaction products are different from those obtained (using the same reactants) from polycyclic aromatic hydrocarbons (PAHs). One such example is the cycloaddition of o-benzyne to C60. It is well established that benzyne adds across one of the rings in naphthalene, anthracene and other PAHs forming the [2+4] cycloaddition product (benzobicyclo[2.2.2.]-octatriene with naphthalene and triptycene with anthracene). However, Hoke et al demonstrated that the only reaction path for o-benzyne with C60 leads to the [2+2] cycloaddition product in which benzyne adds across one of the interpentagonal bonds (forming a cyclobutene ring in the process). Either reaction product results in a loss of aromaticity and distortion of the PAH or fullerene substrate, and in a loss of strain in the benzyne. It is not clear, however, why different products are preferred in these cases. In the current paper, we consider the stability of benzyne-nanotube adducts and the ability of Brenner's potential energy model to describe the structure and stability of these adducts. The Brenner potential has been widely used for describing diamondoid and graphitic carbon. Recently it has also been used for molecular mechanics and molecular dynamics simulations of fullerenes and nanotubes. However, it has not been tested for the case of functionalized fullerenes (especially with highly strained geometries). We use the Brenner potential for our companion nanogear simulations and believe that it should be calibrated to insure that those simulations are physically reasonable. In the present work, Density Functional theory (DFT) calculations are used to determine the preferred geometric structures and energetics for this calibration. The DFT method is a kind of ab initio quantum chemistry method for determining the electronic structure of molecules. For a given basis set expansion, it is comparable in accuracy to the MP2 method (better than Hartree Fock, but less accurate than more extensive electron correlation methods such as MP4 or CCSD). However, for systems with large numbers of basis functions it more efficient than any other methods that include electron correlation effects. In this presentation we show the results of DFT calculations for the reaction of benzyne with naphthalene, C60, and nanotube models. We compare energies for [2+2] and [2+4] cycloaddition products. The preferred products for the naphthalene and C60 reactions have been determined by experiment and, thus, these cases serve as a validation of our quantum chemical approach. We also compare the DFT and Brenner potential results. Finally we can predict the likelihood of reaction between benzyne and nanotubes.
Jaffe, Richard↗
Mid-Field Sonic Boom Extrapolation Methodology
In the design cycle of low-boom airplanes, sonic boom prediction must be accurate and efficient. The classical linear method, Whitham's F-function theory, has been widely applied to predict sonic boom signatures. However, linear theory fails to capture the nonlinear effects created by large civil transport. Computational fluid dynamics (CFD) has been used successfully to predict sonic boom signals at the near and mid fields. Nevertheless, it is computationally expansive in airplane design runs. In the present study, the method of characteristics is used to predict sonic boom signals in an efficient fashion. The governing equations are the axisymmetric Euler's equations with constant enthalpy. Since the method solves Euler's equations, it captures more nonlinear effects than the classical Whitham's F-function technique. Furthermore, the method of characteristics is an efficient marching scheme for initial value problems. In this study, we will first review the current CFD extrapolation technique and the work previously done in sonic boom extrapolation. Then, we will introduce the governing equations and the method of characteristics. Finally, we will show that the present method yields the same accurate results as previous CFD techniques, but with higher efficiency.
Cheung, Samson↗