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 379 records · Page 21

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Revised Kimberlina Groundwater Simulations

Kimberlina Groundwater Qaulity Simulations (Version 3: 2015 October) by Lawrence Livermore National Labratory. Submitted to the EDX team at the National Energy Technology Laboratory for research use. This dataset contains 562 updated groundwater quality simulations based on the Kimberlina site. It is a six-layer model including the aquifers and aquitards inbetween the ground surface and the reservoir. The input files and associated leakage information for these simulations are included in the zip file "nrap_kim_well31_fluxes and inputs.tar.gz". These 562 simulations are the successful results out of 1000 monte carlo simulations, with time length greater than 5 years. Each simulation consists 6 output files, therefore, there are 3372 files (6*562) in total. One single example simulation of full time length (200 year) is also provided in the zip file "single_simulation_200yr_sim0012.tar.gz". These 6 groundwater quality parameters are: - pH (use absolute value) - TDS (mg/kg) - Pressure (Pa) - CO2 Liquid Saturation (fraction range 0.0-1.0) - CO2 Gas Concentration (mol/L) - CO2 Liquid Concentration (mol/kg) Their correponding file names, take the first simulation for example, are: ph0001.ntab (for pH) , tds0001.ntab (for TDS), usnt0001.p.ntab (for pressure), usnt0001.s_liquid.ntab (for liquid saturation), usnt0001.c_co2_gas.ntab (for CO2 gas concentration) and usnt0001.c_co2_liquid.ntab (for CO2 liquid concentration). Each output file contains at least 19 columns (variables). The total number of columns depends on the simulation time length. The column data include the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the jth grid of y-axis 4. k: the kth grid of z-axis 5. element_ref: element type reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The rest (14, 15, 16...) variables are the simulated parameter values over time, take pH as an example, are: 14. ph000y: initial pH value. 15. ph001y: simulated pH value at the end of the 1st year. 16. ph002y: simulated pH value at the end of the 2nd year. ... (and so on)... The zip input files have the input fluxes of brine and CO2 (in kg/s) allocated to model layers 1,2,3 and 5 (layers 4 and 6 are assumed to be no-flow aquitards) through the following nodes: (src1 WELL 7 7 25 25 10 32) #Etchegoin (src2 WELL 7 7 25 25 33 72) #Macoma Chanac (src3 WELL 7 7 25 25 73 83) #Santa Margarita Mclure (src4 WELL 7 7 25 25 84 89) #Fruitvale Round Mountain (aquitard) (src5 WELL 7 7 25 25 90 96) #Olcese Sand (thief zone) (src6 WELL 7 7 25 25 97 99) #tremblor Freeman (aquitard) ex: brn1 is assigned to src 1 from i=(7,7), j=(25,25), k=(10,32) Edited by Ya-Mei Yang at National Energy Technology Laboratory

brine leakage↗

DOE OSTI2021

Computational Fluid Dynamics (CFD) Simulations of Taylor Bubbles in Vertical and Inclined Pipes with Upward and Downward Liquid Flow

Summary Two-phase flow is a common occurrence in pipes of oil and gas developments. Current predictive tools are based on the mechanistic two-fluid model, which requires the use of closure relations to predict integral flow parameters such as liquid holdup (or void fraction) and pressure gradient. However, these closure relations carry the highest uncertainties in the model. In particular, significant discrepancies have been found between experimental data and closure relations for the Taylor bubble velocity in slug flow, which has been determined to strongly affect the mechanistic model predictions (Lizarraga-García 2016). In this work, we study the behavior of Taylor bubbles in vertical and inclined pipes with upward and downward flow using a validated 3D computational fluid dynamics (CFD) approach with level set method implemented in a commercial code. A total of 56 cases are simulated, covering a wide range of fluid properties, pipe diameters, and inclination angles: Eo ∈ [10, 700]; Mo ∈ [1×10–6, 5×103]; ReSL ∈ [–40, 10]; θ ∈ [5°, 90°]. For bubbles in vertical upward flows, the simulated distribution parameter, C0, is successfully compared with an existing model. However, the C0 values of downward and inclined slug flows where the bubble becomes asymmetric are shown to be significantly different from their respective vertical upward flow values, and no current model exists for the fluids simulated here. The main contributions of this work are (1) the relatively large 3D numerical database generated for this type of flow, (2) the study of the asymmetric nature of inclined and some vertical downward slug flows, and (3) the analysis of its impact on the distribution parameter, C0.

Engineering↗

DOE OSTI2022

Witness the Trinity test through Lab artifacts

July 16 marks the 77th anniversary of the Trinity test, conducted in a desert in New Mexico. The test subject, an atomic bomb called The Gadget, was successfully detonated from a 100-foot steel tower. This event marks the commencement of the Atomic Age, a new era where fission capabilities could be employed for national security purposes. Shortly after the Trinity test, two Los Alamos-created atomic weapons were released above Japan, helping to end the world’s bloodiest conflict just weeks later. “Trinity was one of the greatest scientific experiments ever,” said NSRC Senior Historian Alan Carr said. “Los Alamos scientists changed the world forever on that day. Not only was it the dawn of the Atomic Age, but also the beginning of the Lab’s eight decades of cutting-edge science and its national security charge.” To preserve this event, and to continue to learn more about this critical moment in history, the National Security Research Center (NSRC) curates a collection of photographs, films, notes, unclassified artifacts and numerous other materials related to the science of the test. Notably, the collection includes a novel material that formed at the site, trinitite, and artifacts from one of the intriguing scientists present at the test, Enrico Fermi.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

NASA NTRS1971

Longitudinal waves in a perpendicular collisionless plasma shock. III.

This paper considers electrostatic waves in a Vlasov plasma of unmagnetized ions and magnetized electrons undergoing an E x B drift. The linear dispersion relation is solved numerically for an electron temperature approximately equal to the ion temperature. For a fixed ratio of drift velocity to electron thermal velocity, the growth rates of the E x B electron drift instability are smaller, and the waves are stabilized at much smaller values of k.B than in the case where electron temperature is much larger than the ion temperature.

Gary, S. P.↗

NASA NTRS1973

Taylor stability of viscous fluids with application to film boiling

The dispersion relation is evaluated numerically for Taylor waves in a viscous unstable interface with surface tension. The solution takes account of transverse curvature and the numerical evaluations apply to horizontal cylindrical, as well as to plane, interfaces. The result is verified with frequency and wavelength data obtained during film boiling on horizontal wires. A very general empirical correlation is given, en passant, for the vapor blanket thickness during film boiling.

Dhir, V. K.↗

NASA NTRS1977

Effects of finite plasma beta on the lower-hybrid-drift instability

The local dispersion relation for the lower-hybrid-drift (LHD) instability is derived and analyzed, taking into account the finite-beta effects associated with transverse electromagnetic perturbations as well as with resonant and nonresonant electron-orbit modifications due to magnetic-field gradients. The influence of finite-beta effects on the LHD instability is calculated in a fully self-consistent manner for arbitrary values of electron-ion temperature ratio, local beta, cross-field ExB velocity/ion thermal speed ratio, and other plasma parameters. Stability properties are investigated analytically for the case of cold electrons, and the local dispersion relation is solved numerically in the parameter regime of most interest for high-density plasma pinches. The results show that for all parameter regimes studied, the net effect of finite plasma beta is to reduce the maximum growth rate of the LHD instability, although the details can vary, depending on the plasma parameters. Except in the limit where the electron/ion temperature ratio tends to zero, it is found that there is a critical value of plasma beta above which the LHD instability is completely stabilized.

Davidson, R. C.↗

NASA NTRS1980

A generalized vortex lattice method

Several variations of vortex lattice method that are currently available have proved practical and versatile theoretical tools for aerodynamic analysis and design of planar and nonplanar configurations. Success of method is due in great part to relative simplicity of numerical technique involved and to accuracy of results obtained; however, most of available procedures are for subsonic flow applications. VORLAX program was developed to incorporate direct extension of vortex lattice method into supersonic flow regime, thus providing analyst with full flow range capability.

Baker, W. M.↗

NASA NTRS1980

Stability aspects of diverging subsonic flow

The paper describes the computation of two-dimensional, subsonic, diverging internal flows and how they differ from the corresponding converging flows. Such diverging or decelerating flows occur in such obvious places as subsonic diffusers and inlets; however, such flows also occur in supersonic nozzles in the presence of a normal shock. The flow instability and its relation to the numerical method used, boundary conditions, and viscous effects are assessed both analytically and numerically. The inviscid flow is shown to be physically unstable and a poor representation of the true viscous flow.

Cline, M. C.↗

NASA NTRS1983

Further observations on resonance cones in non-Maxwellian plasmas

Results on the angular distribution of the electrostatic potential of a pulsating point charge in a warm magnetized plasma permeated by an electron beam are presented. The theoretical formulation for a finite magnetic field is given, and the solution of the resonance cone dispersion relation is presented. Numerical results on the angular distribution of the potential are shown, and the propagation of waves outside the resonance cones is described. It is demonstrated that with the inclusions of a finite magnetic field, the field patterns of a point charge are qualitatively similar to those obtained for a uniaxial plasma. The Cerenkov radiation occurs at angles much smaller than the cold-cone angle, even with the finite magnetic field. When the beam velocity is well above the thermal velocity of the background electrons, a characteristic wave propagation occurs between the cold-cone angles.

Thiemann, H.↗

NASA NTRS1984

1983 lightning, turbulence, wind shear, and Doppler radar studies at the National Severe Storms Laboratory

As part of continuing research on aviation related weather hazards, numerous experiments were incorporated into the 1983 Spring Observation Program. This year's program was an abbreviated one because of commitments made to the development of the Next Generation Radar (NEXRAD) project. The National Oceanic and Atmospheric Administration's (NOAA) P-3 Orion and the National Aeronautics and Space Administration's (NASA) RB-57B and U-2 were the main aircraft involved in the studies of lightning, wind shear, turbulence, and storm structure. A total of 14 flights were made by these aircraft during the period of May 16 through June 5, 1983. Aircraft instrumentation experiments are described, and resultant data sets available for research are detailed. Aircraft instrumentation and Doppler radar characteristics are detailed.

Lee, J. T.↗

NASA NTRS1986

Theoretical Foundation for Weld Modeling

Differential equations describe physics of tungsten/inert-gas and plasma-arc welding in aluminum. Report collects and describes necessary theoretical foundation upon which numerical welding model is constructed for tungsten/inert gas or plasma-arc welding in aluminum without keyhole. Governing partial differential equations for flow of heat, metal, and current given, together with boundary conditions relevant to welding process. Numerical estimates for relative importance of various phenomena and required properties of 2219 aluminum included

Traugott, S.↗

NASA NTRS1985

Theory of the current-driven ion cyclotron instability in the bottomside ionosphere

A comprehensive treatment of the theory of electrostatic ion cyclotron instability in the collisional bottomside ionosphere is given. The linear dispersion relations for the ion cyclotron instability (electron and ion collisions) are derived. Analytical solutions to the dispersion relation are obtained in three limits: collisionless; weakly collisional; and collisional. In addition, the linear dispersion relation is solved numerically, using parameters typical of the high-latitude bottomside ionosphere. It is pointed out that the linear growth rates in the bottomside collisional regime can be of the same order as the corresponding topside collisionless growth rates. A nonlinear saturation mechanism for the electrostatic ion cyclotron instability is proposed on the basis of Dupree's (1966) resonance broadening theory.

Satyanarayana, P.↗

NASA NTRS1986

Conforming versus non-conforming boundary elements in three-dimensional elastostatics

A critical comparison of two basic formulations in three-dimensional elastostatics, using conforming and nonconforming boundary elements, is presented. The basic structure of the boundary element method is developed. The peculiarities that both types of boundary elements present in relation to the numerical implementation are discussed. Through selected examples, key issues such as the computational advantages and disadvantages of both formulations, mesh discretization and accuracy questions, and optimal location of the collocation nodes in the case of nonconforming elements are addressed. It is shown that conforming elements are able to produce more accurate results than nonconforming ones, with substantial economy in the final size of the system equations.

Manolis, G. D.↗

NASA NTRS1987

Aspects of supercritical wing theory

Concepts and techniques that have proved successful in the computational theory of shockless airfoils are examined. Shockless airfoils have been used in the design of supercritical wings for fuel-efficient aircraft. In transonic aerodynamics, it is has been found that test data agree well with the best results obtained from CFD. This is particularly true for two-dimensional flow, where accurate experiments are difficult to perform and numerical calculations are relatively easy. However, it is important to choose reasonable models in setting up the computations. Particular attention is given to the formulation of the relevant differential equations in ways that have been found to be satisfactory in practice.

Garabedian, P. R.↗

NASA NTRS1988

USSR Space Life Sciences Digest, issue 20

Abstracts of research in the areas of biological rhythms, body fluids, botany, endrocrinology, enzymology, exobiology, genetics, human performance, immunology, life support systems, mathematical modeling, and numerous other topics related to space and life sciences are given.

Hooke, Lydia Razran↗

NASA NTRS1989

Discrete events and solar wind energization

Based on a multiple-magnetic-reconnection picture, an estimation of the energy flux suggests that small-scale EUV exploding events may contribute a significant amount of energy (of order of 100,000 erg/sq cm sec) to solar atmospheric heating and solar-wind acceleration. Most of the dissipated magnetic energy is converted into thermal energy and plasma turbulence. On a related aspect, a numerical study based on the nonlinear one-fluid hydrodynamic equations shows a self-smoothing effect, whereby a multistream structure of the solar wind formed near the sun can be gradually smoothed during its propagation through interplanetary space. This calculation gives support for the possible contribution of discrete energetic events to high-speed solar wind streams.

Yang, W.-H.↗

NASA NTRS1989

A current disruption mechanism in the neutral sheet for triggering substorm expansions

Two main areas were addressed in support of an effort to understand mechanism responsible for the broadband electrostatic noise (BEN) observed in the magnetotail. The first area concerns the generation of BEN in the boundary layer region of the magnetotail whereas the second area concerns the occassional presence of BEN in the neutral sheet region. For the generation of BEN in the boundary layer region, a hybrid simulation code was developed to perform reliable longtime, quiet, highly resolved simulations of field aligned electron and ion beam flow. The result of the simulation shows that broadband emissions cannot be generated by beam-plasma instability if realistic values of the ion beam parameters are used. The waves generated from beam-plasma instability are highly discrete and are of high frequencies. For the plasma sheet boundary layer condition, the wave frequencies are in the kHz range, which is incompatible with the observation that the peak power in BEN occur in the 10's of Hz range. It was found that the BEN characteristics are more consistent with lower hybrid drift instability. For the occasional presence of BEN in the neutral sheet region, a linear analysis of the kinetic cross-field streaming instability appropriate to the neutral sheet condition just prior to onset of substorm expansion was performed. By solving numerically the dispersion relation, it was found that the instability has a growth time comparable to the onset time scale of substorm onset. The excited waves have a mixed polarization in the lower hybrid frequency range. The imposed drift driving the instability corresponds to unmagnetized ions undergoing current sheet acceleration in the presence of a cross-tail electric field. The required electric field strength is in the 10 mV/m range which is well within the observed electric field values detected in the neutral sheet during substorms. This finding can potentially account for the disruption of cross-tail current and its diversion to the ionosphere to form the substorm current wedge. Furthermore, a number of features associated with substorm expansion onset can be understood based on this substorm onset scenario.

Lui, A. T. Y.↗

NASA NTRS1990

The nondeterministic divide

The nondeterministic divide partitions a vector into two non-empty slices by allowing the point of division to be chosen nondeterministically. Support for high-level divide-and-conquer programming provided by the nondeterministic divide is investigated. A diva algorithm is a recursive divide-and-conquer sequential algorithm on one or more vectors of the same range, whose division point for a new pair of recursive calls is chosen nondeterministically before any computation is performed and whose recursive calls are made immediately after the choice of division point; also, access to vector components is only permitted during activations in which the vector parameters have unit length. The notion of diva algorithm is formulated precisely as a diva call, a restricted call on a sequential procedure. Diva calls are proven to be intimately related to associativity. Numerous applications of diva calls are given and strategies are described for translating a diva call into code for a variety of parallel computers. Thus diva algorithms separate logical correctness concerns from implementation concerns.

Charlesworth, Arthur↗