Search NASA⌕ Search

SEARCH · Search NASA

Results for “quantum mechanical simulations”

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 73 records · Page 4

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↗

Photonic Interrogation and Control of Nano Processes

My research activities for the summer of 2003 consisted of two projects: One project was concerned with determining a method for predicting the static and dynamic assembly properties of nano-structures using laser tweezers. The other project was to investigate the generation of Laguerre-Gaussian modes using a spatial light modulator incorporated into an optical tweezers system. Concerning the first project, I initially pursued the approach suggested by my NASA colleague Dr. Art Decker. This approach involved mimicking the model of the structure of atomic nucleus for the assembly of 1 to 100 atoms using allowed quadruple transitions induced by orbital angular momentums of a Laguerre- Gaussian (Doughnut) laser mode. After realizing the inaptness of the nuclear model with the nanostructure model as far as the binding forces and transitions were concerned, I focused on using quantum dot modei. This model was not attuned also for the host lattice influences the electronic structure of the quantum dot. Thus one other option that I decided to pursue was the approach of molecular quantum mechanics. In this approach the nanostructure is treated as a large (10-100 nm) molecule constructed from single element or multi-elements. Subsequent to consultation with Dr. Fred Morales, a chemical engineer at NASA GRC, and Dr. David Ball, a computational chemist at Cleveland State University, I acquired a molecular-quantum computation software, Hyperchem 7.0. This software allows simulation of different molecular structures as far as their static and dynamic behaviors are concerned. The time that I spent on this project was about eight weeks. Once this suitable approach was identified, I realized the need to collaborate with a computational quantum chemist to pursue searching for stable nanostructures in the range of 10-100 nm that we can be assembled using laser tweezers. The second project was about generating laser tweezers that possess orbital angular momentum. As shown, we were able to generate laser tweezers modes of different orbital angular momentum using a spatial light modulator incorporated into a laser tweezers system. The motivation for investigating these types of modes stems from being able to spin particles at high speeds and also to orient two particles in separate traps and then join them together. Also, there has been recent intense interest on fundamental physics research on orbital angular momentum of light. The fact that circularly polarized light may have associated with it angular momentum that relates to the spin of individual photons (spin 0 for the plane polarized light, spin +1 for the right-circularly polarized light and spin -1 for the left-circularly polarized light) was first demonstrated by Beth in 1936. Orbital angular momentum is, however, distinct from spin in that the spin angular momentum of light is intrinsically linked to the behavior of the electric field in the light whereas orbital angular momentum is a consequence of inclined wavefronts. In 1992 L. Allen, et al showed that the Laguerre-Gaussian (LG) modes could possess well-defined orbital angular momentum that can exceed 1 planck's constant, i.e. l plancks constant per photon, where l is the azimuthal index of the mode.

Jassemnejad, Baha↗

Mechanisms for the Formations of the Thymine Under Astrophysical Conditions and Implications for the Origin of Life

Nucleobases are the carriers of the genetic information in ribonucleic acid and deoxyribonucleic acid (DNA) for all life on Earth. Their presence in meteorites clearly indicates that compounds of biological importance can form via non-biological processes in extraterrestrial environments. Recent experimental studies have shown that the pyrimidine-based nucleobases uracil and cytosine can be easily formed from the ultraviolet irradiation of pyrimidine in H2O-rich ice mixtures that simulate astrophysical processes. In contrast, thymine, which is found only in DNA, is more difficult to form under the same experimental conditions, as its formation usually requires a higher photon dose. Earlier quantum chemical studies confirmed that the reaction pathways were favorable provided that several H2O molecules surrounded the reactants. However, the present quantum chemical study shows that the formation of thymine is limited because of the inefficiency of the methylation of pyrimidine and its oxidized derivatives in an H2O ice, as supported by the laboratory studies. Our results constrain the formation of thymine in astrophysical environments and thus the inventory of organic molecules delivered to the early Earth and have implications for the role of thymine and DNA in the origin of life.

Origin↗

Yields of O2(b 1 Sigma g +) from reactions of HO2

The production of O2(b 1 Sigma g +) has been monitored for several reactions of the HO2 radical at 300 K using a discharge-flow apparatus with resonance fluorescence and chemiluminescence detection. In all cases, the resulting quantum efficiencies were found to be less than 0.03. O2(b) was observed when F atoms were added to H2O2 in the gas phase. The signal strengths of O2(b) were proportional to initial concentrations of HO2 formed by the F + H2O2 reaction. Observed /O2(b)/, /HO2/, and /OH/ vs /F/0 were analyzed using a simple three-step mechanism and a more complete computer simulation with 22 reaction steps. The results indicate that the F + HO2 reaction yields O2(b) with an efficiency of (3.6 + or - 1.4) x 10 to the -3rd. Yields from the O + OH2 reaction were less than 0.02, indicating that this reaction cannot be a major source of the O2(b) emission observed in the earth's nightglow.

Keyser, L. F.↗

Development of Large-Scale Magnetic Calorimeter Arrays

We describe performance of large-scale arrays of metallic magnetic calorimeters (MMCs) we are developing to meet requirements of the Lynx X-ray Microcalorimeter (LXM) instrument in the astrophysics mission concept Lynx. We have fabricated prototypes with 55,800 x-ray pixels thermally connected to 5,688 MMC sensors. Subarrays demonstrate three types of pixels, which have different energy and spatial resolution goals in LXM. Pixel pitch is 50 or 25 µm. For two subarrays, use of position sensitive detectors, in a "hydra" configuration with 25 x-ray absorbers per sensor, helps achieve a large focal plane by increasing pixel count relative to readout channels. Since each absorber has a different thermal link to its sensor, it generates a different pulse shape and enables discrimination of pixel position. Superconducting wiring from all sensors was brought out to the perimeter of the overall array using multiple buried metal layers planarized by Chemical Mechanical Planarization to achieve high critical current, low inductance, and high fabrication yield. An automated, algorithmic approach was used to layout the complex wiring pattern. For readout with existing small arrays of Superconducting Quantum Interference (SQUID) amplifiers, 112 selected pixels were connected to wire bond pads. In hydras, we successfully identified 25 pulse shapes using rise-time and pulse height. Pulse shapes were similar to simulations. We measured noise and responsivity to substrate temperature and applied x-ray pulses. While we designed sensors to have high inductance to relative to stray inductance of wiring, we had only mismatched (lower input inductance) SQUIDs available in our current apparatus, giving non-optimized resolution. We compared measured and calculated energy resolutions. We have the required sensitivity to reach < 3 eV FWHM resolution for 6 keV x-rays, required for the LXM Main Array, and the ability to discriminate the 25 different absorber pixels down to energies < 300 eV.

microcalorimeters↗

Modeling the Flow of Rarefied Gases at NASA

At modest temperatures, the thermal energy of atmospheric diatomic gases such as nitrogen is primarily distributed between only translational and rotational energy modes. Furthermore, these energy modes are fully excited such that the specific heat at constant volume is well approximated by the simple expression C(sub v) = 5/2 R. As a result, classical mechanics provides a suitable approximation at such temperatures of the true quantum mechanical behavior of the inter-molecular collisions of such molecules. Using classical mechanics, the transfer of energy between rotational and translation energy modes is studied. The approach of Lordi and Mates is adopted to compute the trajectories and time dependent rotational orientations and energies during the collision of two non-polar diatomic molecules. A Monte-Carlo analysis is performed collecting data from the results of many such simulations in order to estimate the rotational relaxation time. A Graphical Processing Unit (GPU) is employed to improve the performance of the Monte-Carlo analysis. A comparison of the performance of the GPU implementation to an implementation on traditional computer architecture is made. Effects of the assumed inter-molecular potential on the relaxation time are studied. The seminar will also present highlights of computational analyses performed at NASA Johnson Space Center of heat transfer in rarefied gases.

Forrest E. Lumpkin, III↗

Examination of Mars2020 Shock-Layer Conditions via Infrared Laser Absorption Spectroscopy of CO2 and CO

Amid-infrared laser absorption diagnostic was deployed to study a simulated Mars2020 shock layer in the Electric Arc Shock Tube (EAST) facility at NASA’s Ames Research Center. Rapid RF-diplexing techniques enabled quantitative temperature and number density measurements of CO2 and CO with μs-resolution over an incident shock velocity range of 1.39 – 3.75 km/s. Two interband cascade lasers were utilized at 4.17 and 4.19 μm to resolve eight CO2 asymmetric stretch fundamental band (𝜈3) transitions from two different vibrational levels: 0000 (ground state) and 0110 (first excited bending mode). The probed rotational levels span across J” = 58 to J” = 140. Results are compared to DPLR simulations of the shock layer using kinetic mechanisms of Johnston et al. and Cruden et al. At shock velocities below 3.1 km/s, the agreement between the measurements and the Johnston mechanism is typically within 5% for temperature and within 10% for number density. At shock velocities above 3.1 km/s, the CO2 measurement becomes sensitive to a thin boundary layer and corrections of this effect are presented. On test cases with enough energy to dissociate CO2, a quantum cascade laser scanned the P(2, 20), P(0, 31), and P(3, 14) transitions of the CO fundamental band at 4.98 μm. CO formation rate is measured to be close to the Johnston kinetic mechanism at low velocities, and then trending towards the Cruden kinetic mechanism at high velocities. On a few low velocity test cases, rovibrational relaxation of the Martian atmosphere is probed with μs resolution.

Laser Absorption↗

A Quantum Band Model of the nu3 Fundamental of Methanol (CH3OH) and Its Application to Fluorcescence Spectra of Comets

Methanol (CH3OH) radiates efficiently at infrared wavelengths, dominating the C-H stretching region in comets, yet inadequate quantum-mechanical models have imposed limits on the practical use of its emission spectra. Accordingly, we constructed a new line-by-line model for the 3 fundamental band of methanol at 2844 / cm (3.52 micron) and applied it to interpret cometary fluorescence spectra. The new model permits accurate synthesis of line-by-line spectra for a wide range of rotational temperatures, ranging from 10 K to more than 400 K.We validated the model by comparing simulations of CH3OH fluorescent emission with measured spectra of three comets (C/2001 A2 LINEAR, C/2004 Q2 Machholz and 8P/Tuttle) acquired with high-resolution infrared spectrometers at high-altitude sites. The new model accurately describes the complex emission spectrum of the nu3 band, providing distinct rotational temperatures and production rates at greatly improved confidence levels compared with results derived from earlier fluorescence models. The new model reconciles production rates measured at infrared and radio wavelengths in C/2001 A2 (LINEAR). Methanol can now be quantified with unprecedented precision and accuracy in astrophysical sources through high-dispersion spectroscopy at infrared wavelengths

Villanueva, Geronimo L.↗

Machine Phase Fullerene Nanotechnology: 1996

NASA has used exotic materials for spacecraft and experimental aircraft to good effect for many decades. In spite of many advances, transportation to space still costs about $10,000 per pound. Drexler has proposed a hypothetical nanotechnology based on diamond and investigated the properties of such molecular systems. These studies and others suggest enormous potential for aerospace systems. Unfortunately, methods to realize diamonoid nanotechnology are at best highly speculative. Recent computational efforts at NASA Ames Research Center and computation and experiment elsewhere suggest that a nanotechnology of machine phase functionalized fullerenes may be synthetically relatively accessible and of great aerospace interest. Machine phase materials are (hypothetical) materials consisting entirely or in large part of microscopic machines. In a sense, most living matter fits this definition. To begin investigation of fullerene nanotechnology, we used molecular dynamics to study the properties of carbon nanotube based gears and gear/shaft configurations. Experiments on C60 and quantum calculations suggest that benzyne may react with carbon nanotubes to form gear teeth. Han has computationally demonstrated that molecular gears fashioned from (14,0) single-walled carbon nanotubes and benzyne teeth should operate well at 50-100 gigahertz. Results suggest that rotation can be converted to rotating or linear motion, and linear motion may be converted into rotation. Preliminary results suggest that these mechanical systems can be cooled by a helium atmosphere. Furthermore, Deepak has successfully simulated using helical electric fields generated by a laser to power fullerene gears once a positive and negative charge have been added to form a dipole. Even with mechanical motion, cooling, and power; creating a viable nanotechnology requires support structures, computer control, a system architecture, a variety of components, and some approach to manufacture. Additional information is contained within the original extended abstract.

Globus, Al↗

Modeling and simulation of space-based entanglement distribution

In this talk we present a model for analyzing space-based entanglement distribution links. In particular, we analyze deterministic entanglement distribution links relevant to quantum communication applications where all communication nodes must be synchronized with deterministic photon arrival times. The focus is on systems operating within practical limitations set by current technology. The architectures studied in this work are based on all-photonic links with probabilistic (memory-less) entanglement swapping used to signal the success of the entanglement distribution without destroying the entangled state. An end-to-end link model is described which accounts for satellite orbital mechanics as well as the physical effects of entangled photon pair sources, space-to-ground atmospheric channels, two-photon interference, pointing and tracking, and photon counting detectors with and without photon number resolution (this is joint work with Robert Manning, Ian Nemitz, and Bryan Welch).

quantum communication↗

Chemisorption and Diffusion of H on a Graphene Sheet and Single-Wall Carbon Nanotubes

Recent experiments on hydrogen storage in single wall nanotubes and nanotube bundles have reported large fractional weight of stored molecular hydrogen which are not in agreement with theoretical estimates based of simulation of hydrogen storage by physisorption mechanisms. Hydrogen storage in catalytically doped nanotube bundles indicate that atomic H might undergo chemisorption changing the basic nature of the storage mechanism under investigation by many groups. Using a generalized tight-binding molecular dynamics (GTBMD) method for reactive C-H dynamics, we investigate chemisorption and diffusion of atomic H on graphene sheet and C nanotubes. Effective potential energy surfaces (EPS) for chemisorption and diffusion are calculated for graphene sheet and nanotubes of different curvatures. Analysis of the activation barriers and quantum rate constants, computed via wave-packet dynamics method, will be discussed in this presentation.

Srivastava, Deepak↗

Digital Technologies at NASA for Science and Engineering

While scientific and engineering advancements used to rely primarily on theoretical studies and physical experiments, today digital technology enabled by petaflops-scale supercomputers is an equal, if not a greater, contributor to such achievements. In addition, computational modeling and simulation serves as a predictive tool that is not otherwise available. As a result, the use of high performance computing is integral to NASA's work in all mission areas such as space exploration, aeronautics, and scientific discovery. But traditional supercomputing alone is not sufficient for all of the space agency's needs. The success of many NASA missions depends on solving complex computing challenges, some of which are NP-hard (decision theory) if using classical solution methods. Quantum computing promises an unprecedented ability to solve such intractable problems by harnessing quantum mechanical effects such as tunneling, superposition, and entanglement. Another disruptive digital technology is neuromorphic computing that uses brain-inspired lessons to generate new architectures that are much more energy efficient, and capable of massive parallel processing and learning in-situ. Finally, with large amounts of observational and computational data sets, the opportunities of big data and data analytics can be leveraged to enable deep learning and knowledge discovery - it's all a massive digital transformation. This talk will be an overview how NASA utilizes digital technologies for its science and engineering efforts.

Biswas, Rupak↗

Limits on the UV Photodecomposition of Carbonates on Mars

The effect of UV (ultraviolet) light on the stability of calcium carbonate in a simulated martian atmosphere was experimentally investigated. Sample cells containing C-13 labeled calcite were irradiated with a Xe arc lamp in 10 mbar of simulated martian atmosphere and a quadrupole mass spectrometer was used to monitor the headspace gases for the production of (13)CO2. We found no experimental evidence of the UV photodecomposition of calcium carbonate in a simulated martian atmosphere. Extrapolating the lower limit of detection of our experimental system to an upper limit of carbonate decomposition on Mars yields a quantum efficiency of 3.5 x 10(exp -8) molecules/photon over the wavelength interval of 190-390 nm and a maximum UV photodecomposition rate of 1.2 x 10(exp -13) kg m(exp -2) s(exp -1) from a calcite surface. The maximum loss of bulk calcite due to this process would be 2.5 nm yr(exp -1). However, calcite is expected to be thermodynamically stable on the surface of Mars and potential UV photodecomposition reaction mechanisms indicate that while calcium carbonate may decompose under vacuum, it would be stable in a CO2 atmosphere. Given the expected stability of carbonate on Mars and our inability to detect carbonate decomposition, we conclude that it is unlikely that the apparent absence of carbonate on the martian surface is due to UV photo decomposition of calcite in the current environment.

Quinn, Richard↗

Nanoscale Etching and Indentation of Silicon Surfaces with Carbon Nanotubes

The possibility of nanolithography of silicon and germanium surfaces with bare carbon nanotube tips of scanning probe microscopy devices is considered with large scale classical molecular dynamics (MD) simulations employing Tersoff's reactive many-body potential for heteroatomic C/Si/Ge system. Lithography plays a key role in semiconductor manufacturing, and it is expected that future molecular and quantum electronic devices will be fabricated with nanolithographic and nanodeposition techniques. Carbon nanotubes, rolled up sheets of graphene made of carbon, are excellent candidates for use in nanolithography because they are extremely strong along axial direction and yet extremely elastic along radial direction. In the simulations, the interaction of a carbon nanotube tip with silicon surfaces is explored in two regimes. In the first scenario, the nanotubes barely touch the surface, while in the second they are pushed into the surface to make "nano holes". The first - gentle scenario mimics the nanotube-surface chemical reaction induced by the vertical mechanical manipulation of the nanotube. The second -digging - scenario intends to study the indentation profiles. The following results are reported in the two cases. In the first regime, depending on the surface impact site, two major outcomes outcomes are the selective removal of either a single surface atom or a surface dimer off the silicon surface. In the second regime, the indentation of a silicon substrate by the nanotube is observed. Upon the nanotube withdrawal, several surface silicon atoms are adsorbed at the tip of the nanotube causing significant rearrangements of atoms comprising the surface layer of the silicon substrate. The results are explained in terms of relative strength of C-C, C-Si, and Si-Si bonds. The proposed method is very robust and does not require applied voltage between the nanotube tips and the surface. The implications of the reported controllable etching and hole-creating for nanolithography on silicon are discussed in detail.

Dzegilenko, Fedor N.↗

Coupled Structural, Thermal, Phase-change and Electromagnetic Analysis for Superconductors, Volume 2

Described are the theoretical development and computer implementation of reliable and efficient methods for the analysis of coupled mechanical problems that involve the interaction of mechanical, thermal, phase-change and electromag subproblems. The focus application has been the modeling of superconductivity and associated quantum-state phase change phenomena. In support of this objective the work has addressed the following issues: (1) development of variational principles for finite elements, (2) finite element modeling of the electromagnetic problem, (3) coupling of thermel and mechanical effects, and (4) computer implementation and solution of the superconductivity transition problem. The main accomplishments have been: (1) the development of the theory of parametrized and gauged variational principles, (2) the application of those principled to the construction of electromagnetic, thermal and mechanical finite elements, and (3) the coupling of electromagnetic finite elements with thermal and superconducting effects, and (4) the first detailed finite element simulations of bulk superconductors, in particular the Meissner effect and the nature of the normal conducting boundary layer. The theoretical development is described in two volumes. Volume 1 describes mostly formulation specific problems. Volume 2 describes generalization of those formulations.

Felippa, C. A.↗

Coupled Structural, Thermal, Phase-Change and Electromagnetic Analysis for Superconductors

Described are the theoretical development and computer implementation of reliable and efficient methods for the analysis of coupled mechanical problems that involve the interaction of mechanical, thermal, phase-change and electromagnetic subproblems. The focus application has been the modeling of superconductivity and associated quantum-state phase-change phenomena. In support of this objective the work has addressed the following issues: (1) development of variational principles for finite elements, (2) finite element modeling of the electromagnetic problem, (3) coupling of thermal and mechanical effects, and (4) computer implementation and solution of the superconductivity transition problem. The main accomplishments have been: (1) the development of the theory of parametrized and gauged variational principles, (2) the application of those principled to the construction of electromagnetic, thermal and mechanical finite elements, and (3) the coupling of electromagnetic finite elements with thermal and superconducting effects, and (4) the first detailed finite element simulations of bulk superconductors, in particular the Meissner effect and the nature of the normal conducting boundary layer. The theoretical development is described in two volumes. This volume, Volume 1, describes mostly formulations for specific problems. Volume 2 describes generalization of those formulations.

Felippa, C. A.↗

Effects of Atomic-Scale Structure on the Fracture Properties of Amorphous Carbon - Carbon Nanotube Composites

The fracture of carbon materials is a complex process, the understanding of which is critical to the development of next generation high performance materials. While quantum mechanical (QM) calculations are the most accurate way to model fracture, the fracture behavior of many carbon-based composite engineering materials, such as carbon nanotube (CNT) composites, is a multi-scale process that occurs on time and length scales beyond the practical limitations of QM methods. The Reax Force Field (ReaxFF) is capable of predicting mechanical properties involving strong deformation, bond breaking and bond formation in the classical molecular dynamics framework. This has been achieved by adding to the potential energy function a bond-order term that varies continuously with distance. The use of an empirical bond order potential, such as ReaxFF, enables the simulation of failure in molecular systems that are several orders of magnitude larger than would be possible in QM techniques. In this work, the fracture behavior of an amorphous carbon (AC) matrix reinforced with CNTs was modeled using molecular dynamics with the ReaxFF reactive forcefield. Care was taken to select the appropriate simulation parameters, which can be different from those required when using traditional fixed-bond force fields. The effect of CNT arrangement was investigated with three systems: a single-wall nanotube (SWNT) array, a multi-wall nanotube (MWNT) array, and a SWNT bundle system. For each arrangement, covalent bonds are added between the CNTs and AC, with crosslink fractions ranging from 0-25% of the interfacial CNT atoms. The SWNT and MWNT array systems represent ideal cases with evenly spaced CNTs; the SWNT bundle system represents a more realistic case because, in practice, van der Waals interactions lead to the agglomeration of CNTs into bundles. The simulation results will serve as guidance in setting experimental processing conditions to optimize the mechanical properties of CNT composites.

Jensen, Benjamin D.↗

Coupled structural, thermal, phase-change and electromagnetic analysis for superconductors, volume 1

This research program has dealt with the theoretical development and computer implementation of reliable and efficient methods for the analysis of coupled mechanical problems that involve the interaction of mechanical, thermal, phase-change and electromagnetic subproblems. The focus application has been the modeling of superconductivity and associated quantum-state phase-change phenomena. In support of this objective the work has addressed the following issues: (1) development of variational principles for finite elements; (2) finite element modeling of the electromagnetic problem; (3) coupling of thermal and mechanical effects; and (4) computer implementation and solution of the superconductivity transition problem. The research was carried out over the period September 1988 through March 1993. The main accomplishments have been: (1) the development of the theory of parametrized and gauged variational principles; (2) the application of those principled to the construction of electromagnetic, thermal and mechanical finite elements; and (3) the coupling of electromagnetic finite elements with thermal and superconducting effects; and (4) the first detailed finite element simulations of bulk superconductors, in particular the Meissner effect and the nature of the normal conducting boundary layer. The grant has fully supported the thesis work of one doctoral student (James Schuler, who started on January 1989 and completed on January 1993), and partly supported another thesis (Carmelo Militello, who started graduate work on January 1988 completing on August 1991). Twenty-three publications have acknowledged full or part support from this grant, with 16 having appeared in archival journals and 3 in edited books or proceedings.

Felippa, Carlos A.↗