Wide-scale Adoption through Visually-Enhanced HPC codes (WAVE-HPC) for semiconductor R&D
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
In this paper, we demonstrate a class of low-complexity modulation codes satisfying the (d,k) constraint that offer throughput gains over M-PPM on the order of 10-15%, which translate into SNR gains of .4 - .6 dB.
ASME Code Development - Nonmetallics slides to include behavior models, licensing & code, graphite R&D Program, As-Fab'd Properties, Mechanisms and Analysis, and Irradiation. Discussion on where these fit into Irradiation behavior, Ceramic Composites, Design, and Component Failure.
A thermography inversion algorithm has been developed in the open-source Python-based computer code, HYPERION, to calculate the heat flux incident on plasma-facing components (PFCs) in axisymmetric tokamaks. The chosen mesh size at the surface significantly affects the calculated transient heat flux results. The calculated transient heat flux will exceed the real value when the mesh size tends to zero but will underestimate the real value when the mesh size is large. A criterion for determining the appropriate mesh size for the transient heat flux calculation will be discussed. The numerical scheme for HYPERION uses a 2D fully implicit finite-difference approach, allowing temperature-dependent thermal properties of PFC materials. The inversion algorithm is benchmarked against established heat flux calculation codes, TACO and THEODOR, based on thermography data from NSTX and DIII-D respectively. The primary benefits of HYPERION compared to TACO and THEODOR are that it is open-source and it allows for the optimization of mesh thickness along the substrate. The algorithm also accounts for the thermal properties of thin surface layers that characteristically form on PFCs due to plasma-material interactions. The agreement between HYPERION and THEODOR is excellent, as the percent difference between the codes is ~5% on average in the case of the DIII-D data for moderate to high heat flux. Verification tests with TACO show slightly higher average percent differences of 8% and 12%. In using HYPERION to study filaments in heat flux, the initial results indicate that small ELMs filaments significantly broaden the divertor heat flux, and decrease divertor peak flux. Compared to the inter-ELM, the small ELM filaments decrease the divertor peak surface temperature. With intermittent divertor filaments, the divertor heat flux width is comparable with that found in L-mode.
We report the first ultraviolet spectropolarimetry of the peculiar B-type emission-line star with infrared excess HD 45677. The observations were obtained during the 1990 December Astro-l space shuttle mission with the Wisconsin Ultraviolet Photo-Polarimeter Experiment, and cover the spectral range 1400-3220 A. We also present the first optical spectropolarimetry, extending from the atmospheric cutoff to about 7600 A. The observed UV/optical linear polarization displays a strong increase toward shorter wavelengths indicative of scattering by circumstellar dust. The position angle of the intrinsic polarization flips by 90 deg in the near-UV as expected from a bipolar reflection nebula.
We present here a massively parallel implementation of the recently developed CPS(D-3) excitation energy model that is based on cluster perturbation theory. The new algorithm extends the one developed in Baudin et al. [J. Chem. Phys., 150, 134110 (2019)] to leverage multiple nodes and utilize graphical processing units for the acceleration of heavy tensor contractions. Furthermore, we show that the extended algorithm scales efficiently with increasing amounts of computational resources and that the developed code enables CPS(D-3) excitation energy calculations on large molecular systems with a low time-to-solution. More specifically, calculations on systems with over 100 atoms and 1000 basis functions are possible in a few hours of wall clock time. This establishes CPS(D-3) excitation energies as a computationally efficient alternative to those obtained from the coupled-cluster singles and doubles model.
A three dimensional body-fitted coordinate system developed for use in the calculation of inviscid flows over ablated, asymmetric reentry vehicle nosetips is described. Because of the potential geometric asymmetries, no standard coordinate system (e.g., spherical, axisymmetric reference surface-normal) is capable of being closely aligned with the nosetip surface. To generate a 3-D, body-fitted coordinate system an analytic mapping procedure is applied that is conformal within each meridional plane of the nosetip; these transformations are then coupled circumferentially to yield a three dimensional coordinate system. The mappings used are defined in terms of hinge points, which are points selected to approximate the body contours in each meridional plane. The selection of appropriate hinge points was automated to facilitate the use of the resulting nosetip flow field code.
This paper will report on (1) continuation through the 3rd year of a NASA grant for multi-dimensional Stirling CFD code development and validation, (2) continuation through the 3rd and final year of a Department of Energy, Golden Field Office (DOE) regenerator research effort. Results of the NASA multi-D code development effort and the DOE regenerator research efforts will be summarized. Early results and planning for the new regenerator microfabrication contract will also be discussed.
Tests were conducted in a Mach 6 flow to examine the reattachment process of an axisymmetric free shear layer associated with the near wake of a 70 deg. half angle, spherically blunted cone with a cylindrical after body. Model angle of incidence was fixed at 0 deg. and free-stream Reynolds numbers based on body diameter ranged from 0.5 x 10(exp 6) to 4 x 10(exp 6). The sensitivity of wake shear layer transition on reattachment heating was investigated. The present perfect gas study was designed to compliment results obtained previously in facilities capable of producing real gas effects. The instrumented blunted cone model was designed primarily for testing in high enthalpy hypervelocity shock tunnels in both this country and abroad but was amenable for testing in conventional hypersonic blowdown wind tunnels as well. Surface heating rates were inferred from temperature - time histories from coaxial surface thermocouples on the model forebody and thin film resistance gages along the model base and cylindrical after body. General flow feature (bow shock, wake shear layer, and recompression shock) locations were visually identified by schlieren photography. Mean shear layer position and growth were determined from intrusive pitot pressure surveys. In addition, wake surveys with a constant temperature hot-wire anemometer were utilized to qualitatively characterize the state of the shear layer prior to reattachment. Experimental results were compared to laminar perfect gas predictions provided by a 3-D Navier Stokes code (NSHYP). Shear layer impingement on the instrumented cylindrical after body resulted in a localized heating maximum that was 21 to 29 percent of the forebody stagnation point heating. Peak heating resulting from the reattaching shear layer was found to be a factor of 2 higher than laminar predictions, which suggested a transitional shear layer. Schlieren flow visualization and fluctuating voltage time histories and spectra from the hot wire surveys across the shear layer substantiate this observation. The sensitivity of surface heating to forebody roughness was characterized for a reattaching shear layer. For example, at R(sub infinity), d = 4 x 10(exp 6), when the shear layer was transitional, the magnitude of peak heating from shear layer impingement was reduced by approximately 24 percent when transition grit was applied to the forebody. The spatial location of the local peak, however, remained unchanged.
We examine a "two-component" model of the solar wind to see if any of the observed anisotropies of the fields can be explained in light of the need for various quantities, such as the magnetic minimum variance direction, to turn along with the Parker spiral. Previous results used a 3-D MHD spectral code to show that neither Q2D nor slab-wave components will turn their wave vectors in a turning Parker-like field, and that nonlinear interactions between the components are required to reproduce observations. In these new simulations we use higher resolution in both decaying and driven cases, and with and without a turning background field, to see what, if any, conditions lead to variance anisotropies similar to observations. We focus especially on the middle spectral range, and not the energy-containing scales, of the simulation for comparison with the solar wind. Preliminary results have shown that it is very difficult to produce the required variances with a turbulent cascade.
Poloidal field (PF) and central solenoid (CS) coils play a crucial role in sustaining the equilibrium and preserving the shape of highly confined tokamak plasmas. Ensuring that PF coil current and mechanical stress stay within superconducting and structural limitations is an important check in the design assessment. Minimizing the PF coil currents and mechanical stresses influences reliability, cost, and performance. A free-boundary MHD equilibrium code—FreeGS is employed within the fusion reactor design and assessment (FREDA) whole facility modeling (WFM) framework to construct the plasma equilibrium based on the configuration and currents in the PF coils. Here, we present the capability of the FreeGS code to minimize the currents, forces, and electromagnetic stresses on the PF coils by optimizing their number, sizes, structures, and locations while maintaining an MHD stable plasma configuration with a large confinement factor. The workflow is initialized with a configuration of plasma parameters and coils’ locations from the 0-D tokamak build systems code in the FREDA framework. Then, FreeGS is called to calculate the initial equilibrium at the minimum total current in PF coils. Thereafter, FreeGS’s internal optimizer minimizes the currents and hoop and central forces on the PF coils while maintaining the reference equilibrium. Finally, the input configuration is updated with the optimized parameters for equilibria over the ramp-up phase of a burning-plasma operation. FREDA’s whole facility optimization capability, which includes all magnetic field coil systems, blanket, vacuum vessel (VV), first wall, divertor, etc., is under development and out of the scope for this study.
A production 3-D elliptic flow code has been used to calculate non-reacting and reacting flow fields in an experimental mixing section relevant to a rich burn/quick mix/lean burn (RQL) combustion system. A number of test cases have been run to assess the effects of the variation in the number of orifices, mass flow ratio, and rich-zone equivalence ratio on the flow field and mixing rates. The calculated normalized temperature profiles for the non-reacting flow field agree qualitatively well with the normalized conserved variable isopleths for the reacting flow field indicating that non-reacting mixing experiments are appropriate for screening and ranking potential rapid mixing concepts. For a given set of jet momentum-flux ratio, mass flow ratio, and density ratio (J, MR, and DR), the reacting flow calculations show a reduced level of mixing compared to the non-reacting cases. In addition, the rich-zone equivalence ratio has noticeable effect on the mixing flow characteristics for reacting flows.
We solve the Lindblad equation describing the Brownian motion of a Coulombic heavy quark-antiquark pair in a strongly coupled quark-gluon plasma using the highly efficient Monte Carlo wave-function method. The Lindblad equation has been derived in the framework of pNRQCD and fully accounts for the quantum and non-Abelian nature of the system. The hydrodynamics of the plasma is realistically implemented through a 3+1D dissipative hydrodynamics code. We compute the bottomonium nuclear modification factor and compare with the most recent LHC data. The computation does not rely on any free parameter, as it depends on two transport coefficients that have been evaluated independently in lattice QCD. Our final results, which include late-time feed down of excited states, agree well with the available data from LHC 5.02 TeV PbPb collisions.
Laser plasma instabilities (LPI) reduce driver-target coupling, alter implosion symmetry, and therefore can fundamentally limit fusion performance in inertial confinement fusion (ICF). Developing a predictive modeling capability for LPI effects can critically advance the success of the field. We perform vector particle-in-cell simulations of multi-speckled laser beams undergoing stimulated Raman scattering (SRS) at various densities and intensities relevant to mainly indirectly driven and a subset of parameter space for directly driven ICF systems, focusing on the regimes with intensities above the onset of electron trapping. Based on the wavenumber of the SRS daughter electron plasma wave, we identify several regions with underpinning SRS saturation physics: the electron-trapping dominated region with intermediate kλD values, the strong Landau damping region at larger kλD values, and the region where the Langmuir decay instability arises at lower kλD values. We develop a nonlinear SRS reflectivity model that features the base trapping-dominated scaling of (kλD)−4 and its modifications. Electron trapping deforms the initialized electron distribution functions, and we have developed a new δf-Gaussian-mixture algorithm for an accurate characterization of the trapped hot electron population. With this SRS hot electron description, we construct a nonlinear energy deposition model and a hot electron source model—based on a modified Manley–Rowe relation—suitable for including SRS effects as a sub-grid module in a high-fidelity ICF design code.
Recent measurements at the Sandia National Laboratory of the x-ray transmission of iron plasma have inferred opacities much higher than predicted by theory, which casts doubt on modeling of iron x-ray radiative opacity at conditions close to the solar convective zone-radiative zone boundary. An increased radiative opacity of the solar mixture, in particular iron, is a possible explanation for the disagreement in the position of the solar convection zone-radiative zone boundary as measured by helioseismology and predicted by modeling using the most recent photosphere analysis of the elemental composition. In this report we present data from radiation burnthrough experiments, which do not support a large increase in the opacity of iron at conditions close to the base of the solar convection zone and provide a constraint on the possible values of both the mean opacity and the opacity in the x-ray range of the Sandia experiments. The data agree with opacity values from current state-of-the-art opacity modeling using the CASSANDRA opacity code.
Space charge (SC) forces of a circulating beam in a ring have both linear (defocusing) and nonlinear components, due to a nonuniform beam distribution. The linear component of SC forces produces a betatron tune shift, which is the largest for a zero-amplitude particle, while the nonlinear component produces an amplitude-dependent betatron tune spread. These SC effects are responsible for several undesirable phenomena in accelerators: emittance growth, particle losses, beam halo, etc. In this paper, we investigate the possibility to mitigate the distributed SC forces by a thin McMillan lens, providing an axially-symmetric kick, which is qualitatively opposite to the accumulated effect of beam's own SC. Experimentally, the proposed concept can be tested in Fermilab's IOTA ring. A thin McMillan lens can be implemented by a short (70) insertion of an electron beam with a specifically chosen density distribution in transverse directions. In this article, to test if McMillan lenses can reduce the tune spread induced by SC, we make several simulations with a 6-D particle tracking code, Synergia. We choose such beam and lattice parameters that the SC tune spread is roughly 0.5 and the emittance growth due to the half-integer resonance is clearly observed without the SC compensation. Then, we focus on reducing the emittance growth by adjusting the bare betatron tunes using the ring quadrupoles, and reducing the tune spread by the McMillan lenses. The results of reducing a large tune spread (≈ 0.5), reported here, are not perfect, but substantial. There is still room for further investigation. The simulations performed so far in this work indicate that McMillan lenses can cope with an SC tune spread of ≈ 0.1 per lens.
TOUGHREACT V4.12 is a 3-D reactive transport code that handles to non-isothermal, multiphase/multicomponent fluid flow, heat transport, aqueous and gaseous species advection-diffusion, and equilibrium/kinetic water-gas-rock-biological reactions. TOUGHREACT V4.12 is based on TOUGHREACT V3.3, with several new features including: an option to compute activity coefficients using the Pitzer ion interaction model; implementation of ECO2N V2.0 for temperatures up to 300 C; options to simulate forward and reverse osmosis, including 2D Poiseuille flow; addition of EOS8 and revisions for EOS5 and EWASG with H2(g); addition of gas species decay, various new output formats and time step control options; many bug fixes.
Algorithm for constructing random error correcting convolutional codes