Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical model”

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

Dynamic responses of submerged components in advanced reactors: experimental and numerical studies

The seismic design of an advanced nuclear reactor must consider the interaction of vessel internal components with the surrounding coolant: fluid–structure interaction (FSI). Available analytical solutions for FSI of submerged components do not accommodate multiple-component, intense seismic inputs and complex reactor and internal geometries. Physical testing of reactor vessels and internals to inform seismic design is impractical and cost-prohibitive, leaving the use of verified and validated, robust numerical models as the only plausible option for analysis and design. Physical data that could be used for validating such numerical models for multi-component shaking are not available. This article describes an experimental program performed on a 6-degree-of-freedom earthquake simulator to generate data that could support validation of seismic FSI numerical models for submerged components in commercial finite element codes. A scaled model of a base-supported reactor vessel with simplified representations of submerged internals was tested to generate submerged-component response histories for a range of seismic inputs. The generated data were used to validate numerical models in the finite element code LS-DYNA. Numerical models were validated for calculation of hydrodynamic pressure responses on internals, in-water frequencies of internals, and acceleration responses of internals. The generated data and the analysis recommendations could aid engineering analysts designing submerged components and systems for seismic effects.

Engineering↗

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

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↗