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 919 records · Page 51

Search for Dark Matter Satellites Using the FERMI-LAT

Numerical simulations based on the {Lambda}CDM model of cosmology predict a large number of as yet unobserved Galactic dark matter satellites. We report the results of a Large Area Telescope (LAT) search for these satellites via the {gamma}-ray emission expected from the annihilation of weakly interacting massive particle (WIMP) dark matter. Some dark matter satellites are expected to have hard {gamma}-ray spectra, finite angular extents, and a lack of counterparts at other wavelengths. We sought to identify LAT sources with these characteristics, focusing on {gamma}-ray spectra consistent with WIMP annihilation through the b{bar b} channel. We found no viable dark matter satellite candidates using one year of data, and we present a framework for interpreting this result in the context of numerical simulations to constrain the velocity-averaged annihilation cross section for a conventional 100 GeV WIMP annihilating through the b{bar b} channel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Droplet-Gas Interactions in Nanoelectrospray Multiphase Flow

Electrospray (ES) is production of a charged droplet plume that transits from an electrically biased supply capillary to a counter electrode under the influence of an electric field. Drag interactions with the surrounding gas decelerate the droplets while transferring the momentum from droplets to accelerate the surrounding gas, resulting in an induced gas flow. This work seeks to characterize the structure of the resulting gas jet and identify the key mechanisms defining the flow structure. The phenomenon of gas jets produced by momentum transfer from nano-electrospray (nES) plumes is explored with schlieren visualization, thermal anemometry, and numerical simulations. Schlieren visualization experiments provide information on the flow structure in support of simulation predictions, and the hot thermistor anemometry measurements of gas velocities outside the spray demonstrate quantitatively validated simulation results. The study reveals the formation of a moderately high velocity coaxial gas jet within the nES plume and provides insight into the overall flow structure of the induced flow. The multiphase electrohydrodynamic simulations enable numerical experimentation to explore the fundamental physics of coupled droplet-gas transport and the resulting flow structure. The simulations, confirmed by the experiments, reveal gas jetting induced by nES with a narrow (sub hundred micrometers in diameter) core originating from the nES liquid-jet breakup region and a surrounding larger-in-extent zone (several hundred micrometers in diameter) of lower velocity gas flow within the electrospray plume. This behavior is due to nES ejecting a stream of droplets from the tip of a narrow liquid cone-jet, where the combined effect of many small droplets transferring momentum to a confined region of gas yields a narrow but high-velocity gas stream. Furthermore, the micro-jets produced by nES have practical utility for mass spectrometry, 3D printing and fabrication, and thermal management.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

First-Principles Petascale Simulations for Predicting Deflagration to Detonation Transition in Hydrogen-Oxygen Mixtures

Hydrogen has emerged as an important fuel across a range of industries as a means of achieving energy independence and to reduce emissions. DDT and the resulting detonation waves in hydrogen-oxygen can have especially catastrophic consequences in a variety of industrial and energy producing settings related to hydrogen. First-principles numerical simulations of flame acceleration and DDT are required for an in-depth understanding of the phenomena and facilitating design of safe hydrogen systems. The goals of this project were (1) to develop first-principles petascale reactive flow Navier-Stokes simulation code for predicting gaseous high-speed combustion and detonation (HSCD) phenomena and (2) demonstrate feasibility of first-principles simulations of rapid flame acceleration and deflagration-to-detonation transition (DDT) in stoichiometric hydrogen-oxygen mixture (2H 2 + O 2 ). The goals of the project have been accomplished. We have developed a novel numerical simulation code, named HSCD, for performing first-principles direct numerical simulations of high-speed hydrogen combustion. We carried out a series of validating numerical simulations of inert and reactive shock reflection experiments in shock tubes. We then performed a pilot numerical simulation of flame acceleration in a long pipe. The simulation showed the transition of the rapidly accelerating flame into a detonation. The DDT simulations were performed using BG/Q Mira at the Argonne National Laboratory, currently the fourth fastest super-computer in the world.

08 HYDROGEN↗

Grating Polarizers at 170 GHz for ECRH Systems: Low Power Tests and Simulations

The elliptical polarization produced by a pair of sinusoidally profiled miter-bend grating polarizers has been measured experimentally and found to be in excellent agreement with theory. The polarizers were designed to be used at 170 GHz in the miter bends of a 63.5 mm corrugated metallic waveguide for creating the arbitrary elliptically polarized microwave beam needed for electron cyclotron resonance heating of plasma. Using a vector network analyzer to generate a linearly polarized HE11 incident mode, each polarizer was individually tested to measure the rotation α and ellipticity β of the elliptically polarized reflected microwave beam as a function of the grating rotation angle. The grating polarizers were then tested together to measure the elliptical polarization as a function of the angles of the combination of the two polarizers arranged sequentially on the transmission line. The map of the ellipticity versus grating angles agreed very well with numerical simulations using high-frequency structure simulator. Numerical simulations show that up to eight combinations of rotation angles of the two grating polarizers can provide the same ellipticity but with the varying ohmic loss, thus allowing the choice of settings to minimize the ohmic loss, a necessary consideration to minimize heating effects in megawatt power level systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Grating Polarizers at 170 GHz for ECRH Systems: Low Power Tests and Simulations

The elliptical polarization produced by a pair of sinusoidally profiled miter-bend grating polarizers has been measured experimentally and found to be in excellent agreement with theory. The polarizers were designed to be used at 170 GHz in the miter bends of a 63.5 mm corrugated metallic waveguide for creating the arbitrary elliptically polarized microwave beam needed for electron cyclotron resonance heating of plasma. Using a vector network analyzer to generate a linearly polarized HE 11 incident mode, each polarizer was individually tested to measure the rotation α and ellipticity β of the elliptically polarized reflected microwave beam as a function of the grating rotation angle. The grating polarizers were then tested together to measure the elliptical polarization as a function of the angles of the combination of the two polarizers arranged sequentially on the transmission line. The map of the ellipticity versus grating angles agreed very well with numerical simulations using high-frequency structure simulator. Numerical simulations show that up to eight combinations of rotation angles of the two grating polarizers can provide the same ellipticity but with the varying ohmic loss, thus allowing the choice of settings to minimize the ohmic loss, a necessary consideration to minimize heating effects in megawatt power level systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Synergism of the method of characteristics and CAD technology for neutron transport calculation

The method of characteristics (MOC) is a very popular methodology in neutron transport calculation and numerical simulation in recent decades for its unique advantages. One of the key problems determining whether the MOC can be applied in complicated and highly heterogeneous geometry is how to combine an effective geometry processing method with MOC. Most of the existing MOC codes describe the geometry by lines and arcs with extensive input data, such as circles, ellipses, regular polygons and combination of them. Thus they have difficulty in geometry modeling, background meshing and ray tracing for complicated geometry domains. In this study, a new idea making use of a CAD solid modeler MCAM which is a CAD/Image-based Automatic Modeling Program for Neutronics and Radiation Transport developed by FDS Team in China was introduced for geometry modeling and ray tracing of particle transport to remove these geometrical limitations mentioned above. The diamond-difference scheme was applied to MOC to reduce the spatial discretization error of the flat flux approximation in theory. Based on MCAM and MOC, a new MOC code was developed and integrated into SuperMC system, which is a Super Multi-function Computational system for neutronics and radiation simulation. The numerical testing results demonstrated the feasibility and effectiveness of the new idea for geometry treatment in SuperMC. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

INTEGRATION OF PARTICLE-GAS SYSTEMS WITH STIFF MUTUAL DRAG INTERACTION

Numerical simulation of numerous mm/cm-sized particles embedded in a gaseous disk has become an important tool in the study of planet formation and in understanding the dust distribution in observed protoplanetary disks. However, the mutual drag force between the gas and the particles can become so stiff—particularly because of small particles and/or strong local solid concentration—that an explicit integration of this system is computationally formidable. In this work, we consider the integration of the mutual drag force in a system of Eulerian gas and Lagrangian solid particles. Despite the entanglement between the gas and the particles under the particle-mesh construct, we are able to devise a numerical algorithm that effectively decomposes the globally coupled system of equations for the mutual drag force, and makes it possible to integrate this system on a cell-by-cell basis, which considerably reduces the computational task required. We use an analytical solution for the temporal evolution of each cell to relieve the time-step constraint posed by the mutual drag force, as well as to achieve the highest degree of accuracy. To validate our algorithm, we use an extensive suite of benchmarks with known solutions in one, two, and three dimensions, including the linear growth and the nonlinear saturation of the streaming instability. We demonstrate numerical convergence and satisfactory consistency in all cases. Our algorithm can, for example, be applied to model the evolution of the streaming instability with mm/cm-sized pebbles at high mass loading, which has important consequences for the formation scenarios of planetesimals.

79 ASTRONOMY AND ASTROPHYSICS↗

INTEGRATED WORKFLOW MANAGEMENT FOR PARTICLE ACCELERATOR SIMULATION

Supercomputing systems are used for a wide range of computationally demanding tasks in many fields of science and engineering. They play a key role in numerical simulation, in which mathematical models are computed in order to simulate the behavior of physical systems. Scientists and engineers that use supercomputers for numerical simulation often have their productivity limited by the need to manually organize and manage extremely large amounts of data that are often produced and consumed by the software programs run on these systems. Recognizing these limitations, Kitware Inc. (Clifton Park, NY) and SLAC National Accelerator Laboratory (Menlo Park, CA) are developing an advanced software platform that can reduce the cognitive overhead required by knowledge workers when using supercomputers for numerical simulation. Phase I of the project is complete and includes the development of new capabilities for organizing simulation project files, improvements to the user interface and overall usability, and deployment of a “middle tier” server to sit between user desktop machines and supercomputers to offload much of the data management workload. The project also developed prototype software for executing sequences of numerical simulations, and a prototype for migrating supercomputing software to cloud-based computing systems to provide a potential alternative to supercomputers with different logistical and price-to-performance tradeoffs.

Tourtellott, John↗

Development of the T+M coupled flow–geomechanical simulator to describe fracture propagation and coupled flow–thermal–geomechanical processes in tight/shale gas systems

We developed a hydraulic fracturing simulator by coupling a flow simulator to a geomechanics code, namely T+M simulator. Modeling of the vertical fracture development involves continuous updating of the boundary conditions and of the data connectivity, based on the finite element method for geomechanics. The T+M simulator can model the initial fracture development during the hydraulic fracturing operations, after which the domain description changes from single continuum to double or multiple continua in order to rigorously model both flow and geomechanics for fracture-rock matrix systems. The T+H simulator provides two-way coupling between fluid-heat flow and geomechanics, accounting for thermoporomechanics, treats nonlinear permeability and geomechanical moduli explicitly, and dynamically tracks changes in the fracture(s) and in the pore volume. We also fully accounts for leak-off in all directions during hydraulic fracturing. We first validate the T+M simulator, matching numerical solutions with the analytical solutions for poromechanical effects, static fractures, and fracture propagations. Then, from numerical simulation of various cases of the planar fracture propagation, shear failure can limit the vertical fracture propagation of tensile failure, because of leak-off into the reservoirs. Slow injection causes more leak-off, compared with fast injection, when the same amount of fluid is injected. Changes in initial total stress and contributions of shear effective stress to tensile failure can also affect formation of the fractured areas, and the geomechanical responses are still well-posed.

02 PETROLEUM↗

Interpenetration and kinetic effects in converging, high-energy plasma jets

In this work, we report on numerical simulations of laser-driven convergent plasma fusion targets. These “inverted corona” fusion targets are useful for the study of counter-streaming and converging rarefied plasma flows, and previous experiments have demonstrated their potential as neutron sources. The scheme consists of a fuel layer lined along the interior surface of a hollow plastic shell that is laser-ablated and expands inward towards the target center. The plasma streams generated in these targets are initially nearly collisionless as they converge, leading to wide interaction length scales and long interaction time scales as the jets interpenetrate. Such kinetic effects impact mixing of constituent ions - a phenomenon not properly captured by single-fluid hydrodynamic simulations. Here we conduct numerical simulations using two different methods: (1) single-fluid simulations in HYDRA, and (2) kinetic-ion, fluid-electron hybrid particle-in-cell (PIC) simulations in the code Chicago. It is shown that the initially nearly collisionless plasma fronts interpenetrate deeply and lead to broader interaction regions in space and time resulting in significant beam-beam fusion. The two approaches make different, testable predictions for the effect of fuel-layer thickness on neutron yield.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Estimating manufacturing electricity costs by simulating dependence between production parameters

The electricity costs associated with manufacturing generally consist of an energy charge and a demand charge, which are calculated by total kWh × kWh and peak kW × kW, respectively. While these two charges have been extensively studied with production scheduling, dependence between production parameters has not been thoroughly investigated for production energy costs when the production parameters are assumed to be independent. Such dependence, however, can significantly affect the peak kW and resulting energy costs. Therefore, industrial practitioners need a tool to simulate the production power demand and cost to evaluate energy costs for each production plan and for any dependence in the parameters. In this paper, we propose a method for simulating the power demand at the machine and facility levels, focusing on the dependence in milling machine parameters. Using probabilistic techniques with copulas, we quantify the degree of dependence between two milling parameters (processing times and processing amounts) and build a power demand simulator combining discrete event simulation and numerical simulation models. From 35 illustrative examples, we show that a production strategy to adjust dependence can reduce the peak kW and energy costs by more than 15% and 4%, respectively. In addition to dependence, less variability in processing times results in less peak kW and energy costs. Other factors such as a marginal probability distribution and copula types having the same degree of variability have been found to have limited effects on reducing energy costs.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Simulation of laser interaction with ablative plasma and hydrodynamic behavior of laser supported plasma

When an intense laser beam irradiates on a solid target, ambient air ionizes and becomes plasma, while part of the target rises in temperature, melts, vaporizes, ionizes, and yet becomes plasma. A general Godunov finite difference scheme WENO (Weighted Essentially Non-Oscillatory Scheme) with fifth-order accuracy is used to simulate 2-dimensional axis symmetrical laser-supported plasma flow field in the process of laser ablation. The model of the calculation of ionization degree of plasma and the interaction between laser beam and plasma are considered in the simulation. The numerical simulations obtain the profiles of temperature, density, and velocity at different times which show the evolvement of the ablative plasma. The simulated results show that the laser energy is strongly absorbed by plasma on target surface and that the velocity of laser supported detonation (LSD) wave is half of the ideal LSD value derived from Chapman-Jouguet detonation theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical solution of 2D-vector tomography problem using the method of approximate inverse

We propose a numerical solution of reconstruction problem of a two-dimensional vector field in a unit disk from the known values of the longitudinal and transverse ray transforms. The algorithm is based on the method of approximate inverse. Numerical simulations confirm that the proposed method yields good results of reconstruction of vector fields.

97 MATHEMATICS AND COMPUTING↗

Formation of the active medium in high-power repetitively pulsed gas lasers pumped by an electron-beam-controlled discharge

A high-power repetitively pulsed e-beam-controlled discharge CO{sub 2} laser is simulated numerically; the simulation results are compared with experimental data. Optimal sizes and design of electrodes and configuration of the external magnetic field are found, which allow one to introduce no less than 90% electric pump energy into a specified volume of the active medium, including the active volume of a laser with an aperture of 110 × 110 cm. The results obtained can also be used to design other types of highpower gas lasers. (lasers)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Modeling and Analysis of Eavesdropping Attack in 802.11ad mmWave Wireless Networks

Next generation wireless communication networks utilizing 60 GHz millimeter wave (mmWave) frequency bands are expected to achieve multi-gigabit throughput with the use of highly directional phased-array antennas. These directional signal beams provide enhanced security to the legitimate networks due to the increased difficulties of eavesdropping. However, there still exists significant possibility of eavesdropping since 1) the reflections of the signal beam from ambient reflectors enables opportunistic stationary eavesdropping attacks, and; 2) carefully designed beam exploration strategy enables active nomadic eavesdropping attack. This paper discusses eavesdropper attack strategies for 802.11ad mmWave systems and provides the first analytical model to characterize the success possibility of eavesdropping in both opportunistic stationary attacks and active nomadic attacks. We derive the success probability of eavesdropping considering the ambient reflectors in the environment and errors introduced in the beam exploration strategies of the proposed eavesdropping attacker models. We study the success probability for both opportunistic stationary attack scenario and active nomadic attack scenario through numerical simulations. In addition to numerical simulations, we also evaluate the proposed attacker models using an 802.11ad test bed consisting of commercially available off-the-shelf devices.

99 GENERAL AND MISCELLANEOUS↗

Numerical Study of the Voltage Waveform Effect on the Spatiotemporal Characteristics of a Dielectric Barrier Microdischarge in Argon

The effect of the shape of the feeding voltage on the spatiotemporal characteristics of an atmospheric- pressure barrier microdischarge in argon is demonstrated using numerical simulations based on an extended fluid model. Results of simulations performed for sinusoidal and square feeding voltages are analyzed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Prefabricated High-Strength Rebar Systems with High-Performance Concrete for Accelerated Construction of Nuclear Concrete Structures

This project aimed to reduce the field construction times and fabrication costs of reinforced concrete (RC) nuclear structures through the use of: 1) high-strength steel deformed reinforcing bar (rebar); 2) prefabricated high-strength rebar (HSR) assemblies with headed anchorages; and 3) high-strength concrete (HSC). The focus of the project was not to develop new innovations on materials, but rather to fill a major knowledge gap on the effectiveness, code conformity, and viability (i.e., practicality, commercial availability) of existing high-strength materials in nuclear structures, especially in shear walls, their connections/joints, and around penetrations. These advances promote the Nuclear Energy Enabling Technologies Advanced Methods for Manufacturing (NEET-1) goals of accelerating deployment schedules by at least 6 months and reducing component fabrication costs by 20% or more. Toward this vision, the following objectives were identified and completed: Develop limit-benefit and cost-benefit analysis frameworks: The structural and economic limits and benefits of HSR and HSC in nuclear walls were evaluated. Importantly, it was found that modest increases in concrete strength up to 10.0 ksi and rebar yield strength up to 100 ksi result in the most economic designs and improved structural performance (Barbachyn et al., 2017a and Barbachyn et al., 2017b). Also, an industry survey was distributed to determine the cost benefits associated with prefabrication. Importantly, on-site worker-hours are significantly reduced, allowing for significantly accelerated construction schedules (Devine et al. 2018b). Develop an optimization methodology for structural design: An optimization study was deemed not necessary since the trivial solution with maximized use of high-strength materials and prefabrication always resulted in the largest cost and time savings. Instead, the benefits (and limits) of high-strength materials and prefabrication on construction costs and on-site construction time were investigated through a parametric study. Cost-reductions through the use of HSR and HSC were found to be greater than 40%. Also, on-site construction time of nuclear shear walls can be reduced by up to 83% through the use of prefabricated rebar assemblies. Conduct experimental evaluations of nuclear shear walls and wall-foundation joints: A large scale experimental program was completed on shear-critical RC cantilever deep beams. From these experiments, it was demonstrated that increasing rebar yield strength had a greater effect on increasing lateral strength than increasing the concrete compressive strength. Further, using HSR together with HSC resulted in the highest lateral strength and displacement ductility (Devine et al. 2018a). Additionally, large-scale reversed-cyclic lateral load tests were conducted on RC shear walls. It was demonstrated that through the utilization of HSR and HSC, a wall with 55% reduced vertical and horizontal rebar area was able to achieve 91% of the peak lateral strength of a conventional wall using state-of-practice materials with the same geometry. Further, trim reinforcement with headed anchorages around the penetrations proved to be effective. Develop validated numerical simulation models: Numerical modeling was completed in two phases during this project using 3 different finite element software packages. The first phase used existing measured data from previous experimental studies on stocky RC walls to validate simplified predictive models. Models using two different programs were validated based on a database of 38 stocky walls (Barbachyn et al., 2017a). The second phase occurred after experimental testing using a third and more detailed finite element package. Prediction improvements for the initial stiffness, lateral deformations, damage propagation, and failure mechanisms were targeted. The resulting detailed models more accurately predicted the peak lateral strength, initial and secant stiffnesses, lateral deformations, and cracking patterns of the wall test specimens subjected to reversed-cyclic lateral loading as compared to the simplified pretest models. Develop validated design procedures, tools, and criteria: This project evaluated the accuracy of existing closed-form design equations from U.S. building design standards and previous research (Barbachyn et al. 2017a and Devine et al. 2018a). It was shown that the current predictive equations for shear strength capacity need to be revisited as they provide significant scatter and gross overestimates, indicating the potential for unconservative designs. Further, for both the deep beam and shear wall specimens, the numerical models and flexural capacity predictions provided the most accurate strength predictions. Develop cost-effective field construction procedures: To promote the use of prefabricated rebar assemblies, the following construction conditions where prefabrication would be used instead of in-place rebar tying were identified from an industry survey: 1) to save on construction schedule; 2) to improve safety and/or quality control; 3) for areas with heavy rebar congestion; and 4) for structures with significant repetition in rebar layout. All of these conditions directly apply and benefit nuclear construction where construction costs and schedules are extensive, safety and quality are of utmost importance, the reinforcement is extremely dense, and the reinforcement layout is repetitive in many regions of the shear walls. Additionally, an experimental study was conducted to investigate the changes in bar spacing (with respect to code-required rebar spacing and placement tolerances) when a horizontally assembled prefabricated two-dimensional (2D) rebar mat or three-dimensional (3D) rebar cage is lifted (or tripped) to a vertical position. It was found that the bars involved in both the lateral translation of the rebar assemblies (i.e., movement of the assembly in its horizontal position) and the tripping process were most susceptible to large spacing changes, and thus, their positions within the assembly would need to be checked after the placement of the assembly in place within the formwork. Importantly, the spacings of all of the other bars in the assembly were not affected significantly by the lateral translation and tripping process, thereby eliminating the need to recheck their positions within the assembly (Devine et al. 2018b).

36 MATERIALS SCIENCE↗