Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational physics applications”

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 667 records · Page 37

Fundamentals of Plasma Sails Propulsion Concept

The Mini-Magnetospheric Plasma Propulsion (M2P2), originally proposed by Winglee et al. [2000], is based on the two-fluid plasma model and requires a 15-km frontal standoff distance (or 20-km cross-sectional diameter) in order for the magnetic bubble to absorb sufficient momentum from the SW to accelerate a spacecraft to the unprecedented speeds of 50-80 km/s after an acceleration period of about three months. Winglee et al. [2000] derived the above size requirement based on an extrapolation of their simulated results in which a system much smaller than a M2P2 was used (p. 21,074 of their study). We submit, however, that a fluid model has no validity for such a small scale size-even in the region near the plasma source! It is assumed in the MHD fluid model, normally applied to the magnetosphere, that the characteristic scale-size is much greater than the Larmor radius and ion skin depth of the SW. In the case of the M2P2, however, the size of the magnetic bubble is actually less than or, at best, comparable to, the scale of these characteristic parameters and, therefore, a kinetic approach, which addresses the smallscale physical mechanisms involved, must be used. A fully three-dimensional version of the hybrid code is used in our M2P2 (Plasma Sails) studies was originally developed by Delamere et al. [1999]. The M2P2 plasma sail is an excellent application for this hybrid code. The primary advantage of this code is the seamless interface between fluid and kinetic descriptions of the ion populations. A kinetic description is not necessary for the dense inner regions of the magnetic bubble and tremendous computational savings can be realized by treating this dense, magnetized ion population with the fluid description. It is essential, however, that the outer bubble regions be treated kinetically as well as the SW protons. Comparison of full size M2P2 simulation based on 3D MHD and kinetic models show that kinetic treatment introduces much more asymmetry to the considering problem and the possibility of kinetic instabilities development.

Khazanov, G. V.↗

Constrained Kalman Filtering Via Density Function Truncation for Turbofan Engine Health Estimation

Kalman filters are often used to estimate the state variables of a dynamic system. However, in the application of Kalman filters some known signal information is often either ignored or dealt with heuristically. For instance, state variable constraints (which may be based on physical considerations) are often neglected because they do not fit easily into the structure of the Kalman filter. This paper develops an analytic method of incorporating state variable inequality constraints in the Kalman filter. The resultant filter truncates the PDF (probability density function) of the Kalman filter estimate at the known constraints and then computes the constrained filter estimate as the mean of the truncated PDF. The incorporation of state variable constraints increases the computational effort of the filter but significantly improves its estimation accuracy. The improvement is demonstrated via simulation results obtained from a turbofan engine model. The turbofan engine model contains 3 state variables, 11 measurements, and 10 component health parameters. It is also shown that the truncated Kalman filter may be a more accurate way of incorporating inequality constraints than other constrained filters (e.g., the projection approach to constrained filtering).

Simon, Dan↗

Safe Maneuvering Envelope Estimation Based on a Physical Approach

This paper discusses a computationally efficient algorithm for estimating the safe maneuvering envelope of damaged aircraft. The algorithm performs a robust reachability analysis through an optimal control formulation while making use of time scale separation and taking into account uncertainties in the aerodynamic derivatives. This approach differs from others since it is physically inspired. This more transparent approach allows interpreting data in each step, and it is assumed that these physical models based upon flight dynamics theory will therefore facilitate certification for future real life applications.

aerodynamics↗

2040 Vision Study: an Enlargement of Model Based Engineering

Over the last few decades, advances in high-performance computing, new materials characterization methods, and, more recently, an emphasis on integrated computational materials engineering (ICME) and additive manufacturing have been a catalyst for multiscale modeling and simulation-based design of materials and structures in the aerospace industry. As a result, NASA's Transformational Tools and Technology (TTT) Project sponsored a study (performed by a team led by Pratt & Whitney) to define the potential 25-year future state required for integrated multiscale modeling of materials and systems (e.g., load-bearing structures) to accelerate the pace and reduce the expense of innovation in future aerospace and aeronautical systems. This talk will briefly review the findings of this 2040 Vision study (e.g., the 2040 vision state; the required interdependent core technical work areas, Key Element (KE); associated critical gaps and actions to close those gaps; and major recommendations). The study, NASA CR 2018- 219771, envisions the development of a cyber-physical-social ecosystem comprised of experimentally verified and validated computational models, tools, and techniques, along with the associated digital tapestry, that marries two non-mutually exclusive paradigms _ "design of the materials" (material scientist viewpoint) and "design with the materials" (structural analyst viewpoint) _ into a concurrent transformational paradigm that impacts the entire supply chain to enable cost-effective, rapid, and revolutionary design of fit-for-purpose materials, components, and systems. Although the vision focused on aeronautics and space applications, it is believed that other engineering communities (e.g., automotive, biomedical, etc.) can benefit as well from the proposed framework with only minor modifications. Finally, it is TTT's hope and desire that this vision provides the strategic guidance to both public and private research and development decision makers to make the proposed 2040 vision state a reality and thereby provide a significant advancement in the United States global competitiveness.

Arnold, Steven M.↗

Impacts of Aerosol and Cloud Variability on Satellite-Derived All-skies Direct Aerosol Radiative Effects (DARE) over the North and Southeast Atlantic

The Top-of-Atmosphere (TOA) Shortwave (SW) Direct Aerosol Radiative Effect (DARE) in all-sky conditions (i.e., aerosols in clear skies and aerosols above and below all types of clouds) is the global change in upwelling radiative flux due to aerosols. It is one of the strongest indicators of global climate change due to aerosols. SW DARE at TOA depends on the Earth’s surface albedo, cloud fraction, cloud optical properties, and aerosol optical properties, which are all challenging to accurately characterize from space. The overarching goals of our project are to provide state-of-the-art observational all-sky TOA SW DARE, along with guidance on which aerosol and/or cloud properties are the most important to measure, and at which spatio-temporal scales, for accurate DARE observations. We compute all-sky DARE based on state-of-the-art cloud and aerosol retrieval algorithms from CALIOP (Cloud–Aerosol Lidar with Orthogonal Polarization) and MODIS (Moderate Resolution Imaging spectroradiometer) satellite sensors, as well as aerosol intensive properties from MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications, version 2) simulations over three specific regions of the Atlantic Ocean from 2012 to 2016. In this symposium, we present a characterization of the cloud and aerosol optical and physical properties and the observational TOA SW DARE along three individual satellite tracks in the Southeast Atlantic region -- on September 18 and 20, 2016 and August 13, 2017. We then quantify the impact of assuming homogeneous cloud or aerosol fields in space-based TOA SW DARE calculations. And, finally, the resulting space-based all-skies aerosol vertical distribution and DARE calculations are evaluated using remote sensing observations from coincident suborbital flights from the NASA ORACLES (the ObseRvations of Aerosols above CLouds and their intEractionS, ORACLES field campaign) field campaign. The NASA Atmospheric Observing System (AOS) mission addresses the NASA Aerosol, Cloud, Convection and Precipitation (ACCP) designated observables and proposes, as one of its aerosol objectives, to reduce uncertainties in estimates of global mean all-sky SW DARE at TOA. Well characterizing clouds, aerosol vertical distributions, aerosol types and associated all-skies DARE over the Atlantic Ocean will inform the AOS community on where, when, how, and how often the satellite retrievals should be performed to estimate DARE and reduce all-skies DARE uncertainties most accurately. These comprehensive characterizations will also identify the key regions and times when the AOS (or other ACCP-related) suborbital missions should be conducted to evaluate and improve the AOS space-based observations and retrievals.

aerosol↗

2040 Vision Study: NASA’s T3 Implementation Activities

Over the last few decades, advances in high-performance computing, new materials characterization methods, and, more recently, an emphasis on integrated computational materials engineering (ICME) and additive manufacturing have been a catalyst for multiscale modeling and simulation-based design of materials and structures in the aerospace industry. As a result, NASA’s Transformational Tools and Technology (TTT) Project sponsored a study (performed by a team led by Pratt & Whitney) to define the potential 25-year future state required for integrated multiscale modeling of materials and systems (e.g., load-bearing structures) to accelerate the pace and reduce the expense of innovation in future aerospace and aeronautical systems. This talk will briefly review the findings of this 2040 Vision study (e.g., the 2040 vision state; the required interdependent core technical work areas, Key Element (KE); associated critical gaps and actions to close those gaps; and major recommendations) and discuss NASA’s TTT implementation activities, with special emphasis on recent accomplishments. The study, NASA CR 2018- 219771, envisions the development of a cyber-physical-social ecosystem comprised of experimentally verified and validated computational models, tools, and techniques, along with the associated digital tapestry, that marries two non-mutually exclusive paradigms – “design of the materials” (material scientist viewpoint) and “design with the materials” (structural analyst viewpoint) – into a concurrent transformational paradigm that impacts the entire supply chain to enable cost-effective, rapid, and revolutionary design of fit-for-purpose materials, components, and systems. Although the vision focused on aeronautics and space applications, it is believed that other engineering communities (e.g., automotive, biomedical, etc.) can benefit as well from the proposed framework with only minor modifications. Finally, it is TTT’s hope and desire that this vision provides the strategic guidance to both public and private research and development decision makers to make the proposed 2040 vision state a reality and thereby provide a significant advancement in the United State's global competitiveness

Materials↗

Elliptic generation of composite three-dimensional grids about realistic aircraft

An elliptic method for generating composite grids about realistic aircraft is presented. A body-conforming grid is first generated about the entire aircraft by the solution of Poisson's differential equation. This grid has relatively coarse spacing, and it covers the entire physical domain. At boundary surfaces, cell size is controlled and cell skewness is nearly eliminated by inhomogeneous terms, which are found automatically by the program. Certain regions of the grid in which high gradients are expected, and which map into rectangular solids in the computational domain, are then designated for zonal refinement. Spacing in the zonal grids is reduced by adding points with a simple, algebraic scheme. Details of the grid generation method are presented along with results of the present application, a wing-body configuration based on the F-16 fighter aircraft.

Sorenson, R. L.↗

Discontinuous solutions to hyperbolic systems under operator splitting

Two-dimensional systems of linear hyperbolic equations are studied with regard to their behavior under a solution strategy that in alternate time-steps solves exactly the component one-dimensional operators. The initial data is a step function across an oblique discontinuity. The manner in which this discontinuity breaks up under repeated applications of the split operator is analyzed, and it is shown that the split solution will fail to match the true solution in any case where the two operators do not share all their eigenvectors. The special case of the fluid flow equations is analyzed in more detail, and it is shown that arbitrary initial data gives rise to pseudo acoustic waves and a non-physical stationary wave. The implications of these findings for the design of high-resolution computing schemes are discussed.

Roe, P. L.↗

Applications of artificial intelligence to space station: General purpose intelligent sensor interface

This final report describes the accomplishments of the General Purpose Intelligent Sensor Interface task of the Applications of Artificial Intelligence to Space Station grant for the period from October 1, 1987 through September 30, 1988. Portions of the First Biannual Report not revised will not be included but only referenced. The goal is to develop an intelligent sensor system that will simplify the design and development of expert systems using sensors of the physical phenomena as a source of data. This research will concentrate on the integration of image processing sensors and voice processing sensors with a computer designed for expert system development. The result of this research will be the design and documentation of a system in which the user will not need to be an expert in such areas as image processing algorithms, local area networks, image processor hardware selection or interfacing, television camera selection, voice recognition hardware selection, or analog signal processing. The user will be able to access data from video or voice sensors through standard LISP statements without any need to know about the sensor hardware or software.

Mckee, James W.↗

Combining Comparison Functions and Finite Element Approximations in CFD

In a variety of potential flow applications, the modal element method has been shown to significantly reduce the numerical grid, employ a more precise grid termination boundary condition, and give theoretical insight to the flow physics. The method employs eigenfunctions to replace the numerical grid over significant portions of the flow field. Generally, a numerical grid is employed around obstacles with complex geometry while eigenfunctions are applied to regions in the flow field where the boundary conditions can easily be satisfied. To handle a wider class of computational fluid dynamics (CFD) problems, the present paper extends the modal element to include function approximations which do not satisfy the governing differential equation. To accomplish this task, a double modal series approximation and weighted residual constraints are developed to force the comparison functions to satisfy the governing differential equation and to interface properly with the finite element solution. As an example, the method is applied to the problem of potential flow in a channel with two-dimensional cylindrical like obstacles. The calculated flow fields are in excellent agreement with exact analytical solutions.

Baumeister, Kenneth J.↗

Evolutionary Computing Methods for Spectral Retrieval

A methodology for processing spectral images to retrieve information on underlying physical, chemical, and/or biological phenomena is based on evolutionary and related computational methods implemented in software. In a typical case, the solution (the information that one seeks to retrieve) consists of parameters of a mathematical model that represents one or more of the phenomena of interest. The methodology was developed for the initial purpose of retrieving the desired information from spectral image data acquired by remote-sensing instruments aimed at planets (including the Earth). Examples of information desired in such applications include trace gas concentrations, temperature profiles, surface types, day/night fractions, cloud/aerosol fractions, seasons, and viewing angles. The methodology is also potentially useful for retrieving information on chemical and/or biological hazards in terrestrial settings. In this methodology, one utilizes an iterative process that minimizes a fitness function indicative of the degree of dissimilarity between observed and synthetic spectral and angular data. The evolutionary computing methods that lie at the heart of this process yield a population of solutions (sets of the desired parameters) within an accuracy represented by a fitness-function value specified by the user. The evolutionary computing methods (ECM) used in this methodology are Genetic Algorithms and Simulated Annealing, both of which are well-established optimization techniques and have also been described in previous NASA Tech Briefs articles. These are embedded in a conceptual framework, represented in the architecture of the implementing software, that enables automatic retrieval of spectral and angular data and analysis of the retrieved solutions for uniqueness.

Terrile, Richard↗

Elliptic generation of composite three-dimensional grids about realistic aircraft

An elliptic method for generating composite grids about realistic aircraft is presented. A body-conforming grid is first generated about the entire aircraft by the solution of Poisson's differential equation. This grid has relatively coarse spacing, and it covers the entire physical domain. At boundary surfaces, cell size is controlled and cell skewness is nearly eliminated by inhomogeneous terms, which are found automatically by the program. Certain regions of the grid in which high gradients are expected, and which map into rectangular solids in the computational domain, are then designated for zonal refinement. Spacing in the zonal grids is reduced by adding points with a simple, algebraic scheme. Details of the grid-generation method are presented along with results of the present application, a wing/body configuration based on the F-16 fighter aircraft.

Sorenson, Reese L.↗

An application of MSC/NASTRAN in the interdisciplinary analysis of large space-based structures

The Integrated Multidisciplinary Analysis Tool (IMAT), a computer software system developed at NASA Langley to analyze and simulate the dynamics of space-structure/control-system interactions, is described, and its application to the MAST problem (a 60-m truss with fundamental frequency less than 200 mHz and equipped with linear proof-mass actuators, to be deployed from the Space Shuttle as part of COFS-I flight experiment) is demonstrated. Particular attention is given to the IMAT procedures which facilitate the use of the MCS/NASTRAN code to recover physical results from time-domain state-space solutions obtained with an FEM control-design code. Diagrams, drawings, and graphs are provided.

Stockwell, Alan E.↗

ATTRACTOR: Autonomy Teaming and TRAjectories for Complex Trusted Operational Reliability

Autonomous systems (AS) are crucial to realizing the vision of new, complex transportation modes, such as advanced air mobility (AAM) and urban air mobility (UAM). A chief barrier to induction of AS into aviation is insufficient understanding of AS reliability in time-critical and safety-critical environments—an obstacle to certification. ATTRACTOR is aimed at building a basis for certification of classes of autonomous cyber-physical-human systems (CPHS) via establishing metrics and models of trustworthiness and trust in multi-agent team interactions, analyzable trajectories, explainability of computational algorithms (explainable artificial intelligence, or XAI), and persistent modeling and simulation, in the context of missions planning and operation. By “building a basis for certification,” we mean acquiring an understanding of when a system is trustworthy and developing computable means to estimate trustworthiness and trust in order to eventually inform functional requirements that contribute to certification. The outcomes are applicable not just to aviation but to all domains that rely on autonomous systems.

Autonomy↗

Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Model

In-flight icing is an important consideration that affects aircraft design, performance and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of ice accretion testing conducted in the NASA Icing Research Tunnel on a large swept wing section typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small MVD = 25 μm and a large MVD = 110 μm. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. The results for the small MVD cloud conditions were consistent with previous work and showed the changing ice morphology from rime ice at colder temperatures to highly 3D scallop ice in the range of -8.7 ˚C to -3.8 ˚C. The results for the large MVD cloud conditions exhibited some differences in the ice accretion morphology from the small MVD conditions. Scallop-like ice features were observed at total temperature of -3.8˚C. The large MVD ice accretion at lower temperatures tended to be smoother than the corresponding ice accretion for small MVD. The thickness of the main ice shape tended to be larger for the small MVD conditions compared to the corresponding large MVD ice shape while the latter had more ice farther downstream. The measured ice mass was approximately equal between the corresponding small and large MVD ice shapes up to a freezing fraction of 0.6. For higher freezing fraction, the large MVD ice shapes weighed more. Cloud MVD variation from 50 μm to 230 μm while holding the scaling parameters constant resulted in very minor differences in the ice shapes, mostly in the upper and lower surface chordwise extents. Ice volume was computed from 3D scan data and used to calculate a ratio of ice mass to volume. The resulting values were in the range of 240 to 455 Kg/m3. While this is consistent with analogous values previously reported in the literature, more data are needed to determine a specific trend in this ratio as a function of freezing fraction and the corresponding governing physical phenomena. This work has resulted in a significant experimental database of ice accretion for small and large MVD conditions applicable to large-scale swept wings.

Icing↗

Large-Drop Ice Accretion Test Results for a Large Scale Swept Wing Section

In-flight icing is an important consideration that affects aircraft design, performance and safety. Newer regulations combined with increasing demand to reduce fuel burn, emissions and noise are driving a need for improvements in icing simulation capability. To that end, this paper presents the results of ice accretion testing conducted in the NASA Icing Research Tunnel on a large swept wing typical of a modern commercial transport. The model was based upon a section of the Common Research Model wing at the 64% semispan station with a streamwise chord length of 136 in. The test conditions were developed with an icing scaling analysis to generate similar conditions for a small MVD = 25 μm and a large MVD = 110 μm. A series of tests were conducted over a range of total temperature from -23.8 ˚C to -1.4 ˚C with all other conditions held constant. The results for the small MVD cloud conditions were consistent with previous work and showed the changing ice morphology from rime ice at colder temperatures to highly 3D scallop ice in the range of -8.7 ˚C to -3.8 ˚C. The results for the large MVD cloud conditions exhibited some differences in the ice accretion morphology. Scallop-like ice features were observed at total temperature of -3.8˚C. The large MVD ice accretion at lower temperatures tended to be smoother than the corresponding ice accretion for small MVD. The thickness of the main ice shape tended to be larger for the small MVD conditions compared to the corresponding large MVD ice shape while the latter had more ice farther downstream. The measured ice mass was approximately equal between the corresponding small and large MVD ice shapes up to a freezing fraction of 0.6. For higher freezing fraction, the large MVD ice shapes weighed more. Cloud MVD variation from 50 μm to 230 μm while holding the scaling parameters constant resulted in very minor differences in the ice shapes, mostly in the upper and lower surface chordwise extents. Ice volume was computed from 3D scan data and used to calculate a ratio of ice mass to volume. The resulting values were in the range of 240 to 455 Kg/m3. While this is consistent with analogous values previously reported in the literature, more data are needed to determine a specific trend in this ratio as a function of freezing fraction and the corresponding governing physical phenomena. This work has resulted in a significant experimental database of ice accretion for small and large MVD conditions applicable to large-scale swept wings.

Icing↗

Formal Specification and Parametric Verification of the ICAROUS Distributed Merging Protocol for Autonomous Aircraft Systems

ICAROUS is a software architecture that provides highly assured core software modules for building safety-centric autonomous unmanned aircraft applications. One of its core components is the ICAROUS distributed merging (IDM) protocol, which allows for decentralized merging of autonomous aircrafts through a designated intersection. This report presents initial results on formal specification and parametric verification of the IDM protocol. We present the development of a formal, discrete-time specification of the ICAROUS distributed merging protocol in TLA+. The developed TLA+ specification includes an abstracted model of the physical aircraft dynamics, the consensus machinery for leader election and coordination, and the computation of merging schedules. In addition, we present details on a command line tool we developed verimerge, that utilizes the TLC model checker for doing bounded, parametric verification and allows for plotting of these results in 2D parameter spaces. The tool also provides functionality for visualization of concrete protocol behaviors, to aid debugging and understanding. We present preliminary, bounded time verification results for a finite number of aircraft. Limitations of the current techniques and possible future extensions of this work are also discussed.

ICAROUS↗

Focal region fields of distorted reflectors

The problem of the focal region fields scattered by an arbitrary surface reflector under uniform plane wave illumination is solved. The physical optics (PO) approximation is used to calculate the current induced on the reflector. The surface of the reflector is described by a number of triangular domain-wise 5th degree bivariate polynomials. A 2-dimensional Gaussian quadrature is employed to numerically evaluate the integral expressions of the scattered fields. No Freshnel or Fraunhofer zone approximations are made. The relation of the focal fields problem to surface compensation techniques and other applications are mentioned. Several examples of distorted parabolic reflectors are presented. The computer code developed is included, together with instructions on its usage.

Buris, N. E.↗