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 271 records · Page 15

Middle atmosphere project. A semi-spectral numerical model for the large-scale stratospheric circulation

The complete model is a semispectral model in which the longitudinal dependence is represented by expansion in zonal harmonics while the latitude and height dependencies are represented by a finite difference grid. The model is based on the primitive equations in the log pressure coordinate system. The lower boundary of the model domain is set at the 100 mb level (i.e., near the tropopause) and the effects of tropospheric forcing are included in the lower boundary condition. The upper boundary is at approximately 96 km, and the latitudinal extent is either global or hemispheric. The basic differential equations and boundary conditions are outlined. The finite difference equations are described. The initial conditions are discussed and a sample calculation is presented. The FORTRAN code is given in the appendix.

Holton, J. R.↗

A large-scale numerical model of sea ice

The described large-scale sea ice model, which is capable of coupling with atmospheric and oceanic models of comparable resolution, simulates the yearly cycle of ice in both the Northern and Southern Hemispheres. Model results for the yearly cycle of sea ice thickness and extent in both the Arctic and the Antarctic are presented. Horizontally the model resolution is approximately 200 km, while vertically four layers - ice, snow, ocean, and atmosphere - are considered. Thermodynamic processes based on energy balances at the various interfaces and dynamic processes based on wind stress, water stress, Coriolis force, internal ice resistance, and the stress from the tilt of the sea surface are incorporated. It is assumed that the ice within a given grid square is of uniform thickness, although each square has a variable percentage of its area ice free.

Parkinson, C. L.↗

An isentropic and sigma coordinate hybrid numerical model - Model development and some initial tests

A solution is presented for matching boundary conditions across the interface of an isentropic and sigma coordinate hybrid model. A hybrid model based on the flux form of the primitive equations is developed which allows direct vertical exchange between the model domains, satisfies conservation principles with respect to transport processes, and maintains a smooth transition across the interface without need for artificial adjustment or parameterization schemes. The initial hybrid model simulations of a jet streak propagating in a zonal channel are used to test the feasibility of the hybrid model approach. High efficiency of the hybrid model is demonstrated.

Uccellini, L. W.↗

Numerical modeling of baroclinic instability

The Eady model of baroclinic instability applied to a rotating channel with Ekman layers at the top and bottom is extended to a beta-plane channel. Changes in the shape and position of the boundary separating the regions of axially symmetric and wavelike motion are predicted.

Geisler, J. F.↗

The structure and appearance of winds from supercritical accretion disks. I - Numerical models

Equations for the structure and appearance of supercritical accretion disks and the radiation-driven winds which emanate from them are derived and solved by a steady-state hydrodynamic computer code with a relaxation technique used in stellar structure problems. The present model takes into account the mass of the accreting star, the total accretion rate, a generalization of the disk alpha parameter which accounts for heating by processes in addition to viscosity, and the ratio of the total luminosity to the Eddington luminosity. Solutions indicate that for accretion onto a hard-surfaced star, steady, optically thick winds result for even slightly supercritical accretion, and the object will appear as a supergiant star with a high mass loss rate and a nonblackbody spectrum. Winds from black hole accretion disks are expected to depend on the form of the accretion interior to the critical radius, possibly consisting of no ejection at all, a wind similar to that of a hard-surfaced star, or a column of material ejected from a hole in the accretion disk.

Meier, D. L.↗

Middle atmosphere project: A radiative heating and cooling algorithm for a numerical model of the large scale stratospheric circulation

A Curtis matrix is used to compute cooling by the 15 micron and 10 micron bands of carbon dioxide. Escape of radiation to space and exchange the lower boundary are used for the 9.6 micron band of ozone. Voigt line shape, vibrational relaxation, line overlap, and the temperature dependence of line strength distributions and transmission functions are incorporated into the Curtis matrices. The distributions of the atmospheric constituents included in the algorithm, and the method used to compute the Curtis matrices are discussed as well as cooling or heating by the 9.6 micron band of ozone. The FORTRAN programs and subroutines that were developed are described and listed.

Wehrbein, W. M.↗

Compensation for use of monthly-averaged winds in numerical modeling

Ratios R of the monthly averaged wind speeds to the magnitudes of the monthly averaged wind vectors are presented over a 41 x 41 grid covering the southern Ocean and the Antarctic continent. The ratio is found to vary from 1 to over 1000, with an average value of 1.86. These ratios R are relevant for converting from sensible and latent heats calculated with mean monthly data to those calculated with 12 hourly data. The corresponding ratios alpha for wind stress, along with the angle deviations involved, are also presented over the same 41 x 41 grid. The values of alpha generally exceed those for R and average 2.66. Regions in zones of variable wind directions have larger R and alpha ratios, over the ice-covered portions of the southern Ocean averaging 2.74 and 4.35 for R and alpha respectively. Thus adjustments to compensate for the use of mean monthly wind velocities should be stronger for wind stress than for turbulent heats and stronger over ice covered regions than over regions with more persistent wind directions, e.g., those in the belt of mid-latitude westerlies.

Parkinson, C. L.↗

Experimental validation of a numerical model predicting the charging characteristics of Teflon and Kapton under electron beam irradiation

The effect differential charging of spacecraft thermal control surfaces is assessed by studying the dynamics of the charging process. A program to experimentally validate a computer model of the charging process was established. Time resolved measurements of the surface potential were obtained for samples of Kapton and Teflon irradiated with a monoenergetic electron beam. Results indicate that the computer model and experimental measurements agree well and that for Teflon, secondary emission is the governing factor. Experimental data indicate that bulk conductivities play a significant role in the charging of Kapton.

Hazelton, R. C.↗

Interactive initialization of heat flux parameters for numerical models using satellite temperature measurements

Day/night Heat Capacity Mapping Mission image pairs over Kansas and Indiana were used to examine the spatial variation of moisture availability on the mesoscale. Of particular concern was whether patterns of moisture availability and thermal inertia, as well as surface heat fluxes, respond to significant spatial variations in the rainfall pattern. A temperature analysis for the Indiana case is presented.

Carlson, T. N.↗

Numerical models for Halley dust environments

The algebra, parameter values, and results of a reference model for the dust environment of Comet P/Halley are presented. The model employs simplified treatments of the icy conglomerate nucleus and of the hydrodynamic outflow of gas and dust. The model includes a gas production rate derived from observations of Halley's visual magnitude in 1910, gas and dust flux nonuniformly distributed over the nuclear surface, a dust distribution function consistent with large particle inferences from comet antitails, dust density variable with mass, and an approximation to the distribution of particles near the dust envelope. Many parameter values are derived by analogy with the results for other comets. The features are incorporated in a computer program which calculates the concentration, flux, and fluence of dust particles as functions of time, position and particle mass.

Divine, N.↗

Numerical Modeling of Three-Dimensional Confined Flows

A three dimensional confined flow model is presented. The flow field is computed by calculating velocity and enthalpy along a set of streamlines. The finite difference equations are obtained by applying conservation principles to streamtubes constructed around the chosen streamlines. With appropriate substitutions for the body force terms, the approach computes three dimensional magnetohydrodynamic channel flows. A listing of a computer code, based on this approach is presented in FORTRAN IV language. The code computes three dimensional compressible viscous flow through a rectangular duct, with the duct cross section specified along the axis.

Greywall, M. S.↗

Interactive initialization of heat flux parameters for numerical models using satellite temperature measurements

Efforts were made (1) to bring the image processing and boundary layer model operation into a completely interactive mode and (2) to test a method for determining the surface energy budget and surface moisture availability and thermal inertia on a scale appreciably larger than that of the city. A region a few hundred kilometers on a side centered over southern Indiana was examined.

Carlson, T. N.↗

Numerical modeling and optimization of polymer melt processing operations

The application of finite element computer analyses to polymer flows of the type encountered in melt processing operations is described. A code capable of predicting values of fluid veleocity, pressure, shear stress, and temperature at any point within the flow field was developed. As such, is is of value in diagnosing such processing problems as regions of fluid stagnation at which thermal degradation may occur, or regions of excessive shear deformation which lead to thermomechanical damage. It is further able to generate predictions of the forces which must be applied to the melt to achieve the desired flow, and this information is of value in designing processing equipment of optimal efficiency and minimum energy consumption.

Roylance, D.↗

Interactive initialization of heat flux parameters for numerical models using satellite temperature measurements

A method for obtaining patterns of moisture availability (and net evaporation) from satellite infrared measurements employs Carlson's boundary layer model and a variety of image processing routines executed by a minicomputer. To test the method with regard to regional scale moisture analyses, two case studies were chosen because of the availability of HCMM data and because of the presence of a large horizontal gradient in antecedent precipitation and crp moisture index. Results show some correlation in both cases between antecedent precipitation and derived moisture availability. Apparently, regional-scale moisture availability patterns can be determined with some degree of fidelity but the values themselves may be useful only in the relative sense and significant to within plus or minus one category of dryness over a range of 4 or 5 categories between absolutely dry and field saturation. Preliminary results suggest that the derived moisture values correlate best with longer-term precipitation totals, suggesting that the infrared temperatures respond more sensitively to a relatively deep substrate layer.

Carlson, T. N.↗