Search NASA⌕ Search

SEARCH · Search NASA

Results for “consistent discretization”

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 379 records · Page 21

Localized Evaluation for Constructing Discrete Vector Fields

Topological abstractions offer a method to summarize the behavior of vector fields, but computing them robustly can be challenging due to numerical precision issues. One alternative is to represent the vector field using a discrete approach, which constructs a collection of pairs of simplices in the input mesh that satisfies criteria introduced by Forman's discrete Morse theory. While numerous approaches exist to compute pairs in the restricted case of the gradient of a scalar field, state-of-the-art algorithms for the general case of vector fields require expensive optimization procedures. This paper introduces a fast, novel approach for pairing simplices of two-dimensional, triangulated vector fields that do not vary in time. The key insight of our approach is that we can employ a local evaluation, inspired by the approach used to construct a discrete gradient field, where every simplex in a mesh is considered by no more than one of its vertices. Specifically, we observe that for any edge in the input mesh, we can uniquely assign an outward direction of flow. We can further expand this consistent notion of outward flow at each vertex, which corresponds to the concept of a downhill flow in the case of scalar fields. Working with outward flow enables a linear-time algorithm that processes the (outward) neighborhoods of each vertex one-by-one, similar to the approach used for scalar fields. Here, we couple our approach to constructing discrete vector fields with a method to extract, simplify, and visualize topological features. Empirical results on analytic and simulation data demonstrate drastic improvements in running time, produce features similar to the current state-of-the-art, and show the application of simplification to large, complex flows.

97 MATHEMATICS AND COMPUTING↗

Multigroup Thermal Radiation Transport with Tensor Trains

We investigate the application of tensor-train (TT) algorithms to multigroup thermal radiation transport (i.e., photon radiation transport). The TT framework enables simulations at discretizations that might otherwise be computationally infeasible on conventional hardware. We show that solutions to certain multigroup problems possess an intrinsic low-rank structure, which the TT representation leverages effectively. This enables us to solve problems where the discretized solution size exceeds a trillion parameters on a single node. The solver is evaluated on a range of test problems with varying levels of complexity, consistently achieving compression factors greater than 100× and speedups exceeding 2×. We also investigate alternative TT topologies by analyzing the low-rank structure of the merged spatio-spectral core to assess the potential for greater compression. This analysis suggests that compression gains could increase by factors as large as 7. Our results indicate that the low-rank structure of the merged spatio-spectral core captures the spatio-spectral complexity of the solution, largely driven by the opacity structure of the medium. Beyond identifying opportunities for improved compression, this analysis highlights the types of errors that may arise in angle-integrated quantities when exploiting this low-rank structure.

79 ASTRONOMY AND ASTROPHYSICS↗

A 100,000 pulse height analyzer for use in High Energy Cosmic Ray Experiments (HECRE) on high altitude balloon flights

A description is given of a wide dynamic range pulse height analyzer system developed for use on High Energy Cosmic Ray Experiment (HECRE) Balloon Flights. A wide dynamic range of 100,000 is obtained by extending the range of a basic 1024 channel analyzer through the use of multiple ranges and range selection. The system described here contains four 100,000 pulse height analyzers. Each 100,000 pulse height analyzer consists of a group of cordwood welded modules mounted and interconnected on a printed circuit card. Four of these card assemblies, the required clock drive circuitry (discrete components mounted and interconnected on a separate card) and three input-output connectors are interconnected and mounted on the system board.

Cancro, C. A.↗

Penetration of solar irradiances through the atmosphere and plant canopies.

The equation of radiative transfer is applied to an analysis of solar irradiances penetrating into a plant canopy covered by a turbid atmosphere. The method of discrete coordinates is applied to vertically inhomogeneous atmospheres and plant canopies. It is shown that four-point quadrature yields results with an accuracy that is consistent with irradiance measurements. It is believed that the presented computational scheme may have considerable agricultural applications.

Weinman, J. A.↗

Experimental investigation of orthotropic panel flutter at arbitrary yaw angles, and comparison with theory

Flutter characteristics for yaw angles between 15 deg and 90 deg were determined experimentally for two types of corrugation-stiffened panels: those with weak twisting stiffness and those with strong twisting stiffness. By mounting the panels on a remotely controlled turntable, good definition of the flutter boundaries was obtained by rotating the panels into and out of flutter. Flutter tests were conducted at M = 2 and M = 1.6. Before testing, vibration tests and analyses were also performed. The experimental flutter data is compared with flutter theory for orthotropic panels utilizing quasi-steady aerodynamics. Five different corrugated panels were tested consisting of one single skin panel having a length-to-width ratio of 5 on clamped supports and four different square double skin panels on discrete flexible supports. The investigation indicated that flutter speed for corrugated panels is highly dependent on yaw angle. Reasonable flutter correlation between analysis and test was obtained for moderate yaw angles, but extreme sensitivity to structural parameters made the correlation at large yaw angles uncertain.

Shyprykevich, P.↗

Acoustic positioning for space processing experiments

An acoustic positioning system is described that is adaptable to a range of processing chambers and furnace systems. Operation at temperatures exceeding 1000 C is demonstrated in experiments involving the levitation of liquid and solid glass materials up to several ounces in weight. The system consists of a single source of sound that is beamed at a reflecting surface placed a distance away. Stable levitation is achieved at a succession of discrete energy minima contained throughout the volume between the reflector and the sound source. Several specimens can be handled at one time. Metal discs up to 3 inches in diameter can be levitated, solid spheres of dense material up to 0.75 inches diameter, and liquids can be freely suspended in l-g in the form of near-spherical droplets up to 0.25 inch diameter, or flattened liquid discs up to 0.6 inches diameter. Larger specimens may be handled by increasing the size of the sound source or by reducing the sound frequency.

Whymark, R. R.↗

An investigation of coronal active region loop structures using AS&E rocket X-ray images

Simultaneous high spatial resolution observations at 6 cm in soft X-rays, in photospheric magnetograms, and in optical filtergrams were used to compare the most intense sources of centimetric emission in two active regions to coronal loops, sunspots, chromospheric structures, and photospheric magnetic fields. Results show that the majority of the bright microwave components are not associated with sunspots or X-ray emission. A nonthermal mechanism appears necessary to explain the brightest microwave components, discrete regions of continuous particle acceleration may be common in active regions. Studies of the plasma parameters of selected loops imply that the radio emission is consistent with gyro-resonance absorption at the third and fourth harmonic, at least from part of each loop. Results are presented for: (1) X-ray and microwave observations of active regions; (2) comparison of coronal holes observed in soft X-rays and Hel 10830 A spectrosheliograms; and (3) the reappearance of polar coronal holes and the evolution of the solar magnetic field.

Webb, D. F.↗

Exogenic and endogenic albedo and color patterns on Europa

New global and high-resolution multispectral mosaics of Europa have been produced from the Voyager imaging data. Photometric normalizations are based on multiple-image techniques that explicitly account for intrinsic albedo variations through pixel-by-pixel solutions. The exogenic color and albedo pattern on Europa is described by a second-order function of the cosine of the angular distance from the apex of orbital motion. On the basis of this second-order function and of color trends that are different on the leading and trailing hemispheres, the exogenic pattern is interpreted as being due to equilibrium between two dominant processes: (1) impact gardening and (2) magnetospheric interactions, including sulfur-ion implantation and sputtering redistribution. Removal of the model exogenic pattern in the mosaics reveals the endogenic variations, consisting of only two major units: darker (redder) and bright materials. Therefore Europa's visual spectral reflectivity is simple, having one continuous exogenic pattern and two discrete endogenic units.

Mcewen, A. S.↗

Properties of the X-ray emitting gas in early-type galaxies

The properties of the X-ray emitting gas in a sample of 81 E and S0 galaxies observed with the Einstein Observatory are studied. Measured fluxes for 55 of the galaxies and upper limits for 26 of them are reported. An attempt is made to use consistent optical parameters for the galaxies, including a correction to the velocities for the Virgocentric flow. The sample is then used to explore the contribution from discrete sources, the global physical properties of the hot gas, and the implications for heating by supernovae and gravity. Finally, the question of the presence of heavy halos is addressed.

Canizares, Claude R.↗

C IV and Si IV in IUE spectra of normal B8-A0 stars: UV identified Be/Ae stars

Archival IUE high dispersion spectra of 42 B6-A2 stars within 200 pc were surveyed. Five of the program stars show significant C IV and Si IV absorption. All of the stars with detected C IV have v sin i less than or = 190 km/sec. Sharp absorption cores are present in Si II lambda 1533 in 3 of the objects, indicating that these are previously unrecognized shell stars. Three of the stars have variable or asymmetric C IV profiles which are consistent with the C IV and Si IV being produced in stellar winds. One star has C IV in the form of a shortward-shifted discrete absorption component, similar to those observed in Be stars. The data are compared with similar data for Be and B shell stars.

Chiu, Weihsueh A.↗

Studies in astronomical time series analysis. III - Fourier transforms, autocorrelation functions, and cross-correlation functions of unevenly spaced data

This paper develops techniques to evaluate the discrete Fourier transform (DFT), the autocorrelation function (ACF), and the cross-correlation function (CCF) of time series which are not evenly sampled. The series may consist of quantized point data (e.g., yes/no processes such as photon arrival). The DFT, which can be inverted to recover the original data and the sampling, is used to compute correlation functions by means of a procedure which is effectively, but not explicitly, an interpolation. The CCF can be computed for two time series not even sampled at the same set of times. Techniques for removing the distortion of the correlation functions caused by the sampling, determining the value of a constant component to the data, and treating unequally weighted data are also discussed. FORTRAN code for the Fourier transform algorithm and numerical examples of the techniques are given.

Scargle, Jeffrey D.↗

A Hybrid Fiber/Solid-State Regenerative Amplifier with Tunable Pulse Widths for Satellite Laser Ranging

A fiber/solid-state hybrid seeded regenerative amplifier, capable of achieving high output energy with tunable pulse widths, has been developed for satellite laser ranging applications. The regenerative amplifier cavity uses a pair of Nd:YAG zigzag slabs oriented orthogonally to one another in order to make thermal lensing effects symmetrical and simplify optical correction schemes. The seed laser used is a fiber-coupled 1,064-nm narrowband (<0.02 nm) diode laser that is discretely driven in a new short-pulsed mode, enabling continuously tunable seed pulse widths in the 0.2-to-0.4-ns range. The amplifier gain unit consists of a pair of Brewster-cut 6-bounce zigzag Nd:YAG laser slabs, oriented 90deg relative to each other in the amplifier head. This arrangement creates a net-symmetrical thermal lens effect (an opposing singleaxis effect in each slab), and makes thermo-optical corrections simple by optimizing the curvature of the nearest cavity mirror. Each slab is pumped by a single 120-W, pulsed 808-nm laser diode array. In this configuration, the average pump beam distribution in the slabs had a 1-D Gaussian shape, which matches the estimated cavity mode size. A half-wave plate between the slabs reduces losses from Fresnel reflections due to the orthogonal slabs Brewster-cut end faces. Successful "temporal" seeding of the regenerative amplifier cavity results in a cavity Q-switch pulse envelope segmenting into shorter pulses, each having the width of the input seed, and having a uniform temporal separation corresponding to the cavity round-trip time of approx. =10 ns. The pulse energy is allowed to build on successive passes in the regenerative amplifier cavity until a maximum is reached, (when cavity gains and losses are equal), after which the pulse is electro- optically switched out on the next round trip The overall gain of the amplifier is approx. =82 dB (or a factor of 1.26 million). After directing the amplified output through a LBO frequency doubling crystal, approx. = 2.1 W of 532-nm output (>1 mJ) was measured. This corresponds to a nonlinear conversion efficiency of >60%. Furthermore, by pulse pumping this system, a single pulse per laser shot can be created for the SLR (satellite laser ranging) measurement, and this can be ejected into the instrument. This is operated at the precise frequency needed by the measurement, as opposed to commercial short-pulsed, mode-locked systems that need to operate in a continuous fashion, or CW (continuous wave), and create pulses at many MHz. Therefore, this design does not need to throw away or dump 99% of the laser energy to produce what is required; this system can be far smaller, more efficient, cheaper, and readily deployed in the field when packaged efficiently. Finally, by producing custom diode seed pulses electronically, two major advantages over commercial systems are realized: First, this pulse shape is customizable and not affected by the cavity length or gain of the amplifier cavity, and second, it can produce adjustable (selectable) pulse widths by simply adding multiple seed diodes and coupling each into commercial, low-cost fiber-optic combiners.

Coyle, Barry↗

A framework for discrete optimization of stellarator coils

Designing magnets for three-dimensional plasma confinement is a key task for advancing the stellarator as a fusion reactor concept. Stellarator magnets must produce an accurate field while leaving adequate room for other components and being reasonably simple to construct and assemble. In this paper, a framework for coil design and optimization is introduced that enables the attainment of sparse magnet solutions with arbitrary restrictions on where coils may be located. The solution space is formulated as a 'wireframe' consisting of a mesh of interconnected wire segments enclosing the plasma. Two methods are developed for optimizing the current distribution on a wireframe: Regularized Constrained Least Squares, which uses a linear least-squares approach to optimize the currents in each segment, and Greedy Stellarator Coil Optimization, a fully discrete procedure in which loops of current are added to the mesh one by one to achieve the desired magnetic field on the plasma boundary. Examples are presented of solutions obtainable with each method, some of which achieve high field accuracy while obeying spatial constraints that permit easy assembly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Photonuclear cross sections for 197 Au: An update on the gold standard

Cross sections for the 197 Au(γ, n) reaction are broadly used in nuclear physics as a standard for normalizing photonuclear reaction cross-section data at photon beam energies above approximately 8 MeV. In this paper, we report cross-section measurements for the 197 Au(γ, n) 196 Au g+m1 reaction at beam energies from 13 to 31 MeV. Our measurements provide the first cross-section data for this reaction at beam energies above 20 MeV, enabling the use of this reaction as a cross-section standard up to 30 MeV. Also, this work provides first cross-section measurements for the 197 Au(γ, n) 196 Au m2 reaction. In addition, we measured cross-section data for the 197 Au(γ, 3n) 194 Au reaction, which can be used as a cross-section standard above about 25 MeV. These measurements were performed using a new target activation method that is based on the angle-energy correlation of the laser Compton- scattered photon beams at the High Intensity Gamma-ray Source (HIγS). The technique enables measuring photonuclear reaction cross-sections at several discrete beam energies concurrently via a single irradiation on a stack of different targets. Measurements were carried out by irradiating a stack of concentric-ring targets consisting of Au, TiO 2 , Zn, Os, and Au (in order of the γ-ray beam direction). Our data for the 197 Au(γ,n) 196 Au g+m1 reaction in the energy range of 13 to 20 MeV are in good agreement with existing ones measured using monoenergetic γ-ray beams, but differ from data acquired using a bremsstrahlung γ-ray beam. Also, above 18 MeV, our data for the 197 Au(γ,n) 196 Au g+m1 and 197 Au(γ,n) 196 Au m2 reactions differ significantly from the most recent TENDL and JENDL evaluations, suggesting a need to update these data libraries. Furthermore, the TENDL evaluation and existing data are consistent with our data for the 197 Au(γ,3n) reaction, but differ significantly from the JENDL evaluation above 26 MeV.

190 ≤ A ≤ 219↗

Electron polar cap and the boundary of open geomagnetic field lines.

A total of 333 observations of the boundary of the polar access region for electrons (energies greater than 530 keV) provides a comprehensive map of the electron polar cap. The boundary of the electron polar cap, which should occur at the latitude separating open and closed field lines, is consistent with previously reported closed field line limits determined from trapped-particle data. The boundary, which is sharply defined, seems to occur at one of three discrete latitudes. Although the electron flux is generally uniform across the polar cap, a limited region of reduced access is observed about 10% of the time.

Evans, L. C.↗

Comet Halley near-nucleus phenomena in 1986

Computer enhanced high resolution groundbased images taken during Halley's most active phases in 1985-86 exhibit all of the coma structures identified in 1910 photographs. The structures result mostly from discrete dust jet activity producing spirals, arcs, and jagged-edged envelopes. Multiple nucleus and antisunward jet phenomena are also observed. The curvature of the jets is consistent with the location of the Earth with respect to the nucleus equatorial plane. The variation in jet activity during the Vega and Giotto flybys is easily seen on the groundbased images. Images of CN jets show that they usually do not correlate with dust jets. Wide field, medium resolution images of faint outer envelopes are required to augment the smaller field, high resolution images of near-nucleus jets in modeling the jets.

Larson, S.↗

Evaluation of a multi-point method for determining acoustic impedance

An investigation was conducted to explore potential improvements provided by a Multi-Point Method (MPM) over the Standing Wave Method (SWM) and Two-Microphone Method (TMM) for determining acoustic impedance. A wave propagation model was developed to model the standing wave pattern in an impedance tube. The acoustic impedance of a test specimen was calculated from a best fit of this standing wave pattern to pressure measurements obtained along the impedance tube centerline. Three measurement spacing distributions were examined: uniform, random, and selective. Calculated standing wave patterns match the point pressure measurement distributions with good agreement for a reflection factor magnitude range of 0.004 to 0.999. Comparisons of results using 2, 3, 6, and 18 measurement points showed that the most consistent results are obtained when using at least 6 evenly spaced pressure measurements per half-wavelength. Also, data were acquired with broadband noise added to the discrete frequency noise and impedances were calculated using the MPM and TMM algorithms. The results indicate that the MPM will be superior to the TMM in the presence of significant broadband noise levels associated with mean flow.

Jones, Michael G.↗

Surface tension and buoyancy-driven flow in a non-isothermal liquid bridge

The Navier-Stokes-Boussinesq equations governing the transport of momentum, mass and heat in a nonisothermal liquid bridge with a temperature-dependent surface tension are solved using a vorticity-stream-function formulation together with a nonorthogonal coordinate transformation. The equations are discretized using a pseudo-unsteady semi-implicit finite difference scheme and are solved by the ADI method. A Picard-type iteration is adopted which consists of inner and outer iterative processes. The outer iteration is used to update the shape of the free surface. Two schemes have been used for the outer iteration; both use the force balance normal to the free surface as the distinguished boundary condition. The first scheme involves successive approximation by the direct solution of the distinguished boundary condition. The second scheme uses the artificial force imbalance between the fluid pressure, viscous and capillary forces at the free surface which arises when the boundary condition for force balance normal to the surface is not satisfied. This artificial imbalance is then used to change the surface shape until the distinguished boundary condition is satisfied. These schemes have been used to examine a variety of model liquid bridge situations including purely thermocapillary-driven flow situations and mixed thermocapillary- and bouyancy-driven flow.

Zhang, Yiqiang↗