Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Julia”

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 19 records

Julia Programming Language Benchmark Using a Flight Simulation

Julia’s goal to provide scripting language ease-of-coding with compiled language speed is explored. The runtime speed of the relatively new Julia programming language is assessed against other commonly used languages including Python, Java, and C++. An industry-standard missile and rocket simulation, coded in multiple languages, was used as a test bench for runtime speed. All language versions of the simulation, including Julia, were coded to a highly-developed object-oriented simulation architecture tailored specifically for time-domain flight simulation. A “speed-of-coding” second-dimension is plotted against runtime for each language to portray a space that characterizes Julia’s scripting language efficiencies in the context of the other languages. With caveats, Julia runtime speed was found to be in the class of compiled or semi-compiled languages. However, some factors that affect runtime speed at the cost of ease-of-coding are shown. Julia’s built-in functionality for multi-core processing is briefly examined as a means for obtaining even faster runtime speed. The major contribution of this research to the extensive language benchmarking body-of-work is comparing Julia to other mainstream languages using a complex flight simulation as opposed to benchmarking with single algorithms.

Sells, Ray

Experiments with Julia for Astrodynamics Applications

Julia's potential for solving complex astrodynamics problems is studied. Julia is a high-level, new, dynamic programming language with performance approaching C/Fortran and has features like inbuilt parallelism, variable accuracy, integrated numerical libraries and direct C and Fortran interfaces. Two astrodynamics problems are solved in Julia: 1) Lambert's problem, using the vercosine formulation and 2) trajectory integration. Implemented algorithms are compared with C and Fortran based counterparts on key performance parameters (speed, development effort, etc). Using Julia for fast and reliable astrodynamics software development is also discussed.

Arora, Nitin

Julia Programming Language Benchmark Using a Flight Simulation

Julia is a relatively new computer language that aims to reduce the challenge for mathmodelers to develop fast computer tools and simulations. It potentially combines the ease-of-coding feature of scripting languages (like Python) with the performance of compiled languages (like C++). A key question for Julia application to the simulation domain is, “Can Julia, with its obvious coding simplicity, provide runtime speeds comparable to conventional compiled languages for flight simulation?” A unique combination of existing elements can be employed to address the previous question: Extensively documented object-oriented simulation architecture, Industry standard rocket flight simulation, Separate versions (C++, Java, and Python) already benchmarked.

Sells, Ray

Julia Language 1.1 Ephemeris Reader and Gravitational Modeling Program for Solar System Bodies

This paper analyzes the advancements to the Julia Language 1.1 Ephemeris and Physical constants Reader including the addition of gravitational modeling. Originally written in MATLAB, this Julia Language program is intended to be used in for trajectory design. Written in an open-source coding language, this ephemeris reader can output the state of planetary bodies including asteroids as well as other constants such as gravitational parameters. Two primary methods were chosen to calculate the gravitational potentials which include polyhedral modeling and spherical harmonics.

Gray, Brennan

Polarimetric observations of 9 Metis, 15 Eunomia, 89 Julia, and other asteroids.

White-light polarization curves for a number of asteorids are presented. All have well-developed negative branches, suggestive of relatively dark, texturally complex surfaces - possibly regoliths. Geometric albedos, estimated from the positive branch of the polarization curves, show that a significant dispersion in surface albedo exists within the asteroid belt.

Veverka, J.

Performance Enhancement of a Computational Persistent Homology Package

In recent years, persistent homology has become an attractive method for data analysis. It captures topological features, such as connected components, holes, voids, etc., from a point cloud by finding out when these features appear and disappear in the filtration sequence. In this project, we focus on improving the performance of Eirene, a fancy computational persistent homology package. Eirene is a 5000-line opensource software implemented by using the dynamic programming language Julia. We use the Julia profiling tools to identify the performance bottlenecks and develop different methods to manage the bottlenecks, including the parallelization of some time-consuming functions on the multicore/manycore hardware. The empirical results show that the performance can be greatly improved.

Profiling

Tuning the Performance of a Computational Persistent Homology Package

In recent years, persistent homology has become an attractive method for data analysis. It captures topological features, such as connected components, holes, and voids from point cloud data and summarizes the way in which these features appear and disappear in a filtration sequence. In this project, we focus on improving the performanceof Eirene, a computational package for persistent homology. Eirene is a 5000-line open-source software library implemented in the dynamic programming language Julia. We use the Julia profiling tools to identify performance bottlenecks and develop novel methods to manage them, including the parallelization of some time-consuming functions on multicore/manycore hardware. Empirical results show that performance can be greatly improved.

Persistent Homology

Origins of etioporphyrins in sediments - Evidence from stable carbon isotopes

In samples of the Julia Creek and Condor oil shales (Australia, Albian, and early Tertiary, respectively) etioporphyrin III is significantly depleted in C-13 (4 per mil) relative to porphyrins derived from chlorophylls. This isotopic difference suggests a large contribution from some independent source. The haem group found in cytochromes derived from microbial sources is the most likely candidate.

Boreham, Christopher J.

Origin of petroporphyrins. 2. Evidence from stable carbon isotopes

Compared with the carbon-13 isotopic composition of the ubiquitous C32DPEP (DPEP, deoxophylloerythroetioporphyrin) the heavy but equivalent carbon-13 isotopic composition for the porphyrin structures 15(2)-methyl-15,17-ethano-17-nor-H-C30DPEP and 15,17-butano-, 13,15-ethano-13(2),17-propano-, and 13(1)-methyl-13,15-ethano-13(2),17-propanoporphyrin suggests a common precursor, presumably chlorophyll c, for these petroporphyrins isolated from the marine Julia Creek oil shale and the lacustrine Condor oil shale. Similarly, the heavy but variable carbon-13 isotopic composition of 7-nor-H-C31DPEP compared with C32DPEP is consistent with an origin from both chlorophyll b and chlorophyll c3. The equivalent carbon-13 isotopic composition for 13(2)-methyl-C33DPEP compared with C32DPEP suggests a common origin resulting from a weighted average of chlorophyll inputs.

NASA Program Exobiology

Designing the Lunar Regolith Excavation Competition

The project assigned this summer involves designing a lunar regolith mining robotics competition. This process involves consulting several assets available at the Kennedy Space Center. The process involves several steps. The first step is to determine the requirements for the competition. Once these requirements are determined, the dimensions of the playing field are drawn up, first by hand, and then using computer models. After these drawings are tentatively decided upon, the cost of materials must be determined, so as to fit within the allotted budget for the project. The materials are to then be ordered, assembled, broken down, and stored throughout the duration of the competition. We must also design the advertisements and logos for the competition. This is to market and publicize the competition to college level teams. We must also determine the rules for the competition so as to have uniform requirements for all teams. Once these processes are completed, the competition can be finalized and publicized for the public. The contributing parties are Greg Galloway, Robert Mueller, Susan Sawyer, Gloria Murphy, Julia Nething, and Cassandra Liles.

Le, Christopher

Longitudinal Differences of Ionospheric Vertical Density Distribution and Equatorial Electrodynamics

Accurate estimation of global vertical distribution of ionospheric and plasmaspheric density as a function of local time, season, and magnetic activity is required to improve the operation of space-based navigation and communication systems. The vertical density distribution, especially at low and equatorial latitudes, is governed by the equatorial electrodynamics that produces a vertical driving force. The vertical structure of the equatorial density distribution can be observed by using tomographic reconstruction techniques on ground-based global positioning system (GPS) total electron content (TEC). Similarly, the vertical drift, which is one of the driving mechanisms that govern equatorial electrodynamics and strongly affect the structure and dynamics of the ionosphere in the low/midlatitude region, can be estimated using ground magnetometer observations. We present tomographically reconstructed density distribution and the corresponding vertical drifts at two different longitudes: the East African and west South American sectors. Chains of GPS stations in the east African and west South American longitudinal sectors, covering the equatorial anomaly region of meridian approx. 37 deg and 290 deg E, respectively, are used to reconstruct the vertical density distribution. Similarly, magnetometer sites of African Meridian B-field Education and Research (AMBER) and INTERMAGNET for the east African sector and South American Meridional B-field Array (SAMBA) and Low Latitude Ionospheric Sensor Network (LISN) are used to estimate the vertical drift velocity at two distinct longitudes. The comparison between the reconstructed and Jicamarca Incoherent Scatter Radar (ISR) measured density profiles shows excellent agreement, demonstrating the usefulness of tomographic reconstruction technique in providing the vertical density distribution at different longitudes. Similarly, the comparison between magnetometer estimated vertical drift and other independent drift observation, such as from VEFI onboard Communication/Navigation Outage Forecasting System (C/NOFS) satellite and JULIA radar, is equally promising. The observations at different longitudes suggest that the vertical drift velocities and the vertical density distribution have significant longitudinal differences; especially the equatorial anomaly peaks expand to higher latitudes more in American sector than the African sector, indicating that the vertical drift in the American sector is stronger than the African sector.

total electron content

NASA Engineering and Safety Center Technical Bulletin No. 21-02: Genesis Flight Mechanics Simulation

Genesis Flight Mechanics Simulation The NASA Engineering and Safety Center (NESC) consolidated and modernized a suite of legacy flight mechanics simulations, including the Flight Analysis and Simulation Tool (FAST), resulting in Genesis, a generic, multi-vehicle, variable-degree-of-freedom flight mechanics simulation for ascent, aerocapture, entry, descent, and landing (A2EDL) trajectory design. Genesis is more flexible, capable, and performant than FAST. It enables trajectory optimization and interactive trajectory generation. Its interoperability with Copernicus, an exo-atmospheric and interplanetary trajectory design tool, facilitates end-to-end trajectory optimization across all mission phases. Genesis is implemented in Julia, a new language for technical computing that combines the ease of use of scripting languages with the run-time performance of compiled languages.

EDL

Using GPUs and the Parameterization Method for Rapid Search and Refinement of Connections between Tori in Periodically Perturbed Planar Circular Restricted 3-Body Problems

When the planar circular restricted 3-body problem is periodically perturbed, most unstable periodic orbits become invariant tori. However, 2D Poincare ́ sections no longer work to find their manifolds’ intersections; new methods are needed. In this study, we first review a method of restricting the intersection search to only certain manifold subsets. We then implement this search using Julia and OpenCL, representing the manifolds as triangular meshes and gaining a 30x speedup using GPUs. We finally show how to use manifold parametrizations to refine the ap- proximate connections found in the mesh search. We demonstrate the tools on the planar elliptic RTBP.

de la Llave, Rafael

Pterodactyl: Coupled 6-Dof Integration of Guidance and Control Algorithms in Genesis

The NASA-funded Pterodactyl project seeks to advance the current state-of-the-art for entry vehicles by developing novel guidance and control technologies for Deployable Entry Vehicles (DEVs). This paper builds upon the Pterodactyl architecture that employed eight individually articulating flaps with two options for guidance, bank angle modulation with the Fully Numerical Predictor Corrector Entry Guidance technique (FNPEG) and angle of attack and sideslip modulation with FNPEG uncoupled range control (FNPEG URC). These, with a linear quadratic regulator (LQR) controller, had previously been presented in uncoupled 3-DOF trajectories. This work will show results from fully coupled 6-DOF simulations, leveraging recent advancements in trajectory simulation software, namely the Julia-based Genesis package. Preliminary results show good performance for angle of attack and sideslip modulation when using a controller designed at high dynamic pressure conditions.

range targeting

Pterodactyl: 6-DOF Integration of Guidance and Control Algorithms in Genesis

The NASA-funded Pterodactyl project seeks to advance the current state-of-the-art for entry vehicles by developing novel guidance and control technologies for Deployable Entry Vehicles. This paper builds upon the Pterodactyl architecture that employed eight individually articulating flaps with two options for guidance, bank angle modulation with the Fully Numerical Predictor Corrector Entry Guidance (FNPEG) technique and angle of attack and sideslip modulation with FNPEG uncoupled range control. These, with a linear quadratic regulator controller, have previously been presented in separate 3-DOF trajectory simulations, one for translational motion and another for rotational dynamics. This work will show results for integrated 6-DOF simulations, leveraging recent advancements in trajectory simulation tools, namely the Julia-based Genesis package. Results show similar performance for both schemes when compared to the previously presented 3-DOF results once the controllers were adequately tuned. Angle of attack and sideslip modulation necessitated a controller designed at high dynamic pressure conditions, and bank angle modulation necessitated a controller designed at lower dynamic pressure conditions. Notably, angle of attack and sideslip modulation could not achieve the parachute deploy point target Mach number of 2, when only tuning the controller gains, thus future work will be needed to tune guidance specific parameters to achieve acceptable range targeting and guidance command tracking.

DEV

Pterodactyl: 6-DOF Integration of Guidance and Control Algorithms in Genesis

The NASA-funded Pterodactyl project seeks to advance the current state-of-the-art for entry vehicles by developing novel guidance and control technologies for Deployable Entry Vehicles. This paper builds upon the Pterodactyl architecture that employed eight individually articulating flaps with two options for guidance, bank angle modulation with the Fully Numerical Predictor Corrector Entry Guidance (FNPEG) technique and angle of attack and sideslip modulation with FNPEG uncoupled range control. These, with a linear quadratic regulator controller, have previously been presented in separate 3-DOF trajectory simulations, one for translational motion and another for rotational dynamics. This work will show results for integrated 6-DOF simulations, leveraging recent advancements in trajectory simulation tools, namely the Julia-based Genesis flight simulation tool. Results show similar performance for both schemes when compared to the previously presented 3-DOF results once the controllers were adequately tuned. Angle of attack and sideslip angle modulation necessitated a controller designed at high dynamic pressure conditions, and bank angle modulation necessitated a controller designed at lower dynamic pressure conditions. Notably, the angle of attack and sideslip modulation approach could not achieve the parachute deploy point target Mach number of 2, when only tuning the controller gains, so future work will be needed to improve guidance and control tracking.

EDL

Global Carbon Budget 2019

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere – the “global carbon budget” – is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify the five major components of the global carbon budget and their uncertainties. Fossil CO2 emissions (E(FF)) are based on energy statistics and cement production data, while emissions from land use change (E(LUC)), mainly deforestation, are based on land use and land use change data and bookkeeping models. Atmospheric CO2 concentration is measured directly and its growth rate (G(ATM)) is computed from the annual changes in concentration. The ocean CO2 sink (S(OCEAN)) and terrestrial CO2 sink (S(LAND)) are estimated with global process models constrained by observations. The resulting carbon budget imbalance (B(IM)), the difference between the estimated total emissions and the estimated changes in the atmosphere, ocean, and terrestrial biosphere, is a measure of imperfect data and understanding of the contemporary carbon cycle. All uncertainties are reported as ±1σ. For the last decade available (2009–2018), E(FF) was 9.5±0.5 GtC/yr, E(LUC) 1.5±0.7 GtC/yr, G(ATM) 4.9±0.02 GtC/yr (2.3±0.01 ppm/yr), S(OCEAN) 2.5±0.6 GtC/yr, and S(LAND) 3.2±0.6 GtC/yr, with a budget imbalance B(IM) of 0.4 GtC/yr indicating overestimated emissions and/or underestimated sinks. For the year 2018 alone, the growth in E(FF) was about 2.1 % and fossil emissions increased to 10.0±0.5 GtC/yr, reaching 10 GtC/yr for the first time in history, E(LUC) was 1.5±0.7 GtC/yr, for total anthropogenic CO2 emissions of 11.5±0.9 GtC/yr (42.5±3.3 GtCO2). Also for 2018, G(ATM) was 5.1±0.2 GtC/yr (2.4±0.1 ppm/yr), S(OCEAN) was 2.6±0.6 GtC/yr, and S(LAND) was 3.5±0.7 GtC/yr, with a B(IM) of 0.3 GtC. The global atmospheric CO2 concentration reached 407.38±0.1 ppm averaged over 2018. For 2019, preliminary data for the first 6–10 months indicate a reduced growth in E(FF) of +0.6 % (range of −0.2 % to 1.5 %) based on national emissions projections for China, the USA, the EU, and India and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. Overall, the mean and trend in the five components of the global carbon budget are consistently estimated over the period 1959–2018, but discrepancies of up to 1 GtC/yr persist for the representation of semi-decadal variability in CO2 fluxes. A detailed comparison among individual estimates and the introduction of a broad range of observations shows (1) no consensus in the mean and trend in land use change emissions over the last decade, (2) a persistent low agreement between the different methods on the magnitude of the land CO2 flux in the northern extra-tropics, and (3) an apparent underestimation of the CO2 variability by ocean models outside the tropics. This living data update documents changes in the methods and data sets used in this new global carbon budget and the progress in understanding of the global carbon cycle compared with previous publications of this data set (Le Quéré et al., 2018a, b, 2016, 2015a, b, 2014, 2013).

global carbon budget 2019