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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Numerical Simulations of Noisy Variational Quantum Eigensolver Ansatz Circuits

This is a case study of the variational quantum eigensolver (VQE) method using numerical simulations to test the influence of noise on the accuracy of the underlying circuit ansatz. We investigate a computational chemistry application of VQE to calculate the electronic ground state and its energy for Sodium Hydride (NaH), a prototypical two-electron problem. Using a one-parameter ansatz derived from unitary coupled cluster (UCC) theory, we simulate the effects of noise on the energy expectation value and variance with respect to the ansatz parameter. These numerical simulations provide insights into the accuracy of the prepared quantum state and the efficiency of the classical optimizer that iteratively refines the ansatz. We conduct a comparative study between analytical results derived for the UCC ansatz in the absence of noise and the noisy numerical simulation results obtained using an isotropic depolarizing noise model for each gate. We also compare the relative increase in noise on logically equivalent UCC ansatz circuits generated by randomized compiling. Notably, we observe that the intrinsic variance in the energy due to the simplicity of the ansatz itself compares with the noise induced by the bare circuit.

Gowrishankar, Meenambika↗
DOE OSTI2016

NUMERICAL SIMULATIONS OF KELVIN–HELMHOLTZ INSTABILITY: A TWO-DIMENSIONAL PARAMETRIC STUDY

Using two-dimensional simulations, we numerically explore the dependences of Kelvin–Helmholtz (KH) instability upon various physical parameters, including viscosity, the width of the sheared layer, flow speed, and magnetic field strength. In most cases, a multi-vortex phase exists between the initial growth phase and the final single-vortex phase. The parametric study shows that the evolutionary properties, such as phase duration and vortex dynamics, are generally sensitive to these parameters, except in certain regimes. An interesting result is that for supersonic flows, the phase durations and saturation of velocity growth approach constant values asymptotically as the sonic Mach number increases. We confirm that the linear coupling between magnetic field and KH modes is negligible if the magnetic field is weak enough. The morphological behavior suggests that the multi-vortex coalescence might be driven by the underlying wave–wave interaction. Based on these results, we present a preliminary discussion of several events observed in the solar corona. The numerical models need to be further improved to perform a practical diagnostic of the coronal plasma properties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Numerical Simulation Applications in the Design of EGS Collab Experiment 1

The United States Department of Energy, Geothermal Technologies Office (GTO) is funding a collaborative investigation of enhanced geothermal systems (EGS) processes at the meso-scale. This study, referred to as the EGS Collab project, is a unique opportunity for scientists and engineers to investigate the creation of fracture networks and circulation of fluids across those networks under in-situ stress conditions. The EGS Collab project is envisioned to comprise three experiments and the site for the first experiment is on the 4850 Level in phyllite of the Precambrian Poorman formation, at the Sanford Underground Research Facility, located at the former Homestake Gold Mine, in Lead, South Dakota. Principal objectives of the project are to develop a number of intermediate-scale field sites and to conduct well-controlled in situ experiments focused on rock fracture behavior and permeability enhancement. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO’s flagship EGS research effort. As noted, outcomes from the EGS Collab project experiments will serve as benchmarks for computer code verification, but numerical simulation additionally plays an essential role in designing these meso-scale experiments. This paper describes specific numerical simulations supporting the design of Experiment 1, a field test involving hydraulic stimulation of two fractures from notched sections of the injection borehole and fluid circulation between sub-horizontal injection and production boreholes in each fracture individually and collectively, including the circulation of chilled water. Whereas the mine drift allows for accurate and close placement of monitoring instrumentation to the developed fractures, active ventilation in the drift cooled the rock mass within the experimental volume. Numerical simulations were executed to predict seismic events and magnitudes during stimulation, initial fracture orientations for smooth horizontal wellbores, pressure requirements for fracture initiation from notched wellbores, fracture propagation during stimulation between the injection and production boreholes, tracer travel times between the injection and production boreholes, and produced fluid temperatures with chilled water injections, and pressure limits on fluid circulation to avoid fracture growth.

Numerical Simulation, EGS Collab, EGS, Sanford Und↗
DOE OSTI2018

Numerical Simulation Applications in the Design of EGS Collab Experiment 1

The United States Department of Energy, Geothermal Technologies Office (GTO) is funding a collaborative investigation of enhanced geothermal systems (EGS) processes at the meso-scale. This study, referred to as the EGS Collab project, is a unique opportunity for scientists and engineers to investigate the creation of fracture networks and circulation of fluids across those networks under in-situ stress conditions. The EGS Collab project is envisioned to comprise three experiments and the site for the first experiment is on the 4850 Level (4,850 feet below ground surface) in phyllite of the Precambrian Poorman formation, at the Sanford Underground Research Facility, located at the former Homestake Gold Mine, in Lead, South Dakota. Principal objectives of the project are to develop a number of intermediate-scale field sites and to conduct well-controlled in situ experiments focused on rock fracture behavior and permeability enhancement. Data generated during these experiments will be compared against predictions of a suite of computer codes specifically designed to solve problems involving coupled thermal, hydrological, geomechanical, and geochemical processes. Comparisons between experimental and numerical simulation results will provide code developers with direction for improvements and verification of process models, build confidence in the suite of available numerical tools, and ultimately identify critical future development needs for the geothermal modeling community. Moreover, conducting thorough comparisons of models, modelling approaches, measurement approaches and measured data, via the EGS Collab project, will serve to identify techniques that are most likely to succeed at the Frontier Observatory for Research in Geothermal Energy (FORGE), the GTO's flagship EGS research effort. As noted, outcomes from the EGS Collab project experiments will serve as benchmarks for computer code verification, but numerical simulation additionally plays an essential role in designing these meso-scale experiments. This paper describes specific numerical simulations supporting the design of Experiment 1, a field test involving hydraulic stimulation of two fractures from notched sections of the injection borehole and fluid circulation between sub-horizontal injection and production boreholes in each fracture individually and collectively, including the circulation of chilled water. Whereas the mine drift allows for accurate and close placement of monitoring instrumentation to the developed fractures, active ventilation in the drift cooled the rock mass within the experimental volume. Numerical simulations were executed to predict seismic events and magnitudes during stimulation, initial fracture orientations for smooth horizontal wellbores, pressure requirements for fracture initiation from notched wellbores, fracture propagation during stimulation between the injection and production boreholes, tracer travel times between the injection and production boreholes, produced fluid temperatures with chilled water injections, pressure limits on fluid circulation to avoid fracture growth, temperature environment surrounding the 4850 Level drift, and fracture propagation within a stress field altered by drift excavation, ventilation cooling, and dewatering.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Numerical Simulations of Non-Proliferation Experiments

We have implemented a set of numerical simulations to numerically investigate the outcomes of non-proliferation experiments. The high explosive (HE) in the simulations is in the form of a cylinder. The two-dimensional simulations are carried out in the r-z cylindrical coordinate, and the three-dimensional simulation is in the Cartesian coordinate. From the simulation, the cavity formed after detonation stops its expansion after about 0.15 sec. At the time the cavity stops expanding, the cavity shape seems quite spherical although the initial HE is in the shape of cylinder, but the radius in the vertical direction, which is about 19.7m, is a little less than the one in the horizonal direction, which is about 21.5m. Although there is no layer of air between the HE and tuff in real experiments, in addition to the simulations without air, we add a layer of air in one simulation to assess any possible influence of the air on the shape of the cavity.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
NASA NTRS1978

The direct numerical simulations of the turbulent wakes of axisymmetric bodies

Results of direct numerical simulations of turbulence are compared with both laboratory data and self-similarity theory for the case of the turbulent wakes of towed, axisymmetric bodies. In general, the agreement of the simulation results with both the laboratory data and the self-similarity theory is good, although the comparisons are hampered by inadequate procedures for initializing the numerical simulations.

Riley, J. J.↗
NASA NTRS2003

Numerical Simulations of Non-Newtonian Convection in Ice: Application to Europa

Numerical simulations of solid state convection in Europa's ice shell have so far been limited to consideration of Newtonian flow laws, where the viscosity of ice is strongly dependent upon temperature, predicting that a stagnant lid should form at the top (10-40%) of a convecting ice shell. Such large thicknesses seem to contradict estimates of the effective elastic thickness of Europa s ice shell during its geologically active period. Recent laboratory experiments characterize the rheology of ice as the sum of contributions from several temperature and strain rate-dependent creep mechanisms. We present the results of numerical simulations of convection within Europa s ice shell using the finite-element model Citcom, applying the non-Newtonian rheology of grain boundry sliding. Our calculations suggest a shallower brittle/ductile transition and larger interior convective velocities compared to Newtonian rheology. The flow field is time-dependent, with small, localized upwellings and downwellings at the thermal boundary layers that have minimal topographic expression at the surface.

Barr, A. C.↗
NASA NTRS1977

Numerical simulation of spacecraft charging phenomena

A numerical simulation program is being constructed having the following features: (1) infinite circular cylindrical geometry with angle-dependence, (2) inclusion of incident particles, photoelectrons, secondary electrons, backscattered electrons, any gun emissions, and any internal current pathways including surface conductive layers, (3) quasistatic time-dependent iteration, in which sheath potential changes during particle transit times are ignored, (4) use of approximate, locally-dependent space charge density expressions in solving Poisson's equation for sheath potentials, with use of numerical orbit-following to determine surface currents, and (5) incident particle velocity distributions isotropic or beam-like, or some superposition of these. Rationales for each of these features are discussed.

Laframboise, J. G.↗
NASA NTRS1981

Numerical simulation of transonic flows in diffusers

Numerical simulations were made of two-dimensional transonic flows in diffusers, including flow separation induced by a shock or adverse pressure gradient. The mass-averaged, time-dependent, compressible Navier-Stokes equations, simplified by the thin-layer approximation, were solved using MacCormack's hybrid method. The eddy-viscosity formulation was described by the Wilcox-Rubesin's two-equation, k-omega model. Detailed comparison of the computed results with measurements showed good agreement in all cases, including one with massive separation induced by a strong shock. The computation correctly predicted the details of a distinct lambda shock pattern, closely duplicating the configuration observed experimentally in spark-schlieren photographs.

Liou, M.-S.↗
NASA NTRS1983

Numerical simulation of turbulent trailing edge flows

Numerical simulations of the time-dependent, Reynolds-averaged, Navier-Stokes equations, employing a two-equation turbulence model, are presented and compared with measurements from a series of trailing edge experiments at transonic Mach numbers. The test flows include an asymmetric flow with no separation, an asymmetric flow with a small region of separation and a symmetric flow with a large shock-wave induced separated zone. Comparisons are made for mean surface quantities as well as for mean and fluctuating flow-field quantities. For the trailing-edge flows with little or no separation, the solutions correctly predict all the major features of the flow field. Treatment of the viscous-inviscid interaction was found to be important for predicting these test cases. Two-equation eddy-viscosity turbulence models were found to be adequate for these flows. However, for the shock-wave induced separation case, these turbulence models were inadequate to predict this flow field. Modifications of the turbulence model to correct these deficiencies are discussed.

Horstman, C. C.↗
NASA NTRS1987

Numerical simulation of viscous transonic airfoil flows

Numerical simulations of transonic airfoil flows using the Reynolds-averaged Navier-Stokes equations and various turbulence models are presented and compared with experimental data. Three different airfoils were investigated under varying flow conditions ranging from subcritical unseparated flows to supercritical separated flows. The turbulence models investigated consisted of three zero-equation models and one two-equation model. For unseparated flows involving weak viscous-inviscid interactions, the four models were comparable in their agreement with experiment. For separated flows involving strong viscous-inviscid interactions, the nonequilibrium zero-equation model of Johnson and King gave the best overall agreement with experiment.

Coakley, Thomas J.↗
NASA NTRS1987

Numerical simulations of collisions between rotating particles

Numerical simulations of Keplerian systems consisting of 200 mutually colliding rotating particles show that friction and surface irregularity reduce the equilibrium velocity dispersion and transfer some of the energy of random velocities to rotational velocities. The results confirm the theoretical predictions of Salo (1987). Simulations are performed for identical and different particle sizes, and for a power-law distribution of sizes. Taking rotation into account is not found to alter the previously observed Rayleigh distribution of eccentricities and inclinations, while the components of the rotational velocities are found to follow a Gaussian distribution. Application to the rarefied regions of Saturn's rings suggests that friction is able to reduce the equilibrium geometric thickness by about one-half.

Salo, H.↗
NASA NTRS1987

Collisional interchange instability. I - Numerical simulations of intermediate-scale irregularities

Numerical simulations of the generalized Rayleigh-Taylor instability are presented. The model and simulations are applicable to bottomside and topside spread F, unstable barium cloud dynamics, and collisional interchange instability in general. The principal result is that the evolution of the instability tends to an anisotropic state consisting of nearly sinusoidal (quasi-periodic) variation along the effective electric field, and shocklike structures propagating perpendicular to E(eff) along the extrema of the quasi-periodic structures. The spectral properties of the nonlinear state are analyzed using one-dimensional power spectra calculated along spatial trajectories for selected angles to E(eff). In this way a direct comparison to in situ probe data can be made. The inherent anisotropy of the nonlinear state is reflected in major qualitative differences between the spectra taken parallel to and perpendicular to E(eff). The fundamental finding of the present work is that anisotropy in interchange dynamics is much greater than had been previously reported. This strong anisotropy can explain much of the spectral and spatial structural characteristics of both bottomside and topside spread F. In a companion paper a comparison of the simulation results to various in situ data sets is given.

Zargham, S.↗
NASA NTRS1989

Numerical simulations in combustion

This paper reviews numerical simulations in reacting flows in general and combustion phenomena in particular. It is shown that use of implicit schemes and/or adaptive mesh strategies can improve convergence, stability, and accuracy of the solution. Difficulties increase as turbulence and multidimensions are considered, particularly when finite-rate chemistry governs the given combustion problem. Particular attention is given to the areas of solid-propellant combustion dynamics, turbulent diffusion flames, and spray droplet vaporization.

Chung, T. J.↗
NASA NTRS1991

The use of direct numerical simulation data in turbulence modeling

Direct numerical simulations (DNS) of turbulent flows provide a complete data base to develop and to test turbulence models. In this article, the progress made in developing models for the dissipation rate equation is reviewed. New scaling arguments for the various terms in the dissipation rate equation were tested using data from DNS of homogeneous shear flows. Modifications to the epsilon-equation model that take into account near-wall effects were developed using DNS of turbulent channel flows. Testing of new models for flows under mean compression was carried out using data from DNS of isotropically compressed turbulence. In all of these studies the data from the simulations was essential in guiding the model development. The next generation of DNS will be at higher Reynolds numbers, and will undoubtedly lead to improved models for computations of flows of practical interest.

Mansour, N. N.↗
NASA NTRS1990

Numerical simulations of compressible mixing layers

Direct numerical simulations of two-dimensional temporally growing compressible mixing layers are presented. The Kelvin-Helmholtz instability is initially excited by a white-noise perturbation superimposed onto a hyperbolic tangent meanflow profile. The linear regime is studied at low resolution in the case of two flows of equal temperatures, for convective Mach numbers from 0.1 to 1 and for different values of the Reynolds number. At higher resolution, the complete evolution of a two-eddy mixing layer between two flows of different temperatures is simulated at moderate Reynolds number. Similarities and differences between flows of equal convective Mach numbers are discussed.

Normand, Xavier↗
NASA NTRS2003

Towards an Automated Full-Turbofan Engine Numerical Simulation

The objective of this study was to demonstrate the high-fidelity numerical simulation of a modern high-bypass turbofan engine. The simulation utilizes the Numerical Propulsion System Simulation (NPSS) thermodynamic cycle modeling system coupled to a high-fidelity full-engine model represented by a set of coupled three-dimensional computational fluid dynamic (CFD) component models. Boundary conditions from the balanced, steady-state cycle model are used to define component boundary conditions in the full-engine model. Operating characteristics of the three-dimensional component models are integrated into the cycle model via partial performance maps generated automatically from the CFD flow solutions using one-dimensional meanline turbomachinery programs. This paper reports on the progress made towards the full-engine simulation of the GE90-94B engine, highlighting the generation of the high-pressure compressor partial performance map. The ongoing work will provide a system to evaluate the steady and unsteady aerodynamic and mechanical interactions between engine components at design and off-design operating conditions.

Reed, John A.↗