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,081 records · Page 60

Mapping implicit spectral methods to distributed memory architectures

Spectral methods were proven invaluable in numerical simulation of PDEs (Partial Differential Equations), but the frequent global communication required raises a fundamental barrier to their use on highly parallel architectures. To explore this issue, a 3-D implicit spectral method was implemented on an Intel hypercube. Utilization of about 50 percent was achieved on a 32 node iPSC/860 hypercube, for a 64 x 64 x 64 Fourier-spectral grid; finer grids yield higher utilizations. Chebyshev-spectral grids are more problematic, since plane-relaxation based multigrid is required. However, by using a semicoarsening multigrid algorithm, and by relaxing all multigrid levels concurrently, relatively high utilizations were also achieved in this harder case.

Overman, Andrea L.↗

Adsorbent testing and mathematical modeling of a solid amine regenerative CO2 and H2O removal system

The paper examines the design and the construction details of the test bed built for testing a solid-amine-based Regenerable CO2 Removal System (RCRS) built at the NASA/Johnson Space Center for the extended Orbiter missions. The results of tests are presented, including those for the adsorption breakthrough and the adsorption and desorption of CO2 and H2O vapor. A model for predicting the performance of regenerative CO2 and H2O vapor adsorption of the solid amine system under various operating conditions was developed in parallel with the testing of the test stand, using the coefficient of mass transfer calculated from test results. The results of simulations are shown to predict the adsorption performance of the Extended Duration Orbiter test bed fairly well. For the application to the RCRS at various operating conditions the model has to be modified.

Jeng, F. F.↗

Nonlinear evolution of a large-amplitude circularly polarized Alfven wave: High beta

The nonlinear dynamics following saturation of the parametric instabilities of a monochromatic field-aligned large-amplitude circularly polarized Alfven wave is investigated via direct numerical simulation in the case of high plasma beta and no wave dispersion. The magnetohydrodynamic (MHD) code permits nonlinear couplings in the parallel direction to the ambient magnetic field and one perpendicular direction. Compressibility is included in the form of a polytropic equation of state. Turbulent cascades develop after saturation of two coupled oblique three-wave parametric instabilities; one of which is an oblique filamentationlike instability reported earlier. Remnants of the parametric processes, as well as of the original Alfven pump wave, persist during late nonlinear times. Nearly incompressible MHD features such as spectral anisotropies appear as well.

Ghosh, S.↗

Effect of B(sub y) on neutral line ridges and dynamical source ordering

In this paper we study the effect of a uniform magnetic field B(sub y) in the cross-tail y direction on the dynamics and distribution of energetic ions in a current sheet model including a neutral line. Martin and Speiser (1988) have previously shown that a 'ridge' in the velocity space distribution function is a remote neutral line signature in a two-dimensional field without B(sub y). Our results show that as B(sub y) is increased, using nominal tail parameters, there is very little change in the ridge signature for small B(sub y) (up to about 4 times B(sub z)). For intermediate values (up to about 10 times times B(sub Y)) the ridge becomes observable further from the neutral line, while close to the X line the ridge is strongly modified. For large B(sub y) (of the order of B(sub x)) the ridge evolves into the new neutral line signatures which depend strongly on whether the observations point is above or below the current sheet. These results are used to estimate B(sub y) less than 2 nT in the Active Magentospheric Particle Traces Explorers/Ion Release Module (AMPTE/IRM) event modeled by Speiser and Martin (1994), which is consistent with onboard magnetometer measurements. We further show that the cause of the new structures appearing at large B(sub Y) is the same source seperation effect that produced the original ridge and that all these signatures are relatively insensitive to a small parallel electric field. Finally, we point our a dynamical ordering which occurs at large B(sub y): Using high-resolution simulations, we find a common boundary separating regions of positive and negative initial x, y, and z position, as well as initial pitch angle. This boundary is purely dynamical in nature and is independent of modeled source distributions.

Martin, R. F., Jr.↗

Anisotropic turbulence in the solar wind

Solar wind turbulence has been viewed traditionally as composed of parallel propagating ('slab' fluctuations) or otherwise as isotropic turbulence. A variety of recent investigations, reviewed here, indicate that the spectrum may contain a significant admixture of two dimensional fluctuations, having variations mainly perpendicular to the local magnetic field. These indications come from simulations, from the theory of nearly incompressible MHD, from cosmic ray transport studies and from transport theory for solar wind turbulence, as well as from interpretations of direct observations. Thus, solar wind turbulence may be more like bundles of spaghetti than like parallel phase fronts.

Matthaeus, W. H.↗

Computation of High Speed Jet Noise

The objective of this work was to use direct numerical simulation (DNS) techniques to study the physics of noise generation by a high-speed turbulent jet. A Mach 0.9, Reynolds number 3,600 jet was selected because of available experimental data. New numerical methods for generating disturbances at the nozzle and computing far-field sound were developed and reported in the course of this work. Over 25 million mesh points were used in the simulations which ran for over 50,000 timesteps and required over 50,000 processor hours on state-of-the-art parallel computer systems to complete. Figures show a visualization of the jet and sound field, a comparison of the mean flow development with the experiment, a directivity comparison with the experiment, and time spectrum comparison with the experiment. Agreement is seen to be excellent. These are fully document in the attached references. Full details of the work, detailed achievements and conclusions are discussed in appendices, which are copies of publications that resulted from this work. We have studied noise mechanisms in supersonic jets, the refraction of sound by turbulence in subsonic jets, and noise sources in conjunction with a DNS of a Mach 0.9 jet.

Freund, Jonathan B.↗

Parametric Powered-Lift Navier-Stokes Computations

The goal of this work is to enable the computation of large numbers of unsteady high-fidelity flow simulations for a YAV-8B Harrier aircraft in ground effect by improving the solution process and taking advantage of NASA parallel supercomputers. The YAV-8B Harrier aircraft can take off and land vertically, or utilize short runways by directing its four exhaust nozzles toward the ground. Transition to forward flight is achieved by rotating these nozzles into a horizontal position.

Chaderjian, Neal M.↗

Experimental Study of Arcing on High-voltage Solar Arrays

The main obstacle to the implementation of a high-voltage solar array in space is arcing on the conductor-dielectric junctions exposed to the surrounding plasma. One obvious solution to this problem would be the installation of fully encapsulated solar arrays which were not having exposed conductors at all. However, there are many technological difficulties that must be overcome before the employment of fully encapsulated arrays will turn into reality. An alternative solution to raise arc threshold by modifications of conventionally designed solar arrays looks more appealing, at least in the nearest future. A comprehensive study of arc inception mechanism [1-4] suggests that such modifications can be done in the following directions: i) to insulate conductor-dielectric junction from a plasma environment (wrapthrough interconnects); ii) to change a coverglass geometry (overhang); iii) to increase a coverglass thickness; iiii) to outgas areas of conductor-dielectric junctions. The operation of high-voltage array in LEO produces also the parasitic current power drain on the electrical system. Moreover, the current collected from space plasma by solar arrays determines the spacecraft floating potential that is very important for the design of spacecraft and its scientific apparatus. In order to verify the validity of suggested modifications and to measure current collection five different solar array samples have been tested in large vacuum chamber. Each sample (36 silicon based cells) consists of three strings containing 12 cells connected in series. Thus, arc rate and current collection can be measured on every string independently, or on a whole sample when strings are connected in parallel. The heater installed in the chamber provides the possibility to test samples under temperature as high as 80 C that simulates the LEO operational temperature. The experimental setup is described below.

Vayner, Boris↗

An Assessment of GEO Orbital Debris Photometric Properties Derived from Laboratory-Based Measurements

Optical observations of orbital debris offer insights that differ from radar measurements (specifically the size parameter and wavelength regime). For example, time-dependent photometric data yield lightcurves in multiple bandpasses that aid in material identification and possible periodic orientations. This data can also be used to help identify shapes and optical properties at multiple phase angles. Capitalizing on optical data products and applying them to generate a more complete understanding of orbital space objects, is a key objective of NASA s Optical Measurement Program, and a primary driver for creation of the Optical Measurements Center (OMC). The OMC attempts to emulate space-based illumination conditions using equipment and techniques that parallel telescopic observations and source-target-sensor orientations. The OMC uses a 300 Watt Xenon arc lamp as a solar simulator, a CCD camera with Johnson/Bessel colored filters, and a robotic arm to orientate/rotate objects to simulate an object's orbit/rotational period. A high-resolution, high bandwidth (350nm-2500nm) Analytical Spectral Devices (ASD) spectrometer is also employed to baseline various material types. Since observation of GEO targets are generally restricted to the optical regime (due to radar range limitations), analysis of their properties is tailored to those revealed by optical data products. In this connection, much attention has been directed towards understanding the lightcurves of orbital debris with high area-to-mass (A/m) ratios (greater than 0.9 square meters per kilogram). A small population of GEO debris was recently identified, which exhibits the properties of high A/m objects, such as variable eccentricities and inclinations a dynamical characteristic generally resulting from varying solar radiation pressure on high A/m objects. Materials such as multi-layered insulation (MLI) and solar panels are two examples of materials with high area-to mass ratios. Lightcurves for such objects can vary greatly (even for the same object under different illumination conditions). For example, specular reflections from multiple facets of the target surface (e.g. Mylar or Aluminized Kapton) can lead to erratic, orientation-dependent light curves. This paper will investigate published color photometric data for a series of orbital debris targets and compare it to the empirical photometric measurements generated in the OMC. The specific materials investigated (known to exist in GEO) are: an intact piece of MLI, separated layers of MLI, and multiple solar cells materials. Using the data acquired over specific rotational angles through different filters (B, V, R, I), a color index is acquired (B-R, R-I). As a secondary check, the spectrometer is used to define color indexes for the same material. Using these values and their associated lightcurves, this laboratory data is compared to observational data obtained on the 1m telescope of the Astronomical Institute of the University of Bern (AUIB) and the 0.9 m Small and Moderate Aperture Research Telescope System (SMARTS) telescope at Cerro Tololo Inter-American Observatory (CTIO). We will present laboratory generated lightcurves with color indexes of the high A/m materials alongside telescopic data of targets with high A/m values. We will discuss the relationship of laboratory to telescope data in the context of classification of GEO debris objects.

Rodriquez-Cowardin, H.↗

Characterization of Orbital Debris Photometric Properties Derived from Laboratory-Based Measurements

To better characterize and model optical data acquired from ground-based telescopes, the Optical Measurements Center (OMC) at NASA/JSC attempts to emulate illumination conditions seen in space using equipment and techniques that parallel telescopic observations and source-target-sensor orientations. The OMC uses a 75 Watt Xenon arc lamp as a solar simulator, an SBIG CCD camera with standard Johnson/Bessel filters, and a robotic arm to simulate an object's position and rotation. The laboratory uses known shapes, materials suspected to be consistent with the orbital debris population, and three phase angles to best match the lighting conditions of the telescope based data. The fourteen objects studied in the laboratory are fragments or materials acquired through ground-tests of scaled-model satellites/rocket bodies as well as material samples in more/less "flight-ready" condition. All fragments were measured at 10 increments in a full 360 rotation at 6 , 36 , and 60 phase angles. This paper will investigate published color photometric data for a series of orbital debris targets and compare it to the empirical photometric measurements generated in the OMC. Using the data acquired over specific rotational angles through different filters (B, V, R, I), a color index is acquired (B-R, R-I). Using these values and their associated lightcurves, this laboratory data is compared to observational data obtained on the 1 m telescope of the Astronomical Institute of the University of Bern (AUIB), the 0.9 m operated by the Small- and Medium-Aperture Research Telescope System (SMARTS) Consortium and the Curtis-Schmidt 0.6 m Michigan Orbital Debris Space Debris Telescope both located at Cerro Tololo Inter-American Observatory (CTIO). An empirical based optical characterization model will be presented to provide preliminary correlations between laboratory based and telescope-based data in the context of classification of GEO debris objects.

Cowardin, H.↗

Professor Created On-line Biology Laboratory Course

This paper will share the creation, implementation, and modification of an online college level general biology laboratory course offered for non-science majors as a part of a General Education Curriculum. The ability of professors to develop quality online laboratories will address a growing need in Higher Education as more institutions combine course sections and look for suitable alternative course delivery formats due to declining departmental budgets requiring reductions in staffing, equipment, and supplies. Also, there is an equal or greater need for more professors to develop the ability to create online laboratory experiences because many of the currently available online laboratory course packages from publishers do not always adequately parallel on-campus laboratory courses, or are not as aligned with the companion lecture sections. From a variety of scientific simulation and animation web sites, professors can easily identify material that closely fit the specific needs of their courses, instructional environment, and students that they serve. All too often, on-campus laboratory courses in the sciences provide what are termed confirmation experiences that do NOT allow students to experience science as would be carried out by scientists. Creatively developed online laboratory experiences can often provide the type of authentic investigative experiences that are not possible on-campus due to the time constraints of a typical two-hour, once-per-week-meeting laboratory course. In addition, online laboratory courses can address issues related to the need for students to more easily complete missing laboratory assignments, and to have opportunities to extend introductory exercises into more advanced undertakings where a greater sense of scientific discovery can be experienced. Professors are strongly encourages to begin creating online laboratory exercises for their courses, and to consider issues regarding assessment, copyrights, and Intellectual Property concerns.

Bowman, Arthur W.↗

Optical Signature Analysis of Tumbling Rocket Bodies via Laboratory Measurements

The NASA Orbital Debris Program Office has acquired telescopic lightcurve data on massive intact objects, specifically spent rocket bodies, in order to ascertain tumble rates in support of the Active Debris Removal (ADR) task to help remediate the LEO environment. Rotation rates are needed to plan and develop proximity operations for potential future ADR operations. To better characterize and model optical data acquired from ground-based telescopes, the Optical Measurements Center (OMC) at NASA/JSC emulates illumination conditions in space using equipment and techniques that parallel telescopic observations and source-target-sensor orientations. The OMC employs a 75-watt Xenon arc lamp as a solar simulator, an SBIG CCD camera with standard Johnson/Bessel filters, and a robotic arm to simulate an object's position and rotation. The light source is mounted on a rotary arm, allowing access any phase angle between 0 -- 360 degrees. The OMC does not attempt to replicate the rotation rates, but focuses on how an object is rotating as seen from multiple phase angles. The two targets studied are scaled (1:48), SL-8 Cosmos 3M second stages. The first target is painted in the standard government "gray" scheme and the second target is primary white, as used for commercial missions. This paper summarizes results of the two scaled rocket bodies, each rotated about two primary axes: (a) a spin-stabilized rotation and (b) an end-over-end rotation. The two rotation states are being investigated as a basis for possible spin states of rocket bodies, beginning with simple spin states about the two primary axes. The data will be used to create a database of potential spin states for future works to convolve with more complex spin states. The optical signatures will be presented for specific phase angles for each rocket body and shown in conjunction with acquired optical data from multiple telescope sources.

Cowardin, H.↗

Distributed Prognostics based on Structural Model Decomposition

Within systems health management, prognostics focuses on predicting the remaining useful life of a system. In the model-based prognostics paradigm, physics-based models are constructed that describe the operation of a system and how it fails. Such approaches consist of an estimation phase, in which the health state of the system is first identified, and a prediction phase, in which the health state is projected forward in time to determine the end of life. Centralized solutions to these problems are often computationally expensive, do not scale well as the size of the system grows, and introduce a single point of failure. In this paper, we propose a novel distributed model-based prognostics scheme that formally describes how to decompose both the estimation and prediction problems into independent local subproblems whose solutions may be easily composed into a global solution. The decomposition of the prognostics problem is achieved through structural decomposition of the underlying models. The decomposition algorithm creates from the global system model a set of local submodels suitable for prognostics. Independent local estimation and prediction problems are formed based on these local submodels, resulting in a scalable distributed prognostics approach that allows the local subproblems to be solved in parallel, thus offering increases in computational efficiency. Using a centrifugal pump as a case study, we perform a number of simulation-based experiments to demonstrate the distributed approach, compare the performance with a centralized approach, and establish its scalability. Index Terms-model-based prognostics, distributed prognostics, structural model decomposition ABBREVIATIONS

centrifugal pump↗

Designing Liquid Rocket Engine Injectors for Performance, Stability, and Cost

NASA is developing the Space Launch System (SLS) for crewed exploration missions beyond low Earth orbit. Marshall Space Flight Center (MSFC) is designing rocket engines for the SLS Advanced Booster (AB) concepts being developed to replace the Shuttle-derived solid rocket boosters. One AB concept uses large, Rocket-Propellant (RP)-fueled engines that pose significant design challenges. The injectors for these engines require high performance and stable operation while still meeting aggressive cost reduction goals for access to space. Historically, combustion stability problems have been a critical issue for such injector designs. Traditional, empirical injector design tools and methodologies, however, lack the ability to reliably predict complex injector dynamics that often lead to combustion stability. Reliance on these tools alone would likely result in an unaffordable test-fail-fix cycle for injector development. Recently at MSFC, a massively parallel computational fluid dynamics (CFD) program was successfully applied in the SLS AB injector design process. High-fidelity reacting flow simulations were conducted for both single-element and seven-element representations of the full-scale injector. Data from the CFD simulations was then used to significantly augment and improve the empirical design tools, resulting in a high-performance, stable injector design.

Westra, Douglas G.↗

Automated Cooperative Trajectories

This presentation is an overview of the Automated Cooperative Trajectories project. An introduction to the phenomena of wake vortices is given, along with a summary of past research into the possibility of extracting energy from the wake by flying close parallel trajectories. Challenges and barriers to adoption of civilian automatic wake surfing technology are identified. A hardware-in-the-loop simulation is described that will support future research. Finally, a roadmap for future research and technology transition is proposed.

fuel comsumption↗

A Component-Level Model of Automatic Dependent Surveillance - Broadcast (ADS-B)

Automatic Dependent Surveillance – Broadcast (ADS-B) is being employed in numerous peer-to-peer initiatives attempting to expand the capacity of the National Airspace System (NAS) or enable mixed operations of manned and unmanned vehicles. Safety assessments of these initiatives rely, in part, on modeling the accuracy of ADS-B in reporting the position and direction of an ownship and surrounding traffic. Frequently, these initiatives utilize a position uncertainty model that applies a reported ADS-B estimation position uncertainty (EPU) value to a Rayleigh distribution and uses a Gauss-Markov random walk to add error to the ADS-B output of a vehicle. This model of ADS-B state error is easy to implement and apply to numerous problems. However, it has a couple of draw-backs. First, the ADS-B state errors are equally probable in all directions. This is a good assumption in situations where aircraft maneuvering is not constrained. However, in situations where the aircraft maneuvering is constrained such as landing, the error distribution is likely to exhibit directionality and the non-directional model may skew results especially when assessing very low probabilities (e.g., 10(exp -9)) of catastrophic encounters. Second, the model does not account for processing latency in the receiving aircraft. NASA Langley Research Center (LaRC) recently examined the feasibility of decreasing the spacing of aircraft on parallel approaches to runways separated by as little as 700 feet. For Monte-Carlo analysis using a high-fidelity simulation of a large transport, LaRC started with a Gauss-Markov model of ADS-B error but then developed a component level model of ADS-B error to increase the fidelity of results.

Madden, Michael M.↗

Wall-modeled LES of the NASA Juncture Flow Experiment

We perform wall-modeled large eddy simulation (WMLES) of the NASA Juncture Flow experiment (Kegerise & Neuhart, NASA/TM–2019–20286) using an equilibrium wall model and unstructured finite volume solver, to assess its predictive ability for complex geometries. The flow condition simulated corresponds to 5 degrees angle of attack with a freestream Mach number of 0.189, and a Reynolds number based on the mean crank chord of 2.4 million. All simulations neglect the effects of wind tunnel walls, sting and mast present in the experiment. To assess the performance of WMLES on simpler flows with very coarse grids, we simulate flow in a turbulent channel at a friction Reynolds number, Reτ≈2000, and flow through a square duct at anReτ≈1000. The duct flow results using the equilibrium wall model indicate that the skin friction is inaccurate near the corner for grid topologies and resolutions typically used in WMLES, having implications for more complicated flows such as the juncture flow. For the juncture flow simulation, two different approaches are investigated. The first approach consists of a truncated-domain simulation wherein the inflow plane of the computational domain is placed at≈0.2chord length of the wing, and we prescribe the mean flow from a separate Reynolds-averaged Navier-Stokes (RANS) solution along with synthetic turbulence to initiate realistic unsteadiness in the domain. The second approach involves simulating the entire geometry with trips to trigger transition to turbulence. The truncated- and full-domain simulations contain about 62 and 90 million cells, respectively; with 8-10 points per boundary layer thickness in the wall-parallel directions, 14-20 points in the wall-normal direction and a near-wall viscous spacing (∆n+1) of≈100. Preliminary results are encouraging overall in terms of the prediction of wall pressure, wall skin friction, velocity and stresses; but indicate that further work is required to improve the predictions in the separation bubble and wing-fuselage corner regions. The grid resolution used in this study is still fairly coarse, and the results should be interpreted as work in progress.

Prahladh S Iyer↗

Automatic Incorporation of Flaws into Cad and FEA Models Based On Measured NDE Data

Automatic incorporation of non-destructive evaluation (NDE) data into CAD and finite element (FE) models for structural analysis is essential for accurate residual-life and strength analyses of composite structures. Currently, two significant but independent challenges exist for incorporating NDE data in such analyses: interpretation of the NDE signals to locate and characterize the structural defects and the incorporation of that information into structural analysis tools. The first challenge will be addressed through linear inversion of calculated synthetic responses generated by a simulation engine to evaluate flaw size, location, and other characteristics. The second challenge will be addressed by building on prior and parallel defect modeling efforts to integrate suitable flaw models into the CAD/finite-element representation of a specimen. Ultrasonic NDE information is usually represented using C-Scans of time-of-flight or intensity information. Other modalities such as thermography and computed tomography (CT) could also be used to characterize the structure. This information will first be projected into material coordinates and then aligned with the specimen’s finite element mesh. The NDE data can then be incorporated by inserting flaw geometry into the CAD and FE analysis model and automatic modification of local material properties and boundary conditions. This process will enable the use of structural analysis tools even for practical modalities such as ultrasound and thermography which may not directly give accurate-enough geometric information to fully define the flaw. The full flaw characteristics will be inferred using prior expert and modeling knowledge to categorize flaw signatures. Outputs for this project will include demonstrated methods for creating structural analysis geometric models that incorporate flaw characterization information; these models can then be used directly in popular commercial FEA software. In addition, the primary deliverable for the project will include a software that will be able to automatically incorporate different kinds of flaws into an existing FE models. The models will be generated in standard file formats to enable direct import into the analysis software.

Adarsh Krishnamurthy↗