SEARCH · Search NASA

Results for “Dynamic Line Rating”

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 595 records · Page 33

All sources NASA NTRS DOE OSTI

NASA NTRS1981

Kaman 40 kW wind turbine generator - control system dynamics

The generator design incorporates an induction generator for application where a utility line is present and a synchronous generator for standalone applications. A combination of feed forward and feedback control is used to achieve synchronous speed prior to connecting the generator to the load, and to control the power level once the generator is connected. The dynamics of the drive train affect several aspects of the system operation. These were analyzed to arrive at the required shaft stiffness. The rotor parameters that affect the stability of the feedback control loop vary considerably over the wind speed range encountered. Therefore, the controller gain was made a function of wind speed in order to maintain consistent operation over the whole wind speed range. The velocity requirement for the pitch control mechanism is related to the nature of the wind gusts to be encountered, the dynamics of the system, and the acceptable power fluctuations and generator dropout rate. A model was developed that allows the probable dropout rate to be determined from a statistical model of wind gusts and the various system parameters, including the acceptable power fluctuation.

Perley, R.↗

DOE OSTI2023

Inversion-based correction of Double-Torsion (DT) subcritical crack growth tests for crack profile geometry

Because of its simplicity and the ability to produce a stable, slow-propagating crack, the Double-Torsion (DT) method has been used widely for investigating the critical and subcritical propagation of a slow-propagating tensile (mode-I) crack. However, to determine the complex relationship between the crack velocity $\mathcal{v_c}$ vs. the strain energy release rate $\mathscr{G}$ (or the stress intensity factor K) from laboratory measurements, several corrections must be made to account for the impact of sample and crack geometry. Particularly, DT test typically produces a crack with a curved edge profile instead of a straight line, causing the local $\mathcal{v_c}$ and $\mathscr{G}$ vary along the crack front. The experimentally measured $\mathcal{v_c}$ and $\mathscr{G}$ data merely reflect collective, averaged behavior of the crack. This makes inversion for the intrinsic, “true” crack growth kinetics necessary, based upon the knowledge of the crack geometry. Simple and effective correction methods have been proposed and validated for the slow, chemical-reaction-controlled part (Region I) of the $\mathcal{v_c}$-$\mathscr{G}$ curve. However, reliable methods for the highly nonlinear, transport-dominated part (Region II) and its sudden transition to the dynamic propagation part (Region III) are still lacking. Here we propose a method for determining the intrinsic $\mathcal{v_c}$-$\mathscr{G}$ relationship across all three Regions based upon DT test data, using a simple model function and its numerical inversion. The performance of this approach is examined and demonstrated using both synthetic and laboratory data for subcritical crack growth in soda lime glass.

42 ENGINEERING↗

DOE OSTI2023

FLASH: FPGA-Accelerated Smart Switches with GCN Case Study

Some communication switches, e.g., the Mellanox SHArP and those in the IBM BlueGene clusters, are augmented to process packets at the application level with fixed-function collectives. This approach, however, lacks flexibility, which limits their applicability in diverse and dynamic workloads. Recently, a new type of programmable packet processor, which uses high-level languages, e.g., P4, has emerged as possible candidates. P4-based switches, however, fall short in certain applications, including machine learning, where capabilities not currently supported by P4 are needed. These include more complex calculation, such as sparse computation and fused multiply-accumulate, data-intensive floating point operations, data reuse, and significant memory. The problem addressed here is that such a switch augmentation needs to support: a large amount of state, significant flexible compute capability, and ease of programming, all while maintaining full functionality, including ensuring high throughput, and demonstrating utility. In this work, we propose a programmable look-aside-type accelerator that can be embedded into, or attached to, existing communication switch pipelines and that is capable of processing packets at line-rate. The proposed in-switch accelerator is based on mixing an ISA (subset of RISC-V instructions) with dataflow graphs (found in CGRAs). To augment performance, vector instructions are also supported. To facilitate usability, we have developed a complete toolchain to compile user-provided C/C++ codes to appropriate back-end instructions for configuring the accelerator. While this approach is flexible enough to support various workloads, in this paper, we consider Graph Convolutional Networks (GCNs) as a case study. Experimental results show that this approach considerably improves the performance of distributed GCN applications.

Haghi, Pouya↗

DOE OSTI2022

Time dependent supervisory control update with FARM using rolling window

This report describes improvements to the Feasible Actuator Range Modifier (FARM) component of the RAVEN-based HYBRID framework for analysis of Integrated Energy Systems (IES). FARM supports the HERON plug-in that solves the power dispatch problem. The solution to the dispatch problem involves economically optimal dispatches that satisfy limits on production variables and their rates of variation (explicit constraints) as well as process variables tied to the service life of equipment (implicit constraints). FARM serves to validate or confirm that a HERON solution for explicit constraints also satisfies the implicit constraints. FARM-alpha was released by Argonne National Laboratory in January 2021 followed by FARM-Beta in January 2022 with the latter providing increased flexibility for the user. In this report, FARM-Gamma, the latest version of the code, is described. The major improvement is the implementation of a system identification algorithm based on the Dynamic Mode Decomposition with Control (DMDc) coupled with a “Rolling Window” scheme that allows obtaining linear time-varying state-space models. This feature equips FARM with the most accurate approximation of system dynamics, and it relieves the user from the burden of performing an exhaustive off-line characterization of the dynamics. FARM-Gamma capabilities are assessed by solving the power dispatch problem for a representative IES unit. The simulation times corresponding to the different releases are estimated and compared. These values capture the increasing computational burden of the successively higher fidelity state-space models adopted by FARM-Alpha, FARM-Beta and FARM-Gamma. The code implementation provides significant flexibility, i.e., the user can always select the most suitable version of FARM according to the problem to be solved and the available computational resources. It is anticipated that FARM will play a role in addressing several future IES applications. We outline how it can support the coordinated management and safe operation of a nuclear plant coupled to industrial processes to produce hydrogen and synfuels.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Lyα Observations of Comet C/2013 A1 (Siding Spring) Using MAVEN IUVS Echelle

The close approach of comet C/2013 A1 (Siding Spring) to Mars in 2014 October provided a unique opportunity to observe a dynamically new Oort cloud comet with potential for interaction with a planet’s atmosphere. The water-originating hydrogen coma of the comet extended to over 20 million km from the nucleus. Determining the properties of this coma contributes to characterizing the comet’s water content and production rate. The present study analyzes a unique data set of high spectral resolution UV observations of comet C/2013 A1 Siding Spring measured by the Mars Atmosphere and Volatile Evolution spacecraft. The Siding Spring observations capture Lyα emissions from the Martian corona, the interplanetary medium, as well as the cometary H and D reservoirs. The isolated cometary spectra are analyzed to reveal a velocity distribution of H atoms that are consistent with model estimates of H{sub 2}O photodissociated H emissions and of OH photodissociated H emissions, Doppler shifted from the main comet H emission line center by 18 km s{sup −1} and 8 km s{sup −1}, respectively. The variations in comet H brightness with distance from the nucleus are used to constrain cometary water production to a rate of 0.5 × 10{sup 28} molecules s{sup −1} at a time when Siding Spring was at 1.5 au, pre-perihelion.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

Measurement of the Neutron Elastic Electromagnetic Form Factor Ratio at Large Momentum Transfer

Exploring nucleon structure is vital both for understanding its origin and existence as well as for the advancement of the sciences. It helps us answer key questions such as how quark and gluon dynamics create 99% of the nucleon mass. Electron- nucleon scattering has been widely used for precision studies of the nucleon and nuclear structure since the Nobel Prize winning investigations by Robert Hofstadter and collaborators in the 1950s. These studies provide information about the spatial charge and current densities of the nucleon in terms of the electromagnetic form factors. The form factors are functions of four momentum transfer squared (Q2). Extending the electromagnetic form factor measurements to higher Q2 plays a critical role in furthering the understanding of nucleon structure. This motivated the Super BigBite Spectrometer (SBS) program at Jefferson Lab. The open nature of the spectrometers and the direct line of sight from the target to the tracking detector locations in experimental setups such as SBS creates high levels of background at the detectors. This necessitates the use of tracking detectors with high rate capability and good position resolution. Gas Electron Multiplier (GEM) detectors are an excellent choice for tracking detectors in such experiments. Understanding the performance of the GEM detectors is important not just for SBS experiments but also for future high-luminosity experiments. This thesis reports the exploratory results from the measurement of the neutron elastic electromagnetic form factor ratio (Gn E/Gn M) at high momentum transfer. A longitudinally polarized electron beam was scattered off a polarized 3He target, used as an effective polarized neutron target. In this experiment, the polarized 3He target achieved a world record polarization weighted luminosity at a beam current of 45 µA. Double spin asymmetry of the scattered neutron events is used to extract the neutron form factor ratio. Measurements were taken at Q2 = 3.0, 6.8, 9.8 (GeV/c)2. The lowest Q2 measurement is in good agreement with the existing world data, and the higher-Q2 measurements extend the Q2 reach well beyond the existing world data and are expected to remain unmatched for a long time.

Gamage, Vimukthi Haththotuwa [Univ. of Virginia, C↗

NASA NTRS1984

Gas dynamics in the impulsive phase of solar flares. I Thick-target heating by nonthermal electrons

A numerical investigation is carried out of the gas dynamical response of the solar atmosphere to a flare energy input in the form of precipitating nonthermal electrons. Rather than discussing the origin of these electrons, the spectral and temporal characteristics of the injected flux are inferred through a thick-target model of hard X-ray bremsstrahlung production. It is assumed that the electrons spiral about preexisting magnetic field lines, making it possible for a one-dimensional spatial treatment to be performed. It is also assumed that all electron energy losses are due to Coulomb collisions with ambient particles; that is, return-current ohmic effects and collective plasma processes are neglected. The results are contrasted with earlier work on conductive heating of the flare atmosphere. A local temperature peak is seen at a height of approximately 1500 km above the photosphere. This derives from a spatial maximum in the energy deposition rate from an electron beam. It is noted that such a feature is not present in conductively heated models. The associated localized region of high pressure drives material both upward and downward.

Nagai, F.↗

NASA NTRS2013

Sample Analysis at Mars Instrument Simulator

The Sample Analysis at Mars Instrument Simulator (SAMSIM) is a numerical model dedicated to plan and validate operations of the Sample Analysis at Mars (SAM) instrument on the surface of Mars. The SAM instrument suite, currently operating on the Mars Science Laboratory (MSL), is an analytical laboratory designed to investigate the chemical and isotopic composition of the atmosphere and volatiles extracted from solid samples. SAMSIM was developed using Matlab and Simulink libraries of MathWorks Inc. to provide MSL mission planners with accurate predictions of the instrument electrical, thermal, mechanical, and fluid responses to scripted commands. This tool is a first example of a multi-purpose, full-scale numerical modeling of a flight instrument with the purpose of supplementing or even eliminating entirely the need for a hardware engineer model during instrument development and operation. SAMSIM simulates the complex interactions that occur between the instrument Command and Data Handling unit (C&DH) and all subsystems during the execution of experiment sequences. A typical SAM experiment takes many hours to complete and involves hundreds of components. During the simulation, the electrical, mechanical, thermal, and gas dynamics states of each hardware component are accurately modeled and propagated within the simulation environment at faster than real time. This allows the simulation, in just a few minutes, of experiment sequences that takes many hours to execute on the real instrument. The SAMSIM model is divided into five distinct but interacting modules: software, mechanical, thermal, gas flow, and electrical modules. The software module simulates the instrument C&DH by executing a customized version of the instrument flight software in a Matlab environment. The inputs and outputs to this synthetic C&DH are mapped to virtual sensors and command lines that mimic in their structure and connectivity the layout of the instrument harnesses. This module executes, and thus validates, complex command scripts prior to their up-linking to the SAM instrument. As an output, this module generates synthetic data and message logs at a rate that is similar to the actual instrument.

Benna, Mehdi↗

DOE OSTI2014

Heating and dynamics of two flare loop systems observed by AIA and EIS

We investigate heating and evolution of flare loops in a C4.7 two-ribbon flare on 2011 February 13. From Solar Dynamics Observatory/Atmospheric Imaging Assembly (AIA) imaging observations, we can identify two sets of loops. Hinode/EUV Imaging Spectrometer (EIS) spectroscopic observations reveal blueshifts at the feet of both sets of loops. The evolution and dynamics of the two sets are quite different. The first set of loops exhibits blueshifts for about 25 minutes followed by redshifts, while the second set shows stronger blueshifts, which are maintained for about one hour. The UV 1600 observation by AIA also shows that the feet of the second set of loops brighten twice. These suggest that continuous heating may be present in the second set of loops. We use spatially resolved UV light curves to infer heating rates in the few tens of individual loops comprising the two loop systems. With these heating rates, we then compute plasma evolution in these loops with the 'enthalpy-based thermal evolution of loops' model. The results show that, for the first set of loops, the synthetic EUV light curves from the model compare favorably with the observed light curves in six AIA channels and eight EIS spectral lines, and the computed mean enthalpy flow velocities also agree with the Doppler shift measurements by EIS. For the second set of loops modeled with twice-heating, there are some discrepancies between modeled and observed EUV light curves in low-temperature bands, and the model does not fully produce the prolonged blueshift signatures as observed. We discuss possible causes for the discrepancies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Kepler-413B: A slightly misaligned, Neptune-size transiting circumbinary planet

We report the discovery of a transiting, R{sub p} = 4.347 ± 0.099R {sub ⊕}, circumbinary planet (CBP) orbiting the Kepler K+M eclipsing binary (EB) system KIC 12351927 (Kepler-413) every ∼66 days on an eccentric orbit with a{sub p} = 0.355 ± 0.002 AU, e{sub p} = 0.118 ± 0.002. The two stars, with M{sub A} = 0.820 ± 0.015 M {sub ☉}, R{sub A} = 0.776 ± 0.009 R {sub ☉} and M{sub B} = 0.542 ± 0.008 M {sub ☉}, R{sub B} = 0.484 ± 0.024 R {sub ☉}, respectively, revolve around each other every 10.11615 ± 0.00001 days on a nearly circular (e {sub EB} = 0.037 ± 0.002) orbit. The orbital plane of the EB is slightly inclined to the line of sight (i {sub EB} = 87.°33 ± 0.°06), while that of the planet is inclined by ∼2.°5 to the binary plane at the reference epoch. Orbital precession with a period of ∼11 yr causes the inclination of the latter to the sky plane to continuously change. As a result, the planet often fails to transit the primary star at inferior conjunction, causing stretches of hundreds of days with no transits (corresponding to multiple planetary orbital periods). We predict that the next transit will not occur until 2020. The orbital configuration of the system places the planet slightly closer to its host stars than the inner edge of the extended habitable zone. Additionally, the orbital configuration of the system is such that the CBP may experience Cassini State dynamics under the influence of the EB, in which the planet's obliquity precesses with a rate comparable to its orbital precession. Depending on the angular precession frequency of the CBP, it could potentially undergo obliquity fluctuations of dozens of degrees (and complex seasonal cycles) on precession timescales.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

MOSFIRE ABSORPTION LINE SPECTROSCOPY OF z > 2 QUIESCENT GALAXIES: PROBING A PERIOD OF RAPID SIZE GROWTH

Using the MOSFIRE near-infrared multi-slit spectrograph on the Keck 1 Telescope, we have secured high signal-to-noise ratio absorption line spectra for six massive galaxies with redshift 2 < z < 2.5. Five of these galaxies lie on the red sequence and show signatures of passive stellar populations in their rest-frame optical spectra. By fitting broadened spectral templates we have determined stellar velocity dispersions and, with broad-band Hubble Space Telescope and Spitzer photometry and imaging, stellar masses and effective radii. Using this enlarged sample of galaxies, we confirm earlier suggestions that quiescent galaxies at z > 2 have small sizes and large velocity dispersions compared to local galaxies of similar stellar mass. The dynamical masses are in very good agreement with stellar masses (log M {sub *}/M {sub dyn} = –0.02 ± 0.03), although the average stellar-to-dynamical mass ratio is larger than that found at lower redshift (–0.23 ± 0.05). By assuming evolution at fixed velocity dispersion, not only do we confirm a surprisingly rapid rate of size growth but we also consider the necessary evolutionary track on the mass-size plane and find a slope α = dlog R{sub e} /dlog M {sub *} ≳ 2 inconsistent with most numerical simulations of minor mergers. Both results suggest an additional mechanism may be required to explain the size growth of early galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2011

Proposed Framework for Determining Added Mass of Orion Drogue Parachutes

The Crew Exploration Vehicle (CEV) Parachute Assembly System (CPAS) project is executing a program to qualify a parachute system for a next generation human spacecraft. Part of the qualification process involves predicting parachute riser tension during system descent with flight simulations. Human rating the CPAS hardware requires a high degree of confidence in the simulation models used to predict parachute loads. However, uncertainty exists in the heritage added mass models used for loads predictions due to a lack of supporting documentation and data. Even though CPAS anchors flight simulation loads predictions to flight tests, extrapolation of these models outside the test regime carries the risk of producing non-bounding loads. A set of equations based on empirically derived functions of skirt radius is recommended as the simplest and most viable method to test and derive an enhanced added mass model for an inflating parachute. This will increase confidence in the capability to predict parachute loads. The selected equations are based on those published in A Simplified Dynamic Model of Parachute Inflation by Dean Wolf. An Ames 80x120 wind tunnel test campaign is recommended to acquire the reefing line tension and canopy photogrammetric data needed to quantify the terms in the Wolf equations and reduce uncertainties in parachute loads predictions. Once the campaign is completed, the Wolf equations can be used to predict loads in a typical CPAS Drogue Flight test. Comprehensive descriptions of added mass test techniques from the Apollo Era to the current CPAS project are included for reference.

Fraire, Usbaldo, Jr.↗

DOE OSTI2015

A STRONGLY COUPLED REACTOR CORE ISOLATION COOLING SYSTEM MODEL FOR EXTENDED STATION BLACK-OUT ANALYSES

The reactor core isolation cooling (RCIC) system in a boiling water reactor (BWR) provides makeup cooling water to the reactor pressure vessel (RPV) when the main steam lines are isolated and the normal supply of water to the reactor vessel is lost. The RCIC system operates independently of AC power, service air, or external cooling water systems. The only required external energy source is from the battery to maintain the logic circuits to control the opening and/or closure of valves in the RCIC systems in order to control the RPV water level by shutting down the RCIC pump to avoid overfilling the RPV and flooding the steam line to the RCIC turbine. It is generally considered in almost all the existing station black-out accidents (SBO) analyses that loss of the DC power would result in overfilling the steam line and allowing liquid water to flow into the RCIC turbine, where it is assumed that the turbine would then be disabled. This behavior, however, was not observed in the Fukushima Daiichi accidents, where the Unit 2 RCIC functioned without DC power for nearly three days. Therefore, more detailed mechanistic models for RCIC system components are needed to understand the extended SBO for BWRs. As part of the effort to develop the next generation reactor system safety analysis code RELAP-7, we have developed a strongly coupled RCIC system model, which consists of a turbine model, a pump model, a check valve model, a wet well model, and their coupling models. Unlike the traditional SBO simulations where mass flow rates are typically given in the input file through time dependent functions, the real mass flow rates through the turbine and the pump loops in our model are dynamically calculated according to conservation laws and turbine/pump operation curves. A simplified SBO demonstration RELAP-7 model with this RCIC model has been successfully developed. The demonstration model includes the major components for the primary system of a BWR, as well as the safety system components such as the safety relief valve (SRV), the RCIC system, the wet well, and the dry well. The results show reasonable system behaviors while exhibiting rich dynamics such as variable flow rates through RCIC turbine and pump during the SBO transient. The model has the potential to resolve the Fukushima RCIC mystery after adding the off-design two-phase turbine operation model and other additional improvements.

RCIC↗

DOE OSTI2014

MULTI-WAVELENGTH OBSERVATIONS OF COMET C/2011 L4 (PAN-STARRS)

The dynamically new comet C/2011 L4 (Pan-STARRS) is one of the brightest comets observed since the great comet C/1995 O1 (Hale-Bopp). Here, we present our multi-wavelength observations of C/2011 L4 during its in-bound passage to the inner solar system. A strong absorption band of water ice at 2.0 μm was detected in the near-infrared spectra, obtained with the 8 m Gemini-North and 3 m Infrared Telescope Facility Telescopes. The companion 1.5 μm band of water ice, however, was not observed. Spectral modeling shows that the absence of the 1.5 μm feature can be explained by the presence of sub-micron-sized fine ice grains. No gas lines (i.e., CN, HCN, or CO) were observed pre-perihelion in either the optical or the submillimeter. We derived 3σ upper limits for the CN and CO production rates. The comet exhibited a very strong continuum in the optical and its slope seemed to become redder as the comet approached the Sun. Our observations suggest that C/2011 L4 is an unusually dust-rich comet with a dust-to-gas mass ratio >4.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1997

Predictive Modeling in Plasma Reactor and Process Design

Research continues toward the improvement and increased understanding of high-density plasma tools. Such reactor systems are lauded for their independent control of ion flux and energy enabling high etch rates with low ion damage and for their improved ion velocity anisotropy resulting from thin collisionless sheaths and low neutral pressures. Still, with the transition to 300 mm processing, achieving etch uniformity and high etch rates concurrently may be a formidable task for such large diameter wafers for which computational modeling can play an important role in successful reactor and process design. The inductively coupled plasma (ICP) reactor is the focus of the present investigation. The present work attempts to understand the fundamental physical phenomena of such systems through computational modeling. Simulations will be presented using both computational fluid dynamics (CFD) techniques and the direct simulation Monte Carlo (DSMC) method for argon and chlorine discharges. ICP reactors generally operate at pressures on the order of 1 to 10 mTorr. At such low pressures, rarefaction can be significant to the degree that the constitutive relations used in typical CFD techniques become invalid and a particle simulation must be employed. This work will assess the extent to which CFD can be applied and evaluate the degree to which accuracy is lost in prediction of the phenomenon of interest; i.e., etch rate. If the CFD approach is found reasonably accurate and bench-marked with DSMC and experimental results, it has the potential to serve as a design tool due to the rapid time relative to DSMC. The continuum CFD simulation solves the governing equations for plasma flow using a finite difference technique with an implicit Gauss-Seidel Line Relaxation method for time marching toward a converged solution. The equation set consists of mass conservation for each species, separate energy equations for the electrons and heavy species, and momentum equations for the gas. The sheath is modeled by imposing the Bohm velocity to the ions near the walls. The DSMC method simulates each constituent of the gas as a separate species which would be analogous in CFD to employing separate species mass, momentum, and energy equations. All particles including electrons are moved and allowed to collide with one another with the stipulation that the electrons remain tied to the ions consistent with the concept of ambipolar diffusion. The velocities of the electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. These benefits come at a price in terms of computational time and memory. The DSMC and CFD are made as consistent as possible by using similar chemistry and power deposition models. Although the comparison of CFD and DSMC is interesting, the main goal of this work is the increased understanding of high-density plasma flowfields that can then direct improvements in both techniques. This work is unique in the level of the physical models employed in both the DSMC and CFD for high-density plasma reactor applications. For example, the electrons are simulated in the present DSMC work which has not been done before for low temperature plasma processing problems. In the CFD approach, for the first time, the charged particle transport (discharge physics) has been self-consistently coupled to the gas flow and heat transfer.

Hash, D. B.↗

NASA NTRS2004

The Search for Young Planetary Systems And the Evolution of Young Stars

The Space Interferometer Mission (SIM) will provide a census of planetary systems by con- ducting a broad survey of 2,000 stars that will be sensitive to the presence of planets with masses as small as approx. 15 Earth masses (1 Uranus mass) and a deep survey of approx. 250 of the nearest, stars with a mass limit of approx.3 Earth masses. The broad survey will include stars spanning a wide range of ages, spectral types, metallicity, and other important parameters. Within this larger context, the Young Stars and Planets Key Project will study approx. 200 stars with ages from 1 Myr to 100 Myr to understand the formation and dynamical evolution of gas giant planets. The SIM Young Stars and Planets Project will investigate both the frequency of giant planet formation and the early dynamical history of planetary systems. We will gain insight into how common the basic architecture of our solar system is compared with recently discovered systems with close-in giant planets by examining 200 of the nearest (less than 150 pc) and youngest (1-100 Myr) solar-type stars for planets. The sensitivity of the survey for stars located 140 pc away is shown in the planet mass-separation plane. We expect to find anywhere from 10 (assuming that only the presently known fraction of stars. 5-7%, has planets) to 200 (all young stars have planets) planetary systems. W-e have set our sensitivity threshold to ensure the detection of Jupiter-mass planets in the critical orbital range of 1 to 5 AU. These observations, when combined with the results of planetary searches of mature stars, will allow us to test theories of planetary formation and early solar system evolution. By searching for planets around pre-main sequence stars carefully selected to span an age range from 1 to 100 Myr, we will learn a t what epoch and with what frequency giant planets are found at the water-ice snowline where they are expected to form. This will provide insight into the physical mechanisms by which planets form and migrate from their place of birth, and about their survival rate. With these data in hand, we will provide data, for the first time, on such important questions as: What processes affect the formation and dynamical evolution of planets? When and where do planets form? What is initial mass distribution of planetary systems around young stars? How might planets be destroyed? What is the origin of the eccentricity of planetary orbits? What is the origin of the apparent dearth of companion objects between planets and brown dwarfs seen in mature stars? The observational strategy is a compromise between the desire to extend the planetary mass function as low as possible and the essential need to build up sufficient statistics on planetary occurrence. About half of the sample will be used to address the "where" and "when" of planet formation. We will study classical T Tauri stars (cTTs) which have massive accretion disks and post- accretion, weak-lined T Tauri stars (wTTs). Preliminary estimates suggest the sample will consist of approx. 30% cTTs and approx. 70% wTTs, driven in part by the difficulty of making accurate astrometric measurements toward objects with strong variability or prominent disks.

Beichman, Charles A.↗

DOE OSTI2013

STAR FORMATION IN TWO LUMINOUS SPIRAL GALAXIES

We examined star formation in two very luminous (M{sub V} = –22 to –23) Sc-type spiral galaxies, NGC 801 and UGC 2885, using ultra-deep Hα images. We combine these Hα images with UBV and Two-Micron All-Sky Survey JHK images and H I maps to explore the star formation characteristics of disk galaxies at high luminosity. Hα traces star formation in these galaxies to 4-6 disk scale lengths, but the lack of detection of Hα further out is likely due to the loss of Lyman continuum photons. Considering gravitational instabilities alone, we find that the gas and stars in the outer regions are marginally stable in an average sense, but considering dissipative gas and radial and azimuthal forcing, the outer regions are marginally unstable to forming spiral arms. Star formation is taking place in spiral arms, which are regions of locally higher gas densities. Furthermore, we have traced smooth exponential stellar disks over four magnitudes in V-band surface brightness and 4-6 disk scale lengths, in spite of a highly variable gravitational instability parameter. Thus, gravitational instability thresholds do not seem relevant to the stellar disk. One possibility for creating an exponential disk is that the molecular cloud densities and star formation rates have exponential profiles and this fact forces the stellar disk to build up such a profile. Another possibility is that the stellar disk is continuously adjusted to an exponential shape regardless of the star formation profile, for example, through global dynamical processes that scatter stars. However, such scattering processes are only known to operate in spiral systems, in which case they cannot explain the same dilemma of smooth exponential disks observed in dwarf irregular galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2013

THERMAL AND DYNAMICAL PROPERTIES OF GAS ACCRETING ONTO A SUPERMASSIVE BLACK HOLE IN AN ACTIVE GALACTIC NUCLEUS

We study stability of gas accretion in active galactic nuclei (AGNs). Our grid-based simulations cover a radial range from 0.1 to 200 pc, which may enable linking the galactic/cosmological simulations with small-scale black hole (BH) accretion models within a few hundreds of Schwarzschild radii. Here, as in previous studies by our group, we include gas radiative cooling as well as heating by a sub-Eddington X-ray source near the central supermassive BH of 10{sup 8} M{sub Sun }. Our theoretical estimates and simulations show that for the X-ray luminosity, L{sub X} {approx} 0.008 L{sub Edd}, the gas is thermally and convectively unstable within the computational domain. In the simulations, we observe that very tiny fluctuations in an initially smooth, spherically symmetric, accretion flow, grow first linearly and then nonlinearly. Consequently, an initially one-phase flow relatively quickly transitions into a two-phase/cold-hot accretion flow. For L{sub X} = 0.015 L{sub Edd} or higher, the cold clouds continue to accrete but in some regions of the hot phase, the gas starts to move outward. For L{sub X} < 0.015 L{sub Edd}, the cold phase contribution to the total mass accretion rate only moderately dominates over the hot phase contribution. This result might have some consequences for cosmological simulations of the so-called AGN feedback problem. Our simulations confirm the previous results of Barai et al. who used smoothed particle hydrodynamic (SPH) simulations to tackle the same problem. Here, however, because we use a grid-based code to solve equations in one dimension and two dimensions, we are able to follow the gas dynamics at much higher spacial resolution and for longer time compared with the three-dimensional SPH simulations. One of the new features revealed by our simulations is that the cold condensations in the accretion flow initially form long filaments, but at the later times, those filaments may break into smaller clouds advected outward within the hot outflow. Therefore, these simulations may serve as an attractive model for the so-called narrow-line region in AGNs.

79 ASTRONOMY AND ASTROPHYSICS↗