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,459 records · Page 81

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 ↗

Is the bulk mode conversion important in high density helicon plasma?

In a high-density helicon plasma production process, a contribution of Trivelpiece-Gould (TG) wave for surface power deposition is widely accepted. The TG wave can be excited either due to an abrupt density gradient near the plasma edge (surface conversion) or due to linear mode conversion from the helicon wave in a density gradient in the bulk region (bulk mode conversion). By numerically solving the boundary value problem of linear coupling between the helicon and the TG waves in a background with density gradient, we show that the efficiency of the bulk mode conversion strongly depends on the dissipation included in the plasma, and the bulk mode conversion is important when the dissipation is small. Also, by performing FDTD simulation, we show the time evolution of energy flux associated with the helicon and the TG waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A high-efficiency relativistic magnetron with the filled dielectric

Relativistic magnetron (RM) is a popular high power microwave device. Filling the resonant cavities with the low-loss dielectric is a well-proven technology that improves the performance of RM. In order to enhance the power efficiency, a dielectric-filled relativistic magnetron (DFRM) is presented and investigated numerically with KARAT code in this paper. The simulation results indicate that the maximum power efficiency is enhanced from 50.0% in the conventional RM to 56.0% in the DFRM when the applied voltage and the magnetic field are 280 kV and 0.5 T, respectively. Besides, the simulation results are analyzed and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Picosecond lasing in ytterbium fibre laser with nonlinear optical loop mirror: experiment and numerical simulation

The operation regimes of a pulsed all-normal-dispersion polarisation-maintaining fibre laser with a nonlinear optical loop mirror are studied. The use of polarisation-maintaining fibres ensures polarisation and temperature stability of output radiation. The lasing and instability thresholds of the pulsed laser are determined experimentally. A spectral filter placed in the cavity makes it possible to change the centre wavelength of laser radiation within the range 1.02 – 1.05 μm with a spectral full width at half maximum of 2 nm. The average output power is 7 mW, which corresponds to a pulse energy of 0.8 nJ. The autocorrelation function width of the output pulse is 50 ps. The minimum pulse duration achieved after compression by an external pair of diffraction gratings is 1.8 ps. The dynamics of the temporal and spectral parameters of laser pulses is studied using mathematical simulation based on numerical solution of the nonlinear Schrödinger equation. The simulation results coincide with experimental data with a high accuracy. (lasers)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Monitoring of chlorin-based photosensitiser localisation with dual-wavelength fluorescence imaging: numerical simulations

An approach to the estimation of chlorin-based photosensitiser localisation with dual-wavelength fluorescence imaging after topical application or intravenous injection is developed based on numerical simulations. Fluorescence imaging, which is commonly used for the monitoring of photosensitiser distribution within tissue, does not allow for its depth-resolved imaging. Chlorin-based photosensitisers feature two peaks in their absorption spectra corresponding to 402 and 662 nm, which allows for fluorescence excitation at two wavelengths with significantly different light penetration depths into biotissues. The ratio of fluorescence signals detected upon irradiation at different excitation wavelengths is shown to provide the information about photosensitiser localisation. The cases of exponential in-depth decay of the photosensitiser concentration and uniform photosensitiser distribution within the layer at a certain depth within the medium are analysed. These cases simulate topical application of photosensitiser and its intravenous injection and further accumulation in tumour, respectively. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A NEW STELLAR MIXING PROCESS OPERATING BELOW SHELL CONVECTION ZONES FOLLOWING OFF-CENTER IGNITION

During most stages of stellar evolution the nuclear burning of lighter to heavier elements results in a radial composition profile which is stabilizing against buoyant acceleration, with light material residing above heavier material. However, under some circumstances, such as off-center ignition, the composition profile resulting from nuclear burning can be destabilizing and characterized by an outwardly increasing mean molecular weight. The potential for instabilities under these circumstances and the consequences that they may have on stellar structural evolution remain largely unexplored. In this paper we study the development and evolution of instabilities associated with unstable composition gradients in regions that are initially stable according to linear Schwarzschild and Ledoux criteria. In particular, we study the development of turbulent flow under a variety of stellar evolution conditions with multi-dimensional hydrodynamic simulation; the phases studied include the core helium flash in a 1.25 M{sub Sun} star, the core carbon flash in a 9.3 M{sub Sun} star, and oxygen shell burning in a 23 M{sub Sun} star. The results of our simulations reveal a mixing process associated with regions having outwardly increasing mean molecular weight that reside below convection zones. The mixing is not due to overshooting from the convection zone, nor is it due directly to thermohaline mixing which operates on a timescale several orders of magnitude larger than the simulated flows. Instead, the mixing appears to be due to the presence of a wave field induced in the stable layers residing beneath the convection zone which enhances the mixing rate by many orders of magnitude and allows a thermohaline type mixing process to operate on a dynamical, rather than thermal, timescale. The mixing manifests itself in the form of overdense and cold blob-like structures originating from density fluctuations at the lower boundary of convective shell and 'shooting' down into the core. They are enriched with nuclearly processed material, hence leaving behind traces of higher mean molecular weight. In these regions, we find that initially smooth composition gradients steepen into stair-step-like profiles in which homogeneous, mixed regions are separated by composition jumps. These step-like profiles are then seen to evolve by a process of interface migration driven by turbulent entrainment. We discuss our results in terms of related laboratory phenomena and associated theoretical developments. We also discuss the degree to which the simulated mixing rates depend on the numerical resolution, and what future steps can be taken to capture the mixing rates accurately.

79 ASTRONOMY AND ASTROPHYSICS↗

THE CLUSTERING CHARACTERISTICS OF H I-SELECTED GALAXIES FROM THE 40% ALFALFA SURVEY

The 40% Arecibo Legacy Fast ALFA survey catalog ({alpha}.40) of {approx}10,150 H I-selected galaxies is used to analyze the clustering properties of gas-rich galaxies. By employing the Landy-Szalay estimator and a full covariance analysis for the two-point galaxy-galaxy correlation function, we obtain the real-space correlation function and model it as a power law, {xi}(r) = (r/r{sub 0}){sup -{gamma}}, on scales <10 h{sup -1} Mpc. As the largest sample of blindly H I-selected galaxies to date, {alpha}.40 provides detailed understanding of the clustering of this population. We find {gamma} = 1.51 {+-} 0.09 and r{sub 0} = 3.3 + 0.3, -0.2 h{sup -1} Mpc, reinforcing the understanding that gas-rich galaxies represent the most weakly clustered galaxy population known; we also observe a departure from a pure power-law shape at intermediate scales, as predicted in {Lambda}CDM halo occupation distribution models. Furthermore, we measure the bias parameter for the {alpha}.40 galaxy sample and find that H I galaxies are severely antibiased on small scales, but only weakly antibiased on large scales. The robust measurement of the correlation function for gas-rich galaxies obtained via the {alpha}.40 sample constrains models of the distribution of H I in simulated galaxies, and will be employed to better understand the role of gas in environmentally dependent galaxy evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

NEAR THE BOUNDARY OF THE HELIOSPHERE: A FLOW TRANSITION REGION

Since April of 2010, Voyager 1 has been immersed in a region of near zero radial flows, where the solar wind seems to have stopped. The existence of this region contradicts current models that predict that the radial flows will go to zero only at the heliopause. These models, however, do not include the sector region (or include it in a kinematic fashion), where the solar magnetic field periodically reverses polarity. Here we show that the presence of the sector region in the heliosheath, where reconnection occurs, fundamentally alters the flows, giving rise to a Flow Transition Region (FTR), where the flow abruptly turns and the radial velocity becomes near zero or negative. We estimate, based on a simulation, that at the Voyager 1 location, the thickness of the FTR is around 7-11 AU.

79 ASTRONOMY AND ASTROPHYSICS↗

TURBULENCE AND RADIO MINI-HALOS IN THE SLOSHING CORES OF GALAXY CLUSTERS

A number of relaxed, cool-core galaxy clusters exhibit diffuse, steep-spectrum radio sources in their central regions, known as radio mini-halos. It has been proposed that the relativistic electrons responsible for the emission have been reaccelerated by turbulence generated by the sloshing of the cool core gas. We present a high-resolution MHD simulation of gas sloshing in a galaxy cluster coupled with subgrid simulations of relativistic electron acceleration to test this hypothesis. Our simulation shows that the sloshing motions generate turbulence on the order of {delta}v {approx} 50-200 km s{sup -1} on spatial scales of {approx}50-100 kpc and below in the cool core region within the envelope of the sloshing cold fronts, whereas outside the cold fronts, there is negligible turbulence. This turbulence is potentially strong enough to reaccelerate relativistic electron seeds (with initial {gamma} {approx} 100-500) to {gamma} {approx} 10{sup 4} via damping of magnetosonic waves and non-resonant compression. The seed electrons could remain in the cluster from, e.g., past active galactic nucleus activity. In combination with the magnetic field amplification in the core, these electrons then produce diffuse radio synchrotron emission that is coincident with the region bounded by the sloshing cold fronts, as indeed observed in X-rays and the radio. The result holds for different initial spatial distributions of pre-existing relativistic electrons. The power and the steep spectral index ({alpha} Almost-Equal-To 1-2) of the resulting radio emission are consistent with observations of mini-halos, though the theoretical uncertainties of the acceleration mechanisms are high. We also produce simulated maps of inverse-Compton hard X-ray emission from the same population of relativistic electrons.

79 ASTRONOMY AND ASTROPHYSICS↗

PHOTOSPHERIC EMISSION AS THE DOMINANT RADIATION MECHANISM IN LONG-DURATION GAMMA-RAY BURSTS

We present the results of a set of numerical simulations of long-duration gamma-ray burst jets associated with massive, compact stellar progenitors. The simulations extend to large radii and allow us to locate the region in which the peak frequency of the advected radiation is set before the radiation is released at the photosphere. Light curves and spectra are calculated for different viewing angles as well as different progenitor structures and jet properties. We find that the radiation released at the photosphere of matter-dominated jets is able to reproduce the observed Amati and energy-Lorentz factor correlations. Our simulations also predict a correlation between the burst energy and the radiative efficiency of the prompt phase, consistent with observations.

79 ASTRONOMY AND ASTROPHYSICS↗