Analysis of NO Formation During Earth Entry Using the Quasi Classical Trajectory Method
This is a summary of the research carried out by the author while he was a summer intern in the Aerothermodynamics Branch at NASA Ames.
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.
This is a summary of the research carried out by the author while he was a summer intern in the Aerothermodynamics Branch at NASA Ames.
Modern aircraft design methods have produced acceptable designs for large conventional aircraft performance. With revolutionary electronic propulsion technologies fueled by the growth in the small UAS (Unmanned Aerial Systems) industry, these same prediction models are being applied to new smaller, and experimental design concepts requiring a VTOL (Vertical Take Off and Landing) capability for ODM (On Demand Mobility). A 50% sub-scale GL-10 flight model was built and tested to demonstrate the transition from hover to forward flight utilizing DEP (Distributed Electric Propulsion)[1][2]. In 2016 plans were put in place to conduct performance flight testing on the 50% sub-scale GL-10 flight model to support a NASA project called DELIVER (Design Environment for Novel Vertical Lift Vehicles). DELIVER was investigating the feasibility of including smaller and more experimental aircraft configurations into a NASA design tool called NDARC (NASA Design and Analysis of Rotorcraft)[3]. This report covers the performance flight data collected during flight testing of the GL-10 50% sub-scale flight model conducted at Beaver Dam Airpark, VA. Overall the flight test data provides great insight into how well our existing conceptual design tools predict the performance of small scale experimental DEP concepts. Low fidelity conceptual design tools estimated the (L/D)( sub max)of the GL-10 50% sub-scale flight model to be 16. Experimentally measured (L/D)( sub max) for the GL-10 50% scale flight model was 7.2. The aerodynamic performance predicted versus measured highlights the complexity of wing and nacelle interactions which is not currently accounted for in existing low fidelity tools.
Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is applied to the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow which propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillarity to direct the movement of bubbles in space is the fact that surfactant impurities which are unavoidably present in the continuous phase can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature tension gradient, and decreases to near zero the thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have been retarded by the adsorption onto the bubble surface of a surfactant impurity, Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition the rapid kinetic exchange and high bulk concentration maintain a saturated surface with a uniform surface concentrations. This prevents retarding surface tension gradients and keeps the velocity high. In our first report last year, we detailed experimental results which verified the theory of remobilization in ground based experiments in which the steady velocity of rising bubbles was measured in a continuous phase consisting of a glycerol/water mixture containing a polyethylene glycol surfactant C12E6 (CH3(CH2)11(OCH2CH2)6OH). In our report this year, we detail our efforts to describe theoretically the remobilization observed. We construct a model in which a bubble rises steadily by buoyancy in a continuous (Newtonian) viscous fluid containing surfactant with a uniform far field bulk concentration. We account for the effects of inertia as well as viscosity in the flow in the continuous phase caused by the bubble motion (order one Reynolds number), and we assume that the bubble shape remains spherical (viscous and inertial forces are smaller than capillary forces, i e. small Weber and capillary numbers). The surfactant distribution is calculated by solving the mass transfer equations including convection and diffusion in the bulk, and finite kinetic exchange the bulk and the surface. Convective effects dominate diffusive mass transfer in the bulk of the liquid (high Peclet numbers) except in a thin boundary layer near the surface. A finite volume method is used to numerically solve the hydrodynamic and mass transfer equations on a staggered grid which accounts specifically for the thin boundary layer. We present the results of the nondimensional drag as a function of the bulk concentration of surfactant for different rates of kinetic exchange, from which we develop criteria for the concentration necessary to develop a prescribed degree of remobilization. The criteria compare favorably with the experimental results.
Continuing growth in regional and global air travel has resulted in increasing traffic congestion in the air and on the ground. In spite of occasional temporary downturns due to economic recessions and catastrophic events, average growth rates of air travel have remained high since the 1960s. The resulting congestion, which constrains expansion of the air transportation industry, inflicts schedule delays and decreases overall system efficiency, creating a pressing need to develop more efficient methods of air traffic management (ATM). New ATM techniques, procedures, air space automation methods, and decision support tools are being researched and developed for deployment in time frames stretching from the next few years to the year 2020 and beyond. As these methods become more advanced and increase in complexity, the requirements for information generation, sharing and transfer among the relevant entities in the ATM system increase dramatically. However, current aeronautical communications systems will be inadequate to meet the future information transfer demands created by these advanced ATM systems. Therefore, the NASA Glenn Research Center is undertaking research programs to develop communication, methods and key technologies that can meet these future requirements. As part of this process, studies, workshops, testing and experimentation, and research and analysis have established a number of research and technology development needs. The purpose of this paper is to outline the critical research and technology needs that have been identified in these activities, and explain how these needs have been determined.
Isotope harvesting is a method of collecting the long-lived radioisotopes that build up during the operation of ion-beam facilities in a way that is useful for subsequent research. As a demonstration of this method for the collection of a group IV metal at a fragmentation facility, the high-energy Zr 88 secondary beam produced from a 140-MeV/u Mo 92 primary beam at the National Superconducting Cyclotron Laboratory (NSCL) was stopped in a water target. The setup aimed to mimic the aqueous beam dump that will be implemented at the Facility for Rare Isotope Beams (FRIB). The collected Zr 88 and accompanying Y 88 decay daughter were radiochemically extracted from the solution and made into target samples suitable for neutron-capture cross-section measurements. These samples were then irradiated at two reactor facilities, and the Zr 88 average thermal-neutron-capture cross section ( σ T ) and resonance integral ( I ) were determined to be σ T = ( 8.04 ± 0.63 ) × 10 5 b and I = ( 2.53 ± 0.28 ) × 10 6 b. The σ T value agrees well with previous results and I , determined for the first time here, was found to be the largest measured resonance integral by two orders of magnitude. The Y 88 thermal-neutron-capture cross section was determined to be less than 1.8 × 10 4 b. This work demonstrates the steps needed to make cross-section measurements with samples produced via aqueous isotope harvesting.
Five methods for obtaining linear regression fits to bivariate data with unknown or insignificant measurement errors are discussed: ordinary least-squares (OLS) regression of Y on X, OLS regression of X on Y, the bisector of the two OLS lines, orthogonal regression, and 'reduced major-axis' regression. These methods have been used by various researchers in observational astronomy, most importantly in cosmic distance scale applications. Formulas for calculating the slope and intercept coefficients and their uncertainties are given for all the methods, including a new general form of the OLS variance estimates. The accuracy of the formulas was confirmed using numerical simulations. The applicability of the procedures is discussed with respect to their mathematical properties, the nature of the astronomical data under consideration, and the scientific purpose of the regression. It is found that, for problems needing symmetrical treatment of the variables, the OLS bisector performs significantly better than orthogonal or reduced major-axis regression.
Molecular biosignatures are key targets for current, proposed, and future life detection missions. With the high accuracy and low limit of detection (LOD) that new and future instruments will require, decontamination of life detection hardware is necessary to prevent false positives. Lipids are a molecular biosignature of interest, as they are ubiquitous to all life as we know it, can survive unaltered in the geologic record for longer than any other biomolecule (i.e. billions of years), and form through both biotic and abiotic processes. Lipids display origin-diagnostic molecular patterns that can reveal biotic or abiotic synthesis, so finding them and ascertaining their molecular features is important for potentially detecting evidence of life elsewhere. Traditional methods of decontamination, or contamination control (CC), primarily clean hardware through fabrication in sterile (cleanroom) environments, killing microbes, and removing/flushing contaminants off instrument and spacecraft components. However, research suggests that some standard cleaning methods are either unlikely to remove lipid contaminants or are incompatible with life detection instrument materials. To solve this problem, I propose to find, test, and verify a decontamination method that thoroughly cleans instruments by destroying lipid molecules, but is simultaneously compatible with major materials used in these instruments. I will study the effects of traditional CC methods (including Dry Heat Microbial Reduction and Vapor phase Hydrogen Peroxide) and experimental CC methods (Electron Beam Irradiation) on lipid molecules for application to life detection instrumentation. I will then develop a CC plan for a novel lipid detector (ExCALiBR, Extractor for Chemical Analysis of Lipid Biomarkers in Regolith) searching for lipids in either soil or icy world scenarios. This plan will uphold planetary protection regulation requirements and validate experimental analyses of in-situ life detection tests.
Thermocapillary migration is a method for moving bubbles in space in the absence of buoyancy. A temperature gradient is the continuous phase in which a bubble is situated, and the applied gradient impressed on the bubble surface causes one pole of the drop to be cooler than the opposite pole. As the surface tension is a decreasing function of temperature, the cooler pole pulls at the warmer pole, creating a flow that propels the bubble in the direction of the warmer fluid. A major impediment to the practical use of thermocapillary to direct the movement of bubbles in space is the fact that surfactant impurities, which are unavoidably present in the continuous phase, can significantly reduce the migration velocity. A surfactant impurity adsorbed onto the bubble interface is swept to the trailing end of the bubble. When bulk concentrations are low (which is the case with an impurity), diffusion of surfactant to the front end is slow relative to convection, and surfactant collects at the back end of the bubble. Collection at the back lowers the surface tension relative to the front end setting up a reverse tension gradient. (This can also be the case if kinetic desorption of surfactant at the back end of the bubble is much slower than convection.) For buoyancy driven bubble motions in the absence of a thermocapillarity, the tension gradient opposes the surface flow, and reduces the surface and terminal velocities (the interface becomes more solid-like and bubbles translate as solid particles). When thermocapillary forces are present, the reverse tension gradient set up by the surfactant accumulation reduces the temperature-induced tension gradient, and can decrease to near zero the bubble's thermocapillary velocity. The objective of our research is to develop a method for enhancing the thermocapillary migration of bubbles which have be retarded by the adsorption onto the bubble surface of a surfactant impurity. Our remobilization theory proposes to use surfactant molecules which kinetically rapidly exchange between the bulk and the surface and are at high bulk concentrations. Because the remobilizing surfactant is present at much higher concentrations than the impurity, it adsorbs to the bubble surface much faster than the impurity when the bubble is formed, and thereby prevents the impurity from adsorbing onto the surface. In addition, the rapid kinetic exchange and high bulk concentration maintain a saturated surface with uniform surface concentrations. This prevents retarding surface tension gradients and keeps the thermocapillary velocity high. In our reports over the first 2 years, we presented numerical simulations of the bubble motion and surfactant transport which verified theoretically the concept of remobilization, and the development of an apparatus to track and measure the velocity of rising bubbles in a glycerol/water surfactant solution. This year, we detail experimental observations of remobilization. Two polyethylene oxide surfactants were studied, C12E6 (CH3(CH2)11(OCH2)6OH) and C10E8 (CH3(CH2)4(OCH2CH2)8OH). Measurements of the kinetic exchange for these surfactants show that the one with the longer hydrophobe chain C12E6 has a lower rate of kinetic exchange. In addition, this surfactant is much less soluble in the glycerol/water mixture because of the shorter ethoxylate chain. As a result, we found that C12E6 had only a very limited ability to remobilize rising bubbles because of the limited kinetic exchange and reduced solubility. However, C10E8, with its higher solubility and more rapid exchange was found to dramatically remobilize rising bubbles. We also compared our theoretical calculations to the experimental measurements of velocity for both the non-remobilizing and remobilizing surfactants and found excellent agreement. We further observed that for C10E8 at high concentrations, which exceeded the critical micelle concentrations, additional remobilization was measured. In this case the rapid exchange of monomer between micelle and surfactant provides an additional mechanism for maintaining a uniform surface concentrations.
This document summarizes the findings of a strategic planning exercise commissioned by the Weapons Simulation and Computing, Computational Physics (WSC/CP) program at the Lawrence Livermore National Laboratory (LLNL) in FY24. During the year, the committee met with multiple stakeholder communities to gather input, opinions, suggestions and concerns which have been incorporated throughout this document. The key findings from this exercise are summarized: • The development of multiphysics modelling and simulation (mod/sim) codes and software technologies, their deployment on exascale compute platforms, and their broad adoption across the NNSA is a major success of the Advanced Simulation and Computing (ASC) program and the Exascale Computing Project (ECP). Sustained investment in these core technologies is essential. • Today’s state of the art involves running ensembles of O(100K) simulations to perform uncertainty quantification (UQ) and design studies using multiple statistical methods such as Bayesian optimization to understand sensitivities of our models and explore parameterized design spaces. Even with exascale computing, we are practically limited to O(10) parameters in these studies since the number of simulations required to sample the space scales exponentially with the number of design parameters. • The data from these simulation ensembles is increasingly being used to train machine learned (ML) surrogates (or reduced order models, ROMs) which can then be used for optimization or real time design exploration. However, the trained surrogates are still limited in the number of parameters they can represent due to the sampling limitations previously noted. • Augmenting our suite of integrated multiphysics simulation codes, both current and emerging, with the ability to compute gradients (solution derivatives) of arbitrary simulation outputs with respect to (some or all) simulation inputs would be a breakthrough technology, opening the door to a new era of efficient and automated inverse design based on verified and validated mod/sim capabilities. • This capability, which we refer to as differentiable multiphysics codes (DMCs), would revolutionize both UQ and optimization studies by breaking the curse of dimensionality that presently limits our “gradient-free” ensemble based computing approach. A similar breakthrough occurred in the AI/ML community once the ability to compute gradients of arbitrary loss functions using back-propagation became commonplace. Gradient information from the multiphysics codes can also be used to dramatically improve the efficiency and scale of training of ML/ROM surrogates for rapid assessments. • Achieving this in our suite of codes will be a grand challenge, similar to the amount of effort that was required to transition from CPU to GPU computing. It will require buy-in from the entire WSC/CP program and beyond, including all integrated codes, physics and engineering models, third-party library dependencies and performance portability abstractions. It will also require investment in research and development of numerical methods for computing adjoints of coupled physics across multiple adaptively refined moving meshes and of stochastic (Monte Carlo) and mesh free (SPH) methods. • New software and numerical techniques, largely pioneered by the AI/ML community, make this feasible. Chief among these is automatic differentiation (AD), the ability to employ AD at point-wise locations in a physics calculation (instead of traditional black-box approaches) and the ability to perform “back-propagation in time” (or reverse mode AD) for non-linear partial differential equations (PDEs). Fundamentally, the conclusion of this strategic planning exercise is that the time is right to undertake a large scale effort in WSC, centered on the existing integrated codes, to continue the natural evolution of mod/sim in the age of AI/ML. Instead of attempting to replace mod/sim with purely data driven AI/ML models, we believe the key to success is to integrate AI/ML by building on top of the decades of hard-won knowledge and the verified/validated multiphysics modelling capability that is the hallmark of the ASC program.
This paper deals with selected linearized-aerodynamic and vortex-flow methods as applied to aircraft design problems at high subsonic speeds. In particular, the NASA Vortex Lattice and Modified Multhopp methods are the linearized techniques employed, and the suction analogy is used to provide estimates associated with vortex-flow aerodynamics. Many examples are given as to how researchers at Langley have used these methods to design the high subsonic, wing-mean-camber shapes for various configurations such as a supersonic transport, high-aspect-ratio transport, trapezoidal fighter wing, strake wing, tandem wing, joined wing, delta wing, and slender cranked wing. Many of these have been built, tested, and have had their data compared with theory. In addition, a technique for defining efficiently performing strake planforms for use in strake-wing combinations is discussed, and further improvements in wing design are outlined. The latter may be obtained by using higher-ordered linear panel methods as well as nonlinear-transonic methods.
Research activity is described for: ice protection systems, icing instrumentation, experimental methods, analytical modeling for the above, and in flight research. The renewed interest in aircraft icing has come about because of the new need for All-Weather Helicopters and General Aviation aircraft. Because of increased fuel costs, tomorrow's Commercial Transport aircraft will also require new types of ice protection systems and better estimates of the aeropenalties caused by ice on unprotected surfaces. The physics of aircraft icing is very similar to the icing that occurs on ground structures and structures at sea; all involve droplets that freeze on the surfaces because of the cold air. Therefore all icing research groups will benefit greatly by sharing their research information.
Development of probabilistic structural analysis methods for hot engine structures is a major activity at Lewis Research Center. Recent activities have focused on extending the methods to include the combined uncertainties in several factors on structural response. This paper briefly describes recent progress on composite load spectra models, probabilistic finite element structural analysis, and probabilistic strength degradation modeling. Progress is described in terms of fundamental concepts, computer code development, and representative numerical results.
Durability and damage tolerance (D&DT) issues are critical to the development of lighter, safer and more efficient aerospace vehicles. Durability is largely an economic life-cycle design consideration whereas damage tolerance directly addresses the structural airworthiness (safety) of the vehicle. Both D&DT methodologies must address the deleterious effects of changes in material properties and the initiation and growth of damage that may occur during the vehicle s service lifetime. The result of unanticipated D&DT response is often manifested in the form of catastrophic and potentially fatal accidents. As such, durability and damage tolerance requirements must be rigorously addressed for commercial transport aircraft and NASA spacecraft systems. This paper presents an overview of the recent and planned future research in durability and damage tolerance analytical and experimental methods for both metallic and composite aerospace structures at NASA Langley Research Center (LaRC).
A compact and rugged probe based on the phase Doppler method was evaluated as a means for characterizing icing clouds using airborne platforms and for advancing aircraft icing research in large scale wind tunnels. The Phase Doppler Particle Analyzer (PDPA) upon which the new probe was based is now widely recognized as an accurate method for the complete characterization of sprays. The prototype fiber optic-based probe was evaluated in simulated aircraft icing clouds and found to have the qualities essential to providing information that will advance aircraft icing research. Measurement comparisons of the size and velocity distributions made with the standard PDPA and the fiber optic probe were in excellent agreement as were the measurements of number density and liquid water content. Preliminary testing in the NASA Lewis Icing Research Tunnel (IRT) produced reasonable results but revealed some problems with vibration and signal quality at high speeds. The cause of these problems were identified and design changes were proposed to eliminate the shortcomings of the probe.
A shake test and an extensive analysis of results were performed to evaluate the possibility of and the method for dynamically calibrating the Rotor Systems Research Aircraft (RSRA). The RSRA airframe was subjected to known vibratory loads in several degrees of freedom and the responses of many aircraft transducers were recorded. Analysis of the transducer responses using the technique of dynamic force determination showed that the RSRA, when used as a dynamic measurement system, could predict, a posteriori, an excitation force in a single axis to an accuracy of about 5 percent and sometimes better. As the analysis was broadened to include multiple degrees of freedom for the excitation force, the predictive ability of the measurement system degraded to about 20 percent, with the error occasionally reaching 100 percent. The poor performance of the measurement system is explained by the nonlinear response of the RSRA to vibratory forces and the inadequacy of the particular method used in accounting for this nonlinearity.
Detailed temporal measurements of the solar flux at one location are performed. These data are then analyzed and compared to the potential of space measurements which allow one to consider the flux falling on areas of the earth. An important result of the research is that the temporal characteristics of the flux in the presence of a real atmosphere would be difficult to obtain from space and that the variations in the flux can be highly significant in regard to most solar conversion schemes. The detailed results of the research are presented. The instruments developed to separate the direct and scattered solar flux, the computer analysis methods developed, and the results of the research, presented as both graphical and tabular data, are discussed.
Multicopters are becoming one of the more common and popular type of unmanned aircraft systems (UAS) which have both civilian and military applications. One example being the concept of drone deliveries proposed by the distribution company Amazon [1]. The electrical propulsion is considered to have both faster and easier deliveries and also environmental benefits compared to other vehicles that still use fossil fuel. Other examples include surveillance and just simple entertainment. The reason behind their success is often said to be due to their small size, relatively low cost, simple structure and finally simple usage. With an increase in the UAS market comes challenges in terms of security, as both people and other aircrafts could be harmed if not used correctly. Therefore further studies and regulations are needed to ensure that future use of drones, especially in the civilian and public sectors, are safe and efficient. Thorough research has been done on full scale, man or cargo transporting, helicopters so that most parts of flight and performance are fairly well understood. Yet not much of it have been verified for small multicopters. Until today many studies and research projects have been done on the control systems, navigation and aerodynamics of multicopters. Many of the methods used today for building multicopters involve a process of trial an error of what will work well together, and once that is accomplished some structural analysis of the multicopter bodies might be done to verify that the product will be strong enough and have a decent aerodynamic performance. However, not much has been done on the research of the rotor blades, especially in terms of structural stress analyses and ways to ensure that the commonly used parts are indeed safe and follow safety measures. Some producers claim that their propellers indeed have been tested, but again that usually tends towards simple fluid dynamic analyses and even simpler stress analyses. There is no real deflection measurement of said blades and all theory is today based on the theory developed for full scale helicopters. This report thus intends to highlight the problems that come with blade deflection theory and measurements for multicopters. This thesis starts with the introduction and problem formulation where the grounds for what the report contains are laid out. Then a literature review of the research within the area follows where previous and current research and methods are presented. In the background information and theory relevant to the area are presented. After that the methods and set up of the different experimental and computer simulation approaches are depicted. In the results section the obtained results are presented and compared. Lastly there is a discussion where the results and methods are evaluated and analyzed, followed by the most important conclusions.
The purpose of this research was to produce a measurable quantity of radioactive argon-37 through photonuclear production methods. The detection of argon-37, above background levels, uniquely and clearly indicates the underground detonation of a nuclear device.