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,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Polycapillary X-Ray Optics for X-Ray Astronomy

High-energy x-ray spectrometry from astronomical sources is difficult due to low flux rates and high background. Polycapillary x-ray optics can be used to increase the signal-to-noise ratio by focusing photons onto a small energy dispersive detector. Two optics, both designed for x-ray astronomy, have been built and tested as collimating optics with a microfocus x-ray source for x-ray energies up to 70 keV and as focusing optics with parallel x-rays up to 50 keV. Results from these tests show that polycapillary x-ray optics have transmission efficiencies ranging from 5% to 40% and signal gains ranging from 5 to over 100 for energies up to 50 keV. Furthermore the results show that these optics have potential for use up to 100 keV. The experimental results as well as simulations for the measured and for optimized optics will be discussed.

Russell, C. H.↗
NASA NTRS2009

Stress Analysis of Bolted, Segmented Cylindrical Shells Exhibiting Flange Mating-Surface Waviness

Bolted, segmented cylindrical shells are a common structural component in many engineering systems especially for aerospace launch vehicles. Segmented shells are often needed due to limitations of manufacturing capabilities or transportation issues related to very long, large-diameter cylindrical shells. These cylindrical shells typically have a flange or ring welded to opposite ends so that shell segments can be mated together and bolted to form a larger structural system. As the diameter of these shells increases, maintaining strict fabrication tolerances for the flanges to be flat and parallel on a welded structure is an extreme challenge. Local fit-up stresses develop in the structure due to flange mating-surface mismatch (flange waviness). These local stresses need to be considered when predicting a critical initial flaw size. Flange waviness is one contributor to the fit-up stress state. The present paper describes the modeling and analysis effort to simulate fit-up stresses due to flange waviness in a typical bolted, segmented cylindrical shell. Results from parametric studies are presented for various flange mating-surface waviness distributions and amplitudes.

Knight, Norman F., Jr.↗
NASA NTRS2014

Sound Transmission Loss Through a Corrugated-Core Sandwich Panel with Integrated Acoustic Resonators

The goal of this study is to better understand the effect of structurally integrated resonators on the transmission loss of a sandwich panel. The sandwich panel has facesheets over a corrugated core, which creates long aligned chambers that run parallel to the facesheets. When ports are introduced through the facesheet, the long chambers within the core can be used as low-frequency acoustic resonators. By integrating the resonators within the structure they contribute to the static load bearing capability of the panel while also attenuating noise. An analytical model of a panel with embedded resonators is derived and compared with numerical simulations. Predictions show that acoustic resonators can significantly improve the transmission loss of the sandwich panel around the natural frequency of the resonators. In one configuration with 0.813 m long internal chambers, the diffuse field transmission loss is improved by more than 22 dB around 104 Hz. The benefit is achieved with no added mass or volume relative to the baseline structure. The embedded resonators are effective because they radiate sound out-of-phase with the structure. This results in destructive interference, which leads to less transmitted sound power.

Schiller, Noah H.↗
NASA NTRS1987

Fast simulation techniques for switching converters

Techniques for simulating a switching converter are examined. The state equations for the equivalent circuits, which represent the switching converter, are presented and explained. The uses of the Newton-Raphson iteration, low ripple approximation, half-cycle symmetry, and discrete time equations to compute the interval durations are described. An example is presented in which these methods are illustrated by applying them to a parallel-loaded resonant inverter with three equivalent circuits for its continuous mode of operation.

King, Roger J.↗
NASA NTRS1987

Energy buildup in coronal magnetic flux tubes

A time-dependent two-dimensional MHD simulation is used to study the response of the magnetic field in coronal loops to photospheric motion. From an initially uniform field, circular sections of the ends of the loop are slowly rotated to represent the photospheric motion. The evolution of the field and flow is characterized by three phases: (1) a phase of negligible kinetic energy where the current and field are predominantly parallel; (2) a phase where the field twist increases, the axial field at and near the axis increases, and the axial field decreases in two cylindrical regions away from the axis; and (3) a phase in which a significant portion of the field makes several rotations at large radii, with a corresponding reducton in the axial field to a few percent of the initial value.

Steinolfson, R. S.↗
NASA NTRS1991

Numerical simulation of ice growth on a MS-317 swept wing geometry

An effort to develop a 3-D ice accretion modeling method was initiated. This first step towards creation of a complete aircraft icing simulation code builds on previously developed methods for calculating 3-D flow fields and particle trajectories combined with a 2-D ice accretion calculation along coordinate locations corresponding to streamlines. The types of calculations necessary to predict 3-D ice accretion is demonstrated. Results of calculations using 3-D method for a MS-317 swept wing geometry are projected onto a 2-D plane parallel to the free stream direction and compared to experimental results for the same geometry. It is anticipated that many modifications will be made to this approach, however this effort will lay the groundwork for future modeling efforts. Results indicate that rime ice shapes indicate a difficulty in accurately calculating the ice shape in the runback region.

Potapczuk, M. G.↗
NASA NTRS1992

Correlation function ratios and the identification of space plasma instabilities

Wave-particle transport in a collisionless plasma is due to particle scattering by enhanced fluctuations associated with the growth of instabilities. In particular, relatively short wave-length kinetic instabilities are frequently invoked to explain many different types of plasma transport in space. Although there is an extensive theoretical and simulation literature describing the potential applications of many such instabilities, there are only a few cases of clear-cut identification of kinetic modes in space. The research described in this paper uses linear Vlasov dispersion theory to study correlation functions and dimensionless correlation function ratios for fluctuations or instabilities in three space plasma regimes. This research shows that both the compressibility and the parallel compressibility are likely to be useful in distinguishing modes in the magnetosheath as well as in the plasma sheet boundary layer and that helicity remains a useful identifier of electromagnetic ion/ion instabilities in the foreshock.

Gary, S. P.↗
NASA NTRS1992

Modeling deep convection in the Greenland Sea

The development of deep convective events in the high-latitude ocean is studied using a three-dimensional, coupled ice-ocean model. Oceanic mixing is described according to the level 2.5 turbulence closure scheme in which convection occurs in a continuous way, i.e., convective adjustment is not invoked. The model is forced by strong winds and surface cooling. Strong upwelling at the multilyear ice edge and consequent entrainment of warm Atlantic waters into the mixed layer is produced by winds parallel to the ice edge. Concomitant cooling drives deep convection and produces chimneylike structures. Inclusion of a barotropic mean flow over topography to the model provides important preconditioning and selects the location of deep convection. The most efficient preconditioning occurs at locations where the flow ascends a slope. In a stratified environment similar to the Greenland Sea with a 12 m/s wind the model simulations show that localized deep convection takes place after about 10 days to depths of 1000 m.

Hakkinen, S.↗
NASA NTRS2001

Stage Separation Performance Analysis Project

Stage separation process is an important phenomenon in multi-stage launch vehicle operation. The transient flowfield coupled with the multi-body systems is a challenging problem in design analysis. The thermodynamics environment with burning propellants during the upper-stage engine start in the separation processes adds to the complexity of the-entire system. Understanding the underlying flow physics and vehicle dynamics during stage separation is required in designing a multi-stage launch vehicle with good flight performance. A computational fluid dynamics model with the capability to coupling transient multi-body dynamics systems will be a useful tool for simulating the effects of transient flowfield, plume/jet heating and vehicle dynamics. A computational model using generalize mesh system will be used as the basis of this development. The multi-body dynamics system will be solved, by integrating a system of six-degree-of-freedom equations of motion with high accuracy. Multi-body mesh system and their interactions will be modeled using parallel computing algorithms. Adaptive mesh refinement method will also be employed to enhance solution accuracy in the transient process.

Chen, Yen-Sen↗
NASA NTRS2003

[Activities of Bay Area Research Corporation]

During the final year of this effort the HALFSHEL code was converted to work on a fast single processor workstation from it s parallel configuration. This was done because NASA Ames NAS facility stopped supporting space science and we no longer had access to parallel computer time. The single processor version of HALFSHEL was upgraded to address low density cells by using a a 3-D SOR solver to solve the equation Delta central dot E = 0. We then upgraded the ionospheric load packages to provide a multiple species load of the ionosphere out to 1.4 Rm. With these new tools we began to perform a series of simulations to address the major topic of this research effort; determining the loss rate of O(sup +) and O2(sup +) from Mars. The simulations used the nominal Parker spiral field and in one case used a field perpendicular to the solar wind flow. The simulations were performed for three different solar EUV fluxes consistent with the different solar evolutionary states believed to exist before today. The 1 EUV case is the nominal flux of today. The 3 EUV flux is called Epoch 2 and has three times the flux of todays. The 6 EUV case is Epoch 3 and has 6 times the EUV flux of today.

Source record↗
NASA NTRS2011

Introduction to: Atlantic Meridional Overturning Circulation(AMOC)

A striking conclusion of the Intergovernmental Panel on Climate Change 2007 report is the crucial role that the Atlantic Meridional Overturning Circulation (AMOC) may play in anthropogenic climate change. However, these IPCC coupled climate simulations show a broad range of uncertainty in the magnitude and timing of AMOC transport change ranging from none to nearly complete collapse within the 21st century. The potential consequences of large changes in the characteristics of AMOC have motivated the creation in the United States of an interagency program and implementation plan to develop monitoring and prediction capabilities for the AMOC This program parallels the development of substantial monitoring efforts by European, South American and African countries -- notably the UK Rapid and Rapid-Watch programs. The papers contained in this volume are derived from presentations at the First U.S. Atlantic Meridional Overturning Circulation (AMOC) Meeting held 4 - 6 May, 2009 to review the US implementation plan and its coordination with other monitoring activities. The Atlantic Meridional Overturning Circulation consists of multiple components illustrated in an attached figure. Water enters the South Atlantic at upper and intermediate depths through both western and eastern routes (where eddy transport is especially important) and is transported northward across the equator, where it recirculates within the northern subtropical and subpolar gyres. The northern end is defined by the sinking regions of the Nordic Seas and the Labrador Sea where the waters that eventually form the upper and lower branches of North Atlantic Deep Water are conditioned. High surface salinities, the result of high net evaporation in the tropics and subtropics (including the Mediterranean Sea), and presence of regions of the Arctic Ocean that remain ice-free even in winter allow for the rapid cooling and thus densification of surface water. This dense surface water becomes the source of deep water formation in the sinking regions. In addition to transporting mass, the AMOC transports roughly half of the total amount of heat carried northward through the northern subtropics (down the temperature-gradient) by the ocean. In contrast in the Southern Hemisphere AMOC transports heat up-gradient from the cool Circumpolar Current to the warm tropics. Paleoevidence suggests that AMOC heat transport in the two hemispheres has varied over time in ways intimately tied to millennial changes in the Earth's climate. In one example, the abrupt Younger Dryas spell of cold weather over the North Atlantic, which began 13,000 years ago, has generally been linked to a millennial shutdown of the AMOC as a result of massive freshwater discharge from the North American continent. The current AMOC monitoring array consists of a series of instrumented transects located across key passages (see Cunningham et al., 2010 for a recent review). In the Arctic and sub-Arctic, transects cross Fram Strait, Denmark Strait and the Faroe Channel (connecting Greenland, Iceland, and the United Kingdom), as well as the entrance to the Labrador Sea. Further south and extending outwards from the east coast of North America there are a series of monitoring arrays including arrays of the Canadian Atlantic Zone Monitoring Program, deployments of the Rapid Western Atlantic Variability Experiment (WAVE), Line W at 39 N, as well as the Rapid-MOC moored array. The latter spans the entire Atlantic basin along 26.5 N. At tropical latitudes we have the Meridional Overturning Variability Experiment (MOVE) array at 16 N, while in the Southern Hemisphere a corresponding basin-spanning transect is being established at the latitude of Cape of Good Hope, complemented by arrays at Drake Passage.

Hakkinen, Sirpa↗
NASA NTRS2022

Numerical Results for Polarized Light Scattering in a Spherical Atmosphere

We report numerical results for polarized light reflection from the top of a Rayleigh scattering spherical atmosphere with height-dependent single scattering albedo over a dark surface. Michael Mishchenko considered this scenario back in the 1990’s, for a plane-parallel atmosphere of unit optical thickness (OT), for which radiance errors arising from neglecting polarization reaches their highest values. To further extend Mishchenko’s results, we consider a value of OT = 0.25, for which the effect of atmospheric curvature is pronounced. New results are generated using three state-of-the art radiative transfer (RT) codes. These are: the MYSTIC and MCSSA models, which simulate light scattering in a true- spherical atmosphere using Monte Carlo methods; and the discrete ordinate code VLIDORT, operating with a new multiple-scatter spherical correction designed to deliver reasonable approximations to spherical-medium scattering. In this work, we report results for both single and multiple scattering; this will help to support the validation of existing and future polarized spherical RT codes, especially those using approximative methods to deal with sphericity.

S Korkin↗
NASA NTRSDate not supplied

Restoring the Dunes: Using Satellite Imagery to Help Restore Turtle Nesting Grounds In Puerto Rico After Hurricane Maria

The North coast of Puerto Rico sustained heavy damage from Hurricane Maria and other subsequent tropical storms. This damage included erosion of beaches and foredunes. The foredunes are a delicate ecological environment, which provides important nesting grounds for endangered turtles. Local ecologists have been working to restore these dunes with mixed results. Presented here is an evaluation of dune restoration sites, and strategies to improve the success of subsequent restoration efforts. This involves both historical satellite image analysis and field work to evaluate sediment migration pathways. Local ecologists use a technique called Biomimicry to restore dunes. This involves inserting sticks or plants into the ground, simulating the effects of natural vegetation to slow wind speeds and deposit wind-blown sand. Manually planted vegetation can follow to stabilize the trapped sand. While often successful, this technique occasionally fails to accumulate significant amounts of sand. Satellite imagery shows that one of the controlling factors is the orientation of the dune and beach to the prevailing wind. Dunes that are oriented parallel to the prevailing wind are relatively slower-growing due to sediment bypass. Satellite imagery also indicates that the fastest growing dunes are in close proximity to major sediment sources. This is complicated by large dams, which limit the amount of sediment delivered to the coast. Field observation shows that biomimicry matrices are occasionally located too close to the swash zone, where storms over-top the dunes, poisoning stabilizing vegetation with salt water. Another common observation was that many dunes were starting to form behind fallen trees, suggesting that a different type of biomimicry might be successful. Satellite and field observation suggest that strategies to protect nearshore environments and improve dune restoration biomimicry efforts include: 1) orient biomimicry matrices at a higher angle to prevailing winds. 2) Prioritize sites near river mouths, particularly those without upstream dams. 3) Test other orientations of biomimicry, such as horizontal boards to mimic fallen trees. 4) Improve sediment supply by preventing further dam development, and releasing trapped sediment.

Paul Bremner↗
NASA NTRS1989

Upwind-biased, point-implicit relaxation strategies for hypersonic flowfield simulations on supercomputers

An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for three-dimensional, viscous, hypersonic flows in chemical and thermal nonequilibrium is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. The relaxation strategy is well suited for computers employing either vector or parallel architectures, and the relation between computer architecture and algorithm is emphasized. It is also well suited to the numerical solution of the governing equations on unstructured grids. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large. block tri-diagonal systems. A single relaxation step depends only on information from nearest neighbors. Predictions for pressure distributions, surface heating, and aerodynamic coefficients compare well with experimental data for Mach 10 flow over a blunt body. Predictions for the hypersonic flow of air in chemical and thermal nonequilibrium (velocity = 8917 m/s, altitude = 78 km.) over the Aeroassist Flight Experiment (AFE) configuration obtained on a multi-domain grid are discussed.

Gnoffo, Peter A.↗
NASA NTRS1994

Radiative interactions in multi-dimensional chemically reacting flows using Monte Carlo simulations

The Monte Carlo method (MCM) is applied to analyze radiative heat transfer in nongray gases. The nongray model employed is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. The amount and transfer of the emitted radiative energy in a finite volume element within a medium are considered in an exact manner. The spectral correlation between transmittances of two different segments of the same path in a medium makes the statistical relationship different from the conventional relationship, which only provides the non-correlated results for nongray methods is discussed. Validation of the Monte Carlo formulations is conducted by comparing results of this method of other solutions. In order to further establish the validity of the MCM, a relatively simple problem of radiative interactions in laminar parallel plate flows is considered. One-dimensional correlated Monte Carlo formulations are applied to investigate radiative heat transfer. The nongray Monte Carlo solutions are also obtained for the same problem and they also essentially match the available analytical solutions. the exact correlated and non-correlated Monte Carlo formulations are very complicated for multi-dimensional systems. However, by introducing the assumption of an infinitesimal volume element, the approximate correlated and non-correlated formulations are obtained which are much simpler than the exact formulations. Consideration of different problems and comparison of different solutions reveal that the approximate and exact correlated solutions agree very well, and so do the approximate and exact non-correlated solutions. However, the two non-correlated solutions have no physical meaning because they significantly differ from the correlated solutions. An accurate prediction of radiative heat transfer in any nongray and multi-dimensional system is possible by using the approximate correlated formulations. Radiative interactions are investigated in chemically reacting compressible flows of premixed hydrogen and air in an expanding nozzle. The governing equations are based on the fully elliptic Navier-Stokes equations. Chemical reaction mechanisms were described by a finite rate chemistry model. The correlated Monte Carlo method developed earlier was employed to simulate multi-dimensional radiative heat transfer. Results obtained demonstrate that radiative effects on the flowfield are minimal but radiative effects on the wall heat transfer are significant. Extensive parametric studies are conducted to investigate the effects of equivalence ratio, wall temperature, inlet flow temperature, and nozzle size on the radiative and conductive wall fluxes.

Liu, Jiwen↗
NASA NTRS2018

GMI-IPS: Python Processing Software for Aircraft Campaigns

NASA's Atmospheric Tomography Mission (ATom) seeks to understand the impact of anthropogenic air pollution on gases in the Earth's atmosphere. Four flight campaigns are being deployed on a seasonal basis to establish a continuous global-scale data set intended to improve the representation of chemically reactive gases in global atmospheric chemistry models. The Global Modeling Initiative (GMI), is creating chemical transport simulations on a global scale for each of the ATom flight campaigns. To meet the computational demands required to translate the GMI simulation data to grids associated with the flights from the ATom campaigns, the GMI ICARTT Processing Software (GMI-IPS) has been developed and is providing key functionality for data processing and analysis in this ongoing effort. The GMI-IPS is written in Python and provides computational kernels for data interpolation and visualization tasks on GMI simulation data. A key feature of the GMI-IPS, is its ability to read ICARTT files, a text-based file format for airborne instrument data, and extract the required flight information that defines regional and temporal grid parameters associated with an ATom flight. Perhaps most importantly, the GMI-IPS creates ICARTT files containing GMI simulated data, which are used in collaboration with ATom instrument teams and other modeling groups. The initial main task of the GMI-IPS is to interpolate GMI model data to the finer temporal resolution (1-10 seconds) of a given flight. The model data includes basic fields such as temperature and pressure, but the main focus of this effort is to provide species concentrations of chemical gases for ATom flights. The software, which uses parallel computation techniques for data intensive tasks, linearly interpolates each of the model fields to the time resolution of the flight. The temporally interpolated data is then saved to disk, and is used to create additional derived quantities. In order to translate the GMI model data to the spatial grid of the flight path as defined by the pressure, latitude, and longitude points at each flight time record, a weighted average is then calculated from the nearest neighbors in two dimensions (latitude, longitude). Using SciPya's Regular Grid Interpolator, interpolation functions are generated for the GMI model grid and the calculated weighted averages. The flight path points are then extracted from the ATom ICARTT instrument file, and are sent to the multi-dimensional interpolating functions to generate GMI field quantities along the spatial path of the flight. The interpolated field quantities are then written to a ICARTT data file, which is stored for further manipulation. The GMI-IPS is aware of a generic ATom ICARTT header format, containing basic information for all flight campaigns. The GMI-IPS includes logic to edit metadata for the derived field quantities, as well as modify the generic header data such as processing dates and associated instrument files. The ICARTT interpolated data is then appended to the modified header data, and the ICARTT processing is complete for the given flight and ready for collaboration. The output ICARTT data adheres to the ICARTT file format standards V1.1. The visualization component of the GMI-IPS uses Matplotlib extensively and has several functions ranging in complexity. First, it creates a model background curtain for the flight (time versus model eta levels) with the interpolated flight data superimposed on the curtain. Secondly, it creates a time-series plot of the interpolated flight data. Lastly, the visualization component creates averaged 2D model slices (longitude versus latitude) with overlaid flight track circles at key pressure levels. The GMI-IPS consists of a handful of classes and supporting functionality that have been generalized to be compatible with any ICARTT file that adheres to the base class definition. The base class represents a generic ICARTT entry, only defining a single time entry and 3D spatial positioning parameters. Other classes inherit from this base class; several classes for input ICARTT instrument files, which contain the necessary flight positioning information as a basis for data processing, as well as other classes for output ICARTT files, which contain the interpolated model data. Utility classes provide functionality for routine procedures such as: comparing field names among ICARTT files, reading ICARTT entries from a data file and storing them in data structures, and returning a reduced spatial grid based on a collection of ICARTT entries. Although the GMI-IPS is compatible with GMI model data, it can be adapted with reasonable effort for any simulation that creates Hierarchical Data Format (HDF) files. The same can be said of its adaptability to ICARTT files outside of the context of the ATom mission. The GMI-IPS contains just under 30,000 lines of code, eight classes, and a dozen drivers and utility programs. It is maintained with GIT source code management and has been used to deliver processed GMI model data for the ATom campaigns that have taken place to date.

Damon, M. R.↗
NASA NTRS1980

Overseas trip report, CV 990 underflight mission

The scanning microwave radiometer-7 simulator, the ocean temperature scanner, and an imaging scatterometer/altimeter operating at 14 GHz were carried onboard the NASA CV-990 over open oceans, sea ice, and continental ice sheets to gather surface truth information. Data flights were conducted over the Norwegian Sea to map the ocean polar front south and west of Bear Island and to transect several Nimbus-7 footprints in a rectangular pattern parallel to the northern shoreline of Norway. Additional flights were conducted to obtain correlative data on the cryosphere parameters and characteristics of the Greenland ice sheet, and study the frozen lakes near Barrow. The weather conditions and flight path way points for each of the nineteen flights are presented in tables and maps.

Gloersen, P.↗
NASA NTRS1985

Putting the 'super' in supercomputers

Computers used for numerical simulations of physical phenomena, e.g., flowfields, meteorology, structural analysis, etc., replace physical experiments that are too expensive or impossible to perform. The problems considered continually become increasingly more complex and thus demand faster processing times to do all necessary computations. The effects components technologies have on computer speed are leveling off, leaving new architectures and programming as the only currently viable means to upgrade speed. Parallel computations, either in the form of array processors, assembly line processing or multiprocessors are being explored using existing microprocessor technologies. Slower hardware configurations can also be made equivalent to faster supercomputers by economic programming. The availability of rudimentary parallel architecture supercomputers for general industrial use is increasing. Scientific applications continue to drive the development of more sophisticated parallel machines.

Schulbach, C.↗