Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density Functionals”

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 55 records · Page 3

Bond Dissociation Energies of the Tungsten Fluorides and Their Singly-Charged Ions: A Density Functional Survey

The dissociation of WF6 and the related singly-charged cations and anions into the lower fluorides and fluorine atoms has been investigated theoretically using density functional theory (B3LYP) and relativistic effective core potentials, with estimates of spin-orbit effects included using a simple model. The inclusion of spin-orbit is essential for a correct description of the thermochemistry. The total atomization energy of the neutral and anionic WF6 is reproduced to within 25 kcal/mol, but comparison of individual bond dissociation energies with available experimental data shows discrepancies of up to 10 kcal/mol. The results are nevertheless useful to help resolve discrepancies in experimental data and provide estimates of missing data.

Dyall, Kenneth G.↗

Eulerian Mapping Closure Approach for Probability Density Function of Concentration in Shear Flows

The Eulerian mapping closure approach is developed for uncertainty propagation in computational fluid mechanics. The approach is used to study the Probability Density Function (PDF) for the concentration of species advected by a random shear flow. An analytical argument shows that fluctuation of the concentration field at one point in space is non-Gaussian and exhibits stretched exponential form. An Eulerian mapping approach provides an appropriate approximation to both convection and diffusion terms and leads to a closed mapping equation. The results obtained describe the evolution of the initial Gaussian field, which is in agreement with direct numerical simulations.

He, Guowei↗

Modal density function and number of propagating modes in ducts

The question of the number of propagating modes within a small range of mode cut off ratio was raised. The population density of modes were shown to be greatest near cut off and least for the well propagating modes. It was shown that modes of nearly the same cut off ratio behave nearly the same in a sound absorbing duct as well as in the way they propagate to the far. Handling all of the propagating modes individually, they can be grouped into several cut off ratio ranges. It is important to know the modal density function to estimate acoustic power distribution.

Rice, E. J.↗

Linear Scaling Density Functional Calculations with Gaussian Orbitals

Recent advances in linear scaling algorithms that circumvent the computational bottlenecks of large-scale electronic structure simulations make it possible to carry out density functional calculations with Gaussian orbitals on molecules containing more than 1000 atoms and 15000 basis functions using current workstations and personal computers. This paper discusses the recent theoretical developments that have led to these advances and demonstrates in a series of benchmark calculations the present capabilities of state-of-the-art computational quantum chemistry programs for the prediction of molecular structure and properties.

Scuseria, Gustavo E.↗

Investigation of estimators of probability density functions

Four research projects are summarized which include: (1) the generation of random numbers on the IBM 360/44, (2) statistical tests used to check out random number generators, (3) Specht density estimators, and (4) use of estimators of probability density functions in analyzing large amounts of data.

Speed, F. M.↗

Time Dependent Density Functional Theory Calculations of Large Compact PAH Cations: Implications for the Diffuse Interstellar Bands

We investigate the electronic absorption spectra of several maximally pericondensed polycyclic aromatic hydrocarbon radical cations with time dependent density functional theory calculations. We find interesting trends in the vertical excitation energies and oscillator strengths for this series containing pyrene through circumcoronene, the largest species containing more than 50 carbon atoms. We discuss the implications of these new results for the size and structure distribution of the diffuse interstellar band carriers.

Weisman, Jennifer L.↗

Estimated probability density functions for the times between flashes in the storms of 12 September 1975, 26 August 1975, and 13 July 1976

A report is given to supplement the progress report of June 17, 1977. In that progress report gamma, lognormal, and Rayleigh probability density functions were fitted to the times between lightning flashes in the storms of 9/12/75, 8/26/75, and 7/13/76 by the maximum likelihood method. The goodness of fit is checked by the Kolmogoroff-Smirnoff test. Plots of the estimated densities along with normalized histograms are included to provide a visual check on the goodness of fit. The lognormal densities are the most peaked and have the highest tails. This results in the best fit to the normalized histogram in most cases. The Rayleigh densities have too broad and rounded peaks to give good fits. In addition, they have the lowest tails. The gamma densities fall inbetween and give the best fit in a few cases.

Tretter, S. A.↗

Representation of Probability Density Functions from Orbit Determination using the Particle Filter

Statistical orbit determination enables us to obtain estimates of the state and the statistical information of its region of uncertainty. In order to obtain an accurate representation of the probability density function (PDF) that incorporates higher order statistical information, we propose the use of nonlinear estimation methods such as the Particle Filter. The Particle Filter (PF) is capable of providing a PDF representation of the state estimates whose accuracy is dependent on the number of particles or samples used. For this method to be applicable to real case scenarios, we need a way of accurately representing the PDF in a compressed manner with little information loss. Hence we propose using the Independent Component Analysis (ICA) as a non-Gaussian dimensional reduction method that is capable of maintaining higher order statistical information obtained using the PF. Methods such as the Principal Component Analysis (PCA) are based on utilizing up to second order statistics, hence will not suffice in maintaining maximum information content. Both the PCA and the ICA are applied to two scenarios that involve a highly eccentric orbit with a lower apriori uncertainty covariance and a less eccentric orbit with a higher a priori uncertainty covariance, to illustrate the capability of the ICA in relation to the PCA.

Mashiku, Alinda K.↗

Carrier Modulation Via Waveform Probability Density Function

Beyond the classic modes of carrier modulation by varying amplitude (AM), phase (PM), or frequency (FM), we extend the modulation domain of an analog carrier signal to include a class of general modulations which are distinguished by their probability density function histogram. Separate waveform states are easily created by varying the pdf of the transmitted waveform. Individual waveform states are assignable as proxies for digital ONEs or ZEROs. At the receiver, these states are easily detected by accumulating sampled waveform statistics and performing periodic pattern matching, correlation, or statistical filtering. No fundamental natural laws are broken in the detection process. We show how a typical modulation scheme would work in the digital domain and suggest how to build an analog version. We propose that clever variations of the modulating waveform (and thus the histogram) can provide simple steganographic encoding.

Williams, Glenn L.↗

Carrier Modulation Via Waveform Probability Density Function

Beyond the classic modes of carrier modulation by varying amplitude (AM), phase (PM), or frequency (FM), we extend the modulation domain of an analog carrier signal to include a class of general modulations which are distinguished by their probability density function histogram. Separate waveform states are easily created by varying the pdf of the transmitted waveform. Individual waveform states are assignable as proxies for digital one's or zero's. At the receiver, these states are easily detected by accumulating sampled waveform statistics and performing periodic pattern matching, correlation, or statistical filtering. No fundamental physical laws are broken in the detection process. We show how a typical modulation scheme would work in the digital domain and suggest how to build an analog version. We propose that clever variations of the modulating waveform (and thus the histogram) can provide simple steganographic encoding.

Williams, Glenn L.↗

Investigations of turbulent scalar fields using probability density function approach

Scalar fields undergoing random advection have attracted much attention from researchers in both the theoretical and practical sectors. Research interest spans from the study of the small scale structures of turbulent scalar fields to the modeling and simulations of turbulent reacting flows. The probability density function (PDF) method is an effective tool in the study of turbulent scalar fields, especially for those which involve chemical reactions. It has been argued that a one-point, joint PDF approach is the one to choose from among many simulation and closure methods for turbulent combustion and chemically reacting flows based on its practical feasibility in the foreseeable future for multiple reactants. Instead of the multi-point PDF, the joint PDF of a scalar and its gradient which represents the roles of both scalar and scalar diffusion is introduced. A proper closure model for the molecular diffusion term in the PDF equation is investigated. Another direction in this research is to study the mapping closure method that has been recently proposed to deal with the PDF's in turbulent fields. This method seems to have captured the physics correctly when applied to diffusion problems. However, if the turbulent stretching is included, the amplitude mapping has to be supplemented by either adjusting the parameters representing turbulent stretching at each time step or by introducing the coordinate mapping. This technique is still under development and seems to be quite promising. The final objective of this project is to understand some fundamental properties of the turbulent scalar fields and to develop practical numerical schemes that are capable of handling turbulent reacting flows.

Gao, Feng↗

Steady-state probability density function of the phase error for a DPLL with an integrate-and-dump device

The steady-state behavior of a particular type of digital phase-locked loop (DPLL) with an integrate-and-dump circuit following the phase detector is characterized in terms of the probability density function (pdf) of the phase error in the loop. Although the loop is entirely digital from an implementation standpoint, it operates at two extremely different sampling rates. In particular, the combination of a phase detector and an integrate-and-dump circuit operates at a very high rate whereas the loop update rate is very slow by comparison. Because of this dichotomy, the loop can be analyzed by hybrid analog/digital (s/z domain) techniques. The loop is modeled in such a general fashion that previous analyses of the Real-Time Combiner (RTC), Subcarrier Demodulator Assembly (SDA), and Symbol Synchronization Assembly (SSA) fall out as special cases.

Simon, M.↗

Electronic Absorption Spectra of Neutral Perylene (C20H12), Terrylene (C30H16), and Quaterrylene (C40H20) and their Positive and Negative Ions: Ne Matrix-Isolation Spectroscopy and Time Dependent Density Functional Theory Calculations

We present a full experimental and theoretical study of an interesting series of polycyclic aromatic hydrocarbons, the oligorylenes. The absorption spectra of perylene, terrylene and quaterrylene in neutral, cationic and anionic charge states are obtained by matrix-isolation spectroscopy in Ne. The experimental spectra are dominated by a bright state that red shifts with growing molecular size. Excitation energies and state symmetry assignments are supported by calculations using time dependent density functional theory methods. These calculations also provide new insight into the observed trends in oscillator strength and excitation energy for the bright states: the oscillator strength per unit mass of carbon increases along the series.

Halasinski, Thomas M.↗

Constrained Kalman Filtering Via Density Function Truncation for Turbofan Engine Health Estimation

Kalman filters are often used to estimate the state variables of a dynamic system. However, in the application of Kalman filters some known signal information is often either ignored or dealt with heuristically. For instance, state variable constraints (which may be based on physical considerations) are often neglected because they do not fit easily into the structure of the Kalman filter. This paper develops an analytic method of incorporating state variable inequality constraints in the Kalman filter. The resultant filter truncates the PDF (probability density function) of the Kalman filter estimate at the known constraints and then computes the constrained filter estimate as the mean of the truncated PDF. The incorporation of state variable constraints increases the computational effort of the filter but significantly improves its estimation accuracy. The improvement is demonstrated via simulation results obtained from a turbofan engine model. The turbofan engine model contains 3 state variables, 11 measurements, and 10 component health parameters. It is also shown that the truncated Kalman filter may be a more accurate way of incorporating inequality constraints than other constrained filters (e.g., the projection approach to constrained filtering).

Simon, Dan↗

Does Chirality Influence the Stability of Amino Acid –Cu Complexes in the Salt-Induced Peptide Formation Reaction? Insights from Density Functional Theory Calculations

The polymerization of amino acids into peptides and ultimately proteins is critical to evolution of life. Ribosomally-synthesized proteins are homochiral, composed entirely of L-amino acids. Abiotic syntheses of amino acids, however, produce racemic mixtures(equal amounts of L and D enantiomers), meaning that the pool of monomers available for abiotic protein synthesis would be racemic or nearly so. There are no prebiotically plausible mechanisms known to yield homochiral peptides from racemic pools of amino acids. While some mechanisms capable of producing or amplifying enantiomeric excesses have been observed, these excesses are small or not applicable to all biologically relevant amino acids. As a result, a prebiotically-feasible mechanism to create enantiopure pools of amino acids has not been identified. However, it is not clear that these enantio-enriched or enantiopure pools are necessary to produce homochiral peptides. Steric and electronic effects, or a combination of both, may influence how different enantiomers interact with each other, potentially leading to homochiral peptides from racemic mixtures. Here, we use density functional theory (DFT) calculations to assess the stability of homochiral and heterochiral reactive complexes produced during polymerization via the salt-induced peptide formation (SIPF) reaction. Experimental studies of the SIPF reaction have shown that it enables the formation of peptides under diverse environmental conditions, making it a plausible pathway to peptide formation on early Earth. In the SIPF reaction, NaCl acts as a condensation reagent while a divalent metal cation, primarily Cu2+, facilitates polymerization by complexing and activating amino acids forming a monochlorocuprate complex. To assess if hetero-or homochiral monochlorocuprate complexes were energetically favored, we compared the stability of LL, DD, and LD enantiomers of Cu2+–(alanine)2and Cu2+–(valine) 2 complexes. Gaussian 09 density DFT energy-minimization calculations were made for each Cu –amino acid complex in the cis and trans configuration. Models were energy minimized using the BH and HLYP/6-31++G(d,p) level of theory and free energies were compared to determine the most stable configuration.

A C Fox↗

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗

Investigation of Stacking Fault and Antiphase Boundary Segregation in Ni-Based Superalloys Using Density Functional Theory

Ni-base superalloys have a long history of use in jet turbine engines, and efforts to improve their performance in that application are ongoing. It is known that the precipitation of the Ni3Al γ' phase within the disordered FCCγ phase strengthens the overall material. However, in the high-temperature environment found inside a turbine engine during operation, creep can cause the γ' phase to transform to different, weaker phases along stacking faults, leading to a deterioration of performance. In the γ' phase one mode of creep deformation is the formation of stacking fault ribbons, which consist of intrinsic stacking faults further shearing into antiphase boundaries (APBs). It is also known that certain alloying additions exhibit segregation to stacking faults. If segregating elements could be identified which segregate to the intrinsic stacking fault, but not to the APB, the inclusion of such elements could lead to improved creep strength in these alloys. To investigate this possibility, a density functional investigation of the segregation of W, Mo and Cr to both a superlattice intrinsic stacking fault (SISF) and an APB was performed. It was found that W, Mo and Cr all exhibit segregation to the SISF. In contrast, for the APB, Cr was either energy-neutral, or segregates, depending on the presence of additional nearby Cr, or on the specific lattice site upon which it was placed, while Mo and W did not segregate. Because W and Mo segregate to the SISF but not to the APB, the inclusion of these elements could provide a degree of protection against creep-related deterioration.

Good, Brian S.↗

EUPDF: Eulerian Monte Carlo Probability Density Function Solver for Applications With Parallel Computing, Unstructured Grids, and Sprays

The success of any solution methodology used in the study of gas-turbine combustor flows depends a great deal on how well it can model the various complex and rate controlling processes associated with the spray's turbulent transport, mixing, chemical kinetics, evaporation, and spreading rates, as well as convective and radiative heat transfer and other phenomena. The phenomena to be modeled, which are controlled by these processes, often strongly interact with each other at different times and locations. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. The influence of turbulence in a diffusion flame manifests itself in several forms, ranging from the so-called wrinkled, or stretched, flamelets regime to the distributed combustion regime, depending upon how turbulence interacts with various flame scales. Conventional turbulence models have difficulty treating highly nonlinear reaction rates. A solution procedure based on the composition joint probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices (such as extinction, blowoff limits, and emissions predictions) because it can account for nonlinear chemical reaction rates without making approximations. In an attempt to advance the state-of-the-art in multidimensional numerical methods, we at the NASA Lewis Research Center extended our previous work on the PDF method to unstructured grids, parallel computing, and sprays. EUPDF, which was developed by M.S. Raju of Nyma, Inc., was designed to be massively parallel and could easily be coupled with any existing gas-phase and/or spray solvers. EUPDF can use an unstructured mesh with mixed triangular, quadrilateral, and/or tetrahedral elements. The application of the PDF method showed favorable results when applied to several supersonic-diffusion flames and spray flames. The EUPDF source code will be available with the National Combustion Code (NCC) as a complete package.

Raju, M. S.↗