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,441 records · Page 80

Synergistic Effect of Nitrogen and Refractory Material on TiN Formation and Equiaxed Grain Structure of Ferritic Stainless Steel

The effect of nitrogen content on the formation of an equiaxed solidification structure of Fe-16Cr steel was investigated. Moreover, two different kinds of refractory materials, i.e., alumina and magnesia, were employed to control the type of oxide inclusion. The characteristics of TiN(-oxide) inclusions were quantitatively analyzed in both molten steel and solidified samples. When the melting was carried out in the alumina refractory, the grain size continuously decreased with increasing nitrogen content. However, a minimum grain size was observed at a specific nitrogen content (approx. 150 ppm) when the steel was melted in the magnesia refractory. Most of the single TiN particles had a cuboidal shape and fine irregularly shaped particles were located along the grain boundary due to the microsegregation of Ti at the grain boundary during solidification. The type of TiN-oxide hybrid inclusion was strongly affected by the refractory material where Al{sub 2}O{sub 3}-TiN and MgAl{sub 2}O{sub 4}-TiN hybrid-type inclusions were obtained in the alumina and magnesia refractory experiments, respectively. The formation of oxide inclusions was well predicted by thermochemical computations and it was commonly found that oxide particles were initially formed, followed by the nucleation and growth of TiN. When the nitrogen content increased, the number density of TiN linearly increased in the alumina refractory experiments. However, the number of TiN exhibits a maximum at about [N] = 150 ppm, at which a minimum grain size was obtained in the magnesia refractory experiments. Therefore, the larger the number density of TiN, the smaller the primary grain size after solidification. The number density of TiN in the steel melted in the magnesia refractory was greater than that in the steel melted in the alumina refractory at given Ti and N contents, which was due to the lower planar lattice disregistry of MgAl{sub 2}O{sub 4}-TiN interface rather than that of Al{sub 2}O{sub 3}-TiN interface. When ∆T{sub TiN} (= difference between the TiN precipitation temperature and the liquidus of the steel) was 20 K to 40 K, the number density of effective TiN was maximized and thus, the grain size was minimized after solidification. Finally, although most of the TiN particles were smaller than 1 μm in the molten steel samples irrespective of the nitrogen content, TiN particles larger than 10 μm were observed in the solidified samples when the nitrogen content was greater than 150 ppm. The growth of TiN particles during melting and solidification was well predicted by the combinatorial simulation of the ‘Ostwald ripening model’ based on the Lifshitz–Slyozov–Wagner theory in conjunction with the ‘Diffusion controlled model’ using Ohnaka’s microsegregation equation.

36 MATERIALS SCIENCE↗

Removal of Bifilms from Al Melts by Stirring in Unbaffled and Baffled Crucibles

The free surface of aluminum melt, during stirring under different rotation speeds in unbaffled and baffled crucibles, was simulated numerically using a RANS turbulence model coupled with a particle-tracking method. The effect of mechanical stirring of Al melts containing 0.3, 0.7, and 4.5 wt pct Mg in unbaffled and baffled crucibles on the rate of elimination of bifilms from the melt was also investigated using a reduced pressure test. The results showed that the stirring increased the rate of removal of bifilms from the melt due to an increase in the rate of formation of cracks on the bifilms. Using baffles in the crucible accelerated the rate of removal of bifilms by changing the pattern of movement of bifilms in the melt. The critical stirring speeds, above which the oxide layer on the Al melt could submerge into the bulk liquid, were determined by the simulation. It was confirmed by the experimental results to be about 158 and 420 RPM for the unbaffled and baffled crucibles, respectively.

36 MATERIALS SCIENCE↗

A proposal for the realization of universal quantum gates via superconducting qubits inside a cavity

A family of quantum logic gates is proposed via superconducting (SC) qubits coupled to a SC-cavity. The Hamiltonian for SC-charge qubits inside a single mode cavity is considered. Three- and two-qubit operations are generated by applying a classical magnetic field with the flux. Therefore, a number of quantum logic gates are realized. Numerical simulations and calculation of the fidelity are used to prove the success of these operations for these gates. -- Highlights: •A family of quantum logic gates is proposed via SC-qubits coupled to a cavity. •Three- and two-qubit operations are generated via a classical field with the flux. •Numerical simulations and calculation of the fidelity are used to prove the success of these operations for these gates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A mesh partitioning algorithm for preserving spatial locality in arbitrary geometries

Highlights: •An algorithm for partitioning computational meshes is proposed. •The Morton order space-filling curve is modified to achieve improved locality. •A spatial locality metric is defined to compare results with existing approaches. •Results indicate improved performance of the algorithm in complex geometries. -- Abstract: A space-filling curve (SFC) is a proximity preserving linear mapping of any multi-dimensional space and is widely used as a clustering tool. Equi-sized partitioning of an SFC ignores the loss in clustering quality that occurs due to inaccuracies in the mapping. Often, this results in poor locality within partitions, especially for the conceptually simple, Morton order curves. We present a heuristic that improves partition locality in arbitrary geometries by slicing a Morton order curve at points where spatial locality is sacrificed. In addition, we develop algorithms that evenly distribute points to the extent possible while maintaining spatial locality. A metric is defined to estimate relative inter-partition contact as an indicator of communication in parallel computing architectures. Domain partitioning tests have been conducted on geometries relevant to turbulent reactive flow simulations. The results obtained highlight the performance of our method as an unsupervised and computationally inexpensive domain partitioning tool.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Monte Carlo simulation for the electron cascade due to gamma rays in semiconductor radiation detectors

A Monte Carlo code was developed for simulating the electron cascade in radiation detector materials. The electron differential scattering cross sections were derived from measured electron energy-loss and optical spectra, making the method applicable for a wide range of materials. The detector resolution in a simplified model system shows dependence on the bandgap, the plasmon strength and energy, and the valence band width. In principle, these parameters could be optimized to improve detector performance. The intrinsic energy resolution was calculated for three semiconductors: silicon (Si), gallium arsenide (GaAs), and zinc telluride (ZnTe). Setting the ionization thresholds for electrons and holes is identified as a critical issue, as this strongly affects both the average electron-hole pair energy w and the Fano factor F. Using an ionization threshold from impact ionization calculations as an effective bandgap yields pair energies that are well matched to measured values. Fano factors of 0.091 (Si), 0.100 (GaAs), and 0.075 (ZnTe) were calculated. The Fano factor calculated for silicon using this model was lower than some results from past simulations and experiments. This difference could be attributed to problems in simulating inter-band transitions and the scattering of low-energy electrons.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Nonadiabatic electron response in the Hasegawa-Wakatani equations

Tokamak edge turbulence is strongly influenced by parallel electron physics, which relaxes density and potential fluctuations towards electron adiabatic response. Beginning with the paradigmatic Hasegawa-Wakatani equations (HWEs) for resistive tokamak edge turbulence, a unique decomposition of the electric potential (φ) into adiabatic (a) and nonadiabatic (b) portions is derived, based on the requirement that a neither drive nor respond to the parallel current j{sub ∥}. The form of the decomposition clarifies that, at perpendicular scales large relative to the sound radius, the electron adiabatic response controls the nonzonal φ, not the fluctuating density n. Simple energy balance arguments allow one to rigorously bound the ratio of rms nonzonal nonadiabatic fluctuations (b(tilde sign)) relative to adiabatic ones (ã). The role of the vorticity nonlinearity in transferring energy between adiabatic and nonadiabatic fluctuations aids intuitive understanding of self-sustained turbulence in the HWEs. When the normalized parallel resistivity is weak, b(tilde sign) becomes effectively slaved, allowing the reduction to an approximate one-field model that remains valid for strong turbulence. In addition to guiding physical intuition, the one-field reduction should greatly ease further analytical manipulations. Direct numerical simulation of the 2D HWEs confirms the convergence of the asymptotic formula for b(tilde sign)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A comprehensive comparison for simulations of cosmic-ray muons underground

The two leading simulation frameworks used for the simulation of cosmic-ray muons underground are FLUKA and Geant4. There have been in the past various questions raised as to the equivalence of these codes regarding cosmogenically produced neutrons and radioactivity in an underground environment. Many experiments choose one of these frameworks, and because they typically have different geometries or locations, the issues relating to code comparison are compounded. We report on an effort to compare the results of each of these codes in simulations which have simple geometry that is consistent between the two codes. It is seen that in terms of integrated neutron flux and neturon capture statistics the codes agree well in a broad sense. There are, however, differences that will be subject of further study. Comparisons of the simulations to available data are considered and the difficulties of such comparisons are pointed out.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Spiral-like multi-beam emission via transformation electromagnetics

Transformation electromagnetics offers an unconventional approach for the design of novel radiating devices. Here, we propose an electromagnetic metamaterial able to split an isotropic radiation into multiple directive beams. By applying transformations that modify distance and angles, we show how the multiple directive beams can be steered at will. We describe transformation of the metric space and the calculation of the material parameters. Different transformations are proposed for a possible physical realization through the use of engineered artificial metamaterials. Full wave simulations are performed to validate the proposed approach. The idea paves the way to interesting applications in various domains in microwave and optical regimes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Self-generation of dissipative solitons in magnonic quasicrystal active ring resonator

Self-generation of dissipative solitons in the magnonic quasicrystal (MQC) active ring resonator is studied theoretically and experimentally. The developed magnonic crystal has quasiperiodic Fibonacci type structure. Frequency selectivity of the MQC together with the parametric three-wave decay of magnetostatic surface spin wave (MSSW) leads to the dissipative soliton self-generation. The transfer matrix method is used to describe MQC transmission responses. Besides, the model of MQC active ring resonator is suggested. The model includes three coupled differential equations describing the parametric decay of MSSW and two differential equations of linear oscillators describing the frequency selectivity of MQC. Numerical simulation results of dissipative soliton self-generation are in a fair agreement with experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Pressure distribution based optimization of phase-coded acoustical vortices

Based on the acoustic radiation of point source, the physical mechanism of phase-coded acoustical vortices is investigated with formulae derivations of acoustic pressure and vibration velocity. Various factors that affect the optimization of acoustical vortices are analyzed. Numerical simulations of the axial, radial, and circular pressure distributions are performed with different source numbers, frequencies, and axial distances. The results prove that the acoustic pressure of acoustical vortices is linearly proportional to the source number, and lower fluctuations of circular pressure distributions can be produced for more sources. With the increase of source frequency, the acoustic pressure of acoustical vortices increases accordingly with decreased vortex radius. Meanwhile, increased vortex radius with reduced acoustic pressure is also achieved for longer axial distance. With the 6-source experimental system, circular and radial pressure distributions at various frequencies and axial distances have been measured, which have good agreements with the results of numerical simulations. The favorable results of acoustic pressure distributions provide theoretical basis for further studies of acoustical vortices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A methodology to identify and quantify mobility-reducing defects in 4H-silicon carbide power metal-oxide-semiconductor field-effect transistors

In this paper, we present a methodology for the identification and quantification of defects responsible for low channel mobility in 4H-Silicon Carbide (SiC) power metal-oxide-semiconductor field-effect transistors (MOSFETs). To achieve this, we use an algorithm based on 2D-device simulations of a power MOSFET, density functional simulations, and measurement data. Using physical modeling of carrier mobility and interface traps, we reproduce the experimental I-V characteristics of a 4H-SiC doubly implanted MOSFET through drift-diffusion simulation. We extract the position of Fermi level and the occupied trap density as a function of applied bias and temperature. Using these inputs, our algorithm estimates the number of possible trap types, their energy levels, and concentrations at 4H-SiC/SiO{sub 2} interface. Subsequently, we use density functional theory (DFT)-based ab initio simulations to identify the atomic make-up of defects causing these trap levels. We study silicon vacancy and carbon di-interstitial defects in the SiC side of the interface. Our algorithm indicates that the D{sub it} spectrum near the conduction band edge (3.25 eV) is composed of three trap types located at 2.8–2.85 eV, 3.05 eV, and 3.1–3.2 eV, and also calculates their densities. Based on DFT simulations, this work attributes the trap levels very close to the conduction band edge to the C di-interstitial defect.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Drifting solutions with elliptic symmetry for the compressible Navier-Stokes equations with density-dependent viscosity

In this paper, we investigate the analytical solutions of the compressible Navier-Stokes equations with dependent-density viscosity. By using the characteristic method, we successfully obtain a class of drifting solutions with elliptic symmetry for the Navier-Stokes model wherein the velocity components are governed by a generalized Emden dynamical system. In particular, when the viscosity variables are taken the same as Yuen [M. W. Yuen, “Analytical solutions to the Navier-Stokes equations,” J. Math. Phys. 49, 113102 (2008)], our solutions constitute a generalization of that obtained by Yuen. Interestingly, numerical simulations show that the analytical solutions can be used to explain the drifting phenomena of the propagation wave like Tsunamis in oceans.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bias-field controlled phasing and power combination of gyromagnetic nonlinear transmission lines

Gyromagnetic Nonlinear Transmission Lines (NLTLs) generate microwaves through the damped gyromagnetic precession of the magnetic moments in ferrimagnetic material, and are thus utilized as compact, solid-state, frequency agile, high power microwave (HPM) sources. The output frequency of a NLTL can be adjusted by control of the externally applied bias field and incident voltage pulse without physical alteration to the structure of the device. This property provides a frequency tuning capability not seen in many conventional e-beam based HPM sources. The NLTLs developed and tested are mesoband sources capable of generating MW power levels in the L, S, and C bands of the microwave spectrum. For an individual NLTL the output power at a given frequency is determined by several factors including the intrinsic properties of the ferrimagnetic material and the transmission line structure. Hence, if higher power levels are to be achieved, it is necessary to combine the outputs of multiple NLTLs. This can be accomplished in free space using antennas or in a transmission line via a power combiner. Using a bias-field controlled delay, a transient, high voltage, coaxial, three port, power combiner was designed and tested. Experimental results are compared with the results of a transient COMSOL simulation to evaluate combiner performance.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Theoretical and numerical predictions of hypervelocity impact-generated plasma

The hypervelocity impact generated plasmas (HVIGP) in thermodynamic non-equilibrium state were theoretically analyzed, and a physical model was presented to explore the relationship between plasma ionization degree and internal energy of the system by a group of equations including a chemical reaction equilibrium equation, a chemical reaction rate equation, and an energy conservation equation. A series of AUTODYN 3D (a widely used software in dynamic numerical simulations and developed by Century Dynamic Inc.) numerical simulations of the impacts of hypervelocity Al projectile on its targets at different incident angles were performed. The internal energy and the material density obtained from the numerical simulations were then used to calculate the ionization degree and the electron temperature. Based on a self-developed 2D smooth particle hydrodynamic (SPH) code and the theoretical model, the plasmas generated by 6 hypervelocity impacts were directly simulated and their total charges were calculated. The numerical results are in good agreements with the experimental results as well as the empirical formulas, demonstrating that the theoretical model is justified by the AUTODYN 3D and self-developed 2D SPH simulations and applicable to predict HVIGPs. The study is of significance for astrophysical and cosmonautic researches and safety.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantum interference of independently generated telecom-band single photons

We report on high-visibility quantum interference of independently generated telecom O-band (1310 nm) single photons using standard single-mode fibers. The experimental data are shown to agree well with the results of simulations using a comprehensive quantum multimode theory without the need for any fitting parameter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Class of backpropagation techniques for limited-angle reconstruction in microwave tomography

Filtered backpropagation (FBPP) is a well-known technique used in Diffraction Tomography (DT). For accurate reconstruction using FBPP, full 360° angular coverage is necessary. However, it has been shown that using some inherent redundancies in the projection data in a tomographic setup, accurate reconstruction is still possible with 270° coverage which is called the minimal-scan angle range. This can be done by applying weighing functions (or filters) on projection data of the object to eliminate the redundancies and accurately reconstruct the image from 270° coverage. This paper demonstrates procedures to generate many general classes of these weighing filters. These are all equivalent at 270° coverage but vary in performance at lower angular coverages and in presence of noise. This paper does a comparative analysis of different filters when angular coverage is lower than minimal-scan angle of 270°. Simulation studies have been done to find optimum weight filters for sub-minimal angular coverage (<270°)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Optical coherence tomography for nondestructive evaluation of fuel rod degradation

Nuclear power plants regularly inspect fuel rods to ensure safe and reliable operation. Excessive corrosion can cause fuel failures which can have significant repercussions for the plant, including impacts on plant operation, worker exposure to radiation, and the plant's INPO rating. While plants typically inspect for fuel rod corrosion using eddy current techniques, these techniques have known issues with reliability in the presence of tenacious, ferromagnetic crud layers that can deposit during operation, and the nondestructive evaluation (NDE) inspection results can often be in error by a factor of 2 or 3. For this reason, alternative measurement techniques, such as Optical Coherence Tomography (OCT), have been evaluated that are not sensitive to the ferromagnetic nature of the crud. This paper demonstrates that OCT has significant potential to characterize the thickness of crud layers that can deposit on the surfaces of fuel rods during operation. Physical trials have been performed on simulated crud samples, and the resulting data show an apparent correlation between the crud layer thickness and the OCT signal.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Stochastic rotation dynamics for nematic liquid crystals

We introduce a new mesoscopic model for nematic liquid crystals (LCs). We extend the particle-based stochastic rotation dynamics method, which reproduces the Navier-Stokes equation, to anisotropic fluids by including a simplified Ericksen-Leslie formulation of nematodynamics. We verify the applicability of this hybrid model by studying the equilibrium isotropic-nematic phase transition and nonequilibrium problems, such as the dynamics of topological defects and the rheology of sheared LCs. Our simulation results show that this hybrid model captures many essential aspects of LC physics at the mesoscopic scale, while preserving microscopic thermal fluctuations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗