Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical simulation”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

A numerical simulation of seasonal stratospheric climate. II - Energetics

The paper analyzes the energetics involved in the two numerical experiments described in Part I, concerning a longwave radiative transfer model and a simple Newtonian cooling (or heating) model. A three-dimensional quasi-geostrophic model is developed and used over the entire annual cycle to test the climatic responses of the stratospheric circulation to the cited two different longwave radiative heating models. It is shown that the eddy energy parameters (eddy kinetic energy, eddy available potential energy, generation of the latter, and vertical propagation of eddy geopotential energy) undergo annual variations with maximum values in winter and minimum values in summer. The zonal energy parameters exhibit semiannual variations with major maximum values in winter, minor maxima in summer and minima in spring and fall. Analysis of the energy cycle reveals that the vertical propagation of eddy geopotential energy and the generation of eddy available potential energy are the energy sources for energetics in the upper stratosphere. Ways in which the longwave radiative processes may affect the stratospheric energetics are identified.

Chen, T.-C.↗
NASA NTRS1987

Numerical simulation of three-dimensional unsteady vortex flow using a compact vorticity-velocity algorithm

A numerical algorithm is presented which is used to solve the unsteady, fully three-dimensional, incompressible Navier-Stokes equations in vorticity-velocity variables. A discussion of the discrete approximation scheme is presented as well as the solution method used to solve the resulting algebraic set of difference equations. Second order spatial and temporal accuracy is verified through solution comparisons with exact results obtained for steady three-dimensional stagnation point flow and unsteady axisymmetric vortex spin-up. In addition, results are presented for the problem of unsteady bubble-type vortex breakdown with emphasis on internal bubble dynamics and structure.

Gatski, T. B.↗
NASA NTRS1990

Numerical Simulation Of Buckling In Waffle Plants

Accurate results obtained when fillet radii considered. Two reports describe numerical and experimental study of application of PASCO and WAFFLE computer programs to analysis of buckling in integrally machined, biaxially stiffened panel. PASCO (Panal Analysis and Sizing Code) is finite-element stress-and-strain code written for analysis and sizing of uniaxially stiffened panels. WAFFLE program provides comprehensive stress analysis of waffle panel, used to determine bending moments at interfaces.

Yin, Dah N.↗
NASA NTRS1996

A Numerical Simulation of Scattering from One-Dimensional Inhomogeneous Dielectric Random Surfaces

In this paper, an efficient numerical solution for the scattering problem of inhomogeneous dielectric rough surfaces is presented. The inhomogeneous dielectric random surface represents a bare soil surface and is considered to be comprised of a large number of randomly positioned dielectric humps of different sizes, shapes, and dielectric constants above an impedance surface. Clods with nonuniform moisture content and rocks are modeled by inhomogeneous dielectric humps and the underlying smooth wet soil surface is modeled by an impedance surface. In this technique, an efficient numerical solution for the constituent dielectric humps over an impedance surface is obtained using Green's function derived by the exact image theory in conjunction with the method of moments. The scattered field from a sample of the rough surface is obtained by summing the scattered fields from all the individual humps of the surface coherently ignoring the effect of multiple scattering between the humps. The statistical behavior of the scattering coefficient sigma(sup 0) is obtained from the calculation of scattered fields of many different realizations of the surface. Numerical results are presented for several different roughnesses and dielectric constants of the random surfaces. The numerical technique is verified by comparing the numerical solution with the solution based on the small perturbation method and the physical optics model for homogeneous rough surfaces. This technique can be used to study the behavior of scattering coefficient and phase difference statistics of rough soil surfaces for which no analytical solution exists.

Sarabandi, Kamal↗
NASA NTRS1984

Numerical simulation of boundary-layer transition

The transition to turbulence in boundary layers was investigated by direct numerical solution of the nonlinear, three-dimensional, incompressible Navier-Stokes equations in the half-infinite domain over a flat plate. Periodicity was imposed in the streamwise and spanwise directions. A body force was applied to approximate the effect of a nonparallel mean flow. The numerical method was spectra, based on Fourier series and Jacobi polynomials, and used divergence-free basis functions. Extremely rapid convergence was obtained when solving the linear Orr-Sommerfeld equation. The early nonlinear and three-dimensional stages of transition, in a boundary layer disturbed by a vibrating ribbon, were successfully simulated. Excellent qualitative agreement was observed with either experiments or weakly nonlinear theories. In particular, the breakdown pattern was staggered or nonstaggered depending on the disturbance amplitude.

Spalart, P. R.↗
NASA NTRS2013

Combined Experimental and Numerical Simulations of Thermal Barrier Coated Turbine Blades Erosion

A combined experimental and computational study was conducted to investigate the erosion of thermal barrier coated (TBC) blade surfaces by alumina particles ingestion in a single stage turbine. In the experimental investigation, tests of particle surface interactions were performed in specially designed tunnels to determine the erosion rates and particle restitution characteristics under different impact conditions. The experimental results show that the erosion rates increase with increased impingement angle, impact velocity and temperature. In the computational simulations, an Euler-Lagrangian two stage approach is used in obtaining numerical solutions to the three-dimensional compressible Reynolds Averaged Navier-Stokes equations and the particles equations of motion in each blade passage reference frame. User defined functions (UDF) were developed to represent experimentally-based correlations for particle surface interaction models which were employed in the three-dimensional particle trajectory simulations to determine the particle rebound characteristics after each surface impact. The experimentally based erosion UDF model was used to predict the TBC erosion rates on the turbine blade surfaces based on the computed statistical data of the particles impact locations, velocities and angles relative to the blade surface. Computational results are presented for the predicted TBC blade erosion in a single stage commercial APU turbine, for a NASA designed automotive turbine, and for the NASA turbine scaled for modern rotorcraft operating conditions. The erosion patterns in the turbines are discussed for uniform particle ingestion and for particle ingestion concentrated in the inner and outer 5 percent of the stator blade span representing the flow cooling the combustor liner.

Hamed, Awate↗
NASA NTRS1992

The impact of dynamic data assimilation on the numerical simulations of the QE II cyclone and an analysis of the jet streak influencing the precyclogenetic environment

A mesoscale numerical model is combined with a dynamic data assimilation via Newtonian relaxation, or 'nudging', to provide initial conditions for subsequent simulations of the QE II cyclone. Both the nudging technique and the inclusion of supplementary data are shown to have a large positive impact on the simulation of the QE II cyclone during the initial phase of rapid cyclone development. Within the initial development period (from 1200 to 1800 UTC 9 September 1978), the dynamic assimilation of operational and bogus data yields a coherent two-layer divergence pattern that is not well defined in the model run using only the operational data and static initialization. Diagnostic analysis based on the simulations show that the initial development of the QE II storm between 0000 UTC 9 September and 0000 UTC 10 September was embedded within an indirect circulation of an intense 300-hPa jet streak, was related to baroclinic processes extending throughout a deep portion of the troposphere, and was associated with a classic two-layer mass-divergence profile expected for an extratropical cyclone.

Manobianco, John↗
NASA NTRS1991

Numerical simulation of a subtropical squall line over the Taiwan Strait

The present 2D nonhydrostatic, time-dependent numerical cloud model, which is applied to the structure and development of a subtropical squall line observed during the Taiwan Area Mesoscale Experiment, includes (1) a parameterized ice-phase microphysical scheme, (2) longwave and shortwave radiative-transfer processes, and (3) heat and moisture fluxes from the ocean surface. A comparative study is made of the vertical transport of horizontal momentum and of the latent heat release by the simulated subtropical squall system, as well as the heat release of squall systems occurring in other geographic locations.

Tao, Wei-Kuo↗
NASA NTRS1992

An investigation of the information propagation and entropy transport aspects of Stirling machine numerical simulation

Aspects of the information propagation modeling behavior of integral machine computer simulation programs are investigated in terms of a transmission line. In particular, the effects of pressure-linking and temporal integration algorithms on the amplitude ratio and phase angle predictions are compared against experimental and closed-form analytic data. It is concluded that the discretized, first order conservation balances may not be adequate for modeling information propagation effects at characteristic numbers less than about 24. An entropy transport equation suitable for generalized use in Stirling machine simulation is developed. The equation is evaluated by including it in a simulation of an incompressible oscillating flow apparatus designed to demonstrate the effect of flow oscillations on the enhancement of thermal diffusion. Numerical false diffusion is found to be a major factor inhibiting validation of the simulation predictions with experimental and closed-form analytic data. A generalized false diffusion correction algorithm is developed which allows the numerical results to match their analytic counterparts. Under these conditions, the simulation yields entropy predictions which satisfy Clausius' inequality.

Goldberg, Louis F.↗
NASA NTRS1983

Numerical simulation of the atmosphere during a CAT encounter

In an attempt to determine whether clear air turbulence (CAT) which caused a brief, intense turbulence encounter to a jetliner might have been related to gravity waves generated by a line of thunderstorms, local sounding data taken within two hours of the incident have been used as input to a numerical model designed to simulate stratified flow over obstacles. In the model a rigid obstacle at the lower boundary acted as the source of the gravity waves. The model results show a large amplitude disturbance localized over the obstacle at and above flight level, in regions where CAT would be expected. Stratification studies have been conducted, varying the height and shape of the obstacle and the input wind profile. The success of the experiment leads to the conclusion that atmospheric gravity wave simulation can be useful in understanding and possibly avoiding CAT.

Keller, T. L.↗
NASA NTRS1980

Numerical simulation of flow over wavy walls

Spectral methods are used to solve the Navier-Stokes equations numerically for two-dimensional incompressible flow over wavy surfaces. Results obtained for high Reynolds number laminar flows and turbulent flows are compared to experiment and analytical theory. Overall, good agreement has been achieved.

Balasubramanian, R.↗
NASA NTRS1987

Numerical simulation of transition, compressible turbulence, and reacting flows

Some of the recent work at NASA Langley on transition, turbulence, and reacting flows is summarized. Much of this effort has been motivated by outstanding technological problems in high-speed flow. A class of numerical algorithms for these problems has been developed and a variety of physical problems have been simulated. Descriptions are provided of the basic mathematical models, the nature of the numerical methods, and some of the recent simulations.

Zang, T. A.↗
DOE OSTI2024

Numerical simulation of RHIC polarized proton run 17 spin flipper experiments

RHIC nine-magnet spin flipper has been operated successfully during RHIC polarized proton Run 17, with 97% spin flip efficiency achieved. The results show the importance of mirror resonance removal, small spin tune spread, and proper spin flipper driving tune sweep speed. Detailed spin tracking simulations, based on a Lorentz force and Thomas-BMT differential equation numerical solver code for accuracy, have been carried out to understand the experimental results. Agreement within measurement accuracy is obtained at injection energy, 23.8 GeV. It is not as tight at 255 GeV, reasons for that are exposed. These measurements and numerical studies allow to determine the sensitivity of spin-flip efficiency to the dispersion slopes at the two Siberian snakes and to the ac dipole frequency sweep speed. They also provide guidance for future developments at BNL’s electron-ion collider. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Numerical Simulation of Light to Heat Conversion by Plasmonic Nanoheaters

Plasmonic nanoparticles are widely recognized as photothermal conversion agents, i.e., nanotransducers or nanoheaters. Translation of these materials into practical applications requires quantitative analyses of their photothermal conversion efficiencies (η). However, the value of η obtained for different materials is dramatically influenced by the experimental setup and method of calculation. Here, we evaluate the most common methods for estimating η (Roper’s and Wang’s) and compare these with numerical estimates using the simulation software ANSYS. Experiments were performed with colloidal gold nanorod solutions suspended in a hanging droplet irradiated by an 808 nm diode laser and monitored by a thermal camera. The ANSYS simulations accounted for both heating and evaporation, providing η values consistent with the Wang method but higher than the Roper approach. This study details methods for estimating the photothermal efficiency and finds ANSYS to be a robust tool where experimental constraints complicate traditional methods.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
NASA NTRS1992

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

Work performed during the report period is summarized. The first numerical experiment which was performed on the North Carolina Supercomputer Center's CRAY-YMP machine during the second half of FY92 involved a 36 hour simulation of the CCOPE case study. This first coarse-mesh simulation employed the GMASS model with a 178 x 108 x 32 matrix of grid points spaced approximately 24 km apart. The initial data was comprised of the global 2.5 x 2.5 degree analyses as well as all available North American rawinsonde data valid at 0000 UTC 11 July 1981. Highly-smoothed LFM-derived terrain data were utilized so as to determine the mesoscale response of the three-dimensional atmosphere to weak terrain forcing prior to including the observed highly complex terrain of the northern Rocky Mountain region. It was felt that the model should be run with a spectrum of terrain geometries, ranging from observed complex terrain to no terrain at all, to determine how crucial the terrain was in forcing the mesoscale phenomena. Both convection and stratiform (stable) precipitation were not allowed in this simulation so that their relative importance could be determined by inclusion in forth-coming simulations. A full suite of planetary boundary layer forcing was allowed in the simulation, including surface sensible and latent heat fluxes employing the Blakadar PBL formulation. The details of this simulation, which in many ways could be considered the control simulation, including the important synoptic-scale, meso-alpha scale, and meso-beta scale circulations is described. These results are compared to the observations diagnosed by Koch and his colleagues as well as hypotheses set forth in the project proposal for terrain-influences upon the jet stream and their role in the generation of mesoscale wave phenomenon. The fundamental goal of the analyses being the discrimination among background geostrophic adjustment, terrain influences, and shearing instability in the initiation and maintainance of mesoscale internal wave phenomena. Based upon these findings, FY93 plans are discussed. A review of linear theory and theoretical modeling of a geostrophic zonal wind anomaly is included.

Lin, Yuh-Lang↗
NASA NTRS1987

Numerical simulation of a plane turbulent mixing layer, with applications to isothermal, rapid reactions

A hybrid method has been developed for the numerical prediction of turbulent mixing in a spatially-developing, free shear layer. Most significantly, the computation incorporates the effects of large-scale structures, Schmidt number and Reynolds number on mixing, which have been overlooked in the past. In flow field prediction, large-eddy simulation was conducted by a modified 2-D vortex method with subgrid-scale modeling. The predicted mean velocities, shear layer growth rates, Reynolds stresses, and the RMS of longitudinal velocity fluctuations were found to be in good agreement with experiments, although the lateral velocity fluctuations were overpredicted. In scalar transport, the Monte Carlo method was extended to the simulation of the time-dependent pdf transport equation. For the first time, the mixing frequency in Curl's coalescence/dispersion model was estimated by using Broadwell and Breidenthal's theory of micromixing, which involves Schmidt number, Reynolds number and the local vorticity. Numerical tests were performed for a gaseous case and an aqueous case. Evidence that pure freestream fluids are entrained into the layer by large-scale motions was found in the predicted pdf. Mean concentration profiles were found to be insensitive to Schmidt number, while the unmixedness was higher for higher Schmidt number. Applications were made to mixing layers with isothermal, fast reactions. The predicted difference in product thickness of the two cases was in reasonable quantitative agreement with experimental measurements.

Lin, P.↗
NASA NTRS1986

Numerical simulation of the flowfield over ice accretion shapes

The primary goals are directed toward the development of a numerical method for computing flow about ice accretion shapes and determining the influence of these shapes on flow degradation. It is expedient to investigate various aspects of icing independently in order to assess their contribution to the overall icing phenomena. The specific aspects to be examined include the water droplet trajectories with collection efficiencies and phase change on the surface, the flowfield about specified shapes including lift, drag, and heat transfer distribution, and surface roughness effects. The configurations computed were models of ice accretion shapes formed on a circular cylinder in the NASA Lewis Icing Research Tunnel. An existing Navier-Stokes program was modified to compute the flowfield over four shapes (2, 5, and 15 minute models of glaze ice, and a 15 minute accumulation of rime ice).

Scott, J. N.↗
NASA NTRS1992

Numerical simulation and analysis of a prefrontal squall line. I - Observations and basic simulation results. II - Propagation of the squall line as an internal gravity wave

Results of an observational and numerical study of the squall line that occurred on June 17-18, 1978, are presented. This squall line was initially triggered by the strong surface convergence along a cold front and stretched from Illinois to the Texas Panhandle. The squall line was aligned with the surface front during its initial development, but then propagated faster than the front, resulting in a separation of approximately 200 km by 0300 UTC and 300-400 km by 0600 UTC. The Colorado State University Regional Atmospheric Modeling System is used to model the squall-line development and propagation. Results are described from several experiments that tested the sensitivity to the use of the Kuo-type cumulus parameterization scheme and grid-scale microphysical processes. The movement of the squall line in the model is shown to be due to the propagation of a deep tropospheric internal gravity wave in a CISK-like process. The thermal and dynamic perturbations associated with the hypothesized wave are shown to be consistent with internal gravity wave theory.

Cram, Jennifer M.↗