SEARCH · Engineering Papers

Results for “Calculational Margin”

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 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Methanol and hydrogen sulfide in comet P/Halley

The Neutral Mass Spectrometer on the Giotto spacecraft measured the gas and ion composition in the coma of comet P/Halley. A detailed model of the ion chemistry inside the contact surface located at 4660 km is used to interpret the measured ion desnity profiles in the 32 to 35 amu/e mass range. The masses 33 and 35 amu/e are dominated by the protonated methanol and hydrogen sulfide ions CH3OH2(+) and H3S(+). Both profiles are essentially compatible with CH3OH and H2S originating from the nucleus only. The production rates relative to water are Y(CH3OH) = Q(CH3OH)/Q(H2O) = 1.7% and Y(H2S) = 0.41%. Our Y(CH3OH) agrees well with a determination from IR spectra obtained about 6 weeks after the Giotto encounter with P/Halley. In 7 other comets IR and microwave observations give Y(CH3OH) values between about 0.7 and 6%, indicating that the methanol abundance shows a strong variability from comet to comet. In three other comets Y(H2S) values between 0.2 and 0.5% have been reported. In addition to H2S(+), only ions containing minor isotopes of H, C, O and S contribute to mass 34 amu/e (e.g. (34)S(+), (13)CH3OH2(+), CH4DO(+)). These contributions can be calculated from the measured densities of the ions containing the major isotopes and the H2S(+) contribution from the measured H3S(+) density. From mass 34 amu/e we can also derive an upper limit of 1% for the abundance of deuterated methanol. This limit is at most marginally compatible with a direct interstellar origin of the CH3OH in P/Halley as the measured interstellar abundance of deuterated methanol is 1 to 6%.

Eberhardt, P.↗

DOE OSTI2023

A User’s Guide to the PLTEMP/ANL Code

PLTEMP/ANL V4.4 is a program that obtains a steady-state flow and temperature solution, in geometrical extent, for a nuclear reactor core, or a single fuel assembly, including the effect of manufacturing and modelling uncertainties by the means of hot channel factors. It is based on an evolutionary sequence of codes originally used for calculating plate temperatures, hence “PLTEMP”, developed at Argonne National Laboratory over the past 30 years. Fueled and non-fueled regions are modeled. Each fuel assembly consists of one or more plates or tubes separated by coolant channels. The fuel plates may have one to five layers of different materials, each with heat generation. The width of a fuel plate may be divided into multiple longitudinal stripes, each with its own axial power shape. Depending upon the heat transfer method selected, the temperature solution is effectively 2-dimensional or 3-dimensional. The geometry may be either slab or radial, corresponding to fuel assemblies made of a series of flat (or slightly curved) plates, or of nested tubes. A variety of thermal-hydraulic correlations are available to determine safety margins such as onset of nucleate boiling ratio (ONBR), departure from nucleate boiling ratio (DNBR), and onset of flow instability ratio (OFIR). Coolant properties for either light or heavy water are obtained from FORTRAN functions rather than from tables. The code is intended for thermal-hydraulic analysis of research reactor performance in the sub-cooled boiling regime. Both turbulent and laminar flow regimes can be modeled. Options to calculate both forced flow and natural circulation are available. A general search capability for select design and safety parameters is available to greatly reduce the reactor analyst’s time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1977

Analysis of internal flow of J85-13 multistage compressor

Interstage data recorded on a J85-13 engine were used to analyze the internal flow of the compressor. Measured pressures and temperatures were used as input to a streamline analysis program to calculate the velocity diagrams at the inlet and outlet of each blade row. From the velocity diagrams and blade geometry, selected blade-element performance parameters were calculated. From the detailed analysis it is concluded that the compressor is probably hub critical (stall initiates at the hub) in the latter stages for the design speed conditions. As a result, the casing treatment over the blade tips has little or no effect on stall margin at design speed. Radial inlet distortion did not appear to change the flow in the stages that control stall because of the rapid attenuation of the distortion within the compressor.

Hager, R. D.↗

NASA NTRS1998

A Monte Carlo Uncertainty Analysis of Ozone Trend Predictions in a Two Dimensional Model

We use Monte Carlo analysis to estimate the uncertainty in predictions of total O3 trends between 1979 and 1995 made by the Goddard Space Flight Center (GSFC) two-dimensional (2D) model of stratospheric photochemistry and dynamics. The uncertainty is caused by gas-phase chemical reaction rates, photolysis coefficients, and heterogeneous reaction parameters which are model inputs. The uncertainty represents a lower bound to the total model uncertainty assuming the input parameter uncertainties are characterized correctly. Each of the Monte Carlo runs was initialized in 1970 and integrated for 26 model years through the end of 1995. This was repeated 419 times using input parameter sets generated by Latin Hypercube Sampling. The standard deviation (a) of the Monte Carlo ensemble of total 03 trend predictions is used to quantify the model uncertainty. The 34% difference between the model trend in globally and annually averaged total O3 using nominal inputs and atmospheric trends calculated from Nimbus 7 and Meteor 3 total ozone mapping spectrometer (TOMS) version 7 data is less than the 46% calculated 1 (sigma), model uncertainty, so there is no significant difference between the modeled and observed trends. In the northern hemisphere midlatitude spring the modeled and observed total 03 trends differ by more than 1(sigma) but less than 2(sigma), which we refer to as marginal significance. We perform a multiple linear regression analysis of the runs which suggests that only a few of the model reactions contribute significantly to the variance in the model predictions. The lack of significance in these comparisons suggests that they are of questionable use as guides for continuing model development. Large model/measurement differences which are many multiples of the input parameter uncertainty are seen in the meridional gradients of the trend and the peak-to-peak variations in the trends over an annual cycle. These discrepancies unambiguously indicate model formulation problems and provide a measure of model performance which can be used in attempts to improve such models.

Considine, D. B.↗

NASA NTRS2014

Forced Response Analysis of a Fan with Boundary Layer Inlet Distortion

Boundary layer ingesting propulsion systems have the potential to significantly reduce fuel burn for future generations of commercial aircraft, but these systems must be designed to overcome the challenge of high dynamic stresses in fan blades due to forced response. High dynamic stresses can lead to high cycle fatigue failures. High-fidelity computational analysis of the fan aeromechanics is integral to an ongoing effort to design a boundary layer ingesting inlet and fan for a wind-tunnel test. An unsteady flow solution from a Reynoldsaveraged Navier Stokes analysis of a coupled inlet-fan system is used to calculate blade unsteady loading and assess forced response of the fan to distorted inflow. Conducted prior to the mechanical design of a fan, the initial forced response analyses performed in this study provide an early look at the levels of dynamic stresses that are likely to be encountered. For the boundary layer ingesting inlet, the distortion contains strong engine order excitations that act simultaneously. The combined effect of these harmonics was considered in the calculation of the forced response stresses. Together, static and dynamic stresses can provide the information necessary to evaluate whether the blades are likely to fail due to high cycle fatigue. Based on the analyses done, the overspeed condition is likely to result in the smallest stress margin in terms of the mean and alternating stresses. Additional work is ongoing to expand the analyses to off-design conditions, on-resonance conditions, and to include more detailed modeling of the blade structure.

Vibrational stress↗

DOE OSTI2022

The Analysis and Support for LANL Technology Testbed: High Altitude Balloon Gondola and Orion Eagle Gamma Ray Detection Instrument

The Analysis and Support for LANL Technology Testbed: High Altitude Balloon Gondola and Orion Eagle Gamma-Ray Detection Instrument is a project where Los Alamos National Laboratory (LANL), in conjunction with the University of Michigan, are looking into gamma-ray detection in a space-like environment. The purpose of the detection is to work on instruments that can assist with planetary exploration and nuclear deterrents using a High-Altitude Balloon (HAB). For this paper, the main focus is on the mechanical aspects related to the structural analysis of the gondola. It is vital to understand the structural analysis because it helps determine if the gondola can survive the changes in the environment and the different forces or shocks that it will experience when the parachute is deployed and at landing. Additionally, NASA's Balloon Program Office (BPO) has specific requirements to meet to achieve flight certifications. Part of those certifications requires the gondola to simulate different scenarios at specific conditions with various forces applied and orientations. Therefore, it is necessary to adjust the calculations of each method and the force applied to get the right results. Once that is done, the information is inputted into the structural analysis modeling system. The system will aid the mechanical team in assessing the conditions of the gondola at each scenario and determine its margin of safety. After the gondola and instrument are certified, the payload will be able to fly. The HAB and payload are launched from NASA's Columbia Scientific Balloon Facility (CSBF) at Fort Sumner, NM. The official flight was on September 26th, 2021.

42 ENGINEERING↗

NASA NTRS1977

Electron energy distributions in uranium helium mixtures

The high energy portion of the electron energy distribution for mixtures of uranium and helium at 1 atm, 5000 K, and a neutron flux of 2x10 to the 12th power/sq cm-sec have been calculated. The addition of He improves the heat transport characteristics of the plasma and has the feature that the He energy levels lie in the high energy portion of the electron distribution, potentially allowing non maxwellian excitation. It is concluded, however, that the resulting reaction rates are marginal relative to achieving inversion in He.

Makowski, M. A.↗

NASA NTRS1978

The cost of applying current helicopter external noise reduction methods while maintaining realistic vehicle performance

Analytical methods were developed and/or adopted for calculating helicopter component noise, and these methods were incorporated into a unified total vehicle noise calculation model. Analytical methods were also developed for calculating the effects of noise reduction methodology on helicopter design, performance, and cost. These methods were used to calculate changes in noise, design, performance, and cost due to the incorporation of engine and main rotor noise reduction methods. All noise reduction techniques were evaluated in the context of an established mission performance criterion which included consideration of hovering ceiling, forward flight range/speed/payload, and rotor stall margin. The results indicate that small, but meaningful, reductions in helicopter noise can be obtained by treating the turbine engine exhaust duct. Furthermore, these reductions do not result in excessive life cycle cost penalties. Currently available main rotor noise reduction methodology, however, is shown to be inadequate and excessively costly.

Bowes, M. A.↗

NASA NTRS1995

Flight assessment of the onboard propulsion system model for the Performance Seeking Control algorithm on an F-15 aircraft

Performance Seeking Control (PSC), an onboard, adaptive, real-time optimization algorithm, relies upon an onboard propulsion system model. Flight results illustrated propulsion system performance improvements as calculated by the model. These improvements were subject to uncertainty arising from modeling error. Thus to quantify uncertainty in the PSC performance improvements, modeling accuracy must be assessed. A flight test approach to verify PSC-predicted increases in thrust (FNP) and absolute levels of fan stall margin is developed and applied to flight test data. Application of the excess thrust technique shows that increases of FNP agree to within 3 percent of full-scale measurements for most conditions. Accuracy to these levels is significant because uncertainty bands may now be applied to the performance improvements provided by PSC. Assessment of PSC fan stall margin modeling accuracy was completed with analysis of in-flight stall tests. Results indicate that the model overestimates the stall margin by between 5 to 10 percent. Because PSC achieves performance gains by using available stall margin, this overestimation may represent performance improvements to be recovered with increased modeling accuracy. Assessment of thrust and stall margin modeling accuracy provides a critical piece for a comprehensive understanding of PSC's capabilities and limitations.

Orme, John S.↗

NASA NTRS1995

PSC algorithm description

An overview of the performance seeking control (PSC) algorithm and details of the important components of the algorithm are given. The onboard propulsion system models, the linear programming optimization, and engine control interface are described. The PSC algorithm receives input from various computers on the aircraft including the digital flight computer, digital engine control, and electronic inlet control. The PSC algorithm contains compact models of the propulsion system including the inlet, engine, and nozzle. The models compute propulsion system parameters, such as inlet drag and fan stall margin, which are not directly measurable in flight. The compact models also compute sensitivities of the propulsion system parameters to change in control variables. The engine model consists of a linear steady state variable model (SSVM) and a nonlinear model. The SSVM is updated with efficiency factors calculated in the engine model update logic, or Kalman filter. The efficiency factors are used to adjust the SSVM to match the actual engine. The propulsion system models are mathematically integrated to form an overall propulsion system model. The propulsion system model is then optimized using a linear programming optimization scheme. The goal of the optimization is determined from the selected PSC mode of operation. The resulting trims are used to compute a new operating point about which the optimization process is repeated. This process is continued until an overall (global) optimum is reached before applying the trims to the controllers.

Nobbs, Steven G.↗

NASA NTRS2006

Variations in the Sea Ice Edge and the Marginal Ice Zone on Different Spatial Scales as Observed from Different Satellite Sensor

The Marginal sea Ice Zone (MIZ) and the sea ice edge are the most dynamic areas of the sea ice cover. Knowledge of the sea ice edge location is vital for routing shipping in the polar regions. The ice edge is the location of recurrent plankton blooms, and is the habitat for a number of animals, including several which are under severe ecological threat. Polar lows are known to preferentially form along the sea ice edge because of induced atmospheric baroclinicity, and the ice edge is also the location of both vertical and horizontal ocean currents driven by thermal and salinity gradients. Finally, sea ice is both a driver and indicator of climate change and monitoring the position of the ice edge accurately over long time periods enables assessment of the impact of global and regional warming near the poles. Several sensors are currently in orbit that can monitor the sea ice edge. These sensors, though, have different spatial resolutions, different limitations, and different repeat frequencies. Satellite passive microwave sensors can monitor the ice edge on a daily or even twice-daily basis, albeit with low spatial resolution - 25 km for the Special Sensor Microwave Imager (SSM/I) or 12.5 km for the Advanced Microwave Scanning Radiometer (AMSR-E). Although special methods exist that allow the detection of the sea ice edge at a quarter of that nominal resolution (PSSM). Visible and infrared data from the Advanced Very High Resolution Radiometer (AVHRR) and from the Moderate Resolution Imaging Spectroradiometer (MODIS) provide daily coverage at 1 km and 250 m, respectively, but the surface observations me limited to cloud-free periods. The Landsat 7 Enhanced Thematic Mapper (ETM+) has a resolution of 15 to 30 m but is limited to cloud-free periods as well, and does not provide daily coverage. Imagery from Synthetic Aperture Radar (SAR) instruments has resolutions of tens of meters to 100 m, and can be used to distinguish open water and sea ice on the basis of surface and volume scattering characteristics. The Canadian RADARSAT C-band SAR provides data that cover the Arctic Ocean and the MIZ every 3 days. A change-point detection approach was utilized to obtain an ice edge estimate from the RADARSAT data The Quickscat scatterometer provides ice edge information with a resolution of a few kilometers on a near-daily basis. During portions of March and April of 2003 a series of aircraft flights were conducted over the ice edge in the Bering Sea carrying the Polarimetric Scanning Radiometer (PSR), which provides spectral coverage identical with the AMSR-E instrument at a resolution of 500 meters. In this study we investigated these different data sets and analyzed differences in their definition of the sea ice edge and the marginal ice zone and how these differences as well as their individual limitations affect the monitoring of the ice edge dynamics. We also examined how the nature of the sea ice edge, including its location, compactness and shape, changes over the seasons. Our approach was based on calculation of distances between ice edges derived from the satellite and aircraft data sets listed above as well as spectral coherence methods and shape parameters such as tortuosity, curvature, and fractional dimension.

Markus, Thorsten↗

NASA NTRS2019

A Novel Reconfigurable GaN Based Fully Solid-State Microwave Power Module for Communications/Radar Applications

In this paper, we present as a proof-of-concept (POC) the design, integration, and performance of a novel reconfigurable GaN based fully solid-state MPM. The module synergistically integrates diplexers, pre-amplifiers, multistage medium power amplifiers (MPAs), SPDT switches, and CW/Pulsed high power amplifiers (HPAs) with a voltage sequencer, a DC blanking controller, and a low voltage electronic power conditioner. The POC SSMPM operates at both S-band and X-band to serve multiple roles. The SSMPM can be reconfigured to deliver Psat of 39 dBm (8 W CW) at S-band, Psat of 46 dBm (40 W CW) at X-band, and Psat of >50 dBm (>100 W Pulsed) at X-band. Lastly, our link budget calculations indicate that the SSMPM with Psat = 40 W when coupled to a 10 cm X-band transmit antenna on a low Earth orbiting (900 km) satellite can close a 1 Gbps (QPSK) data downlink to a 1 m receive antenna on ground with 3 dB margin.

Simons, Rainee N.↗

DOE OSTI2022

Progress on Demonstration of a MOOSE-Based Coupled Capability for Hot Channel Factors in Fast Reactors

Hot channel factors (HCFs) are computed values that account for the impact on predicted peak fuel, cladding, and coolant temperatures due to uncertainties in the as-built reactor’s material properties and geometry as well as uncertainties due to modeling approximations. Reduction in computed HCF values via reduction or elimination of modeling approximations may translate to significant economic savings if the reactor power can be raised due to the extra temperature margin gained. While limited historical datasets exist for sodium-cooled fast reactors (SFRs), there are no available HCF data for lead-cooled fast reactors (LFRs) outside of work generated previously within NEAMS. The computation of HCFs involves insights from reactor physics, thermal fluids and heat conduction calculations to determine how the peak temperatures respond to various uncertainties in the design. Due to the significant advantages for multi-physics coupling offered by the MOOSE framework, Griffin (MOOSE-based reactor physics code), MOOSE Heat Conduction Module, and Cardinal (MOOSE-wrapped multi-physics application which includes the NekRS thermal fluids code) are being coupled together using the MOOSE MultiApp System to develop a highfidelity multi-physics modeling capability for HCF simulations. This high-fidelity coupling workflow may also be beneficial for other fast reactor applications in the future. In previous work, Griffin and NekRS were individually assessed to ensure the necessary capabilities were in place. This work describes initial efforts to couple the codes (including folding in the MOOSE Heat Conduction Module) and determining the workflow for the perturbed calculations which will leverage the Stochastic Tools Module (STM). To our knowledge, this is the first coupling of Griffin and NekRS as well as the first exploratory use of Stochastic Tools Module for Cardinal. In this report, the neutronics code Griffin, the heat conduction solver in MOOSE, and the MOOSE-wrapped application containing NekRS (Cardinal) are linked together to demonstrate the coupled capability. Griffin and Cardinal are linked dynamically by specifying shared libraries. Different coupling hierarchies are tested for selecting the most appropriate coupling strategy. A coupling scheme is selected based on the efficiency of calculation and ease of data communication. Multiple tests are performed to choose suitable mesh structure, model configurations, scheme setup and boundary conditions to avoid loss of energy due to data interpolation between different modules or weak imposition of fluxes in finite element codes. Computational experiments are performed to study the tolerance control of each type of iteration to avoid false convergence. The coupled capability is demonstrated in both single pin and 7-pin models based on LFR materials and geometry. The study finds that the use of too large a time step size in the heat conduction module can lead to temperature oscillation even though the heat conduction equation does not have a time-derivative kernel, but only the time-dependent boundary condition. A 7-pin model without duct region achieved good convergence in the coupled calculation while a 7 pin model with duct region experienced data communication issues which need to be resolved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1987

Improved techniques for predicting spacecraft power

Radioisotope Thermoelectric Generators (RTGs) are going to supply power for the NASA Galileo and Ulysses spacecraft now scheduled to be launched in 1989 and 1990. The duration of the Galileo mission is expected to be over 8 years. This brings the total RTG lifetime to 13 years. In 13 years, the RTG power drops more than 20 percent leaving a very small power margin over what is consumed by the spacecraft. Thus it is very important to accurately predict the RTG performance and be able to assess the magnitude of errors involved. The paper lists all the error sources involved in the RTG power predictions and describes a statistical method for calculating the tolerance.

Chmielewski, A. B.↗

NASA NTRS1988

Concentration, temperature, and density in a hydrogen-air flame by excimer-induced Raman scattering

Single-pulse, vibrational Raman scattering (VRS) is an attractive laser diagnostic for the study of supersonic hydrogen-air combustion. The VRS technique gives a complete thermodynamic description of the gas mixture at a point in the reacting flow. Single-pulse, vibrational Raman scattering can simultaneously provide independent measurements of density, temperature, and concentration of each major species (H2, H2O, O2 and N2) in a hydrogen/air turbulent combustor. Also the pressure can be calculated using the ideal gas law. However, single-pulse VRS systems in current use for measurement of turbulent combustion have a number of shortcomings when applied to supersonic flows: (1) slow repetition rate (1 to 5 Hz), (2) poor spatial resolution (0.5x0.3x0.3 cu mm), and (3) marginal time resolution. Most of these shortcomings are due to the use of visible wavelength flash-lamp pumped dye lasers. The advent of UV excimer laser allows the possibility of dramatic improvements in the single-pulse, vibrational Raman scattering. The excimer based VRS probe will greatly improve repetition rate (100 to 500 Hz), spatial resolution (0.1x0.1x0.1 cu mm) and time resolution (30ns). These improvements result from the lower divergence of the UV excimer, higher repetition rate, and the increased Raman cross-sections (15 to 20 times higher) at ultra-violet (UV) wavelengths. With this increased capability, single-pulse vibrational Raman scattering promises to be an ideal non-intrusive probe for the study of hypersonic propulsion flows.

Wehrmeyer, Joseph A.↗

NASA NTRS2013

Orion MPCV Touchdown Detection Threshold Development and Testing

A robust method of detecting Orion Multi ]Purpose Crew Vehicle (MPCV) splashdown is necessary to ensure crew and hardware safety during descent and after touchdown. The proposed method uses a triple redundant system to inhibit Reaction Control System (RCS) thruster firings, detach parachute risers from the vehicle, and transition to the post ]landing segment of the Flight Software (FSW). The vehicle crew is the prime input for touchdown detection, followed by an autonomous FSW algorithm, and finally a strictly time based backup timer. RCS thrusters must be inhibited before submersion in water to protect against possible damage due to firing these jets under water. In addition, neglecting to declare touchdown will not allow the vehicle to transition to post ]landing activities such as activating the Crew Module Up ]righting System (CMUS), resulting in possible loss of communication and difficult recovery. A previous AIAA paper gAssessment of an Automated Touchdown Detection Algorithm for the Orion Crew Module h concluded that a strictly Inertial Measurement Unit (IMU) based detection method using an acceleration spike algorithm had the highest safety margins and shortest detection times of other methods considered. That study utilized finite element simulations of vehicle splashdown, generated by LS ]DYNA, which were expanded to a larger set of results using a Kriging surface fit. The study also used the Decelerator Systems Simulation (DSS) to generate flight dynamics during vehicle descent under parachutes. Proto ]type IMU and FSW MATLAB models provided the basis for initial algorithm development and testing. This paper documents an in ]depth trade study, using the same dynamics data and MATLAB simulations as the earlier work, to further develop the acceleration detection method. By studying the combined effects of data rate, filtering on the rotational acceleration correction, data persistence limits and values of acceleration thresholds, an optimal configuration was determined. The lever arm calculation, which removes the centripetal acceleration caused by vehicle rotation, requires that the vehicle angular acceleration be derived from vehicle body rates, necessitating the addition of a 2nd order filter to smooth the data. It was determined that using 200 Hz data directly from the vehicle IMU outperforms the 40 Hz FSW data rate. Data persistence counter values and acceleration thresholds were balanced in order to meet desired safety and performance. The algorithm proved to exhibit ample safety margin against early detection while under parachutes, and adequate performance upon vehicle splashdown. Fall times from algorithm initiation were also studied, and a backup timer length was chosen to provide a large safety margin, yet still trigger detection before CMUS inflation. This timer serves as a backup to the primary acceleration detection method. Additionally, these parameters were tested for safety on actual flight test data, demonstrating expected safety margins.

Daum, Jared↗

NASA NTRS1985

Self-Tuning Adaptive-Controller Using Online Frequency Identification

A real time adaptive controller was designed and tested successfully on a fourth order laboratory dynamic system which features very low structural damping and a noncolocated actuator sensor pair. The controller, implemented in a digital minicomputer, consists of a state estimator, a set of state feedback gains, and a frequency locked loop (FLL) for real time parameter identification. The FLL can detect the closed loop natural frequency of the system being controlled, calculate the mismatch between a plant parameter and its counterpart in the state estimator, and correct the estimator parameter in real time. The adaptation algorithm can correct the controller error and stabilize the system for more than 50% variation in the plant natural frequency, compared with a 10% stability margin in frequency variation for a fixed gain controller having the same performance at the nominal plant condition. After it has locked to the correct plant frequency, the adaptive controller works as well as the fixed gain controller does when there is no parameter mismatch. The very rapid convergence of this adaptive system is demonstrated experimentally, and can also be proven with simple root locus methods.

Chiang, W. W.↗

NASA NTRS2011

AGN Obscuration Through Dusty Infrared Dominated Flows: Radiation-Hydrodynamics Solution for the Wind - 1

We construct a radiation-hydrodynamics model for the obscuring toroidal structure in active galactic nuclei. In this model the obscuration is produced at parsec scale by a dense, dusty wind which is supported by infrared radiation pressure on dust grains. To find the distribution of radiation pressure, we numerically solve the 2D radiation transfer problem in a flux limited diffusion approximation. We iteratively couple the solution with calculations of stationary 1D models for the wind, and obtain the z-component of the velocity. Our results demonstrate that for AGN luminosities greater than 0.1 L(sub edd) external illumination can support a geometrically thick obscuration via outflows driven by infrared radiation pressure. The terminal velocity of marginally Compton-thin models (0.2 < tau(sub T) < 0.6), is comparable to or greater than the escape velocity. In Compton thick models the maximum value of the vertical component of the velocity is lower than the escape velocity, suggesting that a significant part of our torus is in the form of failed wind. The results demonstrate that obscuration via normal or failed infrared-driven winds is a viable option for the AGN torus problem and AGN unification models. Such winds can also provide an important channel for AGN feedback.

Dorodnitsyn, A.↗