SEARCH · Search NASA
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.
Engineering planning document surveys mechanics, utility, and potential improvements of codes which control noise-generated errors introduced into radio transmission of binary information. Chronological development of certain codes and a comprehensive bibliography are included.
This poster is about the project that I did this year as a CCI intern. The project is code to control a waveform generator as well as send custom waveform types.
A transient, semi-coupled, multi-dimensional thermal and flow analysis methodology was developed to predict the thermal/gas dynamic conditions in the field joint region of the Redesigned Solid Rocket Motor (RSRM). Transient temperature response, pressure history, and combustion gas flow rate (within the field joint region), were of principle interest, in the course of this study. The thermal environment in the field joint was modeled using SINDA, a finite difference based thermal network analyzer. The combustion gas flow boundary condition was generated using the FLAP code; this code performs a transient, lumped-parameter, control volume analysis to solve the mass, momentum, and energy conservation equations. The FLAP computer code was modified to account for erosion of the NBR insulation material, following ignition. An independent grid sensitivity study was conducted to determine an appropriate grid distribution near the wall. The predicted results, obtained using an optimum grid distribution and computer generated flow boundary condition, were compared with subscale test data.
An explicit-implicit and an implicit two-dimensional Navier-Stokes code along with various grid generation capabilities were developed. A series of classical benckmark cases were simulated using these codes.
This study supports the NASA Langley sponsored project aimed at determining the viability of using Euler technology for preliminary design use. The primary objective of this study was to assess the accuracy and efficiency of the Boeing, St. Louis unstructured grid flow field analysis system, consisting of the MACGS grid generation and NASTD flow solver codes. Euler solutions about the Aero Configuration/Weapons Fighter Technology (ACWFT) 1204 aircraft configuration were generated. Several variations of the geometry were investigated including a standard wing, cambered wing, deflected elevon, and deflected body flap. A wide range of flow conditions, most of which were in the non-linear regimes of the flight envelope, including variations in speed (subsonic, transonic, supersonic), angles of attack, and sideslip were investigated. Several flowfield non-linearities were present in these solutions including shock waves, vortical flows and the resulting interactions. The accuracy of this method was evaluated by comparing solutions with test data and Navier-Stokes solutions. The ability to accurately predict lateral-directional characteristics and control effectiveness was investigated by computing solutions with sideslip, and with deflected control surfaces. Problem set up times and computational resource requirements were documented and used to evaluate the efficiency of this approach for use in the fast paced preliminary design environment.
Program invariants can help software developers identify program properties that must be preserved as the software evolves, however, formulating correct invariants can be challenging. In this work, we introduce iDiscovery, a technique which leverages symbolic execution to improve the quality of dynamically discovered invariants computed by Daikon. Candidate invariants generated by Daikon are synthesized into assertions and instrumented onto the program. The instrumented code is executed symbolically to generate new test cases that are fed back to Daikon to help further re ne the set of candidate invariants. This feedback loop is executed until a x-point is reached. To mitigate the cost of symbolic execution, we present optimizations to prune the symbolic state space and to reduce the complexity of the generated path conditions. We also leverage recent advances in constraint solution reuse techniques to avoid computing results for the same constraints across iterations. Experimental results show that iDiscovery converges to a set of higher quality invariants compared to the initial set of candidate invariants in a small number of iterations.
3D heat flux calculations at the lower outer divertor plate of SPARC using the HEAT code show that 3D fields generated from error field correction coils can lead to enhanced peak heat fluxes up to 15 times larger compared to the axisymmetric case. Previously employed to simulate axisymmetric heat flux on 3D plasma facing components, the HEAT code can now predict 3D heat flux generated by non-axisymmetric plasmas. This is achieved via a new HEAT module which leverages the 3D field line tracing capabilities of MAFOT starting from an M3D-C1 (MHD resistive code) perturbed equilibria. The resulting heat flux is assigned using the magnetic footprint and the heat flux layer model, an extension of the 2D heat flux model also known as the Eich, to 3D non-axisymmetric plasmas. For SPARC, the new capabilities of HEAT are used to calculate the 3D heat loads resulting from n = 1 perturbation fields (with n indicating the toroidal periodicity) applied through a toroidal array of six picture frame coils with different amplitude. The comparison with the unperturbed case shows significant changes in shape and intensity of the heat flux profile. The results show that the application of n = 1 3D field leads to a localized enhancement of the heat flux peak, influenced by the wetted area impacted by the magnetic footprint, and the appearance of a secondary heat flux peak, whose intensity depends on amplitude of the applied 3D field and toroidal location.
A method is presented for calculating the locations and sizes of delaminations which occur in a rectangular, fiber reinforced composite plate subjected to nonpenetrating (low velocity) impact of a solid object. The plate may be simply supported or clamped along its edges. In-plane loads or in-plane strains may be imposed on the plate during the impact. The method includes two steps. First, the stresses and strains in the plate are calculated by a three-dimensional, transient finite element method using 8-node brick elements with incompatible modes. Second, the locations, lengths, and widths of delaminations inside the plate are predicted by means of a proposed failure criterion, which is based on the concept of dimensional analysis. The finite element method and the failure criterion were implemented by a computer code which can be used to calculate the impactor position and velocity, the displacements of the plate, the stresses and strains inside the plate during the impact, and the locations and dimensions of the delaminations after the impact. Parametric studies were performed to illustrate the information which can be generated by the computer code.
Enhanced geothermal systems (EGS) involve fracturing low permeability material to establish well connectivity and then injecting and circulating fluid into the fractured subsurface for geothermal power production. Changes in fracture aperture from contraction of the cooling matrix rock may alter network connectivity and risk thermal short-circuiting. Thermo-hydro-mechanical (THM) models are a useful tool to study these processes. However, as fracture networks are complex, and data may be limited, fracture networks in THM models are often stochastically generated. Given reliance on stochastic fracture networks and THM modeling to represent the subsurface and assess productivity of EGS, increased understanding of the influence of such statistically derived fracture networks on flow and heat transport in THM models is needed. Here, a new fracture process model is developed in the reactive transport code PFLOTRAN to stochastically generate fracture families and simulate changes in fracture aperture over time due to temperature changes of the rock matrix. Sixty-four different fracture networks ranging from well to poorly-connected, are modeled in PFLOTRAN with and without mechanical processes (THM vs TH). Results indicate that for well-connected fracture networks, thermal short-circuiting is less of a concern due to the abundance of available alternative flowpaths. For poorly-connected fracture networks, inclusion of mechanical processes showed steep thermal drawdown coincident with increase in fracture aperture along developing colder flowpaths, demonstrating the risk of thermal short-circuiting. Simulations with additional, larger fractures engineered to establish connectivity in a poorly-fractured subsurface, indicate that while stochastic variation of fracture orientation of the background network had limited influence, such variation in the engineered fractures significantly affected flow and heat transport.
Improvements in microscopy have enabled the acquisition of data at a scale that is difficult to process manually, making automated machine learning approaches to analyzing experimental images essential. In this project, we developed and applied machine learning (ML) workflows for atomic resolution scanning transmission electron microscopy (STEM) images. This development included improving both methodology as well as generating user-friendly codes. We developed machine learning architectures which, after training, automatically identify the location and types of defects throughout a material. We used these data to produce class-averaged images of 2D atomic coordinates with up to 0.3 pm precision, uncovering the structure and oscillations of long-range strain fields around point defects in WSe 2-2x Te 2x . We also resolved a long-standing problem in this field in the training of ML models, a lack of labeled experimental data, by developing a cycle-GAN that transformed simulated-generated labeled data into labeled data indistinguishable from experiment and therefore suitable for training. This removed the remaining parts of the ML data processing workflow where human intervention was still critical and therefore a bottleneck to working at scale. Codes have been developed and released for this full machine learning workflow. ML approaches to partially automate STEM acquisition were also developed. Finally we applied ML and other advanced data processing methods to several materials science problems in two-dimensional materials, including studying the evolution of hyperuniformity with defect concentration in WSe2, understanding phase transformations in transition metal dichalcogenides during in-situ heating in the STEM, and exploring how 2D interfaces transform from twisted into aligned structures.
A methodology for modeling nonlinear unsteady aerodynamic responses, for subsequent use in aeroservoelastic analysis and design, using the Volterra-Wiener theory of nonlinear systems is presented. The methodology is extended to predict nonlinear unsteady aerodynamic responses of arbitrary frequency. The Volterra-Wiener theory uses multidimensional convolution integrals to predict the response of nonlinear systems to arbitrary inputs. The CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) code is used to generate linear and nonlinear unit impulse responses that correspond to each of the integrals for a rectangular wing with a NACA 0012 section with pitch and plunge degrees of freedom. The computed kernels then are used to predict linear and nonlinear unsteady aerodynamic responses via convolution and compared to responses obtained using the CAP-TSD code directly. The results indicate that the approach can be used to predict linear unsteady aerodynamic responses exactly for any input amplitude or frequency at a significant cost savings. Convolution of the nonlinear terms results in nonlinear unsteady aerodynamic responses that compare reasonably well with those computed using the CAP-TSD code directly but at significant computational cost savings.
This paper presents an alternate approach for the generation of volumetric grids for supersonic and hypersonic flows about complex configurations. The method uses parametric two dimensional block face grid definition within the framework of GRIDGEN2D. The incorporation of face decomposition reduces complex surfaces to simple shapes. These simple shapes are combined to obtain the final face definition. The advantages of this method include the reduction of overall grid generation time through the use of vectorized computer code, the elimination of the need to generate matching block faces, and the implementation of simplified boundary conditions. A simple axisymmetric grid is used to illustrate this method. In addition, volume grids for two complex configurations, the Langley Lifting Body (HL-20) and the Space Shuttle Orbiter, are shown.
Plasma instabilities excited in collisionless shocks are responsible for particle acceleration. We have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electronpositron jet propagating into an unmagnetized electron-positron plasma. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic-like shock structure. In the leading shock, electron density increases by a factor of about 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. These magnetic fields contribute to the electrons transverse deflection behind the shock. We calculate the radiation from deflected electrons in the turbulent magnetic fields. Radiation from electrons near the trailing shock will be variable due to fluctuations of density and electromagnetic fields. The properties of this radiation may be important for rapid variability in relativistic jets such as AGN jets and blazars.
Langley Research Center has developed a novel external path sensor which monitors water vapor along an optical path between an airplane window and reflective material on the plane's engine. An infrared tunable diode laser is wavelength modulated across a water vapor absorption line at a frequency f. The 2f and DC signals are measured by a detector mounted adjacent to the laser. The 2f/DC ratio depends on the amount of wavelength modulation, the water vapor absorption line being observed, and the temperature, pressure, and water vapor content of the atmosphere. The present work concerns efforts to quantify the contributions of these factors and to derive a method for extracting the water vapor mixing ratio from the measurements. A 3 m cell was fabricated in order to perform laboratory tests of the sensor. Measurements of 2f/DC were made for a series of pressures and modulation amplitudes. During my 1994 faculty fellowship, a computer program was created which allowed 2f/DC to be calculated for any combination of the variables which effect it. This code was used to generate 2f/DC values for the conditions measured in the laboratory. The experimental and theoretical values agreed to within a few percent. As a result, the laser modulation amplitude can now be set in the field by comparing the response of the instrument to the calculated response as a function of modulation amplitude. Once the validity of the computer code was established, it was used to investigate possible candidate absorption lines. 2f/DC values were calculated for pressures, temperatures, and water vapor mixing ratios expected to be encountered in future missions. The results have been incorporated into a database which will be used to select the best line for a particular mission. The database will also be used to select a retrieval technique. For examples under some circumstances there is little temperature dependence in 2f/DC so temperature can be neglected. In other cases, there is a dependence with temperature for a particular pressure, requiring a more complicated retrieval algorithm. Future experimental work is necessary to test agreement with the theoretical values over a range of temperatures and mixing ratios. Additionally, retrieval algorithms for forthcoming missions must be incorporated into the software package which controls the instrument.
The technique of Pressure-Sensitive Paint (PSP) is commonly used in the aerospace industry to measure surface pressures on the model of launch vehicles and airplanes in the wind tunnel test. Recent research has demonstrated that Unsteady Pressure-Sensitive Paint (uPSP) can be an essential tool for the assessment of the unsteady, aerodynamic phenomena. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. This paper describes the spectral analysis of integrated pressures on patches of the scale model of the Space Launch System (SLS) Block 1 cargo vehicle with the uPSP measurements, which were collected in the Ascent Transient Aerodynamics Test (ATAT) with the Unitary Plan Wind Tunnel 11-by-11-foot Transonic Wind Tunnel in September 2019 at NASA Ames Research Center. The patches are defined with x station values, indicating the position along length of the SLS vehicle, and azimuth angles of the scale model. For each patch, the polygons are determined from the surface cells of the grid of the model, clipped with the edges of the patch, and each of the polygons is divided into triangles. The integrated pressure of the patch is determined as the ratio of the sum of the forces on the triangles over the sum of the areas of the triangles. For each run of the test, the Cross Power Spectral Density (CPSD) and magnitude squared coherence are computed from the time series of the integrated pressures on the patches. The pressure integration is coded in C++ and the spectral analysis is coded in MATLAB. The results were generated with the execution of the compiled C++ and MATLAB codes in parallel on the NASA Pleiades supercomputer. The results of pressure integration and spectral analysis are presented in this paper, and the data consistency of the test is also demonstrated. Funding for this research was provided by the NASA Aerosciences Evaluation and Test Capabilities Project.
Our MSAD-funded research project is to provide numerical modeling support for the VIBES (Vibration Isolation Box Experiment System) which is an IML2 flight experiment being built by the Japanese research team of Dr. H. Azuma of the Japanese National Aerospace Laboratory. During this reporting period, the following have been accomplished: A semi-consistent mass finite element projection algorithm for 2D and 3D Boussinesq flows has been implemented on Sun, HP And Cray Platforms. The algorithm has better phase speed accuracy than similar finite difference or lumped mass finite element algorithms, an attribute which is essential for addressing realistic g-jitter effects as well as convectively-dominated transient systems. The projection algorithm has been benchmarked against solutions generated via the commercial code FIDAP. The algorithm appears to be accurate as well as computationally efficient. Optimization and potential parallelization studies are underway. Our implementation to date has focused on execution of the basic algorithm with at most a concern for vectorization. The initial time-varying gravity Boussinesq flow simulation is being set up. The mesh is being designed and the input file is being generated. Some preliminary 'small mesh' cases will be attempted on our HP9000/735 while our request to MSAD for supercomputing resources is being addressed. The Japanese research team for VIBES was visited, the current set up and status of the physical experiment was obtained and ongoing E-Mail communication link was established.
The Lunar Atmosphere Dust Environment Explorer (LADEE) mission will orbit the moon in order to measure the density, composition and time variability of the lunar dust environment. The ground-side and onboard flight software for the mission is being developed using a Model-Based Software methodology. In this technique, models of the spacecraft and flight software are developed in a graphical dynamics modeling package. Flight Software requirements are prototyped and refined using the simulated models. After the model is shown to work as desired in this simulation framework, C-code software is automatically generated from the models. The generated software is then tested in real time Processor-in-the-Loop and Hardware-in-the-Loop test beds. Travelling Road Show test beds were used for early integration tests with payloads and other subsystems. Traditional techniques for verifying computational sciences models are used to characterize the spacecraft simulation. A lightweight set of formal methods analysis, static analysis, formal inspection and code coverage analyses are utilized to further reduce defects in the onboard flight software artifacts. These techniques are applied early and often in the development process, iteratively increasing the capabilities of the software and the fidelity of the vehicle models and test beds.
The Lunar Atmosphere Dust Environment Explorer (LADEE) mission orbited the moon in order to measure the density, composition, and time variability of the lunar dust environment. The successful mission launched September 7, 2013 and was de-orbited and impacted the moon's surface on April 17, 2014. The ground-side and onboard flight software for the mission was developed using a “Model-Based Software Engineering” (MBSE) methodology combined with strong reuse of Government and Commercial Off-The Shelf (G/COTS) components. Models of the spacecraft and flight software were developed in a graphical dynamics modeling package. Flight Software requirements were prototyped and refined using the simulated models. After the model was shown to work as desired in the simulation framework, C-code software was automatically generated from the models. The auto-generated software was then tested in real-time Processor-in-the-Loop and Hardware-in-the-Loop test beds. “Traveling Road Show” test beds were used for early integration tests with payloads and other subsystems. Traditional techniques for verifying computational sciences models were used to characterize the spacecraft simulation. A lightweight set of formal methods analysis, static analysis, formal inspection, and code coverage analyses were utilized to further reduce defects in the onboard flight software artifacts. These techniques were applied early and often in the development process, iteratively increasing the capabilities of software and fidelity of vehicle models and test beds.