Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “COMPUTER 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,117 records · Page 62

LABORATORY MEASUREMENTS OF WHITE DWARF PHOTOSPHERIC SPECTRAL LINES: Hβ

We spectroscopically measure multiple hydrogen Balmer line profiles from laboratory plasmas to investigate the theoretical line profiles used in white dwarf (WD) atmosphere models. X-ray radiation produced at the Z Pulsed Power Facility at Sandia National Laboratories initiates plasma formation in a hydrogen-filled gas cell, replicating WD photospheric conditions. Here we present time-resolved measurements of Hβ and fit this line using different theoretical line profiles to diagnose electron density, n{sub e}, and n = 2 level population, n{sub 2}. Aided by synthetic tests, we characterize the validity of our diagnostic method for this experimental platform. During a single experiment, we infer a continuous range of electron densities increasing from n{sub e} ∼ 4 to ∼30 × 10{sup 16} cm{sup −3} throughout a 120-ns evolution of our plasma. Also, we observe n{sub 2} to be initially elevated with respect to local thermodynamic equilibrium (LTE); it then equilibrates within ∼55 ns to become consistent with LTE. This supports our electron-temperature determination of T{sub e} ∼ 1.3 eV (∼15,000 K) after this time. At n{sub e} ≳ 10{sup 17} cm{sup −3}, we find that computer-simulation-based line-profile calculations provide better fits (lower reduced χ{sup 2}) than the line profiles currently used in the WD astronomy community. The inferred conditions, however, are in good quantitative agreement. This work establishes an experimental foundation for the future investigation of relative shapes and strengths between different hydrogen Balmer lines.

79 ASTRONOMY AND ASTROPHYSICS↗

Holographic Imaging of Evolving Laser-Plasma Structures

In the 1870s, English photographer Eadweard Muybridge captured motion pictures within one cycle of a horse’s gallop, which settled a hotly debated question of his time by showing that the horse became temporarily airborne. In the 1940s, Manhattan project photographer Berlin Brixner captured a nuclear blast at a million frames per second, and resolved a dispute about the explosion’s shape and speed. In this project, we developed methods to capture detailed motion pictures of evolving, light-velocity objects created by a laser pulse propagating through matter. These objects include electron density waves used to accelerate charged particles, laser-induced refractive index changes used for micromachining, and ionization tracks used for atmospheric chemical analysis, guide star creation and ranging. Our “movies”, like Muybridge’s and Brixner’s, are obtained in one shot, since the laser-created objects of interest are insufficiently repeatable for accurate stroboscopic imaging. Our high-speed photographs have begun to resolve controversies about how laser-created objects form and evolve, questions that previously could be addressed only by intensive computer simulations based on estimated initial conditions. Resolving such questions helps develop better tabletop particle accelerators, atmospheric ranging devices and many other applications of laser-matter interactions. Our photographic methods all begin by splitting one or more “probe” pulses from the laser pulse that creates the light-speed object. A probe illuminates the object and obtains information about its structure without altering it. We developed three single-shot visualization methods that differ in how the probes interact with the object of interest or are recorded. (1) Frequency-Domain Holography (FDH). In FDH, there are 2 probes, like “object” and “reference” beams in conventional holography. Our “object” probe surrounds the light-speed object, like a fleas swarming around a sprinting animal. The object modifies the probe, imprinting information about its structure. Meanwhile, our “reference” probe co-propagates ahead of the object, free of its influence. After the interaction, object and reference combine to record a hologram. For technical reasons, our recording device is a spectrometer (a frequency-measuring device), hence the name “frequency-domain” holography. We read the hologram electronically to obtain a “snapshot” of the object’s average structure as it transits the medium. Our published work shows numerous snapshots of electron density waves (“laser wakes”) in ionized gas (“plasma”), analogous to a water wake behind a boat. Such waves are the basis of tabletop particle accelerators, in which charged particles surf on the light-speed wave, gaining energy. Comparing our snapshots to computer simulations deepens understanding of laser wakes. FDH takes snapshots of objects that are quasi-static --- i.e. like Muybridge’s horse standing still on a treadmill. If the object changes shape, FDH images blur, as when a subject moves while a camera shutter is open. Many laser-generated objects of interest do evolve as they propagate. To overcome this limit of FDH, we developed .... (2) Frequency-Domain Tomography (FDT). In FDT, 5 to 10 probe pulses are fired simultaneously across the object’s path at different angles, like a crossfire of bullets. The object imprints a “streaked” record of its evolution on each probe, which we record as in FDH, then recover a multi-frame “movie” of the object’s evolving structure using algorithms of computerized tomography. When propagation distance exceeds a few millimeters, reconstructed FDT images distort. This is because the lenses that image probes to detector have limited depth of field, like cameras that cannot focus simultaneously on both nearby and distant objects. But some laser-generated objects of interest propagate over meters. For these applications we developed … (3) Multi-Object-Plane Phase-Contrast Imaging (MOP-PCI). In MOP-PCI, we image FDT-like probes to the detector from multiple “object planes” --- like recording an event simultaneously with several cameras, some focused on nearby, others on distant, objects. To increase sensitivity, we exploit a phase-contrast imaging technique developed by Dutch Nobel laureate Fritz Zernike in the 1930s. Using MOP-PCI we recorded single-shot movies of laser pulse tracks through more than 10 cm of air. We plan to record images of meter-long tracks of electron bunches propagating through plasma in an experiment at the Stanford Linear Accelerator Center (SLAC). This will help SLAC scientists understand, optimize and scale small plasma-based particle accelerators that have applications in medicine, industry, materials science and high-energy physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structure and magnetic properties of irradiated Fe/Fe oxide core-shell nanoclusters

A cluster deposition method was used to produce a film of loosely aggregated particles of Fe-Fe3O4 core-shell nanoclusters with an 8 nm iron core size and 2 nm oxide shell thickness. The film of particles on a silicon substrate was irradiated with 5.5 MeV Si2+ ions to a fluence of 1016 cm-2 near room temperature, and computer simulations based on the SRIM (Stopping and Range of Ions in Matter) code show that the implanted Si species stops near the film-substrate interface. The ion irradiation creates a structural change in the film with corresponding chemical and magnetic changes. X-ray diffraction shows that the core size and chemistry stay the same but the shell becomes FeO that grows to a thickness of 17 nm. Helium ion microscopy shows that the previously separate particles have densified into a nearly continuous film. Major loop magnetic hysteresis measurements show a decrease in saturation magnetization that we attribute to the presence of the antiferromagnetic (AFM) FeO shell. First-order reversal curve measurements on the irradiated film performed with a vibrating sample magnetometer show that the AFM shell prevents the particles from interacting magnetically, leading to low coercivity from the iron core and little bias field from the core interactions. These results, and others reported previously on different compositions (Fe3O4 or FeO+Fe3N nanoclusters), show that the ion irradiation behavior of nanocluster films such as these depends strongly on the initial nanostructure and chemistry.

Environmental Molecular Sciences Laboratory↗

Validation of Kestrel IDDES Simulations for SLS Transition Analysis

Complex computational simulations are needed to support the Space Launch System (SLS) program, and the fidelity of the computational results must be defended. More specifically, a number of databases are produced for the transition phase of flight, which occurs after the rocket clears the tower but before reaching transonic speeds. In an effort to reduce computational uncertainty, many of the computational parameters were altered to determine the sensitivity of the results to the value of the parameter and then updating the best practice procedures. The baseline routines were developed over many years through the maturation of the SLS program, and this paper delivers a detailed discussion of these baseline results. This section is followed by demonstration of perturbing some of the more significant parameters, including time step and turbulence model, and concluded with a summary of the herein determined best practices for such analysis.

CFD↗

Long period nodal motion of sun synchronous orbits

An approximative model is formulated for assessing these perturbations that significantly affect long term modal motion of sun synchronous orbits. Computer simulations with several independent computer programs consider zonal and tesseral gravitational harmonics, third body gravitational disturbances induced by the sun and the moon, and atmospheric drag. A pendulum model consisting of evenzonal harmonics through order 4 and solar gravity dominated nodal motion approximation. This pendulum motion results from solar gravity inducing an inclination oscillation which couples into the nodal precession induced by the earth's oblateness. The pendulum model correlated well with simulations observed flight data.

Duck, K. I.↗

A novel electron gun for inline MRI-linac configurations

Purpose: This work introduces a new electron gun geometry capable of robust functioning in the presence of a high strength external magnetic field for axisymmetric magnetic resonance imaging (MRI)-linac configurations. This allows an inline MRI-linac to operate without the need to isolate the linear accelerator (linac) using a magnetic shield. This MRI-linac integration approach not only leaves the magnet homogeneity unchanged but also provides the linac flexibility to move along the magnet axis of symmetry if the source to target distance needs to be adjusted. Methods: Simple electron gun geometry modifications of a Varian 600C electron gun are considered and solved in the presence of an external magnetic field in order to determine a set of design principles for the new geometry. Based on these results, a new gun geometry is proposed and optimized in the fringe field of a 0.5 T open bore MRI magnet (GE Signa SP). A computer model for the 6 MeV Varian 600C linac is used to determine the capture efficiency of the new electron gun-linac system in the presence of the fringe field of the same MRI scanner. The behavior of the new electron gun plus the linac system is also studied in the fringe fields of two other magnets, a 1.0 T prototype open bore magnet and a 1.5 T GE Conquest scanner. Results: Simple geometrical modifications of the original electron gun geometry do not provide feasible solutions. However, these tests show that a smaller transverse cathode diameter with a flat surface and a slightly larger anode diameter could alleviate the current loss due to beam interactions with the anode in the presence of magnetic fields. Based on these findings, an initial geometry resembling a parallel plate capacitor with a hole in the anode is proposed. The optimization procedure finds a cathode-anode distance of 5 mm, a focusing electrode angle of 5°, and an anode drift tube length of 17.1 mm. Also, the linac can be displaced with ±15 cm along the axis of the 0.5 T magnet without capture efficiency reduction below the experimental value in zero field. In this range of linac displacements, the electron beam generated by the new gun geometry is more effectively injected into the linac in the presence of an external magnetic field, resulting in approximately 20% increase of the target current compared to the original gun geometry behavior at zero field. The new gun geometry can generate and accelerate electron beams in external magnetic fields without current loss for fields higher than 0.11 T. The new electron-gun geometry is robust enough to function in the fringe fields of the other two magnets with a target current loss of no more than 16% with respect to the current obtained with no external magnetic fields. Conclusions: In this work, a specially designed electron gun was presented which can operate in the presence of axisymmetric strong magnetic fringe fields of MRI magnets. Computer simulations show that the electron gun can produce high quality beams which can be injected into a straight through linac such as Varian 600C and accelerated with more efficiency in the presence of the external magnetic fields. Also, the new configuration allows linac displacements along the magnet axis in a range equal to the diameter of the imaging spherical volume of the magnet under consideration. The new electron gun-linac system can function in the fringe field of a MRI magnet if the field strength at the cathode position is higher than 0.11 T. The capture efficiency of the linac depends on the magnetic field strength and the field gradient. The higher the gradient the better the capture efficiency. The capture efficiency does not degrade more than 16%.

43 PARTICLE ACCELERATORS↗

Development of an x-ray Talbot-Lau moiré deflectometer for fast density profile measurements of dense plasmas generated by beam-target interactions

We report that for the first time an x-ray Talbot-Lau moiré deflectometer is being developed that will use a flash tube source and fast detector for dynamic density gradient measurements. In Talbot-Lau moiré deflectometry, an x-ray grating makes an image of itself on a second grating (the Talbot effect) to produce a moiré pattern on a detector. The test object is placed between these gratings, with variations in index of refraction changing the pattern. A third grating in front of an incoherent x-ray source produces an array of coherent sources. With a 150 kV x-ray flash tube as the source, the gratings are placed in a glancing angle setup for performance at ~60 keV. The detector is a gated CCD with a fast scintillator for x-ray conversion. This diagnostic, designed for the Dual-Axis Radiographic Hydrodynamic Test facility (DARHT) at Los Alamos National Laboratory, measures the density profile of dense plasma plumes ejected from beam-target interactions. DARHT has two high-current, pulsed, inductive linear electron accelerators with bremsstrahlung targets at the end of each beam line to create 2-D radiographic images of hydrodynamic tests. One multi-pulse accelerator has up to four beam pulses striking the same target within 2 μs. Computer simulations that model target evolution and ejected material between pulses are used to design these targets for optimal radiographic performance; the x-ray deflectometer will directly measure density gradients in the ejected plumes and provide the first experimental constraints to these models. During the first year, currently underway, the diagnostic systems are being designed. In year two, the flash tube and fast detector will be deployed at DARHT for radiographic imaging while the deflectometer is built and tested on the bench with a continuous source. In year three, the fast deflectometer will be installed on DARHT and density measurements will be performed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Review of Discrete Element Method (DEM) Particle Shapes and Size Distributions for Lunar Soil

As part of ongoing efforts to develop models of lunar soil mechanics, this report reviews two topics that are important to discrete element method (DEM) modeling the behavior of soils (such as lunar soils): (1) methods of modeling particle shapes and (2) analytical representations of particle size distribution. The choice of particle shape complexity is driven primarily by opposing tradeoffs with total number of particles, computer memory, and total simulation computer processing time. The choice is also dependent on available DEM software capabilities. For example, PFC2D/PFC3D and EDEM support clustering of spheres; MIMES incorporates superquadric particle shapes; and BLOKS3D provides polyhedra shapes. Most commercial and custom DEM software supports some type of complex particle shape beyond the standard sphere. Convex polyhedra, clusters of spheres and single parametric particle shapes such as the ellipsoid, polyellipsoid, and superquadric, are all motivated by the desire to introduce asymmetry into the particle shape, as well as edges and corners, in order to better simulate actual granular particle shapes and behavior. An empirical particle size distribution (PSD) formula is shown to fit desert sand data from Bagnold. Particle size data of JSC-1a obtained from a fine particle analyzer at the NASA Kennedy Space Center is also fitted to a similar empirical PSD function.

Lane, John E.↗

Preface to the special volume on the third Sandia Fracture Challenge

The mounting reliance on computational simulations to predict all aspects of the lifecycle of a mechanical system, from fabrication to failure, has prompted the mechanics community to selfassess its abilities to perform those predictions. Benchmark problems in mechanics that compare simulations that use different computational approaches with experiments have sprung up lately, including the NIST AM-Bench looking at additively manufactured (AM) materials (https://www.nist.gov/ambench),the Contact-Mechanics Challenge (Miiser, 2017) considering adhesion between two nominally flat surfaces, Numisheet providing semiannual benchmarking activities in sheet metal forming (http://numisheet2018.org),and the Sandia Fracture Challenge (SFC) (Boyce, 2014 and Boyce, 2016) investigating ductile failure. The previous SFCs have shown that progress has been made in computations of ductile failure, but improvements still can be made, hence the third Sandia Fracture Challenge (SFC3), the subject of this Special Volume. The most recent installment of SFC is building on previous successes and tackling the difficult problem of fracture in an AM 316L stainless steel structure.

42 ENGINEERING↗

Program optimizations: The interplay between power, performance, and energy

Practical considerations for future supercomputer designs will impose limits on both instantaneous power consumption and total energy consumption. Working within these constraints while providing the maximum possible performance, application developers will need to optimize their code for speed alongside power and energy concerns. This paper analyzes the effectiveness of several code optimizations including loop fusion, data structure transformations, and global allocations. A per component measurement and analysis of different architectures is performed, enabling the examination of code optimizations on different compute subsystems. Using an explicit hydrodynamics proxy application from the U.S. Department of Energy, LULESH, we show how code optimizations impact different computational phases of the simulation. This provides insight for simulation developers into the best optimizations to use during particular simulation compute phases when optimizing code for future supercomputing platforms. Here, we examine and contrast both x86 and Blue Gene architectures with respect to these optimizations.

97 MATHEMATICS AND COMPUTING↗

Computational Molecular Simulation of the Oxidative Adsorption of Ferrous Iron at the Hematite (001)-Water Interface

The interaction of Fe(II) with ferric oxide/oxyhydroxide phases is central to the biogeochemical redox chemistry of iron. Molecular simulation techniques were employed to determine the mechanisms and quantify the rates of Fe(II) oxidative adsorption at the hematite (001)-water interface. Molecular dynamics potential of mean force calculations of Fe(II) adsorbing on the hematite surface revealed the presence of three free energy minima corresponding to Fe(II) adsorbed in an outersphere complex, a monodentate innersphere complex, and a tridentate innersphere complex. The free energy barrier for adsorption from the outersphere position to the monodentate innersphere site was calculated to be similar to the activation enthalpy for water exchange around aqueous Fe(II). Adsorption at both innersphere sites was predicted to be unfavorable unless accompanied by release of protons. Molecular dynamics umbrella sampling simulations and ab initio cluster calculations were performed to determine the rates of electron transfer from Fe(II) adsorbed as an innersphere and outersphere complex. The electron transfer rates were calculated to range from 10^-4 to 10^2 s-1, depending on the adsorption site and the potential parameter set, and were generally slower than those obtained in the bulk hematite lattice. The most reliable estimate of the rate of electron transfer from Fe(II) adsorbed as an outersphere complex to lattice Fe(III) was commensurate with the rate of adsorption as an innersphere complex suggesting that adsorption does not necessarily need to precede oxidation.

Environmental Molecular Sciences Laboratory↗

Behavior of nanoporous thermal barrier coatings under cyclic thermal loading. Computer-aided simulation

The work presents the development of new algorithms for calculating the fraction of thermal energy dissipated during the irradiation on the inner surface of pores. On the basis of these algorithms, the simulation of heat transfer in three-layered systems was carried out taking into account the dissipation of thermal energy in specimens having different porosity. We have performed quantitative estimates of the portion of dissipating thermal energy and its influence on the distribution of thermal stresses in thermal barrier coating systems. It was demonstrated that the presence of pores with a large internal surface area in the intermediate layer material enables two-fold decrease of the internal thermal stresses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗