Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Numerical Simulation of Flow Through an Artificial Heart

A solution procedure was developed that solves the unsteady, incompressible Navier-Stokes equations, and was used to numerically simulate viscous incompressible flow through a model of the Pennsylvania State artificial heart. The solution algorithm is based on the artificial compressibility method, and uses flux-difference splitting to upwind the convective terms; a line-relaxation scheme is used to solve the equations. The time-accuracy of the method is obtained by iteratively solving the equations at each physical time step. The artificial heart geometry involves a piston-type action with a moving solid wall. A single H-grid is fit inside the heart chamber. The grid is continuously compressed and expanded with a constant number of grid points to accommodate the moving piston. The computational domain ends at the valve openings where nonreflective boundary conditions based on the method of characteristics are applied. Although a number of simplifing assumptions were made regarding the geometry, the computational results agreed reasonably well with an experimental picture. The computer time requirements for this flow simulation, however, are quite extensive. Computational study of this type of geometry would benefit greatly from improvements in computer hardware speed and algorithm efficiency enhancements.

Rogers, Stuart E.↗
NASA NTRS1990

A time-accurate adaptive grid method and the numerical simulation of a shock-vortex interaction

A time accurate, general purpose, adaptive grid method is developed that is suitable for multidimensional steady and unsteady numerical simulations. The grid point movement is performed in a manner that generates smooth grids which resolve the severe solution gradients and the sharp transitions in the solution gradients. The temporal coupling of the adaptive grid and the PDE solver is performed with a grid prediction correction method that is simple to implement and ensures the time accuracy of the grid. Time accurate solutions of the 2-D Euler equations for an unsteady shock vortex interaction demonstrate the ability of the adaptive method to accurately adapt the grid to multiple solution features.

Bockelie, Michael J.↗
NASA NTRS1993

Numerical simulation of plasma processes driven by transverse ion heating

The plasma processes driven by transverse ion heating in a diverging flux tube are investigated with numerical simulation. The heating is found to drive a host of plasma processes, in addition to the well-known phenomenon of ion conics. The downward electric field near the reverse shock generates a doublestreaming situation consisting of two upflowing ion populations with different average flow velocities. The electric field in the reverse shock region is modulated by the ion-ion instability driven by the multistreaming ions. The oscillating fields in this region have the possibility of heating electrons. These results from the simulations are compared with results from a previous study based on a hydrodynamical model. Effects of spatial resolutions provided by simulations on the evolution of the plasma are discussed.

Singh, Nagendra↗
DOE OSTI2016

Three-dimensional direct numerical simulation of turbulent lean premixed methane combustion with detailed kinetics

The interaction of maintained homogeneous isotropic turbulence with lean premixed methane flames is investigated using direct numerical simulation with detailed chemistry. The conditions are chosen to be close to those found in atmospheric laboratory experiments. As the Karlovitz number is increased from 1 to 36, the preheat zone becomes thickened, while the reaction zone remains largely unaffected. A negative correlation of fuel consumption with mean flame surface curvature is observed. With increasing turbulence intensity, the chemical composition in the preheat zone tends towards that of an idealised unity Lewis number flame, which we argue is the onset of the transition to distributed burning, and the response of the various chemical species is shown to fall into broad classes. Smaller-scale simulations are used to isolate the specific role of species diffusion at high turbulent intensities. Diffusion of atomic hydrogen is shown to be related to the observed curvature correlations, but does not have significant consequential impact on the thickening of the preheat zone. It is also shown that susceptibility of the preheat zone to thickening by turbulence is related to the 'global' Lewis number (the Lewis number of the deficient reactant); higher global Lewis number flames tend to be more prone to thickening.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

A direct numerical simulation of Jet A flame kernel quenching

The safe operation of aeronautical engines requires an understanding of flame ignition, propagation and extinction. In this study, direct numerical simulations are performed using a 29 species reduced chemical mechanism for jet fuel surrogate Jet A to understand the flame quenching process. Here, initially laminar spherical flames of varying sizes and equivalence ratios are subject to an identical periodic domain of decaying and isotropic high intensity turbulence with a turbulent Reynolds number of 2400. All cases become quenched, except for the larger kernel with lower Karlovitz number. An analysis of the flame structure shows broadened preheat zone, flame shortening on the product side, differential species diffusion and partial fuel pyrolysis in the fresh mixture. Two extinction mechanisms are identified arising from flame shortening and high flame stretch. Flame shortening occurs due to turbulence-chemistry interactions that resemble the flame–flame interaction in a laminar counterflow reactant-to-reactant configuration, which contorts and breaks up the ignition kernel. Flame stretch is a local effect that attenuates the heat release rate and causes the flame to retreat towards the product mixtures, similar to what has been observed for reactant-to-product laminar counterflow flames. Chemical explosive mode analysis was also performed to quantify the flame structure and local combustion mode. The diffusion–reaction balance in pinched-off flame islands favors extinction of these smaller structures, while auto-ignition modes are observed within the flame kernel after fresh mixture is engulfed and preheated in the product kernel. Statistics of the density-weighted displacement speed conditional on local combustion mode indicates strong correlation between the local extinction mode and negative displacement speed. The local balance between diffusion and reaction ultimately determines the propensity for local extinction in both laminar and turbulent flames, the extent of which has an impact on global flame propagation.

42 ENGINEERING↗
NASA NTRS1986

Numerical simulations of auroral plasma processes. I - Ion beams and conics. II - Electric fields

The characteristics of ion beams and conics as seen in two-dimensional numerical simulations in which the plasmas are driven by current sheets of a finite thickness are described. It is shown that the most energetic ions in the simulations have pitch angles near 90 deg, implying a large perpendicular acceleration of the ions. Large perpendicular electric fields, similar to those measured in electrostatic shocks, are mainly confined near the edges of the current sheets, where E-perpendicular is much larger than E-parallel. Double layers with E-perpendicular about equal to, and spatially separated from, E-parallel form in the interior of wide current sheets. Such double layers have upward electric fields in the region of the upward current. Downward electric fields develop outside the current sheet.

Singh, Nagendra↗
DOE OSTI2021

Propagation of Input Uncertainties in Numerical Simulations of Laser Powder Bed Fusion

Laser powder bed fusion has the potential of redefining state-of-the-art processing and production methods, but defect formation and inconsistent build quality have limited the implementation of this process. Numerical models are widely used to study this process and predict the formation of these defects. Presently, the uncertainties of model input parameters and thermophysical properties used by these numerical simulations have not been investigated. In the present study, the uncertainty in these input parameters and material properties are quantified for laser powder bed fusion, with and without a simulated powder bed, to determine their influence on the predictive accuracy of an experimentally validated numerical model. Accounting for all possible sources of uncertainty quickly becomes computationally expensive on account of the curse of dimensionality. Uncertainty in laser absorption, solid, and liquid specific heat of the metal were found to have the largest effect on model prediction reliability with or without the use of a powder bed. Results also illustrate that accounting for these three uncertain parameters still captures the majority of model prediction uncertainty. Furthermore, the methodology of this study may be used to understand the uncertainty in as-built microstructure through propagation to microstructure prediction models, or applied under processing conditions where high Péclet numbers are observed and the thermal convection and fluid flow within the molten pool are substantial.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1982

A numerical simulation of the magnetospheric gate model for the X-ray bursters

A Lagrangian, stellar interior computer code which assumes radiation transport by diffusion is used to numerically simulate an X-ray burst occurring on a rapidly rotating, 1 solar mass neutron star accreting material from a companion in noncontinuous, episodic fashion. The blob of gas from the companion extends from the surface of the neutron star to a radius of 185 km, is optically thick, has mass of 10 to the -11th solar masses, and is initially at rest with respect to the surface. The gas initially has a temperature of about 10 million K and a surface luminosity of 0.1 that of the sun. The infall results in a burst which lasts about 0.1 second and reaches a peak luminosity and effective temperature of 240,000 suns and nine million K. The burst is followed by a phase of oscillations with a period of 0.2 second. It is found that the burst is too cool and rapid to resemble normal bursts, although the time scale is in agreement with some observations.

Starrfield, S.↗
NASA NTRS1983

Numerical simulations of loops heated to solar flare temperatures. I - Gasdynamics. II - X-ray and UV spectroscopy

The NRL's Dynamic Flux Tube Model is used to numerically simulate the dynamic response of a coronal magnetic loop to an energy input of the order encountered in solar flares. The coronal plasma is heated by the deposition of flare energy at the top of the loop to more than 10 million K, yielding a conduction front that moves toward the chromosphere, where the plasma is heated by the large downward conductive flux and ablates upward to the coronal part of the loop at velocities of a few hundred km/sec. The conduction front simultaneously produces chromospheric ablation and compresses the material ahead of it. With the aid of compressional instabilities, the compressed plasma grows throughout the flare heating phase, presenting a possible source of the flare optical continuum emission which is correlated with soft X-ray radiation. The observational consequences of rapidly heated loop gas dynamic processes are discussed. In the second part of this presentation, the dynamical calculation results previously obtained are used to predict the spectral line intensities, profiles and wavelengths of several X-ray lines and the UV line of Fe XXI at 1354.1 A. Three different viewing orientations of the loop are considered.

Cheng, C.-C.↗
NASA NTRS1996

Numerical Simulation of Forced and Free-to-Roll Delta-Wing Motions

The three-dimensional, Reynolds-averaged, Navier-Stokes (RANS) equations are used to numerically simulate nonsteady vortical flow about a 65-deg sweep delta wing at 30-deg angle of attack. Two large-amplitude, high-rate, forced-roll motions, and a damped free-to-roll motion are presented. The free-to-roll motion is computed by coupling the time-dependent RANS equations to the flight dynamic equation of motion. The computed results are in good agreement with the forces, moments, and roll-angle time histories. Vortex breakdown is present in each case. Significant time lags in the vortex breakdown motions relative to the body motions strongly influence the dynamic forces and moments.

Chaderjian, Neal M.↗
DOE OSTI2025

Mixing by internal gravity waves in stars: assessing numerical simulations against theory

ABSTRACT Here we present a study of radial chemical mixing in non-rotating massive main-sequence stars driven by internal gravity waves (IGWs), based on multidimensional hydrodynamical simulations with the fully compressible code MUSIC. We examine two proposed mechanisms of material mixing in stars by IGWs that are commonly quoted, relating to thermal diffusion and sub-wavelength shearing. Thermal diffusion provides a non-restorative effect to the waves, leaving material displaced from its previous equilibrium, while shearing arising within the waves drives weak localized flows, mixing the fluid there. Using IGW spectra from the simulations, we evaluate theoretical predictions of mixing rates due to these mechanisms. We show, for $20\, \mathrm{M}_\odot$ main-sequence stars, that neither of these mechanisms are likely to create mixing sufficient to correct inaccuracies in current stellar evolution models. Furthermore, we compare these predictions to results obtained from Lagrangian tracer particles, following a method recently used for global simulations of stellar interiors to measure mixing by IGWs in their radiative zones. We demonstrate that tracer particle methods face significant numerical challenges in measuring the small diffusion coefficients predicted by the aforementioned theories, for which they are prone to yielding artificially enhanced coefficients. Diffusion coefficients based on such methods are currently used with stellar evolution codes for asteroseismic studies, but should be viewed with caution. Finally, in a case where tracer particles do not suffer from numerical artefacts, we suggest that a diffusion model is not suitable for time-scales typically considered by 2D numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Direct numerical simulation of hot jets

The ultimate motivation of this work is to investigate the stability of two dimensional heated jets and its implications for aerodynamic sound generation from data obtained with direct numerical simulations (DNS). As pointed out in our last report, these flows undergo two types of instabilities, convective or absolute, depending on their temperature. We also described the limits of earlier experimental and theoretical studies and explained why a numerical investigation could give us new insight into the physics of these instabilities. The aeroacoustical interest of these flows was also underlined. In order to reach this goal, we first need to succeed in the DNS of heated jets. Our past efforts have been focused on this issue which encountered several difficulties. Our numerical difficulties are directly related to the physical problem we want to investigate since these absolutely or almost absolutely unstable flows are by definition very sensitive to the smallest disturbances and are very likely to reach nonlinear saturation through a numerical feedback mechanism. As a result, it is very difficult to compute a steady laminar solution using a spatial DNS. A steady state was reached only for strongly co-flowed jets, but these flows are almost equivalent to two independent mixing layers. Thus they are far from absolute instability and have much lower growth rates.

Jacob, Marc C.↗
DOE OSTI2019

Monitoring of chlorin-based photosensitiser localisation with dual-wavelength fluorescence imaging: numerical simulations

An approach to the estimation of chlorin-based photosensitiser localisation with dual-wavelength fluorescence imaging after topical application or intravenous injection is developed based on numerical simulations. Fluorescence imaging, which is commonly used for the monitoring of photosensitiser distribution within tissue, does not allow for its depth-resolved imaging. Chlorin-based photosensitisers feature two peaks in their absorption spectra corresponding to 402 and 662 nm, which allows for fluorescence excitation at two wavelengths with significantly different light penetration depths into biotissues. The ratio of fluorescence signals detected upon irradiation at different excitation wavelengths is shown to provide the information about photosensitiser localisation. The cases of exponential in-depth decay of the photosensitiser concentration and uniform photosensitiser distribution within the layer at a certain depth within the medium are analysed. These cases simulate topical application of photosensitiser and its intravenous injection and further accumulation in tumour, respectively. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Ephemeral antibubbles: Spatiotemporal evolution from direct numerical simulations

Antibubbles, which consist of a shell of a low-density fluid inside a high-density fluid, have several promising applications. We show, via extensive direct numerical simulations (DNSs), in both two and three dimensions, that the spatiotemporal evolution of antibubbles can be described naturally by the coupled Cahn-Hilliard-Navier-Stokes (CHNS) equations for a binary fluid. Our DNSs capture elegantly the gravity-induced thinning and breakup of an antibubble via the time evolution of the Cahn-Hilliard scalar-order-parameter field $\phi$, which varies continuously across interfaces, so we do not have to enforce complicated boundary conditions at the moving antibubble interfaces. To ensure that our results are robust, we supplement our CHNS simulations with sharp-interface volume-of-fluid DNSs. We track the thickness of the antibubble and calculate the dependence of the lifetime of an antibubble on several parameters; we show that our DNS results agree with various experimental results; in particular, the velocity with which the arms of the antibubble retract after breakup scales as σ 1/2 , where σ is the surface tension.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1987

Analysis of subgrid-scale eddy viscosity with use of results from direct numerical simulations

Without resort to any modeling, subgrid-scale eddy viscosity is computed from the results of high-resolution (64 cubed and 128 cubed grid points) direct numerical simulations of three-dimensional homogeneous isotropic decaying turbulence. In these simulations the eddy viscosity peaks sharply at the cutoff wavenumber, in rough agreement with the results of Kraichnan (1976). In addition, in the low-wavenumber range the eddy viscosity may be negative, contrary to the generally accepted concept of subgrid-scale eddy viscosity. Some possible explanations of this behavior are discussed.

Domaradzki, J. Andrzej↗
NASA NTRS1987

Numerical simulations of magnetospheric plasmas - Review and quadrennial report to the IUGG

This article provides a broad survey of U.S. progress during the quadrennium 1983-1986 in the category 'numerical simulations of magnetospheric plasmas'. There has been a substantial increase of activity in this area during this period. Simulations have been instrumental in providing valuable insights into large scale dynamic phenomena, nonlinear effects, and complex kinetic phenomena in a wide variety of subject areas, including shocks and double layers, ionosphere-magnetosphere coupling phenomena, and important microphysical processes such as broadband electrostatic noise. The methodology of computer simulation has also been advanced during this quadrennium. Vlasov algorithms have been improved; hybrid codes in 2 and 3D have been developed and applied to magnetospheric problems; and complex problems have been subjected with increasing frequency to a multipronged attack in which several types of simulation models, each designed to accurately model phenomena within a particular range of temporal or spatial scales, are employed synergistically.

Palmadesso, Peter J.↗
NASA NTRS1987

Numerical simulation of solar wind and magnetospheric phenomena

The nature of the nonlinear evolution of Kelvin-Helmholtz instability in the presence of sheared magnetic fields was investigated via numerical simulation. Models of the magnetosheath-magnetopause boundary in earth's tail and stream interaction regions in the inner heliosphere were studied. The development of the instability is influenced strongly by the orientation and magnitude of the magnetic field. Large vortical structures that resemble observations in the earth's tail can form while other cases generate turbulent spectra that provide insight into the generation of Alfven turbulence in the solar wind.

Goldstein, Melvyn L.↗
NASA NTRS1989

Numerical simulation of a three-dimensional turbulent boundary layer

The effects of transverse strain on an initially two-dimensional turbulent boundary layer are studied in a direct numerical simulation of a planar channel flow with impulsively started transverse pressure gradient. Consistent with experiments in three-dimensional boundary layers, the simulation shows a drop in the Reynolds shear stress with increasing transverse strain. Also, the directions of the Reynolds shear stress vector and the mean velocity gradient vector were found to differ. In addition, the simulation shows a drop in the turbulent kinetic energy. The terms in the budget of the Reynolds stress were computed. In general, the individual terms do not show a significant change. However, an increase in dissipation rate rather than reduction in total production appears to be responsible for the drop in turbulent kinetic energy.

Moin, P.↗