SEARCH · Engineering Papers

Results for “generalized linear model”

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 1,027 records · Page 57

All sources NASA NTRS DOE OSTI

NASA NTRS2014

Applying Motion Constraints Based on Test Data

MSC ADAMS is a simulation software that is used to analyze multibody dynamics. Using user subroutines, it is possible to apply motion constraints to the rigid bodies so that they match the motion profile collected from test data. This presentation describes the process of taking test data and passing it to ADAMS using user subroutines, and uses the Morpheus free-flight 4 test as an example of motion data used for this purpose. Morpheus is the name of a prototype lander vehicle built by NASA that serves as a test bed for various experimental technologies (see backup slides for details) MSC.ADAMS"TM" is used to play back telemetry data (vehicle orientation and position) from each test as the inputs to a 6-DoF general motion constraint (details in backup slides) The MSC.ADAMS"TM" playback simulations allow engineers to examine and analyze flight trajectory as well as observe vehicle motion from any angle and at any playback speed. This facilitates the development of robust and stable control algorithms, increasing reliability and reducing development costs of this developmental engine The simulation also incorporates a 3D model of the artificial hazard field, allowing engineers to visualize and measure performance of the developmental autonomous landing and hazard avoidance technology ADAMS is a multibody dynamics solver. It uses forces, constraints, and mass properties to numerically integrate equations of motion. The ADAMS solver will ask the motion subroutine for position, velocity, and acceleration values at various time steps. Those values must be continuous over the whole time domain. Each degree of freedom in the telemetry data can be examined separately; however, linear interpolation of the telemetry data is invalid, since there will be discontinuities in velocity and acceleration.

Burlone, Michael↗

NASA NTRS2014

Simulation of VSPT Experimental Cascade Under High and Low Free-Stream Turbulence Conditions

Variable-Speed Power Turbines (VSPT) for rotorcraft applications operate at low Reynolds number and over a wide range in incidence associated with shaft speed change. A comprehensive linear cascade data set obtained includes the effects of Reynolds number, free-stream turbulence and incidence is available and this paper concerns itself with the presentation and numerical simulation of conditions resulting in a selected set of those data. As such, post-dictions of blade pressure loading, total-pressure loss and exit flow angles under conditions of high and low turbulence intensity for a single Reynolds number are presented. Analyses are performed with the three-equation turbulence models of Walters-Leylek and Walters and Cokljat. Transition, loading, total-pressure loss and exit angle variations are presented and comparisons are made with experimental data as available. It is concluded that at the low freestream turbulence conditions the Walters-Cokljat model is better suited to predictions while for high freestream conditions the two models generate similar predications that are generally satisfactory.

low pressure trubine↗

NASA NTRS2013

Spring-Summer Temperatures Since AD 1780 Reconstructed from Stable Oxygen Isotope Ratios in White Spruce Tree-Rings from the Mackenzie Delta, Northwestern Canada

High-latitude delta(exp 18)O archives deriving from meteoric water (e.g., tree-rings and ice-cores) can provide valuable information on past temperature variability, but stationarity of temperature signals in these archives depends on the stability of moisture source/trajectory and precipitation seasonality, both of which can be affected by atmospheric circulation changes. A tree-ring delta(exp 18)O record (AD 1780-2003) from the Mackenzie Delta is evaluated as a temperature proxy based on linear regression diagnostics. The primary source of moisture for this region is the North Pacific and, thus, North Pacific atmospheric circulation variability could potentially affect the tree-ring delta(exp 18)O-temperature signal. Over the instrumental period (AD 1892-2003), tree-ring delta(exp 18)O explained 29% of interannual variability in April-July minimum temperatures, and the explained variability increases substantially at lower-frequencies. A split-period calibration/verification analysis found the delta(exp 18)O-temperature relation was time-stable, which supported a temperature reconstruction back to AD 1780. The stability of the delta(exp 18)O-temperature signal indirectly implies the study region is insensitive to North Pacific circulation effects, since North Pacific circulation was not constant over the calibration period. Simulations from the NASA-GISS ModelE isotope-enabled general circulation model confirm that meteoric delta(exp 18)O and precipitation seasonality in the study region are likely insensitive to North Pacific circulation effects, highlighting the paleoclimatic value of tree-ring and possibly other delta(exp 18)O records from this region. Our delta(exp 18)O-based temperature reconstruction is the first of its kind in northwestern North America, and one of few worldwide, and provides a long-term context for evaluating recent climate warming in the Mackenzie Delta region.

Canada↗

NASA NTRS1998

Making Intercontinental Burning Influences Quantitative: Lessons on Transport and Geochemistry From the 3-D Simulation of CO in PEM-T A and TRACE-A

We report on two aspects of the simulation of global transport of plumes originating from subtropical biomass fires. We use of meteorological assimilation (MM5) at 2-degree resolution with a Grail cloud parameterization and a Blackadar-based planetary boundary layer parameterization. Ames's GRACES model provides emissions, transport, and an appropriate level of simulated chemical transformation. We have worked with passive-tracer CO or linear chemistry. This is appropriate since we find major work to be done in evaluating CO source strengths and transport mechanisms before chemical integrations could be meaningful. First, we present mechanisms by which CO and other pollutants are introduced into the free troposphere, and are then transported with hate dilution from approx. 0 to approx. 180 degrees longitude. One principal conduit for these plumes is the vernal subtropical jet; however the plumes appear at various altitudes and latitudes as they influenced by frontal motions and (most likely) radiative processes. A common, repeated pattern of transport has pollutant plumes arriving in the distant Pacific Ocean from Africa and South America at 25 degrees south and 14 km altitude. Following this, there is then a general appearance of pollution at extending down to 5 kin at more equatorial (10 S) and polar latitudes (to 45 S). Second, we evaluate the quantitative success of our simulation. (Such success requires efforts considerably beyond trajectory analyses, and is necessary for our community to claim an understanding of the effects of biomass burning on global atmospheric chemistry and the planet's trend in oxidizing capacity.) We find that we simulate most pollution episodes sampled by Glenn Sachse's CO instrument and the Blake hydrocarbon analyses during PENT A. We will present our current ideas on why our general levels appear satisfactory when the observations are within 20 ppb of background levels, but substantially miss the variability associated with the most intense biomass burning plumes. As time allows, we will speculate on the origins of a class of plumes that appear in the mid-troposphere, but are missing from our simulations.

Chatfield, Robert↗

NASA NTRS2010

Efficient Web Services Policy Combination

Large-scale Web security systems usually involve cooperation between domains with non-identical policies. The network management and Web communication software used by the different organizations presents a stumbling block. Many of the tools used by the various divisions do not have the ability to communicate network management data with each other. At best, this means that manual human intervention into the communication protocols used at various network routers and endpoints is required. Developing practical, sound, and automated ways to compose policies to bridge these differences is a long-standing problem. One of the key subtleties is the need to deal with inconsistencies and defaults where one organization proposes a rule on a particular feature, and another has a different rule or expresses no rule. A general approach is to assign priorities to rules and observe the rules with the highest priorities when there are conflicts. The present methods have inherent inefficiency, which heavily restrict their practical applications. A new, efficient algorithm combines policies utilized for Web services. The method is based on an algorithm that allows an automatic and scalable composition of security policies between multiple organizations. It is based on defeasible policy composition, a promising approach for finding conflicts and resolving priorities between rules. In the general case, policy negotiation is an intractable problem. A promising method, suggested in the literature, is when policies are represented in defeasible logic, and composition is based on rules for non-monotonic inference. In this system, policy writers construct metapolicies describing both the policy that they wish to enforce and annotations describing their composition preferences. These annotations can indicate whether certain policy assertions are required by the policy writer or, if not, under what circumstances the policy writer is willing to compromise and allow other assertions to take precedence. Meta-policies are specified in defeasible logic, a computationally efficient non-monotonic logic developed to model human reasoning. One drawback of this method is that at one point the algorithm starts an exhaustive search of all subsets of the set of conclusions of a defeasible theory. Although the propositional defeasible logic has linear complexity, the set of conclusions here may be large, especially in real-life practical cases. This phenomenon leads to an inefficient exponential explosion of complexity. The current process of getting a Web security policy from combination of two meta-policies consists of two steps. The first is generating a new meta-policy that is a composition of the input meta-policies, and the second is mapping the meta-policy onto a security policy. The new algorithm avoids the exhaustive search in the current algorithm, and provides a security policy that matches all requirements of the involved metapolicies.

Vatan, Farrokh↗

NASA NTRS2013

NIR Imaging Spectroscopy of the Inner Few Arcseconds of NGC 4151 with OSIRIS at Keck

We present H- and K-band data from the inner arcsecond of the Seyfert 1.5 galaxy NGC 4151 obtained with the adaptive optics assisted near-infrared imaging field spectrograph OSIRIS at the Keck Observatory. The angular resolution is about a few parsecs on-site and thus competes easily with optical images taken previously with the Hubble Space Telescope. We present the morphology and dynamics of most species detected but focus on the morphology and dynamics of the narrow line region (as traced by emission of [FeII]λ1.644 μm), the interplay between plasma ejected from the nucleus (as traced by 21 cm continuum radio data) and hot H2 gas and characterize the detected nuclear HeIλ2.058 μm absorption feature as a narrow absorption line (NAL) phenomenon. Emission from the narrow line region (NLR) as traced by [FeII] reveals a biconical morphology and we compare the measured dynamics in the [FeII] emission line with models proposing acceleration of gas in the NLR and simple ejection of gas into the NLR. In the inner 2.5 arcseconds the acceleration model reveals a better fit to our data than the ejection model.We also see evidence that the jet very locally enhances emission in [FeII] at certain positions in our field-of-view such that we were able to distinct the kinematics of these clouds from clouds generally accelerated in the NLR. Further, the radio jet is aligned with the bicone surface rather than the bicone axis such that we assume that the jet is not the dominant mechanism responsible for driving the kinematics of clouds in the NLR. The hot H2 gas is thermal with a temperature of about 1700 K. We observe a remarkable correlation between individual H2 clouds at systemic velocity with the 21 cm continuum radio jet. We propose that the radio jet is at least partially embedded in the galactic disk of NGC 4151 such that deviations from a linear radio structure are invoked by interactions of jet plasma with H2 clouds that are moving into the path of the jet because of rotation of the galactic disk of NGC 4151. Additionally, we observe a correlation of the jet as traced by the radio data, with gas as traced in Brγ and H2, at velocities between systemic and +/- 200 km/s at several locations along the path of the jet. The HeIλ2.058 μm line in NGC 4151 appears in emission with a blueshifted absorption component from an outflow. The emission (absorption) component has a velocity offset of 10 km/s (-280 km/s) with a Gaussian (Lorentzian) full-width (half-width) at half maximum of 160 km/s (440 km/s). The absorption component remains spatially unresolved and its kinematic measures differ from that of UV resonance absorption lines. From the amount of absorption we derive a lower limit of the HeI 2S column density of 1 × 10(exp 14) cm−2 with a covering factor along the line-of-sight of C(sub los) approximately equal to 0.1.

K-band data↗

NASA NTRS2022

SatNet: A Benchmark for Satellite Scheduling Optimization

Satellites provide essential services such as networking and weather tracking, and the number of near-earth and deep space satellites are expected to grow rapidly in the coming years. Communications with terrestrial ground stations is one of the critical functionalities of any space mission. Satellite scheduling is a problem that has been scientifically investigated since the 1970s. A central aspect of this problem is the need to consider resource contention and satellite visibility constraints as they require line of sight. Due to the combinatorial nature of the problem, prior solutions such as linear programs and evolutionary algorithms require extensive compute capabilities to output a feasible schedule for each scenario. Machine learning based scheduling can provide an alternative solution by training a model with historical data and generating a schedule quickly with model inference. We present SatNet, a benchmark for satellite scheduling optimization based on historical data from the NASA Deep Space Network. We propose formulation of the satellite scheduling problem as a Markov Decision Process and use reinforcement learning (RL) policies to generate schedules. The nature of constraints imposed by SatNet differ from other combinatorial optimization problems such as vehicle routing studied in prior literature. Our initial results indicate that RL is an alternative optimization approach that can generate candidate solutions of comparable quality to existing state-of-the-practice results. However, we also find that RL policies overfit to the training dataset and do not generalize well to new data, thereby necessitating continued research on reusable and generalizable agents.

Wilson, Brian↗

NASA NTRS1986

The single scattering phase functions of Jupiter's clouds

The determination of the single scattering phase functions of Jupiter's clouds and a thin upper haze by Tomasko et al. was refined and extended to seven latitudes in blue and red light. The phase function is well-constrained by the Pioneer 10 and 11 photometric data sets. Multiple scattering models were computed to match the limb darkening at each latitude at up to 15 phase angles from 12 deg to 151 deg. Ground-based observations were used for absolute calibration and to extend the data to lower phase angles. The phase functions were parameterized using the double Henyey-Greenstein function. The three Henyey-Greenstein parameters and the single scattering albedo were determined using a non-linear least squares method for the haze and the clouds below. The phase functions derived for the northen zone and belt are remarkably similar to the phase functions of the corresponding regions in the south, with most of the differences in brightness of the northern and southern features resulting from minor differences in single scattering albedo. Analysis of the Equatorial Region is complicated by the presence of numerous small features, but the phase function required is generally similar to that seen in the more homogeneous regions. Details of the phase functions of the haze and clouds are presented, and the differences between the cloud phase functions at low and high latitudes in red and blue light are discussed.

Doose, L. R.↗

NASA NTRS2014

A Technique for Mapping Characteristic Lengths to Preserve Energy Dissipated via Strain Softening in a Multiscale Analysis

It is often advantageous to account for the microstructure of the material directly using multiscale modeling. For computational tractability, an idealized repeating unit cell (RUC) is used to capture all of the pertinent features of the microstructure. Typically, the RUC is dimensionless and depends only on the relative volume fractions of the different phases in the material. This works well for non-linear and inelastic behavior exhibiting a positive-definite constitutive response. Although, once the material exhibits strain softening, or localization, a mesh objective failure theories, such as smeared fracture theories, nodal and element enrichment theories (XFEM), cohesive elements or virtual crack closure technique (VCCT), can be utilized at the microscale, but the dimensions of the RUC must then be defined. One major challenge in multiscale progressive damage modeling is relating the characteristic lengths across the scales in order to preserve the energy that is dissipated via localization at the microscale. If there is no effort to relate the size of the macroscale element to the microscale RUC, then the energy that is dissipated will remain mesh dependent at the macroscale, even if it is regularized at the microscale. Here, a technique for mapping characteristic lengths across the scales is proposed. The RUC will be modeled using the generalized method of cells (GMC) micromechanics theory, and local failure in the matrix constituent subcells will be modeled using the crack band theory. The subcell characteristic lengths used in the crack band calculations will be mapped to the macroscale finite element in order to regularize the local energy in a manner consistent with the global length scale. Examples will be provided with and without the regularization, and they will be compared to a baseline case where the size and shape of the element and RUC are coincident (ensuring energy is preserved across the scales).

softening↗

NASA NTRS1990

Fully elliptic incompressible flow calculations on regular grid by a new pressure substitution method

A new method is presented for the solution of incompressible flow in generalized coordinates. This method is based on the substitution of the pressure weighted form of the momentum equations into the continuity equation. The algorithm is rigorously derived and a Fourier analysis is used to assess its suitability to act as an error smoother. Linear stability analysis results indicate that the performance of the new pressure substitution method (PSM) and the pressure correction method (PCM) is about the same at low Reynolds numbers, with no significant pressure gradient. At high Reynolds numbers the PSM shows much faster convergence. Likewise prediction of various flows indicate that the PSM has better accuracy for high Reynolds number flows with significant pressure gradients. Since most practical aerodynamic flows have significant pressure gradients, the PSM seems to be attractive for such flows. Solutions for both laminar and turbulent flow are compared with the experimental data. A two-equation low Reynolds number turbulence model is used to resolve the turbulent flowfield.

Hobson, G. V.↗

NASA NTRS2001

Extending unbiased stereology of brain ultrastructure to three-dimensional volumes

OBJECTIVE: Analysis of brain ultrastructure is needed to reveal how neurons communicate with one another via synapses and how disease processes alter this communication. In the past, such analyses have usually been based on single or paired sections obtained by electron microscopy. Reconstruction from multiple serial sections provides a much needed, richer representation of the three-dimensional organization of the brain. This paper introduces a new reconstruction system and new methods for analyzing in three dimensions the location and ultrastructure of neuronal components, such as synapses, which are distributed non-randomly throughout the brain. DESIGN AND MEASUREMENTS: Volumes are reconstructed by defining transformations that align the entire area of adjacent sections. Whole-field alignment requires rotation, translation, skew, scaling, and second-order nonlinear deformations. Such transformations are implemented by a linear combination of bivariate polynomials. Computer software for generating transformations based on user input is described. Stereological techniques for assessing structural distributions in reconstructed volumes are the unbiased bricking, disector, unbiased ratio, and per-length counting techniques. A new general method, the fractional counter, is also described. This unbiased technique relies on the counting of fractions of objects contained in a test volume. A volume of brain tissue from stratum radiatum of hippocampal area CA1 is reconstructed and analyzed for synaptic density to demonstrate and compare the techniques. RESULTS AND CONCLUSIONS: Reconstruction makes practicable volume-oriented analysis of ultrastructure using such techniques as the unbiased bricking and fractional counter methods. These analysis methods are less sensitive to the section-to-section variations in counts and section thickness, factors that contribute to the inaccuracy of other stereological methods. In addition, volume reconstruction facilitates visualization and modeling of structures and analysis of three-dimensional relationships such as synaptic connectivity.

NASA Discipline Neuroscience↗

NASA NTRS2007

Using EIGER for Antenna Design and Analysis

EIGER (Electromagnetic Interactions GenERalized) is a frequency-domain electromagnetics software package that is built upon a flexible framework, designed using object-oriented techniques. The analysis methods used include moment method solutions of integral equations, finite element solutions of partial differential equations, and combinations thereof. The framework design permits new analysis techniques (boundary conditions, Green#s functions, etc.) to be added to the software suite with a sensible effort. The code has been designed to execute (in serial or parallel) on a wide variety of platforms from Intel-based PCs and Unix-based workstations. Recently, new potential integration scheme s that avoid singularity extraction techniques have been added for integral equation analysis. These new integration schemes are required for facilitating the use of higher-order elements and basis functions. Higher-order elements are better able to model geometrical curvature using fewer elements than when using linear elements. Higher-order basis functions are beneficial for simulating structures with rapidly varying fields or currents. Results presented here will demonstrate curren t and future capabilities of EIGER with respect to analysis of installed antenna system performance in support of NASA#s mission of exploration. Examples include antenna coupling within an enclosed environment and antenna analysis on electrically large manned space vehicles.

Champagne, Nathan J.↗

NASA NTRS2010

Notes on SAW Tag Interrogation Techniques

We consider the problem of interrogating a single SAW RFID tag with a known ID and known range in the presence of multiple interfering tags under the following assumptions: (1) The RF propagation environment is well approximated as a simple delay channel with geometric power-decay constant alpha >/= 2. (2) The interfering tag IDs are unknown but well approximated as independent, identically distributed random samples from a probability distribution of tag ID waveforms with known second-order properties, and the tag of interest is drawn independently from the same distribution. (3) The ranges of the interfering tags are unknown but well approximated as independent, identically distributed realizations of a random variable rho with a known probability distribution f(sub rho) , and the tag ranges are independent of the tag ID waveforms. In particular, we model the tag waveforms as random impulse responses from a wide-sense-stationary, uncorrelated-scattering (WSSUS) fading channel with known bandwidth and scattering function. A brief discussion of the properties of such channels and the notation used to describe them in this document is given in the Appendix. Under these assumptions, we derive the expression for the output signal-to-noise ratio (SNR) for an arbitrary combination of transmitted interrogation signal and linear receiver filter. Based on this expression, we derive the optimal interrogator configuration (i.e., transmitted signal/receiver filter combination) in the two extreme noise/interference regimes, i.e., noise-limited and interference-limited, under the additional assumption that the coherence bandwidth of the tags is much smaller than the total tag bandwidth. Finally, we evaluate the performance of both optimal interrogators over a broad range of operating scenarios using both numerical simulation based on the assumed model and Monte Carlo simulation based on a small sample of measured tag waveforms. The performance evaluation results not only provide guidelines for proper interrogator design, but also provide some insight on the validity of the assumed signal model. It should be noted that the assumption that the impulse response of the tag of interest is known precisely implies that the temperature and range of the tag are also known precisely, which is generally not the case in practice. However, analyzing interrogator performance under this simplifying assumption is much more straightforward and still provides a great deal of insight into the nature of the problem.

Barton, Richard J.↗

NASA NTRS2000

The Use of Ferrofluids to Model Materials Processing (MSFC Center Director's Discretionary Fund)

Many crystals grown in space have structural flaws believed to result from convective motions during the growth phase. The character of these instabilities is not well understood but is associated with thermal and solutal density variations near the solidification interface in the presence of residual gravity and g-jitter. To study these instabilities in a separate, controlled space experiment, a concentration gradient would first have to be artificially established in a timely manner as an initial condition. This is generally difficult to accomplish in a microgravity environment because the momentum of the fluid injected into a test cell tends to swirl around and mix in the absence of a restoring force. The use of magnetic fields to control the motion and position of liquids has received recent, growing interest. The possibility of using the force exerted by a non-uniform magnetic field on a ferrofluid to not only achieve fluid manipulation but also to actively control fluid motion makes it an attractive candidate for space applications. This paper describes a technique for quickly establishing a linear or exponential fluid concentration gradient using a magnetic field in place of gravity to stabilize the deployment. Also discussed is a photometric technique for measuring the concentration profile using light attenuation. Although any range of concentrations can be realized, photometric constraints impose some limitations on measurements. Results of the ground-based experiments indicate that the species distribution is within 3 percent of the predicted value.

Leslie, F.↗

NASA NTRS1994

Physical Retrievals of Over-Ocean Rain Rate from Multichannel Microwave Imagery: Theoretical Characteristics of Normalized Polarization and Scattering Indices - Part 1

Microwave rain rate retrieval algorithms have most often been formulated in terms of the raw brightness temperatures observed by one or more channels of a satellite radiometer. Taken individually, single-channel brightness temperatures generally represent a near-arbitrary combination of positive contributions due to liquid water emission and negative contributions due to scattering by ice and/or visibility of the radiometrically cold ocean surface. Unfortunately, for a given rain rate, emission by liquid water below the freezing level and scattering by ice particles above the freezing level are rather loosely coupled in both a physical and statistical sense. Furthermore, microwave brightness temperatures may vary significantly (approx. 30-70 K) in response to geophysical parameters other than liquid water and precipitation. Because of these complications, physical algorithms which attempt to directly invert observed brightness temperatures have typically relied on the iterative adjustment of detailed micro-physical profiles or cloud models, guided by explicit forward microwave radiative transfer calculations. In support of an effort to develop a significantly simpler and more efficient inversion-type rain rate algorithm, the physical information content of two linear transformations of single-frequency, dual-polarization brightness temperatures is studied: the normalized polarization difference P of Petty and Katsaros (1990, 1992), which is intended as a measure of footprint-averaged rain cloud transmittance for a given frequency; and a scattering index S (similar to the polarization corrected temperature of Spencer et al.,1989) which is sensitive almost exclusively to ice. A reverse Monte Carlo radiative transfer model is used to elucidate the qualitative response of these physically distinct single-frequency indices to idealized 3-dimensional rain clouds and to demonstrate their advantages over raw brightness temperatures both as stand-alone indices of precipitation activity and as primary variables in physical, multichannel rain rate retrieval schemes. As a byproduct of the present analysis, it is shown that conventional plane-parallel analyses of the well-known foot-print-filling problem for emission-based algorithms may in some cases give seriously misleading results.

Petty, G. W.↗

NASA NTRSDate not supplied

An Investigation Into Transecting Satellites in Future Space Traffic Management Scenarios

The number of satellites expected to populate the near-Earth space environment is set to dramatically increase in the coming decade as new large constellations are approved and deployed. Current strategies for deploying new batches of these satellites often involve launching into an initial orbit, and then performing apogee raising maneuvers to reach a target altitude. Similarly, end-of-life planning for these constellation satellites can consist of de-orbit burns that lower perigee to permit disposal via re-entry. Both the raising and de-orbiting maneuvers can result in the individual satellites traveling in transecting orbits that have the potential to cross other spacecraft trajectories. While individual large constellations may be able to coexist in separate altitude and inclination bands, having thousands of satellites moving between these bands as new satellites are added and old satellites are removed could pose additional collision risks. Similar concerns have been raised regarding the impact that large numbers of university-class CubeSats might have in terms of their overall collision risk, especially as these satellites typically do not have propulsion systems for active maneuvering. To assess the impact that transecting satellites might have on future space traffic management strategies, this study explored a variety of realistic future scenarios using a high-fidelity simulation tool. The model can simulate the orbit of tens of thousands of resident space objects (RSOs) simultaneously, to include active satellites, debris, rocket bodies, or even future hypothetical satellite constellations, using a realistic force model that incorporates non-spherical gravity, atmospheric drag, and solar radiation pressure, as well as station-keeping. The simulation can be customized to accommodate different methods of calculating the probability of collision, as well as the process for determining probability ellipsoids and screening volumes. This makes it possible to replicate, and compare, different processes used by different spacecraft operators and space situational awareness (SSA) providers. As the model is run forward in time, each conjunction event is recorded, allowing for the analysis of statistics and meta-data related to these events, providing insight into the nature and frequency of potential collisions, such as whether are they active or passive objects, what size are the two objects, and who owns the objects (if known). This information makes it possible to characterize how changes to the status quo affect the number and type of conjunctions that occur, as well as the distributional effects on various types of satellite operators. To assess the general risks that transecting satellites might pose for hypothetical future space object environments, approximately 60,000 new large constellation satellites were considered, in addition to the existing catalog of approximately 7800 known resident space objects (RSOs), over a simulation period of one year. The results indicate that the future space environment will introduce a non-linear increase in conjunction events as the number of RSOs also increase. This will require adjustments to spacecraft fuel budgets in order to conduct the avoidance maneuvers necessary to minimize collision risk, both for existing and new satellites. This increase is due in large part to the higher density of RSOs and the overlap between some constellation orbits. Current catalog objects were shown to require three times more ∆V for collision avoidance (CA) maneuvers in the simulated future environment, and some constellation spacecraft were estimated to devote the majority of their annual ∆V to CA. The impact of small satellites was found to be proportional for the current space environment, and actually decreased in terms of percentage for the future scenario, suggesting that small satellites do not pose an outsized collision risk. Lastly, transecting satellites were found to contribute thousands of additional conjunctions outside of their operational orbit, and may require up to an additional 5% in CA maneuver fuel allocation.

conjunction assessment↗

NASA NTRS2010

Analyzing the Effect of Intraseasonal Meteorological Variability and Land Cover on Aerosol-Cloud Interactions During the Amazonian Biomass Burning Season

High resolution aerosol, cloud, water vapor, and atmospheric profile data from the Moderate Resolution Imaging Spectroradiometer (MODIS) are utilized to examine the impact of aerosols on clouds during the Amazonian biomass burning season in Rondnia, Brazil. It is found that increasing background column water vapor (CWV) throughout this transition season between the Amazon dry and wet seasons exerts a strong effect on cloud properties. As a result, aerosol-cloud correlations should be stratified by column water vapor to achieve a more accurate assessment of the effect of aerosols on clouds. Previous studies ignored the systematic changes to meteorological factors during the transition season, leading to possible misinterpretation of their results. Cloud fraction is shown generally to increase with aerosol optical depth (AOD) for both low and high values of column water vapor, whereas the relationship between cloud optical depth (COD) and AOD exhibits a different relationship. COD increases with AOD until AOD approx. 0.25 due to the first indirect (microphysical) effect. At higher values of AOD, COD is found to decrease with increasing AOD, which may be due to: (1) the inhibition of cloud development by absorbing aerosols (radiative effect) and/or (2) a retrieval artifact in which the measured reflectance in the visible is less than expected from a cloud top either from the darkening of clouds through the addition of carbonaceous biomass burning aerosols or subpixel dark surface contamination in the measured cloud reflectance. If (1) is a contributing mechanism, as we suspect, then a linear relationship between the indirect effect and increasing AOD, assumed in a majority of GCMs, is inaccurate since these models do not include treatment of aerosol absorption in and around clouds. The effect of aerosols on both column water vapor and clouds over varying land surface types is also analyzed. The study finds that the difference in column water vapor between forest and pasture is not correlated with aerosol loading, supporting the assumption that temporal variation of column water vapor is primarily a function of the larger-scale meteorology. However, a difference in the response of cloud fraction to increasing AOD is observed between forest and pasture. This suggests that dissimilarities between other meteorological factors, such as atmospheric stability, are likely to have an impact on aerosol-cloud correlations between different land-cover types.

TenHoeve, J. E.↗

NASA NTRS2025

Generalized Reference Targeting for Spaceflight

For spaceflight programs to achieve some of the aggressive exploration initiatives such as visiting and landing on other celestial bodies, rendezvousing with other orbiting vehicles, and ultimately returning crew safely to Earth, an assortment of targeting algorithms to compute the necessary burns to strategically maneuver a spacecraft to a variety of destinations are required. Numerous examples exist, but rather than creating and implementing multiple targeting solutions, is it possible to have a general targeting model that can accommodate a variety of applications? Originally motivated for mission design and analysis purposes, this paper outlines a generalized reference targeting algorithm for spaceflight that may also have applications in on-orbit flight software. It accommodates arbitrary flight dynamics and both impulsive and finite burns for either absolute or relative targeting applications. It also allows for an arbitrary number of targeting design parameters such as multiple discrete correction burns or finite thrust parameters to satisfy numerous combinations of targeting constraints that can have fixed or variable time epochs. Given a reference trajectory, this targeting technique provides a general framework to quickly solve an assortment of targeting problems that may be well-defined, over-determined, or under-determined while naturally producing metrics providing insight into the controllability for a given problem formulation. Due to the derivation, speed, and accuracy of the algorithm, it lends to supporting rapid linear covariance analysis and robust trajectory design applications for a variety of flight phases such as rendezvous and docking, cislunar transfer, interplanetary flight, orbit maintenance, de-orbit, and powered descent and landing.

Targeting↗