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,351 records · Page 75

Numerical results on the transcendence of constants involving pi, e, and Euler's constant

The existence of simple polynomial equations (integer relations) for the constants e/pi, e + pi, log pi, gamma (Euler's constant), e exp gamma, gamma/e, gamma/pi, and log gamma is investigated by means of numerical computations. The recursive form of the Ferguson-Fourcade algorithm (Ferguson and Fourcade, 1979; Ferguson, 1986 and 1987) is implemented on the Cray-2 supercomputer at NASA Ames, applying multiprecision techniques similar to those described by Bailey (1988) except that FFTs are used instead of dual-prime-modulus transforms for multiplication. It is shown that none of the constants has an integer relation of degree eight or less with coefficients of Euclidean norm 10 to the 9th or less.

Bailey, David H.↗

Flame balls - Past, present and future

This paper discusses analytical and numerical work that has been carried out in order to understand flame-balls and related phenomena that have been observed in microgravity experiments. The importance of heat losses is identified, whether from conduction, convection, or radiation. Accurate numerical simulations for hydrogen-air mixtures with radiation losses reveal a flammability limit of 3.5 percent of hydrogen by volume, a value close to the experimental one. The important role of stability analyses is emphasized, with particular attention to the role of three-dimensional instabilities in explaining unsteady spheroidal flames and flame-strings, objects that are observed in the experiments. We speculate that the dynamics of flame-strings is affected by in-depth radiation absorption for mixtures containing SF6, and report on some preliminary calculations in which this phenomenon is accounted for.

Buckmaster, J.↗

Evaluation of a pulse control law for flexible spacecraft

The following analytical and experimental studies were conducted: (1) A simple algorithm was developed to suppress the structural vibrations of 3-dimensional distributed parameter systems, subjected to interface motion and/or directly applied forces. The algorithm is designed to cope with structural oscillations superposed on top of rigid-body motion: a situation identical to that encountered by the SCOLE components. A significant feature of the method is that only local measurements of the structural displacements and velocities relative to the moving frame of reference are needed. (2) A numerical simulation study was conducted on a simple linear finite element model of a cantilevered plate which was subjected to test excitations consisting of impulsive base motion and of nonstationary wide-band random excitation applied at its root. In each situation, the aim was to suppress the vibrations of the plate relative to the moving base. (3) A small mechanical model resembling an aircraft wing was designed and fabricated to investigate the control algorithm under realistic laboratory conditions.

Source record↗

MEANDERING SHALLOW ATMOSPHERIC JET AS A MODEL OF SATURN'S NORTH-POLAR HEXAGON

The Voyager flybys of Saturn in 1980–1981 revealed a circumpolar Hexagon at ∼78° north planetographic latitude that has persisted for over 30 Earth years, more than one Saturn year, and has been observed by ground-based telescopes, Hubble Space Telescope and multiple instruments on board the Cassini orbiter. Its average phase speed is very slow with respect to the System III rotation rate, defined by the primary periodicity in the Saturn Kilometric Radiation during the Voyager era. Cloud tracking wind measurements reveal the presence of a prograde jet-stream whose path traces the Hexagon’s shape. Previous numerical models have produced large-amplitude, n = 6, wavy structures with westward intrinsic phase propagation (relative to the jet). However, the observed net phase speed has proven to be more difficult to achieve. Here we present numerical simulations showing that instabilities in shallow jets can equilibrate as meanders closely resembling the observed morphology and phase speed of Saturn’s northern Hexagon. We also find that the winds at the bottom of the model are as important as the winds at the cloud level in matching the observed Hexagon’s characteristics.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamic Pressure Drag on Rising Buoyant Thermals in a Neutrally Stable Environment

Abstract This study examines dynamic pressure drag on rising dry buoyant thermals. A theoretical expression for drag coefficient C d as a function of several other nondimensional parameters governing thermal dynamics is derived based on combining the thermal momentum budget with the similarity theory of Scorer. Using values for these nondimensional parameters from previous studies, the theory suggests drag on thermals is small relative to that on solid spheres in laminar or turbulent flow. Two sets of numerical simulations of thermals in an unstratified, neutrally stable environment using an LES configuration of the Cloud Model 1 (CM1) are analyzed. One set has a relatively low effective Reynolds number R e and the other has an order-of-magnitude-higher R e ; these produce laminar and turbulent resolved-scale flows, respectively. Consistent with the theoretical C d , the magnitude of drag is small in all simulations. However, whereas the laminar thermals have C d ≈ 0.01, the turbulent thermals have weakly negative drag ( C d ≈ −0.1). This difference is explained by the laminar thermals having near vertical symmetry but the turbulent thermals exhibiting considerable vertical asymmetry of their azimuthally averaged flows. In the laminar thermals, buoyancy rapidly becomes concentrated around the main centers of rotation located along the horizontal central axis, leading to expansion of thermals via baroclinic vorticity generation but doing little to break vertical symmetry of the flow. Vertical asymmetry of the azimuthally averaged flow of turbulent thermals is attributed mainly to small-scale resolved eddies that are concentrated in the upper part of the thermals.

Meteorology & Atmospheric Sciences↗

A numerical and experimental study on weld lines formation and strength in extrusion

In this work, we present a study on numerical and experimental approaches aiming to improve our understanding about the relation between the flow conditions and the strength of weld lines in extruded profiles. For this purpose, a prototype extrusion die that comprises a movable spider leg is used to produce tapes under different flow conditions, by varying the location of the spider leg and the flow rate. The same case studies are numerically studied, considering a viscoelastic Giesekus constitutive model. The numerical results are then investigated to identify potential relationships between the flow conditions, in which the weld lines are formed, with the strength of the weld line experimentally characterized. The results obtained allow identifying a clear correlation between the flow conditions and the weld‐line strength. POLYM. ENG. SCI., 58:249–260, 2018. © 2017 Society of Plastics Engineers

Ferrás, Luís L.↗

The Relative Contribution to Heavy Metals Production from Binary Neutron Star Mergers and Neutron Star–Black Hole Mergers

The origin of the heavy elements in the universe is not fully determined. Neutron star–black hole (NSBH) and binary neutron star (BNS) mergers may both produce heavy elements via rapid neutron-capture (r-process). We use the recent detection of gravitational waves from NSBHs, improved measurements of the neutron star equation of state (EoS), and the most modern numerical simulations of ejected material from binary collisions to measure the relative contribution of NSBHs and BNSs to the production of heavy elements. As the amount of r-process ejecta depends on the mass and spin distribution of the compact objects, as well as on the EoS of the neutron stars, we consider various models for these quantities, informed by gravitational-wave and pulsar data. We find that in most scenarios, BNSs have produced more r-process elements than NSBHs over the past 2.5 billion years. If black holes have preferentially small spins, BNSs can produce at least twice the amount of r-process elements than NSBHs. If black hole spins are small and there is a dearth of low-mass (<5M ⊙ ) black holes within NSBH binaries, BNSs can account for the near totality of the r-process elements from binaries. For NSBH to produce a large fraction of r-process elements, black holes in NSBHs must have small masses and large aligned spins, which is disfavored by current data.

79 ASTRONOMY AND ASTROPHYSICS↗

Linear and nonlinear acoustic wave propagation in the atmosphere

The investigation of the acoustic wave propagation theory and numerical implementation for the situation of an isothermal atmosphere is described. A one-dimensional model to validate an asymptotic theory and a 3-D situation to relate to a realistic situation are considered. In addition, nonlinear wave propagation and the numerical treatment are included. It is known that the gravitational effects play a crucial role in the low frequency acoustic wave propagation. They propagate large distances and, as such, the numerical treatment of those problems become difficult in terms of posing boundary conditions which are valid for all frequencies.

Hariharan, S. I.↗

Low Reynolds number kappa-epsilon modeling with the aid of direct simulation data

The constant C(sub mu) and the near-wall damping function f(sub mu) in the eddy-viscosity relation of the kappa-epsilon model are evaluated from direct numerical simulation (DNS) data for developed channel and boundary layer flow at two Reynolds numbers each. Various existing f(sub mu) model functions are compared with the DNS data, and a new function is fitted to the high-Reynolds-number channel flow data. The epsilon-budget is computed for the fully developed channel flow. The relative magnitude of the terms in the epsilon-equation is analyzed with the aid of scaling arguments, and the parameter governing this magnitude is established. Models for the sum of all source and sink terms in the epsilon-equation are tested against the DNS data, and an improved model is proposed.

Rodi, W.↗

Low Reynolds number k-epsilon modelling with the aid of direct simulation data

The constant C sub mu and the near-wall damping function f sub mu in the eddy-viscosity relation of the k-epsilon model are evaluated from direct numerical simulation (DNS) data for developed channel and boundary layer flow at two Reynolds numbers each. Various existing f sub mu model functions are compared with the DNS data, and a new function is fitted to the high-Reynolds-number channel flow data. The epsilon-budget is computed for the fully developed channel flow. The relative magnitude of the terms in the epsilon-equation is analyzed with the aid of scaling arguments, and the parameter governing this magnitude is established. Models for the sum of all source and sink terms in the epsilon-equation are tested against the DNS data, and an improved model is proposed.

Rodi, W.↗

Small-scale fluctuations of the X-ray background

Fluctuations of the diffuse X-ray background are investigated using the Einstein Observatory Deep Survey data. Small angular-scale variations of the numbers of counts in the energy range 0.8-3.5 keV in the central 32 x 32 sq arcmin of the IPC are analyzed by means of the two-point correlation function. It is shown that detector modulation contributes substantially to the observed fluctuations. Instrumental effects are estimated and subtracted from the raw observations. The residual fluctuations are used to estimate the log N-log S relation of the X-ray sources below the Einstein sensitivity limit. Taking into account numerous uncertainties involved in this analysis, it is shown that the amplitude of fluctuations observed in the Deep Survey data is consistent with a relatively small number of sources producing the X-ray background.

Soltan, A. M.↗

Advanced Numerical Model for Irradiated Concrete

In this report, we establish a numerical model for concrete exposed to irradiation to address these three critical points. The model accounts for creep in the cement paste and its coupling with damage, temperature and relative humidity. The shift in failure mode with the loading rate is also properly represented. The numerical model for creep has been validated and calibrated against different experiments in the literature [Wittmann, 1970, Le Roy, 1995]. Results from a simplified model are shown to showcase the ability of numerical homogenization to simulate irradiation effects in concrete. In future works, the complete model will be applied to the analysis of the irradiation experiments of Elleuch et al. [1972] and Kelly et al. [1969]. This requires a careful examination of the experimental environmental conditions as in both cases certain critical information are missing, including the relative humidity history. A sensitivity analysis will be conducted to provide lower and upper bounds of the concrete expansion under irradiation, and check if the scatter in the simulated results matches the one found in experiments. The numerical and experimental results will be compared in terms of expansion and loss of mechanical stiffness and strength. Both effects should be captured accordingly by the model to validate it. Once the model has been validated on these two experiments, it can be applied to simulate concrete from nuclear power plants. To do so, the materials used in these concrete must be as well characterized as possible. The main parameters required are the mechanical properties of each constituent in the concrete (aggregates, cement paste), namely the elastic modulus, the creep properties, the tensile and compressive strength, the thermal expansion coefficient, and the drying shrinkage. These can be either measured experimentally, estimated from the initial composition in the case of cement paste, or back-calculated from mechanical tests on concrete. If some are unknown, a sensitivity analysis must be carried out to provide lower and upper bounds of the material behaviour. Finally, the model can be used as a basis to formulate a macroscopic material model for concrete subject to irradiation, which later can be used in structural analyses to estimate the structural impact of irradiation on nuclear power plants.

36 MATERIALS SCIENCE↗

Characterizing Mesoscale Cellular Convection in Marine Cold Air Outbreaks With a Machine Learning Approach

Abstract During marine cold‐air outbreaks (MCAOs), when cold polar air moves over warmer ocean, a well‐recognized cloud pattern develops, with open or closed mesoscale cellular convection (MCC) at larger fetch over open water. The Cold‐Air Outbreaks in the Marine Boundary Layer Experiment provided a comprehensive set of ground‐based in situ and remote sensing observations of MCAOs at a coastal location in northern Norway. MCAO periods that unambiguously exhibit open or closed MCC are determined. Individual cells observed with a profiling Ka‐band radar are identified using a watershed segmentation method. Using self‐organizing maps (SOMs), these cells are then objectively classified based on the variability in their vertical structure. The SOM nodes contain some information about the location of the cell transect relative to the center of the MCC. This adds classification noise, requiring numerous cell transects to isolate cell dynamical information. The SOM‐based classification shows that comparatively intense convection occurs only in open MCC. This convection undergoes an apparent lifecycle. Developing cells are associated with stronger updrafts, large spectrum width, larger amounts of liquid water, lower surface precipitation rates, and lower cloud tops than mature and weakening cells. The weakening of these cells is associated with the development of precipitation‐induced cold pools. The SOM classification also reveals less intense convection, with a similar lifecycle. More stratiform vertical cloud structures with weak vertical motions are common during closed MCC periods and are separated into precipitating and non‐precipitating stratiform cores. Convection is observed only occasionally in the closed MCC environment.

Meteorology & Atmospheric Sciences↗

On the combustion of n-butyl acetate droplets

This paper presents an experimental and numerical study of the combustion of isolated n‑butyl acetate droplets in the standard atmosphere. Numerical simulations are reported using a model that incorporates unsteady gas and liquid transport, variable properties, and radiation. Three skeletal mechanisms of n‑butyl acetate, derived from a large detailed mechanism comprised of 819 species and 52,698 reactions, were used in the numerical simulations to evaluate the influence of the kinetic mechanism on burning. The reduced mechanisms comprised 212 species and 5413 reactions, 157 species and 3089 reactions, and 105 species and 1035 reactions. The numerical model did not include soot formation, though qualitatively mild sooting was noted only for droplets larger than 0.7 mm. The numerical predictions were in good agreement with experimental measurements of droplet and flame diameters. Flame extinction was numerically predicted which was attributed to a decrease of the characteristic diffusion time relative to the chemical time as droplet burned. Effects of initial droplet diameter on the evolution of maximum gas temperature (Tmax) and peak mole fractions of CO2 and CO are also examined numerically.

09 BIOMASS FUELS↗

Numerical methods for the simulation of complex multi-body flows with applications for the integrated Space Shuttle vehicle

This project forms part of the long term computational effort to simulate the time dependent flow over the integrated Space Shuttle vehicle (orbiter, solid rocket boosters (SRB's), external tank (ET), and attach hardware) during its ascent mode for various nominal and abort flight conditions. Due to the limitations of experimental data such as wind tunnel wall effects and the difficulty of safely obtaining valid flight data, numerical simulations are undertaken to supplement the existing data base. This data can then be used to predict the aerodynamic behavior over a wide range of flight conditions. Existing computational results show relatively good overall comparison with experiments but further refinement is required to reduce numerical errors and to obtain finer agreements over a larger parameter space. One of the important goals of this project is to obtain better comparisons between numerical simulations and experiments. In the simulations performed so far, the geometry has been simplified in various ways to reduce the complexity so that useful results can be obtained in a reasonable time frame due to limitations in computer resources. In this project, the finer details of the major components of the Space Shuttle are modeled better by including more complexity in the geometry definition. Smaller components not included in early Space Shuttle simulations will now be modeled and gridded.

Chan, William M.↗

Uncertainty Analysis of Decomposing Polyurethane Foam

Sensitivity/uncertainty analyses are necessary to determine where to allocate resources for improved predictions in support of our nation's nuclear safety mission. Yet, sensitivity/uncertainty analyses are not commonly performed on complex combustion models because the calculations are time consuming, CPU intensive, nontrivial exercises that can lead to deceptive results. To illustrate these ideas, a variety of sensitivity/uncertainty analyses were used to determine the uncertainty associated with thermal decomposition of polyurethane foam exposed to high radiative flux boundary conditions. The polyurethane used in this study is a rigid closed-cell foam used as an encapsulant. Related polyurethane binders such as Estane are used in many energetic materials of interest to the JANNAF community. The complex, finite element foam decomposition model used in this study has 25 input parameters that include chemistry, polymer structure, and thermophysical properties. The response variable was selected as the steady-state decomposition front velocity calculated as the derivative of the decomposition front location versus time. An analytical mean value sensitivity/uncertainty (MV) analysis was used to determine the standard deviation by taking numerical derivatives of the response variable with respect to each of the 25 input parameters. Since the response variable is also a derivative, the standard deviation was essentially determined from a second derivative that was extremely sensitive to numerical noise. To minimize the numerical noise, 50-micrometer element dimensions and approximately 1-msec time steps were required to obtain stable uncertainty results. As an alternative method to determine the uncertainty and sensitivity in the decomposition front velocity, surrogate response surfaces were generated for use with a constrained Latin Hypercube Sampling (LHS) technique. Two surrogate response surfaces were investigated: 1) a linear surrogate response surface (LIN) and 2) a quadratic response surface (QUAD). The LHS techniques do not require derivatives of the response variable and are subsequently relatively insensitive to numerical noise. To compare the LIN and QUAD methods to the MV method, a direct LHS analysis (DLHS) was performed using the full grid and timestep resolved finite element model. The surrogate response models (LIN and QUAD) are shown to give acceptable values of the mean and standard deviation when compared to the fully converged DLHS model.

Hobbs, Michael L.↗

Atmospheric Diabatic Heating in Different Weather States and the General Circulation

Analysis of multiple global satellite products identifies distinctive weather states of the atmosphere from the mesoscale pattern of cloud properties and quantifies the associated diabatic heating/cooling by radiative flux divergence, precipitation, and surface sensible heat flux. The results show that the forcing for the atmospheric general circulation is a very dynamic process, varying strongly at weather space-time scales, comprising relatively infrequent, strong heating events by ''stormy'' weather and more nearly continuous, weak cooling by ''fair'' weather. Such behavior undercuts the value of analyses of time-averaged energy exchanges in observations or numerical models. It is proposed that an analysis of the joint time-related variations of the global weather states and the general circulation on weather space-time scales might be used to establish useful ''feedback like'' relationships between cloud processes and the large-scale circulation.

Cooling↗