Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel simulation”

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 1,621 records · Page 90

Bounce-averaged Hamiltonian for charged particles in an axisymmetric but nondipolar model magnetosphere

In order to facilitate bounce-averaged guiding center simulations of geomagnetically trapped particles, we express the kinetic energy of a particle with magnetic coordinates (L,phi) as an analytic function of the first two adiabatic invariants (M,J) and the L value of the field line. The magnetic field model is axisymmetric, consisting of a dipolar vector-B field plus a uniform southward magnetic field parallel to the dipole moment mu(sub E). This model magnetosphere is surrounded by a circular equatorial neutral line whose radius b is an adjustable parameter. Our formulation provides a computationally efficient method for tracing the bounce-averaged adiabatic motion (conserving all three invariants) and nonadiabatic transport (violating the third invariant while conserving the first two invariants) of geomagnetically trapped particles in the model magnetosphere.

Schulz, Michael↗

Genetically Engineered Microelectronic Infrared Filters

A genetic algorithm is used for design of infrared filters and in the understanding of the material structure of a resonant tunneling diode. These two components are examples of microdevices and nanodevices that can be numerically simulated using fundamental mathematical and physical models. Because the number of parameters that can be used in the design of one of these devices is large, and because experimental exploration of the design space is unfeasible, reliable software models integrated with global optimization methods are examined The genetic algorithm and engineering design codes have been implemented on massively parallel computers to exploit their high performance. Design results are presented for the infrared filter showing new and optimized device design. Results for nanodevices are presented in a companion paper at this workshop.

Cwik, Tom↗

Accelerated Simulation of Air Pollution Using NVIDIA RAPIDS

Atmospheric chemistry models are a central tool to study and forecast the impact of air pollution on the environment, vegetation, and human health. However, the numerical simulation of chemical kinetics is computationally expensive due to the stiffness of the system of ordinary differential equations that describes atmospheric chemistry. Here we present an alternative approach to the computation of atmospheric chemistry based on machine learning. Our training data set is produced using the NASA Goddard Earth Observing System (GEOS) model with GEOS-Chem chemistry, run on the NASA Center for Climate Simulation (NCCS) Discover supercomputing cluster on 384 Intel Xeon Haswell cores. This model spends more than 50% of total run time on solving atmospheric chemistry. The data set contains as input features the air pollution concentrations before solving the differential equations, together with some key physical parameters such as temperature and sun intensity. As target variables we define the air pollution concentrations after solving the differential equations. Using Dask-cuDF and Dask-XGBoost on the NVIDIA RAPIDS platform on 8 Tesla V100 GPUs, we generate from this training set gradient boosted decision tree models that can reproduce the simulation of chemical kinetics. We do this on the NCCS Advanced Data Analytics Platform (ADAPT) science cloud environment. Our application takes full advantage of recent advances in Dask-XGBoost, such as multi-node and multi-GPU scaling for distributed training with large data sets. The increase in training data size enabled by this is critical to capture the full range of chemical environments encountered across the globe and all annual seasons.The boosted tree models offer good predictability and show many of the features of the full chemistry reference simulation. Further improvements can be achieved through mass balance considerations and by accounting for error correlations. We incorporate the boosted tree models into the GEOS reference model using XGBoost's C API. This enables a seamless integration of the GPU trained models into GEOS-Chem, which is written in Fortran and optimized for use in a massively parallel CPU environment. We show the benefits of this approach and discuss the potential speedup of this machine learning accelerated atmospheric chemistry model.

Keller, Christoph A.↗

Liquid Oxygen/Liquid Methane Test Results of the RS-18 Lunar Ascent Engine at Simulated Altitude Conditions at NASA White Sands Test Facility

Tests were conducted with the RS-18 rocket engine using liquid oxygen (LO2) and liquid methane (LCH4) propellants under simulated altitude conditions at NASA Johnson Space Center White Sands Test Facility (WSTF). This project is part of NASA's Propulsion and Cryogenics Advanced Development (PCAD) project. "Green" propellants, such as LO2/LCH4, offer savings in both performance and safety over equivalently sized hypergolic propulsion systems in spacecraft applications such as ascent engines or service module engines. Altitude simulation was achieved using the WSTF Large Altitude Simulation System, which provided altitude conditions equivalent up to ~122,000 ft (~37 km). For specific impulse calculations, engine thrust and propellant mass flow rates were measured. LO2 flow ranged from 5.9 - 9.5 lbm/sec (2.7 - 4.3 kg/sec), and LCH4 flow varied from 3.0 - 4.4 lbm/sec (1.4 - 2.0 kg/sec) during the RS-18 hot-fire test series. Propellant flow rate was measured using a coriolis mass-flow meter and compared with a serial turbine-style flow meter. Results showed a significant performance measurement difference during ignition startup due to two-phase flow effects. Subsequent cold-flow testing demonstrated that the propellant manifolds must be adequately flushed in order for the coriolis flow meters to give accurate data. The coriolis flow meters were later shown to provide accurate steady-state data, but the turbine flow meter data should be used in transient phases of operation. Thrust was measured using three load cells in parallel, which also provides the capability to calculate thrust vector alignment. Ignition was demonstrated using a gaseous oxygen/methane spark torch igniter. Test objectives for the RS-18 project are 1) conduct a shakedown of the test stand for LO2/methane lunar ascent engines, 2) obtain vacuum ignition data for the torch and pyrotechnic igniters, and 3) obtain nozzle kinetics data to anchor two-dimensional kinetics codes. All of these objectives were met with the RS-18 data and additional testing data from subsequent LO2/methane test programs in 2009 which included the first simulated-altitude pyrotechnic ignition demonstration of LO2/methane.

Melcher, John C., IV↗

Three methods of presenting flight vector information in a head-up display during simulated STOL approaches

A simulator study was conducted to determine the usefulness of adding flight path vector symbology to a head-up display designed to improve glide-slope tracking performance during steep 7.5 deg visual approaches in STOL aircraft. All displays included a fixed attitude symbol, a pitch- and roll-stabilized horizon bar, and a glide-slope reference bar parallel to and 7.5 deg below the horizon bar. The displays differed with respect to the flight-path marker (FPM) symbol: display 1 had no FPM symbol; display 2 had an air-referenced FPM, and display 3 had a ground-referenced FPM. No differences between displays 1 and 2 were found on any of the performance measures. Display 3 was found to decrease height error in the early part of the approach and to reduce descent rate variation over the entire approach. Two measures of workload did not indicate any differences between the displays.

Dwyer, J. H., III↗

Calculation of secondary-electron escape currents from inclined-spacecraft surfaces in a magnetic field

In low Earth orbit, the geomagnetic field B(vector) is strong enough that secondary electrons emitted from spacecraft surfaces have an average gyroradius much smaller than typical dimensions of large spacecraft. This implies that escape of secondaries will be strongly inhibited on surfaces which are nearly parallel to B(vector), even if a repelling electric field exists outside them. This effect is likely to make an important contribution to the current balance and hence the equilibrium potential of such surfaces, making high voltage charging of them more likely. Numerically calculated escaping secondary electron fluxes are presented for these conditions. For use in numerical spacecraft charging simulations, an analytic curve fit to these results is given which is accurate to within 3% of the emitted current.

Laframboise, J. G.↗

Transient aerodynamic forces on a fighter model during simulated approach and landing with thrust reversers

Previous wind tunnel tests of fighter configurations have shown that thrust reverser jets can induce large, unsteady aerodynamic forces and moments during operation in ground proximity. This is a concern for STOL configurations using partial reversing to spoil the thrust while keeping the engine output near military (MIL) power during landing approach. A novel test technique to simulate approach and landing was developed under a cooperative Northrop/NASA/USAF program. The NASA LaRC Vortex Research Facility was used for the experiments in which a 7-percent F-18 model was moved horizontally at speeds of up to 100 feet per second over a ramp simulating an aircraft to ground rate of closure similar to a no-flare STOL approach and landing. This paper presents an analysis of data showing the effect of reverser jet orientation and jet dynamic pressure ratio on the transient forces for different angles of attack, and flap and horizontal tail deflection. It was found, for reverser jets acting parallel to the plane of symmetry, that the jets interacted strongly with the ground, starting approximately half a span above the ground board. Unsteady rolling moment transients, large enough to cause the probable upset of an aircraft, and strong normal force and pitching moment transients were measured. For jets directed 40 degrees outboard, the transients were similar to the jet-off case, implying only minor interaction.

Humphreys, A. P.↗

The interaction of collisionless shocks in astrophysical plasmas

The interaction between collisionless shocks, such as may arise in a number of situations in space plasmas, is examined by means of hybrid numerical simulations. The production of energetic particles by shock collisions is investigated. When the shocks are quasi-perpendicular, ions with energies greater than 10E(0) are produced, where E(0) is the kinetic energy of the unshocked ions. The effect is optimized when the two colliding shocks have roughly equal strengths. As the shocks become closer to quasi-parallel, the number and energy of the accelerated ions increases. Energies in excess of 30E(0) are obtained in the quasi-parallel regime, with a small fall-off when the shocks are exactly parallel. For quasi-parallel collisions, the ion energization can be attributed to the interaction of hot ions with both strong electromagnetic waves and the motional electric field located at the shocks. These results are applied to shock collisions in the distant solar wind, at planetary bow shocks, as well as to the impulsive phase of solar flares.

Cargill, P. J.↗

Extended ecosystem signatures with application to Eos synergism requirements

The primary objective is to define the advantages of synergistically combining optical and microwave remote sensing measurements for the determination of biophysical properties important in ecosystem modeling. This objective was approached in a stepwise fashion starting with ground-based observations of controlled agricultural and orchard canopies and progressing to airborne observations of more natural forest ecosystems. This observational program is complemented by a parallel effort to model the visible reflectance and microwave scattering properties of composite vegetation canopies. The goals of the modeling studies are to verify our basic understanding of the sensor-scene interaction physics and to provide the basis for development of inverse models optimized for retrieval of key biophysical properties. These retrieval algorithms can then be used to simulate the expected performance of various aspects of Eos including the need for simultaneous SAR and HIRIS observations or justification for other (non-synchronous) relative timing constraints and the frequency, polarization, and angle of incidence requirements for accurate biophysical parameter extractions. This program completed a very successful series of truck-mounted experiments, made remarkable progress in development and validation of optical reflectance and microwave scattering models for vegetation, extended the scattering models to accommodate discontinuous and periodic canopies, developed inversion approaches for surface and canopy properties, and disseminated these results widely through symposia and journal publications. In addition, the third generation of the computer code for the microwave scattering models was provided to a number of other US, Canadian, Australian, and European investigators who are currently presenting and publishing results using the MIMICS research code.

Ulaby, Fawwaz T.↗

Laboratory Astrophysics using a Microcalorimeter and Bragg Crystal Spectrometer on an Electron Beam Ion Trap

During the past year we have been preparing our new microcalorimeter system for permanent delivery to the NIST EBIT. Unfortunately, there have been delays due to technical difficulties in the fabrication of the two-stage adiabatic demagnetization refrigerator and in the life expectancy of the thin windows used for internal thermal baffling of the infrared radiation. These problems have been solved and we are completing tests of the entire system and it will be set up at NIST during the first week of May. Several photos of the new system are shown in Figures 1A and 1B. This microcalorimeter spectrometer only requires helium refills every three days (as opposed to every 24 hours) and it will hold a temperature! of 65 mK for up to 48 hours (as opposed to 8 hours). Consequently, the efficiency of data acquisition will improve dramatically. In parallel we have published a paper that reviews our previous work (Takacs et al. 2003), especially on Fe XVII, in the context of recent measurements by other groups. This paper is included. We highlight a recent measurement of a broad band spectrum of Fe in Figure 2 that simultaneously includes L and K radiation. It is compared with the simulated spectrum of the Perseus Cluster that one could expect to obtain with a microcalorimeter in the focus of a grazing incidence telescope such as the one being designed for Constellation X. Both the charge state distributions and the relative intensity ratios of the emission lines within the particular charge state are very similar in the two spectra. This further demonstrates the importance and relevance of the laboratory measurements in predicting the components of cosmic spectra.

Brinton, John C.↗

Recent Improvements to the LAURA and HARA Codes

This paper describes recent improvements to the LAURA and HARA codes. LAURA is a CFD code for aerothermodynamics, and HARA evaluates the shock-layer radiation that provides the radiative source term for the flowfield energy equations and radiative heating to a surface. The next release of LAURA and HARA includes a variety of new capabilities. These new capabilities include an automated uncertainty quantification workflow for radiative heat transfer, options for specifying surface roughness and turbulent transition location in the algebraic turbulence models, and improved grid and solution interpolation techniques. Additionally, the computational efficiency of both LAURA and HARA have been improved. Optimization of the MPI communication routines in LAURA are shown to improve the parallel efficiency of the primary flow when running with multiple processes per block, and recent optimization of HARA leverage graphics processing unit (GPU) acceleration in the radiation calculations. Using GPU acceleration of HARA is shown to decrease the cost of the radiation line-of-sight calculation by approximately one order of magnitude for a 10.5 km/s Earth entry simulation.

LAURA HARA CFD 5.6↗

The design and implementation of CRT displays in the TCV real-time simulation

The design and application of computer graphics to the Terminal Configured Vehicle (TCV) program were described. A Boeing 737-100 series aircraft was modified with a second flight deck and several computers installed in the passenger cabin. One of the elements in support of the TCV program is a sophisticated simulation system developed to duplicate the operation of the aft flight deck. This facility consists of an aft flight deck simulator, equipped with realistic flight instrumentation, a CDC 6600 computer, and an Adage graphics terminal; this terminal presents to the simulator pilot displays similar to those used on the aircraft with equivalent man-machine interactions. These two displays form the primary flight instrumentation for the pilot and are dynamic images depicting critical flight information. The graphics terminal is a high speed interactive refresh-type graphics system. To support the cockpit display, two remote CRT's were wired in parallel with two of the Adage scopes.

Leavitt, J. B.↗

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

The primary objective of this research is to extend current capabilities of Large Eddy Simulations (LES) and Direct Numerical Simulations (DNS) for the computational analyses of high speed reacting flows. Our efforts in the first two years of this research have been concentrated on a priori investigations of single-point Probability Density Function (PDF) methods for providing subgrid closures in reacting turbulent flows. In the efforts initiated in the third year, our primary focus has been on performing actual LES by means of PDF methods. The approach is based on assumed PDF methods and we have performed extensive analysis of turbulent reacting flows by means of LES. This includes simulations of both three-dimensional (3D) isotropic compressible flows and two-dimensional reacting planar mixing layers. In addition to these LES analyses, some work is in progress to assess the extent of validity of our assumed PDF methods. This assessment is done by making detailed companions with recent laboratory data in predicting the rate of reactant conversion in parallel reacting shear flows. This report provides a summary of our achievements for the first six months of the third year of this program.

Givi, P.↗

A Study on the Effect of Communication Performance on Message-Passing Parallel Programs: Methodology and Case Studies

From a source-program perspective, the performance achieved on distributed/parallel systems is governed by the underlying message-passing library overhead and the network capabilities of the architecture. Studying the impact of changes in these features on the source-program. can have a significant influence in the development of next-generation system designs. In this paper we introduce a simple and robust tool that can be used for this purpose. This tool is based on event-driven simulation of programs that generates a new set of trace events - that preserves causality and partial order - corresponding to the expected execution of the program in the simulated environment. Trace events can be visualized and source-level profile information can be used to pin-point locations of program which are most significantly affected with changing system parameters in the simulated environment. We present a number of examples from the NAS benchmark suite, executed on the Intel Paragon and iPSC/860 that are used to identify and expose performance bottlenecks with varying system parameters. Specific aspects of the system that significantly effect these benchmarks are presented and discussed,

Sarukkai, Sekhar R.↗

Electromagnetic ion beam instabilities - Comparison of oneand two-dimensional simulations

Numerical simulations of electromagnetic instabilities, driven by a cool tenuous ion beam propagating along an ambient magnetic field, have been conducted in one and two spatial dimensions. The calculations employ particle ions, fluid electrons, and a predictor corrector scheme for solving the electromagnetic field in two dimensions that is described in some detail. While the principal features of the one-dimensional calculations (which reproduce previous work) are retained, the two-dimensional simulations show some reduction of the overall level of the magnetic field fluctuations. Enhancement of the heating of the beam ions at the expense of the core ions also occurs in the case where the beam density is sufficiently large that the right-hand nonresonant instability dominates. Implications of the results for modeling of the ion foreshock and quasi-parallel shocks are discussed.

Winske, D.↗

MIT research in telerobotics

Ongoing MIT research in telerobotics (vehicles capable of some autonomous sensing and manipulating, having some remote supervisory control by people) and teleoperation (vehicles for sensing and manipulating which are fully controlled remotely by people) is discussed. The current efforts mix human and artificial intelligence/control. The idea of adjustable impedance at either end of pure master-slave teleoperation, and simultaneous coordinated control of teleoperator/telerobotic systems which have more than six degrees of freedom (e.g., a combined vehicle and arm, each with five or six DOF) are discussed. A new cable-controlled parallel link arm which offers many advantages over conventional arms for space is briefly described. Predictor displays to compensate for time delay in teleoperator loops, the use of state estimation to help human control decisions in space, and ongoing research in supervisory command language are covered. Finally, efforts to build a human flyable real-time dynamic computer-graphic telerobot simulator are described. These projects represent most, but not all, of the telerobotics research in our laboratory, supported by JPL, NASA Ames and NOAA.

Sheridan, T. B.↗

Adaptive Strategies for Controls of Flexible Arms

An adaptive controller for a modern manipulator has been designed based on asymptotical stability via the Lyapunov criterion with the output error between the system and a reference model used as the actuating control signal. Computer simulations were carried out to test the design. The combination of the adaptive controller and a system vibration and mode shape estimator show that the flexible arm should move along a pre-defined trajectory with high-speed motion and fast vibration setting time. An existing computer-controlled prototype two link manipulator, RALF (Robotic Arm, Large Flexible), with a parallel mechanism driven by hydraulic actuators was used to verify the mathematical analysis. The experimental results illustrate that assumed modes found from finite element techniques can be used to derive the equations of motion with acceptable accuracy. The robust adaptive (modal) control is implemented to compensate for unmodelled modes and nonlinearities and is compared with the joint feedback control in additional experiments. Preliminary results show promise for the experimental control algorithm.

Yuan, Bau-San↗

Parallel arrays of Josephson junctions for submillimeter local oscillators

In this paper we discuss the influence of the DC biasing circuit on operation of parallel biased quasioptical Josephson junction oscillator arrays. Because of nonuniform distribution of the DC biasing current along the length of the bias lines, there is a nonuniform distribution of magnetic flux in superconducting loops connecting every two junctions of the array. These DC self-field effects determine the state of the array. We present analysis and time-domain numerical simulations of these states for four biasing configurations. We find conditions for the in-phase states with maximum power output. We compare arrays with small and large inductances and determine the low inductance limit for nearly-in-phase array operation. We show how arrays can be steered in H-plane using the externally applied DC magnetic field.

Pance, Aleksandar↗