SEARCH · Engineering Papers
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.
Current declarative systems are limited to a small number of types of states and resources that they can represent. Thus software is usually hand-crafted to meet the needs of the state validation-estimation-projection for real domains where estimation is not sufficient.
The instability and subsequent atomization of a viscous liquid jet emanated into a high-pressure gaseous surrounding is studied both computationally and experimentally. Liquid water issued into nitrogen gas at elevated pressures is used to simulate the flow conditions in a coaxial shear injector element relevant to liquid propellant rocket engines. The theoretical analysis is based on a simplified mathematical formulation of the continuity and momentum equations in their conservative form. Numerical solutions of the governing equations subject to appropriate initial and boundary conditions are obtained via a robust finite difference scheme. The computations yield real-time evolution and subsequent breakup characteristics of the liquid jet. The experimental investigation utilizes a digital imaging technique to measure resultant drop sizes. Data were collected for liquid Reynolds number between 2,500 and 25,000, aerodynamic Weber number range of 50-500 and ambient gas pressures from 150 to 1200 psia. Comparison of the model predictions and experimental data for drop sizes at gas pressures of 150 and 300 psia reveal satisfactory agreement particularly for lower values of investigated Weber number. The present model is intended as a component of a practical tool to facilitate design and optimization of coaxial shear atomizers.
The NASA Kepler mission recently announced over 1200 exoplanet candidates. While some are common Hot Jupiters, a large number are Neptune size and smaller, transit depths suggest sizes down to the radius of Earth. The Kepler project has a fairly high confidence that most of these candidates are real exoplanets. Many analysis steps and lessons learned from Kepler light curves are used during the vetting process. This talk will cover some new results in the areas of stellar variability, solar systems with multiple planets, and how transit-like signatures are vetted for false positives, especially those indicative of small planets.
Dynamic environments in which objectives and environmental features change with respect to time pose a difficult problem with regards to planning optimal paths through these environments. Path planning methods are typically computationally expensive, and are often difficult to implement in real time if system objectives are changed. This computational problem is compounded when multiple agents are present in the system, as the state and action space grows exponentially with the number of agents in the system. In this work, we use cooperative coevolutionary algorithms in order to develop policies which control agent motion in a dynamic multiagent unmanned aerial system environment such that goals and perceptions change, while ensuring safety constraints are not violated. Rather than replanning new paths when the environment changes, we develop a policy which can map the new environmental features to a trajectory for the agent while ensuring safe and reliable operation, while providing 92% of the theoretically optimal performance.
During inactive phases of Madden-Julian oscillation (MJO), there are plenty of deep but small convective systems and far fewer deep and large ones. During active phases of MJO, a manifestation of an increase in the occurrence of large and deep cloud clusters results from an amplification of large-scale motions by stronger convective heating. This study is designed to quantitatively examine the roles of small and large cloud clusters during the MJO life cycle. We analyze the cloud object data from Aqua CERES observations for tropical deep convective (DC) and cirrostratus (CS) cloud object types according to the real-time multivariate MJO index. The cloud object is a contiguous region of the earth with a single dominant cloud-system type. The size distributions, defined as the footprint numbers as a function of cloud object diameters, for particular MJO phases depart greatly from the combined (8-phase) distribution at large cloud-object diameters due to the reduced/increased numbers of cloud objects related to changes in the large-scale environments. The medium diameter corresponding to the combined distribution is determined and used to partition all cloud objects into "small" and "large" groups of a particular phase. The two groups corresponding to the combined distribution have nearly equal numbers of footprints. The medium diameters are 502 km for DC and 310 km for cirrostratus. The range of the variation between two extreme phases (typically, the most active and depressed phases) for the small group is 6-11% in terms of the numbers of cloud objects and the total footprint numbers. The corresponding range for the large group is 19-44%. In terms of the probability density functions of radiative and cloud physical properties, there are virtually no differences between the MJO phases for the small group, but there are significant differences for the large groups for both DC and CS types. These results suggest that the intreseasonal variation signals reside at the large cloud clusters while the small cloud clusters represent the background noises resulting from various types of the tropical waves with different wavenumbers and propagation directions/speeds.
NASA's Earth Observing System Data and Information System (EOSDIS) has been in operation since August 1994, and serving a diverse user community around the world with Earth science data from satellites, aircraft, field campaigns and research investigations. The ESDIS Project, responsible for EOSDIS is a Network Member of the International Council for Sciences (ICSU) World Data System (WDS). Nine of the 12 Distributed Active Archive Centers (DAACs), which are part of EOSDIS, are Regular Members of the ICSUWDS. This poster presents the EOSDIS mission objectives, key characteristics of the DAACs that make them world class Earth science data centers, successes, challenges and best practices of EOSDIS focusing on the years 2014-2016, and illustrates some highlights of accomplishments of EOSDIS. The highlights include: high customer satisfaction, growing archive and distribution volumes, exponential growth in number of products distributed to users around the world, unified metadata model and common metadata repository, flexibility provided to uses by supporting data transformations to suit their applications, near-real-time capabilities to support various operational and research applications, and full resolution image browse capabilities to help users select data of interest. The poster also illustrates how the ESDIS Project is actively involved in several US and international data system organizations.
The OpenSE Cookbook is an open-sourced collection of patterns, procedures, and best practices targeted for systems engineers who seek guidance on applying model-based and executable systems engineering (MBSE) using SysML. Its content has emerged from the system level modeling effort on the European Framework Program 6 (FP6) and the Thirty Meter Telescope (TMT). The TMT MBSE approach applied the Executable Systems Engineering Method (ESEM) and the open-source Engineering Environment (OpenMBEE) to specify, analyze, and verify requirements of TMT’s Alignment and Phasing System (APS) and the Narrow Field Infrared Adaptive Optics System (NFIRAOS). In these applications, implicit dependencies are made explicit in a formal model through the use of ESEM, OpenMBEE, and SysML modeling constructs. The value proposition for applying this MBSE approach was to establish precise requirements and fine-grained traceability to system designs, and to verify key requirements beginning early in development. The integration of ESEM and the OpenMBEE tooling infrastructure (providing linked-data and web-operability) is a significant added value for the MBSE approach. The APS is responsible for the overall pre-adaptive optics wavefront quality, using starlight to measure wavefront errors and align the TMT optics. In the formally integrated and executable SysML model, simulations are performed to analyze the impact of changed requirements and verify specified constraints for various operational scenarios. The APS team used several modeling patterns to capture information such as the requirements, the operational scenarios, involved subsystems and their interaction points, the estimated or required time durations, and the mass and power consumption. Adaptive optics systems are designed to sense real-time atmospheric turbulence and correct the telescope’s optical beam to remove its effect. The system model for the adaptive optics operational modes was developed to capture sequence behaviors and operational scenarios to run Monte-Carlo simulations for verifying acquisition time, observing efficiency, and operational behavior requirements. The model is particularly useful for investigating the effect of parallelization, identifying interface issues, and re-ordering sequence acquisition tasks. A former version of the Cookbook (which is now updated to MBSE challenges, goals, and lessons learned) included modeling guidelines and conventions for all system aspects, hierarchy levels, and views, which were developed during for the Active Phasing Experiment (APE), an opto-mechatronical system technology demonstrator for the Extremely Large Telescope (ELT). The Cookbook utilizes the above mentioned system models as real-world case-studies to demonstrate and document the applications of the recipes, providing also instructional examples and addressing the available tooling support. The Cookbook is accompanied by a number of SysML models and aodel libraries which facilitate model authoring and maintenance. The Cookbook covers the different aspects of Systems Engineering such as management of Requirements, Design (behavior and structure), Interfaces, Interdisciplinary Integration, Analysis, Trade Studies, and Technical Resources. This paper presents the background, motivation, architecture, and highlights some key content of the Cookbook. For example, interface management, error budget management, requirements verification, Monte Carlo driven analysis, and timing analysis of operational scenarios. The paper discusses how the capabilities of OpenMBEE contributed significantly to the adoption of executable systems engineering.
Contemporary advancements in low-cost automation and computation, reduced barrier to entry in developing artificial intelligence/machine learning (AI/ML), and increased ability to represent complex materials in digital form have led to a number of accelerated materials discovery platforms. However, many of these approaches operate with completely rigid vertical integration in an isolated feedback loop using limited modalities. In order to make a substantial impact on discovering new energy materials, AI-driven experiments must operate collaboratively with each other and researchers and over multiple measurement modalities. Herein, we describe the potential for an “internet of things” approach to self-driving enterprise beamlines that merges core information technologies, robotics, and multi-modal AI. The approach will enable full utility of light sources, collaborate effectively with other remote materials acceleration platforms, and help stride toward the world’s energy future.
Conversion-type electrode materials for lithium-ion batteries experience significant structural changes during the first discharge–charge cycle, where a single particle is taken apart into a number of nanoparticles. This structural evolution may affect the following lithium insertion reactions; however, how lithiation occurs in pre-cycled electrode materials is elusive. In this work, in situ transmission electron microscopy was employed to see the lithium-induced structural and chemical evolutions in pre-cycled nickel oxide as a model system. The introduction of lithium ions induced the evolution of metallic nickel, with volume expansion as a result of a conversion reaction. After pre-cycling, the phase evolutions occurred in two separate areas almost at the same time. This is different from the first lithiation, where the phase change takes place successively, with a boundary dividing the reacted and unreacted areas. Structural changes were restricted at the areas having large amount of fluorine, implying the residuals from the decomposition of electrolytes may have hindered the electrochemical reactions. This work provides insights into phase and chemical evolutions in pre-cycled conversion-type materials, which govern electrochemical properties during operation.
In this innovation, a digital downconverter has been created that produces a large (16 or greater) number of output channels of smaller bandwidths. Additionally, this design has the flexibility to tune each channel independently to anywhere in the input bandwidth to cover a wide range of output bandwidths (from 32 MHz down to 1 kHz). Both the flexibility in channel frequency selection and the more than four orders of magnitude range in output bandwidths (decimation rates from 32 to 640,000) presented significant challenges to be solved. The solution involved breaking the digital downconversion process into a two-stage process. The first stage is a 2 oversampled filter bank that divides the whole input bandwidth as a real input signal into seven overlapping, contiguous channels represented with complex samples. Using the symmetry of the sine and cosine functions in a similar way to that of an FFT (fast Fourier transform), this downconversion is very efficient and gives seven channels fixed in frequency. An arbitrary number of smaller bandwidth channels can be formed from second-stage downconverters placed after the first stage of downconversion. Because of the overlapping of the first stage, there is no gap in coverage of the entire input bandwidth. The input to any of the second-stage downconverting channels has a multiplexer that chooses one of the seven wideband channels from the first stage. These second-stage downconverters take up fewer resources because they operate at lower bandwidths than doing the entire downconversion process from the input bandwidth for each independent channel. These second-stage downconverters are each independent with fine frequency control tuning, providing extreme flexibility in positioning the center frequency of a downconverted channel. Finally, these second-stage downconverters have flexible decimation factors over four orders of magnitude The algorithm was developed to run in an FPGA (field programmable gate array) at input data sampling rates of up to 1,280 MHz. The current implementation takes a 1,280-MHz real input, and first breaks it up into seven 160-MHz complex channels, each spaced 80 MHz apart. The eighth channel at baseband was not required for this implementation, and led to more optimization. Afterwards, 16 second stage narrow band channels with independently tunable center frequencies and bandwidth settings are implemented A future implementation in a larger Xilinx FPGA will hold up to 32 independent second-stage channels.
Agile development methods, which have gone from outlier to mainstream in software development, are poised to expand into all aspects of space mission development. Modern software development operates on a principle of continuous deployment, where progress is verified not with conventional metrics, but with a continuous build, available to key stakeholders, enabling direct examination of the state of the code base, and assessment of progress through demonstration of capability. Delivery times are measured in weeks, not months. Stakeholders are part of the process on an ongoing basis. The cost of change is comparatively low and requirements, which often are not precisely defined at the start of a project, may be iteratively refined in a series of agile development cycles. Agile methods are compatible with traditional system engineering methods and may be tailored to the space operations environment. The low cost of change and iterative development cycles of agile enable requirements to be defined as outcomes and constraints, with design details to be refined during the development cycle. We are now at a point where agile methods may be extended beyond software, to Mission Systems, including the Mission Operations System and the Ground Data System. For NASA’s VIPER Lunar Rover Mission, scheduled to land at a lunar pole in late 2023, we are developing the Mission System using agile methods. As in agile software, where the measure of progress is working code, in agile mission system development, the measure of capability is what we can demonstrate. Demonstrations over presentations. We demonstrate mission system capability using simulations. The concept of operations, from commanding, to driving the rover, to how we downlink images for evaluation for a near-real time command cycle, will be tested and proven in simulation, years before we begin the traditional simulation cycle for training. “Say it then simulate it.” We develop and refine our designs using simulations, with an emphasis on new components of the system that are not well known early. For example, the required duration of a mission planning cycle for a lunar surface asset such as VIPER, that operates twenty-four hours a day, seven days a week, with continuous communications and a unique set of constraints based on the physics of the lunar poles and the line of site to Earth, is a unique problem in mission planning that is unlikely to be solved in a series of meetings. A small number of requirements specifying the outcomes may serve as the jumping off point to an agile development cycle, with demonstration in simulations. We have already demonstrated this process with simulations of rover driver decision time. VIPER is driven using near-real time command and control to waypoints. The driver decision time between waypoints is a fundamental enabling unit of productivity to accomplish the mission timeline. We have validated driver decision time in simulations of rover driving at the lunar South Pole, using the prototype mission tools for driving, command and control. The capability to develop and refine designs using simulations as part of agile Mission System development cycle changes the nature of team interactions, creating a focus on doing, rather than analyzing and documenting. Waterfall development cycles were, in part, a product of the significant cost of change in the early days of spaceflight. When the cost of change is high, it is vital to get your requirements right at the outset, because the system will be built to those specifications, and, when change is expensive, you better get it right early. However, modern technology has greatly lowered the cost of change, enabling iterative, rapid development cycles, in which key operations concepts may be tested and refined during development. Extending agile development to the Mission System for VIPER is a significant step in moving agile development methods for space operations beyond software, to the Mission System.
Implicit methods developed by Beam and Warming (1978 and Briley and McDonald (1977) make it possible to overcome the hyperbolic stiffness of the conservative compressible Navier-Stokes equations in the fine wall region computational mesh for high Reynold's number flow. Certain difficulties related to the use of these methods could be overcome by employing an approach reported by Roe (1981). In the present investigation Roe's conceptual framework has been adopted for constructing globally conservative finite difference methods. A globally conservative upwind finite difference method (CSCM) consisting of both implicit interior point and boundary point equations is constructed from a new characteristics based flux difference splitting. It is found that the employed upwind eigenvector split scheme which combines fully coupled implicit interior point and boundary point approximations has the desired properties of robust stability and accuracy.
The current definition of Coordinated Universal Time (UTC) dates from 1972. The duration of a UTC second is defined in terms of the frequency of a hyperfine transition in the ground state of cesium. This standard frequency is realized in a number of different laboratories using ensembles of commercial cesium clocks and a few primary frequency standards. The data from all of these devices are transmitted periodically to the Bureau International des Poids et Mesures (BIPM) in Sevres, France, where they are combined in a statistical procedure to produce International Atomic Time (TAI). The time of this scale is adjusted as needed ('coordinated') by adding or dropping integer seconds so as to keep it within plus or minus 0.9 s of UT1, a time scale based on the observation of the transit times of stars and corrected for the predicted seasonal variations in these observations. When the leap seconds are included into TAI, the result is called UTC. The difference between TAI and UTC is therefore an exact integer number of seconds. This difference is currently 29 s and will become 30 s at 0 UTC on 1 January 1996.
Local heat transfer coefficients were measured along the midchord of a three-times-size turbine vane airfoil in a static cascade operated at roon temperature over a range of Reynolds numbers. The test surface consisted of a composite of commercially available materials: a Mylar sheet with a layer of cholestric liquid crystals, which change color with temperature, and a heater made of a polyester sheet coated with vapor-deposited gold, which produces uniform heat flux. After the initial selection and calibration of the composite sheet, accurate, quantitative, and continuous heat transfer coefficients were mapped over the airfoil surface. Tests were conducted at two free-stream turbulence intensities: 0.6 percent, which is typical of wind tunnels; and 10 percent, which is typical of real engine conditions. In addition to a smooth airfoil, the effects of local leading-edge sand roughness were also examined for a value greater than the critical roughness. The local heat transfer coefficients are presented for both free-stream turbulence intensities for inlet Reynolds numbers from 1.20 to 5.55 x 10 to the 5th power. Comparisons are also made with analytical values of heat transfer coefficients obtained from the STAN5 boundary layer code.
Local heat transfer coefficients were measured along the midchord of a three-times-size turbine vane airfoil in a static cascade operated at room temperature over a range of Reynolds numbers. The test surface consisted of a composite of commercially available materials: a Mylar sheet with a layer of cholestric liquid crystals, which change color with temperature, and a heater made of a polyester sheet coated with vapor-deposited gold, which produces uniform heat flux. After the initial selection and calibration of the composite sheet, accurate, quantitative, and continuous heat transfer coefficients were mapped over the airfoil surface. Tests were conducted at two free-stream turbulence intensities: 0.6 percent, which is typical of wind tunnels; and 10 percent, which is typical of real engine conditions. In addition to a smooth airfoil, the effects of local leading-edge sand roughness were also examined for a value greater than the critical roughness. The local heat transfer coefficients are presented for both free-stream turbulence intensities for inlet Reynolds numbers from 1.20 to 5.55 x 10 to the 5th power. Comparisons are also made with analytical values of heat transfer coefficients obtained from the STAN5 boundary layer code.
We present a new methodology for computing sensitivities in evolutionary systems using a model-driven low-rank approximation. To this end, we formulate a variational principle that seeks to minimize the distance between the time derivative of the reduced approximation and sensitivity dynamics. The first order optimality condition of the variational principle leads to a system of closed form evolution equations for an orthonormal basis and corresponding sensitivity coefficients. This approach allows for the computation of sensitivities with respect to a large number of parameters in an accurate and tractable manner by extracting correlations between different sensitivities on the fly. The presented method requires solving forward evolution equations, sidestepping the restrictions imposed by the forward/backward workflow of adjoint sensitivities. For example, the presented method, unlike the adjoint equation, does not impose any input/output load and can be used in applications in which real-time sensitivities are of interest. We demonstrate the utility of the method for three test cases: (1) computing sensitivity with respect to model parameters in the Rössler system, (2) computing sensitivity with respect to an infinite-dimensional forcing parameter in the chaotic Kuramoto--Sivashinsky equation, and (3) computing sensitivity with respect to reaction parameters for species transport in a turbulent reacting flow.
Two-phase water-steam flow conditions are frequently encountered in many engineering applications, including geothermal reservoirs. Although routine calculations are based on the multiphase Darcy’s law, the role of the topology of the flowing phases at the pore-scale is usually neglected in the estimation of relative permeabilities. Instead, the latter are frequently computed using empirical models like the Corey correlation. In this work, we first apply the model for relative permeabilities based on pore-scale flow regimes developed by Picchi and Battiato (2019), Relative permeability scaling from pore-scale flow regimes, Water Resour. Res. 55, 3215–3233, to scenarios typical of geothermal reservoirs and then extend it by deriving the scaling laws for the transmissibilities and the thermal properties as a function of temperature. First, we discuss the scaling behavior of normalized relative permeabilities in terms of viscosity ratio and capillary number of water-steam systems and, then, we provide a validation of the model against experimental data available in the literature. The model captures well the data trends collected in real 3D porous media. These results suggest that water-steam relative permeabilities follow the same scaling behavior of gas-liquid systems where the non-wetting phase is much less viscous than the wetting phase. Finally, we investigate the impact that relative permeabilities have on heat transfer rates at two-phase flow conditions and the scaling of mass and energy transmissibility and thermal properties of the mixture. An estimation of the exergy carried by the two-phase water-steam mixture is also included.
Transition metal complexes capable of photoinduced spin crossover have been greatly investigated because of their potential to enable ultrafast optical control of information processing. Yet, any real application of photoswitchable molecules requires that spin crossover be paired with additional functionality such as long-range magnetic order. Important advances combining these functions are notably reported for a number of bimetallic Prussian Blue analogues; however, to date, PBA-based magnetic photoswitches can only operate below 150 K due to loss of magnetic order. In contrast, cobalt ferrite is a ferrimagnetic semiconductor with a Curie temperature of 790 K and extremely favorable magnetic properties by comparison to state-of-the-art PBAs. The mixed valence electronic structure of cobalt ferrite is reminiscent of cobalt–iron PBA, which is a well-known photoswitch. To investigate the potential for photoswitching in this material, we employ transient XUV spectroscopy to probe charge and spin dynamics with element-specific resolution on the femtosecond time scale. Results show that 400 nm light excites a metal-to-metal charge transfer transition, which drives the crossover of high-spin Co 2+ to low-spin Co 3+ with a time constant of 405 ± 29 fs and an internal quantum efficiency of unity. This result establishes the existence of efficient photoswitching in a new class of robust ferrimagnetic spinel ferrites.