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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Numerical simulation of the reactive flow in advanced (HSR) combustors using KIVA-2

Recent work has been done with the goal of establishing ultralow emission aircraft gas turbine combustors. A significant portion of the effort is the development of three dimensional computational combustor models. The KIVA-II computer code which is based on the Implicit Continuous Eulerian Difference mesh Arbitrary Lagrangian Eulerian (ICED-ALE) numerical scheme is one of the codes selected by NASA to achieve these goals. This report involves a simulation of jet injection through slanted slots within the Rich burn/Quick quench/Lean burn (RQL) baseline experimental rig. The RQL combustor distinguishes three regions of combustion. This work specifically focuses on modeling the quick quench mixer region in which secondary injection air is introduced radially through 12 equally spaced slots around the mixer circumference. Steady state solutions are achieved with modifications to the KIVA-II program. Work currently underway will evaluate thermal mixing as a function of injection air velocity and angle of inclination of the slots.

Winowich, Nicholas S.↗
NASA NTRS1998

Numerical Simulation of Combustion and Extinction of a Solid Cylinder in Low-Speed Cross Flow

The combustion and extinction behavior of a diffusion flame around a solid fuel cylinder (PMMA) in low-speed forced flow in zero gravity was studied numerically using a quasi-steady gas phase model. This model includes two-dimensional continuity, full Navier Stokes' momentum, energy, and species equations with a one-step overall chemical reaction and second-order finite-rate Arrhenius kinetics. Surface radiation and Arrhenius pyrolysis kinetics are included on the solid fuel surface description and a parameter Phi, representing the percentage of gas-phase conductive heat flux going into the solid, is introduced into the interfacial energy balance boundary condition to complete the description for the quasi-steady gas-phase system. The model was solved numerically using a body-fitted coordinate transformation and the SIMPLE algorithm. The effects of varying freestream velocity and Phi were studied. These parameters have a significant effect on the flame structure and extinction limits. Two flame modes were identified: envelope flame and wake flame. Two kinds of flammability limits were found: quenching at low-flow speeds due to radiative loss and blow-off at high flow speeds due to insufficient gas residence time. A flammability map was constructed showing the existence of maximum Phi above which the solid is not flammable at any freestream velocity.

Tien, J. S.↗
NASA NTRS2001

Numerical Simulation of Two Phase Flows

Two phase flows can be found in broad situations in nature, biology, and industry devices and can involve diverse and complex mechanisms. While the physical models may be specific for certain situations, the mathematical formulation and numerical treatment for solving the governing equations can be general. Hence, we will require information concerning each individual phase as needed in a single phase. but also the interactions between them. These interaction terms, however, pose additional numerical challenges because they are beyond the basis that we use to construct modern numerical schemes, namely the hyperbolicity of equations. Moreover, due to disparate differences in time scales, fluid compressibility and nonlinearity become acute, further complicating the numerical procedures. In this paper, we will show the ideas and procedure how the AUSM-family schemes are extended for solving two phase flows problems. Specifically, both phases are assumed in thermodynamic equilibrium, namely, the time scales involved in phase interactions are extremely short in comparison with those in fluid speeds and pressure fluctuations. Details of the numerical formulation and issues involved are discussed and the effectiveness of the method are demonstrated for several industrial examples.

Liou, Meng-Sing↗
NASA NTRS1975

Numerical simulation of turbulent jet noise, part 1

Flow characteristics, such as quadrupole moments are examined in order to study generation of aerodynamic noise. The mean flow quantities are set in accordance with experimental data and the incompressible Navier-Stokes are solved numerically. Isolated downstream sections of a turbulent jet are modelled separately with the mean flow characteristics held constant in time. The flows are allowed to evolve until the fluctuating velocity components reach a statistically steady state. Cross section contour plots of the velocity components and the quadrupole moments at three different downstream positions are presented.

Metcalfe, R. W.↗
NASA NTRS1985

Numerical simulation of a turbulent flame stabilized behind a rearward-facing step

Flow of combustible mixtures in a plane channel past a smooth contraction followed by an abrupt expansion, in a typical dump combustor configuration, is modeled by a two-dimensional numerical technique based on the random vortex method. Both the inert and the reacting case are considered. In the latter, the flame is treated as an interface, self-advancing at a prescribed normal burning speed, while the dynamic effects of expansion due to the exothermicity of combustion are expressed by volumetric source lines delineated by its front. Solutions are shown to be in satisfactory agreement with experimental results, especially with respect to global properties such as the average velocity profiles and the reattachment length. The stochastic turbulent velocity components manifest interesting differences, especially near the walls where three-dimensional effects of turbulence are expected to be of importance.

Hsiao, C. C.↗
NASA NTRS1998

Numerical Simulation of THM Growth of CdTe in Presence of Rotating Magnetic Fields

The influence of rotating magnetic fields (RMF) on the flow pattern and compositional uniformity in the solution zone of a traveling heater method (THM) system for growth of CdTe is numerically investigated. The analysis is conducted at the 10(exp -6) and 10(exp -1) g(sub 0) as representative of space and ground processing conditions. It is shown that under microgravity conditions, application of RMF can be used to overwhelm residual buoyancy-induced convection and to control the uniformity of solution-zone composition at the growth front without appreciable modification of the growth interface shape. At high gravity levels, RMF is found not to be able to completely dominate buoyancy-induced convection. In this regime, for the range of field strengths studied, RMF is found to result in a) complex flow structures in the solution zone, b) enhancement of compositional non-uniformities at the growth front, and c) increased convexity of the growth interface. A scaling analysis of convection in the solution zone is used to generate a non-dimensional map delineating the RMF- and gravity-dominated flow regimes.

Ghaddar, Chahid K.↗
NASA NTRS1988

Numerical simulation of the vortical flow over a round-edged double-delta wing

Simulations of three-dimensional vortical flows over a thin double-delta wing with an aspect ratio of 2.05 have been performed. Steady-state solutions to the unsteady incompressible Navier-Stokes equations are obtained using a new implicit upwind-relaxation finite-difference scheme. The method is second-order accurate spatially and naturally dissipative. Numerical results indicate that key features of both vortical interaction and vortex breakdown are successfully simulated. The computed lift coefficients and lateral trajectories of the vortical cores are in good agreement with the experimental data.

Hsu, C.-H.↗
DOE OSTI2015

Numerical Simulation of Solidification of Colloids Inside a Differentially Heated Cavity

Development of the solid–liquid interface, distribution of the particle concentration field, as well as the development of thermosolutal convection during solidification of colloidal suspensions in a differentially heated cavity are investigated. The numerical model is based on the one-fluid mixture approach combined with the single-domain enthalpy porosity model for phase change, and it is implemented in FLUENT software package. The linear dependence of the liquidus and solidus temperatures with the concentration of the nanoparticles was assumed. A colloidal suspension consisting of water and copper or alumina nanoparticles were considered. In the current investigation, the nanoparticle size selected was 5 and 2 nm. The suspension was solidified unidirectionally inside a square differentially heated cavity that was cooled from the left side. It was found that the solid–liquid interface changed its morphology from a planar shape to a dendritic one as the solidification process proceeds in time, due to the constitutional supercooling that resulted from the increased concentration of particles at the solid–liquid interface rejected from the crystalline phase. Initially, the flow consisted of two vortices rotating in opposite directions. However, at later times, only one counter clockwise rotating cell survived. Changing the material of the particle to alumina resulted in crystallized phase with a higher concentration of particles. If it is compared to that of the solid phase resulted from freezing the copper–water colloidal suspension. Furthermore, decreasing the segregation coefficient destabilizes the solid–liquid interface and increases the intensity of the convection cell with respect to that of the case of no particle rejection. At slow freezing rates, the resulting crystal phase consisted of lower particle content compared to the case of higher freezing rate.

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

Three-dimensional numerical simulation of near-surface flows over the Martian north polar cap

Measurements made by Viking Lander VL-2 (48 N) have shown that the near-surface wind and temperature regime on Mars displays striking similarities to terrestrial counterparts. The diurnal radiative cycle is responsible for establishment of a well-defined thermal circulation in which downslope (Katabatic) flows prevail during the nighttime hours and weak upslope (anabatic) conditions prevail during the daytime. Previous work has indicated that the slope flows are much like those found on Earth, particularly the Katabatic winds, which show striking similarities to drainage flows observed over Antarctica. The low-level wind regime appears to be an important factor in the scouring of the martian landscape. The north polar cap shows evidence of eolian features such as dunes, frost streaks, and grooves from Viking imagery. The direction of the prevailing wind can in cases be inferred from the eolian features. We examine the thermally induced flows that result from the radiative heating and cooling of the martian north polar region using a comprehensive three-dimensional atmospheric mesoscale numerical model. The same model has been used previously for simulation of Antarctic Katabatic winds. The model equations are written in terrain-following coordinates to allow for irregular terrain; prognostic equations include the flux forms of the horizontal momentum equations, temperature, continuity. A surface energy budget equation is also incorporated in which the surface temperature is determined. Explicit parameterization of both terrestrial (longwave) and solar (shortwave) radiation is included. Turbulent transfer of heat and momentum in the martian atmosphere is assumed to follow the similarity expressions in the surface boundary layer on Earth. The terrain heights for the martian north polar region have been obtained from the U.S. Geological Survey map and digitized onto a 57x57 grid with a spacing of 75 km. The resulting terrain map is shown in Fig. 1. The vertical grid consists of 15 levels (delta equals 0.998, 0.99, 0.98, 0.97, 0.96, 0.94, 0.92, 0.90, 0.85, 0.775, 0.70, 0.60, 0.50, 0.30, 0.10). The high resolution in the lower atmosphere is necessary to capture details of the boundary layer flows. The lowest level corresponds to a height of approximately 20 m above the ground, the seconds level 100 m.

Parish, Thomas R.↗
DOE OSTI2015

Investigation of the effect of aggregates' morphology on concrete creep properties by numerical simulations

Prestress losses due to creep of concrete is a matter of interest for long-term operations of nuclear power plants containment buildings. Experimental studies by Granger (1995) have shown that concretes with similar formulations have different creep behaviors. The aim of this paper is to numerically investigate the effect of size distribution and shape of elastic inclusions on the long-term creep of concrete. Several microstructures with prescribed size distribution and spherical or polyhedral shape of inclusions are generated. By using the 3D numerical homogenization procedure for viscoelastic microstructures proposed by Šmilauer and Bažant (2010), it is shown that the size distribution and shape of inclusions have no measurable influence on the overall creep behavior. Moreover, a mean-field estimate provides close predictions. An Interfacial Transition Zone was introduced according to the model of Nadeau (2003). It is shown that this feature of concrete's microstructure can explain differences between creep behaviors.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Numerical simulations of chaotic and complex spatiotemporal patterns in fractional reaction–diffusion systems

The generalized fractional reaction–diffusion equations which exist in the form of noninteger order partial differential equations have now found wide application for illustrating important and useful physical phenomena, such as subdiffusive and superdiffusive scenarios. The space fractional derivatives are defined in the Riesz sense on the intervals 0 < α < 1 and 1 < α ≤ 2. We propose robust numerical techniques, such as a spectral representation of the fractional Laplacian operator in conjunction with the exponential time differencing method, and present the equivalent relationship between the Riesz fractional derivative and fractional Laplacian operator. We apply these techniques to numerically solve a range of chaotic processes, such as the Chua’s equations, Rössler system, Lorenz and Lorenz-type systems. Simulation results revealed various complex and spatiotemporal chaos, spiral chaos, intermittent chaos and spots patterns in two-dimensional space fractional reaction–diffusion problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Numerical Simulation and Experimental Investigation of Nitrogen Transfer Mechanism from Gas to Liquid Steel During Pressurized Electroslag Remelting Process

Since the gas–slag–metal three-phase transfer mechanism of nitrogen is not suitable to account for the nitrogen pickup during pressurized electroslag remelting process, the laboratory experiments have been carried out to study the transfer mechanism of nitrogen using different nitrogen partial pressures and electrode immersion depths. Then, two possible transfer mechanisms of nitrogen have been proposed, and a 2D transient numerical model has been developed to investigate the behavior of nitrogen based on the new proposed mechanisms. The simulated results are compared with the experiment to validate the feasibility of new proposed mechanism. The results show that the nitrogen content in ingot increases with the increasing nitrogen partial pressure and decreasing electrode immersion depth. At the current experimental condition, the electrode immersion depth seems to play a more important role in the nitrogen pickup than the nitrogen partial pressure. Due to the limits of thermodynamic and kinetic conditions, the nitrogen pickup is negligible via the transfer from gas to liquid metal through molten slag whether nitrogen chemically or physically dissolves into slag. The mechanism that nitrogen directly reacts with the partially exposed liquid metal film under the electrode tip due to the fluctuation of gas/slag interface could reasonably account for the experimental results. Furthermore, based on the mechanism, the predicted variation tendency of nitrogen content in ingot is in agreement with the experiment.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Numerical simulation of free shear flows: Towards a predictive computational aeroacoustics capability

Implicit and explicit spatial differencing techniques with fourth order accuracy have been developed. The implicit technique is based on the Pade compact scheme. A Dispersion Relation Preserving concept has been incorporated into both of the numerical schemes. Two dimensional Euler computation of a spatially-developing free shear flow, with and without external excitation, has been performed to demonstrate the capability of numerical schemes developed. Results are in good agreement with theory and experimental observation regarding the growth rate of fluctuating velocity, the convective velocity, and the vortex-pairing process.

Soh, W. Y.↗
NASA NTRS1999

Numerical Simulation of Complex Turbomachinery Flows

An unsteady, multiblock, Reynolds Averaged Navier Stokes solver based on Runge-Kutta scheme and Pseudo-time step for turbo-machinery applications was developed. The code was validated and assessed against analytical and experimental data. It was used to study a variety of physical mechanisms of unsteady, three-dimensional, turbulent, transitional, and cooling flows in compressors and turbines. Flow over a cylinder has been used to study effects of numerical aspects on accuracy of prediction of wake decay and transition, and to modify K-epsilon models. The following simulations have been performed: (a) Unsteady flow in a compressor cascade: Three low Reynolds number turbulence models have been assessed and data compared with Euler/boundary layer predictions. Major flow features associated with wake induced transition were predicted and studied; (b) Nozzle wake-rotor interaction in a turbine: Results compared to LDV data in design and off-design conditions, and cause and effect of unsteady flow in turbine rotors were analyzed; (c) Flow in the low-pressure turbine: Assessed capability of the code to predict transitional, attached and separated flows at a wide range of low Reynolds numbers and inlet freestream turbulence intensity. Several turbulence and transition models have been employed and comparisons made to experiments; (d) leading edge film cooling at compound angle: Comparisons were made with experiments, and the flow physics of the associated vortical structures were studied; and (e) Tip leakage flow in a turbine. The physics of the secondary flow in a rotor was studied and sources of loss identified.

Chernobrovkin, A. A.↗
DOE OSTI2019

Numerical Simulation and Validation of Gas and Molten Metal Flows in Close-Coupled Gas Atomization

Molten metal atomization in close-coupled gas atomization dies can operate between two limiting conditions, jetting and filming, together with several complex mechanisms: liquid-gas drafting, downward/upward shearing, melt bouncing, etc. Liquid jet deformation depends on flow and geometric parameters, such as liquid Reynolds, liquid Weber, and gas Mach numbers, as well as gas jet apex angle and melt tube tip extension and aspect ratio, among others. Understanding their effect is of importance for the metal powder making industry. Numerical gas atomization studies can provide approximated flow information and consider a wide range of conditions, beyond experimental reach. Here, 3D high-resolution simulations employing a 5-equation compressible flow model coupled with the volume-of-fluid method are compared with experiments, for liquid Weber number in the range of 1–30 and liquid Reynolds number below 10,000. This validation explores the predicting capabilities of the numerical model.

Hernandez, F.↗
NASA NTRS1993

Evolution of planetesimals. I - Dynamics: Relaxation in a thin disk. II - Numerical simulations

The study examines the effects of density inhomogeneity and differential rotation as well as inelastic collisions on the dynamical evolution of planetesimals. Consideration is given to a three-step analysis: the dynamical evolution of the planetesimals, collisions and mass accumulation, and interaction with gas. It is shown that the velocity dispersion of a cold system of planetesimals increases rapidly due to elastic gravitational scattering. When the dispersion in the epicycle amplitude becomes comparable to the planetesimals' Roche radius, energy is transferred from the systematic Keplerian shear to the dispersive motion. With a numerical N-body scheme, gravitational scattering and physical collisions among a system of planetesimals is simulated. It is shown that dynamical equilibrium is attained with a velocity dispersion comparable to the surface escape velocity of those planetesimals which contribute most of the system mass.

Palmer, P. L.↗
DOE OSTI2015

NUMERICAL SIMULATION OF FAST-MODE MAGNETOSONIC WAVES EXCITED BY PLASMOID EJECTIONS IN THE SOLAR CORONA

The Atmospheric Imaging Assembly instrument on board the Solar Dynamics Observatory has directly imaged the fast-propagating magnetosonic waves (FMWs) successively propagating outward along coronal magnetic funnels. In this study we perform a numerical investigation of the excitation of FMWs in the interchange reconnection scenario, with footpoint shearing flow being used to energize the system and drive the reconnection. The modeling results show that as a result of magnetic reconnection, the plasma in the current sheet is heated up by Joule dissipation to ∼10 MK and is ejected rapidly, developing the hot outflows. Meanwhile, the current sheet is torn into plasmoids, which are shot quickly both upward and downward. When the plasmoids reach the outflow regions, they impact and collide with the ambient magnetic field there, which consecutively launches FMWs. The FMWs propagate outward divergently away from the impact regions, with a phase speed of the Alfvén speed of ∼1000 km s{sup –1}. In the k – ω diagram of the Fourier wave power, the FMWs display a broad frequency distribution with a straight ridge that represents the dispersion relation. With the WKB approximation, at the distance of 15 Mm from the wave source region, we estimate the energy flux of FMWs to be E ∼ 7.0 × 10{sup 6} erg cm{sup –2} s{sup –1}, which is ∼50 times smaller than the energy flux related to the tube-channeled reconnection outflow. These simulation results indicate that energetically and dynamically the outflow is far more important than the waves.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1992

Assessment of modern methods in numerical simulations of high speed flows

Results of extensive studies on CFD algorithms for 2D inviscid flows in Cartesian and body fitted coordinates geometries are reviewed. These studies represent part of an ongoing investigation of combustion instabilities involving the interactions of high-speed nonlinear acoustic waves. Four numerical methods for the treatment of high speed flows are compared, namely, Roe-Sweby TVD, Yee symmetric TVD; Osher-Chakravarthy TVD; and the Colella's multi-dimensional Godunov method.

Pindera, M. Z.↗