Results for “Extreme value analysis” — Page 7

DOE OSTI2025

Examining Clustered MCSs and Their Precipitation Significance Over Global Land MCS Hotspots

Multiple independent mesoscale convective systems (MCSs) can occur closely in space and time, subjecting common geographic regions to multiple precipitation events over short periods. Motivated by an earlier study highlighting significant hydrological impacts of clustered MCSs in the US, we extend the analysis to seven global continental MCS hotspots using a 20‐year global MCS data set. Our findings reveal that (a) 40%–80% of MCSs occur in clusters, depending on the region; (b) clustered MCSs exhibit larger rainfall areas and higher mean areal rainfall, with the highest values observed in Argentina and India; and (c) clustered MCSs contribute increasing fractions to more extreme 24‐ and 72‐hr precipitation events across all regions. Our results underscore the critical role of clustered MCSs in driving hydrological extremes over global land regions and highlight the need to better understand how their formation mechanisms differ from non‐clustered events across the globe.

Hu, Huancui [Pacific Northwest National Laboratory↗

DOE OSTI2024

DiffESM: Conditional Emulation of Temperature and Precipitation in Earth System Models With 3D Diffusion Models

Earth system models (ESMs) are essential for understanding the interaction between human activities and the Earth's climate. However, the computational demands of ESMs often limit the number of simulations that can be run, hindering the robust analysis of risks associated with extreme weather events. While low-cost climate emulators have emerged as an alternative to emulate ESMs and enable rapid analysis of future climate, many of these emulators only provide output on at most a monthly frequency. This temporal resolution is insufficient for analyzing events that require daily characterization, such as heat waves or heavy precipitation. We propose using diffusion models, a class of generative deep learning models, to effectively downscale ESM output from a monthly to a daily frequency. Trained on a handful of ESM realizations, reflecting a wide range of radiative forcings, our DiffESM model takes monthly mean precipitation or temperature as input, and is capable of producing daily values with statistical characteristics close to ESM output. Combined with a low-cost emulator providing monthly means, this approach requires only a small fraction of the computational resources needed to run a large ensemble. We evaluate model behavior using a number of extreme metrics, showing that DiffESM closely matches the spatio-temporal behavior of the ESM output it emulates in terms of the frequency and spatial characteristics of phenomena such as heat waves, dry spells, or rainfall intensity.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

How Flood Hazards in a Warming Climate Could Be Amplified by Changes in Spatiotemporal Patterns and Mechanisms of Water Available for Runoff

Prior research on climate change impacts on flooding has primarily focused on changes in extreme rainfall magnitudes, often neglecting snow processes and spatiotemporal storm patterns, such as hyetograph shapes and areal reduction factors (ARFs). This study examines projected changes in extreme water available for runoff (W) events in two snow-dominated basins in the western United States: the Yakima River Basin (YRB) in Washington State and the Walker River Basin (WRB) spanning the California-Nevada border. We analyze changes in W magnitudes, mechanisms, hyetograph shapes, and ARFs, and study their compounded impacts on flood hazard. Our findings suggest increased extreme W magnitudes across a large portion of the basins, with steeper or flatter hyetographs, and higher ARF values under the future climate. These changes are driven by a shift from seasonal snowmelt to more rain-on-snow events at higher elevations and by increased rainfall at lower elevations. We then use a single event-based rainfall-runoff model to estimate flood hazard changes based on extreme W magnitudes, hyetograph shapes, ARFs, and their compounded impacts. Our analysis reveals that focusing solely on the magnitude of changes in extreme W can significantly underestimate future flood hazards and uncertainties. Ignoring future changes in spatiotemporal patterns can underestimate future flood hazards by 63% and underestimate the uncertainty in future flood events by 18% in the WRB. These results underscore the necessity of incorporating spatiotemporal dynamics into future flood hazard assessments to provide a more accurate evaluation of potential impacts.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2023

Adaptive Independent Verification and Validation (IV&V) Reduces Risk of Software Impacting Safety in Artemis Missions

The National Aeronautics and Space Administration (NASA) is asking more of its human spaceflight programs than ever before through the collective Artemis Missions. The NASA Independent Verification and Validation (IV&V) Program contributes to NASA’s human spaceflight goals by providing IV&V services for NASA’s critical spacecraft and ground software. The IV&V Program is tasked with providing assurance from both individual and integrated mission software perspectives. The Artemis IV&V organization is actively supporting six distinct development efforts: Orion, the Space Launch System (SLS), Exploration Ground Systems (EGS), Mission Control Center (MCC), the Lunar Gateway, and the Human Landing System (HLS), representing a wide diversity of developer organizations, management structures, and development approaches. With much of this extremely complex flight and ground software being essential to human safety both on the ground and in space, Artemis IV&V is likewise challenged to provide more value-added assurance to future Artemis missions within a constrained budget. To meet this challenge, Artemis IV&V employs a variety of novel and evolving “Adaptive IV&V” approaches for planning and executing IV&V analysis to increase both the efficiency and effectiveness of the IV&V Program’s assurance activities, and to address the difficulties imposed by assuring software for a large, highly integrated, multi-mission enterprise managed and executed by physically and organizationally distinct programs. Instilling agile principles like iterative planning cycles, self-organizing teams, and regular retrospectives, into IV&V planning and execution has led to a more rapid turnaround of a minimum viable assurance product and allowed for increased alignment of assurance activities with development progress. Adopting an assurance case methodology has led to greater consistency and clearer communication of assurance design and provided a foundation for long-term maintenance of assurance plans, products, and results across missions. The IV&V-developed Assurance / Safety Case Analytical Network (A-SCAN) framework and tool has enabled the quantification and tracking of system/software risk and confidence. These confidence measures provide a means to repeatedly express the impact of planned and completed assurance work and the remaining residual risk. Applied as part of a “Follow-the-Risk” organizational ethos, this allows consistent rightsizing of analysis rigor and intensity commensurate with the perceived risk of defects, as well as appropriate targeting of the highest risk areas of the software to find safety issues before they can manifest. Finally, the development of the IV&V Advanced Risk Reduction Integrated Software Test and Operations Tri-program Lightweight Environment (ARRISTOTLE), an integrated software-only simulation of Orion, SLS, and EGS systems, has made it possible to independently test integrated pad and flight scenarios and inject faults to observe how the Artemis multi-program, mission software behaves in degraded modes and in response to hazards. These adaptive IV&V investments have enabled Artemis IV&V to become more efficient and effective in IV&V planning and execution and respond more readily to changes in the risk landscape, increasing the breadth and depth of risk reduction possible within the available resources. Residual risk tracking allows IV&V to communicate more effectively with stakeholders, both internal and external at all levels, and inform key decision-making personnel. This evolving assurance design approach provides IV&V surety that work is performed in the highest risk, most value-added areas of the software, to keep our astronauts and ground crews safe and ensure mission success.

Gerek A Whitman↗

NASA NTRS2023

Adaptive Independent Verification and Validation (IV&V) Reduces Risk of Software Impacting Safety in Artemis Missions

The National Aeronautics and Space Administration (NASA) is asking more of its human spaceflight programs than ever before through the collective Artemis Missions. The NASA Independent Verification and Validation (IV&V) Program contributes to NASA’s human spaceflight goals by providing IV&V services for NASA’s critical spacecraft and ground software. The IV&V Program is tasked with providing assurance from both individual and integrated mission software perspectives. The Artemis IV&V organization is actively supporting six distinct development efforts: Orion, the Space Launch System (SLS), Exploration Ground Systems (EGS), Mission Control Center (MCC), the Lunar Gateway, and the Human Landing System (HLS), representing a wide diversity of developer organizations, management structures, and development approaches. With much of this extremely complex flight and ground software being essential to human safety both on the ground and in space, Artemis IV&V is likewise challenged to provide more value-added assurance to future Artemis missions within a constrained budget. To meet this challenge, Artemis IV&V employs a variety of novel and evolving “Adaptive IV&V” approaches for planning and executing IV&V analysis to increase both the efficiency and effectiveness of the IV&V Program’s assurance activities, and to address the difficulties imposed by assuring software for a large, highly integrated, multi-mission enterprise managed and executed by physically and organizationally distinct programs. Instilling agile principles like iterative planning cycles, self-organizing teams, and regular retrospectives, into IV&V planning and execution has led to a more rapid turnaround of a minimum viable assurance product and allowed for increased alignment of assurance activities with development progress. Adopting an assurance case methodology has led to greater consistency and clearer communication of assurance design and provided a foundation for long-term maintenance of assurance plans, products, and results across missions. The IV&V-developed Assurance / Safety Case Analytical Network (A-SCAN) framework and tool has enabled the quantification and tracking of system/software risk and confidence. These confidence measures provide a means to repeatedly express the impact of planned and completed assurance work and the remaining residual risk. Applied as part of a “Follow-the-Risk” organizational ethos, this allows consistent rightsizing of analysis rigor and intensity commensurate with the perceived risk of defects, as well as appropriate targeting of the highest risk areas of the software to find safety issues before they can manifest. Finally, the development of the IV&V Advanced Risk Reduction Integrated Software Test and Operations Tri-program Lightweight Environment (ARRISTOTLE), an integrated software-only simulation of Orion, SLS, and EGS systems, has made it possible to independently test integrated pad and flight scenarios and inject faults to observe how the Artemis multi-program, mission software behaves in degraded modes and in response to hazards. These adaptive IV&V investments have enabled Artemis IV&V to become more efficient and effective in IV&V planning and execution and respond more readily to changes in the risk landscape, increasing the breadth and depth of risk reduction possible within the available resources. Residual risk tracking allows IV&V to communicate more effectively with stakeholders, both internal and external at all levels, and inform key decision-making personnel. This evolving assurance design approach provides IV&V surety that work is performed in the highest risk, most value-added areas of the software, to keep our astronauts and ground crews safe and ensure mission success.

Gerek Whitman↗

NASA NTRS1994

Zero potential vorticity envelopes for the zonal-mean velocity of the Venus/Titan atmospheres

The diagnostic analysis of numerical simulations of the Venus/Titan wind regime reveals an overlooked constraint upon the latitudinal structure of their zonal-mean angular momentum. The numerical experiments, as well as the limited planetary observations, are approximately consistent with the hypothesis that within the latitudes bounded by the wind maxima the total Ertel potential vorticity associated with the zonal-mean motion is approximately well mixed with respect to the neutral equatorial value for a stable circulation. The implied latitudinal profile of angular momentum is of the form M equal to or less than M(sub e)(cos lambda)(exp 2/Ri), where lambda is the latitude and Ri the local Richardson number, generally intermediate between the two extremes of uniform angular momentum (Ri approaches infinity) and uniform angular velocity (Ri = 1). The full range of angular momentum profile variation appears to be realized within the observed meridional - vertical structure of the Venus atmosphere, at least crudely approaching the implied relationship between stratification and zonal velocity there. While not itself indicative of a particular eddy mechanism or specific to atmospheric superrotation, the zero potential vorticity (ZPV) constraint represents a limiting bound for the eddy - mean flow adjustment of a neutrally stable baroclinic circulation and may be usefully applied to the diagnostic analysis of future remote sounding and in situ measurements from planetary spacecraft.

Allison, Michael↗

NASA NTRS1986

Solar irradiance variations from photometry of active regions

The Extreme Limb Photometer (ELP) has been used to measure the irradiance fluctuation of the sun due to selected active regions. Forty-five active regions that were completely scanned at various disk positions are included in the analysis. The contribution of these active regions to a global solar irradiance fluctuation has been correlated with photometric sunspot and facular indices (PSI and PFI) using published values of sunspot and calcium plage areas. The measured ELP fluctuations are converted to a global brightness fluctuation, Delta B/B. The sunspot component of Delta B/B correlates with PSI with r = 0.95. The facular component of Delta B/B correlates with PFI with r = 0.72. The expression for PFI is important to the question of energy balance between sunspots and faculae and the results presented here are not incompatible with energy balance between the two phenomena; that is the energy deficit of sunspots may be balanced by the energy excess of faculae.

Chapman, G. A.↗

NASA NTRS2011

Martian Atmospheric Modeling of Scale Factors for MarsGRAM 2005 and the MAVEN Project

For spacecraft missions to Mars, especially the navigation of Martian orbiters and landers, an extensive knowledge of the Martian atmosphere is extremely important. The generally-accepted NASA standard for modeling (MarsGRAM), which was developed at Marshall Space Flight Center. MarsGRAM is useful for task such as aerobraking, performance analysis and operations planning for aerobraking, entry descent and landing, and aerocapture. Unfortunately, the densities for the Martian atmosphere in MarsGRAM are based on table look-up and not on an analytical algorithm. Also, these values can vary drastically from the densities actually experienced by the spacecraft. This does not have much of an impact on simple integrations but drastically affects its usefulness in other applications, especially those in navigation. For example, the navigation team for the Mars Atmosphere Volatile Environment (MAVEN) Project uses MarsGRAM to target the desired atmospheric density for the orbiter's pariapse passage, its closet approach to the planet. After the satellite's passage through pariapsis the computed density is compared to the MarsGRAM model and a scale factor is assigned to the model to account for the difference. Therefore, large variations in the atmosphere from the model can cause unexpected deviations from the spacecraft's planned trajectory. In order to account for this, an analytic stochastic model of the scale factor's behavior is desired. The development of this model will allow for the MAVEN navigation team to determine the probability of various Martian atmospheric variations and their effects on the spacecraft.

Martian atmosphere↗

NASA NTRS1985

Cosmic ray charge and energy spectrum measurements using a new large area Cerenkov x dE/dx telescope

In September, 1981, a new 0.5 square meter ster cosmic ray telescope was flown to study the charge composition and energy spectrum of cosmic ray nuclei between 0.3 and 4 GeV/nuc. A high resolution Cerenkov counter, and three dE/dx measuring scintillation counters, including two position scintillators were contained in the telescope used for the charge and energy spectrum measurements. The analysis procedures did not require any large charge or energy dependent corrections, and absolute fluxes could be obtained to an accuracy approximately 5%. The spectral measurements made in 1981, at a time of extreme solar modulation, could be compared with measurements with a similar telescope made by our group in 1977, at a time of minimum modulation and can be used to derive absolute intensity values for the HEAO measurements made in 1979 to 80. Using both data sets precise energy spectra and abundance ratios can be derived over the entire energy range from 0.3 to greater than 15 GeV/nuc.

Webber, W. R.↗

NASA NTRS1986

An analysis on the shapes of ultraviolet extinction curves. I - The 2175 A bump

IUE data are used to investigate the properties of the 2175 A interstellar extinction bump toward 45 reddened Milky Way stars. Using analytic fitting functions, parameters are derived which describe the central position, FWHM, and strength of the bump. The bump position is very stable, with an extreme deviation which is still significantly larger than the measurement uncertainty, indicating true variability. The bump FWHM has a large range of intrinsic values, 0.77/micron to 1.25/micron. The normalized bump strength varies by more than a factor of three. There are no convincing correlations among the three bump parameters, but the width of the bump is strongly correlated with the dust grain environments. Dense quiescent regions yield broad bumps, while diffuse interstellar medium and regions of recent star formation yield narrower bumps. The absence of any correlation between bump width and central position argues strongly against graphite grains in the size range prescribed by the Mathis, Rumpl, and Nordsieck (1977) model as carriers of the bump.

Fitzpatrick, E. L.↗

NASA NTRS1996

Analysis of Cumulus Solar Irradiance Reflectance (CSIR) Events

Clouds are extremely important with regard to the transfer of solar radiation at the earth's surface. This study investigates Cumulus Solar Irradiance Reflection (CSIR) using ground-based pyranometers. CSIR events are short-term increases in solar radiation observed at the surface as a result of reflection off the sides of convective clouds. When sun-cloud observer geometry is favorable, these occurrences produce characteristic spikes in the pyranometer traces and solar irradiance values may exceed expected clear-sky values. Ultraviolet CSIR events were investigated during the summer of 1995 using Yankee Environmental Systems UVA-1 and UVB-1 pyranometers. Observed data were compared to clear-sky curves which were generated using a third degree polynomial best-fit line technique. Periods during which the observed data exceeded this clear-sky curve were identified as CSIR events. The magnitude of a CSIR event was determined by two different quantitative calculations. The MAC (magnitude above clear-sky) is an absolute measure of the difference between the observed and clear-sky irradiances. Maximum MAC values of 3.4 Wm(exp -2) and 0.069 Wm(exp -2) were observed at the UV-A and UV-B wavelengths, respectively. The second calculation determined the percentage above clear-sky (PAC) which indicated the relative magnitude of a CSIR event. Maximum UV-A and UV-B PAC magnitudes of 10.1% and 7.8%, respectively, were observed during the study. Also of interest was the duration of the CSIR events which is a function of sun-cloud-sensor geometry and the speed of cloud propagation over the measuring site. In both the UV-A and UV-B wavelengths, significant CSIR durations of up to 30 minutes were observed.

Laird, John L.↗

NASA NTRS1992

The Space Transportation System summer environment on launch pad

This paper describes a 2D flow and thermalanalysis to determine the solar effect on the Space Shuttle launch components subsequent to the external tank (ET) loading operation in extremely hot conditions. An existing CFD code Parabolic Hyperbolic or Elliptical Numerical Integration Code Series was used in the study. The analysis was done for a 2D slice between planes perpendicular to the longitudinal axis of the STS and passing through the lower portions of the Redesigned Solid Rocket Motors (RSRMs), the ET, and the wing of the Orbiter. The results are presented as local and average values of the heat transfer coefficient, and the Nusselt number, and the surface temperature around the RSRMs and the ET. Solar heating effects increased the surface temperatures of the RSRMs by 9-11 F. Higher prelaunch surface temperatures measured on the east and west RSRMs (in the inboard region between the RSRMs and the ET) during 19 most recent launches of the STS are correlated as a function of the ambient temperature.

Ahmad, R. A.↗

NASA NTRS2000

Hyperpolarized (129)Xe T (1) in oxygenated and deoxygenated blood

The viability of the new technique of hyperpolarized (129)Xe MRI (HypX-MRI) for imaging organs other than the lungs depends on whether the spin-lattice relaxation time, T(1), of (129)Xe is sufficiently long in the blood. In previous experiments by the authors, the T(1) was found to be strongly dependent upon the oxygenation of the blood, with T(1) increasing from about 3 s in deoxygenated samples to about 10 s in oxygenated samples. Contrarily, Tseng et al. (J. Magn. Reson. 1997; 126: 79-86) reported extremely long T(1) values deduced from an indirect experiment in which hyperpolarized (129)Xe was used to create a 'blood-foam'. They found that oxygenation decreased T(1). Pivotal to their experiment is the continual and rapid exchange of hyperpolarized (129)Xe between the gas phase (within blood-foam bubbles) and the dissolved phase (in the skin of the bubbles); this necessitated a complicated analysis to extract the T(1) of (129)Xe in blood. In the present study, the experimental design minimizes gas exchange after the initial bolus of hyperpolarized (129)Xe has been bubbled through the sample. This study confirms that oxygenation increases the T(1) of (129)Xe in blood, from about 4 s in freshly drawn venous blood, to about 13 s in blood oxygenated to arterial levels, and also shifts the red blood cell resonance to higher frequency. Copyright 2000 John Wiley & Sons, Ltd. Abbreviations used BOLD blood oxygen level dependent NOE nuclear overhouses effect PO(2) oxygen partial pressure RBC red blood cells RF radio frequency SNR signal-to-noise ratio.

NASA Discipline Life Support Systems↗

NASA NTRS2023

Reconstructing PM 2.5 Data Record for the Kathmandu Valley Using a Machine Learning Model

This paper presents a method for reconstructing the historical hourly concentrations of particulate matter 2.5 (PM2.5) over the Kathmandu Valley from 1980 to the present. The method uses a machine learning model that is trained using PM2.5 readings from US Embassy (Phora Durbar) as a ground truth, and the meteorological data from Modern-Era Retrospective Analysis for Research and Applications v2 (MERRA2) as input. The Extreme Gradient Boosting (XGBoost) model acquires a credible 10-fold cross-validation (CV) score of ~83.4%, an r2-score of ~84%, a Root Mean Square Error (RMSE) of ~15.82 µg/m3, and a Mean Absolute Error (MAE) of ~10.27 µg/m3. Further demonstrating the model's applicability to years other than those for which truth values are unavailable, the multiple cross-test with an unseen data set offered r2-scores for 2018, 2019, and 2020 ranging from 56% to 67%. The model-predicted data agrees with true values and indicates that MERRA2 underestimates PM2.5 over the region. It strongly agrees with ground-based evidence showing substantially higher mass concentrations in the dry pre- and post-monsoon seasons than in the monsoon months. It also shows a strong anti-correlation between PM2.5 concentration and humidity. The results also demonstrate that none of the years fulfilled the annual mean air quality index (AQI) standards set by the World Health Organization (WHO).

machine learning↗

NASA NTRS1975

A determination of the lunar moment of inertia

An estimate of the second zonal coefficient of the spherical-harmonic representation of the lunar gravitational field has been obtained from an analysis of particular orbital-element variations of the Explorer 35 and Explorer 49 spacecraft. Data from these spacecraft were used because the orbital configurations resulted in variations of the longitude of periapse and node which were, to first order, dependent only on the even zonal harmonics. The data time span for each satellite was extremely long: 2138 days for Explorer 35 and 230 days for Explorer 49. The value of the harmonic coefficient is determined and used to obtain a value of the lunar moment of inertia.

Gapcynski, J. P.↗

NASA NTRS1923

A Study of Taking off and Landing an Airplane

This report covers the results of an investigation carried on at the Langley Memorial Aeronautical Laboratory of the National Advisory Committee for Aeronautics for the purpose of discussing the various methods of effecting the take-off and the landing of an airplane, and to make a direct analysis of the control movements, the accelerations, and air speeds during these the maneuvers. The recording instruments developed at the laboratory were used in this test and the records obtained by them were made the basis for a comparative study of the two extreme methods of taking off (the tail-high and tail-low methods) and of various types of landings. The records should be of considerable value to a student pilot in enabling him to visualize the movements of the controls and the consequent effect upon the air speed and acceleration.

Carroll, T↗

DOE OSTI2024

Improving dairy manure hydrolysis and acidification through microbial community restructuring by adaptation to hyperthermophilic conditions

Dairy manure (DM) contributes significantly to greenhouse gas emissions and ecosystem degradation, yet its resistance to biodegradation hinders widespread bioprocessing applications. Lignocellulosic materials in DM pose a particular challenge because of their recalcitrance. Bioprocessing under hyperthermophilic (≥70 °C) conditions potentially offers an advantage over traditional fermentation temperatures due to enhanced activity of enzymes and the kinetics of enzymatic reactions. This can lead to a higher conversion rate and a greater extent of biomass hydrolysis and acidification. To test the validity of this hypothesis, the current study evaluated the efficacy of anaerobic hydrolysis and acidogenic fermentation of DM under mesophilic, thermophilic, and hyperthermophilic conditions. All inocula were adapted to corresponding temperatures but were derived from the same mesophilic source. Hyperthermophilic conditions resulted in superior DM hydrolysis efficiency (53%) compared to mesophilic (34%) and thermophilic (42%) conditions. The hyperthermophilic environment was particularly favorable to the decomposition of crude proteins and hemicellulose, which were reduced by 64% and 54%, respectively. Furthermore, hyperthermophilic fermentation also yielded the highest volatile fatty acid (VFA) production rate of 460 mg/L/day during the first four days, representing improvements of 50% and 90% over mesophilic and thermophilic conditions. In part, this was attributed to the enhanced production of branched-chain VFAs, including an increase of 6–10% in isobutyric acid and 12–13% in isovaleric acid. At hyperthermophilic conditions, however, there was no accumulation of VFAs during the days 5–8 of fermentation, which could be due to acetate conversion by the syntrophic acetate-oxidizing bacteria. A considerable gain in hydrolysis efficiency and VFA production rate were accompanied by a reduction in microbial diversity, which suggests that hyperthermophilic temperature is a favorable environment for the selection of organisms with enhanced DM hydrolysis and fermentation capabilities. A significantly increased relative abundance of xylanolytic Caldicoprobacter (23% of population) and proteolytic Thermovirga (9% of population) could be the major contributors to improved decomposition of hemicellulose and protein. As revealed by the techno-economic analysis, acidogenic fermentation of DM at 70 °C and a retention period of 4 days provides the greatest positive net present value, highest internal rate of return of 9.2%, and shortest investment payback period of 9 years. Furthermore, this study demonstrates that hyperthermophilic conditions enable superior deconstruction and bioconversion of lignocellulose-containing biomass into VFAs under reduced retention times, offering a promising approach for improving DM management and generating bioproducts.

09 BIOMASS FUELS↗

NASA NTRS2002

Formation Flying Satellite Control Around the L2 Sun-Earth Libration Point

This paper discusses the development of a linear control algorithm for formations in the vicinity of the L2 sun-Earth libration point. The development of a simplified extended Kalman filter is included as well. Simulations are created for the analysis of the stationkeeping and various formation maneuvers of the Stellar Imager mission. The simulations provide tracking error, estimation error, and control effort results. For formation maneuvering, the formation spacecraft track to within 4 meters of their desired position and within 1.5 millimeters per second of their desired zero velocity. The filter, with few exceptions, keeps the estimation errors within their three-sigma values. Without noise, the controller performs extremely well, with the formation spacecraft tracking to within several micrometers. Each spacecraft uses around 1 to 2 grams of propellant per maneuver, depending on the circumstances.

Hamilton, Nicholas H.↗