Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation Use Case”

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 145 records · Page 8

Application of the ENS3DAE Euler/Navier-Stokes Aeroelastic method

This paper discusses recent applications of the ENS3DAE computational aeroelasticity method. In particular, it describes aeroelastic and unsteady aerodynamics calculations performed on wings with trailing edge control surfaces. These simulations include the investigation of control reversal for a structurally flexible wing with a deflected control surface, and a static and dynamic analysis of a rigid wing with an oscillating control surface. The two sets of calculations were performed independently on different wings using different grid topologies. The control reversal simulation represents an inviscid Euler static aeroelastic analysis of a thin wing with a rectangular planform. The geometry of this wing makes it suitable for computations using more approximate, inviscid aerodynamics methods. Thus, the results of the present Euler computations are compared with numerical data generated by a validated computational aeroelasticity code which uses a simpler aerodynamic formulation. The second illustrated case involves the simulation of a significantly more complex flowfield and the static and dynamic analyses of this geometry were performed using the viscous Navier-Stokes equation option in ENS3DAE. Results of both the steady and unsteady calculations on this wing are compared with existing experimental data.

Schuster, David M.↗

Formationkeeping of Spacecraft via Differential Drag

This paper examines the use of differential drag in the formationkeeping of spacecraft. In many future space missions one satellite will be required to fly in a specific position with respect to another satellite; this action is referred to as formationkeeping. In this study, differential drag is the difference in drag between the two satellites. Reasons to use differential drag as an actuator for formationkeeping include the avoidance of jet plume impingement effects on closely spaced satellites and possible fuel savings. The equations of relative motion between the two satellites are derived and a mathematical transformation is made to reduce the formationkeeping problem to the simultaneous solution of a double integrator and a harmonic oscillator. A two-part control law is developed that simultaneously and dependently solves the double oscillator and harmonic integrator. Results are presented of eight test cases being driven to a target position; two different simulations are used. The validity of assumptions made in the derivation of the control law is examined in the comparison of similar test cases run through different simulations. The control law developed can drive a satellite from an initial position to a target position and maintain the satellite at that location.

Leonard, Carolina Lee↗

Case studies in Aircraft Parameter Identification

Three case studies in Aircraft Parameter Identification using simulated data for X-22 VTOL aircraft and flight test data for HL-10 and M2/F3 lifting bodies. After a brief discussion of the previous techniques and their limitations, a technique based on the Maximum Likelihood criterion is described. The problems of identifiability and uniqueness in determining the Stability and Control derivatives from flight test data are discussed and several methods for alleviating these problems are presented. The flight test data is analyzed in several different ways for obtaining physically meaningful estimates for the aircraft parameters.

Mehra, R. K.↗

Some issues related to simulation of the tracking and communications computer network

The Communications Performance and Integration branch of the Tracking and Communications Division has an ongoing involvement in the simulation of its flight hardware for Space Station Freedom. Specifically, the communication process between central processor(s) and orbital replaceable units (ORU's) is simulated with varying degrees of fidelity. The results of investigations into three aspects of this simulation effort are given. The most general area involves the use of computer assisted software engineering (CASE) tools for this particular simulation. The second area of interest is simulation methods for systems of mixed hardware and software. The final area investigated is the application of simulation methods to one of the proposed computer network protocols for space station, specifically IEEE 802.4.

Lacovara, Robert C.↗

Ku-band interferences for the Space Station multiple access link

The multiple access (MA) system, an essential part of the United States Space Station, will be used to provide simultaneous communication services, voice and command, telemetry and video, between the Space Station and all the different coorbiting user vehicles. The system is designed to operate in the Ku-band frequencies between 14 GHz and 14.89 GHz. Within this band, there are other users that could pose serious interference potential. The extent of the interference and its impact on the system operations are analyzed using a Monte Carlo simulation. Worst-case analysis indicated that certain links are very susceptible to interference. Mitigation techniques involving both hardware changes and operation constraints are presented. With the modification of some of the design and implementation of some of the mitigation techniques, the interference levels were found to be only a few percent and to occur in a burst manner.

Loh, Y. C.↗

Updated Results of Deterministic Reconfigurable Control Design for the X-33 Vehicle

In the event of a control surface failure, the purpose of a reconfigurable control system is to redistribute the control effort among the remaining working surfaces such that satisfactory stability and performance are retained. An Off-line Nonlinear General Constrained Optimization approach was used for the reconfigurable X-33 control design method. Three examples of failure are shown using a high fidelity 6 DOF simulation (case 1: ascent with a left body flap jammed at 25 deg.; case 2: entry with a right inboard elevon jam at 25 deg. and case 3: landing (TAEM) (Terminal Area Energy Management) with a left rudder jam at -30 deg.) Failure comparisons between responses with the nominal controller and reconfigurable controllers show the benefits of reconfiguration. Single jam aerosurface failures were considered, and failure detection and identification is considered accomplished in the actuator controller. The X-33 flight control system will incorporate reconfigurable flight control in the baseline system.

Cotting, M. Christopher↗

Coupled Model Simulation of Snowfall Events Over the Black Hills

Although many long-term simulations of snow accumulation and oblation have been made using stand-alone land surface models and surface models coupled with GCMs, less research has focused on short-term event simulations. Actually, accurate event simulations of snow-related processes are the basis for successful long-term simulation. Three advantages of event simulations of snowfall and snow melting are availability of: (1) intensive observation data from field experiments for validation; (2) more physically-realistic precipitation schemes for use in atmospheric models to simulate snowfall; and (3) a more detailed analysis of the snow melting processes. In addition to the complexities of snow related processes themselves, terrain-induced effects on snowfall/snow melting make simulations of snow events more difficult. Climatological observations indicate that terrain features such as the Black Hills of South Dakota and Wyoming can exert important effects on snow accumulation and snow oblation processes. One of the primary effects is that the orography causes forced uplift of airflow and causes atmospheric waves to form both upwind and downwind of it. Airflow often splits around the obstacle, converging on the lee side. This convergence may lead to precipitation enhancement. It also provides an elevated heat and moisture source that enhances atmospheric instability. During the period of April 5-May 5, 1999, the Upper Missouri River Basin Pilot Project (UMRBPP) made intensive observations on precipitation events occurring in the Black Hills. Two moderate snowfall events were captured during the period. The resulting high temporal and spatial resolution data provides opportunities to investigate terrain effects on snowfall amount, distribution, and melting. Successful simulation of snowfall amount, distribution, and evolution using atmospheric models is important to subsequent modeling of snow melting using snow sub-models in land surface schemes. In this paper, a coupled model system, consisting of an atmosphere model (ARPS) and a land-surface model (revised NCAR LSM), is used to simulate one of these UMRBPP cases.

Wang, Jianzhong↗

Full Envelope Reconfigurable Control Design for the X-33 Vehicle

In the event of a control surface failure, the purpose of a reconfigurable control system is to redistribute the control effort among the remaining working surfaces such that satisfactory stability and performance are retained. An Off-line Nonlinear General Constrained Optimization (ONCO) approach was used for the reconfigurable X-33 control design method. Three example failures are shown using a high fidelity 6 DOF simulation (case I ascent with a left body flap jammed at 25 deg.; case 2 entry with a right inboard elevon jam at 25 deg.; and case 3, landing (TAEM) with a left rudder jam at -30 deg.) Failure comparisons between responses with the nominal controller and reconfigurable controllers show the benefits of reconfiguration. Single jam aerosurface failures were considered, and failure detection and identification is considered accomplished in the actuator controller. The X-33 flight control system will incorporate reconfigurable flight control in the baseline system.

Cotting, M. Christopher↗

Aeroacoustic Simulation of a Nose Landing Gear in an Open Jet Facility Using FUN3D

Numerical simulations have been performed for a partially-dressed, cavity-closed nose landing gear configuration that was tested in NASA Langley s closed-wall Basic Aerodynamic Research Tunnel (BART) and in the University of Florida s open-jet acoustic facility known as UFAFF. The unstructured-grid flow solver, FUN3D, developed at NASA Langley Research center is used to compute the unsteady flow field for this configuration. A hybrid Reynolds-averaged Navier-Stokes/large eddy simulation (RANS/LES) turbulence model is used for these computations. Time-averaged and instantaneous solutions compare favorably with the measured data. Unsteady flowfield data obtained from the FUN3D code are used as input to a Ffowcs Williams-Hawkings noise propagation code to compute the sound pressure levels at microphones placed in the farfield. Significant improvement in predicted noise levels is obtained when the flowfield data from the open jet UFAFF simulations is used as compared to the case using flowfield data from the closed-wall BART configuration.

Vatsa, Veer N.↗

Current-State Constrained Filter Bank for Wald Testing of Spacecraft Conjunctions

We propose a filter bank consisting of an ordinary current-state extended Kalman filter, and two similar but constrained filters: one is constrained by a null hypothesis that the miss distance between two conjuncting spacecraft is inside their combined hard body radius at the predicted time of closest approach, and one is constrained by an alternative complementary hypothesis. The unconstrained filter is the basis of an initial screening for close approaches of interest. Once the initial screening detects a possibly risky conjunction, the unconstrained filter also governs measurement editing for all three filters, and predicts the time of closest approach. The constrained filters operate only when conjunctions of interest occur. The computed likelihoods of the innovations of the two constrained filters form a ratio for a Wald sequential probability ratio test. The Wald test guides risk mitigation maneuver decisions based on explicit false alarm and missed detection criteria. Since only current-state Kalman filtering is required to compute the innovations for the likelihood ratio, the present approach does not require the mapping of probability density forward to the time of closest approach. Instead, the hard-body constraint manifold is mapped to the filter update time by applying a sigma-point transformation to a projection function. Although many projectors are available, we choose one based on Lambert-style differential correction of the current-state velocity. We have tested our method using a scenario based on the Magnetospheric Multi-Scale mission, scheduled for launch in late 2014. This mission involves formation flight in highly elliptical orbits of four spinning spacecraft equipped with antennas extending 120 meters tip-to-tip. Eccentricities range from 0.82 to 0.91, and close approaches generally occur in the vicinity of perigee, where rapid changes in geometry may occur. Testing the method using two 12,000-case Monte Carlo simulations, we found the method achieved a missed detection rate of 0.1%, and a false alarm rate of 2%.

Carpenter, J. Russell↗

Qualitative Fault Isolation of Hybrid Systems: A Structural Model Decomposition-Based Approach

Quick and robust fault diagnosis is critical to ensuring safe operation of complex engineering systems. A large number of techniques are available to provide fault diagnosis in systems with continuous dynamics. However, many systems in aerospace and industrial environments are best represented as hybrid systems that consist of discrete behavioral modes, each with its own continuous dynamics. These hybrid dynamics make the on-line fault diagnosis task computationally more complex due to the large number of possible system modes and the existence of autonomous mode transitions. This paper presents a qualitative fault isolation framework for hybrid systems based on structural model decomposition. The fault isolation is performed by analyzing the qualitative information of the residual deviations. However, in hybrid systems this process becomes complex due to possible existence of observation delays, which can cause observed deviations to be inconsistent with the expected deviations for the current mode in the system. The great advantage of structural model decomposition is that (i) it allows to design residuals that respond to only a subset of the faults, and (ii) every time a mode change occurs, only a subset of the residuals will need to be reconfigured, thus reducing the complexity of the reasoning process for isolation purposes. To demonstrate and test the validity of our approach, we use an electric circuit simulation as the case study.

Hybrid Systems Diagnosis↗

Parasitic Interference in Nulling Interferometry

Nulling interferometry aims to detect faint objects close to bright stars. Its principle is to produce a destructive interference along the line of sight so that the stellar flux is rejected, while the flux of the off-axis source can be transmitted. In practice, various instrumental perturbations can degrade the nulling performance. Any imperfection in phase, amplitude or polarization produces a spurious flux that leaks to the interferometer output and corrupts the transmitted off-axis flux. One of these instrumental perturbations is the crosstalk phenomenon, which occurs because of multiple parasitic reflections inside transmitting optics, and/or diffraction effects related to beam propagation along finite size optics. It can include a crosstalk of a beam with itself, and a mutual crosstalk between different beams. This can create a parasitic interference pattern, which degrades the intrinsic transmission map - or intensity response - of the interferometer. In this context, we describe how this instrumental effect impairs the performance of a Bracewell interferometer. A simple formalism is developed to derive the corresponding modified intensity response of the interferometer, as a function of the two parameters of interest: the crosstalk level (or contamination rate) and the phase shift between the primary and secondary - parasitic - beams. We then apply our mathematical approach to a few scientific cases, both analytically and using the GENIESIM simulation software, adapted to handle coherent crosstalk. Our results show that a coherent crosstalk level of about 1 per cent implies a 20 per cent drop of the signal-to-noise ratio at most. Careful attention should thus be paid to reduce the crosstalk level inside an interferometric instrument and ensure an instrumental stability that provides the necessary sensitivity through calibration procedures.

circumstellar matter↗

Reduced Order Modeling for Transonic Aeroservoelastic Control Law Development

As aircraft become more flexible, aeroelastic considerations become increasingly important and complex, particularly for transonic flight where nonlinearities in the flow render linear analysis tools less effective. In order to analyze these aeroelastic interactions between the fluid and the structure efficiently, reduced order models (ROMs) are sometimes generated from and used in place of computational fluid dynamics solutions. In this paper, several aerodynamic ROMs are generated and coupled with structural models to form aeroelastic ROMs. The aerodynamic ROMs generated here include the effects of control surface motion. Hence, the aeroelastic ROMs presented here are appropriate for use in aeroservoelastic applications and are intended to be used for aeroservoelastic control law development. These ROMs are used to simulate a number of test cases with and without control surface involvement. Results show that several of the ROMs generated in the paper are able to predict results similar to solutions of higher-order computational methods.

Waite, Josiah M.↗

Computational Aerodynamic Analysis in Support of the CRM Tail Cone Thruster Configuration Wind Tunnel Test

NASA’s Advanced Air Transport Technology (AATT) project is breaching the boundaries of aircraft design in pursuit of eco-friendly solutions that are compatible with urban noise comfort levels. Boundary layer ingesting (BLI) propulsion systems promise to reduce fuel burn with additional potential benefits in noise reduction. Type-II BLI systems of the STARC-ABL type are the subject of a test campaign planned for fiscal year 2022 in the National Transonic Facility (NTF), for which a CRM-based model with a retrofitted tail cone thruster (TCT) has been designed. The present work is a precursor to the NTF test, where the 240 cases planned for the experiment were simulated using the Launch, Ascent and Vehicle Aerodynamics (LAVA)computational framework. Solution sensitivity to angle of attack, engine operating conditions, and the presence of the supporting structure (sting) in the wind tunnel test are analyzed in the extensive dataset. The main flow features contributing to the inlet distortion are identified as the vertical tail wake, wing downwash and fuselage upsweep vortices, with the latter two experiencing the greatest sensitivity to angle of attack. Finally, results from an inlet-guide-vane(IGV) design/integration study are presented. The LAVA team and the turbomachinery design team at NASA’s Ames and Glenn Research Centers are collaborating in an effort to reduce flow distortion upstream of the fan by means of integrating an IGV system into the CRM-TCT model. A significant improvement in flow distortion metrics has been achieved since the initial design iteration. Results employing an actuator zone model with realistic radially-varying thrust profiles to simulate first-order fan effects within LAVA are presented.

AATT↗

A Radiative Transfer Simulator for PACE: Theory and Applications

A radiative transfer simulator was developed to compute the synthetic data of all three instruments onboard NASA’s Plankton Aerosol, Cloud, ocean Ecosystem (PACE) observatory, at the top of the atmosphere (TOA). The instrument suite includes the ocean color instrument (OCI), the HyperAngular Rainbow Polarimeter 2 (HARP2), and the Spectro-Polarimeter for Planetary Exploration 1 (SPEXone). The PACE simulator is wrapped around a monochromatic radiative transfer model based on the successive order of scattering (RTSOS), which accounts for atmosphere and ocean coupling, polarization, and gas absorption. Inelastic scattering, including Raman scattering from pure ocean water, fluorescence due to chlorophyll, and colored dissolved organic matter (CDOM), is also simulated. This PACE simulator can be used to explore the sensitivity of the hyperspectral and polarized reflectance of the Earth system with tunable atmosphere and ocean parameters, which include aerosol and cloud number concentration, refractive indices, and size distribution, ocean particle microphysical parameters, and solar and sensor-viewing geometry. The PACE simulator is used to study two important case studies. One is the impact of the significant uncertainty in pure ocean water absorption coefficient to the radiance field in the ultraviolet (UV) spectral region, which can be as much as 6%. The other is the influence of different amounts of brown carbon aerosols and CDOM on the polarized radiance field at TOA. The percentage variation of the radiance field due to CDOM is mostly for wavelengths smaller than 600 nm, while brown aerosols affect the whole spectrum from 350 to 890 nm, primarily due to covaried soot aerosols. Both case studies are important for aerosol and ocean color remote sensing and have not been previously reported in the literature.

PACE↗

Validated Reduced Computational Methods for Realistic RDC Injection Modeling

This paper outlines a simplified computational method without chemistry for single and multi-injector performance for rotating detonation combustion (RDC) environment as well as relevant performance metrics. For non-reactive cases, an unsteady pressure profile from a Method of Characteristics (MOC) simulation was used to simulate an RDC combustion wave travelling around the annulus. In reactive cases, a pre-mixed stoichiometric mixture of hydrogen and air is injected into a 2D unwrapped combustion chamber with a specified injector area. The performance is evaluated based on several performance metrics including unsteady discharge coefficient, average mass flux coming into the combustor, and percentage time obstructed. Finally, the computational method was validated through an experimental RDC test setup with a single element injector in a hydrogen-air detonation driver environment.

Propulsion↗

Probabilistic Modeling of a Three-Stage Human Landing System Architecture

Space Policy Directive-1 has led to NASA partnerships with commercial entities on procurement which includes the development of the Human Landing System (HLS) [1]. With the goal of delivering human crew to the lunar surface by 2024, system uncertainties become an important obstacle to the maturation of multiple new, driving technologies and mission concepts of the HLS program. As unmitigated uncertainties have previously led to failed development programs, these risks and their impacts must be understood and handled to ensure program success [2]. Sources of uncertainty include novel engine designs and configurations, increased reliance on cryogenic fluid management(CFM), and refueling technologies—which propagate as high-level performance metrics such as overall propellant mass and engine performance. Also, the occurrence of operational uncertainties—e.g. launch conditions or need to abort during the mission—can cause cascading effects on the rest of the mission that are difficult to definitively quantify, and are outside the scope of control. These concrete examples and other occurrences can be categorized as either epistemic or aleatory uncertainties.Epistemic uncertainty arises due to a lack of knowledge and can be alleviated with design and program maturation. Aleatory uncertainty is due to the inherent randomness of the system and cannot be directly reduced, unlike epistemic uncertainty. Robust design and probabilistic methods can compensate for aleatory effects. A taxonomy of uncertainty is referred to for this work [3]. In this paper, a probabilistic methodology to handle uncertainties has been demonstrated on a three-element HLS concept [1, 4], which allows tracking of current best estimates of the concept and assessment of concept design robustness against uncertainties. A sample case has been completed for this abstract, and an expansion on the methodology will be included in the final paper. This methodology has two key parts: first, the creation of a dynamic architecture model of a three-element HLS concept; and second, its use with surrogate modeling and range estimating techniques to capture and propagate uncertainties. This abstract will cover the basics of the approach used, and further details and justifications will be in the final paper.The mission profile associated with this three-element concept (Fig 1) was modeled as a set of mission events that facilitated mass changes, idles, or spacecraft maneuvers. The mission profile scope starts with each element’s NRHO orbit insertion and aggregation and ends at post-sortie rendezvous with Orion. More detail on the mission profile will be in the final paper. The DYnamic Rocket EQuation Tool (DYREQT), a space systems synthesis and sizing framework used by NASA, was used as the physics framework to model the HLS architecture for applying the probabilistic methodology [5, 6]. Specifically, a parametric representation of the lander, ascent, and transfer elements and the mission profile of each element was established, with vehicle and mission parameters available as inputs to allow for a dynamic model. Each vehicle stage was modeled with high-level performance metrics, using Isp and propellant mass fraction (PMF) to remain parametric. For the probabilistic analysis, uncertainties of interest within the HLS concept were enumerated and represented as parameters within the DYREQT model as inputs for vehicle stages or mission profile events. These parameters were frozen at their nominal values for the purposes of baselining architecture performance and sizing the vehicle appropriately based on reference documentation [1]. Range estimating—a probabilistic method that combines Monte Carlo sampling, focus on critical parameters, and heuristics to assess risk and opportunities—is traditionally used with Mass Equipment Lists (MELs), but has been adapted with operational parameters as well as vehicle parameters in theDYREQT model to capture mission uncertainty alongside vehicle uncertainty [7, 3]. This method was selected due to its application and insight on a system from a bottom-up perspective, independence from historical rules of thumb, and ability to generate sensitivities based on design decisions and uncertainties. As a sample case for the abstract, the boiloff rates of the vehicle elements and the loiter times during the mission (simulating launch time variations and changing window of opportunities) were used with range estimating to provide preliminary results. To perform the range estimation portion of this methodology (depicted in Fig. 3, further details in final paper), the DYREQT model was sampled using a Design of Experiments (DoE) to efficiently explore the architecture design space with respect to the sample set of uncertainty parameters; 5,000 cases via Latin Hypercube Sampling were computed on the DYREQT architecture model. Then, the results were used to create surrogate models, multivariate regressions that can visualize hypercube trends in the design space, of the architecture with respect to the uncertainty parameters. Range estimating was applied to the surrogates instead of the actual models, which saves computational expense due to the bulk of cases needed for the Monte Carlo simulation as part of range estimating. Uncertainty parameters were sampled independently from triangular distributions using the DoE ranges as ‘min’ and ‘max’, and the nominal value as ‘most likely’. Based engineering intuition, some uncertainty parameters are correlated—e.g. if the main propellant has a high boil-off rate, the oxidizer should follow suit as both are related to CFM technology.While a Monte Carlo simulation samples all inputs as independent, the results would show model correlations; thus, it is efficient to sample the inputs as correlated. Using a correlation matrix constructed for the uncertainty parameters, previously independent samples were transformed to perform a Correlated Monte Carlo. A table for the DoE ranges and probability distribution parameters is shown in Table 1, and more details on Correlated Monte Carlo Simulations will be discussed in the final paper. The model’s resulting DoE showed that multivariate polynomial equations fit via least squares method captured its behavior accurately for the sample case. For the Correlated Monte Carlo Simulation, a positive correlation between fuel and oxidizer boiloff rates was used as a demonstration. 10,000 cases were computed with the surrogates and the launched masses for each vehicle element was collated. The results can be displayed in a probability density function (PDF), showing the impact of the uncertainty parameters chosen. Integrating the PDFs will yield a cumulative distribution function (CDF) that shows the cumulative probability of a given value on the x-axis. For the sample case, the elements’ launch mass margin was calculated and represented in as CDFs, as a demonstrated representation of figures of merit for the HLS concept. For the lander and ascent elements, the NRHO mass insertion limit is 16t; the transfer element has a limit of 30t [1]. It can be seen with Figure 2 that this probabilistic methodology can provide insight into mass margin with respect to the uncertainties being modeled. Currently, the results show that the lander (descent) vehicle element has the most restrictive design space; it is the only element to show a 10% probability of negative margin. Further analysis on the Monte Carlo results will show sensitivities for driving constraints and parameters for architecture feasibility, which can lead to establishing potential mission rules.The combination of range estimating with a parametric architecture model for HLS demonstrated the capability of this probabilistic methodology in a sample case. As the HLS development progresses, this methodology has the potential for keeping current best estimates of architecture performance for awarded concepts due to the flexibility in DYREQT’s modeling framework and its parametric nature. Concept maturation and increased epistemic knowledge can be injected into the model probabilistic modeling, and thus continue to track probability of mission success.

Stephanie Y Zhu↗

Estimation of Timing Performance with Pulsed Laser Links for Small Satellite PNT

Satellite constellations are growing in size and in utility for applications as diverse as providing reliable, low-latency internet service to rural areas orEarth science missions.This expansion of satellite constellations brings attention to the need for responsive and reliable satellite communication.Current research in the PrecisionSpace Systems Laboratory at the University of Florida focuses on using pulsed laser communications as the method of relaying time signals instead of the radio frequency (RF)signals.The upcoming CLICK (CubeSat Laser Infrared CrosslinK) missions (Massachusetts Institute of Technology, University of Florida, NASA Ames Research Center) will demonstrate a spacecraft-to-spacecraft laser communication link and time-transfer.For pre-flight technology assessment of missions such as CLICK, a simulation of a spacecraft constellation in Earth orbit has been created to estimate clock synchronization and precision orbit determination based on measured instrumentation performance.We have developed a novel MATLAB-based numerical simulator to model spacecraft-to-spacecraft laser time-transfer and estimate the offset between the spacecraft clocks over time.This simulation includes timing errors associated with laser pulse generation and detection, as well as non-Gaussian clock drift models.The two on-board clocks modeled are a cesium-based Chip-Scale Atomic Clock (CSAC) and a rubidium-based Miniature Atomic Clock (MAC), both produced by Microchip.An example case of two spacecraft in a circular, low-Earth orbit receiving GPS position, GPS timing, and laser pulse time-of-flight measurements is simulated.The positions and velocities of the two spacecraft at a reference epoch and the constant clock model coefficients are estimated.Polynomial models of different orders were used as clock models. The effect of clock model order on the root-mean-square (RMS) of the clock error is apparent in the case using GPS and lasing measurements, showing that the clock model improves with increasing clock model order. Results compare the estimated clock model of a mission operation that only uses GPS measurements and a mission operation that uses both GPS and laser pulse time-of-flight measurements between spacecraft referenced to their on-board CSACs or on-board MACs. Including lasing measurements reduces the RMS clock model error by approximately 80% of the RMS of the cases with only GPS measurements.This simulation tool can be used to optimize the lasing operations schedule based on mission timing performance objectives. During this presentation, I will discuss the details of the simulation and clock model estimation, as well as discuss the results.I will compare the results of incorporating lasing timing measurements to using only GPS measurements and the outcomes of using different orders of polynomials.

Timing Performance↗