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,171 records · Page 65

All-optical XNOR gate based on 2D photonic-crystal ring resonators

A novel all-optical XNOR gate is proposed, which combines the nonlinear Kerr effect with photonic-crystal ring resonators (PCRRs). The total size of the proposed optical XNOR gate based on photonic crystals with a square lattice of silicon rods is equal to 35 × 21 μm. The proposed structure has a bandgap in the range from 0.32 to 0.44. To confirm the operation and feasibility of the overall system use is made of analytical and numerical simulation using the dimensional finite difference time domain (FDTD) and plane wave expansion (PWE) methods. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Contactless method for studying temperature within the active element of a multidisk cryogenic amplifier

A new original method has been developed and experimentally implemented, allowing temperature fields to be contactlessly measured in the pump region of active elements in high-power-diode-pumped laser amplifiers, including those operating at cryogenic temperatures. The presence of a temperature gradient of ∼57 K mm{sup −1} along the pump beam axis at the centre of the active element of the laser amplification unit operating at cryogenic temperatures with a pulse repetition rate up to 1 kHz is simulated and experimentally confirmed. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Formula for the ionisation rate of an atom or ion in a strong electromagnetic field for numerical simulation

A formula is proposed for calculating the rate of field ionisation of an atom or ion, taking into account both tunnel ionisation and barrier suppression ionisation. Compared with the previous formula proposed in 2018, it more accurately describes the transition region between both ionisations and is determined mainly by two parameters: the ionisation potential of an atom or ion and the amplitude of an external electric field. This makes the presented formula suitable for use in numerical packages simulating the interaction of high-power laser radiation with matter by the particle-in-cell (PIC) method. (paper)

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

PARTICLE-IN-CELL SIMULATIONS OF PARTICLE ENERGIZATION VIA SHOCK DRIFT ACCELERATION FROM LOW MACH NUMBER QUASI-PERPENDICULAR SHOCKS IN SOLAR FLARES

Low Mach number, high beta fast mode shocks can occur in the magnetic reconnection outflows of solar flares. These shocks, which occur above flare loop tops, may provide the electron energization responsible for some of the observed hard X-rays and contemporaneous radio emission. Here we present new two-dimensional particle-in-cell simulations of low Mach number/high beta quasi-perpendicular shocks. The simulations show that electrons above a certain energy threshold experience shock-drift-acceleration. The transition energy between the thermal and non-thermal spectrum and the spectral index from the simulations are consistent with some of the X-ray spectra from RHESSI in the energy regime of E {approx}< 40 {approx} 100 keV. Plasma instabilities associated with the shock structure such as the modified-two-stream and the electron whistler instabilities are identified using numerical solutions of the kinetic dispersion relations. We also show that the results from PIC simulations with reduced ion/electron mass ratio can be scaled to those with the realistic mass ratio.

79 ASTRONOMY AND ASTROPHYSICS↗

MAGNETOHYDRODYNAMIC SIMULATIONS OF INTERPLANETARY CORONAL MASS EJECTIONS

We describe a new MHD model for the propagation of interplanetary coronal mass ejections (ICMEs) in the solar wind. Accurately following the propagation of ICMEs is important for determining space weather conditions. Our model solves the MHD equations in spherical coordinates from a lower boundary above the critical point to Earth and beyond. On this spherical surface, we prescribe the magnetic field, velocity, density, and temperature calculated typically directly from a coronal MHD model as time-dependent boundary conditions. However, any model that can provide such quantities either in the inertial or rotating frame of the Sun is suitable. We present two validations of the technique employed in our new model and a more realistic simulation of the propagation of an ICME from the Sun to Earth.

79 ASTRONOMY AND ASTROPHYSICS↗

Statistical mechanics of collisionless orbits. IV. Distribution of angular momentum

It has been shown in previous work that DARKexp, which is a theoretically derived, maximum entropy, one shape parameter model for isotropic collisionless systems, provides very good fits to simulated and observed dark matter halos. Specifically, it fits the energy distribution, N(E), and the density profiles, including the central cusp. Here, we extend DARKexp N(E) to include the distribution in angular momentum, L {sup 2}, for spherically symmetric systems. First, we argue, based on theoretical, semi-analytical, and simulation results, that while dark matter halos are relaxed in energy, they are not nearly as relaxed in angular momentum, which precludes using maximum entropy to uniquely derive N(E, L {sup 2}). Instead, we require that when integrating N(E, L {sup 2}) over squared angular momenta one retrieves the DARKexp N(E). Starting with a general expression for N(E, L {sup 2}) we show how the distribution of particles in L {sup 2} is related to the shape of the velocity distribution function, VDF, and velocity anisotropy profile, β(r). We then demonstrate that astrophysically realistic halos, as judged by the VDF shape and β(r), must have linear or convex distributions in L {sup 2}, for each separate energy bin. The distribution in energy of the most bound particles must be nearly flat, and become more tilted in favor of radial orbits for less bound particles. These results are consistent with numerical simulations and represent an important step toward deriving the full distribution function for spherically symmetric dark matter halos.

79 ASTRONOMY AND ASTROPHYSICS↗

Three-dimensional MHD simulation of the Caltech plasma jet experiment: first results

Magnetic fields are believed to play an essential role in astrophysical jets with observations suggesting the presence of helical magnetic fields. Here, we present three-dimensional (3D) ideal MHD simulations of the Caltech plasma jet experiment using a magnetic tower scenario as the baseline model. Magnetic fields consist of an initially localized dipole-like poloidal component and a toroidal component that is continuously being injected into the domain. This flux injection mimics the poloidal currents driven by the anode-cathode voltage drop in the experiment. The injected toroidal field stretches the poloidal fields to large distances, while forming a collimated jet along with several other key features. Detailed comparisons between 3D MHD simulations and experimental measurements provide a comprehensive description of the interplay among magnetic force, pressure, and flow effects. In particular, we delineate both the jet structure and the transition process that converts the injected magnetic energy to other forms. With suitably chosen parameters that are derived from experiments, the jet in the simulation agrees quantitatively with the experimental jet in terms of magnetic/kinetic/inertial energy, total poloidal current, voltage, jet radius, and jet propagation velocity. Specifically, the jet velocity in the simulation is proportional to the poloidal current divided by the square root of the jet density, in agreement with both the experiment and analytical theory. This work provides a new and quantitative method for relating experiments, numerical simulations, and astrophysical observation, and demonstrates the possibility of using terrestrial laboratory experiments to study astrophysical jets.

79 ASTRONOMY AND ASTROPHYSICS↗

OPACITY BROADENING OF {sup 13}CO LINEWIDTHS AND ITS EFFECT ON THE VARIANCE-SONIC MACH NUMBER RELATION

We study how the estimation of the sonic Mach number (M{sub s} ) from {sup 13}CO linewidths relates to the actual three-dimensional sonic Mach number. For this purpose we analyze MHD simulations that include post-processing to take radiative transfer effects into account. As expected, we find very good agreement between the linewidth estimated sonic Mach number and the actual sonic Mach number of the simulations for optically thin tracers. However, we find that opacity broadening causes M{sub s} to be overestimated by a factor of ≈1.16-1.3 when calculated from optically thick {sup 13}CO lines. We also find that there is a dependence on the magnetic field: super-Alfvénic turbulence shows increased line broadening compared with sub-Alfvénic turbulence for all values of optical depth for supersonic turbulence. Our results have implications for the observationally derived sonic Mach number-density standard deviation (σ{sub ρ/(ρ)}) relationship, σ{sub ρ/〈ρ〉}{sup 2}=b{sup 2}M{sub s}{sup 2}, and the related column density standard deviation (σ {sub N/(N)}) sonic Mach number relationship. In particular, we find that the parameter b, as an indicator of solenoidal versus compressive driving, will be underestimated as a result of opacity broadening. We compare the σ {sub N/(N)}-M{sub s} relation derived from synthetic dust extinction maps and {sup 13}CO linewidths with recent observational studies and find that solenoidally driven MHD turbulence simulations have values of σ {sub N/(N)}which are lower than real molecular clouds. This may be due to the influence of self-gravity which should be included in simulations of molecular cloud dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗

Microstructure of vapor deposited coatings on curved substrates

Thermal barrier coating systems consisting of a metallic bond coat and ceramic over layer are widely used to extend the life of gas turbine engine components. They are applied using either high-vacuum physical vapor deposition techniques in which vapor atoms rarely experience scattering collisions during propagation to a substrate, or by gas jet assisted (low-vacuum) vapor deposition techniques that utilize scattering from streamlines to enable non-line-of-sight deposition. Both approaches require substrate motion to coat a substrate of complex shape. Here, direct simulation Monte Carlo and kinetic Monte Carlo simulation methods are combined to simulate the deposition of a nickel coating over the concave and convex surfaces of a model airfoil, and the simulation results are compared with those from experimental depositions. The simulation method successfully predicted variations in coating thickness, columnar growth angle, and porosity during both stationary and substrate rotated deposition. It was then used to investigate a wide range of vapor deposition conditions spanning high-vacuum physical vapor deposition to low-vacuum gas jet assisted vapor deposition. The average coating thickness was found to increase initially with gas pressure reaching a maximum at a chamber pressure of 8–10 Pa, but the best coating thickness uniformity was achieved under high vacuum deposition conditions. However, high vacuum conditions increased the variation in the coatings pore volume fraction over the surface of the airfoil. The simulation approach was combined with an optimization algorithm and used to investigate novel deposition concepts to tailor the local coating thickness.

42 ENGINEERING↗

SU-E-T-08: A Beam Source Model for Monte Carlo Simulations of a Double-Scattering Proton Beam Delivery System Using Beam Current Modulation

Purpose: To reconstruct phase-space information upstream of patient specific collimators for Monte Carlo simulations using only radiotherapy planning system data. Methods: The proton energies are calculated based on residual ranges, e.g., sum of prescribed ranges in a patient and SSD. The Kapchinskij and Vladimirskij (KV) distribution was applied to sample proton’s x-y positions and momentum direction and the beam shape was assumed to be a circle. Free parameters, e.g., the initial energy spread and the emittance of KV distribution were estimated from the benchmarking with commissioning data in a commercial treatment planning system for an operational proton therapy center. The number of histories, which defines the height of individual pristine Bragg peaks (BP) of Spread-out Bragg peak (SOBP), are weighted based on beam current modulation and a correction factor is applied to take into account the fluence reduction as the residual range decreases due to the rotation of the range modulator wheel. The timedependent behaviors, e.g., the changes of the residual range and histories per a pristine BP, are realized by utilizing TOPAS (Tool for Particle Simulation). Results: Benchmarking simulations for selected SOBPs ranging 7.5 cm to 15.5 cm matched within 2 mm in range and up to 5 mm in SOBP width against measurement data in water phantom. We found this model tends to underestimate entrance dose by about 5 % in comparison to measurement. This was attributed to the situation that the energy distribution used in the model was limited in its granularity at the limit of single energy spectrum for the narrow angle modulator steps used in the proximal pull back region of the SOBPs. Conclusion: Within these limitations the source modeling method proved itself an acceptable alternative of a full treatment head simulation when the machine geometry and materials information are not available.

60 APPLIED LIFE SCIENCES↗

SU-D-BRC-04: Development of Proton Tissue Equivalent Materials for Calibration and Dosimetry Studies

Purpose: To develop new proton tissue equivalent materials (PTEM), urethane and fiberglass based, for proton therapy calibration and dosimetry studies. Existing tissue equivalent plastics are applicable only for x-rays because they focus on matching mass attenuation coefficients. This study aims to create new plastics that match mass stopping powers for proton therapy applications instead. Methods: New PTEMs were constructed using urethane and fiberglass resin materials for soft, fat, bone, and lung tissue. The stoichiometric analysis method was first used to determine the elemental composition of each unknown constituent. New initial formulae were then developed for each of the 4 PTEMs using the new elemental compositions and various additives. Samples of each plastic were then created and exposed to a well defined proton beam at the UF Health Proton Therapy Institute (UFHPTI) to validate its mass stopping power. Results: The stoichiometric analysis method revealed the elemental composition of the 3 components used in creating the PTEMs. These urethane and fiberglass based resins were combined with additives such as calcium carbonate, aluminum hydroxide, and phenolic micro spheres to achieve the desired mass stopping powers and densities. Validation at the UFHPTI revealed adjustments had to be made to the formulae, but the plastics eventually had the desired properties after a few iterations. The mass stopping power, density, and Hounsfield Unit of each of the 4 PTEMs were within acceptable tolerances. Conclusion: Four proton tissue equivalent plastics were developed: soft, fat, bone, and lung tissue. These plastics match each of the corresponding tissue’s mass stopping power, density, and Hounsfield Unit, meaning they are truly tissue equivalent for proton therapy applications. They can now be used to calibrate proton therapy treatment planning systems, improve range uncertainties, validate proton therapy Monte Carlo simulations, and assess in-field and out-of-field organ doses.

60 APPLIED LIFE SCIENCES↗

SU-F-T-146: Comparing Monte Carlo Simulations with Commissioning Beam Data for Mevion S250 Proton Therapy System

Purpose: To evaluate the Monte Carlo simulated beam data with the measured commissioning data for the Mevion S250 proton therapy system. Method: The Mevion S250 proton therapy system utilizes a passive double scattering technique with a unique gantry mounted superconducting accelerator and offers effective proton therapy in a compact design concept. The field shaping system (FSS) includes first scattering foil, range modulator wheel (RMW), second scattering foil and post absorber and offers two field sizes and a total of 24 treatment options from proton range of 5 cm to 32 cm. The treatment nozzle was modeled in detail using TOPAS (TOolkit for PArticle Simulation) Monte Carlo code. The timing feathers of the moving modulator wheels were also implemented to generate the Spread Out Bragg Peak (SOBP). The simulation results including pristine Bragg Peak, SOBP and dose profiles were compared with the data measured during beam commissioning. Results: The comparison between the measured data and the simulation data show excellent agreement. For pristine proton Bragg Peaks, the simulated proton range (depth of distal 90%) values agreed well with the measured range values within 1 mm accuracy. The differences of the distal falloffs (depth from distal 80% to 20%) were also found to be less than 1 mm between the simulations and measurements. For the SOBP, the widths of modulation (depth of proximal 95% to distal 90%) were also found to agree with the measurement within 1 mm. The flatness of the simulated and measured lateral profiles was found to be 0.6 % and 1.1 %, respectively. Conclusion: The agreement between simulations and measurements demonstrate that TOPAS could be used as a viable platform to proton therapy applications. The matched simulation results offer a great tool and open opportunity for variety of applications.

60 APPLIED LIFE SCIENCES↗

WE-DE-207B-01: Optimization for Contrast-Enhanced Spectral Mammography Based On Photon-Counting Detectors

Purpose: To investigate the feasibility of optimizing the imaging parameters for contrast-enhanced spectral mammography based on Si strip photon-counting detectors. Methods: A computer simulation model using polyenergetic spectra from a tungsten anode x-ray tube and a Si-based photon-counting detector was evaluated for contrast-enhanced spectral mammography. The simulation traces the emission of photons from the x-ray source, attenuation through the breast and subsequent absorption in the detector. The breast was modeled as a mixture of adipose and mammary gland tissues with a breast density of 30%. A 4 mm iodine signal with a concentration of 4 mg/ml was used to simulate the enhancement of a lesion. Quantum efficiency of the detector was calculated based on the effective attenuation length in the Si strips. The figure-of-merit (FOM), which was defined as the decomposed iodine signal-to-noise ratio (SNR) with respect to the square root of the mean glandular dose (MGD), was chosen to optimize the imaging parameters, in terms of beam energy, splitting energy, and pre-filtrations for breast of various thicknesses and densities. Results: The optimal imaging parameters, which lead to the highest FOM, were found at a beam energy of 45 kVp with a splitting energy at 34 keV for an averaged breast thickness of 4 cm with a standard 0.75 mm Al pre-filtration. The optimal tube voltage varied slightly from 46 to 44 kVp as the breast thickness increases from 2 to 8 cm. The optimal tube voltage decreased to 42 kVp when the Al pre-filtration was increased to 3 mm. Conclusion: This simulation study predicted the optimal imaging parameters for application of photon-counting spectral mammography to contrast-enhanced imaging. The simulation results laid the ground work for future phantom and clinical studies. Grant funding from Philips Medical Systems.

60 APPLIED LIFE SCIENCES↗

Kink model of state switching in quasi-one-dimensional nanosystems

A method for the analytical calculation of the kinetics of state switching induced by external field in quasi-one-dimensional systems has been developed. This method is based on the Kolmogorov approach, which was proposed to describe statistical crystallization. As applied to magnetic systems, it allows one to calculate magnetization curves and hysteresis loops. The results of analytical calculation are confirmed by stochastic Monte Carlo simulation, which also makes it possible to find characteristics of mixed domain states formed as a result of the multiple cycling of the external field amplitude.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Double-Crystal Rocking Curve Simulation Using 2D Spectral Angular Diagrams of X-Ray Radiation

A new approach to numerical simulation of double-crystal rocking curves is proposed. This approach is based on the use of experimental spectral angular diagrams of X-ray intensity distribution. Special calculation algorithms, which take into account the instrumental function of X-ray diffractometer and possible effects of dispersion and Bragg reflection asymmetry, have been developed and applied. A specific feature of the proposed approach is the possibility of visualizing the 2D spectral angular diagram of X-ray beam after its interaction with each element of the scheme. The approach makes it possible to perform calculations for a wide range of radiation sources (from an X-ray tube with any anode to a synchrotron radiation source) and X-ray optical elements (slits and monochromators). A comparison of simulation results and experimental data for a Si(110) crystal sample has confirmed adequacy of the proposed approach and its applicability for simulating diffraction patterns recorded in real experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Accretion and orbital inspiral in gas-assisted supermassive black hole binary mergers

Many galaxies are expected to harbor binary supermassive black holes (SMBHs) in their centers. Their interaction with the surrounding gas results in the accretion and exchange of angular momentum via tidal torques, facilitating binary inspiral. Here, we explore the non-trivial coupling between these two processes and analyze how the global properties of externally supplied circumbinary disks depend on the binary accretion rate. By formulating our results in terms of the angular momentum flux driven by internal stresses, we come up with a very simple classification of the possible global disk structures, which differ from the standard constant M-dot accretion disk solution. The suppression of accretion by the binary tides, leading to a significant mass accumulation in the inner disk, accelerates binary inspiral. We show that once the disk region strongly perturbed by the viscously transmitted tidal torque exceeds the binary semimajor axis, the binary can merge in less than its mass-doubling time due to accretion. Thus, unlike the inspirals driven by stellar scattering, the gas-assisted merger can occur even if the binary is embedded in a relatively low-mass disk (lower than its own mass). This is important for resolving the “last parsec” problem for SMBH binaries and understanding powerful gravitational wave sources in the universe. We argue that the enhancement of accretion by the binary found in some recent simulations cannot persist for a long time and should not affect the long-term orbital inspiral. We also review existing simulations of SMBH binary–disk coupling and propose a numerical setup which is particularly well suited to verifying our theoretical predictions.

79 ASTRONOMY AND ASTROPHYSICS↗

Excitation and Depletion of the Asteroid Belt in the Early Instability Scenario

Containing only a few percentages of the mass of the moon, the current asteroid belt is around three to four orders of magnitude smaller than its primordial mass inferred from disk models. Yet dynamical studies have shown that the asteroid belt could not have been depleted by more than about an order of magnitude over the past ∼4 Gyr. The remainder of the mass loss must have taken place during an earlier phase of the solar system’s evolution. An orbital instability in the outer solar system occurring during the process of terrestrial planet formation can reproduce the broad characteristics of the inner solar system. Here, we test the viability of this model within the constraints of the main belt’s low present-day mass and orbital structure. Although previous studies modeled asteroids as massless test particles because of limited computing power, our work uses graphics processing unit acceleration to model a fully self-gravitating asteroid belt. We find that depletion in the main belt is related to the giant planets’ exact evolution within the orbital instability. Simulations that produce the closest matches to the giant planets’ current orbits deplete the main belt by two to three orders of magnitude. These simulated asteroid belts are also good matches to the actual asteroid belt in terms of their radial mixing and broad orbital structure.

79 ASTRONOMY AND ASTROPHYSICS↗