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,045 records · Page 58

Cell module and fuel conditioner

Measurements of stack height changes with temperature and cell material characteristics were made. Stack 559 was assembled and components were fabricated for 560, 561, and 562. Stack 425 was transferred from the parallel DOE program and installed in the OS/IES simulation loop for mechanical and electrical testing. Construction and preliminary checkout of the 2 kW test facility was completed and design and procurement of the 8 kW test facility was initiated. The fuel conditioning subsystem design continued to evolve and the state points for the current design were calculated at full and part load conditions. Steam reforming catalyst activity tests were essentially completed and aging tests and CO shift converter tests were initiated.

Hoover, D. Q., Jr.↗

Monte Carlo simulation of charged particles upstream of the earth's bow shock

A simple model of the earth's bow shock is presented and compared with ISEE-1 observations of diffuse ions upstream of the shock. Ion behavior is modeled in a test particle simulation using Monte Carlo techniques for a plane parallel shock at which cold solar wind ions are heated sufficiently to allow diffusion back upstream. Fitting of the model to the observations indicates that the diffusion coefficient must increase faster than linearly with velocity, and provides an upper limit to the acceleration efficiency defined as the proportion of incoming solar wind ions accelerated and observed upstream. Results suggest that astrophysical shocks are capable of accelerating ions directly from the thermal plasma without the need for a distinct preenergizing injection stage.

Ellison, D. C.↗

Quasi non-intrusive sampling and analysis of gases associated with the boundary layer on the tethered satellite and similar supersonic and hypersonic research vehicles

The effect that one candidate inlet, an assembly of capillary openings in a thin glass plate (a multichannel plate), might have on the overall sampling process is considered. The flow characteristics of the plate, under a variety of conditions of external pressure and mass flow, is evaluated. A review of capillary flow theory is presented with some development of the pertinent equations. The predicted mass flow will be compared to that determined perimentally to assess the effect that mass flow through one capillary might have upon a neighboring capillary. Mass spectrometric measurements of mixtures flowing through the multichannel plate (MCP) is also considered. In the first part of the experiments, the flow was in a direction normal to the surface of the plate. The experimental gases were Argon and mixtures of carbon dioxide in air. Ongoing experiments are discussed which are conducted with the flow parallel to the surface of the plate, a simulation of the kind of flow that a flight vehicle would experience.

Fishel, C.↗

Numerical experiments on the stability of controlled boundary layers

Nonlinear simulations are presented for instability and transition in parallel water boundary layers subjected to pressure gradient, suction, or heating control. In the nonlinear regime, finite amplitude, 2-D Tollmein-Schlichting waves grow faster than is predicted by linear theory. Moreover, this discrepancy is greatest in the case of heating control. Likewise, heating control is found to be the least effective in delaying secondary instabilities of both the fundamental and subharmonic type. Flow field details (including temperature profiles) are presented for both the uncontrolled boundary layer and the heated boundary layer.

Zang, Thomas A.↗

A mathematical model of the structure and evolution of small scale discrete auroral arcs

A three dimensional fluid model which includes the dispersive effect of electron inertia is used to study the nonlinear macroscopic plasma dynamics of small scale discrete auroral arcs within the auroral acceleration zone and ionosphere. The motion of the Alfven wave source relative to the magnetospheric and ionospheric plasma forms an oblique Alfven wave which is reflected from the topside ionosphere by the negative density gradient. The superposition of the incident and reflected wave can be described by a steady state analytical solution of the model equations with the appropriate boundary conditions. This two dimensional discrete auroral arc equilibrium provides a simple explanation of auroral acceleration associated with the parallel electric field. Three dimensional fully nonlinear numerical simulations indicate that the equilibrium arc configuration evolves three dimensionally through collisionless tearing and reconnection of the current layer. The interaction of the perturbed flow and the transverse magnetic field produces complex transverse structure that may be the origin of the folds and curls observed to be associated with small scale discrete arcs.

Seyler, C. E.↗

On the wall-normal velocity of the compressible boundary-layer equations

Numerical methods for the compressible boundary-layer equations are facilitated by transformation from the physical (x,y) plane to a computational (xi,eta) plane in which the evolution of the flow is 'slow' in the time-like xi direction. The commonly used Levy-Lees transformation results in a computationally well-behaved problem for a wide class of non-similar boundary-layer flows, but it complicates interpretation of the solution in physical space. Specifically, the transformation is inherently nonlinear, and the physical wall-normal velocity is transformed out of the problem and is not readily recovered. In light of recent research which shows mean-flow non-parallelism to significantly influence the stability of high-speed compressible flows, the contribution of the wall-normal velocity in the analysis of stability should not be routinely neglected. Conventional methods extract the wall-normal velocity in physical space from the continuity equation, using finite-difference techniques and interpolation procedures. The present spectrally-accurate method extracts the wall-normal velocity directly from the transformation itself, without interpolation, leaving the continuity equation free as a check on the quality of the solution. The present method for recovering wall-normal velocity, when used in conjunction with a highly-accurate spectral collocation method for solving the compressible boundary-layer equations, results in a discrete solution which is extraordinarily smooth and accurate, and which satisfies the continuity equation nearly to machine precision. These qualities make the method well suited to the computation of the non-parallel mean flows needed by spatial direct numerical simulations (DNS) and parabolized stability equation (PSE) approaches to the analysis of stability.

Pruett, C. David↗

Propagation characteristics of waves upstream and downstream of quasi-parallel shocks

The propagation characteristics of waves upstream and downstream of quasi-parallel shocks are investigated by using 2D hybrid simulations. At low Alfven Mach numbers, M(A) below about 2, the shock is initially associated with upstream phase-standing whistlers. At later times, backstreaming ions excite longer-wavelength whistlers via the right-hand resonant ion/ion instability. These waves propagate along the magnetic field at a group velocity no smaller than the upstream flow speed, so that the waves remain in the upstream region. At higher MA (above about 3), these waves are convected back into the shock, causing its reformation and downstream perturbations. Shock transmitted waves mode-convert into Alfven/ion-cyclotron waves which have a wave vector along the shock normal (pointing upstream) and convect downstream. The 2D simulation results confirm our earlier suggestion that the upstream waves should be field aligned, and that their convection into the downstream is associated with linear mode conversion into the Alfven/ion-cyclotron branch.

Krauss-Varban, D.↗

The differential emission measure of nested hot and cool magnetic loops

The detailed thermal structure of the magnetized solar transition region, as measured by its differential emission measure (DEM(T)), is poorly known. Building on the fact that the solar surface is strongly magnetized and thereby structured, proposals have been made that envision a significant lower-temperature contribution to the energy balance from (ion) heat flux across an arcade of different temperature loops. In this paper, we describe a self-consistent 2D MHD simulation, which includes the full thermal effects of parallel stability and anisotropic conduction, of a nested-loop model of the thermal and magnetic structure of the transition region. We then demonstrate that the predicted DEM agrees with observations in the conceptually elusive T less than 10 exp 5 K regime.

Van Hoven, G.↗

First CLIPS Conference Proceedings, volume 2

The topics of volume 2 of First CLIPS Conference are associated with following applications: quality control; intelligent data bases and networks; Space Station Freedom; Space Shuttle and satellite; user interface; artificial neural systems and fuzzy logic; parallel and distributed processing; enchancements to CLIPS; aerospace; simulation and defense; advisory systems and tutors; and intelligent control.

Source record↗

Drop Breakup in Fixed Bed Flows as Model Stochastic Flow Fields

We examine drop breakup in a class of stochastic flow fields as a model for the flow through fixed fiber beds and to elucidate the general mechanisms whereby drops breakup in disordered, Lagrangian unsteady flows. Our study consists of two parallel streams of investigation. First, large scale numerical simulations of drop breakup in a class of anisotropic Gaussian fields will be presented. These fields are generated spectrally and have been shown in a previous publication to be exact representations of the flow in a dilute disordered bed of fibers if close interactions between the fibers and the drops are dynamically unimportant. In these simulations the drop shape is represented by second and third order small deformation theories which have been shown to be excellent for the prediction of drop breakup in steady strong flows. We show via these simulations that the mechanisms of drop breakup in these flows are quite different than in steady flows. The predominant mechanism of breakup appears to be very short lived twist breakups. Moreover, the occurrence of breakup events is poorly predicted by either the strength of the local flow in which the drop finds itself at breakup, or the degree of deformation that the drop achieves prior to breakup. It is suggested that a correlation function of both is necessary to be predictive of breakup events. In the second part of our research experiments are presented where the drop deformation and breakup in PDMS/polyisobutylene emulsions is considered. We consider very dilute emulsions such that coalescence is unimportant. The flows considered are simple shear and the flow through fixed fiber beds. Turbidity, small angle light scattering, dichroism and microscopy are used to interrogate the drop deformation process in both flows. It is demonstrated that breakup at very low capillary numbers occurs in both flows but larger drop deformation occurs in the fixed bed flow. Moreover, it is witnessed that breakup in the bed occurs continuously during flow and apparently with uniform probability through the bed length. The drop deformations witnessed in our experiments are larger than those predicted by the numerical simulations, and future plans to investigate these differences are discussed.

Shaqfeh, Eric S. G.↗

A New Approach to Parallel Interference Cancellation for CDMA

This paper introduces an improved nonlinear parallel interference cancellation scheme that significantly reduces the degrading effect of user interference with implementation complexity linear in the number of users. The scheme operates on the fact that parallel processing simultaneously removes from each user a part of the interference produced by the remaining users accessing the channel the amount being proportional to their reliability. The parallel processing can be done in multiple stages. Simulation results are given for a multitude of different situations, in particular those cases for which the analysis is too complex.

interference cancellation multiuser communication ↗

High Resolution Aerospace Applications using the NASA Columbia Supercomputer

This paper focuses on the parallel performance of two high-performance aerodynamic simulation packages on the newly installed NASA Columbia supercomputer. These packages include both a high-fidelity, unstructured, Reynolds-averaged Navier-Stokes solver, and a fully-automated inviscid flow package for cut-cell Cartesian grids. The complementary combination of these two simulation codes enables high-fidelity characterization of aerospace vehicle design performance over the entire flight envelope through extensive parametric analysis and detailed simulation of critical regions of the flight envelope. Both packages. are industrial-level codes designed for complex geometry and incorpor.ats. CuStomized multigrid solution algorithms. The performance of these codes on Columbia is examined using both MPI and OpenMP and using both the NUMAlink and InfiniBand interconnect fabrics. Numerical results demonstrate good scalability on up to 2016 CPUs using the NUMAIink4 interconnect, with measured computational rates in the vicinity of 3 TFLOP/s, while InfiniBand showed some performance degradation at high CPU counts, particularly with multigrid. Nonetheless, the results are encouraging enough to indicate that larger test cases using combined MPI/OpenMP communication should scale well on even more processors.

Mavriplis, Dimitri J.↗

Efficient Parallel Engineering Computing on Linux Workstations

A C software module has been developed that creates lightweight processes (LWPs) dynamically to achieve parallel computing performance in a variety of engineering simulation and analysis applications to support NASA and DoD project tasks. The required interface between the module and the application it supports is simple, minimal and almost completely transparent to the user applications, and it can achieve nearly ideal computing speed-up on multi-CPU engineering workstations of all operating system platforms. The module can be integrated into an existing application (C, C++, Fortran and others) either as part of a compiled module or as a dynamically linked library (DLL).

Lou, John Z.↗

A 4-D Climatology (1979-2009) of the Monthly Tropospheric Aerosol Optical Depth Distribution over the Mediterranean Region from a Comparative Evaluation and Blending of Remote Sensing and Model Products

Since the 1980s several spaceborne sensors have been used to retrieve the aerosol optical depth (AOD) over the Mediterranean region. In parallel, AOD climatologies coming from different numerical model simulations are now also available, permitting to distinguish the contribution of several aerosol types to the total AOD. In this work, we perform a comparative analysis of this unique multiyear database in terms of total AOD and of its apportionment by the five main aerosol types (soil dust, seasalt, sulfate, black and organic carbon). We use 9 different satellite-derived monthly AOD products: NOAA/AVHRR, SeaWiFS (2 products), TERRA/MISR, TERRA/MODIS, AQUA/MODIS, ENVISAT/MERIS, PARASOL/POLDER and MSG/SEVIRI, as well as 3 more historical datasets: NIMBUS7/CZCS, TOMS (onboard NIMBUS7 and Earth- Probe) and METEOSAT/MVIRI. Monthly model datasets include the aerosol climatology from Tegen et al. (1997), the climate-chemistry models LMDz-OR-INCA and RegCM-4, the multi-model mean coming from the ACCMIP exercise, and the reanalyses GEMS and MACC. Ground-based Level- 2 AERONET AOD observations from 47 stations around the basin are used here to evaluate the model and satellite data. The sensor MODIS (on AQUA and TERRA) has the best average AOD scores over this region, showing a relevant spatio-temporal variability and highlighting high dust loads over Northern Africa and the sea (spring and summer), and sulfate aerosols over continental Europe (summer). The comparison also shows limitations of certain datasets (especially MERIS and SeaWiFS standard products). Models reproduce the main patterns of the AOD variability over the basin. The MACC reanalysis is the closest to AERONET data, but appears to underestimate dust over Northern Africa, where RegCM-4 is found closer to MODIS thanks to its interactive scheme for dust emissions. The vertical dimension is also investigated using the CALIOP instrument. This study confirms differences of vertical distribution between dust aerosols showing a large vertical spread, and other continental and marine aerosols which are confined in the boundary layer. From this compilation, we propose a 4-D blended product from model and satellite data, consisting in monthly time series of 3-D aerosol distribution at a 50 km horizontal resolution over the Euro-Mediterranean marine and continental region for the 2003-2009 period. The product is based on the total AOD from AQUA/MODIS, apportioned into sulfates, black and organic carbon from the MACC reanalysis, and into dust and sea-salt aerosols from RegCM-4 simulations, which are distributed vertically based on CALIOP climatology.We extend the 2003-2009 reconstruction to the past up to 1979 using the 2003-2009 average and applying the decreasing trend in sulfate aerosols from LMDz-OR-INCA, whose AOD trends over Europe and the Mediterranean are median among the ACCMIP models. Finally optical properties of the different aerosol types in this region are proposed from Mie calculations so that this reconstruction can be included in regional climate models for aerosol radiative forcing and aerosol-climate studies.

troposphere↗

Resource Prospector Propulsion System Cold Flow Testing

Resource Prospector (RP) is a NASA mission being led by NASA Ames Research Center with current plans to deliver a scientific payload package aboard a rover to the lunar surface. As part of an early risk reduction activity, Marshall Space Flight Center (MSFC) and Johnson Space Flight Center (JSC) have jointly developed a government-version concept of a lunar lander for the mission. The spacecraft consists of two parts, the lander and the rover which carries the scientific instruments. The lander holds the rover during launch, cruise, and landing on the surface. Following terminal descent and landing the lander portion of the spacecraft become dormant after the rover embarks on the science mission. The lander will be equipped with a propulsion system for lunar descent and landing, as well as trajectory correction and attitude control maneuvers during transit to the moon. Hypergolic propellants monomethyl hydrazine and nitrogen tetroxide will be used to fuel sixteen 70-lbf descent thrusters and twelve 5-lbf attitude control thrusters. A total of four metal-diaphragm tanks, two per propellant, will be used along with a high-pressure composite-overwrapped pressure vessel for the helium pressurant gas. Many of the major propulsion system components are heritage missile hardware obtained by NASA from the Air Force. In parallel with the flight system design activities, a simulated propulsion system based on flight drawings was built for conducting a series of water flow tests to characterize the transient fluid flow of the propulsion system feed lines and to verify the critical operation modes such as system priming, waterhammer, and crucial mission duty cycles. The primary objective of the cold flow testing was to simulate the RP propulsion system fluid flow operation through water flow testing and to obtain data for anchoring analytical models. The models will be used to predict the transient and steady state flow behaviors in the actual flight operations. All design and build efforts, including the analytical modeling, have been performed. The cold flow testing of the propulsion system was set up and conducted at a NASA MSFC test facility. All testing was completed in the summer of 2014, and this paper documents the results of that testing and the associated fluid system modeling efforts.

Williams, Hunter↗

A New Code SORD for Simulation of Polarized Light Scattering in the Earth Atmosphere

We report a new publicly available radiative transfer (RT) code for numerical simulation of polarized light scattering in plane-parallel atmosphere of the Earth. Using 44 benchmark tests, we prove high accuracy of the new RT code, SORD (Successive ORDers of scattering). We describe capabilities of SORD and show run time for each test on two different machines. At present, SORD is supposed to work as part of the Aerosol Robotic NETwork (AERONET) inversion algorithm. For natural integration with the AERONET software, SORD is coded in Fortran 90/95. The code is available by email request from the corresponding (first) author or from ftp://climate1.gsfc.nasa.gov/skorkin/SORD/.

AERONET↗

New Features of the NEQAIR Radiation Code

The longest-lived code for predicting shock layer radiation, NEQAIR, is now in its 5th decade of service. Substantial changes to the code have been made over the previous decade, the most recent report of which was at the 5th Workshop on Radiation in High Temperature Gases in 2014, for the version referred to as NEQAIR14. This paper will review some of the improvements made to the NEQAIR code since then, which is now at v15.2. Some of these features are discussed briefly below. NEQAIR15 and subsequent versions have enabled parallel evaluation of multiple lines of sight. This is accomplished by utilizing the HDF5 file format and placing multiple lines into a single file, LOS.h5, which is used for both input and output. This approach enables straightforward parallel execution both over the number of lines of sight and the number of points per line. For large problems, runtime reduces linearly with the number of nodes deployed since each line is processed independently by a subset of MPI ranks. Three applications of the multi-line solver are discussed. The first has to do with performing loosely coupled radiation-flowfield solutions. In this case the computed absorption and emission coefficients are used to evaluate the total energy absorbed or emitted at each point, allowing evaluation of the volumetric source term in the flowfield. The second computation is for obtaining heat flux from nonuniform flows, which require integration over spherical co-ordinates. These are of particular interest for evaluating radiation on the vehicle backshell. This 3D option improves the angular integration scheme and allows adaptive line selection that together reduce the number of lines required by about an order of magnitude. The final application is for remote observation, which is essentially the 3D integration problem over a small solid angle. For all three of these computations, data can be stored in the HDF5 file which allows a NEQAIR run to be restarted when it times out, or to add atmospheric absorption or instrument scan functions. An additional level of parallelism is enabled in NEQAIR15.2 using GPU routines. The GPU parallelism has realized up to 8x speed-up when running on a single core but diminishes as CPU parallelism is increased. For running multi-line simulations, it may be easier to reserve a large number of CPU nodes than to obtain the number of GPU nodes required for similar performance. A GUI, known as NEQTPY, allows for reading and creating input files, running NEQAIR, and displaying results. A significant feature of NEQTPY is the ability to perform spectral fits to data. The fits can operate on a single line spectrum (radiance vs. wavelength) or a 3D input file with multiple columns of data. Other new features include improved constants, additional species, more detailed non-Boltzmann modelling, advanced user controls, the ability to read and calculate spectra from HITRAN datafiles, photodissociation and photoionization cross-sections. A “fast” automatic grid option may reduce the size and time of spectral calculations while still maintaining good accuracy for total heat flux.

Brett A Cruden↗

Editorial: Flux Rope Interaction With the Ambient Corona: From Jets to CMEs

Eruptive events within the Sun’s corona occur across a broad range of scales, from abundant small-scale jets to highly energetic coronal mass ejections (CMEs) (e.g., Webb and Howard, 2012; Raouafi et al., 2016; Kumar et al., 2021). Flux ropes have been understood to be a fundamental constituent of CMEs for many years, but it has only been more recently that their role in smaller eruptive events has become more appreciated (e.g., Sterling et al., 2015; Wyper et al., 2017). The key to understanding the differing morphology and nature of eruptions on these vastly differing scales is to understand the nature of the interaction between the flux ropes involved and the magnetic field of the surrounding corona. This Research Topic invited perspectives on (and examples of) flux rope eruptions across this broad range of scales, from formation to ejection, with the aim of highlighting commonalities and differences to aid in ultimately building a common framework for their understanding. Starting at the smallest scales, Schmieder presented a historically ordered review of jets and surges. The parallel improvements in instrumentation, physical insight and numerical simulations over time were highlighted along with the ways in which improvements in one area led to developments in another. Overall, jets of all forms are presented as the result of magnetic reconnection between two magnetic domains, with the reconnection involving a null point, bald patch separatrix or quasi-separatrix layer. In particular, jets with observed rotational motions often involve the transfer of twist from a small-scale filament or flux rope by reconnection across the separatrix.

solar corona and wind↗