Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Cluster computing”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Real-Time Kennedy Space Center and Cape Canaveral Air Force Station High-Resolution Model Implementation and Verification

Customer: NASA's Launch Services Program (LSP), Ground Systems Development and Operations (GSDO), and Space Launch System (SLS) programs. NASA's LSP, GSDO, SLS and other programs at Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS) use the daily and weekly weather forecasts issued by the 45th Weather Squadron (45 WS) as decision tools for their day-to-day and launch operations on the Eastern Range (ER). For example, to determine if they need to limit activities such as vehicle transport to the launch pad, protect people, structures or exposed launch vehicles given a threat of severe weather, or reschedule other critical operations. The 45 WS uses numerical weather prediction models as a guide for these weather forecasts, particularly the Air Force Weather Agency (AFWA) 1.67 kilometer Weather Research and Forecasting (WRF) model. Considering the 45 WS forecasters' and Launch Weather Officers' (LWO) extensive use of the AFWA model, the 45 WS proposed a task at the September 2013 Applied Meteorology Unit (AMU) Tasking Meeting requesting the AMU verify this model. Due to the lack of archived model data available from AFWA, verification is not yet possible. Instead, the AMU proposed to implement and verify the performance of an ER version of the AMU high-resolution WRF Environmental Modeling System (EMS) model (Watson 2013) in real-time. The tasking group agreed to this proposal; therefore the AMU implemented the WRF-EMS model on the second of two NASA AMU modeling clusters. The model was set up with a triple-nested grid configuration over KSC/CCAFS based on previous AMU work (Watson 2013). The outer domain (D01) has 12-kilometer grid spacing, the middle domain (D02) has 4-kilometer grid spacing, and the inner domain (D03) has 1.33-kilometer grid spacing. The model runs a 12-hour forecast every hour, D01 and D02 domain outputs are available once an hour and D03 is every 15 minutes during the forecast period. The AMU assessed the WRF-EMS 1.33-kilometer domain model performance for the 2014 warm season (May-September). Verification statistics were computed using the Model Evaluation Tools, which compared the model forecasts to observations. The mean error values were close to 0 and the root mean square error values were less than 1.8 for mean sea-level pressure (millibars), temperature (degrees Kelvin), dewpoint temperature (degrees Kelvin), and wind speed (per millisecond), all very small differences between the forecast and observations considering the normal magnitudes of the parameters. The precipitation forecast verification results showed consistent under-forecasting of the precipitation object size. This could be an artifact of calculating the statistics for each hour rather than for the entire 12-hour period. The AMU will continue to generate verification statistics for the 1.33-kilometer WRF-EMS domain as data become available in future cool and warm seasons. More data will produce more robust statistics and reveal a more accurate assessment of model performance. Once the formal task was complete, the AMU conducted additional work to better understand the wind direction results. The results were stratified diurnally and by wind speed to determine what effects the stratifications would have on the model wind direction verification statistics. The results are summarized in the addendum at the end of this report. In addition to verifying the model's performance, the AMU also made the output available in the Advanced Weather Interactive Processing System II (AWIPS II). This allows the 45 WS and AMU staff to customize the model output display on the AMU and Range Weather Operations AWIPS II client computers and conduct real-time subjective analyses. In the future, the AMU will implement an updated version of the WRF-EMS model that incorporates local data assimilation. This model will also run in real-time and be made available in AWIPS II.

Numerical Weather Prediction↗
NASA NTRSDate not supplied

Communicating Metrics of Land Surface Temperature Variability Using Multi-sensor Machine Learning

Land surface temperature (LST) is a key climate observable used to detect changes in the Earth’s surface energy budget that influence carbon and water cycles. Land surface temperature exhibits strong diurnal variability, which geostationary satellites can observe at scale thanks to their temporal resolution. Due to anthropogenic climate and land use changes, the surface energy balance has been considerably modified and may be described by changes in diurnal temperature range and extremes. Using high performance computing and datasets from the NASA Earth Exchange, we exploit co-located, co-temporal observations from low-earth orbit (LEO) and geostationary (GEO) sensors to develop a deep learning-based method for LEO-to-GEO algorithm emulation. Our model is trained to predict MODIS Terra LST from GOES-16 thermal bands and achieves validation error <2K. Application of the model to unseen times of day (observed by MODIS Aqua) and a new GEO sensor (Himawari-8) observing an unseen spatial domain, demonstrate the generalization of the deep learning model across space, time and spectra. Further, time series clustering approaches are examined with the objective of identifying key indicators of change in diurnal cycling and extremes on a continental scale. Communicating LST variability observed by geostationary satellites can have impacts in multiple disciplines, from understanding of snow, vegetation and soil dynamics, to recognizing trends in heat events relevant to human health.

Kate Duffy↗
NASA NTRS2006

Hera - The HEASARC's New Data Analysis Service

Hera is the new computer service provided by the HEASARC at the NASA Goddard Space Flight Center that enables qualified student and professional astronomical researchers to immediately begin analyzing scientific data from high-energy astrophysics missions. All the necessary resources needed to do the data analysis are freely provided by Hera, including: * the latest version of the hundreds of scientific analysis programs in the HEASARC's HEASOFT package, as well as most of the programs in the Chandra CIAO package and the XMM-Newton SAS package. * high speed access to the terabytes of data in the HEASARC's high energy astrophysics Browse data archive. * a cluster of fast Linw workstations to run the software * ample local disk space to temporarily store the data and results. Some of the many features and different modes of using Hera are illustrated in this poster presentation.

Pence, William↗
NASA NTRS1980

Documentation of program AFTBDY to generate coordinate system for 3D after body using body fitted curvilinear coordinates, part 1

The computer program AFTBDY generates a body fitted curvilinear coordinate system for a wedge curved after body. This wedge curved after body is being used in an experimental program. The coordinate system generated by AFTBDY is used to solve 3D compressible N.S. equations. The coordinate system in the physical plane is a cartesian x,y,z system, whereas, in the transformed plane a rectangular xi, eta, zeta system is used. The coordinate system generated is such that in the transformed plane coordinate spacing in the xi, eta, zeta direction is constant and equal to unity. The physical plane coordinate lines in the different regions are clustered heavily or sparsely depending on the regions where physical quantities to be solved for by the N.S. equations have high or low gradients. The coordinate distribution in the physical plane is such that x stays constant in eta and zeta direction, whereas, z stays constant in xi and eta direction. The desired distribution in x and z is input to the program. Consequently, only the y-coordinate is solved for by the program AFTBDY.

Kumar, D.↗
NASA NTRS2005

New Directions in Space Operations Services in Support of Interplanetary Exploration

To gain access to the necessary operational processes and data in support of NASA's Lunar/Mars Exploration Initiative, new services, adequate levels of computing cycles and access to myriad forms of data must be provided to onboard spacecraft and ground based personnel/systems (earth, lunar and Martian) to enable interplanetary exploration by humans. These systems, cycles and access to vast amounts of development, test and operational data will be required to provide a new level of services not currently available to existing spacecraft, on board crews and other operational personnel. Although current voice, video and data systems in support of current space based operations has been adequate, new highly reliable and autonomous processes and services will be necessary for future space exploration activities. These services will range from the more mundane voice in LEO to voice in interplanetary travel which because of the high latencies will require new voice processes and standards. New services, like component failure predictions based on data mining of significant quantities of data, located at disparate locations, will be required. 3D or holographic representation of onboard components, systems or family members will greatly improve maintenance, operations and service restoration not to mention crew morale. Current operational systems and standards, like the Internet Protocol, will not able to provide the level of service required end to end from an end point on the Martian surface like a scientific instrument to a researcher at a university. Ground operations whether earth, lunar or Martian and in flight operations to the moon and especially to Mars will require significant autonomy that will require access to highly reliable processing capabilities, data storage based on network storage technologies. Significant processing cycles will be needed onboard but could be borrowed from other locations either ground based or onboard other spacecraft. Reliability will be a key factor with onboard and distributed backup processing an absolutely necessary requirement. Current cluster processing/Grid technologies may provide the basis for providing these services. An overview of existing services, future services that will be required and the technologies and standards required to be developed will be presented. The purpose of this paper will be to initiate a technological roadmap, albeit at a high level, of current voice, video, data and network technologies and standards (which show promise for adaptation or evolution) to what technologies and standards need to be redefined, adjusted or areas where new ones require development. The roadmap should begin the differentiation between non manned and manned processes/services where applicable. The paper will be based in part on the activities of the CCSDS Monitor and Control working group which is beginning the process of standardization of the these processes. Another element of the paper will be based on an analysis of current technologies supporting space flight processes and services at JSC, MSFC, GSFC and to a lesser extent at KSC. Work being accomplished in areas such as Grid computing, data mining and network storage at ARC, IBM and the University of Alabama at Huntsville will be researched and analyzed.

Bradford, Robert N.↗
NASA NTRS1975

Wildlife management by habitat units: A preliminary plan of action

Procedures for yielding vegetation type maps were developed using LANDSAT data and a computer assisted classification analysis (LARSYS) to assist in managing populations of wildlife species by defined area units. Ground cover in Travis County, Texas was classified on two occasions using a modified version of the unsupervised approach to classification. The first classification produced a total of 17 classes. Examination revealed that further grouping was justified. A second analysis produced 10 classes which were displayed on printouts which were later color-coded. The final classification was 82 percent accurate. While the classification map appeared to satisfactorily depict the existing vegetation, two classes were determined to contain significant error. The major sources of error could have been eliminated by stratifying cluster sites more closely among previously mapped soil associations that are identified with particular plant associations and by precisely defining class nomenclature using established criteria early in the analysis.

Frentress, C. D.↗
NASA NTRS1983

Directivity of planar array feeds for satellite reflector applications

Array-fed reflector antennas are used extensively in today's contour and multiple beam satellite antennas. To determine the directivity of these antennas theoretically, the total radiated power of the array feed must be accurately computed. In this paper, a closed-form expression for the radiated power is obtained for arrays with the well-known (cos theta q-type element patterns. The formulation is general, and takes into account polarization, nonsymmetric E- and H-plane element patterns, nonuniform element spacings and arbitrary complex excitation coefficients. Selective numerical data are presented to demonstrate the usefulness of this result. Comparisons are made with the available results obtained using direct numerical integration techniques, and with other available data based on less general formulations. Excellent agreement is observed for all cases. In particular, data are presented on the directivity of seven-element cluster feeds used in multiple beam designs, and an array feed for producing a contour beam covering the United States Eastern Time Zone (ETZ).

Rahmat-Samii, Y.↗
NASA NTRS2002

Constitutive Relation in Transitional Granular Flows

To study the constitutive behavior of granular materials, the presence of gravity is detrimental. Although empirical relations have been obtained for engineering designs to control granular flows on Earth, it is not known how well these Earthbound relations can be used in another gravity field. Fundamental understanding must be derived to reliably design for granular flows in space exploration. There are two extremes of granular flows of which significant amount of knowledge is available. One deals with a dense and quasi-static situation where the deformation rate nearly vanishes. The other deals with dilute and rapidly fluctuating grain velocities where particle inertia dominates. This project, funded by the NASA Microgravity Fluid Physics Program, aims to study this transitional regime via physical experiments and computer simulations. A conceptual model has been established as described below. There are two natural time scales in a granular flow. One is the travel time between two consecutive collisions and the other is the duration of a collision contact. At a very low shear-rate, the shear-induced particle velocity is low. Hence the travel time between collisions is longer than the contact time between colliding particles. Binary collisions prevail. As the shear-rate increases, the traveling time between collisions reduces and the probability of multiple collisions goes up. These particle groups disperse shortly after and new groups form. When shear-rate is further increased, clusters grow in size due to an increasing chance for free particles to join before groups have the time to disperse. The maximum cluster size may depend on the global concentration and material properties. As the solid concentration approaches zero, the cluster size goes to one particle diameter. The maximum possible cluster size under any condition is the container size, provided that the shear flow is inside a container. The critical shear-rates that dictate the initiation of the multiple contacts, and the size and lifetime of the collision clusters, are functions of the concentration also. A 'regime' theory has been proposed by Babic et al. This theory suggested that both the solid concentration, C, and the non-dimensional shear-rate, B, are important in determining the regimes of the granular constitutive law. Additional information is included in the original extended abstract.

Shen, Hayley H.↗
NASA NTRS1982

Methods of editing cloud and atmospheric layer affected pixels from satellite data

Practical methods of computer screening cloud-contaminated pixels from data of various satellite systems are proposed. Examples are given of the location of clouds and representative landscape features in HCMM spectral space of reflectance (VIS) vs emission (IR). Methods of screening out cloud affected HCMM are discussed. The character of subvisible absorbing-emitting atmospheric layers (subvisible cirrus or SCi) in HCMM data is considered and radiosonde soundings are examined in relation to the presence of SCi. The statistical characteristics of multispectral meteorological satellite data in clear and SCi affected areas are discussed. Examples in TIROS-N and NOAA-7 data from several states and Mexico are presented. The VIS-IR cluster screening method for removing clouds is applied to a 262, 144 pixel HCMM scene from south Texas and northeast Mexico. The SCi that remain after cluster screening are sited out by applying a statistically determined IR limit.

Nixon, P. R.↗
NASA NTRS2003

Iron catalyst chemistry in modeling a high-pressure carbon monoxide nanotube reactor

The high-pressure carbon monoxide (HiPco) technique for producing single-wall carbon nanotubes (SWNTs) is analyzed with the use of a chemical reaction model coupled with flow properties calculated along streamlines, calculated by the FLUENT code for pure carbon monoxide. Cold iron pentacarbonyl, diluted in CO at about 30 atmospheres, is injected into a conical mixing zone, where hot CO is also introduced via three jets at 30 degrees with respect to the axis. Hot CO decomposes the Fe(CO)5 to release atomic Fe. Then iron nucleates and forms clusters that catalyze the formation of SWNTs by a disproportionation reaction (Boudouard) of CO on Fe-containing clusters. Alternative nucleation rates are estimated from the theory of hard sphere collision dynamics with an activation energy barrier. The rate coefficient for carbon nanotube growth is estimated from activation energies in the literature. The calculated growth was found be about an order of magnitude greater than measured, regardless of the nucleation rate. A study of cluster formation in an incubation zone prior to injection into the reactor shows that direct dimer formation from Fe atoms is not as important as formation via an exchange reaction of Fe with CO in FeCO.

Evaluation Studies↗
NASA NTRS1992

Research in computational fluid dynamics and analysis of algorithms

Recently, higher-order compact schemes have seen increasing use in the DNS (Direct Numerical Simulations) of the Navier-Stokes equations. Although they do not have the spatial resolution of spectral methods, they offer significant increases in accuracy over conventional second order methods. They can be used on any smooth grid, and do not have an overly restrictive CFL dependence as compared with the O(N(exp -2)) CFL dependence observed in Chebyshev spectral methods on finite domains. In addition, they are generally more robust and less costly than spectral methods. The issue of the relative cost of higher-order schemes (accuracy weighted against physical and numerical cost) is a far more complex issue, depending ultimately on what features of the solution are sought and how accurately they must be resolved. In any event, the further development of the underlying stability theory of these schemes is important. The approach of devising suitable boundary clusters and then testing them with various stability techniques (such as finding the norm) is entirely the wrong approach when dealing with high-order methods. Very seldom are high-order boundary closures stable, making them difficult to isolate. An alternative approach is to begin with a norm which satisfies all the stability criteria for the hyperbolic system, and look for the boundary closure forms which will match the norm exactly. This method was used recently by Strand to isolate stable boundary closure schemes for the explicit central fourth- and sixth-order schemes. The norm used was an energy norm mimicking the norm for the differential equations. Further research should be devoted to BC for high order schemes in order to make sure that the results obtained are reliable. The compact fourth order and sixth order finite difference scheme had been incorporated into a code to simulate flow past circular cylinders. This code will serve as a verification of the full spectral codes. A detailed stability analysis by Carpenter (from the fluid Mechanics Division) and Gottlieb gave analytic conditions for stability as well as asymptotic stability. This had been incorporated in the code in form of stable boundary conditions. Effects of the cylinder rotations had been studied. The results differ from the known theoretical results. We are in the middle of analyzing the results. A detailed analysis of the effects of the heating of the cylinder on the shedding frequency had been studied using the above schemes. It has been found that the shedding frequency decreases when the wire was heated. Experimental work is being carried out to affirm this result.

Gottlieb, David↗
NASA NTRS2012

NASA Center for Climate Simulation (NCCS) Presentation

The NASA Center for Climate Simulation (NCCS) offers integrated supercomputing, visualization, and data interaction technologies to enhance NASA's weather and climate prediction capabilities. It serves hundreds of users at NASA Goddard Space Flight Center, as well as other NASA centers, laboratories, and universities across the US. Over the past year, NCCS has continued expanding its data-centric computing environment to meet the increasingly data-intensive challenges of climate science. We doubled our Discover supercomputer's peak performance to more than 800 teraflops by adding 7,680 Intel Xeon Sandy Bridge processor-cores and most recently 240 Intel Xeon Phi Many Integrated Core (MIG) co-processors. A supercomputing-class analysis system named Dali gives users rapid access to their data on Discover and high-performance software including the Ultra-scale Visualization Climate Data Analysis Tools (UV-CDAT), with interfaces from user desktops and a 17- by 6-foot visualization wall. NCCS also is exploring highly efficient climate data services and management with a new MapReduce/Hadoop cluster while augmenting its data distribution to the science community. Using NCCS resources, NASA completed its modeling contributions to the Intergovernmental Panel on Climate Change (IPCG) Fifth Assessment Report this summer as part of the ongoing Coupled Modellntercomparison Project Phase 5 (CMIP5). Ensembles of simulations run on Discover reached back to the year 1000 to test model accuracy and projected climate change through the year 2300 based on four different scenarios of greenhouse gases, aerosols, and land use. The data resulting from several thousand IPCC/CMIP5 simulations, as well as a variety of other simulation, reanalysis, and observationdatasets, are available to scientists and decision makers through an enhanced NCCS Earth System Grid Federation Gateway. Worldwide downloads have totaled over 110 terabytes of data.

Webster, William P.↗
NASA NTRS2001

Coupled-Flow Simulation of HP-LP Turbines Has Resulted in Significant Fuel Savings

Our objective was to create a high-fidelity Navier-Stokes computer simulation of the flow through the turbines of a modern high-bypass-ratio turbofan engine. The simulation would have to capture the aerodynamic interactions between closely coupled high- and low-pressure turbines. A computer simulation of the flow in the GE90 turbofan engine's high-pressure (HP) and low-pressure (LP) turbines was created at GE Aircraft Engines under contract with the NASA Glenn Research Center. The three-dimensional steady-state computer simulation was performed using Glenn's average-passage approach named APNASA. The areas upstream and downstream of each blade row mutually interact with each other during engine operation. The embedded blade row operating conditions are modeled since the average passage equations in APNASA actively include the effects of the adjacent blade rows. The turbine airfoils, platforms, and casing are actively cooled by compressor bleed air. Hot gas leaks around the tips of rotors through labyrinth seals. The flow exiting the high work HP turbines is partially transonic and, therefore, has a strong shock system in the transition region. The simulation was done using 121 processors of a Silicon Graphics Origin 2000 (NAS 02K) cluster at the NASA Ames Research Center, with a parallel efficiency of 87 percent in 15 hr. The typical average-passage analysis mesh size per blade row was 280 by 45 by 55, or approx.700,000 grid points. The total number of blade rows was 18 for a combined HP and LP turbine system including the struts in the transition duct and exit guide vane, which contain 12.6 million grid points. Design cycle turnaround time requirements ran typically from 24 to 48 hr of wall clock time. The number of iterations for convergence was 10,000 at 8.03x10(exp -5) sec/iteration/grid point (NAS O2K). Parallel processing by up to 40 processors is required to meet the design cycle time constraints. This is the first-ever flow simulation of an HP and LP turbine. In addition, it includes the struts in the transition duct and exit guide vanes.

Veres, Joseph P.↗
NASA NTRS2011

Design of Refractory Metal Heat Pipe Life Test Environment Chamber, Cooling System, and Radio Frequency Heating System

A series of 16 Mo-44.5%Re alloy/sodium heat pipes will be experimentally tested to examine heat pipe aging. To support this evaluation, an environmental test chamber and a number of auxiliary subsystems are required. These subsystems include radio frequency (RF) power supplies/inductive coils, recirculation water coolant loops, and chamber gas conditioning. The heat pipes will be grouped, based on like power and gas mixture requirements, into three clusters of five units each, configured in a pentagonal arrangement. The highest powered heat pipe will be tested separately. Test chamber atmospheric purity is targeted at <0.3 ppb oxygen at an approximate operating pressure of 76 torr (.1.5 psia), maintained by active purification (oxygen level is comparable to a 10(exp -6) torr environment). Treated water will be used in two independent cooling circuits to remove .85 kW. One circuit will service the RF hardware while the other will maintain the heat pipe calorimetry. Initial procedures for the startup and operation of support systems have been identified. Each of these subsystems is outfitted with a variety of instrumentation, integrated with distributed real-time controllers and computers. A local area network provides communication between all devices. This data and control network continuously monitors the health of the test hardware, providing warning indicators followed by automatic shutdown should potentially damaging conditions develop. During hardware construction, a number of checkout tests.many making use of stainless steel prototype heat pipes that are already fabricated.will be required to verify operation.

Martin, J. J.↗
NASA NTRS1979

The evolution of rotating stars. III - Predicted surface rotation velocities for stars which conserve total angular momentum

Predicted surface rotation velocities for Population I stars at 10, 7, 5, 3, and 1.5 solar masses are presented. The surface velocities were computed for angular momentum with no radial redistribution, complete redistribution, and partial redistribution as predicted by consideration of circulation currents in rotating stars. Near the main sequence, rotational effects can reduce the moment of inertia of a star, so nonrotating models underestimate the expected velocities for evolving stars. On the red giant branch, angular momentum redistribution reduces the surface velocity by a factor of 2 or more, relative to the velocity expected for no radial redistribution. This removes the discrepancy between predicted and observed rotation rates for the K giants and makes it unlikely that these stars lose significant amounts of angular momentum by stellar winds. Calculations indicate that improved observations of the red giants in the Hyades cluster can be used to determine how angular momentum is redistributed by convection

Endal, A. S.↗
NASA NTRS1988

Synthetic Stromgren photometry for F dwarf stars

Recent synthetic spectrum and color calculations for cool dwarf star models are tested by comparison with observation. The accuracy of the computed dependence of the thermal colors B-V and b-y on effective temperature is examined, and H-beta indices are presented and compared with observed values. The accuracy of the predictions of the Stromgren uvby system metal-abundance indicator m1 and luminosity indicator c1 are tested. A new calibration of the c1, b-y diagram in terms of absolute magnitudes is given, making use of recent calculations of stellar isochrones. Observations of very metal-poor subdwarfs are used to study the accuracy of the isochrones. The c1, b-y diagram of the subdwarfs is compared with that of the turnoff-region stars in the very metal-poor globular cluster NGC 6397.

Bell, R. A.↗
NASA NTRS2000

Design and Evaluation of a New Boundary-Layer Rake for Flight Testing

A new boundary-layer rake has been designed and built for flight testing on the NASA Dryden Flight Research Center F-15B/Flight Test Fixture. A feature unique to this rake is its curved body, which allows pitot tubes to be more densely clustered in the near-wall region than conventional rakes allow. This curved rake design has a complex three-dimensional shape that requires innovative solid-modeling and machining techniques. Finite-element stress analysis of the new design shows high factors of safety. The rake has passed a ground test in which random vibration measuring 12 g rms was applied for 20 min in each of the three normal directions. Aerodynamic evaluation of the rake has been conducted in the NASA Glenn Research Center 8 x 6 Supersonic Wind Tunnel at Mach 0-2. The pitot pressures from the new rake agree with conventional rake data over the range of Mach numbers tested. The boundary-layer profiles computed from the rake data have been shown to have the standard logarithmic-law profile. Skin friction values computed from the rake data using the Clauser plot method agree with the Preston tube results and the van Driest II compressible skin friction correlation to approximately +/-5 percent.

Bui, Trong T.↗
NASA NTRS2002

On the Formation of Massive Stars

We calculate numerically the collapse of slowly rotating, nonmagnetic, massive molecular clumps of masses 30,60, and 120 Stellar Mass, which conceivably could lead to the formation of massive stars. Because radiative acceleration on dust grains plays a critical role in the clump's dynamical evolution, we have improved the module for continuum radiation transfer in an existing two-dimensional (axial symmetry assumed) radiation hydrodynamic code. In particular, rather than using "gray" dust opacities and "gray" radiation transfer, we calculate the dust's wavelength-dependent absorption and emission simultaneously with the radiation density at each wavelength and the equilibrium temperatures of three grain components: amorphous carbon particles. silicates, and " dirty ice " -coated silicates. Because our simulations cannot spatially resolve the innermost regions of the molecular clump, however, we cannot distinguish between the formation of a dense central cluster or a single massive object. Furthermore, we cannot exclude significant mass loss from the central object(s) that may interact with the inflow into the central grid cell. Thus, with our basic assumption that all material in the innermost grid cell accretes onto a single object. we are able to provide only an upper limit to the mass of stars that could possibly be formed. We introduce a semianalytical scheme for augmenting existing evolutionary tracks of pre-main-sequence protostars by including the effects of accretion. By considering an open outermost boundary, an arbitrary amount of material could, in principal, be accreted onto this central star. However, for the three cases considered (30, 60, and 120 Stellar Mass originally within the computation grid), radiation acceleration limited the final masses to 3 1.6, 33.6, and 42.9 Stellar Mass, respectively, for wavelength-dependent radiation transfer and to 19.1, 20.1, and 22.9 Stellar Mass. for the corresponding simulations with gray radiation transfer. Our calculations demonstrate that massive stars can in principle be formed via accretion through a disk. The accretion rate onto the central source increases rapidly after one initial free-fall time and decreases monotonically afterward. By enhancing the nonisotropic character of the radiation field, the accretion disk reduces the effects of radiative acceleration in the radial direction - a process we call the "flashlight effect." The flashlight effect is further amplified in our case by including the effects of frequency-dependent radiation transfer. We conclude with the warning that a careful treatment of radiation transfer is a mandatory requirement for realistic simulations of the formation of massive stars.

Yorke, Harold W.↗