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 685 records · Page 38
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Numerical simulation of lava flows: Applications to the terrestrial planets

Lava flows are the visible expression of the extrusion of volcanic materials on a variety of planetary surfaces. A computer program described by Ishihara et al. appears to be well suited for application to different environments, and we have undertaken tests to evaluate their approach. Our results are somewhat mixed; the program does reproduce reasonable lava flow behavior in many situations, but we have encountered some conditions common to planetary environments for which the current program is inadequate. Here we present our initial efforts to identify the 'parameter space' for reasonable numerical simulations of lava flows.

Zimbelman, James R.↗
DOE OSTI2018

Numerical Simulation and Industrial Experimental Research on the Coherent Jet with “CH{sub 4} + N{sub 2}” Mixed Fuel Gas

Coherent jet technology is widely used in the electric arc furnace (EAF) steelmaking process to deliver more energy and momentum into the molten steel bath. Meanwhile, the characteristics of a coherent jet using pure CH{sub 4} as the fuel gas have been well investigated in previous studies. To reduce the consumption of CH{sub 4}, coherent jet technology using “CH{sub 4} + N{sub 2}” mixed fuel gas instead of pure CH{sub 4} was proposed and studied in detail by numerical simulation in the present work. The Eddy Dissipation Concept model, which has detailed chemical kinetic mechanisms, was adopted to model the fuel gas combustion reactions. Experimental measurements were carried out to validate the accuracy of the computational model. The present study shows that the jet characteristics of the main oxygen improve along with the increase of the CH{sub 4} ratio in fuel gas and with the increase of the flow rate of fuel gas. When the CH{sub 4} ratio in the fuel gas is 25 pct, the fuel gas flow rate only has a limited influence on the jet characteristics, unlike the rest of the fuel gas compositions, because a high N{sub 2} proportion deteriorates the combustion performance and leads to severe incomplete combustion. Moreover, a false potential core phenomenon was observed and explained in the present study. Based on the average values, the jet length of a coherent jet with 75 pct CH{sub 4} can achieve 89.8 pct of that with 100 pct CH{sub 4}. Finally, an industrial experiment was carried out on a commercial 100t EAF using coherent jet with 75 pct CH{sub 4}, showing that the average CH{sub 4} consumption was reduced from 3.84 to 3.05 Nm{sup 3} t{sup −1} under the premise of no obvious changes in the other production indexes.

42 ENGINEERING↗
NASA NTRS1990

An Improved K-Epsilon Model for Near-Wall Turbulence and Comparison with Direct Numerical Simulation

An improved k-epsilon model for low Reynolds number turbulence near a wall is presented. The near-wall asymptotic behavior of the eddy viscosity and the pressure transport term in the turbulent kinetic energy equation is analyzed. Based on this analysis, a modified eddy viscosity model, having correct near-wall behavior, is suggested, and a model for the pressure transport term in the k-equation is proposed. In addition, a modeled dissipation rate equation is reformulated. Fully developed channel flows were used for model testing. The calculations using various k-epsilon models are compared with direct numerical simulations. The results show that the present k-epsilon model performs well in predicting the behavior of near-wall turbulence. Significant improvement over previous k-epsilon models is obtained.

Shih, T. H.↗
DOE OSTI2015

Mechanism of formation of wiggly compaction bands in porous sandstone: 2. Numerical simulation using discrete element method

Wiggly compaction bands in porous aeolian s andstone vary from chevron shape to wavy shape to nearly straight. In some outcrops these variations occur along a single band. Here, a bonded close-packed discrete element model is used to investigate what mechanical properties control the formation of wiggly compaction bands (CBs). To simulate the volumetric yielding failure of porous sandstone, a discrete element shrinks when the force state of one of its bonds reaches the yielding cap defined by the failure force and the aspect ratio (k) of the yielding ellipse. A Matlab code “MatDEM 3D ” has been developed on the basis of this enhanced discrete element method. Mechanical parameters of elements are chosen according to the elastic properties and the strengths of porous sandstone. In numerical simulations, the failure angle between the band segment and maximum principle stress decreases from 90° to approximately 45° as k increases from 0.5 to 2, and compaction bands vary from straight to chevron shape. With increasing strain, subsequent compaction occurs inside or beside compacted elements, which leads to further compaction and thickening of bands. The simulations indicate that a greater yielding stress promotes chevron CBs, and a greater cement strength promotes straight CBs. Combined with the microscopic analysis introduced in the companion paper, we conclude that the shape of wiggly CBs is controlled by the mechanical properties of sandstone, including the aspect ratio of the yielding ellipse, the cri tical yielding stress, and the cement strength, w hich are d etermined primarily by petrophysical attributes, e.g., grain sorting, porosity, and cementation.

58 GEOSCIENCES↗
DOE OSTI2018

Numerical Simulation and Model Analyses of In-situ Biogeochemical Remediation of Chromium at the LANL - 18306

Environmental management problems associated with contamination of groundwater resources with heavy metals (such as chromium) typically require complex biogeochemical model simulations and analysis of in situ remediation. Frequently, the in situ remediation is the only viable treatment option. In these cases, a three-dimensional reactive transport model of relevant biogeochemical processes is invaluable. Recently, we have developed a novel unique biogeochemical simulator called CHROTRAN (http://chrotran.lanl.gov; Hansen et al. 2017, Vesselinov et al. 2018a) capable to represent the full three-dimensional dynamics of processes associated with in situ biogeochemical remediation of heavy metals in heterogeneous aquifers. CHROTRAN accounts for spatial and temporal transients of (1) the heavy metal to be remediated, (2) introduced amendments (e.g., through well injection), and (3) biomass growth and decay. The introduced amendments can include (1) an electron donor (e.g., molasses), (2) a nontoxic-conservative bio-inhibitor (e.g., ethanol), and (3) a biocide (e.g., dithionite). In addition, direct abiotic reduction by donor-metal interaction as well as donor-driven biomass growth and bio-reduction are modeled in CHROTRAN. Other critical processes, such as donor sorption, bio-fouling and biomass death are also modeled. CHROTRAN can simulate the groundwater flow dynamics introduced due a series of injection and extraction wells and account for reduction of aquifer permeability near the wells due bio clogging. CHROTRAN is based on the existing PFLOTRAN (http://pflotran.org; Lichtner et al, 2017a) code framework, taking advantage of the reaction interface allowing for implementation of complex model features and leverage other facets of PFLOTRAN, such as multi-component geochemistry and high-performance computing. Here, we present numerical simulations and model analyses related to the Los Alamos National Laboratory (LANL) chromium site. The analyses are performed by coupling of CHROTRAN with MADS (Model Analysis and Decision Support; http://mads.lanl.gov; Vesselinov et al., 2018b). We present analyses related to uncertainty quantification of model predictions, sensitivity of model parameters and remedial decision support. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2014

Numerical simulation of flow regions in red mud separation thickener's feedwell by analysis of residence-time distribution

The residence-time distribution (RTD) and the compartment model were applied to characterizing the flow regions in red mud separation thickener’s feedwells. Combined with the experimental work, validated mathematical model as well as three-dimensional computational fluid dynamics (CFD) model was established to analyze the flow regions of feedwells on an industrial scale. The concept of RTD, although a well-known method for the characterization of mixing behavior in conventional mixers and reactors, is still a novel measure for the characterization of mixing in feedwells. Numerical simulation results show that the inlet feed rate and the aspect ratio of feedwells are the most critical parameters which affect the RTD of feedwell. Further simulation experiments were then carried out. Under the optimal operation conditions, the volume fraction of dead zone can reduce by10.8% and an increasement of mixing flow volume fraction by 6.5% is also observed. There is an optimum feed inlet rate depending on the feedwell design. The CFD model in conjunction with the RTD analysis then can be used as an effective tool in the design, evaluation and optimization of thickener feedwell in the red mud separation.

42 ENGINEERING↗
DOE OSTI2016

Numerical simulations of early-stage dynamics of electron bunches emitted from plasmonic photocathodes

High-brightness electron sources are a key ingredient to the development of compact accelerator-based light sources. The electron sources are commonly based on (linear) photoemission process where a laser pulse with proper wavelength impinges on the surface of a metallic or semiconductor photocathode. Very recently, the use of plasmonic cathodes - cathodes with a nano-patterned surface - have demonstrated great enhancement in quantum efficiencies (Li et al., 2013 [1]). Alternatively, this type of photocathodes could support the formation of structured beams composed of transversely-separated beamlets. In this paper we discuss numerical simulations of the early-stage beam dynamics of the emission process from plasmonic cathodes carried out using the Warp (Friedman et al., 2014 [2]) framework. Furthermore, the model is used to investigate the properties of beams emitted from these photocathode and subsequently combined with particle-in-cell simulations to explore the imaging of cathode pattern after acceleration in a radiofrequency gun.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Numerical Simulation of Electron Magnetohydrodynamics with Landau-quantized Electrons in Magnetar Crusts

Abstract In magnetar crusts, magnetic fields are sufficiently strong to confine electrons into a small to moderate number of quantized Landau levels. This can have a dramatic effect on the crust's thermodynamic properties, generating field-dependent de Haas–van Alphen oscillations. We previously argued that the large-amplitude oscillations of the magnetic susceptibility could enhance the ohmic dissipation of the magnetic field by continuously generating small-scale, rapidly dissipating field features. This could be important to magnetar field evolution and contribute to their observed higher temperatures. To study this, we performed quasi-3D numerical simulations of electron MHD in a representative volume of neutron star crust matter, for the first time including the magnetization and magnetic susceptibility resulting from Landau quantization. We find that the potential enhancement in the ohmic dissipation rate due to this effect can be a factor ∼3 for temperatures of the order of 10 8 K, and ∼4.5 for temperatures of the order of 5 × 10 7 K, depending on the magnetic field configuration. The nonlinear Hall term is crucial to this amplification: without it, the magnetic field decay is only enhanced by a factor ≲2 even at 5 × 10 7 K. These effects generate a high wavenumber plateau in the magnetic energy spectrum associated with the small-scale de Haas–van Alphen oscillations. Our results suggest that this mechanism could help explain the magnetar heating problem, though due to the effect's temperature-dependence, full magneto-thermal evolution simulations in a realistic stellar model are needed to judge whether it is viable explanation.

Rau, Peter B. (ORCID:0000000152209277)↗
NASA NTRS1989

Stereo image visualization of numerically simulated turbulence

Stereo imaging techniques are used in the current study to analyze the results of a direct Navier-Stokes simulation of a turbulent boundary layer. Coherent structures embedded within the unsteady, highly three-dimensional turbulent flow fields are significantly enhanced through the use of stereo computer graphic renderings. In the paper, the geometrical foundations of two different methods for generating stereo images are reviewed. Sufficient detail is provided to enable readers to code stereo display algorithms from scratch. The effects of varying the geometric parameters of the stereo algorithm are demonstrated. Sample image pairs from animations of the turbulence numerical simulation are presented in color. Stereo imagery, especially when used with time-evolving data, provides new insight into the physical processes responsible for the generation an maintenance of turbulence within boundary layers.

Robinson, Stephen K.↗
DOE OSTI2021

Unraveling the implications of finite specimen size on the interpretation of dynamic experiments for polycrystalline aluminum through direct numerical simulations

Normal and Pressure-shear plate impact (NPI and PSPI) experiments are popular experimental techniques for studying the mean-field macroscopic behavior of polycrystalline metals under high-rate dynamic loading. However, since both configurations rely upon geometry for subjecting the specimen to high strain rates, these experiments often involve a limited specimen size. Moreover, because of the inherent heterogeneities present within polycrystalline metals, it is difficult to ascertain if the size of the specimen and/or regions where measurements are made are sufficiently large for making representative inferences about the mean-field macroscopic properties from single-point velocity measurements. In the present study, we quantify the expected measurement variability on observable point measurements in NPI and PSPI experiments by carrying out direct numerical simulations (DNS) of statistically representative polycrystalline microstructures subjected to dynamic compression and compression-shear loading. In particular, we consider the role of specific material heterogeneities (e.g. the grain-to-grain difference in size, crystallographic orientation) on dispersion in the normal and transverse particle velocity records and on local fluctuations in key state variables (e.g. velocity, accumulated plastic strain) by incorporating these effects directly into a representative synthetic microstructure geometry and crystalline description of pure polycrystalline aluminum. The form of the present study is a large parametric investigation, consisting of ten ensembles of one hundred simulations. Each of the thousand simulations reflects a randomly realized synthetic microstructure in one of five cases of decreasing average grain size for the two loading configurations. Our analysis of the DNS results demonstrates that for both of these experimental configurations, the grain size directly correlates with the coefficient of variation (CV) in simulated point measurements, showing a convergent decrease in CV to zero (i.e. particle velocity record approaches the mean-field value) with decreasing grain size. Remarkably, the magnitude of variations in the particle velocity record is shown to be largest where the deviatoric stresses are most significant. In the case of NPI, this occurs at the elastic and plastic wavefront, whereas, in the case of PSPI, the magnitude of fluctuations are approximately constant throughout the experimental window time. The reasoning for the scatter in particle velocity due to the heterogeneous microstructure is demonstrated to be dependent on the mechanisms for accommodating deformation and on the interaction of reflection waves generated at sites of heterogeneities occurring at the scale of grains. Lastly, we develop a power-law description for the magnitude of scattering versus characteristic length, which provides a statistical framework for assessing the required number of grains per characteristic specimen dimension for minimizing scatter within these two experimental configurations (NPI, PSPI).

36 MATERIALS SCIENCE↗
DOE OSTI2022

A direct numerical simulation study of the dilution tolerance of propane combustion under spark-ignition engine conditions

Modern spark ignition internal combustion (IC) engines rely on highly diluted fuel-air mixtures to achieve high brake thermal efficiencies. To support this, new engine designs have introduced high stroke-to-bore ratios and cylinder head designs that promote high tumble flow and turbulence intensities. However, mixture dilution through exhaust gas recirculation (EGR) is limited by combustion instabilities manifested in the form of cycle-to-cycle variability. Propane has been observed to have superior EGR dilution tolerance than gasoline, which makes it a very competitive low-carbon fuel for the new IC engines without sacrificing efficiency. Two-dimensional direct numerical simulations (DNS) are performed with detailed chemistry to study and contrast the effect of turbulence intensity and dilution on propane and iso-octane premixed flames at high pressure conditions similar to those in-cylinder. A new reduced mechanism for propane consisting of 53 transported species and 17 quasi-steady state species is developed based on a previously published mechanism and used in these simulations. Three levels of turbulence intensity and two levels of exhaust gas dilution are chosen based on conditions relevant to IC engine operation. The DNS results are analyzed based on the evolution of the flame surface area and the statistics of its driving terms, which are found to be similar for both fuels when there is no dilution but considerably different under high dilution. The analysis of the DNS data provides fundamental insights into the underlying mechanisms for improved stability under dilution.

42 ENGINEERING↗
DOE OSTI2020

Failure investigation of nuclear grade POCO graphite target in high energy neutrino physics through numerical simulation

In many high energy physics neutrino beamlines, targets are made up of isotropic POCO graphite grade for production of neutrinos for high energy particle physics research and are bombarded with highly energetic pulsed proton beam. The pulsed proton beam with a small beam size creates thermal stress waves as well as radiation damage such as displacement damage, void formation, swelling, and gas formation amongst others. As a result, after a few years of operation one of the longest operated targets appears to have undergone bulk swelling and fractured along the beam centerline. A complex interaction of dynamic loading due to beam, material degradation due to irradiation are thought to have caused such a failure. Here, we present a numerical simulation to capture the combined effects of swelling and the dynamic effect of thermal stress waves due to beam heating, to explain the failure of target. An empirical formula has been developed to take into account swelling as a function of temperature and proton fluence and has been implemented in a commercial finite element code. Extensive X-ray diffraction (XRD) were performed on graphite samples from the failed target to understand the lattice parameter changes due to irradiation and build a valid empirical model. Simulation results show that swelling plays a major role in elevating the stress state in the material exceeding the failure strength while the pulsed beam heating also introduces fatigue loading that would explain brittle fracture emanating from inside of the material and propagating outward. Such empirical formulations will help in identifying the threshold values of critical parameters to extend the service life of future multi-megawatt neutrino targets and targets in other accelerator environments.

43 PARTICLE ACCELERATORS↗
DOE OSTI2017

High-stress study of bioinspired multifunctional PEDOT:PSS/nanoclay nanocomposites using AFM, SEM and numerical simulation

Bioinspired design has been central in the development of hierarchical nanocomposites. Particularly, the nacre-mimetic brick-and-mortar structure has shown excellent mechanical properties, as well as gas-barrier properties and optical transparency. Along with these intrinsic properties, the layered structure has also been utilized in sensing devices. Here we extend the multifunctionality of nacre-mimetics by designing an optically transparent and electron conductive coating based on PEDOT:PSS and nanoclays Laponite RD and Cloisite Na + . We carry out extensive characterization of the nanocomposite using transmittance spectra (trans-parency), conductive atomic force microscopy (conductivity), contact-resonance force microscopy (mechanical properties), and SEM combined with a variety of stress-strain AFM experiments and AFM numerical simulations (internal structure). We further study the nanoclay’s response to the application of pressure with multifrequency AFM and conductive AFM, whereby increases and decreases in conductivity can occur for the Laponite RD composites. We offer a possible mechanism to explain the changes in conductivity by modeling the coating as a 1-dimensional multibarrier potential for electron transport, and show that conductivity can change when the separation between the barriers changes under the application of pressure, and that the direction of the change depends on the energy of the electrons. We did not observe changes in conductivity under the application of pressure with AFM for the Cloisite Na+ nanocomposite, which has a large platelet size compared with the AFM probe diameter. No pressure-induced changes in conductivity were observed in the clay-free polymer either.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1996

Visualizing Time-Varying Phenomena In Numerical Simulations Of Unsteady Flows

Streamlines, contour lines, vector plots, and volume slices (cutting planes) are commonly used for flow visualization. These techniques are sometimes referred to as instantaneous flow visualization techniques because calculations are based on an instant of the flowfield in time. Although instantaneous flow visualization techniques are effective for depicting phenomena in steady flows,they sometimes do not adequately depict time-varying phenomena in unsteady flows. Streaklines and timelines are effective visualization techniques for depicting vortex shedding, vortex breakdown, and shock waves in unsteady flows. These techniques are examples of time-dependent flow visualization techniques, which are based on many instants of the flowfields in time. This paper describes the algorithms for computing streaklines and timelines. Using numerically simulated unsteady flows, streaklines and timelines are compared with streamlines, contour lines, and vector plots. It is shown that streaklines and timelines reveal vortex shedding and vortex breakdown more clearly than instantaneous flow visualization techniques.

Lane, David A.↗
NASA NTRS2013

Towards Bridging the Gaps in Holistic Transition Prediction via Numerical Simulations

The economic and environmental benefits of laminar flow technology via reduced fuel burn of subsonic and supersonic aircraft cannot be realized without minimizing the uncertainty in drag prediction in general and transition prediction in particular. Transition research under NASA's Aeronautical Sciences Project seeks to develop a validated set of variable fidelity prediction tools with known strengths and limitations, so as to enable "sufficiently" accurate transition prediction and practical transition control for future vehicle concepts. This paper provides a summary of selected research activities targeting the current gaps in high-fidelity transition prediction, specifically those related to the receptivity and laminar breakdown phases of crossflow induced transition in a subsonic swept-wing boundary layer. The results of direct numerical simulations are used to obtain an enhanced understanding of the laminar breakdown region as well as to validate reduced order prediction methods.

Choudhari, Meelan M.↗
DOE OSTI2014

Numerical simulations and experimental measurements of steel and ice impacts on concrete for acoustic interrogation of delaminations in bridge decks

Delaminations in bridge decks typically result from corrosion of the top mat of reinforcing steel, which leads to a localized separation of the concrete cover from the underlying concrete. Because delaminations cannot be detected using visual inspection, rapid, large-area interrogation methods are desired to characterize bridge decks without disruption to traffic, without the subjectivity inherent in existing methods, and with increased inspector safety. To this end, disposable impactors such as water droplets or ice chips can be dropped using automatic dispensers onto concrete surfaces to excite mechanical vibrations while acoustic responses can be recorded using air-coupled microphones. In this work, numerical simulations are used to characterize the flexural response of a model concrete bridge deck subject to both steel and ice impactors, and the results are compared with similar experiments performed in the laboratory on a partially delaminated concrete bridge deck slab. The simulations offer greater understanding of the kinetics of impacts and the responses of materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1982

A numerical simulation of an atmospheric vortex street

A simplified mesoscale mixed-layer model is utilized to study the flow past an ocean island during a case of cold air outbreak over the Kuroshio Current. Using data taken during AMTEX '75, the governing equations are integrated in time to simulate the development of a Karman vortex street downstream of the island of Cheju-do. The surface layer, inversion layer, and overlying stable layer are all parametrized in the model. Very weak sinusoidal waviness but no vortex street was found in the case of the island. For a mixed layer barely deep enough to allow no part of the mountain on the island to protrude, there was more pronounced waviness but still no vortices. This result may not hold for boundary layers with poor vertical mixing. A numerical simulation of the vortex street yielded good results in comparison. The speed ratio, spacing ratio, and Strouhal number agreed rather closely with observations.

Ruscher, P. H.↗
DOE OSTI2021

Evaluation of Energy Storage Potential of Unconventional Shale Reservoirs Using Numerical Simulation of Cyclic Gas Injection

Compressed air energy storage (CAES) stores energy as compressed air in underground formations, typically salt dome caverns. When electricity demand grows, the compressed air is released through a turbine to produce electricity. CAES in the US is limited to one plant built in 1991, due in part to the inherent risk and uncertainty of developing subsurface storage reservoirs. As an alternative to CAES, we propose using some of the hundreds of thousands of hydraulically fractured horizontal wells to store energy as compressed natural gas in unconventional shale reservoirs. To store energy, produced or “sales” natural gas is injected back into the formation using excess electricity and is later produced through an expander to generate electricity. To evaluate this concept, we performed numerical simulations of cyclic natural gas injection into unconventional shale reservoirs using cmg-gem commercial reservoir modeling software. We tested short-term (diurnal) and long-term (seasonal) energy storage potential by modeling well injection and production gas flowrates as a function of bottom-hole pressure. First, we developed a conceptual model of a single fracture stage in an unconventional shale reservoir to characterize reservoir behavior during cyclic injection and production. Next, we modeled cyclic injection in the Marcellus shale gas play using published data. Results indicate that Marcellus unconventional shale reservoirs could support both short- and long-term energy storage at capacities of 100–1000 kWe per well. The results indicate that energy storage in unconventional shale gas wells may be feasible and warrants further investigation.

25 ENERGY STORAGE↗