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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Direct Numerical Simulation Database of High-Speed Flow over Parameterized Curved Walls

This study presents a direct numerical simulation (DNS) database of high-speed turbulent boundary layers (TBLs) subject to pressure gradients due to parametrically varied backward-facing and forward-facing wall curvatures, with an inflow Mach number of 4.9 and a friction Reynolds number of [Formula: see text] immediately before the onset of wall curvature. The Mach and Reynolds numbers are significantly higher than those reported in the literature for the DNS of pressure-gradient TBLs. The flow conditions and baseline wall geometries are representative of experiments in the high-speed blowdown wind tunnel at the National Aerothermochemistry Laboratory at Texas A&M University. The wall steepness of the baseline geometry for both the backward-facing and forward-facing walls was systematically varied to cause attached, incipiently separated, and fully separated flows. Precomputed flow statistics, including turbulent kinetic energy budgets, are available on the website of the Turbulence Modeling Resource of the NASA Langley Research Center, allowing other investigators to query any property of interest.

Engineering↗
DOE OSTI2018

Numerical Simulation of Ion Transport in a Nano-Electrospray Ion Source at Atmospheric Pressure

Understanding ion transport properties from the ion source to the mass spectrometer (MS) is essential for optimizing device performance. Numerical simulation helps in understanding of ion transport properties and, furthermore, facilitates instrument design. In contrast to previously reported numerical studies, ion transport simulations in a continuous injection mode whilst considering realistic space-charge effects have been carried out. The flow field was solved using Reynolds-averaged Navier-Stokes (RANS) equations, and a particle-in-cell (PIC) method was applied to solve a time-dependent electric field with local charge density. A series of ion transport simulations were carried out at different cone gas flow rates, ion source currents, and capillary voltages. A force evaluation analysis reveals that the electric force, the drag force, and the Brownian force are the three dominant forces acting on the ions. Both the experimental and simulation results indicate that cone gas flow rates of ≤250 slph (standard liter per hour) are important for high ion transmission efficiency, as higher cone gas flow rates reduce the ion signal significantly. The simulation results also show that the ion transmission efficiency reduces exponentially with an increased ion source current. Additionally, the ion loss due to space-charge effects has been found to be predominant at a higher ion source current, a lower capillary voltage, and a stronger cone gas counterflow. The interaction of the ion driving force, ion opposing force, and ion dispersion is discussed to illustrate ion transport mechanism in the ion source at atmospheric pressure. .

42 ENGINEERING↗
DOE OSTI2015

Development of numerical simulation system for thermal-hydraulic analysis in fuel assembly of sodium-cooled fast reactor

A numerical simulation system, which consists of a deformation analysis program and three kinds of thermal-hydraulics analysis programs, is being developed in Japan Atomic Energy Agency in order to offer methodologies to clarify thermal-hydraulic phenomena in fuel assemblies of sodium-cooled fast reactors under various operating conditions. This paper gives the outline of the system and its applications to fuel assembly analyses as a validation study.

42 ENGINEERING↗
DOE OSTI2020

Numerical simulation of an atmospheric pressure plasma jet with coaxial shielding gas

A two-dimensional ( r, z ) numerical simulation of the discharge characteristics of an atmospheric pressure plasma jet (APPJ), with coaxial shielding gas, was performed. The helium working gas flowed in a central capillary tube, engulfing a needle electrode powered by 13.7 MHz radio frequency sinusoidal voltage. The N 2 shielding gas flowed in the annular space of a coaxial tube. These gases emerged, in laminar flow, in a 78%N 2 -21%O 2 -1%Ar dry air ambient. The characteristics of the APPJ with shielding gas were compared to those of the APPJ without shielding gas. The nitrogen shielding gas hindered the diffusion of oxygen and argon from the ambient air into the helium jet. With the shielding gas present, more nitrogen penetrated into the helium core, causing a shorter plasma ‘plume’. The flow rates of the working and shielding gas, critically affected the gas temperature, and in turn the discharge characteristics. For a He flow of 2 standard liters per minute (slm), switching on the nitrogen shielding gas flow (at 4.5 slm) reduced the on-axis O 2 and Ar mole fractions from 3.9 × 10 - 4 to 6.8 × 10 - 5 and from 1.9 × 10 - 5 to 3.3 × 10 - 6 , respectively, at an axial distance of 3 mm downstream of the nozzle. The radial profiles of the mole fractions of the ambient gases were monotonically and strongly decreasing towards the system axis, for short axial distances from the nozzle (~1 mm), but became progressively flatter at longer distances from the nozzle (3 mm and 5 mm). Simulation predictions captured the salient features of experimental data of ambient species mole fractions in the plasma jet, and the 706 nm optical emission intensity profiles of the He 3 3 S excited state.

Physics↗
DOE OSTI2015

Numerical simulations of stable explosive-emission center plasma expansion in vacuum

We present the results of our numerical simulations of the stable cathode-flare expansion stage in the vacuum-breakdown switching period. We used the model of two-fluid hydrodynamics with corrections that allowed it to be applied to a consistent description of the electron and ion transport in plasma and the electron transport in a plasma-anode vacuum gap without imposing any joining conditions at the expanding plasma boundary. We have analyzed the dynamics of the density, directed velocity, and temperature profiles for both components and the influence of various factors on the plasma-cloud expansion velocity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Direct Numerical Simulation of hydrogen combustion at auto-ignitive conditions: Ignition, stability and turbulent reaction-front velocity

Direct Numerical Simulations (DNS) are performed to investigate the process of spontaneous ignition of hydrogen flames at laminar, turbulent, adiabatic and non-adiabatic conditions. Mixtures of hydrogen and vitiated air at temperatures representing gas-turbine reheat combustion are considered. Adiabatic spontaneous ignition processes are investigated first, providing a quantitative characterization of stable and unstable flames. Results indicate that, in hydrogen reheat combustion, compressibility effects play a key role in flame stability and that unstable ignition and combustion are consistently encountered for reactant temperatures close to the mixture’s characteristic crossover temperature. Furthermore, it is also found that the characterization of the adiabatic processes is also valid in the presence of non-adiabaticity due to wall heat-loss. Finally, a quantitative characterization of the instantaneous fuel consumption rate within the reaction front is obtained and of its ability, at auto-ignitive conditions, to advance against the approaching turbulent flow of the reactants, for a range of different turbulence intensities, temperatures and pressure levels.

42 ENGINEERING↗
DOE OSTI2018

Direct numerical simulation of mechanical response in synthetic additively manufactured microstructures

Additive manufacturing (AM) processes for metals can yield as-built microstructures that vary significantly from their cast or wrought counterparts. These microstructural variations can in turn, have profound effects on the properties of a component. Here, a modeling methodology is presented to investigate microstructurally-influenced mechanical response in additively manufactured structures via direct numeral simulation. Three-dimensional, synthetic voxelized microstructures are generated by kinetic Monte Carlo (kMC) additive manufacturing process simulations performed at four scan speeds to create a thin-wall cylindrical geometry notionally constructed using a concentric-pathed directed energy deposition AM process. The kMC simulations utilize a steady-state molten pool geometry that is held constant throughout the study. Resultant microstructures are mapped onto a highly-refined conformal finite-element mesh of a part geometry. A grain-scale anisotropic crystal elasticity model is then used to represent the constitutive response of each grain. Here, the response of the structure subjected to relatively simple load conditions is studied in order to provide understanding of both the influence of AM processing on microstructure as well as the microstructure's influence on the macroscale mechanical response.

42 ENGINEERING↗
DOE OSTI2021

Numerical simulation of NSLS-II fast orbit feedback system

We discuss recent work on the construction of a numerical simulation of the fast orbit feedback system at NSLS-II. The simulation operates in the time domain and also includes the spatial domain, i.e. all beam position monitors (BPMs) and correctors. It can accept inputs of the real beam orbit and the measured orbit response matrix. It can also add errors to every stage of the calculation. We present methods to verify the simulation results by comparing the simulation results with experimental data collected at NSLS-II. The results are in very good agreement. The effects of errors from BPMs and correctors, as well as the amplitude of excitation, on the feedback performance, are also explored. This simulation can be used to predict the behavior and performance of the system for a future upgrade.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Numerical simulation to study the transient self focusing of laser beam in plasma

In this paper, we present the numerical simulation for the coupled system of equations governing the dynamics of laser and Ion Acoustic Wave (IAW) in a collisionless plasma, when the coupling between the waves is through ponderomotive non-linearity. The nonlinear evolution of the laser beam is studied when the pump laser is perturbed by a periodic perturbation. By changing the perturbation wave number, we have studied its effect on the nonlinear evolution pattern of laser beam. In order to have a physical insight into the nonlinear dynamics of laser beam evolution in time and space, we have studied the laser and IAW spectra containing spatial harmonics. The magnitude of these harmonics changes with time and leads to time dependent localization of laser beam in spatial domain. The nonlinear dynamics of this localization is investigated in detail by using simulation and a semi-analytical model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Gaussian variational ansatz in the problem of anomalous sea waves: Comparison with direct numerical simulation

The nonlinear dynamics of an obliquely oriented wave packet on a sea surface is analyzed analytically and numerically for various initial parameters of the packet in relation to the problem of the so-called rogue waves. Within the Gaussian variational ansatz applied to the corresponding (1+2)-dimensional hyperbolic nonlinear Schrödinger equation (NLSE), a simplified Lagrangian system of differential equations is derived that describes the evolution of the coefficients of the real and imaginary quadratic forms appearing in the Gaussian. This model provides a semi-quantitative description of the process of nonlinear spatiotemporal focusing, which is one of the most probable mechanisms of rogue wave formation in random wave fields. The system of equations is integrated in quadratures, which allows one to better understand the qualitative differences between linear and nonlinear focusing regimes of a wave packet. Predictions of the Gaussian model are compared with the results of direct numerical simulation of fully nonlinear long-crested waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Numerical Simulation of Systems of Shear Bands in Ductile Metal with Inclusions

These are slides for a presentation on numerical simulation of systems of shear bands in ductile metal with inclusions, performed at Los Alamos National Laboratory. The conclusions are the following: A shear band is much thinner (0.1 μm) than the spacing between inclusions (100μm). Therefore fully resolved simulation is not feasible; asymptotic analysis allows modeling of a shear band as a velocity discontinuity within a moderately sized cell; formation criterion (critical strain) is determined by rate-dependent viscoplastic flow rule; inclusions cause shear bands to form at smaller strain; under expansion, shear bands form at the same strain than under compression; this can be applied to crystal plasticity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Direct numerical simulation of an auto-ignitive turbulent flame in a stratified dimethyl-ether (DME)/air mixture

Direct numerical simulations (DNS) of a turbulent premixed flame in an auto-ignitive dimethyl-ether (DME)/air mixture exhibiting two-staged ignition were conducted at elevated pressure and temperature using a statistically stationary planar configuration. Three sets of conditions with an identical turbulence intensity, but different stratification/turbulence length scales and correlations between temperature (T) and equivalence ratio (f) fields were simulated to study the dynamics of turbulent premixed flame propagation at auto-ignitive conditions. The influence of different stratification length scales and T-f correlations on turbulent flame speed was examined by comparing the flame statistics obtained from 2D simulations against laminar 1D solution. A global analysis of the front propagation speed was carried out for determining the effect of turbulent wrinkling and the role of upstream auto-ignition on the overall burning rate.

Desai, Swapnil S.↗
DOE OSTI2020

Direct numerical simulations of turbulent reacting flows with shock waves and stiff chemistry using many-core/GPU acceleration

Compressible reacting flows may display sharp spatial variation related to shocks, contact discontinuities or reactive zones embedded within relatively smooth regions. The presence of such phenomena emphasizes the relevance of shock-capturing schemes such as the weighted essentially non-oscillatory (WENO) scheme as an essential ingredient of the numerical solver. However, these schemes are complex and have more computational cost than the simple high-order compact or non-compact schemes. In this paper, we present the implementation of a seventh-order, minimally-dissipative mapped WENO (WENO7M) scheme in a newly developed direct numerical simulation (DNS) code called KAUST Adaptive Reactive Flows Solver (KARFS). In order to make efficient use of the computer resources and reduce the solution time, without compromising the resolution requirement, the WENO routines are accelerated via graphics processing unit (GPU) computation. The performance characteristics and scalability of the code are studied using different grid sizes and block decomposition. Furthermore, the performance portability of KARFS is demonstrated on a variety of architectures including NVIDIA Tesla P100 GPUs and NVIDIA Kepler K20X GPUs. In addition, the capability and potential of the newly implemented WENO7M scheme in KARFS to perform DNS of compressible flows is also demonstrated with model problems involving shocks, isotropic turbulence, detonations and flame propagation into a stratified mixture with complex chemical kinetics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A Numerical Simulation of a Single Shock-Accelerated Particle

Particle drag models, which capture macro viscous and pressure effects, have been developed over the years for various flow regimes to enable cost effective simulations of particle-laden flows. The relatively recent derivation by Maxey and Riley has provided an exact equation of motion for spherical particles in a flow field based on the continuum assumption. Many models that have been simplified from these equations have provided reasonable approximations; however, the sensitivity of particle-laden flows to particle drag requires a very accurate model to simulate. To develop such a model, a 2D axisymmetric Navier-Stokes direct numerical simulation of a single particle in a transient, shock-driven flow field was conducted in the hydrocode FLAG. FLAG’s capability to run arbitrary Lagrangian-Eulerian hydrodynamics coupled with solid mechanic models makes it an ideal code to capture the physics of the flow field around and in the particle as it is shock-accelerated – a challenging regime to study. The goal of this work is twofold: to provide a validation for FLAG’s Navier-Stokes and heat diffusion solutions, and to provide a rationale for recent experimental particle drag measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Direct numerical simulations in solid mechanics for quantifying the macroscale effects of microstructure and material model-form error

Two fundamental approximations in macroscale solid-mechanics modeling are (1) the assumption of scale separation in homogenization theory and (2) the use of a macroscopic plasticity material model that represents, in a mean sense, the multitude of inelastic processes occurring at the microscale. With the goal of quantifying the errors induced by these approximations on engineering quantities of interest, we perform a set of direct numerical simulations (DNS) in which polycrystalline microstructures are embedded throughout a macroscale structure. The largest simulations model over 50,000 grains. The microstructure is idealized using a randomly close-packed Voronoi tessellation in which each polyhedral Voronoi cell represents a grain. An face centered cubic crystal-plasticity model is used to model the mechanical response of each grain. The overall grain structure is equiaxed, and each grain is randomly oriented with no overall texture. The detailed results from the DNS simulations are compared to results obtained from conventional macroscale simulations that use homogeneous isotropic plasticity models. The macroscale plasticity models are calibrated using a representative volume element of the idealized microstructure. Furthermore, we envision that DNS modeling will be used to gain new insights into the mechanics of material deformation and failure.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Edge flame structure in a turbulent lifted flame: A direct numerical simulation study

This paper presents here a statistical analysis of edge flames in a turbulent lifted flame using direct numerical simulation (DNS). To investigate the dynamics of edge flames, a theoretical framework describing the edge-flame propagation velocity as a function of propagation velocities of mixture-fraction and product-mass fraction iso-surfaces at the flame base is used. The correlations between these propagation velocities and several other variables are then studied, including iso-surface curvatures, iso-surface orientations, strain rates, scalar dissipation rate and gradients of product mass fraction. The contribution of these parameters to the overall behaviour of the edge flame is also investigated using conditional averaging on two-dimensional spatial locations at the flame base. The analysis reveals that the tangential and normal strain rates in addition to the curvatures and scalar dissipation rates have significant contributions to the overall behaviour of the edge flame. The elliptical motion of the flame base described in our earlier study is extended to provide a clearer picture of how these various parameters affect the large fluctuations of edge-flame velocity observed at the flame base.

42 ENGINEERING↗
DOE OSTI2024

Direct numerical simulations of three-component Rayleigh–Taylor mixing and an improved model for multicomponent reacting mixtures

We present direct numerical simulations of a three-layer Rayleigh–Taylor instability (RTI) problem with a configuration based on the experiments of Suchandra & Ranjan ( J. Fluid Mech. , vol. 974, 2023, A35) and Jacobs & Dalziel ( J. Fluid Mech. , vol. 542, 2005, pp. 251–279). The problem consists of a layer of light fluid between two layers of heavy fluid with an Atwood number of 0.3. These simulations are first validated through comparison with available experimental data. The validated simulations are then utilized to analyse statistics in this three-component flow. First, length scales are examined utilizing spectra and two-point spatial correlations of velocity and species concentration fluctuations. Next, joint probability density functions (p.d.f.s) of species concentration are compared against several model p.d.f.s representing generalizations of the bivariate beta distribution. Notably, the joint p.d.f.s do not appear to be accurately described by a Dirichlet distribution, indicating the marginal distributions do not conform to a beta distribution. Finally, similarity of the present configuration to three-component mixing found in inertial confinement fusion (ICF) applications is exploited to develop and validate an improved model for the impact of multicomponent mixing on thermonuclear (TN) reaction rates. A single time instant from the present simulations is chosen for a TN burn calculation under the hypothetical assumption of ICF materials and temperatures. Total TN output from this second calculation is then compared against the prediction of the improved model. The new model is found to accurately predict TN reaction rates in both premixed and non-premixed configurations.

42 ENGINEERING↗