Search NASA⌕ Search

SEARCH · Search NASA

Results for “Parallel in Time”

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 901 records · Page 50

Space Technology 5 Multipoint Observations of Temporal and Spatial Variability of Field-Aligned Currents

Space Technology 5 (ST5) is a constellation mission consisting of three microsatellites. It provides the first multipoint magnetic field measurements in low Earth orbit, which enables us to separate spatial and temporal variations. In this paper, we present a study of the temporal variability of field-aligned currents using the ST5 data. We examine the field-aligned current observations during and after a geomagnetic storm and compare the magnetic field profiles at the three spacecraft. The multipoint data demonstrate that mesoscale current structures, commonly embedded within large-scale current sheets, are very dynamic with highly variable current density and/or polarity in approx.10 min time scales. On the other hand, the data also show that the time scales for the currents to be relatively stable are approx.1 min for mesoscale currents and approx.10 min for large-scale currents. These temporal features are very likely associated with dynamic variations of their charge carriers (mainly electrons) as they respond to the variations of the parallel electric field in auroral acceleration region. The characteristic time scales for the temporal variability of mesoscale field-aligned currents are found to be consistent with those of auroral parallel electric field.

Le, G.↗

Parallel Implementation of Power System Dynamic Simulation

Dynamic simulation of power system transient stability is important for planning, monitoring, operation, and control of electrical power systems. However, modeling the system dynamics and network involves the computationally intensive time-domain solution of numerous differential and algebraic equations (DAE). This results in a transient stability implementation that may not maintain the real-time constraints of an online security assessment. This paper presents a parallel implementation of the dynamic simulation on a high-performance computing (HPC) platform using parallel simulation algorithms and computation architectures. It enables the simulation to run even faster than real time, enabling the “look-ahead” capability of upcoming stability problems in the power grid.

dynamic simulation↗

Ultrahigh-order Maxwell solver with extreme scalability for electromagnetic PIC simulations of plasmas

The advent of massively parallel supercomputers, with their distributed-memory technology using many processing units, has favored the development of highly-scalable local low-order solvers at the expense of harder-to-scale global very high-order spectral methods. Indeed, FFT-based methods, which were very popular on shared memory computers, have been largely replaced by finite-difference (FD) methods for the solution of many problems, including plasmas simulations with electromagnetic Particle-In-Cell methods. For some problems, such as the modeling of so-called “plasma mirrors” for the generation of high-energy particles and ultra-short radiations, we have shown that the inaccuracies of standard FD-based PIC methods prevent the modeling on present supercomputers at sufficient accuracy. Here, we demonstrate here that a new method, based on the use of local FFTs, enables ultrahigh-order accuracy with unprecedented scalability, and thus for the first time the accurate modeling of plasma mirrors in 3D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A laser sounder for U.S. Navy helicopters

The design and operating principles of the laser sounder developed for use in airborne coastal-zone bathymetric surveys with a U.S. Navy P-3 aircraft are described and illustrated with diagrams. The sounder components are listed and characterized, including the Nd:YAG transmitter (operating at 532 nm for bottom ranging and 1.064 microns for surface ranging), the scanning transceiver, the 10 x 6-inch-aperture controlled-FOV receiver/digitizer, the constant-fraction discriminator, the time-to-digital converter, the navigation and data-recording subsystems, and the parallel distributed processing computer (comprising a data collection and control system and a real-time processing system). Consideration is also given to the phase-I and phase-II data-reduction software being developed to process the approximately 228 million soundings to be obtained annually. The sounder can be used day or night to sound clear ocean water up to 20 m deep.

Harris, Michael M.↗

Computational Methods for HSCT-Inlet Controls/CFD Interdisciplinary Research

A program aimed at facilitating the use of computational fluid dynamics (CFD) simulations by the controls discipline is presented. The objective is to reduce the development time and cost for propulsion system controls by using CFD simulations to obtain high-fidelity system models for control design and as numerical test beds for control system testing and validation. An interdisciplinary team has been formed to develop analytical and computational tools in three discipline areas: controls, CFD, and computational technology. The controls effort has focused on specifying requirements for an interface between the controls specialist and CFD simulations and a new method for extracting linear, reduced-order control models from CFD simulations. Existing CFD codes are being modified to permit time accurate execution and provide realistic boundary conditions for controls studies. Parallel processing and distributed computing techniques, along with existing system integration software, are being used to reduce CFD execution times and to support the development of an integrated analysis/design system. This paper describes: the initial application for the technology being developed, the high speed civil transport (HSCT) inlet control problem; activities being pursued in each discipline area; and a prototype analysis/design system in place for interactive operation and visualization of a time-accurate HSCT-inlet simulation.

Cole, Gary L.↗

Automated Generation of Message-Passing Programs: An Evaluation Using CAPTools

Scientists at NASA Ames Research Center have been developing computational aeroscience applications on highly parallel architectures over the past ten years. During that same time period, a steady transition of hardware and system software also occurred, forcing us to expend great efforts into migrating and re-coding our applications. As applications and machine architectures become increasingly complex, the cost and time required for this process will become prohibitive. In this paper, we present the first set of results in our evaluation of interactive parallelization tools. In particular, we evaluate CAPTool's ability to parallelize computational aeroscience applications. CAPTools was tested on serial versions of the NAS Parallel Benchmarks and ARC3D, a computational fluid dynamics application, on two platforms: the SGI Origin 2000 and the Cray T3E. This evaluation includes performance, amount of user interaction required, limitations and portability. Based on these results, a discussion on the feasibility of computer aided parallelization of aerospace applications is presented along with suggestions for future work.

Hribar, Michelle R.↗

Automated Generation of Message-Passing Programs: An Evaluation Using CAPTools

Scientists at NASA Ames Research Center have been developing computational aeroscience applications on highly parallel architectures over the past ten years. During that same time period, a steady transition of hardware and system software also occurred, forcing us to expend great efforts into migrating and re-coding our applications. As applications and machine architectures become increasingly complex, the cost and time required for this process will become prohibitive. In this paper, we present the first set of results in our evaluation of interactive parallelization tools. In particular, we evaluate CAPTool's ability to parallelize computational aeroscience applications. CAPTools was tested on serial versions of the NAS Parallel Benchmarks and ARC3D, a computational fluid dynamics application, on two platforms: the SGI Origin 2000 and the Cray T3E. This evaluation includes performance, amount of user interaction required, limitations and portability. Based on these results, a discussion on the feasibility of computer aided parallelization of aerospace applications is presented along with suggestions for future work.

Hribar, Michelle R.↗

An efficient parallel implementation of explicit multirate Runge–Kutta schemes for discontinuous Galerkin computations

Although explicit time integration schemes require small computational efforts per time step, their efficiency is severely restricted by their stability limits. Indeed, the multi-scale nature of some physical processes combined with highly unstructured meshes can lead some elements to impose a severely small stable time step for a global problem. Multirate methods offer a way to increase the global efficiency by gathering grid cells in appropriate groups under local stability conditions. These methods are well suited to the discontinuous Galerkin framework. The parallelization of the multirate strategy is challenging because grid cells have different workloads. The computational cost is different for each sub-time step depending on the elements involved and a classical partitioning strategy is not adequate any more. In this paper, we propose a solution that makes use of multi-constraint mesh partitioning. It tends to minimize the inter-processor communications, while ensuring that the workload is almost equally shared by every computer core at every stage of the algorithm. Particular attention is given to the simplicity of the parallel multirate algorithm while minimizing computational and communication overheads. Our implementation makes use of the MeTiS library for mesh partitioning and the Message Passing Interface for inter-processor communication. Performance analyses for two and three-dimensional practical applications confirm that multirate methods preserve important computational advantages of explicit methods up to a significant number of processors.

97 MATHEMATICS AND COMPUTING↗

LEWICE droplet trajectory calculations on a parallel computer

A parallel computer implementation (128 processors) of LEWICE, a NASA Lewis code used to predict the time-dependent ice accretion process for two-dimensional aerodynamic bodies of simple geometries, is described. Two-dimensional parallel droplet trajectory calculations are performed to demonstrate the potential benefits of applying parallel processing to ice accretion analysis. Parallel performance is evaluated as a function of the number of trajectories and the number of processors. For comparison, similar trajectory calculations are performed on single-processor Cray computers, and the best parallel results are found to be 33 and 23 times faster, respectively, than those of the Cray XMP and YMP.

Caruso, Steven C.↗

Stochastic density functional theory

Linear-scaling implementations of density functional theory (DFT) reach their intended efficiency regime only when applied to systems having a physical size larger than the range of their Kohn–Sham density matrix (DM). This causes a problem since many types of large systems of interest have a rather broad DM range and are therefore not amenable to analysis using DFT methods. For this reason, the recently proposed stochastic DFT (sDFT), avoiding exhaustive DM evaluations, is emerging as an attractive alternative linear-scaling approach. This review develops a general formulation of sDFT in terms of a (non)orthogonal basis representation and offers an analysis of the statistical errors (SEs) involved in the calculation. Using a new Gaussian-type basis-set implementation of sDFT, applied to water clusters and silicon nanocrystals, it demonstrates and explains how the standard deviation and the bias depend on the sampling rate and the system size in various types of calculations. We also develop a basis-set embedded-fragments theory, demonstrating its utility for reducing the SEs for energy, density of states and nuclear force calculations. Finally, we discuss the algorithmic complexity of sDFT, showing it has CPU wall-time linear-scaling. Lastly, the method parallelizes well over distributed processors with good scalability and therefore may find use in the upcoming exascale computing architectures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A generalized parallel replica dynamics

Metastability is a common obstacle to performing long molecular dynamics simulations. Many numerical methods have been proposed to overcome it. One method is parallel replica dynamics, which relies on the rapid convergence of the underlying stochastic process to a quasi-stationary distribution. Two requirements for applying parallel replica dynamics are knowledge of the time scale on which the process converges to the quasi-stationary distribution and a mechanism for generating samples from this distribution. By combining a Fleming–Viot particle system with convergence diagnostics to simultaneously identify when the process converges while also generating samples, we can address both points. This variation on the algorithm is illustrated with various numerical examples, including those with entropic barriers and the 2D Lennard-Jones cluster of seven atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ultra-Fast Vertical Ordering of Lamellar Block Copolymer Films on Unmodified Substrates

To utilize the full potential of block copolymer (BCP) thin films for use in technological devices ranging from ion conduction membranes, transistors to nanowire array antennas, rapid self- assembly of lamellar block copolymers (l-BCPs) with vertically oriented lamellar domains on a variety of unmodified substrates is needed. l-BCPs have an inherently larger interfacial area for transport compared to their cylindrical counterpart. Our observations demonstrate that the as-cast weakly ordered vertically oriented state of l-BCP films of polystyrene-block-poly(methyl methacrylate) (PS-b-PMMA) from directional evaporation of select solvents, act as “seed templates” for their ultra-fast evolution (~30 s) into well-ordered vertically oriented nanostructures, using a thermal gradient-based Cold Zone Annealing (CZA) technique. Furthermore, vertical lamellae are obtained on unmodified substrates, Quartz and Kapton, and the kinetics of l-BCP ordering is much faster by CZA as compared to the isotropic oven annealing. The rapid ordering kinetics of vertical l-BCPs is tested and found applicable to different molecular weights and film thicknesses ranging from 20 nm to 480 nm, which ultimately flip over to their equilibrium parallel morphology at upper limits of annealing times. This rapid ordering strategy for vertical orientation of l-BCPs using roll-to-roll compatible CZA would be highly relevant for fundamental studies of interfacial transport as well as for industrial applications from membranes to nanowires.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Quantum Dot Nanobiosensor for Rapid Detection of Botulinum Neurotoxin Serotype E

Botulinum neurotoxins (BoNTs) are potent toxins produced by Clostridium bacteria that are responsible for the illness botulism and are listed as bioterrorism agents. BoNT serotype E (BoNT/E) is one of four BoNT serotypes that cause human botulism and is the second most frequent cause of foodborne botulism. Rapid detection and discrimination of BoNT serotypes implicated in human disease are critical for ensuring timely treatment of patients and identifying sources of toxins, but there have been few reported detection methods for BoNT/E and even fewer methods usable for BoNT serotyping. We report a nanobiosensor based on Forster resonance energy transfer (FRET) between semiconductor nanocrystals (quantum dots, QDs) and dark quencher-labeled peptide probes to detect biologically active BoNT/E in aqueous media. The peptide probes contain a specific cleavage site for active BoNT/E. QD photoluminescence, which changes intensity due to FRET when the peptide probe is cleaved, was used to indicate toxin presence and quantity. The detection of a BoNT/E light chain (LcE) and holotoxin was observed within 3 h. The limits of detection were 0.02 and 2 ng/mL for LcE and holotoxin, respectively. The nanobiosensor shows good specificity toward the target in tests with nontarget BoNT serotypes. The high sensitivity, simple operation, short detection time, and ability to be used in parallel with probes developed for other BoNT serotypes indicate that the nanobiosensor will be useful for rapid BoNT/E detection and serotype discrimination in food analysis.

botulinum neurotoxin E↗

Understanding resistances in capacitive deionization devices

The desalination characteristics of capacitive deionization (CDI) are significantly impacted by the cell resistance. Here we apply electrochemical impedance spectroscopy (EIS) to analyze the resistive properties of a CDI device, which include EIS measurements on a complete assembled cell (in two wire mode) and EIS measurements with additional contacts to decouple electrode/electrolyte and contact impedances (four wire mode). These measurements shed light on the interface between current collector and electrode as well as the internal capacitive and resistive elements of the cell. With two-wire and four-wire EIS measurements, we find that the often observed high-frequency arc(s) in the impedance spectra can be due to contact effects and/or an internal ionic-transfer resistance that is locally in parallel with macropore wall capacitance. In some cases, the two associated arcs in the Nyquist plot are overlapping and present simultaneously, making their identification in the EIS spectra difficult unless both two-wire and four-wire measurements are made. Despite causing a large high-frequency arc in the Nyquist plot, we find that the apparent internal ionic-transfer resistance is negligible compared to the resistances of the electrode, separator, and external electronic components. Here, an investigation of the contact impedance between the electrode and current collector reveals a parallel RC arc that can grow over time if an improper material is used (e.g., titanium). In our system, the electrode ionic resistance dominates the total Ohmic power dissipation, which has implications for scale-up of CDI systems and future improvements of cell design.

36 MATERIALS SCIENCE↗

A brightness exceeding simulated Langmuir limit

When an excitation of the first lens determines a beam is parallel beam, a brightness that is 100 times higher than Langmuir limit is measured experimentally, where Langmuir limits are estimated using a simulated axial cathode current density which is simulated based on a measured emission current. The measured brightness is comparable to Langmuir limit, when the lens excitation is such that an image position is slightly shorter than a lens position. Previously measured values of brightness for cathode apical radii of curvature 20, 60, 120, 240, and 480 μm were 8.7, 5.3, 3.3, 2.4, and 3.9 times higher than their corresponding Langmuir limits, respectively, in this experiment, the lens excitation was such that the lens and the image positions were 180 mm and 400 mm, respectively. From these measured brightness for three different lens excitation conditions, it is concluded that the brightness depends on the first lens excitation. For the electron gun operated in a space charge limited condition, some of the electrons emitted from the cathode are returned to the cathode without having crossed a virtual cathode. Therefore, method that assumes a Langmuir limit defining method using a Maxwellian distribution of electron velocities may need to be revised. For the condition in which the values of the exceeding the Langmuir limit are measured, the simulated trajectories of electrons that are emitted from the cathode do not cross the optical axis at the crossover, thus the law of sines may not be valid for high brightness electron beam systems.

43 PARTICLE ACCELERATORS↗

Design parameters and commissioning of vertical inserts used for testing the XFEL superconducting cavities

The European XFEL is a new research facility currently under construction at DESY in the Hamburg area in Germany. From 2015 on, it will generate extremely intense X-ray flashes that will be used by researchers from all over the world. The superconducting XFEL linear accelerator consists of 100 accelerator modules with more than 800 RF-cavities inside. The accelerator modules, superconducting magnets and cavities will be tested in the accelerator module test facility (AMTF). This paper gives an overview of the design parameters and the commissioning of the vertical insert, used in two cryostats (XATC) of the AMTF-hall. The Insert serves as a holder for 4 nine-cell cavities. This gives the possibility to cool down 4 cavities to 2K in parallel and, consequently, to reduce the testing time. The following RF measurement, selected as quality check, will be done separately for each cavity. Afterwards the cavities will be warmed up again and will be sent to the accelerator module assembly.

43 PARTICLE ACCELERATORS↗

Nonlinear simulations of GAEs in NSTX-U

A set of nonlinear simulations has been performed in order to study the nonlinear evolution of unstable global Alfvén eigenmodes in the National Spherical Torus Experiment-Upgrade (NSTX-U). Results of the single toroidal mode number, n, simulations are compared with a full nonlinear simulation (all toroidal harmonics included). In single-n simulations, the conservation of two integrals of motion of a particle in a cyclotron resonance with a monochromatic wave is demonstrated, resulting in a one-dimensional evolution of the particle distribution in (E,μ,pϕ) phase-space. Nonlinear simulations (both single-n and full nonlinear) show a significant redistribution of the resonant fast ions, especially in the pitch parameter. Thus, the changes in the resonant particle's parallel and perpendicular energies can be several times larger than the total particle energy change, with only a small fraction transferring into the excitation of the mode itself. This implies that even a relatively small amplitude mode can significantly modify the beam distribution in the resonant region. For the NSTX-U case considered, the single-n simulation results are close to full nonlinear simulation only for the most unstable mode, in which case the saturation amplitudes and changes in the fast ion distribution are comparable. In contrast, peak amplitudes of subdominant modes in all-n simulations are smaller by a factor of 3–10 compared to single-n runs due to the flattening of the beam ion distribution by the fastest growing mode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗