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,243 records · Page 69

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↗

Rayleigh-type parametric chemical oscillation

We consider a nonlinear chemical dynamical system of two phase space variables in a stable steady state. When the system is driven by a time-dependent sinusoidal forcing of a suitable scaling parameter at a frequency twice the output frequency and the strength of perturbation exceeds a threshold, the system undergoes sustained Rayleigh-type periodic oscillation, wellknown for parametric oscillation in pipe organs and distinct from the usual forced quasiperiodic oscillation of a damped nonlinear system where the system is oscillatory even in absence of any external forcing. Our theoretical analysis of the parametric chemical oscillation is corroborated by full numerical simulation of two well known models of chemical dynamics, chlorite-iodine-malonic acid and iodine-clock reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab initio implementation of quantum trajectory mean-field approach and dynamical simulation of the N{sub 2}CO photodissociation

In this work, the recently introduced quantum trajectory mean-field (QTMF) approach is implemented and employed to explore photodissociation dynamics of diazirinone (N{sub 2}CO), which are based on the high-level ab initio calculation. For comparison, the photodissociation process has been simulated as well with the fewest-switches surface hopping (FSSH) and the ab initio multiple spawning (AIMS) methods. Overall, the dynamical behavior predicted by the three methods is consistent. The N{sub 2}CO photodissociation at λ > 335 nm is an ultrafast process and the two C—N bonds are broken in a stepwise way, giving birth to CO and N{sub 2} as the final products in the ground state. Meanwhile, some noticeable differences were found in the QTMF, FSSH, and AIMS simulated time constants for fission of the C—N bonds, excited-state lifetime, and nonadiabatic transition ratios in different intersection regions. These have been discussed in detail. The present study provides a clear evidence that direct ab initio QTMF approach is one of the reliable tools for simulating nonadiabatic dynamics processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Research on calorimeter for high-power microwave measurements

Based on measurement of the volume increment of polar liquid that is a result of heating by absorbed microwave energy, two types of calorimeters with coaxial capacitive probes for measurement of high-power microwave energy are designed in this paper. The first is an “inline” calorimeter, which is placed as an absorbing load at the end of the output waveguide, and the second is an “offline” calorimeter that is placed 20 cm away from the radiation horn of the high-power microwave generator. Ethanol and high density polyethylene are used as the absorbing and housing materials, respectively. Results from both simulations and a “cold test” on a 9.3 GHz klystron show that the “inline” calorimeter has a measurement range of more than 100 J and an energy absorption coefficient of 93%, while the experimental results on a 9.3 GHz relativistic backward-wave oscillator show that the device’s power capacity is approximately 0.9 GW. The same experiments were also carried out for the “offline” calorimeter, and the results indicate that it can be used to eliminate the effects of the shock of the solenoid on the measurement curves and that the device has a higher power capacity of 2.5 GW. The results of the numerical simulations, the “cold tests,” and the experiments show good agreement.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Computational modeling of the effect of external electron injection into a direct-current microdischarge

Electron emission from the electrode surface plays an important role in determining the structure of a direct-current microdischarge. Here we have developed a computational model of a direct-current microdischarge to study the effect of external electron injection from the cathode surface into the discharge to manipulate its properties. The model provides a self-consistent, multi-species, multi-temperature fluid representation of the plasma. A microdischarge with a metal-insulator-metal configuration is chosen for this study. The effect of external electron injection on the structure and properties of the microdischarge is described. The transient behavior of the microdischarge during the electron injection is examined. The nonlinearities in the dynamics of the plasma result in a large increase of conduction current after active electron injection. For the conditions simulated a switching time of ∼100 ns from a low-current to high-current discharge state is realized.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reformulations of the Yang-Mills theory toward quark confinement and mass gap

We propose the reformulations of the SU (N) Yang-Mills theory toward quark confinement and mass gap. In fact, we have given a new framework for reformulating the SU (N) Yang-Mills theory using new field variables. This includes the preceding works given by Cho, Faddeev and Niemi, as a special case called the maximal option in our reformulations. The advantage of our reformulations is that the original non-Abelian gauge field variables can be changed into the new field variables such that one of them called the restricted field gives the dominant contribution to quark confinement in the gauge-independent way. Our reformulations can be combined with the SU (N) extension of the Diakonov-Petrov version of the non-Abelian Stokes theorem for the Wilson loop operator to give a gauge-invariant definition for the magnetic monopole in the SU (N) Yang-Mills theory without the scalar field. In the so-called minimal option, especially, the restricted field is non-Abelian and involves the non-Abelian magnetic monopole with the stability group U (N− 1). This suggests the non-Abelian dual superconductivity picture for quark confinement. This should be compared with the maximal option: the restricted field is Abelian and involves only the Abelian magnetic monopoles with the stability group U(1){sup N−1}, just like the Abelian projection. We give some applications of this reformulation, e.g., the stability for the homogeneous chromomagnetic condensation of the Savvidy type, the large N treatment for deriving the dimensional transmutation and understanding the mass gap, and also the numerical simulations on a lattice which are given by Dr. Shibata in a subsequent talk.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The propagation of particulate gravity currents in a V-shaped triangular cross section channel: Lock-release experiments and shallow-water numerical simulations

In this study, we investigate the motion of particulate gravity currents in a horizontal V-shaped channel. The particulate currents consisted of particles whose size varied between 0 and 100 μm but whose mean size increased. Particles were poorly sorted as the variance of the grain size distributions varied between 50 and 200. While the phases of propagation of homogeneous currents in such a geometry have been studied in the literature, this study considers the effects of the grain size on the propagation. The distance of propagation and front velocity of full-depth high-Reynolds-number lock-release experiments and shallow-water equation simulations were analyzed as the mean grain size of the initial particle distributions, defined by mass, was increased from 19 to 58 μm. Similar to the homogeneous currents, three consecutive phases of the front velocity could be identified but their characteristics and extent depend on the particle size. The initial phase, in particular, depends on a dimensionless settling number β that is defined as the ratio of two characteristic time scales, the propagation time x{sub 0}/U, where U is the scale for the front speed and x{sub 0} the lock length, and the settling time h{sub 0}/v{sub s}, where v{sub s} is the scale for the settling velocity and h{sub 0} the initial height of the current. For dimensionless settling numbers less than 0.001, the initial phase is characterized by a constant velocity for over about 6-7 lock lengths that is alike the initial slumping phase of perfectly constant velocity of the homogeneous currents. For dimensionless settling numbers greater than 0.001 and less than 0.015, the initial phase is no longer characterized by a constant velocity but an almost constant velocity for over about a similar 6-7 lock lengths. For dimensionless settling numbers greater than 0.015, however, as such, this phase is no longer seen. This initial phase is followed by a continuous decrease of the front advance, which results from the sedimentation of the particles. Unlike the homogeneous currents, this phase is a non-self-similar propagation. This phase is ended by a viscosity-dominated phase appearing to vary as ∼t{sup 1/7}. The good agreement between the front advance of the experiments and shallow-water equation simulations demonstrates that the mean size by mass is a fairly good proxy of poorly sorted particles.

42 ENGINEERING↗

Collisional transport across the magnetic field in drift-fluid models

Drift ordered fluid models are widely applied in studies of low-frequency turbulence in the edge and scrape-off layer regions of magnetically confined plasmas. Here, we show how collisional transport across the magnetic field is self-consistently incorporated into drift-fluid models without altering the drift-fluid energy integral. We demonstrate that the inclusion of collisional transport in drift-fluid models gives rise to diffusion of particle density, momentum, and pressures in drift-fluid turbulence models and, thereby, obviates the customary use of artificial diffusion in turbulence simulations. We further derive a computationally efficient, two-dimensional model, which can be time integrated for several turbulence de-correlation times using only limited computational resources. The model describes interchange turbulence in a two-dimensional plane perpendicular to the magnetic field located at the outboard midplane of a tokamak. The model domain has two regions modeling open and closed field lines. The model employs a computational expedient model for collisional transport. Numerical simulations show good agreement between the full and the simplified model for collisional transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

MHD simulation of the Bastille day event

We describe a time-dependent, thermodynamic, three-dimensional MHD simulation of the July 14, 2000 coronal mass ejection (CME) and flare. The simulation starts with a background corona developed using an MDI-derived magnetic map for the boundary condition. Flux ropes using the modified Titov-Demoulin (TDm) model are used to energize the pre-event active region, which is then destabilized by photospheric flows that cancel flux near the polarity inversion line. More than 10{sup 33} ergs are impulsively released in the simulated eruption, driving a CME at 1500 km/s, close to the observed speed of 1700km/s. The post-flare emission in the simulation is morphologically similar to the observed post-flare loops. The resulting flux rope that propagates to 1 AU is similar in character to the flux rope observed at 1 AU, but the simulated ICME center passes 15° north of Earth.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Calibration of imaging plates to electrons between 40 and 180 MeV

This paper presents the response calibration of Imaging Plates (IPs) for electrons in the 40-180 MeV range using laser-accelerated electrons at Laboratoire d’Optique Appliquée (LOA), Palaiseau, France. In the calibration process, the energy spectrum and charge of electron beams are measured by an independent system composed of a magnetic spectrometer and a Lanex scintillator screen used as a calibrated reference detector. It is possible to insert IPs of different types or stacks of IPs in this spectrometer in order to detect dispersed electrons simultaneously. The response values are inferred from the signal on the IPs, due to an appropriate charge calibration of the reference detector. The effect of thin layers of tungsten in front and/or behind IPs is studied in detail. GEANT4 simulations are used in order to analyze our measurements.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗