Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computerized 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 1,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Ab Initio and Monte Carlo Approaches for the Magnetocaloric Effect in BaMnO{sub 3} Oxide Perovskite

The self-consistent ab initio calculations, based on density functional theory approach (DFT) and using full potential linearised augmented plane wave (FLAPW) method, have been used to investigate both electronic and magnetic properties of the BaMnO{sub 3} perovskite. Spin-polarised calculations, including the spin-orbit interaction, are used to determine the energy of the ferromagnetic (FM) and antiferromagnetic (AFM) states of BaMnO{sub 3} perovskite. Obtained data from ab initio calculations are used as input for the Monte Carlo simulations to compute other magnetic parameters. Magnetisation, specific heat and magnetic entropy change have been given using the Monte Carlo simulations. The adiabatic temperature change, transition temperature and relative cooling power have been established.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2017

Modeling and Simulation of the Off-gas in an Electric Arc Furnace

The following paper describes an approach to process modeling and simulation of the gas phase in an electric arc furnace (EAF). The work presented represents the continuation of research by Logar, Dovžan, and Škrjanc on modeling the heat and mass transfer and the thermochemistry in an EAF. Due to the lack of off-gas measurements, Logar et al. modeled a simplified gas phase under consideration of five gas components and simplified chemical reactions. The off-gas is one of the main continuously measurable EAF process values and the off-gas flow represents a heat loss up to 30 pct of the entire EAF energy input. Therefore, gas phase modeling offers further development opportunities for future EAF optimization. This paper presents the enhancement of the previous EAF gas phase modeling by the consideration of additional gas components and a more detailed heat and mass transfer modeling. In order to avoid the increase of simulation time due to more complex modeling, the EAF model has been newly implemented to use an efficient numerical solver for ordinary differential equations. Compared to the original model, the chemical components H{sub 2}, H{sub 2}O, and CH{sub 4} are included in the gas phase and equilibrium reactions are implemented. The results show high levels of similarity between the measured operational data from an industrial scale EAF and the theoretical data from the simulation within a reasonable simulation time. In the future, the dynamic EAF model will be applicable for on- and offline optimizations, e.g., to analyze alternative input materials and mode of operations.

42 ENGINEERING↗
DOE OSTI2018

Motion of Solid Grains During Magnetic Field-Assisted Directional Solidification

In this paper, we report the visible evidence for thermoelectric magnetic forces (TEMFs) during magnetic field-assisted directional solidification, and their potential to control the motion of solid grains(dendrite fragments or equiaxed grains). These motions are observed by means of synchrotron X-ray radiography and compared with analytic calculations for a spherical particle’s motion driven only by TEMFs, which confirms that the observed solid grain motions are the combined result of the TEMFs and gravity. We also carried out corresponding 3D numerical simulations to validate the calculations and further prove our conclusion that TEMF acts on the solid grain and affects its motion trajectory.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Numerical Simulation of Recirculating Flow and Physical Model of Slag–Metal Behavior in an RH Reactor: Application to Desulfurization

Computational fluid dynamics (CFD) techniques and a 1:7.5 physical model of an RH degasser have been used to evaluate the flow of gas and metal inside an RH reactor for vacuum degassing of liquid steel. The effect of gas injection on the gas spatial distribution, steel circulation rate and flow field inside the ladle, snorkels and vacuum chamber have been assessed. N-pentane oil was employed to evaluate the average residence time as well as the slag droplet size distribution. The predicted radial gas distribution and liquid circulation rate have been validated against experimental data from a physical model. The results with incorporation of the virtual mass force coefficient of 0.25 and the turbulence dispersion force showed better predictions of gas distribution in the up-snorkel as well as circulation rate. Full-scale simulations were performed, and the predicted circulation rate was significantly affected by the argon bubble expansion. Data from these simulations were used to analyze the degree of desulfurization performed by the addition of desulfurizing agents inside the vacuum chamber. A model of the kinetics of desulfurization based on the results from the physical model and CFD simulation and on slag dispersion inside liquid steel yields degrees of desulfurization similar to the industrial trials reported in the literature.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Investigation on Cracking Behavior and Development of a Fracture Model of Ti-47Al-2Nb-2Cr Alloy During Hot Deformation

When alloys are strained to a certain limit under the plastic deformation, ductile fracture occurs. Whether the required deformation can be put into effect before fracture appears has become a serious concern. A combination of compression tests, ductile fracture criteria as well as finite element simulation (FEM) were employed in this study to establish the fracture initiation prediction model of Ti-47Al-2Nb-2Cr alloy. To obtain the critical fracture strain, a technique was proposed to calculate the data for deformed cylindrical samples through isothermal compression tests in the temperature range between 1223 and 1423 K with strain rates ranging from 0.001 to 1 s{sup −1}. Five commonly used ductile fracture criteria were compared by FE simulation of compressing tests. It was found that the Brozzo criterion was found to be suitable and accordingly a failure criterion, considering temperature and strain ratio was established. Hot working parameters (critical fracture value and critical strain) in deformation could be obtained through the finite element simulation, revealing that this newly established failure criterion could accurately predict cracking initiation in compressing deformation of Ti-47Al-2Nb-2Cr alloy ingot.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Almost periodic oscillations for delay impulsive stochastic Nicholson’s blowflies timescale model

In this paper, we employ the new concepts of piecewise almost periodic stochastic process on a new periodic time scale and study a type of delay impulsive stochastic Nicholson’s blowflies timescale model under shifts operators. Some sufficient conditions for the existence and exponential stability of piecewise mean-square almost periodic solutions to the model are established. In addition, we provide numerical examples and computer simulations to demonstrate the feasibility of the obtained results. The methods and definitions introduced in the paper can be applied to study other types of stochastic dynamic models involving q-difference models and more.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

Artifacts in X-ray Dark-Field Tomography

Grating-based x-ray imaging provides three principle kinds of contrast: absorption, phase, and dark-field. Due to the availability of tomographic reconstruction algorithms for the dark-field contrast, it is now possible to take advantage of quantitative scatter information. However, the published algorithm is based on several assumptions that might be violated in reality. We use numerical simulations in order to identify artifacts in the reconstructions, which is crucial for the interpretation of experimental data.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2011

Efficient generation of fast ions from surface modulated nanostructure targets irradiated by high intensity short-pulse lasers

It's shown that the imposition of sub-laser wavelength relief structures on the surface of mass-limited-targets results into several folds higher short-pulse laser absorption, and consequently the efficient generation of fast ions. The optimum relief parameters for enhanced short-pulse laser absorption and higher ion acceleration are estimated numerically by particle-in-cell simulations and then corroborated by analytical scalings. The stability of the pre-imposed surface modulation during the laser pulse foil interaction is also examined.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Microwave diagnostics of femtosecond laser-generated plasma filaments

We present a simple non-intrusive experimental method allowing a complete single shot temporal measurement of laser produced plasma filament conductivity. The method is based on filament interaction with low intensity microwave radiation in a rectangular waveguide. The suggested diagnostics allow a complete single shot temporal analysis of filament plasma decay with resolution better than 0.3 ns and high spatial resolution along the filament. The experimental results are compared to numerical simulations, and an initial electron density of 7 x 10{sup 16 }cm{sup -3} and decay time of 3 ns are obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Molecular Dynamics Simulation of Binary Fluid in a Nanochannel

This paper presents the results from a molecular dynamics simulation of binary fluid (mixture of argon and krypton) in the nanochannel flow. The computational software LAMMPS is used for carrying out the molecular dynamics simulations. Binary fluids of argon and krypton with varying concentration of atom species were taken for two densities 0.65 and 0.45. The fluid flow takes place between two parallel plates and is bounded by horizontal walls in one direction and periodic boundary conditions are imposed in the other two directions. To drive the flow, a constant force is applied in one direction. Each fluid atom interacts with other fluid atoms and wall atoms through Week-Chandler-Anderson (WCA) potential. The velocity profile has been looked at for three nanochannel widths i.e for 12{sigma}, 14{sigma} and 16{sigma} and also for the different concentration of two species. The velocity profile of the binary fluid predicted by the simulations agrees with the quadratic shape of the analytical solution of a Poiseuille flow in continuum theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2011

Demixing-stimulated lane formation in binary complex plasma

Recently lane formation and phase separation have been reported for experiments with binary complex plasmas in the PK3-Plus laboratory onboard the International Space Station (ISS). Positive non-additivity of particle interactions is known to stimulate phase separation (demixing), but its effect on lane formation is unknown. In this work, we used Langevin dynamics (LD) simulation to probe the role of non-additivity interactions on lane formation. The competition between laning and demixing leads to thicker lanes. Analysis based on anisotropic scaling indices reveals a crossover from normal laning mode to a demixing-stimulated laning mode. Extensive numerical simulations enabled us to identify a critical value of the non-additivity parameter {Delta} for the crossover.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Stray RF Power Estimates From EC Exploitation During ITER Plasma Operations

The EC H and CD system of ITER tokamak is an essential tool for all the phases of ITER operation. Different levels of EC power are required through all the plasma discharge: up to 6.7 MW for assisting the breakdown and burn through, up to 20 MW for current drive and saw-teeth control from the equatorial launcher and up to 20 MW for NTM stabilization from the upper launchers. The assistance to breakdown and burn through is characterized by a very low (if not negligible) RF power absorption by the plasma. A significant level of stray radiation may also arise from partial absorption due to non-optimal plasma parameters and/or wrong injected polarization. The stray power radiated in the vacuum chamber is estimated as a first step toward mitigating potential harmful consequence to in-vessel structures and diagnostics. Power loading of the chamber walls (peek power and average power density for straight beam propagation in the empty chamber) and diffuse stray radiation effects are simulated to infer suitable strategies to avoid damage to first wall and to microwave sensitive components.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

Effects of the Spatial Extent of Multiple Harmonic Layers

An analytic model for single particle motion in the presence of a wave field and multiple cyclotron harmonics is developed and investigated. The model suggests that even in the absence of Doppler broadening, cyclotron harmonic layers have finite spatial extent. This allows for particles to interact with more than one harmonic layer simultaneously, provided the layers are tightly packed. The latter phenomenon is investigated in the context of the model using symplectic mapping techniques. Then the model behavior is compared with numerical simulations of neutral beam particle trajectories in NSTX using the full-orbit code SPIRAL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2012

Regions for Brillouin seed pulse growth in relativistic laser-plasma interaction

Parametric plasma processes received renewed interest in the context of generating ultraintense and ultrashort laser pulses up to the exawatt-zetawatt regime. For Brillouin scattering and seed pulse amplification at high intensities, the strong coupling regime is of special interest. The intensity-driven low-frequency modes depend on the amplitudes of the laser fields. It is investigated here how these modes develop in the relativistic regime. Then, a unified treatment of Raman and Brillouin processes becomes necessary. Assuming circular polarization, it is shown that with increasing intensity an overlap of the originally different Raman, Brillouin, and modulational instability branches occurs. Numerical simulations with a linearized Maxwell-fluid code confirm the analytically predicted behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Simulation of ion beam transport through the 400 Kv ion implanter at Michigan Ion Beam Laboratory

The Michigan Ion Beam Laboratory houses a 400 kV ion implanter. An application that simulates the ion beam trajectories through the implanter from the ion source to the target was developed using the SIMION Registered-Sign code. The goals were to have a tool to develop an intuitive understanding of abstract physics phenomena and diagnose ion trajectories. Using this application, new implanter users of different fields in science quickly understand how the machine works and quickly learn to operate it. In this article we describe the implanter simulation application and compare the parameters of the implanter components obtained from the simulations with the measured ones. The overall agreement between the simulated and measured values of magnetic fields and electric potentials is {approx}10%.

43 PARTICLE ACCELERATORS↗
DOE OSTI2013

Nonlinear interaction of kinetic Alfven wave and whistler: Turbulent spectra and anisotropic scaling

In this work, we are presenting the excitation of oblique propagating whistler wave as a consequence of nonlinear interaction between whistler wave and kinetic Alfven wave (KAW) in intermediate beta plasmas. Numerical simulation has been done to study the transient evolution of magnetic field structures of KAW when the nonlinearity arises due to ponderomotive effects by taking the adiabatic response of the background density. Weak oblique propagating whistler signals in these nonlinear plasma density filaments (produced by KAW localization) get amplified. The spectral indices of the power spectrum at different times are calculated with given initial conditions of the simulations. Anisotropic scaling laws for KAW and whistlers are presented. The relevance of the present investigation to solar wind turbulence and its acceleration is also pointed out.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

2-dimensional ion velocity distributions measured by laser-induced fluorescence above a radio-frequency biased silicon wafer

The dynamics of ions traversing sheaths in low temperature plasmas are important to the formation of the ion energy distribution incident onto surfaces during microelectronics fabrication. Ion dynamics have been measured using laser-induced fluorescence (LIF) in the sheath above a 30 cm diameter, 2.2 MHz-biased silicon wafer in a commercial inductively coupled plasma processing reactor. The velocity distribution of argon ions was measured at thousands of positions above and radially along the surface of the wafer by utilizing a planar laser sheet from a pulsed, tunable dye laser. Velocities were measured both parallel and perpendicular to the wafer over an energy range of 0.4–600 eV. The resulting fluorescence was recorded using a fast CCD camera, which provided resolution of 0.4 mm in space and 30 ns in time. Data were taken at eight different phases during the 2.2 MHz cycle. The ion velocity distributions (IVDs) in the sheath were found to be spatially non-uniform near the edge of the wafer and phase-dependent as a function of height. Several cm above the wafer the IVD is Maxwellian and independent of phase. Experimental results were compared with simulations. The experimental time-averaged ion energy distribution function as a function of height compare favorably with results from the computer model.

36 MATERIALS SCIENCE↗
DOE OSTI2013

GIOVE, a shallow laboratory Ge-spectrometer with 100 μBq/kg sensitivity

A new germanium gamma spectrometer called GIOVE (Germanium spectrometer with Inner and Outer Veto) has been set up at the underground/shallow laboratory (15 m w.e.) of MPI-K. Its double plastic scintillator veto system and neutron moderation interlayer lower the background by more than one order of magnitude compared to the other existing spectrometer at this facility. The integral (40-2700 keV) background rate of about 290 counts (day kg){sup −1} is just a factor 4 to 8 above that of the GeMPI spectrometers operated at LNGS (3800 m w.e.) and thus proves that even under shallow overburden sub mBq/kg sensitivities are achievable. Extended material screening and neutron attenuation studies preceded the final design of the spectrometer. The technical realization of the spectrometer is described in detail with special emphasis on the inner veto system. For its optimisation a simulation model was developed for light collection on small low activity PMT’s under various geometrical conditions. Radon suppression is accomplished by employing a gas tight sample container and a nitrogen flushed glove-box system with an airlock. The active volume of the crystal was modelled by absorption scanning measurements and Monte Carlo simulations. The complete shield is implemented in a Geant4 based simulation framework.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗