Results for “Numerical Relativity” — Page 50

DOE OSTI2021

PETSc Users Manual (Revision 3.15)

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

PETSc/TAO Users Manual: Revision 3.18

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Separation of CO 2 from Flue Gas and Potential for Geologic Sequestration

The objectives of this study were to review various methods reported in the literature for the separation and geologic sequestration of carbon dioxide and evaluate the potential of TVA fossil fuel-burning plant locations for onsite geologic sequestration of CO 2 from stack emissions. Several conventional and nonconventional technologies for the separation of CO 2 from flue gas, including absorption, adsorption, cryogenic distillation, membranes, hydrate formation and dissociation, and ammonia carbonation, have been reviewed in terms of separation mechanisms, flow diagrams, and costs. Most of the technologies that have been reviewed are still at the research and development stage. Critical information needed to assess and compare these technologies is still lacking. In addition, information on some of the technologies that have been tested at a pilot or industrial scale has not been fully disclosed in the open literature. Because of this lack of data, it is difficult to make a critical assessment of each of the separation technologies. Based on limited information, it was concluded that the most promising methods are membrane separation and the Mitsubishi process for chemical absorption. Both processes involve separating CO 2 at high temperature, minimizing the cost for cooling prior to separation. Physical and chemical geologic formations of CO 2 were also reviewed. It was concluded that due to the geologic time scale of CO 2 sequestration periods, relatively safe conditions, general proximity to CO 2 sources, and extensive knowledge of underground conditions, sequestration of CO 2 in underground aquifers and coal beds is a very promising method of mitigating greenhouse gas emissions. The cost is predicted to be relatively low and the suitable sites are numerous for this application, with many of these sites located close to the plants.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2023

PETSc/TAO Users Manual (Rev. 3.19)

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication. PETSc/TAO includes a large suite of parallel linear solvers, nonlinear solvers, time integrators, and opti mization that may be used in application codes written in Fortran, C, C++, and Python (via petsc4py; see Getting Started). PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than packages such as MATLAB or a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or lldb, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than “rolling them” yourself. For many tasks a package such as MATLAB is often the best tool; PETSc is not intended for the classes of problems for which effective MATLAB code can be written. There are several packages, built on PETSc, that may satisfy your needs without requiring directly using PETSc. We recommend reviewing these packages functionality before starting to code directly with PETSc. PETSc can be used to provide a “MPI parallel linear solver” in an otherwise sequential, or OpenMP parallel code. This approach cannot provide extremely large improvements in the application time by utilizing large numbers of MPI processes but can still improve the performance. Certainly all parts of a previously sequential code need not be parallelized but the matrix generation portion must be parallelized to expect true scalability to large numbers of MPI processes. See PCMPI for details on how to utilize the PETSc MPI linear solver server. Since PETSc is under continued development, small changes in usage and calling sequences of routines will occur. PETSc has been supported for twenty-five years; see mailing list information on our website for information on contacting support.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

PETSc/TAO Users Manual (Rev. 3.20)

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Late Holocene evolution of a coupled, mud-dominated delta plain–chenier plain system, coastal Louisiana, USA

Abstract. Major deltas and their adjacent coastal plains are commonly linked by means of coast-parallel fluxes of water, sediment, and nutrients. Observations of the evolution of these interlinked systems over centennial to millennial timescales are essential to understand the interaction between point sources of sediment discharge (i.e. deltaic distributaries) and adjacent coastal plains across large spatial (i.e. hundreds of kilometres) scales. This information is needed to constrain future generations of numerical models to predict coastal evolution in relation to climate change and other human activities. Here we examine the coastal plain (Chenier Plain, CP) adjacent to the Mississippi River delta, one of the world's largest deltas. We use a refined chronology based on 22 new optically stimulated luminescence and 22 new radiocarbon ages to test the hypothesis that cyclic Mississippi subdelta shifting has influenced the evolution of the adjacent CP. We show that over the past 3 kyr, accumulation rates in the CP were generally 0–1 Mt yr −1 . However, between 1.2 and 0.5 ka, when the Mississippi River shifted to a position more proximal to the CP, these rates increased to 2.9 ±1.1 Mt yr −1 or 0.5–1.5 % of the total sediment load of the Mississippi River. We conclude that CP evolution during the past 3 kyr was partly a direct consequence of shifting subdeltas, in addition to changing regional sediment sources and modest rates of relative sea-level (RSL) rise. The RSL history of the CP during this time period was constrained by new limiting data points from the base of overwash deposits associated with the cheniers. These findings have implications for Mississippi River sediment diversions that are currently being planned to restore portions of this vulnerable coast. Only if such diversions are located in the western portion of the Mississippi Delta plain could they potentially contribute to sustaining the CP shoreline. Our findings highlight the importance of a better understanding of mud-dominated shorelines that are often associated with major deltas, in light of the enormous investments in coastal management and restoration that will likely be made around the globe, now and especially later during this century.

Hijma, Marc P.↗

NASA NTRS1975

Optimum quantisers for a Gaussian input probability density and for the magnitude-error distortion measure

The parameters of non-uniform and uniform quantizers up to ten bits of quantization, optimum for a Gaussian input probability and for the magnitude-error distortion criterion are computed. Optimum quantizers must be understood as quantizers with minimum distortion. The numerical method used for the optimization converges relatively rapidly. The comparison between optimum non-uniform quantizers and optimum uniform quantizers is made.

Demaret, J. C.↗

NASA NTRS1979

Global magnetohydrodynamic simulation of the two-dimensional magnetosphere

The time-dependent magnetohydrodynamic interaction of the solar wind with a two-dimensional dipole magnetic field has been simulated using a novel Lagrangian particle type of MHD code that can treat local low density or vacuum regions without numerical instability. This enables one to simulate the time-dependent magnetic tail. When the solar wind field is southward, a magnetic field line topology consistent with Dungey's model emerges in steady state. The tail, however, is short, and the x-points are only slightly shifted from their vacuum locations, because of strong numerical resistivity. Different configurations resulting from different relative orientations of the solar wind magnetic field and dipole axis are also presented. While the magnetic field is relatively steady, the density and flow in the magnetosheath are turbulent, as are the bow shock and magnetopause; the Kelvin-Helmholtz instability may account for these phenomena. We also model a 'substorm' as the passage of a rotational discontinuity in the solar wind over the dipole. Both 90 and 180 deg shifts to a southward solar wind field cause a violent readjustment of the magnetic tail which eventually settles down to the Dungey configuration.

Leboeuf, J. N.↗

NASA NTRS1980

Viscous compressible flow about blunt bodies using a numerically generated orthogonal coordinate system

A numerical solution to the Navier-Stokes equations was obtained for blunt axisymmetric entry bodies of arbitrary shape in supersonic flow. These equations are solved on a finite difference mesh obtained from a simple numerical technique which generates orthogonal coordinates between arbitrary boundaries. The governing equations are solved in time dependent form using Stetter's improved stability three step predictor corrector method. For the present application, the metric coefficients were obtained numerically using fourth order accurate, finite difference relations and proved to be totally reliable for the highly stretched mesh used to resolve the thin viscous boundary layer. Solutions are obtained for a range of blunt body nose shapes including concavities.

Graves, R. A., Jr.↗

NASA NTRS1976

Spectroscopy of small cometary particles

Knowledge of chemical composition derived from the spectroscopy of small cometary meteoroids when they enter the Earth's atmosphere at high velocity and become visible as meteors is briefly summarized. For statistically reliable results it is necessary to have a large number of observations and this requires the photography of relatively faint meteors, which are considerably more numerous than the bright fireballs. The improvement in observational techniques and cameras has made it possible to extend the range of spectrographic data from cometary particles of mass between a kilogram and a gram down to those of only 0.01 or 0.001 grams in mass. The corresponding increase in the number-value of the statistics was from less than 100 data points to several 1000.

Millman, P. M.↗

NASA NTRS1981

Physics of auroral arc formation; Proceedings of the Chapman Conference on Formation of Auroral Arcs, Fairbanks, AK, July 21-25, 1980

The results of satellite, balloon-borne, and ground-based observations of auroral phenomena are discussed in terms of the morphology of auroral arcs and the behavior and characteristics of auroral electrons and ions. Further attention is given to auroral electric fields and field aligned currents, and to the development of models of auroral potential structures and energization of auroral particles. A simulation of space plasma phenomena is presented, along with numerical simulations of auroral potential structures and related phenomena, including the V-potential double layers and auroral arc deformations, and a simulation of auroral arcs. Plasma waves observed on auroral field lines and in laboratory conditions are reported, and theoretical studies of waves and turbulence in auroral plasmas are described.

Akasofu, S.-I.↗

NASA NTRS1982

Modulation of energetic particle fluxes by a mixed mode of transverse and compressional waves

Modulation characteristics of particle fluxes in the presence of a mixed mode of compressional and transverse magnetic waves at hydromagnetic frequencies are investigated through kinetic perturbation of the distribution function. The magnetospheric medium where the particles are modulated contains both the magnetic and pressure gradients. The modulation features are found to be strongly dependent on the energy and pitch angle of the particles. Drifting particles can resonate with waves whose phase velocities are close to their drift velocities. When this occurs, the modulation amplitudes become significantly large and large phase shifts will occur. It is pointed out that resonance is important for particles with mid pitch angles (40-70 deg). The phase shift between the particle modulations and the magnetic field oscillations are strongly controlled by the combined effects of transverse and compressional wave components and/or the occurrence of drift resonance. Numerical calculations are performed using the dispersion relation of drift mirror Alfven waves as an example of waves with both compressional and transverse components.

Lin, C. S.↗

NASA NTRS1984

Optimization of fringe-type laser anemometers for turbine engine component testing

The fringe type laser anemometer is analyzed using the Cramer-Rao bound for the variance of the estimate of the Doppler frequency as a figure of merit. Mie scattering theory is used to calculate the Doppler signal wherein both the amplitude and phase of the scattered light are taken into account. The noise from wall scatter is calculated using the wall bidirectional reflectivity and the irradiance of the incident beams. A procedure is described to determine the optimum aperture mask for the probe volume located a given distance from a wall. The expected performance of counter type processors is also discussed in relation to the Cramer-Rao bound. Numerical examples are presented for a coaxial backscatter anemometer.

Seasholtz, R. G.↗

NASA NTRS1984

Optimization of fringe-type laser anemometers for turbine engine component testing

The fringe type laser anemometer is analyzed using the Cramer-Rao bound for the variance of the estimate of the Doppler frequency as a figure of merit. Mie scattering theory is used to calculate the Doppler signal wherein both the amplitude and phase of the scattered light are taken into account. The noise from wall scatter is calculated using the wall bidirectional reflectivity and the irradiance of the incident beams. A procedure is described to determine the optimum aperture mask for the probe volume located a given distance from a wall. The expected performance of counter type processors is also discussed in relation to the Cramer-Rao bound. Numerical examples are presented for a coaxial backscatter anemometer. Previously announced in STAR as N84-25019

Seasholtz, R. G.↗

NASA NTRS1984

Some features of auroral electric fields as seen in 2D numerical simulations

Results of 2D plasma simulations are presented and related to auroral observations. The formation of V-shaped potentials is studied with a 2 1/2 dimensional electrostatic particle-in-cell code for a magnetized plasma. It is shown that amplitudes for perpendicular electric fields are larger than for parallel electric fields, and for Te less than 100 eV, the amplitudes are comparable to the electric fields associated with the electrostatic shocks observed from the S3-3 satellite. The excitation of electrostatic ion-cyclotron EIC waves which occurs in the region below the parallel potential drop is discussed. In auroral plasmas EIC waves are observed above the V-shaped double layers in association with ion beams and field-aligned currents. The results also show that oppositely directed electric fields in the center and at the edges of the simulation region produce oppositely directed currents. Precipitating auroral ions in association with electron inverted-V events are seen by the DMSP-F6 satellite.

Thiemann, H.↗

NASA NTRS1986

NASA-Chinese Aeronautical Establishment (CAE) Symposium

Several topics relative to combustion research are discussed. A numerical study of combustion processes in afterburners; the modeling of turbulent, reactive flow; gas turbine research; modeling of dilution jet flow fields; and chemical shock tubes as tools for studying high-temperature chemical kinetics are among the topics covered.

Source record↗

NASA NTRS1986

Calculation of external-internal flow fields for mixed-compression inlets

Supersonic inlet flows with mixed external-internal compressions were computed using a combined implicit-explicit (Beam-Warming-Steger/MacCormack) method for solving the three-dimensional unsteady, compressible Navier-Stokes equations in conservation form. Numerical calculations were made of various flows related to such inlet operations as the shock-wave intersections, subsonic spillage around the cowl lip, and inlet started versus unstarted conditions. Some of the computed results were compared with wind tunnel data.

Chyu, W. J.↗

NASA NTRS1987

Calculation of external-internal flow fields for mixed-compression inlets

Supersonic inlet flows with mixed external-internal compressions were computed using a combined implicit-explicit (Beam-Warming-Steger/MacCormack) method for solving the three-dimensional unsteady, compressible Navier-Stokes equations in conservation form. Numerical calculations were made of various flows related to such inlet operations as the shock-wave intersections, subsonic spillage around the cowl lip, and inlet started versus unstarted conditions. Some of the computed results were compared with wind tunnel data.

Chyu, W. J.↗