SEARCH · Engineering Papers

Results for “continuum scale modeling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Time-resolved spectrophotometry of the emission lines in the galactic X-ray source H2252-035

The galactic X-ray source H2252-035 has been the object of a concentrated observational effort, since Griffiths et al. (1980) have identified it with a cataclysmic variable-like spectrum. Since the observed 859 s periodicity corresponds to the beat frequency between the X-ray and orbital periods, it has been interpreted as related to a reprocessing of the X-ray flux into optical emission in some region of the system which is rotating with the orbital period and heated by beamed X-ray emission from the compact star. In connection with the present investigation, H2252-035 has been observed with a 120 inch (3 m) telescope and image dissector scanner on 1980 September 9/10 and on 1981 September 20/21. At a time resolution of approximately 2 minutes, four bright emission lines showed short time-scale flickering and a pronounced orbital period variation in phase with the continuum modulation. A two-component model for the line emission region is proposed.

Clarke, J. T.↗

NASA NTRS2023

Remote Sensing of Dissolved Organic Carbon (DOC) Stocks, Fluxes and Transformations Along the Land-Ocean Aquatic Continuum: Advances, Challenges, and Opportunities

From stream headwaters to the deep ocean, dissolved organic carbon (DOC) is ubiquitous across the world’s aquatic systems, where it fulfills important biogeochemical and ecological functions. Importantly, DOC is a dynamic component of the carbon cycle and represents a global reservoir (mass of carbon) equivalent in size to the atmospheric CO 2 pool. The mobilization, transport, and transformations of DOC along the land-ocean aquatic continuum biogeochemically connect terrestrial landscapes, freshwater and marine ecosystems, and the atmosphere. As a result, DOC budgets are vulnerable to climate- and human-driven changes. Having the ability to quantify DOC stocks, fluxes, and transformations consistently and accurately in aquatic systems is therefore critical to our understanding of DOC’s functions in the environment and of their susceptibility to change. In recent years, optical techniques and synoptic earth observations have opened new doors to achieve this objective on a global scale. In this review, we explore how in situ measurements, earth observations, and the models they inform have facilitated the large-scale quantification of DOC stocks, fluxes, and transformations along the land-ocean aquatic continuum and how they have helped shape our current understanding of the global DOC cycle. The review also explores the limitations and challenges we still face today, and discusses how upcoming and future sensors are expected to provide new insights into our understanding of the global DOC cycle.

Dissolved organic carbon↗

NASA NTRS2006

Precision Fe K-Alpha and Fe K-Beta Line Spectroscopy of the Seyfert 1.9 Galaxy NGC 2992 with Suzaku

We present detailed time-averaged X-ray spectroscopy in the 0.5-10 keV band of the Seyfert 1.9 galaxy NGC 2992 with the Suzaku X-ray Imaging Spectrometers (XIS). The source had a factor approximately 3 higher 2-10 keV flux (approximately 1.2 x l0(exp -11) erg per square cm per s) than the historical minimum and a factor approximately 7 less than the historical maximum. The XIS spectrum of NGC 2992 can be described by several components. There is a primary continuum, modeled as a power-law with a photon index of Gamma = 1.57(sup +0.06) (sup -0.03) that is obscured by a Compton-thin absorber with a column density of 8.01(sup +0.6) (sup -0.5)x l0 (exp 21) per square cm. . There is another, weaker, unabsorbed power-law component (modeled with the same slope as the primary), that is likely to be due to the primary continuum being electron-scattered into our line-of-sight by a region extended on a scale of hundreds of parsecs. We measure the Thomson depth of the scattering zone to be Tau = 0.072 +/- 0.021. An optically-thin thermal continuum emission component, which probably originates in the same extended region, is included in the model and yields a temperature and luminosity of KT = 0.656(sup +0.088) (sup -0.0.61) keV and approximately 1.2 +/- 0.4 x l0 (exp 40) erg per s respectively. We detect an Fe K emission complex which we model with broad and narrow lines and we show that the intensities of the two components are decoupled at a confidence level > 3 sigma. The broad Fe K alpha line has an equivalent width of 118(sup +32) (sup -61) eV and could originate in an accretion disk (with inclination angle greater than approximately 30 deg) around the putative central black hole. The narrow Fe K alpha line has an equivalent width of 1632(sup +47) (sup -26) eV and is unresolved (FWHM < 4630 km per s) and likely originates in distant matter. The absolute flux in the narrow line implies that the column density out of the line-of-sight could be much higher than measured in the line-of-sight, and that the mean (historically-averaged) continuum luminosity responsible for forming the line could be a factor of several higher than that measured from the data. We also detect the Fe K Beta line (corresponding to the narrow Fe K alpha line) with a high signal-to-noise ratio and describe a new robust method to constrain the ionization state of Fe responsible for the Fe K alpha and Fe K Beta lines that does not require any knowledge of possible gravitational and Doppler energy shifts affecting the line energies. For the distant line-emitting matter (e. g. the putative obscuring torus) we deduce that the predominant ionization state is lower than Fe VIII (at 99% confidence), conservatively taking into account residual calibration uncertainties in the XIS energy scale and theoretical and experimental uncertainties in the Fe K fluorescent line energies. From the limits on a possible Compton-reflection continuum it is likely that the narrow Fe K alpha and Fe K Beta lines originate in a Compton-thin structure.

Yaqoob, Tahir↗

NASA NTRS1983

Modeling global structural damping in trusses using simple continuum models

Truss beams with members having viscous damping are modeled as continuum Timoshenko beams. Procedures for deriving the equivalent beam stiffnesses and damping are presented. The global damping for the continuum beam is explicitly expressed in terms of the damping coefficients of the individual truss members. The continuum beam model is used to study transient vibration problems and the solutions are compared well with the full scale finite element solutions. The gradient method is used for parameter estimations in conjunction with the Timoshenko beam model. It is shown that the Timoshenko beam model can be easily updated with measured data and the updated model can yield very accurate transient solutions.

Sun, C. T.↗

NASA NTRS2002

Thermomechanical Response Variability of Stiffened Composite Panels

A significant numerical simulation capability now exists for studying the various phenomena associated with the response, failure, and performance of multilayered composite panels and shells subjected to combined pressure, mechanical, and thermal loads. The phenomena involved cover a wide range of length scales from local to global structural response. The modeling approaches used for multilayered panels include micromechanical models, three-dimensional continuum models, quasi-three-dimensional models, and two-dimensional plate and shell models. Within each category a number of models with several levels of sophistication has evolved. The four categories are described in review papers. Despite the extensive literature cited in the afore-mentioned references, only a few studies have been reported on the effects of stiffness discontinuities, such as those associated with an abrupt stiffener termination or dropped plies, on the response of composite panel. Stiffener termination is often necessary in composite aerospace structures to satisfy detailed design requirements, and therefore an understanding and a prediction of its effect on the response and failure of composite panels are desirable. Such a prediction must take into account the fact that current measurement technology does not allow the accurate determination of the material parameters that are used in the analytical models.

Noor, Ahmed K.↗

NASA NTRS2018

Plasticity and Damage Modeling of Stress Asymmetry and Dynamic Behavior of AFS Additive Manufactured Aluminum Alloy 2219

The Solid State Additive Manufacturing (AM) process referred as MELD that fabricated the samples in this study, provides a new path for repairing, coating, joining and additive manufacturing metals and metal matrix composites. This research will be the first application of a physics-based microstructure dependent internal state variable (ISV) plasticity and damage material model to capture the mechanical response of an AM Aluminum Alloy (AA) 2219 via the MELD process. In this research, a microstructure-based internal state variable (ISV) plasticity-damage model was used to capture the mechanical behavior of AFS 2219 aluminum alloy. Aeroprobe Corporation, creator and patent holder for the MELD process, fabricated the material by pushing a solid filler rod of AA2219-T861 material through a hollow rotating tool onto an AA2219 T851 plate substrate. As feedstock, solid or powder precursor metals are pushed through a nonconsumable rotating cylindrical tool. Herein, added layers are deposited and metallurgically bonded to substrate material or previously deposited layers by the heat generated from the rotating tool through plastic deformation of the filler material. Once a layer has been added, the tool height increases, and starts the deposition of the next layer. This process results in beneficial properties such as grain refinement, homogenization and reduced porosity (fully dense). This process will experience temperatures similar to those in the weld nugget zone (WNZ) in friction stir welding (FSW), ranging from 0.6-0.9 Tm, with Tm being the melting point of the material. MELD is highly scalable with AA deposition rates reaching over 1000 cm3/hr, which allows for MELD being used for repairs, coatings, and building components. A motivating factor driving the research for physics-based history dependent material modeling of MELD components is the ability to accurately capture the stress-state and strain rate dependence in the material caused by variations in material microstructure from the MELD processing of new or repaired components. The ISV model incorporates microstructural content and is consistent with continuum level kinematics, kinetics, and thermodynamics. These features allow the ISV model to capture large deformations at the structural scale using the kinematic and isotropic hardening, while microscale damage is obtained from the microstructural features. The benefits of the ISV model arise from the inclusion of structure-property relationships identified from microstructural characterization and experimentation. The Bauschinger effect (BE) is an important concept, vital in the accurate prediction of cyclic stress-strain response of ductile materials such as metals. The ISV model has been successfully used to capture the behavior and damage, and the BE of different aluminum alloys and steels. The ISV model uses kinematic and isotropic hardening to help capture deformations of the material at the macro scale. To understand this hardening relationship, calculating the kinematic and isotropic hardening relationship in the material is warranted for a high-fidelity model. Electron Backscattered Diffraction (EBSD) was used to characterize the as-fabricated microstructure, where a fully-dense equiaxed grain morphology with average grain size of 2.5 μm was observed. Microhardness mapping of the as-built structures, monotonic tension and compression experiments at both quasi-static (0.001/s) strain rates, tension-followed-by-compression and compression-followed-by-tension experiments were performed to obtain the set of plasticity and damage constants necessary to capture strain rate and stress state behavior of this additive material. To calibrate the plasticity-damage model, a single set of constants were determined to capture the different stress states the MELD AA2219. One set of the constants was determined from experimental true stress-strain curves for the tension and compression data. Additionally, microstructural information and data from the open literature were used as the other model constants. This research is a first of its kind for AFS AA2219, includes correlating the ISV model to the monotonic experimental results that capture the isotropic and kinematic plasticity mechanical response.

Rivera, O. G.↗

NASA NTRS2018

Validation of a Mesoscale Fiber Kinking Model Through Test and Analysis of Double Edge Notch Compression Specimens

Tests and analyses were conducted on a series of geometrically-scaled double edge notch compression specimens to validate the capability of the fiber-kinking model in the NASA continuum damage mechanics code for progressive damage analysis, CompDam. The tests and analyses focused on the fiber-kinking damage mechanism, which is of critical importance for predicting accurately longitudinal compression failure in carbon fiber reinforced polymer laminates. The fiber-kinking model in CompDam was augmented with a new capability to represent a hardening response once a kink band has fully formed, enabling prediction of band broadening. Correlation of the test and analysis results showed agreement in stiffness and strength with less than 10% error. The analysis predicts the same sequence of events leading to ultimate failure that was found in the test. The overall excellent correlation in terms of stiffness, strength, and failure process validates the capability of the model for predicting longitudinal compression failure in notched laminates with cross-ply layups.

Bergan, Andrew C.↗

NASA NTRS1978

A comparison of synthetic and measured solar continuum intensities and limb darkening coefficients

The paper compares recent center-limb studies of optical and IR continuum intensities (over the wavelength range 0.3-2.4 microns, with the data tabulated in the form of second- and fifth-order polynomial fits to empirical limb darkening curves) by Pierce and Slaughter (1977) and by Pierce et al. (1977) with predictions based on several current photospheric thermal structure models, using calibrated specific intensity measurements by Labs and Neckel (1968, 1970) to set an absolute temperature scale. The thermal structure in best quantitative agreement with the measurements of Pierce and Slaughter and Pierce et al. is the semiempirical model of Vernazza, Avrett and Loeser (1976). The temperatures in this model must be scaled upward, however, by a factor of 1.015 plus or minus 0.0005 to be consistent with the Labs and Neckel absolute calibration of continuum high points in the 0.40-0.65 micron region.

Ayres, T. R.↗

NASA NTRS1980

A mathematical model of post-instability in continuum mechanics

The post-instability can appear in the form of failure of hyperbolicity or in the form of cascade instability (stretching of vorticity in an inviscid fluid or instability of the Navier-Stokes equations). All the cases are characterized by an unlimited decrease of the scale of the motions, in the course of which the derivatives of the corresponding functions tend to infinity, though the functions themselves remain finite. Such an instability shows that the corresponding state of a continuum cannot be properly simulated by smooth functions. Hence, the original model must be corrected by giving up the requirement about differentiability and enlarging the class of functions. Applications to a model of an inviscid fluid and to a model of turbulence are discussed.

Zak, M.↗

NASA NTRS1988

Temperature determinations of hot DA white dwarfs using IUE continuum fluxes

Effective temperatures of 15 DA white dwarfs hotter than 20,000 K were derived from low-dispersion far ultraviolet spectra obtained with IUE. The analysis was carried out by comparing the observed far ultraviolet fluxes with model fluxes scaled to the V-band flux. Accurate calibration of the IUE spectra is critical for this analysis. Observations at all epochs were corrected to the 1980 IUE calibration using the time-dependent corrections of Bohlin (1988). Taking advantage of the smooth and well-defined continuum fluxes provided by DA white dwarfs, seven white dwarfs for which accurate, independent temperature determinations were made from line profile analyses were used to improve the accuracy of the IUE flux calibration. The correction to the original calibration is as great as 20 percent in individual 5 A wavelength bins, while the average over the IUE wavelength range is 5 percent. The final calibration correction and the temperatures for the hot white dwarfs are presented.

Finley, David↗

NASA NTRS2005

Multi-Scale Modeling of Cross-Linked Nanotube Materials

The effect of cross-linking single-walled carbon nanotubes on the Young's modulus of a nanotube-reinforced composite is modeled with a multi-scale method. The Young's modulus is predicted as a function of nanotube volume fraction and cross-link density. In this method, the constitutive properties of molecular representative volume elements are determined using molecular dynamics simulation and equivalent-continuum modeling. The Young's modulus is subsequently calculated for cross-linked nanotubes in a matrix which consists of the unreacted cross-linking agent. Two different cross-linking agents are used in this study, one that is short and rigid (Molecule A), and one that is long and flexible (Molecule B). Direct comparisons between the predicted elastic constants are made for the models in which the nanotubes are either covalently bonded or not chemically bonded to the cross-linking agent. At a nanotube volume fraction of 10%, the Young's modulus of Material A is not affected by nanotube crosslinking, while the Young's modulus of Material B is reduced by 64% when the nanotubes are cross-linked relative to the non-cross-linked material with the same matrix.

Frankland, S. J. V.↗

NASA NTRS2012

Integrated Modeling of Aerosol, Cloud, Precipitation and Land Processes at Satellite-Resolved Scales

In this talk, I will present recent results from a project led at NASA/GSFC, in collaboration with NASA/MSFC and JHU, focused on the development and application of an observation-driven integrated modeling system that represents aerosol, cloud, precipitation and land processes at satellite-resolved scales. The project, known as the NASA Unified WRF (NU-WRF), is funded by NASA's Modeling and Analysis Program, and leverages prior investments from the Air Force Weather Agency and NASA's Earth Science Technology Office (ESTO). We define "satellite-resolved" scales as being within a typical mesoscale atmospheric modeling grid (roughly 1-25 km), although this work is designed to bridge the continuum between local (microscale), regional (mesoscale) and global (synoptic) processes. NU-WRF is a superset of the standard NCAR Advanced Research WRF model, achieved by fully integrating the GSFC Land Information System (LIS, already coupled to WRF), the WRF/Chem enabled version of the Goddard Chemistry Aerosols Radiation Transport (GOCART) model, the Goddard Satellite Data Simulation Unit (SDSU), and boundary/initial condition preprocessors for MERRA and GEOS-5 into a single software release (with source code available by agreement with NASA/GSFC). I will show examples where the full coupling between aerosol, cloud, precipitation and land processes is critical for predicting local, regional, and global water and energy cycles, including some high-impact phenomena such as floods, hurricanes, mesoscale convective systems, droughts, and monsoons.

Peters-Lidard, Christa↗

NASA NTRS2021

From Pseudostreamer Jets to Coronal Mass Ejections:: Observations of the Breakout Continuum

The magnetic breakout model, in which reconnection in the corona leads to destabilization of a filament channel, explains numerous features of eruptive solar events, from small-scale jets to global-scale coronal mass ejections (CMEs). The underlying multipolar topology, pre-eruption activities, and sequence of magnetic-reconnection onsets (first breakout, then flare) of many observed fast CMEs/eruptive flares are fully consistent with the model. Recently, we have demonstrated that most observed coronal-hole jets in fan/spine topologies also are induced by breakout reconnection at the null point above a filament channel (with or without a filament).For these two types of eruptions occurring in similar topologies, the key question is, why do some events generate jets while others form CMEs? We focused on the initiation of eruptions in large bright points/small active regions that were located in coronal holes and clearly exhibited null-point (fan/spine) topologies: such configurations are referred to as pseudostreamers. We analyzed and compared SDO/AIA,SOHO/LASCO, and RHESSI observations of three events. Our analysis of the events revealed two new observable signatures of breakout reconnection prior to the explosive jet/CME outflows and flare onset: coronal dimming and the opening-up of field lines above the breakout current sheet. Most key properties were similar among the selected erupting structures, thereby eliminating region size, photospheric field strength, magnetic configuration, and pre-eruptive evolution as discriminating factors between jets and CMEs. We consider the factors that contribute to the different types of dynamic behavior, and conclude that the determining factor is the ratio of the magnetic free energy associated with the filament channel compared to the energy associated with the overlying flux inside and outside the pseudostreamer dome. Our results imply that a pseudostreamer eruption must be fast in order to produce an interplanetary CME.

Sun↗

NASA NTRS1993

Models of the Cartwheel ring galaxy: Spokes and starbursts

Recent observations of this famous ring galaxy, including optical and near-infrared CCD surface photometry, and VLA radio continuum and 21 cm line mapping (Higdon 1992b, in prep.), have inspired a renewed modeling effort. Toomre's (1978, in The Large-scale Structure of the Universe, eds. Longair and Einasto) series of restricted three-body simulations demonstrated how the multiple rings could be produced in a nearly head-on galaxy collision. New models with a halo-dominated potential based on the 21 cm rotation curve are able to reproduce such details as the spacing between rings, ring widths, offset of the nucleus, and several kinematical features, thus providing strong support for the collisional theory. The new observations have shown there are little or no old stars in Cartwheel; it may consist almost entirely of gas and stars produced as a result of compression in the ring wave. To model this process Smooth Particle Hydrodynamics (SPH) simulations of the Cartwheel disk have been performed. Fixed gravitational potentials were used to represent the Cartwheel and a roughly 30 percent mass collision partner. The interaction dynamics was treated as in the usual restricted three-body approximation, and the effects of local self-gravity between disk particles were calculated. We are particularly interested in testing the theory that enhanced star formation in waves is the result of gravitational instability in the compressed region (see e.g. Kennicutt 1989, ApJ 344, 685). The gas surface density in a number of simulations was initialized to a value slightly below the threshold for local gravitational instability throughout most of the disk. The first ring wave produces relatively modest compressions (a factor of order a few), triggering instability in a narrow range of wavelengths. Self-gravity in the disk is calculated over a comparable range of scales. Simulations were run with isothermal, adiabatic, and adiabatic with radiative cooling characterized by a relatively short timescale. The isothermal approximation is good except in the vicinity of the strong second (inner) ring, and several snapshots from one case are shown in the figure below. Flocculent spiral segments are present before the collision, and these are compressed into dense knots in the ring wave. These knots are likely to be sites of vigorous star formation. In the strong rarefaction behind the outer ring most of the knots are radially stretched and sheared, giving rise to spoke-like features. A few dense knots are evidently very tightly bound, because they retain their coherence and are stretched relatively little through the rarefaction. This is in accord with evidence for continuing star formation in some spokes (Marcum, Appleton and Higdon 1992). The number and spacing of spokes is a direct function of the scale of the gravitational instability in the disk. Thus, the gravitational instability theory, together with the hypothesis that massive stars are only formed in dense knots of gas, can account for most of the distinct morphology of the Cartwheel.

Struck-Marcell, Curtis↗

NASA NTRS1990

The temperature scale of hot DA white dwarfs - Temperatures from far-ultraviolet continuum fluxes

Observed far-UV to visible flux ratios of all DA white dwarfs hotter than about 25,000 K are compared with ratios predicted by models in order to derive effective temperatures for the stars. All suitable IUE spectra of seven hot DA white dwarfs for which accurate temperatures based on hydrogen line profiles were available are used to derive a flux correction to the 1980 IUE absolute calibration. Application of the corrections to the IUE spectra results in derived temperature which are consistent with the line profile temperatures.

Finley, David S.↗

NASA NTRS1983

The Spectral Properties of Gamma-ray Bursts: a Review of Recent Development

Developments in the spectroscopy of gamma ray bursts (GRB) are reviewed. The general question of the validity of the spectral results, particularly with regard to features in the spectrum, is discussed. Confirmations of these spectral features are summarized. Results from the KONUS experiments on Venera 13 and 14 are reviewed. The status of models of the continuum spectrum is summarized. A number of different radiation mechanisms appear capable of fitting the data. These include thermal bremsstrahlung, thermal synchrotron and inverse Compton. Rapid variability of the spectra shape on time scales 76] 0.25 sec. was reported. The characteristic energy of the spectrum was observed to vary over nearly an order of magnitude during individual events. A strong correlation between spectral hardness and luminosity was found. Low-energy (50 keV) absorption features and high-energy (400 keV) emission features continue to appear in GRB spectra. Understanding the origin of these lines in the context of the existing continuum models remains a difficult problem.

Teegarden, B. J.↗

NASA NTRSDate not supplied

Updates on the Predictive Materials Modeling Software Tools

Updates on NASA‘s efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. The PMM effort is part of the Entry Systems Modeling (ESM) project under NASA’s Game Changing Development (GCD) program. To reduce the need for extensive testing and accelerate the design cycle process, ESM is developing simulation and modeling tools that enable the characterization of the properties of thermal protection materials and their response to extremely hot plasma. The Porous Microstructure Analysis (PuMA) software has been developed to compute effective material properties and perform material response simulations on digitized microstructures of porous media. PuMA is able to import three-dimensional digital images obtained from X-ray microtomography or to generate artificial microstructures that mimic real materials. PuMA also provides a module for interactive 3D visualizations. Version 3, which was recently released as open-source, includes modules to compute simple morphological properties such as porosity, volume fractions, pore diameter, and specific surface area. Additional capabilities include the determination of effective thermal and electrical conductivity (both radiative and solid conduction - including the ability to simulate local anisotropy for the latter); effective diffusivity and tortuosity from the continuum to the rarefied regime; techniques to determine the local material orientation, as well as mechanical properties (elasticity coefficients), and permeability. Computed properties are then used to inform a macro-scale material response model, such as those implemented in the Porous material Analysis Toolbox based on OpenFOAM (PATO) software developed within ESM. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale, assuming Local Thermal Equilibrium. Recent efforts include the development of a mechanical erosion model as well as a unified model allowing an intrinsic coupling between fluid and material. Comparison to flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI] and Mars 2020 MEDLI2) is critical in order to validate these computational tools. Examples of ablative material response using the code will be presented, including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrated the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains through the use of massively parallel computations.

material modeling↗

NASA NTRS1997

Computational Modeling as a Design Tool in Microelectronics Manufacturing

Plans to introduce pilot lines or fabs for 300 mm processing are in progress. The IC technology is simultaneously moving towards 0.25/0.18 micron. The convergence of these two trends places unprecedented stringent demands on processes and equipments. More than ever, computational modeling is called upon to play a complementary role in equipment and process design. The pace in hardware/process development needs a matching pace in software development: an aggressive move towards developing "virtual reactors" is desirable and essential to reduce design cycle and costs. This goal has three elements: reactor scale model, feature level model, and database of physical/chemical properties. With these elements coupled, the complete model should function as a design aid in a CAD environment. This talk would aim at the description of various elements. At the reactor level, continuum, DSMC(or particle) and hybrid models will be discussed and compared using examples of plasma and thermal process simulations. In microtopography evolution, approaches such as level set methods compete with conventional geometric models. Regardless of the approach, the reliance on empricism is to be eliminated through coupling to reactor model and computational surface science. This coupling poses challenging issues of orders of magnitude variation in length and time scales. Finally, database development has fallen behind; current situation is rapidly aggravated by the ever newer chemistries emerging to meet process metrics. The virtual reactor would be a useless concept without an accompanying reliable database that consists of: thermal reaction pathways and rate constants, electron-molecule cross sections, thermochemical properties, transport properties, and finally, surface data on the interaction of radicals, atoms and ions with various surfaces. Large scale computational chemistry efforts are critical as experiments alone cannot meet database needs due to the difficulties associated with such controlled experiments and costs.

Meyyappan, Meyya↗