NONLINEAR MAPPING OF UNCERTAINTIES: A DIFFERENTIAL ALGEBRAIC APPROACH
|
|
- Laurence Kelly
- 5 years ago
- Views:
Transcription
1 NONLINEAR MAPPING OF UNCERTAINTIES: A DIFFERENTIAL ALGEBRAIC APPROACH R. Armellin 1, P. DiLizia 1, F. Bernelli-Zazzera 1, and M. Berz 2 1 Dipartimento di Ingegneria Aerospaziale, Politecnico di Milano, Via La masa 34, 2156 Milano, Italy 2 Department o Physics and Astronomy, Michigan State University, East Lansing, Michigan, 48824, USA ABSTRACT A method or the nonlinear propagation o uncertainties in celestial mechanics based on dierential algebra is presented. The arbitrary order Taylor expansion o the low o ordinary dierential equations with respect to the initial condition delivered by dierential algebra is exploited to implement an accurate and computationally eicient Monte Carlo algorithm, in which thousands o pointwise integrations are substituted by polynomial evaluations. The algorithm is applied to study the close encounter o asteroid Apophis with our planet in 229. To this aim, we irst compute the high order Taylor expansion o Apophis close encounter distance rom the Earth by means o map inversion and composition; then we run the proposed Monte Carlo algorithm to perorm the statistical analysis. Key words: Uncertainties Propagation; Monte Carlo Simulation; Apophis Close Encounter; Dierential Algebra. The propagation o uncertainties in orbital mechanics is usually addressed by linear propagation models [1, 7, 17] or nonlinear Monte Carlo simulations [14]. The main advantage o the linear methods is the simpliication o the problem, but their accuracy drops o or highly nonlinear systems and/or long time propagations. On the other hand, Monte Carlo simulations provide true trajectory statistics, but are computationally intensive. The tools currently used or the robust detection and prediction o planetary encounters and potential impacts o Near Earth Objects (NEO) are based on these techinques [8, 9, 15], and thus suer the same limitations. The eect o the coordinate system on the propagated statistics is analyzed by Junkins et al. [18, 19] and used to develop an alternative approach to orbit uncertainty propagation. However, this method is based on a linear assumption and thus cannot map nonlinearities. An alternate way to analyze trajectory statistics by incorporating higher-order Taylor series terms that describe localized nonlinear motion is proposed by Park and Scheeres [2]. Their approach is based on proving the integral invariance o the probability density unction via solutions o the Fokker Planck equations or diusionless systems, and by combining this result with the nonlinear state propagation to derive an analytic representation o the nonlinear uncertainty propagation. As a result, the method enables the nonlinear mapping o Gaussian statistics, bypassing the drawbacks o Monte Carlo simulations. However, it is limited to systems derived rom a single potential. 1. INTRODUCTION Dierential algebraic (DA) techniques are proposed as a valuable tool to develop an alternative approach to tackle the previous tasks. Dierential algebra supplies the tools to compute the derivatives o unctions within a computer environment [4, 5, 6]. More speciically, by substituting the classical implementation o real algebra with the implementation o a new algebra o Taylor polynomials, any unction o n variables is expanded into its Taylor polynomial up to an arbitrary order k. This has an important consequence when the numerical integra-
2 tion o an ordinary dierential equation (ODE) is perormed by means o an arbitrary integration scheme. Any explicit integration scheme is based on algebraic operations, involving the evaluation o the ODE right hand side at several integration points. Thereore, starting rom the DA representation o the initial condition and carrying out all the evaluations in the DA ramework, the low o an ODE is obtained at each step as its Taylor expansion in the initial condition [1]. The availability o such high order expansions is exploited when problems with uncertain initial conditions have to be analyzed. As the accuracy o the Taylor expansion can be kept arbitrarily high by adjusting the expansion order, the approach o classical Monte Carlo simulations can be enhanced by replacing thousands o integrations with evaluations o the Taylor expansion o the low. As a result, the computational time reduces considerably without any signiicant loss in accuracy. The algorithm is applied to the prediction o Apophis planetary encounter and potential impact, taking into account its measurement uncertainties. The availability o high order maps in space and time, and intrinsic tools or their inversion, are exploited to reduce the computation o the close encounter distance (CED) rom the Earth o all the asteroids belonging to the initial uncertainty cloud (commonly reerred to as virtual asteroids [16]) to the simple evaluation o polynomials. Similar techniques exploiting high order Taylor expansions o the low o ODE and their inverses obtained with DA techniques have already been eiciently utilized in beam physics. Two noticeable applications are the reconstruction o trajectories in particle spectrographs together with the reconstructive correction o residual aberrations [2], and the end-to-end simulations o ragment separators [11]. This paper presents an application to celestial mechanics. 2. DA INTEGRATION OF APOPHIS DY- NAMICS 2.1. Dynamical Models The study o the motion o a NEO in the Solar System with an accuracy suicient to predict collisions requires the inclusion o various relativistic corrections to the well-known Newtonian orces based on the Kepler s orce law. Speciically, the ull equation o motion in the Solar System including the relevant relativistic eects are considered. It is assumed that the object we are integrating is aected by the gravitational attraction o k bodies, but has no gravitational eect on them; i.e., we are adopting the restricted (n + 1)-body problem approximation. The positions, velocities, and accelerations o the k bodies are considered as given values, computed by cubic spline interpolations o data retrieved rom HORIZONS Web-Interace ( These interpolations are necessary as in the DA ramework the use o external code is not permitted. The cubic splines are built so as to keep the maximum error with respect to HORIZONS ephemerides o the order o 1 9 AU and 1 1 AU/day or bodies position and velocity, respectively (see [3] or details). In our integrations k includes the Sun, planets, the Moon, Ceres, Pallas, and Vesta. For planets with moons, with the exception o the Earth, the center o mass o the system is considered. The dynamical model is written in the J2. Ecliptic reerence rame and is commonly reerred to as Standard Dynamical Model [12]. To improve the integration accuracy the dynamics are scaled by Earth semimajor axis and Sun gravitational parameter (i.e., a E = 1 and µ S = Gm s = 1). We must mention that, to obtain a ull understanding o the dynamics o a body in the Solar System, other eects should be taken into account, such as: the orces due to other natural satellites and asteroids, the J 2 (and higher order harmonics o the potential) eect o the Earth and other bodies, Yarkosvsky and solar radiation pressure eects [12]. When the asteroid approaches the Earth, a dierent set o ODE are integrated to avoid cancellation errors associated to repetitive subtraction o Apophis and Earth s position vectors occurring across the lyby pericenter. In this case the equations o motion are written in the J2. Earth- Centered Inertial reerence rame. The same gravitational bodies o the heliocentric phase are considered, whereas relativistic corrections are neglected as their eect during a ast close encounter is negligible. In this phase the dynamics are scaled by the radius o the Earth and by the Earth gravitational parameter (i.e., R E = 1 and µ E = Gm E = 1).
3 2.2. Flow Expansion The high order expansion o the low o ODE can be straightorwardly obtained by evaluating any explicit numerical integration scheme within the DA ramework. The results presented here are obtained by applying a DA-based 8-th order Runge Kutta Fehlberg (RKF78) scheme with absolute and relative tolerance o The integration window is June 18, 29 to April 16, 229, being April 13, 229, the day in which the close approach occurs. The nominal initial state and the associated σ o Apophis, expressed in equinoctial variables p = (a, P 1, P 2, Q 1, Q 2, l), are taken rom the Near Earth Object Dynamic Site (newton.dm.unipi.it/neodys) and summarized in Table 1. Apophis initial condition is initialized y [AU] Earth Apophis x [AU] Figure 1. Apophis heliocentric phase trajectory. 2 x 1-3 Apophis Moon Table 1. Apophis equinoctial variable at 3456 MJD2 (June 18, 29) and associated σ values. Units: a [AU] and l [deg]. y [AU] 1 Nom Value a P P Q Q l as DA variables [p ] = p + 3σδp, where 3σ is used as a scaling actor. These variables are converted into cartesian coordinates using the relations given in Battin [1], evaluated in the DA ramework and then numerically propagated. Note that the solution o the Kepler equation, required or the computation o the eccentric longitude, is carried out by applying the DA-algorithm introduced in Di Lizia et al. [1]. The nominal heliocentric trajectories o Apophis and the Earth are shown in Fig. 1 by the solid and dotted lines, respectively. Figure 2 shows a zoom o Apophis close approach with the Earth in the geocentric reerence rame. It is worth mentioning that the maximum norm o the dierence between the computed trajectory and Apophis HORIZONS ephemerides σ x [AU] x 1-3 Figure 2. Apophis close encounter trajectory. is less than AU. The mismatch is due to all dierent initial conditions, dynamical and ephemeris model, and integration scheme. An analysis on the accuracy o the low expansion is mandatory beore introducing the DAbased Monte Carlo algorithm. Figures 3 and 4 show the maximum position and velocity error o the Taylor representation o the low at the corners o the initial set, with respect to the pointwise integration o the same points. Initial widths o 3, 6, and 9 σ and expansion orders rom 1 to 5 are considered. The expansion error decreases when higher expansion orders are selected and when smaller uncertain sets are considered. The errors tend to decrease exponentially with the expansion order, until reaching a lower limit o approximately [AU] on position and [AU/day] on velocity. It is worth notic-
4 Expansion error on position [AU] σ 6 σ 9 σ Expansion error on velocity [AU/day] σ 6 σ 9 σ Expansion order Figure 3. Accuracy o the Taylor expansion o the low corresponding to dierent expansion orders and initial uncertainty sets: position error Expansion order Figure 4. Accuracy o the Taylor expansion o the low corresponding to dierent expansion orders uncertainty sets: velocity error. ing that a ith order expansion guarantees a gain o approximately three order o magnitude in the low representation with respect to linear methods. This gain can be crucial when impact probability and/or resonant returns are studied. The igure clearly shows that Taylor polynomial accuracy is a unction o both the expansion order and domain width. The drawback or obtaining the Taylor expansion o the low with respect to the initial condition is the computational time to perorm a single integration, as shown in Fig. 5. In this igure the ratio between the computational time o a k-th order DA integration and a single pointwise integration is illustrated, underlining that a 5-th order integration is approximately eight times slower. On the other hand, the availability o the low expansion enables the development o a computationally eicient Monte Carlo method, as described in the next section. 3. DA-BASED MONTE CARLO Within the dynamical models adopted and the chosen integration scheme, the asteroid reerence solution has a close encounter distance rom the Earth center o mass o km at epoch MJD2. In order to evaluate the possibility o an Earth s impact it is necessary to accurately propagate the statistics o the asteroid. The accurate computation o statistics in nonlinear dynamical systems oten relies on Monte Carlo simulations. The algorithmic low o a Monte Carlo simulation is: Ratio o computational time Expansion order Figure 5. DA integration computational time compared to single pointwise integration. 1. Generate random samples based on the statistical distribution o the uncertainty to be propagated. 2. Run a pointwise integration o each sample in the ully nonlinear dynamics. 3. Perorm the statistical analysis o the results. There are three critical disadvantages when using this approach: convergence o the statistics usually requires a large number o sample trajectories to be propagated,
5 the simulation needs to be repeated or dierent initial distributions, it does not provide the user with analytical inormation, useul or additional analyses. These problems aect both the computational burden associated to the Monte Carlo simulation and its validity or dierent statistics [2]. The previous drawbacks become dramatic when thousands o long-term integrations are required, as or the analysis o possible NEO close encounter with the Earth [16]. DA integration delivers an arbitrary order Taylor expansion o the low o the ODE, which is analytic. Furthermore, the accuracy o the map expansion can be controlled by acting on the expansion order. For these reasons, it is possible to substitute the thousands o pointwise integrations required or classical Monte Carlo simulations with an equal number o map evaluations, i.e. ast polynomials evaluations. The resulting DA-based Monte Carlo simulation can be summarized as: 1. Perorm a single DA integration selecting the expansion order according to the demanded accuracy. 2. Generate random samples based on the statistical distribution o the uncertainty to be propagated. 3. Evaluate the low expansion map or all the samples, requiring only ast polynomial evaluations. 4. Perorm the statistical analysis o the results. The ratio between the computational time o a DA-based Monte Carlo simulation and its pointwise counterpart is given by t n + n s t e n s t, (1) where t n, t e, and t are the computational times o a k-th order DA integration, a low map evaluation, and a pointwise integration, respectively; and n s the number o samples o the Monte Carlo simulation. The computational cost o a Taylor map evaluation, depends on the expansion order, but it is negligible compared to a pointwise integration. For this reason, expression (1) can be approximated by m n s, in which m is the ratio between the computational time o a k-th order DA integration and a pointwise integration (see Fig. 5). The value o m strongly depends on the expansion order, but it is ew orders o magnitude smaller than the number o samples required or a good representation o the statistics. For this reason, the ratio m n s is small, proving that the proposed DA-based Monte Carlo simulation is computationally eicient. As an example, in Sec. 5, Fig. 1 will show that the computational time is reduced by a actor o at least 1 or a typical sample size o 1 virtual asteroids. In case new statistics need to be propagated, it is not necessary to perorm an additional DA integration as only steps 2 4 are required. Furthermore, i the statistical analysis is perormed or a dierent inal time, the possibility o obtaining Taylor expansions with respect to the inal time can be exploited (see Sec. 6.1). Moreover, as the low expansion is analytical, an analytic ramework is delivered. In conclusion, all the major drawbacks o a classical Monte Carlo approach are circumvented. These properties are better highlighted in Sect. 6 by applying the algorithm to the study o Apophis close encounter with the Earth in CED ALGORITHM Let us suppose the close approach o the nominal asteroid occurs at the epoch t, and consider the integration o the asteroid dynamics rom t = t to t = t. Initialize the initial state and the inal integration epoch as DA variables; i.e., [x ] = x + δx [t ] = t + δt, (2) where x is the initial condition corresponding to the nominal asteroid. Using the DA-based RKF78 integrator obtain the map [x ] = x + M x (δx, δt ). (3) The map (3) is the k-th order Taylor expansion o the low with respect to the initial condition
6 and the inal epoch about their nominal values x and t. Based on a mere DA-based computation, the inal solution x can be used to compute the Taylor expansion o distance rom the Earth [d ] = d + M d (δx, δt ). (4) More speciically, map (4) describes how the distance varies depending on the virtual asteroid and the inal integration epoch. Using the derivation operator available in the DA ramework, the Taylor expansion o the derivative d = d(d )/dt can be obtained [d ] = d + M d (δx, δt ). (5) The constant part o the map (5), d, is the derivative o the distance rom the Earth o the nominal solution at its close encounter; i.e., at CED epoch. Consequently, this is a stationary point or the nominal solution, and d =. Then, the map (5) reduces to δd = M d (δx, δt ), (6) in which we omit the bracket operator or the sake o a simpler notation. Consider the map ( ) ( ) ( ) δd Md δx =, (7) δx I x δt which is built by concatenating M d with the identity map or δx. Map (7) can now be inverted to obtain ( ) ( ) 1 ( ) δx Md δd =. (8) δt I x δx This is a ull nonlinear map inversion that is obtained by applying the algorithm illustrated in [5]. This algorithm reduces the map inversion problem to the solution o an equivalent ixed point problem, which can be solved with a ixed amount o eort in the DA setting. Map (4) is then concatenated to the identity map or δt to obtain ( ) ( ) ( ) d Md δx =. (9) δt I t δt Map (9) can now be composed with map (8) to obtain ( ) ( ) ( ) 1 ( ) d Md Md δd =, δt I t I x δx (1) which relates d and δt to the displacements o the derivative o the inal distance δd and o the state vector o the virtual asteroid δx i rom their values. As or the reerence value d =, the necessary condition or CED computation is Substituting into (1) yields ( d δt ) δd =. (11) ( ) ( Md Md = I t I x ) 1 ( δx ). (12) Eventually, map (12) delivers the desired explicit relation between the CED (d ) and the epoch at which it is reached (t + δt ) with the displacement δx i in terms o Taylor polynomials. Given any virtual asteroid belonging to the initial set (which corresponds to a speciic value o the displacement δx i ), the simple evaluation o the polynomials in (12) delivers the CED and the epoch at which it is reached. In this way, the problem highlighted by Milani et al. [15] is solved. 5. CED STATISTICAL ANALYSIS The DA-based Monte Carlo simulation introduced in Sec. 3 is run on map (12) to perorm the nonlinear mapping o the initial uncertainties on the CED. More speciically, 1 virtual asteroids are generated with a normal random distribution with mean value and standard deviation as in Table 1. For each sample, the displacement with respect to the nominal initial conditions is computed and map (12) is evaluated to obtain its CED and the associated epoch. The result is reported in Fig. 6 in terms o probability distribution or the CED. The analysis o the results shows that the mean value is km with a standard deviation o km, thus the possibility o having an Earth impact in 229 is ruled out. For the same virtual asteroids, map (12) is also evaluated to obtain the close encounter epochs. The result is presented in Fig. 7 in terms o the probability distribution o the displacement δt rom the nominal epoch t. The maximum displacement is o the order o 3 s. In Fig. 8 the
7 Probability % CED [km] x 1 4 Figure 6. Monte Carlo analysis o virtual asteroids CED. CED [km] 4.1 x δ t * [s] Figure 8. CED vs δt or 1 virtual asteroids. Probability % Distance rom the Earth [km] 4.1 x CED map evaluation Actual trajectory δ t * [s] Time [MJD2] x 1 4 Figure 7. Monte Carlo analysis o virtual asteroids CED epochs. Figure 9. Accuracy o the CED algorithm: virtual asteroids actual trajectories. 1 virtual asteroids are plotted in the CEDδt plane. For the sake o completeness, an accuracy analysis o the results is presented in Fig. 9. Ten virtual asteroids are randomly selected rom the initial set. For each virtual asteroid, the minimum distance and the corresponding epoch, resulting rom map (12), are reported in the igure. Then, a pointwise integration o the motion o each asteroid is perormed to obtain the proile o Earth s distances shown in the dotted lines. Although the accuracy on the identiication o the epoch o the close encounter is not clearly visualized, due to the very little displacement in t, it is clearly shown that the algorithm is able to accurately identiy the CED values o the resulting trajectories. Figure 1 concludes the analysis by showing the ratio o the computational time between the proposed DA-based Monte Carlo simulation and its pointwise counterpart as a unc- tion o the expansion order when 1 virtual asteroids are considered. It is apparent that the drawback o the higher computational cost required by a DA integration is rewarded by the signiicant time saving achieved by substituting 1 pointwise integrations with the same number o polynomials evaluations. 6. CONCLUSIONS The paper introduced a Monte Carlo simulation based on the high order Taylor expansion o the low o ODE, enabled by the use o dierential algebra. Being based on the replacement o pointwise integrations with ast evaluation o polynomials, the proposed algorithm guarantees signiicant computational time savings. The ac-
8 .2 ACKNOWLEDGMENTS Ratio o computational time % This work was inspired by [2] and was partially conducted under the ESA Ariadna study NEO Encounter 229. The authors are grateul to Dario Izzo or having encouraged the work and to Jon Giorgini or his precious comments Expansion order REFERENCES Figure 1. Percentage o computational time required by a DA-based Monte Carlo run versus a classical Monte Carlo simulation or 1 virtual asteroids. curacy o the algorithm can be suitably tuned by varying the low expansion order. Furthermore, the availability o analytic Taylor expansions and the use o DA embedded tools as map inversion, composition, and derivation allows the user to compute arbitrary order maps o the quantities on which the statistical analysis is perormed; thus, the algorithm is not limited to the low o ODE. More speciically, a technique or the automatic computation o both CED and CED epochs or all the virtual asteroids belonging to the initial uncertainty cloud has been developed. The eiciency and eectiveness o the methods are proven by applying them to the analysis o Apophis close encounter with the Earth occurring in April 229. In particular, it is shown that the nonlinear mapping o uncertainties can be perormed or any complex and arbitrary dynamics, even when long-term integrations are required; a ith order expansion increases the accuracy o the computation o the CED by approximately two orders o magnitude with respect to classical linear methods; the expansion in time allows or the proper identiication o the CED epoch or all the virtual asteroids. As an additional result, the occurrence o an impact with the Earth in April 229 can be ruled out. [1] Battin, R.H.: An Introduction to the Mathematics and Methods o Astrodynamics. AIAA Education Series, New York, [2] Berz, M., Joh, K., Nolen, J.A., Sherrill, B.M., and Zeller, A.F.: Reconstructive correction o aberration in nuclear particle spectrographs. Physical Review C, 47, pp , (1993). [3] Bernelli-Zazzera, F., Berz, M., Lavagna, M., Makino, K., Armellin, R., Di Lizia, P., Jagasia, R., and Topputo, F.: NEO Encounter 229: Orbital Prediction via Dierential Algebra and Taylor Models. Final Report, Ariadna id: 8-433, Contract No. 2271/6/NL/HI (29). [4] Berz, M.: Dierential Algebraic Techniques. Entry in Handbook o Accelerator Physics and Engineering, World Scientiic, New York, (1999a). [5] Berz, M.: Modern Map Methods in Particle Beam Physics. Academic Press, (1999b). [6] Berz, M., and Makino, K.: COSY INFIN- ITY version 9 reerence manual. MSU Report MSUHEP-683, Michigan State University, East Lansing, MI 48824, pp 1-84, (26). [7] Crassidis, J.L., and Junkins, J.L.: Optimal Estimation o Dynamics Systems. CRC Press LLC, Boca Raton, FL, pp , (24). [8] Chesley, S.R., and Milani, A.: An automatic Earth-asteroid collision monitoring system. Bull. Am. Astron. Soc., 32, pp 682, (2). [9] Chodas, P.W., and Yeomans, D.K.: Predicting close approaches and estimating impact probabilities or near-earth projects. AAS/AIAA Astrodynamics Specialists Conerence, Girdwood, Alaska, (1999).
9 [1] Di Lizia, P.: Robust Space Trajectory and Space System Design using Dierential Algebra. Ph.D. dissertation, Politecnico di Milano (28). [11] Erdelyi, B., Bandura, L., Nolen, J., and Manikonda, S.: Code development or Next- Generation High-Intensity Large Acceptance Fragment Separators, Proceedings o PAC7, Albuquerque, New Mexico, USA, (27). [12] Giorgini, J.D., Benner, L.A.M., Ostro, S.J., Nolan, M.C., and Busch, M.W.: Predicting the Earth encounters o (99942) Apophis, Icarus 193, pp 1 19, (28). [13] Hoekens, J., Berz, M., and Makino, K.: Controlling the Wrapping Eect in the Solution o ODEs or Asteroids, Reliable Computing, 8, pp 21-41, (23). [14] Maybeck, P.S., Stochastic Models, Estimation, and Control. Academic Press, New York, 2, pp , (1982). [15] Milani, A., Chesley, S.R., and Valsecchi, G.B.: Asteroid close encounters with Earth: Risk assessment. Planet. Space Sci., 48, pp , (2). [16] Milani, A. and Chesley, S.R. and Chodas, P.W. and Valsecchi, G.B.: Asteroid close approaches: analysis and potential impact detection. Asteroids III, pp 89 11, (22). [17] Montenbruck, O., and Gill, E.: Satellite Orbits. 2nd ed., SpringerVerlag, New York, pp , (21). [18] Junkins, J., Akella, M., and Alriend, K.: Non-Gaussian Error Propagation in Orbit Mechanics. Journal o the Astronautical Sciences, 44, pp , (1996). [19] Junkins, J., and Singla, P.: How Nonlinear Is It? A Tutorial on Nonlinearity o Orbit and Attitude Dynamics. Journal o the Astronautical Sciences, 52, pp 7 6 (24). [2] Park, R., and Scheeres, D.: Nonlinear Mapping o Gaussian Statistics: Theory and Applications to Spacecrat Trajectory Design. Journal o guidance, Control, and Dynamics, 29, pp , (26). [21] Park, R.S. and Scheeres, D.J.: Nonlinear Semi-Analytic Methods or Trajectory Estimation. Journal o Guidance, Control and Dynamics, 3, , (27). [22] Seidelmann, P.K.: Explanatory Supplement to the Astronomical Almanac. University Science Books, Mill Valley, Caliornia, (1992).
Apophis Encounter 2029: Differential Algebra and Taylor Model Approaches. 1 st IAA Planetary Defense Conference: Protecting Earth from Asteroids
Apophis Encounter 9: Differential Algebra and Taylor Model Approaches 1 st IAA Planetary Defense Conference: Protecting Earth from Asteroids 7- April 9 Granada, Spain R. Armellin (1), P. Di Lizia (1),
More informationValidated integration of Solar System Dynamics. 1 st IAA Planetary Defense Conference: Protecting Earth from Asteroids
Validated integration of Solar System Dynamics 1 st IAA Planetary Defense Conference: Protecting Earth from Asteroids 27-3 April 29 Granada, Spain P. Di Lizia (1), R. Armellin (1), F. Bernelli Zazzera
More informationVerified High-Order Optimal Control in Space Flight Dynamics
Verified High-Order Optimal Control in Space Flight Dynamics R. Armellin, P. Di Lizia, F. Bernelli-Zazzera K. Makino and M. Berz Fourth International Workshop on Taylor Methods Boca Raton, December 16
More informationRelevant Applications of Differential Algebra in Astrodynamics
AstroNet-II International Final Conference Tossa de Mar 15-19 June 2015 Relevant Applications of Differential Algebra in Astrodynamics Pierluigi Di Lizia in collaboration with: Roberto Armellin, Alessandro
More informationFrom Particle Accelerators to Celestial Dynamics: An Introduction to Differential Algebra
AstroNet II Summer School 2014 Friday, 27 June 2014 Zielona Góra (Poland) From Particle Accelerators to Celestial Dynamics: An Introduction to Differential Algebra 1, Pierluigi di Lizia 1, Roberto Armellin
More informationTwo-Point Boundary Value Problem and Optimal Feedback Control based on Differential Algebra
Two-Point Boundary Value Problem and Optimal Feedback Control based on Differential Algebra Politecnico di Milano Department of Aerospace Engineering Milan, Italy Taylor Methods and Computer Assisted Proofs
More informationTrajectory Optimisation under Uncertainties
Trajectory Optimisation under Uncertainties Final Report Authors: F. Bernelli-Zazzera, M. Lavagna, R. Armellin, P. Di Lizia, A. Morselli Affiliation: Department of Aerospace Engineering, Politecnico di
More informationGlobal Optimization of Impulsive Interplanetary Transfers
Global Optimization of Impulsive Interplanetary Transfers R. Armellin, Dipartimento di Ingegneria Aerospaziale, Politecnico di Milano Taylor Methods and Computer Assisted Proofs Barcelona, June, 3 7, 2008
More informationRigorous Global Optimization of Impulsive Space Trajectories
Rigorous Global Optimization of Impulsive Space Trajectories P. Di Lizia, R. Armellin, M. Lavagna K. Makino, M. Berz Fourth International Workshop on Taylor Methods Boca Raton, December 16 19, 2006 Motivation
More informationCOSY INFINITY version 8
Nuclear Instruments and Methods in Physics Research A 427 (1999) 338}343 COSY INFINITY version 8 Kyoko Makino*, Martin Berz Department of Physics and National Superconducting Cyclotron Laboratory, Michigan
More informationCOSY INFINITY Version 7
COSY INFINITY Version 7 Kyoko Makino and Martin Berz Department of Physics and Astronomy and National Superconducting Cyclotron Laboratory Michigan State University, East Lansing, MI 48824 Abstract An
More informationHigh-Order Representation of Poincaré Maps
High-Order Representation of Poincaré Maps Johannes Grote, Martin Berz, and Kyoko Makino Department of Physics and Astronomy, Michigan State University, East Lansing, MI, USA. {grotejoh,berz,makino}@msu.edu
More informationGlobal Trajectory Optimisation: Can We Prune the Solution Space when Considering Deep Space Maneuvers?
Global Trajectory Optimisation: Can We Prune the Solution Space when Considering Deep Space Maneuvers? Final Report Authors: F. Bernelli-Zazzera, M. Lavagna, R. Armellin, P. Di Lizia, F. Topputo Affiliation:
More informationReceived: 30 July 2017; Accepted: 29 September 2017; Published: 8 October 2017
mathematics Article Least-Squares Solution o Linear Dierential Equations Daniele Mortari ID Aerospace Engineering, Texas A&M University, College Station, TX 77843, USA; mortari@tamu.edu; Tel.: +1-979-845-734
More informationAccurate numerical orbit propagation using Polynomial Algebra Computational Engine PACE. ISSFD 2015 congress, Germany. Dated: September 14, 2015
Accurate numerical orbit propagation using Polynomial Algebra Computational Engine PACE Emmanuel Bignon (), Pierre mercier (), Vincent Azzopardi (), Romain Pinède () ISSFD 205 congress, Germany Dated:
More informationPRECISION ZEM/ZEV FEEDBACK GUIDANCE ALGORITHM UTILIZING VINTI S ANALYTIC SOLUTION OF PERTURBED KEPLER PROBLEM
AAS 16-345 PRECISION ZEM/ZEV FEEDBACK GUIDANCE ALGORITHM UTILIZING VINTI S ANALYTIC SOLUTION OF PERTURBED KEPLER PROBLEM Jaemyung Ahn, * Yanning Guo, and Bong Wie A new implementation o a zero-eort-miss/zero-eort-velocity
More informationThe achievable limits of operational modal analysis. * Siu-Kui Au 1)
The achievable limits o operational modal analysis * Siu-Kui Au 1) 1) Center or Engineering Dynamics and Institute or Risk and Uncertainty, University o Liverpool, Liverpool L69 3GH, United Kingdom 1)
More informationThird Body Perturbation
Third Body Perturbation p. 1/30 Third Body Perturbation Modeling the Space Environment Manuel Ruiz Delgado European Masters in Aeronautics and Space E.T.S.I. Aeronáuticos Universidad Politécnica de Madrid
More informationINVARIANT RELATIVE SATELLITE MOTION. Marco Sabatini Dip. Ingegneria Aerospaziale, Università di Roma La Sapienza
AC Workshop on Innovative Concepts. ESA-ESEC 8-9 January 8. INVARIAN RELAIVE SAELLIE MOION Marco Sabatini Dip. Ingegneria Aerospaziale, Università di Roma La Sapienza marcosabatini@hotmail.it Dario Izzo
More informationA Study of Six Near-Earth Asteroids Summary Introduction
A Study of Six Near-Earth Asteroids John Junkins 1, Puneet Singla 1, Daniele Mortari 1, William Bottke 2, Daniel Durda 2 Summary We consider here 6 Earth-orbit-crossing asteroids (ECAs) as possible targets
More informationAn Introduction to Differential Algebra
An Introduction to Differential Algebra Alexander Wittig1, P. Di Lizia, R. Armellin, et al. 1 ESA Advanced Concepts Team (TEC-SF) SRL, Milan Dinamica Outline 1 Overview Five Views of Differential Algebra
More informationAPPENDIX 1 ERROR ESTIMATION
1 APPENDIX 1 ERROR ESTIMATION Measurements are always subject to some uncertainties no matter how modern and expensive equipment is used or how careully the measurements are perormed These uncertainties
More informationImpact Risk Assessment of a Fragmented Asteroid in Earth Resonant Orbits
Impact Risk Assessment of a Fragmented Asteroid in Earth Resonant Orbits George Vardaxis and Bong Wie Iowa State University, Ames, IA 50011-2271, USA Taking advantage of direct numerical simulation, analytic
More informationUncertainties Management for Orbit Propagation using Taylor Differential Algebra
Uncertainties Management for Orbit Propagation using aylor Differential Algebra Vincent Azzopardi hales ervices, pace Flight Dynamics department, oulouse, France Institut upérieur de l Aéronautique et
More informationIMPACT RISK ANALYSIS OF NEAR-EARTH ASTEROIDS WITH MULTIPLE SUCCESSIVE EARTH ENCOUNTERS
(Preprint) AAS 16-267 IMPACT RISK ANALYSIS OF NEAR-EARTH ASTEROIDS WITH MULTIPLE SUCCESSIVE EARTH ENCOUNTERS George Vardaxis and Bong Wie INTRODUCTION Accurate estimation of the impact risk associated
More informationCOVARIANCE DETERMINATION, PROPAGATION AND INTERPOLATION TECHNIQUES FOR SPACE SURVEILLANCE. European Space Surveillance Conference 7-9 June 2011
COVARIANCE DETERMINATION, PROPAGATION AND INTERPOLATION TECHNIQUES FOR SPACE SURVEILLANCE European Space Surveillance Conference 7-9 June 2011 Pablo García (1), Diego Escobar (1), Alberto Águeda (1), Francisco
More informationNumerical analysis of the compliance of interplanetary CubeSats with planetary protection requirements
Numerical analysis of the compliance of interplanetary CubeSats with planetary protection requirements Francesca Letizia, Camilla Colombo University of Southampton 5th Interplanetary CubeSat Workshop Session
More informationProbabilistic Optimisation applied to Spacecraft Rendezvous on Keplerian Orbits
Probabilistic Optimisation applied to pacecrat Rendezvous on Keplerian Orbits Grégory aive a, Massimiliano Vasile b a Université de Liège, Faculté des ciences Appliquées, Belgium b Dipartimento di Ingegneria
More informationProper initial conditions for long-term integrations of the solar system
Proper initial conditions for long-term integrations of the solar system M. Arminjon Laboratoire Sols, Solides, Structures [CNRS / Université Joseph Fourier / Institut National Polytechnique de Grenoble]
More informationFeedback Optimal Control of Low-thrust Orbit Transfer in Central Gravity Field
Vol. 4, No. 4, 23 Feedback Optimal Control of Low-thrust Orbit Transfer in Central Gravity Field Ashraf H. Owis Department of Astronomy, Space and Meteorology, Faculty of Science, Cairo University Department
More informationEFFECTS OF THE WIGGLER ON THE HEFEI LIGHT SOURCE STORAGE RING
International Journal of Modern Physics A Vol. 24, No. 5 (2009) 1057 1067 c World Scientific Publishing Company EFFECTS OF THE WIGGLER ON THE HEFEI LIGHT SOURCE STORAGE RING HE ZHANG and MARTIN BERZ Department
More informationINTERPLANETARY AND LUNAR TRANSFERS USING LIBRATION POINTS
INTERPLANETARY AND LUNAR TRANSFERS USING LIBRATION POINTS Francesco Topputo (), Massimiliano Vasile () and Franco Bernelli-Zazzera () () Dipartimento di Ingegneria Aerospaziale, Politecnico di Milano,
More informationDIFFERENTIAL ALGEBRA SPACE TOOLBOX FOR NONLINEAR UNCERTAINTY PROPAGATION IN SPACE DYNAMICS
DIFFERENTIAL ALGEBRA SPACE TOOLBOX FOR NONLINEAR UNCERTAINTY PROPAGATION IN SPACE DYNAMICS M. Rasotto, A. Morselli, A. Wittig, M. Massari, P. Di Lizia, R. Armellin, C. Y. Valles and G. Ortega Dinamica
More informationA fast method for estimation of the impact probability of near-earth objects
doi:101093/mnras/stu2097 A fast method for estimation of the impact probability of near-earth objects D E Vavilov and Yu D Medvedev Institute of Applied Astronomy of Russian Academy of Science, StPetersburg,
More informationThe Ascent Trajectory Optimization of Two-Stage-To-Orbit Aerospace Plane Based on Pseudospectral Method
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 00 (014) 000 000 www.elsevier.com/locate/procedia APISAT014, 014 Asia-Paciic International Symposium on Aerospace Technology,
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 5. Numerical Methods Gaëtan Kerschen Space Structures & Systems Lab (S3L) Why Different Propagators? Analytic propagation: Better understanding of the perturbing forces. Useful
More informationA posteriori reading of Virtual Impactors impact probability
A posteriori reading of Virtual Impactors impact probability Germano D Abramo Istituto di Astrofisica Spaziale e Fisica Cosmica, Area di Ricerca CNR Tor Vergata, Roma, Italy E mail: dabramo@rm.iasf.cnr.it
More informationaddition, multiplication, and dierentiation on the space of innitely often dierentiable functions to a suitable smaller space that can be described on
Modern Map Methods for Charged Particle Optics Martin Berz Department of Physics and National Superconducting Cyclotron Laboratory Michigan State University East Lansing, MI 48824 Published in \Nuclear
More informationOn High-Rate Cryptographic Compression Functions
On High-Rate Cryptographic Compression Functions Richard Ostertág and Martin Stanek Department o Computer Science Faculty o Mathematics, Physics and Inormatics Comenius University Mlynská dolina, 842 48
More informationROBUST STABILITY AND PERFORMANCE ANALYSIS OF UNSTABLE PROCESS WITH DEAD TIME USING Mu SYNTHESIS
ROBUST STABILITY AND PERFORMANCE ANALYSIS OF UNSTABLE PROCESS WITH DEAD TIME USING Mu SYNTHESIS I. Thirunavukkarasu 1, V. I. George 1, G. Saravana Kumar 1 and A. Ramakalyan 2 1 Department o Instrumentation
More informationAUTOMATIC MINIMUM PRINCIPLE FORMULATION FOR LOW THRUST OPTIMAL CONTROL IN ORBIT TRANSFERS USING COMPLEX NUMBERS
AUOMAIC MINIMUM PRINCIPLE FORMULAION FOR LOW HRUS OPIMAL CONROL IN ORBI RANSFERS USING COMPLEX NUMBERS hierry Dargent hales Alenia Space, 1 Bd du Midi, BP 99, 6322 Cannes La Bocca Cedex France, E-mail:
More informationA study upon Eris. I. Describing and characterizing the orbit of Eris around the Sun. I. Breda 1
Astronomy & Astrophysics manuscript no. Eris c ESO 2013 March 27, 2013 A study upon Eris I. Describing and characterizing the orbit of Eris around the Sun I. Breda 1 Faculty of Sciences (FCUP), University
More informationEarth Risk Corridor Computations for 2011 AG 5 on 5 February Introduction
Introduction The near-earth object (NEO) designated 2011 AG 5 is currently associated with a relatively high linearized Earth impact probability of 0.001557 (odds of 1-in-642.3) during its predicted 5
More informationAPPLYING FLYWHEEL ENERGY STORAGE TO SOLAR ELECTRIC ORBITAL TRANSFERS THESIS. Presented to the Faculty of the Graduate School of Engineering
AFIT/GA/ENY/97D-02 APPLYING FLYWHEEL ENERGY STORAGE TO SOLAR ELECTRIC ORBITAL TRANSFERS THESIS Presented to the Faculty o the Graduate School o Engineering o the Air Force Institute o Technology Air University
More informationScattered Data Approximation of Noisy Data via Iterated Moving Least Squares
Scattered Data Approximation o Noisy Data via Iterated Moving Least Squares Gregory E. Fasshauer and Jack G. Zhang Abstract. In this paper we ocus on two methods or multivariate approximation problems
More informationLeast-Squares Spectral Analysis Theory Summary
Least-Squares Spectral Analysis Theory Summary Reerence: Mtamakaya, J. D. (2012). Assessment o Atmospheric Pressure Loading on the International GNSS REPRO1 Solutions Periodic Signatures. Ph.D. dissertation,
More informationBayesian Technique for Reducing Uncertainty in Fatigue Failure Model
9IDM- Bayesian Technique or Reducing Uncertainty in Fatigue Failure Model Sriram Pattabhiraman and Nam H. Kim University o Florida, Gainesville, FL, 36 Copyright 8 SAE International ABSTRACT In this paper,
More informationScattering of Solitons of Modified KdV Equation with Self-consistent Sources
Commun. Theor. Phys. Beijing, China 49 8 pp. 89 84 c Chinese Physical Society Vol. 49, No. 4, April 5, 8 Scattering o Solitons o Modiied KdV Equation with Sel-consistent Sources ZHANG Da-Jun and WU Hua
More informationWhen Does the Uncertainty Become Non-Gaussian. Kyle T. Alfriend 1 Texas A&M University Inkwan Park 2 Texas A&M University
When Does the Uncertainty Become Non-Gaussian Kyle T. Alfriend Texas A&M University Inkwan Park 2 Texas A&M University ABSTRACT The orbit state covariance is used in the conjunction assessment/probability
More informationA high order method for orbital conjunctions analysis: Monte Carlo collision probability computation
A high order method for orbital conjunctions analysis: Monte Carlo collision probability computation Alessandro Morselli 1, Dipartimento di Scienze e Tecnologie Aerospaziali, Politecnico di Milano, Via
More informationEnd-Of-Life Disposal Concepts for Lagrange-Point and Highly Elliptical Orbit Missions
End-Of-Life Disposal Concepts for Lagrange-Point and Highly Elliptical Orbit Missions Executive summary of the main study and the study extension Version 1.0 12 February 2015 ESA/ESOC contract No. 4000107624/13/F/MOS
More informationCALCULATION OF POSITION AND VELOCITY OF GLONASS SATELLITE BASED ON ANALYTICAL THEORY OF MOTION
ARTIFICIAL SATELLITES, Vol. 50, No. 3 2015 DOI: 10.1515/arsa-2015-0008 CALCULATION OF POSITION AND VELOCITY OF GLONASS SATELLITE BASED ON ANALYTICAL THEORY OF MOTION W. Góral, B. Skorupa AGH University
More informationBézier Description of Space Trajectories
Bézier Description of Space Trajectories Francesco de Dilectis, Daniele Mortari, Texas A&M University, College Station, Texas and Renato Zanetti NASA Jonhson Space Center, Houston, Texas I. Introduction
More informationimin...
Pulsar Timing For a detailed look at pulsar timing and other pulsar observing techniques, see the Handbook of Pulsar Astronomy by Duncan Lorimer and Michael Kramer. Pulsars are intrinsically interesting
More informationPubl. Astron. Obs. Belgrade No. 90 (2010), DYNAMICAL CHARACTERISTICS OF HUNGARIA ASTEROIDS 1. INTRODUCTION
Publ. Astron. Obs. Belgrade No. 9 (21), 11-18 Invited Lecture DYNAMICAL CHARACTERISTICS OF HUNGARIA ASTEROIDS Z. KNEŽEVIĆ1, B. NOVAKOVIĆ2, A. MILANI 3 1 Astronomical Observatory, Volgina 7, 116 Belgrade,
More informationAn Analysis of N-Body Trajectory Propagation. Senior Project. In Partial Fulfillment. of the Requirements for the Degree
An Analysis of N-Body Trajectory Propagation Senior Project In Partial Fulfillment of the Requirements for the Degree Bachelor of Science in Aerospace Engineering by Emerson Frees June, 2011 An Analysis
More informationCelestial Mechanics I. Introduction Kepler s Laws
Celestial Mechanics I Introduction Kepler s Laws Goals of the Course The student will be able to provide a detailed account of fundamental celestial mechanics The student will learn to perform detailed
More informationChapter 6 Reliability-based design and code developments
Chapter 6 Reliability-based design and code developments 6. General Reliability technology has become a powerul tool or the design engineer and is widely employed in practice. Structural reliability analysis
More informationExtending the Patched-Conic Approximation to the Restricted Four-Body Problem
Monografías de la Real Academia de Ciencias de Zaragoza 3, 133 146, (6). Extending the Patched-Conic Approximation to the Restricted Four-Body Problem Thomas R. Reppert Department of Aerospace and Ocean
More informationSteered Spacecraft Deployment Using Interspacecraft Coulomb Forces
Steered Spacecrat Deployment Using Interspacecrat Coulomb Forces Gordon G. Parker, Lyon B. King and Hanspeter Schaub Abstract Recent work has shown that Coulomb orces can be used to maintain ixed-shape
More informationNEOPROP: A NEO PROPAGATOR FOR SPACE SITUATIONAL AWARENESS FOR THE 6 TH IAASS CONFERENCE
NEOPROP: A NEO PROPAGATOR FOR SPACE SITUATIONAL AWARENESS FOR THE 6 TH IAASS CONFERENCE Valentino Zuccarelli (1), David Bancelin (2),Sven Weikert (3), William Thuillot (4), Daniel Hestroffer (5), Celia
More informationExperimental Analysis of Low Earth Orbit Satellites due to Atmospheric Perturbations
Experimental Analysis of Low Earth Orbit Satellites due to Atmospheric Perturbations Aman Saluja #1, Manish Bansal #2, M Raja #3, Mohd Maaz #4 #Aerospace Department, University of Petroleum and Energy
More informationOBSERVER/KALMAN AND SUBSPACE IDENTIFICATION OF THE UBC BENCHMARK STRUCTURAL MODEL
OBSERVER/KALMAN AND SUBSPACE IDENTIFICATION OF THE UBC BENCHMARK STRUCTURAL MODEL Dionisio Bernal, Burcu Gunes Associate Proessor, Graduate Student Department o Civil and Environmental Engineering, 7 Snell
More informationAnalytical Non-Linear Uncertainty Propagation: Theory And Implementation
Analytical Non-Linear Uncertainty Propagation: Theory And Implementation Kohei Fujimoto 1 Texas A&M University, College Station, TX, 77843 Daniel J. Scheeres 2 The University of Colorado-Boulder, Boulder,
More informationorbits Moon, Planets Spacecrafts Calculating the and by Dr. Shiu-Sing TONG
A Science Enrichment Programme for Secondary 3-4 Students : Teaching and Learning Resources the and Spacecrafts orbits Moon, Planets Calculating the 171 of by Dr. Shiu-Sing TONG 172 Calculating the orbits
More informationTheory and equations for Craters from Impacts and Explosions. Keith A. Holsapple
Theory and equations or Craters rom Impacts and Explosions Keith A. Holsapple 1. Dimensionless Forms As the primary example, the volume V o a crater ormed by a given impact can be expected to depend on
More informationSolar System dynamics: recent developments
Solar System dynamics: recent developments G. B. Valsecchi, IASF-Roma, INAF, Roma (Italy) Recent developments: planets Dynamics of the planets Slow chaos in planetary motions. Chaotic spins of the inner
More informationMISS DISTANCE GENERALIZED VARIANCE NON-CENTRAL CHI DISTRIBUTION. Ken Chan ABSTRACT
MISS DISTANCE GENERALIZED VARIANCE NON-CENTRAL CI DISTRIBUTION en Chan E-Mail: ChanAerospace@verizon.net ABSTRACT In many current practical applications, the probability o collision is oten considered
More informationCelestial Mechanics III. Time and reference frames Orbital elements Calculation of ephemerides Orbit determination
Celestial Mechanics III Time and reference frames Orbital elements Calculation of ephemerides Orbit determination Orbital position versus time: The choice of units Gravitational constant: SI units ([m],[kg],[s])
More informationRoberto s Notes on Differential Calculus Chapter 8: Graphical analysis Section 1. Extreme points
Roberto s Notes on Dierential Calculus Chapter 8: Graphical analysis Section 1 Extreme points What you need to know already: How to solve basic algebraic and trigonometric equations. All basic techniques
More informationBenefits of a Geosynchronous Orbit (GEO) Observation Point for Orbit Determination
Benefits of a Geosynchronous Orbit (GEO) Observation Point for Orbit Determination Ray Byrne, Michael Griesmeyer, Ron Schmidt, Jeff Shaddix, and Dave Bodette Sandia National Laboratories ABSTRACT Determining
More informationState Vector Reference Body Transformations Throughout The Solar System
Introduction One of the more mind-boggling programming challenges encountered during interplanetary trajectory design software development is how to transform a specified inertial position and velocity
More informationGEO Collisional Risk Assessment Based on Analysis of NASA-WISE Data and Modeling
GEO Collisional Risk Assessment Based on Analysis o NASA-WISE Data and Modeling Jeremy Murray Krezan, Samantha Howard, Phan D. Dao, Derek Surka AFRL Space Vehicles Directorate, Applied Technoy Associates
More informationRobust Residual Selection for Fault Detection
Robust Residual Selection or Fault Detection Hamed Khorasgani*, Daniel E Jung**, Gautam Biswas*, Erik Frisk**, and Mattias Krysander** Abstract A number o residual generation methods have been developed
More informationAIM RS: Radio Science Investigation with AIM
Prepared by: University of Bologna Ref. number: ALMARS012016 Version: 1.0 Date: 08/03/2017 PROPOSAL TO ESA FOR AIM RS Radio Science Investigation with AIM ITT Reference: Partners: Radio Science and Planetary
More informationRESOLUTION MSC.362(92) (Adopted on 14 June 2013) REVISED RECOMMENDATION ON A STANDARD METHOD FOR EVALUATING CROSS-FLOODING ARRANGEMENTS
(Adopted on 4 June 203) (Adopted on 4 June 203) ANNEX 8 (Adopted on 4 June 203) MSC 92/26/Add. Annex 8, page THE MARITIME SAFETY COMMITTEE, RECALLING Article 28(b) o the Convention on the International
More informationAn Alternative Poincaré Section for Steady-State Responses and Bifurcations of a Duffing-Van der Pol Oscillator
An Alternative Poincaré Section or Steady-State Responses and Biurcations o a Duing-Van der Pol Oscillator Jang-Der Jeng, Yuan Kang *, Yeon-Pun Chang Department o Mechanical Engineering, National United
More informationRELATIVISTIC EFFECTS O N SV CLOCKS DUE TO ORBIT CHANGES, AND DUE TO EARTH S OBLATENESS
33rdAnnual Precise Time and Time Interval (PTTI)Meeting RELATIVISTIC EFFECTS O N SV CLOCKS DUE TO ORBIT CHANGES, AND DUE TO EARTH S OBLATENESS Neil Ashby Department o Physics, UCB 390 University o Colorado,
More informationVerified Global Optimization with Taylor Model based Range Bounders
Verified Global Optimization with Taylor Model based Range Bounders KYOKO MAKINO and MARTIN BERZ Department of Physics and Astronomy Michigan State University East Lansing, MI 88 USA makino@msu.edu, berz@msu.edu
More informationGuidance, Navigation, and Control Conference and Exhibit August 5 8, 2002/Monterey, CA
AIAA Guidance, Navigation, and Control Conerence and Exhibit 5-8 August 22, Monterey, Caliornia AIAA 22-53 AIAA 22-53 Tracking and pointing o target by a Biocal Relay Mirror Spacecrat using attitude control
More informationPDC2015 Frascati, Roma, Italy IAA-PDC-15-P-84 GLOBAL IMPACT DISTRIBUTION OF ASTEROIDS AND AFFECTED POPULATION. Kingdom,
PDC2015 Frascati, Roma, Italy Planetary Defense Recent Progress & Plans NEO Discovery NEO Characterization Mitigation Techniques & Missions Impact Effects that Inform Warning, Mitigation & Costs Consequence
More informationSupplementary material for Continuous-action planning for discounted infinite-horizon nonlinear optimal control with Lipschitz values
Supplementary material or Continuous-action planning or discounted ininite-horizon nonlinear optimal control with Lipschitz values List o main notations x, X, u, U state, state space, action, action space,
More informationEstimation of Sample Reactivity Worth with Differential Operator Sampling Method
Progress in NUCLEAR SCIENCE and TECHNOLOGY, Vol. 2, pp.842-850 (2011) ARTICLE Estimation o Sample Reactivity Worth with Dierential Operator Sampling Method Yasunobu NAGAYA and Takamasa MORI Japan Atomic
More informationModel for the spacetime evolution of 200A GeV Au-Au collisions
PHYSICAL REVIEW C 70, 021903(R) (2004) Model or the spacetime evolution o 200A GeV Au-Au collisions Thorsten Renk Department o Physics, Duke University, P.O. Box 90305, Durham, North Carolina 27708, USA
More informationSimulation and Optimization of the Tevatron Accelerator
Simulation and Optimization of the Tevatron Accelerator Pavel Snopok, Carol Johnstone, and Martin Berz 3 Fermi National Accelerator Laboratory, Michigan State University, snopok@pa.msu.edu Fermi National
More informationarxiv: v1 [astro-ph.ep] 21 Jan 2019
On the predictability of Galileo disposal orbits David J. Gondelach Roberto Armellin Alexander Wittig arxiv:1901.06947v1 [astro-ph.ep] 21 Jan 2019 Abstract The end-of-life disposal of Galileo satellites
More informationFeasibility of a Multi-Pass Thomson Scattering System with Confocal Spherical Mirrors
Plasma and Fusion Research: Letters Volume 5, 044 200) Feasibility o a Multi-Pass Thomson Scattering System with Conocal Spherical Mirrors Junichi HIRATSUKA, Akira EJIRI, Yuichi TAKASE and Takashi YAMAGUCHI
More informationTwo-step self-tuning phase-shifting interferometry
Two-step sel-tuning phase-shiting intererometry J. Vargas, 1,* J. Antonio Quiroga, T. Belenguer, 1 M. Servín, 3 J. C. Estrada 3 1 Laboratorio de Instrumentación Espacial, Instituto Nacional de Técnica
More informationIntroduction to Simulation - Lecture 2. Equation Formulation Methods. Jacob White. Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy
Introduction to Simulation - Lecture Equation Formulation Methods Jacob White Thanks to Deepak Ramaswamy, Michal Rewienski, and Karen Veroy Outline Formulating Equations rom Schematics Struts and Joints
More informationProbabilistic Model of Error in Fixed-Point Arithmetic Gaussian Pyramid
Probabilistic Model o Error in Fixed-Point Arithmetic Gaussian Pyramid Antoine Méler John A. Ruiz-Hernandez James L. Crowley INRIA Grenoble - Rhône-Alpes 655 avenue de l Europe 38 334 Saint Ismier Cedex
More informationNew method for two-point nonuniformity correction of microbolometer detectors
10 th International Conerence on Quantitative InraRed Thermography July 27-30, 2010, Québec (Canada) New method or two-point nonuniormity correction o microbolometer detectors by R. Olbrycht*, B. Wiecek*,
More information2 INTRODUCTION 3 ORBIT DETERMINATION
Investigation of Properties and Characteristics of High-Area-to-Mass Ratio Objects Based on Examples of Optical Observation Data of Space Debris Objects in GEO-like Orbits Carolin Früh Astronomical Institute,
More informationGENERALIZATIONS AND APPLICATIONS OF THE LAGRANGE IMPLICIT FUNCTION THEOREM
AAS 8-3 GENERALIZATIONS AND APPLICATIONS OF THE LAGRANGE IMPLICIT FUNCTION THEOREM John L. Junkins *, James D. Turner, Manoranjan Majji Dedicated to F. Landis Markley INTRODUCTION The Implicit unction
More information2.6 Two-dimensional continuous interpolation 3: Kriging - introduction to geostatistics. References - geostatistics. References geostatistics (cntd.
.6 Two-dimensional continuous interpolation 3: Kriging - introduction to geostatistics Spline interpolation was originally developed or image processing. In GIS, it is mainly used in visualization o spatial
More informationAnalytical Non-Linear Uncertainty Propagation: Theory And Implementation
Analytical Non-Linear Uncertainty Propagation: Theory And Implementation Kohei Fujimoto 1), Daniel J. Scheeres 2), and K. Terry Alfriend 1) May 13, 2013 1) Texas A&M University 2) The University of Colorado
More informationLong-Term Evolution of High Earth Orbits: Effects of Direct Solar Radiation Pressure and Comparison of Trajectory Propagators
Long-Term Evolution of High Earth Orbits: Effects of Direct Solar Radiation Pressure and Comparison of Trajectory Propagators by L. Anselmo and C. Pardini (Luciano.Anselmo@isti.cnr.it & Carmen.Pardini@isti.cnr.it)
More informationRADIALLY ADAPTIVE EVALUATION OF THE SPHERICAL HARMONIC GRAVITY SERIES FOR NUMERICAL ORBITAL PROPAGATION
AAS 15-440 RADIALLY ADAPTIVE EVALUATION OF THE SPHERICAL HARMONIC GRAVITY SERIES FOR NUMERICAL ORBITAL PROPAGATION Austin B. Probe, * Brent Macomber,* Julie I. Read,* Robyn M. Woollands,* and John L. Junkins
More informationHM NAUTICAL ALMANAC OFFICE NAO TECHNICAL N OTE
Council for the Central Laboratory of the Research Councils HM NAUTICAL ALMANAC OFFICE NAO TECHNICAL N OTE No. 73 2003 June Approximate ephemerides and phenomena of natural satellites by D.B. Taylor and
More informationSpace Surveillance using Star Trackers. Part I: Simulations
AAS 06-231 Space Surveillance using Star Trackers. Part I: Simulations Iohan Ettouati, Daniele Mortari, and Thomas Pollock Texas A&M University, College Station, Texas 77843-3141 Abstract This paper presents
More informationNUMERICAL INTEGRATION OF A SATELLITE ORBIT WITH KS TRANSFORMATION
NUMERICAL INTEGRATION OF A SATELLITE ORBIT WITH KS TRANSFORMATION Jhonathan Murcia Piñeros 1 Maxime Koffi 1 Helio Koiti Kuga 1, 2 jhonathan.pineros@inpe.br koffi.maxime@yahoo.com helio.kuga@inpe.br 1 INPE
More information