SEARCH · Engineering Papers

Results for “Direct 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,369 records · Page 76

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Spectral multigrid methods for the solution of homogeneous turbulence problems

New three-dimensional spectral multigrid algorithms are analyzed and implemented to solve the variable coefficient Helmholtz equation. Periodicity is assumed in all three directions which leads to a Fourier collocation representation. Convergence rates are theoretically predicted and confirmed through numerical tests. Residual averaging results in a spectral radius of 0.2 for the variable coefficient Poisson equation. In general, non-stationary Richardson must be used for the Helmholtz equation. The algorithms developed are applied to the large-eddy simulation of incompressible isotropic turbulence.

Erlebacher, G.↗

NASA NTRS1989

Time-dependent computational studies of flames in microgravity

The research performed at the Center for Reactive Flow and Dynamical Systems in the Laboratory for Computational Physics and Fluid Dynamics, at the Naval Research Laboratory, in support of the NASA Microgravity Science and Applications Program is described. The primary focus was on investigating fundamental questions concerning the propagation and extinction of premixed flames in Earth gravity and in microgravity environments. The approach was to use detailed time-dependent, multispecies, numerical models as tools to simulate flames in different gravity environments. The models include a detailed chemical kinetics mechanism consisting of elementary reactions among the eight reactive species involved in hydrogen combustion, coupled to algorithms for convection, thermal conduction, viscosity, molecular and thermal diffusion, and external forces. The external force, gravity, can be put in any direction relative to flame propagation and can have a range of values. A combination of one-dimensional and two-dimensional simulations was used to investigate the effects of curvature and dilution on ignition and propagation of flames, to help resolve fundamental questions on the existence of flammability limits when there are no external losses or buoyancy forces in the system, to understand the mechanism leading to cellular instability, and to study the effects of gravity on the transition to cellular structure. A flame in a microgravity environment can be extinguished without external losses, and the mechanism leading to cellular structure is not preferential diffusion but a thermo-diffusive instability. The simulations have also lead to a better understanding of the interactions between buoyancy forces and the processes leading to thermo-diffusive instability.

Oran, Elaine S.↗

NASA NTRS1994

Performance issues for iterative solvers in device simulation

Due to memory limitations, iterative methods have become the method of choice for large scale semiconductor device simulation. However, it is well known that these methods still suffer from reliability problems. The linear systems which appear in numerical simulation of semiconductor devices are notoriously ill-conditioned. In order to produce robust algorithms for practical problems, careful attention must be given to many implementation issues. This paper concentrates on strategies for developing robust preconditioners. In addition, effective data structures and convergence check issues are also discussed. These algorithms are compared with a standard direct sparse matrix solver on a variety of problems.

Fan, Qing↗

NASA NTRS1984

Ballistic transport in Saturn's rings - An analytic theory

Ejecta from impacts of micrometeoroids on Saturn's ring particles will, in most cases, remain in orbit about Saturn and eventually be reaccreted by the rings, possibly at a different radial location. The resulting mass transport has been suggested as the cause of some of the features observed in Saturn's rings. Previous attempts to model this transport have used numerical simulations which have not included the effects of the angular momentum transport coincident with mass transport. An analytical model for ballistic mass transport in Saturn's rings is developed. The model includes the effects of angular momentum advection and shows that the net material movement due to angular momentum advection is comparable to that caused by direct ballistic mass transport.

Lissauer, J. J.↗

NASA NTRS1988

Traveling waves in large-aspect-ratio thermosolutal convection

Numerical simulations of two-dimensional thermosolutal convection in a long box of aspect ratio d = 16, with free but impenetrable boundaries, reveal solutions in the form of traveling, standing, and modulated waves. These nonlinear solutions are closely related to the theoretical predictions for an unbounded layer. An unexpected feature is that the traveling waves reverse their direction of propagation in an episodic manner; this is not the case for the modulated waves. The calculations are related to recent experiments on convection in binary fluid mixtures.

Deane, A. E.↗

NASA NTRS1988

Analysis of sea ice dynamics

The ongoing work has established the basis for using multiyear sea ice concentrations from SMMR passive microwave for studies of largescale advection and convergence/divergence of the Arctic sea ice pack. Comparisons were made with numerical model simulations and buoy data showing qualitative agreement on daily to interannual time scales. Analysis of the 7-year SMMR data set shows significant interannual variations in the total area of multiyear ice. The scientific objective is to investigate the dynamics, mass balance, and interannual variability of the Arctic sea ice pack. The research emphasizes the direct application of sea ice parameters derived from passive microwave data (SMMR and SSMI) and collaborative studies using a sea ice dynamics model. The possible causes of observed interannual variations in the multiyear ice area are being examined. The relative effects of variations in the large scale advection and convergence/divergence within the ice pack on a regional and seasonal basis are investigated. The effects of anomolous atmospheric forcings are being examined, including the long-lived effects of synoptic events and monthly variations in the mean geostrophic winds. Estimates to be made will include the amount of new ice production within the ice pack during winter and the amount of ice exported from the pack.

Zwally, J.↗

NASA NTRS1991

Advanced electromagnetic methods for aerospace vehicles

The Advanced Helicopter Electromagnetics is centered on issues that advance technology related to helicopter electromagnetics. Progress was made on three major topics: composite materials; precipitation static corona discharge; and antenna technology. In composite materials, the research has focused on the measurements of their electrical properties, and the modeling of material discontinuities and their effect on the radiation pattern of antennas mounted on or near material surfaces. The electrical properties were used to model antenna performance when mounted on composite materials. Since helicopter platforms include several antenna systems at VHF and UHF bands, measuring techniques are being explored that can be used to measure the properties at these bands. The effort on corona discharge and precipitation static was directed toward the development of a new two dimensional Voltage Finite Difference Time Domain computer program. Results indicate the feasibility of using potentials for simulating electromagnetic problems in the cases where potentials become primary sources. In antenna technology the focus was on Polarization Diverse Conformal Microstrip Antennas, Cavity Backed Slot Antennas, and Varactor Tuned Circular Patch Antennas. Numerical codes were developed for the analysis of two probe fed rectangular and circular microstrip patch antennas fed by resistive and reactive power divider networks.

Balanis, Constantine A.↗

NASA NTRS2020

Impacts of Water Clarity Variability on Temperature and Biogeochemistry in the Chesapeake Bay

Estuarine water clarity depends on the concentrations of aquatic constituents, such as colored dissolved organic matter, phytoplankton, inorganic suspended solids, and detritus, which are influenced by variations in riverine inputs. These constituents directly affect temperature because when water is opaque, sunlight heats a shallower layer of the water compared to when it is clear. Despite the importance of accurately predicting temperature variability, many numerical modeling studies do not adequately account for this key process. In this study, we quantify the effect of water clarity on heating by comparing two simulations of a hydrodynamic-biogeochemical model of the Chesapeake Bay for the years 2001–2005, in which (1) water clarity is constant in space and time for the computation of solar heating, compared to (2) a simulation where water clarity varies with modeled concentrations of light-attenuating materials. In the variable water clarity simulation, the water is more opaque, particularly in the northern region of the Bay. This decrease in water clarity reduces the total heat, phytoplankton, and nitrate throughout the Bay. During the spring and summer months, surface temperatures in the northern Bay are warmer by 0.1 °C and bottom temperatures are colder by 0.2 °C in the variable light attenuation simulation. Warmer surface temperatures encourage phytoplankton growth and nutrient uptake near the head of the Bay, and fewer nutrients are transported downstream. These impacts are greater during higher river flow years, when differences in temperature, nutrients, phytoplankton, and zooplankton extend further seaward compared to other years. This study demonstrates the consequences of utilizing different light calculations for estuarine heating and biogeochemistry.

water clarity↗

NASA NTRS2003

Modified Involute Helical Gears: Computerized Design, Simulation of Meshing and Stress Analysis

The computerized design, methods for generation, simulation of meshing, and enhanced stress analysis of modified involute helical gears is presented. The approaches proposed for modification of conventional involute helical gears are based on conjugation of double-crowned pinion with a conventional helical involute gear. Double-crowning of the pinion means deviation of cross-profile from an involute one and deviation in longitudinal direction from a helicoid surface. Using the method developed, the pinion-gear tooth surfaces are in point-contact, the bearing contact is localized and oriented longitudinally, and edge contact is avoided. Also, the influence of errors of alignment on the shift of bearing contact, vibration, and noise are reduced substantially. The theory developed is illustrated with numerical examples that confirm the advantages of the gear drives of the modified geometry in comparison with conventional helical involute gears.

COMPUTER AIDED DESIGN↗

NASA NTRS1995

Coronal heating by the resonant absorption of Alfven waves: Wavenumber scaling laws.

The importance of global modes in coronal loop heating is well established. In the present work the scaling of the global-mode resonant heating rate with the perturbation wavenumbers is studied with the numerical solution of the linearized time-dependent MHD equations for a full compressible, low-beta, resistive plasma using an implicit integration scheme. The numerical simulations demonstrate that the dissipation on inhomogeneties in the background Alfven speed occurs in narrow resonant layer with the highest heating rate at the global-mode frequency. The global-mode heating rate H (sub r) was found to scale as H (Sub r) approximately k (sub y) (exp 1.03) when k (sub z) = 0.1, and as H (sub r) approximately k (sub y) (exp -1.93) when k (sub z) = 0.75, where k (sub y) and k (sub z) are the wavenumbers in the perpendicular and parallel to the magnetic field directions, respectively, while the dependence of H (sub r) on k (sub z) is more complex. The quality factor Q of the MHD resonance cavity scales as Q approximately k (sub y) (exp -1.8) for k (sub z) = 0.75 and as Q approximately k (sub y) (exp -1.46) for k (sub z) = 0.1. The numerically determined heating rate scaling, the global-mode fequency, and the quality factor are in good agreement with the analytical linear theory. The magnitude of the perturbed velocities was found to decrease with k (sub y). Assuming typical coronal loop parameters (B (sub 0) = 100-200 G, upsilon (sub A) = 2000-4000 km/s), the Alfven waves can supply the required heating to a low-Q loops.

Ofman, L.↗

NASA NTRS2023

CFD Modeling of Bi-Directional PMD inside Cryogenic Propellant Tanks Onboard Parabolic Flights

Future cryogenic propulsion systems will require efficient methods with which to transfer cryogenic propellants from a depot storage tank to a customer receiver tank to minimize cost and maximize reusability. The Reduced Gravity Cryogenic Transfer project is currently developing advanced cryogenic fluid management technology and developing and validating new numerical models for three phases of transfer: line chilldown, tank chilldown, and tank fill. Additionally, multiple liquid nitrogen (LN 2 ) parabolic flight transfer rigs are being designed by universities and NASA to investigate the gravitational sensitivities that exist in these three technologies. In order to maximize the collection of low-g data during flights, it is required to extract as much (LN 2 as possible from the supply tank, despite variable gravity levels. The purpose of this paper is to present computational fluid dynamics (CFD) volume of fluid simulations of (LN 2 behavior in the supply tank onboard parabolic flights to validate the optimal design of a bi-directional propellant management device (PMD) using the commercial software FLOW-3D. A parametric study is conducted on the effects of gravity level, fill level, pore size, open area, thickness, and type of baffle on PMD performance. Based on results, the PMD as designed exceeds the targeted expulsion efficiency.

Jason Hartwig↗

NASA NTRS2011

Impacts of Aerosol Direct Effects on the South Asian Climate: Assessment of Radiative Feedback Processes Using Model Simulations and Satellite/Surface Measurements

Current assessment of aerosol radiative effect is hindered by our incomplete knowledge of aerosol optical properties, especially absorption, and our current inability to quantify physical and microphysical processes. In this research, we investigate direct aerosol radiative effect over heavy aerosol loading areas (e.g., Indo-Gangetic Plains, South/East Asia) and its feedbacks on the South Asian climate during the pre-monsoon season (March-June) using the Purdue Regional Climate Model (PRCM) with prescribed aerosol data derived by the NASA Goddard Earth Observing System Model (GEOS-5). Our modeling domain covers South and East Asia (60-140E and 0-50N) with spatial resolutions of 45 km in horizontal and 28 layers in vertical. The model is integrated from 15 February to 30 June 2008 continuously without nudging (i.e., only forced by initial/boundary conditions). Two numerical experiments are conducted with and without the aerosol-radiation effects. Both simulations are successful in reproducing the synoptic patterns on seasonal-to-interannual time scales and capturing a pre-monsoon feature of the northward rainfall propagation over Indian region in early June which shown in Tropical Rainfall Measuring Mission (TRMM) observation. Preliminary result suggests aerosol-radiation interactions mainly alter surface-atmosphere energetics and further result in an adjustment of the vertical temperature distribution in lower atmosphere (below 700 hPa). The modifications of temperature and associated rainfall and circulation feedbacks on the regional climate will be discussed in the presentation.

Wang, Sheng-Hsiang↗

NASA NTRS1994

A forward-reverse shock pair in the solar wind driven by over-expanison of a coronal mass ejection: Ulysses observations

A previously unidentified type of solar wind forward-reverse shock pair has been observed by Ulysses at 4.64 AU and S32.5 deg. In contrast to most solar wind forward-reverse shock pairs, which are driven by the speed difference between fast solar wind plasma and slower plasma ahead, this particular shock pair was driven purely by the over-expansion of a coronal mass ejection (CME) in transit from the Sun. A simple numerical simulation indicates that the over-expansion was a result of a high initial internal plasma and magnetic field pressure within the CME. The CME observed at 4.64 AU had the internal field structure of a magnetic flux rope. This event was associated with a solar disturbance in which new magnetic loops formed in the corona almost directly beneath Ulysses approximately 11 days earlier. This association suggests that the flux rope was created as a result of reconnection between the 'legs' of neighboring magnetic loops within the rising CME.

Gosling, J. T.↗

NASA NTRS1996

Internal Passage Heat Transfer Prediction Using Multiblock Grids and a Kappa-Omega Turbulence Model

Numerical simulations of the three-dimensional flow and heat transfer in a rectangular duct with a 180 C bend were performed. Results are presented for Reynolds numbers of 17,000 and 37,000 and for aspect ratios of 0.5 and I.O. A kappa-omega turbulence model with no reference to distance to a wall is used. Direct comparison between single block and multiblock grid calculations are made. Heat transfer and velocity distributions are compared to available literature with good agreement. The multi-block grid system is seen to produce more accurate results compared to a single-block grid with the same number of cells.

Rigby, David L.↗

NASA NTRS2001

Global MHD Models of the Solar Corona

Global magnetohydrodynamic (MHD) models of the solar corona are computationally intensive, numerically complex simulations that have produced important new results over the past few years. After a brief overview of how these models usually work, I will address three topics: (1) How these models are now routinely used to predict the morphology of the corona and analyze Earth and space-based remote observations of the Sun; (2) The direct application of these models to the analysis of physical processes in the corona and chromosphere and to the interpretation of in situ solar wind observations; and (3) The use of results from global models to validate the approximations used to make detailed studies of physical processes in the corona that are not otherwise possible using the global models themselves.

Suess, S. T.↗

NASA NTRS2002

Numerical Simulations for Landing Gear Noise Generation and Radiation

Aerodynamic noise from a landing gear in a uniform flow is computed using the Ffowcs Williams -Hawkings (FW-H) equation. The time accurate flow data on the surface is obtained using a finite volume flow solver on an unstructured and. The Ffowcs Williams-Hawkings equation is solved using surface integrals over the landing gear surface and over a permeable surface away from the landing gear. Two geometric configurations are tested in order to assess the impact of two lateral struts on the sound level and directivity in the far-field. Predictions from the Ffowcs Williams-Hawkings code are compared with direct calculations by the flow solver at several observer locations inside the computational domain. The permeable Ffowcs Williams-Hawkings surface predictions match those of the flow solver in the near-field. Far-field noise calculations coincide for both integration surfaces. The increase in drag observed between the two landing gear configurations is reflected in the sound pressure level and directivity mainly in the streamwise direction.

Morris, Philip J.↗

NASA NTRS2016

Overview of the CHarring Ablator Response (CHAR) Code

An overview of the capabilities of the CHarring Ablator Response (CHAR) code is presented. CHAR is a one-, two-, and three-dimensional unstructured continuous Galerkin finite-element heat conduction and ablation solver with both direct and inverse modes. Additionally, CHAR includes a coupled linear thermoelastic solver for determination of internal stresses induced from the temperature field and surface loading. Background on the development process, governing equations, material models, discretization techniques, and numerical methods is provided. Special focus is put on the available boundary conditions including thermochemical ablation and contact interfaces, and example simulations are included. Finally, a discussion of ongoing development efforts is presented.

Amar, Adam J.↗

NASA NTRS1973

Numerical calculation of the three dimensional transonic flow over a yawed wing.

Results are presented of calculations of the three dimensional steady transonic flow over a finite yawed wing. The full potential flow equation is solved in a transformed coordinate system which permits the boundary conditions to be satisfied exactly. The correct differential properties are enforced by rotating the difference scheme to conform with the flow direction, and fast convergence is assured by simulating a time dependent equation designed to settle quickly to a steady state. Computed lift drag ratios are consistent with the results of wind tunnel tests of a yawed wing conducted by R. T. Jones (1972).-

Jameson, A.↗