Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical force fields”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

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↗
NASA NTRS2022

Building and Breaking Carbon Composites with REACTER

Carbon-based composites have become indispensable materials in aerospace and other high-performance applications, yet obtaining a detailed, nanoscale understanding of their morphology and failure mechanisms using only experimental methods remains a difficult challenge. REACTER is a versatile computational modeling tool for atomistic molecular dynamics simulations designed to model chemical reactions at the speed and length scales of classical force fields. In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites: carbon nanotube (CNT) composites and carbon fiber reinforced polymers (CFRP). A network of CNTs was grown dynamically using the new ‘create atoms’ feature of REACTER. The CNT filler was embedded into a polyarylacetylene (PAA) matrix by simulated in situ polymerization to obtain the final composite model. To generate the second class of carbon composite, fully carbonized (graphitic) carbon fiber morphologies were created by the method of Desai et al. [1], but using the advanced reaction constraints framework of REACTER. Two fiber models were created, representing a circular carbon fiber core and a flat surface, and similarly infiltrated with resin to obtain the final CFRP structure. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third order DFT-based tight-binding simulations, via a reaction constraint on the total potential energy of the involved atoms.

Polymer↗
NASA NTRSDate not supplied

REACTER: A Versatile Tool for Large-Scale Reactive Molecular Dynamics

Accurately describing reactive events over long length and time scales remains a grand challenge of computational materials science. REACTER is a general protocol for modeling chemical reactions using classical force fields, and is implemented in the popular molecular dynamics software LAMMPS. REACTER has a growing user base and has been used as a model-building tool for a variety of materials, including thermoplastics, thermosets, glassy materials and composites. The method has also been applied to accelerated modeling of reversible chemical reactions, such as the formation of electrochemical components for batteries. Recently, the REACTER protocol has received some major upgrades to enhance its ability to predict when reactions occur and to make it easier to use. Force field parameters can now be automatically assigned to newly created bonds, angles and other interactions. Advanced reaction constraints have also been added, including an Arrhenius constraint to enforce an effective activation energy, a root-mean-square-deviation option for complex geometrical constraints, and a constraint based on the total potential energy of the atoms involved in a reactive site. This potential energy constraint allows for the accurate reproduction of DFT-based tight-binding (DFTB3) predicted bond dissociation curves for mechanically induced bond breaking.

polymer simulations, molecular dynamics↗
NASA NTRS2002

Effects of Structural Deformation and Tube Chirality on Electronic Conductance of Carbon Nanotubes

A combination of large scale classical force-field (UFF), density functional theory (DFT), and tight-binding Green's function transport calculations is used to study the electronic properties of carbon nanotubes under the twist, bending, and atomic force microscope (AFM)-tip deformation. We found that in agreement with experiment a significant change in electronic conductance can be induced by AFM-tip deformation of metallic zigzag tubes and by twist deformation of armchair tubes. The effect is explained in terms of bandstructure change under deformation.

Svizhenko, Alexei↗
NASA NTRSDate not supplied

REACTER: A Heuristic Method for Reactive Molecular Dynamics

REACTER (www.reacter.org) is a heuristic protocol that enables the simulation of complex reactions using atomistic molecular dynamics (MD) with a fixed-valence force field. Incorporating reactions into classical MD with this approach allows modeling of reactive systems over greatly-increased time scales, enabling systems to be modeled with MD that would not otherwise be feasible. One or more competing multi-step reactions or series of reactions can be invoked simultaneously. Special treatment can be applied to neighboring atoms to relax high energy configurations while the simulation progresses. The original version of REACTER, which was implemented in the open-source LAMMPS simulation package as fix bond/react, was only available for serial simulations. This work describes the expansion of fix bond/react for use in parallel simulations, as well as the addition of various new options, including deletion of reaction by-products, reversible reactions, and custom reaction constraints. These new capabilities are demonstrated through large-scale simulations (200,000+ atoms) of the polymerization of polystyrene and nylon 6,6. The morphologies of both polymers are analyzed after reaching >99% extent of polymerization. Finally, the newly-added reversible reactions feature is demonstrated by rupturing these highly-entangled systems under uniaxial strain by defining a chain scission reaction.

polymer simulations, molecular dynamics↗
NASA NTRS1973

Study of liquid-solid transition for materials processing in space. Analysis of field effects on dense liquid materials

The behavior of dense liquids near the solidification point under the influence of magnetic fields or near-zero gravity conditions were analyzed within the framework of existing liquid state models and classical field theory. The dynamic body forces which are an effect of the magnetic field on liquids are discussed. It is concluded that both magnetic fields and low gravity conditions produce an increase in the diffusion coefficient which result in an increased growth rate of crystals, and time varying magnetic fields induce eddy currents in the liquid which produce body forces and tend to disrupt convection.

Miller, R. I.↗
NASA NTRS1996

Magnetothermal Convection in Nonconducting Diamagnetic and Paramagnetic Fluids

Nonuniform magnetic fields exert a magnetic body force on electrically nonconducting classical fluids. These include paramagnetic fluids such as gaseous and liquid oxygen and diamagnetic fluids such as helium. Recent experiments show that this force can overwhelm the force of gravity even at the surface of the earth; it can levitate liquids and gases, quench candle flames, block gas flows, and suppress heat transport. Thermal gradients render the magnetic force nonuniform through the temperature-dependent magnetic susceptibility. These thermal gradients can therefore drive magnetic convection analogous to buoyancy-driven convection. This magnetothermal convection can overwhelm convection driven by gravitational buoyancy in terrestrial experiments. The objectives of the proposed ground-based theoretical study are (a) to supply the magnetothermohydrodynamic theory necessary to understand these recent experiments and (b) to explore the consequences of nonuniform magnetic fields in microgravity. Even the linear theory for the onset of magnetothermal convection is lacking in the literature. We intend to supply the linear and nonlinear theory based on the thermohydrodynamic equations supplemented by the magnetic body force. We intend to investigate the effect of magnetic fields on gas blockage and heat transport in microgravity. Since magnetic fields provide a means of creating arbitrary, controllable body force distributions, we intend to investigate the possibility of using magnetic fields to position and control fluids in microgravity. We also intend to investigate the possibility of creating stationary terrestrial microgravity environments by using the magnetic force to effectively cancel gravity. These investigations may aid in the design of space-based heat-transfer, combustion, and human-life-support equipment.

Edwards, Boyd F.↗
NASA NTRS1978

Expansion and broadening of coronal loop transients - A theoretical explanation

Consequences are examined of the assumption that an observed coronal loop transient is a twisted rope of magnetic-field lines expanding and broadening in the background coronal plasma and magnetic field. It is shown that the expansion can be accounted for by the azimuthal component of the field; the observed broadening of the loop as it moves outward can be accounted for by the longitudinal component of the field. In order to have a net outward force and at the same time avoid a classical pinch (sausage) instability, the two components of the field must satisfy a certain inequality. It is predicted that, as the loop rises, the width (h) of its top portion should vary proportionally with distance (R) from the sun's center. This is in good agreement with measurements that show h is proportional to the 0.8 power of R. The prediction that the radius of curvature of the top portion of the loop should be proportional to R differs from the measured variation. The difference could be accounted for by a drag due to the background coronal field that flattens the loop's top.

Mouschovias, T. CH.↗
NASA NTRS1998

Computational Investigation of Helical Traveling Wave Tube Transverse RF Field Forces

In a previous study using a fully three-dimensional (3D) helical slow-wave circuit cold- test model it was found, contrary to classical helical circuit analyses, that transverse FF electric fields have significant amplitudes compared with the longitudinal component. The RF fields obtained using this helical cold-test model have been scaled to correspond to those of an actual TWT. At the output of the tube, RF field forces reach 61%, 26% and 132% for radial, azimuthal and longitudinal components, respectively, compared to radial space charge forces indicating the importance of considering them in the design of electron beam focusing.

Kory, Carol L.↗
NASA NTRS1990

On the numerical computation of nonlinear force-free magnetic fields

An algorithm has been developed to extrapolate nonlinear force-free magnetic fields from the photosphere, given the proper boundary conditions. This paper presents the results of this work, describing the mathematical formalism that was developed, the numerical techniques employed, and comments on the stability criteria and accuracy developed for these numerical schemes. An analytical solution is used for a benchmark test; the results show that the computational accuracy for the case of a nonlinear force-free magnetic field was on the order of a few percent (less than 5 percent). This newly developed scheme was applied to analyze a solar vector magnetogram, and the results were compared with the results deduced from the classical potential field method. The comparison shows that additional physical features of the vector magnetogram were revealed in the nonlinear force-free case.

Wu, S. T.↗
NASA NTRS2015

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.↗
NASA NTRS1989

Potential theory of radiation

A theoretical method is being developed by which the structure of a radiation field can be predicted by a radiation potential theory, similar to a classical potential theory. The introduction of a scalar potential is justified on the grounds that the spectral intensity vector is irrotational. The vector is also solenoidal in the limits of a radiation field in complete radiative equilibrium or in a vacuum. This method provides an exact, elliptic type equation that will upgrade the accuracy and the efficiency of the current CFD programs required for the prediction of radiation and flow fields. A number of interesting results emerge from the present study. First, a steady state radiation field exhibits an optically modulated inverse square law distribution character. Secondly, the unsteady radiation field is structured with two conjugate scalar potentials. Each is governed by a Klein-Gordon equation with a frictional force and a restoring force. This steady potential field structure and the propagation of radiation potentials are consistent with the well known results of classical electromagnetic theory. The extension of the radiation potential theory for spray combustion and hypersonic flow is also recommended.

Chiu, Huei-Huang↗
NASA NTRS1989

Multilayer theory for delamination analysis of a composite curved bar subjected to end forces and end moments

A composite test specimen in the shape of a semicircular curved bar subjected to bending offers an excellent stress field for studying the open-mode delamination behavior of laminated composite materials. This is because the open-mode delamination nucleates at the midspan of the curved bar. The classical anisotropic elasticity theory was used to construct a multilayer theory for the calculations of the stress and deformation fields induced in the multilayered composite semicircular curved bar subjected to end forces and end moments. The radial location and intensity of the open-mode delamination stress were calculated and were compared with the results obtained from the anisotropic continuum theory and from the finite element method. The multilayer theory gave more accurate predictions of the location and the intensity of the open-mode delamination stress than those calculated from the anisotropic continuum theory.

Ko, William L.↗
NASA NTRS1989

Multilayer theory for delamination analysis of a composite curved bar subjected to end forces and end moments

A composite test specimen in the shape of a semicircular curved bar subjected to bending offers an excellent stress field for studying the open-mode delamination behavior of laminated composite materials. This is because the open-mode delamination nucleates at the midspan of the curved bar. The classical anisotropic elasticity theory was used to construct a 'multilayer' theory for the calculations of the stress and deformation fields induced in the multilayered composite semicircular curved bar subjected to end forces and end moments. The radial location and intensity of the open-mode delamination stress were calculated and were compared with the results obtained from the anisotropic continuum theory and from the finite element method. The multilayer theory gave more accurate predictions of the location and the intensity of the open-mode delamination stress than those calculated from the anisotropic continuum theory.

Ko, William L.↗
NASA NTRS2017

Modeling Chemical Reactions in Classical Molecular Dynamics Simulations

An algorithm capable of incorporating multi-step reaction mechanisms into atomistic molecular dynamics (MD) simulations using traditional fixed valence force fields is proposed and implemented within the framework of LAMMPS (Large-scale Atomic Molecular Massively Parallel Simulator). This extension, referred to as fix bond/react, enables bonding topology modifications during a running MD simulation using pre- and post-reaction bonding templates to carry out a pre-specified reaction. Candidate reactants are first identified by interatomic separation, followed by the application of a generalized topology matching algorithm to confirm they match the pre-reaction template. This is followed by a topology conversion to match the post-reaction template and a dynamic relaxation to minimize high energy configurations. Two case studies, the condensation polymerization of nylon 6,6 and the formation of a highly-crosslinked epoxy, are simulated to demonstrate the robustness, stability, and speed of the algorithm. Improvements which could increase its utility are discussed.

Jacob R Gissinger↗
NASA NTRS1990

Studies of the thermal and optical responses of H atoms in solid H2

It was the goal of this reserch project to model both the storage of energy in solid hydrogen in the form of atoms and the conversion of this stored energy into other forms of useful energy. The basic ideas of rocket propulsion originate in classical physics and they remain unchanged. To escape a strong gravitational field, the 'burn time' must be minimized but in negligible force fields, the burn time is unimportant and only the relative masses of rocket to fuel determine a specific exhaust velocity. It is in this latter case that low mass fuels such as hydrogen become very important. The burning of hydrogen in oxygen is a 'benchmark' fuel today providing a specific impulse of 400 seconds or better. More exotic fuels will be needed for many of the interesting explorations of the future but they still must have large energy releases per unit mass. It is in this context that propulsion based on hydrogen atom recombination receives attention and these studies will serve as engineering guides.

Gaines, James R.↗
NASA NTRS1986

Simulation of the thermospheric tides by use of the NCAR thermospheric general circulation model

Numerical calculations of the thermospheric tidal winds and temperatures at equinox are presented. The calculations were made using the National Center for Atmospheric Research Thermospheric General Circulation Model (TGCM), which includes the effects of viscosity, conductivity, diffusion, ion drag, winds, and temperature gradients. The thermospheric diurnal and semidiurnal tides are excited in situ by solar heating and by ion-neutral momentum coupling. The semidiurnal tidal calculations also include the effects of upward propagating waves generated by heating in the lower atmosphere. This semidiurnal propagating component is modeled by use of the classical tidal perturbations as lower boundary conditions. The model is tuned by adjusting the propagating tidal forcing term until calculated semidiurnal wind and temperature fields best approximate incoherent scatter observations. The tidal TGCM results are consistent with previous theoretical work and successfully reproduce high altitude temperature and meridional velocity data, but they give significantly lower magnitudes for velocities and temperatures near 160 km than are seen by observations.

Fensen, C. G.↗
NASA NTRS1990

Random field estimation approach to robot dynamics

The difference equations of Kalman filtering and smoothing recursively factor and invert the covariance of the output of a linear state-space system driven by a white-noise process. Here it is shown that similar recursive techniques factor and invert the inertia matrix of a multibody robot system. The random field models are based on the assumption that all of the inertial (D'Alembert) forces in the system are represented by a spatially distributed white-noise model. They are easier to describe than the models based on classical mechanics, which typically require extensive derivation and manipulation of equations of motion for complex mechanical systems. With the spatially random models, more primitive locally specified computations result in a global collective system behavior equivalent to that obtained with deterministic models. The primary goal of applying random field estimation is to provide a concise analytical foundation for solving robot control and motion planning problems.

Rodriguez, Guillermo↗