Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical modelling”

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,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Verification of thermal analysis codes for modeling solid rocket nozzles

One of the objectives of the Solid Propulsion Integrity Program (SPIP) at Marshall Space Flight Center (MSFC) is development of thermal analysis codes capable of accurately predicting the temperature field, pore pressure field and the surface recession experienced by decomposing polymers which are used as thermal barriers in solid rocket nozzles. The objective of this study is to provide means for verifications of thermal analysis codes developed for modeling of flow and heat transfer in solid rocket nozzles. In order to meet the stated objective, a test facility was designed and constructed for measurement of the transient temperature field in a sample composite subjected to a constant heat flux boundary condition. The heating was provided via a steel thin-foil with a thickness of 0.025 mm. The designed electrical circuit can provide a heating rate of 1800 W. The heater was sandwiched between two identical samples, and thus ensure equal power distribution between them. The samples were fitted with Type K thermocouples, and the exact location of the thermocouples were determined via X-rays. The experiments were modeled via a one-dimensional code (UT1D) as a conduction and phase change heat transfer process. Since the pyrolysis gas flow was in the direction normal to the heat flow, the numerical model could not account for the convection cooling effect of the pyrolysis gas flow. Therefore, the predicted values in the decomposition zone are considered to be an upper estimate of the temperature. From the analysis of the experimental and the numerical results the following are concluded: (1) The virgin and char specific heat data for FM 5055 as reported by SoRI can not be used to obtain any reasonable agreement between the measured temperatures and the predictions. However, use of virgin and char specific heat data given in Acurex report produced good agreement for most of the measured temperatures. (2) Constant heat flux heating process can produce a much higher heating rate than the radiative heating process. The results show that heating rates of 125 C/s (225 F/s) can be achieved. (3) The electrical resistance of the FM 5055 samples were about 150 Omega in the virgin state, and decreased to about 35 Omega when charred. A reliable scheme must be developed to electrically insulate the composite from the heater foil in order to prevent any current leakage into the sample which can result in volumetric heating of the char zone of the composite.

Keyhani, M.↗
NASA NTRS1993

Convection equation modeling: A non-iterative direct matrix solution algorithm for use with SINDA

The determination of the boundary conditions for a component-level analysis, applying discrete finite element and finite difference modeling techniques often requires an analysis of complex coupled phenomenon that cannot be described algebraically. For example, an analysis of the temperature field of a coldplate surface with an integral fluid loop requires a solution to the parabolic heat equation and also requires the boundary conditions that describe the local fluid temperature. However, the local fluid temperature is described by a convection equation that can only be solved with the knowledge of the locally-coupled coldplate temperatures. Generally speaking, it is not computationally efficient, and sometimes, not even possible to perform a direct, coupled phenomenon analysis of the component-level and boundary condition models within a single analysis code. An alternative is to perform a disjoint analysis, but transmit the necessary information between models during the simulation to provide an indirect coupling. For this approach to be effective, the component-level model retains full detail while the boundary condition model is simplified to provide a fast, first-order prediction of the phenomenon in question. Specifically for the present study, the coldplate structure is analyzed with a discrete, numerical model (SINDA) while the fluid loop convection equation is analyzed with a discrete, analytical model (direct matrix solution). This indirect coupling allows a satisfactory prediction of the boundary condition, while not subjugating the overall computational efficiency of the component-level analysis. In the present study a discussion of the complete analysis of the derivation and direct matrix solution algorithm of the convection equation is presented. Discretization is analyzed and discussed to extend of solution accuracy, stability and computation speed. Case studies considering a pulsed and harmonic inlet disturbance to the fluid loop are analyzed to assist in the discussion of numerical dissipation and accuracy. In addition, the issues of code melding or integration with standard class solvers such as SINDA are discussed to advise the user of the potential problems to be encountered.

Schrage, Dean S.↗
DOE OSTI2020

Reversible structural transformations in supercooled liquid water from 135 to 245 K

Water has many anomalous properties compared to “simple” liquids, and these anomalies are typically enhanced in supercooled water.1-3 While numerous models have been proposed, including the liquid-liquid critical point,4,5 the singularity-free scenario,6 and the stability limit conjecture,1 a molecular-level understanding remains elusive. The main difficulty in determining which, if any, of these models is correct is the limited amount of data in the relevant temperature and pressure ranges. For water at ambient pressures, which is the focus of this work, data is largely missing from 160 – 232 K (“No Man’s Land”) due to rapid crystallization.2,3 Whether rapid crystallization is just an experimental obstacle, or a fundamental problem signaling the inability of water to thermally equilibrate prior to crystallization is also a major unanswered question.5,7 Here, we investigate the structural transformations of transiently-heated, supercooled water with nanosecond time resolution using infrared vibrational spectroscopy. The experiments demonstrate three key results. First, water’s structure relaxes from its initial configuration to a “steady-state” configuration prior to the onset of crystallization over a wide temperature range. Second, water’s steady-state structure can be reproduced by a linear combination of two, temperature-independent structures that correspond to a “high-temperature liquid” and a “low-temperature liquid.” Third, the observed structural changes are reversible over the full temperature range. Taken together, these results show that supercooled water can equilibrate prior to crystallization for temperatures from the homogeneous nucleation temperature, TH ~232 K,3 down to the glass transition temperature (Tg ~ 136 K). Second, the results provide support for the hypothesis that supercooled water can be described as a mixture of two, structurally-distinct, interconvertible liquids from 135 K to 245 K.5,8-18

Kringle, Loni M.↗
DOE OSTI2022

Volumetric Rendering on Wavelet-Based Adaptive Grid

Numerical modeling of physical phenomena frequently involves processes across a wide range of spatial and temporal scales. In the last two decades, the advancements in wavelet-based numerical methodologies to solve partial differential equations, combined with the unique properties of wavelet analysis to resolve localized structures of the solution on dynamically adaptive computational meshes, make it feasible to perform large-scale numerical simulations of a variety of physical systems on a dynamically adaptive computational mesh that changes both in space and time. Volumetric visualization of the solution is an essential part of scientific computing, yet the existing volumetric visualization techniques do not take full advantage of multi-resolution wavelet analysis and are not fully tailored for visualization of a compressed solution on the wavelet-based adaptive computational mesh. Our objective is to explore the alternatives for the visualization of time-dependent data on space-time varying adaptive mesh using volume rendering while capitalizing on the available sparse data representation. Two alternative formulations are explored. The first one is based on volumetric ray casting of multi-scale datasets in wavelet space. Rather than working with the wavelets at the finest possible resolution, a partial inverse wavelet transform is performed as a preprocessing step to obtain scaling functions on a uniform grid at a user-prescribed resolution. As a result, a solution in physical space is represented by a superposition of scaling functions on a coarse regular grid and wavelets on an adaptive mesh. An efficient and accurate ray casting algorithm is based just on these coarse scaling functions. Additional details are added during the ray tracing by taking an appropriate number of wavelets into account based on support overlap with the interpolation point, wavelet coefficient magnitude, and other characteristics, such as opacity accumulation (front to back ordering) and deviation from frontal viewing direction. The second approach is based on complementing of wavelet-based adaptive mesh to the traditional Adaptive Mesh Refinement (AMR) mesh. Both algorithms are illustrated and compared to the existing volume visualization software for Rayleigh-Benard thermal convection and electron density data sets in terms of rendering time and visual quality for different data compression of both wavelet-based and AMR adaptive meshes.

Vezolainen, Alexei V.↗
NASA NTRS1992

Mechanical energy flow models of rods and beams

It has been proposed that the flow of mechanical energy through a structural/acoustic system may be modeled in a manner similar to that of flow of thermal energy/in a heat conduction problem. If this hypothesis is true, it would result in relatively efficient numerical models of structure-borne energy in large built-up structures. Fewer parameters are required to approximate the energy solution than are required to model the characteristic wave behavior of structural vibration by using traditional displacement formulations. The energy flow hypothesis is tested in this investigation for both longitudinal vibration in rods and transverse flexural vibrations of beams. The rod is shown to behave approximately according to the thermal energy flow analogy. However, the beam solutions behave significantly differently than predicted by the thermal analogy unless locally-space-averaged energy and power are considered. Several techniques for coupling dissimilar rods and beams are also discussed. Illustrations of the solution accuracy of the methods are included.

Wohlever, J. C.↗
NASA NTRS1994

The r-process and neutrino-heated supernova ejecta

As a neutron star is formed by the collapse of the iron core of a massive star, its Kelvin-Helmholtz evolution is characterized by the release of gravitational binding energy as neutrinos. The interaction of these neutrinos with heated material above the neutron star generates a hot bubble in an atmosphere that is nearly in hydrostatic equilibrium and heated, after approximately 10 s, to an entropy of S/N(sub AS)k greater than or approximately = 400. The neutron-to-proton ratio for material moving outward through this bubble is set by the balance between neutrino and antineutrino capture on nucleons. Because the electron antineutrino spectrum at this time is hotter than the electron neutrino spectrum, the bubble is neutron-rich (0.38 less than or approximately = Y(sub e) less than or approximately = 0.47). Previous work using a schematic model has shown that these conditions are well suited to the production of heavy elements by the r-process. In this paper we have advanced the numerical modeling of a 20 solar mass 'delayed' supernova explosion to the point that we can follow the detailed evolution of material moving through the bubble at the late times appropiate to r-process nucleosynthesis. The supernova model predicts a final kinetic energy for the ejecta of 1.5 x 10(exp 51) ergs and leaves behind a remnant with a baryon mass of 1.50 solar mass (and a gravitational mass of 1.445 solar mass). We follow the thermodynamic and compositional evolution of 40 trajectories in rho(t), T(t), Y(sub e)(t) for a logarithmic grid of mass elements for the last approximately = 0.03 solar mass to be ejected by the proto-neutron star down to the last less than 10(exp -6) solar mass of material expelled at up to approximately = 18 s after core collapse. We find that an excellent fit to the solar r-process abundance distribution is obtained with no adjustable parameters in the nucleosynthesis calculations. Moreover, the abundances are produced in the quantities required to account for the present Galactic abundances. However, at earlier times, this one-dimensional model ejects too much material with entropies S/N(sub A)k approximately 50 and Y(sub e) approximately 0.46. This leads to an acceptable over production of N = 50 nuclei, particularly Sr-88, Y-89, and Zr-90, relative to their solar abundances. We speculate on various means to avoid the early overproduction and/or ejection of N = 50 isotonic nuclei while still producing and ejecting the correct amount of r-process material.

Woosley, S. E.↗
NASA NTRS2014

Concerning the Charging of an Exploration Craft on and near a Small Asteroid

Introduction: An object immersed in an airless plasma environment will experience a natural process of surface charging in order to acheieve current balance, or zero net electric current to the object. It has been shown in recent computer simulations that the small-body plasma environment is very complex [1], considering effects of photoemission, topography, and formation of a plasma wake. For this work we consider an exploration craft (or astronaut) immersed within a plasma environment near an asteroid, which exhibits widely varying solar wind and photoelectric particle fluxes and continuously evolving illumination conditions. Objective: We aim to determine how an explo-ration craft or astronaut suit accumulates charge while located in the "nightside" asteroid wake where the particle fluxes are reduced, and in the dayside near-surface photoelectron sheath, by combining an object charging model [2] with kinetic simulations of a near-asteroid plasma environment [1]. We consider an astronaut floating near the asteroid while not in contact with the surface, as well as an astronaut moving along the surface using their hands/gloves to crawl along. Results: The modeling results suggest that remediation of triboelectric charge via accumulation of plasma currents is an important factor to consider when designing future NEA mission infrastructure, especially if repeated and frequent contact with the surface is planned. In shadowed regions such as the location shown in Fig. 1a, the plasma currents are so low (and the effective charge-remediation timescale so long, e.g. minutes to hours) that repeated contact with the surface tribocharges the glove in an uncontrollable fashion, as shown for two representative electron temperatures in Fig. 2a. The resulting buildup of significant negative charge would eventually initiate some other "current of last resort" [4] such as transport of positively-charged dust, field-emission from the glove, or significant alteration of environmental ion currents within the wake. In contrast, the few-meters-thick dayside photoelectron sheath in which the astronaut of Fig. 1b is immersed in is so rich in electrons (and hence so electrically conductive) that accumulated tribocharge dissipates almost instantaneously (e.g. in less than a ms) as shown in Fig. 2b. As our model astronaut orbits the NEA they would experience plasma currents and associated charge re-mediation times spanning many orders of magnitude, and the fusion between our numerical models provides a detailed understanding of the charging hazards possibly associated with contact-based NEA exploration.

Asteroid↗
NASA NTRS2002

Meteorology, Macrophysics, Microphysics, Microwaves, and Mesoscale Modeling of Mediterranean Mountain Storms: The M8 Laboratory

Comprehensive understanding of the microphysical nature of Mediterranean storms can be accomplished by a combination of in situ meteorological data analysis and radar-passive microwave data analysis, effectively integrated with numerical modeling studies at various scales, from synoptic scale down through the mesoscale, the cloud macrophysical scale, and ultimately the cloud microphysical scale. The microphysical properties of and their controls on severe storms are intrinsically related to meteorological processes under which storms have evolved, processes which eventually select and control the dominant microphysical properties themselves. This involves intense convective development, stratiform decay, orographic lifting, and sloped frontal lifting processes, as well as the associated vertical motions and thermodynamical instabilities governing physical processes that affect details of the size distributions and fall rates of the various types of hydrometeors found within the storm environment. Insofar as hazardous Mediterranean storms, highlighted in this study by three mountain storms producing damaging floods in northern Italy between 1992 and 2000, developing a comprehensive microphysical interpretation requires an understanding of the multiple phases of storm evolution and the heterogeneous nature of precipitation fields within a storm domain. This involves convective development, stratiform transition and decay, orographic lifting, and sloped frontal lifting processes. This also involves vertical motions and thermodynamical instabilities governing physical processes that determine details of the liquid/ice water contents, size disi:ributions, and fall rates of the various modes of hydrometeors found within hazardous storm environments.

Starr, David O.↗
NASA NTRS1993

Global dynamic modeling of a transmission system

The work performed on global dynamic simulation and noise correlation of gear transmission systems at the University of Akron is outlined. The objective is to develop a comprehensive procedure to simulate the dynamics of the gear transmission system coupled with the effects of gear box vibrations. The developed numerical model is benchmarked with results from experimental tests at NASA Lewis Research Center. The modal synthesis approach is used to develop the global transient vibration analysis procedure used in the model. Modal dynamic characteristics of the rotor-gear-bearing system are calculated by the matrix transfer method while those of the gear box are evaluated by the finite element method (NASTRAN). A three-dimensional, axial-lateral coupled bearing model is used to couple the rotor vibrations with the gear box motion. The vibrations between the individual rotor systems are coupled through the nonlinear gear mesh interactions. The global equations of motion are solved in modal coordinates and the transient vibration of the system is evaluated by a variable time-stepping integration scheme. The relationship between housing vibration and resulting noise of the gear transmission system is generated by linear transfer functions using experimental data. A nonlinear relationship of the noise components to the fundamental mesh frequency is developed using the hypercoherence function. The numerically simulated vibrations and predicted noise of the gear transmission system are compared with the experimental results from the gear noise test rig at NASA Lewis Research Center. Results of the comparison indicate that the global dynamic model developed can accurately simulate the dynamics of a gear transmission system.

Choy, F. K.↗
DOE OSTI2023

Automated power-following control for nuclear thermal propulsion startup and shutdown using MOOSE-based applications

This paper describes an investigation of various automated control strategies applied to a full-core multiphysics Griffin/Bison/RELAP-7 model of a prototypical nuclear thermal propulsion system. In all cases, control is achieved by actuating control drums based on the demanded power and predicted quantities from the numerical model. One key finding is that hybrid proportional integral derivative controllers – a novel type of controller that uses both power and reactivity predicted signals – can demonstrate a level of performance rivaling that of period-generated control. The former requires parameter tuning, while the latter mostly necessitates providing reactivity coefficients and temperature rates of change, which could be very challenging to accurately measure in real-time. In addition, decay heat plays an important role in determining cooling requirements during cooldown phases. A decay heat model, accounting for burn time and throttling, was derived, and is incorporated within the model to simulate the steady-state and shutdown phases and satisfyingly follow the power demand. However, temperature overshoots and non-constant specific impulse during throttling will necessitate further improvements.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2019

Secondary Pore Network Design Optimization Analytical Diffusion Model for Lithium Ion Battery

Fast ion transport along the electrode thickness is required for Lithium ion battery (LIB) to sustain continuous high-rate discharge or charge currents, which is a desired feature for the fast charging of electric vehicle batteries. Advance in LIB electrode manufacturing allows to control, to a certain extent, the electrode microstructure architecture in order to improve the macroscopic performances of these devices. For instance, several groups have experimentally tested architectures tailored specifically to increase the through plane ionic diffusion in regards for fast charge application by adding channels or groove-like large pores oriented along the electrode thickness. These voids do not replace the existing pores of the electrode microstructure but consist in a second, coarser, pore network, called the secondary pore network (SPN). The model quantifies the impact of SPN on ionic diffusion with a composite electrode through a simple analytical approach, considering various pore channel geometries and comparing to standard electrodes with identical gravimetric and identical volumetric specific theoretical capacity. The goal is to find optimal SPN design parameters according to three optimization objectives aiming to balance the improved overall through-plane, thanks to the coarse aligned channels, and the degraded in-plane diffusion, due to the porous matrix densification required to maintain gravimetric and volumetric specific theoretical capacity. The analytical model must be considered as a pre-screening tool to quickly investigate large design space, that can be refined later with more complex electrochemical CPU-expensive numerical models.

Usseglio Viretta, Francois↗
NASA NTRS1984

On the Kelvin-Helmholtz instabilities of supersonic shear layers

Analytical consideration is given to a Kelvin-Helmholtz instability in a finite width supersonic shear layer with linear and sinusoidal velocity profiles in the transition region. An analogy is explored between the pressure perturbation and the bound and scattering wave states of the one-dimensional Schroedinger equation. A numerical model is obtained for generating contour maps of the real and imaginary parts of the perturbation frequency in the (A,B)-plane. A comparison is made between the numerical results and analytical expressions derived for very small Mach numbers. Two unstable regimes are characterized in the (A,B)-plane. The method is concluded useful for both laboratory and astrophysics problems.

Choudhury, S. R.↗
NASA NTRS1985

Atmospheric electric field and current configurations in the vicinity of mountains

A number of investigations have been conducted regarding the electrical distortion produced by the earth's orography. Hays and Roble (1979) utilized their global model of atmospheric electricity to study the effect of large-scale orographic features on the currents and fields of the global circuit. The present paper is concerned with an extension of the previous work, taking into account an application of model calculations to orographic features with different configurations and an examination of the electric mapping of these features to ionospheric heights. A two-dimensional quasi-static numerical model of atmospheric electricity is employed. The model contains a detailed electrical conductivity profile. The model region extends from the surface to 100 km and includes the equalization layer located above approximately 70 km. The obtained results show that the electric field and current configurations above mountains depend upon the curvature of the mountain slopes, on the width of the mountain, and on the columnar resistance above the mountain (or mountain height).

Tzur, I.↗
NASA NTRS1990

Low-temperature Galactic fountains

Analytic and one- and two-dimensional hydrodynamic calculations are performed to determine whether or not the Galactic fountain model can explain observations of neutral gas in the Galactic halo. Galactic halo gasdynamics are examined, including thermal conduction, Galactic rotation, and Galactic fountains. Results are presented from numerical modeling of low-temperature Galactic fountains. It is found that the best reproduction of the observations is provided by a model near the transonic regime with a temperature and density at the base of the fountain of 300,000 K and 0.001/cu cm, respectively. The cooling time for the hot gas is 33 million yr and the time for newly formed clouds to return to the disk is 47 million yr, neglecting drag forces.

Houck, John C.↗
NASA NTRS1991

Modeled performance of monolithic, 3-terminal InP/Ga(0.47)In(0.53)As concentrator solar cells as a function of temperature and concentration ratio

Using measured device parameters from a monolithic three-terminal In/GaInAs (0.75 eV) tandem concentrator solar cell, a numerical model is constructed that calculates efficiency as a function of temperature and concentration ratio. The device measurements indicate that the series resistance in the InP top cell severely limits the maximum efficiency at high concentration ratios. Results from the model in which a single-junction InP concentrator solar cell that has a lower resistance by a factor of six was substituted for the top cell show that peak 20 C air-mass-zero (AM0) efficiencies should approach 30 percent at concentration ratios greater than 100. At 80 C, this tandem cell should exceed 24 percent AM0 efficiency.

Osterwald, C. R.↗
NASA NTRS1993

Advanced adaptive computational methods for Navier-Stokes simulations in rotorcraft aerodynamics

A phase 2 research and development effort was conducted in area transonic, compressible, inviscid flows with an ultimate goal of numerically modeling complex flows inherent in advanced helicopter blade designs. The algorithms and methodologies therefore are classified as adaptive methods, which are error estimation techniques for approximating the local numerical error, and automatically refine or unrefine the mesh so as to deliver a given level of accuracy. The result is a scheme which attempts to produce the best possible results with the least number of grid points, degrees of freedom, and operations. These types of schemes automatically locate and resolve shocks, shear layers, and other flow details to an accuracy level specified by the user of the code. The phase 1 work involved a feasibility study of h-adaptive methods for steady viscous flows, with emphasis on accurate simulation of vortex initiation, migration, and interaction. Phase 2 effort focused on extending these algorithms and methodologies to a three-dimensional topology.

Stowers, S. T.↗
NASA NTRS1996

A New Axisymmetric MHD Model of the Interaction of the Solar Wind with Venus

A new two-dimensional axisymmetric MHD model is used to study the interaction of the solar wind with Venus under conditions where the interplanetary field is approximately aligned with the solar wind velocity. This numerical model solves the MHD transport equations for density, velocity, pressure, and magnetic field on an adaptively refined, unstructured grid system. This use of an adaptive grid allows high spatial resolution in regions of large density/velocity gradients and yet can be run on a workstation. The actual grid sizes vary from about 0.06 R(sub v) near the bowshock to 2 R(sub v) in the unperturbed solar wind. The results of the calculations are compared with observed magnetic field values obtained from the magnetometer on the Pioneer Venus Orbiter, at a time when the angle between the solar wind velocity vector and the interplanetary magnetic field (IMF) was only 7.6 deg. Good qualitative agreement between the observed and calculated field behavior is found. The overall results suggest that the induced magnetotail disappears when the IMF is radial for an extended time period and implies that it weakens when the field rotated through a near-radial orientation.

DeZeeuw, Darren L.↗
NASA NTRSDate not supplied

Liquid Hydrogen Tank Chill and No-Vent Fill Prediction using Computational Fluid Dynamics

Cryogenic tank chill and fill is an important cryogenic fluid management (CFM) technology that supports and enables many of NASA’s long-duration space missions. For no-vent fill, the receiver tank pressure remains below the supply tank pressure during the entire duration of the fill so that the tank does not require venting. This is especially advantageous for tank fill operations in low gravity where the position of the liquid is not always known and venting the tank may cause loss of propellant by venting liquid. In lieu of expensive tests conducted on-orbit, accurate computational models capable of predicting receiver tank pressure during cryogenic propellant tank fill may be used to reduce system and propellant mass as well as mission risk. However, these numerical models must be validated or anchored to test data. This study presents a computational fluid dynamics (CFD) model with conjugate heat transfer that is used to predict a liquid hydrogen tank chill and no-vent fill ground test conducted at Lewis Research Center (now Glenn Research Center) in 1991. The specific test case chosen for model validation implemented an upward-facing jet near the bottom of a cylindrical 34 liter tank. CFD predictions are compared to experimental measurements of tank pressure, fill level, fluid temperatures, and wall temperatures. The CFD results show reasonable agreement to the test data but overpredict the pressure collapse near the end of the fill despite the liquid jet penetrating the liquid-vapor interface. Several sensitivity studies are considered due to notable uncertainties in the experiment.

Computational Fluid Dynamics↗