Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical models”

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 505 records · Page 28

A Numerical and Experimental Investigation of Flow Separation Control over a Wall-Mounted Hump Model

Numerical simulations and wind tunnel experiments were performed to investigate flow separation and its control over a wall-mounted hump model. Three-dimensional, unsteady fluid dynamic simulations were used to capture the relevant flow physics of unsteady flow separation and its control via active and passive flow control methods. Surface pressure measurements and surface oilflow visualization were obtained in the experiments. Current measurements documented a considerably thinner incoming boundary layer, but similar pressure distribution, compared to the reference data reported in the CFD (Computational Fluid Dynamics) Validation (CFDVAL 2004) Experiments. Consistent with the data available in the literature, numerical simulations predicted a longer separation bubble, which manifested itself as a shift in the pressure distributions. Streamwise vortices generated by passive vortex generators increased the suction pressure and provided substantial pressure recovery. Numerical simulations of vortex generators qualitatively predicted the effects of control, and revealed patterns of attached/separated flow regions. On the other hand, active flow control using steady discrete jets reduced, but was not able to eliminate, flow separation for the tested configuration. Although the major disagreement between the numerical and experimental results was in the pressure recovery region due to the inaccurate prediction of separation bubbles, the relative pressure recoveries from the baselines were comparable. Surface oilflow visualization and simulated surface streamlines were qualitatively in good agreement revealing the key flow features.

Koklu, Mehti↗

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

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

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↗

Finite-Strain Fractional-Order Viscoelastic (FOV) Material Models and Numerical Methods for Solving Them

Fraction-order viscoelastic (FOV) material models have been proposed and studied in 1D since the 1930's, and were extended into three dimensions in the 1970's under the assumption of infinitesimal straining. It was not until 1997 that Drozdov introduced the first finite-strain FOV constitutive equations. In our presentation, we shall continue in this tradition by extending the standard, FOV, fluid and solid, material models introduced in 1971 by Caputo and Mainardi into 3D constitutive formula applicable for finite-strain analyses. To achieve this, we generalize both the convected and co-rotational derivatives of tensor fields to fractional order. This is accomplished by defining them first as body tensor fields and then mapping them into space as objective Cartesian tensor fields. Constitutive equations are constructed using both variants for fractional rate, and their responses are contrasted in simple shear. After five years of research and development, we now possess a basic suite of numerical tools necessary to study finite-strain FOV constitutive equations and their iterative refinement into a mature collection of material models. Numerical methods still need to be developed for efficiently solving fraction al-order integrals, derivatives, and differential equations in a finite element setting where such constitutive formulae would need to be solved at each Gauss point in each element of a finite model, which can number into the millions in today's analysis.

Freed, Alan D.↗