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,693 records · Page 94

Evaporating Spray in Supersonic Streams Including Turbulence Effects

Evaporating spray plays an important role in spray combustion processes. This paper describes the development of a new finite-conductivity evaporation model, based on the two-temperature film theory, for two-phase numerical simulation using Eulerian-Lagrangian method. The model is a natural extension of the T-blob/T-TAB atomization/spray model which supplies the turbulence characteristics for estimating effective thermal diffusivity within the droplet phase. Both one-way and two-way coupled calculations were performed to investigate the performance of this model. Validation results indicate the superiority of the finite-conductivity model in low speed parallel flow evaporating sprays. High speed cross flow spray results indicate the effectiveness of the T-blob/T-TAB model and point to the needed improvements in high speed evaporating spray modeling.

Balasubramanyam, M. S.↗

Assessment of Polarization Effect on Efficiency of Levenberg-Marquardt Algorithm in Case of Thin Atmosphere Over Black Surface

The Levenberg-Marquardt algorithm [1, 2] provides a numerical iterative solution to the problem of minimization of a function over a space of its parameters. In our work, the Levenberg-Marquardt algorithm retrieves optical parameters of a thin (single scattering) plane parallel atmosphere irradiated by collimated infinitely wide monochromatic beam of light. Black ground surface is assumed. Computational accuracy, sensitivity to the initial guess and the presence of noise in the signal, and other properties of the algorithm are investigated in scalar (using intensity only) and vector (including polarization) modes. We consider an atmosphere that contains a mixture of coarse and fine fractions. Following [3], the fractions are simulated using Henyey-Greenstein model. Though not realistic, this assumption is very convenient for tests [4, p.354]. In our case it yields analytical evaluation of Jacobian matrix. Assuming the MISR geometry of observation [5] as an example, the average scattering cosines and the ratio of coarse and fine fractions, the atmosphere optical depth, and the single scattering albedo, are the five parameters to be determined numerically. In our implementation of the algorithm, the system of five linear equations is solved using the fast Cramer s rule [6]. A simple subroutine developed by the authors, makes the algorithm independent from external libraries. All Fortran 90/95 codes discussed in the presentation will be available immediately after the meeting from sergey.v.korkin@nasa.gov by request.

Korkin, S.↗

Lower-Hybrid Wave Structures and Interactions With Electrons Observed in Magnetotail Reconnection Diffusion Regions

We investigate waves close to the lower-hybrid frequency in 12 magnetotail reconnection electron diffusion region (EDR) events with guide field levels of near-zero to 30%. In about half of the events, the wave vector has a small component along the current sheet normal, consistent with known lower-hybrid drift wave properties, but the perpendicular magnetic field fluctuations can be comparable or greater than the parallel component, a feature unique to the waves inside and adjacent to EDRs. Another new wave property is that the wave vector has a significant component along the current sheet normal in some events and completely along the normal for one event. In 1/4of the events, the ∇∙𝑷(sub 𝑒) term has a significant contribution to the wave electric field, possibly a feature of lower-hybrid waves more likely to exist in the diffusion region than further away from the X-line. Electron temperature variations are correlated with the wave potential, due to wave electric field acceleration and crossings at the corrugated separatrix region with different amounts of mixing between reconnection inflowing and outflowing populations. The latter also leads to the anti-correlation between parallel and perpendicular temperature components. Using four-spacecraft measurements, the magnetic field line twisting is demonstrated by the correlated fluctuations in (∇×𝑽(sub 𝐸×𝐵))||and (∇×𝐁)||. The lower-hybrid wave in the EDR of weak guide field reconnection may be generated near separatrices and penetrate to the mid-plane or locally generated, and the latter possibility is beyond the prediction of previous reconnection simulations.

Magnetic reconnection↗

The tail lobe ion spectrometer

A two-dimensional kinetic 'tail lobe ion spectrometer' model for the transport of ionospheric ions from polar cleft ionosphere into the tail lobes is developed to semiquantitatively simulate the behavior of the observed O(+) ion streams in the magnetotail lobes. The consequences of the present model include: (1) the increase in velocity of the O(+) streams away from the tail midplane; (2) the existence of 'tongues' of O(+) density from the cleft ionosphere into the tail lobes whose distribution configuration depends on the convection electric field and the source thermal energies; and (3) very low parallel thermal energies and speeds in the lobes which are correlated with the parallel ion bulk velocities. These consequences are reasonably consistent with the trends in lobe streams observed by recent spacecraft.

Horwitz, J. L.↗

Applications considerations in the system design of highly concurrent multiprocessors

A flow model processor approach to parallel processing is described, using very-high-performance individual processors, high-speed circuit switched interconnection networks, and a high-speed synchronization capability to minimize the effect of the inherently serial portions of applications on performance. Design studies related to the determination of the number of processors, the memory organization, and the structure of the networks used to interconnect the processor and memory resources are discussed. Simulations indicate that applications centered on the large shared data memory should be able to sustain over 500 million floating point operations per second.

Lundstrom, Stephen F.↗

Composite load spectra for select space propulsion structural components

The objective of this program is to develop generic load models with multiple levels of progressive sophistication to simulate the composite load spectra that are induced in space propulsion system components, representative of Space Shuttle Main Engines (SSME), such as transfer ducts, turbine blades, and liquid oxygen (LOX) posts and system ducting. These models will be developed using two independent approaches. The first approach consists of using state-of-the-art probabilistic methods to describe the individual loading conditions and combinations of these loading conditions to synthesize the composite load spectra simulation. The methodology required to combine the various individual load simulation models (hot-gas dynamic, vibrations, instantaneous position, centrifugal field, etc.) into composite load spectra simulation models will be developed under this program. A computer code incorporating the various individual and composite load spectra models will be developed to construct the specific load model desired. The second approach, which is covered under the options portion of the contract, will consist of developing coupled models for composite load spectra simulation which combine the (deterministic) models for composite load dynamic, acoustic, high-pressure and high rotational speed, etc., load simulation using statistically varying coefficients. These coefficients will then be determined using advanced probabilistic simulation methods with and without strategically selected experimental data. This report covers the efforts of the third year of the contract. The overall program status is that the turbine blade loads have been completed and implemented. The transfer duct loads are defined and are being implemented. The thermal loads for all components are defined and coding is being developed. A dynamic pressure load model is under development. The parallel work on the probabilistic methodology is essentially completed. The overall effort is being integrated in an expert system code specifically developed for this project.

Newell, J. F.↗

Direct Simulation of Evolution and Control of Nonlinear Instabilities in Attachment-Line Boundary Layers

The unsteady, incompressible Navier-Stokes equations are used for the direct numerical simulation (DNS) of spatially evolving disturbances in a three-dimensional (3-D) attachment-line boundary layer. Two-dimensional (2-D) disturbances are introduced either by forcing at the in ow or by harmonic-source generators at the wall; 3-D disturbances are introduced by harmonic-source generators at the wall. The DNS results are in good agreement with both 2-D non-parallel theory (for small-amplitude disturbances) and weakly nonlinear theory (for finite-amplitude disturbances), which validates the two theories. The 2-D DNS results indicate that nonlinear disturbance growth occurs near branch II of the neutral stability curve; however, steady suction can be used to stabilize this disturbance growth. For 3-D instabilities that are generated o the attachment line, spreading both toward and away from the attachment line causes energy transfer to the attachment-line and downstream instabilities; suction stabilizes these instabilities. Furthermore, 3-D instabilities are more stable than 2-D or quasi-2-D instabilities.

Joslin, Ronald D.↗

Cross-Scale Coupling in the Inner Magnetosphere

Magnetosphere-ionosphere (MI) coupling has interested scientists for decades and, in spite of experimental and theoretical research efforts, is still one of the least well-known dynamic processes in space plasma. The reason for this is that the numerous physical processes associated with MI coupling occur over multiple spatial lengths and temporal scales. One typical example of MI coupling is small- and large-scale ring current (RC) electrodynamic coupling. In this talk, we will address the two primary issues of RC electrodynamic coupling: (1) RC self-consistent coupling with electromagnetic ion cyclotron (EMIC) waves (small-scale electrodynamic coupling) and (2) RC self-consistent MI coupling that includes calculation of the magnetospheric electric field (large-scale electrodynamic coupling). We also will emphasize the role of the heavy ions in the number of wave-particle interaction magnetospheric processes. In particular, we will discuss some of the experimental and theoretical studies that have investigated the role of the heavy ions (mainly He(+) and O(+)) in generation and propagation of electromagnetic ion cyclotron waves and their contribution to the heating of magnetospheric electrons and ions. The more recent studies have also shown that the heavy ions can greatly contribute to a generation of lower hybrid waves, ring current precipitation phenomena, and the overall energy redistribution in the inner magnetosphere. Using newly developed 2.5-dimensional particle-in-cell simulations, we study the energization and nonlinear coupling of different plasma waves in the presence of the heavy ions. We have shown that the high frequency wave modes critically depend on the heavy ion density and irrespective of the driven wave modes, both the light and heavy ions undergo significant transverse acceleration. But for the large heavy-ion densities, even the electrons are significantly accelerated in the parallel direction by the waves below the LH frequency.

Khazanov, G. V.↗

The Role of the Heavy Ions in the Wave Magnetospheric Phenomena

This talk will emphasize the role of the heavy ions in the number of wave-particle interaction magnetospheric processes. In particular, we will discuss some of the experimental and theoretical studies that have investigated the role of the heavy ions (mainly He(+) and O(+)) in generation and propagation of electromagnetic ion cyclotron waves and their contribution to the heating of magnetospheric electrons and ions. The more recent studies have also shown that the heavy ions can greatly contribute to a generation of lower hybrid waves, ring current precipitation phenomena, and the overall energy redistribution in the inner magnetosphere. Using newly developed 2.5-dimensional particle-in-cell simulations, we study the energization and nonlinear coupling of different plasma waves in the presence of the heavy ions. We have shown that the high frequency wave modes critically depend on the heavy ion density and irrespective of the driven wave modes, both the light and heavy ions undergo significant transverse acceleration. But for the large heavy-ion densities, even the electrons are significantly accelerated in the parallel direction by the waves below the LH frequency.

Khazanov, G. V.↗

Assessment of the Impacts of ACLS on the ISS Life Support System using Dynamic Simulations in V-HAB

The Advanced Closed Loop System (ACLS) is currently under development by Airbus Defense and Space and is slated for launch to the International Space Station (ISS) in 2017. The addition of new hardware into an already complex system such as the ISS life support system (LSS) always poses operational risks. It is therefore important to understand the impacts ACLS will have on the existing systems to ensure smooth operations for the ISS. This analysis can be done by using dynamic computer simulations and one possible tool for such a simulation is Virtual Habitat (V-HAB). Based on Matlab (Registered Trademark) V-HAB has been under development at the Institute of Astronautics of the Technical University Munich (TUM) since 2006 and in the past has been successfully used to simulate the ISS life support systems. The existing V-HAB ISS simulation model treated the interior volume of the space station as one large ideally-stirred container. This model was improved to allow the calculation of the atmospheric composition inside the individual modules of the ISS by splitting it into ten distinct volumes. The virtual volumes are connected by a simulation of the inter-module ventilation flows. This allows for a combined simulation of the LSS hardware and the atmospheric composition aboard the ISS. A dynamic model of ACLS is added to the ISS simulation and different operating modes for both ACLS and the existing ISS life support systems are studied to determine the impacts of ACLS on the rest of the system. The results suggest that the US, Russian and ACLS CO2 systems can operate at the same time without impeding each other. Furthermore, based on the results of this analysis, the US and ACLS Sabatier systems can be operated in parallel as well to achieve the highest possible CO2 recycling together with a low CO2 concentration.

Puetz, Daniel↗

On the Origin of the Photospheric Magnetic Field

This article presents results that challenge the paradigms that (1) the convection zone is the source of the radial magnetic field in the photosphere and (2) that coronal currents are neutralized from the perspective of the photosphere. We demonstrate, using a new analysis tool applied to simulations and observations, that bare or partially dressed current channels are supported by the solar corona and that fingerprints of these coronal current systems can be detected in the photosphere. These coronal current channels can be a significant source of the radial component of the magnetic field in the photosphere. The roots of these coronal current channels in the photosphere are the source of the magnetic field component parallel to the polarity inversion line in active region NOAA 12673. These analyses and observations transform our theoretical understanding of coronal evolution and argue for a reexamination of the present paradigm in which the convection zone is the sole source of the photospheric magnetic field.

magnetic fields↗

Thermal ions at Saturn - Plasma parameters and implications

Plasma data from the PLS experiment on Voyagers 1 and 2 are analyzed at the Saturn encounters. Measured ion currents are simulated using the full response function of the PLS instrument assuming plasma distributions are isotropic Maxwellians. The plasma velocity, densities, and temperatures which best simulate the observations are determined. Saturn's magnetosphere is found to subcorotate outside of L = 5.5 for Voyager 1 and L = 7 for Voyager 2. Large radial and vertical flows are observed, especially in the outer magnetosphere and on the nightside. Ion temperatures increase with L to a maximum at L = 10, outside of which they vary rapidly. A temperature anisotropy with T-perpendicular greater than T-parallel is necessary to account for the observed temperatures and densities. The ions present are H(+) and O(+); there is no evidence of heavier ions or higher charge states. Calculations of N(L-squared) show an O(+) peak at L = 10 for both encounters, indicating the major O(+) source is at L = 10.

Richardson, J. D.↗

Waves associated with quasi-parallel shocks: Generation, mode conversion and implications

Ions that are energized at quasi-parallel collisionless shocks and move back upstream generate low-frequency waves, largely on the fast/magnetosonic branch. At sufficient Mach number, the waves are convected back into the shock, lead to shock re-formation, and are mode-converted into downstream (magnetosheath) Alfvenic turbulence. Other waves are generated more locally at the interface of the incoming solar wind and the partially thermalized plasma. This paper reviews how recent simulation studies of collisionless shocks in conjunction with linear kinetic theory and proper wave diagnostics have aided in our understanding of the upstream and magnetosheath waves.

Krauss-Varban, D.↗

A Probabilistic Framework for the Validation and Certification of Computer Simulations

The paper presents a methodology for quantifying, propagating, and managing the uncertainty in the data required to initialize computer simulations of complex phenomena. The purpose of the methodology is to permit the quantitative assessment of a certification level to be associated with the predictions from the simulations, as well as the design of a data acquisition strategy to achieve a target level of certification. The value of a methodology that can address the above issues is obvious, specially in light of the trend in the availability of computational resources, as well as the trend in sensor technology. These two trends make it possible to probe physical phenomena both with physical sensors, as well as with complex models, at previously inconceivable levels. With these new abilities arises the need to develop the knowledge to integrate the information from sensors and computer simulations. This is achieved in the present work by tracing both activities back to a level of abstraction that highlights their commonalities, thus allowing them to be manipulated in a mathematically consistent fashion. In particular, the mathematical theory underlying computer simulations has long been associated with partial differential equations and functional analysis concepts such as Hilbert spares and orthogonal projections. By relying on a probabilistic framework for the modeling of data, a Hilbert space framework emerges that permits the modeling of coefficients in the governing equations as random variables, or equivalently, as elements in a Hilbert space. This permits the development of an approximation theory for probabilistic problems that parallels that of deterministic approximation theory. According to this formalism, the solution of the problem is identified by its projection on a basis in the Hilbert space of random variables, as opposed to more traditional techniques where the solution is approximated by its first or second-order statistics. The present representation, in addition to capturing significantly more information than the traditional approach, facilitates the linkage between different interacting stochastic systems as is typically observed in real-life situations.

Ghanem, Roger↗

Pulsar Emission Geometry and Accelerating Field Strength

The high-quality Fermi LAT observations of gamma-ray pulsars have opened a new window to understanding the generation mechanisms of high-energy emission from these systems, The high statistics allow for careful modeling of the light curve features as well as for phase resolved spectral modeling. We modeled the LAT light curves of the Vela and CTA I pulsars with simulated high-energy light curves generated from geometrical representations of the outer gap and slot gap emission models. within the vacuum retarded dipole and force-free fields. A Markov Chain Monte Carlo maximum likelihood method was used to explore the phase space of the magnetic inclination angle, viewing angle. maximum emission radius, and gap width. We also used the measured spectral cutoff energies to estimate the accelerating parallel electric field dependence on radius. under the assumptions that the high-energy emission is dominated by curvature radiation and the geometry (radius of emission and minimum radius of curvature of the magnetic field lines) is determined by the best fitting light curves for each model. We find that light curves from the vacuum field more closely match the observed light curves and multiwavelength constraints, and that the calculated parallel electric field can place additional constraints on the emission geometry

DeCesar, Megan E.↗

Experiment Description and Results for Arrival Operations Using Interval Management with Spacing to Parallel Dependent Runways (IMSPiDR)

The predicted increase in the number of commercial aircraft operations creates a need for improved operational efficiency. Two areas believed to offer increases in aircraft efficiency are optimized profile descents and dependent parallel runway operations. Using Flight deck Interval Management (FIM) software and procedures during these operations, flight crews can achieve by the runway threshold an interval assigned by air traffic control (ATC) behind the preceding aircraft that maximizes runway throughput while minimizing additional fuel consumption and pilot workload. This document describes an experiment where 24 pilots flew arrivals into the Dallas Fort-Worth terminal environment using one of three simulators at NASA?s Langley Research Center. Results indicate that pilots delivered their aircraft to the runway threshold within +/- 3.5 seconds of their assigned time interval, and reported low workload levels. In general, pilots found the FIM concept, procedures, speeds, and interface acceptable. Analysis of the time error and FIM speed changes as a function of arrival stream position suggest the spacing algorithm generates stable behavior while in the presence of continuous (wind) or impulse (offset) error. Concerns reported included multiple speed changes within a short time period, and an airspeed increase followed shortly by an airspeed decrease.

Baxley, Brian T.↗

Method and apparatus for measuring minority carrier lifetimes and bulk diffusion length in P-N junction solar cells

Carrier lifetimes and bulk diffusion length are qualitatively measured as a means for qualification of a P-N junction photovoltaic solar cell. High frequency (blue) monochromatic light pulses and low-frequency (red) monochromatic light pulses were alternately applied to the cell while it was irradiated by light from a solar simulator, and synchronously displaying the derivative of the output voltage of the cell on an oscilloscope. The output voltage is a measure of the lifetimes of the minority carriers (holes) in the diffused N layer and majority carriers (electrons) in the bulk P material, and of the diffusion length of the bulk silicon. By connecting a reference cell in this manner with a test cell to be tested in reverse parallel, the display of a test cell that matches the reference cell will be a substantially zero output.

Vonroos, O. H.↗

Aerothermal loads analysis for high speed flow over a quilted surface configuration

Attention is given to hypersonic laminar flow over a quilted surface configuration that simulates an array of Space Shuttle Thermal Protection System panels bowed in a spherical shape as a result of thermal gradient through the panel thickness. Pressure and heating loads to the surface are determined. The flow field over the configuration was mathematically modeled by means of time-dependent, three-dimensional conservation of mass, momentum, and energy equations. A boundary mapping technique was then used to obtain a rectangular, parallel piped computational domain, and an explicit MacCormack (1972) explicit time-split predictor corrector finite difference algorithm was used to obtain steady state solutions. Total integrated heating loads vary linearly with bowed height when this value does not exceed the local boundary layer thickness.

Olsen, G. C.↗