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 577 records · Page 32

Lehigh University Industrial Assessment Center (Final Technical Report for the Period 2016 to 2021)

This is the final technical report summarizing the activities Industrial Assessment Center at Lehigh University sponsored by the Department of Energy for the period September 1, 2016 till September 30, 2021. During this period, the Industrial Assessment Center at Lehigh University was successful in training numerous energy engineers of the future, and helped many manufacturing plants in New Jersey and Eastern Pennsylvania. The center and its activities helped save a great deal of energy, related fuel costs and along the way mitigated tons of CO2 emissions. The report covers many of the numerical details related to the center activities. The Industrial Center at Lehigh University is one of the most successful ones and this is made evident by the statistics included as part of this report and best center award bestowed upon the center.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

National plans for aircraft icing and improved aircraft icing forecasts and associated warning services

Recently, the United States has increased its activities related to aircraft icing in numerous fields: ice phobics, revised characterization of icing conditions, instrument development/evaluation, de-ice/anti-ice devices, simulated supercooled clouds, computer simulation and flight tests. The Federal Coordinator for Meteorology is involved in two efforts, one a National Plan on Aircraft Icing and the other a plan for Improved Aircraft Icing Forecasts and Associated Warning Services. These two plans will provide an approved structure for future U.S. activities related to aircraft icing. The recommended activities will significantly improve the position of government agencies to perform mandated activities and to enable U.S. manufacturers to be competitive in the world market.

Pass, Ralph P.↗

Numerical experiments on the climatic sensitivity of an atmospheric hydrologic cycle

It is shown for an intermediate numerical model that fractional cloudiness and relative humidity decrease with increasing temperature. The fractional cloudiness decreases at a rate about 1 per deg K. This occurs in spite of an increase in the evaporation, water transport, condensation, precipitation and cloud water content with increasing temperature. These results are quite similar to those found from models with more highly parameterized clouds, notably the NCAR model. The fractional cloudiness in this model is measured by the fractional coverage of total cloud water and the fractional coverage of positive condensation, in addition to the relative humidity. It is also shown that some of the characteristics of a temperate climate can be simulated in an intermediate numerical model with periodic, antisymmetric and symmetric boundary conditions on an f plane. Intermediate models of this sort may therefore be useful to investigate general questions about the earth's hydrologic cycle on climatic space and time scales

Roads, J. O.↗

Comet-solar wind interaction through ion-proton beam instability

The acceleration of cometary ions in the distant cometary tail by an electromagnetic beam instability, which is caused by the relative drift velocity between the cometary ions and solar protons is studied. The linear phase of the instability is analyzed by numerical evaluation of the Vlasov dispersion relation. A particle simulation is used to study the nonlinear phase of the instability. Right-hand polarized cyclotron waves are made unstable by the presence of ion beams. The cometary tail ions are accelerated mainly during the reduced-growth phase of the instability. The cyclotron waves cause momentum transfer from solar wind protons to cometary tail ions. The acceleration of cometary ions through the electromagnetic beam instability is found to be 500-5000 cm/s-squared.

Price, Channon P.↗

Analysis of the NASA White Sands Test Facility (WSTF) Test System for Friction-Ignition of Metallic Materials

Friction is a known ignition source for metals in oxygen-enriched atmospheres. The test system developed by the NASA White Sands Test Facility in response to ASTM G-94 has been used successfully to determine the relative ignition from friction of numerous metallic materials and metallic materials pairs. These results have been ranked in terms of a pressure-velocity product (PV) as measured under the prescribed test conditions. A high value of 4.1(exp 8) watts per square meter for Inconel MA 754 is used to imply resistance to friction ignition, whereas a low value of 1.04(exp 8) watts per square meter for stainless steel 304 is taken as indicating material susceptible to friction ignition. No attempt has been made to relate PV values to other material properties. This work reports the analysis of the WSTF friction-ignition test system for producing fundamental properties of metallic materials relating to ignition through friction. Three materials, aluminum, titanium, and nickel were tested in the WSTF frictional ignition instrument system under atmospheres of oxygen or nitrogen. Test conditions were modified to reach a steady state of operation, that is applied, the force was reduced and the rotational speed was reduced. Additional temperature measurements were made on the stator sample. The aluminum immediately galled on contact (reproducible) and the test was stopped. Titanium immediately ignited as a result of non-uniform contact of the stator and rotor. This was reproducible. A portion of the stator sampled burned, but the test continued. Temperature measurements on the stator were used to validate the mathematical model used for estimating the interface (stator/rotor) temperature. These interface temperature measurements and the associate thermal flux into the stator were used to distinguish material-phase transitions, chemical reaction, and mechanical work. The mechanical work was used to analyze surface asperities in the materials and to estimate a coefficient of fiction. The coefficient of fiction was analyzed in terms of material properties that is, hardness, Young's modulus and elasticity/plasticity of the material.

Shoffstall, Michael S.↗

Verification of MOOSE/Bison's Heat Conduction Solver Using Combined Spatiotemporal Convergence Analysis

Bison is a computational physics code that uses the finite element method to model the thermo-mechanical response of nuclear fuel. Since Bison is used to inform high-consequence decisions, it is important that its computational results are reliable and predictive. One important step in assessing the reliability and predictive capabilities of a simulation tool is the verification process, which quantifies numerical errors in a discrete solution relative to the exact solution of the mathematical model. One step in the verification process—called code verification—ensures that the implemented numerical algorithm is a faithful representation of the underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. In this paper, the code verification process is applied to spatiotemporal heat conduction problems in Bison. Simultaneous refinement of the discretization in space and time is employed to reveal any potential mistakes in the numerical algorithms for the interactions between the spatial and temporal components of the solution. For each verification problem, the correct spatial and temporal order of accuracy is demonstrated for both first- and second-order accurate finite elements and a variety of time-integration schemes. Furthermore, these results provide strong evidence that the Bison numerical algorithm for solving spatiotemporal problems reliably represents the underlying mathematical model in MOOSE. The selected test problems can also be used in other simulation tools that numerically solve for conduction or diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Dynamical friction from scalar dark matter in the relativistic regime

Light bosonic scalars (e.g., axions) may form clouds around black holes via superradiant instabilities or via accretion if they form some component of the dark matter. It has been suggested that their presence may lead to a distinctive dephasing of the gravitational wave signal when a small compact object spirals into a larger black hole. Motivated by this, we have studied numerically the dynamical friction force on a black hole moving at relativistic velocities in a background scalar field with an asymptotically homogeneous energy density. Here we show that the relativistic scaling is analogous to that found for supersonic collisional fluids, assuming an approximate expression for the pressure correction which depends on the velocity and scalar mass. While we focus on a complex scalar field, our results confirm the expectation that real scalars would exert a force which oscillates between positive and negative values in time with a frequency set by the scalar mass. The complex field describes the time averaged value of this force, but in a real scalar, the rapid force oscillations could, in principle, leave an imprint on the trajectory. The approximation we obtain can be used to inform estimates of dephasing in the final stages of an extreme mass ratio inspiral.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Doppler echo evaluation of pulmonary venous-left atrial pressure gradients: human and numerical model studies

The simplified Bernoulli equation relates fluid convective energy derived from flow velocities to a pressure gradient and is commonly used in clinical echocardiography to determine pressure differences across stenotic orifices. Its application to pulmonary venous flow has not been described in humans. Twelve patients undergoing cardiac surgery had simultaneous high-fidelity pulmonary venous and left atrial pressure measurements and pulmonary venous pulsed Doppler echocardiography performed. Convective gradients for the systolic (S), diastolic (D), and atrial reversal (AR) phases of pulmonary venous flow were determined using the simplified Bernoulli equation and correlated with measured actual pressure differences. A linear relationship was observed between the convective (y) and actual (x) pressure differences for the S (y = 0.23x + 0.0074, r = 0.82) and D (y = 0.22x + 0.092, r = 0.81) waves, but not for the AR wave (y = 0. 030x + 0.13, r = 0.10). Numerical modeling resulted in similar slopes for the S (y = 0.200x - 0.127, r = 0.97), D (y = 0.247x - 0. 354, r = 0.99), and AR (y = 0.087x - 0.083, r = 0.96) waves. Consistent with numerical modeling, the convective term strongly correlates with but significantly underestimates actual gradient because of large inertial forces.

Non-NASA Center↗

Sample processor for chemical analysis

An apparatus is provided which can process numerous samples that must be chemically analyzed by the application of fluids such as liquid reagents, solvents and purge gases, as well as the application of dumps for receiving the applied fluid after they pass across the sample, in a manner that permits numerous samples to be processed in a relatively short time and with minimal manpower. The processor includes a rotor which can hold numerous cartridges containing inert or adsorbent material for holding samples, and a pair of stators on opposite sides of the rotor. The stators form stations spaced along the path of the cartridges which lie in the rotor, and each station can include an aperture in one stator through which a fluid can be applied to a cartridge resting at that station, and an aperture in the other stator which can receive the fluid which has passed through the cartridge. The stators are sealed to the ends of the cartridges lying on the rotor, to thereby isolate the stations from one another.

Boettger, Heinz G.↗

Verification of numerical models for seismic fluid-structure interaction analysis of internal components in liquid-filled advanced reactors

Earthquake shaking of a liquid-filled advanced reactor induces fluid-structure interaction (FSI) between the reactor vessel (tank), its internal components, and the contained liquid. Seismic design, qualification, and risk assessment of such reactors must consider fluid-structure responses, their geometries and support conditions, and three-directional seismic inputs, all of which require numerical simulations. Numerical models used for analysis of safety-related nuclear equipment must be verified and validated. Here, this paper verifies numerical models of submerged components using analytical solutions. Seismic FSI analysis of the numerical models is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Prior analytical solutions are reworked in this paper, and calculation errors are identified and corrected. The solutions address frequencies of two concentric cylindrical pipes filled with liquid: the inner pipe is a submerged component, and the outer pipe contains the liquid. A unitless frequency coefficient is defined and presented for pipes of different materials, filled with different liquids, and with a range of dimensions for application to advanced reactors. The numerical models are verified here by comparing the lateral frequencies of submerged components with those calculated using the corrected analytical solutions. Recommendations for verification of numerical models of internal components in advanced reactors for seismic FSI analysis are provided. Although the unitless frequency coefficients and verification procedures are developed for application to advanced nuclear reactors, they are broadly applicable to FSI analysis of submerged components in liquid-filled vessels such as storage tanks, boilers, and steam generators.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Ultrasonic/Sonic Rotary-Hammer Drills

Ultrasonic/sonic rotary-hammer drill (USRoHD) is a recent addition to the collection of apparatuses based on ultrasonic/sonic drill corer (USDC). As described below, the USRoHD has several features, not present in a basic USDC, that increase efficiency and provide some redundancy against partial failure. USDCs and related apparatuses were conceived for boring into, and/or acquiring samples of, rock or other hard, brittle materials of geological interest. They have been described in numerous previous NASA Tech Briefs articles. To recapitulate: A USDC can be characterized as a lightweight, lowpower, piezoelectrically driven jackhammer in which ultrasonic and sonic vibrations are generated and coupled to a tool bit. A basic USDC includes a piezoelectric stack, an ultrasonic transducer horn connected to the stack, a free mass ( free in the sense that it can bounce axially a short distance between hard stops on the horn and the bit), and a tool bit. The piezoelectric stack creates ultrasonic vibrations that are mechanically amplified by the horn. The bouncing of the free mass between the hard stops generates the sonic vibrations. The combination of ultrasonic and sonic vibrations gives rise to a hammering action (and a resulting chiseling action at the tip of the tool bit) that is more effective for drilling than is the microhammering action of ultrasonic vibrations alone. The hammering and chiseling actions are so effective that unlike in conventional twist drilling, little applied axial force is needed to make the apparatus advance into the material of interest. There are numerous potential applications for USDCs and related apparatuses in geological exploration on Earth and on remote planets. In early USDC experiments, it was observed that accumulation of cuttings in a drilled hole causes the rate of penetration of the USDC to decrease steeply with depth, and that the rate of penetration can be increased by removing the cuttings. The USRoHD concept provides for removal of cuttings in the same manner as that of a twist drill: An USRoHD includes a USDC and a motor with gearhead (see figure). The USDC provides the bit hammering and the motor provides the bit rotation. Like a twist drill bit, the shank of the tool bit of the USRoHD is fluted. As in the operation of a twist drill, the rotation of the fluted drill bit removes cuttings from the drilled hole. The USRoHD tool bit is tipped with a replaceable crown having cutting teeth on its front surface. The teeth are shaped to promote fracturing of the rock face through a combination of hammering and rotation of the tool bit. Helical channels on the outer cylindrical surface of the crown serve as a continuation of the fluted surface of the shank, helping to remove cuttings. In the event of a failure of the USDC, the USRoHD can continue to operate with reduced efficiency as a twist drill. Similarly, in the event of a failure of the gearmotor, the USRoHD can continue to operate with reduced efficiency as a USDC.

Badescu, Mircea↗

Effect of Cyclo-Pentane Impurities on the Autoignition Reactivity and Properties of a Gasoline Surrogate Fuel

Surrogate fuels that reproduce the characteristics of full-boiling range fuels are key tools to enable numerical simulations of fuel-related processes and ensure reproducibility of experiments by eliminating batch-to-batch variability. Within the PACE initiative, a surrogate fuel for regular-grade E10 (10%vol ethanol) gasoline representative of a U.S. market gasoline, termed PACE-20, was developed and adopted as baseline fuel for the consortium. Although extensive testing demonstrated that PACE-20 replicates the properties and combustion behavior of the full-boiling range gasoline, several concerns arose regarding the purity level required for the species that compose PACE-20. This is particularly important for cyclo-pentane, since commercial-grade cyclo-pentane typically shows 60%–85% purity. Here, in the present work, the effects of the purity level of cyclo-pentane on the properties and combustion characteristics of PACE-20 were studied. Chemical kinetic simulations were performed to predict the effects of cyclo-pentane impurities on the properties, octane rating, and autoignition reactivity under homogeneous charge compression-ignition conditions of PACE-20. From the numerical results, cyclo-pentane with 85% purity or higher is required to reasonably match both the research octane number and motor octane number of the target gasoline. Finally, homogeneous charge compression-ignition engine simulations show that impurities have only a modest effect on reactivity at naturally aspirated conditions, but cyclo-pentane purity is critical to properly replicate the pressure dependency of the reactivity.

33 ADVANCED PROPULSION SYSTEMS↗

Review of some vortex relations

The evaluation of the circulation from numerical solutions of the momentum and energy equations is discussed for incompressible and compressible flows. It is shown how artificial damping directly influences the time ratio of change of the circulation.

Krause, E.↗

Review of some vortex relations

The evaluation of the circulation from numerical solutions of the momentum and energy equations is discussed for incompressible and compressible flows. It is shown how artificial damping directly influences the time rate of change of the circulation.

Krause, E.↗

Relating Satellite-Based Contrail Detection to NWA Output

Contrail-induced cloud cover has been shown to be a significant factor in regional climate change over the United States of America. As air traffic increases, the potential for globally significant impacts also rises. To better understand and predict these potential climatic effects, it is necessary to develop models that can accurately represent contrail properties based on ambient atmospheric variables including temperature, relative humidity and winds. Several high-resolution numerical weather analyses (NWA) including the 20-km Rapid Update Cycle (RUC) and the University of Oklahoma Center for Analysis and Prediction of Storms Advanced Regional Prediction System (ARPS) can provide the temperature, humidity and wind information necessary to diagnose contrail formation. One outstanding problem that must be addressed to achieve a realistic simulation of contrails is the large uncertainties in upper tropospheric relative humidity (UTH) in numerical weather analyses. Current numerical weather analyses tend to underestimate UTH due to large dry biases in the balloon soundings used to construct the analyses. To evaluate each of the models, we match several months of contrail properties derived from satellite and surface observations to the NWA-derived humidity, vertical velocity, wind shear and atmospheric stability. The relationships between several contrail properties (including optical depth and areal contrail coverage) and the NWA-derived statistics will be analyzed to determine under which atmospheric conditions widespread contrail outbreaks are favored.

Duda, David P.↗

Numerical Arc Segmentation Algorithm for a Radio Conference-NASARC, Version 2.0: User's Manual

The information contained in the NASARC (Version 2.0) Technical Manual (NASA TM-100160) and the NASARC (Version 2.0) User's Manual (NASA TM-100161) relates to the state of the Numerical Arc Segmentation Algorithm for a Radio Conference (NASARC) software development through October 16, 1987. The technical manual describes the NASARC concept and the algorithms which are used to implement it. The User's Manual provides information on computer system considerations, installation instructions, description of input files, and program operation instructions. Significant revisions have been incorporated in the Version 2.0 software over prior versions. These revisions have enhanced the modeling capabilities of the NASARC procedure while greatly reducing the computer run time and memory requirements. Array dimensions within the software have been structured to fit into the currently available 6-megabyte memory capacity of the International Frequency Registration Board (IFRB) computer facility. A piecewise approach to predetermined arc generation in NASARC (Version 2.0) allows worldwide scenarios to be accommodated within these memory constraints while at the same time reducing computer run time.

Whyte, Wayne A., Jr.↗

Ameliorating the Courant-Friedrichs-Lewy condition in spherical coordinates: A double FFT filter method for general relativistic MHD in dynamical spacetimes

Numerical simulations of merging compact objects and their remnants form the theoretical foundation for gravitational wave and multimessenger astronomy. While Cartesian-coordinate-based adaptive mesh refinement is commonly used for simulations, spherical-like coordinates are more suitable for nearly spherical remnants and azimuthal flows due to lower numerical dissipation in the evolution of fluid angular momentum, as well as requiring fewer numbers of computational cells. However, the use of spherical coordinates to numerically solve hyperbolic partial differential equations can result in severe Courant-Friedrichs-Lewy (CFL) stability condition time step limitations, which can make simulations prohibitively expensive. This paper addresses this issue for the numerical solution of coupled spacetime and general relativistic magnetohydrodynamics evolutions by introducing a double fast Fourier transform (FFT) filter and implementing it within the fully message passing interface (mpi)-parallelized sphericalnr framework in the einstein toolkit. In conclusion, we demonstrate the effectiveness and robustness of the filtering algorithm by applying it to a number of challenging code tests, and show that it passes these tests effectively, demonstrating convergence while also increasing the time step significantly compared to unfiltered simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Advanced Computational Aeroacoustics Methods for Fan Noise Prediction

Direct computation of fan noise is presently not possible. One of the major difficulties is the geometrical complexity of the problem. In the case of fan noise, the blade geometry is critical to the loading on the blade and hence the intensity of the radiated noise. The precise geometry must be incorporated into the computation. In computational fluid dynamics (CFD), there are two general ways to handle problems with complex geometry. One way is to use unstructured grids. The other is to use body fitted overset grids. In the overset grid method, accurate data transfer is of utmost importance. For acoustic computation, it is not clear that the currently used data transfer methods are sufficiently accurate as not to contaminate the very small amplitude acoustic disturbances. In CFD, low order schemes are, invariably, used in conjunction with unstructured grids. However, low order schemes are known to be numerically dispersive and dissipative. dissipative errors are extremely undesirable for acoustic wave problems. The objective of this project is to develop a high order unstructured grid Dispersion-Relation-Preserving (DRP) scheme. would minimize numerical dispersion and dissipation errors. contains the results of the funded portion of the project. scheme on an unstructured grid has been developed. constructed in the wave number space. The characteristics of the scheme can be improved by the inclusion of additional constraints. Stability of the scheme has been investigated. Stability can be improved by adopting the upwinding strategy.

Envia, Edmane↗