Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “explicit time integration”

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 271 records · Page 15

An Open Computing Infrastructure that Facilitates Integrated Product and Process Development from a Decision-Based Perspective

Computer applications for design have evolved rapidly over the past several decades, and significant payoffs are being achieved by organizations through reductions in design cycle times. These applications are overwhelmed by the requirements imposed during complex, open engineering systems design. Organizations are faced with a number of different methodologies, numerous legacy disciplinary tools, and a very large amount of data. Yet they are also faced with few interdisciplinary tools for design collaboration or methods for achieving the revolutionary product designs required to maintain a competitive advantage in the future. These organizations are looking for a software infrastructure that integrates current corporate design practices with newer simulation and solution techniques. Such an infrastructure must be robust to changes in both corporate needs and enabling technologies. In addition, this infrastructure must be user-friendly, modular and scalable. This need is the motivation for the research described in this dissertation. The research is focused on the development of an open computing infrastructure that facilitates product and process design. In addition, this research explicitly deals with human interactions during design through a model that focuses on the role of a designer as that of decision-maker. The research perspective here is taken from that of design as a discipline with a focus on Decision-Based Design, Theory of Languages, Information Science, and Integration Technology. Given this background, a Model of IPPD is developed and implemented along the lines of a traditional experimental procedure: with the steps of establishing context, formalizing a theory, building an apparatus, conducting an experiment, reviewing results, and providing recommendations. Based on this Model, Design Processes and Specification can be explored in a structured and implementable architecture. An architecture for exploring design called DREAMS (Developing Robust Engineering Analysis Models and Specifications) has been developed which supports the activities of both meta-design and actual design execution. This is accomplished through a systematic process which is comprised of the stages of Formulation, Translation, and Evaluation. During this process, elements from a Design Specification are integrated into Design Processes. In addition, a software infrastructure was developed and is called IMAGE (Intelligent Multidisciplinary Aircraft Generation Environment). This represents a virtual apparatus in the Design Experiment conducted in this research. IMAGE is an innovative architecture because it explicitly supports design-related activities. This is accomplished through a GUI driven and Agent-based implementation of DREAMS. A HSCT design has been adopted from the Framework for Interdisciplinary Design Optimization (FIDO) and is implemented in IMAGE. This problem shows how Design Processes and Specification interact in a design system. In addition, the problem utilizes two different solution models concurrently: optimal and satisfying. The satisfying model allows for more design flexibility and allows a designer to maintain design freedom. As a result of following this experimental procedure, this infrastructure is an open system that it is robust to changes in both corporate needs and computer technologies. The development of this infrastructure leads to a number of significant intellectual contributions: 1) A new approach to implementing IPPD with the aid of a computer; 2) A formal Design Experiment; 3) A combined Process and Specification architecture that is language-based; 4) An infrastructure for exploring design; 5) An integration strategy for implementing computer resources; and 6) A seamless modeling language. The need for these contributions is emphasized by the demand by industry and government agencies for the development of these technologies.

Hale, Mark A.↗

The OpenSE Cookbook: A Practical, Recipe Based Collection of Patterns, Procedures, and Best Practices for Executable Systems Engineering for the Thirty Meter Telescope

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.

Brower, Eric↗

A Study of the Response of Deep Tropical Clouds to Mesoscale Processes: Sensitivities to Microphysics, Radiation, and Surface Fluxes - Part 2

The Goddard Cumulus Ensemble (GCE) model is used to examine the sensitivities of surface fluxes, explicit radiation, and ice microphysical processes on multi-day simulations of deep tropical convection over the Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE). The simulations incorporate large-scale advective temperature and moisture forcing, as well as large-scale momentum, that are updated every time step on a periodic lateral boundary grid. This study shows that when surface fluxes are eliminated, the mean atmosphere is much cooler and drier, convection and CAPE are much weaker, precipitation is less, and cloud coverage in stratiform regions much greater. Surface fluxes using the TOGA COARE flux algorithm are weaker than with the aerodynamic formulation, but closer to the observed fluxes. In addition, similar trends noted above for the case without surface fluxes are produced for the TOGA flux case, albeit to a much lesser extent. The elimination of explicit shortwave and longwave radiation is found to have only minimal effects on the mean thermodynamics, convection, and precipitation. However explicit radiation does have a significant impact on cloud temperatures and structure above 200 mb and on the overall mean vertical circulation. The removal of ice processes produces major changes in the structure of the cloud. Much of the liquid water is transported aloft and into anvils above the melting layer (600 mb), leaving narrow, but intense bands of rainfall in convective regions. The elimination of melting processes leads to greater hydrometeor mass below the melting layer, and produces a much warmer and moister boundary layer, leading to a greater mean CAPE. Finally, the elimination of the graupel species has only a small impact on mean total precipitation, thermodynamics, and dynamics of the simulation, but does produce much greater snow mass just above the melting layer. Some of these results differ from previous CRM studies of tropical systems, which is likely due to the type of simulated system, total time integration, and model setup.

Johnson, Daniel↗

Mixed models and reduction method for dynamic analysis of anisotropic shells

A time-domain computational procedure is presented for predicting the dynamic response of laminated anisotropic shells. The two key elements of the procedure are: (1) use of mixed finite element models having independent interpolation (shape) functions for stress resultants and generalized displacements for the spatial discretization of the shell, with the stress resultants allowed to be discontinuous at interelement boundaries; and (2) use of a dynamic reduction method, with the global approximation vectors consisting of the static solution and an orthogonal set of Lanczos vectors. The dynamic reduction is accomplished by means of successive application of the finite element method and the classical Rayleigh-Ritz technique. The finite element method is first used to generate the global approximation vectors. Then the Rayleigh-Ritz technique is used to generate a reduced system of ordinary differential equations in the amplitudes of these modes. The temporal integration of the reduced differential equations is performed by using an explicit half-station central difference scheme (Leap-frog method). The effectiveness of the proposed procedure is demonstrated by means of a numerical example and its advantages over reduction methods used with the displacement formulation are discussed.

Noor, A. K.↗

Large Eddy Simulation (LES) of Particle-Laden Temporal Mixing Layers

High-fidelity models of plume-regolith interaction are difficult to develop because of the widely disparate flow conditions that exist in this process. The gas in the core of a rocket plume can often be modeled as a time-dependent, high-temperature, turbulent, reacting continuum flow. However, due to the vacuum conditions on the lunar surface, the mean molecular path in the outer parts of the plume is too long for the continuum assumption to remain valid. Molecular methods are better suited to model this region of the flow. Finally, granular and multiphase flow models must be employed to describe the dust and debris that are displaced from the surface, as well as how a crater is formed in the regolith. At present, standard commercial CFD (computational fluid dynamics) software is not capable of coupling each of these flow regimes to provide an accurate representation of this flow process, necessitating the development of custom software. This software solves the fluid-flow-governing equations in an Eulerian framework, coupled with the particle transport equations that are solved in a Lagrangian framework. It uses a fourth-order explicit Runge-Kutta scheme for temporal integration, an eighth-order central finite differencing scheme for spatial discretization. The non-linear terms in the governing equations are recast in cubic skew symmetric form to reduce aliasing error. The second derivative viscous terms are computed using eighth-order narrow stencils that provide better diffusion for the highest resolved wave numbers. A fourth-order Lagrange interpolation procedure is used to obtain gas-phase variable values at the particle locations.

Bellan, Josette↗

Adaptive Data-based Predictive Control for Short Take-off and Landing (STOL) Aircraft

Data-based Predictive Control is an emerging control method that stems from Model Predictive Control (MPC). MPC computes current control action based on a prediction of the system output a number of time steps into the future and is generally derived from a known model of the system. Data-based predictive control has the advantage of deriving predictive models and controller gains from input-output data. Thus, a controller can be designed from the outputs of complex simulation code or a physical system where no explicit model exists. If the output data happens to be corrupted by periodic disturbances, the designed controller will also have the built-in ability to reject these disturbances without the need to know them. When data-based predictive control is implemented online, it becomes a version of adaptive control. The characteristics of adaptive data-based predictive control are particularly appropriate for the control of nonlinear and time-varying systems, such as Short Take-off and Landing (STOL) aircraft. STOL is a capability of interest to NASA because conceptual Cruise Efficient Short Take-off and Landing (CESTOL) transport aircraft offer the ability to reduce congestion in the terminal area by utilizing existing shorter runways at airports, as well as to lower community noise by flying steep approach and climb-out patterns that reduce the noise footprint of the aircraft. In this study, adaptive data-based predictive control is implemented as an integrated flight-propulsion controller for the outer-loop control of a CESTOL-type aircraft. Results show that the controller successfully tracks velocity while attempting to maintain a constant flight path angle, using longitudinal command, thrust and flap setting as the control inputs.

Barlow, Jonathan Spencer↗

Hypersonic vehicle model and control law development using H(infinity) and micron synthesis

The control system design for a Single Stage To Orbit (SSTO) air breathing vehicle will be central to a successful mission because a precise ascent trajectory will preserve narrow payload margins. The air breathing propulsion system requires the vehicle to fly roughly halfway around the Earth through atmospheric turbulence. The turbulence, the high sensitivity of the propulsion system to inlet flow conditions, the relatively large uncertainty of the parameters characterizing the vehicle, and continuous acceleration make the problem especially challenging. Adequate stability margins must be provided without sacrificing payload mass since payload margins are critical. Therefore, a multivariable control theory capable of explicitly including both uncertainty and performance is needed. The H(infinity) controller in general provides good robustness but can result in conservative solutions for practical problems involving structured uncertainty. Structured singular value mu framework for analysis and synthesis is potentially much less conservative and hence more appropriate for problems with tight margins. An SSTO control system requires: highly accurate tracking of velocity and altitude commands while limiting angle-of-attack oscillations, minimized control power usage, and a stabilized vehicle when atmospheric turbulence and system uncertainty are present. The controller designs using H(infinity) and mu-synthesis procedures were compared. An integrated flight/propulsion dynamic mathematical model of a conical accelerator vehicle was linearized as the vehicle accelerated through Mach 8. Vehicle acceleration through the selected flight condition gives rise to parametric variation that was modeled as a structured uncertainty. The mu-analysis approach was used in the frequency domain to conduct controller analysis and was confirmed by time history plots. Results demonstrate the inherent advantages of the mu framework for this class of problems.

Gregory, Irene M.↗

Line Emission from an Accretion Disk Around a Rotating Black Hole: Toward a Measurement of Frame Dragging

Line emission from an accretion disk and a corotating hot spot about a rotating black hole are considered for possible signatures of the frame-dragging effect. We explicitly compare integrated line profiles from a geometrically thin disk about a Schwarzschild and an extreme Kerr black hole, and show that the line profile differences are small if the inner radius of the disk is near or above the Schwarzschild stable-orbit limit of radius 6GM/sq c. However, if the inner disk radius extends below this limit, as is Possible in the extreme Kerr spacetime, then differences can become significant, especially if the disk emissivity is stronger near the inner regions. We demonstrate that the first three moments of a line profile define a three-dimensional space in which the presence of material at small radii becomes quantitatively evident in broad classes of disk models. In the context of the simple, thin disk paradigm, this moment-mapping scheme suggests formally that the iron line detected by the Advanced Satellite,for Cosmology and Astrophysics mission from MCG --6-30-15 (Tanaka et al.) is approximately 3 times more likely to originate from a disk about a rotating black hole than from a Schwarzschild system. A statistically significant detection of black hole rotation in this way may be achieved after only modest improvements in the quality of data. We also consider light curves and frequency shifts in line emission as a function of time for corotating hot spots in extreme Kerr and Schwarzschild geometries. The frequency-shift profile is a valuable measure of orbital parameters and might possibly be used to detect frame dragging even at radii approaching 6GM/sq c if the inclination angle of the orbital plane is large. The light curve from a hot spot shows differences as well, although these too are pronounced only at large inclination angles.

Bromley, Benjamin C.↗

Anthropometric Measurement Procedures at the Anthropometry and Biomechanics Facility

The Anthropometry and Biomechanics Facility (ABF) at NASA Johnson Space Center (JSC) is the centralized source of anthropometry data at NASA for use in requirement development, engineering design, and verification of human-system integrations such as spacesuits and vehicles. To ensure the collection of consistent, accurate, and valid anthropometry data, ABF has defined anthropometric measurements and anatomical landmarks which are generally based on the Anthropometric Survey for US Army Personnel (ANSUR), but with additional measurements which have been uniquely designed for spacesuit and hardware accommodation purposes. Unlike ANSUR, NASA anthropometry measurements primarily rely on 3D body scanning, a method which allows for explicitly defining the process of collecting measurements and the quality control methods used to ensure measurement reliability. In addition to ABF’s anthropometry measurement definitions and procedures, this presentation will also compare specifics of different measurement techniques ABF uses, for example, laser scanning versus manual measurement methods (e.g., tape measure or anthropometer-based), and the benefits and considerations of each method. The evolution of ABF’s measurement process will also be addressed, as new spaceflight programs have been introduced which have required unique measurements and the need to capture body shape changes over time was identified to be critical for spacesuit fit in some instances.

Garima Gupta↗

Anthropometric Measurement Procedures at the Anthropometry and Biomechanics Facility

The Anthropometry and Biomechanics Facility (ABF) at NASA Johnson Space Center (JSC) is the centralized source of anthropometry data at NASA for use in requirement development, engineering design, and verification of human-system integrations such as spacesuits and vehicles. To ensure the collection of consistent, accurate, and valid anthropometry data, ABF has defined anthropometric measurements and anatomical landmarks which are generally based on the Anthropometric Survey for US Army Personnel (ANSUR), but with additional measurements which have been uniquely designed for spacesuit and hardware accommodation purposes. Unlike ANSUR, NASA anthropometry measurements primarily rely on 3D body scanning, a method which allows for explicitly defining the process of collecting measurements and the quality control methods used to ensure measurement reliability. In addition to ABF’s anthropometry measurement definitions and procedures, this presentation will also compare specifics of different measurement techniques ABF uses, for example, laser scanning versus manual measurement methods (e.g., tape measure or anthropometer-based), and the benefits and considerations of each method. The evolution of ABF’s measurement process will also be addressed, as new spaceflight programs have been introduced which have required unique measurements and the need to capture body shape changes over time was identified to be critical for spacesuit fit in some instances.

Garima Gupta↗

Mapping National Forest Aboveground Biomass in Mexico By Integrating GEDI and Landsat Times Series Data

Mexico is one of the countries with great potential for the UN's Reducing Emissions from Deforestation and Forest Degradation (REDD+) program, a key nature-based solution for the forest sector. To monitor carbon stock changes, there is a growing demand for unbiased Monitoring Reporting Verification (MRV) systems to facilitate effective forest management and climate change mitigation strategies. Remote sensing-based national aboveground biomass density (AGBD) estimation over Mexico is scarce and often limited to one-time static mapping, leading to spatiotemporal inconsistency in inputs. As an effort under NASA's Carbon Monitoring System (CMS) program, we have developed a remote sensing-based approach to create consistent historical AGBD maps of Mexico using multi-stream remote sensing data, including spaceborne lidar GEDI and long-term Landsat time series, as well as topographic information. We employ the continuous change detection and classification (CCDC) algorithm for temporal modeling of Landsat surface reflectance, followed by the inference of forest AGBD using a random forest machine learning algorithm with the temporal information of land surface dynamics extracted by the CCDC as input. GEDI provides unprecedented forest structure and AGBD sampling datasets for model training and validation practices. In this presentation, we share the progress made in developing a spatially explicit mapping of historical AGBD changes associated with land surface changes and post-disturbance landscapes.

Taejin Park↗

Some considerations on velocity vector accuracy in dust trajectory analysis

The relative contributions of comets and asteroids to the reservoir of dust in the interplanetary medium is not known. There are direct observations of dust released from comets and there is evidence to associate the IRAS dust bands with possible collisions of asteroids in the main belt. A means towards sorting out the parent sources has been proposed in the establishment of a dust collector in orbit about the Earth. The purpose of such a facility would be to collect not only cosmic dust particles intact but also the state vectors, as they arrive at the detector, the idea being that one may combine analytical laboratory analysis of the physics and chemistry of the captured particles with orbital data in order to help distinguish between bodies and identify parent bodies. The theoretical study of dust particle orbits in the solar system takes on greatly more importance if we use collected trajectory data. The orbital motion of dust when radiation and forces alone are acting is well understood. When gravitational forces due to the planets are included, the motion can become quite complex. In order to characterize the orbits of particles as they crossed the Earth's orbits, a study of the long-time dust orbital evolution was undertaken. We have considered various parameters associated with these dust orbits to see if one may in a general way discriminate between particles evolved from comets and asteroids. We proceed in this study as we have done previously. That is, we considered the dust particles as ideal black bodies, of density 1 gm/cc, spherical, with radii 10-100 microns. Particles of this size are affected by radiation forces, photon pressure, and Poynting-Robertson drag. Account was also taken of solar wind drag, which amounts to about 30 percent of the Poynting-Robertson drag negligible. The gravitational forces due to the planets are included, unlike in our previous study; the planetary orbits are those of true n-body interaction so that the possibility of secular resonance is included. Our method was to calculate explicitly by a numerical procedure the orbits of dust particles after they left their parent bodies. The motion is determined numerically with the implicit Runge-Kutta integrator using Gauss-Radau spacings.

Jackson, A. A.↗

Current capabilities for simulating the extreme distortion of thin structures subjected to severe impacts

The explicit transient dynamics technology in use today for simulating the impact and subsequent transient dynamic response of a structure has its origins in the 'hydrocodes' dating back to the late 1940's. The growth in capability in explicit transient dynamics technology parallels the growth in speed and size of digital computers. Computer software for simulating the explicit transient dynamic response of a structure is characterized by algorithms that use a large number of small steps. In explicit transient dynamics software there is a significant emphasis on speed and simplicity. The finite element technology used to generate the spatial discretization of a structure is based on a compromise between completeness of the representation for the physical processes modelled and speed in execution. That is, since it is expected in every calculation that the deformation will be finite and the material will be strained beyond the elastic range, the geometry and the associated gradient operators must be reconstructed, as well as complex stress-strain models evaluated at every time step. As a result, finite elements derived for explicit transient dynamics software use the simplest and barest constructions possible for computational efficiency while retaining an essential representation of the physical behavior. The best example of this technology is the four-node bending quadrilateral derived by Belytschko, Lin and Tsay. Today, the speed, memory capacity and availability of computer hardware allows a number of the previously used algorithms to be 'improved.' That is, it is possible with today's computing hardware to modify many of the standard algorithms to improve their representation of the physical process at the expense of added complexity and computational effort. The purpose is to review a number of these algorithms and identify the improvements possible. In many instances, both the older, faster version of the algorithm and the improved and somewhat slower version of the algorithm are found implemented together in software. Specifically, the following seven algorithmic items are examined: the invariant time derivatives of stress used in material models expressed in rate form; incremental objectivity and strain used in the numerical integration of the material models; the use of one-point element integration versus mean quadrature; shell elements used to represent the behavior of thin structural components; beam elements based on stress-resultant plasticity versus cross-section integration; the fidelity of elastic-plastic material models in their representation of ductile metals; and the use of Courant subcycling to reduce computational effort.

Key, Samuel W.↗

Direct Fault Tolerant RLV Altitude Control: A Singular Perturbation Approach

In this paper, we present a direct fault tolerant control (DFTC) technique, where by "direct" we mean that no explicit fault identification is used. The technique will be presented for the attitude controller (autopilot) for a reusable launch vehicle (RLV), although in principle it can be applied to many other applications. Any partial or complete failure of control actuators and effectors will be inferred from saturation of one or more commanded control signals generated by the controller. The saturation causes a reduction in the effective gain, or bandwidth of the feedback loop, which can be modeled as an increase in singular perturbation in the loop. In order to maintain stability, the bandwidth of the nominal (reduced-order) system will be reduced proportionally according to the singular perturbation theory. The presented DFTC technique automatically handles momentary saturations and integrator windup caused by excessive disturbances, guidance command or dispersions under normal vehicle conditions. For multi-input, multi-output (MIMO) systems with redundant control effectors, such as the RLV attitude control system, an algorithm is presented for determining the direction of bandwidth cutback using the method of minimum-time optimal control with constrained control in order to maintain the best performance that is possible with the reduced control authority. Other bandwidth cutback logic, such as one that preserves the commanded direction of the bandwidth or favors a preferred direction when the commanded direction cannot be achieved, is also discussed. In this extended abstract, a simplistic example is proved to demonstrate the idea. In the final paper, test results on the high fidelity 6-DOF X-33 model with severe dispersions will be presented.

Zhu, J. J.↗

A Parallel Compact Multi-Dimensional Numerical Algorithm with Aeroacoustics Applications

In this study we propose a novel method to parallelize high-order compact numerical algorithms for the solution of three-dimensional PDEs (Partial Differential Equations) in a space-time domain. For this numerical integration most of the computer time is spent in computation of spatial derivatives at each stage of the Runge-Kutta temporal update. The most efficient direct method to compute spatial derivatives on a serial computer is a version of Gaussian elimination for narrow linear banded systems known as the Thomas algorithm. In a straightforward pipelined implementation of the Thomas algorithm processors are idle due to the forward and backward recurrences of the Thomas algorithm. To utilize processors during this time, we propose to use them for either non-local data independent computations, solving lines in the next spatial direction, or local data-dependent computations by the Runge-Kutta method. To achieve this goal, control of processor communication and computations by a static schedule is adopted. Thus, our parallel code is driven by a communication and computation schedule instead of the usual "creative, programming" approach. The obtained parallelization speed-up of the novel algorithm is about twice as much as that for the standard pipelined algorithm and close to that for the explicit DRP algorithm.

Povitsky, Alex↗

An Interoperability Concept for Detect and Avoid and Collision Avoidance Systems: Results from a Human-In-The-Loop Simulation

The integration of Unmanned Aircraft Systems (UAS) into the National Airspace System (NAS) poses a variety of technical challenges to UAS developers and aviation regulators. In response to growing demand for access to civil airspace in the United States, the Federal Aviation Administration (FAA) has produced a roadmap identifying key areas requiring further research and development. One such technical challenge is the development of a "detect and avoid" system (DAA) capable of providing a means of compliance with the "see and avoid" requirement in manned aviation. The purpose of the DAA system is to support the pilot, situated at a ground control station (GCS), in maintaining "DAA well clear" of nearby aircraft through the use of GCS displays and alerts. In addition to its primary function of aiding the pilot in maintaining DAA well clear, the DAA system must also safely interoperate with existing NAS systems and operations, such as the airspace management procedures of air traffic controllers (ATC) and Collision Avoidance (CA) systems currently in use by manned aircraft, namely the Traffic Alert and Collision Avoidance System (TCAS II). It is anticipated that many UAS architectures will integrate both a DAA system and a TCAS II. It is therefore necessary to explicitly study the integration of DAA and TCAS II alerting structures and maneuver guidance formats to ensure that pilots understand the appropriate type and urgency of their response to the various alerts. This paper presents a concept of interoperability for the two systems. The concept was developed with the goal of avoiding any negative impact on the performance level of TCAS II while retaining a DAA system that still effectively enables pilots to maintain DAA well clear. The interoperability concept described in the paper focuses primarily on facilitating the transition from a late-stage DAA encounter (where a loss of DAA well clear is imminent) to a TCAS II Corrective Resolution Advisory (RA), which requires pilot compliance within five seconds of its issuance. The interoperability concept was presented to 10 participants (6 active UAS pilots and 4 active commercial pilots) in a medium-fidelity, human-in-the-loop simulation designed to stress different aspects of the DAA and TCAS II systems. Pilots' ability to maintain separation, their rate of compliance and response times using the interoperability concept are reported. Results indicated that pilots exhibited comprehension of, and appropriate prioritization within, the DAA-TCAS II combined alert structure. Pilots demonstrated a high rate of compliance with TCAS II RAs and were also seen to respond to corrective RAs within the five second requirement established for manned aircraft. The DAA system presented under test was also shown to be effective in supporting pilots' ability to maintain DAA well clear in the overwhelming majority of cases in which pilots had sufficient time to respond.

Rorie, R. Conrad↗

Euler and Navier-Stokes solutions for hypersonic flows

An upwind finite-element technique that uses cell-centered quantities and implicit and/or explicit time marching has been developed for computing hypersonic laminar viscous flows using adaptive unstructured grids in two and three dimensions. A perfect gas model as well as an equilibrium air model is implemented for solving high-speed flows. A first-order basic scheme and a higher-order flux-corrected transport (FCT) scheme have been implemented. This technique has been used to predict 'Type III and IV' shock interactions on a cylinder in two dimensions and a swept cylinder in three dimensions, with a view to determine the pressure and heating rate augmentation caused by an impinging shock on the leading edge of a cowl lip of an engine inlet. The predictions of wall pressure and heating rates compare very well with experimental data. The flow features are very distinctly captured with a sequence of adaptively-generated grids. Three-dimensional corner flow, typically encountered in engine inlets due to compression of the flow by ramps in the walls, is also modeled. This procedure is the first step in developing an integrated fluid, thermal, structural analysis capability for hypersonic flight vehicles like the National Aero-Space Plane.

Thareja, Rajiv R.↗

The mesoscale forcing of a midlatitude upper-tropospheric jet streak by a simulated convective system. 2: Kinetic energy and resolution analysis

A kinetic energy (KE) analysis of the forcing of a mesoscale upper-tropospheric jet streak by organized diabatic processes within the simulated convective system (SCS) that was discussed in Part 1 is presented in this study. The relative contributions of the ageostrophic components of motion to the generation of KE of the convectively generated jet streak are compared, along with the KE generation by the rotational (nondivergent) and irrotational (divergent) mass transport. The sensitivity of the numerical simulations of SCS development to resolution is also briefly examined. Analysis within isentropic coordinates provides for an explicit determination of the influence of the diabatic processes on the generation of KE. The upper-level production of specific KE is due predominatly to the inertial advective ageostrophic component (IAD), and as such represents the primary process through which the KE of the convectively generated jet streak is realized. A secondary contribution by the inertial diabatic (IDI) term is observed. Partitioning the KE generation into its rotational and irrotational components reveals that the latter, which is directly linked to the diabatic heating within the SCS through isentropic continuity requirements, is the ultimate source of KE generation as the global area integral of generation by the rotational component vanishes. Comparison with an identical dry simulation reveals that the net generation of KE must be attributed to latent heating. Both the IAD and IDI ageostrophic components play important roles in this regard. Examination of results from simulations conducted at several resolutions supports the previous findings in that the effects of diabatic processes and ageostrophic motion on KE generation remain consistent. Resolution does impact the location and timing of SCS development, a result that has important implications in forecasting the onset of convection that develops from evolution of the large-scale flow and moisture transport. Marked differences are observed in the momentum field aloft subsequent to the life cycle of the SCS in the 1 deg, 30-level base case (MP130) simulation discussed in Part 1 versus its 2 deg counterparts in that the MP130 simulation with higher spatial resolution contains from 14% to 30% more total KE.

Wolf, Bart J.↗