Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “global shape modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 325 records · Page 18

The NASA Micro Pulse Lidar Network (MPLNET): Introduction of the New Version 3 Release

The NASA Micro-Pulse Lidar Network (MPLNET) is a global federated network of polarized Micro-Pulse Lidar (MPL) systems running continuously. MPLNET began in 2000, and there have been over 70 sites deployed worldwide, with 24 sites currently active and a few more planned over the next year. Seven of the long-term sites have 10+ years of data, and many more have 5+ years. Most sites are co-located with AERONET providing joint data on column and vertically resolved aerosol and cloud information. This presentation will introduce our new Version 3 MPLNET data. All sites in the network now feature eye-safe polarized backscatter MPL instruments, providing information on attenuated backscatter and particle shape. In addition to change with our signal data, we have an enhanced cloud product suite, a new PBL height product, and inclusion of the new AERONET lunar aerosol optical depth into MPLNET aerosol retrievals. A new quality flag process will be used to better describe all data products. Finally, a new data portal will provide near-real-time (NRT) access to all data products, including new quality assured NRT L1.5 products. Custom products developed for model specific applications will also be provided.

Welton, Ellsworth J.↗

Evidence for a Stratified Upper Mantle Preserved within the South Pole - Aitken Basin

The evolution and compositional structure of the lunar mantle has been extensively modeled but insufficiently constrained by observations. Here, we identify and characterize mantle materials exposed by the Moon's largest impact basin to better understand the composition, stratigraphy, and evolution of the upper mantle. The vast South Pole-Aitken Basin (SPA) exhibits a broad, crescent-shaped thorium and potassium distribution. These incompatible elements are predicted to be concentrated in the dregs of the lunar magma ocean during end-stage crystallization. Through consideration of basin formation models convolved with subsequent geologic evolution, we demonstrate that the distribution and implied stratigraphy of Th- and K-bearing materials across SPA are consistent with an upper mantle ejecta origin. The most pristine exposures of these materials are confined to northwest SPA and also exhibit elevated Ti and Fe (relative to the farside highlands) in association with a gabbronoritic mineralogy. This is consistent with latestage magma ocean assemblages predicted by petrologic models. In contrast, SPA impact melt derived from greater depths is associated with a low-Ca pyroxene-dominated assemblage. Together, these compositional patterns are evidence for a stratified ancient upper mantle. Importantly, the incompatible-element-enriched, ilmenite-bearing ferroan gabbronoritic cumulates evidently had not participated in gravitational overturn at the time of SPA formation. Contrary to recent hypotheses invoking nearside sequestration of incompatible elements to explain hemispherical differences in crustal building and volcanic resurfacing, it follows that incompatible elements were globally distributed in the magma ocean at the time of SPA formation.

Lunar dichotomy↗

Pinatubo Aerosol Evolution: Using Composite Data Sets to Build the Global- To Micro-Scale Picture and Assess Consistency of Different Measurements

This paper brings together experimental evidence required to build realistic models of the global evolution of physical, chemical, and optical properties of the aerosol resulting from the 1991 Pinatubo volcanic eruption. Such models are needed to compute the effects of the aerosol on atmospheric chemistry, dynamics, radiation, and temperature. Whereas there is now a large and crowing body of post-Pinatubo measurements by a variety of techniques, some results are in conflict, and a self-consistent, unified picture is needed, along with an assessment of remaining uncertainties. This paper examines, data from photometers, radiometers, impactors, optical counter/sizers, and lidars operated on the ground, aircraft, balloons, and spacecraft. Example data sources include: (1) Tracking sunphotometers and lidars at Mauna Loa Observatory (MLO) and on the DC-8. (2) Particle spectrometers and wire impactors on the ER-2 and DC-8. (3) Dustsondes (particle counter/sizers on balloons). and (3) SAGE II, SAM II, AVHRR, CLAES, and ISAMS sensors on a variety of satellites. We assess the mutual consistency of these disparate data sets and recommend 'consensus' properties and uncertainties in the process of developing a composite data set. Recommended properties include the spatial and temporal evolution of particle chemical composition, shape, wavelength-and temperature-dependent refractive index, size distribution, and optical depth spectra. Supporting references are cited and representative data shown.

Russell, Philip B.↗

Intercomparison of numerical models of flaring coronal loops

The proposed Benchmark Problem consists of an infinitesimal magnetic flux tube containing a low-beta plasma. The field strength is assumed to be so large that the plasma can move only along the flux tube, whose shape remains invariant with time (i.e., the fluid motion is essentially one-dimensional). The flux tube cross section is taken to be constant over its entire length. In planar view the flux tube has a semi-circular shape, symmetric about its midpoint s = s sub max and intersecting the chromosphere-corona interface (CCI) perpendicularly at each foot point. The arc length from the loop apex to the CCI is 10,000 km. The flux tube extends an additional 2000 km below the CCI to include the chromosphere, which initially has a uniform temperature of 8000 K. The temperature at the top of the loop was fixed initially at 2 X 1 million K. The plasma is assumed to be a perfect gas (gamma = 5/3), consisting of pure hydrogen which is considered to be fully ionized at all temperatures. For simplicity, moreover, the electron and ion temperatures are taken to be everywhere equal at all times (corresponding to an artificially enhanced electron-ion collisional coupling). While there was more-or-less unanimous agreement as to certain global properties of the system behavior (peak temperature reached, thermal-wave time scales, etc.), no two groups could claim satisfactory accord when a more detailed comparison of solutions was attempted.

Kopp, R. A.↗

Modeling South Pacific Ice-Ocean Interactions in the Global Climate System

The objective of this project has been to improve the modeling of interactions between large Antarctic ice shelves and adjacent regions of the Southern Ocean. Our larger goal is to gain a better understanding of the extent to which the ocean controls ice shelf attrition, thereby influencing the size and dynamics of the Antarctic Ice Sheet. Melting and freezing under ice shelves also impacts seawater properties, regional upwelling and sinking and the larger-scale ocean circulation. Modifying an isopycnal coordinate general circulation model for use in sub-ice shelf cavities, we found that the abrupt change in water column thickness at an ice shelf front does not form a strong barrier to buoyancy-driven circulation across the front. Outflow along the ice shelf base, driven by melting of the thickest ice, is balanced by deep inflow. Substantial effort was focused on the Filchner-Ronne cavity, where other models have been applied and time-series records are available from instruments suspended beneath the ice. A model comparison indicated that observed changes in the production of High Salinity Shelf Water could have a major impact on circulation within the cavity. This water propagates into the cavity with an asymmetric seasonal signal that has similar phasing and shape in the model and observations, and can be related to winter production at the sea surface. Even remote parts of the sub-ice shelf cavity are impacted by external forcing on sub-annual time scales. This shows that cavity circulations and products, and therefore cavity shape, will respond to interannual variability in sea ice production and longer-term climate change. The isopycnal model gives generally lower net melt rates than have been obtained from other models and oceanographic data, perhaps due to its boundary layer formulation, or the lack of tidal forcing. Work continues on a manuscript describing the Ross cavity results.

Holland, David M.↗

Global Magnetohydrodynamic Modeling of the Solar Corona

The coronal magnetic field defines the structure of the solar corona, the position of the heliospheric current sheet, the regions of fast and slow solar wind, and the most likely sites of coronal mass ejections. There are few measurements of the magnetic fields in the corona, but the line-of-sight component of the global magnetic fields in the photosphere have been routinely measured for many years (for example, at Stanford's Wilcox Solar Observatory, and at the National Solar Observatory at Kitt Peak). The SOI/MDI instrument is now providing high-resolution full-disk magnetograms several times a day. Understanding the large-scale structure of the solar corona and inner heliosphere requires accurately mapping the measured photospheric magnetic field into the corona and outward. Ideally, a model should not only extrapolate the magnetic field, but should self-consistently reconstruct both the plasma and magnetic fields in the corona and solar wind. Support from our NASA SR&T contract has allowed us to develop three-dimensional magnetohydrodynamic (MHD) computations of the solar corona that incorporate observed photospheric magnetic fields into the boundary conditions. These calculations not only describe the magnetic field in the corona and interplanetary spice, but also predict the plasma properties as well. Our computations thus far have been successful in reproducing many aspects of both coronal and interplanetary data, including the structure of the streamer belt, the location of coronal hole boundaries, and the position and shape of the heliospheric current sheet. The most widely used technique for extrapolating the photospheric magnetic field into the corona and heliosphere are potential field models, such as the potential field source-surface model (PFSS),and the potential field current-sheet (PFCS) model

Linker, Jon A.↗

Towards Soil Moisture Profile Estimation in the Root Zone Using L- and P-Band Radiometer Observations: A Coherent Modelling Approach

Precision irrigation management and crop water stress assessment rely on accurate estimation of root zone soil moisture. However, only the top 5cm soil moisture can be estimated using the two current passive microwave satellite missions, Soil Moisture and Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP), which operate at L-band (wavelength of ~21cm). Since the contributing depth of the soil to brightness temperature increases with observation wavelength, it is expected that a P-band (wavelength of ~40cm) radiometer could potentially provide soil moisture information from deeper layers of the soil profile. Moreover, by combining both L- and P- bands, it is hypothesized that the soil moisture profile can be estimated even beyond their individual observation depths. The aim of this study was to demonstrate the potential of combined L-band and P-band radiometer observations to estimate the soil moisture profile under flat bare soil using a stratified coherent forward model. Brightness temperature observations at L-band and P-band from a tower based experimental site across a dry (April 2019) and a wet (March 2020) period, covering different soil moisture profile shapes, were used in this study. Results from an initial synthetic study showed that the performance of a combined L-band and P-band approach was better than the performance of using either band individually, with an average depth over which reliable soil moisture profile information could be estimated (i.e. with a target root mean square error (RMSE) of less than 0.04 m3/m3) being 20cm for linear and 15cm for second-order polynomial functions. Other functions were also tested but found to have a poorer performance. Applying the method to the tower-based brightness temperature achieved an average estimation depth of 28cm (20cm) and 5cm (5cm) during the dry and wet periods respectively when using a second-order polynomial (linear) function. These findings highlight the opportunity of a satellite mission with L-band and P-band observations to accurately estimate the soil moisture profile to as deep as 30cm globally.

Foad Brakhasi↗

Object-Based Comparison of Data-Driven and Physics-Driven Satellite Estimates of Extreme Rainfall

The Global Precipitation Measurement (GPM) constellation of spaceborne sensors provides a variety of direct and indirect measurements of precipitation processes. Such observations can be employed to derive spatially and temporally consistent gridded precipitation estimates either via data-driven retrieval algorithms or by assimilation into physically based numerical weather models. We compare the data-driven Integrated Multisatellite Retrievals for GPM (IMERG) and the assimilation-enabled NASA-Unified Weather Research and Forecasting (NU-WRF) model against Stage IV reference precipitation for four major extreme rainfall events in the southeastern United States using an object-based analysis framework that decomposes gridded precipitation fields into storm objects. As an alternative to conventional ‘‘grid-by-grid analysis,’’ the object-based approach provides a promising way to diagnose spatial properties of storms, trace them through space and time, and connect their accuracy to storm types and input data sources. The evolution of two tropical cyclones are generally captured by IMERG and NU-WRF, while the less organized spatial patterns of two mesoscale convective systems pose challenges for both. NU-WRF rain rates are generally more accurate, while IMERG better captures storm location and shape. Both show higher skill in detecting large, intense storms compared to smaller, weaker storms. IMERG’s accuracy depends on the input microwave and infrared data sources; NU-WRF does not appear to exhibit this dependence. Findings highlight that an object-oriented view can provide deeper insights into satellite precipitation performance and that the satellite precipitation community should further explore the potential for ‘‘hybrid’’ data-driven and physics-driven estimates in order to make optimal usage of satellite observations.

extreme events↗

A Hyperspectral Inversion Framework for Estimating Absorbing Inherent Optical Properties and Biogeochemical Parameters in Inland and Coastal Waters

The simultaneous remote estimation of biogeochemical parameters (BPs) and inherent optical properties (IOPs) from hyperspectral satellite imagery of globally distributed optically distinct inland and coastal waters is a complex, unsolved, non-unique inverse problem. To tackle this problem, we leverage a machine-learning model termed Mixture Density Networks (MDNs). MDNs outperform operational algorithms by calculating the covariance between the simultaneously estimated products. We train the MDNs on a large ( N = 8237) dataset of co-aligned, in situ measured, hyperspectral remote sensing reflectance (R rs ), BPs, and absorbing IOPs from globally representative optically distinct inland and coastal waters. The estimated IOPs include absorption due to phytoplankton (a ph ), chromophoric dissolved organic matter (a cdom ), and non-algal particles (a nap ). The estimated BPs include chlorophyll-a, total suspended solids, and phycocyanin (PC). MDNs dramatically reduce uncertainty in the retrievals, relative to operational algorithms, when using a 50/50 dataset split, where the MDNs are trained on a randomly selected half of the in situ dataset and validated on the other half. Our model is shown to have higher, or equivalent, generalization performance than the calculated operational algorithms available for all BPs and IOPs (except PC) via a leave-one-out cross-validation assessment. The MDNs are sensitive to uncertainties in the hyperspectral satellite R rs , resulting from instrument noise and atmospheric correction; there is a difference of ~37.4–62.8% (using median symmetric accuracy) between the MDNs’ estimates derived from co-located satellite-derived R rs and in situ R rs . Of the IOPs, a cdom and a nap are less sensitive to uncertainties in hyperspectral satellite imagery relative to a ph , with remote estimates of a ph exhibiting incorrect spectral shape and magnitude relative to in situ measured IOPs. Despite the uncertainties in satellite derived R rs , the spatial distributions of BPs and IOPs in MDN-derived product maps of Lake Erie and the Curonian Lagoon, based on imagery taken with the Hyperspectral Imager for the Coastal Ocean (HICO) and PRecursore Iper-Spettrale della Missione Applicativa (PRISMA), are confirmed via co-aligned in situ measurements and agree with the literature’s understanding of these well-studied regions. The consistency and accuracy of the model on HICO and PRISMA imagery, despite radiometric uncertainties, demonstrate its applicability to future hyperspectral missions, such as the Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission, where the simultaneous estimation model will serve as a key part of phytoplankton community composition analysis.

Ryan E. O'Shea↗

Fomalhaut's Debris Disk and Planet: Constraining the Mass of Formalhaut B from Disk Morphology

Following the optical imaging of exoplanet candidate Fomalhaut b (Fom b), we present a numerical model of how Fomalhaut's debris disk is gravitationally shaped by a single interior planet. The model is simple, adaptable to other debris disks, and can be extended to accommodate multiple planets. If Fom b is the dominant perturber of the belt, then to produce the observed disk morphology it must have a mass M(sub pl) < 3M(sub J), an orbital semimajor axis a(sub pl) > 101.5AU, and an orbital eccentricity e(sub pl) = 0.11 - 0.13. These conclusions are independent of Fom b's photometry. To not disrupt the disk, a greater mass for Fom b demands a smaller orbit farther removed from the disk; thus, future astrometric measurement of Fom b's orbit, combined with our model of planet-disk interaction, can be used to determine the mass more precisely. The inner edge of the debris disk at a approximately equals 133AU lies at the periphery of Fom b's chaotic zone, and the mean disk eccentricity of e approximately equals 0.11 is secularly forced by the planet, supporting predictions made prior to the discovery of Fom b. However, previous mass constraints based on disk morphology rely on several oversimplifications. We explain why our constraint is more reliable. It is based on a global model of the disk that is not restricted to the planet's chaotic zone boundary. Moreover, we screen disk parent bodies for dynamical stability over the system age of approximately 100 Myr, and model them separately from their dust grain progeny; the latter's orbits are strongly affected by radiation pressure and their lifetimes are limited to approximately 0.1 Myr by destructive grain-grain collisions. The single planet model predicts that planet and disk orbits be apsidally aligned. Fomalhaut b's nominal space velocity does not bear this out, but the astrometric uncertainties are difficult to quantify. Even if the apsidal misalignment proves real, our calculated upper mass limit of 3 M(sub J) still holds. Parent bodies are evacuated from mean-motion resonances with Fom b; these empty resonances are akin to the Kirkwood gaps opened by Jupiter. The belt contains at least 3M(sub Earth) of solids that are grinding down to dust, their velocity dispersions stirred so strongly by Fom b that collisions are destructive. Such a large mass in solids is consistent with Fom b having formed in situ.

Chiang, E.↗

Airborne Lidar Measurements of Atmospheric Pressure Made Using the Oxygen A-Band

Accurate measurements of greenhouse gas mixing ratios on a global scale are currently needed to gain a better understanding of climate change and its possible impact on our planet. In order to remotely measure greenhouse gas concentrations in the atmosphere with regard to dry air, the air number density in the atmosphere is also needed in deriving the greenhouse gas concentrations. Since oxygen is stable and uniformly mixed in the atmosphere at 20.95%, the measurement of an oxygen absorption in the atmosphere can be used to infer the dry air density and used to calculate the dry air mixing ratio of a greenhouse gas, such as carbon dioxide or methane. OUT technique of measuring Oxygen uses integrated path differential absorption (IPDA) with an Erbium Doped Fiber Amplifier (EDF A) laser system and single photon counting module (SPCM). It measures the absorbance of several on- and off-line wavelengths tuned to an O2 absorption line in the A-band at 764.7 nm. The choice of wavelengths allows us to maximize the pressure sensitivity using the trough between two absorptions in the Oxygen A-band. Our retrieval algorithm uses ancillary meteorological and aircraft altitude information to fit the experimentally obtained lidar O2 line shapes to a model atmosphere and derives the pressure from the profiles of the two lines. We have demonstrated O2 measurements from the ground and from an airborne platform. In this paper we will report on our airborne measurements during our 2011 campaign for the ASCENDS program.

Riris, Haris↗

Informing Algorithms from Ground Validation: The GPM Combined Algorithm

Ground validation (GV) for the Global Precipitation Measurement (GPM) mission encompasses in-situ (e.g., gauge, disdrometer) measurements, ground radar products, and comprehensive datasets from dedicated airborne field campaigns. These datasets are used for direct validation of the precipitation products from GPM as well as to inform assumptions used by the algorithms that produce these products. This presentation will focus on current and potential such uses of GV data in the GPM combined radar-radiometer algorithm.The GPM combined algorithm, by virtue of using data from the Dual-frequency Precipitation Radar (DPR) and GPM Microwave Imager (GMI), is the most well-constrained instantaneous precipitation product from GPM. It also plays an important in role the passive microwave algorithms as a basis for the construction of brightness temperature-precipitation profile databases and is a calibrator for the multi- satellite gridded product IMERG. Since both radar and radiometer data are used as observational inputs, and even with these data the retrievals are underconstrained, microphysical properties of the hydrometeor profile that are relevant over the range of GMI and DPR wavelengths and incidence angles need to be assumed by the algorithm forward models. Such properties include the normalized intercept parameter (N w ) and shape parameter () of the particle size distribution, ice size-density relationship, and particle size-aspect ratio relationships. The sub-beam variability must also be prescribed in order to accurately simulate the observed radar reflectivity profiles without introducing significant biases. The sensitivity of the algorithm to these parameters, along with their mean values and variability, will be discussed. Examples of future directions include refinement with new data (particularly the assumptions related to ice), and seeking relationships between assumptions and observable or environmental data.

Munchak, Joe.↗

Transient Creep of a Composite Lower Crust: A Polymineralic Basis for Rapidly Evolving Postseismic Deformation Modes - 2

Postseismic horizontal strain and displacement following the June 28, 1992, Landers, California, earthquake (M(sub W) 7.3) is broad scale and cannot be explained solely by delayed afterslip located at the rupturing fault trace. Both the observed strain at Pifion Flat Observatory (PFO) and observed Global Positioning System receiver velocities evolve rapidly after the Landers-Big Bear earthquake sequence. The observed exponential decay of these motions, with timescales of 4-34 days, may reflect a soft creep rheology in the lower crust and brittle-ductile transition zone or even within the seismogenic crust itself. Here a simple model of a two-dimensional screw dislocation in a layered Maxwell viscoelastic Earth is employed in conjunction with a composite rheology to demonstrate that the short timescale transient response modes (approx. = 4-34 days) are consistent with the behavior of a biviscous lower crust. The lowest viscosity of this system is derivable from laboratory experimental data on the long-term creep of natural quarztites, and the highest viscosity is consistent with isostasy-related lower crustal flow in a continental extensional tectonic environment. The model predicts significant stress relaxation at the base of the seismogenic crust. Near the base of the seismogenic zone, and about 4 km away from the mainshock, the rate of predicted relaxation is of the order of 0.01 MPa/ d during the first 20 days of postseismic flow. Oblate spheroidal inclusions at 5% concentration levels that are both aligned and fairly flat in shape and that have a viscosity of 3-4 x 10(exp 15) Pa s are consistent with both the amplitude and decay time of horizontal crustal strain observed at PFO after the Landers mainshock. It is speculated that the structures exposed in cross sections and in seismic reflection profiles of the lower crust that have mylonitic associations are, in part, the cause of such rapid postseismic evolution in southeastern California. Unmylonitized quartz-rich rock at sufficiently elevated temperatures could also contribute to the rapid decay modes.

Ivins, Erik R.↗

Geochemical Mapping of 4 Vesta Begins

By December, the NASA Dawn spacecraft will have descended to a low altitude mapping orbit (LAMO), where the Gamma Ray and Neutron Detector (GRaND) will acquire global mapping data for up to four months. Measurements by GRaND will help answer elusive questions about how Vesta differentiated and the nature of processes that shaped Vesta s surface. The data will be analyzed to determine the abundances of Mg, Si, Fe, K, Th, and H at a spatial resolution of roughly 300 km full-width-at-half-maximum from a 465 km radius orbit. Thermal and fast neutron counting data will be analyzed to determine the neutron macroscopic absorption cross section and average atomic mass, providing constraints on additional elements, such as Ca and Al. GRaND will quantify the elemental composition of coarse spatial units identified by Dawn s Framing Camera (FC) and the Visible & Infrared Spectrometer (VIR). In addition, GRaND will map the mixing ratio of compositional end members selected from the howardite, eucrite and diogenite (HED) meteorites, determine the relative proportions of plagioclase and mafic minerals, and search for compositions that are absent or under-represented in the meteorite collection. While it is generally thought that Vesta s crust on a regional scale should be well-represented by linear mixing of HED whole-rock compositions, there are hints that Vesta may be more diverse than implied by this model. For example, the discovery of K-rich impact glasses in howardites suggests that K-rich rocks may be present on a portion of Vesta s surface, and the analysis of diogenites indicates considerable variability in the magmatic processes that formed them. The chemical composition of materials within Vesta s south polar structure may provide further clues to how it formed. An impact might have exposed mantle and lower crustal materials, which should have a distinctive compositional signature. We present the analysis of data acquired by GRaND from cruise through the descent to LAMO, including GRaND s sensitivity to different elements and geochemical processes.

Prettyman, Thomas H.↗

A Curved, Elastostatic Boundary Element for Plane Anisotropic Structures

The plane-stress equations of linear elasticity are used in conjunction with those of the boundary element method to develop a novel curved, quadratic boundary element applicable to structures composed of anisotropic materials in a state of plane stress or plane strain. The curved boundary element is developed to solve two-dimensional, elastostatic problems of arbitrary shape, connectivity, and material type. As a result of the anisotropy, complex variables are employed in the fundamental solution derivations for a concentrated unit-magnitude force in an infinite elastic anisotropic medium. Once known, the fundamental solutions are evaluated numerically by using the known displacement and traction boundary values in an integral formulation with Gaussian quadrature. All the integral equations of the boundary element method are evaluated using one of two methods: either regular Gaussian quadrature or a combination of regular and logarithmic Gaussian quadrature. The regular Gaussian quadrature is used to evaluate most of the integrals along the boundary, and the combined scheme is employed for integrals that are singular. Individual element contributions are assembled into the global matrices of the standard boundary element method, manipulated to form a system of linear equations, and the resulting system is solved. The interior displacements and stresses are found through a separate set of auxiliary equations that are derived using an Airy-type stress function in terms of complex variables. The capabilities and accuracy of this method are demonstrated for a laminated-composite plate with a central, elliptical cutout that is subjected to uniform tension along one of the straight edges of the plate. Comparison of the boundary element results for this problem with corresponding results from an analytical model show a difference of less than 1%.

Smeltzer, Stanley S.↗

X-Ray Spectra from MHD Simulations of Accreting Black Holes

We present the results of a new global radiation transport code coupled to a general relativistic magneto-hydrodynamic simulation of an accreting, nonrotating black hole. For the first time, we are able to explain from first principles in a self-consistent way the X-ray spectra observed from stellar-mass black holes, including a thermal peak, Compton reflection hump, power-law tail, and broad iron line. Varying only the mass accretion rate, we are able to reproduce the low/hard, steep power-law, and thermal-dominant states seen in most galactic black hole sources. The temperature in the corona is T(sub e) ~ 10 keV in a boundary layer near the disk and rises smoothly to T(sub e) greater than or approximately 100 keV in low-density regions far above the disk. Even as the disk's reflection edge varies from the horizon out to approximately equal to 6M as the accretion rate decreases, we find that the shape of the Fe Ka line is remarkably constant. This is because photons emitted from the plunging region are strongly beamed into the horizon and never reach the observer. We have also carried out a basic timing analysis of the spectra and find that the fractional variability increases with photon energy and viewer inclination angle, consistent with the coronal hot spot model for X-ray fluctuations.

Schnittman, Jeremy D.↗

Trajectory Engineering with Modular Patched Conics for Entry Systems and TPS (TEMPEST)

Brief Presenter Biography (35 word limit): Bohdan Wesely is an Aerospace Engineer in the Entry Systems and Technology Division at Ames. He has worked on a variety of projects for NASA including integrated TPS (thermal protection system) flight hardware deliveries and testing services for commercial partners. Introduction: TEMPEST is a new trajectory analysis framework that is designed to fill the gap between dedicated flight mechanics tools and aerothermal and TPS sizing tools. The project started as an SJSU master’s thesis and has since evolved into a general conceptual design tool capable of studying a wide variety of entry problems. Development is ongoing in the Entry Systems and Technology Division at NASA ARC. Why TEMPEST: Space missions involving entry into a planetary atmosphere involve a series of unique requirements across multiple disciplines. Whether it is traditional entry descent and landing (EDL), or aerocapture, the vehicle must navigate to its target landing location or orbit state, and the TPS must protect the payload during entry. The design process typically involves iterative handoffs between various flight mechanics, flow solver, and material response level tools. During the early conceptual phase, a wide variety of feasible trajectories are simulated in a Monte Carlo scenario which broadly satisfy the mission or landing requirements. Next, computational fluid dynamics (CFD), direct simulation Monte Carlo (DSMC), and other flow solver analyses are performed at various key trajectory points to generate an aero-database, heating and TPS design requirements also emerge at this stage. At this point, with updated aerodynamics from the various flow solvers, trajectories can be re-run, this in turn can change the required freestream conditions for the CFD tools, and as a project progresses, these analyses converge, and uncertainty is reduced. However, there is always a “hand-off” occurring between two inherently coupled phenomena. Analysis Description: One of the goals with TEMPEST is to use a variety of first principles estimation methods coupled with an atmosphere model to predict vehicle aerothermodynamics across the entire flight regime while propagating a 3 or 6 degree of freedom (DoF) trajectory. Aerodynamics methods include modified Newtonian, Maxwell and Cercignani- Lampis-Lord (CLL) for continuum, transitional, and free molecular flow regimes. Aerothermodynamics include boundary layer and reference enthalpy methods, and Mutation++ for non-equilibrium chemistry modeling. TEMPEST is also capable of stitching multiple trajectory segments together to study mission scenarios like multi-pass aerocapture and aero-gravity assists. Most of the program is implemented in MATLAB using modern system objects, it relies on several C++ shared libraries for supporting tools like Gmsh, the Global Reference Atmospheric Model (GRAM), and Mutation++. The various first principles aerothermal estimation methods are discretized across either a structured axisymmetric panel mesh or an unstructured tri-mesh generated from an open-source tool such as Gmsh, this allows solutions on the same mesh to be compared across tools such as CB-Aero. CFD Coupling. A physics-aware, gaussian process CFD anchoring scheme is proposed to adjust the various first principles methods as a CFD database is populated. One goal for this anchoring module is to inform the project where CFD should be run. Full knowledge of the entire trajectory, atmosphere, and aerothermodynamics allows for easier identification of high sensitivity areas and uncertainty quantification. While the first principles effects are well known and proven accurate in existing tools such as CB- Aero and Cart3D, a physics aware CFD anchoring scheme increases tool credibility across a project lifecycle. Material Response Modeling. Correct TPS sizing is critical for optimizing mass for science payloads and ensuring mission success. The process typically involves a thermal analysis along the trajectory with surface heating environments as a boundary condition. Several design constraints are maximum bondline temperature and maximum recession with various margining techniques. The material response tool FIAT, developed out of NASA Ames, is currently being integrated into the TEMPEST environment. TPS recession, shape change, mass loss, and mass property alteration are all factors that can perturb an entry trajectory. For missions like Mars 2020, recession was minimal and was safely handled separately as a post process. For missions such as Jupiter Galileo with a high TPS mass fraction or asteroid entries, recession plays a major role. The proposed fully coupled scheme is to use an epoch-based approach where the trajectory integration is halted after a recession threshold, the energy balance and FIAT are solved at each panel, the mesh, aerodynamics, and mass properties are updated, and the trajectory continues. Several computational tradeoffs have been made during the development of TEMPEST to limit the cost of a single trajectory and preserve its utility as a conceptual, rapid iteration tool. Conclusion: Development of TEMPEST is ongoing and the project is still in its infancy. This talk aims to showcase its unique capabilities to support future NASA entry systems missions.

Bohdan O Wesely↗

Exploring the Utility of the Planned CYGNSS Mission for Investigating the Initiation and Development of the Madden-Julian Oscillation

CYGNSS is a planned constellation consisting of multiple micro-satellites that leverage the Global Positioning System (GPS) to provide rapidly updated, high resolution (approx. 15-50 km, approx. 4 h) surface wind speeds (via bi-static scatterometry) over the tropical oceans in any weather condition, including heavy rainfall. The approach of the work to be presented at this conference is to utilize a limited-domain, cloud-system resolving model (Weather Research and Forecasting or WRF) and its attendant data assimilation scheme (Three-Dimensional Variational Assimilation or 3DVAR) to investigate the utility of the CYGNSS mission for helping characterize key convectiveto- mesoscale processes - such as surface evaporation, moisture advection and convergence, and upscale development of precipitation systems - that help drive the initiation and development of the Madden-Julian Oscillation (MJO) in the equatorial Indian Ocean. The proposed work will focus on three scientific objectives. Objective 1 is to produce a high-resolution surface wind dataset resolution (approx. 0.5 h, approx. 1-4 km) for multiple MJO onsets using WRF-assimilated winds and other data from the DYNAmics of the MJO (DYNAMO) field campaign, which took place during October 2011 - March 2012. Objective 2 is to study the variability of surface winds during MJO onsets at temporal and spatial scales of finer resolution than future CYGNSS data. The goal is to understand how sub-CYGNSS-resolution processes will shape the observations made by the satellite constellation. Objective 3 is to ingest simulated CYGNSS data into the WRF model in order to perform observing system simulation experiments (OSSEs). These will be used to test and quantify the potential beneficial effects provided by CYGNSS, particularly for characterizing the physical processes driving convective organization and upscale development during the initiation and development of the MJO. The proposed research is ideal for answering important questions about the CYGNSS mission, such as the representativeness of surface wind retrievals in the context of the complex airflow processes that occur during heavy precipitation, as well as the tradeoffs in retrieval accuracy that result from finer spatial resolution of the CYGNSS winds versus increased errors/noisiness in those data. Research plans and initial progress toward these objectives will be presented.

Lang, Timothy↗