Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TIME SCALE”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Numerical Prediction of Cyclic Variability in a Spark Ignition Engine Using a Parallel Large Eddy Simulation Approach

Cycle-to-cycle variability (CCV) is detrimental to IC engine operation and can lead to partial burn, misfire, and knock. Predicting CCV numerically is extremely challenging due to two key reasons. First, high-fidelity methods such as large eddy simulation (LES) are required to accurately resolve the in-cylinder turbulent flow field both spatially and temporally. Second, CCV is experienced over long timescales and hence the simulations need to be performed for hundreds of consecutive cycles. Ameen et al. (2017, "Parallel Methodology to Capture Cyclic Variability in Motored Engines,"Int. J. Engine Res., 18(4), pp. 366-377.) developed a parallel perturbation model (PPM) approach to dissociate this long time-scale problem into several shorter time-scale problems. This strategy was demonstrated for motored engine and it was shown that the mean and variance of the in-cylinder flow field was captured reasonably well by this approach. In the present study, this PPM approach is extended to simulate the CCV in a fired port-fuel injected (PFI) spark ignition (SI) engine. The predictions from this approach are also shown to be similar to the consecutive LES cycles. It is shown that the parallel approach is able to predict the coefficient of variation (COV) of the in-cylinder pressure and burn raterelated parameters with sufficient accuracy, and is also able to predict the qualitative trends in CCV with changing operating conditions. It is shown that this new approach is able to give accurate predictions of the CCV in fired engines in less than one-tenth of the time required for the conventional approach of simulating consecutive engine cycles.

Ameen, Muhsin M.↗
NASA NTRS1996

Paleobiology, community ecology, and scales of ecological pattern

The fossil record provides a wealth of data on the role of regional processes and historical events in shaping biological communities over a variety of time scales. The Quaternary record with its evidence of repeated climatic change shows that both terrestrial and marine species shifted independently rather than as cohesive assemblages over scales of thousands of years. Larger scale patterns also show a strong individualistic component to taxon dynamics; assemblage stability, when it occurs, is difficult to separate from shared responses to low rates of environmental change. Nevertheless, the fossil record does suggest that some biotic interactions influence large-scale ecological and evolutionary patterns, albeit in more diffuse and protracted fashions than those generally studied by community ecologists. These include: (1) the resistance by incumbents to the establishment of new or invading taxa, with episodes of explosive diversification often appearing contingent on the removal of incumbents at extinction events; (2) steady states of within-habitat and global diversity at longer time scales (10(7)-l0(8) yr), despite enormous turnover of taxa; and (3) morphological and biogeographic responses to increased intensities of predation and substratum disturbance over similarly long time scales. The behavior of species and communities over the array of temporal and spatial scales in the fossil record takes on additional significance for framing conservation strategies, and for understanding recovery of species, lineages, and communities from environmental changes.

NASA Discipline Exobiology↗
NASA NTRS1981

Orbital inclination of Iapetus and the rotation of the Laplacian plane

This paper explores the possibility that the orbit of Iapetus, with its relatively large inclination but small eccentricity, was generated by a rapid dispersal of a gaseous circumplanetary disk, assumed to be the progenitor of the satellite system. The orientation of the local Laplacian plane is shown to be a sensitive function of the disk's structure. Modification of the disk on a time scale comparable to the precession of the orbit's nodal line can produce a large inclination from one that is initially zero, while leaving the eccentricity unchanged. This time scale is of the same order of magnitude as the viscous evolution time scale for a fully turbulent disk. Hence, Iapetus need not be a captured satellite to account for its curious orbital signature.

Ward, W. R.↗
DOE OSTI2026

Dynamic Nanoscale Spatial Heterogeneity in a Perovskite-to-Brownmillerite Topotactic Phase Transformation

Phase transitions are omnipresent in modern condensed matter physics and its applications. In solids, first-order phase transformations typically occur by nucleation and growth under nonequilibrium conditions. Under constant external conditions, e.g., constant annealing temperature and pressure, the nucleation and growth dynamics are often thought of as spatially and temporally independent. Here, in situ Bragg X-ray photon correlation spectroscopy (XPCS) reveals nanoscale spatial and dynamical heterogeneity in the perovskite-to-brownmillerite topotactic phase transformation in La 0.7 Sr 0.3 CoO 3 thin films annealed under constant reducing conditions over a time span of multiple hours. Specifically, a time scale associated with domain growth remains stable, with a corresponding domain wall speed of v d = 6 ± 0.5 × 10 –4 nm/s (2 ± 0.2 nm/h), while a slower time scale, associated with temperature-driven depinning of domains, leads to accelerating dynamics with time scales following an aging power law with exponent −2.2 ± 0.5. This experiment demonstrates that Bragg XPCS is a powerful tool to study nanoscale dynamics in structural phase transformations, with the ability to extract quantitative average values related to nanodomain motion in situ. Furthermore, the results are relevant for phase engineering of phase-change devices, as they show that nanoscale dynamics, linked to domain and domain-wall motion, can continuously evolve and speed up with time, even hours after the initiation of the phase transformation, with potential repercussions on electrical performance.

X-ray photon correlation spectroscopy↗
DOE OSTI2016

Supramolecular Double-Helix Formation by Diastereoisomeric Conformations of Configurationally Enantiomeric Macrocycles

Solid-state superstructures, resulting from assemblies programmed by homochirality, are attracting considerable attention. In addition, artificial double-helical architectures are being investigated, especially in relation to the ways in which homochiral small molecules can be induced to yield helical forms as a result of chiral induction. In this work, we report the highly specific self-assembly upon crystallization of a double-helical superstructure from an enantiopure macrocyclic dimer which adopts two diastereoisomeric conformations in a molar ratio of 1.5:1 in dimethyl sulfoxide. These two conformational diastereoisomers self-organize—and self-sort—in the crystalline phase in equimolar proportions to form two single-handed helices which are complementary to each other, giving rise to the assembly of a double helix that is stabilized by intermolecular [C–H···O] and π–π stacking interactions. The observed self-sorting phenomenon occurs on going from a mixed-solvent system containing two equilibrating conformational diastereoisomers, presumably present in unequal molar proportions, into the solid state. The diastereoisomeric conformations are captured upon crystallization in a 1:1 molar ratio in the double-helical superstructure, whose handedness is dictated by the choice of the enantiomeric macrocyclic dimer. The interconversion of the two conformational diastereoisomers derived from each configurationally enantiomeric macrocycle was investigated in CD 3 SOCD 3 solution by variable-temperature 1 H NMR spectroscopy (VT NMR) and circular dichroism (VT CD). The merging of the resonances for the protons corresponding to the two diastereoisomers at a range of coalescence temperatures in the VT NMR spectra and occurrence of the isosbestic points in the VT CD spectra indicate that the two diastereoisomers are interconverting slowly in solution on the 1 H NMR time scale but rapidly on the laboratory time scale. To the best of our knowledge, the self-assembly of such solid-state superstructures from two conformational diastereoisomers of a homochiral macrocycle is a rare, if not unique, occurrence.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Formation of turing patterns in strongly magnetized electric discharges

Abstract Pattern formation and self-organization in many biological and non-biological systems can be explained through Turing’s activator-inhibitor model. Here we show how this model can be employed to describe the formation of filamentary structures in a low-pressure electric discharge exposed to a strong magnetic field. Theoretical investigation reveals that the fluid equations describing a magnetized plasma can be rearranged to take the mathematical form of Turing’s activator-inhibitor model. Numerical simulations based on the equations derived from this approach could reproduce the various patterns observed in the experiments. Also, it is shown that a density imbalance between electrons and ions exists in the bulk of the magnetized plasma that generates an electric field structure transverse to the applied magnetic field. This electric field is responsible for the stability of the filamentary patterns in the magnetized plasma over time scales much longer than the characteristic time scales of the electric discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Real-time power management technique for microgrid with flexible boundaries

In order to diminish the impacts brought by high penetration of renewable energy on the reliability of distribution systems, some distribution networks (e.g. Chattanooga electric power board system) have deployed smart switches (SSs) to island some areas to mitigate outage losses. However, due to intermittency and sharply changing rate of renewable energy, it is likely to experience insufficient or excessive power for islanded areas. Therefore, a microgrid controller featured with flexible boundaries is proposed. With proposed microgrid controller, the microgrid can not only shrink or expand its boundaries according to current renewable energy supply, but also disconnect/connect to the main grid with a designated SS. Furthermore, to ensure the microgrid controller could obtain suitable boundaries on the time scale of seconds, a real-time power management technique with alternative generating algorithm is designed to generate all possible alternative boundaries and choose the optimal one, which is scalable to any topology. In addition, in order to maintain state of charge of batteries within a desirable range, anti-overcharge/discharge strategies are designed. Four comprehensive experiments verify that the implementation of the microgrid controllers can realise flexible boundaries and deal with sharply changing rate of renewable generation or load on the time scale of seconds.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2011

Fluctuation driven transport and stationary profiles

Transport equations for particles and energy can be derived when the fluctuations conserve adiabatic invariants. The transport equations determine both stationary density and pressure profiles and the direction of the turbulence-driven fluxes which can be inward or outward. An inward turbulent pinch is predicted which creates stationary profiles and reverses direction depending on the density and temperature gradients. The transport fluxes are independent of the underlying drive that leads to plasma turbulence. For low frequency turbulence, the formulation remains correct when the collisional time scale is faster than the confinement time scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Ponderomotive force driven density modifications parallel to B on the LAPD

The ponderomotive force has previously been identified as a possible driver of observed density modifications close to radio frequency actuators during operation. This nonlinear force redistributes density in regions of gradients in the magnitude of an oscillating electric field and describes the influence of the fast time scale RF wave dynamics on slow time scale plasma transport. Depletion of the saturation current (a proxy for the density) measured at the Large Plasma Device (LAPD) was 30–35% during ion cyclotron range of frequencies operation. A coupled 1D plasma transport and cold plasma frequency domain wave solver was developed to self-consistently describe ponderomotive effects and was used to compare with results obtained from the LAPD experiment. The scaled current density driver for the wave model yielded an RF B field in close agreement with two components of the experimental data. However, the 1D parallel model did not accurately reproduce the amplitude or spatial distribution observed in experimental measurements of By. Within the limitations of the 1D model, initial simulation results showed that the ponderomotive force depleted up to 8% for high power (1 MW) and around 1% for the experimental power of 120 kW. This could suggest that the ponderomotive force is not the main driver of density modification for the LAPD experiments presented in this paper. Higher fidelity tools of at least 2D will be required to give a more realistic description of the RF E fields and the effect of the ponderomotive force on the LAPD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Explaining Cold-Pulse Dynamics in Tokamak Plasmas Using Local Turbulent Transport Models

A long-standing enigma in plasma transport has been resolved by modeling of cold-pulse experiments conducted on the Alcator C-Mod tokamak. Controlled edge cooling of fusion plasmas triggers core electron heating on time scales faster than an energy confinement time, which has long been interpreted as strong evidence of nonlocal transport. Here, this Letter shows that the steady-state profiles, the cold-pulse rise time, and disappearance at higher density as measured in these experiments are successfully captured by a recent local quasilinear turbulent transport model, demonstrating that the existence of nonlocal transport phenomena is not necessary for explaining the behavior and time scales of cold-pulse experiments in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Putting multipoint kinetics to the test: preliminary analysis of the MUSE-4 experiment with APOLLO3

Modeling neutron kinetics is a challenging task and over the years, several approaches have been developed offering different tradeoffs between precision and cost. Probably the most popular one is point kinetics (PK), which, despite being very cheap, provides good results for many applications. One of the limitations of PK is that it is not able to model effects that happen before the prompt jump time scale. In practice, applications, where those time scales are interesting, are rare and those for which experimental data are available are even rarer. One of those cases is the MUSE-4 experiment for which the focus is on the very fast behavior and where we have no thermal feedback. As we have shown in a recent paper, a simple extension of PK, known as multipoint kinetics (MPK), can overcome this limitation. In this paper, which is intended as a preliminary study, we show how a simple MPK model can qualitatively reproduce the fast behavior that we cannot model with PK. Moreover, we will show how to exploit its low cost and compact form to model also slow effects as the precursors built up, behaviours that are not affordable to model with more expensive techniques, such as kinetic Monte Carlo simulations. For the periodic pulse transients with no feedbacks we will also provide an analytical solution in closed form. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2010

Fundamental Science-Based Simulation of Nuclear Waste Forms

This report presents a hierarchical multiscale modeling scheme based on two-way information exchange. To account for all essential phenomena in waste forms over geological time scales, the models have to span length scales from nanometer to kilometer and time scales from picoseconds to millenia. A single model cannot cover this wide range and a multi-scale approach that integrates a number of different at-scale models is called for. The approach outlined here involves integration of quantum mechanical calculations, classical molecular dynamics simulations, kinetic Monte Carlo and phase field methods at the mesoscale, and continuum models. The ultimate aim is to provide science-based input in the form of constitutive equations to integrated codes. The atomistic component of this scheme is demonstrated in the promising waste form xenotime. Density functional theory calculations have yielded valuable information about defect formation energies. This data can be used to develop interatomic potentials for molecular dynamics simulations of radiation damage. Potentials developed in the present work show a good match for the equilibrium lattice constants, elastic constants and thermal expansion of xenotime. In novel waste forms, such as xenotime, a considerable amount of data needed to validate the models is not available. Integration of multiscale modeling with experimental work is essential to generate missing data needed to validate the modeling scheme and the individual models. Density functional theory can also be used to fill knowledge gaps. Key challenges lie in the areas of uncertainty quantification, verification and validation, which must be performed at each level of the multiscale model and across scales. The approach used to exchange information between different levels must also be rigorously validated. The outlook for multiscale modeling of wasteforms is quite promising.

Continuum Model↗
DOE OSTI2017

The Calibration System of the E989 Experiment at Fermilab

The muon anomaly aµ is one of the most precise quantity known in physics experimentally and theoretically. The high level of accuracy permits to use the measurement of aµ as a test of the Standard Model comparing with the theoretical calculation. After the impressive result obtained at Brookhaven National Laboratory in 2001 with a total accuracy of 0.54 ppm, a new experiment E989 is under construction at Fermilab, motivated by the diff of aexp SM µ − aµ ∼ 3σ. The purpose of the E989 experiment is a fourfold reduction of the error, with a goal of 0.14 ppm, improving both the systematic and statistical uncertainty. With the use of the Fermilab beam complex a statistic of × 21 with respect to BNL will be reached in almost 2 years of data taking improving the statistical uncertainty to 0.1 ppm. Improvement on the systematic error involves the measurement technique of ωa and ωp, the anomalous precession frequency of the muon and the Larmor precession frequency of the proton respectively. The measurement of ωp involves the magnetic field measurement and improvements on this sector related to the uniformity of the field should reduce the systematic uncertainty with respect to BNL from 170 ppb to 70 ppb. A reduction from 180 ppb to 70 ppb is also required for the measurement of ωa; new DAQ, a faster electronics and new detectors and calibration system will be implemented with respect to E821 to reach this goal. In particular the laser calibration system will reduce the systematic error due to gain fl of the photodetectors from 0.12 to 0.02 ppm. The 0.02 ppm limit on systematic requires a system with a stability of 10−4 on short time scale (700 µs) while on longer time scale the stability is at the percent level. The 10−4 stability level required is almost an order of magnitude better than the existing laser calibration system in particle physics, making the calibration system a very challenging item. In addition to the high level of stability a particular environment, due to the presence of a 14 m diameter storage ring, a highly uniform magnetic field and the detector distribution around the storage ring, set specific guidelines and constraints. This thesis will focus on the final design of the Laser Calibration System developed for the E989 experiment. Chapter 1 introduces the subject of the anomalous magnetic moment of the muon; chapter 2 presents previous measurement of g-2, while chapter 3 discusses the Standard Model prediction and possible new physics scenario. Chapter 4 describes the E989 experiment. In this chapter will be described the experimental technique and also will be presented the experimental apparatus focusing on the improvements necessary to reduce the statistical and systematic errors. The main item of the thesis is discussed in the last two chapters: chapter 5 is focused on the Laser Calibration system while chapter 6 describes the Test Beam performed at the Beam Test Facility of Laboratori Nazionali di Frascati from the 29th February to the 7th March as a final test for the full calibrations system. An introduction explain the physics motivation of the system and the diff t devices implemented. In the final chapter the setup used will be described and some of the results obtained will be presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Explaining Cold-Pulse Dynamics in Tokamak Plasmas Using Local Turbulent Transport Models

A long-standing enigma in plasma transport has been resolved by modeling of cold-pulse experiments conducted on the Alcator C-Mod tokamak. Controlled edge cooling of fusion plasmas triggers core electron heating on time scales faster than an energy confinement time, which has long been interpreted as strong evidence of nonlocal transport. This Letter shows that the steady-state profiles, the cold-pulse rise time, and disappearance at higher density as measured in these experiments are successfully captured by a recent local quasilinear turbulent transport model, demonstrating that the existence of nonlocal transport phenomena is not necessary for explaining the behavior and time scales of cold-pulse experiments in tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1999

Dynamics and Structure of Weakly-Strained Flames In Normal- and Micro-Gravity

Strained laminar flames have been systematically studied, as the understanding of their structure and dynamic behavior is of relevance to turbulent combustion. Most of these studies have been conducted in opposed-jet, stagnation-type flow configurations. Flame studies in stagnation flows also allow for the determination of fundamental flame properties such as laminar flame speeds and extinction states that can be also conveniently modeled in detail. Studies at high strain rates are important in quantifying and understanding the response of vigorously-burning flames under conditions in which the transport time scales become comparable to the chemical time scales. Studies of weakly-strained flames can be of particular interest for all stoichiometries. For example, the laminar flame speeds for any equivalence ratio, phi, can be accurately determined by using the counterflow technique only if measurements are obtained at very low strain rates. Furthermore, near-limit flames can be only stabilized by weak strain rates. Previous studies have shown that weakly-burning flames are particularly sensitive to chain mechanisms, thermal radiation, and unsteadiness. The stabilization and study of weakly-strained flames is complicated by the presence of buoyancy that can render the flames unstable to the point of extinction. Such instabilities are caused either by the induced natural convection or the fact that higher density fluid can find itself on top of a lower density fluid. Thus, the use of microgravity (mu-g) becomes essential in order to provide meaningful insight into this important combustion regime. In view of the foregoing considerations, the main objectives of the program are to: (1) Experimentally determine the laminar flame speed at near-zero strain rates; (2) Experimentally determine the extinction limits of near-limit flames; (3) Experimentally determine the response of near-limit flames to unsteadiness and heat loss; (4) Introduce Digital Particle Image Velocimetry (DPIV) in microgravity; (5) Conduct detailed numerical simulations of the experiments; and (6) Provide physical insight into the underlying physico-chemical mechanisms.

Zhang, Hai↗
NASA NTRS1998

Four Years of EUVE Observations of the Bursting Gamma-Ray Emitting Blazar Markarian 421

We present spectral and timing analysis of all the data collected by the Extreme Ultraviolet Explorer Satellite (EUVE) for the bright, nearby BL Lacertae object, Markarian 421, during the four-year period 1994-1997. During these years Mrk 421 has been observed by EUVE 4 times with the Deep-Survey/Spectrograph and 2 times with the imaging telescopes for a total of approximately 1.4 7nillions seconds. From 1993 to 1996 three very bright gamma ray flares were also detected by the Whipple observatory. In 1994 Mrk 421 was observed simultaneously by EUVE (Apr 2-12) and IUE one month before the TeV flare; of the 2 EUVF, observations (Feb 4-7 and Apr 25-May 13) of 1995, the second was part of a multiwavelength campaign that mapped the evolution of the TeV flare. In 1996 we observed Mrk 421 twice simultaneously with XTE: one immediately before (Apr 17-30) and another one (May 10-11) right after the May 7 1996 TeV flare. And finally in 1997 from Feb. 7 to Feb. ll. The total light curve seems to be smoothly varying on the long time-scale while on a shorter time-scale there is evidence of an EUVE flare well correlated to the TeV energy 1995 flare. We have analysed the three spectral data set in an homogenous way using the appropriate calibration data for the off-axis observations and our analysis confirms the presence of the absorption features around approximately 70A, in the entire 1995 dataset and possibly in the 1996 data set. We also show the first power spectrum analysis of the Mrk 421 EUVE lightcurves and a comparison with the power spectra predicted by current theoretical models.

Cagnoni, Ilaria↗
NASA NTRS2000

Structure of the Circumnuclear Region of Seyfert 2 Galaxies Revealed by RXTE Hard X-Ray Observations of NGC 4945

NGC 4945 is one of the brightest Se.yfert galaxies on the sky at 100 keV, but is completely absorbed below 10 keV, implying an optical depth of the absorber to electron scattering of a few; its absorption column is probably the largest which still allows a direct view of the nucleus at hard X-ray energies. Our observations of it with the Rossi X-ray Timing Explorer (RXTE) satellite confirm the large absorption, which for a simple phenomenological fit using an absorber with Solar abundances implies a column of 4.5(sup 0.4, sub -0.4) x 10(exp 24) /sq cm. Using a a more realistic scenario (requiring Monte Carlo modeling of the scattering), we infer the optical depth to Thomson scattering of approximately 2.4. If such a scattering medium were to subtend a large solid angle from the nucleus, it should smear out any intrinsic hard X-ray variability on time scales shorter than the light travel time through it. The rapid (with a time scale of approximately a day) hard X-ray variability of NGC 4945 we observed with the RXTE implies that the bulk of the extreme absorption in this object does not originate in a parsec-size, geometrically thick molecular torus. Limits on the amount of scattered flux require that the optically thick material on parsec scales must be rather geometrically thin, subtending a half-angle < 10 deg. This is only marginally consistent with the recent determinations of the obscuring column in hard X-rays, where only a quarter of Seyfert 2s have columns which are optically thick, and presents a problem in accounting for the Cosmic X-ray Background primarily with AGN possessing the geometry as that inferred by us. The small solid angle of the obscuring material, together with the black hole mass (of approximately 1.4 x 10(exp 6) solar mass) from megamaser measurements. allows a robust determination of the source luminosity, which in turn implies that the source radiates at approximately 10% of the Eddington limit.

Madejski, G.↗