Engineering Papers⌕ Search

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.

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

At least 433 records · Page 24

Nonstandard Finite Difference Schemes: Relations Between Time and Space Step-Sizes in Numerical Schemes for PDE's That Follow from Positivity Condition

A large class of physical phenomena can be modeled by evolution and wave type Partial Differential Equations (PDE). Few of these equations have known explicit exact solutions. Finite-difference techniques are a popular method for constructing discrete representations of these equations for the purpose of numerical integration. However, the solutions to the difference equations often contain so called numerical instabilities; these are solutions to the difference equations that do not correspond to any solution of the PDE's. For explicit schemes, the elimination of this behavior requires functional relations to exist between the time and space steps-sizes. We show that such functional relations can be obtained for certain PDE's by use of a positivity condition. The PDE's studied are the Burgers, Fisher, and linearized Euler equations.

Mickens, Ronald E.↗

Sound diffraction around screens and wedges for arbitrary point source locations

A development is presented for the Green's function for a point source in the vicinity of a rigid wedge. The diffraction contributions to the Green's function for arbitrary source and listener location is expressed in a form which can be readily evaluated using the Laguerre technique for numerical integration. The present approach offers the advantages of efficient numerical evaluation and of relatively straightforward reduction to well-known analytical approximations in limiting cases. Comparisons with previously obtained experimental and numerical results obtained by Ambaud and Bergassoli (1972) are presented. The comparison with the experimental results is excellent; the advantages of the present numerical technique, vis a vis that of Ambaud and Bergassoli, are pointed out.

Hadden, W. J., Jr.↗

An Overview of the Null-Field Method. II: Convergence and Numerical Stability

In this paper we provide an analysis of the convergence and numerical stability of the null-field method with discrete sources. We show that (i) if the null-field scheme is numerically stable then we can decide whether or not convergence can be achieved; (ii) if the null-field scheme is numerically unstable then we cannot draw any conclusion about the convergence issue; and (iii) the numerical stability is closely related to the property of a tangential system of radiating discrete sources to form a Riesz basis. Our numerical analysis indicates that for prolate spheroids and localized vector spherical wave functions, the null-field scheme is numerically unstable (this system of vector functions does not form a Riesz basis), while for distributed vector spherical wave functions, the numerical instability is not so pronounced (this system of discrete sources almost possesses the property of being a Riesz basis). We also describe an analytical method for computing the surface integrals in the framework of the conventional null-field method with localized vector spherical wave functions which increases the stability of the numerical scheme.

Null-field method↗

Simulation of Propellant Loading System Senior Design Implement in Computer Algorithm

Propellant loading from the Storage Tank to the External Tank is one of the very important and time consuming pre-launch ground operations for the launch vehicle. The propellant loading system is a complex integrated system involving many physical components such as the storage tank filled with cryogenic fluid at a very low temperature, the long pipe line connecting the storage tank with the external tank, the external tank along with the flare stack, and vent systems for releasing the excess fuel. Some of the very important parameters useful for design purpose are the prediction of pre-chill time, loading time, amount of fuel lost, the maximum pressure rise etc. The physics involved for mathematical modeling is quite complex due to the fact the process is unsteady, there is phase change as some of the fuel changes from liquid to gas state, then conjugate heat transfer in the pipe walls as well as between solid-to-fluid region. The simulation is very tedious and time consuming too. So overall, this is a complex system and the objective of the work is student's involvement and work in the parametric study and optimization of numerical modeling towards the design of such system. The students have to first become familiar and understand the physical process, the related mathematics and the numerical algorithm. The work involves exploring (i) improved algorithm to make the transient simulation computationally effective (reduced CPU time) and (ii) Parametric study to evaluate design parameters by changing the operational conditions

Bandyopadhyay, Alak↗

Nonlinear evolution of quadratic gravity in 3 + 1 dimensions

We present a numerically stable system of ( 3 + 1 ) evolution equations for the nonlinear gravitational dynamics of quadratic-curvature corrections to General Relativity (quadratic gravity). We also report on the numerical implementation of these evolution equations. We recover a well-known linear instability and gather evidence that—aside from said instability—quadratic gravity exhibits a physically stable Ricci-flat subsector. In particular, we demonstrate that Teukolsky-wave perturbations of a Schwarzschild black hole as well as a full binary inspiral (evolved up to merger) remain Ricci flat throughout evolution. Furthermore, this suggests that, at least in vacuum, classical quadratic gravity can mimic General Relativity, even in the fully nonlinear strong-gravity regime.

79 ASTRONOMY AND ASTROPHYSICS↗

Morrow1 Model - No Capillary Pressure

Farnsworth Unit (FWU) CO2EOR Eclipse compositional model: Morrow1 model uses the Morrow1 relative permeability curve with no capillary pressure assigned.

capillary pressure↗

Morrow1 Model - with Capillary Pressure

Farnsworth Unit (FWU) CO2EOR Eclipse compositional model: Morrow1 model uses the Morrow1 relative permeability curve with capillary pressure assigned homogeneously

Capillary Pressure↗

MATBOX, an Open-Source Microstructure Analysis Toolbox for Meshing, Generation, Segmentation, and Characterization of 3D Heterogenous Volumes

Battery performance is strongly correlated with electrode microstructural properties. To account for its impact, lithium-ion battery (LIB) models either abstract the microstructural heterogeneity of composite electrodes using effective macroscopic properties (macro- or meso- scale models) or directly solve the system of equations on the microstructure geometry or mesh (microstructure-scale models). Therefore, to be adequate, both families of models require information from the microstructure geometry, which can be provided by the numerical tool presented in this work. MATBOX is a MATLAB open-source application [1] developed by NREL for performing various microstructure-related tasks including microstructure numerical generation, image filtering and microstructure segmentation, microstructure characterization and correlation, visualization, and microstructure meshing. MATBOX was originally developed for the analysis of LIB electrode microstructures; however, the algorithms provided by the toolbox are widely applicable to other heterogeneous materials. The toolbox provides a user-friendly experience thanks to a Graphical-User Interface, requires no coding by the user, and is well documented. This presentation will illustrate various MATBOX features for the characterization of a LIB electrode, including a fully automated Representative Volume Element (RVE) analysis, the numerical generation of complex 'virtual' microstructure, including dual-layer electrodes and carbon-binder additive phase, and the meshing of a complex NMC/graphite full cell microstructure suitable for 3D finite-element modeling. Other modules (segmentation, visualization, and correlation) will be briefly presented. Thanks to its modular, open-source approach, MATBOX can easily incorporate third-party algorithms to eventually build a standard in the field that will benefit the whole scientific community. Effective diffusion coefficient [2], additive phase numerical generation [3], and meshing [4] third-party algorithms have been already integrated in the toolbox with more to come.

DIRECT ENERGY CONVERSION,MATHEMATICS AND COMPUTING↗

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗

Capillary Pressure to Relative Permeability: Task 7.3.4 Relative Permeability Analysis Final Report

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability relationship can be determined. Relative permeability is a numerical description of the interplay between two or more fluids and the porous media they flow through. It is a critical parameter for various tools used to characterized subsurface multiphase flow systems, such as numerical simulation for oil and gas development, carbon sequestration, and groundwater contamination remediation. Therefore, it is critical to get a good statistic distribution of relative permeability in the porous media under study. Empirical relationships for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility required to match laboratory derive relative permeability relationships. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability relationships. Existing laboratory-measured relative permeability data along with mercury intrusion capillary (MICP) data coupled with interfacial tension and contact angle measurements were used to determine the efficacy of this approach to relative permeability curve creation. The relative permeability relationships determined with this method were fit to existing laboratory data to elucidate common fitting parameters that were then used to create relative permeability relationships from MICP data that does not have an associated laboratory-measured relative permeability relationship.

58 GEOSCIENCES↗

Numerical formulation for a higher order plane finite dynamic element

The paper describes the development of an eight-node plane rectangular finite dynamic element and presents detailed descriptions of the associated numerical formulation involving the higher order dynamic correction terms pertaining to the related stiffness and inertia matrices. Numerical test results of free vibration analyses are presented in detail for the newly developed eight-node element and also the corresponding four-node element in order to make a clear comparison of the relative efficiencies of the corresponding finite element and dynamic element procedures. Such results indicate a superior pattern of solution convergence of the presently developed dynamic element.

Gupta, K. K.↗

A Review and Reflections on the Sun-Climate Connection

The field of sun-climate is beset with an extraordinary number of numerical correlations attempting to relate various periodicities of solar activity with changes in the Earth's weather and climate. Signatures representing climatological variability have been sought for cycles as short as the solar 28-day rotational period up to Milankovich periods of thousands of years, although a majority of correlations have concentrated on the 11-year sunspot and 22-year Hale double sunspot cycles. For the shorter term, parameters including temperature, pressure, winds storm tracks, rainfall, and water levels in rivers and lakes, etc. have been correlated with solar variability. For longer periods, it has been necessary to seek more indirect evidence in ice cores, tree rings, and geologic deep sea cores. Other atmospheric parameters relating to atmospheric electricity and the global electric circuit have also been correlated in similar fashion. Unfortunately, few, if any, of this wide spectrum of numerical correlations have been associated with any viable physical explanation, making most studies in the field an exercise in numerical statistics. More recently, a few suggestions for plausible coupling processes have begun to appear. These, coupled with new and stronger correlations involving selective binning of climatological data sets have injected new life and hope to this field. An overview is given of the historical past and current perspectives, to evaluate possible avenues for defining physical linking processes in the future.

Goldberg, Richard A.↗

Hydrostratigraphic Region 4 Model - with HSU3-4Pc

Farnsworth Unit (FWU) CO2EOR Eclipse compositional model: Hydrostratigraphic Region 4 model uses the Morrow2 (UNOCAL) relative permeability curve with HSU3-4 capillary pressure assigned homogeneously

Capillary Pressure↗

Hydrostratigraphic Region 4 Model - with HSUPc

Farnsworth Unit (FWU) CO2EOR Eclipse compositional model: Hydrostratigraphic Region 4 model uses the Morrow2 (UNOCAL) relative permeability curve with capillary pressure assigned heterogeneously by hydrostratigraphic unit

Capillary Pressure↗

Pc-to-RP Method

These excel spreadsheets were used to develop the two-phase and three-phase methods for converting capillary pressure to relative permeability. Data was collected by the SWP and literature sources cited in the spreadsheets.

Capillary Pressure↗