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.

SourceResearch connectionFrom yearThrough year

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

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Particles climbing along a vertically vibrating tube: numerical simulation using the Discrete Element Method (DEM)

It has been reported experimentally that granular particles can climb along a vertically vibrating tube partially inserted inside a granular silo. Here, we use the Discrete Element Method (DEM) available in the Multiphase Flow with Interphase eXchanges (MFIX) code to investigate this phenomenon. By tracking the movement of individual particles, the climbing mechanism was illustrated and analyzed. The numerical results show that a sufficiently high vibration strength is needed to form a low solids volume fraction region inside the lower end of the vibrating tube, a dense region in the middle of the tube, and to bring the particles outside from the top layers down to fill in the void. The results also show that particle compaction in the middle section of the tube is the main cause of the climbing. Consequently, varying parameters which influence the compacted region, such as the restitution coefficient, change the climbing height.

42 ENGINEERING↗

NASA NTRS1990

Numerical simulation of solidification in a horizontal cylindrical annulus charged with an aqueous salt solution

A computational study of solidification of a binary Na2CO3 solution in a horizontal cylindrical annulus is performed using a continuum formulation with a control-volume based, finite-difference scheme. The initial conditions were selected to facilitate the study of counter thermal and solutal convection, accompanied by extensive mushy region growth. Numerical results are compared with experimental data with mixed success. Qualitative agreement is obtained for the overall solidification process and associated physical phenomena. However, the plume thickness calculated for the solutally-driven convective upflow is substantially smaller than the observed value. Evolution of double-diffusive layers is predicted, but over a time scale much smaller than that observed experimentally. Good agreement is obtained between predicted and measured results for solid growth, but the mushy region thickness is significantly overpredicted.

Neilson, D. G.↗

DOE OSTI2019

Hot-spot induced mild ignition: Numerical simulation and scaling analysis

Mild ignition is a 'non-ideal' autoignition process observed at some conditions in rapid compression machine, shock tube and flow reactor experiments where flames and/or reaction fronts initiate due to non-uniformities within the test mixture. These flames/fronts can consume the entire mixture, or compression-heat the unburned gas, causing bulk ignition at times earlier than in homogeneous environments. This study numerically investigates mild ignition features using computational fluid dynamics techniques under a scenario where a small hot-spot is employed to initiate the process in a spherically-symmetric geometry. Only subsonic regimes are considered. The intent is to develop greater fundamental understandings of competing phenomena, and formulate a robust framework to quantify the perturbative behavior. Syngas blends (CO/H 2 = 80/20) are employed as the fuel, mixed with air at lean fuel loadings of phi = 0.2 and 0.5, with the initial temperatures and pressures covering 900-1120K, and 1.5-15 bar, respectively. Hot-spots ranging from 50 to 2,000 mu m in diameter are utilized with temperatures elevated from the surroundings by 1 to 90%, i.e., T hot /T 0 =1.01-1.90. Scaling arguments are derived to analyze competing processes of: (a) quenching vs. ignition, (b) extinction vs. propagation, and (c) flame consumption/forced autoignition vs. homogeneous ignition. Here, this work provides new insight into non-uniform autoignition phenomena, and discusses opportunities and challenges associated with mitigating it. In particular, scaling arguments indicate that experimental parameters could have potential to either suppress hot-spot initiated flame formation (e.g., via bath gases with high thermal diffusivity), or avert flame compression of the end gas (e.g., via lean/diluted mixtures, or bath gases with high heat capacity).

42 ENGINEERING↗

NASA NTRS1992

Numerical simulation of entry flow over blunt swept-wing planes

Viscous, heat-conducting flow with chemical and vibrational relaxation processes of the constituent gases surrounding winged spacecraft is considered in the continuum regime. The Navier-Stokes equations are appended by additional vibrational energy and species rate equations and supplemented by the equations of state and the phenomenological laws based on mixture rules or collisional cross sections. Numerical convective flux can be obtained from several forms of one-dimensional Riemann solver, with or without entropy correction. High-order accuracy is obtained from two types of reconstructive interpolation. A number of explicit and implicit numerical schemes have been implemented as a means to yield converged solutions. Both shock-fitting, finite-difference and shock-capturing, finite-volume techniques have been tested for configurations such as a sphere, double ellipsoid, blunt-edge delta wing, a European Hermes vehicle, and the U.S. Shuttle Orbiter. The shock-fitting code provides excellent results only for simple configurations, whereas the shock-capturing code leads to overall satisfying solutions for complex geometries.

Li, C. P.↗

NASA NTRS1988

Numerical simulation of hypersonic flow over an aeroassist flight experiment vehicle

Solutions of wind-tunnel and entry-flight flow around the vehicle are obtained from the Navier-Stokes equations coupled with the chemical species continuity equations if needed. The time-iterative method employs several techniques: shock fitting, chemistry-split ADI and an algebraic grid in conformal spherical-polar space. Sensitivities of the results to numerical parameters and to frozen, equilibrium and finite rate reactions are investigated in the forebody computation. Quantitative results are obtained for the shock layer and the near wake for the entire vehicle corresponding to both ground test and flight conditions. Complex flow characteristics are analyzed on the basis of the complete flowfield over the aerobrake and simplified afterbodies. The method is stable and cost effective, and has yielded shock locations and wall pressure distributions which are in good agreement with wind-tunnel data.

Li, C. P.↗

DOE OSTI2015

Numerical simulation of nonlinear processes in a beam-plasma system

In the present paper we consider the efficiency of the electromagnetic radiation generation due to various nonlinear processes in the beam-plasma system. The beam and plasma parameters were chosen close to the parameters in the experiment on the GOL-3 facility (BINP SB RAS). The model of the collisionless plasma is described by system of the Vlasov-Maxwell equations with periodic boundary conditions. The parallel numerical algorithm is based on the particles-in-cell method (PIC) with mixed Euler-Lagrangian domain decomposition. Various scenarios of nonlinear evolution in the beam-plasma system under the influence of an external magnetic field in case of a low density beam were studied. The energy transfer from one unstable mode to the others modes was observed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1993

Numerical simulation of rarefied flow through a slit. I - Direct simulation Monte Carlo results

The pressure-driven flow of a rarefied monatomic gas through a two-dimensional slit is simulated using the direct simulation Monte Carlo technique. Of particular interest is the change in flow field structure as pressure ratio and Knudsen number are varied. Comparisons are made to quantify the limits of validity of free-molecular theory and approximate, nearly free-molecular iterative methods. Also addressed is the sensitivity of the numerical solutions to grid structure and boundary conditions. The free-molecular theory is found to predict quantitative flow field properties (e.g., centerline velocities or downstream flux profiles) reasonably well for large finite Knudsen number with the error dependent on the pressure ratio. The nearly free-molecular corrections are shown to have limited range of applicability. A previously derived parameter is found to correlate total mass flux well as a function of pressure ratio and Knudsen number over a large portion of the transitional regime.

Wadsworth, D. C.↗

DOE OSTI2018

Numerical simulation of the geometrical-optics reduction of CE2 and comparisons to quasilinear dynamics

Zonal flows have been observed to appear spontaneously from turbulence in a number of physical settings. A complete theory for their behavior is still lacking. Recently, a number of studies have investigated the dynamics of zonal flows using quasilinear (QL) theories and the statistical framework of a second-order cumulant expansion (CE2). A geometrical-optics (GO) reduction of CE2, derived under an assumption of separation of scales between the fluctuations and the zonal flow, is studied here numerically. The reduced model, CE2-GO, has a similar phase-space mathematical structure to the traditional wave-kinetic equation, but that wave-kinetic equation has been shown to fail to preserve enstrophy conservation and to exhibit an ultraviolet catastrophe. CE2-GO, in contrast, preserves nonlinear conservation of both energy and enstrophy. We show here how to retain these conservation properties in a pseudospectral simulation of CE2-GO. We then present nonlinear simulations of CE2-GO and compare with direct simulations of quasilinear (QL) dynamics. Furthermore, we find that CE2-GO retains some similarities to QL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1994

Conditional statistics in a turbulent premixed flame derived from direct numerical simulation

The objective of this paper is to briefly introduce conditional moment closure (CMC) methods for premixed systems and to derive the transport equation for the conditional species mass fraction conditioned on the progress variable based on the enthalpy. Our statistical analysis will be based on the 3-D DNS database of Trouve and Poinsot available at the Center for Turbulence Research. The initial conditions and characteristics (turbulence, thermo-diffusive properties) as well as the numerical method utilized in the DNS of Trouve and Poinsot are presented, and some details concerning our statistical analysis are also given. From the analysis of DNS results, the effects of the position in the flame brush, of the Damkoehler and Lewis numbers on the conditional mean scalar dissipation, and conditional mean velocity are presented and discussed. Information concerning unconditional turbulent fluxes are also presented. The anomaly found in previous studies of counter-gradient diffusion for the turbulent flux of the progress variable is investigated.

Mantel, Thierry↗

NASA NTRS1988

Numerical simulation of axisymmetric turbulent flow in combustors and diffusors

A method for predicting turbulent flow in combustors and diffusers is developed. The Navier-Stokes equations, incorporating a turbulence kappa-epsilon model equation, were solved in a nonorthogonal curvilinear coordinate system. The solution applied the finite volume method to discretize the differential equations and utilized the SIMPLE algorithm iteratively to solve the differenced equations. A zonal grid method, wherein the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. However, treatment of the zonal boundaries required special handling. Boundary overlap and interpolating techniques were used and an adjustment of the flow variables was required to assure conservation of mass, momentum and energy fluxes. The numerical accuracy was assessed using different finite differencing methods, i.e., hybrid, quadratic upwind and skew upwind, to represent the convection terms. Flows in different geometries of combustors and diffusers were simulated and results compared with experimental data and good agreement was obtained.

Yung, Chain Nan↗

DOE OSTI2017

Numerical Simulations of Supernova Remnant Evolution in a Cloudy Interstellar Medium

The mixed morphology class of supernova remnants has centrally peaked X-ray emission along with a shell-like morphology in radio emission. White and Long proposed that these remnants are evolving in a cloudy medium wherein the clouds are evaporated via thermal conduction once being overrun by the expanding shock. Their analytical model made detailed predictions regarding temperature, density, and emission profiles as well as shock evolution. We present numerical hydrodynamical models in 2D and 3D including thermal conduction, testing the White and Long model and presenting results for the evolution and emission from remnants evolving in a cloudy medium. We find that, while certain general results of the White and Long model hold, such as the way the remnants expand and the flattening of the X-ray surface brightness distribution, in detail there are substantial differences. In particular we find that the X-ray luminosity is dominated by emission from shocked cloud gas early on, leading to a bright peak, which then declines and flattens as evaporation becomes more important. In addition, the effects of thermal conduction on the intercloud gas, which is not included in the White and Long model, are important and lead to further flattening of the X-ray brightness profile as well as lower X-ray emission temperatures.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1987

Numerical Simulation For Supersonic Inlets

Flows calculated for realistic engine-inlet conditions. Computer code LAPIN, large-perturbation inlet, developed to analyze large-perturbation, transient-flow fields in supersonic inlets. Robust, quick-running code capable of solving unsteady quasi-one-dimensional, inviscid-flow problems in mixed subsonic and supersonic regimes for inlets. Approach based upon quasi-one-dimensional, inviscid, unsteady formulation including engineering models of unstart/restart, bleed, bypass, and geometrical effects. Numerical solution of governing time-dependent equations of motion accomplished through shock-capturing, finite-difference algorithm. Program written in FORTRAN IV.

Varner, M. O.↗

NASA NTRS1989

Numerical simulation of three-dimensional transonic turbulent projectile aerodynamics by TVD schemes

The two-dimensional symmetric TVD scheme proposed by Yee has been extended to and investigated for three-dimensional thin-layer Navier-Stokes simulation of complex aerodynamic problems. An existing three-dimensional Navier-stokes code based on the beam and warming algorithm is modified to provide an option of using the TVD algorithm and the flow problem considered is a transonic turbulent flow past a projectile with sting at ten-degree angle of attack. Numerical experiments conducted for three flow cases, free-stream Mach numbers of 0.91, 0.96 and 1.20 show that the symmetric TVD algorithm can provide surface pressure distribution in excellent agreement with measured data; moreover, the rate of convergence to attain a steady state solution is about two times faster than the original beam and warming algorithm.

Shiau, Nae-Haur↗

NASA NTRS2010

Numerical Simulation of Vitiation Effects on a Hydrogen-Fueled Dual-Mode Scramjet

The Wind-US computational fluid dynamics (CFD) flow solver was used to simulate dual-mode direct-connect ramjet/scramjet engine flowpath tests conducted in the University of Virginia (UVa) Supersonic Combustion Facility (SCF). The objective was to develop a computational capability within Wind-US to aid current hypersonic research and provide insight to flow as well as chemistry details that are not resolved by instruments available. Computational results are compared with experimental data to validate the accuracy of the numerical modeling. These results include two fuel-off non-reacting and eight fuel-on reacting cases with different equivalence ratios, split between one set with a clean (non-vitiated) air supply and the other set with a vitiated air supply (12 percent H2O vapor). The Peters and Rogg hydrogen-air chemical kinetics model was selected for the scramjet simulations. A limited sensitivity study was done to investigate the choice of turbulence model and inviscid flux scheme and led to the selection of the k-epsilon model and Harten, Lax and van Leer (for contact waves) (HLLC) scheme for general use. Simulation results show reasonably good agreement with experimental data and the overall vitiation effects were captured.

Vyas, Manan A.↗

NASA NTRS1988

A two-dimensional numerical simulation of a supersonic, chemically reacting mixing layer

Research has been undertaken to achieve an improved understanding of physical phenomena present when a supersonic flow undergoes chemical reaction. A detailed understanding of supersonic reacting flows is necessary to successfully develop advanced propulsion systems now planned for use late in this century and beyond. In order to explore such flows, a study was begun to create appropriate physical models for describing supersonic combustion, and to develop accurate and efficient numerical techniques for solving the governing equations that result from these models. From this work, two computer programs were written to study reacting flows. Both programs were constructed to consider the multicomponent diffusion and convection of important chemical species, the finite rate reaction of these species, and the resulting interaction of the fluid mechanics and the chemistry. The first program employed a finite difference scheme for integrating the governing equations, whereas the second used a hybrid Chebyshev pseudospectral technique for improved accuracy.

Drummond, J. Philip↗

NASA NTRS2003

Numerical Simulation And Experimental Investigation Of The Lift-Off And Blowout Of Enclosed Laminar Flames

Enclosed flames are found in many industrial applications such as power plants, gas-turbine combustors and jet engine afterburners. A better understanding of the burner stability limits can lead to development of combustion systems that extend the lean and rich limits of combustor operations. This paper reports a fundamental study of the stability limits of co-flow laminar jet diffusion flames. A numerical study was conducted that used an adaptive mesh refinement scheme in the calculation. Experiments were conducted in two test rigs with two different fuels and diluted with three inert species. The numerical stability limits were compared with microgravity experimental data. Additional normal-gravity experimental results were also presented.

Venuturmilli, Rajasekhar↗

NASA NTRS2013

On the Theory and Numerical Simulation of Cohesive Crack Propagation with Application to Fiber-Reinforced Composites

The phenomenon of crack propagation is among the predominant modes of failure in many natural and engineering structures, often leading to severe loss of structural integrity and catastrophic failure. Thus, the ability to understand and a priori simulate the evolution of this failure mode has been one of the cornerstones of applied mechanics and structural engineering and is broadly referred to as "fracture mechanics." The work reported herein focuses on extending this understanding, in the context of through-thickness crack propagation in cohesive materials, through the development of a continuum-level multiscale numerical framework, which represents cracks as displacement discontinuities across a surface of zero measure. This report presents the relevant theory, mathematical framework, numerical modeling, and experimental investigations of through-thickness crack propagation in fiber-reinforced composites using the Variational Multiscale Cohesive Method (VMCM) developed by the authors.

Rudraraju, Siva Shankar↗

DOE OSTI2020

Numerical simulations of colliding jets in an external wind: application to 3C 75

ABSTRACT The radio galaxy 3C 75 is remarkable because it contains a pair of radio-loud active galaxies, each of which produces a two-sided jet, with the jet beams appearing to collide and merge to the west of the galaxies. Motivated by 3C 75, we have conducted three-dimensional hydrodynamic simulations of jet collisions. We have extended previous studies by modelling the physical properties of the cluster atmosphere, including an external wind, and using realistic jet powers obtained from observational data. We are able to produce a morphology similar to that of 3C 75. The simulations imply that direct contact between the bulk jet flows on the west of the source is required to produce a morphology consistent with 3C 75. We quantify how the merging jets decelerate, how the wind deflects the jets and cocoons, the entrainment of intra-cluster material into the cocoons, the cocoon energetics, and how the jet interactions generate enstrophy. By comparing simulations of pairs of two-sided jets with those of single two-sided sources, we determine how the interaction between two bipolar jets changes their evolution. The unprecedented sensitivity and angular resolution of upcoming observatories will lead to the detection of many more complex sources at high redshift, where interacting jets are expected to be more numerous. The morphology of these complex sources can provide significant insight into the conditions in their environments.

Musoke, G.↗