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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

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.↗
NASA NTRS1992

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.↗
NASA NTRS1980

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.↗
NASA NTRS2013

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.↗
NASA NTRS2004

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.↗
NASA NTRS1992

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↗
NASA NTRS1987

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.↗
NASA NTRS2009

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↗
NASA NTRS1987

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↗
DOE OSTI2021

Numerical Modeling of CO2 Sequestration within a Five-Spot Well Pattern in the Morrow B Sandstone of the Farnsworth Hydrocarbon Field: Comparison of the TOUGHREACT, STOMP-EOR, and GEM Simulators

The objectives of this study were (1) to assess the fate and impact of CO2 injected into the Morrow B Sandstone in the Farnsworth Unit (FWU) through numerical non-isothermal reactive transport modeling, and (2) to compare the performance of three major reactive solute transport simulators, TOUGHREACT, STOMP-EOR, and GEM, under the same input conditions. The models were based on a quarter of a five-spot well pattern where CO2 was injected on a water-alternating-gas schedule for the first 25 years of the 1000 year simulation. The reservoir pore fluid consisted of water with or without petroleum. The results of the models have numerous broad similarities, such as the pattern of reservoir cooling caused by the injected fluids, a large initial pH drop followed by gradual pH neutralization, the long-term persistence of an immiscible CO2 gas phase, the continuous dissolution of calcite, very small decreases in porosity, and the increasing importance over time of carbonate mineral CO2 sequestration. The models differed in their predicted fluid pressure evolutions; amounts of mineral precipitation and dissolution; and distribution of CO2 among immiscible gas, petroleum, formation water, and carbonate minerals. The results of the study show the usefulness of numerical simulations in identifying broad patterns of behavior associated with CO2 injection, but also point to significant uncertainties in the numerical values of many model output parameters.

STOMP↗
NASA NTRS2007

Numerical Modeling of Propellant Boiloff in Cryogenic Storage Tank

This Technical Memorandum (TM) describes the thermal modeling effort undertaken at Marshall Space Flight Center to support the Cryogenic Test Laboratory at Kennedy Space Center (KSC) for a study of insulation materials for cryogenic tanks in order to reduce propellant boiloff during long-term storage. The Generalized Fluid System Simulation program has been used to model boiloff in 1,000-L demonstration tanks built for testing the thermal performance of glass bubbles and perlite insulation. Numerical predictions of boiloff rate and ullage temperature have been compared with the measured data from the testing of demonstration tanks. A satisfactory comparison between measured and predicted data has been observed for both liquid nitrogen and hydrogen tests. Based on the experience gained with the modeling of the demonstration tanks, a numerical model of the liquid hydrogen storage tank at launch complex 39 at KSC was built. The predicted boiloff rate of hydrogen has been found to be in good agreement with observed field data. This TM describes three different models that have been developed during this period of study (March 2005 to June 2006), comparisons with test data, and results of parametric studies.

Majumdar, A. K.↗
DOE OSTI2022

Numerical modeling of ion transport and adsorption in porous media: A pore-scale study for capacitive deionization desalination

In this report a pore-scale model is presented to simulate the dynamic ion transport and adsorption processes in porous electrodes used for capacitive deionization (CDI). The Stokes equation governing water flow is solved using the lattice Boltzmann method and Nernst-Planck equation describing ion transport is solved using the finite volume method. The ion adsorption process is considered at the surface of carbon electrodes. After validation against analytical solutions and published results, the model is used to study the coupled water flow, ion transport and adsorption in both two-dimensional and three-dimensional porous CDI electrodes at the pore scale. The effect of electrode microstructure, electrical potential and flow velocity on the adsorption processes is quantitatively investigated, and the relative importance of various parameters is determined. The presented model can be a powerful numerical tool to quantitatively analyze ion transport and adsorption in porous electrodes, and may provide useful information for the design and optimization of CDI electrodes and operating conditions for desired desalination efficiency, water throughput and cost.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1985

Current status of one- and two-dimensional numerical models: Successes and limitations

The capabilities of one and two-dimensional numerical solar cell modeling programs (SCAP1D and SCAP2D) are described. The occasions when a two-dimensional model is required are discussed. The application of the models to design, analysis, and prediction are presented along with a discussion of problem areas for solar cell modeling.

Schwartz, R. J.↗
NASA NTRS1992

Analytic and numerical models of wave-CISK with conditional heating

Wave-CISK with conditional heating is investigated in the equatorial zonal-height plane by analytic and numerical techniques. For two- and three-level models, previous results are extended to give additional evidence that the most unstable mode has a single wet region of infinitesimal width. A three-level model has qualitatively similar behavior as the two-level model except that propagating solutions are possible due to coalescence of internal vertical modes. Phase speeds with conditional heating are found to be slightly greater than those for unconditional heating. The structure has one circulation cell in the vertical and is asymmetric in longitude with stronger motion on the leading edge. Growth rate is inversely proportional to the width of the single wet region. That width can be limited by second-order diffusion. A general integral relationship between growth rate, viscosity, phase speed, and heating is derived. The main conclusion is that the linear wave-CISK catastrophe is modified by conditional heating but not eliminated. The preferred mode of instability has one wet region, but it occurs on the smallest possible scale. It is likely that numerical models that use conditional heating are sensitive to resolution, especially for the commonly used spectral truncations, unless there is sufficiently strong damping at the smallest scales.

Crum, Francis X.↗
DOE OSTI2023

Numerical modeling of hydrogen mixing in a direct-injection engine fueled with gaseous hydrogen

Hydrogen is considered as one of the most promising options to achieve effective decarbonization of the energy and transportation sectors. As such, it has recently been receiving increasing attention because of its promising potential as an energy carrier for advanced energy and propulsion systems. With a focus on internal combustion engines, direct injection (DI) of gaseous hydrogen during the compression stroke offers great potential for high engine efficiency and specific power while reducing the risk of backfiring and pre-ignition issues. Therefore, many experimental and numerical efforts have recently been dedicated to understanding the physical and chemical behaviors of hydrogen in engine during mixing and combustion. This study focuses on computational fluid dynamics (CFD) modeling of the hydrogen DI process in a hydrogen optical research engine. Under the conditions studied, gaseous hydrogen is injected into the combustion chamber via a centrally located single-hole injector at a pressure of 100 bar. Two configurations, namely low-and high-tumble, are investigated to understand the impact of different in-cylinder flow patterns on the fuel-air mixture preparation. Simulations are carried out using the commercial CFD software CONVERGE. Here, the in-cylinder turbulence is modeled with an unsteady Reynolds-averaged Navier-Stokes (URANS) formulation closed by the renormalization group (RNG) k-ε model. Several numerical methods and model constants, including but not limited to turbulent Schmidt number, are evaluated. The numerical results are systematically compared against experimental measurements of velocity and hydrogen concentration fields on the vertical center plane to assess the performance of the CFD model, unveil the physics of hydrogen mixing, and establish best practices for modeling hydrogen DI under relatively high injection pressure conditions.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1989

Numerical modeling of a vortex stabilized arcjet thruster

A numerical method to solve the equations governing a vortex stabilized arcjet thruster is being developed. The model will allow the effects of swirling flow and geometry on arcjet thruster performance to be determined. The propellant flow equations, i.e., the axisymmetric, thin layer, Navier-Stokes equations, are solved using a Gauss-Siedel line-relaxation procedure. An implicit FTCS method is used to solve the electromagnetic field equations. A grid generation scheme was developed for an arbitrary arcjet geometry. The model will allow the radial and axial components of velocity and current distributions to be determined from a region upstream of the cathode, through the constrictor, to the exit plane of the nozzle.

Keith, Theo G., Jr.↗