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 631 records · Page 35

A numerical circulation model with topography for the Martian Southern Hemisphere

A quasi-geostrophic numerical model, including friction, radiation, and the observed planetary topography, is applied to the general circulation of the Martian atmosphere in the Southern Hemisphere at latitudes south of about 35 deg. Near equilibrium weather systems developed after about 5 model days. To avoid violating the quasi-geostrophic approximation, only 0.8 of the already smoothed relief was employed. Weather systems and velocity fields are strikingly tied to topography. A 2mb middle latitude jet stream is found of remarkably terrestrial aspect. Highest surface velocities, both horizontal and vertical, are predicted in western Hellas Planitia and eastern Argyre Planitia, which are observed to be preferred sites of origin of major Martian dust storms. Mean horizontal velocities and vertical velocities are found just above the surface velocity boundary layer.

Mass, C.↗

A numerical circulation model with topography for the Martian southern hemisphere

A two-level numerical model taking account of the observed planetary topography is applied to the general circulation of the Martian atmosphere in the southern hemisphere. The physics of this model is based on the assumption of a hydrostatic quasi-geostrophic ideal-gas atmosphere composed entirely of carbon dioxide; radiative transfer among the atmosphere, space, and the ground is included along with surface friction. The calculations indicate the development of near-equilibrium weather systems after about five model days, weather systems and velocity fields whose positions appear to be connected with topography, and a midlatitude jet stream strikingly similar to the terrestrial case. The results also suggest that the vertical and horizontal wind components in western Hellas and eastern Argyre produce copious amounts of suspendible dust each year, that both major and minor dust storms should be a common occurrence on Mars, and that the entire south circumpolar zone appears to be a dust sink in winter.

Mass, C.↗

Equivalent beam modeling using numerical reduction techniques

Numerical procedures that can accomplish model reductions for space trusses were developed. Three techniques are presented that can be implemented using current capabilities within NASTRAN. The proposed techniques accomplish their model reductions numerically through use of NASTRAN structural analyses and as such are termed numerical in contrast to the previously developed analytical techniques. Numerical procedures are developed that permit reductions of large truss models containing full modeling detail of the truss and its joints. Three techniques are presented that accomplish these model reductions with various levels of structural accuracy. These numerical techniques are designated as equivalent beam, truss element reduction, and post-assembly reduction methods. These techniques are discussed in detail.

Chapman, J. M.↗

A numerical study of process complexity in permafrost dominated regions

Numerical modeling of permafrost dynamics requires adequate representation of atmospheric and surface processes, a reasonable parameter estimation strategy, and site-specific model development. The three main research objectives of the study are: (i) to propose a novel methodology that determines the required level of surface process complexity of permafrost models by conducting parameter sensitivity and calibration, (ii) to design and compare three numerical models of increasing surface process complexity, and (iii) to calibrate and validate the numerical models at the Yakou catchment on the Qinghai-Tibet Plateau as an exemplary study site. The calibration was carried out by coupling the Advanced Terrestrial Simulator (numerical model) and PEST (calibration tool). Simulation results showed that (i) A simple numerical model that considers only subsurface processes can simulate active layer development with the same accuracy as other more complex models that include surface processes. (ii) Peat and mineral soil layer permeability, Van Genuchten alpha, and porosity are highly sensitive. (iii) Liquid precipitation aids in increasing the rate of permafrost degradation. (iv) Deposition of snow insulated the subsurface during the thaw initiation period. We have developed and released an integrated code that couples the numerical software ATS to the calibration software PEST. The numerical model can be further used to determine the impacts of climate change on permafrost degradation.

Calibration↗

Middle Atmosphere Dynamics with Gravity Wave Interactions in the Numerical Spectral Model: Tides and Planetary Waves

As Lindzen (1981) had shown, small-scale gravity waves (GW) produce the observed reversals of the zonal-mean circulation and temperature variations in the upper mesosphere. The waves also play a major role in modulating and amplifying the diurnal tides (DT) (e.g., Waltersheid, 1981; Fritts and Vincent, 1987; Fritts, 1995a). We summarize here the modeling studies with the mechanistic numerical spectral model (NSM) with Doppler spread parameterization for GW (Hines, 1997a, b), which describes in the middle atmosphere: (a) migrating and non-migrating DT, (b) planetary waves (PW), and (c) global-scale inertio gravity waves. Numerical experiments are discussed that illuminate the influence of GW filtering and nonlinear interactions between DT, PW, and zonal mean variations. Keywords: Theoretical modeling, Middle atmosphere dynamics, Gravity wave interactions, Migrating and non-migrating tides, Planetary waves, Global-scale inertio gravity waves.

Mayr, Hans G.↗

Impact of turbulence modeling on numerical accuracy and efficiency of compressible flow simulations

Discussed is the numerical implementation of turbulence models used in viscous compressible flow simulations and their performance described with respect to numerical accuracy, efficiency and stability. The basic approach utilizes the Reynolds-averaged compressible Navier-Stokes equations in which the Reynolds stresses and heat fluxes are mathematically modeled by suitable turbulence models. The turbulence models investigated include zero-, one-, and two-equation eddy viscosity models. The flow fields investigated include theoretical and supersonic flows about two-dimensional and axisymmetric bodies. Discussions concerning the numerical implementation of models include differencing procedures and boundary conditions used to assume numerical stability and accuracy. Numerical performance is also evaluated by comparing computations with experimental results.

Coakley, Thomas J.↗

Thermal Analysis and Design of Multi-layer Insulation for Re-entry Aerodynamic Heating

The combined radiation/conduction heat transfer in high-temperature multi-layer insulations was modeled using a finite volume numerical model. The numerical model was validated by comparison with steady-state effective thermal conductivity measurements, and by transient thermal tests simulating re-entry aerodynamic heating conditions. A design of experiments technique was used to investigate optimum design of multi-layer insulations for re-entry aerodynamic heating. It was found that use of 2 mm foil spacing and locating the foils near the hot boundary with the top foil 2 mm away from the hot boundary resulted in the most effective insulation design. A 76.2 mm thick multi-layer insulation using 1, 4, or 16 foils resulted in 2.9, 7.2, or 22.2 percent mass per unit area savings compared to a fibrous insulation sample at the same thickness, respectively.

Daryabeigi, Kamran↗

A numerical cloud model for the support of laboratory experimentation

A numerical cloud model is presented which can describe the evolution of a cloud starting from moist aerosol-laden air through the diffusional growth regime. The model is designed for the direct support of cloud chamber laboratory experimentation, i.e., experiment preparation, real-time control and data analysis. In the model the thermodynamics is uncoupled from the droplet growth processes. Analytic solutions for the cloud droplet growth equations are developed which can be applied in most laboratory situations. The model is applied to a variety of representative experiments.

Hagen, D. E.↗

Modeling of Chill Down in Cryogenic Transfer Lines

A numerical model to predict chill down in cryogenic transfer lines has been developed. Three chill down cases using hydrogen as the working fluid are solved: 1) a simplified model amenable to analytical solution, 2) a realistic model of superheated vapor flow, and 3) a realistic model of initially subcooled liquid flow. The first case compares a numerical model with an analytical solution with very good agreement between the two. Additionally, the analytical solution provides a convenient way to look at parametric effects on the chill down. The second and third cases are numerical models which provide temperature histories of the fluid and solid tube wall during chill down as well as several other quantities of interest such as pressure and mass flow rate. Of great interest is the ability to predict accurate values of chill down time (the time required to achieve steady-state cryogenic flow). The models predict that a 26 in. long, 3/16 in. ID aluminum tube has a shorter chill down time (approx. equal to 100 sec) and uses less hydrogen with superheated vapor flow than with initially subcooled liquid flow (greater than 200 sec for chill down).

Cross, Matthew F.↗

Tests of a procedure for inserting satellite radiance measurements into a numerical circulation model.

Description of experiments designed to test the sensitivity of a new method for retrieving vertical temperature profiles from satellite radiance measurements to the initial guess of temperature distribution. Emphasis is placed on initial guesses of the kind that would be provided by a numerical model of the atmospheric circulation, with the objective of developing a procedure for direct use of satellite radiation measurements as input to a numerical forecasting scheme.

Hogan, J. S.↗

Computational general relativistic force-free electrodynamics: II. Characterization of numerical diffusivity

Scientific codes are an indispensable link between theory and experiment; in (astro-)plasma physics, such numerical tools are one window into the universe’s most extreme flows of energy. The discretization of Maxwell’s equations – needed to make highly magnetized (astro)physical plasma amenable to its numerical modeling – introduces numerical diffusion. It acts as a source of dissipation independent of the system’s physical constituents. Understanding the numerical diffusion of scientific codes is the key to classifying their reliability. It gives specific limits in which the results of numerical experiments are physical. We aim at quantifying and characterizing the numerical diffusion properties of our recently developed numerical tool for the simulation of general relativistic force-free electrodynamics by calibrating and comparing it with other strategies found in the literature. Our code correctly models smooth waves of highly magnetized plasma. We evaluate the limits of general relativistic force-free electrodynamics in the context of current sheets and tearing mode instabilities. We identify that the current parallel to the magnetic field (j ∥ ), in combination with the breakdown of general relativistic force-free electrodynamics across current sheets, impairs the physical modeling of resistive instabilities. We find that at least eight numerical cells per characteristic size of interest (e.g., the wavelength in plasma waves or the transverse width of a current sheet) are needed to find consistency between resistivity of numerical and of physical origins. High-order discretization of the force-free current allows us to provide almost ideal orders of convergence for (smooth) plasma wave dynamics. The physical modeling of resistive layers requires suitable current prescriptions or a sub-grid modeling for the evolution of j ∥ .

79 ASTRONOMY AND ASTROPHYSICS↗

Estimation of surface heat and moisture fluxes over a prairie grassland. IV - Impact of satellite remote sensing of slow canopy variables on performance of a hybrid biosphere model

Numerical experiments are conducted using the Ex-BATS model of Crosson and Weng (1993), which is an adaptation the Dickinson (1983, 1984) and Dickinson et al. (1986) biosphere model BATS. The purpose of these experiments is the assessment of the Ex-BATS performance when using remotely sensed data for the estimation of three key canopy variables retrieved from NOAA-AVHRR measurements: the total surface albedo, the leaf area index (LAI), and the nondiurnally varying component of stomatal resistance, r(s). The results of the simulations, which cover the entire FIFE 1987 time period, show that the satellite retrievals of r(s) are only 20 to 30 percent less accurate than the idealized results of the control experiment. The performance of the model which used satellite retrieval of the surface albedo and LAI was essentially equivalent to the hypothetical version.

Crosson, William L.↗

Wave Energy Converter Power Take-Off Modeling and Validation From Experimental Bench Tests

This article describes the implementation of a new numerical model of the power take-off system installed in the Monterey Bay Aquarium Research Institute wave energy converter, a device developed to provide power to various oceanic research missions. The simultaneous presence of hydraulic, pneumatic, and electrical subsystems in the power take-off system represents a significant challenge in forging an accurate model able to replicate the main dynamic characteristics of the system. The validation of the new numerical model is addressed by comparing simulations with the measurements obtained during a series of bench tests. Data from the bench tests show good agreement with the numerical model. The validated model provides deeper insights into the complex nonlinear dynamics of the power take-off system and will support further performance improvements in the future.

16 TIDAL AND WAVE POWER↗

Prospect of Using Numerical Dynamo Model for Prediction of Geomagnetic Secular Variation

Modeling of the Earth's core has reached a level of maturity to where the incorporation of observations into the simulations through data assimilation has become feasible. Data assimilation is a method by which observations of a system are combined with a model output (or forecast) to obtain a best guess of the state of the system, called the analysis. The analysis is then used as an initial condition for the next forecast. By doing assimilation, not only we shall be able to predict partially secular variation of the core field, we could also use observations to further our understanding of dynamical states in the Earth's core. One of the first steps in the development of an assimilation system is a comparison between the observations and the model solution. The highly turbulent nature of core dynamics, along with the absence of any regular external forcing and constraint (which occurs in atmospheric dynamics, for example) means that short time comparisons (approx. 1000 years) cannot be made between model and observations. In order to make sensible comparisons, a direct insertion assimilation method has been implemented. In this approach, magnetic field observations at the Earth's surface have been substituted into the numerical model, such that the ratio of the multiple components and the dipole component from observation is adjusted at the core-mantle boundary and extended to the interior of the core, while the total magnetic energy remains unchanged. This adjusted magnetic field is then used as the initial field for a new simulation. In this way, a time tugged simulation is created which can then be compared directly with observations. We present numerical solutions with and without data insertion and discuss their implications for the development of a more rigorous assimilation system.

Kuang, Weijia↗

Laser-Processed Condensing Heat Exchanger for Space Applications

In sizing a condensing heat exchanger for space life support systems, there is a limit to allowable pressure drop along with a required vapor recuperation rate. This sizing is performed by using known heat transfer and pressure drop correlations in a numerical model. Dimpling the surface of an air-liquid channeled heat exchanger increases the convective heat transfer rate however creates difficulties within the numerical model as the flat plate correlations are no longer accurate. To increase the fidelity of the numerical model used for sizing the heat exchanger Computational Fluid Dynamics (CFD) with an SST k-ω turbulence model was performed on a periodic portion of the heat exchanger. Through CFD modeling a heat flux and pressure drop was measured across the heat exchanger channel, this information was used to create a localized Nusselt number that was incorporated into the numerical model along with a more accurate pressure drop due to the dimpled surface. Additionally, through modeling the heat exchanger after manufacturing and assembly in Computer Aided Drafting (CAD) software there was an increase in the measurement accuracy of geometric features which was incorporated into the numerical model for further increase in fidelity.

heat exchanger↗

Meso-beta scale numerical simulation studies of terrain-induced jet streak mass and momentum perturbations

An in-depth analysis of observed gravity waves and their relationship to precipitation bands over the Montana mesonetwork during the 11-12 July 1981 CCOPE case study indicated two episodes of coherent waves. While geostrophic adjustment, shearing instability, and terrain were all implicated separately or in combination as possible wave generation mechanisms, the lack of upper-air data within the wave genesis region made it difficult to define the genesis processes from observations alone. The first part of this paper, 3D Numerical Modeling Studies of Terrain-Induced Mass/Momentum Perturbations, employs a mesoscale numerical model to help diagnose the intricate early wave generation mechanisms during the first observed gravity wave episode. The meso-beta scale numerical model is used to study various simulations of the role of multiple geostrophic adjustment processes in focusing a region for gravity wave genesis. The second part of this paper, Linear Theory and Theoretical Modeling, investigates the response of non-resting rotating homogeneous and continuously stratified Boussinesq models of the terrestrial atmosphere to temporally impulsive and uniformly propagating three-dimensional localized zonal momentum sources representative of midlatitude jet streaks. The methods of linear perturbation theory applied to the potential vorticity (PV) and wave field equations are used to study the geostrophic adjustment dynamics. The total zonal and meridional wind perturbations are separated into geostrophic and ageostrophic components in order to define and follow the evolution of both the primary and secondary mesocirculations accompanying midlatitude jetogenesis forced by geostrophic adjustment processes. This problem is addressed to help fill the gap in understanding the dynamics and structure of mesoscale inertia-gravity waves forced by geostrophic adjustment processes in simple two-dimensional quiescent current systems and those produced by mesoscale numerical models simulating the orographic and diabatic perturbation of three-dimensional quasi-geostrophically balanced synoptic scale jet streaks associated with complex baroclinic severe storm producing environments.

Lin, Yuh-Lang↗

A Numerical Turbulence Model for Multiphase Flows in the Protoplanetary Nebula

It is thought that planets form from solid particles in a flattened, rotating, 99% gaseous nebula. These grains gradually coagulate into millimeter-to-meter sized aggregates which settle toward the midplane of the nebula. It is widely believed that the resulting dense layer eventually becomes gravitationally unstable and collapses into 'planetesimals.' A new numerical model is presented to simulate the predominant processes (gravitation, vertical convection, and shear-driven turbulence) during the stage while the particulate material is still dispersed about the midplane of the nebula. In our previous work, particles were assumed to be spheres of a single radius; in the present work, particles are spheres of different radii. Results indicate that neither a broad nor a narrow distribution of particle sizes is likely to become gravitationally unstable.

Champney, Joelle M.↗

Investigation of the Seasonal Freeze/Thaw Cycle of Soils in the GAPP Regions

The objective of this project is to understand the current and future impacts and feedbacks of seasonally frozen ground on the regional and global climate system and the ecosystem as a whole. More specifically, we proposed to: investigate the dynamics of the seasonally freezing and thawing processes and their relations to current climatic conditions through in situ data analysis and numerical modeling; develop a comprehensive algorithm to detect surface soil freeze/thaw status for snow-free land surface using passive microwave remote sensing data and for snow-covered land surface using numerical modeling, and validate the algorithm and numerical model using available ground-based measurements; investigate seasonal and inter-annual variations of frozen soils in the GAPP regions using data generated from the new validated algorithm and their relations to the environmental factors such as air temperature, snow cover, surface morphology, soil type and soil moisture; generate a frozen soil data set which describes the timing, duration, thickness, and area extent of seasonally frozen ground for the period from 1978 to present. The proposed research will be accomplished through data analysis, remote sensing, and numerical modeling.

FREEZE-THAW TESTS↗