Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 217 records · Page 12

Preliminary Estimation of Black Carbon Deposition from Nepal Climate Observatory-Pyramid Data and Its Possible Impact on Snow Albedo Changes Over Himalayan Glaciers During the Pre-Monsoon Season

The possible minimal range of reduction in snow surface albedo due to dry deposition of black carbon (BC) in the pre-monsoon period (March-May) was estimated as a lower bound together with the estimation of its accuracy, based on atmospheric observations at the Nepal Climate Observatory-Pyramid (NCO-P) sited at 5079 m a.s.l. in the Himalayan region. We estimated a total BC deposition rate of 2.89 g m-2 day-1 providing a total deposition of 266 micrograms/ square m for March-May at the site, based on a calculation with a minimal deposition velocity of 1.0 10(exp -4) m/s with atmospheric data of equivalent BC concentration. Main BC size at NCO-P site was determined as 103.1-669.8 nm by correlation analysis between equivalent BC concentration and particulate size distribution in the atmosphere. We also estimated BC deposition from the size distribution data and found that 8.7% of the estimated dry deposition corresponds to the estimated BC deposition from equivalent BC concentration data. If all the BC is deposited uniformly on the top 2-cm pure snow, the corresponding BC concentration is 26.0-68.2 microgram/kg assuming snow density variations of 195-512 kg/ cubic m of Yala Glacier close to NCO-P site. Such a concentration of BC in snow could result in 2.0-5.2% albedo reductions. From a simple numerical calculations and if assuming these albedo reductions continue throughout the year, this would lead to a runoff increases of 70-204 mm of water drainage equivalent of 11.6-33.9% of the annual discharge of a typical Tibetan glacier. Our estimates of BC concentration in snow surface for pre-monsoon season can be considered comparable to those at similar altitude in the Himalayan region, where glaciers and perpetual snow region starts in the vicinity of NCO-P. Our estimates from only BC are likely to represent a lower bound for snow albedo reductions, since a fixed slower deposition velocity was used and atmospheric wind and turbulence effects, snow aging, dust deposition, and snow albedo feedbacks were not considered. This study represents the first investigation about BC deposition on snow from atmospheric aerosol data in Himalayas and related albedo effect is especially the first track at the southern slope of Himalayas.

Yasunari, T. J.↗

Aqua's First 10 Years: An Overview

NASA's Aqua spacecraft was launched at 2:55 a.m. on May 4, 2002, from Vandenberg Air Force Base in California, into a near-polar, sun-synchronous orbit at an altitude of 705 km. Aqua carries six Earth-observing instruments to collect data on water in all its forms (liquid, vapor, and solid) and on a wide variety of additional Earth system variables (Parkinson 2003). The design lifetime for Aqua's prime mission was 6 years, and Aqua is now well into its extended mission, approaching 10 years of successful operations. The Aqua data have been used for hundreds of scientific studies and continue to be used for scientific discovery and numerous practical applications.

Parkinson, Claire L.↗

Earth Science Data Analysis in the Era of Big Data

Anyone with even a cursory interest in information technology cannot help but recognize that "Big Data" is one of the most fashionable catchphrases of late. From accurate voice and facial recognition, language translation, and airfare prediction and comparison, to monitoring the real-time spread of flu, Big Data techniques have been applied to many seemingly intractable problems with spectacular successes. They appear to be a rewarding way to approach many currently unsolved problems. Few fields of research can claim a longer history with problems involving voluminous data than Earth science. The problems we are facing today with our Earth's future are more complex and carry potentially graver consequences than the examples given above. How has our climate changed? Beside natural variations, what is causing these changes? What are the processes involved and through what mechanisms are these connected? How will they impact life as we know it? In attempts to answer these questions, we have resorted to observations and numerical simulations with ever-finer resolutions, which continue to feed the "data deluge." Plausibly, many Earth scientists are wondering: How will Big Data technologies benefit Earth science research? As an example from the global water cycle, one subdomain among many in Earth science, how would these technologies accelerate the analysis of decades of global precipitation to ascertain the changes in its characteristics, to validate these changes in predictive climate models, and to infer the implications of these changes to ecosystems, economies, and public health? Earth science researchers need a viable way to harness the power of Big Data technologies to analyze large volumes and varieties of data with velocity and veracity. Beyond providing speedy data analysis capabilities, Big Data technologies can also play a crucial, albeit indirect, role in boosting scientific productivity by facilitating effective collaboration within an analysis environment. To illustrate the effects of combining a Big Data technology with an effective means of collaboration, we relate the (fictitious) experience of an early-career Earth science researcher a few years beyond the present, interlaced and contrasted with reminiscences of its recent past (i.e., the present).

Kuo, K.-S.↗

Is the Sky Really Falling? An Overview of Orbital Debris

Orbital debris has been a prominent topic for a while, even before the movie Gravity came out. An anti-satellite test and a collision with an operational satellite both produced large highly-publicized debris clouds within recent years. While large objects like abandoned satellites and rocket bodies may be the most recognizable and identifiable concerns, a majority of the daily threat comes from the much more numerous smaller particles. In fact, small particle penetration continues to rank among the leading risks for manned space missions to the International Space Station and beyond. How much 'stuff' is up there, where did it come from, what harm can it do, and what is being done about it? These questions and more will be discussed.

Orbital Debris↗

Mini-Stamp as a Micro-Display for At-A-Glance Subsystem Information for DSN Links

Operators of the Deep Space Network (DSN) attend to numerous tasks with the overall goal of providing continuous support for the world's deep space missions. This high-stakes operations environment requires operators to understand the state of the Deep Space Network and predict what will happen next. Under the Follow-the-Sun initiative which requires remote operations of the highly complex telecommunications equipment, operators will need to remain aware of the state of the entire network rather than just their own facility, and transitioning fluidly between periods of low activity and periods of high demand. I designed a micro-display for operators to see, at a glance, the state of a Deep Space Network support including its subsystems. Using in-depth user-centered and participatory design techniques to identify information requirements, I designed what I called a Postage Stamp (NTR-49720) for individual operators to be able to maintain awareness of their own assigned supports. However, under Follow the Sun, operators must remain aware of all supports. The area occupied by the Postage Stamp must shrink to allow operators to see the state of the entire system, e.g., via a Big Board posted prominently in the operations room. Micro-displays are tools for mental model re-alignment, helping operators to keep their mental models of how the system works and behaves aligned with the changing state of the complex system. Data-driven micro-displays such as the Postage Stamp and Mini-Stamp display information about the system in a consistent way. Like a traffic light, the format of the micro-display never changes: the operator always knows where to look to find a specific piece of information. The Mini-Stamp always looks like the Mini-Stamp, and all of its data fields always lie in the same place on the micro-display. Real-time data flows through the Mini-Stamp to provide information to the operator.

Holloway, Alexandra↗

International VLBI Service for Geodesy and Astrometry 2017+2018 Biennial Report

Since its creation in 1999, the International VLBI Service for Geodesy and Astrometry (IVS) has documented its progress and current status in the form of annual or biennial reports. The first sixteen years were recorded in the form of annual reports, while starting with the years 2015+2016 the rhythm was changed to two years. Hence, the 2017+2018 Biennial Report constitutes the second installment of the two-yearly reporting scheme, documenting the work of the IVS components for the calendar years 2017 and 2018. As a testament to its usefulness, the general structure of the reporting has remained stable over the years: the individual components of the IVS contributed short reports describing their numerous activities, progress, and future plans. Without the continued input from the VLBI groups of the international geodetic and astrometric community, this publication could not be compiled and the IVS itself would not be able to flourish. So, once again many thanks to all IVS components who contributed to this Biennial Report.

Armstrong, Kyla L.↗

International VLBI Service for Geodesy and Astrometry 2021+2022 Biennial Report

Since its creation in 1999, the International VLBI Service for Geodesy and Astrometry (IVS) has documented its progress and current status in the form of annual or biennial reports. The first sixteen years were recorded in the form of annual reports. Starting in 2015/2016 the cadence was changed to two years; hence, this Biennial Report marks eight years of using the two-yearly reporting scheme. The general structure of the reporting has remained stable over the years: the individual components of the IVS contributed short reports describing their numerous activities, progress, and future plans. Without the continued input from the VLBI groups of the international geodetic and astrometric community, this publication could not be compiled and the IVS itself would not be able to flourish. So, once again many thanks to all IVS components who contributed to this Biennial Report

Space Geodesy↗

Asymptotic continuation method for trajectory optimization

A continuation method is applied to a singular perturbation parameter to obtain a new numerical method for computing optimal trajectories. This method allows one to use simply calculated reduced order approximations as starting solutions and continue the perturbation parameter until the optimal full-order solution is obtained. The method does not require the calculation of higher order correction terms nor does it require the perturbation parameter to be small - thus, it has potentially superior convergence properties compared to conventional asymptotic expansions when the perturbation parameter is large. A simple trajectory optimization problem is considered to illustrate the method.

Washburn, R. B., Jr.↗

Shifting the closed-loop spectrum in the optimal linear quadratic regulator problem for hereditary systems

In the optimal linear quadratic regulator problem for finite dimensional systems, the method known as an alpha-shift can be used to produce a closed-loop system whose spectrum lies to the left of some specified vertical line; that is, a closed-loop system with a prescribed degree of stability. This paper treats the extension of the alpha-shift to hereditary systems. As infinite dimensions, the shift can be accomplished by adding alpha times the identity to the open-loop semigroup generator and then solving an optimal regulator problem. However, this approach does not work with a new approximation scheme for hereditary control problems recently developed by Kappel and Salamon. Since this scheme is among the best to date for the numerical solution of the linear regulator problem for hereditary systems, an alternative method for shifting the closed-loop spectrum is needed. An alpha-shift technique that can be used with the Kappel-Salamon approximation scheme is developed. Both the continuous-time and discrete-time problems are considered. A numerical example which demonstrates the feasibility of the method is included.

Gibson, J. S.↗

Shifting the closed-loop spectrum in the optimal linear quadratic regulator problem for hereditary systems

In the optimal linear quadratic regulator problem for finite dimensional systems, the method known as an alpha-shift can be used to produce a closed-loop system whose spectrum lies to the left of some specified vertical line; that is, a closed-loop system with a prescribed degree of stability. This paper treats the extension of the alpha-shift to hereditary systems. As infinite dimensions, the shift can be accomplished by adding alpha times the identity to the open-loop semigroup generator and then solving an optimal regulator problem. However, this approach does not work with a new approximation scheme for hereditary control problems recently developed by Kappel and Salamon. Since this scheme is among the best to date for the numerical solution of the linear regulator problem for hereditary systems, an alternative method for shifting the closed-loop spectrum is needed. An alpha-shift technique that can be used with the Kappel-Salamon approximation scheme is developed. Both the continuous-time and discrete-time problems are considered. A numerical example which demonstrates the feasibility of the method is included.

Gibson, J. S.↗

Numerical algorithms for computations of feedback laws arising in control of flexible systems

Several continuous models will be examined, which describe flexible structures with boundary or point control/observation. Issues related to the computation of feedback laws are examined (particularly stabilizing feedbacks) with sensors and actuators located either on the boundary or at specific point locations of the structure. One of the main difficulties is due to the great sensitivity of the system (hyperbolic systems with unbounded control actions), with respect to perturbations caused either by uncertainty of the model or by the errors introduced in implementing numerical algorithms. Thus, special care must be taken in the choice of the appropriate numerical schemes which eventually lead to implementable finite dimensional solutions. Finite dimensional algorithms are constructed on a basis of a priority analysis of the properties of the original, continuous (infinite diversional) systems with the following criteria in mind: (1) convergence and stability of the algorithms and (2) robustness (reasonable insensitivity with respect to the unknown parameters of the systems). Examples with mixed finite element methods and spectral methods are provided.

Lasiecka, Irena↗

Computation of Transonic Flows Using Potential Methods

The proposed paper will describe the state of the art associated with numerical solution of the full or exact velocity potential equation for solving transonic, external-aerodynamic flows. The presentation will begin with a review of the literature emphasizing research activities of the past decade. Next, the various forms of the full or exact velocity potential equation, the equation's corresponding mathematical characteristics, and the derivation assumptions will be presented and described in detail. Impact of the derivation assumptions on simulation accuracy, especially with respect to shock wave capture, will be presented and discussed relative to the more complete Euler or Navier-Stokes formulations. The technical presentation will continue with a description of recently developed full potential numerical approach characteristics. This description will include governing equation nondimensionalization, physical-to-computational-domain mapping procedures, a limited description of grid generation requirements, the spatial discretization scheme, numerical implementation of boundary conditions, and the iteration scheme. The next portion of the presentation will present and discuss numerical results for several two- and three-dimensional aerodynamic applications. Included in the results section will be a discussion and demonstration of a typical grid refinement analysis for determining spatial convergence of the numerical solution and level of solution accuracy. Computer timings for a variety of full potential applications will be compared and contrasted with similar results for the Euler equation formulation. Finally. the presentation will end with concluding remarks and recommendations for future work.

Hoist, Terry L.↗

Plasma expansion characteristics of ionized clouds in the ionosphere - Macroscopic formulation

A macroscopic plasma expansion model, based on a numerical solution of the type-dependent nonlinear coupled continuity and momentum equations for background O(+) ions and several released ion species, was used to examine plasma expansion characteristics of ionized clouds in the ionosphere. Information is obtained on Ba(+), Li(+), and Ba(+) - Li(+) clouds; different cloud sizes (Gaussian half widths of 0.1, 0.5, and 1 km); cloud/background ion density ratios covering two orders of magnitude; electron/ion temperature ratios of 1, 5, and 10; and several cloud-background relative velocities along B (0, 1, 2, 4, 6, and 8 km/sec). The macroscopic expansion features obtained were found to be in general agreement with those obtained from the small-scale numerical simulations.

Schunk, R. W.↗

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗

Synthesis of offset dual reflector antennas transforming a given feed illumination pattern into a specified aperture distribution

The problem of transforming a given primary feed pattern into a desired aperture field distribution through two reflections by an offset dual reflector system is investigated using the concepts of geometrical optics. A numerically rigorous solution for the reflector surfaces is developed which realizes an exact aperture phase distribution and an aperture amplitude distribution that is accurate to within an arbitrarily small numerical tolerance. However, this procedure does not always yield a smooth solution, i.e., the reflector surfaces thus realized may not be continuous or their slopes may vary too rapidly. In the event of nonexistence of a numerically rigorous smooth solution, an approximate solution that enforces the smoothness of the reflector surfaces can be obtained. In the approximate solution, only the requirement for the aperture amplitude distribution is relaxed, and the condition on the aperture phase distribution is continued to be satisfied exactly.

Mittra, R.↗

Staffing the ISS Control Centers: Lessons Learned from Long-Duration Human Space Flight

The International Space Station (ISS) has been in operation with a permanent human presence in space for over five years, and plans for continued operations stretch ten years into the future. Ground control and support operations are, likewise, a 15-year enterprise. This long-term, 24-hour per day, 7 day per week support has presented numerous challenges in the areas of ground crew training, initial and continued certification, and console staffing. The Mission Control Center in Houston, Texas and the Payload Operations Center in Huntsville, Alabama have both tackled these challenges, with similar, yet distinct, approaches. This paper describes the evolution of the staffing and training policies of both control centers in a chronological progression. The relative merits and shortcomings of the various policies employed are discussed and a summary of "lessons learned" is presented. Finally, recommendations are made as best practices for future long-term space missions.

Olsen, Carrie D.↗

A combined algorithm for minimum time slewing of flexible spacecraft

The use of Pontryagin's Maximum Principle for the large-angle slewing of large flexible structures usually results in the so-called two-point boundary-value problem (TPBVP), in which many requirements (e.g., minimum time, small flexible amplitude, and limited control powers, etc.) must be satisfied simultaneously. The successful solution of this problem depends largely on the use of an efficient numerical computational algorithm. There are many candidate algorithms available for his problem (e.g., quasilinearization, gradient, and shooting, etc.). In this paper, a proposed algorithm, which combines the quasilinearization method with a time shortening technique and a shooting method, is applied to the minimum-time, three-dimensional, and large-angle maneuver of flexible spacecraft, particularly the orbiting Spacecraft Control Laboratory Experiment (SCOLE) configuration. Theoretically, the nonlinear TPBVP can be solved only through the shooting method to find the 'exact' switching times for the bang-bang controls. However, computationally, a suitable guess for the missing initial costates is crucial because the convergence range of the unknown initial costates is usually narrow, especially for systems with high dimensions and when a multi-bang-bang control strategy is needed. On the other hand, the problems of near minimum time attitude maneuver of general rigid spacecraft and fast slewing of flexible spacecraft have been examined by the authors through a numerical approach based on the quasilinearization algorithm with a time shortening technique. Computational results have demonstrated its broad convergence range and insensitivity to initial costate choices. Consequently, a combined approach is naturally suggested here to solve the minimum time slewing problem. That is, in the computational process, the quasilinearization method is used first to obtain a near minimum time solution. Then, the acquired converged initial costates from the quasilinearization approach are transformed (tailored) to and used as the initial costate guess for starting the shooting method. Finally, the shooting method takes over the remaining calculations until the minimum-time solution converges. The nonlinear equations of motion of the SCOLE are formulated by using Lagrange's equations, with the mast modeled as a continuous beam subject to three-dimensional deformations. The numerical results will be presented and some related computational issues will also be discussed.

Bainum, P. M.↗

Numerical simulation of backward step and jet exhaust flows

Preparatory to a continuing study of efficient, high resolution flow structure capture on patched mesh systems, turbulent flowfield solutions for a tactical missile base flow problem are computed on two first level composite grids of differing topologies. One mesh is of the wraparound character and the other is a step mesh on the Cartesian topology. A simplified Baldwin-Lomax algebraic eddy viscosity model is applied in the wake region. Results for the two meshes are compared with each other and with other computations reported in the literature. The totality of results suggest both mesh topology and resolution are important factors in the accuracy of predictive results as well as turbulence modeling. Related to interests in high altitude laminar flows bluff bodies such as the AOTV, the paper also presents results from a solution from laminar flow over a backward facing step and makes comparison with the available classical experiment of Chapman, et al.

Lombard, C. K.↗