Search NASASearch

SEARCH · Search NASA

Results for “numerical models”

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 163 records · Page 9

A study of circumferentially-heated and block-heated heat pipes. I - Experimental analysis and generalized analytical prediction of capillary limits. II - Three-dimensional numerical modeling as a conjugate problem

The wall and centerline vapor temperatures and heat output are determined experimentally for a low-temperature copper-water heat pipe under uniform circumferential heating and block heating. The time required to reach a vapor temperature of 60 C from an initial ambient temperature of 21 C is determined for both modes of heating. The experimental capillary limit of the heat pipe is compared to the generalized capillary limits for block-heated pipes over a range of vapor temperatures. A three-dimensional numerical model is then developed for determining the temperature, pressure, and velocity distributions in the entire domain of a circumferentially heated and a block-heated pipe. The problem is formulated as a conjugate problem, without the assumption of a uniform vapor temperature. The predictions of the model are found to be in good agreement with the experimental data.

Schmalhofer, Joseph

Numerical Modeling of Pulse Detonation Rocket Engine Gasdynamics and Performance

Pulse detonation engines (PDB) have generated considerable research interest in recent years as a chemical propulsion system potentially offering improved performance and reduced complexity compared to conventional gas turbines and rocket engines. The detonative mode of combustion employed by these devices offers a theoretical thermodynamic advantage over the constant-pressure deflagrative combustion mode used in conventional engines. However, the unsteady blowdown process intrinsic to all pulse detonation devices has made realistic estimates of the actual propulsive performance of PDES problematic. The recent review article by Kailasanath highlights some of the progress that has been made in comparing the available experimental measurements with analytical and numerical models.

Morris, C. I.

Numerical Modeling & Optimization of the iProTech Pitching Inertial Pump (PIP) Wave Energy Converter (WEC) (CRADA Final Report)

This project represents a continuation of the collaboration between iProTech and NLR to simulate, optimize and design the iProTech Pitching Inertial Pump (PIP) device. The objectives of this TEAMER project are twofold: 1. Refining the physical characteristics of the existing iProTech PIP WEC-Sim model to enhance the model’s fidelity and include controllable components. Key model enhancements target the inclusion of Coulomb friction, the introduction of a controllable bypass valve, and the replacement of traditional check valves with advanced motorized ones. 2. Exploring traditional and advanced control algorithms. From traditional methods like latching control to cutting-edge reinforcement learning (RL) algorithms, the goal is to ensure the PIP device's adaptability and optimal performance across a range of ocean conditions. NLR is tasked with augmenting the WEC-Sim model and implementing the control algorithms, culminating in performance comparison analyses. iProTech will update their existing 3D models, advise on model improvements, and determine crucial system metrics. WEC-Sim, developed in MATLAB/SIMULINK with Simscape Multibody, is the main piece of software that will be used in this project. Coupled with the MATLAB RL Toolbox, it offers a robust platform for in-depth simulation and optimization of the iProTech PIP device. Building on previous work to explore the PIP design space and optimize its geometry, mass distribution, center of gravity and other key parameters, this project aims to refine iProTech’s existing numerical models and develop effective control algorithms that can seamlessly integrate into their future hardware testing campaigns.

16 TIDAL AND WAVE POWER

Axisymmetric Numerical Modeling of Pulse Detonation Rocket Engines

Pulse detonation rocket engines (PDREs) have generated research interest in recent years as a chemical propulsion system potentially offering improved performance and reduced complexity compared to conventional rocket engines. The detonative mode of combustion employed by these devices offers a thermodynamic advantage over the constant-pressure deflagrative combustion mode used in conventional rocket engines and gas turbines. However, while this theoretical advantage has spurred considerable interest in building PDRE devices, the unsteady blowdown process intrinsic to the PDRE has made realistic estimates of the actual propulsive performance problematic. The recent review article by Kailasanath highlights some of the progress that has been made in comparing the available experimental measurements with analytical and numerical models. In recent work by the author, a quasi-one-dimensional, finite rate chemistry CFD model was utilized to study the gasdynamics and performance characteristics of PDREs over a range of blowdown pressure ratios from 1-1000. Models of this type are computationally inexpensive, and enable first-order parametric studies of the effect of several nozzle and extension geometries on PDRE performance over a wide range of conditions. However, the quasi-one-dimensional approach is limited in that it cannot properly capture the multidimensional blast wave and flow expansion downstream of the PDRE, nor can it resolve nozzle flow separation if present. Moreover, the previous work was limited to single-pulse calculations. In this paper, an axisymmetric finite rate chemistry model is described and utilized to study these issues in greater detail. Example Mach number contour plots showing the multidimensional blast wave and nozzle exhaust plume are shown. The performance results are compared with the quasi-one-dimensional results from the previous paper. Both Euler and Navier-Stokes solutions are calculated in order to determine the effect of viscous effects in the nozzle flowfield. Additionally, comparisons of the model results to performance data from CalTech, as well as experimental flowfield measurements from Stanford University, are also reported.

Morris, Christopher I.

First-Principle Dynamic Electro-Thermal Numerical Model of a Scanning Radiometer for Earth Radiation Budget Applications

Low Earth Observing instruments that are used to monitor the incoming solar and outgoing long wave radiation have been a crucial part of studying the Earth’s radiation budget for the past three decades. These instruments go through several robust design phases followed by vigorous ground calibration campaigns to set their baseline characterization spectrally, spatially, temporally and radiometrically. The knowledge from building and calibrating these instruments has aided in technology advancements and the need for developing more accurate instruments has increased. In order to understand the on-ground instrument performance, NASA Langley Research Center has partnered with the Thermal Radiation Group of Virginia Tech to develop a first-principle, dynamic, electrothermal, numerical model of a scanning radiometer that can be used to enhance the interpretation of an Earth radiation budget-like instrument on orbit. This paper will summarize the current efforts of developing this high-fidelity end-to-end model and also highlight how it can be applied to an Earth radiation budget instrument.

Ashraf, Anum Rauf Bark

Numerical Modeling of No Vent Filling of a Cryogenic Tank

This paper presents a multi-node finite volume model of the chilling and filling of a cryogenic tank using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code. Vented Chill and No Vent Fill (VCNVF) tests were conducted at Marshall Space Flight Center where a flight-like tank was filled with liquid nitrogen from a supply tank. In the VCNVF tests, the tank was partially chilled while the vent valve was open. After the partial chilling, the vent valve was closed and the tank was filled without any venting. An integrated numerical model of the test set up was developed. The model included the transfer line from the supply tank, the target tank with the spray nozzle and solid walls, and the discharge line with the vent valve. The tank was discretized into multiple fluid nodes and branches to represent the ullage and liquid nitrogen and multiple solid nodes to represent the tank wall and structure. The heat transfer between solid to fluid was calculated from pool boiling correlations which include film, transition, and nucleate boiling, as well as natural convection during pre- and post-boiling. The model also accounts for the condensation of vapor in the tank when it comes in contact with the liquid spray. The predicted pressure, resident mass, wall and ullage temperature in the tank were compared with the test data.

Alok Majumdar

Numerical Modeling of No Vent Filling of a Cryogenic Tank

This paper presents a multi-node finite volume model of the chilling and filling of a cryogenic tank using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code. Vented Chill and No Vent Fill (VCNVF) tests were conducted at Marshall Space Flight Center where a flight-like tank was filled with liquid nitrogen from a supply tank. In the VCNVF tests, the tank was partially chilled while the vent valve was open. After the partial chilling, the vent valve was closed and the tank was filled without any venting. An integrated numerical model of the test set up was developed. The model included the transfer line from the supply tank, the target tank with the spray nozzle and solid walls, and the discharge line with the vent valve. The tank was discretized into multiple fluid nodes and branches to represent the ullage and liquid nitrogen and multiple solid nodes to represent the tank wall and structure. The heat transfer between solid to fluid was calculated from pool boiling correlations which include film, transition, and nucleate boiling, as well as natural convection during pre- and post-boiling. The model also accounts for the condensation of vapor in the tank when it comes in contact with the liquid spray. The predicted pressure, resident mass, wall and ullage temperature in the tank were compared with the test data.

Alok Majumdar

Numerical Modeling of Ocular Dysfunction in Space

Upon introduction to microgravity, the near-loss of hydrostatic pressure causes a marked cephalic (headward) shift of fluid in an astronaut's body. The fluid shift, along with other factors of spaceflight, induces a cascade of interdependent physiological responses which occur at varying time scales. Long-duration missions carry an increased risk for the development of the Visual Impairment and Intracranial Pressure (VIIP) syndrome, a spectrum of ophthalmic changes including posterior globe flattening, choroidal folds, distension of the optic nerve sheath, kinking of the optic nerve and potentially permanent degradation of visual function. In the cases of VIIP found to date, the initial onset of symptoms occurred after several weeks to several months of spaceflight, by which time the gross bodily fluid distribution is well established. We are developing a suite of numerical models to simulate the effects of fluid shift on the cardiovascular, central nervous and ocular systems. These models calculate the modified mean volumes, flow rates and pressures that are characteristic of the altered quasi-homeostatic state in microgravity, including intracranial and intraocular pressures. The results of the lumped models provide initial and boundary data to a 3D finite element biomechanics simulation of the globe, optic nerve head and retrobulbar subarachnoid space. The integrated set of models will be used to investigate the evolution of the biomechanical stress state in the ocular tissues due to long-term exposure to microgravity.

intraocular pressure

Numerical Modeling of No Vent Filling of a Cryogenic Tank with Thermo-dynamic Vent System Assisted Injector

This paper presents a multi-node finite volume model of No Vent Filling (NVF) of a cryogenic tank with Thermo-dynamic Vent System (TVS) assisted injector, using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code. NVF tests were conducted in a CRYOgenic Orbital Testbed (CRYOTE) tank which was filled by liquid nitrogen using a TVS assisted injector with vent valve closed during the entire filling process. In TVS assisted injector, the liquid flow splits into two streams: one stream is routed to a Joule-Thomson (J-T) orifice where the flow immediately flashes from liquid to vapor or two phase mixture because the downstream to J-T leg is maintained at vacuum level.; the other stream is injected into the tank after being cooled by cold vapor of J-T leg in a heat exchanger. The flow through J-T leg is also used to cool the outer metal matrix of the injector which in turn cools the vapor in the tank ullage. The cooling of vapor in the ullage by cold injector surface reduces ullage pressure that allows liquid to enter and fill the tank. An integrated numerical model of the test set up was developed. The model included the dual lines of TVS assisted injector, tank and tank wall. The tank was discretized into multiple fluid nodes and branches to represent the ullage and liquid nitrogen and multiple solid nodes to represent the tank wall and structure. The heat transfer between solid to fluid was calculated from pool boiling correlations which include film, transition, and nucleate boiling, as well as natural convection during pre- and post-boiling. The model also accounts for the condensation of vapor at ullage-injector interface and when it comes in contact with the liquid spray. The predicted pressure, resident mass, wall temperature in the tank were compared with the test data.

Alok Majumdar

Numerical models of mantle convection

An overview of numerical methods describing the structure and dynamics of the mantle is presented with attention given to novel 3D modeling techniques. The paper reviews 3D spherical and Cartesian models for constant viscosity emphasizing the assumptions regarding style of convection, time dependence, and implications for the mantle. Similarly treated are 3D Cartesian models with temperature-dependent viscosities, and briefly examined are models that are based on compressibility, nonlinear viscosity, or plates. Extensive illustrations are presented detailing: (1) temperature variations from models of 3D thermal convection in spherical shells; (2) thermal anomalies in equatorial cross sections; and (3) temperature variations in a spherical shell heated from within. The discussion relates the numerical results of the models with real mantle-convection events, and the simulations are shown to yield increasingly realistic representations of material behavior.

Schubert, G.

Experimental, physical and numerical modeling of lunar regolith and lunar regolith structures

An extensive series of laboratory strength and deformation experiments have been performed on a lunar regolith simulant. Results of these experiments are compared to results from experiments on real lunar regolith from the Apollo and Luna missions to illustrate the suitability of this material in capturing the engineering properties lunar regolith. In addition, these results are used to calibrate a constitutive model used to describe its stress-strain behavior. This model, in conjunction with numerical analysis techniques, is used to predict the response i.e. material parameters) of lunar simulant under 1/6-g and low confining stress conditions. These tools are also used to predict the displacement response of a lunar soil embankment structure used to cover a first generation human habitat module, which might be used to accommodate the first astronauts revisiting the moon. These predictions are compared to physical models of this structure, which are tested in a geotechnical centrifuge in order to satisfy scaling relationships between prototype and model.

Perkins, Steven W.

Thermal Pollution Mathematical Model. Volume 3: User's Manual for One-Dimensional Numerical Model for the Seasonal Thermocline

A user's manual for a one dimensional thermal model to predict the temperature profiles of a deep body of water for any number of annual cycles is presented. The model is essentially a set of partial differential equations which are solved by finite difference methods using a high speed digital computer. The model features the effects of area change with depth, nonlinear interaction of wind generated turbulence and buoyancy, adsorption of radiative heat flux below the surface, thermal discharges, and the effects of vertical convection caused by discharge. The main assumption in the formulation is horizontal homogeneity. The environmental impact of thermal discharges from power plants is emphasized. Although the model is applicable to most lakes, a specific site (Lake Keowee, S.C.) application is described in detail. The programs are written in FORTRAN 5.

Lee, S. S.

Numerical Modeling of Ablation Heat Transfer

A unique numerical method has been developed for solving one-dimensional ablation heat transfer problems. This paper provides a comprehensive description of the method, along with detailed derivations of the governing equations. This methodology supports solutions for traditional ablation modeling including such effects as heat transfer, material decomposition, pyrolysis gas permeation and heat exchange, and thermochemical surface erosion. The numerical scheme utilizes a control-volume approach with a variable grid to account for surface movement. This method directly supports implementation of nontraditional models such as material swelling and mechanical erosion, extending capabilities for modeling complex ablation phenomena. Verifications of the numerical implementation are provided using analytical solutions, code comparisons, and the method of manufactured solutions. These verifications are used to demonstrate solution accuracy and proper error convergence rates. A simple demonstration of a mechanical erosion (spallation) model is also provided to illustrate the unique capabilities of the method.

Ewing, Mark E.

Numerical Modeling of Glaciers in Martian Paleoclimates

Numerous geologic features suggest the presence of ice flow on the surface of mars. These features include lobate debris aprons, concentric crater fill, and lineated valley fill. The lateral extent of these features can range from 100 meters to over 20 km. Previous work has demonstrated that these features could not have formed in current Martian conditions. It has long been speculated that changes in Mars orbital properties, namely its obliquity, eccentricity, and argument of perihelion, can result in dramatic changes to climate. Recent climate model studies have shown that at periods of increased obliquity north polar water ice is mobilized southward and deposited at low ad mid latitudes. Mid latitude accumulation of ice would provide the necessary conditions for rock glaciers to form. A time-marching, finite element glacier model is used to demonstrate the ability of ice and ice-rock mixtures to flow under Martian paleoclimate conditions. Input to this model is constrained by the NASA Ames Mars General Circulation Model (MGCM).

Colaprete, A.

Numerical modeling of transonic juncture flow

A numerical investigation of the interaction between a wind tunnel sidewall boundary layer and a thin low-aspect-ratio wing has been performed for transonic speeds and flight Reynolds numbers. A three-dimensional Navier-Stokes code was applied to calculate the flowfields. The results indicated that the sidewall boundary layer had a strong influence on the flowfield around the wing. The computed wing pressure distributions showed vast improvements over previous free-air computations, and were in excellent agreement with experimental data. The low momentum of the sidewall boundary layer resulted in higher pressures in the juncture region, which decreased the favorable spanwise pressure gradient. This significantly decreased the spanwise migration of the wing boundary layer. Weak vortices were predicted in both the upper and lower surface juncture regions. These vortices are believed to have been generated by lateral skewing of the streamlines in the approaching boundary layer.

Milholen, William E., II

A global time-dependent model of thunderstorm electricity. I - Mathematical properties of the physical and numerical models

A time-dependent model is introduced that can be used to simulate the interaction of a thunderstorm with its global electrical environment. The model solves the continuity equation of the Maxwell current, which is assumed to be composed of the conduction, displacement, and source currents. Boundary conditions which can be used in conjunction with the continuity equation to form a well-posed initial-boundary value problem are determined. Properties of various components of solutions of the initial-boundary value problem are analytically determined. The results indicate that the problem has two time scales, one determined by the background electrical conductivity and the other by the time variation of the source function. A numerical method for obtaining quantitative results is introduced, and its properties are studied. Some simulation results on the evolution of the displacement and conduction currents during the electrification of a storm are presented.

Browning, G. L.

Combined Modeling of Acceleration, Transport, and Hydrodynamic Response in Solar Flares: The Numerical Model - 1

Acceleration and transport of high-energy particles and fluid dynamics of atmospheric plasma are interrelated aspects of solar flares, but for convenience and simplicity they were artificially separated in the past. We present here self consistently combined Fokker-Planck modeling of particles and hydrodynamic simulation of flare plasma. Energetic electrons are modeled with the Stanford unified code of acceleration, transport, and radiation, while plasma is modeled with the Naval Research Laboratory flux tube code. We calculated the collisional heating rate directly from the particle transport code, which is more accurate than those in previous studies based on approximate analytical solutions. We repeated the simulation of Mariska et al. with an injection of power law, downward-beamed electrons using the new heating rate. For this case, a -10% difference was found from their old result. We also used a more realistic spectrum of injected electrons provided by the stochastic acceleration model, which has a smooth transition from a quasi-thermal background at low energies to a non thermal tail at high energies. The inclusion of low-energy electrons results in relatively more heating in the corona (versus chromosphere) and thus a larger downward heat conduction flux. The interplay of electron heating, conduction, and radiative loss leads to stronger chromospheric evaporation than obtained in previous studies, which had a deficit in low-energy electrons due to an arbitrarily assumed low-energy cutoff. The energy and spatial distributions of energetic electrons and bremsstrahlung photons bear signatures of the changing density distribution caused by chromospheric evaporation. In particular, the density jump at the evaporation front gives rise to enhanced emission, which, in principle, can be imaged by X-ray telescopes. This model can be applied to investigate a variety of high-energy processes in solar, space, and astrophysical plasmas.

STOCHASTIC ACCELERATION MODEL

Use of simulated satellite radiances from a mesoscale numerical model to understand kinematic and dynamic processes

Techniques by which mesoscale numerical weather prediction model output and radiative transfer codes are combined to simulate the radiance fields that a given passive temperature/moisture satellite sensor would see if viewing the evolving model atmosphere are introduced. The goals are to diagnose the dynamical atmospheric processes responsible for recurring patterns in observed satellite radiance fields, and to develop techniques to anticipate the ability of satellite sensor systems to depict atmospheric structures and provide information useful for numerical weather prediction (NWP). The concept of linking radiative transfer and dynamical NWP codes is demonstrated with time sequences of simulated radiance imagery in the 24 TIROS vertical sounder channels derived from model integrations for March 6, 1982.

Kalb, Michael