Search NASASearch

SEARCH · Search NASA

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 19 records

Numerical modeling of diffusive-convective physical vapor transport in cylindrical vertical ampoules

Diffusive-convective physical vapor transport (PVT) in cylindrical, vertical ampoules of aspect ratio (length/radius) between 0.5 and 10 was modeled numerically. The transport of a crystal forming component through an inert component that undergoes zero net transport was considered. Systems were treated in which: (1) with unequal molecular weight of the components and with temperature gradients typically employed in PVT, convective flow arises dominantly from solutal density gradients; and (2) with equal molecular weight of the components, convective flow can arise only from thermal expansion. It was found that, due to the diffusion-induced horizontal density gradients, buoyancy-driven convective flows are superimposed on the diffusive-advective fluxes without threshold. Net recirculation sets in adjacent to the growing interface, in contrast to the corresponding monocomponent situation where marginally stable convective modes fill the whole fluid space. Depending on the orientation of the main transport direction with respect to gravity, convection can either reduce or enhance the diffusion-induced radial concentration gradients. Significant enhancement of the net transport rate was found to occur only when the whole vapor space between source and growing crystal is filled by a convective recirculation roll. Solutal and thermal convection results are similar; yet for quantitative discussions, thermal and solutal Rayleigh numbers are not interchangeable in contrast to convective situations that lack net mass transport across the fluid space.

Markham, B. L.

Thermal Pollution Mathematical Model: Verification of One-Dimensional Numerical Model at Lake Keowee - Volume 2

A one dimensional model for studying the thermal dynamics of cooling lakes was developed and verified. The model is essentially a set of partial differential equations which are solved by finite difference methods. The model includes the effects of variation of area with depth, surface heating due to solar radiation absorbed at the upper layer, and internal heating due to the transmission of solar radiation to the sub-surface layers. The exchange of mechanical energy between the lake and the atmosphere is included through the coupling of thermal diffusivity and wind speed. The effects of discharge and intake by power plants are also included. The numerical model was calibrated by applying it to Cayuga Lake. The model was then verified through a long term simulation using Lake Keowee data base. The comparison between measured and predicted vertical temperature profiles for the nine years is good. The physical limnology of Lake Keowee is presented through a set of graphical representations of the measured data base.

Lee, S. S.

A mesoscale sixth-order numerical modelling system

A numerical simulation system is currently under development for NASA which is intended to improve the modeling of subsynoptic and mesoscale adjustments associated with cyclogenesis, severe storm development and atmospheric transport processes. The model utilizes a standard hydrostatic sigma-p coordinate primitive equation set, with x,y-space differencing accurate to eighth order. A three-step dynamic initialization procedure is employed between the analysis of real-time data and grid interpolation. Results of an 18-hour simulation during which synoptic scale cyclogenesis, subsynoptic scale jet streak adjustments, mesoscale convergence zones and tornadic storms were observed have shown the present model to have the potential for simulating the fine-scale structure of features associated with cyclogenesis and intense squall-line development. The mesoscale model was also found to produce less truncation than the NWS LFM model, although a frictionless version of the mesoscale model somewhat overdeepens and overaccelerates features.

Kaplan, M. L.

Interactive applications of satellite observations and mesoscale numerical models

The impact of numerical weather prediction (NWP) and satellite meteorology on operational weather forecasting has become overwhelming in the past few years. The paper looks toward the merger of these technologies in making short range 6-18 h forecasts through the use of mesoscale NWP models. A short-range (2-18 h) mesoscale forecast system envisioned for the near future is described that includes four components: hydrodynamic numerical models, large-scale and mesoscale; satellites, polar orbiting for high latitudes and geostationary for low latitudes; mesoclimatology, derived in large part from satellite data; and special-purpose simple models and empirical relations. It is important that the components of the forecast system be developed in parallel rather than in series if the system is to be completed within five years. There is enough evidence to substantiate the revolution in the mesoscale weather prediction in the next five years.

Kreitzberg, C. W.

Supernova remnant revolution in an inhomogeneous medium. I - Numerical models

The first numerical simulations of supernova remnant evolution in an inhomogeneous gas are presented. Evolution in the lowest density substrate (the intercloud) is assumed to be spherically symmetric with a large intercloud filling factor and many dense regions (clouds) within the remnant; however, mass momentum and energy transfer between cloud and intercloud are included and the position and morphology of individual clouds tracked. Evolution is considered in several different models of the interstellar medium, both those in which the intercloud gas is diffuse (0.001 to 0.01/cu cm) and those in which it is relatively dense (n approximately 0.3/cu cm) under a variety of assumptions about the efficiency of thermal evaporation from the clouds into the intercloud medium.

Cowie, L. L.

Numerical modeling of quasi-static coronal loops. I - Uniform energy input

A quasi-static numerical model for coronal loops is considered for the case of a uniform energy input per unit volume into the loops. A line dipole model is used to represent the loop magnetic field, and the variations in loop cross section observed in X-ray photographs are parameterized by the ratio between the cross-sectional areas at the loop apex and base. The results of numerical modeling indicate that for an area ratio greater than unity, increases in the area ratio of a loop with a given length and apex area cause a general rise in electron density and a fall in the temperature gradient, leading to large increases in the differential emission factor at high temperatures. The differential function obtained is significantly different from that predicted by analytical models; however, analytical predictions for the temperature-electron density relations are comparable to numerical results. It is also concluded that even a symmetrical loop may have a maximum temperature away from the apex.

Vesecky, J. F.

Numerical model of the moon's rotation

The differential equations for the Euler angles that define the moon's orientation with respect to an inertial coordinate frame have been integrated numerically. This numerical model of the moon's rotation, when used in conjunction with current lunar-orbit and earth-rotation models, fits lunar laser ranging observations over a five-year period within 28 cm (rms). The present model is also compared with the numerical model of Williams (1975) and Eckhardt's (1981) semianalytical model.

Cappalo, R. J.

Numerical modeling of two-dimensional confined flows

A numerical model of two-dimensional confined flows is presented. The flow in the duct is partitioned into finite streams. The difference equations are then obtained by applying conservation principles directly to the individual streams. A listing of a computer code based on this approach in FORTRAN 4 language is presented. The code computes two dimensional compressible turbulent flows in ducts when the duct area along the flow is specified and the pressure gradient is unknown.

Greywall, M. S.

A numerical model for low-latitude ionospheric TEC

A numerical model of total electron content (TEC) at low latitudes has been constructed from the Faraday rotation data recorded at a number of locations in India during the period Oct. 1975-July 1976 when the geostationary satellite ATS-6 was located at 35 deg E. In all, 36 coefficients for each season are required to represent the model. The model can be used for satellite tracking systems in estimating quickly the range, range rate and angular refraction errors. The present model is an improvement over the previously constructed model based on orbiting satellite data.

Sethia, G.

A numerical model study of long-term planetary wave predictability

A numerical model was developed to study planetary wave behavior over periods up to 4 months under the influences of stochastic and diffusive forcing by synoptic scale motions and of sea surface temperature anomalies. The model is suitable for Monte Carlo experiments which provide ensemble mean and variance predictions and permit the response to heating anomalies to be partly separated from random variability. The model incorporates the linear-balanced dynamical equations for wind at two tropospheric levels and temperature at one level. The model atmosphere is forced by semiempirical January heating linearly dependent on local air temperature. Topography and the lower boundary layer are represented by vertical motion and drag functions, respectively, of the lower level wind. The linear phase speeds of Rossby waves are modified to reflect the effect of the stratosphere.

Wobus, R. L.

Verification and transfer of thermal pollution model. Volume 5: Verification of 2-dimensional numerical model

The six-volume report: describes the theory of a three dimensional (3-D) mathematical thermal discharge model and a related one dimensional (1-D) model, includes model verification at two sites, and provides a separate user's manual for each model. The 3-D model has two forms: free surface and rigid lid. The former, verified at Anclote Anchorate (FL), allows a free air/water interface and is suited for significant surface wave heights compared to mean water depth; e.g., estuaries and coastal regions. The latter, verified at Lake Keowee (SC), is suited for small surface wave heights compared to depth (e.g., natural or man-made inland lakes) because surface elevation has been removed as a parameter. These models allow computation of time dependent velocity and temperature fields for given initial conditions and time-varying boundary conditions.

Lee, S. S.

Numerical Modeling of the Transient Chilldown Process of a Cryogenic Propellant Transfer Line

Before cryogenic fuel depots can be fully realized, efficient methods with which to chill down the spacecraft transfer line and receiver tank are required. This paper presents numerical modeling of the chilldown of a liquid hydrogen tank-to-tank propellant transfer line using the Generalized Fluid System Simulation Program (GFSSP). To compare with data from recently concluded turbulent LH2 chill down experiments, seven different cases were run across a range of inlet liquid temperatures and mass flow rates. Both trickle and pulse chill down methods were simulated. The GFSSP model qualitatively matches external skin mounted temperature readings, but large differences are shown between measured and predicted internal stream temperatures. Discrepancies are attributed to the simplified model correlation used to compute two-phase flow boiling heat transfer. Flow visualization from testing shows that the initial bottoming out of skin mounted sensors corresponds to annular flow, but that considerable time is required for the stream sensor to achieve steady state as the system moves through annular, churn, and bubbly flow. The GFSSP model does adequately well in tracking trends in the data but further work is needed to refine the two-phase flow modeling to better match observed test data.

Temperature Sensors

A mesoscale numerical model and the development of a severe storm prediction system

The use of a mesoscale numerical model for predicting preferred zones of severe storm development is analyzed. A 60 consecutive day real-time test of the prediction system during the spring of 1978 proved useful in determining the problems and potentialities of such a system. A case study of severe storm development from this test period is described and compared to the model forecast fields.

Kaplan, M. L.

A survey of numerical models for wind prediction

A literature review is presented of the work done in the numerical modeling of wind flows. Pertinent computational techniques are described, as well as the necessary assumptions used to simplify the governing equations. A steady state model is outlined, based on the data obtained at the Deep Space Communications complex at Goldstone, California.

Schonfeld, D.

Low-level circulations induced by an upper-level jet in a fine mesh numerical model

Attention is given to varying the initial conditions in the Penn State three-dimensional numerical model in the development of the low-level jet over the Ohio Valley during the 10 May 1973 severe weather event. The results of two forecasts are summarized: (1) the control forecast which interpolated the National Meteorological Center's operational analysis to the mesoscale model grid; and (2) modifying the control forecast by inserting artificial soundings and surface observations. Enhancements are made on the basis of subsynoptic scale analyses by Uccellini and surface data. The results indicate that the low-level jet was in part induced by pressure changes associated with deceleration and horizontal divergence of air in the exit region of the upper-level jet. The enhancement of the initial moisture field using surface reports of clouds and precipitation has a substantial effect on precipitation.

Burkhart, R. P.

Project Fog Drops 5. Task 1: A numerical model of advection fog. Task 2: Recommendations for simplified individual zero-gravity cloud physics experiments

A two-dimensional numerical model was used to investigate the formation of marine advection fog. The model predicts the evolution of potential temperature, horizontal wind, water vapor content, and liquid water content in a vertical cross section of the atmosphere as determined by vertical turbulent transfer and horizontal advection, as well as radiative cooling and drop sedimentation. The model is designed to simulate the formation, development, or dissipation of advection fog in response to transfer of heat and moisture between the atmosphere and the surface as driven by advection over horizontal discontinuities in the surface temperature. Results from numerical simulations of advection fog formation are discussed with reference to observations of marine fog. A survey of candidate fog or cloud microphysics experiments which might be performed in the low gravity environment of a shuttle-type spacecraft in presented. Recommendations are given for relatively simple experiments which are relevent to fog modification problems.

Rogers, C. W.

The multiple nesting of mesoscale, submesoscale, and nonhydrostatic microscale numerical models for a case of mesocyclogenesis and severe storm development

A fundamental approach involving the integration of two different numerical models over several different spatial mesh lengths is considered. Model I is a hydrostatic, incompressible, and adiabatic primitive (Newtonian) equation set with no boundary layer forcing, i.e., surface heat or moisture flux or external frictional stresses. Typically this model is initialized from conventional radiosonde data and run for a 12-15 hour forecast period over a 40 km horizontal mesh. If phase I of the tornado developmental processes develops, the data is 'frozen' in space and time, interpolated to a 20 km mesh, and integrated for a shorter time period (3 hours) with a shorter time step. Model II is a nonhydrostatic compressible 18-level complement of model I. If model II indicates the development of the intensification of phase V dynamical processes, the data is frozen and interpolated to finer mesh lengths until the tornadic circulation system is approximated

Kaplan, M. L.

Random element method for numerical modeling of diffusional processes

The random element method is a generalization of the random vortex method that was developed for the numerical modeling of momentum transport processes as expressed in terms of the Navier-Stokes equations. The method is based on the concept that random walk, as exemplified by Brownian motion, is the stochastic manifestation of diffusional processes. The algorithm based on this method is grid-free and does not require the diffusion equation to be discritized over a mesh, it is thus devoid of numerical diffusion associated with finite difference methods. Moreover, the algorithm is self-adaptive in space and explicit in time, resulting in an improved numerical resolution of gradients as well as a simple and efficient computational procedure. The method is applied here to an assortment of problems of diffusion of momentum and energy in one-dimension as well as heat conduction in two-dimensions in order to assess its validity and accuracy. The numerical solutions obtained are found to be in good agreement with exact solution except for a statistical error introduced by using a finite number of elements, the error can be reduced by increasing the number of elements or by using ensemble averaging over a number of solutions.

Ghoniem, A. F.