An Adams Guy Does the Runge-Kutta
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.
Explore the source record for details and available documents.
This paper describes an architecture for an autonomous deep space tracking station (DS-T).
A transonic unsteady aerodynamic and aeroelastic code called CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) was developed for application to realistic aircraft configurations. It permits the calculation of steady and unsteady flows about complete aircraft configurations for aeroelastic analysis of the flutter critical transonic speed range. The CAP-TSD code uses a time accurate approximate factorization algorithm for solution of the unsteady transonic small disturbance potential equation. An overview is given of the CAP-TSD code development effort along with recent algorithm modifications which are listed and discussed. Calculations are presented for several configurations including the General Dynamics 1/9th scale F-16C aircraft model to evaluate the algorithm and hence the reliability of the CAP-TSD code in general. Calculations are also presented for a flutter analysis of a 45 deg sweptback wing which agree well with the experimental data. Descriptions are presented of the CAP-TSD code and algorithm details along with results and comparisons which demonstrate the stability, accuracy, efficiency, and utility of CAP-TSD.
A double serially concatenated code with two interleavers consists of the cascade of an outer encoder, an interleaver permuting the outer codeword bits, a middle encoder, another interleaver permuting the middle codeword bits and an inner encoder whose input words are the permuted middle codewords.
We discuss theoretical and numerical aspects of gyrokinetics as a Lagrangian field theory when the field perturbation is introduced into the symplectic part. A consequence is that the field equations and particle equations of motion in general depend on the time derivatives of the field. The most well-known example is when the parallel vector potential is introduced as a perturbation, where a time derivative of the field arises only in the equations of motion, so an explicit equation for the fields may still be written. We will consider the conceptually more problematic case where the time-dependent fields appear in both the field equations and equations of motion, but where the additional term in the field equations is formally small. The conceptual issues were described by Burby (J. Plasma Phys., vol. 82 (3), 2016, 905820304): these terms lead to apparent additional degrees of freedom to the problem, so that the electric field now requires an initial condition, which is not required in low-frequency (Darwin) Vlasov–Maxwell equations. Also, the small terms in the Euler–Lagrange equations are a singular perturbation, and these two issues are interlinked. For well-behaved problems the apparent additional degrees of freedom are spurious, and the physically relevant solution may be directly identified. Because we needed to assume that the system is well behaved for small perturbations when deriving gyrokinetic theory, we must continue to assume that when solving it, and the physical solutions are thus the regular ones. The spurious nature of the singular degrees of freedom may also be seen by changing coordinate systems so the varying field appears only in the Hamiltonian. We then describe how methods appropriate for singular perturbation theory may be used to solve these asymptotic equations numerically. We then describe a proof-of-principle implementation of these methods for an electrostatic strong-flow gyrokinetic system; two basic test cases are presented to illustrate code functionality.
We use the global gyrokinetic stellarator code XGC-S to look at turbulence physics in geometry of the proposed quasi-axisymmetric stellarator, QUASAR. In this article, the XGC-S code is first verified for nonlinear simulations of ion temperature gradient-driven turbulence by a successful benchmark with the established nonlinear tokamak code XGC1 in a CYCLONE-like tokamak case. We then apply the XGC-S code to QUASAR for nonlinear turbulence physics. We find that the highly poloidally localized linear mode gives way to a more evenly distributed nonlinear turbulent structure. Furthermore, the radial extent of the mode increases and occupies a large fraction of the volume, despite the initially localized temperature profile, as the profiles relax.
Inertial confinement fusion (ICF) and high-energy density (HED) physics experiments experience complicated forcing for instability growth and mix due to the ubiquitous presence of multiple shocks interacting with perturbations on multiple material interfaces. One common driver of instability growth is successive shocks from the same direction. However, there is a severe lack of analytic work and modeling validation for same-sided successive shocks since they are extremely difficult to achieve with conventional (non-HED) drivers. Successive shocks access a large instability parameter space; idealized fluid theory [K. O. Mikaelian, Phys. Rev. A 31, 410 (1985)] predicts 15 different interface evolution scenarios for a sinusoidal perturbation. Growth becomes more complex for multi-mode, compressible HED systems. The Mshock campaign is the first experiment in any fluid regime to probe a wide portion of successive shock parameter space. This is enabled by our development of a hybrid direct/indirect drive platform capable of creating independently controllable successive shocks on the National Ignition Facility. These experiments have delivered the first data capable of rigorously challenging our models and their ability to accurately capture Richtmyer–Meshkov growth under successive shocks. Single-mode and two-mode experiments have successfully demonstrated the ability to access and control the various growth scenarios of the shocked interface, including re-inversion, freeze out, and continued growth. Simulations and theoretical modeling are shown to accurately capture the experimental observations in the linear growth phase, giving us confidence in our ICF/HED design codes.
This report reviews the modeling and simulation capabilities of Argonne National Laboratory’s DASSH code that is used in present reactor analysis activities. These capabilities will be used to establish the set of verification tasks necessary to verify DASSH for use on commercial projects. A similar approach was taken for the PERSENT, REBUS and DIF3D software packages. The DASSH program is a thermal analysis code designed to rapidly allow a reactor design engineer to obtain flow rates requirements that satisfy peak temperature constraints in the domain. DASSH is a follow-on development to the SE2-ANL software and SUPERENERGY-2 software that it is based upon. DASSH was designed to account for both neutron and gamma heating and is inherently connected to the GAMSOR part of the ARC suite of fast reactor analysis software. SE2-ANL is a developed piece of software from the 1980s while DASSH is a modern implementation with notable improvements in geometry handling. The most important upgrade of DASSH relative to SE2-ANL is that it can analyze multiple time points in a single run where SE2-ANL can only treat a single time point. This allows the user to understand the impact of and search the flow distribution for the entire operational period of a reactor design considering pressure drop, peak coolant and fuel temperatures, and thermal striping. DASSH has three input paths that have to be verified. The first input path builds the geometry and power distribution based upon the DIF3D model but ignores the gamma heating aspects of the problem. The second input path also builds the geometry from the DIF3D model but it takes the neutron and gamma heating distributions from GAMSOR. The third input path is to take the geometry and power distribution directly from user input (i.e. not coupled to DIF3D or GAMSOR). DASSH also has many built in correlations for material properties along with a user defined specification of the fuel, structure, and coolant properties. There are correlations for flow split, mixing, pressure drop, and heat transfer coefficients (subchannel rather than a direct methodology). In total, verification of DASSH will require an extensive testing to cover all possible user features of the software.
It is shown what choice there is in assigning output digits to transitions of binary rate 1/n code trellis so that the latter will correspond to a convolutional code. A new upper bound on free distance of rate 1/n convolutional codes is also derived, and the results obtained are used to determine the length of the largest input sequence that can conceivably result in an output whose weight is equal to the free distance of a code of rate 1/2.
Wind tunnel tests were conducted to determine the subsonic longitudinal aerodynamic characteristics of lifting configuration consisting of a 60 deg delta main wing with two smaller 60 deg delta wings (called sub-wings) attached underneath. The test was designed to determine the effects on lift, drag, and pitching moment due to various placement of the subwings in relation to the main wing. Test results indicate the increasing vertical separation between the main wing and the sub-wings produced the most significant results; a 23.1% increase in maximum lift coefficient, a reduction in drag coefficient at high lift coefficients, and an increase in longitudinal stability. Lateral separation of the sub-wings produced no significant changes. Placement of the sub-wings rearward increases the initial lift curve slope and maximum lift coefficient and also increase the longitudinal stability. Results of a computer study using a vortex lattice code supported the experimental conclusions.
The performance of a direct-detection optical communication system in which the laser transmitter sends short optical pulses of selected nonoverlapping center frequencies is analysed. This modulation format, in which a single light pulse is sent in one of M time slots at one of N optical center frequencies, is referred to as color coded optical pulse position modulation (CCPPM). The optimum energy-efficiency of this system, as measured by the channel capacity in nats per photon, exceeds that of ordinary optical pulse position modulation which uses a pulsed laser of fixed optical frequency. Reliable communication at optimal energy efficiency is easily achieved through the use of modest block length Reed-Solomon codes with the code words represented as CCPPM symbols.
The applicability of and the limits on the applicability of 11 near wall similarity laws characterizing three-dimensional turbulent boundary layer flows were determined. A direct force sensing local wall shear stress meter was used in both pressure-driven and shear-driven three-dimensional turbulent boundary layers, together with extensive mean velocity field and wall pressure field data. This resulted in a relatively large number of graphical comparisons of the predictive ability of 10 of these 11 similarity models relative to measured data over a wide range of flow conditions. Documentation of a complex, separated three-dimensional turbulent flow as a standard test case for evaluating the predictive ability of numerical codes solving such flows is presented.
The power of the STAGS family of programs was greatly enhanced. Members of the family include STAGS-C1 and RRSYS. As a result of improvements implemented, it is now possible to address the full collapse of a structural system, up to and beyond critical points where its resistance to the applied loads vanishes or suddenly changes. This also includes the important class of problems where a multiplicity of solutions exists at a given point (bifurcation), and where until now no solution could be obtained along any alternate (secondary) load path with any standard production finite element code.
The implementation of a partitioned numerical algorithm for determining the dynamic response of coupled structure/controller/estimator finite-dimensional systems is reviewed. The partitioned approach leads to a set of coupled first and second-order linear differential equations which are numerically integrated with extrapolation and implicit step methods. The present software implementation, ACSIS, utilizes parallel processing techniques at various levels to optimize performance on a shared-memory concurrent/vector processing system. A general procedure for the design of controller and filter gains is also implemented, which utilizes the vibration characteristics of the structure to be solved. Also presented are: example problems; a user's guide to the software; the procedures and algorithm scripts; a stability analysis for the algorithm; and the source code for the parallel implementation.
The objective of this research paper was to link two computer codes, PDCOMP (for Progressive Damage Analysis for Laminated Composites) and 3DIMPACT (for the prediction of the extent of delaminations in laminated composites resulting from point impact loads), in order to predict impact damage by taking into account local damage and material degradation and to estimate residual stiffness of composites after impact. The resulting graphs and analysis versus test results are presented along with the conclusive results of the codes' performances.
The VASIMR space propulsion development effort relies on a high power (greater than 10kW) helicon source to produce a dense flowing plasma (H, D and He) target for ion cyclotron resonance (ICR) acceleration of the ions. Subsequent expansion in an expanding magnetic field (magnetic nozzle) converts ion lunetic energy to directed momentum. This plasma source must have critical features to enable an effective propulsion device. First, it must ionize most of the input neutral flux of gas, thus producing a plasma stream with a high degree of ionization for application of ICR power. This avoids propellant waste and potential power losses due to charge exchange. Next, the plasma stream must flow into a region of high magnetic field (approximately 0.5 T) for efficient ICR acceleration. Third, the ratio of input power to plasma flux must be low, providing an energy per ion-electron pair approaching 100 eV. Lastly, the source must be robust and capable of very long life-times (years). In our helicon experiment (VX-10) we have measured a ratio of input gas to plasma flux near 100%. The plasma flows from the helicon region (B approximately 0.1 T) into a region with a peak magnetic field of 0.8 T. The energy input per ion-electron pair has been measured at 300 plus or minus 100 eV. Recent results at Oak Ridge National Laboratory (ORNL) show an enhanced efficiency mode of operation with a high power density, over 5 kW in a 5 cm diameter tube. Our helicon is presently 9 cm in diameter and operates up to 3.5 kW of input power. An upgrade to a power level of 10 kW is underway. Much of our recent work has been with a Boswell double-saddle antenna design. We are also converting the antenna design to a helical type. With these modifications, we anticipate an improvement in the ionization efficiency. This paper presents the results from scaling the helicon in the VX-10 device from 3.5 to 10 kW. We also compare the operation with a double-saddle to a helical antenna design. Finally, we discuss modeling of these configurations using ORNL's EMIR code.
The Coronagraph Performance Error Budget (CPEB) tool automates many of the key steps required to evaluate the scattered starlight contrast in the dark hole of a space-based coronagraph. The tool uses a Code V prescription of the optical train, and uses MATLAB programs to call ray-trace code that generates linear beam-walk and aberration sensitivity matrices for motions of the optical elements and line-of-sight pointing, with and without controlled fine-steering mirrors (FSMs). The sensitivity matrices are imported by macros into Excel 2007, where the error budget is evaluated. The user specifies the particular optics of interest, and chooses the quality of each optic from a predefined set of PSDs. The spreadsheet creates a nominal set of thermal and jitter motions, and combines that with the sensitivity matrices to generate an error budget for the system. CPEB also contains a combination of form and ActiveX controls with Visual Basic for Applications code to allow for user interaction in which the user can perform trade studies such as changing engineering requirements, and identifying and isolating stringent requirements. It contains summary tables and graphics that can be instantly used for reporting results in view graphs. The entire process to obtain a coronagraphic telescope performance error budget has been automated into three stages: conversion of optical prescription from Zemax or Code V to MACOS (in-house optical modeling and analysis tool), a linear models process, and an error budget tool process. The first process was improved by developing a MATLAB package based on the Class Constructor Method with a number of user-defined functions that allow the user to modify the MACOS optical prescription. The second process was modified by creating a MATLAB package that contains user-defined functions that automate the process. The user interfaces with the process by utilizing an initialization file where the user defines the parameters of the linear model computations. Other than this, the process is fully automated. The third process was developed based on the Terrestrial Planet Finder coronagraph Error Budget Tool, but was fully automated by using VBA code, form, and ActiveX controls.