Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Code of Accounts”

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 919 records · Page 51

Symmetry reduction of tensor networks in many-body theory: I. Automated symbolic evaluation of SU(2) algebra

Abstract The ongoing progress in (nuclear) many-body theory is accompanied by an ever-rising increase in complexity of the underlying formalisms used to solve the stationary Schrödinger equation. The associated working equations at play in state-of-the-art ab initio nuclear many-body methods can be analytically reduced with respect to angular-momentum, i.e. SU (2), quantum numbers whenever they are effectively employed in a symmetry-restricted context. The corresponding procedure constitutes a tedious and error-prone but yet an integral part of the implementation of those many-body frameworks. Indeed, this symmetry reduction is a key step to advance modern simulations to higher accuracy since the use of symmetry-adapted tensors can decrease the computational complexity by orders of magnitude. While attempts have been made in the past to automate the (anti-) commutation rules linked to Fermionic and Bosonic algebras at play in the derivation of the working equations, there is no systematic account to achieve the same goal for their symmetry reduction. In this work, the first version of an automated tool performing graph-theory-based angular-momentum reduction is presented. Taking the symmetry-unrestricted expressions of a generic tensor network as an input, the code provides their angular-momentum-reduced form in an error-safe way in a matter of seconds. Several state-of-the-art many-body methods serve as examples to demonstrate the generality of the approach and to highlight the potential impact on the many-body community.

Tichai, A.↗

An experimental study of pressures on 60 deg Delta wings with leading edge vortex flaps

An experimental study was conducted in the Virginia Tech Stability Wind Tunnel to determine surface pressures over a 60 deg sweep delta wing with three vortex flap designs. Extensive pressure data was collected to provide a base data set for comparison with computational design codes and to allow a better understanding of the flow over vortex flaps. The results indicated that vortex flaps can be designed which will contain the leading edge vortex with no spillage onto the wing upper surface. However, the tests also showed that flaps designed without accounting for flap thickness will not be optimum and the result can be oversized flaps, early flap vortex reattachment and a second separation and vortex at the wing/flap hinge line.

Marchman, J. F., III↗

A hybrid method for prediction of propeller performance

An analytical/computational code has been developed to predict the performance of advanced turboprop propellers. The method of solution is based on strip theory and conservation of linear and angular momenta applied to a control volume enclosing the propeller and its far wake. An iterative procedure is used to account for the induced axial and rotational velocities. Sectional aerodynamic coefficients are explicitly accounted for by thin airfoil theory and refinement of the results is obtained using a two-dimensional vortex panel method and Euler equations. The results obtained compare well with wind-tunnel data obtained for advanced propellers, such as the Hamilton Standard SR-2 and SR-7.

Takallu, M. A.↗

Transient Vibration Prediction for Rotors on Ball Bearings Using Load-dependent Non-linear Bearing Stiffness

Rolling-element bearing forces vary nonlinearly with bearing deflection. Thus an accurate rotordynamic transient analysis requires bearing forces to be determined at each step of the transient solution. Analyses have been carried out to show the effect of accurate bearing transient forces (accounting for non-linear speed and load dependent bearing stiffness) as compared to conventional use of average rolling-element bearing stiffness. Bearing forces were calculated by COBRA-AHS (Computer Optimized Ball and Roller Bearing Analysis - Advanced High Speed) and supplied to the rotordynamics code ARDS (Analysis of Rotor Dynamic Systems) for accurate simulation of rotor transient behavior. COBRA-AHS is a fast-running 5 degree-of-freedom computer code able to calculate high speed rolling-element bearing load-displacement data for radial and angular contact ball bearings and also for cylindrical and tapered roller beatings. Results show that use of nonlinear bearing characteristics is essential for accurate prediction of rotordynamic behavior.

Fleming, David P.↗

Code for Analyzing and Designing Spacecraft Power System Radiators

GPHRAD is a computer code for analysis and design of disk or circular-sector heat-rejecting radiators for spacecraft power systems. A specific application is for Stirling-cycle/linear-alternator electric-power systems coupled to radioisotope general-purpose heat sources. GPHRAD affords capabilities and options to account for thermophysical properties (thermal conductivity, density) of either metal-alloy or composite radiator materials.

Juhasz, Albert↗

Implementing Badhwar-O'Neill Galactic Cosmic Ray Model for the Analysis of Space Radiation Exposure

For the analysis of radiation risks to astronauts and planning exploratory space missions, accurate energy spectrum of galactic cosmic radiation (GCR) is necessary. Characterization of the ionizing radiation environment is challenging because the interplanetary plasma and radiation fields are modulated by solar disturbances and the radiation doses received by astronauts in interplanetary space are likewise influenced. A model of the Badhwar‐O'Neill 2011 (BO11) GCR environment, which is represented by GCR deceleration potential theta, has been derived by utilizing all of the GCR measurements from balloons, satellites, and the newer NASA Advanced Composition Explorer (ACE). In the BO11 model, the solar modulation level is derived from the mean international sunspot numbers with time‐delay, which has been calibrated with actual flight instrument measurements to produce better GCR flux data fit during solar minima. GCR fluxes provided by the BO11 model were compared with various spacecraft measurements at 1 AU, and further comparisons were made for the tissue equivalent proportional counters measurements at low Earth orbits using the high‐charge and energy transport (HZETRN) code and various GCR models. For the comparison of the absorbed dose and dose equivalent calculations with the measurements by Radiation Assessment Detector (RAD) at Gale crater on Mars, the intensities and energies of GCR entering the heliosphere were calculated by using the BO11 model, which accounts for time‐dependent attenuation of the local interstellar spectrum of each element. The BO11 model, which has emphasized for the last 24 solar minima, showed in relatively good agreement with the RAD data for the first 200 sols, but it was resulted in to be less well during near the solar maximum of solar cycle 24 due to subtleties in the changing heliospheric conditions. By performing the error analysis of the BO11 model and the optimization in reducing overall uncertainty, the resultant BO13 model corrects the fit at solar maxima as well as being accurate at solar minima. The BO13 model is implemented to the NASA Space Cancer Risk model for the assessment of radiation risks. Overall cumulative probability distribution of solar modulation parameters represents the percentile rank of the average interplanetary GCR environment, and the probabilistic radiation risks can be assessed for various levels of GCR environment to support mission design and operational planning for future manned space exploration missions.

Kim, Myung-Hee Y.↗

Magnetosphere-Ionosphere Coupling of Precipitated Electrons in Diffuse Aurora Driven by Time Domain Structures

Recent theoretical studies and the spacecraft conjugate observations between the Time History of Events and Macroscale Interactions during Substorms (THEMIS) mission and the low-altitude Enhanced Polar Outflow Probe (e-POP) spacecraft demonstrated the connections between broadband electrostatic fluctuations, the so called time domain structures (TDSs), and electron precipitation in the region of diffuse aurora. In this letter, we used the SuperThermal Electron Transport (STET) code to implement these theoretical and experimental results and analyze magnetosphere-ionosphere energy interplay of the precipitated electrons that are driven by TDSs. To put TDSs electron scattering processes in the context with other wave scattering activities, ECH and whistler chorus waves are taken into account. Integrated electron energy fluxes in the diffuse aurora driven by TDSs are comparable to those driven by whistler waves. We find that including Magnetosphere-Ionosphere coupling processes increases the integrated electron energy fluxes by more than two times than only driven by TDSs.

diffuse aurora↗

Large Eddy simulation of turbulence: A subgrid scale model including shear, vorticity, rotation, and buoyancy

The Reynolds numbers that characterize geophysical and astrophysical turbulence (Re approximately equals 10(exp 8) for the planetary boundary layer and Re approximately equals 10(exp 14) for the Sun's interior) are too large to allow a direct numerical simulation (DNS) of the fundamental Navier-Stokes and temperature equations. In fact, the spatial number of grid points N approximately Re(exp 9/4) exceeds the computational capability of today's supercomputers. Alternative treatments are the ensemble-time average approach, and/or the volume average approach. Since the first method (Reynolds stress approach) is largely analytical, the resulting turbulence equations entail manageable computational requirements and can thus be linked to a stellar evolutionary code or, in the geophysical case, to general circulation models. In the volume average approach, one carries out a large eddy simulation (LES) which resolves numerically the largest scales, while the unresolved scales must be treated theoretically with a subgrid scale model (SGS). Contrary to the ensemble average approach, the LES+SGS approach has considerable computational requirements. Even if this prevents (for the time being) a LES+SGS model to be linked to stellar or geophysical codes, it is still of the greatest relevance as an 'experimental tool' to be used, inter alia, to improve the parameterizations needed in the ensemble average approach. Such a methodology has been successfully adopted in studies of the convective planetary boundary layer. Experienc e with the LES+SGS approach from different fields has shown that its reliability depends on the healthiness of the SGS model for numerical stability as well as for physical completeness. At present, the most widely used SGS model, the Smagorinsky model, accounts for the effect of the shear induced by the large resolved scales on the unresolved scales but does not account for the effects of buoyancy, anisotropy, rotation, and stable stratification. The latter phenomenon, which affects both geophysical and astrophysical turbulence (e.g., oceanic structure and convective overshooting in stars), has been singularly difficult to account for in turbulence modeling. For example, the widely used model of Deardorff has not been confirmed by recent LES results. As of today, there is no SGS model capable of incorporating buoyancy, rotation, shear, anistropy, and stable stratification (gravity waves). In this paper, we construct such a model which we call CM (complete model). We also present a hierarchy of simpler algebraic models (called AM) of varying complexity. Finally, we present a set of models which are simplified even further (called SM), the simplest of which is the Smagorinsky-Lilly model. The incorporation of these models into the presently available LES codes should begin with the SM, to be followed by the AM and finally by the CM.

Canuto, V. M.↗

Adaptive mesh refinement for steady flows in Nek5000

Adaptive mesh refinement is performed in the framework of the spectral element method augmented by approaches to error estimation and control. The h-refinement technique is used for adapting the mesh, where selected grid elements are split by a quadtree (2D) or octree (3D) structure. Continuity between parent-child elements is enforced by high-order interpolation of the solution across the common faces. Parallel mesh partitioning and grid management respectively, are taken care of by the external libraries ParMETIS and p4est. Two methods are considered for estimating and controlling the error of the solution. The first error estimate is local and based on the spectral properties of the solution on each element. This method gives a local measure of the L-2-norm of the solution over the entire computational domain. The second error estimate uses the dual-weighted residuals method - it is based on and takes into account both the local properties of the solution and the global dependence of the error in the solution via an adjoint problem. The objective of this second approach is to optimize the computation of a given functional of physical interest. The simulations are performed by using the code Nek5000 and three steady-state test cases are studied: a two-dimensional lid-driven cavity at Re = 7, 500, a two-dimensional flow past a cylinder at Re = 40, and a three-dimensional lid-driven cavity at Re = 2,000 with a moving lid tilted by an angle of 30 degrees. Finally, the efficiency of both error estimators is compared in terms of refinement patterns and accuracy on the functional of interest. In the case of the adjoint error estimators, the trend on the error of the functional is shown to be correctly represented up to a multiplicative constant.

42 ENGINEERING↗

Enskog Transport Calculations with Magpie

he revised Enskog theory for mixtures has been implemented, as a library, in the LANL institutional thermochemical code Magpie. This report provides the necessary physical background and several representative examples of transport calculations with this library. The particular attention is paid to how the contributions of non-translational degrees of freedom (molecular vibrations, rotations, electronic excitations) to the energy transport can be consistently accounted for via the so called Eucken correction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Flexible Tool for Beam Induced Background Simulations at a Muon Collider

A Muon Collider represents a very interesting possibility for a future machine to explore the energy frontier in particle physics.However, to reach the needed luminosity, beam intensities of the order of $10^{9-11}$ muons per bunch are needed. In this context, the Beam Induced Background must be taken into account for its effects on the magnets and on the detector. Several optimizations can be conceived with the aim to mitigate them. In this view, it is of crucial importance to develop a flexible tool that allows to easily reconstruct the machine geometry in a Monte Carlo code, allowing to simulate in detail the interaction of muon decay products in the machine, while being able to change the machine optics itself to find the best configuration.In this contribution, a possible approach to such a purpose is presented, based on FLUKA for the Monte Carlo simulation and on LineBuilder for the geometry reconstruction.First results for the 1.5TeV machine optics developed by the MAP collaboration are discussed in this paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Analysis of coolant flow and heat transfer in the SSME HPOTP Number 4 bearing assembly

The PHOENICS code has been applied to simulate the flow of liquid oxygen through the number 4 ball bearing assembly of the Space Shuttle Main Engine High-Pressure Oxidizer Turbopump. The Body-Fitted Coordinate system capability of PHOENICS was utilized to create a geometrically accurate model. Effects accounted for by the model include the rotation rate of the calculation domain, viscous heating in the liquid oxygen and two-phase effects due to LOX boiling. A separation program was used to account for conduction within a ball bearing element which was coupled with the PHOENICS flow and heat transfer analysis. Solutions have been obtained for the velocity and temperature fields within the LOX coolant and the temperature within the ball element. The predicted ball temperatures indicate that the coolant is boiling in the region near the ball surface.

Owens, S. F.↗

Assessment and Improvement of RANS-based Transition Models based on Experimental Data of the Common Research Model with Natural Laminar Flow

Transition models based on auxiliary transport equations augmenting the Reynolds-averaged Navier-Stokes (RANS) framework often rely upon the correlations that were derived from a limited number of low-speed experiments and do not account for all of the transition mechanisms and/or their variation with the significant flow parameters. Available data from a recent experiment in the National Transonic Facility at the NASA Langley Research Center are used to assess the current transition modeling capability in NASA's OVERFLOW 2.2o code for a swept wing configuration at transonic cruise conditions. Specifically, the OVERFLOW solutions are used together with detailed stability analysis of the boundary layer flow over the new Common Research Model with Natural Laminar Flow (CRM-NLF) to evaluate the accuracy and the robustness of the transport-equation-based transition models, with the goal of proposing improvements that would help to strengthen the physical basis of these models for the important class of flows involving the combined effects of crossflow and flow compressibility. Results highlight the significant underprediction of the laminar flow extent within the inboard region of the wing, wherein the onset of transition may be attributed to a gradual amplification of Tollmien-Schlichting instabilities.

Boundary layer transition↗

Application of the NASA Multiscale Analysis Tool: Multiscale Integration and Interoperability

In order to demonstrate NASMAT’s multiscale operability, a series of illustrative examples will be presented that focus on the application of NASMAT to practical problems. First, the multiscale integration and data recursion is demonstrated by performing a multiscale analysis using only built-in micromechanics methods. NASMAT’s integration is then highlighted by running a multiscale analysis where an external finite element software calls NASMAT. In this case, at each integration point within the finite element model, a local NASMAT analysis is performed to account for failure behavior at the constituent scale. In a similar example, an external program is called from within NASMAT. This case would be relevant for a user wanting to implement an outside micromechanics technique. A combination of these examples is then presented to further illustrate the code’s flexibility when interfacing with outside codes in a multiscale framework. For all examples, data is presented using a custom-developed visualization tool. Additional potential use cases are also addressed. Finally, the plan for upcoming features and added capabilities is discussed.

NASMAT↗

A Model fot the Sources of the Slow Solar Wind

Models for the origin of the slow solar wind must account for two seemingly contradictory observations: the slow wind has the composition of the closed-field corona, implying that it originates from the continuous opening and closing of flux at the boundary between open and closed field. On the other hand, the slow wind also has large angular width, up to approx.60deg, suggesting that its source extends far from the open-closed boundary. We propose a model that can explain both observations. The key idea is that the source of the slow wind at the Sun is a network of narrow (possibly singular) open-field corridors that map to a web of separatrices and quasi-separatrix layers in the heliosphere. We compute analytically the topology of an open-field corridor and show that it produces a quasi-separatrix layer in the heliosphere that extends to angles far from the heliospheric current sheet. We then use an MHD code and MDI/SOHO observations of the photospheric magnetic field to calculate numerically, with high spatial resolution, the quasi-steady solar wind, and magnetic field for a time period preceding the 2008 August 1 total solar eclipse. Our numerical results imply that, at least for this time period, a web of separatrices (which we term an S-web) forms with sufficient density and extent in the heliosphere to account for the observed properties of the slow wind. We discuss the implications of our S-web model for the structure and dynamics of the corona and heliosphere and propose further tests of the model. Key words: solar wind - Sun: corona - Sun: magnetic topology

Antiochos, S. K.↗

A Model for the Sources of the Slow Solar Wind

Models for the origin of the slow solar wind must account for two seemingly contradictory observations: The slow wind has the composition of the closed-field corona, implying that it originates from the continuous opening and closing of flux at the boundary between open and closed field. On the other hand, the slow wind has large angular width, up to approximately 60 degrees, suggesting that its source extends far from the open-closed boundary. We propose a model that can explain both observations. The key idea is that the source of the slow wind at the Sun is a network of narrow (possibly singular) open-field corridors that map to a web of separatrices and quasi-separatrix layers in the heliosphere. We compute analytically the topology of an open-field corridor and show that it produces a quasi-separatrix layer in the heliosphere that extends to angles far front the heliospheric current sheet. We then use an MHD code and MIDI/SOHO observations of the photospheric magnetic field to calculate numerically, with high spatial resolution, the quasi-steady solar wind and magnetic field for a time period preceding the August 1, 2008 total solar eclipse. Our numerical results imply that, at least for this time period, a web of separatrices (which we term an S-web) forms with sufficient density and extent in the heliosphere to account for the observed properties of the slow wind. We discuss the implications of our S-web model for the structure and dynamics of the corona and heliosphere, and propose further tests of the model.

Antiochos, Spiro K.↗

Modelling type 1 quasar colours in the era of Rubin and Euclid

We construct a parametric SED model which is able to reproduce the average observed SDSS–UKIDSS–WISE quasar colours to within one-tenth of a magnitude across a wide range of redshift (0 < z < 5) and luminosity (–22 > M i > –29). This model is shown to provide accurate predictions for the colours of known quasars which are less luminous than those used to calibrate the model parameters, and also those at higher redshifts z > 5. Using a single parameter, the model encapsulates an up-to-date understanding of the intra-population variance in the rest-frame ultraviolet and optical emission lines of luminous quasars. At fixed redshift, there are systematic changes in the average quasar colours with apparent i-band magnitude, which we find to be well explained by the contribution from the host galaxy and our parametrization of the emission-line properties. By including redshift as an additional free parameter, the model could be used to provide photometric redshifts for individual objects. For the population as a whole we find that the average emission line and host-galaxy contributions can be well described by simple functions of luminosity which account for the observed changes in the average quasar colours across 18.1 < i AB < 21.5. Finally, we use these trends to provide predictions for quasar colours at the luminosities and redshifts which will be probed by the Rubin Observatory LSST and ESA-Euclid wide survey. The model code is applicable to a wide range of upcoming photometric and spectroscopic surveys, and is made publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

PORFLOW 6.42.9 Testing and Verification Document

This report is a continuation of the series of PORFLOW QA documents including (Aleman, 2007; Whiteside, 2007, 2008, 2010a,c,b, 2016b,a, 2017a,b). In this report we have updated the test cases where necessary to account for new features in PORFLOW and to ensure those features perform as described in the documentation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗