Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “D CODES”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Boundary condition effects on runaway electron mitigation coil modeling for the SPARC and DIII-D tokamaks

Abstract Extended-MHD modeling of planned Runaway Electron Mitigation Coils (REMC) for SPARC and DIII-D is performed with the NIMROD code. A coil has been designed for each machine, with the two differing in shape and location, but both having n = 1 symmetry (with n the toroidal mode number). Compared to previous modeling efforts, three improvements are made to the simulations boundary conditions. First a resistive wall model is used in place of an ideal wall. Second, the ThinCurr code is used to compute the time-dependent 3D fields used as magnetic boundary conditions for the simulations. Third, the simulation boundary is moved from the first-wall location to the Vacuum Vessel (VV), which extends the boundary past the location of the internal REMC. To remove the 3D coil from the simulation domain, an equivalent set of 3D fields is calculated at the VV boundary that produce approximately the same field distribution at the last closed flux surface assuming vacuum between the two. Each of these three boundary condition improvements leads to an improvement in the predicted performance of the REMC for both machines. The resistive wall alone primarily effects the resonance of the coil with the plasma after the TQ, affecting the q-profile evolution in the SPARC modeling, and allowing the applied spectrum to be modified in response to the plasma in the DIII-D modeling. The movement of the simulation boundary has the most significant effect on the RE confinement overall, including in the early stages, particularly for a DIII-D inner wall limited equilibrium, where the RE loss fraction increases from 90% to > 99%, with SPARC RE losses also occurring much earlier when the boundary is placed at the VV.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

Low-Rate Turbo Codes for Deep-Space Communications

It is shown how turbo codes and decoders can be used to improve the coding gain for deep-space communications, while decreasing the decoding complexity with respect to the large constraint length convolutional codes currently in use. Similar code constructions were used to build multiple-encoder turbo codes. This generalizes the turbo decoding concept to a truly distributed decoding system.

turbo codes deep-space communications coding decod↗
NASA NTRS2011

Towards Simulating a Realistic Planetary Seismic Wavefield: The Contribution of the Megaregolith and Low-Velocity Waveguides

Lunar seismograms are distinctly different from their terrestrial counterparts. The Apollo lunar seismometers recorded moonquakes without distinct P- or S-wave arrivals; instead waves arrive as a diffuse coda that decays over several hours making the identification of body waves difficult. The unusual character of the lunar seismic wavefield is generally tied to properties of the megaregolith: it consists of highly fractured and broken crustal rock, the result of extensive bombardment of the Moon. The megaregolith extends several kilometers into the lunar crust, possibly into the mantle in some regions, and is covered by a thin coating of fine-scale dust. These materials possess very low seismic velocities that strongly scatter the seismic wavefield at high frequencies. Directly modeling the effects of the megaregolith to simulate an accurate lunar seismic wavefield is a challenging computational problem, owing to the inherent 3-D nature of the problem and the high frequencies (greater than 1 Hz) required. Here we focus on modeling the long duration code, studying the effects of the low velocities found in the megaregolith. We produce synthetic seismograms using 1-D slowness integration methodologies, GEMINI and reflectivity, and a 3-D Cartesian finite difference code, Wave Propagation Program, to study the effect of thin layers of low velocity on the surface of a planet. These codes allow us generate seismograms with dominant frequencies of approximately 1 Hz. For background lunar seismic structure we explore several models, including the recent model of Weber et al., Science, 2011. We also investigate variations in megaregolithic thickness, velocity, attenuation, and seismogram frequency content. Our results are compared to the Apollo seismic dataset, using both a cross correlation technique and integrated envelope approach to investigate coda decay. We find our new high frequency results strongly support the hypothesis that the long duration of the lunar seismic codes is generated by the presence of the low velocity megaregolith, and that the diffuse arrivals are a combination of scattered energy and multiple reverberations within this layer. The 3-D modeling indicates the extreme surface topography of the Moon adds only a small contribution to scattering effects, though local geology may play a larger role. We also study the effects of the megaregolith on core reflected and converted phases and other body waves. Our analysis indicates detection of core interacting arrivals with a polarization filter technique is robust and lends the possibility of detecting other body waves from the Moon.

Schmerr, Nicholas C.↗
DOE OSTI2024

Upgrade of the Lyman-alpha diagnostic system on DIII-D for main chamber edge neutral studies

The LLAMA (Lyman Alpha Measurement Apparatus) pinhole camera diagnostic had previously been deployed on DIII-D to measure radial profiles of the Lyman-α (Ly-α) deuterium neutral line brightness across the plasma boundary in the lower chamber to infer neutral deuterium density and ionization rate profiles. This system has recently been upgraded with a new diagnostic head, named ALPACA, that also encloses two pinhole cameras and duplicates the LLAMA views in the upper chamber. Similar to LLAMA, ALPACA provides two times 20 lines of sight, viewing the plasma edge on the inboard and outboard sides with a radial resolution of ~2.5 cm (FWHM) and an effective time resolution of ~1 ms that allows for the investigation of inter-ELM dynamics. The extended Ly-α system provides better coverage to study neutrals in experiments with various plasma shapes utilizing both the upper and lower divertors. Furthermore, post-campaign calibration of the LLAMA diagnostic has successfully been demonstrated for the first time. This was facilitated by various upgrades to the calibration set-up and detailed measurements of the emissivity distribution of the Ly-α calibration source using a pinhole collimator. It was found that the sensitivity of the inboard LLAMA pinhole camera was reduced by a factor of 2.0 ± 0.2 over the course of six months of plasma operation in 2021. In conclusion, the upgraded Ly-α system, equipped with improved absolute calibration, will provide key input for neutral fueling and pedestal particle transport studies and for 2D edge transport code validation on the DIII-D tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Gyrokinetic analysis of inter-edge localized mode transport mechanisms in a DIII-D pedestal

In this study, gyrokinetic simulations are used to study pedestal fluctuations for DIII-D discharge 174082 using the GENE code. Nonlinear local simulations indicate that electron heat flux has contributions from electron temperature gradient-driven transport but at levels insufficient to satisfy power balance. We show that microtearing modes (MTM) and neoclassical transport are likely to account for the remaining observed energy losses in the electron and ion channels, respectively. The MTM instabilities found in the simulations are consistent with the high-frequency fluctuations identified in the magnetic fluctuation data from Mirnov coils. The fluctuation data in this discharge also exhibit a low-frequency band of fluctuations. By modifying the equilibrium profiles and plasma β, simulations produce MHD modes, which may be responsible for these observed low-frequency fluctuations. We compare several metrics involving ratios of fluctuation amplitudes and transport quantities for both MTMs and MHD modes. This analysis suggests that the available data are consistent with the simultaneous activity of both MHD modes and MTMs provided that the former is limited largely to the particle transport channel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Dispersive shell pellet injection modeling and validation for DIII-D disruption mitigation

Dispersive shell pellet injection for disruption mitigation on the DIII-D tokamak is modeled with the 3D MHD code NIMROD. The simulations attempt to reproduce as closely as possible a pellet velocity scan that was performed on the DIII-D tokamak [Hollmann, et al., Phys. Rev. Lett. 122, 65001 (2019)]. This scan involved injecting a hollow diamond shell filled with boron dust that is released in the core of the plasma. Trends that are reproduced in the simulation include improved mitigation, reduced plasma current spike amplitude, and increased likelihood of hot-tail runaway electron production with increased pellet speed. The simulations provide additional physics insight into these findings. Limitations of the present model and requirements for a predictive model are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

The development and application of the self-adaptive grid code, SAGE

The multidimensional self-adaptive grid code, SAGE, has proven to be a flexible and useful tool in the solution of complex flow problems. Both 2- and 3-D examples given in this report show the code to be reliable and to substantially improve flowfield solutions. Since the adaptive procedure is a marching scheme the code is extremely fast and uses insignificant CPU time compared to the corresponding flow solver. The SAGE program is also machine and flow solver independent. Significant effort was made to simplify user interaction, though some parameters still need to be chosen with care. It is also difficult to tell when the adaption process has provided its best possible solution. This is particularly true if no experimental data are available or if there is a lack of theoretical understanding of the flow. Another difficulty occurs if local features are important but missing in the original grid; the adaption to this solution will not result in any improvement, and only grid refinement can result in an improved solution. These are complex issues that need to be explored within the context of each specific problem.

Davies, Carol B.↗
NASA NTRS2009

Simulation of Relativistic Shocks and Associated Radiation from Turbulent Magnetic Fields

Plasma instabilities (e.g., Buneman, Weibel and other two-stream instabilities) excited in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a new 3-D relativistic particle-in-cell code, we have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. The simulation has been performed using a long simulation system in order to study the nonlinear stages of the Weibel instability, the particle acceleration mechanism, and the shock structure. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic (HD) like shock structure. In the leading shock, electron density increases by a factor of <_ 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. We discuss the possible implication of our simulation results within the AGN and GRB context. We have calculated the time evolution of the spectrum from two electrons propagating in a uniform parallel magnetic field to verify the technique. The same technique will be used to calculate radiation from accelerated electrons (positrons) in turbulent magnetic fields generated by Weibel instability.

Nishikawa, K.↗
NASA NTRS1994

Error control techniques for satellite and space communications

The unequal error protection capabilities of convolutional and trellis codes are studied. In certain environments, a discrepancy in the amount of error protection placed on different information bits is desirable. Examples of environments which have data of varying importance are a number of speech coding algorithms, packet switched networks, multi-user systems, embedded coding systems, and high definition television. Encoders which provide more than one level of error protection to information bits are called unequal error protection (UEP) codes. In this work, the effective free distance vector, d, is defined as an alternative to the free distance as a primary performance parameter for UEP convolutional and trellis encoders. For a given (n, k), convolutional encoder, G, the effective free distance vector is defined as the k-dimensional vector d = (d(sub 0), d(sub 1), ..., d(sub k-1)), where d(sub j), the j(exp th) effective free distance, is the lowest Hamming weight among all code sequences that are generated by input sequences with at least one '1' in the j(exp th) position. It is shown that, although the free distance for a code is unique to the code and independent of the encoder realization, the effective distance vector is dependent on the encoder realization.

Costello, Daniel J., Jr.↗
NASA NTRS1992

Solution of nonlinear flow equations for complex aerodynamic shapes

Solution-adaptive CFD codes based on unstructured methods for 3-D complex geometries in subsonic to supersonic regimes were investigated, and the computed solution data were analyzed in conjunction with experimental data obtained from wind tunnel measurements in order to assess and validate the predictability of the code. Specifically, the FELISA code was assessed and improved in cooperation with NASA Langley and Imperial College, Swansea, U.K.

Djomehri, M. Jahed↗
DOE OSTI2021

Modelling dust transport in DIII-D with DTOKS-Upgrade

Comprehensive upgrades to the dust transport code Dust in TOKamaks (DTOKS) that extend the plasma-dust interaction model are presented and compared with recent measurements of dust transport in DIII-D. Simulations incorporating variation in physical properties of graphite dust with temperature and size in a stationary plasma background suggest a substantial decrease in lifetimes due principally to thermal expansion. The trajectories of 53 dust grains identified from analysis of visible camera data taken across two similar shots were used to measure the dust particle velocity distributions. Dust tracks terminated mostly at the outer divertor strike point having a mean observation time of 2.1 ± 0.4 ms. Stochastic modelling of 200 graphite dust particles in the DIII-D tokamak performed with DTOKS-U using plasma simulations generated by OEDGE found similar behaviour, with particles ablating rapidly after acquiring a positive charge in the region close to the outer strike point, creating an acute source of neutral carbon atoms. The simulated mean lifetime, 11 ± 2 ms, showed approximate agreement with experimental observation when corrected by accounting for dust visibility and ignoring the longest trajectories 1.5 ± 0.2 ms. Synthetic diagnostic data generated from coupling the results of DTOKS-U with the visualisation software Calcam offers a powerful new tool for validation of simulations and predictive calculations of dust dynamics.

Physics↗
DOE OSTI2022

A smeared crack modeling framework accommodating multi-directional fracture at finite strains

A generic smeared crack modeling framework predicated on the deformation gradient decomposition (DGD) approach is proposed for use in dynamic fracture problems at finite strains, accommodating failure along multiple mutually orthogonal fracture planes embedded within an independently defined bulk material model. Within this constitutive framework, the traction equilibrium conditions imposed at each failure surface are used to determine the associated crack displacements stored as internal state variables. In general, the enforcement of interfacial equilibrium entails the implicit solution of a non-linear system of equations within the constitutive update procedure. However, if inertial effects arising due to the relative motion of the fractured material are incorporated within the model, the traction equilibrium conditions are shown to give rise to corresponding dynamic equations of motion governing the time-evolution of the crack opening displacements. For dynamic problems, an explicit time-integration procedure is devised to efficiently update the material state, subject to a set of internal frictionless contact constraints to prevent material inversion. Finally, the efficacy of the proposed modeling framework is investigated through several benchmark dynamic fracture problems run within the explicit finite element code DYNA3D.

42 ENGINEERING↗
DOE OSTI2020

Prediction of tissue optical properties using the Monte Carlo modeling of photon transport in turbid media and integrating spheres

Monte Carlo methods are an established technique for simulating light transport in biological tissue. Integrating spheres make experimental measurements of the reflectance and transmittance of a sample straightforward and inexpensive. This work presents an extension to existing Monte Carlo photon transport methods to simulate integrating sphere experiments. Crosstalk between spheres in dual-sphere experiments is accounted for in the method. Analytical models, previous works on Monte Carlo photon transport, and experimental measurements of a synthetic tissue phantom validate this method. We present two approaches for using this method to back-calculate the optical properties of samples. Experimental and simulation uncertainties are propagated through both methods. Both back-calculation methods find the optical properties of a sample accurately and precisely. Our model is implemented in standard Python 3 and CUDA C++ [J. Nickolls, I. Buck, M. Garland, and K. Skadron, ACM Queue 6 , 40 ( 2008 ) ] and is publicly available in Code 1.

Cook, Patrick D. (ORCID:0000000279345428)↗
NASA NTRS1967

The organization of skilled response.

Skilled response organization, discussing stimulus coherence, tracking task, spatial and temporal coherence, secondary task, sequence length and task coding

STIMULUS↗
NASA NTRS1980

A supersonic three-dimensional code for flow over blunt bodies: Program documentation and test cases

The use of a computer code for the calculation of steady, supersonic, three dimensional, inviscid flow over blunt bodies is illustrated. Input and output are given and explained for two cases: a pointed code of 20 deg half angle at 15 deg angle of attack in a free stream with M sub infinite = 7, and a cone-ogive-cylinder at 10 deg angle of attack with M sub infinite = 2.86. A source listing of the computer code is provided.

Chaussee, D. S.↗