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.

SourceResearch connectionFrom yearThrough year

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

At least 217 records · Page 12

All sources NASA NTRS DOE OSTI

NASA NTRS1996

Numerical modelling of the circulation of superrotating atmospheres: Venus and Titan

Atmospheric superrotation is known to exist on Venus where the cloud cover, at 70 km above the surface, rotates about 60 times faster than the solid planet. Superrotation is strongly suspected on Titan on the basis of numerical studies and observations. The numerical simulation of superrotation is described. It is shown how the general circulation models, originally developed for terrestrial meteorology and climatology, can be used to infer the circulation in such atmospheres. The emphasis is on the case of Titan for which a specific general circulation model was developed.

Hourding, F.↗

NASA NTRS1992

Numerical models for the collapse and fragmentation of centrally condensed molecular cloud cores

The gravitational collapse and fragmentation of centrally condensed molecular cloud cores are investigated using a new hydrodynamical code. The numerical scheme is second-order accurate and uses explicit finite difference methods to advance the fluid variables on a 3D Cartesian grid. Two initial power-law density profiles, rho varies as r exp -1 and rho varies as r exp -2, are considered, as well as two initial density perturbations in the azimuthal coordinate theta, rho-i - rho(1 + a cos 2theta) where a = 0.1 and 0.5. Fragmentation is found to be possible in these centrally condensed cores if the initial conditions also include differential rotation. Models which collapse with initial uniform rotation do not produce fragments. If molecular cloud cores are indeed centrally condensed, as suggested by observations of star-forming regions and by studies of ambipolar diffusion, then differential rotation may be a mechanism for producing binary protostars during gravitational collapse.

Myhill, Elizabeth A.↗

NASA NTRS1990

Two-dimensional atmospheric transport and chemistry model - Numerical experiments with a new advection algorithm

Extensive testing of the advective scheme proposed by Prather (1986) has been carried out in support of the California Institute of Technology-Jet Propulsion Laboratory two-dimensional model of the middle atmosphere. The original scheme is generalized to include higher-order moments. In addition, it is shown how well the scheme works in the presence of chemistry as well as eddy diffusion. Six types of numerical experiments including simple clock motion and pure advection in two dimensions have been investigated in detail. By comparison with analytic solutions, it is shown that the new algorithm can faithfully preserve concentration profiles, has essentially no numerical diffusion, and is superior to a typical fourth-order finite difference scheme.

Shia, Run-Lie↗

NASA NTRS2001

Numerical Modeling of Liquid-Vapor Phase Change

We implemented a two- and three-dimensional finite difference/front tracking technique to solve liquid-vapor phase change problems. The mathematical and the numerical features of the method were explained in great detail in our previous reports, Briefly, we used a single formula representation which incorporated jump conditions into the governing equations. The interfacial terms were distributed as singular terms using delta functions so that the governing equations would be the same as conventional conservation equations away from the interface and in the vicinity of the interface they would provide correct jump conditions. We used a fixed staggered grid to discretize these equations and an unstructured grid to explicitly track the front. While in two dimensions the front was simply a connection of small line segments, in three dimensions it was represented by a connection of small triangular elements. The equations were written in conservative forms and during the course of computations we used regriding to control the size of the elements of the unstructured grid. Moreover, we implemented a coalescence in two dimensions which allowed the merging of different fronts or two segments of the same front when they were sufficiently close. We used our code to study thermocapillary migration of bubbles, burst of bubbles at a free surface, buoyancy-driven interactions of bubbles, evaporation of drops, rapid evaporation of an interface, planar solidification of an undercooled melt, dendritic solidification, and a host of other problems cited in the reference.

Esmaeeli, Asghar↗

NASA NTRS1996

Numerical Modeling of Spray Combustion with an Unstructured-Grid Method

The present unstructured-grid method follows strictly the basic finite volume forms of the conservation laws of the governing equations for the entire flow domain. High-order spatially accurate formulation has been employed for the numerical solutions of the Navier-Stokes equations. A two-equation k-epsilon turbulence model is also incorporated in the unstructured-grid solver. The convergence of the resulted linear algebraic equation is accelerated with preconditioned Conjugate Gradient method. A statistical spray combustion model has been incorporated into the present unstructured-grid solver. In this model, spray is represented by discrete particles, rather than by continuous distributions. A finite number of computational particles are used to predict a sample of total population of particles. Particle trajectories are integrated using their momentum and motion equations and particles exchange mass, momentum and energy with the gas within the computational cell in which they are located. The interaction calculations are performed simultaneously and eliminate global iteration for the two-phase momentum exchange. A transient spray flame in a high pressure combustion chamber is predicted and then the solution of liquid-fuel combusting flow with a rotating cup atomizer is presented and compared with the experimental data. The major conclusion of this investigation is that the unstructured-grid method can be employed to study very complicated flow fields of turbulent spray combustion. Grid adaptation can be easily achieved in any flow domain such as droplet evaporation and combustion zone. Future applications of the present model can be found in the full three-dimensional study of flow fields of gas turbine and liquid propulsion engine combustion chambers with multi-injectors.

Shang, H. M.↗

NASA NTRS1989

Numerical modeling of turbulent supersonic reacting coaxial jets

The paper considers the mixing and subsequent combustion within turbulent reacting shear layers. A computer program was developed to solve the axisymmetric Reynolds averaged Navier-Stokes equations. The numerical method integrates the Reynolds averaged Navier-Stokes equations using a finite volume approach while advancing the solution forward in time using a Runge-Kutta scheme. Three separate flowfields are investigated and it is found that no single turbulence model considered could accurately predict the degree of mixing for all three cases.

Eklund, Dean R.↗

NASA NTRS1995

Numerical modeling of flame-balls in fuel-air mixtures

At low gravity, when buoyancy effects are small, flame-balls can be generated. These are stationary spherical structures whose existence appears to require a near-limit mixture, a small Lewis number and heat losses from radiation. It is our goal to combine computational modeling with existing experimental and theoretical studies (NASA) of these structures so that an improved understanding of flammability limits and near-limit phenomena will occur. The question of flammability limits is of fundamental importance and has long been examined. It is of great practical importance to predict, from first principles, a limit mixture strength that agrees with experimental values for the configuration at hand. Flame-balls provide an excellent configuration in which convective losses can be eliminated and the resulting stable solutions are produced from a diffusive, reactive and radiative balance. Although analytical modeling provides convincing evidence that the key physical ingredients of flame-balls have been identified, quantitative confirmation can only come from detailed numerical simulations. Our goal is to predict theoretically the mass fractions of the species and the temperature as functions of the independent coordinate r.

Smooke, Mitchell D.↗

NASA NTRS1993

Numerical Modeling of Anti-icing Systems and Comparison to Test Results on a NACA 0012 Airfoil

A series of experimental tests were conducted in the NASA Lewis IRT on an electro-thermally heated NACA 0012 airfoil. Quantitative comparisons between the experimental results and those predicted by a computer simulation code were made to assess the validity of a recently developed anti-icing model. An infrared camera was utilized to scan the instantaneous temperature contours of the skin surface. Despite some experimental difficulties, good agreement between the numerical predictions and the experiment results were generally obtained for the surface temperature and the possibility for each runback to freeze. Some recommendations were given for an efficient operation of a thermal anti-icing system.

Al-Khalil, Kamel M.↗

NASA NTRS1993

Numerical modeling of anti-icing systems and comparison to test results on a NACA 0012 airfoil

A series of experimental tests were conducted in the NASA Lewis IRT on an electro-thermally heated NACA 0012 airfoil. Quantitative comparisons between the experimental results and those predicted by a computer simulation code were made to assess the validity of a recently developed anti-icing model. An infrared camera was utilized to scan the instantaneous temperature contours of the skin surface. Despite some experimental difficulties, good agreement between the numerical predictions and the experimental results were generally obtained for the surface temperature and the possibility for the runback to freeze. Some recommendations were given for an efficient operation of a thermal anti-icing system.

Al-Khalil, Kamel M.↗

NASA NTRS1982

Numerical model of a two-dimensional, non-plane transient magnetohydrodynamic flow

The equations describing two-dimensional three-component magnetohydrodynamic (MHD) transient flows are formulated for a system of spherical coordinates. With the numerical code based on Implicit Continuous Fluid Eulerian (ICE) scheme, MHD flows resulting from a sudden energy release in a stratified medium are examined. Because of the inclusion of out-of-plane components of velocity and magnetic fields, MHD transverse waves are observed in addition to fast, slow and entropy waves. Numerical results for compressible MHD shocks are found in satisfactory agreement with the theoretical predictions.

Han, S. M.↗

NASA NTRS1993

Analytical and Numerical Modeling of Flame-Balls in Hydrogen-Air Mixtures

Flame-balls are stationary spherical premixed flames observed in certain near-limit mixtures. It is believed that radiative heat losses are an important stabilizing influence. Numerical solutions of flame balls are constructed for hydrogen-air mixtures using an accurate description of the chemical kinetics, diffusive transport, and radiation losses. A lean limit equivalence ratio of 0.0866 is predicted and a rich limit of 2,828. For any equivalence ratio between the two limits there are two solutions. One is characterized by a small flame, incomplete reactant consumption, and negligible radiation losses. The other by a large flame, complete consumption of one of the reactants, and significant radiation losses. The maximum temperature varies between 1200 and 900 K as the two solution branches are traversed. Much of our discussion is a reprise and modification of previously published analytical results, for these provide physical insight into the nature of the solutions, and suggest that a portion of the large flame branch near the lean limit is stable and so corresponds to observable flames.

Buckmaster, J.↗

NASA NTRS2005

Numerical Modelling of the Deep Impact Mission Experiment

NASA s Deep Impact Mission (launched January 2005) will provide, for the first time ever, insights into the interior of a comet (Tempel 1) by shooting a approx.370 kg projectile onto the surface of a comets nucleus. Although it is usually assumed that comets consist of a very porous mixture of water ice and rock, little is known about the internal structure and in particular the constitutive material properties of a comet. It is therefore difficult to predict the dimensions of the excavated crater. Estimates of the crater size are based on laboratory experiments of impacts into various target compositions of different densities and porosities using appropriate scaling laws; they range between 10 s of meters up to ~250 m in diameter [1]. The size of the crater depends mainly on the physical process(es) that govern formation: Smaller sizes are expected if (1) strength, rather than gravity, limits crater growth; and, perhaps even more crucially, if (2) internal energy losses by pore-space collapse reduce the coupling efficiency (compaction craters). To investigate the effect of pore space collapse and strength of the target we conducted a suite of numerical experiments and implemented a novel approach for modeling porosity and the compaction of pores in hydrocode calculations.

Wuennemann, K.↗

NASA NTRS2008

Numerically Modeling the Erosion of Lunar Soil by Rocket Exhaust Plumes

In preparation for the Apollo program, Leonard Roberts of the NASA Langley Research Center developed a remarkable analytical theory that predicts the blowing of lunar soil and dust beneath a rocket exhaust plume. Roberts assumed that the erosion rate was determined by the excess shear stress in the gas (the amount of shear stress greater than what causes grains to roll). The acceleration of particles to their final velocity in the gas consumes a portion of the shear stress. The erosion rate continues to increase until the excess shear stress is exactly consumed, thus determining the erosion rate. Roberts calculated the largest and smallest particles that could be eroded based on forces at the particle scale, but the erosion rate equation assumed that only one particle size existed in the soil. He assumed that particle ejection angles were determined entirely by the shape of the terrain, which acts like a ballistic ramp, with the particle aerodynamics being negligible. The predicted erosion rate and the upper limit of particle size appeared to be within an order of magnitude of small-scale terrestrial experiments but could not be tested more quantitatively at the time. The lower limit of particle size and the predictions of ejection angle were not tested. We observed in the Apollo landing videos that the ejection angles of particles streaming out from individual craters were time-varying and correlated to the Lunar Module thrust, thus implying that particle aerodynamics dominate. We modified Roberts theory in two ways. First, we used ad hoc the ejection angles measured in the Apollo landing videos, in lieu of developing a more sophisticated method. Second, we integrated Roberts equations over the lunar-particle size distribution and obtained a compact expression that could be implemented in a numerical code. We also added a material damage model that predicts the number and size of divots which the impinging particles will cause in hardware surrounding the landing rocket. Then, we performed a long-range ballistics analysis for the ejected particulates.

Source record↗

NASA NTRS1988

Mesoscale spiral vortex embedded within a Lake Michigan snow squall band - High resolution satellite observations and numerical model simulations

It is known that Great Lakes snow squall convection occurs in a variety of different modes depending on various factors such as air-water temperature contrast, boundary-layer wind shear, and geostrophic wind direction. An exceptional and often neglected source of data for mesoscale cloud studies is the ultrahigh resolution multispectral data produced by Landsat satellites. On October 19, 1972, a clearly defined spiral vortex was noted in a Landsat-1 image near the southern end of Lake Michigan during an exceptionally early cold air outbreak over a still very warm lake. In a numerical simulation using a three-dimensional Eulerian hydrostatic primitive equation mesoscale model with an initially uniform wind field, a definite analog to the observed vortex was generated. This suggests that intense surface heating can be a principal cause in the development of a low-level mesoscale vortex.

Lyons, Walter A.↗

NASA NTRS1986

Optimization methods and silicon solar cell numerical models

The goal of this project is the development of an optimization algorithm for use with a solar cell model. It is possible to simultaneously vary design variables such as impurity concentrations, front junction depth, back junctions depth, and cell thickness to maximize the predicted cell efficiency. An optimization algorithm has been developed and interfaced with the Solar Cell Analysis Program in 1 Dimension (SCAPID). SCAPID uses finite difference methods to solve the differential equations which, along with several relations from the physics of semiconductors, describe mathematically the operation of a solar cell. A major obstacle is that the numerical methods used in SCAPID require a significant amount of computer time, and during an optimization the model is called iteratively until the design variables converge to the value associated with the maximum efficiency. This problem has been alleviated by designing an optimization code specifically for use with numerically intensive simulations, to reduce the number of times the efficiency has to be calculated to achieve convergence to the optimal solution. Adapting SCAPID so that it could be called iteratively by the optimization code provided another means of reducing the cpu time required to complete an optimization. Instead of calculating the entire I-V curve, as is usually done in SCAPID, only the efficiency is calculated (maximum power voltage and current) and the solution from previous calculations is used to initiate the next solution.

Girardini, K.↗

NASA NTRS2003

Numerical Modeling of Unsteady Thermofluid Dynamics in Cryogenic Systems

A finite volume based network analysis procedure has been applied to model unsteady flow without and with heat transfer. Liquid has been modeled as compressible fluid where the compressibility factor is computed from the equation of state for a real fluid. The modeling approach recognizes that the pressure oscillation is linked with the variation of the compressibility factor; therefore, the speed of sound does not explicitly appear in the governing equations. The numerical results of chilldown process also suggest that the flow and heat transfer are strongly coupled. This is evident by observing that the mass flow rate during 90-second chilldown process increases by factor of ten.

Majumdar, Alok↗

NASA NTRS1993

Numerical models of jet disruption in cluster cooling flows

We present a coherent picture for the formation of the observed diverse radio morphological structures in dominant cluster galaxies based on the jet Mach number. Realistic, supersonic, steady-state cooling flow atmospheres are evolved numerically and then used as the ambient medium through which jets of various properties are propagated. Low Mach number jets effectively stagnate due to the ram pressure of the cooling flow atmosphere while medium Mach number jets become unstable and disrupt in the cooling flow to form amorphous structures. High Mach number jets manage to avoid disruption and are able to propagate through the cooling flow.

Loken, Chris↗

NASA NTRS2005

Shock Melting of Permafrost on Mars: Water Ice Multiphase Equation of State for Numerical Modeling and Its Testing

The presence of water/ice/brine in upper layers of Martian crust affects many processes of impact cratering. Modeling of these effects promises better understanding of Martian cratering records. We present here the new ANEOS-based multiphase equation of state for water/ice constructed for usage in hydrocodes and first numerical experiments on permafrost shock melting. Preliminary results show that due to multiple shock compression of ice inclusions in rocks the entropy jump in shocked ice is smaller than in pure ice for the same shock pressure. Hence previous estimates of ice melting during impact cratering on Mars should be re-evaluated. Additional information is included in the original extended abstract.

Ivanov, B. A.↗