Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accelerated computing”

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,153 records · Page 64

Validation of a Flexible Aircraft Take-Off and Landing Analysis (FATOLA)

Modifications to improve the analytical simulation capabilities of a multi-degree-of-freedom flexible aircraft take-off and landing analysis (FATOLA) computer program are discussed. The FATOLA program was used to simulate the landing behavior of a stiff body X-24B reentry research vehicle and of a flexible body supersonic cruise YF-12A research airplane. The analytical results were compared with flight test data, and correlations of vehicle motions, attitudes, forces, and accelerations during the landing impact and rollout were good. For the YF-12A airplane, airframe flexibility was found to be important for nose gear loading. Based upon the correlation study presented, the versatility and validity of the FATOLA program for the study of landing dynamics of aircraft are confirmed.

Carden, H. D.↗

Neutron Star Structure in the Presence of Scalar Fields

Motivated by the possible presence of scalar fields on astrophysical scales, suggested by the apparent acceleration of the universe implied by the supernovae surveys, we present models of neutron star structure including the contribution of a (massless) scalar field to the stress energy momentum tensor, in addition to that made by the normal matter. To that end we solve the coupled Einstein -- scalar field -- hydrostatic balance equations to compute the effect of the presence of the scalar field on the neutron star structure. We find that the presence of the scalar field does change the structure of the neutron star, especially in cases of strong coupling between the scalar field and the matter density. We present the neutron star radius as a function of the matter--scalar field coupling constant for different values of the neutron star central density. The presence of the scalar field affects both the maximum neutron star mass and Its radius, the latter increasing with the value of the above coupling constant. We also compute particle and photon geodesics in the geometry of these neutron stars as well as to the geometry of black holes with different values of the scalar field. Our results may be testable with timing observations of accreting neutron stars.

Kazanas, Demosthenes↗

Computational Design of Eutectic Molten Salt Mixtures: What Can Thermodynamic Models Do?

In the search for efficient energy storage battery technologies, designing stable electrolytes has been a long-standing challenge. Electrolytes based on molten salt eutectics are known for their stability with minimum parasitic reactions when compared to their widely used organic counterparts. However, the operating temperatures of these molten salt electrolyte-based batteries are dictated by the melting point of the eutectic mixtures. Design and high throughput screening of low melting temperature eutectic molten salt mixtures have been hindered by the lack of computational models. In this work, we develop thermodynamic models to predict the eutectic points of several molten salt mixtures. The framework of the COSMO-SAC model is used for the predictions and is compared with experimental data and other thermodynamic approaches. Rapid thermodynamics-based approaches, as shown in this study, can accelerate the discovery of new materials, complementing experimental techniques.

Ashwin Ravichandran↗

Development of Fast Algorithms Using Recursion, Nesting and Iterations for Computational Electromagnetics

In the first phase of our work, we have concentrated on laying the foundation to develop fast algorithms, including the use of recursive structure like the recursive aggregate interaction matrix algorithm (RAIMA), the nested equivalence principle algorithm (NEPAL), the ray-propagation fast multipole algorithm (RPFMA), and the multi-level fast multipole algorithm (MLFMA). We have also investigated the use of curvilinear patches to build a basic method of moments code where these acceleration techniques can be used later. In the second phase, which is mainly reported on here, we have concentrated on implementing three-dimensional NEPAL on a massively parallel machine, the Connection Machine CM-5, and have been able to obtain some 3D scattering results. In order to understand the parallelization of codes on the Connection Machine, we have also studied the parallelization of 3D finite-difference time-domain (FDTD) code with PML material absorbing boundary condition (ABC). We found that simple algorithms like the FDTD with material ABC can be parallelized very well allowing us to solve within a minute a problem of over a million nodes. In addition, we have studied the use of the fast multipole method and the ray-propagation fast multipole algorithm to expedite matrix-vector multiplication in a conjugate-gradient solution to integral equations of scattering. We find that these methods are faster than LU decomposition for one incident angle, but are slower than LU decomposition when many incident angles are needed as in the monostatic RCS calculations.

Chew, W. C.↗

Role of High-End Computing in Meeting NASA's Science and Engineering Challenges

Two years ago, NASA was on the verge of dramatically increasing its HEC capability and capacity. With the 10,240-processor supercomputer, Columbia, now in production for 18 months, HEC has an even greater impact within the Agency and extending to partner institutions. Advanced science and engineering simulations in space exploration, shuttle operations, Earth sciences, and fundamental aeronautics research are occurring on Columbia, demonstrating its ability to accelerate NASA s exploration vision. This talk describes how the integrated production environment fostered at the NASA Advanced Supercomputing (NAS) facility at Ames Research Center is accelerating scientific discovery, achieving parametric analyses of multiple scenarios, and enhancing safety for NASA missions. We focus on Columbia s impact on two key engineering and science disciplines: Aerospace, and Climate. We also discuss future mission challenges and plans for NASA s next-generation HEC environment.

Biswas, Rupak↗

The solution of radiative transfer problems in molecular bands without the LTE assumption by accelerated lambda iteration methods

An iterative method based on the use of approximate transfer operators, which was designed initially to solve multilevel NLTE line formation problems in stellar atmospheres, is adapted and applied to the solution of the NLTE molecular band radiative transfer in planetary atmospheres. The matrices to be constructed and inverted are much smaller than those used in the traditional Curtis matrix technique, which makes possible the treatment of more realistic problems using relatively small computers. This technique converges much more rapidly than straightforward iteration between the transfer equation and the equations of statistical equilibrium. A test application of this new technique to the solution of NLTE radiative transfer problems for optically thick and thin bands (the 4.3 micron CO2 band in the Venusian atmosphere and the 4.7 and 2.3 micron CO bands in the earth's atmosphere) is described.

Kutepov, A. A.↗

A single axis study of flight simulator kinematics by difference techniques

The kinematic parameters of position, velocity, and acceleration of a flight simulator may be calculated by knowing the distance between two or more points on an axis and the time the simulator takes to traverse the space between each set of points. These parameters are calculated through the use of difference techniques. Given the true kinematic response of the simulator to computer generated commands, the entire motion system loop may be calibrated, and system operability verified.

Webster, L. D.↗

DE 1 observations of theta aurora plasma source regions and Birkeland current charge carriers

Detailed analyses of the DE 1 high-altitude plasma instrument electron and ion data have been performed for four passes during which theta auroras were observed. The data indicate that the theta auroras occur on what appear to be closed field lines with particle signatures and plasma parameters that are quite similar to those of the magnetospheric boundary plasma sheet. The field-aligned currents computed from particle fluxes in the energy range 18-13 keV above the theta auroras are observed to be generally downward on the dawnside of the arcs with a narrower region of larger (higher density) upward currents on the duskside of the arcs. These currents are carried predominantly by field-aligned beams of accelerated cold electrons. Of particualr interest in regions of upward field-aligned current are downward electron beams at energies less than the inferred potential drop above the spacecraft.

Menietti, J. D.↗

Efficient control with an order (n) recursive inversion of the Jacobian for an n-link serial manipulator

A method is presented to order (n) recursively invert the Jacobian matrix for a serial n-link manipulator. By realizing that the Jacobian relationship involving accelerations is a special case for the forward dynamics problem, it is possible to obtain the Jacobian inverse by directly applying Rodriguez's (1987) recursive forward dynamics algorithm. In this technique, the n-link robot equation is formulated as a spatially recursive algorithm in the form of a filtering and smoothing problem. To compute the Jacobian inverse via this algorithm, the n-link manipulator is modeled with only a unit end-point mass and zero applied torques.

Meldrum, D. R.↗

Full field flow visualization and computer-aided velocity measurements in a bank of cylinders in a wind tunnel

The full field flow tracking (FFFT) method that is presented in this paper uses a laser-generated, mechanically strobed planar sheet of light, a low luminosity TV camera coupled with a long distance microscope, and a computer-controlled video recorder to study nonintrusively and qualitatively the flow structures in a bank of cylinders that are placed in a wind tunnel. This setup simulates an upscale version of the geometry of internal cooling passageways characteristic of small air-cooled radial turbines. The qualitative images supplied by the FFFT system are processed by means of a computer-integrated image quantification (CIIQ) method into quantitative information, trajectories and velocities, that describe the flow upstream of and within the bank of cylinders. The tracking method is Lagrangian in concept, and permits identification and tracking of the same particle, thus facilitating construction of time dependent trajectories and the calculation of true velocities and accelerations. The error analysis evaluates the accuracy with which the seed particles follow the flow and the errors incurred during the quantitative processing of the raw data derived from the FFFT/CIIQ method.

Braun, M. J.↗

Multigrid convergence of an implicit symmetric relaxation scheme

The multigrid method has been applied to an existing three-dimensional compressible Euler solver to accelerate the convergence of the implicit symmetric relaxation scheme. This lower-upper symmetric Gauss-Seidel implicit scheme is shown to be an effective multigrid driver in three dimensions. A grid refinement study is performed including the effects of large cell aspect ratio meshes. Performance figures of the present multigrid code on Cray computers including the new C90 are presented. A reduction of three orders of magnitude in the residual for a three-dimensional transonic inviscid flow using 920 k grid points is obtained in less than 4 min on a Cray C90.

Yoon, Seokkwan↗

Analysis of Factors Affecting the Performance of RLV Thrust Cell Liners

The reusable launch vehicle (RLV) thrust cell liner, or thrust chamber, is a critical component of the Space Shuttle Main Engine (SSME). It is designed to operate in some of the most severe conditions seen in engineering practice. This requirement, in conjunction with experimentally observed 'dog-house' failure modes characterized by bulging and thinning of the cooling channel wall, provides the motivation to study the factors that influence RLV thrust cell liner performance. Factors or parameters believed to be directly related to the observed characteristic deformation modes leading to failure under in-service loading conditions are identified, and subsequently investigated using the cylindrical version of the higher-order theory for functionally graded materials in conjunction with the Robinson's unified viscoplasticity theory and the power-law creep model for modeling the response of the liner s constituents. Configurations are analyzed in which specific modifications in cooling channel wall thickness or constituent materials are made to determine the influence of these parameters on the deformations resulting in the observed failure modes in the outer walls of the cooling channel. The application of thermal barrier coatings and functional grading are also investigated within this context. Comparison of the higher-order theory results based on the Robinson and power-law creep model predictions has demonstrated that, using the available material parameters, the power-law creep model predicts more precisely the experimentally observed deformation leading to the 'dog-house' failure mode for multiple short cycles, while also providing much improved computational efficiency. However, for a single long cycle, both models predict virtually identical deformations. Increasing the power-law creep model coefficients produces appreciable deformations after just one long cycle that would normally be obtained after multiple cycles, thereby enhancing the efficiency of the analysis. This provides a basis for the development of an accelerated modeling procedure to further characterize dog-house deformation modes in RLV thrust cell liners. Additionally, the results presented herein have demonstrated that the mechanism responsible for deformation leading to 'dog-house' failure modes is driven by pressure, creep/relaxation and geometric effects.

Arnold, Steven M.↗

Lighting Condition Analysis for Mars Moon Phobos

A manned mission to Phobos may be an important precursor and catalyst for the human exploration of Mars, as it will fully demonstrate the technologies for a successful Mars mission. A comprehensive understanding of Phobos' environment such as lighting condition and gravitational acceleration are essential to the mission success. The lighting condition is one of many critical factors for landing zone selection, vehicle power subsystem design, and surface mobility vehicle path planning. Due to the orbital characteristic of Phobos, the lighting condition will change dramatically from one Martian season to another. This study uses high fidelity computer simulation to investigate the lighting conditions, specifically the solar radiation flux over the surface, on Phobos. Ephemeris data from the Jet Propulsion Laboratory (JPL) DE405 model was used to model the state of the Sun, the Earth, and Mars. An occultation model was developed to simulate Phobos' self-shadowing and its solar eclipses by Mars. The propagated Phobos' state was compared with data from JPL's Horizon system to ensure the accuracy of the result. Results for Phobos lighting condition over one Martian year are presented in this paper, which include length of solar eclipse, average solar radiation intensity, surface exposure time, total maximum solar energy, and total surface solar energy (constrained by incident angle). The results show that Phobos' solar eclipse time changes throughout the Martian year with the maximum eclipse time occurring during the Martian spring and fall equinox and no solar eclipse during the Martian summer and winter solstice. Solar radiation intensity is close to minimum at the summer solstice and close to maximum at the winter solstice. Total surface exposure time is longer near the north pole and around the anti- Mars point. Total maximum solar energy is larger around the anti-Mars point. Total surface solar energy is higher around the anti-Mars point near the equator. The results from this study and others like it will be important in determining landing site selection, vehicle system design and mission operations for the human exploration of Phobos and subsequently Mars.

Li, Zu Qun↗

A New Cell-Centered Implicit Numerical Scheme for Ions in the 2-D Axisymmetric Code Hall2de

We present a new algorithm in the Hall2De code to simulate the ion hydrodynamics in the acceleration channel and near plume regions of Hall-effect thrusters. This implementation constitutes an upgrade of the capabilities built in the Hall2De code. The equations of mass conservation and momentum for unmagnetized ions are solved using a conservative, finite-volume, cell-centered scheme on a magnetic-field-aligned grid. Major computational savings are achieved by making use of an implicit predictor/multi-corrector algorithm for time evolution. Inaccuracies in the prediction of the motion of low-energy ions in the near plume in hydrodynamics approaches are addressed by implementing a multi-fluid algorithm that tracks ions of different energies separately. A wide range of comparisons with measurements are performed to validate the new ion algorithms. Several numerical experiments with the location and value of the anomalous collision frequency are also presented. Differences in the plasma properties in the near-plume between the single fluid and multi-fluid approaches are discussed. We complete our validation by comparing predicted erosion rates at the channel walls of the thruster with measurements. Erosion rates predicted by the plasma properties obtained from simulations replicate accurately measured rates of erosion within the uncertainty range of the sputtering models employed.

electric propulsion↗

Using Big Data Technologies with Earth Science Data in HDF5: HDF5 Scalable Solutions

HDF5 (Hierarchical Data Format 5) is open-source, high-performance software that consists of an abstract data model, library, and fileformat used for storing and managing extremely large and/or complex data collections. NASA Earth Observing System (EOS) Data and Information Systems use HDF5 as an archival format to store remote sensing data from EOS satellites. HDF5 is also used to store other types of Geoscience and Strophysical data, e.g., seismic data and data from Low-Frequency Array (LOFAR) radio telescopes. Data stored in HDF5 has reached tens of petabytes and is growing at an accelerated rate.With the growing amout of HDF5 Earth Science data to analyze and process, scientists need to adopt big data technologies including new storage paradigms such as cloud and object storage. To run models and perform data analysis they also need to utilizied efficient and diverse ways to access data, from high-performance computing's (HPC) Message Passing Interface (MPI) I/O and deep memory hierarchies (DMH) to non-HPC frameworks such as Apache Hadoop, Spark, and Drill. The HDF Group continually works to enable usage of big data technologies in HDF software.

Knox, Larry↗

The 'convection surge' mechanism of ion acceleration during substorms

To explain the presence of the so-called 'bounce-phase-bunched' ion distributions (eV range to about 20 keV) observed in the earth's geosynchronous magnetosphere in association with substorm expansion phases, Quinn and Southwood proposed the 'convection surge' mechanism of ion energization. This mechanism is associated with a sudden earthward displacement of curved field lines resulting from the short-lived application of an intense, east-west electric field. A quantitative computer model of the mechanism has been constructed here. It is found that the mechanism easily generates the bounce-phase-bunched ion distributions in question. Additionally, the mechanism can generate dramatically field-aligned distributions. It is proposed that a key general characteristic of geosynchronous ion pitch angle distributions, field-aligned low energies, and field-perpendicular high energies is a signature of this energization mechanism. The presence of field-aligned ions cannot be presumed to be a signature of recent ionospheric extraction.

Mauk, B. H.↗

A Dose of Reality: Radiation Analysis for Realistic Human Spacecraft

INTRODUCTION As with most computational analyses, a tradeoff exists between problem complexity, resource availability and response accuracy when modeling radiation transport from the source to a detector. The largest amount of analyst time for setting up an analysis is often spent ensuring that any simplifications made have minimal impact on the results. The vehicle shield geometry of interest is typically simplified from the original CAD design in order to reduce computation time, but this simplification requires the analyst to "re-draw" the geometry with a limited set of volumes in order to accommodate a specific radiation transport software package. The resulting low-fidelity geometry model cannot be shared with or compared to other radiation transport software packages, and the process can be error prone with increased model complexity. The work presented here demonstrates the use of the DAGMC (Direct Accelerated Geometry for Monte Carlo) Toolkit from the University of Wisconsin, to model the impacts of several space radiation sources on a CAD drawing of the US Lab module. METHODS The DAGMC toolkit workflow begins with the export of an existing CAD geometry from the native CAD to the ACIS format. The ACIS format file is then cleaned using SpaceClaim to remove small holes and component overlaps. Metadata is then assigned to the cleaned geometry file using CUBIT/Trelis from csimsoft (Registered Trademark). The DAGMC plugin script removes duplicate shared surfaces, facets the geometry to a specified tolerance, and ensures that the faceted geometry is water tight. This step also writes the material and scoring information to a standard input file format that the analyst can alter as desired prior to running the radiation transport program. The scoring results can be transformed, via python script, into a 3D format that is viewable in a standard graphics program. RESULTS The CAD model of the US Lab module of the International Space Station, inclusive of all the racks and components, was simplified to remove holes and volume overlaps. Problematic features within the drawing were also removed or repaired to prevent runtime issues. The cleaned drawing was then run through the DAGMC workflow to prepare for analysis. Pilot tests modeling transport of 1GeV proton and 800MeV/A oxygen sources show that reasonable results are converged upon in an acceptable amount of overall computation time from drawing preparation to data analysis. The FLUKA radiation transport code will next be used to model both a GCR and a trapped radiation source. These results will then be compared with measurements that have been made by the radiation instrumentation deployed inside the US Lab module. DISCUSSION Early analyses have indicated that the DAGMC workflow is a promising toolkit for running vehicle geometries of interest to NASA through multiple radiation transport codes. In addition, recent work has shown that a realistic human phantom, provided via a subcontract with the University of Florida, can be placed inside any vehicle geometry for a combinatorial analysis. This added functionality gives the user the ability to score various parameters at the organ level, and the results can then be used as input for cancer risk models.

Barzilla, J. E.↗

Ride comfort control in large flexible aircraft

The problem of ameliorating the discomfort of passengers on a large air transport subject to flight disturbances is examined. The longitudinal dynamics of the aircraft, including effects of body flexing, are developed in terms of linear, constant coefficient differential equations in state variables. A cost functional, penalizing the rigid body displacements and flexure accelerations over the surface of the aircraft is formulated as a quadratic form. The resulting control problem, to minimize the cost subject to the state equation constraints, is of a class whose solutions are well known. The feedback gains for the optimal controller are calculated digitally, and the resulting autopilot is simulated on an analog computer and its performance evaluated.

Warren, M. E.↗