Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “test vector”

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,657 records · Page 92

Attention-Augmented Parametric Kernel Graph Neural Network (APKGNN) for Node Classification

We present a new graph neural network, the Attention-based Parametric-Kernel augmented Graph Neural Network (APKGNN), developed for node classification tasks. Despite extensive work on modeling multi-faceted relationships between connected nodes of a graph, the effect of attention on edge features mapped to relationships has not yet been analyzed through learning representation. This study derives such an attention vector by first calculating node features corresponding to endpoints of an edge and then aggregating these with extracted local intrinsic patches of a given graph to generate augmented local patch vectors. This process uses a parametric kernel based on Gaussian mixture models (GMMs) to embed local neighborhoods of the graph in local patches. The patch vectors then convolve with the above node features to produce an updated node representation. We show that this new learning representation (APKGNN) achieves higher node classification accuracy on tasks - both standard benchmarks (Cora, PubMed, Citeseer) and new experimental short text corpora where nodes correspond to text documents and words. This implementation of the GNN convolution layer outperforms state-of-the-art (SOTA) algorithms, achieving higher training, validation, and test accuracy by a significant margin on three standard benchmark data sets under both SOTA experimental settings and those for new testbeds.

Bose, Avishek↗

The peculiar velocities of rich clusters in the hot and cold dark matter scenarios

We present the results of a study of the peculiar velocities of rich clusters of galaxies. The peculiar motion of rich clusters in various cosmological scenarios is of interest for a number of reasons. Observationally, one can measure the peculiar motion of clusters to greater distances than galaxies because cluster peculiar motions can be determined to greater accuracy. One can also test the slope of distance indicator relations using clusters to see if galaxy properties vary with environment. We have used N-body simulations to measure the amplitude and rms cluster peculiar velocity as a function of bias parameter in the hot and cold dark matter scenarios. In addition to measuring the mean and rms peculiar velocity of clusters in the two models, we determined whether the peculiar velocity vector of a given cluster is well aligned with the gravity vector due to all the particles in the simulation and the gravity vector due to the particles present only in the clusters. We have investigated the peculiar velocities of rich clusters of galaxies in the cold dark matter and hot dark matter galaxy formation scenarios. We have derived peculiar velocities and associated errors for the scenarios using four values of the bias parameter ranging from b = 1 to b = 2.5. The growth of the mean peculiar velocity with scale factor has been determined and compared to that predicted by linear theory. In addition, we have compared the orientation of force and velocity in these simulations to see if a program such as that proposed by Bertschinger and Dekel (1989) for elliptical galaxy peculiar motions can be applied to clusters. The method they describe enables one to recover the density field from large scale redshift distance samples. The method makes it possible to do this when only radial velocities are known by assuming that the velocity field is curl free. Our analysis suggests that this program if applied to clusters is only realizable for models with a low value of the bias parameter, i.e., models in which the peculiar velocities of clusters are large enough that the errors do not render the analysis impracticable.

Rhee, George F.↗

Bringing "Scientific Expeditions" Into the Schools

Two new technologies, the FASTexpedition and Remote FAST, have been developed that provide remote, 3D, high resolution, dynamic, interactive viewing of scientific data (such as simulations or measurements of fluid dynamics). The FASTexpedition permits one to access scientific data from the World Wide Web, take guided expeditions through the data, and continue with self controlled expeditions through the data. Remote FAST permits collaborators at remote sites to simultaneously view an analysis of scientific data being controlled by one of the collaborators. Control can be transferred between sites. These technologies are now being used for remote collaboration in joint university, industry, and NASA projects in computational fluid dynamics (CFD) and wind tunnel testing. Also, NASA Ames Research Center has initiated a project to make scientific data and guided expeditions through the data available as FASTexpeditions on the World Wide Web for educational purposes. Previously, remote visualiZation of dynamic data was done using video format (transmitting pixel information) such as video conferencing or MPEG movies on the Internet. The concept for this new technology is to send the raw data (e.g., grids, vectors, and scalars) along with viewing scripts over the Internet and have the pixels generated by a visualization tool running on the viewer's local workstation. The visualization tool that is currently used is FAST (Flow Analysis Software Toolkit). The advantages of this new technology over using video format are: 1. The visual is much higher in resolution (1280xl024 pixels with 24 bits of color) than typical video format transmitted over the network. 2. The form of the visualization can be controlled interactively (because the viewer is interactively controlling the visualization tool running on his workstation). 3. A rich variety of guided expeditions through the data can be included easily. 4. A capability is provided for other sites to see a visual analysis of one site as the analysis is interactively performed. Control of the analysis can be passed from site to site. 5. The scenes can be viewed in 3D using stereo vision. 6. The network bandwidth used for the visualization using this new technology is much smaller than when using video format. (The measured peak bandwidth used was 1 Kbit/sec whereas the measured bandwidth for a small video picture was 500 Kbits/sec.)

Watson, Val↗

Anisotropic Solution Adaptive Unstructured Grid Generation Using AFLR

An existing volume grid generation procedure, AFLR3, was successfully modified to generate anisotropic tetrahedral elements using a directional metric transformation defined at source nodes. The procedure can be coupled with a solver and an error estimator as part of an overall anisotropic solution adaptation methodology. It is suitable for use with an error estimator based on an adjoint, optimization, sensitivity derivative, or related approach. This offers many advantages, including more efficient point placement along with robust and efficient error estimation. It also serves as a framework for true grid optimization wherein error estimation and computational resources can be used as cost functions to determine the optimal point distribution. Within AFLR3 the metric transformation is implemented using a set of transformation vectors and associated aspect ratios. The modified overall procedure is presented along with details of the anisotropic transformation implementation. Multiple two-and three-dimensional examples are also presented that demonstrate the capability of the modified AFLR procedure to generate anisotropic elements using a set of source nodes with anisotropic transformation metrics. The example cases presented use moderate levels of anisotropy and result in usable element quality. Future testing with various flow solvers and methods for obtaining transformation metric information is needed to determine practical limits and evaluate the efficacy of the overall approach.

Marcum, David L.↗

Properties, Propagation, and Excitation of EMIC Waves Properties, Propagation, and Excitation of EMIC Waves

Electromagnetic ion cyclotron (EMIC) waves (0.1-5 Hz) play an important role in particle dynamics in the Earth's magnetosphere. EMIC waves are preferentially excited in regions where hot anisotropic ions and cold dense plasma populations spatially overlap. While the generation region of EMIC waves is usually on or near the magnetic equatorial plane in the inner magnetosphere, EMIC waves have both equatorial and off-equator source regions on the dayside in the compressed outer magnetosphere. Using field and plasma measurements from the Magnetospheric Multiscale (MMS) mission, we perform a case study of EMIC waves and associated local plasma conditions observed on 19 October 2015. From 0315 to 0810 UT, before crossing the magnetopause into the magnetosheath, all four MMS spacecraft detected long-lasting He(exp +)-band EMIC wave emissions around local noon (MLT = 12.7 - 14.0) at high L-shells (L = 8.8 - 15.2) and low magnetic latitudes (MLAT = -21.8deg - -30.3deg). Energetic (greater than 1 keV) and anisotropic ions were present throughout this event that was in the recovery phase of a weak geomagnetic storm (min. Dst = -48 nT at 1000 UT on 18 October 2015). The testing of linear theory suggests that the EMIC waves were excited locally. Although the wave event is dominated by small normal angles, its polarization is mixed with right- and left-handedness and its propagation is bi-directional with regard to the background magnetic field. The short inter-spacecraft distances (as low as ~15 km) of the MMS mission make it possible to accurately determine the k vector of the waves using the phase difference technique. Preliminary analysis finds that the k vector magnitude, phase speed, and wavelength of the 0.3-Hz wave packet at 0453:55 UT are 0.005 km(exp -1), 372.9 km/s, and 1242.9 km, respectively.

Zhang, Jichun↗

A Study of a Magnetic Cloud Propagating through Large-Amplitude Alfvén Waves

We discuss Wind observations of a long and slow magnetic cloud (MC) propagating through large-amplitude Alfvén waves (LAAWs). The MC axis has a strong component along GSE X, as also confirmed by a Grad-Shafranov reconstruction. It is overtaking the solar wind at a speed roughly equal to the upstream Alfvén speed, leading to a weak shock wave âLij17 hrs. ahead. We give evidence to show that the nominal sheath region is populated by LAAWs: (i) a well-defined deHoffman-Teller frame in which there is excellent correlation between the field and flow vectors, (ii) constant field and total pressure, and (iii) an Alfvén ratio (i.e. ratio of kinetic-to-magnetic energy of the fluctuations) near unity at frequencies much lower than the ion cyclotron frequency in the spacecraft frame. In the region where the LAAWs approach the MCâAZs front boundary there are field and flow discontinuities. At the first, magnetic reconnection is taking place, as deduced from a stress balance test (Walen test). This severs connection of some field lines to the Sun and the solar wind strahl disappears. There follows a aLij2-hour interval where the magnetic field strength is diminished while pressure balance is maintained. Here the bi-directionality of the suprathermal electron flows is intermittently disrupted. This interval ends with a slow expansion fan downstream of which there is a dropout of halo electrons just inside the front boundary of the MC. This study illustrates an untypical case of a slow MC interacting with LAWWs in the slow solar wind.

C. J. Farrugia↗

Requirement Discovery Using Embedded Knowledge Graph With ChatGPT

The field of Advanced Air Mobility (AAM) is witnessing a transformation with innovations such as electric aircraft and increasingly automated airspace operations. Within AAM, the Urban Air Mobility (UAM) con-cept focuses on providing air-taxi services in densely populated urban areas. This research introduces the utilization of Large Language Models (LLMs), such as OpenAI's GPT-4, to enhance the UAM Requirement discovery process. This study explores two distinct approaches to leverage LLMs in the context of UAM Requirement discovery. The first approach evaluates the LLM's ability to provide responses without relying on additional outside systems, such as a relational or graph database. Instead, a vector store provides relevant information to the LLM based on the user’s question, a process known as Retrieval Augmented Generation (RAG). The second approach integrates the LLM with a graph database. The LLM acts as an intermediary between the user and the graph database, translating user questions into cypher queries for the database and database responses into human-readable answers for the user. Our team implemented and tested both solutions to analyze require-ments within a UAM dataset. This paper will talk about our approaches, implementations, and findings related to both approaches.

systems engineering↗

Requirement Discovery Using Embedded Knowledge Graph with ChatGPT

The field of Advanced Air Mobility (AAM) is witnessing a transformation with innovations such as electric aircraft and increasingly automated airspace operations. Within AAM, the Urban Air Mobility (UAM) concept focuses on providing air-taxi services in densely populated urban areas. This research introduces the utilization of Large Language Models (LLMs), such as OpenAI's GPT-4, to enhance the UAM Requirement discovery process. This study explores two distinct approaches to leverage LLMs in the context of UAM Requirement discovery. The first approach evaluates the LLM's ability to provide responses without relying on additional outside systems, such as a relational or graph database. Instead, a vector store provides relevant information to the LLM based on the user’s question, a process known as Retrieval Augmented Generation (RAG). The second approach integrates the LLM with a graph database. The LLM acts as an intermediary between the user and the graph database, translating user questions into cypher queries for the database and database responses into human-readable answers for the user. Our team implemented and tested both solutions to analyze requirements within a UAM dataset. This paper will talk about our approaches, implementations, and findings related to both approaches.

systems engineering↗

Requirement Discovery Using Embedded Knowledge Graph With ChatGPT - Poster

The field of Advanced Air Mobility (AAM) is witnessing a transformation with innovations such as electric aircraft and increasingly automated airspace operations. Within AAM, the Urban Air Mobility (UAM) con-cept focuses on providing air-taxi services in densely populated urban areas. This research introduces the utilization of Large Language Models (LLMs), such as OpenAI's GPT-4, to enhance the UAM Requirement discovery process. This study explores two distinct approaches to leverage LLMs in the context of UAM Requirement discovery. The first approach evaluates the LLM's ability to provide responses without relying on additional outside systems, such as a relational or graph database. Instead, a vector store provides relevant information to the LLM based on the user’s question, a process known as Retrieval Augmented Generation (RAG). The second approach integrates the LLM with a graph database. The LLM acts as an intermediary between the user and the graph database, translating user questions into cypher queries for the database and database responses into human-readable answers for the user. Our team implemented and tested both solutions to analyze require-ments within a UAM dataset. This paper will talk about our approaches, implementations, and findings related to both approaches.

systems engineering↗

Roofline Analysis in the Intel® Advisor to Deliver Optimized Performance for applications on Intel® Xeon Phi™ Processor

In this session we show, in two case studies, how the roofline feature of Intel Advisor has been utilized to optimize the performance of kernels of the XGC1 and PICSAR codes in preparation for Intel Knights Landing architecture. The impact of the implemented optimizations and the benefits of using the automatic roofline feature of Intel Advisor to study performance of large applications will be presented. This demonstrates an effective optimization strategy that has enabled these science applications to achieve up to 4.6 times speed-up and prepare for future exascale architectures. # Goal/Relevance of Session The roofline model [1,2] is a powerful tool for analyzing the performance of applications with respect to the theoretical peak achievable on a given computer architecture. It allows one to graphically represent the performance of an application in terms of operational intensity, i.e. the ratio of flops performed and bytes moved from memory in order to guide optimization efforts. Given the scale and complexity of modern science applications, it can often be a tedious task for the user to perform the analysis on the level of functions or loops to identify where performance gains can be made. With new Intel tools, it is now possible to automate this task, as well as base the estimates of peak performance on measurements rather than vendor specifications. The goal of this session is to demonstrate how the roofline feature of Intel Advisor can be used to balance memory vs. computation related optimization efforts and effectively identify performance bottlenecks. A series of typical optimization techniques: cache blocking, structure refactoring, data alignment, and vectorization illustrated by the kernel cases will be addressed. # Description of the codes ## XGC1 The XGC1 code [3] is a magnetic fusion Particle-In-Cell code that uses an unstructured mesh for its Poisson solver that allows it to accurately resolve the edge plasma of a magnetic fusion device. After recent optimizations to its collision kernel [4], most of the computing time is spent in the electron push (pushe) kernel, where these optimization efforts have been focused. The kernel code scaled well with MPI+OpenMP but had almost no automatic compiler vectorization, in part due to indirect memory addresses and in part due to low trip counts of low-level loops that would be candidates for vectorization. Particle blocking and sorting have been implemented to increase trip counts of low-level loops and improve memory locality, and OpenMP directives have been added to vectorize compute-intensive loops that were identified by Advisor. The optimizations have improved the performance of the pushe kernel 2x on Haswell processors and 1.7x on KNL. The KNL node-for-node performance has been brought to within 30% of a NERSC Cori phase I Haswell node and we expect to bridge this gap by reducing the memory footprint of compute intensive routines to improve cache reuse. ## PICSAR is a Fortran/Python high-performance Particle-In-Cell library targeting at MIC architectures first designed to be coupled with the PIC code WARP for the simulation of laser-matter interaction and particle accelerators. PICSAR also contains a FORTRAN stand-alone kernel for performance studies and benchmarks. A MPI domain decomposition is used between NUMA domains and a tile decomposition (cache-blocking) handled by OpenMP has been added for shared-memory parallelism and better cache management. The so-called current deposition and field gathering steps that compose the PIC time loop constitute major hotspots that have been rewritten to enable more efficient vectorization. Particle communications between tiles and MPI domain has been merged and parallelized. All considered, these improvements provide speedups of 3.1 for order 1 and 4.6 for order 3 interpolation shape factors on KNL configured in SNC4 quadrant flat mode. Performance is similar between a node of cori phase 1 and KNL at order 1 and better on KNL by a factor 1.6 at order 3 with the considered test case (homogeneous thermal plasma).

97 MATHEMATICS AND COMPUTING↗

Stability of a tachyon braneworld

Within the braneworld paradigm the tachyonic scalar field has been used to generate models that attempt to solve some of the open problems that physics faces nowadays, both in cosmology and high energy physics as well. When these field configurations are produced by the interplay of higher dimensional warped gravity with some matter content, braneworld models must prove to be stable under the whole set of small fluctuations of the gravitational and matter fields background, among other consistency tests. Here we present a complete proof of the stability under scalar perturbations of tachyonic thick braneworlds with an embedded maximally symmetric 4D space-time, revealing its physical consistency. This family of models contains a recently reported tachyonic de Sitter thick braneworld which possesses a series of appealing properties. These features encompass complete regularity, asymptotic flatness (instead of being asymptotically dS or AdS) even when it contains a negative bulk cosmological constant, a relevant 3-brane with dS metric which naturally arises from the full set of field equations of the 5D background (it is not imposed), qualitatively describing the inflationary epochs of our Universe, and a graviton spectrum with a single zero mode bound state that accounts for the 4D graviton localised on the brane and is separated from the continuum of Kaluza-Klein massive graviton excitations by a mass gap. The presence of this mass gap in the graviton spectrum makes the extra-dimensional corrections to Newton’s law decay exponentially. Gauge vector fields with a single massless bound state in its mass spectrum are also localised on this braneworld model a fact that allows us to recover the Coulomb’s law of our 4D world. All these properties of the above referred tachyonic braneworld together with the positive stability analysis provided in this work, constitute a firm step towards the construction of realistic cosmological models within the braneworld paradigm.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

LncRNA NKILA upregulation mediates oxygen glucose deprivation/re-oxygenation-induced neuronal cell death by inhibiting NF-κB signaling

Highlights: • OGDR induces NKILA expression to inhibit NF-κB activation in neuronal cells. • NKILA shRNA attenuates OGDR-induced NF-κB inhibition and neuronal cell death. • NKILA overexpression augments OGDR-induced neuronal cell death. • miR-103/107 downregulation mediates NKILA upregulation by OGDR in neuronal cells. Oxygen glucose deprivation (OGD)/re-oxygenation (OGDR) induces severe injury to neuronal cells. The expression and potential function of NKILA (NF-κB Interacting LncRNA) in OGDR-treated neuronal cells were tested in this study. We show that OGDR induced NKILA upregulation to in-activate NF-κB signaling in SH-SY5Y cells and primary murine hippocampal neurons. Conversely, shRNA-mediated NKILA silencing almost reversed OGDR-induced NF-κB inhibition. OGDR-induced neuronal cell viability reduction, apoptosis and necrosis were largely attenuated by NKILA shRNA as well. Conversely, ectopic overexpression of NKILA by a lentiviral vector enhanced OGDR-induced SH-SY5Y cell death. For the mechanism study, we show that OGDR downregulated miR-103 and miR-107 to induce NKILA upregulation in neuronal cells. Transfection of miR-103 mimic or miR-107 mimic almost reversed OGDR-induced NKILA upregulation, NF-κB in-activation and SH-SY5Y cell death. Taken together, OGDR induces NKILA upregulation to in-activate NF-κB signaling, which mediates subsequent neuronal cell death. NKILA could be a novel therapeutic target of ischemic neuronal injury.

60 APPLIED LIFE SCIENCES↗

Constraining $\mathcal{f(R)}$ gravity with a $\mathcal{k}$-cut cosmic shear analysis of the Hyper Suprime-Cam first-year data

Here, using Subaru Hyper Suprime-Cam (HSC) year 1 data, we perform the first k-cut cosmic shear analysis constraining both Λ CDM and f(R) Hu-Sawicki modified gravity. To generate the f(R) cosmic shear theory vector, we use the matter power spectrum emulator trained on COLA (COmoving Lagrangian Acceleration) simulations [Phys. Rev. D 103, 123525 (2021)]. The k-cut method is used to significantly down-weight sensitivity to small scale (k > 1h Mpc -1 ) modes in the matter power spectrum where the emulator is less accurate, while simultaneously ensuring our results are robust to baryonic feedback model uncertainty. We have also developed a test to ensure that the effects of poorly modeled small scales are nulled as intended. For Λ CDM we find S 8 =σ 8 (Ω m /0.3) 0.5 =0.789$^{+0.039}_{-0.022}$, while the constraints on the f(R) modified gravity parameters are prior dominated. In the future, the k-cut method could be used to constrain a large number of theories of gravity where computational limitations make it infeasible to model the matter power spectrum down to extremely small scales.

79 ASTRONOMY AND ASTROPHYSICS↗

Forecasting longitudinal changes in oropharyngeal tumor morphology throughout the course of head and neck radiation therapy

Purpose: To create models that forecast longitudinal trends in changing tumor morphology and to evaluate and compare their predictive potential throughout the course of radiation therapy. Methods: Two morphology feature vectors were used to describe 35 gross tumor volumes (GTVs) throughout the course of intensity-modulated radiation therapy for oropharyngeal tumors. The feature vectors comprised the coordinates of the GTV centroids and a description of GTV shape using either interlandmark distances or a spherical harmonic decomposition of these distances. The change in the morphology feature vector observed at 33 time points throughout the course of treatment was described using static, linear, and mean models. Models were adjusted at 0, 1, 2, 3, or 5 different time points (adjustment points) to improve prediction accuracy. The potential of these models to forecast GTV morphology was evaluated using leave-one-out cross-validation, and the accuracy of the models was compared using Wilcoxon signed-rank tests. Results: Adding a single adjustment point to the static model without any adjustment points decreased the median error in forecasting the position of GTV surface landmarks by the largest amount (1.2 mm). Additional adjustment points further decreased the forecast error by about 0.4 mm each. Selection of the linear model decreased the forecast error for both the distance-based and spherical harmonic morphology descriptors (0.2 mm), while the mean model decreased the forecast error for the distance-based descriptor only (0.2 mm). The magnitude and statistical significance of these improvements decreased with each additional adjustment point, and the effect from model selection was not as large as that from adding the initial points. Conclusions: The authors present models that anticipate longitudinal changes in tumor morphology using various models and model adjustment schemes. The accuracy of these models depended on their form, and the utility of these models includes the characterization of patient-specific response with implications for treatment management and research study design.

60 APPLIED LIFE SCIENCES↗

TMD program at JLab

Measurements of multiplicities and asymmetries of multiparticle final states in wide kinematics, including dihadrons and vector mesons, will be crucial for a separation of different structure functions, and different contributions to specific structure functions, used in phenomenological studies of Semi-Inclusive DIS (SIDIS). Measurements of various beam Single Spin Asymmetries are discussed as a possible tool to separate different kinematical regions as well as different dynamical contributions. Precision measurements in JLab, including detailed studies of the $Q^2$ dependences of observables, would also allow testing the impact of several theoretical assumptions used in phenomenological studies of generalized PDFs (GPDs, TMDs) , also providing a mechanism for validation of the extraction frameworks, which is critical for the proper evaluation of systematic uncertainties.

Avagyan, Harutyun↗

Vector thrust induced lift effects for several ejector exhaust locations on a V/STOL wind tunnel model at forward speed

The results and analysis of aerodynamic force data obtained from a small scale model of a V/STOL research vehicle in a low speed wind tunnel are presented. The analysis of the data includes the evaluation of aerodynamic-propulsive lift performance when operating twin ejector nozzles with thrust deflected. Three different types of thrust deflector systems were examined: 90 deg downward deflected nozzle, 90 deg slotted nozzle with boundary layer control, and an externally blown flap configuration. Several nozzle locations were tested, including over and underwing positions. The interference lift of the nacelle and model due to jet exhaust thrust is compared and results show that 90 deg turned nozzles located over the wing (near the trailing edge) produce the largest interference lift increment for an untrimmed aircraft, and that the slotted nozzle located under the wing near the trailing edge (in conjunction with a BLC flap) gives a comparable interference lift in the trimmed condition. The externally blown flap nozzle produced the least interference lift and significantly less total lift due to jet thrust effects.

Sharon, A. D.↗

Improved definition of crustal magnetic anomalies for MAGSAT data

The routine correction of MAGSAT vector magnetometer data for external field effects such as the ring current and the daily variation by filtering long wavelength harmonics from the data is described. Separation of fields due to low altitude sources from those caused by high altitude sources is affected by means of dual harmonic expansions in the solution of Dirichlet's problem. This regression/harmonic filter procedure is applied on an orbit by orbit basis, and initial tests on MAGSAT data from orbit 1176 show reduction in external field residuals by 24.33 nT RMS in the horizontal component, and 10.95 nT RMS in the radial component.

Brown, R. D.↗

Predictor-Corrector Entry Guidance for Reusable Launch Vehicles

An online entry guidance algorithm has been developed using a predictor-corrector approach. The algorithm is designed for the Reusable Launch Vehicle (RLV) and is demonstrated by using, the X-33 model. The objective of the design is to handle widely dispersed entry conditions and deliver the vehicle at the Terminal Area Energy Management (TAEM) interface box within an acceptable tolerance and without violating any of the vehicle physical constraints. Combination of several control variables is used in testing the performance and computational requirement of the algorithm. The control variables are the bank angle, angle-of-attack and the time for roll reversal. The bank angle and angle-of-attack profiles are the nominal profiles plus the perturbations in each direction. The initial guess of the bank profile is a 45 degrees bank angle with reversal at 360 seconds from liftoff. A six-element state vector is propagated to the TAEM interface box through the integration of the equations of motion (EOM). Altitude, heading and range errors are computed between the desired and the achieved state at the TAEM interface. These errors are used to correct the initial guess of the control variables. This process is repeated until the errors meet an acceptable level at the TAEM interface. Several numerical optimization methods are used to evaluate the convergent property of the predictor-predictor methodology. Successful results are demonstrated using the X-33 model.

Youssef, Hussein↗