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 91 records · Page 5

The computation of steady 3-D separated flows over aerodynamic bodies at incidence and yaw

This paper describes the implementation of a general purpose 3-D NS code and its application to simulated 3-D separated vortical flows over aerodynamic bodies. The thin-layer Reynolds-averaged NS equations are solved by an implicit approximate factorization scheme. The pencil data structure enables the code to run on very fine grids using only limited incore memories. Solutions of a low subsonic flow over an inclined ellipsoid are compared with experimental data to validate the code. Transonic flows over a yawed elliptical wing at incidence are computed and separations occurred at different yaw angles are discussed.

Pulliam, T. H.↗

An MHD Code for the Study of Magnetic Structures in the Solar Wind

We have developed a 2.5D MHD code designed to study how the solar wind influences the evolution of transient events in the solar corona and inner heliosphere. The code includes thermal conduction, coronal heating and radiative cooling. Thermal conduction is assumed to be magnetic field-aligned in the inner corona and transitions to a collisionless formulation in the outer corona. We have developed a stable method to handle field-aligned conduction around magnetic null points. The inner boundary is placed in the upper transition region, and the mass flux across the boundary is determined from 1D field-aligned characteristics and a 'radiative energy balance' condition. The 2.5D nature of this code makes it ideal for parameter studies not yet possible with 3D codes. We have made this code publicly available as a tool for the community. To this end we have developed a graphical interface to aid in the selection of appropriate options and a graphical interface that can process and visualize the data produced by the simulation. As an example, we show a simulation of a dipole field stretched into a helmet streamer by the solar wind. Plasmoids periodically erupt from the streamer, and we perform a parameter study of how the frequency and location of these eruptions changed in response to different levels of coronal heating. As a further example, we show the solar wind stretching a compact multi-polar flux system. This flux system will be used to study breakout coronal mass ejections in the presence of the solar wind.

corona↗

Gauging nexus between topological and fracton phases

Coupled layer constructions are a valuable tool for capturing the universal properties of certain interacting quantum phases of matter in terms of the simpler data that characterizes the underlying layers. In the study of fracton phases, the X-Cube model in 3+1D can be realized via such a construction by starting with a stack of 2+1D Toric Codes and turning on a coupling which condenses a composite "particle-string" object. In a recent work [Phys. Rev. B 112, 125124 (2025)], we have demonstrated that in fact, the particle-string can be viewed as a symmetry defect of a topological 1-form symmetry. In this paper, we study the result of gauging this symmetry in depth. We unveil a rich gauging web relating the X-Cube model to symmetry protected topological (SPT) phases protected by a mix of subsystem and higher-form symmetries, subsystem symmetry fractionalization in the 3+1D Toric Code, and non-trivial extensions of topological symmetries by subsystem symmetries. Here, our work emphasizes the importance of topological symmetries in non-topological, geometric phases of matter.

Anyons↗

New multi-level codes over GF(q)

Set partitioning to multi-dimensional signal spaces over GF(q), particularly GF sup q-1(q) and GF sup q (q), and show how to construct both multi-level block codes and multi-level trellis codes over GF(q). Two classes of multi-level (n, k, d) block codes over GF(q) with block length n, number of information symbols k, and minimum distance d sub min greater than or = d, are presented. These two classes of codes use Reed-Solomon codes as component codes. They can be easily decoded as block length q-1 Reed-Solomon codes or block length q or q + 1 extended Reed-Solomon codes using multi-stage decoding. Many of these codes have larger distances than comparable q-ary block codes, as component codes. Low rate q-ary convolutional codes, work error correcting convolutional codes, and binary-to-q-ary convolutional codes can also be used to construct multi-level trellis codes over GF(q) or binary-to-q-ary trellis codes, some of which have better performance than the above block codes. All of the new codes have simple decoding algorithms based on hard decision multi-stage decoding.

Wu, Jiantian↗

New multilevel codes over GF(q)

Set partitioning to multi-dimensional signal spaces over GF(q), particularly GF sup q-1(q) and GF sup q (q), and show how to construct both multi-level block codes and multi-level trellis codes over GF(q). Two classes of multi-level (n, k, d) block codes over GF(q) with block length n, number of information symbols k, and minimum distance d sub min greater than or = d, are presented. These two classes of codes use Reed-Solomon codes as component codes. They can be easily decoded as block length q-1 Reed-Solomon codes or block length q or q + 1 extended Reed-Solomon codes using multi-stage decoding. Many of these codes have larger distances than comparable q-ary block codes, as component codes. Low rate q-ary convolutional codes, work error correcting convolutional codes, and binary-to-q-ary convolutional codes can also be used to construct multi-level trellis codes over GF(q) or binary-to-q-ary trellis codes, some of which have better performance than the above block codes. All of the new codes have simple decoding algorithms based on hard decision multi-stage decoding.

Wu, Jiantian↗

Performance comparison of combined ECC/RLL codes

In this paper, we present a performance comparison of several combined error correcting/run-lenth limited (ECC/RLL) codes created by concatenating a convolutional code with a run-length limited code. In each case, encoding and decoding are accomplished using a single trellis based on the combined code. Half of the codes under investigation use conventionally (d,k) run-length limited codes, where d is the minimum and k is the maximum allowable run of 0's between 1's. The other half of the combined codes use a special class of (d,k) codes known as distance preserving codes. These codes have the property that pairwise Hamming distances out of the (d,k) encoder are at least as large as the corresponding distances into the encoder (i.e., the codes preserve distance). Thus a combined code, created using a convolutional code concatenated with a distance preserving (d,k) code, will have a free distance (dfree) no smaller than the free distance of the original convolutional code. It should be noted that this does not hold if the (d,k) code was not distance preserving. A computer simulation is used to compare the performance of these two types of codes over the binary symmetric channel for various (d,k) constraints, rates, free distances, and numbers of states. Of particular interest for magnetic recording applications are codes with run-length constraints (1,3), (1,7), and (2,7).

French, C.↗

3-D GRMHD Simulations of Disk-Jet Coupling and Emission

We have performed a fully 3-D GRMHD simulation of jet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The simulation results show that a bipolar jet is initially created. At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface of the thickened (torus-like) disk. This evolution of disk-jet coupling suggests that the jet fades with a thickened accretion disk. Recently we have developed two new codes: 3 -D GRMHD: RelAtivIStic magnetoHydrodynamica1 sImulatioN (RAISHIN) code constructed by modern high-resolution shock-capturing (HRSC) techniques and 3-D GRPIC code. We have calculated free-free and synchrotron emission from the disks and jet/outflows obtained from our GRMHD simulations using a fully covariant radiative transfer formulation.

Nishikawa, K.-I.↗

Coupled neutronic-thermal-hydraulic simulations of the European SFR core

Within the European SFR (Sodium-cooled Fast Reactor) - Safety Measures Assessment and Research Tools (ESFR-SMART) project, steady-state coupled simulation of the ESFR core has been performed using several core analysis packages, with the objective of quantifying the coupling effect. Focus is on the fuel Doppler effect and coolant expansion effect. Standalone neutronics calculations in TRACE/PARCS (PSI), DYN3D (HZDR) and WIMS (Jacobs) showed superb agreement with the reference Serpent power distribution (root mean square - rms discrepancies of 1.3%, 1.5% and 0.7% respectively). Results for COUNTHER (CIEMAT) were also in reasonable agreement, but with a somewhat higher discrepancy of 3.7%. Temperature distributions from thermal-hydraulic calculations were also compared and are found to be in good agreement. The effect of Doppler and coolant density feedback on core power distribution was predicted by TRACE/PARCS, DYN3D and WIMS to be between 0.4% and 0.8% rms difference in assembly powers. Reactivity coefficients for perturbations in the inlet temperature, flow rate and core power were shown to be negative for these three codes, with values of roughly -0.5 pcm/C. degree, -0.3 pcm/C. degree and -3.5 pcm/C. degree respectively. A preliminary investigation of differential thermal expansion effects indicates that this may have a significant effect on core power distribution of a few %, greater than anticipated a priori and may warrant inclusion in coupled core analysis to ensure the accurate calculation of power distributions. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

LINFLUX-AE: A Turbomachinery Aeroelastic Code Based on a 3-D Linearized Euler Solver

This report describes the development and validation of LINFLUX-AE, a turbomachinery aeroelastic code based on the linearized unsteady 3-D Euler solver, LINFLUX. A helical fan with flat plate geometry is selected as the test case for numerical validation. The steady solution required by LINFLUX is obtained from the nonlinear Euler/Navier Stokes solver TURBO-AE. The report briefly describes the salient features of LINFLUX and the details of the aeroelastic extension. The aeroelastic formulation is based on a modal approach. An eigenvalue formulation is used for flutter analysis. The unsteady aerodynamic forces required for flutter are obtained by running LINFLUX for each mode, interblade phase angle and frequency of interest. The unsteady aerodynamic forces for forced response analysis are obtained from LINFLUX for the prescribed excitation, interblade phase angle, and frequency. The forced response amplitude is calculated from the modal summation of the generalized displacements. The unsteady pressures, work done per cycle, eigenvalues and forced response amplitudes obtained from LINFLUX are compared with those obtained from LINSUB, TURBO-AE, ASTROP2, and ANSYS.

Reddy, T. S. R.↗

Study of Two-Dimensional Compressible Non-Acoustic Modeling of Stirling Machine Type Components

A two-dimensional (2-D) computer code was developed for modeling enclosed volumes of gas with oscillating boundaries. such as Stirling machine components. An existing 2-D incompressible flow computer code, CAST, was used as the starting point for the project. CAST was modified to use the compressible non-acoustic Navier-Stokes equations to model an enclosed volume including an oscillating piston. The devices modeled have low Mach numbers and are sufficiently small that the time required for acoustics to propagate across them is negligible. Therefore. acoustics were excluded to enable more time efficient computation. Background information about the project is presented. The compressible non-acoustic flow assumptions are discussed. The governing equations used in the model are presented in transport equation format. A brief description is given of the numerical methods used. Comparisons of code predictions with experimental data are then discussed.

Roy C. Tew, Jr↗

Physics of Boundaries and their Interactions in Space Plasmas

This report describes the work done by SciberNet, Inc. during the month of October. We are working on the further refinement of the model used in our large-scale hybrid simulations of the magnetopause. Specifically, we are experimenting with several ways of modeling the effects of cold magnetospheric ions into our simulations. In addition, we are preparing two presentations for the upcoming Fall AGU highlighting the results of these simulations. We have also made progress in our development of a new kinetic linear code which we are using to study the linear properties of the Kelvin-Helmholtz instability at the magnetopause. We have extended the code from the electrostatic limit to the fully electromagnetic regime and are currently in the process of debugging and testing the code. Finally, we have made several test runs with our 2-D hybrid code for the magnetopause. The inflow-outflow boundary conditions are working properly. However, there are issues related to the setup and evolution of the original equilibrium that we are still trying to resolve. Finally, we are preparing several presentations for the upcoming Fall AGU.

Omidi, Nojan↗

Ablator Modeling: Why Not Much Has Changed Over the Past 45 plus Years

In the late 1950s, the earliest models describing the thermal response of ablative materials were based on the heat of ablation concept, which is an empirical approach that was reasonable for the types of materials of interest at that time. In the early-mid 60s the models were expanded to include pyrolysis since organic resin composites became the TPS materials of interest. However, surface recession was still predominantly modeled via empirical correlation. The development of the 1-D CMA finite difference code in the mid-late 60s introduced the thermochemical ablation approach for gas/surface interactions. Since that time investigators have developed finite volume and finite element codes, in 1-D, 2-D and 3-D, but the basic modeling has not evolved significantly. Models describing internal gas pressure due to pyrolysis, particle impact erosion, in-depth radiant transport, etc., have been added to address specific problems, but the fundamental modeling has not evolved. The reasons for this stagnation, as viewed by the author, will be described.

Laub, B.↗

Publicly Available Numerical Codes for Modeling the X-ray and Microwave Emissions from Solar and Stellar Activity

We have posted numerical codes on the Web for modeling the bremsstrahlung x-ray emission and the a gyrosynchrotron radio emission from solar and stellar activity. In addition to radiation codes, steady-state and time-dependent Fokker-Planck codes are provided for computing the distribution and evolution of accelerated electrons. A 1-D hydrodynamics code computes the response of the stellar atmosphere (chromospheric evaporation). A code for modeling gamma-ray line spectra is also available. On-line documentation is provided for each code. These codes have been developed for modeling results from the High Energy Solar Spectroscopic Imager (HESSI) along related microwave observations of solar flares. Comprehensive codes for modeling images and spectra of solar flares are under development. The posted codes can be obtained on NASA/Goddard's HESSI Web Site at http://hesperia.gsfc.nasa.gov/hessi/modelware.htm. This work is supported in part by the NASA Sun-Earth Connection Program.

Holman, Gordon D.↗

Simulations of plasma current induced by toroidal field ramping down on tokamaks

In tokamak experiments, the BT ramping down can drive off-axis and parallel inductive current. This approach leads to a decrease in q95 and an increase in normalized beta, βN. The off-axis inductive current also broadens the current profile. Typically, the 1.5D transport code is used to simulate the time evolution of the plasma current profile, which is based on the flux-surface-averaged Faraday’s law. The ONETWO code is one of those transport codes. However, this code cannot simulate the situation of the evolving toroidal field BT. This article proposes a new time-dependent model to take into account the BT ramping down situation. A modified formula of flux-surface-averaged Faraday’s law was derived to consider the effect of the BT ramping rate on the current evolution, and it was implemented in the ONETWO code. Then, the modified ONETWO code was used to simulate the BT ramping down experiment on DIII-D. The simulation result of the plasma current evolution with BT ramping down shows a broader current profile with smaller ohmic current induced by the poloidal field, compared to that without BT ramping down.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Wall interference assessment and corrections

Wind tunnel wall interference assessment and correction (WIAC) concepts, applications, and typical results are discussed in terms of several nonlinear transonic codes and one panel method code developed for and being implemented at NASA-Langley. Contrasts between 2-D and 3-D transonic testing factors which affect WIAC procedures are illustrated using airfoil data from the 0.3 m Transonic Cryogenic Tunnel and Pathfinder 1 data from the National Transonic Facility. Initial results from the 3-D WIAC codes are encouraging; research on and implementation of WIAC concepts continue.

Newman, P. A.↗

SAM-ML: Integrating data-driven closure with nuclear system code SAM for improved modeling capability

Advanced reactors often involve complicated thermal-fluid (T-F) phenomena. Modeling such phenomena with the traditional one-dimensional (1-D) system code is a challenging task. The System Analysis Module (SAM), a modern nuclear system code, has developed a coarse mesh multi-dimensional (multi-D) flow model to capture the spatial effect of T-F phenomena in advanced reactors. As a coarse mesh solver, constitutive relations are required for SAM's multi-D model for unresolved fine-scale physics, such as turbulence. Here this work presents a novel approach that integrates neural networks as data-driven closure for SAM's multi-D flow model. The data-driven closure is trained with fine-resolution data to ensure its accuracy while maintaining a coarse mesh setup to ensure its efficiency and consistency with SAM. We demonstrate the applicability of this SAM-ML capability in an open volume thermal stratification problem, where a neural network model serves as the eddy viscosity closure. A customized interface between the neural network and SAM is developed to ensure flexible and efficient data exchange. The SAM-ML results demonstrate superior performance compared to SAM's built-in zero-equation eddy viscosity closure. The case study shows that although the generalization capability of the data-driven closure still needs to be improved for different transient case or different geometric setup, SAM -ML demonstrates good potential for challenging simulation problems with improved accuracy and computational efficiency.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Nonlinear Waves and Particle Heating in the Auroral Ionosphere

In the two year period, the researchers have developed innovative 2-D simulation codes for modeling the interaction between electron beams and plasma waves and particles in the auroral ionosphere. These new simulations include kinetic wave particles as well as wave-wave interactions near and above approximately 500 km. The principal numerical model consists of partial-differential equations which evolve large amplitude magnetized 2-D Langmuir waves self-consistently with the electron distribution function. Progress was also made in developing 2-D Particle in Cell (PIC) codes.

Goldman, Martin↗