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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2012

Numerical simulation of current-free double layers created in a helicon plasma device

Two-dimensional simulations reveal that when radially confined source plasma with magnetized electrons and unmagnetized ions expands into diverging magnetic field B, a current-free double layer (CFDL) embedded in a conical density structure forms, as experimentally measured in the Australian helicon plasma device (HPD). The magnetized electrons follow the diverging B while the unmagnetized ions tend to flow directly downstream of the source, resulting in a radial electric field (E{sub Up-Tack }) structure, which couples the ion and electron flows. Ions are transversely (radially) accelerated by E{sub Up-Tack} on the high potential side of the double layer in the CFDL. The accelerated ions are trapped near the conical surface, where E{sub Up-Tack} reverses direction. The potential structure of the CFDL is U-shaped and the plasma density is enhanced on the conical surface. The plasma density is severely depleted downstream of the parallel potential drop ({phi}{sub Double-Vertical-Line Double-Vertical-Line o}) in the CFDL; the density depletion and the potential drop are related by quasi-neutrality condition, including the divergence in the magnetic field and in the plasma flow in the conical structure. The potential and density structures, the CFDL spatial size, its electric field strengths and the electron and ion velocities and energy distributions in the CFDL are found to be in good agreements with those measured in the Australian experiment. The applicability of our results to measured axial potential profiles in magnetic nozzle experiments in HPDs is discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Coalescence kinetics in surfactant stabilized emulsions: Evolution equations from direct numerical simulations

Lattice Boltzmann simulations were used to study the coalescence kinetics in emulsions with amphiphilic surfactant, under neutrally buoyant conditions, and with a significant kinematic viscosity contrast between the phases (emulating water in oil emulsions). The 3D simulation domain was large enough (256 3rd power -- 10 7th power grid points) to obtain good statistics with droplet numbers ranging from a few thousand at early times to a few hundred near equilibrium. Increased surfactant contents slowed down the coalescence rate between droplets due to the Gibbs-Marangoni effect, and the coalescence was driven by a quasi-turbulent velocity field. The kinetic energy decayed at a relatively slow rate at early times, due to conversion of interfacial energy to kinetic energy in the flow during coalescence. Phenomenological, coupled differential equations for the mean droplet diameter D(t) and the number density nd(t) were obtained from the simulation data and from film draining theories. Local (in time) power law exponents for the growth of the mean diameter (and for the concomitant decrease of nd) were established in terms of the instantaneous values of the kinetic energy, coalescence probability, Gibbs elasticity, and interfacial area. The model studies indicated that true power laws for the growth of the droplet size and decrease of the number of droplets with time may not be justified, since the exponents derived using the phenomenological model were time dependent. In contrast to earlier simulation results for symmetric blends with surfactant, we found no evidence for stretched logarithmic scaling of the formD -- [ln (ct)]a for the morphology length, or exponential scalings associated with arrested growth, on the basis of the phenomenological model.

58 GEOSCIENCES↗
DOE OSTI2014

Numerical simulation of ion charge breeding in electron beam ion source

The Electron Beam Ion Source particle-in-cell code (EBIS-PIC) tracks ions in an EBIS electron beam while updating electric potential self-consistently and atomic processes by the Monte Carlo method. Recent improvements to the code are reported in this paper. The ionization module has been improved by using experimental ionization energies and shell effects. The acceptance of injected ions and the emittance of extracted ion beam are calculated by extending EBIS-PIC to the beam line transport region. An EBIS-PIC simulation is performed for a Cs charge-breeding experiment at BNL. The charge state distribution agrees well with experiments, and additional simulation results of radial profiles and velocity space distributions of the trapped ions are presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2013

Numerical simulation of temperature field, microstructure evolution and mechanical properties of HSS during hot stamping

The hot stamping of boron steels is widely used to produce ultra high strength automobile components without any spring back. The ultra high strength of final products is attributed to the fully martensitic microstructure that is obtained through the simultaneous forming and quenching of the hot blanks after austenization. In the present study, a mathematical model incorporating both heat transfer and the transformation of austenite is presented. A FORTRAN program based on finite element technique has been developed which permits the temperature distribution and microstructure evolution of high strength steel during hot stamping process. Two empirical diffusion-dependent transformation models under isothermal conditions were employed respectively, and the prediction capability on mechanical properties of the models were compared with the hot stamping experiment of an automobile B-pillar part.

36 MATERIALS SCIENCE↗
DOE OSTI2014

A numerical simulation study of gallium-phosphide/silicon heterojunction passivated emitter and rear solar cells

The performance of passivated emitter and rear (PERC) solar cells made of p-type Si wafers is often limited by recombination in the phosphorus-doped emitter. To overcome this limitation, a realistic PERC solar cell is simulated, whereby the conventional phosphorus-doped emitter is replaced by a thin, crystalline gallium phosphide (GaP) layer. The resulting GaP/Si PERC cell is compared to Si PERC cells, which have (i) a standard POCl{sub 3} diffused emitter, (ii) a solid-state diffused emitter, or (iii) a high efficiency ion-implanted emitter. The maximum efficiencies for these realistic PERC cells are between 20.5% and 21.2% for the phosphorus-doped emitters (i)–(iii), and up to 21.6% for the GaP emitter. The major advantage of this GaP hetero-emitter is a significantly reduced recombination loss, resulting in a higher V{sub oc}. This is so because the high valence band offset between GaP and Si acts as a nearly ideal minority carrier blocker. This effect is comparable to amorphous Si. However, the GaP layer can be contacted with metal fingers like crystalline Si, so no conductive oxide is necessary. Compared to the conventional PERC structure, the GaP/Si PERC cell requires a lower Si base doping density, which reduces the impact of the boron-oxygen complexes. Despite the lower base doping, fewer rear local contacts are necessary. This is so because the GaP emitter shows reduced recombination, leading to a higher minority electron density in the base and, in turn, to a higher base conductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Numerical simulation of the sustaining discharge in radio frequency hollow cathode discharge in argon

In this paper, a two-dimensional fluid model was developed to study the radio frequency (RF) hollow cathode discharge (HCD) in argon at 1 Torr. The evolutions of the particle density distribution and the ionization rate distribution in RF HCD at 13.56 MHz indicate that the discharge mainly occurs inside the hollow cathode. The spatio-temporal distributions of the ionization rate and the power deposition within the hollow cathode imply that sheath oscillation heating is the primary mechanism to sustain the RF HCD, whereas secondary electron emission plays a negligible role. However, as driving frequency decreases, secondary electron heating becomes a dominant mechanism to sustain the discharge in RF hollow cathode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Note: Numerical simulation and experimental validation of accelerating voltage formation for a pulsed electron accelerator

This paper describes the development of a computation model of a pulsed voltage generator for a repetitive electron accelerator. The model is based on a principle circuit of the generator, supplemented with the parasitics elements of the construction. Verification of the principle model was achieved by comparison of simulation with experimental results, where reasonable agreement was demonstrated for a wide range of generator load resistance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2014

Numerical simulation: Toward the design of high-efficiency planar perovskite solar cells

Organo-metal halide perovskite solar cells based on planar architecture have been reported to achieve remarkably high power conversion efficiency (PCE, >16%), rendering them highly competitive to the conventional silicon based solar cells. A thorough understanding of the role of each component in solar cells and their effects as a whole is still required for further improvement in PCE. In this work, the planar heterojunction-based perovskite solar cells were simulated with the program AMPS (analysis of microelectronic and photonic structures)-1D. Simulation results revealed a great dependence of PCE on the thickness and defect density of the perovskite layer. Meanwhile, parameters including the work function of the back contact as well as the hole mobility and acceptor density in hole transport materials were identified to significantly influence the performance of the device. Strikingly, an efficiency over 20% was obtained under the moderate simulation conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Numerical simulations used for a validity check on the laser induced photo-detachment diagnostic method in electronegative plasmas

Laser photo-detachment is used as a method to measure or determine the negative ion density and temperature in electronegative plasmas. In essence, the method consists of producing an electropositive channel (negative ion free region) via pulsed laser photo-detachment within an electronegative plasma bulk. Electrostatic probes placed in this channel measure the change in the electron density. A second pulse might be used to track the negative ion recovery. From this, the negative ion density and temperature can be determined. We study the formation and relaxation of the electropositive channel via a two-dimensional Particle-In-Cell/Mote Carlo collision model. The simulation is mainly carried out in a Hydrogen plasma with an electronegativity of α = 1, with a parametric study for α up to 20. The temporal and spatial evolution of the plasma potential and the electron densities shows the formation of a double layer (DL) confining the photo-detached electrons within the electropositive channel. This DL evolves into two fronts that move in the opposite directions inside and outside of the laser spot region. As a consequence, within the laser spot region, the background and photo-detached electron energy distribution function relaxes/thermalizes via collisionless effects such as Fermi acceleration and Landau damping. Moreover, the simulations show that collisional effects and the DL electric field strength might play a non-negligible role in the negative ion recovery within the laser spot region, leading to a two-temperature negative ion distribution. The latter result might have important effects in the determination of the negative ion density and temperature from laser photo detachment diagnostic.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Numerical simulation of heat and mass transfer in unsteady nanofluid between two orthogonally moving porous coaxial disks

The paper deals with the study of heat and mass transfer in an unsteady viscous incompressible water-based nanofluid (containing Titanium dioxide nanoparticles) between two orthogonally moving porous coaxial disks with suction. A combination of iterative (successive over relaxation) and a direct method is employed for solving the sparse systems of linear algebraic equations arising from the FD discretization of the linearized self similar ODEs. It has been noticed that the rate of mass transfer at the disks decreases with the permeability Reynolds number whether the disks are approaching or receding. The findings of the present investigation may be beneficial for the electronic industry in maintaining the electronic components under effective and safe operational conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Numerical simulation of gas dynamics and heat exchange tasks in fuel assemblies of the nuclear reactors

This report presents a PC-based program for solution gas dynamics and heat exchange mathematical tasks in fuel assemblies of the fast-neutron nuclear reactors. A fuel assembly consisting of bulk heat-generating elements, which are integrated together by the system of supply and pressure manifolds, is examined. Spherical heat-generating microelements, which contain nuclear fuel, are pulled into the heat-generating elements. Gaseous coolant proceed from supply manifolds to heat-generating elements, where it withdraws the nuclear reaction heat and assembles in pressure manifolds.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

3D numerical simulation of the long range propagation of acoustical shock waves through a heterogeneous and moving medium

Many situations involve the propagation of acoustical shock waves through flows. Natural sources such as lightning, volcano explosions, or meteoroid atmospheric entries, emit loud, low frequency, and impulsive sound that is influenced by atmospheric wind and turbulence. The sonic boom produced by a supersonic aircraft and explosion noises are examples of intense anthropogenic sources in the atmosphere. The Buzz-Saw-Noise produced by turbo-engine fan blades rotating at supersonic speed also propagates in a fast flow within the engine nacelle. Simulating these situations is challenging, given the 3D nature of the problem, the long range propagation distances relative to the central wavelength, the strongly nonlinear behavior of shocks associated to a wide-band spectrum, and finally the key role of the flow motion. With this in view, the so-called FLHOWARD (acronym for FLow and Heterogeneous One-Way Approximation for Resolution of Diffraction) method is presented with three-dimensional applications. A scalar nonlinear wave equation is established in the framework of atmospheric applications, assuming weak heterogeneities and a slow wind. It takes into account diffraction, absorption and relaxation properties of the atmosphere, quadratic nonlinearities including weak shock waves, heterogeneities of the medium in sound speed and density, and presence of a flow (assuming a mean stratified wind and 3D turbulent ? flow fluctuations of smaller amplitude). This equation is solved in the framework of the one-way method. A split-step technique allows the splitting of the non-linear wave equation into simpler equations, each corresponding to a physical effect. Each sub-equation is solved using an analytical method if possible, and finite-differences otherwise. Nonlinear effects are solved in the time domain, and others in the frequency domain. Homogeneous diffraction is handled by means of the angular spectrum method. Ground is assumed perfectly flat and rigid. Due to the 3D aspect, the code was massively parallelized using the single program, multiple data paradigm with the Message Passing Interfaces (MPI) for distributed memory architectures. This allows us to handle problems in the order of a thousand billion mesh points in the four dimensions (3 dimensions of space plus time). The validity of the method has been thoroughly evaluated on many cases with known solutions: linear piston, scattering of plane wave by a heterogeneous sphere, propagation in a waveguide with a shear flow, scattering by a finite amplitude vortex and nonlinear propagation in a thermoviscous medium. This validation process allows for a detailed assessment of the advantages and limitations of the method. Finally, applications to atmospheric propagation of shock waves will be presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Numerical simulations of the first operational conditions of the negative ion test facility SPIDER

In view of the realization of the negative ion beam injectors for ITER, a test facility, named SPIDER, is under construction in Padova (Italy) to study and optimize production and extraction of negative ions. The present paper is devoted to the analysis of the expected first operations of SPIDER in terms of single-beamlet and multiple-beamlet simulations of the hydrogen beam optics in various operational conditions. The effectiveness of the methods adopted to compensate for the magnetic deflection of the particles is also assessed. Indications for a sequence of the experimental activities are obtained.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Numerical simulation and experiment on split tungsten carbide cylinder of high pressure apparatus

A new high pressure device with a split cylinder was investigated on the basis of the belt-type apparatus. The belt-type die is subjected to excessive tangential tensile stress and the tungsten carbide cylinder is easily damaged in the running process. Taking into account the operating conditions and material properties of the tungsten carbide cylinder, it is divided into 6 blocks to eliminate the tangential tensile stress. We studied two forms of the split type: radial split and tangential split. Simulation results indicate that the split cylinder has more uniform stress distribution and smaller equivalent stress compared with the belt-type cylinder. The inner wall of the tangential split cylinder is in the situation that compressive stress is distributed in the axial, radial, and tangential directions. It is similar to the condition of hydrostatic pressure, and it is the best condition for tungsten carbide materials. The experimental results also verify that the tangential split die can bear the highest chamber pressure. Therefore, the tangential split structure can increase the pressure bearing capacity significantly.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗