Engineering Papers⌕ Search

Engineering topics

Bauschlicher, Charles W., Jr.

Publications and source records attributed to Bauschlicher, Charles W., Jr..

At least 19 records

Multiscale Modeling of Carbon/Phenolic Composite Thermal Protection Materials: Atomistic to Effective Properties

Next generation ablative thermal protection systems are expected to consist of 3D woven composite architectures. It is well known that composites can be tailored to achieve desired mechanical and thermal properties in various directions and thus can be made fit-for-purpose if the proper combination of constituent materials and microstructures can be realized. In the present work, the first, multiscale, atomistically-informed, computational analysis of mechanical and thermal properties of a present day - Carbon/Phenolic composite Thermal Protection System (TPS) material is conducted. Model results are compared to measured in-plane and out-of-plane mechanical and thermal properties to validate the computational approach. Results indicate that given sufficient microstructural fidelity, along with lowerscale, constituent properties derived from molecular dynamics simulations, accurate composite level (effective) thermo-elastic properties can be obtained. This suggests that next generation TPS properties can be accurately estimated via atomistically informed multiscale analysis.

Materials Science↗

Physics-Based Modeling of Meteor Entry and Breakup

A new research effort at NASA Ames Research Center has been initiated in Planetary Defense, which integrates the disciplines of planetary science, atmospheric entry physics, and physics-based risk assessment. This paper describes work within the new program and is focused on meteor entry and breakup. Over the last six decades significant effort was expended in the US and in Europe to understand meteor entry including ablation, fragmentation and airburst (if any) for various types of meteors ranging from stony to iron spectral types. These efforts have produced primarily empirical mathematical models based on observations. Weaknesses of these models, apart from their empiricism, are reliance on idealized shapes (spheres, cylinders, etc.) and simplified models for thermal response of meteoritic materials to aerodynamic and radiative heating. Furthermore, the fragmentation and energy release of meteors (airburst) is poorly understood. On the other hand, flight of human-made atmospheric entry capsules is well understood. The capsules and their requisite heat shields are designed and margined to survive entry. However, the highest speed Earth entry for capsules is 13 kms (Stardust). Furthermore, Earth entry capsules have never exceeded diameters of 5 m, nor have their peak aerothermal environments exceeded 0.3 atm and 1 kWcm2. The aims of the current work are: (i) to define the aerothermal environments for objects with entry velocities from 13 to 20 kms; (ii) to explore various hypotheses of fragmentation and airburst of stony meteors in the near term; (iii) to explore the possibility of performing relevant ground-based tests to verify candidate hypotheses; and (iv) to quantify the energy released in airbursts. The results of the new simulations will be used to anchor said risk assessment analyses.With these aims in mind, state-of-the-art entry capsule design tools are being extended for meteor entries. We describe: (i) applications of current simulation tools to spherical geometries of diameters ranging from 1 to 100 m for an entry velocity of 20 kms and stagnation pressures ranging from 1 to 100 atm; (ii) the influence of shape and departure of heating environment predictions from those for a simple spherical geometry; (iii) assessment of thermal response models for silica subject to intense radiation; and (iv) results for porosity-driven gross fragmentation of meteors, idealized as a collection of smaller objects. Lessons learned from these simulations will be used to help understand the Chelyabinsk meteor entry up to its first point of fragmentation.

Atmospheric Entry↗

Physics-Based Modeling of Meteor Entry and Breakup

A new research effort at NASA Ames Research Center has been initiated in Planetary Defense, which integrates the disciplines of planetary science, atmospheric entry physics, and physics-based risk assessment. This paper describes work within the new program and is focused on meteor entry and breakup.Over the last six decades significant effort was expended in the US and in Europe to understand meteor entry including ablation, fragmentation and airburst (if any) for various types of meteors ranging from stony to iron spectral types. These efforts have produced primarily empirical mathematical models based on observations. Weaknesses of these models, apart from their empiricism, are reliance on idealized shapes (spheres, cylinders, etc.) and simplified models for thermal response of meteoritic materials to aerodynamic and radiative heating. Furthermore, the fragmentation and energy release of meteors (airburst) is poorly understood.On the other hand, flight of human-made atmospheric entry capsules is well understood. The capsules and their requisite heatshields are designed and margined to survive entry. However, the highest speed Earth entry for capsules is 13 kms (Stardust). Furthermore, Earth entry capsules have never exceeded diameters of 5 m, nor have their peak aerothermal environments exceeded 0.3 atm and 1 kW/sq cm. The aims of the current work are: (i) to define the aerothermal environments for objects with entry velocities from 13 to 20 kms; (ii) to explore various hypotheses of fragmentation and airburst of stony meteors in the near term; (iii) to explore the possibility of performing relevant ground-based tests to verify candidate hypotheses; and (iv) to quantify the energy released in airbursts. The results of the new simulations will be used to anchor said risk assessment analyses. With these aims in mind, state-of-the-art entry capsule design tools are being extended for meteor entries. We describe: (i) applications of current simulation tools to spherical geometries of diameters ranging from 1 to 100 m for an entry velocity of 20 kms and stagnation pressures ranging from 1 to 100 atm; (ii) the influence of shape and departure of heating environment predictions from those for a simple spherical geometry; (iii) assessment of thermal response models for silica subject to intense radiation; and (iv) results for porosity-driven gross fragmentation of meteors, idealized as a collection of smaller objects. Lessons learned from these simulations will be used to help understand the Chelyabinsk meteor entry up to its first point of fragmentation.

Near-Earth Asteroids↗

Investigation of Structure and Transport in Li-Doped Ionic Liquid Electrolytes: [pyr14][TFSI], [pyr13][FSI], [EMIM][BF4]

Ionic liquid electrolytes have been proposed as a means of improving the safety and cycling behavior of advanced lithium batteries; however, the properties of these electrolytes under high lithium doping are poorly understood. Here, we employ both polarizable molecular dynamics simulation and experiment to investigate the structure, thermodynamics and transport of three potential electrolytes, N-methyl-N-butylpyrrolidinium bis(trifluoromethylsufonyl)imide ([pyr14][TFSI]), N- methyl-N-propylpyrrolidinium bis(fluorosufonyl)imide ([pyr13][FSI]), and 1-ethyl-3-- methylimidazolium boron tetrafluoride ([EMIM][BF4]), as a function of Li (-) salt concentration and temperature. Structurally, Li(+) is shown to be solvated by three anion neighbors in [pyr14][TFSI] and four anion neighbors in both [pyr13][FSI] and [EMIM][BF4], and at all levels of xLi we find the presence of lithium aggregates. Furthermore, the computed density, diffusion, viscosity, and ionic conductivity show excellent agreement with experimental data. While the diffusion and viscosity exhibit a systematic decrease and increase, respectively, with increasing xLi, the contribution of Li(+) to ionic conductivity increases until reaching a saturation doping level of xLi 0.10. Comparatively, the Li(+) conductivity of [pyr14][TFSI] is an order of magnitude lower than that of the other liquids, which range between 0.1-0.3 mScm. The differences in Li(+) transport are reflected in the residence times of Li(+) with the anions, which are revealed to be much larger for [pyr14][TFSI] (up to 100 ns at the highest doping levels) than in either [EMIM][BF4] or [pyr13][FSI]. Finally, we comment on the relative kinetics of Li(+) transport in each liquid and we present strong evidence for transport through anion exchange (hopping) as opposed to the net motion of Li(+) with its solvation shell (vehicular).

molecular dynamics↗

Investigation of Structure and Transport in Li-Doped Ionic Liquid Electrolytes: [pyr14][TFSI], [pyr13][FSI], [EMIM][BF4]

Ionic liquid electrolytes have been proposed as a means of improving the safety and cycling behavior of advanced lithium batteries; however, the properties of these electrolytes under high lithium doping are poorly understood. Here, we employ both polarizable molecular dynamics simulation and experiment to investigate the structure, thermodynamics and transport of three potential electrolytes, N-methyl-N-butylpyrrolidinium bis(trifluoromethylsufonyl)imide ([pyr14][TFSI]), N- methyl-N-propylpyrrolidinium bis(fluorosufonyl)imide ([pyr13][FSI]), and 1-ethyl-3-- methylimidazolium boron tetrafluoride ([EMIM][BF4]), as a function of Li (-) salt concentration and temperature. Structurally, Li(+) is shown to be solvated by three anion neighbors in [pyr14][TFSI] and four anion neighbors in both [pyr13][FSI] and [EMIM][BF4], and at all levels of xLi we find the presence of lithium aggregates. Furthermore, the computed density, diffusion, viscosity, and ionic conductivity show excellent agreement with experimental data. While the diffusion and viscosity exhibit a systematic decrease and increase, respectively, with increasing xLi, the contribution of Li(+) to ionic conductivity increases until reaching a saturation doping level of xLi 0.10. Comparatively, the Li(+) conductivity of [pyr14][TFSI] is an order of magnitude lower than that of the other liquids, which range between 0.1-0.3 mScm. The differences in Li(+) transport are reflected in the residence times of Li(+) with the anions, which are revealed to be much larger for [pyr14][TFSI] (up to 100 ns at the highest doping levels) than in either [EMIM][BF4] or [pyr13][FSI]. Finally, we comment on the relative kinetics of Li(+) transport in each liquid and we present strong evidence for transport through anion exchange (hopping) as opposed to the net motion of Li(+) with its solvation shell (vehicular).

ionic liquids↗

Multiscale Modeling of Grain Boundaries in ZrB2: Structure, Energetics, and Thermal Resistance

A combination of ab initio, atomistic and finite element methods (FEM) were used to investigate the structures, energetics and lattice thermal conductance of grain boundaries for the ultra high temperature ceramic ZrB2. Atomic models of idealized boundaries were relaxed using density functional theory. Information about bonding across the interfaces was determined from the electron localization function. The Kapitza conductance of larger scale versions of the boundary models were computed using non-equilibrium molecular dynamics. The interfacial thermal parameters together with single crystal thermal conductivities were used as parameters in microstructural computations. FEM meshes were constructed on top of microstructural images. From these computations, the effective thermal conductivity of the polycrystalline structure was determined.

Lawson, John W.↗

Lattice Thermal Conductivity of Ultra High Temperature Ceramics ZrB2 and HfB2 from Atomistic Simulations

Atomistic Green-Kubo simulations are performed to evaluate the lattice thermal conductivity for single crystals of the ultra high temperature ceramics ZrB2 and HfB2 for a range of temperatures. Recently developed interatomic potentials are used for these simulations. Heat current correlation functions show rapid oscillations which can be identified with mixed metal-Boron optical phonon modes. Agreement with available experimental data is good.

Lawson, John W.↗

Application of Density Functional Theory to Systems Containing Metal Atoms

The accuracy of density functional theory (DFT) for problems involving metal atoms is considered. The DFT results are compared with experiment as well as results obtained using the coupled cluster approach. The comparisons include geometries, frequencies, and bond energies. The systems considered include MO2, M(OH)+n, MNO+, and MCO+2. The DFT works well for frequencies and geometries, even in case with symmetry breaking; however, some examples have been found where the symmetry breaking is quite severe and the DFT methods do not work well. The calculation of bond energies is more difficult and examples of successes as well as failures of DFT will be given.

Bauschlicher, Charles W., Jr.↗

The Infrared Spectrum of Matrix Isolated Aminoacetonitrile: A Precursor to the Amino Acid Glycine

We present infrared (IR) spectral data from matrix isolation experiments and density functional theory calculations on the pre-biologically interesting molecule aminoacetonitrile, a precursor to glycine. We find that this nitrile has an unusually weak nitrile (C=N) stretch in the infrared, in contrast to expectations based on measurements and models of other nitriles under astrophysical conditions. The absence of an observable nitrile absorption feature in the infrared will make the IR search for this molecule considerably more difficult, and will raise estimates of upper limits on nitriles in interstellar and outer Solar System ices. This is also of relevance to assessing the formation routes of the amino acid glycine, since aminoacetonitrile is the putative precursor to glycine via the Strecker synthesis, the mechanism postulated to have produced the amino acids in meteorites.

Bernstein, Max P.↗

Deuterated PAHs in Space

The cosmic deuterium to hydrogen (D/H) ratio is of key importance from a cosmological and stellar evolution perspective since deuterium originates from big-bang nucleosynthesis and is destroyed by stellar thermonuclear reactions. Further, from the interstellar perspective, the galactic distribution of deuterium and the D/H ratio among various molecular species also traces interstellar chemical evolution. Over the past few decades, radio observations have enabled the study of a handful of small, deuterated interstellar species. However, the number of deuterated species detected and environments probed are limited, raising issues of selection effects that hamper generalization and applications to other environments. Infrared spectroscopy of the interstellar medium offers a distinct advantage in this regard as the extent of deuteration of entire chemical families, rather than one species, can be probed. These observations require spaceborne telescopes because the molecular vibrations involving D which produce the strongest IR bands fall in spectral regions which are obscured by terrestrial CO2 absorption. Here we report the tentative detection of the C-D stretching vibration from deuterated interstellar polycyclic aromatic hydrocarbons in the Orion nebula. Since the PAH emission features are widespread and probe many different types of cosmic environments, follow up observations of deuterated PAHs will provide fundamental, far reaching new insight and perspective into galactic and extragalactic processes.

Peeters, Els↗

High Coverages of Hydrogen on a (10,0) Carbon Nanotube

The binding energy of H to a (10,0) carbon nanotube is calculated at 24, 50, and 100% coverage. Several different bonding configurations are considered for the 50% coverage case. Using the ONIOM (our own n-layered integrated molecular orbital and molecular mechanics) approach, the average C-H bond energy for the most stable 50% coverage and for the 100% coverage are 57.3 and 38.6 kcal/mol, respectively. Considering the size of the bond energy of H2, these values suggest that it will be difficult to achieve 100% atomic H coverage on a (10,0) nanotube.

Bauschlicher, Charles W., Jr.↗

On Interpreting the Photoelectron Spectra of MgO

The (sup 2)Sigma(+) and (sup 2)Pi states of MgO(-) and the (sup 1)Sigma(+), (sup 1)Pi, and (sup 3)Pi states of MgO are studied using the averaged coupled-pair functional (ACPF) approach. The computed spectroscopic constants are in good agreement with the available experimental data. The computed Franck-Condon factors and photodetachment overlaps are compared with experiment.

Bauschlicher, Charles W., Jr.↗

Closed-Shell Polycyclic Aromatic Hydrocarbon Cations: A New Category of Interstellar Polycyclic Aromatic Hydrocarbons

Density functional theory has been employed to calculate the harmonic frequencies and intensities of a range of polycyclic aromatic hydrocarbon (PAH) cations that explore both size and electronic structure effects of the infrared spectroscopic properties of these species. The sample extends the size range of PAH species considered to more than 50 carbon atoms and includes several representatives from each of two heretofore unexplored categories of PAH cations: (1) fully benzenoid PAH cations whose carbon skeleton is composed of an odd number of carbon atoms (C(sub odd) PAHs); and (2) protonated PAH cations (HPAH+). Unlike the radical electronic structures of the PAH cations that have been the subject of previous theoretical and experimental work, the species in these two classes have a closed-shell electronic configuration. The calculated spectra of circumcoronene, C54H18 in both neutral and (radical) cationic form are also reported and compared with those of the other species. Overall, the C(sub odd) PAHs spectra are dominated by strong CC stretching modes near 1600 cm(exp -1) and display spectra that are remarkably insensitive to molecular size. The HPAH+ species evince a more complex spectrum consistent with the added contributions of aliphatic modes and their generally lower symmetry. Finally, for both classes of closed-shell cations, the intensity of the aromatic CH stretching modes is found to increase with molecular size far out of proportion with the number of CH groups, approaching a value more typical of neutral PAHs for the largest species studied.

Hudgins, Douglas M.↗

The A6Sigma+ - X6Sigma+ Transition of CrH, Einstein Coefficients and an Improved Description of the A State

The spectrum of CrH has been reinvestigated in the 9000-15000/cm region using the Fourier transform spectrometer of the National Solar Observatory. The 1-0 and 1-1 bands of the A6Sigma+ - X6Sigma+ transition have been measured and improved spectroscopic constants have been determined. A value for the 2-0 band origin has been obtained from the band head using estimated spectroscopic constants. These data provide a set of much improved equilibrium vibrational and rotational constants for the A6Sigma+ state. An accurate description of the A-X transition has been obtained using a multi-reference configuration interaction approach. The inclusion of both scalar relativity and Cr 3s3p correlation are required to obtain a good description of both states. The ab initio computed Einstein coefficients and radiative lifetimes are reported.

Bauschlicher, Charles W., Jr.↗

Is the Lamb shift chemically significant?

The contribution of the Lamb shift to the atomization energies of some prototype molecules, BF3, AlF3, and GaF3, is estimated by a perturbation procedure. It is found to be in the range of 3-5% of the one-electron scalar relativistic contribution to the atomization energy. The maximum absolute value is 0.2 kcal/mol for GaF3. These sample calculations indicate that the Lamb shift is probably small enough to be neglected for energetics of molecules containing light atoms if the target accuracy is 1 kcal/mol, but for higher accuracy calculations and for molecules containing heavy elements it must be considered.

Dyall, Kenneth G.↗

A New Interpretation of the CrO-Photoelectron Detachment Experiment

CASSCF/MRCI calculations are performed on CrO(-) and CrO. The ground state of CrO(-) is identified as 4-II, not 6-Epsilon(+), as recently suggested. We assign the experimental detachment features to the seven one-electron processes arising from the two ion states.

Bauschlicher, Charles W., Jr.↗

A Theoretical Investigation of the Infrared Spectroscopic Properties of Closed-Shell Polycyclic Aromatic Hydrocarbon Cations

Density functional theory has been employed to calculate the harmonic frequencies and intensities of a range of PAH cations which explore both size and electronic structure effects on the infrared spectroscopic of these species. The sample extends the size range of PAH species considered to more than 50 carbon atoms and includes several representatives from each of two heretofore unexplored categories of PAH cations: (1) fully benzenoid PAH cations whose carbon skeleton is composed of an odd number of carbon atoms and (2) protonated PAH cations. Unlike the radical electronic structures of the PAH cations that have been the subject of previous theoretical and experimental work, the species in these two classes have a closed-shell electronic configuration. The calculated spectra of circumcoronene, C54H18, in both neutral and (radical) cationic form are also reported and compared to those of the other species. Closed-shell species are inherently less reactive than radical (or open-shell) cations and are known to play a role in combustion chemistry. Since interstellar PAHs are typically exposed to abundant atomic hydrogen and are thought to originate under pseudo-combustion conditions in carbon-rich circumstellar shells, such species may represent an important component of the interstellar PAH population. Furthermore, species larger than 50 carbon atoms are more representative of the size of the typical interstellar PAH. Overall, as has been the case for previous studies of PAH radical cations, the general pattern of band positions and intensities are consistent with that of the interstellar infrared emission spectrum. In addition, the spectra of closed-shell and open-shell cations are found to converge with increasing molecular size and are found to be relatively similar for species containing about 50 carbon atoms.

Hudgins, Douglas M.↗

The Reactions of Polycyclic Aromatic Hydrocarbons with OH

The OH radical adds to naphthalene and naphthalene cation without a barrier. For the neutrals, the most favorable path for this intermediate is the loss of the OH, and the next most favorable option is the loss of an H atom to form the alcohol. For the cation, the most favorable path appears to be a hydrogen migration followed by the loss of a hydrogen to form the alcohol. The OH at carbon atom 1 is energetically most favorable for both the initial complex and final product. This is true for both the neutrals and cations.

Ricca, Alessandra↗