Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “F 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 1,531 records · Page 85

A simulation study of cold electron heating by loss cone instabilities

Results are presented of a computer simulation study of electron cyclotron harmonic waves; an electrostatic finite size particle simulation code is used. The initial electron velocity distribution function is modeled as a hot ring perpendicular to the ambient magnetic field and a cold Maxwellian component. Nonlinear cyclotron resonance is found to be the saturation mechanism of the unstable waves and the heating mechanism of the cold electrons. The results can be related to certain ionospheric phenomena.

Ashour-Abdalla, M.↗

Quasi-complementary codes - A new technique for MST radar sounding

The binary phase complementary code set has been used in stratospheric radar sounding. The quasi-complementary code technique presented here has nearly complementary properties and offers significant reduction of side lobes caused by imperfections in the radar transmitter. The new technique and its application are described; a simulation of the performance of the two techniques when certain imperfections are present in the transmitter is presented. The results are compared with the actual performance of the Arecibo 430-MHz radar transmitter.

Sulzer, M. P.↗

Paneling methods with vorticity effects and corrections for nonlinear compressibility

Supersonic panel methods and axisymmetric body-modeling singularity methods are presently combined with corrections for nonlinear flow phenomena to a complete missile, its airbreathing inlets, and wing-body combinations. The computer code LRCDM2 is used as an illustrative example of the methods in question. Attention is given to a preliminary method which employs panels to estimate additive drag and lift acting on supersonic rectangular inlets, as well as to the method used to correct off-body flowfields for the presence of a shock. Examples of missile applications of these methods with the appropriate nonlinear corrections are presented.

Dillenius, Marnix F. E.↗

Development and application of a reverse Monte Carlo radiative transfer code for rocket plume base heating

A reverse Monte Carlo radiative transfer code to predict rocket plume base heating is presented. In this technique rays representing the radiation propagation are traced backwards in time from the receiving surface to the point of emission in the plume. This increases the computational efficiency relative to the forward Monte Carlo technique when calculating the radiation reaching a specific point, as only the rays that strike the receiving point are considered.

Everson, John↗

Physics of the Isotopic Dependence of Galactic Cosmic Ray Fluence Behind Shielding

For over 25 years, NASA has supported the development of space radiation transport models for shielding applications. The NASA space radiation transport model now predicts dose and dose equivalent in Earth and Mars orbit to an accuracy of plus or minus 20%. However, because larger errors may occur in particle fluence predictions, there is interest in further assessments and improvements in NASA's space radiation transport model. In this paper, we consider the effects of the isotopic composition of the primary galactic cosmic rays (GCR) and the isotopic dependence of nuclear fragmentation cross-sections on the solution to transport models used for shielding studies. Satellite measurements are used to describe the isotopic composition of the GCR. Using NASA's quantum multiple-scattering theory of nuclear fragmentation (QMSFRG) and high-charge and energy (HZETRN) transport code, we study the effect of the isotopic dependence of the primary GCR composition and secondary nuclei on shielding calculations. The QMSFRG is shown to accurately describe the iso-spin dependence of nuclear fragmentation. The principal finding of this study is that large errors (plus or minus 100%) will occur in the mass-fluence spectra when comparing transport models that use a complete isotope grid (approximately 170 ions) to ones that use a reduced isotope grid, for example the 59 ion-grid used in the HZETRN code in the past, however less significant errors (less than 20%) occur in the elemental-fluence spectra. Because a complete isotope grid is readily handled on small computer workstations and is needed for several applications studying GCR propagation and scattering, it is recommended that they be used for future GCR studies.

Cucinotta, Francis A.↗

VLTI-MATISSE Chromatic Aperture-Synthesis Imaging of η Carinae’s Stellar Wind Across the Brα Line: Periastron Passage Observations in February 2020

Context: Eta Carinae is a highly eccentric, massive binary system (semimajor axis∼15.5 au) with powerful stellar winds and a phase-dependent wind-wind collision (WWC) zone. The primary star, ηCar A, is a luminous blue variable (LBV); the secondary, η Car B, is a Wolf-Rayet or O star with a faster but less dense wind. Aperture-synthesis imaging allows us to study the mass loss from the enigmatic LBV η Car. Understanding LBVs is a crucial step toward improving our knowledge about massive stars and their evolution. Aims: Our aim is to study the intensity distribution and kinematics of η Car’s WWC zone. Methods: Using the VLTI-MATISSE mid-infrared interferometry instrument, we perform Brα imaging of η Car’s distorted wind. Results: We present the first VLTI-MATISSE aperture-synthesis images of η Car A’s stellar wind in several spectral channels distributed across the Brα 4.052 μm line (spectral resolving power R∼960). Our observations were performed close to periastron passage in February 2020 (orbital phase∼14.0022). The reconstructed iso-velocity images show the dependence of the primary stellar wind on wavelength or line-of-sight (LOS) velocity with a spatial resolution of 6 mas (∼14 au). The radius of the faintest outer wind regions is∼26 mas (∼60 au). At several negative LOS velocities, the primary stellar wind is less extended to the northwest than in other directions. This asymmetry is most likely caused by the WWC. Therefore, we see both the velocity field of the undisturbed primary wind and the WWC cavity. In continuum spectral channels, the primary star wind is more compact than in line channels. A fit of the observed continuum visibilities with the visibilities of a stellar wind CMFGEN model (CMFGEN is an atmosphere code developed to model the spectra of a variety of objects) provides a full width at half maximum fit diameter of the primary stellar wind of 2.84±0.06 mas (6.54±0.14 au). We compare the derived intensity distributions with the CMFGEN stellar wind model and hydrodynamic WWC models.

G Weigelt↗

SOFE2023 Presentation: "Tritium Safety and Sensitivity Analysis of Tritium Extraction eXperiment (TEX) PbLi Loop Using MELCOR-TMAP"

The efficient extraction of tritium from the breeder blanket (BB) in a fusion reactor is vital for maintaining the tritium economy for a sustainable and self-sufficient fusion power plant. The vacuum permeator (VP) design is one of several promising technology platforms for extracting tritium from the BB with the advantages of being capable of operating continuously at high temperatures with the aim of achieving high tritium extraction efficiencies. The Tritium Extraction eXperiment (TEX), under construction at the Safety and Tritium Applied Research (STAR) Facility at Idaho National Laboratory (INL), is designed to test the viability of VP technology by introducing atomic hydrogen (H, D, and T) from a source permeator (SP) into a forced convection lead lithium eutectic (PbLi) loop and determining the hydrogen extraction rates from the vanadium VP. Presented here is a computational analysis of the hydrogen inventory and permeation behavior in TEX using the hydrogen transport and accident analysis code MELCOR-TMAP developed at INL. The results of a sensitivity analysis based on MELCOR-TMAP simulations of hydrogen extraction efficiencies from the TEX test section using varying PbLi flow rates, hydrogen source rates from the SP, and material property uncertainties are also reported. Operation scenarios for TEX in consideration of the analyses based on results from MELCOR-TMAP are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

GeantV: Results from the Prototype of Concurrent Vector Particle Transport Simulation in HEP

Full detector simulation was among the largest CPU consumers in all CERN experiment software stacks for the first two runs of the Large Hadron Collider. In the early 2010s, it was projected that simulation demands would scale linearly with increasing luminosity, with only partial compensation from increasing computing resources. The extension of fast simulation approaches to cover more use cases that represent a larger fraction of the simulation budget is only part of the solution, because of intrinsic precision limitations. The remainder corresponds to speeding up the simulation software by several factors, which is not achievable by just applying simple optimizations to the current code base. In this context, the GeantV R&D project was launched, aiming to redesign the legacy particle transport code in order to benefit from features of fine-grained parallelism, including vectorization and increased locality of both instruction and data. This paper provides an extensive presentation of the results and achievements of this R&D project, as well as the conclusions and lessons learned from the beta version prototype.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

hPIC2: A hardware-accelerated, hybrid particle-in-cell code for dynamic plasma-material interactions

The exascale era of high performance computing promises to bring the field of computational plasma physics ever closer to the goal of accurate multiscale modeling. Such computers will rely on hardware acceleration to offload work to dedicated components, notably general-purpose graphics processing units (GPUs). However, devices from different manufacturers require software to be written with different parallel programming models, greatly increasing the code maintenance burden of applications designed to perform on more than one such device. hPIC2 is a hybrid plasma simulation code developed with the Kokkos performance portability framework to target the architectures that will drive exascale computing for the foreseeable future. As a hybrid simulation code, hPIC2 investigates the simultaneous use of various plasma models on the same domain, at the same time. hPIC2 also optionally couples to RustBCA, which accurately models ion-material interactions using the binary collision approximation (BCA) method. In conclusion, hPIC2 therefore achieves scalable performance on a variety of computing architectures when simulating complex and diverse plasmas, particularly near plasma-material interfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Study of Alfvén eigenmode stability in Quasi-Poloidal Stellarator (QPS) plasma using a Landau closure model

The aim of this study is to analyze the linear stability of Alfvén eigenmodes (AE) in the QPS device heated by a tangential neutral beam injector (NBI). The analysis is performed using the gyro-fluid code FAR3d, that solves the reduced MHD equations for the thermal plasma coupled with moments of the kinetic equation for the energetic particles (EP). The AE stability is calculated in several operational regimes of the tangential NBI: EP β between 0.001 and 0.1, EP energy between 12 and 180 keV and different radial locations of the beam. The analysis is performed for vacuum and finite β equilibria as well as QPS configurations with two and three periods. The EP β threshold in the vacuum case is 0.001 and the AE frequency is lower as the energy of the EP population decreases. Toroidal Alfvén eigenmodes with f = 80–120 kHz and elliptical AE between f = 120–350 kHz are triggered between the middle-outer plasma region ($r/a \gt 0.5$). The AE stability improves in the simulations with finite β equilibria and three period configurations with respect to the vacuum case with two periods because the continuum gaps are slender, leading to a higher threshold of the EP β, above 0.03 for the AEs triggered by the helical mode families. Helical effects are not strong enough to destabilize Helical Alfvén eigenmodes, the AEs with the largest growth rates are triggered by the n = 1 and n = 2 toroidal families.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Program Verification for Extreme-Scale Applications (Final Scientific/Technical Report)

The project seeks to develop tools and techniques to help software developers verify the correctness and accuracy of their code. The target domain is scientific software of the kind widely used and developed in the Department of Energy research community, with a particular focus on "extreme scale" programs - those that are expected to involve possibly millions of parallel threads of execution. The report covers the University of Delaware contribution to the collaborative project. The project had a number of successful outcomes, especially regarding the development and extension of the CIVL software verification framework. CIVL is a verification tool for C or Fortran programs that use MPI, OpenMP, CUDA, and/or Pthreads for parallelization. CIVL went through vast improvements and extensions, and was successfully applied to a number of challenging codes. It found a subtle bug in the Devito PDE framework. It was able to verify the functional correctness of a conjugate gradient solver using a novel probabilistic technique with vanishingly small chance of error. CIVL was used very successfully in verification competitions, and the CIVL solutions to the competition challenges were presented and published. CIVL was also used successfully by other researchers on a computational chemistry kernel.

97 MATHEMATICS AND COMPUTING↗

An advanced imaging communication system for planetary exploration

An advanced imaging communication system (AICS) for planetary exploration is presented. The system offers 'end-to-end' information rate improvements of 3 to 5 times over existing systems in addition to extensive user flexibility to adapt his rate/fidelity priorities to fit a particular mission. AICS contains two major system elements. The first is a concatenated Reed-Solomon/Viterbi coded channel. This provides a powerful, yet practical, solution to the usual 'error vulnerability' problem associated with compressed data. The second major element is an extremely adaptive image data compression algorithm called RM2. The details of this algorithm as presently simulated are discussed in considerable detail. Used in conjunction with the virtually error free performance of the Reed-Solomon/Viterbi channel the stated AICS advantages are obtained.

Rice, R. F.↗

Numerical solution techniques for unsteady transonic aerodynamics problems

Basic concepts of finite difference solution techniques for unsteady transonic flows are presented. The hierarchy of mathematical forumulations that approximate the Navier-Stokes equations are reviewed. The basic concepts involved in constructing numerical algorthms to solve these formulations are given. Semi-implicit and implicit schemes are constructed and analyzed. The discussion focuses primarily on techniques for solving the low frequency transonic small disturbance equation. This is the simplest formulation that contains the essence of inviscid unsteady transonic flow physics. The low frequency formulation is emphasized here because codes based on this theory can be run in minutes of processor time on currently available computers. Furthermore, numerical techniques involved in solving this simple formulation also apply to the more complicated formulations. Extensions to these formulations are briefly described. An indication of the present capability for solving unsteady transonic flows is provided. Important areas of future research for the advancement of computational unsteady transonic aerodynamics are described.

Ballhaus, W. F.↗

Postbuckling behavior of selected flat stiffened graphite-epoxy panels loaded in compression

Results of an experimental study of the postbuckling behavior of selected flat stiffened graphite-epoxy panels loaded in compression are presented. The postbuckling response and failure characteristics of undamaged panels and panels damaged by low-speed impact are described. Each panel had four equally-spaced I-spaced stiffeners and 16- or 24-ply quasi-isotropic skins. Panels with three different stiffener spacings were tested. Some undamaged specimens supported as much as three times their initial buckling load before failing. Failure of all panels initiated in a skin-stiffener interface region. Analytical results obtained from a nonlinear general shell finite element analysis computer code correlate well with typical postbuckling test results up to failure. The analytical modeling detail necessary to predict accurately the response of a panel is described. Test results show that low-speed impact damage can reduce the postbuckling strength of a stiffened panel and that the skin-stiffener interface region is more sensitive to impact damage than the skin midway between stiffeners.

Starnes, J. H., Jr.↗

The application of a low-order panel method - program VSAERO to powerplant and airframe flow studies

Results from the application of VSAERO, a low-order panel method, to three practical aircraft configurations are presented. The Grumman 698-411 tilt-nacelle V/STOL model is analyzed with particular emphasis on the inlet pressures and the nacelle/fuselage interference effects. Excellent correlation with experiment is reported for the inlet pressure ratio and the inlet operational boundaries. Analysis of an inlet designed for a tilt-rotor/nacelle aircraft is presented. The code was used in a design environment for this configuration to determine an inlet geometry that maintained attached flow for three design flight conditions: hover, hover-transition and cruise. VSAERO is also used to examine the prop-slipstream induced loading for the Langley prop-fan configuration. The versatility and economy of the aerodynamic modeling program, VSAERO, is demonstrated.

Strash, D. J.↗

Codeless GPS Applications to Multi-Path: CGAMP

Cordless Global Positioning System (GPS) Applications to Multi-Path (CGAMP) is meeting the challenge of exploiting the L-band signals from the Global Positioning System (GPS) satellites for the measurement of the impulse response of radio transmission channels over space-Earth paths. This approach was originally suggested by E. K. Smith and has been pursued by J. Lemmon, without an affordable implementation being identifiable. In addition to the high cost of a suitable P code correlating GPS receiver, there is also the major impediment of the often announced Department of Defense policy of selective availability/anti-spoof (SA/AS) that clouds reliable access to the wideband (20 MHz) P channel of the GPS signals without cryptographic access. A technique proposed by MacDoran utilizes codeless methods for exploiting the P channel signals implemented by the use of a pair of antennas and cross correlation signal detection.

Macdoran, P. F.↗

On the symbolic manipulation and code generation for elasto-plastic material matrices

A computerized procedure for symbolic manipulations and FORTRAN code generation of an elasto-plastic material matrix for finite element applications is presented. Special emphasis is placed on expression simplifications during intermediate derivations, optimal code generation, and interface with the main program. A systematic procedure is outlined to avoid redundant algebraic manipulations. Symbolic expressions of the derived material stiffness matrix are automatically converted to RATFOR code which is then translated into FORTRAN statements through a preprocessor. To minimize the interface problem with the main program, a template file is prepared so that the translated FORTRAN statements can be merged into the file to form a subroutine (or a submodule). Three constitutive models; namely, von Mises plasticity, Drucker-Prager model, and a concrete plasticity model, are used as illustrative examples.

Chang, T. Y.↗

Modal Test/Analysis Correlation of Space Station Structures Using Nonlinear Sensitivity

The modal correlation problem is formulated as a constrained optimization problem for validation of finite element models (FEM's). For large-scale structural applications, a pragmatic procedure for substructuring, model verification, and system integration is described to achieve effective modal correlation. The space station substructure FEM's are reduced using Lanczos vectors and integrated into a system FEM using Craig-Bampton component modal synthesis. The optimization code is interfaced with MSC/NASTRAN to solve the problem of modal test/analysis correlation; that is, the problem of validating FEM's for launch and on-orbit coupled loads analysis against experimentally observed frequencies and mode shapes. An iterative perturbation algorithm is derived and implemented to update nonlinear sensitivity (derivatives of eigenvalues and eigenvectors) during optimizer iterations, which reduced the number of finite element analyses.

Gupta, Viney K.↗