Search NASA⌕ Search

SEARCH · Search NASA

Results for “MOLECULAR DYNAMICS METHOD”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Effective Forces Between Colloidal Particles

Colloidal suspensions have proven to be excellent model systems for the study of condensed matter and its phase behavior. Many of the properties of colloidal suspensions can be investigated with a systematic variation of the characteristics of the systems and, in addition, the energy, length and time scales associated with them allow for experimental probing of otherwise inaccessible regimes. The latter property also makes colloidal systems vulnerable to external influences such as gravity. Experiments performed in micro-ravity by Chaikin and Russell have been invaluable in extracting the true behavior of the systems without an external field. Weitz and Pusey intend to use mixtures of colloidal particles with additives such as polymers to induce aggregation and form weak, tenuous, highly disordered fractal structures that would be stable in the absence of gravitational forces. When dispersed in a polarizable medium, colloidal particles can ionize, emitting counterions into the solution. The standard interaction potential in these charged colloidal suspensions was first obtained by Derjaguin, Landau, Verwey and Overbeek. The DLVO potential is obtained in the mean-field linearized Poisson-Boltzmann approximation and thus has limited applicability. For more precise calculations, we have used ab initio density functional theory. In our model, colloidal particles are charged hard spheres, the counterions are described by a continuum density field and the solvent is treated as a homogeneous medium with a specified dielectric constant. We calculate the effective forces between charged colloidal particles by integrating over the solvent and counterion degrees of freedom, taking into account the direct interactions between the particles as well as particle-counterion, counterion-counterion Coulomb, counterion entropic and correlation contributions. We obtain the effective interaction potential between charged colloidal particles in different configurations. We evaluate two- and three-body forces in the bulk as well as study the influence of soft walls. We qualitatively explain the effects of the walls on the forces and demonstrate that many-body effects are negligible in our system. With adjustments in the parameters, the DLVO pair-potential can describe the results quantitatively. Besides electrostatic interactions, entropic depletion effects that arise from (hard-core) exclusion play an important role in determining the behavior of multi-component colloidal suspensions. A standard theory for depletion forces is due to Asakura and Oosawa and is based on the ideal gas approximation. To go beyond this approximation, we have studied entropic forces in molecular dynamics simulations of systems of hard spheres (the effects of the solvent have been ignored). The effective depletion forces for these systems can be found either from equilibrium distribution functions or from direct momentum transfer calculations. Our results obtained by either method show qualitative differences from the Asakura-Oosawa forces, indicating a longer range, higher value at contact and most importantly a more complicated structure, comprising of several maxima and minima. Our calculations include the determination of effective forces between two spheres, a hard sphere and a wall, and the behavior of a hard sphere near a step-edge and a corner. We also demonstrate that such entropic forces do not necessarily satisfy pairwise additivity.

Tehver, Riina↗
NASA NTRS2000

Modeling Molecules

The molecule modeling method known as Multibody Order (N) Dynamics, or MBO(N)D, was developed by Moldyn, Inc. at Goddard Space Flight Center through funding provided by the SBIR program. The software can model the dynamics of molecules through technology which stimulates low-frequency molecular motions and properties, such as movements among a molecule's constituent parts. With MBO(N)D, a molecule is substructured into a set of interconnected rigid and flexible bodies. These bodies replace the computation burden of mapping individual atoms. Moldyn's technology cuts computation time while increasing accuracy. The MBO(N)D technology is available as Insight II 97.0 from Molecular Simulations, Inc. Currently the technology is used to account for forces on spacecraft parts and to perform molecular analyses for pharmaceutical purposes. It permits the solution of molecular dynamics problems on a moderate workstation, as opposed to on a supercomputer.

Source record↗
NASA NTRSDate not supplied

Effect of Vacuum on Force Response of an Ultrasonic Penetrator

Introduction: The Apollo astronauts encountered higher than expected resistances when interacting with the lunar soil via the Apollo Lunar Surface Drill (ALSD) and the trenching tool. Reducing the force required to move tools or other mechanical components through regolith will impact many steps of the resource extraction process. Force reduction has been achieved in soil materials by imparting vibration to tooling interfaces such as a vibratory farming cultivator, a percussive scoop, and ultrasonically resonant penetrators. Vibration-assisted tools in granular media reduce interaction forces by fluidizing a volume around the tool, allowing the tool to progress through a dynamic (fluid) medium instead of a static (solid) medium. This work seeks to quantify ultrasonic vibration’s effect on the force response of a penetrator in lunar soil simulant in vacuum sufficient to be within the molecular flow regime of any disturbed gases. Methods: A custom vacuum chamber setup, CUBEvac, was designed and built to facilitate penetration testing in a high vacuum environment, for comparison to penetration behavior in ambient terrestrial environment. A two-stage pumping system (Agilent Triscroll 600 roughing pump, Agilent VHS-6 oil diffusion pump) reached chamber pressures of about 5x10-6 Torr with regolith simulant in place. Figure 1 is a schematic of the heart of the assembly (note the penetration drive mechanisms above the chamber feedthrough and the regolith simulant sample in the bottom are not shown). The penetration actuation stack was comprised of a stepper motor driving a lead screw to move the ultra-sonic probe vertically inside the chamber. Motion was coordinated with an Arduino Uno. GRC-3 lunar simulant was used for this set of experiments. Samples were prepared in a four-liter stain-less steel, cylindrical pot with an internal diameter of approximately 15.56 cm (6.125 in) and a depth of 19.37 cm (7 5/8 in) for testing. The soil was baked out prior to compaction preparation as a measure to reduce soil moisture which interfered with pump down capacity. The soil was not baked again if it was removed from the vacuum chamber, prepped, and immediately returned to the vacuum chamber for pump down. The soil was compacted using a 60 Hz vibration table with a surcharge of 34 kg place on top of the soil in the container. Prepared soil samples weighed approximately 6.5 kg (bulk density 1.895-1.934 g/cm3). Two probe end effectors were tested: A cone penetrometer (static only) with a nominal diameter of 12.7 mm (0.5 in) and a nominal height of 28.6 mm (1.125 in); and a vibrating cylindrical probe measuring 12.7 mm in diameter and 50.8 mm in effective length from the tip (Figure 2). The cylindrical probe vibrated resonantly at 20 kHz with an amplitude 23 μm. The cone penetration tests were conducted to assess potential soil behavior differences in vacuum. The cylinder probe tests were conducted as the primary subject of this investigation to assess force response in vacuum. For each test, a regolith simulant sample was load-ed and compacted in the chamber, which was then evacuated for roughly 18 hours to reach the lowest possible pressure (approximately 5x10-6 Torr for most tests). The probe was then moved to about 10 mm above the soil surface before being pushed to a depth of 50 mm for the cylinder probe tests and to a depth of 100 mm for the cone penetrometer tests, both at a speed of 2 mm/s. The simulant samples were prepared the same for all tests. Ideally, they would respond consistently to probe penetration under ambient and vacuum conditions. This was evaluated by measuring the resistance of representative prepared simulant beds with a standard cone penetrometer in both environments. Results and Discussion: The resistance of the simulant samples in the vacuum tests was consistently lower than in the ambient tests as determined by the cone penetration tests. Thus, the ambient and vacuum results cannot be compared directly; work is underway to de-confound and better correlate the data. Still, figures 3-6 show that probe penetration forces are lower overall in the vacuum environment. In both environments, resonant vibration of the probe provides two useful effects: It reduces the probe penetration force and smooths the force-depth curve, significantly reducing local maxima. These effects have implications for various potential applications, such as astronaut hand-tools, where benefits (reducing astronaut effort) outweigh the cost of the additional energy re-quired to generate vibration. These results demonstrate that resonantly vibrating tools can meaningfully reduce the penetration force required for excavation, probing, and drilling tools in simulated lunar regolith deposits under vacuum levels approaching those that will be experience on the Moon’s surface. Lunar-gravity, ambient environment tests are scheduled soon. The effects of realistic temperatures and temperature gradients and deeper vacuum remain to be tested.

E Rezich↗
NASA NTRS2021

Molecular Dynamics Investigation of the Structural and Mechanical Properties of Off-Stoichiometric Epoxy Resins

Molecular dynamic (MD) simulations were performed to compute the mechanical properties of off-stoichiometric epoxy resins as a function of hardener/epoxy mixture ratio (r). Properties were characterized in relation to their microscopic structures. Such resins have been used recently for adhesive-free bonding of large-scale composite structures using the co-curing-ply method. In this process, two partially precured composite panels with hardener-poor (HP) off-stoichiometric resins are coupled with ply(ies) of complementary hardener-rich (HR) formulations and then cured simultaneously. This bonding process has the potential to produce reliable and certifiable composite joints without the need for additional fasteners, which are often required for many conventional bonding methods because even small amounts of contamination can cause a weak bond. The reflow and mixing of the HP/HR resin in this bonding process result in a joint with no discernable interface that should not be susceptible to surface contamination. However, incomplete mixing of the two offset resins may result in chemical heterogeneity of the cured polymeric joint. Thus, different r values may be obtained across the joint. Classical MD simulations were performed to compute the Young’s modulus of polymers with different r values and correlate their properties to network structures. High stiffness was associated with molecular packing due to chemical crosslinking, leading to a single network structure. Moreover, the networks became denser as the ratio approached the stoichiometric value r = 1. Thus, the r = 1 systems were single clusters, with high stiffness, high molecular weight, and a high degree of crosslinking. Structural properties such as radius of gyration and mean square displacement were determined to investigate the variation in the stiffness with respect to r. This MD simulation study was validated with experimental measurements.

Chang Woon Jang↗
NASA NTRS1988

A TEOM (tm) particulate monitor for comet dust, near Earth space, and planetary atmospheres

Scientific missions to comets, near earth space, and planetary atmospheres require particulate and mass accumulation instrumentation for both scientific and navigation purposes. The Rupprecht & Patashnick tapered element oscillating microbalance can accurately measure both mass flux and mass distribution of particulates over a wide range of particle sizes and loadings. Individual particles of milligram size down to a few picograms can be resolved and counted, and the accumulation of smaller particles or molecular deposition can be accurately measured using the sensors perfected and toughened under this contract. No other sensor has the dynamic range or sensitivity attained by these picogram direct mass measurement sensors. The purpose of this contract was to develop and implement reliable and repeatable manufacturing methods; build and test prototype sensors; and outline a quality control program. A dust 'thrower' was to be designed and built, and used to verify performance. Characterization and improvement of the optical motion detection system and drive feedback circuitry was to be undertaken, with emphasis on reliability, low noise, and low power consumption. All the goals of the contract were met or exceeded. An automated glass puller was built and used to make repeatable tapered elements. Materials and assembly methods were standardized, and controllers and calibrated fixtures were developed and used in all phases of preparing, coating and assembling the sensors. Quality control and reliability resulted from the use of calibrated manufacturing equipment with measurable working parameters. Thermal and vibration testing of completed prototypes showed low temperature sensitivity and high vibration tolerance. An electrostatic dust thrower was used in vacuum to throw particles from 2 x 10(exp 6) g to 7 x 10(exp -12) g in size. Using long averaging times, particles as small as 0.7 to 4 x 10(exp 11) g were weighted to resolutions in the 5 to 9 x 10(exp -13) g range. The drive circuit and optics systems were developed beyond what was anticipated in the contract, and are now virtually flight prototypes. There is already commercial interest in the developed capability of measuring picogram mass losses and gains. One area is contamination and outgassing research, both measuring picogram losses from samples and collecting products of outgassing.

Source record↗
NASA NTRS2004

Computational studies of Ras and PI3K

Until recently, experimental techniques in molecular cell biology have been the primary means to investigate biological risk upon space radiation. However, computational modeling provides an alternative theoretical approach, which utilizes various computational tools to simulate proteins, nucleotides, and their interactions. In this study, we are focused on using molecular mechanics (MM) and molecular dynamics (MD) to study the mechanism of protein-protein binding and to estimate the binding free energy between proteins. Ras is a key element in a variety of cell processes, and its activation of phosphoinositide 3-kinase (PI3K) is important for survival of transformed cells. Different computational approaches for this particular study are presented to calculate the solvation energies and binding free energies of H-Ras and PI3K. The goal of this study is to establish computational methods to investigate the roles of different proteins played in the cellular responses to space radiation, including modification of protein function through gene mutation, and to support the studies in molecular cell biology and theoretical kinetics models for our risk assessment project.

Ren, Lei↗
NASA NTRS2020

Hierarchical Coupling of Molecular Dynamics and Micromechanics to Predict the Elastic Properties of Three-Phase and Four-Phase Silicon Carbide Composites

The results obtained from previously conducted molecular dynamics analysis of silicon carbide (-SiC (6H, 4H, & 2H-SiC), -SiC (3C SiC)), silicon and boron nitride, were utilized as inputs in the MAC/GMC micromechanics software to model and evaluate the elastic properties of three-phase SiC/BN/SiC and four-phase SiC/BN/Si/SiC composites. This method of analysis eliminates the need for back-calculation of the apparent properties of the base constituents from the measured ceramic matrix composites properties. The multiscale models are validated against the available data in literature.

Aluko, Olanrewaju↗
NASA NTRS2005

Optimal Sampling of a Reaction Coordinate in Molecular Dynamics

Estimating how free energy changes with the state of a system is a central goal in applications of statistical mechanics to problems of chemical or biological interest. From these free energy changes it is possible, for example, to establish which states of the system are stable, what are their probabilities and how the equilibria between these states are influenced by external conditions. Free energies are also of great utility in determining kinetics of transitions between different states. A variety of methods have been developed to compute free energies of condensed phase systems. Here, I will focus on one class of methods - those that allow for calculating free energy changes along one or several generalized coordinates in the system, often called reaction coordinates or order parameters . Considering that in almost all cases of practical interest a significant computational effort is required to determine free energy changes along such coordinates it is hardly surprising that efficiencies of different methods are of great concern. In most cases, the main difficulty is associated with its shape along the reaction coordinate. If the free energy changes markedly along this coordinate Boltzmann sampling of its different values becomes highly non-uniform. This, in turn, may have considerable, detrimental effect on the performance of many methods for calculating free energies.

Pohorille, Andrew↗
NASA NTRS2012

Large Eddy Simulation (LES) of Particle-Laden Temporal Mixing Layers

High-fidelity models of plume-regolith interaction are difficult to develop because of the widely disparate flow conditions that exist in this process. The gas in the core of a rocket plume can often be modeled as a time-dependent, high-temperature, turbulent, reacting continuum flow. However, due to the vacuum conditions on the lunar surface, the mean molecular path in the outer parts of the plume is too long for the continuum assumption to remain valid. Molecular methods are better suited to model this region of the flow. Finally, granular and multiphase flow models must be employed to describe the dust and debris that are displaced from the surface, as well as how a crater is formed in the regolith. At present, standard commercial CFD (computational fluid dynamics) software is not capable of coupling each of these flow regimes to provide an accurate representation of this flow process, necessitating the development of custom software. This software solves the fluid-flow-governing equations in an Eulerian framework, coupled with the particle transport equations that are solved in a Lagrangian framework. It uses a fourth-order explicit Runge-Kutta scheme for temporal integration, an eighth-order central finite differencing scheme for spatial discretization. The non-linear terms in the governing equations are recast in cubic skew symmetric form to reduce aliasing error. The second derivative viscous terms are computed using eighth-order narrow stencils that provide better diffusion for the highest resolved wave numbers. A fourth-order Lagrange interpolation procedure is used to obtain gas-phase variable values at the particle locations.

Bellan, Josette↗
NASA NTRS2005

Laser-Based Flowfield Imaging in a Lean Premixed Prevaporized Sector Combustor

OH and fuel planar laser-induced fluorescence (PLIF) is used qualitatively in this study to observe the flame structure resultant from different fuel injector dome configurations within the 3-cup sector combustor test rig. The fluorescence images are compared with some computational fluid dynamics (CFD) results. Interferences in obtaining OH fluorescence signals due to the emission of other species are assessed. NO PLIF images are presented and compared to gas analysis results. The comparison shows that PLIF NO can be an excellent method for measuring NO in the flame. Additionally, we present flow visualization of the molecular species C2.

Hicks, Yolanda R.↗
NASA NTRS1973

Role of pressure transients in the detection and identification of lunar surface gas sources.

The dynamic behavior of neutral gases emitted from lunar surface sources is investigated. Pressure transients to be expected from several types of possible lunar surface sources are constructed. An approximate method is developed to estimate the source-sensor separation distance and the amount of effluent release. The case for which the source location and effluent molecular weight can be determined from ancillary data is discussed. It is shown that, if such data are available, pressure transients detected by either Apollo lunar surface experiment package ion spectrometers or cold cathode gages can be used to determine the release rate and the total amount of effluent release from the source.

Hall, F. G.↗
NASA NTRS2002

Modeling of Switching and Hysteresis in Molecular Transport

The conventional way of modeling current transport in two and three terminal molecular devices could be inadequate for certain cases involving switching and hysteresis. Here we present an alternate approach. Contrary to the regular way where applied bias directly modulates the conducting energy levels of the molecule, our method introduces a nonlinear potential energy surface varying with the applied bias as a control parameter. A time-dynamics is also introduced properly accounting for switching and hysteresis behavior. Although the model is phenomenological at this stage, we believe any detailed model would contain similar descriptions at its core.

Samanta, Manoj P.↗
NASA NTRS1990

Thermal instabilities in protogalactic clouds

The means by which a protogalaxy can fragment to form the first generation of stars and globular clusters remains an important problem in astrophysics. Gravitational instabilities grow on timescales too long to drive fragmentation before the background density grows by many orders of magnitude (see Murray and Lin 1989a, and references therein). Thermal instability provides a much more likely mechanism. After its initial collapse, a protogalactic cloud is expected to be shock heated to its virial temperature approx. 10(exp 6) K. Cooling by H and He+ below 10(exp 6) K has a negative slope, so that the cloud is subject to strong thermal instabilities. Density enhancements may then grow rapidly, fragmenting the protogalaxy as it cools to lower temperatures. The role of dynamical effects upon the growth of perturbations is considered here. The method used is similar to that used in Murray and Lin (1989a; see also the Erratum to appear September 15), which examined the growth of thermal instabilities with a one-dimensional Lagrangian hydrodynamics code, written for spherical symmetry. Perturbed regions therefore take the form of shells. The dynamical variables are integrated explicitly, while the temperature, ionization fraction, and molecular fraction are integrated implicitly, and account is taken for non-equilibrium values of these quantities.

Murray, S. D.↗
NASA NTRS1999

Molecular Modeling of the Poling of Piezoelectric Polyimides

The computational method described in this paper allows the calculation of the dielectric relaxation strength of an amorphous polymer based solely upon its chemical structure. The 4,4' oxydiphthalic anhydride (ODPA) dianhydride and bis-aminophenoxybenzene (APB) diamine based polyimides, (beta-CN) APB-ODPA and APB-ODPA were studied. Amorphous cells were constructed and then poled using molecular dynamics. Dielectric relaxation strengths of Delta(epsilon) = 17.8 for (beta-CN) APB-ODPA and Delta(epsilon) = 7.7 for APB-ODPA were predicted. These values are in excellent agreement with the experimental values. It was found that both the pendant nitrile dipole and the backbone anhydride residue dipole make significant contributions to the polyimides dielectric response. Specifically, it was shown that the difference in the magnitude of the dielectric relaxations is directly attributable to the nitrile dipole. The size of the relaxations indicate an absence of cooperative dipolar motions, The model was used to explain these results in terms of the average orientation of the nitrile and anhydride dipoles to within 51 deg. and 63 deg., respectively, of the applied electric field.

Young, J. A.↗
NASA NTRS1980

Modeling of a solar-pumped iodine laser

The direct conversion in space of solar radiation into laser radiation for power transmission to earth, satellites, or deep space probes shows promise as a reasonably simple technology and may have cost advantage in deployment and greater reliability compared to other methods of space power generation and transmission. The main candidates for solar pumping are the gas dynamic, photochemical, and direct photoexcited lasers. Here consideration is given to the photochemical reaction of alkyliodides which predominantly excite the I(2P1/2) state which then lases at 1.315 microns. The iodine ground state is eventually lost to reconstituting the gas or in the formation of molecular iodine. The rates at which the gas is required to be recycled through the laser system are modest. The side exposure at 100-fold solar concentration of a 100-m long tube with a 1 sq m cross section is estimated to provide 20 kW of continuous laser output. Scaling laws and optimum operating conditions of this system are discussed.

Wilson, J. W.↗
NASA NTRS1982

A numerical solution of the supersonic flow over a rearward facing step with transverse non-reacting hydrogen injection

The present investigation represents an application of computational fluid dynamics to a problem associated with the flow in the combustor region of a supersonic combustion ramjet engine (scramjet). The governing equations are considered, taking into account the Navier-Stokes equations, a molecular viscosity calculation, the molecular thermal conductivity, molecular diffusion, and a turbulence model. The employed numerical solution is patterned after the explicit, time-dependent, unsplit, predictor-corrector, finite-difference method given by MacCormack (1969). The calculation is concerned with the supersonic flow over a rearward-facing step with transverse H2 injection at conditions germane to the combustor region of a scramjet engine. The H2 jet acts as an effective body which essentially shields the primary flow from the rearward-facing step, thus substantially changing the wave pattern in the primary flow.

Berman, H. A.↗
NASA NTRS1943

A Method for Removal of Bakelite-Impregnated Wire Strain Gages

An increasing interest is being shown in the use of wire-type strain gages to measure static and dynamic stresses in aircraft-engine parts. Bakelite cement has been found satisfactory as a bonding agent for attaching the strain gages to machine parts that must operate at elevated temperatures. On many occasions, it is desired to remove Bakelite-cemented gages from the test parts for the purpose of replacing faulty gages or of returning the parts to service after strain measurements have been completed. Removal of the gages by means of scraping without prior treatment is very unsatisfactory because it is tedious and almost invariably damages the finished surface. Various solvents have been tried, but all attempts in this direction have been unsuccessful inasmuch as Bakelite cement, when properly baked, forms a polymer of very high molecular weight that resists the action of solvents. This report presents a gas-flame method of removal that is rapid and does not inure the structural part.

Kemp, Richard H.↗
NASA NTRS1990

Plume flowfield analysis of the shuttle primary Reaction Control System (RCS) rocket engine

A solution was generated for the physical properties of the Shuttle RCS 4000 N (900 lb) rocket engine exhaust plume flowfield. The modeled exhaust gas consists of the five most abundant molecular species, H2, N2, H2O, CO, and CO2. The solution is for a bare RCS engine firing into a vacuum; the only additional hardware surface in the flowfield is a cylinder (=engine mount) which coincides with the nozzle lip outer corner at X = 0, extends to the flowfield outer boundary at X = -137 m and is coaxial with the negative symmetry axis. Continuum gas dynamic methods and the Direct Simulation Monte Carlo (DSMC) method were combined in an iterative procedure to produce a selfconsistent solution. Continuum methods were used in the RCS nozzle and in the plume as far as the P = 0.03 breakdown contour; the DSMC method was used downstream of this continuum flow boundary. The DSMC flowfield extends beyond 100 m from the nozzle exit and thus the solution includes the farfield flow properties, but substantial information is developed on lip flow dynamics and thus results are also presented for the flow properties in the vicinity of the nozzle lip.

Hueser, J. E.↗