Coupling High Fidelity Body Modeling with Non- Keplerian Dynamics to Design AIM-MASCOT-2 Landing Trajectories on Didymos Binary Asteroid

Similar documents
LOTNAV: A LOW-THRUST INTERPLANETARY NAVIGATION TOOL

Forward-Integration Riccati-Based Feedback Control for Spacecraft Rendezvous Maneuvers on Elliptic Orbits

Image Processing Adam Finkelstein & Tim Weyrich Princeton University

Introductory Thoughts

1 Course Notes in Introductory Physics Jeffrey Seguritan

COUPLING HIGH FIDELITY BODY MODELING WITH NON-KEPLERIAN DYNAMICS TO DESIGN AIM-MASCOT-2 LANDING TRAJECTORIES ON DIDYMOS BINARY ASTEROID

Eric Klein and Ning Sa

Introduction: A Generalized approach for computing the trajectories associated with the Newtonian N Body Problem

NGSS High School Physics Domain Model

Phys101 Final Code: 1 Term: 132 Wednesday, May 21, 2014 Page: 1

Faculty of Engineering and Department of Physics Engineering Physics 131 Midterm Examination February 27, 2006; 7:00 pm 8:30 pm

The influence of a semi-infinite atmosphere on solar oscillations

Technology, Dhauj, Faridabad Technology, Dhauj, Faridabad

Trajectory Design For A Visible Geosynchronous Earth Imager

The ultra-high energy cosmic rays image of Virgo A

3D FE Modeling Simulation of Cold Rotary Forging with Double Symmetry Rolls X. H. Han 1, a, L. Hua 1, b, Y. M. Zhao 1, c

THE INTERMEDIATE EXPERIMENTAL VEHICLE IN FLIGHT EXPERIMENTATION. FROM THE DESIGN TO THE FLIGHT.

THERMAL-VACUUM VERSUS THERMAL- ATMOSPHERIC TESTS OF ELECTRONIC ASSEMBLIES

Uncertainties in TRP Measurements Due to Finite Range Lengths

The MAVEN Mission: Exploring Mars Climate History. Bruce Jakosky LASP / University of Colorado

Study Group Report: Plate-fin Heat Exchangers: AEA Technology

Computational modeling techniques

ASTRODYNAMICS. o o o. Early Space Exploration. Kepler's Laws. Nicolaus Copernicus ( ) Placed Sun at center of solar system

Numerical Simulation of the Thermal Resposne Test Within the Comsol Multiphysics Environment

EXPERIMENTAL STUDY ON DISCHARGE COEFFICIENT OF OUTFLOW OPENING FOR PREDICTING CROSS-VENTILATION FLOW RATE

Fabrication Thermal Test. Methodology for a Safe Cask Thermal Performance

Aircraft Performance - Drag

Computational modeling techniques

Solution to HW14 Fall-2002

Statistics, Numerical Models and Ensembles

Comparison of hybrid ensemble-4dvar with EnKF and 4DVar for regional-scale data assimilation

Fall 2013 Physics 172 Recitation 3 Momentum and Springs

THE ASTEROID IMPACT MISSION: CONSOLIDATED MISSION ANALYSIS AND SCIENTIFIC PAYLOAD OPERATIONS AT BINARY ASTEROID DIDYMOS

14. Which shows the direction of the centripetal force acting on a mass spun in a vertical circle?

14. Which shows the direction of the centripetal force acting on a mass spun in a vertical circle?

Main Menu. SEG Houston 2009 International Exposition and Annual Meeting. Summary

Performance Bounds for Detect and Avoid Signal Sensing

Ice Accretion Prediction Code

DESIGN OPTIMIZATION OF HIGH-LIFT CONFIGURATIONS USING A VISCOUS ADJOINT-BASED METHOD

Math 302 Learning Objectives

Free Vibrations of Catenary Risers with Internal Fluid

Computational modeling techniques

Distributions, spatial statistics and a Bayesian perspective

The Kullback-Leibler Kernel as a Framework for Discriminant and Localized Representations for Visual Recognition

Increasing Heat Transfer in Microchannels with Surface Acoustic Waves*

ELECTRON CYCLOTRON HEATING OF AN ANISOTROPIC PLASMA. December 4, PLP No. 322

Kinematic transformation of mechanical behavior Neville Hogan

Kinetic Complexity in Divertor Modeling

AP CHEMISTRY CHAPTER 6 NOTES THERMOCHEMISTRY

Professional Development. Implementing the NGSS: High School Physics

Perfrmance f Sensitizing Rules n Shewhart Cntrl Charts with Autcrrelated Data Key Wrds: Autregressive, Mving Average, Runs Tests, Shewhart Cntrl Chart

Numerical model about hygro-lock effect in the mechano-sorptive behavior

Curvature Effects on Thermal Buckling Load of DWCNT Under Axial Compression Force

MANIPAL INSTITUTE OF TECHNOLOGY

Q x = cos 1 30 = 53.1 South

FIELD QUALITY IN ACCELERATOR MAGNETS

Part 3 Introduction to statistical classification techniques

Multiple Source Multiple. using Network Coding

POLARISATION VISUAL PHYSICS ONLINE. View video on polarisation of light

Simple Models of Foundation-Soil Interactions

New Angles-only Algorithms for Initial Orbit Determination

NUROP CONGRESS PAPER CHINESE PINYIN TO CHINESE CHARACTER CONVERSION

Module 4: General Formulation of Electric Circuit Theory

AP Physics Laboratory #4.1: Projectile Launcher

Optical property of Few-Mode Fiber with Non-uniform Refractive Index for Cylindrical Vector Beam Generation

ES201 - Examination 2 Winter Adams and Richards NAME BOX NUMBER

Chaotic behavior of the piccolo

Effect of Paramagnetic Ions in Aqueous Solution for Precision Measurement of Proton Gyromagnetic Ratio

Application of Moment of Momentum Theorem in Tie the Knot with A Hoop and A Chain

Pattern Recognition 2014 Support Vector Machines

EE247B/ME218: Introduction to MEMS Design Lecture 7m1: Lithography, Etching, & Doping CTN 2/6/18

Robust trajectory design using invariant manifolds. Application to (65803) Didymos

Checking the resolved resonance region in EXFOR database

Guaranteeing Reliability with Vibration Simulation and Testing. Dr. Nathan Blattau

Chapter 4. Unsteady State Conduction

Kinetics of Particles. Chapter 3

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System

CESAR Science Case The differential rotation of the Sun and its Chromosphere. Introduction. Material that is necessary during the laboratory

Chapter 3: Cluster Analysis

A Study on Pullout Strength of Cast-in-place Anchor bolt in Concrete under High Temperature

Chapter 2. Kinematics in One Dimension. Kinematics deals with the concepts that are needed to describe motion.

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System

and the Doppler frequency rate f R , can be related to the coefficients of this polynomial. The relationships are:

SEISMIC STABILITY ANALYSIS OF FOOTING ADJACENT TO SLOPES BY SLIP LINE METHOD

arxiv:hep-ph/ v1 2 Jun 1995

Yeu-Sheng Paul Shiue, Ph.D 薛宇盛 Professor and Chair Mechanical Engineering Department Christian Brothers University 650 East Parkway South Memphis, TN

Phys102 Final-061 Zero Version Coordinator: Nasser Wednesday, January 24, 2007 Page: 1

Q1. In figure 1, Q = 60 µc, q = 20 µc, a = 3.0 m, and b = 4.0 m. Calculate the total electric force on q due to the other 2 charges.

Phys. 344 Ch 7 Lecture 8 Fri., April. 10 th,

Conceptual Dynamics SDC. An Interactive Text and Workbook. Kirstie Plantenberg Richard Hill. Better Textbooks. Lower Prices.

Work, Energy, and Power

Reinforcement Learning" CMPSCI 383 Nov 29, 2011!

Chapter 2 GAUSS LAW Recommended Problems:

Evaluation of the Humboldt GeoGauge TM on Dry Cohesionless Silica Sand in a Cubical Test Bin. Lary R. Lenke, P.E. Matt Grush R. Gordon McKeen, P.E.

Lecture 2: Supervised vs. unsupervised learning, bias-variance tradeoff

You need to be able to define the following terms and answer basic questions about them:

Internal vs. external validity. External validity. This section is based on Stock and Watson s Chapter 9.

Determining the Accuracy of Modal Parameter Estimation Methods

Examiner: Dr. Mohamed Elsharnoby Time: 180 min. Attempt all the following questions Solve the following five questions, and assume any missing data

Direct Monte Carlo Simulation of Time- Dependent Problems

Transcription:

6 th Internatinal Cnference n Astrdynamics Tls and Techniques (ICATT) Darmstadt 14 17 March 2016 Cupling High Fidelity Bdy Mdeling with Nn- Keplerian Dynamics t Design AIM-MASCOT-2 Landing Trajectries n Didyms Binary Asterid Fabi Ferrari*, Michèle Lavagna*, Ian Carnelli * Plitecnic di Milan, Department f Aerspace Science and Technlgy Eurpean Space Agency, General Study Prgram Office, ESA HQ

Cntent 2 Intrductin: AIM and MASCOT-2 landing Dynamical envirnment arund Didyms system Gravity mdel Estimate f dynamical perturbatins Design f MASCOT-2 ballistic descent Requirements Assumptins Escape velcities in the three-bdy system Sensitivity analyses Successful landing prbability Landing dispersin at rest Fabi Ferrari

Asterid Impact Missin (AIM) 3 65803 Didyms will transit near Earth (less than 0.1 AU) in late 2022 [ Didyms = Didymain + Didymn ] Credits: ESA AIDA - Asterid Impact & Deflectin Assessment AIM - Asterid Impact Missin (ESA) Gals: Study f binary system Deplyment f a lander (MASCOT-2) Deplyment f cubesats (COPINS) DART - Duble Asterid Redirectin Test (NASA) Gals: High velcity kinetic impact Fabi Ferrari

MASCOT-2 landing highlights 4 MASCOT-2 shall land n the smaller asterid (Didymn) MASCOT-2 is a passive prbe N actuatrs fr trajectry cntrl N actuatrs r devices fr anchring t the surface Extremely lw gravity field n Didymn s surface MASCOT-2 can bunce but it shall stay n the surface MASCOT-2 release pint shall be safe frm AIM spacecraft pint f view Fabi Ferrari

Didyms system 5 1.18 km Credits: ESA Diameter [m] Mass [kg] Didymain 775 5.2 e11 Didymn 163 4.8 e9 Tw rders f magnitude lwer than 67/P C-G (Rsetta/Philae) Fabi Ferrari

Didyms system 6 1.18 km Credits: ESA Diameter [m] Mass [kg] Didymain 775 5.2 e11 Didymn 163 4.8 e9 Fur rders f magnitude lwer than 67/P C-G (Rsetta/Philae) + Didymain s perturbatin Fabi Ferrari

Didyms system Gravity mdel 7 1.18 km MASCOT-2 dynamics Three-bdy system L* 1.18 km µ 9.2 e-3 Credits: ESA T 11.92 h Diameter [m] Mass [kg] Didymain 775 5.2 e11 Didymn 163 4.8 e9 CR3BP t mdel the mtin f the asterids center f mass Fabi Ferrari

Didyms system Gravity mdel 8 CR3BP t mdel the mtin f the asterids center f mass Gravity surce Didymain (sphere) + Didymn (sphere) Didymain (sphere) + Didymn (ellipsid) Didymain (plyhedrn) + Didymn (ellipsid) Dynamical mdel S1+S2 (CR3BP) S1+E2 P1+E2 S1+E2 P1+E2 Fabi Ferrari

Didyms system Dynamical perturbatins 9 Relevant effects at 1 km [m/s 2 ] Acceleratin due t Didyms system s gravity e-5 Acceleratin due t SRP n MASCOT-2 lander e-8 Perturbatin due t Sun s gravity (third bdy) e-11 Perturbatin due t Earth s gravity (third bdy) e-13 Nt relevant fr MASCOT-2 scenari Irregularities in Didyms gravity field are the mst relevant perturbatins N relevant differences fund between S1+E2 and P1+E2 mdeling Fabi Ferrari

MASCOT-2 landing design Requirements 10 MASCOT-2 shall land n the smaller asterid (Didymn) MASCOT-2 is a passive prbe N actuatrs fr trajectry cntrl N actuatrs r devices fr anchring t the surface Extremely lw gravity field n Didymn s surface MASCOT-2 can bunce but it shall stay n the surface MASCOT-2 release pint shall be safe frm AIM spacecraft pint f view Ballistic descent with n escape after buncing v aftertd < v esc Fabi Ferrari

Escape velcity frm Didymn s surface Theretical limits 11 Dynamical mdel v esc [cm/s] S2 7.70 S1+S2 4.95 S1+E2 4.58 P1+E2 4.57 Values t escape Didymn s SOI (r neighbrhd) It depends n latitude they depend n latitude and lngitude The escape velcity is the tuch dwn velcity fr a pure ballistic landing frm the SOI v TD v esc Minimum energy t land / escape: thrugh L1 pint Fabi Ferrari

Tuch-dwn velcity Mnte Carl backwards integratin 12 Energy v TD < v L1 v L1 < v TD < v L2 v TD > v L2 Minimum v td [cm/s] Thrugh L1 neck (v L1 ) Thrugh L2 neck (v L2 ) S1+S2 4.95 5.23 S1+E2 4.58 5.11 P1+E2 4.57 5.11 Fabi Ferrari

MASCOT-2 landing design Trajectry selectin 13 Requirements Ballistic descent with n escape after buncing v aftertd < v esc Safety: release far frm Didymn cnstraints the dynamics t v TD v esc Design strategy Keep v TD as lw as pssible Trajectries with lw v TD are thse passing thrugh the L1 r L2 neck with lw energy (L2 neck is cnsidered fr safer release) Maniflds assciated t L2 pint Fabi Ferrari

MASCOT-2 landing design Example f landing slutin 14 Release frm 200 m altitude Manifld High time f flight Nt rbust after release dispersin Fabi Ferrari

MASCOT-2 landing design Example f landing slutin 15 Release frm 200 m altitude Near the manifld Lw time f flight Rbust after release dispersin Fabi Ferrari

MASCOT-2 landing design Assumptins 16 BALLISTIC DESCENT: MASCOT-2 dispersin at release Navigatin errr Release mechanism errr 5 deg (1σ) half cne angle arund nminal release directin 0.5 cm/s (1σ) uncertainty in release velcity BOUNCING DYNAMICS: velcity after tuch dwn Directin accrding t sil inclinatin Unifrm distributin in azimuth Gaussian distributin in elevatin: 90 ± 70 deg (µ ± 3σ) Nrm f velcity dumped accrding t the restitutin cefficient η = v aftertd v befretd = 0.9 0.6 = 0.54 asterid s sil lander s structure Fabi Ferrari

MASCOT-2 landing design Escape prbability 17 Frm release up t escape r stay n Didymn s surface Mnte Carl simulatin (200000 cases fr each release pint) S1+E2 mdel used (equivalent t P1+E2) Release frm 200 m altitude Release frm 250 m altitude Fabi Ferrari

MASCOT-2 landing design Escape prbability 18 Frm release up t escape r stay n Didymn s surface Mnte Carl simulatin (200000 cases fr each release pint) S1+E2 mdel used (equivalent t P1+E2) Release altitude [m] Escaped trajectries [%] After release After tuch dwn TOT 100 1.14 0.00 1.14 150 1.26 0.00 1.27 200 3.18 0.01 3.19 250 5.79 0.04 5.83 300 6.03 0.32 6.35 Fabi Ferrari

MASCOT-2 landing design Landing dispersin at rest 19 Frm release up t rest n Didymn s surface: utging vertical velcity lwer than 0.5 cm/s Mnte Carl simulatin (200000 cases fr each release pint) S1+E2 mdel used (equivalent t P1+E2) Release altitude [m] Dispersin at rest (µ ± 3σ) Latitude [deg] Lngitude [deg] Tf [h] 100 0.2 ± 32.5 23.3 ± 87.9 1.81 ± 0.95 150 0.0 ± 31.5 20.1 ± 66.7 2.19 ± 1.06 200 0.1 ± 32.9 19.7 ± 60.5 2.50 ± 1.21 250-0.1 ± 36.6 20.5 ± 64.8 2.77 ± 1.41 300 0.1 ± 44.3 19.5 ± 86.0 2.95 ± 1.35 Fabi Ferrari

MASCOT-2 landing design Landing dispersin at rest: release frm 200 m 20 Release altitude [m] Dispersin at rest (µ ± 3σ) Latitude [deg] Lngitude [deg] Tf [h] 200 0.1 ± 32.9 19.7 ± 60.5 2.50 ± 1.21 Fabi Ferrari

MASCOT-2 landing design Landing dispersin at rest: release frm 200 m 21 Release altitude [m] Dispersin at rest (µ ± 3σ) Latitude [deg] Lngitude [deg] Tf [h] 200 0.1 ± 32.9 19.7 ± 60.5 2.50 ± 1.21 Fabi Ferrari

MASCOT-2 landing design Cnclusin 22 Summary The main perturbing actin is due t the uncertainty n Didyms/Didymn gravity field Higher successful landing prbability fr lwer release pints Landing dispersin at rest is cnfined within a certain latitude-lngitude band Gd results with current assumptins n release dispersin, sil uncertainty and restitutin cefficient Release altitude [m] Successful prbability [%] 100 98.86 150 98.73 200 96.81 250 94.17 300 93.65 Fabi Ferrari

6 th Internatinal Cnference n Astrdynamics Tls and Techniques (ICATT) Darmstadt 14 17 March 2016 Cupling High Fidelity Bdy Mdeling with Nn- Keplerian Dynamics t Design AIM-MASCOT-2 Landing Trajectries n Didyms Binary Asterid Fabi Ferrari*, Michèle Lavagna*, Ian Carnelli * Plitecnic di Milan, Department f Aerspace Science and Technlgy Eurpean Space Agency, General Study Prgram Office, ESA HQ