Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Local time stepping”

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 343 records · Page 19

Solar System Escape Architecture for Revolutionary Science

Our most distance spacecraft, Voyager 1 and Voyager 2, recently reached the boundary of the solar system known as the Heliopause. Beyond this boundary lies interstellar space, and, at tremendous distance, the stars. The Voyager spacecraft have shown us that the Heliopause is a dynamic, rapidly changing environment, varying greatly as the fields and effects emanating from the sun are met by a host of forces pushing in from our nearby galactic neighborhood. This is the region we aim to explore. Specifically we aim to develop a mission architecture and spacecraft concept capable of reaching the Heliopause (approximately100 AU (Astronomical Units)) region in multiple directions of interest, with a variety of scientific instrument suites, within a reasonable timeframe (about 15 years.) This report details our Phase 1 work including science goal definition, trajectory planning, technology research, mission planning, instrument selection, spacecraft design, and more. We present an architecture employing a 250-meter by 250-meter solar sail with a 175 kilogram spacecraft. This spacecraft could carry a variety of instrument suites depending on the destination, time, and other factors. The spin-stabilized spacecraft would powered by a small radioisotope power system (about 20 watts), makes use of an optical communication system, and carries 3 extensible booms for instrument accommodation as visible in the sketch and rendering below. We conclude that it would be reasonable to implement the architecture described herein in the 2020-2030 timeframe based on likely near-term technological and material advances. We are confident that cruise times of around 15 years could be achieved to science-rich destinations. This report details the approach we used to come to these conclusions and identifies some of the steps along the way. We utilized a variety of methodologies involving our full capability at JPL (Jet Propulsion Laboratory) - from small point studies to large collaborative engineering processes. This paper also describes potential industry partnerships and costing information, and includes the methods used to account for our project’s significant deviation from any prior mission in the NASA costing tool database. Outreach and public engagement is an important part of NASA activities. We engaged in several different forms of this outreach and describe them in this report. We conclude with some open questions and a summary of activities we hope to engage in with Phase 2 funding. The entire team would like to thank the NIAC (NASA Innovative Advanced Concepts) program for the opportunity to investigate this fascinating and challenging concept. The Voyager program gives us tantalizing clues as to the nature and behavior of the Heliopause – the boundary between the sun’s influence and the interstellar medium. This information comes from forty-year old instruments designed to study the outer planets. A targeted Heliopause investigation would give insight into the formation of the solar system, the role of the sun in the local interstellar neighborhood, and contribute to human exploration planning by helping to predict periods of low galactic cosmic ray (GCR) penetration into the inner solar system.

Spacecraft↗

Estimating Resolution Lengths of Hybrid Turbulence Models

A two-stage procedure has been devised for estimating the spatial resolution achievable in the simulation of a given flow on a given computational grid by a computational fluid dynamics (CFD) code that incorporates a hybrid model of turbulence. The hybrid models to which this procedure is especially relevant are those of the Reynolds-averaged Navier-Stokes (RANS) and the partial-averaged Navier-Stokes (PANS) approaches. This procedure represents the first step toward adding variable-resolution turbulence-modeling capabilities to CFD codes as part of a continuing effort to increase the accuracy and robustness of CFD simulations of unsteady flows. Some background information is prerequisite to a meaningful summary of the procedure. Among experts in CFD, it is well known that combination of the Reynolds-averaged Navier-Stokes (RANS) approach and eddy-viscosity turbulence models offers limited capability for simulating unsteady and complex flows. The RANS approach includes an assumption that most of the energy in a given flow is modeled through turbulence-transport equations and is resolved in a computational grid used to simulate the flow. RANS also overpredicts eddy viscosity, thereby yielding excessive damping of unsteady motion. The eddy viscosity attains an unphysically large value because of unresolved scales, and suppresses most temporal and spatial fluctuations in the resolved flow field. One approach used to overcome this deficiency is to provide a mechanism for the RANS equations to resolve motion only on the largest scales and to use a hybrid model to represent effects at smaller scales. The RANS approach involves the use of a standard two-equation turbulence model in which the effect of turbulence is summarized by a viscosity that is a function of (1) the time-averaged kinetic- energy density (k) associated with the local fluctuating (turbulent) component of flow and (2) the time-averaged rate of dissipation of the turbulent-kinetic- energy density ( ). In PANS, which was developed to overcome the grid dependency associated with other hybrid turbulence models (including that of RANS), the standard two-equation turbulence model is replaced by another two-equation model in which one solves for the previously unresolved k and , which are now allowed to vary in space and time.

Abdol-Hamid, Khaled S.↗

Dynamics of Liquids in Edges and Corners (DYLCO): IML-2 Experiment for the BDPU

Knowledge of the behavior of fluids possessing free surfaces is important to many fluid systems, particularly in space, where the normally subtle effects of surface wettability play a more dramatic and often surprising role. DYLCO for the IML-2 mission was proposed as a simple experiment to probe the particular behavior of capillary surfaces in containers of irregular cross section. Temperature control was utilized to vary the fluid-solid contact angle, a questionable thermodynamic parameter of the system, small changes in which can dramatically influence the configuration, stability, and flow of a capillary surface. Container shapes, test fluid, and temperature ranges were selected for observing both local changes in interface curvature as well as a global change in fluid orientation due to a critical wetting phenomenon. The experiment hardware performed beyond what was expected and fluid interfaces could be readily digitized post flight to show the dependence of the interface curvature on temperature. For each of the containers tested surfaces were observed which did not satisfy the classic equations for the prediction of interface shape with constant contact angle boundary condition. This is explained by the presence of contact angle hysteresis arising from expansion and contraction of the liquid during the heating and cooling steps of the test procedure. More importantly, surfaces exceeding the critical surface curvature required for critical wetting were measured, yet no wetting was observed. These findings are indeed curious and pose key questions concerning the role of hysteresis for this critical wetting phenomena. The stability of such surfaces was determined numerically and it is shown that stability is enhance (reduced) when a surface is in its 'advancing' ('receding') state, The analysis shows complete instability as the critical wetting condition is reached. The case of ideal dynamic wetting is addressed analytically in detail with results of significant flow characteristics presented in closed form. The solutions indicate a square root of T dependence of the capillary 'rise' rate which is corroborated by drop tower tests. The analysis clearly shows that infinite time is necessary for surfaces to reorient at the critical wetting transition.

Langbein, D.↗

Convective flow effects on protein crystal growth

The long-term stability of the interferometric setup for the monitoring of protein morphologies has been improved. Growth or dissolution of a crystal on a 100 A scale can now be clearly distinguished from dimensional changes occurring within the optical path of the interferometer. This capability of simultaneously monitoring the local interfacial displacement at several widely-spaced positions on the crystal surface with high local depth resolution, has already yielded novel results. We found with lysozyme that (1) the normal growth rate is oscillatory, and (2) the mean growth step density is greater at the periphery of a facet than in its center. The repartitioning of Na(+) and Cl(-) ions between lysozyme solutions and crystals was studied for a wide range of crystallization conditions. A nucleation-growth-repartitioning model was developed to interpret the large body of data in a unified way. The results strongly suggests that (1) the ion to lysozyme ratio in the crystal depends mostly on kinetic rather than crystallographic parameters, and (2) lysozyme crystals possess a salt-rich core with a diameter on the order of 10 microns. The computational model for diffusive-convective transport in protein crystallization (see the First Report) has been applied to a realistic growth cell geometry, taking into account the findings of the above repartitioning studies. These results show that some elements of a moving boundary problem must be incorporated into the model in order to obtain a more realistic description. Our experimental setup for light scattering investigations of aggregation and nucleation in protein solutions has been extensively tested. Scattering intensity measurements with a true Rayleigh scatterer produced systematically increased forward scattering, indicating problems with glare. These have been resolved. Preliminary measurements with supersaturated lysozyme solutions revealed that the scatterers grow with time. Work has begun on a computer program for the unified evaluation of simultaneously obtained, multi-angle static and dynamic light scattering data.

Rosenberger, Franz↗

The Quantum Approximation Optimization Algorithm for MaxCut: A Fermionic View

Farhi et al. recently proposed a class of quantum algorithms, the Quantum Approximate Optimization Algorithm (QAOA), for approximately solving combinatorial optimization problems. A level-p QAOA circuit consists of steps in which a classical Hamiltonian, derived from the cost function, is applied followed by a mixing Hamiltonian. The 2p times for which these two Hamiltonians are applied are the parameters of the algorithm. As p increases, however, the parameter search space grows quickly. The success of the QAOA approach will depend, in part, on finding effective parameter-setting strategies. Here, we analytically and numerically study parameter setting for QAOA applied to MAXCUT. For level-1 QAOA, we derive an analytical expression for a general graph. In principle, expressions for higher p could be derived, but the number of terms quickly becomes prohibitive. For a special case of MAXCUT, the Ring of Disagrees, or the 1D antiferromagnetic ring, we provide an analysis for arbitrarily high level. Using a Fermionic representation, the evolution of the system under QAOA translates into quantum optimal control of an ensemble of independent spins. This treatment enables us to obtain analytical expressions for the performance of QAOA for any p. It also greatly simplifies numerical search for the optimal values of the parameters. By exploring symmetries, we identify a lower-dimensional sub-manifold of interest; the search effort can be accordingly reduced. This analysis also explains an observed symmetry in the optimal parameter values. Further, we numerically investigate the parameter landscape and show that it is a simple one in the sense of having no local optima.

quantum algorithm↗

Adjoint Techniques for Topology Optimization of Structures Under Damage Conditions

The objective of this cooperative agreement was to seek computationally efficient ways to optimize aerospace structures subject to damage tolerance criteria. Optimization was to involve sizing as well as topology optimization. The work was done in collaboration with Steve Scotti, Chauncey Wu and Joanne Walsh at the NASA Langley Research Center. Computation of constraint sensitivity is normally the most time-consuming step of an optimization procedure. The cooperative work first focused on this issue and implemented the adjoint method of sensitivity computation (Haftka and Gurdal, 1992) in an optimization code (runstream) written in Engineering Analysis Language (EAL). The method was implemented both for bar and plate elements including buckling sensitivity for the latter. Lumping of constraints was investigated as a means to reduce the computational cost. Adjoint sensitivity computation was developed and implemented for lumped stress and buckling constraints. Cost of the direct method and the adjoint method was compared for various structures with and without lumping. The results were reported in two papers (Akgun et al., 1998a and 1999). It is desirable to optimize topology of an aerospace structure subject to a large number of damage scenarios so that a damage tolerant structure is obtained. Including damage scenarios in the design procedure is critical in order to avoid large mass penalties at later stages (Haftka et al., 1983). A common method for topology optimization is that of compliance minimization (Bendsoe, 1995) which has not been used for damage tolerant design. In the present work, topology optimization is treated as a conventional problem aiming to minimize the weight subject to stress constraints. Multiple damage configurations (scenarios) are considered. Each configuration has its own structural stiffness matrix and, normally, requires factoring of the matrix and solution of the system of equations. Damage that is expected to be tolerated is local and represents a small change in the stiffness matrix compared to the baseline (undamaged) structure. The exact solution to a slightly modified set of equations can be obtained from the baseline solution economically without actually solving the modified system.. Shennan-Morrison-Woodbury (SMW) formulas are matrix update formulas that allow this (Akgun et al., 1998b). SMW formulas were therefore used here to compute adjoint displacements for sensitivity computation and structural displacements in damaged configurations.

Akgun, Mehmet A.↗

Channel Wall Nozzle Manufacturing Technology Advancements for Liquid Rocket Engines

A regeneratively-cooled nozzle is a critical component for expansion of the hot gases to enable high temperature and performance liquid rocket engines systems. Channel wall nozzles are a design solution used across the propulsion industry as a simplified method to fabricate the nozzle structure with internal coolant passages. The scale and complexity of the channel wall nozzle (CWN) design is challenging to fabricate leading to extended lead times and higher costs. Some of these challenges include: 1) unique and high temperature materials, 2) Tight tolerances during manufacturing and assembly to contain high pressure propellants, 3) thin-walled features to maintain adequate wall temperatures, and 4) Unique manufacturing process operations and tooling. The United States (U.S.) National Aeronautics and Space Administration (NASA) along with U.S. specialty manufacturing vendors are maturing modern fabrication techniques to reduce complexity and decrease costs associated with channel wall nozzle manufacturing technology. Additive Manufacturing (AM) is one of the key technology advancements being evaluated for channel wall nozzles. Much of additive manufacturing for propulsion components has focused on powder bed fusion, but the scale is not yet feasible for application to large scale nozzles. NASA is evolving directed energy deposition (DED) techniques for nozzles including arc-based deposition, blown powder deposition, and Laser Wire Direct Closeout (LWDC). There are different approaches being considered for fabrication of the nozzle and each of these DED processes offer unique process steps for rapid fabrication. The arc-based and blown powder deposition techniques are being used for the forming of the CWN liner. A variety of materials are being demonstrated including Inconel 625, Haynes 230, JBK-75, and NASA HR-1. The blown powder DED process is also being demonstrated for forming an integral channel nozzle in a single operation in similar materials. The LWDC process is a method for closing out the channels within the liner and forming the structural jacket using a localized laser wire deposition technique. Identical materials mentioned above have been used for this process in addition to bimetallic closeout (C-18150–SS347, and C-18150–Inconel 625). NASA has completed process development, material characterization, and hot-fire testing on a variety of these channel wall nozzle fabrication techniques. This paper will present an overview of the various processes and materials being evaluated and the results from the hot-fire testing. Future development and technology focus areas will also be discussed relative to channel wall nozzle manufacturing.

Gradl, Paul R.↗

Technology Advancements for Channel Wall Nozzle Manufacturing In Liquid Rocket Engines

A regeneratively-cooled or dump-cooled nozzle is a critical component for expansion of hot gases to enable high temperature and performance in liquid rocket engines systems. Regeneratively-cooled channel wall nozzles are a design solution used across the propulsion industry as a simplified method to fabricate the nozzle structure with internal coolant passages. The scale and complexity of the channel wall nozzle (CWN) design can be challenging to fabricate which results in extended lead times and higher costs. Some of these challenges include: 1) Unique and high temperature materials, 2) Tight tolerances on large parts during manufacturing and assembly to contain high pressure propellants, 3) Thin-walled features to maintain adequate wall temperatures, and 4) Unique manufacturing process operations and complex tooling. The United States (U.S.) National Aeronautics and Space Administration (NASA) and U.S. specialty manufacturing vendors are maturing modern fabrication techniques to reduce complexity and decrease costs associated with channel wall nozzle manufacturing technology. Additive Manufacturing (AM) is one of the key technology advancements under evaluation for channel wall nozzles. Much of additive manufacturing for propulsion components has focused on laser powder bed fusion (L-PBF), but the scale is not yet feasible for application to large scale nozzles. NASA is evolving directed energy deposition (DED) techniques for nozzles including arc-based deposition, blown powder deposition, and Laser Wire Direct Closeout (LWDC). There are different approaches being considered for fabrication of the nozzle, and each of these DED processes offer unique process steps for rapid fabrication. The arc-based and blown powder deposition techniques are used for the forming of the CWN liner. A variety of materials are being demonstrated including Inconel 625, Haynes 230, JBK-75, and NASA HR-1. The blown powder DED process is also being demonstrated for forming an integral channel nozzle in a single operation in similar materials. The LWDC process is a method for closing out the channels within the liner and forming the structural jacket using a localized laser wire deposition technique. Identical materials mentioned above have been used for this process in addition to bimetallic closeout (C-18150–SS347, and C-18150–Inconel 625). NASA has completed process development, material characterization, and hot-fire testing on a variety of these channel wall nozzle fabrication technique. This publication presents an overview of the various channel wall nozzle manufacturing processes and materials under evaluation including results from the hot-fire testing. Future development and technology focus areas is also discussed relative to channel wall nozzle manufacturing.

Additive manufacturing↗

Design Criteria for X-CRV Honeycomb Panels: A Preliminary Study

The objective of this project is to perform the first step in developing structural design criteria for composite sandwich panels that are to be used in the aeroshell of the crew return vehicle (X-CRV). The preliminary concept includes a simplified method for assessing the allowable strength in the laminate material. Ultimately, it is intended that the design criteria be extended to address the global response of the vehicle. This task will require execution of a test program as outlined in the recommendation section of this report. The aeroshell of the X-CRV is comprised of composite sandwich panels consisting of fiberite face sheets and a phenolic honeycomb core. The function of the crew return vehicle is to enable the safe return of injured or ill crewpersons from space station, the evacuation of crew in case of emergency or the return of crew if an orbiter is not available. A significant objective of the X-CRV project is to demonstrate that this vehicle can be designed, built and operated at lower cost and at a significantly faster development time. Development time can be reduced by driving out issues in both structural design and manufacturing concurrently. This means that structural design and analysis progresses in conjunction with manufacturing and testing. Preliminary tests results on laminate coupons are presented in the report. Based on these results a method for detection material failure in the material is presented. In the long term, extrapolation of coupon data to large scale structures may be inadequate. Test coupons used to develop failure criteria at the material scale are typically small when compared to the overall structure. Their inherent small size indicates that the material failure criteria can be used to predict localized failure of the structure, however, it can not be used to predict failure for all failure modes. Some failure modes occur only when the structure or one of its sub-components are studied as a whole. Conversely, localized failure may not indicate failure of the structure as a whole and the amount of reserve capacity, if any, should be assessed. To develop a complete design criteria experimental studies of the sandwich panel are needed. Only then can a conservative and accurate design criteria be developed. This criteria should include effects of flaws and defects, and environmental factors such as temperature and moisture. Preliminary results presented in this report suggest that a simplified analysis can be used to predict the strength of a laminate. Testing for environmental effects have yet to be included in this work. The so called 'rogue flaw test' appears to be a promising method for assessing the effect of a defect in a laminate. This method fits in quite well with the philosophy of achieving a damage tolerant design.

Caccese, Vincent↗

Configuration and Evaluation of a Dual-Doppler 3-D Wind Field System

Current LSP, GSDO, and SLS space vehicle operations are halted when wind speeds from specific directions exceed defined thresholds and when lightning is a threat. Strong winds and lightning are difficult parameters for the 45th Weather Squadron (45 WS) to forecast, yet are important in the protection of customer vehicle operations and the personnel that conduct them. A display of the low-level horizontal wind field to reveal areas of high winds or convergence would be a valuable tool for forecasters in assessing the timing of high winds, or convection initiation and subsequent lightning occurrence. This is especially important for areas where no weather observation platforms exist. Developing a dual-Doppler radar capability would provide such a display to assist forecasters in predicting high winds and convection initiation. The wind fields can also be used to initialize a local mesoscale numerical weather prediction model to help improve the model forecast winds, convection initiation, and other phenomena. The 45 WS and NWS MLB tasked the Applied Meteorology Unit (AMU) to develop a dual- Doppler wind field display using data from the 45th Space Wing radar, known as the Weather Surveillance Radar (WSR), NWS MLB Weather Surveillance Radar 1988 Doppler (KMLB), and the Orlando International Airport Terminal Doppler Weather Radar (KMCO). They also stipulated that the software used should be freely available. The AMU evaluated two software packages and, with concurrence from NWS MLB and the 45 WS, chose the Warning Decision Support System-Integrated Information (WDSS-II). The AMU collected data from two significant weather cases: a tornadic event on 14 April 2013 and a severe wind and hail event on 12 February 2014. For the 14 April case, the data were from WSR and KMLB. For the 12 February case, the data were from KMCO and KMLB. The AMU installed WDSS-II on a Linux PC, then processed and quality controlled the radar data for display and analysis using WDSS-II tools. Because of issues with de-aliasing the WSR velocity field, the AMU did not use data from this radar in this study and only analyzed the 12 February case. Merging the data to create the dual-Doppler analysis involved several steps. The AMU used instructions from the WDSS-II website and discussion forum to determine the correct tools to use for the analysis, and was successful in creating a merged reflectivity field, which was critical to the success of creating a merged velocity field. However, the AMU was unable to create a merged velocity field. The AMU researched the WDSS-II forum for discussions on similar issues, asked questions on the forum, and tested different options and values in the merger tool with no success. Developing a dual-Doppler wind field was the main goal of this task, but that was not accomplished. It could be an issue of not using the correct options or the correct value for the options used, or there could be issues with the radar data. There is a follow-on AMU task to install the operational version of WDSS-II in the NWS MLB office. This will provide more opportunities to try different options and input values in order to create a merged wind field from KMCO and KMLB.

Doppler radar↗

Machine Learning Techniques in Optimal Design

Many important applications can be formalized as constrained optimization tasks. For example, we are studying the engineering domain of two-dimensional (2-D) structural design. In this task, the goal is to design a structure of minimum weight that bears a set of loads. A solution to a design problem in which there is a single load (L) and two stationary support points (S1 and S2) consists of four members, E1, E2, E3, and E4 that connect the load to the support points is discussed. In principle, optimal solutions to problems of this kind can be found by numerical optimization techniques. However, in practice [Vanderplaats, 1984] these methods are slow and they can produce different local solutions whose quality (ratio to the global optimum) varies with the choice of starting points. Hence, their applicability to real-world problems is severely restricted. To overcome these limitations, we propose to augment numerical optimization by first performing a symbolic compilation stage to produce: (a) objective functions that are faster to evaluate and that depend less on the choice of the starting point and (b) selection rules that associate problem instances to a set of recommended solutions. These goals are accomplished by successive specializations of the problem class and of the associated objective functions. In the end, this process reduces the problem to a collection of independent functions that are fast to evaluate, that can be differentiated symbolically, and that represent smaller regions of the overall search space. However, the specialization process can produce a large number of sub-problems. This is overcome by deriving inductively selection rules which associate problems to small sets of specialized independent sub-problems. Each set of candidate solutions is chosen to minimize a cost function which expresses the tradeoff between the quality of the solution that can be obtained from the sub-problem and the time it takes to produce it. The overall solution to the problem, is then obtained by solving in parallel each of the sub-problems in the set and computing the one with the minimum cost. In addition to speeding up the optimization process, our use of learning methods also relieves the expert from the burden of identifying rules that exactly pinpoint optimal candidate sub-problems. In real engineering tasks it is usually too costly to the engineers to derive such rules. Therefore, this paper also contributes to a further step towards the solution of the knowledge acquisition bottleneck [Feigenbaum, 1977] which has somewhat impaired the construction of rulebased expert systems.

Cerbone, Giuseppe↗

AC and DC Fault Management for Megawatt Electrified Aircraft Electrical Powertrains Task 3: Lifetime and Reliability of Electrical Insulators

This research project was a collaborative investigation between researchers at the RTX Technology Research Center (RTRC) and the University of Texas at Austin and made a significant contribution to enabling electric aircraft. The transport of electric power between the points of generation and use requires power cables. These cables must be smaller, lighter and provide a more predictable life than power cables used in stationary applications. Consequently, this investigation provided heretofore unavailable information supporting the safety and reliability of smaller lighter power cables for electrified aircraft. In addition, the research identified key additional engineering data needed to support quantitative reliability assessments. Important advances included: • Demonstrated that at least one manufacturer can make a novel, smaller, lighter power cable that is free from serious defects. • Developed and published an appropriate analytical construct to describe the life of this novel cable. This is a necessary step for use in aviation where the understanding of remaining life is critical. • Demonstrated thermal-mechanical aging that suggested 1000+ flights before the thermal-mechanical processes produced defects large enough that the defect growth was accelerated electrically. • Showed that electrical aging took place at two rates. The first possibly lasting weeks to months and the second possibly days to weeks. If robust, this provides a good diagnostic for cable replacement. • Demonstrated that the traditional electrical testing of cable materials using manufactured voids can be misleading due to the size of the voids. Emerging laser drilling technology permitted demonstration that the physics of failure in realistically small voids is different from that in the unrealistically large voids used in earlier research, which is very important for high-quality, high-performance, small aircraft cables. Although this project represents a significant contribution to the specifics of cable aging in the aircraft environment, important additional research remains to be completed, including: • Non-uniform thermal cycling by applying the heat from the center conductor to maximize thermal stress next to the core area where the electric gradient is the strongest. This builds on the uniform thermal cycling that has been completed. • The augmentation of the thermal-mechanical failure rate by electrical processes. Better understanding of these time constants strongly affects the ability to predict life. • Termination design: Terminations provide not only electrical reflection potential, but a location for a series arc fault and an area where ozone can diffuse into the center conductor and negatively affect cable insulation. • The abrasion and ozone resistance of the cable jacket. • Pressure cycling as an accelerant of thermal, mechanical, and/or electrical aging. • Possible methods for online PD detection and offline PD localization

model↗

Analysis of Aircraft Clusters to Measure Sector-Independent Airspace Congestion

The Distributed Air/Ground Traffic Management (DAG-TM) concept of operations* permits appropriately equipped aircraft to conduct Free Maneuvering operations. These independent aircraft have the freedom to optimize their trajectories in real time according to user preferences; however, they also take on the responsibility to separate themselves from other aircraft while conforming to any local Traffic Flow Management (TFM) constraints imposed by the air traffic service provider (ATSP). Examples of local-TFM constraints include temporal constraints such as a required time of arrival (RTA), as well as spatial constraints such as regions of convective weather, special use airspace, and congested airspace. Under current operations, congested airspace typically refers to a sector(s) that cannot accept additional aircraft due to controller workload limitations; hence Dynamic Density (a metric that is indicative of controller workload) can be used to quantify airspace congestion. However, for Free Maneuvering operations under DAG-TM, an additional metric is needed to quantify the airspace congestion problem from the perspective of independent aircraft. Such a metric would enable the ATSP to prevent independent aircraft from entering any local areas of congestion in which the flight deck based systems and procedures may not be able to ensure separation. This new metric, called Gaggle Density, offers the ATSP a mode of control to regulate normal operations and to ensure safety and stability during rare-normal or off-normal situations (e.g., system failures). It may be difficult to certify Free Maneuvering systems for unrestricted operations, but it may be easier to certify systems and procedures for specified levels of Gaggle Density that could be monitored by the ATSP, and maintained through relatively minor flow-rate (RTA type) restrictions. Since flight deck based separation assurance is airspace independent, the challenge is to measure congestion independent of sector boundaries. Figure 1 , reproduced from Ref. 1, depicts an example traffic situation. When the situation is analyzed by sector boundaries (left side of figure), a Dynamic Density metric would identify excessive congestion in the central sector. When the same traffic situation is analyzed independent of sector boundaries (right side of figure), a Gaggle Density metric would identify congestion in two dynamically defined areas covering portions of several sectors. The first step towards measuring airspace-independent congestion is to identify aircraft clusters, i.e., groups of closely spaced aircraft. The objective of this work is to develop techniques to detect and classify clusters of aircraft.

Bilimoria, Karl D.↗

Application of Simulated Annealing and Related Algorithms to TWTA Design

Simulated Annealing (SA) is a stochastic optimization algorithm used to search for global minima in complex design surfaces where exhaustive searches are not computationally feasible. The algorithm is derived by simulating the annealing process, whereby a solid is heated to a liquid state and then cooled slowly to reach thermodynamic equilibrium at each temperature. The idea is that atoms in the solid continually bond and re-bond at various quantum energy levels, and with sufficient cooling time they will rearrange at the minimum energy state to form a perfect crystal. The distribution of energy levels is given by the Boltzmann distribution: as temperature drops, the probability of the presence of high-energy bonds decreases. In searching for an optimal design, local minima and discontinuities are often present in a design surface. SA presents a distinct advantage over other optimization algorithms in its ability to escape from these local minima. Just as high-energy atomic configurations are visited in the actual annealing process in order to eventually reach the minimum energy state, in SA highly non-optimal configurations are visited in order to find otherwise inaccessible global minima. The SA algorithm produces a Markov chain of points in the design space at each temperature, with a monotonically decreasing temperature. A random point is started upon, and the objective function is evaluated at that point. A stochastic perturbation is then made to the parameters of the point to arrive at a proposed new point in the design space, at which the objection function is evaluated as well. If the change in objective function values (Delta)E is negative, the proposed new point is accepted. If (Delta)E is positive, the proposed new point is accepted according to the Metropolis criterion: rho((Delta)f) = exp((-Delta)E/T), where T is the temperature for the current Markov chain. The process then repeats for the remainder of the Markov chain, after which the temperature is decremented and the process repeats. Eventually (and hopefully), a near-globally optimal solution is attained as T approaches zero. Several exciting variants of SA have recently emerged, including Discrete-State Simulated Annealing (DSSA) and Simulated Tempering (ST). The DSSA algorithm takes the thermodynamic analogy one step further by categorizing objective function evaluations into discrete states. In doing so, many of the case-specific problems associated with fine-tuning the SA algorithm can be avoided; for example, theoretical approximations for the initial and final temperature can be derived independently of the case. In this manner, DSSA provides a scheme that is more robust with respect to widely differing design surfaces. ST differs from SA in that the temperature T becomes an additional random variable in the optimization. The system is also kept in equilibrium as the temperature changes, as opposed to the system being driven out of equilibrium as temperature changes in SA. ST is designed to overcome obstacles in design surfaces where numerous local minima are separated by high barriers. These algorithms are incorporated into the optimal design of the traveling-wave tube amplifier (TWTA). The area under scrutiny is the collector, in which it would be ideal to use negative potential to decelerate the spent electron beam to zero kinetic energy just as it reaches the collector surface. In reality this is not plausible due to a number of physical limitations, including repulsion and differing levels of kinetic energy among individual electrons. Instead, the collector is designed with multiple stages depressed below ground potential. The design of this multiple-stage collector is the optimization problem of interest. One remaining problem in SA and DSSA is the difficulty in determining when equilibrium has been reached so that the current Markov chain can be terminated. It has been suggested in recent literature that simulating the thermodynamic properties opecific heat, entropy, and internal energy from the Boltzmann distribution can provide good indicators of having reached equilibrium at a certain temperature. These properties are tested for their efficacy and implemented in SA and DSSA code with respect to TWTA collector optimization.

Radke, Eric M.↗

Large Scale Turbulent Structures in Supersonic Jets

Jet noise is a major concern in the design of commercial aircraft. Studies by various researchers suggest that aerodynamic noise is a major contributor to jet noise. Some of these studies indicate that most of the aerodynamic jet noise due to turbulent mixing occurs when there is a rapid variation in turbulent structure, i.e. rapidly growing or decaying vortices. The objective of this research was to simulate a compressible round jet to study the non-linear evolution of vortices and the resulting acoustic radiations. In particular, to understand the effect of turbulence structure on the noise. An ideal technique to study this problem is Direct Numerical Simulations(DNS), because it provides precise control on the initial and boundary conditions that lead to the turbulent structures studied. It also provides complete 3-dimensional time dependent data. Since the dynamics of a temporally evolving jet are not greatly different from those, of a spatially evolving jet, a temporal jet problem was solved, using periodicity ill the direction of the jet axis. This enables the application of Fourier spectral methods in the streamwise direction. Physically this means that turbulent structures in the jet are repeated in successive downstream cells instead of being gradually modified downstream into a jet plume. The DNS jet simulation helps us understand the various turbulent scales and mechanisms of turbulence generation in the evolution of a compressible round jet. These accurate flow solutions will be used in future research to estimate near-field acoustic radiation by computing the total outward flux across a surface and determine how it is related to the evolution of the turbulent solutions. Furthermore, these simulations allow us to investigate the sensitivity of acoustic radiations to inlet/boundary conditions, with possible application to active noise suppression. In addition, the data generated can be used to compute various turbulence quantities such as mean velocities, turbulent stresses, etc. which will aid in turbulence modeling. This report will be presented in two chapters. The first chapter describes some work on the linear stability of a supersonic round jet and the implications of this for the jet noise problem. The second chapter is an extensive discussion of numerical work using the spectral method which we use to solve the compressible Navier-Stokes equations to study turbulent jet flows. The method uses Fourier expansions in the azimuthal and streamwise direction and a 1-D B-spline basis representation in the radial direction. The B-spline basis is locally supported and this ensures block diagonal matrix equations which can be solved in O(N) steps. This is a modification of a boundary layer code developed by Robert Moser. A very accurate highly resolved Direct Numerical Simulation (DNS) of a turbulent jet flow is produced.

Rao, Ram Mohan↗

Large Scale Turbulent Structures in Supersonic Jets

Jet noise is a major concern in the design of commercial aircraft. Studies by various researchers suggest that aerodynamic noise is a major contributor to jet noise. Some of these studies indicate that most of the aerodynamic jet noise due to turbulent mixing occurs when there is a rapid variation in turbulent structure, i.e. rapidly growing or decaying vortices. The objective of this research was to simulate a compressible round jet to study the non-linear evolution of vortices and the resulting acoustic radiations. In particular, to understand the effect of turbulence structure on the noise. An ideal technique to study this problem is Direct Numerical Simulations (DNS), because it provides precise control on the initial and boundary conditions that lead to the turbulent structures studied. It also provides complete 3-dimensional time dependent data. Since the dynamics of a temporally evolving jet are not greatly different from those of a spatially evolving jet, a temporal jet problem was solved, using periodicity in the direction of the jet axis. This enables the application of Fourier spectral methods in the streamwise direction. Physically this means that turbulent structures in the jet are repeated in successive downstream cells instead of being gradually modified downstream into a jet plume. The DNS jet simulation helps us understand the various turbulent scales and mechanisms of turbulence generation in the evolution of a compressible round jet. These accurate flow solutions will be used in future research to estimate near-field acoustic radiation by computing the total outward flux across a surface and determine how it is related to the evolution of the turbulent solutions. Furthermore, these simulations allow us to investigate the sensitivity of acoustic radiations to inlet/boundary conditions, with possible appli(,a- tion to active noise suppression. In addition, the data generated can be used to compute, various turbulence quantities such as mean velocities, turbulent stresses, etc. which will aid in turbulence modeling. This report will be presented in two chapters. The first chapter describes some work on the linear stability of a supersonic round jet and the implications of this for the jet noise problem. The second chapter is an extensive discussion of numerical work using the spectral method which we use to solve the compressible Navier-Stokes equations to study turbulent jet flows. The method uses Fourier expansions in the azimuthal and streamwise direction and a 1-D B-spline basis representation in the radial direction. The B-spline basis is locally supported and this ensures block diagonal matrix equations which can be solved in O(N) steps. This is a modification of a boundary layer code developed by Robert Moser. A very accurate highly resolved DNS of a turbulent jet flow is produced.

Rao, Ram Mohan↗

Diamond Heat-Spreader for Submillimeter-Wave Frequency Multipliers

The planar GaAs Shottky diode frequency multiplier is a critical technology for the local oscillator (LO) for submillimeter- wave heterodyne receivers due to low mass, tenability, long lifetime, and room-temperature operation. The use of a W-band (75-100 GHz) power amplifier followed by a frequency multiplier is the most common for submillimeter-wave sources. Its greatest challenge is to provide enough input power to the LO for instruments onboard future planetary missions. Recently, JPL produced 800 mW at 92.5 GHz by combining four MMICs in parallel in a balanced configuration. As more power at W-band is available to the multipliers, their power-handling capability be comes more important. High operating temperatures can lead to degradation of conversion efficiency or catastrophic failure. The goal of this innovation is to reduce the thermal resistance by attaching diamond film as a heat-spreader on the backside of multipliers to improve their power-handling capability. Polycrystalline diamond is deposited by hot-filament chemical vapor deposition (CVD). This diamond film acts as a heat-spreader to both the existing 250- and 300-GHz triplers, and has a high thermal conductivity (1,000-1,200 W/mK). It is approximately 2.5 times greater than copper (401 W/mK) and 20 times greater than GaAs (46 W/mK). It is an electrical insulator (resistivity approx. equals 10(exp 15) Ohms-cm), and has a low relative dielectric constant of 5.7. Diamond heat-spreaders reduce by at least 200 C at 250 mW of input power, compared to the tripler without diamond, according to thermal simulation. This superior thermal management provides a 100-percent increase in power-handling capability. For example, with this innovation, 40-mW output power has been achieved from a 250-GHz tripler at 350-mW input power, while the previous triplers, without diamond, suffered catastrophic failures. This breakthrough provides a stepping-stone for frequency multipliers-based LO up to 3 THz. The future work for this design is to apply the high output power from both the 250 and 300 GHz to multiple chains in order to generate milliwatts at 2.3 THz. Using the first generation of results for this innovation, 40 mW of output power were produced from a 240-GHz tripler at 350-mW input power, and 27- mW output power was produced from a 300-GHz tripler at 408-mW input power. This is two times higher than the current state-of-the-art output power capability. A finite-element thermal simulation also shows that 30-microns thick diamond dropped the temperature of the anodes by at least 200 C.

Lin, Robert H.↗

The IOAG Recommendations on Spacecraft Emergency Cross Support

In 2014, the Inter-agency Operations Advisory Group (IOAG) chartered a multi-agency team effort to study how to best handle spacecraft emergency as part of cross support among network assets. The original intent is to put in place for the first time a process and guidelines to make emergency support as part of a permanent cross-support capability between space agencies. From the perspective of space communications service providers, a few key issues, common to all agencies, concerning the Spacecraft Emergency Cross Support (SECS) have been explored. They are: context of emergency support, provision of emergency support under a cross-support agreement, provision of emergency support with no cross-support agreement, legal and liability issues, response time, support priority, services provided, sustaining emergency support capabilities, and charge for the support. The positions of the various participating agencies with respect to these issues have been collected, analyzed, and finally harmonized to form the recommended IOAG positions. A central challenge most communications service providers are facing is that since the spacecraft emergency is an unplanned critical event, it typically requires fast response to the emergency call, hence lacking an international standard process for the operational interfaces seems to exacerbate the difficulty in providing SECS. For reducing the response time, i.e. from the time of accepting a request for SECS to the readiness for support, it is recommended that a “cross-support emergency system” be established by IOAG member agencies. Along with it, the IOAG core services, just-in-time ground communications line, cross support service management (CSSM), and standard operations procedures (SOP) for operational interfaces form the basic foundation of the “cross support emergency system”. Of the above, new to the cross support conducted thus far in the IOAG community is the concept of the SOP specifically for the interfaces between the service provider and service user during the SECS. Use cases, salient features, and definition of the key operational activities/tasks relevant to the interfaces are addressed by the effort. Underpinning such a system is the availability of the RF license granted by the local authority, at national and/or regional level, for a given ground station to communicate with and track the spacecraft in emergency mode at the uplink and downlink frequencies assigned to that spacecraft. That means it is critical for the IOAG member agencies to obtain a priori all-band licenses (for the entire X-band or S-band) for some, if not all, of their ground stations that are most capable of or likely to provide SECS. It is also recommended that certain prior arrangements be made with the relevant local licensing authorities for a process that will allow expedited authorization to transmit/receive signals to/from the declared spacecraft over the declared ground stations specifically and solely for the emergency case. Our analysis of the SECS has also uncovered a few fundamental programmatic issues. Recognizing any IOAG positions reached on these issues do not necessarily lead to any binding authority, it is recommended that they, along with those key attributes of the cross support emergency system, be explicitly stated as multi-agency guidelines to guide the implementation and provision of the SECS. The paper will present the results of the working group on Spacecraft Emergency Cross Support and, in particular, its findings, products and recommendations. It will also identify the next steps of this work that may include the production of some international standards as an extension of this work, or the process to make this emergency system usable by other spacecraft operators beyond the space agencies.

Tai, Wallace S.↗