Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Dynamic Simulation”

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 811 records · Page 45

Computational Analysis of the Combustion Processes in an Axisymmetric, RBCC Flowpath

Computational fluid dynamic simulations have been used to study the combustion processes within an axisymmetric, RBCC flowpath. Two distinct operating modes have been analyzed to date, including the independent ramjet stream (IRS) cycle and the supersonic combustion ramjet (scramJet) cycle. The IRS cycle investigation examined the influence of fuel-air ratio, fuel distribution, and rocket chamber pressure upon the combustion physics and thermal choke characteristics. Results indicate that adjustment of the amount and radial distribution of fuel can control the thermal choke point. The secondary massflow rate was very sensitive to the fuel-air ratio and the rocket chamber pressure. The scramjet investigation examined the influence of fuel-air ratio and fuel injection schedule upon combustion performance estimates. An analysis of the mesh-dependence of these calculations was presented. Jet penetration data was extracted from the three-dimensional simulations and compared favorably with experimental correlations of similar flows. Results indicate that combustion efficiency was very sensitive to the fuel schedule.

Steffen, Christopher J., Jr.↗

Generic Wing-Body Aerodynamics Data Base

The wing-body aerodynamics data base consists of a series of CFD (Computational Fluid Dynamics) simulations about a generic wing body configuration consisting of a ogive-circular-cylinder fuselage and a simple symmetric wing mid-mounted on the fuselage. Solutions have been obtained for Nonlinear Potential (P), Euler (E) and Navier-Stokes (N) solvers over a range of subsonic and transonic Mach numbers and angles of attack. In addition, each solution has been computed on a series of grids, coarse, medium and fine to permit an assessment of grid refinement errors.

Holst, Terry L.↗

Visualizing Flow Over Parametric Surfaces Using Line Integral Convolution

Line Integral Convolution (LIC) is a powerful technique for imaging and animating vector fields. We extend the LIC paradigm in three ways: (1) The existing technique is limited to vector fields over a regular Cartesian grid. We extend it to vector fields over parametric surfaces, such as those found in curvilinear grids, used in computational fluid dynamics simulations; (2) Periodic motion filters can be used to animate the flow visualization. When the flow lies on a parametric surface, however, the motion appears misleading. We explain why this problem arises and show how to adjust the LIC algorithm to handle it; (3) We introduce a technique to visualize vector magnitudes as well as vector direction. Cabral and Leedom have suggested a method for variable-speed animation, which is based on varying the frequency of the filter function. We develop a different technique based on kernel phase shifts which we have found to show substantially better results. Our implementation of these algorithms utilizes texture-mapping hardware to run in real time, which allows them to be included in interactive applications.

Forssell, Lisa↗

Automation of a N-S S and C Database Generation for the Harrier in Ground Effect

A method of automating the generation of a time-dependent, Navier-Stokes static stability and control database for the Harrier aircraft in ground effect is outlined. Reusable, lightweight components arc described which allow different facets of the computational fluid dynamic simulation process to utilize a consistent interface to a remote database. These components also allow changes and customizations to easily be facilitated into the solution process to enhance performance, without relying upon third-party support. An analysis of the multi-level parallel solver OVERFLOW-MLP is presented, and the results indicate that it is feasible to utilize large numbers of processors (= 100) even with a grid system with relatively small number of cells (= 10(exp 6)). A more detailed discussion of the simulation process, as well as refined data for the scaling of the OVERFLOW-MLP flow solver will be included in the full paper.

Murman, Scott M.↗

Spaceborne Microwave Instrument for High Resolution Remote Sensing of the Earth's Surface Using a Large-Aperture Mesh Antenna

This report describes a two-year study of a large-aperture, lightweight, deployable mesh antenna system for radiometer and radar remote sensing of the Earth from space. The study focused specifically on an instrument to measure ocean salinity and Soil moisture. Measurements of ocean salinity and soil moisture are of critical . importance in improving knowledge and prediction of key ocean and land surface processes, but are not currently obtainable from space. A mission using this instrument would be the first demonstration of deployable mesh antenna technology for remote sensing and could lead to potential applications in other remote sensing disciplines that require high spatial resolution measurements. The study concept features a rotating 6-m-diameter deployable mesh antenna, with radiometer and radar sensors, to measure microwave emission and backscatter from the Earth's surface. The sensors operate at L and S bands, with multiple polarizations and a constant look angle, scanning across a wide swath. The study included detailed analyses of science requirements, reflector and feedhorn design and performance, microwave emissivity measurements of mesh samples, design and test of lightweight radar electronic~., launch vehicle accommodations, rotational dynamics simulations, and an analysis of attitude control issues associated with the antenna and spacecraft, The goal of the study was to advance the technology readiness of the overall concept to a level appropriate for an Earth science emission.

Njoku, E.↗

elVis: An Interactive System For Visualization of Unsteady Fluid Flow

ElVis is a prototype system with allows for the interactive visualization of unsteady fluid flow. The increasing computational power applied to fluid dynamics simulations presents the enormous challenge to the visualization system designer to apply a wide range of technologies to the analysis process with ever increasing demands on performance. Visualization of the results of unsteady fluid flow simulations presents the challenge of exploring very large and complex data sets. Since exploration is a trial and error process, it is of utmost importance that the time required to execute a trial (i.e., create a visualization) be at a minimum in order to provide real time interaction.

Gerald-Yamasaki, Michael↗

Structure and Dynamics of the Quiet Solar Chromosphere

The grant supported research on the structure of the quiet, nonmagnetic chromosphere and on wave excitation and propagation in both the nonmagnetic chromosphere and the magnetic network. The work on the structure of the chromosphere culminated in the recognition that between two competing views of the solar chromosphere, older models by Avrett and collaborators (referred to as VAL) and the newer, dynamical model by Carlsson & Stein (referred to as CS), the clear decision is in favor of the older models, and this in spite of the evident lack of physics, which does not include wave motion and oscillations. The contrast between the static VAL models and the dynamical CS model can be stated most succinctly by comparing the temperature variation implied by the VAL models and the temperature fluctuations of the CS model, which are, respectively, of the order of 10% for the VAL model (at heights where hydrogen is 50% ionized) and a factor of 10 (at the upper boundary of their chromospheric model). The huge fluctuations of the CS model have never been observed, whereas the smaller temperature variations of the VAL models are consistent with ground-based and space-based observations. While it should be obvious which model describes the Sun and which one fails, the case is far from settled in the minds of solar physicists. Thus, much educational work remains to be done and, of course, more research to develop arguments that make the case more convincing. The research on waves and oscillations has been based on a unified theory of excitation of acoustic waves in the field-free atmosphere and of transverse and longitudinal waves in magnetic flux tubes located in the magnetic network by noting, first, that impulsive excitation of all these waves in gravitationally stratified media leads to oscillations at the respective cutoff frequencies and, second, that the observed oscillation frequencies in the nonmagnetic and magnetic parts of the chromosphere match corresponding cutoff frequencies in the upper solar photosphere. The dynamical simulations by Carlsson and Stein have been most instructive and of fundamental importance for understanding wave propagation in a stratified medium by their "flaws", the most important of which is an intensity excess at the H2v emission peak in the H line of Ca II, which is surprising since the observed intensity should have been an upper limit to the simulated intensity. The only plausible explanation for a predicted intensity that is higher than observed is that energy is spread horizontally in upward propagation in the Sun, but not in the plane-wave modeling of CS (as well as by almost everyone else). Investigations of the horizontal size of the region disturbed by the upward-propagating shock in the acoustic-wave propagation implies that the waves in H2v bright-point oscillations emanate from a point source with a diameter corresponding to the width of an intergranular lane, about 100 km, and reach a size of about 4000 km in the upper layers of the chromosphere. Linear, analytic modeling of waves emanating from a point source in a stratified atmosphere shows that the upward-expanding propagation channel does not have sharp boundaries and that the shape of the wave front depends on the order of the wave behind the initial pulse. Otherwise, the behavior of the linear waves resembles that of the nonlinear shock waves observed in the Sun. Research that needs to be done to firm up the conclusions reached above concerns the numerical simulations of nonlinear waves and oscillations in a three-dimensional stratified atmosphere with impulsive excitation, and observations linking directly the horizontal size of the disturbed area in upward propagation to individual waves.

Kalkofen, Wolfgang↗

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↗

Comparison of Oxygen Gauche Effects in Poly(Oxyethylene) and Poly(ethylene terephtylene) Based on Quantum Chemistry Calculations

The so-called oxygen gauche effect in poly(oxyethylene) (POE) and its model molecules such as 1,2-dimethoxyethane (DME) and diglyme (CH3OC2H4OC2H4OCH3) is manifested in the preference for gauche C-C bond conformations over trans. This has also been observed for poly(ethylene terephthalate) (PET). Our previous quantum chemistry calculations demonstrated that the large C-C gauche population in DME is due, in part, to a low-lying tg +/- g+ conformer that exhibits a substantial 1,5 CH ... O attraction. New calculations will be described that demonstrate the accuracy of the original quantum chemistry calculations. In addition, an extension of this work to model molecules for PET will be presented. It is seen that the C-C gauche preference is much stronger in 1,2 diacetoxyethane than in DME. In addition, there exist low-lying tg +/- g+/- and g+/-g+/-g+/- conformers that exhibit 1,5 CH ... O attractions involving the carbonyl oxygens. It is expected that the -O-C-C-O- torsional properties will be quite different in these two polymers. The quantum chemistry results are used to parameterize rotational isomeric states models (RIS) and force fields for molecular dynamics simulations of these polymers.

Jaffe, Richard↗

Origin of Water on Mars

Dynamical simulations suggest that the Earth's water budget was delivered primarily from the asteroid belt, in the form of large planetary embryos. The same simulations present a very different picture for Mars its water came from a mixture of cometary and small asteroidal bodies. Additional information is contained in the original extended abstract.

Lunine, J. I.↗

Observation and Modeling of Single Wall Carbon Nanotube Bend Junctions

Single wall carbon nanotube (SWNT) bends, with diameters from approx. 1.0 to 2.5 nm and bend angles from 18 deg. to 34 deg., are observed in catalytic decomposition of hydrocarbons at 600 - 1200 C. An algorithm using molecular dynamics simulation (MD) techniques is developed to model these structures that are considered to be SWNT junctions formed by topological defects (i.e. pentagon-heptagon pairs). The algorithm is used to predict the tube helicities and defect configurations for bend junctions using the observed tube diameters and bend angles. The number and arrangement of the defects at the junction interfaces are found to depend on the tube helicities and bend angle. The structural and energetic calculations using the Brenner potential show a number of stable junction configurations for each bend angle with the 34 deg. bends being more stable than the others. Tight binding calculations for local density of state (LDOS) and transmission coefficients are carried out to investigate electrical properties of the bend junctions.

Han, Jie↗

Electronic Properties of Carbon Nanotubes and Junctions

Metallic and semiconducting Single Wall Carbon Nanotubes (CNT) have recently been characterized using scanning tunneling microscopy (STM) and the manipulation of individual CNT has been demonstrated. These developments make the prospect of using CNT as molecular wires and possibly as electronic devices an even more interesting one. We have been modeling various electronic properties such as the density of states and the transmission coefficient of CNT wires and junctions. These studies involve first calculating the stability of junctions using molecular dynamics simulations and then calculating the electronic properties using a pi-electron tight binding Hamiltonian. We have developed the expertise to calculate the electronic properties of both finite-sized CNT and CNT systems with semi-infinite boundary conditions. In this poster, we will present an overview of some of our results. The electronic application of CNT that is most promising at this time is their use as molecular wires. The conductance can however be greatly reduced because of reflection due to defects and contacts. We have modeled the transmission through CNT in the presence of two types of defects: weak uniform disorder and strong isolated scatterers. We find that the conductance is affected in significantly different manners due to these defects Junctions of CNT have also been imaged using STM. This makes it essential to derive rules for the formation of junctions between tubes of different chirality, study their relative energies and electronic properties. We have generalized the rules for connecting two different CNT and have calculated the transmission and density of states through CNT junctions. Metallic and semiconducting CNT can be joined to form a stable junction and their current versus voltage characteristics are asymmetric. CNT are deformed by the application of external forces including interactions with a substrate or other CNT. In many experiments, these deformation are expected to occur naturally. We will present some preliminary results of our calculations of the modification of CNT electronic properties as a result of deformations.

Anantram, M. P.↗

How Human Factors Drove the Design and Implementation of the Virtual Windtunnel

This viewgraph presentation describes decisions made at the NASA Ames Research Center during its development of a virtual windtunnel to assist humans in visualizing computational fluid dynamics simulations. User requirements for the system include the simulation of vortical structure, pressure distribution, and overall sense of flow. In addition, the system needs to support a variety of interfaces, excluding head mounts, and needs to use an object oriented approach. Direct manipulation of the system is most user-friendly when limited to only grab and point gestures. 'Visualization control tools' (vtools) improve the realism of the system. The necessary object oriented programming is in C++ and openGL. Users interact with objects called 'tools', which include vtools and tools which control the virtual wind tunnel environment. Other objects include data objects accessed by the visualizations. Visualizations can be added to the system by the user. The presentation includes a discussion of run-time architecture, and issues related to computation and implementation.

Bryson, Steve↗

Nanoscale Etching and Indentation of Silicon(001) Surface with Carbon Nanotube Tips

The possibility of nanoscale etching and indentation of Si(001)(2x1) surface by (8,0) and (10,10) carbon nanotube tips is demonstrated, for the first time, by classical molecular dynamics simulations employing Tersoff's many-body potential for a mixed C/Si/Ge system. In the nanotube tip barely touching the surface scenario atomistic etching is observed, where as in the nanoindentation scenario nanotube tip penetrates the surface without much hindrance. The results are explained in terms of the relative strength of C-C, C-Si, and Si-Si bonds.

Dzegilenko, Fendor N.↗

Nanomechanics of Carbon and CxByNz Nanotubes: Via a Quantum Molecular Dynamics Method

Nanomechanics of single-wall C, BN and BC$_3$ and B doped C nanotubes under axial compression and tension are investigated through a generalized tight-binding molecular dynamics (GTBMD) and {\it ab-initio} electronic structure methods. The dynamic strength of BN, BC$_3$ and B doped C nanotubes for small axial strain are comparable to each other. The main difference is in the critical strain at which structural collapse occurs. For example, even a shallow doping with B lowers the value of critical strain for C nanotubes. The critical strain for BN nanotube is found to be more than that for the similar C nanotube. Once the structural collapse starts to occur we find that carbon nanotubes irreversibly go into plastic deformation regime via the formation of tetrahedral (four-fold coordinated) bonds at the location of sharp pinches or kinks. This finding is considerably different from the classical MD (molecular dynamics) simulation results known so far. The energetics and electronic densities of states of the collapsed structures, investigated with {\it ab-initio) methods, will also be discussed.

Srivastava, Deepak↗

Extreme Velocity Wind Sensor

This presentation provides an overview of the development of new hurricane wind sensor (Extreme Velocity Wind Sensor) for the Kennedy Space Center (KSC) which is designed to withstand winds of up to three hundred miles an hour. The proposed Extreme Velocity Wind Sensor contains no moveable components that would be exposed to extreme wind conditions. Topics covered include: need for new hurricane wind sensor, conceptual design, software applications, computational fluid dynamic simulations of design concept, preliminary performance tests, and project status.

Perotti, Jose↗

Motion of Drops on Surfaces with Wettability Gradients

A liquid drop present on a solid surface can move because of a gradient in wettability along the surface, as manifested by a gradient in the contact angle. The contact angle at a given point on the contact line between a solid and a liquid in a gaseous medium is the angle between the tangent planes to the liquid and the solid surfaces at that point and is measured within the liquid side, by convention. The motion of the drop occurs in the direction of increasing wettability. The cause of the motion is the net force exerted on the drop by the solid surface because of the variation of the contact angle around the periphery. This force causes acceleration of an initially stationary drop, and leads to its motion in the direction of decreasing contact angle. The nature of the motion is determined by the balance between the motivating force and the resisting hydrodynamic force from the solid surface and the surrounding gaseous medium. A wettability gradient can be chemically induced as shown by Chaudhury and Whitesides who provided unambiguous experimental evidence that drops can move in such gradients. The phenomenon can be important in heat transfer applications in low gravity, such as when condensation occurs on a surface. Daniel et al have demonstrated that the velocity of a drop on a surface due to a wettability gradient in the presence of condensation can be more than two orders of magnitude larger than that observed in the absence of condensation. In the present research program, we have begun to study the motion of a drop in a wettability gradient systematically using a model system. Our initial efforts will be restricted to a system in which no condensation occurs. The experiments are performed as follows. First, a rectangular strip of approximate dimensions 10 x 20 mm is cut out of a silicon wafer. The strip is cleaned thoroughly and its surface is exposed to the vapor from an alkylchlorosilane for a period lasting between one and two minutes inside a desiccator. This is done using an approximate line source of the vapor in the form of a string soaked in the alkylchlorosilane. Ordinarily, many fluids, including water, wet the surface of silicon quite well. This means that the contact angle is small. But the silanized surface resists wetting, with contact angles that are as large as 100 degs. Therefore, a gradient of wettability is formed on the silicon surface. The region near the string is highly hydrophobic, and the contact angle decreases gradually toward a small value at the hydrophilic end away from this region. The change in wettability occurs over a distance of several mm. The strip is placed on a platform within a Plexiglas cell. Drops of a suitable liquid are introduced on top of the strip near the hydrophobic end. An optical system attached to a video camera is trained on the drop so that images of the moving drop can be captured on videotape for subsequent analysis. We have performed preliminary experiments with water as well as ethylene glycol drops. Results from these experiments will be presented in the poster. Future plans include the refinement of the experimental system so as to permit images to be recorded from the side as well as the top, and the conduct of a systematic study in which the drop size is varied over a good range. Experiments will be conducted with different fluids so as to obtain the largest possible range of suitably defined Reynolds and Capillary numbers. Also, an effort will be initiated on theoretical modeling of this motion. The challenges in the development of the theoretical description lie in the proper analysis of the region in the vicinity of the contact line, as well as in the free boundary nature of the problem. It is known that continuum models assuming the no slip condition all the way to the contact line fail by predicting that the stress on the solid surface becomes singular as the contact line is approached. One approach for dealing with this issue has been to relax the no-slip boundary condition using the Navier model. Molecular dynamics simulations of the contact line region show that for a non-polar liquid on a solid surface, the no-slip boundary condition is in fact incorrect near the contact line. Furthermore, the same simulations also show that the usual relationship between stress and the rate of deformation breaks down in the vicinity of the contact line. In developing continuum theoretical models of the system, we shall accommodate this knowledge to the extent possible.

Subramanian, R. Shankar↗