Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “distributed parallelization”

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,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Microwave Kinetic Inductance Detector with Selective Polarization Coupling

A conventional low-noise detector requires a technique to both absorb incident power and convert it to an electrical signal at cryogenic temperatures. This innovation combines low-noise detector and readout functionality into one device while maintaining high absorption, controlled polarization sensitivity, and broadband detection capability. The resulting far-infrared detectors can be read out with a simple approach, which is compact and minimizes thermal loading. The proposed microwave kinetic inductance detector (MKID) consists of three basic elements. The first is the absorptive section in which the incident power is coupled to a superconducting resonator at far-infrared frequency above its superconducting critical frequency (where superconductor becomes normal conductor). This absorber's shape effectively absorbs signals in the desired polarization state and is resonant at the radio frequency (RF) used for readout of the device. Control over the metal film used in the absorber allows realization of structures with either a 50% broadband or 100% resonance absorptance over a 30% fractional bandwidth. The second element is a microwave resonator - which is realized from the thin metal films used to make the absorber as transmission lines - whose resonance frequency changes due to a variation in its kinetic inductance. The resonator's kinetic inductance is a function of the power absorbed by the device. A low-loss dielectric (mono-crystalline silicon) is used in a parallel-plate transmission line structure to realize the desired superconducting resonators. There is negligible coupling among the adjacent elements used to define the polarization sensitivity of each detector. The final component of the device is a microwave transmission line, which is coupled to the resonator, and allows detection of changes in resonance frequency for each detector in the focal plane array. The spiral shape of the detector's absorber allows incident power with two polarizations to couple to the detector equally. A stepped impedance resonator was used that allows the incident power absorbed in the detecting membrane area to be uniformly distributed in the detector's transmission line at the RF readout frequency. This maximizes the sensitivity of the detector. The signal is read out via a frequency multiplexing technique that requires a minimum number of interface transmission lines for readout. This reduces the packaging complexity and coupling to the device's thermal environment.

Wollack, Edward↗
NASA NTRS2018

A Novel Multi-Spacecraft Interplanetary Global Trajectory Optimization Transcription

As the frontier of space exploration continues to advance, so does the design complexity of future interplanetary missions. One avenue of this increasing complexity includes a class of designs known as "Distributed Spacecraft Missions"; missions where multiple spacecraft coordinate to perform shared objectives. Current approaches for the global trajectory optimization of these Multi-Vehicle Missions (MVMs) are prone to shortcomings including laborious iterative design, considerable human-in-the-loop effort, treatment of the multi-vehicle problem as multiple separate trajectory optimization subproblems (resulting in suboptimal solutions where the whole is less than the sum of its parts), and poor handling of coordination objectives and constraints. There are only a handful of software platforms in existence capable of fully-automated, rapid, interplanetary mission and systems global optimization including the Parallel Global Multiobjective Optimizer (PaGMO), the Gravity Assisted Low-thrust Local Optimization Program (GALLOP), and the Evolutionary Mission Trajectory Generator (EMTG). However, none of these tools is capable of performing such tasks for MVM designs. The work outlined in this paper lays the groundwork for a technique to begin addressing these shortcomings. We present a fully-automated technique which frames interplanetary MVMs as Multi-Objective, Multi-Agent Hybrid Optimal Control Problems (MOMA HOCP). First, the basic functionality of this technique is validated on the single-vehicle problem of reproducing the Cassini interplanetary cruise.

Interplanetary↗
NASA NTRS2018

A Novel Multi-Spacecraft Interplanetary Global Trajectory Optimization Transcription

As the frontier of space exploration continues to advance, so does the design complexity of future interplanetary missions. One avenue of this increasing complexity includes a class of designs known as ``Distributed Spacecraft Missions"; missions where multiple spacecraft coordinate to perform shared objectives. Current approaches for the global trajectory optimization of these Multi-Vehicle Missions (MVMs) are prone to shortcomings including laborious iterative design, considerable human-in-the-loop effort, treatment of the multi-vehicle problem as multiple separate trajectory optimization subproblems (resulting in suboptimal solutions where the whole is less than the sum of its parts), and poor handling of coordination objectives and constraints. There are only a handful of software platforms in existence capable of fully-automated, rapid, interplanetary mission and systems global optimization including the Parallel Global Multiobjective Optimizer (PaGMO), the Gravity Assisted Low-thrust Local Optimization Program (GALLOP), and the Evolutionary Mission Trajectory Generator (EMTG). However, none of these tools is capable of performing such tasks for MVM designs. The work outlined in this paper lays the groundwork for a technique to begin addressing these shortcomings. We present a fully-automated technique which frames interplanetary MVMs as Multi-Objective, Multi-Agent Hybrid Optimal Control Problems (MOMA HOCP). First, the basic functionality of this technique is validated on the single-vehicle problem of reproducing the Cassini interplanetary cruise.

Napier, Sean W.↗
NASA NTRSDate not supplied

NASA Langley FUN3D Analyses in Support of the 1st AIAA Stability and Control Prediction Workshop

This work summarizes the results of FUN3D analyses conducted for the 1st AIAA Stability and Control Workshop on behalf of participants from the NASA Langley Research Center. The workshop was created to establish best practices for the prediction of stability and control derivatives using computational fluid dynamics and assess the limitations of these methods when those best practices are applied. The inaugural workshop considered the ONERA version of the NASA/Boeing Common Research Model, which includes the wing, body, horizontal tail, and a vertical tail designed by ONERA. Wind tunnel data at small sideslip angles remain unpublished and served as ‘blind’ data for computational comparisons. The present research generated workshop test case data using the NASA FUN3D code, which is a parallelized, unstructured, node-based, finite-volume discretization, Reynolds-averaged Navier-Stokes flow solver. Steady- state numerical simulations were conducted for workshop test cases investigating the following: grid convergence, Mach number effect on static stability, wind tunnel sting increments, static stability-derivative calculations, and a sideslip angle sweep. Results were generated for two series of unstructured, mixed-element grids, one set provided by the workshop and another set created using the HeldenMesh grid generation software. The results provided include total- and component-level breakdowns of the force and moment coefficients, in addition to sectional pressure distributions for the wing and tail components for comparisons to wind tunnel data.

CFD↗
NASA NTRS1997

Mid-Infrared Observational and Theoretical Studies of Star Formation and Early Solar Systems

The first 2 years of this program were used to make mid-IR observations of regions of star formation in the Orion nebula with the UCSD mid-IR camera at the UCSD/University of Minnesota telescope at Mt. Lemmon. These observations attempted to make the first systematic study of an extended region, known to have newly forming stars, and expected to have complex mid-IR emission. We discovered, to our surprise, that most of the thermal emission originated from extended sources rather than from point sources. This interesting observation made the analysis of the data much more complex, since the chop/nod procedures used at these wavelengths produce a differential measurement of the emission in one region compared to that in the adjacent region. Disentangling complex extended emission in such a situation is very difficult. In parallel with this work we were also observing comets in the thermal infrared, the other component of the original proposal. Some spectacular data on the comet Swift-Tuttle was acquired and published. A changing jet structure observed over a 2 week period is described. The rotation period of the comet can be measured at 66 hours. The size of the nucleus can also be estimated (at 30 km) from the observed excess flux from the nucleus. These data have lead to the development of models describing the action of dust particles of differing sizes and composition leaving the nucleus. The spatial distribution of the predicted IR emission has been compared to the observed jet structures, leading to estimates of both particles sizes, relative amounts of silicate vs organic grains, and the amounts of dust emitted in the jets vs isotopic emission.

Jones, Barbara↗
NASA NTRS1975

Experimental and analytical investigation of 0 G condensation in a mechanical refrigeration system application

Basic equations of momentum and energy are presented and discussed with respect to heat transfer and pressure drop for forced flow condensation in horizontal tubes under 1-g and 0-g conditions. Some experimental results are presented for condensing refrigerant-12 in a system of three parallel-connected quartz tubes (3-mm inside diameter, G = 1.037 to 3.456 x 105 lbm/hr-sq. ft). From high speed photographs, measurements were obtained of film thickness, phase velocities, disturbance wavelengths, and flow regimes and their transitions. Based upon these measurements various dimensionless force ratios (flow and instability parameters) were calculated. Under 0-g conditions a uniformly thick redistribution of liquid condensate about the tube walls was found to result in a lowered heat transfer coefficient as compared with 1-g conditions, based upon fundamental heat transfer theory. A model is proposed that takes into account the difference in heat transfer due to condensate distribution under 1-g and 0-g conditions.

Keshock, E. G.↗
NASA NTRS1984

Severe storm initiation and development from satellite infrared imagery and Rawinsonde data

The geographical distribution of potential temperatures, mixing ratio, and streamlines of flow patterns at 850, 700, and 500 mb heights are used to understand the prestorm convection and the horizontal convergence of moisture. From the analysis of 21 tornadoes the following conclusions are reached: (1) Strong horizontal convergence of moisture appeared at the 850, 700, and 500 mb levels in the area 12 hours before the storm formation; (2) An abundantly moist atmosphere below 3 km (700 mb) becomes convectively unstable during the time period between 12 and 24 hours before the initiation of the severe storms; (3) Strong winds veering with height with direction parallel to the movement of a dryline, surface fronts, etc; (4) During a 36-hour period, a tropopause height in the areas of interest is lowest at the time of tornadic cloud formation; (5) A train of gravity waves is detected before and during the cloud formation period. Rapid-scan infrared imagery provides near real-time information on the life cycle of the storm which can be summarized as follows: (1) Enhanced convection produced an overshooting cloud top penetrating above the tropopause, making the mass density of the overshooting cloud much greater than the mass density of the surrounding air; (2) The overshooting cloud top collapsed at the end of the mature stage of the cloud development; (3) The tornado touchdown followed the collapse of the overshooting cloud top.

Hung, R. J.↗
NASA NTRS2018

Evaluation of X-Ray Reflectors by Optical Diffraction Patterns

Performance of X-ray reflectors affects that of X-ray mirrors. Modern X-ray mirrors have thousands of reflectors to gain large effective area. Evaluation of the reflectors is an important process in production of the mirrors. A diffraction pattern dominates reflector image when the parallel optical beam illuminates the reflector along its optical axis because the reflectors are used at grazing incident angles of around 1 deg and their effective width are 1–10 mm. A diffraction pattern from the entire reflector surface can be acquired at once with the aid of a lens. The diffraction pattern holds information of the surface profiles of the reflectors. To quantitatively evaluate the reflectors with the diffraction pattern, we created a diffraction pattern model by the wave optics with the ideal surface profile and fitted it to data. As a result, a correlation between fitting residual and the normal vector distribution of the surface profile was found. With our method, the reflectors can be evaluated and sorted out more efficiently.

Hayashi, Takayuki↗
NASA NTRS2000

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.

Tryggvason, G.↗
NASA NTRS1982

On the structure of the Io Torus

Since the plasma flow near Io is reduced, neutrals originating in charge exchange are not energetic enough to leave the Jovian system and are therefore distributed over an extensive region, as indicated by the sodium cloud. New ions subsequently created in the distributed neutral atomic cloud as a result of charge exchange or electron impact ionization are picked up by the corotating magnetic field, and the radial current driven by the pickup process cannot close in the Io torus, so that it must instead be connected to the planetary ionosphere by field-aligned currents. These field-aligned currents will flow away from the equator at the outer edge of the neutral cloud, and towards it at the inner edge. It is found that the Jovian ionospheric photoelectrons also cannot supply the current flowing away from the equator, so that torus ions accelerated by a parallel electric field could be involved. The parallel potential drop is large enough to push the torus into the Jovian atmosphere, explaining both the sharp, discontinuous change of flux tube content and ion temperature at L equals 5.6, as well as the generation of Auroral-type hiss at that point.

Goertz, C. K.↗
NASA NTRS1985

Heating of thermal helium in the equatorial magnetosphere - A simulation study

The heating of He(+) ions observed by the GEOS satellite at L of about seven starting from an initial anisotropic proton distribution is studied. A simulation of electromagnetic ion cyclotron waves in a plasma consisting of hot anisotropic protons, a dominant thermal hydrogen plasma, and a minority species of helium ions is utilized. It is shown that the thermal ions are heated to 100 times their initial temperature due to the growth of large-amplitude ion cyclotron waves. He(+) ions are heated mainly in the perpendicular direction, and they are heated more than cold H(+) ions. The He(+) ions are heated in a two-step process. First, they are set into oscillations by the growing wave until some of the ions reach a parallel velocity on the order of the resonant velocity; then, strong heating occurs.

Omura, Y.↗
NASA NTRS1987

Particle simulation of plasmas on the massively parallel processor

Particle simulations, in which collective phenomena in plasmas are studied by following the self consistent motions of many discrete particles, involve several highly repetitive sets of calculations that are readily adaptable to SIMD parallel processing. A fully electromagnetic, relativistic plasma simulation for the massively parallel processor is described. The particle motions are followed in 2 1/2 dimensions on a 128 x 128 grid, with periodic boundary conditions. The two dimensional simulation space is mapped directly onto the processor network; a Fast Fourier Transform is used to solve the field equations. Particle data are stored according to an Eulerian scheme, i.e., the information associated with each particle is moved from one local memory to another as the particle moves across the spatial grid. The method is applied to the study of the nonlinear development of the whistler instability in a magnetospheric plasma model, with an anisotropic electron temperature. The wave distribution function is included as a new diagnostic to allow simulation results to be compared with satellite observations.

Gledhill, I. M. A.↗
NASA NTRS1992

Modeling, simulation, and high-autonomy control of a Martian oxygen production plant

Progress on a project for the development of a high-autonomy intelligent command and control architecture for process plants used to produce oxygen from local planetary resources is reported. A distributed command and control architecture is being developed and implemented so that an oxygen production plant, or other equipment, can be reliably commanded and controlled over an extended time period in a high-autonomy mode with high-level task-oriented teleoperation from one or several remote locations. During the reporting period, progress was made at all levels of the architecture. At the remote site, several remote observers can now participate in monitoring the plant. At the local site, a command and control center was introduced for increased flexibility, reliability, and robustness. The local control architecture was enhanced to control multiple tubes in parallel, and was refined for increased robustness. The simulation model was enhanced to full dynamics descriptions.

Schooley, L. C.↗
NASA NTRS1993

Solar flare model atmospheres

Solar flare model atmospheres computed under the assumption of energetic equilibrium in the chromosphere are presented. The models use a static, one-dimensional plane parallel geometry and are designed within a physically self-consistent coronal loop. Assumed flare heating mechanisms include collisions from a flux of non-thermal electrons and x-ray heating of the chromosphere by the corona. The heating by energetic electrons accounts explicitly for variations of the ionized fraction with depth in the atmosphere. X-ray heating of the chromosphere by the corona incorporates a flare loop geometry by approximating distant portions of the loop with a series of point sources, while treating the loop leg closest to the chromospheric footpoint in the plane-parallel approximation. Coronal flare heating leads to increased heat conduction, chromospheric evaporation and subsequent changes in coronal pressure; these effects are included self-consistently in the models. Cooling in the chromosphere is computed in detail for the important optically thick HI, CaII and MgII transitions using the non-LTE prescription in the program MULTI. Hydrogen ionization rates from x-ray photo-ionization and collisional ionization by non-thermal electrons are included explicitly in the rate equations. The models are computed in the 'impulsive' and 'equilibrium' limits, and in a set of intermediate 'evolving' states. The impulsive atmospheres have the density distribution frozen in pre-flare configuration, while the equilibrium models assume the entire atmosphere is in hydrostatic and energetic equilibrium. The evolving atmospheres represent intermediate stages where hydrostatic equilibrium has been established in the chromosphere and corona, but the corona is not yet in energetic equilibrium with the flare heating source. Thus, for example, chromospheric evaporation is still in the process of occurring.

Hawley, Suzanne L.↗
NASA NTRS2022

4. Parallels in Communication and Navigation Technology and Natural Phenomenon

This chapter represents an introspective search for hidden connections between technology and nature and inadvertent influence on biomimicry – especially as it pertains to links to electromagnetics in navigation and communication systems. Several devices, phenomena, and approaches to problem solving are considered and an attempt is made to identify linkages between experiences, observables and the way we think about and create new technology.

Advanced Air Mobility↗
NASA NTRS1986

A distributable APSE

A distributed Ada program library is a key element in a distributed Ada Program Support Environment (APSE). To implement this successfully, the program library universe as defined by the Ada Reference Manual must be broken up into independently manageable pieces. This in turn requires the support of a distributed database system, as well as a mechanism for identifying compilation units, linkable subprograms, and Ada types in a decentralized way, to avoid falling victim to the bottlenecks of a global database and/or global unique-identifier manager. It was found that the ability to decentralize Ada program library activity is a major advantage in the management of large Ada programs. Currently, there are 18 resource-catalog revision sets, each in its own Host Interface (HIF) partition, plus 18 partitions for testing each of these, plus 11 partitions for the top-level compiler/linker/program library manager components. Compiling and other development work can proceed in parallel in each of these partitions, without suffering the performance bottlenecks of global locks or global unique-identifier generation.

Taft, Tucker, S.↗
NASA NTRS2003

Matrix Cracking in 3D Orthogonal Melt-Infiltrated SiC/SiC Composites with Various Z-Fiber Types

The occurrence of matrix cracks in melt-infiltrated SiC/SiC composites with a 3D orthogonal architecture was determined at room temperature for specimens tested in tension oriented in the X-direction (parallel to Z-bundle weave direction) and Y-direction (perpendicular to Z-bundle weave direction) and Y-direction (perpendicular to Z-bundle weave direction). The fiber-types were Sylramic and Sylramic-IBN in the X and Y-directions and lower modulus ZMI, T300, and rayon in the Z-direction. Acoustic emission (AE) was used to monitor the matrix cracking activity. For Y-direction composites, the AE data was used to determine the exact (+/- 0.25 mm) location where matrix cracks occurred in the 3D orthogonal architecture. This enabled the determination of the stress-dependent matrix crack distributions for small but repeatable matrix rich 'unidirectional' and the matrix poor 'cross-ply' regions within the architecture. It was found that matrix cracking initiated at very low stresses (approx. 40 MPa) in the 'unidirectional' regions for the largest z-direction fiber tow composites. Decreasing the size of the z-fiber bundle, increased the stress for matrix cracking in the 'unidirectional' regions. Matrix cracking in the 'cross-ply' regions always occurred at higher stresses than in 'unidirectional' regions, and the stress-dependent matrix crack distribution of the 'cross-ply' regions was always over a wider stress-range than the 'unidirectional' regions. For composites tested in the X-direction, a lower elastic modulus and a narrower and lower stress-range for matrix cracking were observed compared to composites tested in the Y-direction.

Morscher, Gregory N.↗
NASA NTRS2005

Particle acceleration magnetic field generation, and emission in Relativistic pair jets

Plasma waves and their associated instabilities (e.g., the Buneman instability, two-streaming instability, and the Weibel instability) are responsible for particle acceleration in relativistic pair jets. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic pair jet propagating through a pair plasma. Simulations show that the Weibel instability created in the collisionless shock accelerates particles perpendicular and parallel to the jet propagation direction. Simulation results show that this instability generates and amplifies highly nonuniform, small-scale magnetic fields, which contribute to the electron's transverse deflection behind the jet head. The "jitter' I radiation from deflected electrons can have different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants. The growth rate of the Weibel instability and the resulting particle acceleration depend on the magnetic field strength and orientation, and on the initial particle distribution function. In this presentation we explore some of the dependencies of the Weibel instability and resulting particle acceleration on the magnetic field strength and orientation, and the particle distribution function.

Nishikawa, K.-I.↗