Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical Relativity”

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,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Dispersion characteristics of three-dimensional dielectric-loaded grating for terahertz Smith-Purcell radiation

In this paper, a dielectric-loaded grating for Smith-Purcell device is proposed. The three-dimensional (3D) analytical theory for hot dispersion relation is obtained by using field matched method, which is solved by numerical simulations. The first and second order growth rates for the proposal model are analyzed, which is obtained by expanding hot dispersion equation at the operating point. The results show that the dispersion can be effectively weakened by introducing dielectric-loaded grating, in which the cutoff frequency is affected by the grating thickness. The dispersion curve becomes flatter and shifts towards lower frequency at the optimum grating parameters. The 3D particle-in-cell (PIC) simulation is also performed and the results are in good agreement with theoretical calculations. Comparing the first order growth rate with the second one, it reveals that the discrepancy is small when electron beam parameters are selected with small values. Otherwise, the discrepancy is large and cannot be ignored. To accurately describe the process of beam-wave interaction, the second order growth rate is necessary to apply.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Gravitational Self-force Errors of Poisson Solvers on Adaptively Refined Meshes

An error in the gravitational force that the source of gravity induces on itself (a self-force error) violates both the conservation of linear momentum and the conservation of energy. If such errors are present in a self-gravitating system and are not sufficiently random to average out, the obtained numerical solution will become progressively more unphysical with time: the system will acquire or lose momentum and energy due to numerical effects. In this paper, we demonstrate how self-force errors can arise in the case where self-gravity is solved on an adaptively refined mesh when the refinement is nonuniform. Here, we provide the analytical expression for the self-force error and numerical examples that demonstrate such self-force errors in idealized settings. We also show how these errors can be corrected to an arbitrary order by straightforward addition of correction terms at the refinement boundaries.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1992

Numerical simulations of the electrodynamic interactions between the Tethered-Satellite-System and space plasma

The first Tethered-Satellite-System (TSS-1), scheduled for a flight in late 1992, is expected to provide relevant information related to the concept of generating an emf in a 20-km-long (or longer) conducting wire. This paper presents numerical simulations of the electrodynamic interactions between the TSS system and space plasma, using a 2D and 3D models of the system. The 2D case code simulates the motion of a long cylinder past a plasma, which is composed of electrons and H(+) ions. The system is solved by allowing the plasma to flow past the cylinder with an imposed magnetic field. The more complex 3D case is considered to study the dynamics in great detail. Results of 2D simulation show that the interaction of a satellite with plasma flowing perpendicularly to the magnetic field results in an enhancement in the current collection.

Vashi, Bharat I.↗
NASA NTRS1984

A numerical weather prediction system designed to simulate atmospheric downburst phenomena

It is pointed out that an increase in the understanding of weather-related aircraft accidents can save hundreds of human lives and million of dollars. A better understanding regarding the interaction between aircraft operation and severe weather conditions can be obtained with the aid of flight simulator facilities. It is shown that numerical weather modeling is one of the most precise and cost-effective inputs for flight simulators in the long run. A comprehensive weather modeling system is being developed for the simulation of different scales of atmospheric phenomena. The modeling system utilizes two numerical weather models, including the Mesoscale Atmospheric Simulation system, and the Terminal Area Simulation System.

Chuang, S.↗
DOE OSTI2021

S-matrix bootstrap in 3+1 dimensions: regularization and dual convex problem

The S-matrix bootstrap maps out the space of S-matrices allowed by analyticity, crossing, unitarity, and other constraints. For the 2 → 2 scattering matrix S 2→2 such space is an infinite dimensional convex space whose boundary can be determined by maximizing linear functionals. On the boundary interesting theories can be found, many times at vertices of the space. Here we consider 3 + 1 dimensional theories and focus on the equivalent dual convex minimization problem that provides strict upper bounds for the regularized primal problem and has interesting practical and physical advantages over the primal problem. Its variables are dual partial waves k ℓ (s) that are free variables, namely they do not have to obey any crossing, unitarity or other constraints. Nevertheless they are directly related to the partial waves f ℓ (s), for which all crossing, unitarity and symmetry properties result from the minimization. Numerically, it requires only a few dual partial waves, much as one wants to possibly match experimental results. We consider the case of scalar fields which is related to pion physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Computation of resistive instabilities by matched asymptotic expansions

Here, we present a method for determining the linear resistive magnetohydrodynamic (MHD) stability of an axisymmetric toroidal plasma, based on the method of matched asymptotic expansions. The plasma is partitioned into a set of ideal MHD outer regions, connected through resistive MHD inner regions about singular layers where q = m/n, with m and n toroidal mode numbers, respectively, and q the safety factor. The outer regions satisfy the ideal MHD equations with zero-frequency, which are identical to the Euler-Lagrange equations for minimizing the potential energy delta W. The solutions to these equations go to infinity at the singular surfaces. The inner regions satisfy the equations of motion of resistive MHD with a finite eigenvalue, resolving the singularity. Both outer and inner regions are solved numerically by newly developed singular Galerkin methods, using specialized basis functions. These solutions are matched asymptotically, providing a complex dispersion relation which is solved for global eigenvalues and eigenfunctions in full toroidal geometry. The dispersion relation may have multiple complex unstable roots, which are found by advanced root-finding methods. These methods are much faster and more robust than the previous numerical methods. The new methods are applicable to more challenging high-pressure and strongly shaped plasma equilibria and generalizable to more realistic inner region dynamics. In the thermonuclear regime, where the outer and inner regions overlap, they are also much faster and more accurate than the straight-through methods, which treat the resistive MHD equations in the whole plasma volume.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS1972

Weldability, strength, and high temperature stability of chemically vapor deposited tungsten

Three types of CVD tungsten (fluoride-produced, chloride-produced and the combination of the two which is termed duplex) were evaluated to determine their weldability, high temperature strength and structural stability during 5000 hour exposure to temperatures of 1540 C and 1700 C. Each type of CVD tungsten could be successfully electron beam welded but the results for the chloride product were not as satisfactory as those of the other two materials. The high temperature strength behavior of the three materials did not differ greatly. However a large difference was noted for the grain growth behavior of the two basic CVD tungsten materials. Fluoride tungsten was found to be relatively stable while for the most part the grain size of chloride tungsten increased appreciably. The examination of freshly fractured surfaces with a scanning electron microscope revealed numerous bubbles in the fluoride material following its exposure to 1700 C for 5000 hours. Less severe thermal treatments produced relatively few bubbles in this material. Only at certain locations within the chloride material associated with the interruption of tungsten were bubbles noted.

Bryant, W. A.↗
DOE OSTI2017

The Radiative Efficiency and Spectra of Slowly Accreting Black Holes from Two-temperature GRRMHD Simulations

In this paper, we present axisymmetric numerical simulations of radiatively inefficient accretion flows onto black holes combining general relativity, magnetohydrodynamics, self-consistent electron thermodynamics, and frequency-dependent radiation transport. We investigate a range of accretion rates up to ${10}^{-5}\,{\dot{M}}_{\mathrm{Edd}}$ onto a ${10}^{8}\,{M}_{\odot }$ black hole with spin ${a}_{\star }=0.5$. We report on averaged flow thermodynamics as a function of accretion rate. We present the spectra of outgoing radiation and find that it varies strongly with accretion rate, from synchrotron-dominated in the radio at low $\dot{M}$ to inverse-Compton-dominated at our highest $\dot{M}$. In contrast to canonical analytic models, we find that by $\dot{M}\approx {10}^{-5}\,{\dot{M}}_{\mathrm{Edd}}$, the flow approaches $\sim 1 \% $ radiative efficiency, with much of the radiation due to inverse-Compton scattering off Coulomb-heated electrons far from the black hole. Finally, these results have broad implications for modeling of accreting black holes across a large fraction of the accretion rates realized in observed systems.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2004

Topographic Evidence for Eruptive Style Changes and Magma Evolution of Small Plains-style Volcanoes on Earth and Mars

Topographic profiles and surface characteristics of small (5 - 25 km diameter) plains-style shield volcanoes on the eastern Snake River Plain (ESRP) provide a method to evaluate eruptive processes and magmatic evolution on Martian volcanic plains. The ESRP is an ideal place to observe Mars-like volcanic features where hundreds of small monogenetic basaltic shields dominate the volcanic-sedimentary depositional sequence, and numerous planetary analogues are evident: coalescent mafic shields, hydromagmatic explosive eruptions, the interaction of lava flows with surficial water and glacial ice, and abundant eolian sand and loess. Single flows cannot be correlated over great distances, and are spatially restricted. These relations are useful for planetary exploration when inferring volcanic evolutionary patterns in lava plains represented by numerous eruptive vents. High spatial resolution imagery and digital topographic data for Mars from MOC, MOLA, and THEMIS is allowing for improvements in the level of detail of stratigraphic mapping of fields of small (< 25 km in diameter) volcanoes as well as studies of the morphological characteristics of individual volcanoes. In order to compare Mars and Earth volcanic features, elevation data from U.S.G.S. 10-meter digital elevation models (DEMs) and high-precision GPS field measurements are used in this study to generate approx. 20m spacing topographic profiles from which slope and surface morphology can be extracted. Average ESRP flank and crater slopes are calculated using 100 - 200 m spacing for optimum comparison to MOLA data, and to reduce the effects of surface irregularities.

Hughes, S. S.↗
DOE OSTI2020

Comparison of Numerical Methods for Modeling the Wave Field Effects Generated by Individual Wave Energy Converters and Multiple Converter Wave Farms

This numerical study compares the wave field generated by the spectral wave action balance code, SNL-SWAN, to the linear-wave boundary-element method (BEM) code, WAMIT. The objective of this study is to assess the performance of SNL-SWAN for modeling wave field effects produced by individual wave energy converters (WECs) and wave farms comprising multiple WECs by comparing results from SNL-SWAN with those produced by the BEM code WAMIT. BEM codes better model the physics of wave-body interactions and thus simulate a more accurate near-field wave field than spectral codes. In SNL-SWAN, the wave field’s energy extraction is modeled parametrically based on the WEC’s power curve. The comparison between SNL-SWAN and WAMIT is made over a range of incident wave conditions, including short-, medium-, and long-wavelength waves with various amounts of directional spreading, and for three WEC archetypes: a point absorber (PA), a pitching flap (PF) terminator, and a hinged raft (HR) attenuator. Individual WECs and wave farms of five WECs in various configuration were studied with qualitative comparisons made of wave height and spectra at specific locations, and quantitative comparisons of the wave fields over circular arcs around the WECs as a function of radial distance. Results from this numerical study demonstrate that in the near-field, the difference between SNL-SWAN and WAMIT is relatively large (between 20% and 50%), but in the far-field from the array the differences are minimal (between 1% and 5%). The resultant wave field generated by the two different numerical approaches is highly dependent on parameters such as: directional wave spreading, wave reflection or scattering, and the WEC’s power curve.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Towards the scale-up production of cellulose nanocomposites using melt processing: A critical review on structure-processing-property relationships

Cellulose nanomaterials (CNMs) naturally exist in plant biomass. The success of extraction of CNMs opened up a new era of using plant biomass for innovative industrial applications. Because CNMs are abundant, renewable, biodegradable, transparent, light weight and low in cost, they are ideal materials for large volume applications such as packaging, automotive, building and infrastructure. In many potential application areas, CNM-enable products appear in a composite form, mostly polymer composites. Additionally, the industrial-scale manufacturing of CNM/thermoplastic composites remains as a set of unsolved problems for academia and industry. A prime challenge in applications is the nanoscale dispersion of CNMs in thermoplastic matrices during melt processing. Both bench-scale and pilot-scale studies have been conducted to solve the dispersion issue of CNMs. In this article, research related to the dispersion of CNMs in thermoplastic matrices during melt processing were critically reviewed. All research papers were classified into three groups: chemically-aided dispersion, physically-aided dispersion and mechanically-aided dispersion. Numerous factors affect the CNM dispersion and the mechanical performance of its nanocomposites. There are material-related factors, including CNM types and forms, polymer matrices, surface modification, coupling agents, etc. Furthermore, extrusion processing parameters also play a significant role, covering screw rotation speed, extrusion barrel temperature settings and screw design. In addition, the material-related factors interact with the processing-related factors. Understanding all factors and their interactions are important for moving CNM nanocomposites research a step further towards industrial-relevant production, which is the final ambitious goal of this manuscript.

42 ENGINEERING↗
NASA NTRS1990

Energy and angular momentum transfer in binary galaxies

The authors numerically studied tidal effects of a massive perturber on a satellite galaxy. The model consists of a spherical satellite galaxy and a point mass perturber and the encounter is non-penetrating. A wide range of density ratios and eccentricities of the relative orbits have been considered. The disruption of the satellite galaxy has been observed when the numerical value of the fractional change in the energy is greater than two. The changes in the energy and angular momentum show smooth variation in the case of unbound orbits and irregular variation in the bound orbit cases. It is shown that for a constant pericentral distance, increasing the density ratio decreases the tidal effects; and for a given density ratio an increase in the eccentricity decreases the tidal effects.

Namboodiri, P. M. S.↗
DOE OSTI2022

E-transit-bench: simulation platform for analyzing electric public transit bus fleet operations

When electrified transit systems make grid aware choices, improved social welfare is achieved by reducing grid stress, reducing system loss, and minimizing power quality issues. Electrifying transit fleet has numerous challenges like non availability of buses during charging, varying charging costs and so on, that are related the electric grid behavior. However, transit systems do not have access to the information about the co-evolution of the grid's power flow and therefore cannot account for the power grid's needs in its day-to-day operation. In this paper we propose a framework of transportation-grid co-simulation, analyzing the spatio-temporal interaction between the transit operations with electric buses and the power distribution grid. Real-world data for a day's traffic from Chattanooga city's transit system is simulated in SUMO and integrated with a realistic distribution grid simulation (using GridLAB-D) to understand the grid impact due to transit electrification. Charging information is obtained from the transportation simulation to feed into grid simulation to assess the impact of charging. We also discuss the impact to the grid with higher degree of transit electrification that further necessitates such an integrated transportation-grid co-simulation to operate the integrated system optimally. Our future work includes extending the platform for optimizing the charging and trip assignment operations.

Sen, Rishav↗
DOE OSTI2020

Direct numerical simulations of turbulent reacting flows with shock waves and stiff chemistry using many-core/GPU acceleration

Compressible reacting flows may display sharp spatial variation related to shocks, contact discontinuities or reactive zones embedded within relatively smooth regions. The presence of such phenomena emphasizes the relevance of shock-capturing schemes such as the weighted essentially non-oscillatory (WENO) scheme as an essential ingredient of the numerical solver. However, these schemes are complex and have more computational cost than the simple high-order compact or non-compact schemes. In this paper, we present the implementation of a seventh-order, minimally-dissipative mapped WENO (WENO7M) scheme in a newly developed direct numerical simulation (DNS) code called KAUST Adaptive Reactive Flows Solver (KARFS). In order to make efficient use of the computer resources and reduce the solution time, without compromising the resolution requirement, the WENO routines are accelerated via graphics processing unit (GPU) computation. The performance characteristics and scalability of the code are studied using different grid sizes and block decomposition. Furthermore, the performance portability of KARFS is demonstrated on a variety of architectures including NVIDIA Tesla P100 GPUs and NVIDIA Kepler K20X GPUs. In addition, the capability and potential of the newly implemented WENO7M scheme in KARFS to perform DNS of compressible flows is also demonstrated with model problems involving shocks, isotropic turbulence, detonations and flame propagation into a stratified mixture with complex chemical kinetics.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Fatigue, pilot deviations and time of day

The relationships between pilot fatigue, pilot deviations, reported incidents, and time of day are examined. A sample of 200 Aviation Safety Reporting System (ASRS) reports were analyzed from 1985 and 200 reports from 1987, plus 100 reports from late 1987 and early 1988 that were selected because of possible association with fatigue. The FAA pilot deviation data and incident data were analyzed in relation to denominator data that summarized the hourly operations (landings and takeoffs of scheduled flights) at major U.S. airports. Using as numerators FAA data on pilot deviations and incidents reported to the FAA, the rates by time of day were calculated. Pilot age was also analyzed in relation to the time of day, phase of flight, and type of incident.

Baker, Susan P.↗
NASA NTRS2005

The QBO as Potential Amplifier and Conduit to Lower Altitudes of Solar Cycle Influence

In several papers, the solar cycle (SC) effect in the lower atmosphere has been linked to the Quasi-biennial Oscillation (QBO), which is generated primarily by small-scale gravity waves. Salby and Callaghan (2000) analyzed the observed zonal winds of the QBO over more than 40 years and found that it contains a relatively large SC signature at 20 km. Following up on an earlier 2D study with our global-scale Numerical Spectral Model (NSM), we discuss here a 3D study with the QBO under the influence of the SC. For a SC period of 10 years, the amplitude of the relative variations of radiative forcing is taken to vary for simplicity from 0.2% at the surface to 2% at 50 km to 20% at 100 km and above. Covering a limited time span of 40 years, this model produces in the lower stratosphere a relatively large modulation of the QBO, which appears to be related to the SC and is in qualitative agreement with the observations. Some of the energy in the QBO, confined to low latitudes primarily, is redistributed globally by the meridional circulation and planetary waves presumably, so that a measurable SC modulation is generated in the tropospheric temperatures of the polar regions. Further studies are needed, (1) to determine whether the effect is real and prevails in more extensive simulations and whether the results are robust when shorter integration steps are employed, and (2) to explore the mechanism(s) that may ample the apparent SC influence of the UV radiation extending into the lower atmosphere. Quasi-decadal oscillations, generated internally by the QBO interacting with the seasonal cycles, may interfere with or aid the SC effect.

Mayr, Hans G.↗
DOE OSTI2018

Numerical Tools for Improved Convergence of Meshfree Peridynamic Discretizations

Peridynamic models have been employed to simulate a broad range of engineering applications concerning material failure and damage, with the majority of these simulations using a meshfree discretization. This chapter reviews that meshfree discretization, related issues present in peridynamic convergence studies, and possible remedies proposed in the literature. In particular, we discuss two numerical tools, partial-volume algorithms and influence functions, to improve the convergence behavior of numerical solutions in peridynamics. Numerical studies in this chapter involve static and dynamic simulations for linear elastic state-based peridynamic problems.

Seleson, Pablo D.↗
DOE OSTI2024

Two-loop master integrals for leading-color $$ pp\to t\overline{t}H $$ amplitudes with a light-quark loop

Abstract We compute the two-loop master integrals for leading-color QCD scattering amplitudes including a closed light-quark loop in$$ t\overline{t}H $$ t t ¯ H production at hadron colliders. Exploiting numerical evaluations in modular arithmetic, we construct a basis of master integrals satisfying a system of differential equations inϵ-factorized form. We present the analytic form of the differential equations in terms of a minimal set of differential one-forms. We explore properties of the function space of analytic solutions to the differential equations in terms of iterative integrals which can be exploited for studying the analytic form of related scattering amplitudes. Finally, we solve the differential equations using generalized series expansions to numerically evaluate the master integrals in physical phase space. As the first computation of a set of two-loop seven-scale master integrals, our results provide valuable input for analytic studies of scattering amplitudes in processes involving massive particles and a large number of kinematic scales.

Physics↗