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 631 records · Page 35

From Weakly Chaotic Dynamics to Deterministic Subdiffusion via Copula Modeling

Copula modeling consists in finding a probabilistic distribution, called copula, whereby its coupling with the marginal distributions of a set of random variables produces their joint distribution. The present work aims to use this technique to connect the statistical distributions of weakly chaotic dynamics and deterministic subdiffusion. More precisely, we decompose the jumps distribution of Geisel–Thomae map into a bivariate one and determine the marginal and copula distributions respectively by infinite ergodic theory and statistical inference techniques. We verify therefore that the characteristic tail distribution of subdiffusion is an extreme value copula coupling Mittag–Leffler distributions. We also present a method to calculate the exact copula and joint distributions in the case where weakly chaotic dynamics and deterministic subdiffusion statistical distributions are already known. Numerical simulations and consistency with the dynamical aspects of the map support our results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Full quantum mechanical analysis of atomic three-grating Mach–Zehnder interferometry

Atomic three-grating Mach–Zehnder interferometry constitutes an important tool to probe fundamental aspects of the quantum theory. There is, however, a remarkable gap in the literature between the oversimplified models and robust numerical simulations considered to describe the corresponding experiments. Consequently, the former usually lead to paradoxical scenarios, such as the wave–particle dual behavior of atoms, while the latter make difficult the data analysis in simple terms. Here these issues are tackled by means of a simple grating working model consisting of evenly-spaced Gaussian slits. As is shown, this model suffices to explore and explain such experiments both analytically and numerically, giving a good account of the full atomic journey inside the interferometer, and hence contributing to make less mystic the physics involved. More specifically, it provides a clear and unambiguous picture of the wavefront splitting that takes place inside the interferometer, illustrating how the momentum along each emerging diffraction order is well defined even though the wave function itself still displays a rather complex shape. To this end, the local transverse momentum is also introduced in this context as a reliable analytical tool. The splitting, apart from being a key issue to understand atomic Mach–Zehnder interferometry, also demonstrates at a fundamental level how wave and particle aspects are always present in the experiment, without incurring in any contradiction or interpretive paradox. On the other hand, at a practical level, the generality and versatility of the model and methodology presented, makes them suitable to attack analogous problems in a simple manner after a convenient tuning. - Highlights: • A simple model is proposed to analyze experiments based on atomic Mach–Zehnder interferometry. • The model can be easily handled both analytically and computationally. • A theoretical analysis based on the combination of the position and momentum representations is considered. • Wave and particle aspects are shown to coexist within the same experiment, thus removing the old wave-corpuscle dichotomy. • A good agreement between numerical simulations and experimental data is found without appealing to best-fit procedures.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Stability of a modified Peaceman–Rachford method for the paraxial Helmholtz equation on adaptive grids

This study concerns the asymptotic stability of an eikonal, or ray, transformation based Peaceman–Rachford splitting method for solving the paraxial Helmholtz equation with high wave numbers. Arbitrary nonuniform grids are considered in transverse and beam propagation directions. The differential equation targeted has been used for modeling propagations of high intensity laser pulses over a long distance without diffractions. Self-focusing of high intensity beams may be balanced with the de-focusing effect of created ionized plasma channel in the situation, and applications of grid adaptations are frequently essential. It is shown rigorously that the fully discretized oscillation-free decomposition method on arbitrary adaptive grids is asymptotically stable with a stability index one. Simulation experiments are carried out to illustrate our concern and conclusions.

97 MATHEMATICS AND COMPUTING↗

A dynamical polynomial chaos approach for long-time evolution of SPDEs

We propose a Dynamical generalized Polynomial Chaos (DgPC) method to solve time-dependent stochastic partial differential equations (SPDEs) with white noise forcing. The long-time simulation of SPDE solutions by Polynomial Chaos (PC) methods is notoriously difficult as the dimension of the stochastic variables increases linearly with time. Exploiting the Markovian property of white noise, DgPC implements a restart procedure that allows us to expand solutions at future times in terms of orthogonal polynomials of the measure describing the solution at a given time and the future white noise. The dimension of the representation is kept minimal by application of a Karhunen–Loeve (KL) expansion. Using frequent restarts and low degree polynomials on sparse multi-index sets, the method allows us to perform long time simulations, including the calculation of invariant measures for systems which possess one. We apply the method to the numerical simulation of stochastic Burgers and Navier–Stokes equations with white noise forcing. Our method also allows us to incorporate time-independent random coefficients such as a random viscosity. We propose several numerical simulations and show that the algorithm compares favorably with standard Monte Carlo methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Imaging of free carriers in semiconductors via optical feedback in terahertz quantum cascade lasers

To monitor the density of photo-generated charge carriers on a semiconductor surface, we demonstrate a detectorless imaging system based on the analysis of the optical feedback in terahertz quantum cascade lasers. Photo-excited free electron carriers are created in high resistivity n-type silicon wafers via low power (≅40 mW/cm{sup 2}) continuous wave pump laser in the near infrared spectral range. A spatial light modulator allows to directly reconfigure and control the photo-patterned intensity and the associated free-carrier density distribution. The experimental results are in good agreement with the numerical simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Experimentally driven atomistic model of 1,2 polybutadiene

We present an efficient method of combining wide angle neutron scattering data with detailed atomistic models, allowing us to perform a quantitative and qualitative mapping of the organisation of the chain conformation in both glass and liquid phases. The structural refinement method presented in this work is based on the exploitation of the intrachain features of the diffraction pattern and its intimate linkage with atomistic models by the use of internal coordinates for bond lengths, valence angles, and torsion rotations. Atomic connectivity is defined through these coordinates that are in turn assigned by pre-defined probability distributions, thus allowing for the models in question to be built stochastically. Incremental variation of these coordinates allows for the construction of models that minimise the differences between the observed and calculated structure factors. We present a series of neutron scattering data of 1,2 polybutadiene at the region 120–400 K. Analysis of the experimental data yields bond lengths for Cî—¸C and C î—» C of 1.54 Å and 1.35 Å, respectively. Valence angles of the backbone were found to be at 112° and the torsion distributions are characterised by five rotational states, a three-fold trans-skew± for the backbone and gauche± for the vinyl group. Rotational states of the vinyl group were found to be equally populated, indicating a largely atactic chan. The two backbone torsion angles exhibit different behaviour with respect to temperature of their trans population, with one of them adopting an almost all trans sequence. Consequently, the resulting configuration leads to a rather persistent chain, something indicated by the value of the characteristic ratio extrapolated from the model. We compare our results with theoretical predictions, computer simulations, RIS models and previously reported experimental results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Comparison of bifurcation dynamics of turbulent transport models for the L-H transition

In more than three decades, a large amount of models and mechanisms have been proposed to describe a very beneficial feature of magnetically confined fusion plasmas: the L-H transition. Bifurcation theory can be used to compare these different models based on their dynamical transition structure. In this paper, we employ bifurcation theory to distinguish two fundamentally different descriptions of the interaction between turbulence levels and sheared flows. The analytic bifurcation analysis characterises the parameter space structure of the transition dynamics. Herewith, in these models three dynamically different types of transitions are characterised, sharp transitions, oscillatory transitions, and smooth transitions. One of the two models has a very robust transition structure and is therefore likely to be more accurate for such a robust phenomenon as the L-H transition. The other model needs more fine-tuning to get non-oscillatory transitions. These conclusions from the analytic bifurcation analysis are confirmed by dedicated numerical simulations, with the newly developed code Bifurcator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Determination of hydrogen/deuterium ratio with neutron measurements on MAST

On MAST, compressional Alfvén eigenmodes can be destabilized by the presence of a sufficiently large population of energetic particles in the plasma. This dependence was studied in a series of very similar discharges in which increasing amounts of hydrogen were puffed into a deuterium plasma. A simple method to estimate the isotopic ratio n{sub H}/n{sub D} using neutron emission measurements is here described. The inferred isotopic ratio ranged from 0.0 to 0.6 and no experimental indication of changes in radial profile of n{sub H}/n{sub D} were observed. These findings are confirmed by TRANSP/NUBEAM simulations of the neutron emission.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Non-contact direct measurement of the magnetocaloric effect in thin samples

An experimental setup, based on a non-contact temperature sensor, is proposed to directly measure the magnetocaloric effect of samples few micrometers thick. The measurement of the adiabatic temperature change of foils and ribbons is fundamental to design innovative devices based on magnetocaloric thin materials or micro-structuring bulk samples. The reliability of the proposed setup is demonstrated by comparing the measurements performed on a bulk gadolinium sample with the results obtained by an experimental setup based on a Cernox bare chip thermoresistance and by in-field differential scanning calorimetry. We show that this technique can measure the adiabatic temperature variation on gadolinium sheets as thin as 27 μm. Heat transfer simulations are added to describe the capability of the presented technique.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Suppressing light reflection from polycrystalline silicon thin films through surface texturing and silver nanostructures

This work demonstrates a novel method combining ion implantation and silver nanostructures for suppressing light reflection from polycrystalline silicon thin films. Samples were implanted with 20-keV hydrogen ions to a dose of 10¹⁷/cm², and some of them received an additional argon ion implant to a dose of 5×10¹⁵ /cm² at an energy between 30 and 300 keV. Compared to the case with a single H implant, the processing involved both H and Ar implants and post-implantation annealing has created a much higher degree of surface texturing, leading to a more dramatic reduction of light reflection from polycrystalline Si films over a broadband range between 300 and 1200 nm, e.g., optical reflection from the air/Si interface in the AM1.5 sunlight condition decreasing from ~30% with an untextured surface to below 5% for a highly textured surface after post-implantation annealing at 1000°C. Formation of Ag nanostructures on these ion beam processed surfaces further reduces light reflection, and surface texturing is expected to have the benefit of diminishing light absorption losses within large-size (>100 nm) Ag nanoparticles, yielding an increased light trapping efficiency within Si as opposed to the case with Ag nanostructures on a smooth surface. A discussion of the effects of surface textures and Ag nanoparticles on light trapping within Si thin films is also presented with the aid of computer simulations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Illusion thermodynamics: A camouflage technique changing an object into another one with arbitrary cross section

The previously reported magical thermal devices, such as the thermal invisible cloak and the thermal concentrator, are generalized into one general case named here thermal illusion device. The thermal illusion device is displayed by the design of a thermal reshaper which can reshape an arbitrary thermal object into another one with arbitrary cross section. General expressions of the material parameters for the thermal reshaper are derived unambiguously to greatly facilitate the design of general thermal illusion device. We believe that this work will broaden the current research and pave a path to the thermal invisibility. Numerical simulations show good agreement with the analytical results of the thermal illusion device.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Large amplitude oscillation of magnetization in spin-torque oscillator stabilized by field-like torque

Oscillation frequency of spin torque oscillator with a perpendicularly magnetized free layer and an in-plane magnetized pinned layer is theoretically investigated by taking into account the field-like torque. It is shown that the field-like torque plays an important role in finding the balance between the energy supplied by the spin torque and the dissipation due to the damping, which results in a steady precession. The validity of the developed theory is confirmed by performing numerical simulations based on the Landau-Lifshitz-Gilbert equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Exact milestoning

A new theory and an exact computer algorithm for calculating kinetics and thermodynamic properties of a particle system are described. The algorithm avoids trapping in metastable states, which are typical challenges for Molecular Dynamics (MD) simulations on rough energy landscapes. It is based on the division of the full space into Voronoi cells. Prior knowledge or coarse sampling of space points provides the centers of the Voronoi cells. Short time trajectories are computed between the boundaries of the cells that we call milestones and are used to determine fluxes at the milestones. The flux function, an essential component of the new theory, provides a complete description of the statistical mechanics of the system at the resolution of the milestones. We illustrate the accuracy and efficiency of the exact Milestoning approach by comparing numerical results obtained on a model system using exact Milestoning with the results of long trajectories and with a solution of the corresponding Fokker-Planck equation. The theory uses an equation that resembles the approximate Milestoning method that was introduced in 2004 [A. K. Faradjian and R. Elber, J. Chem. Phys. 120(23), 10880-10889 (2004)]. However, the current formulation is exact and is still significantly more efficient than straightforward MD simulations on the system studied.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Damage imaging using Lamb waves for SHM applications

2-D ultrasonic arrays, due to their beam-steering capability and all azimuth angle coverage are a very promising tool for the inspection of plate-like structures using Lamb waves (LW). Contrary to the classical linear phased arrays (PAs) the 2D arrays enable unequivocal defect localization and they are even capable of mode selectivity of the received LWs . Recently, it has been shown that multistatic synthetic focusing (SF) algorithms applied for 2D arrays are much more effective than the classical phase array mode commonly used in NDT. The multistatic SF assumes multiple transmissions of elements in a transmitting aperture and off-line processing of the data acquired by a receiving aperture. In the simplest implementation of the technique, only a single multiplexed input and a number of output channels are required, which results in significant hardware simplification compared with the PA systems. On the one hand implementation of the multistatic SF to 2D arrays creates additional degrees of freedom during the design of the array topology, which complicates the array design process. On the other hand, it enables designing sparse arrays with performance similar to that of the fully populated dense arrays. In this paper we present a general systematic approach to the design and optimization of imaging systems based on the 2D array operating in the multistatic mode. We start from presenting principles of the SF schemes applied to LW imaging. Then, we outline the coarray concept and demonstrate how it can be used for reducing number of elements of the 2D arrays. Finally, efficient tools for the investigation and experimental verification of the designed 2D array prototypes are presented. The first step in the investigation is theoretical evaluation performed using frequency-dependent structure transfer function (STF), which enables approximate simulation of an array excited with a tone-burst in a dispersive medium. Finally, we show how scanning laser vibrometer, sensing waves in multiple points corresponding to the locations of the 2D receiving array elements, can be used as a tool for rapid experimental verification of the developed topologies. The presented methods are discussed in terms of the beampatterns and sparse versions of the fully populated array topologies are be presented. The effect of apodization applied to the array elements is also investigated. Both simulated and experimental results are included.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Molecular dynamics simulations of soliton-like structures in a dusty plasma medium

The existence and evolution of soliton-like structures in a dusty plasma medium are investigated in a first principles approach using molecular dynamic (MD) simulations of particles interacting via a Yukawa potential. These localized structures are found to exist in both weakly and strongly coupled regimes with their structures becoming sharper as the correlation effects between the dust particles get stronger. A surprising result, compared to fluid simulations, is the existence of rarefactive soliton-like structures in our non-dissipative system, a feature that arises from the charge conjugation symmetry property of the Yukawa fluid. Our simulation findings closely resemble many diverse experimental results reported in the past.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation of the 2-dimensional Drude’s model using molecular dynamics method

In this paper, we reported the results of the simulation of the electronic conduction in solids. The simulation is based on the Drude’s models by applying molecular dynamics (MD) method, which uses the fifth-order predictor-corrector algorithm. A formula of the electrical conductivity as a function of lattice length and ion diameter τ(L, d) cand be obtained empirically based on the simulation results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Measurements of pion production in eA with the CLAS detector

Preliminary results on semi-inclusive charged pion production in eA collisions at E{sub beam}=5 GeV/c{sup 2} are presented. The data were collected using the CLAS detector, which is a multipurpose, large acceptance, magnetic spectrometer located in Hall B at the Thomas Jefferson National Accelerator Facility. Selected distributions in W, Q{sup 2}, p{sub π}, and θ{sub π} are shown for data produced using deuterium, carbon, iron and lead targets. Preliminary comparisons with data simulated using the GENIE generator are made. The motivation for this work is to provide distributions useful for tuning the hadronic final state interaction models in order to reduce uncertainties in the results extracted from current and next-generation neutrino oscillation experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Study of plastic strain localization mechanisms caused by nonequilibrium transitions in mesodefect ensembles under high-speed loading

The behavior of specimens dynamically loaded during split Hopkinson (Kolsky) bar tests in a regime close to simple shear conditions was studied. The lateral surface of the specimens was investigated in-situ using a high-speed infrared camera CEDIP Silver 450M. The temperature field distribution obtained at different time allowed one to trace the evolution of plastic strain localization. The process of target perforation involving plug formation and ejection was examined using a high-speed infrared camera and a VISAR velocity measurement system. The microstructure of tested specimens was analyzed using an optical interferometer-profiler and a scanning electron microscope. The development of plastic shear instability regions has been simulated numerically.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗