Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Vector Processing”

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 1,171 records · Page 65

CHESS 2025: Waveform LiDAR data from NEON AOP surveys

This dataset provides Level 1 (L1) full-waveform light detection and ranging (LiDAR) data collected for the 2025 Colorado Headwaters Ecological Spectroscopy Study (CHESS). These data were acquired to enable characterization of vegetation structure and other three-dimensional features of the land surface, and to evaluate structural changes that may have occurred between a prior LiDAR acquisition in 2018 and the 2025 overflight. Waveform LiDAR data can provide more detailed information about objects on the ground than discrete point clouds typically do, and they are often used for granular target segmentation and characterization of subcanopy vegetation. The data were acquired over three study domains in the Upper Gunnison river basin: the upper East River watershed (CRBU); Almont Triangle and Taylor Canyon (ALMO); and Upper Taylor River watershed (UPTA) between 2025-06-13 and 2025-07-15. LiDAR data were acquired using the Optech Galaxy Prime Airborne LiDAR Terrain Mapper onboard the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP). These are the primary waveform LiDAR data delivered by NEON and are provided per flightline in compressed Pulsewaves format, an open-source binary file standard. A Pulsewaves object comprises a two files: a pulse (.pls) file, which stores the geographic origin, outgoing vector, and metadata for every laser pulse emitted by the scanner, and a wave file (.wvs), which stores the sequential amplitude samples of the outgoing pulse and the returning signals. The files are published here in their compressed forms (.plz, .wvz). All waveform data were processed following the theoretical workflow described in the NEON L0-to-L1 Waveform LiDAR Algorithm Theoretical Basis Document (Krause and Goulden 2022a); however, the Pulsewaves output format differs from a legacy format described in that document. Waveform amplitude samples are recorded at 1 nanosecond intervals. All coordinates are provided in meters. Horizontal coordinates are referenced in Universal Transverse Mercator (UTM) zone 13N and the World Geodetic System (WGS) 1984 ensemble datum. Elevations are referenced to Geoid12A. Waveform data for the UPTA survey area were collected without incident and the published records are complete. However, both the ALMO and CRBU collections experienced issues that resulted in incomplete data for those areas. On collection day 2018-06-16 a hardware failure caused the waveform digitizer to lose data from the eastern edge of the ALMO site (Figure 22). The waveform data for flightlines 2–20 could not be extracted from the digitizer, and the data proved unrecoverable. As a result, a portion of the site does not have coverage with waveform data. Although no hardware failure was observed during collection over the CRBU area, final waveform files generated by vendor software contained only ~25% of the expected number of return pulses. After discovery, NEON initiated troubleshooting with the vendor. The root cause of the data ablation had not been identified at the time of publication. Additional data will be published in an update to this package if further recovery proves successful. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgement: Field and remote-sensing data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). This work was also supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗

Model error estimation for large space systems

In-flight estimation of large structure model errors may have to be carried out in order to detect inevitable deficiencies in large structure controller/estimator models. These error estimates can most efficiently be obtained by the minimization of a quadratic functional of the model errors and on the subsequent analysis of the resulting optimal model error estimates. An integral operator approach to estimation leads to a geometrical interpretation of the model error estimation process. One of the significant insights gained with this interpretation is that the actual but unknown model errors can be decomposed as the sum of two distinct components that are orthogonal in some sense. One of these components is a so-called minimal error vector that retains many of the significant dynamics of the actual errors. The basic ideas in the model error estimation approach are first set forth with a two-dimensional analogy that has most of the essential features of the general estimation problem. The generalized results are then established and their application to a reference large structure model illustrated.

Rodriguez, G.↗

A Modified Sequence-to-point HVAC Load Disaggregation Algorithm

This paper presents a modified sequence-to-point (S2P) algorithm for disaggregating the heat, ventilation, and air conditioning (HVAC) load from the total building electricity consumption. The original S2P model is convolutional neural network (CNN) based, which uses load profiles as inputs. We propose three modifications. First, the input convolution layer is changed from 1D to 2D so that normalized temperature profiles are also used inputs to the S2P model. Second, a drop-out layer is added to improve adaptability and generalizability so that the model trained in one area can be transferred to other geographical areas without labelled HVAC data. Third, a fine-tuning process is proposed for areas with a small amount of labelled HVAC data so that the pre-trained S2P model can be fine-tuned to achieve higher disaggregation accuracy (i.e., better transferability) in other areas. The model is first trained and tested using smart meter and sub-metered HVAC data collected in Austin, Texas. Then, the trained model is tested on two other areas: Boulder, Colorado and San Diego, California. Simulation results show that the proposed modified S2P algorithm outperforms the original S2P model and the support-vector machine based approach in accuracy, adaptability, and transferability.

Ye, Kai↗

Latent map Gaussian processes for mixed variable metamodeling

Gaussian processes (GPs) are ubiquitously used in sciences and engineering as metamodels. Standard GPs, however, can only handle numerical or quantitative variables. Here we introduce latent map Gaussian processes (LMGPs) that inherit the attractive properties of GPs and are also applicable to mixed data which have both quantitative and qualitative inputs. The core idea behind LMGPs is to learn a continuous, low-dimensional latent space or manifold which encodes all qualitative inputs. To learn this manifold, we first assign a unique prior vector representation to each combination of qualitative inputs. We then use a low-rank linear map to project these priors on a manifold that characterizes the posterior representations. As the posteriors are quantitative, they can be directly used in any standard correlation function such as the Gaussian or Matern. Hence, the optimal map and the corresponding manifold, along with other hyperparameters of the correlation function, can be systematically learned via maximum likelihood estimation. Through a wide range of analytic and real-world examples, we demonstrate the advantages of LMGPs over state-of-the-art methods in terms of accuracy and versatility. In particular, we show that LMGPs can handle variable-length inputs, have an explainable neural network interpretation, and provide insights into how qualitative inputs affect the response or interact with each other. We also employ LMGPs in Bayesian optimization and illustrate that they can discover optimal compound compositions more efficiently than conventional methods that convert compositions to qualitative variables via manual featurization.

42 ENGINEERING↗

Elimination of secular terms from the differential equations for the elements of perturbed two-body motion

In 1961, Sperling linearized and regularized the differential equations of motion of the two-body problem by changing the independent variable from time to fictitious time by Sundman's transformation (r = dt/ds) and by embedding the two-body energy integral and the Laplace vector. In 1968, Burdet developed a perturbation theory which was uniformly valid for all types of orbits using a variation of parameters approach on the elements which appeared in Sperling's equations for the two-body solution. In 1973, Bond and Hanssen improved Burdet's set of differential equations by embedding the total energy (which is a constant when the potential function is explicitly dependent upon time.) The Jacobian constant was used as an element to replace the total energy in a reformulation of the differential equations of motion. In the process, another element which is proportional to a component of the angular momentum was introduced. Recently trajectories computed during numerical studies of atmospheric entry from circular orbits and low thrust beginning in near-circular orbits exhibited numerical instability when solved by the method of Bond and Gottlieb (1989) for long time intervals. It was found that this instability was due to secular terms which appear on the righthand sides of the differential equations of some of the elements. In this paper, this instability is removed by the introduction of another vector integral called the delta integral (which replaces the Laplace Vector) and another scalar integral which removes the secular terms. The introduction of these integrals requires a new derivation of the differential equations for most of the elements. For this rederivation, the Lagrange method of variation of parameters is used, making the development more concise. Numerical examples of this improvement are presented.

Bond, Victor R.↗

An Overview of the NASA ADM-Aeolus Cal/Val Test Flight Campaign

Lidars are uniquely capable of collecting high precision and high spatio-temporal observations that have been used for atmospheric process studies from the ground, aircraft, and space for over 30 years. The Atmospheric Dynamics Mission Aeolus (ADM-Aeolus), the first-ever spaceborne Doppler wind lidar, was developed by the European Space Agency and launched in August 2018. The ADM-Aeolus Atmospheric Laser Doppler Instrument (ALADIN) observes profiles of the component of the wind vector and aerosol optical properties along the instrument’s line-of-sight direction on a global scale. Two airborne lidar systems have been developed at NASA Langley Research Center in recent years that collect measurements in support of several NASA Earth Science Division focus areas. The Doppler Aerosol WiNd (DAWN) lidar measures vertical profiles of horizontal wind speed and direction. The High Altitude Lidar Observatory (HALO) measure high resolution profiles of atmospheric water vapor (WV) and aerosol and cloud optical properties. Though new space-based missions with advanced capabilities such as Aeolus are continuously being conceived and developed, there are limitations in terms of spatial and vertical detail and measurement precision that can be accomplished from space, thus airborne remote sensing observations like those from DAWN and HALO are required to fill these observational gaps as well as to calibrate and validate space-based measurements. Over a two-week period in April 2019 during their Aeolus Cal/Val Test Flight campaign, NASA conducted five research flights over the Eastern Pacific Ocean with the DC-8 aircraft to 1) demonstrate the DAWN Doppler wind lidar, HALO high spectral resolution lidar aerosol/cloud measurement and differential absorption lidar WV measurement capabilities across a range of atmospheric conditions, 2) demonstrate Aeolus Cal/Val flight strategies and compare DAWN and HALO measurements with Aeolus to gain an initial perspective of Aeolus performance, and 3) demonstrate how atmospheric dynamic processes can be resolved and better understood through simultaneous observations of wind, WV, and aerosol profile observations, coupled with numerical model and other remote sensing observations. This paper provides a brief description of the DAWN and HALO instruments, discusses the synergistic observations collected across a wide range of atmospheric conditions sampled during the DC-8 flights, and a summary of the validation of DAWN, HALO, and Aeolus observations and comparisons.

Kristopher M. Bedka↗

Observation of electroweak production of two jets in association with an isolated photon and missing transverse momentum, and search for a Higgs boson decaying into invisible particles at 13 TeV with the ATLAS detector

This paper presents a measurement of the electroweak production of two jets in association with a Zγ pair, with the Z boson decaying into two neutrinos. It also presents a search for invisible or partially invisible decays of a Higgs boson with a mass of 125 GeV produced through vector-boson fusion with a photon in the final state. These results use data from LHC proton–proton collisions at $\sqrt{s}$ = 13 TeV collected with the ATLAS detector and corresponding to an integrated luminosity of 139 fb -1 . The event signature, shared by all benchmark processes considered for the measurements and searches, is characterized by a significant amount of unbalanced transverse momentum and a photon in the final state, in addition to a pair of forward jets. Electroweak $Zγ$ production in association with two jets is observed in this final state with a significance of 5.2 (5.1 expected) standard deviations. The measured fiducial cross-section for this process is 1.31 ± 0.29 fb. An observed (expected) upper limit of 0.37 (0.34$^{+0.15}_{-0.10}$) at 95% confidence level is set on the branching ratio of a 125 GeV Higgs boson to invisible particles, assuming the Standard Model production cross-section. The signature is also interpreted in the context of decays of a Higgs boson into a photon and a dark photon. An observed (expected) 95% CL upper limit on the branching ratio for this decay is set at 0.018 (0.017$^{+0.007}_{-0.005}$), assuming the Standard Model production cross-section for a 125 GeV Higgs boson.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

MAGSAT science investigations

A preliminary magnetization model for the U.S. was produced from MAGSAT data. The double grid processing method resulted in a 100 km source spacing. While some spurious features are present in the map, much of the detail appears to be real. The determination of the optional source spacing for inversion of MAGSAT delta B data is in progress. A fine attitude vector data set was obtained for comparison with the theoretical vector field associated with the magnetization model.

Source record↗

A demonstration of real-time connected element interferometry for spacecraft navigation

Connected element interferometry is a technique of observing a celestial radio source at two spatially separated antennas, and then interfering the received signals to extract the relative phase of the signal at the two antennas. The high precision of the resulting phase delay data type can provide an accurate determination of the angular position of the radio source relative to the baseline vector between the two stations. A connected element interferometer on a 21-km baseline between two antennas at the Deep Space Network's Goldstone, CA tracking complex is developed. Fiber optic links are used to transmit the data at 112 Mbit/sec to a common site for processing. A real-time correlator to process these data in real-time is implemented. The architecture of the system is described, and observational data is presented to characterize the potential performance of such a system. The real-time processing capability offers potential advantages in terms of increased reliability and improved delivery of navigational data for time-critical operations. Angular accuracies of 50-100 nrad are achievable on this baseline.

Edwards, C.↗

Development of a Robust star identification technique for use in attitude determination of the ACE spacecraft

The Advanced Composition Explorer (ACE) spacecraft is designed to fly in a spin-stabilized attitude. The spacecraft will carry two attitude sensors - a digital fine Sun sensor and a charge coupled device (CCD) star tracker - to allow ground-based determination of the spacecraft attitude and spin rate. Part of the processing that must be performed on the CCD star tracker data is the star identification. Star data received from the spacecraft must be matched with star information in the SKYMAP catalog to determine exactly which stars the sensor is tracking. This information, along with the Sun vector measured by the Sun sensor, is used to determine the spacecraft attitude. Several existing star identification (star ID) systems were examined to determine whether they could be modified for use on the ACE mission. Star ID systems which exist for three-axis stabilized spacecraft tend to be complex in nature and many require fairly good knowledge of the spacecraft attitude, making their use for ACE excessive. Star ID systems used for spinners carrying traditional slit star sensors would have to be modified to model the CCD star tracker. The ACE star ID algorithm must also be robust, in that it will be able to correctly identify stars even though the attitude is not known to a high degree of accuracy, and must be very efficient to allow real-time star identification. The paper presents the star ID algorithm that was developed for ACE. Results from prototype testing are also presented to demonstrate the efficiency, accuracy, and robustness of the algorithm.

Woodard, Mark↗

SUMI - The Solar Ultraviolet Magnetograph Investigation

Solar physics has been successful in characterizing the full vector magnetic field in the photosphere, where the ratio of gas pressure to magnetic pressure (Beta) is >1. However, at higher levels in the atmosphere, where Beta <<1 and flares and CMEs are believed to be triggered, observations are difficult, severely limiting the understanding of these processes. In response to this situation, we are developing SUMI (the Solar Ultraviolet Magnetograph Investigation) a unique instrument designed to measure the circular and linear polarization of upper chromospheric Mg II lines (280 nm) and circular polarization of transition region C IV lines (155 nm). To date the telescope mirrors have been built, tested and coated with dielectric stacks designed to reflect only the wavelengths of interest. We have also developed a unique UV polarimeter and completed the design of a high-resolution spectrograph that uses dual toroidal varied- line-space (TVLS) gratings. Incorporating measurements of those components developed so far, the revised estimate of the system throughput exceeds our original estimate by more than an order of magnitude. A sounding rocket flight is anticipated in 2006. Our objectives and progress are detailed in this presentation.

Porter, J. G.↗

Mission and Navigation Design for the 2009 Mars Science Laboratory Mission

NASA s Mars Science Laboratory mission will launch the next mobile science laboratory to Mars in the fall of 2009 with arrival at Mars occurring in the summer of 2010. A heat shield, parachute, and rocket-powered descent stage, including a sky crane, will be used to land the rover safely on the surface of Mars. The direction of the atmospheric entry vehicle lift vector will be controlled by a hypersonic entry guidance algorithm to compensate for entry trajectory errors and counteract atmospheric and aerodynamic dispersions. The key challenges for mission design are (1) develop a launch/arrival strategy that provides communications coverage during the Entry, Descent, and Landing phase either from an X-band direct-to-Earth link or from a Ultra High Frequency link to the Mars Reconnaissance Orbiter for landing latitudes between 30 deg North and 30 deg South, while satisfying mission constraints on Earth departure energy and Mars atmospheric entry speed, and (2) generate Earth-departure targets for the Atlas V-541 launch vehicle for the specified launch/arrival strategy. The launch/arrival strategy employs a 30-day baseline launch period and a 27-day extended launch period with varying arrival dates at Mars. The key challenges for navigation design are (1) deliver the spacecraft to the atmospheric entry interface point (Mars radius of 3522.2 km) with an inertial entry flight path angle error of +/- 0.20 deg (3 sigma), (2) provide knowledge of the entry state vector accurate to +/- 2.8 km (3 sigma) in position and +/- 2.0 m/s (3 sigma) in velocity for initializing the entry guidance algorithm, and (3) ensure a 99% probability of successful delivery at Mars with respect to available cruise stage propellant. Orbit determination is accomplished via ground processing of multiple complimentary radiometric data types: Doppler, range, and Delta-Differential One-way Ranging (a Very Long Baseline Interferometry measurement). The navigation strategy makes use of up to five interplanetary trajectory correction maneuvers to achieve entry targeting requirements. The requirements for cruise propellant usage and atmospheric entry targeting and knowledge are met with ample margins.

Mission design↗

The Updated GEO Population for ORDEM 3.1

The limited availability of data for satellite fragmentations and debris in the geosynchronous orbit (GEO) region creates challenges to building accurate models for the orbital debris environment at such altitudes. Updated methods to properly incorporate and extrapolate measurement data have become a cornerstone of the GEO component in the newest version of the NASA Orbital Debris Engineering Model (ORDEM), ORDEM 3.1. For the GEO region, the Space Surveillance Network (SSN) catalog provides coverage down to a limit of approximately 1 m. A more statistically complete representation of the GEO population for smaller objects, which can pose a high risk to operational spacecraft, is thus dependent on dedicated observations by instruments optimized to observe debris smaller than the SSN cataloging threshold. For ORDEM 3.1, optical data from the Michigan Orbital DEbris Survey Telescope (MODEST) provided the input for building the GEO population down to approximately 30 cm (converting absolute magnitude to size). For smaller sizes, the size distribution of debris in the MODEST dataset was extrapolated down to 10 cm, and orbital parameters were estimated based on the orbits of the larger objects. When compared to previous versions of the model, significant improvements were made to the process of building the GEO population in ORDEM 3.1, both in the assessment of fragmentation debris in the data and assignment of orbital elements within the model. A so-called “debris ring filter,” based on a range of angles between an orbit’s angular momentum vector and that of the stable Laplace plane, was applied to the data to reduce biases from non- GEO objects, such as objects in a GEO-transfer orbit. In addition, a new approach was implemented to assign noncircular mean motions and eccentricities to the fragmentation debris observed by MODEST because the short observation window (5 min) in GEO limits orbit resolution to a circular orbit assumption for assigning orbital parameters. For ORDEM 3.1, non-circular orbital elements were assigned using relationships that were identified between mean motion and the angle between the orbit plane and the stable Laplace plane, as well as between mean motion and eccentricity, based on breakup clouds modeled by the NASA Standard Breakup Model. This approach has yielded a high-fidelity GEO model that has been validated with data from more recent MODEST observation campaigns.

Manis, A.↗

Space Environment Factors Affecting the Performance of International Space Station Materials: The First Two Years of Flight Operations

In this paper, the natural and induced space environment factors affecting materials performance on ISS are described in some detail. The emphasis will be on ISS flight experience and the more significant design and development issues of the last two years. The intent is to identify and document the set of space environment factors, affecting materials, that are producing the largest impacts on the ISS flight hardware verification and acceptance process and on ISS flight operations. Orbital inclination (S1.6 ) and altitude (nominal3S0 km to 400 km altitude) determine the set of natural environment factors affecting the functional life of materials and subsystems on ISS. ISS operates in the F2 region of Earth's ionosphere in well-defined fluxes of atomic oxygen, other ionospheric plasma species, and solar UV, VUV, and x-ray radiation, as well as galactic cosmic rays, trapped radiation, and solar cosmic rays (1,2). The high latitude orbital environment also exposes external surfaces to significantly less well-defined or predictable fluxes of higher energy trapped electrons and auroral electrons (3 ,4). The micrometeoroid and orbital debris environment is an important determinant of spacecraft design and operations in any orbital inclination. Environment factors induced by ISS flight operations include ram-wake effects, magnetic induction voltages arising from flight through Earth's magnetic field, hypergolic thruster plume impingement from proximity operations of visiting vehicles, materials outgassing, venting and dumping of fluids, ISS thruster operations, as well as specific electrical power system interactions with the ionospheric plasma (S-7). ISS must fly in a very limited number of approved flight attitudes leading to location specific environmental exposures and extreme local thermal environments (8). ISS is a large vehicle and produces a deep wake structure from which both ionospheric plasma and neutrals (atomic oxygen) are largely excluded (9-11). At high latitude, the ISS wake may produce a spacecraft charging environment similar to that experienced by the DMSP and Freja satellites (800 to 100 km altitude polar orbits), especially during geo-magnetic disturbances (12-14). ISS is also subject to magnetic induction voltages (VxB L) on conducting structure, a result of high velocity flight through Earth's magnetic field. The magnitude of the magnetic induction voltage varies with location on ISS, as well as the relative orientation of the vehicle velocity vector and planetary magnetic field vector, leading to maximum induction voltages at high latitude (15). The space environment factors, natural and induced, that have had the largest impact on pre-launch ISS flight hardware verification and flight operations during the first two years of ISS flight operations are listed below and grouped according to the physical and chemical processes driving their interaction with ISS materials.

Koontz, Steven L.↗

Electron Scattering by Low-frequency Whistler Waves at Earth’s Bow Shock

Electrons are accelerated to nonthermal energies at shocks in space and astrophysical environments. While shock drift acceleration (SDA) has been considered a key process of electron acceleration at Earth’s bow shock, it has also been recognized that SDA needs to be combined with an additional stochastic process to explain the observed power-law energy spectra. Here, we show mildly energetic (∼0.5 keV) electrons are locally scattered (and accelerated while being confined) by magnetosonic-whistler waves within the shock transition layer, especially when the shock angle is large (θ(sub Bn) approximately equal or greater than 70°). When measured by the Magnetospheric Multiscale mission at a high cadence, ∼0.5 keV electron flux increased exponentially in the shock transition layer. However, the flux profile was not entirely smooth and the fluctuation showed temporal/spectral association with large-amplitude (δB/B ~ 0.3), low-frequency (approximately equal or less than 0.1 Ω(sub ce) where Ω(sub ce) is the cyclotron frequency), obliquely propagating (θ(sub kB) ~ 30°–60°, where θ(sub kB) is the angle between the wave vector and background magnetic field) whistler waves, indicating that the particles were interacting with the waves. Particle simulations demonstrate that, although linear cyclotron resonances with ∼0.5 keV electrons are unlikely due to the obliquity and low frequencies of the waves, the electrons are still scattered beyond 90° pitch angle by (1) resonant mirroring (transit-time damping), (2) non-resonant mirroring, and (3) subharmonic cyclotron resonances. Such coupled nonlinear scattering processes are likely to provide the stochasticity needed to explain the power-law formation.

Oka, M.↗

Improvement of Automated POST Case Success Rate Using Support Vector Machines

During early conceptual design of complex systems, concept down selection can have a large impact upon program life-cycle cost. Therefore, any concepts selected during early design will inherently commit program costs and affect the overall probability of program success. For this reason it is important to consider as large a design space as possible in order to better inform the down selection process. For conceptual design of launch vehicles, trajectory analysis and optimization often presents the largest obstacle to evaluating large trade spaces. This is due to the sensitivity of the trajectory discipline to changes in all other aspects of the vehicle design. Small deltas in the performance of other subsystems can result in relatively large fluctuations in the ascent trajectory because the solution space is non-linear and multi-modal. In order to help capture large design spaces for new launch vehicles, the authors have performed previous work seeking to automate the execution of the industry standard tool, Program to Optimize Simulated Trajectories (POST). This work initially focused on implementation of analyst heuristics to enable closure of cases in an automated fashion, with the goal of applying the concepts of design of experiments (DOE) and surrogate modeling to enable near instantaneous throughput of vehicle cases.3 As noted in [4] work was then completed to improve the DOE process by utilizing a graph theory based approach to connect similar design points.

Zwack, Mathew R.↗

Improvement of Automated POST Case Success Rate Using Support Vector Machines

During early conceptual design of complex systems, concept down selection can have a large impact upon program life-cycle cost. Therefore, any concepts selected during early design will inherently commit program costs and affect the overall probability of program success. For this reason it is important to consider as large a design space as possible in order to better inform the down selection process. For conceptual design of launch vehicles, trajectory analysis and optimization often presents the largest obstacle to evaluating large trade spaces. This is due to the sensitivity of the trajectory discipline to changes in all other aspects of the vehicle design. Small deltas in the performance of other subsystems can result in relatively large fluctuations in the ascent trajectory because the solution space is non-linear and multi-modal [1]. In order to help capture large design spaces for new launch vehicles, the authors have performed previous work seeking to automate the execution of the industry standard tool, Program to Optimize Simulated Trajectories (POST). This work initially focused on implementation of analyst heuristics to enable closure of cases in an automated fashion, with the goal of applying the concepts of design of experiments (DOE) and surrogate modeling to enable near instantaneous throughput of vehicle cases [2]. Additional work was then completed to improve the DOE process by utilizing a graph theory based approach to connect similar design points [3]. The conclusion of the previous work illustrated the utility of the graph theory approach for completing a DOE through POST. However, this approach was still dependent upon the use of random repetitions to generate seed points for the graph. As noted in [3], only 8% of these random repetitions resulted in converged trajectories. This ultimately affects the ability of the random reps method to confidently approach the global optima for a given vehicle case in a reasonable amount of time. With only an 8% pass rate, tens or hundreds of thousands of reps may be needed to be confident that the best repetition is at least close to the global optima. However, typical design study time constraints require that fewer repetitions be attempted, sometimes resulting in seed points that have only a handful of successful completions. If a small number of successful repetitions are used to generate a seed point, the graph method may inherit some inaccuracies as it chains DOE cases from the non-global-optimal seed points. This creates inherent noise in the graph data, which can limit the accuracy of the resulting surrogate models. For this reason, the goal of this work is to improve the seed point generation method and ultimately the accuracy of the resulting POST surrogate model. The work focuses on increasing the case pass rate for seed point generation.

Zwack, Matthew R.↗

A benchmark dataset for Hydrogen Combustion

The generation of reference data for deep learning models is challenging for reactive systems, and more so for combustion reactions due to the extreme conditions that create radical species and alternative spin states during the combustion process. Here, we extend intrinsic reaction coordinate (IRC) calculations with ab initio MD simulations and normal mode displacement calculations to more extensively cover the potential energy surface for 19 reaction channels for hydrogen combustion. A total of ~290,000 potential energies and ~1,270,000 nuclear force vectors are evaluated with a high quality range-separated hybrid density functional, ωB97X-V, to construct the reference data set, including transition state ensembles, for the deep learning models to study hydrogen combustion reaction.

08 HYDROGEN↗