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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

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

Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Recent simulations show that the Weibel instability created by relativistic pair jets is responsible for particle (electron, positron, and ion) acceleration. Using a 3-D relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic jet propagating through an ambient plasma with and without initial magnetic fields. The growth rates of the Weibel instability depends on the distribution of pair jets. The Weibel instability created in the collisionless shock accelerates particles perpendicular and parallel to the jet propagation direction. This instability is also responsible for generating and amplifying highly nonuniform, small-scale magnetic fields, which contribute to the electron s transverse deflection behind the jet head. The jitter radiation from deflected electrons has 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.

Nishikawa, K.-I.↗
NASA NTRS2013

Comparison of Precision of Biomass Estimates in Regional Field Sample Surveys and Airborne LiDAR-Assisted Surveys in Hedmark County, Norway

Airborne scanning LiDAR (Light Detection and Ranging) has emerged as a promising tool to provide auxiliary data for sample surveys aiming at estimation of above-ground tree biomass (AGB), with potential applications in REDD forest monitoring. For larger geographical regions such as counties, states or nations, it is not feasible to collect airborne LiDAR data continuously ("wall-to-wall") over the entire area of interest. Two-stage cluster survey designs have therefore been demonstrated by which LiDAR data are collected along selected individual flight-lines treated as clusters and with ground plots sampled along these LiDAR swaths. Recently, analytical AGB estimators and associated variance estimators that quantify the sampling variability have been proposed. Empirical studies employing these estimators have shown a seemingly equal or even larger uncertainty of the AGB estimates obtained with extensive use of LiDAR data to support the estimation as compared to pure field-based estimates employing estimators appropriate under simple random sampling (SRS). However, comparison of uncertainty estimates under SRS and sophisticated two-stage designs is complicated by large differences in the designs and assumptions. In this study, probability-based principles to estimation and inference were followed. We assumed designs of a field sample and a LiDAR-assisted survey of Hedmark County (HC) (27,390 km2), Norway, considered to be more comparable than those assumed in previous studies. The field sample consisted of 659 systematically distributed National Forest Inventory (NFI) plots and the airborne scanning LiDAR data were collected along 53 parallel flight-lines flown over the NFI plots. We compared AGB estimates based on the field survey only assuming SRS against corresponding estimates assuming two-phase (double) sampling with LiDAR and employing model-assisted estimators. We also compared AGB estimates based on the field survey only assuming two-stage sampling (the NFI plots being grouped in clusters) against corresponding estimates assuming two-stage sampling with the LiDAR and employing model-assisted estimators. For each of the two comparisons, the standard errors of the AGB estimates were consistently lower for the LiDAR-assisted designs. The overall reduction of the standard errors in the LiDAR-assisted estimation was around 40-60% compared to the pure field survey. We conclude that the previously proposed two-stage model-assisted estimators are inappropriate for surveys with unequal lengths of the LiDAR flight-lines and new estimators are needed. Some options for design of LiDAR-assisted sample surveys under REDD are also discussed, which capitalize on the flexibility offered when the field survey is designed as an integrated part of the overall survey design as opposed to previous LiDAR-assisted sample surveys in the boreal and temperate zones which have been restricted by the current design of an existing NFI.

Comparison↗
NASA NTRS2017

Lower Hybrid Frequency Range Waves Generated by Ion Polarization Drift Due to Electromagnetic Ion Cyclotron Waves: Analysis of an Event Observed by the Van Allen Probe B

We analyze a wave event that occurred near noon between 07:03 and 07:08 UT on 23 February 2014 detected by the Van Allen Probes B spacecraft, where waves in the lower hybrid frequency range (LHFR) and electromagnetic ion cyclotron (EMIC) waves are observed to be highly correlated, with Pearson correlation coefficient of approximately 0.86. We assume that the correlation is the result of LHFR wave generation by the ions polarization drift in the electric field of the EMIC waves. To check this assumption the drift velocities of electrons and H+, He+, and O+ ions in the measured EMIC wave electric field were modeled. Then the LHFR wave linear instantaneous growth rates for plasma with these changing drift velocities and different plasma compositions were calculated. The time distribution of these growth rates, their frequency distribution, and the frequency dependence of the ratio of the LHFR wave power spectral density (PSD)parallel and perpendicular to the ambient magnetic eld to the total PSD were found. These characteristics of the growth rates were compared with the corresponding characteristics of the observed LHFR activity. Reasonable agreement between these features and the strong correlation between EMIC and LHFR energy densities support the assumption that the LHFR wave generation can be caused by the ions polarization drift in the electric field of an EMIC wave.

Van Allen Probes↗
NASA NTRS2025

An Efficient GPU-Accelerated Multi-Source Global Fit Pipeline for LISA Data Analysis

The large-scale analysis task of deciphering gravitational wave signals in the LISA data stream will be difficult, requiring a large amount of computational resources and extensive development of computational methods. Its high dimensionality, multiple model types, and complicated noise profile require a global fit to all parameters and input models simultaneously. In this work, we detail our global fit algorithm, called “Erebor,” designed to accomplish this challenging task. It is capable of analysing current state-of-the-art datasets and then growing into the future as more pieces of the pipeline are completed and added. We describe our pipeline strategy, the algorithmic setup, and the results from our analysis of the LDC2A Sangria dataset, which contains Massive Black Hole Binaries, compact Galactic Binaries, and a parameterized noise spectrum whose parameters are unknown to the user. The Erebor algorithm includes three unique and very useful contributions: GPU acceleration for enhanced computational efficiency; ensemble MCMC sampling with multiple MCMC walkers per temperature for better mixing and parallelized sample creation; and special online updates to reversible-jump (or trans-dimensional) sampling distributions to ensure sampler mixing and accurate initial estimates for detectable sources in the data. We recover posterior distributions for all 15 (6) of the injected MBHBs in the LDC2A training (hidden) dataset. We catalog ∼12000 Galactic Binaries (∼8000 as high confidence detections) for both the training and hidden datasets. All of the sources and their posterior distributions are provided in publicly available catalogs.

LISA global fit↗
NASA NTRS2026

Influence of Processing Parameters on the Mechanical Properties of 3D Printed Borosilicate Particulate Reinforced Polymer Composites

Emerging composite materials are expanding the potential of additive manufacturing and enabling applications previously restricted by traditional manufacturing methods. The multi-phase nature of these composite materials combined with the complex in-ternal geometry of additively manufactured parts have enabled unique behavior, and potentially new applications. Additionally, these materials can be pyrolyzed to create dense metal, ceramic, and glass parts with geometries typically not achievable by tra-ditional processes. Additive manufacturing of borosilicate glass-based systems can open new applications in nuclear engineering, astronomy, and bone regrowth therapy. To elucidate the process-parameter relationship of borosilicate-polylactic acid (PLA) composites, mechanical testing was conducted and compared with a pure polylactic acid polymer baseline. Test specimens were fabricated by fused-filament fabrication with minimal post-processing. Yield strength, ultimate strength, and elastic modulus were calculated from stress-strain curves. Optical and scanning electron microscopy were conducted to observe the specimen microstructure before and after testing. The highest compressive yield strength for the composite was 28.22 MPa, and the highest compressive yield strength for PLA was 49.30 MPa. Print orientation was found to benefit the composite material but have a detrimental effect on the pure matrix material. An elastic modulus of 2.66 GPa was recorded for the borosilicate-PLA composite at 100% infill, 1 shell wall, and layer lines parallel to compression axis. Microscopy revealed that lower modulus composite specimens had the particulates re-distributed within the matrix. Tensile testing was done according to a polymer testing standard, which caused difficulties obtaining consistent fracture within the gauge length.

mechanical testing↗
NASA NTRS2001

Design and Implementation of a Parallel Multivariate Ensemble Kalman Filter for the Poseidon Ocean General Circulation Model

A multivariate ensemble Kalman filter (MvEnKF) implemented on a massively parallel computer architecture has been implemented for the Poseidon ocean circulation model and tested with a Pacific Basin model configuration. There are about two million prognostic state-vector variables. Parallelism for the data assimilation step is achieved by regionalization of the background-error covariances that are calculated from the phase-space distribution of the ensemble. Each processing element (PE) collects elements of a matrix measurement functional from nearby PEs. To avoid the introduction of spurious long-range covariances associated with finite ensemble sizes, the background-error covariances are given compact support by means of a Hadamard (element by element) product with a three-dimensional canonical correlation function. The methodology and the MvEnKF configuration are discussed. It is shown that the regionalization of the background covariances; has a negligible impact on the quality of the analyses. The parallel algorithm is very efficient for large numbers of observations but does not scale well beyond 100 PEs at the current model resolution. On a platform with distributed memory, memory rather than speed is the limiting factor.

Keppenne, Christian L.↗
NASA NTRS1993

The Generation of Oblique Magnetosonic Waves

One of the outstanding issues regarding the excitation of magnetosonic waves has been the observational evidence that obliquely propagating waves are dominant at comets and planetary foreshocks despite the predictions of linear theory that maximum growth occurs at parallel propagation. To address this issue, we have conducted a detailed linear theory using a beam-ring distribution function. The results have shown that such distribution functions are associated with four separate instabilities. Two of these instabilities are similar to the right hand resonant ion/ion and the non-resonant instabilities which are also present in the case of a field aligned beam. The other two instabilities are associated with the presence of ring part of the distribution function. One of these instabilities excites magnetosonic waves with maximum growth in the oblique directions. The other excites Alfven waves with maximum growth in the oblique directions. The importance of the former instability is that it may explain the oblique nature of the magnetosonic waves observed at planetary foreshocks and comets. In order to understand the nonlinear properties of these various instabilities, we have also conducted 2-D hybrid (particle ions, fluid electrons) simulations. We have found that when the beam density is sufficiently large so that the non-resonant instability has the largest growth rate, these waves dominate the initial wave power in the system. At later times, however, the obliquely propagating magnetosonic waves become dominant. Another important finding was that when the two instabilities which excite magnetosonic waves have the largest growth rates, the system is dominated by the obliquely propagating magnetosonic waves. This is despite the fact that the largest growth rate for one of the instabilities occurs in the parallel direction. The exact cause of this is not currently understood and is under investigation.

Source record↗
NASA NTRS1977

Unsteady supersonic aerodynamic theory for interfering surfaces by the method of potential gradient

A generalized solution of the hyperbolic wave equation was further developed to relate the velocity components at a field point to the potential gradient distribution in the dependence domain. Singular integrals were evaluated in closed form, with numerical integration methods for more complex but analytic functions. Idealization of the lifting surfaces by trapezoidal elements with two sides parallel to the streamlines is computationally efficient. Streamwise integrals were performed analytically, and spanwise integrals were neccessary only on element leading and trailing sides. All integrands vanish on the Mach cone. Pressure distribution on a double delta wing and generalized aerodynamic coefficients for three AGARD planforms were calculated and compared with available results.

Jones, W. P.↗
NASA NTRS2018

Electron Bulk Acceleration and Thermalization at Earth's Quasiperpendicular Bow Shock

Electron heating at Earth's quasiperpendicular bow shock has been surmised to be due to the combined effects of a quasistatic electric potential and scattering through wave-particle interaction. Here we report the observation of electron distribution functions indicating a new electron heating process occurring at the leading edge of the shock front. Incident solar wind electrons are accelerated parallel to the magnetic field toward downstream, reaching an electron-ion relative drift speed exceeding the electron thermal speed. The bulk acceleration is associated with an electric field pulse embedded in a whistler-mode wave. The high electron-ion relative drift is relaxed primarily through a nonlinear current-driven instability. The relaxed distributions contain a beam traveling toward the shock as a remnant of the accelerated electrons. Similar distribution functions prevail throughout the shock transition layer, suggesting that the observedacceleration and thermalization is essential to the cross-shock electron heating.

Chen, L.-J.↗
NASA NTRS1976

Test particle propagation in magnetostatic turbulence. 1. Failure of the diffusion approximation

The equation which governs the quasi-linear approximation to the ensemble and gyro-phase averaged one-body probability distribution function is constructed from first principles. This derived equation is subjected to a thorough investigation in order to calculate the possible limitations of the quasi-linear approximation. It is shown that the reduction of this equation to a standard diffusion equation in the Markovian limit can be accomplished through the application of the adiabatic approximation. A numerical solution of the standard diffusion equation in the Markovian limit is obtained for the narrow parallel beam injection. Comparison of the diabatic and adiabatic results explicitly demonstrates the failure of the Markovian description of the probability distribution function. Through the use of a linear time-scale extension the failure of the adiabatic approximation, which leads to the Markovian limit, is shown to be due to mixing of the relaxation and interaction time scales in the presence of the strong mean field.

Klimas, A. J.↗
NASA NTRS1988

Proton beam generation of oblique whistler waves

It is known that ion beams are capable of generating whistler waves that propagate parallel to the mean magnetic field. Such waves may have been observed both upstream of the earth's bow shock and in the vicinity of comets. Previous analyses are extended to include propagation oblique to the mean magnetic field. The instability is generated by the perpendicular component of free energy in the ions, which can arise either via a temperature anisotropy or via a gyrating distribution. In the former case, the generation of whistler waves is confined to a fairly narrow cone of propagation directions centered about parallel propagation; in the latter case, the maximum growth of the instability can occur at fairly large obliquities (theta equal to about 50 deg).

Wong, H. K.↗
NASA NTRS2022

Whistler Waves Generated by Nongyrotropic and Gyrotropic Electron Beams During Asymmetric Guide Field Reconnection

Using a two-dimensional particle-in-cell simulation of asymmetric reconnection with a guide field whose strength is 0.3 times the reconnecting magnetic field, we study electron distribution functions and wave intensities in the diffusion region, focusing on the electron diffusion region (EDR). Wave activities with frequencies below the electron cyclotron frequency are observed, and these are whistler waves propagating almost anti-parallel to the magnetic field. The waves are concentrated near the magnetospheric separatrix away from the X line, but the wave activity also spreads through the EDR near the X line. The reconnection outflows are asymmetric in the outflow direction in the magnetospheric side, and the wave intensity is stronger in the side of the faster electron outflow. We study the whistler waves using the fast Fourier transform, analyses of electron velocity distribution functions, and the dispersion solver calculation. Along the magnetospheric separatrix in the stronger outflow side, highly anisotropic electron beams exist with super-Alfvénic drift speeds. The dispersion analysis shows that there are two modes: a temperature anisotropy mode and a beam mode. Outside the EDR, the whistler wave intensity is highest near the separatrix, but the wave intensity decreases if we move away from the separatrix toward the magnetic neutral line because of the increase in the electron population near zero parallel velocity. In the EDR, in the velocity plane perpendicular to the magnetic field, ring/crescent electron distribution functions are observed. Near the X-line, the wave power is enhanced where nongyrotropic electrons contribute to increase the perpendicular temperature anisotropy.

magnetic reconnection↗
NASA NTRS2006

LDCM Grid Prototype (LGP)

The LGP successfully demonstrated that grid technology could be used to create a collaboration among research scientists, their science development machines, and distributed data to create a science production system in a nationally distributed environment. Grid technology provides a low cost and effective method of enabling production of science products by the science community. To demonstrate this, the LGP partnered with NASA GSFC scientists and used their existing science algorithms to generate virtual Landsat-like data products using distributed data resources. LGP created 48 output composite scenes with 4 input scenes each for a total of 192 scienes processed in parallel. The demonstration took 12 hours, which beat the requirement by almost 50 percent, well within the LDCM requirement to process 250 scenes per day. The LGP project also showed the successful use of workflow tools to automate the processing. Investing in this technology has led to funding for a ROSES ACCESS proposal. The proposal intends to enable an expert science user to produce products from a number of similar distributed instrument data sets using the Land Cover Change Community-based Processing and Analysis System (LC-ComPS) Toolbox. The LC-ComPS Toolbox is a collection of science algorithms that enable the generation of data with ground resolution on the order of Landsat-class instruments.

Weinstein, Beth↗
NASA NTRS1994

Equation solvers for distributed-memory computers

A large number of scientific and engineering problems require the rapid solution of large systems of simultaneous equations. The performance of parallel computers in this area now dwarfs traditional vector computers by nearly an order of magnitude. This talk describes the major issues involved in parallel equation solvers with particular emphasis on the Intel Paragon, IBM SP-1 and SP-2 processors.

Storaasli, Olaf O.↗
NASA NTRS1990

Development of pitch angle anisotropy and velocity diffusion of pickup ion shell distribution by solar wind turbulence

The evolution of pitch angle anisotropy in initially isotropic spherical shell distribution of pickup ions subjected to a continuous influence of weak or moderate intrinsic solar wind turbulence is investigated using a physical model which assumes that the wave of the solar wind turbulence is propagating mainly in the direction parallel to the ambient magnetic field. It is shown that, because of the pitch angle dependence of the velocity diffusion process, a significant pitch angle anisotropy of pickup ion shell distribution develops in the early stage of the diffusion process. Although it is smeared out later on, the result can be of significant importance, especially if the pickup ion density is large, because the pitch angle anisotropy can excite collective instabilities and increase the wave level, thus accelerating the diffusion process.

Yoon, Peter H.↗
NASA NTRS1997

Parallelizing Navier-Stokes Computations on a Variety of Architectural Platforms

We study the computational, communication, and scalability characteristics of a Computational Fluid Dynamics application, which solves the time accurate flow field of a jet using the compressible Navier-Stokes equations, on a variety of parallel architectural platforms. The platforms chosen for this study are a cluster of workstations (the LACE experimental testbed at NASA Lewis), a shared memory multiprocessor (the Cray YMP), distributed memory multiprocessors with different topologies-the IBM SP and the Cray T3D. We investigate the impact of various networks, connecting the cluster of workstations, on the performance of the application and the overheads induced by popular message passing libraries used for parallelization. The work also highlights the importance of matching the memory bandwidth to the processor speed for good single processor performance. By studying the performance of an application on a variety of architectures, we are able to point out the strengths and weaknesses of each of the example computing platforms.

Jayasimha, D. N.↗
NASA NTRS1971

Development of MHD flow in the entrance region of a channel.

A general method has been developed for determining the developing magnetohydrodynamic flow in the entrance region of a parallel-plate channel. Closed form solutions are obtained. Numerical results are presented for the development of velocity distribution, pressure drop, and magnetohydrodynamic entrance length at different Hartmann numbers. At large Hartmann numbers, the magnetic field effect becomes dominant and the magnetohydrodynamic entrance length becomes negligible. Also, the pressure drop is enhanced by the magnetic field. Comparison of results with direct numerical analysis shows excellent agreement.-

Wu, S. T.↗
NASA NTRS1975

Sketch of a unifying auroral theory

On the basis of field and particle observations, it is suggested that a bright auroral display is a part of a magnetosphere-ionosphere current system which is fed by a charge-separation process in the outer magnetosphere (or the solar wind). The upward magnetic-field-aligned current is flowing out of the display, carried mainly by down-flowing electrons from the hot-particle populations in the outer magnetosphere (the ambient cold electrons being depleted at high altitudes). As a result of the magnetic mirroring of these downflowing current carriers, a large potential drop is set up along the magnetic field, increasing both the number flux and the kinetic energy of the precipitating electrons. It is found that this simple basic model, when combined with wave-particle interactions, may be able to explain a highly diversified selection of auroral particle observations. It may thus be possible to explain both inverted-V events and auroral rays in terms of a static parallel electric field, and the electric field may be compatible with a strongly variable pitch-angle distribution of the precipitating electrons, including distributions peaked at 90 deg as well as 0 deg. This model may also provide a simple explanation of the simultaneous precipitation of electrons and collimated positive ions.

Lennartsson, W.↗