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 1,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Numerical study of heat transfer enhancement in laminar and turbulent flow in enhanced tubes with uniform and non-uniform thermal boundary conditions

Various types of tube enhancements have been proposed in order to increase heat transfer while minimizing any associated increases in pressure drop. This paper presents a numerical study of two different types of enhanced tubes: the relatively common twisted tape insert tube, and a novel tube with airfoil shaped pins on the inner surface. Each tube is studied in both the laminar and turbulent flow regime. In addition to the typical constant temperature and constant flux boundary conditions, a non-constant heat flux boundary condition is applied on one half of the tube, a situation which may be encountered when absorbing solar radiation. In general, each type of enhanced tube offers a thermal performance factor greater than unity at certain Re. In addition, the non-uniform heat flux boundary condition appears to function similarly as the uniform heat flux boundary condition.

Sandlin, M↗
NASA NTRS1989

Higher order dispersion in the propagation of a gravity wave packet

To the first order of approximation, the complex amplitude of a wave packet in an anisotropic and dispersive medium is convected with the group of velocity. However, a gravity wave is a vector wave. Its wave packet must be formed by superposition of various wave numbers with corresponding frequencies, as is the case for scalar waves, and additionally by superposing many eigenmodes which also depend on the wave number. To represent the vector wave packet self-consistently, it is found that a gradient term must be included in the expansion. For a Guassian wave packet, this gradient term is shown to have important implications on the velocity vector as represented by its hodograph. Numerical results show that the hodograph is influenced by the location of the relative position of interest from the center of a Gaussian pulse. Higher order expansion shows that an initial Gaussian wave packet will retain its Gaussian shape as it propagates, but the pulse will spread in all directions with its major axis undergoing a rotation. Numerical results indicate that these higher order dispersive effects may be marginally observable in the atmosphere.

Yeh, K. C.↗
NASA NTRS1993

Evaluation of a distributed catchment scale water balance model

The validity of some of the simplifying assumptions in a conceptual water balance model is investigated by comparing simulation results from the conceptual model with simulation results from a three-dimensional physically based numerical model and with field observations. We examine, in particular, assumptions and simplifications related to water table dynamics, vertical soil moisture and pressure head distributions, and subsurface flow contributions to stream discharge. The conceptual model relies on a topographic index to predict saturation excess runoff and on Philip's infiltration equation to predict infiltration excess runoff. The numerical model solves the three-dimensional Richards equation describing flow in variably saturated porous media, and handles seepage face boundaries, infiltration excess and saturation excess runoff production, and soil driven and atmosphere driven surface fluxes. The study catchments (a 7.2 sq km catchment and a 0.64 sq km subcatchment) are located in the North Appalachian ridge and valley region of eastern Pennsylvania. Hydrologic data collected during the MACHYDRO 90 field experiment are used to calibrate the models and to evaluate simulation results. It is found that water table dynamics as predicted by the conceptual model are close to the observations in a shallow water well and therefore, that a linear relationship between a topographic index and the local water table depth is found to be a reasonable assumption for catchment scale modeling. However, the hydraulic equilibrium assumption is not valid for the upper 100 cm layer of the unsaturated zone and a conceptual model that incorporates a root zone is suggested. Furthermore, theoretical subsurface flow characteristics from the conceptual model are found to be different from field observations, numerical simulation results, and theoretical baseflow recession characteristics based on Boussinesq's groundwater equation.

Troch, Peter A.↗
DOE OSTI2022

A Robust Numerical Treatment of Solid-Phase Diffusion in Pseudo Two-Dimensional Lithium-Ion Battery Models

Solid-phase diffusion in active materials of lithium-ion batteries significantly affects charging and safety-related behavior of lithium-ion batteries. Therefore, it is essential to develop an efficient and robust numerical algorithm for solving solid-phase diffusion equations in physics-based battery models. In this work, we discuss the origins of numerical instabilities that can occur when solving the solid-phase diffusion equations using iterative methods. Then, in order to resolve such issues, we propose a simple numerical treatment to the surface flux term of discretized solid-phase diffusion equations. To demonstrate its numerical robustness, the proposed method is implemented into a pseudo two-dimensional (P2D) physics-based battery model and simulations are conducted at wide ranges of operating conditions. Even with extremely poor initial guesses for the Li+ concentrations of the active materials, computations using the proposed method do not diverge and the their computational speeds are comparable to those with conventional initial guesses. Comprehensive tests of the proposed method are also performed with a dynamic current profile based on US06 driving profile and a multi-stage charging profile with very high initial C-rate (12C).

battery modeling↗
DOE OSTI2020

Simulation Analysis of Geometry and Material Effects for Dropkinson Bar

The reported research is motivated by the need to address a key issue affecting the Dropkinson bar apparatus. This unresolved issue is the interference of the stress wave reflected from the bar-beam boundary with the measurement of the stress-strain response of a material tested in the apparatus. The purpose of the wave beam that is currently connected to the bar is to dissipate the stress wave, but the portion of the wave reflected from the bar-beam boundary is still significant. First, we focused on understanding which parameters affect the reflected wave's arrival time at a strain gauge. Specifically, we used finite-element numerical simulations with the Sierra/SM module to study the effects of various bar-beam connection fixities, alternative wave beam materials, and alternative geometries of the Dropkinson bar system based on a monolithic design. The conclusion of this study is that a partial reflection always occurs at the bar-beam boundary (or, for a monolithic design, at a point where the bar geometry changes). Therefore, given a fixed total length of the bar, it is impossible to increase the reflected wave's arrival time by any significant amount. After reaching this conclusion, we focused instead on trying to minimize the energy of the reflected stress wave circulating up and down through the bar over a relatively long period of time (10 ms). Once again, we used numerical simulations with the Sierra/SM module to investigate the effects of various bar-beam connection fixities, alternative wave beam materials, and parameters of an asymmetric monolithic design of the bar-and-beam system. This study demonstrated that various parameters can significantly affect the energy of the wave reflections, with the difference between best and worst configurations being about one order of magnitude in terms of energy. Based on the obtained results, we conclude with concrete takeaways for Dropkinson bar users and propose potential directions for future research and optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Challenges and Opportunities in Magnetospheric Space Weather Prediction

Abstract Space Weather is the study of the dynamics of the coupled Solar‐Terrestrial environment, as these dynamics impact technological systems and human activity. This paper reviews a selection of the advances, challenges, and new opportunities for magnetospheric space weather, and while the focus is on specific phenomena, many other aspects of space weather will have similar challenges and needs. As a scientific field with direct applications, the field of space weather is partly driven by imperatives from both policy and operations. We provide an introduction to some of the context in which the field exists and discuss how this might shape future developments and norms within the space weather enterprise. We briefly examine benchmarking, as a policy and operationally driven activity, as it provides immediate societal relevance and an opportunity to stretch scientific understanding. As numerical space weather prediction now becomes routine, and exascale computing is in the near future, we identify challenges relating to computational expense and big data, capturing and accounting for uncertainties, and specification of boundary conditions. Here, as with the observations supporting numerical space weather prediction, the key challenge lies in extending the lead time of predictions. We also discuss the role of data, particularly in regard to model validation and empirical modeling. Due to the growing societal impact of space weather, we also examine the relationships between space weather and its terrestrial counterpart and look at the importance of continuous evaluation, monitoring progress in predictive capability, and communication with researchers, forecasters, and end users.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

Optimization and acceleration guidance of flight trajectories in a windshear

Guidance strategies for near-optimum performance in a wind shear are examined. The takeoff problem is considered with reference to flight in a vertical plane; the presence of a downdraft is assumed. Trajectories for optimum performance in a wind shear are determined for different wind shear models and intensities. Numerical experiments with the optimum control approach lead to the conclusion that, for weak to moderate shear/downdraft combinations, the optimal trajectory is characterized by a monotonic climb, and for severe shear/downdraft combinations, it is characterized by an initial climb, followed by a nearly horizontal flight, followed by renewed climbing after the aircraft has passed through the shear region. An acceleration guidance scheme based on relative acceleration is presented in both analytical form and feedback control form. Numerical results with this scheme result in trajectories close to the optimum and considerably superior to those arising from alternative guidance schemes.

Miele, A.↗
NASA NTRS1986

Spacecraft charging in the auroral plasma: Progress toward understanding the physical effects involved

The work is presented in four parts. First, main differences between the plasma environments in geostationary orbit and low polar orbit with regard to high-voltage charging situations are reviewed. Next, results are presented from a calculation of secondary-electron escape currents from negatively-charged spacecraft surfaces having various orientations relative to the local magnetic-field direction. It is shown that for finite ranges of combinations of electric and magnetic field directions, secondary-electron escape is completely suppressed and therefore cannot help to discharge the spacecraft. In such circumstances, secondary electrons may travel distances many times their gyroradii before reimpacting, and this may produce greatly increased secondary-electron surface currents. Thirdly, a simple rough estimate of the required conditions for high-voltage auroral-zone charging is developed. The results suggest that for any given spacecraft, surface potentials are likely to depend more strongly on the ratio of ambient flux of high-energy electrons to that of all ions, than on any other environmental parameter. Finally, preliminary results are presented of numerical simulation work directed toward testing this hypothesis. Numerical instabilities encountered in doing this simulation work probably are closely related to physical sensitivities inherent in the physics of the ion wake behind the spacecraft, and especially to beam-like constituents of the ion population in the wake.

Laframboise, J. G.↗
DOE OSTI2020

Intrinsic Elastic Anisotropy of Westerly Granite Observed by Ultrasound Measurements, Microstructural Investigations, and Neutron Diffraction

Westerly granite (WG) has been generally accepted as an isotropic homogeneous rock. Here, we return to WG and observe significant elastic anisotropy using multidirectional ultrasonic sounding on spherical samples at pressures up to 400 MPa. Thermal treatment of WG leads to formation of microcracks that reduce elastic wave velocities and increase its elastic anisotropy. The 3D distribution of P-wave velocities at low pressure is close to orthorhombic symmetry. Application of hydrostatic pressure closes most of thermally induced microcracks and decreases elastic anisotropy of WG, but at high pressure the anisotropy is practically reversed compared to low pressure: maximum P-wave velocity direction at low pressures is near minimum velocity direction at high pressure and vice versa. To understand this effect, microstructures of the rock were investigated by optical and scanning electron microscopy. Preferred orientations of four major rock-forming minerals—quartz, orthoclase, plagioclase, and biotite—were measured by time-of-flight neutron diffraction, which confirms significant crystal alignment. All these data were used to numerically model anisotropic elastic properties of WG. It is shown that WG possesses weak intrinsic elastic anisotropy related mainly to the preferred orientation of feldspars formed during igneous crystallization. Observed microcracks are mostly related to the cleavage planes of feldspars and biotite, and thus also demonstrate preferred orientation. Higher preheating temperatures produce larger quantity of longer microcracks. A numerical model shows that these microcracks act against the weak intrinsic elastic anisotropy of WG, and define the elastic anisotropy at low pressures.

58 GEOSCIENCES↗
DOE OSTI2021

Inclusive prompt photon-jet correlations as a probe of gluon saturation in electron-nucleus scattering at small x

We compute the differential cross-section for inclusive prompt photon+quark production in deeply inelastic scattering of electrons off nuclei at small x (e + A DIS) in the framework of the Color Glass Condensate effective field theory. The result is expressed as a convolution of the leading order (in the strong coupling α s ) impact factor for the process and universal dipole matrix elements, in the limit of hard photon transverse momentum relative to the nuclear saturation scale Q s,A (x). We perform a numerical study of this process for the kinematics of the Electron-Ion Collider (EIC), exploring in particular the azimuthal angle correlations between the final state photon and quark. We observe a systematic suppression and broadening pattern of the back-to-back peak in the relative azimuthal angle distribution, as the saturation scale is increased by replacing proton targets with gold nuclei. Our results suggest that photon+jet final states in inclusive e + A DIS at high energies are in general a promising channel for exploring gluon saturation that is complementary to inclusive and diffractive dijet production. They also provide a sensitive empirical test of the universality of dipole matrix elements when compared to identical measurements in proton-nucleus collisions. However because photon+jet correlations at small x in EIC kinematics require jet reconstruction at small k ⊥ , it will be important to study their feasibility relative to photon-hadron correlations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2001

Searching the ASRS Database Using QUORUM Keyword Search, Phrase Search, Phrase Generation, and Phrase Discovery

To support Search Requests and Quick Responses at the Aviation Safety Reporting System (ASRS), four new QUORUM methods have been developed: keyword search, phrase search, phrase generation, and phrase discovery. These methods build upon the core QUORUM methods of text analysis, modeling, and relevance-ranking. QUORUM keyword search retrieves ASRS incident narratives that contain one or more user-specified keywords in typical or selected contexts, and ranks the narratives on their relevance to the keywords in context. QUORUM phrase search retrieves narratives that contain one or more user-specified phrases, and ranks the narratives on their relevance to the phrases. QUORUM phrase generation produces a list of phrases from the ASRS database that contain a user-specified word or phrase. QUORUM phrase discovery finds phrases that are related to topics of interest. Phrase generation and phrase discovery are particularly useful for finding query phrases for input to QUORUM phrase search. The presentation of the new QUORUM methods includes: a brief review of the underlying core QUORUM methods; an overview of the new methods; numerous, concrete examples of ASRS database searches using the new methods; discussion of related methods; and, in the appendices, detailed descriptions of the new methods.

McGreevy, Michael W.↗
NASA NTRS1987

Numerical arc segmentation algorithm for a radio conference-NASARC (version 2.0) technical manual

The information contained in the NASARC (Version 2.0) Technical Manual (NASA TM-100160) and NASARC (Version 2.0) User's Manual (NASA TM-100161) relates to the state of NASARC software development through October 16, 1987. The Technical Manual describes the Numerical Arc Segmentation Algorithm for a Radio Conference (NASARC) concept and the algorithms used to implement the concept. The User's Manual provides information on computer system considerations, installation instructions, description of input files, and program operating instructions. Significant revisions have been incorporated in the Version 2.0 software. 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 within 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 effecting an overall reduction in computer run time.

Whyte, Wayne A., Jr.↗
NASA NTRS1992

A numerical field experiment approach for determining probabilities of microburst intensity

Several investigators had determined that some atmospheric parameters were related to the formation and severity of microbursts. For example, Caracena pointed out the relationship between a dry adiabatic lapse rate and microbursts in 'The crash of Delta Flight 191 at Dallas-Fort Worth international airport'. These early investigations led to the idea that numeric modeling of microbursts with varying atmospheric parameters might define 'signatures' that could lead to determining the probability of microburst intensity. The idea was that, by using already available sensors (such as static air temperature, pressure altitude, and radar reflectivity) onboard an aircraft, a reliable prediction of microburst existence and intensity could be formed. Such data could be used to create an 'expert meteorologist' using either artificial intelligence or other techniques that could be used in either reactive or look-ahead systems to vary sensitivity thresholds and coordinate the inputs from different detecting systems. To this end, Honeywell contracted to have the microburst simulations run. The questions to be addressed were the following: using the sensor set available to the aircraft (e.g. temperature, radar reflectivity, etc.), can we calculate the probability that (1) a microburst could be formed? and (2) that the resultant winds would be of sufficient magnitude to threaten the aircraft? Over a two year period, a data set of 1800 microburst simulations was accumulated. Verification of the microburst simulation was obtained using the results of other independent researchers and actual comparison to microburst events in Orlando and Denver. Some of the results from the simulation have already been incorporated into Honeywell's Windshear Detection and Guidance System with excellent results. Various aspects of this investigation are presented in viewgraph form.

Droegemeier, Kelvin K.↗
NASA NTRS2010

The Writhe of Helical Structures in the Solar Corona

Context. Helicity is a fundamental property of magnetic fields, conserved in ideal MHD. In flux rope topology, it consists of twist and writhe helicity. Despite the common occurrence of helical structures in the solar atmosphere, little is known about how their shape relates to the writhe, which fraction of helicity is contained in writhe, and how much helicity is exchanged between twist and writhe when they erupt. Aims. Here we perform a quantitative investigation of these questions relevant for coronal flux ropes. Methods. The decomposition of the writhe of a curve into local and nonlocal components greatly facilitates its computation. We use it to study the relation between writhe and projected S shape of helical curves and to measure writhe and twist in numerical simulations of flux rope instabilities. The results are discussed with regard to filament eruptions and coronal mass ejections (CMEs). Results. (1) We demonstrate that the relation between writhe and projected S shape is not unique in principle, but that the ambiguity does not affect low-lying structures, thus supporting the established empirical rule which associates stable forward (reverse) S shaped structures low in the corona with positive (negative) helicity. (2) Kink-unstable erupting flux ropes are found to transform a far smaller fraction of their twist helicity into writhe helicity than often assumed. (3) Confined flux rope eruptions tend to show stronger writhe at low heights than ejective eruptions (CMEs). This argues against suggestions that the writhing facilitates the rise of the rope through the overlying field. (4) Erupting filaments which are S shaped already before the eruption and keep the sign of their axis writhe (which is expected if field of one chirality dominates the source volume of the eruption), must reverse their S shape in the course of the rise. Implications for the occurrence of the helical kink instability in such events are discussed.

WRITHE↗
DOE OSTI2021

Multi-fidelity modeling to predict the rheological properties of fiber suspensions

Unveiling the rheological properties of fiber suspensions is of paramount interest to many industrial applications like biofuel production. The 3D numerical simulations of the suspension of fibers are often computationally expensive and time-consuming. Machine learning methods such as neural networks can simplify the prediction of rheological behavior; however, they require a relatively large training data set. Multi -fidelity models, which combine high-fidelity data from numerical simulations and less expensive lower fidelity data from resources such as simplified physical equations, can lead to optimized predictions. Here, we focus on a neural network with two levels of fidelity, i.e., high and low fidelity networks. To produce high-fidelity data, we perform direct numerical simulations to model the fibers as one-dimensional inextensible slender bodies that obey the Euler- Bernoulli beam equation. The Navier-Stokes equations govern the suspended fluid, and an immersed boundary method is used to couple the fluid and solid motion. The low-fidelity data is produced by using constitutive equations. Noticeable improvements have been observed in the accuracy of predicting the rheological behavior when a multi-fidelity network is used compared to the single-fidelity network.

Boodaghidizaji, Miad↗
NASA NTRS2005

The QBO as Potential Amplifier and Conduit to Lower Altitudes of Solar Cycle Influence

The solar cycle (SC) effect in the lower atmosphere has been linked observationally to the Quasi-biennial Oscillation (QBO), which is generated primarily by small-scale gravity waves. Salby and Callaghan analyzed the QBO observations covering more than 40 years and found that it contains a relatively large SC signature at 20 km. Following up on a 2D study with our Numerical Spectral Model (NSM), we discuss here a 3D study in which we simulated the QBO under the influence of the SC. For a SC period of 10 years, the amplitude of the relative variations of radiative forcing is taken to vary from 0.2% at the surface to 2% at 50 km to 20% at 100 km and above. This model produces in the lower stratosphere a relatively large modulation of the QBO, which appears to be related to the SC and is in qualitative agreement with the observations. Further studies are needed, (1) to determine whether the effect is real and the results are robust and (2) to explore the mechanism(s) that may amplify the SC effect. Quasi-decadal oscillations, generated internally by the QBO interacting with the seasonal cycles, may interfere with or aid the SC effect.

Mayr, Hans G.↗
NASA NTRS1985

Implicit dual control for general stochastic systems

A new implicit dual technique is presented for stochastic adaptive control synthesis. In this technique, denoted as the method of utility costs (MOUC), a new control policy is derived from a known starting control policy by using a single iteration in time-varying policy space. An important result is that the new control policy yields a corresponding improvement in adaptive performance relative to the known starting policy. Since the starting policy can be chosen arbitrarily, the new policy has potentially superior performance relative to any existing design. Unlike other implicit dual designs, the MOUC can be applied to systems having constrained inputs, discontinuities, discrete valued disturbances, as well as in situations where the variances of the associated process statistics are large. A numerical example is given in which a controller is synthesized to provide improved performance relative to the Open-Loop Feedback policy of Dreyfus.

Bayard, D. S.↗