SEARCH · Engineering Papers

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Numerical simulation of the helium gas spin-up channel performance of the relativity gyroscope

The dependence of the spin-up system efficiency on each geometrical parameter of the spin-up channel and the exhaust passage of the Gravity Probe-B (GPB) is individually investigated. The spin-up model is coded into a computer program which simulates the spin-up process. Numerical results reveal optimal combinations of the geometrical parameters for the ultimate spin-up performance. Comparisons are also made between the numerical results and experimental data. The experimental leakage rate can only be reached when the gap between the channel lip and the rotor surface increases beyond physical limit. The computed rotating frequency is roughly twice as high as the measured ones although the spin-up torques fairly match.

Karr, Gerald R.↗

NASA NTRS1975

Study of viscous flow about airfoils by the integro-differential method

An integro-differential method was used for numerically solving unsteady incompressible viscous flow problems. A computer program was prepared to solve the problem of an impulsively started 9% thick symmetric Joukowski airfoil at an angle of attack of 15 deg and a Reynolds number of 1000. Some of the results obtained for this problem were discussed and compared with related work completed previously. Two numerical procedures were used, an Alternating Direction Implicit (ADI) method and a Successive Line Relaxation (SLR) method. Generally, the ADI solution agrees well with the SLR solution and with previous results are stations away from the trailing edge. At the trailing edge station, the ADI solution differs substantially from previous results, while the vorticity profiles obtained from the SLR method there are in good qualitative agreement with previous results.

Wu, J. C.↗

NASA NTRS2015

Tetrahedral-Mesh Simulation of Turbulent Flows with the Space-Time Conservative Schemes

Direct numerical simulations of turbulent flows are predominantly carried out using structured, hexahedral meshes despite decades of development in unstructured mesh methods. Tetrahedral meshes offer ease of mesh generation around complex geometries and the potential of an orientation free grid that would provide un-biased small-scale dissipation and more accurate intermediate scale solutions. However, due to the lack of consistent multi-dimensional numerical formulations in conventional schemes for triangular and tetrahedral meshes at the cell interfaces, numerical issues exist when flow discontinuities or stagnation regions are present. The space-time conservative conservation element solution element (CESE) method - due to its Riemann-solver-free shock capturing capabilities, non-dissipative baseline schemes, and flux conservation in time as well as space - has the potential to more accurately simulate turbulent flows using unstructured tetrahedral meshes. To pave the way towards accurate simulation of shock/turbulent boundary-layer interaction, a series of wave and shock interaction benchmark problems that increase in complexity, are computed in this paper with triangular/tetrahedral meshes. Preliminary computations for the normal shock/turbulence interactions are carried out with a relatively coarse mesh, by direct numerical simulations standards, in order to assess other effects such as boundary conditions and the necessity of a buffer domain. The results indicate that qualitative agreement with previous studies can be obtained for flows where, strong shocks co-exist along with unsteady waves that display a broad range of scales, with a relatively compact computational domain and less stringent requirements for grid clustering near the shock. With the space-time conservation properties, stable solutions without any spurious wave reflections can be obtained without a need for buffer domains near the outflow/farfield boundaries. Computational results for the isotropic turbulent flow decay, at a relatively high turbulent Mach number, show a nicely behaved spectral decay rate for medium to high wave numbers. The high-order CESE schemes offer very robust solutions even with the presence of strong shocks or widespread shocklets. The explicit formulation in conjunction with a close to unity theoretical upper Courant number bound has the potential to offer an efficient numerical framework for general compressible turbulent flow simulations with unstructured meshes.

Chang, Chau-Lyan↗

DOE OSTI2022

Dispersion relation and instability for an anisotropic nonuniform flowing plasma

A generalized formula for wave instability is developed for an anisotropic nonuniform plasma with finite flows and temperatures. Six-moment fluid equations are solved to give the analytic expression for wave instability in arbitrarily nonuniform plasmas. The analytic formula explicitly states the dependence of wave instability on the nonuniformities of number density, flow velocity, and anisotropic or isotropic pressure. The accuracy of the formalism is verified by a numerical calculation of implicit dispersion relations in complex Fourier space. The analysis shows that nonuniformity plays a critical role in plasma instability, while the flow velocity and anisotropic pressures determine the growth rate of the instability. Lastly, the instability diagram and associated instability criterion for anisotropy-driven instability are introduced as applications of the formalism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1992

Ground-based measurements of galactic cosmic ray fragmentation in shielding

The mean free path for nuclear interactions of galactic cosmic-rays is comparable to shielding and tissue thicknesses present in human interplanetary exploration, resulting in a significant fraction of nuclear reaction products at depth. In order to characterize the radiation field, the energy spectrum, the angular distribution, and the multiplicity of each type of secondary particles must also be known as a function of depth. Reactions can take place anywhere in a thick absorber; therefore, it is necessary to know these quantities as a function of particle energy for all particles produced. HZE transport methods are used to predict the radiation field; they are dependent on models of the interaction of man-made systems with the space environment to an even greater extent than methods used for other types of radiation. Hence, there is a major need to validate these transport codes by comparison with experimental data. The most cost-effective method of validation is a comparison with ground-based experimental measurements. A research program to provide such validation measurements using neon, iron and other accelerated heavy ion beams will be discussed and illustrated using results from ongoing experiments and their comparison with current transport codes. The extent to which physical measurements yield radiobiological predictions will be discussed.

NASA Discipline Radiation Health↗

NASA NTRS2010

Black Hole Mergers, Gravitational Waves, and Multi-Messenger Astronomy

The final merger of two black holes is expected to be the strongest source of gravitational waves for both ground-based detectors such as LIGO and VIRGO, as well as the space-based LISA. Since the merger takes place in the regime of strong dynamical gravity, computing the resulting gravitational waveforms requires solving the full Einstein equations of general relativity on a computer. Although numerical codes designed to simulate black hole mergers were plagued for many years by a host of instabilities, recent breakthroughs have conquered these problems and opened up this field dramatically. This talk will focus on the resulting gold rush of new results that is revealing the dynamics and waveforms of binary black hole mergers, and their applications in gravitational wave detection, astrophysics, and testing general relativity.

Centrella, Joan M.↗

DOE OSTI2022

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↗

NASA NTRS1988

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.↗

NASA NTRS1978

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.↗

NASA NTRS1988

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.↗

NASA NTRS2000

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.↗

DOE OSTI2022

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↗

DOE OSTI2021

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↗

NASA NTRS2000

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↗

NASA NTRS1980

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.↗

DOE OSTI2021

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↗

NASA NTRS2010

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↗

DOE OSTI2024

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↗