Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 1,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

An Approach for Simulating Transport of Nonspherical, Sublimating Particles in Rarefied Gas Flows via the Direct Simulation Monte Carlo Method

This paper extends the particle Green’s function approach of Gallis et al. [Phys. Fluids 13, 3482–3492 (2001)] to compute the force and heat transfer on nonspherical, sublimating particles in rarefied gas flows. Green’s functions are derived for rods and plates, which are representative of ice particle shapes, and a numerical procedure is developed for computing the Green’s functions for a general discretized convex shape. Additionally, a sublimation model is developed to simulate the evolution of ice particles in rarefied gas flows. These models are implemented in a direct simulation Monte Carlo code and demonstrated for the case of ice-particle-laden hypersonic flow over a ramp.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1990

Fluid dynamics and low gravity effects of chemical vapor deposition

Based on the comparison between experimental data and numerical results for the growth of GaAs from TMGa, it was shown that 3D simulations are necessary to simulate rectangular CVD reactors even when operated under subcritical (Ra) conditions. The important points found are summarized in the three attached reprints. The experimental studies of mixed convection in horizontal channels have shown three regimes of high Ra (22,220) number flows. At Re = 18.5, the rolls develop very quickly, significantly modulating the axial velocity even before it reaches the beginning of the hot plate. A few centimeters downstream, the velocities become asymmetric about the vertical centerplane and at x = 12 cm, become unsteady. These asymmetries were predicted theoretically, but experimental evidence has not been published prior to this work. At Re = 36, the axial velocity is only slightly modified at x = 0. Although the flow remains steady and symmetric about the vertical centerplane, there is a small spatial oscillation in the velocities over the length of the channel. The period of this oscillation was around 5 cm. At Re = 54, the longitudinal rolls developed smoothly over a length of 30 cm, with no asymmetries, unsteadiness, or spatial oscillations. Comparison of numerical simulations of these flows to experiments has revealed the importance and difficulty of setting proper thermal boundary conditions on the sidewalls. Calculated flows and experimentally measured flows showed very similar profiles, but at different axial locations, with the rolls developing more rapidly in the experiments. This is directly attributable to partially conducting sidewalls of the apparatus being hotter in the entrance section than the adiabatic walls of the simulations. A thorough comparison of the experimental data and numerical results for a variety of sidewall boundary conditions is in preparation.

Nyce, Thomas A.↗
DOE OSTI2017

Bi-resonant klynac

Here in this work, we present the first analysis of a new device called the bi-resonant klynac, which is a combined klystron and linac. In a bi-resonant klynac, all RF cells, except for the cell that acts as the input for the klystron section, belong to a single resonant circuit. This resonant coupling configuration leads to increased operational stability and can tolerate significant temperature variations. In this paper, a basic analysis of this device is presented, including discussions of how it operates and of the advantages of resonantly coupling the RF generation directly to the linac. We additionally describe the approach used to numerically model the klynac and we include detailed simulations of a 50-kV, 10-A klynac that produces a 1-MeV, 0.1-A output beam. In conclusion, this type of device may be especially useful for situations where an electron beam is needed at low cost.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

PSCF+: An Extended and Improved Open-Source Software Package for Polymer Self-Consistent Field Calculations

This work introduces PSCF+, a freely available open‑source software tool for modeling how block copolymers self‑assemble into ordered nanoscale structures. The software allows researchers to study a wide range of polymer models and interaction types, making it possible to directly compare theoretical predictions with molecular and simulation‑based studies in a consistent way. PSCF+ incorporates several recently developed numerical methods that significantly reduce computer memory requirements and speed up calculations, enabling efficient use of modern GPU hardware. Here, it also provides automated tools to calculate free‑energy landscapes and phase boundaries, which are essential for understanding and predicting material behavior. Together, these advances make PSCF+ a powerful and flexible platform for studying polymer self‑assembly and support the design of new polymeric materials for scientific and technological applications.

Chemical calculations↗
DOE OSTI2024

Misaligned Wind‐Waves Behind Atmospheric Cold Fronts

Abstract Atmospheric fronts embedded in extratropical cyclones are high‐impact weather phenomena, contributing significantly to mid‐latitude winter precipitation. The three vital characteristics of the atmospheric fronts, high wind speeds, abrupt change in wind direction, and rapid translation, force the induced surface waves to be misaligned with winds exclusively behind the cold fronts. The effects of the misaligned waves under atmospheric cold fronts on air‐sea fluxes remain undocumented. Using the multi‐year in situ near‐surface observations and direct covariance flux measurements from the Pioneer Array off the coast of New England, we find that the majority of the passing cold fronts generate misaligned waves behind the cold front. Once generated, the waves remain misaligned, on average, for about 8 hr. The parameterized effect of misaligned waves in a fully coupled model significantly increases the roughness length (185%), drag coefficient (19%), and air‐sea momentum flux (11%). The increased surface drag reduces the wind speeds in the surface layer. The upward turbulent heat flux is weakly decreased by the misaligned waves because of the decrease in temperature and humidity scaling parameters being greater than the increase in friction velocity. The misaligned wave effect is not accurately represented in a commonly used wave‐based bulk flux algorithm. Yet, considering this effect in the current formulation improves the overall accuracy of parameterized momentum flux estimates. The results imply that better representing a directional wind‐wave coupling in the bulk formula of the numerical models may help improve the air‐sea interaction simulations under the passing atmospheric fronts in the mid‐latitudes.

Oceanography↗
DOE OSTI2024

Towards Construction of a Novel Nanometer-Resolution MeV-STEM for Imaging Thick Frozen Biological Samples

Driven by life-science applications, a mega-electron-volt Scanning Transmission Electron Microscope (MeV-STEM) has been proposed here to image thick frozen biological samples as a conventional Transmission Electron Microscope (TEM) may not be suitable to image samples thicker than 300–500 nm and various volume electron microscopy (EM) techniques either suffering from low resolution, or low speed. The high penetration of inelastic scattering signals of MeV electrons could make the MeV-STEM an appropriate microscope for biological samples as thick as 10 μm or more with a nanoscale resolution, considering the effect of electron energy, beam broadening, and low-dose limit on resolution. The best resolution is inversely related to the sample thickness and changes from 6 nm to 24 nm when the sample thickness increases from 1 μm to 10 μm. To achieve such a resolution in STEM, the imaging electrons must be focused on the specimen with a nm size and an mrad semi-convergence angle. This requires an electron beam emittance of a few picometers, which is ~1000 times smaller than the presently achieved nm emittance, in conjunction with less than 10 -4 energy spread and 1 nA current. We numerically simulated two different approaches that are potentially applicable to build a compact MeV-STEM instrument: (1) DC (Direct Current) gun, aperture, superconducting radio-frequency (SRF) cavities, and STEM column; (2) SRF gun, aperture, SRF cavities, and STEM column. Beam dynamic simulations show promising results, which meet the needs of an MeV-STEM, a few-picometer emittance, less than 10 -4 energy spread, and 0.1–1 nA current from both options. Also, we designed a compact STEM column based on permanent quadrupole quintuplet, not only to demagnify the beam size from 1 μm at the source point to 2 nm at the specimen but also to provide the freedom of changing the magnifications at the specimen and a scanning system to raster the electron beam across the sample with a step size of 2 nm and the repetition rate of 1 MHz. This makes it possible to build a compact MeV-STEM and use it to study thick, large-volume samples in cell biology.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Measuring Black Hole Light Echoes with Very Long Baseline Interferometry

Light passing near a black hole can follow multiple paths from an emission source to an observer due to strong gravitational lensing. Photons following different paths take different amounts of time to reach the observer, which produces an echo signature in the image. The characteristic echo delay is determined primarily by the mass of the black hole, but it is also influenced by the black hole spin and inclination to the observer. In the Kerr geometry, echo images are demagnified, rotated, and sheared copies of the direct image and lie within a restricted region of the image. Echo images have exponentially suppressed flux, and temporal correlations within the flow make it challenging to directly detect light echoes from the total light curve. In this Letter, we propose a novel method to search for light echoes by correlating the total light curve with the interferometric signal at high spatial frequencies, which is a proxy for indirect emission. We explore the viability of our method using numerical general relativistic magnetohydrodynamic simulations of a near-face-on accretion system scaled to M87-like parameters. We demonstrate that our method can be used to directly infer the echo delay period in simulated data. An echo detection would be clear evidence that we have captured photons that have circled the black hole, and a high-fidelity echo measurement would provide an independent measure of fundamental black hole parameters. Our results suggest that detecting echoes may be achievable through interferometric observations with a modest space-based very long baseline interferometry mission.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1985

Simulation of Multistage Turbine Flows

The numerical simulation of turbine flows serves to enhance the understanding of the flow phenomena within multistage turbomachinery components. The direct benefit of this activity is improved modeling capability, which can be used to improve component efficiency and durability. A hierarchy of equations was formulated to assess the difficulty in analyzing the flow field within multistage turbomachinery components. The Navier-Stokes equations provides the most complete description. The simplest description is given by a set of equations that govern the quasi-one-dimensional flow. The number of unknowns to be solved for increases monotonically above the number of equations. The development of the additional set of equations needed to mathematically close the system of equations forms the closure problem associated with that level of description. For the Navier-Stokes equation there is no closure problem. For the quasi-one-dimensional equation set random flow fluctuations, unsteady fluctuations, nonaxisymmetric flow variations, and hub-to-shroud variations on the quasi-one-dimensional flow must be accounted for.

Celestina, M. L.↗
NASA NTRS1989

Computation of rarefied hypersonic flows

Numerical techniques for the simulation of hypersonic flows of rarefied gases are examined in an analytical review. The direct-simulation Monte Carlo (DSMC) method developed to interpret measurement data obtained by the Space Shuttle in the SUMS project is described; the fundamental limitations of the DMSC approach are discussed; and modifications to improve the physical plausibility of DSMC predictions are proposed. Particular attention is given to the use of a double-peaked molecular distribution function for the internal flow in the SUMS probe, a downstream vacuum-reservoir condition as a simplifying assumption, and an explicit forward-time centered-space differencing scheme for the discretization of the SUMS problem. Typical simulation results are presented in extensive graphs and briefly characterized.

Cheng, Sin-I↗
NASA NTRS1993

Near wake structure for a generic ASTV configuration

Results of a numerical study are presented for hypersonic low-density flow about a 70-deg blunt cone using the direct simulation Monte Carlo method. Particular emphasis is given to the near wake flow and its sensitivity to rarefaction and other parametric variations. The flow conditions simulated are attainable in existing low-density hypersonic wind tunnels; that is, Mach 20 nitrogen flow encompassing freestream Knudsen numbers of 0.03 to 0.001. A stable vortex forms in the near wake at and below a freestream Knudsen number of 0.01 and the size of the vortex increases with decreasing freestream Knudsen number. The base region of the flow remains in thermal nonequilibrium for all cases. There is no formation of a lip separation shock or a distinct wake shock at these rarefied conditions.

Dogra, Virendra K.↗
NASA NTRS1999

Bubble Dynamics on a Heated Surface

In this work, we study steady and oscillatory thermocapillary and natural convective flows generated by a bubble on a heated solid surface. The interaction between gas and vapor bubbles with the surrounding fluid is of interest for both space and ground-based processing. A combined numerical-experimental approach is adopted here. The temperature field is visualized using Mach-Zehnder and/or Wollaston Prism Interferometry and the flow field is observed by a laser sheet flow visualization technique. A finite element numerical model is developed which solves the transient two-dimensional continuity, momentum, and energy equations and includes the effects of temperature-dependent surface tension and bubble surface deformation. Below the critical Marangoni number, the steady state low-g and 1-g temperature and velocity fields predicted by the finite element model are in excellent agreement with both the visualization experiments in our laboratory and recently published experimental results in the literature. Above the critical Marangoni number, the model predicts an oscillatory flow which is also closely confirmed by experiments. It is shown that the dynamics of the oscillatory flow are directly controlled by the thermal and hydrodynamic interactions brought about by combined natural and thermocapillary convection. Therefore, as numerical simulations show, there are considerable differences between the 1-g and low-g temperature and flow fields at both low and high Marangoni numbers. This has serious implications for both materials processing and fluid management in space.

Kassemi, M.↗
NASA NTRS1994

Regional-scale flows in mountainous terrain. Part 1: A numerical and observational comparison

This study uses observed data and a numerical simulation to examine the generation of thermally driven flows across the Colorado mountain barrier on meso-beta to meso-alpha scales. The observations were collected from remote surface observing systems at exposed mountaintop locations throughout the state of Colorado, over the summers of 1984-88, as part of the Rocky Mountain Peaks Experiment (ROMPEX). The data show the development of a recurrent circulation system across the Colorado mountain barrier, operating on a diurnal timescale. From the observations, the basic structure of the flow system appears as a daytime inflow toward the highest terrain, and a nocturnal outflow away from it. However, when examined in detail, the flow system exhibits more unusual behavior, especially west of the barrier crest. Here, winds in the early evening are occasionally observed to onset abruptly from an easterly direction, generally counter to the upper-level winds. Observations from ROMPEX for 26 August 1985 are used to provide comparison data for a numerical simulation with the Regional Atmospheric Modeling System (RAMS). This three-dimensional case study experiment is initialized with data from the National Meteorological Center and incorporates two-way interactive grid nesting. From the observed data and case study simulation, four distinct phases of the regional-scale circulation system have been identified. In the development phase, a deep mountain-plains solenoid is generated through terrain heating along the Front Range. This circulation system transforms in the late afternoon transition phase into a westward-propagating density current (WPDC). The third phase, called the 'density-current propagation phase,' occurs as the WPDC moves westward across the mountains, leaving in its wake strong southeasterly flow at the mountaintop level. This current appears to be the cause of the peculiar easterly component winds found in the ROMPEX mountaintop observations along the western slope. In the final late-night adjustment phase, the WPDC dissipates near the western edge of the Colorado mountains and steady southerly flow evolves over the high mountain terrain. This southerly flow is the steady response to the differential heating that develops between the low-lying plains and the intermountain region.

Bossert, James E.↗
NASA NTRS1982

Scattering and collapse of Langmuir waves driven by a weak electron beam

Wave-wave interactions are examined, particularly for significance in relation to Langmuir turbulence associated with type III solar radio emissions. The Zakharov equations are considered as two coupled nonlinear partial differential equations for the low-frequency density and the Langmuir envelope, and the electron and ion temperatures are regarded as equal. The Zhakharov equations are solved in one dimension, and display no collapse of a driven broadband pump, indicating that scattering instabilities are important in the saturation of a beam-plasma instability. A two-dimensional numerical simulation is developed to interpret the role of parametric instabilities, and it is concluded that the distortion of a wave packet is due to scattering instability, and not a change in the physics of direct collapse due to an artificial numerical damping. Further discussions are presented for comparisons of models of Langmuir turbulence in type III bursts.

Hafizi, B.↗
NASA NTRS1998

3-D Modeling of Directional Solidification of a Non-Dilute Alloy with Temperature and Concentration Fields Coupling via Materials Properties Dependence and via Double Diffusive Convection

Numerical simulation of the HgCdTe growth by the vertical Bridgman method was performed using FIDAP finite element code. Double-diffusive melt convection is analyzed, as the primary factor at controls inhomogeneity of the solidified material. Temperature and concentration fields in the model are also coupled via material properties, such as thermal and solutal expansion coefficients with the dependence on both temperature and concentration, and melting temperature evaluation from pseudobinary CdTe-HgTe phase diagram. Experimental measurements were used to obtain temperature boundary conditions. Parametric study of the melt convection dependence on the gravity conditions was undertaken. It was found, that the maximum convection velocity in the melt can be reduced under certain conditions. Optimal conditions to obtain a near flat solidified interface are discussed. The predicted interface shape is in agreement with one obtained experimentally by quenching. The results of 3-D calculations are compared with previous 2- D findings. A video film featuring 3-D melt convection will be presented.

Bune, Andris V.↗
NASA NTRS1989

Numerical investigation of the control of the secondary instability process in boundary layers

A numerical model has been developed for investigating boundary layer transition control for a flat plate boundary layer. Active control of a periodically forced boundary layer in an incompressible fluid is studied using surface heating techniques. The spatially evolving boundary layer is simulated. The Navier-Stokes and energy equations are integrated using a fully implicit finite difference/spectral method. Temperature perturbations are introduced locally along finite heater strips to directly attenuate the instability waves in the flow. A feedback control loop is employed in which a downstream sensor is used to monitor wall shear stress fluctuations. Active control of small amplitude two-dimensional and three-dimensional disturbances is numerically simulated. With proper phase control, in-phase reinforcement and out-of-phase attenuation are demonstrated. A receptivity study of the localized temperature perturbations is made. It is shown that narrow heater strips are more receptive in that they maximize the amplitude level of the disturbances in the flow. Active control of the early stages of the fundamental breakdown process is also numerically simulated. Control is achieved with either two-dimensional or three-dimensional control inputs.

Kral, L. D.↗
NASA NTRS1985

Evolution of a vortical structure associated with the bursting event in a channel flow

The temporal evolution of a horseshoe vortex structure in a channel is investigated by a numerical simulation. A spectral numerical method is employed to integrate the time-dependent, three-dimensional Navier-Stokes equations. The initial vortical structure is obtained by applying a conditional sampling technique to a data base generated from a direct simulation of a turbulent channel flow. The evolution of this vortical structure under the influence of the self-induced motion and the mean shear is presented. It is shown that the initial sheet-like vortical structure rolls up into a vortex tube as it is convected downstream. Turbulence structures associated with the vortex stretching are investigated. Production of vorticity stretching is high inside the vortex legs, although it is also substantial in the tip region and above the legs. High Reynolds shear stress is produced near the tip of the vortex.

Kim, J.↗
NASA NTRS2002

Annual Research Briefs, 2002: Center for Turbulence Research

Turbulent combustion remains the largest component of the CTR's core program. This program and several related activities at CTR are supported by NASA's Ultra Efficient Engine Technology Program. It is also intimately connected with the Department of Energy's ASCI program at Stanford which develops the technology for numerical simulation of realistic aircraft engines using state of the art massively parallel computers. In combustion modeling the attention has been directed to the modeling of higher levels of complexity such as spray dynamics, radiation and soot formation. Major aircraft engine manufacturers have shown considerable interest in this program; in particular, a significant active collaboration exists between CTR and the Pratt & Whitney Corporation. CTR's combustion program is essentially based on the large-eddy simulation technique, LES, which is actively being pursued at CTR for this and many other applications. Important accomplishments in LES included simulations with three-dimensional filters which result in grid independent calculations (that is why we call it "true" LES), and the development of the methodology for integration of LES and Reynolds Averaged computations. Optimization techniques are being studied and used for the important problem of wall boundary conditions for LES as well as for optimal shape design for aeroacoustic and aerodynamic performance gains.

ULTRA EFFICIENT ENGINE TECHNOLOGY↗
NASA NTRS1990

Numerical simulation of spatially-evolving instability control in plane channel flow

The applicability of active control by periodic suction-blowing in spatially evolving plane Poiseuille flow is investigated by the direct simulations of the two-dimensional, incompressible Navier-Stokes equations. All the computations were performed for Reynolds number (Re = 7500) which is linearly unstable for this flow. The result reveal that significant reductions in perturbation amplitudes can be obtained by a proper choice of the control wave amplitude and phase even for large disturbance amplitudes. The upstream influence of the control wave is also investigated.

Danabasoglu, G.↗