Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “V 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,045 records · Page 58

A Verification and Validation Approach for COMSOL Multiphysics to Support High Flux Isotope Reactor (HFIR)

Over the last several decades, many reactors have successfully been converted from high enriched uranium (HEU) to low enriched uranium (LEU) fuels in United States in support of its global non-proliferation objectives. Of the reactors slated for conversion, five high-performance research reactors (HPRRs) remain. The high-flux isotope reactor (HFIR) at Oak Ridge National Laboratory (ORNL) is one of the five. Conversion of HFIR requires the qualification of a new fuel. To aid in the conversion process as well as improve safety margins, COMSOL Multiphysics was chosen to support and eventually supplement the steady-state heat-transfer code (SSHTC). In order to use COMSOL as a supplementary tool, i.e. one used for safety basis calculations, the code needs to undergo a verification and validation process. Verification and validation are important processes for the use of any software; however, physics-based solvers have a unique challenge, in that the validation is much more rigorous and its domain is only applicable to the general scope of the problem. This document details the verification and validation of COMOSL for analysis of HFIR. This report is intended for two different users: (1) current and future users of COMSOL for HFIR and (2) for individuals interested in assessing the scope of the validation of COMSOL for HFIR. This report is divided into six sections. The first of which details the process of verification and validation. Section one provides an overview of key concepts from ASME V&V 10, 20, and 40 and how these concepts are used to provide a validation of any multiphysics software. Section two covers the approach needed in order to validate multiphysics software for a given HFIR analysis. Parts three through five split up the key concepts highlighted in section two. The final section provides concluding remarks which demonstrate the scope and limitations of the validation study and how these may be expanded and improved.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Explainable Machine Learning for Functional Data

Black-box machine learning models are recognized as useful tools for prediction applications, but the algorithmic complexity of some models causes interpretation challenges. Explainability methods have been proposed to provide insight into these models, but there is little research focused on supervised modeling with functional data inputs. We argue that, especially in applications of high consequence, it is important to explicitly model the functional dependence in a black-box analysis to not obscure or misrepresent patterns in explanations. As such, we propose the V ariable importance E xplainable E lastic S hape A nalysis (VEESA) pipeline for training supervised machine learning models with functional inputs. The pipeline is an analysis process that includes the data preprocessing, modeling, and post-hoc explanations. The preprocessing is done using elastic functional principal components analysis, which accounts for vertical and horizontal variability in functional data and, ultimately, allows for explanations in the original data space that identify the important functional variability without bias due to correlated variables. Here, we demonstrate the pipeline on two high-consequence applications: explosives classification for national security and inkjet printer identification in forensic science. The applications exhibit the VEESA pipeline’s ability to provide an understanding of the characteristics of the functional data useful for prediction. Code for implementing the pipeline is available in the veesa R package (and supplemental python code).

Elastic Shape Analysis↗

A study and evaluation of image analysis techniques applied to remotely sensed data

An analysis of phenomena causing nonlinearities in the transformation from Landsat multispectral scanner coordinates to ground coordinates is presented. Experimental results comparing rms errors at ground control points indicated a slight improvement when a nonlinear (8-parameter) transformation was used instead of an affine (6-parameter) transformation. Using a preliminary ground truth map of a test site in Alabama covering the Mobile Bay area and six Landsat images of the same scene, several classification methods were assessed. A methodology was developed for automatic change detection using classification/cluster maps. A coding scheme was employed for generation of change depiction maps indicating specific types of changes. Inter- and intraseasonal data of the Mobile Bay test area were compared to illustrate the method. A beginning was made in the study of data compression by applying a Karhunen-Loeve transform technique to a small section of the test data set. The second part of the report provides a formal documentation of the several programs developed for the analysis and assessments presented.

Atkinson, R. J.↗

A full potential inverse method based on a density linearization scheme for wing design

A mixed analysis inverse procedure based on the full potential equation in conservation form was developed to recontour a given base wing to produce density linearization scheme in applying the pressure boundary condition in terms of the velocity potential. The FL030 finite volume analysis code was modified to include the inverse option. The new surface shape information, associated with the modified pressure boundary condition, is calculated at a constant span station based on a mass flux integration. The inverse method is shown to recover the original shape when the analysis pressure is not altered. Inverse calculations for weakening of a strong shock system and for a laminar flow control (LFC) pressure distribution are presented. Two methods for a trailing edge closure model are proposed for further study.

Shankar, V.↗

An experimental investigation of the structural dynamics of a torsionally soft rotor in vacuum

An extensive data base of structural dynamic characteristics has been generated from an experimental program conducted on a torsionally soft two-bladed model helicopter rotor system. Measurements of vibratory strains for five modes of vibration were made at twenty-one locations on the two blades at speeds varying from 0 to 1000 RPM and for several combinations of precone, droop and flexure stiffness. Tests were conducted in vacuum under carefully controlled conditions using a unique excitation device with a system of piezoelectric crystals bonded to the blade surface near the root. Frequencies, strain mode shapes and dampings are extracted from the time histories and can be used to validate structural dynamics codes. The dynamics of the system are such that there is a clear tendency for the first torsion and second flap modes to couple within the speed range considered. Strain mode shapes vary significantly with speed and configuration. This feature is important in the calcualtion of aeroelastic instabilities. The tension axis tests confirmed that the modulus-weighted centroid for the nonhomogeneous airfoil is slightly off the mass centroid and validated previous static tests done to determine location of the tension axis.

Srinivasan, A. V.↗

Development of a unified constitutive model for an isotropic nickel base superalloy Rene 80

Accurate analysis of stress-strain behavior is of critical importance in the evaluation of life capabilities of hot section turbine engine components such as turbine blades and vanes. The constitutive equations used in the finite element analysis of such components must be capable of modeling a variety of complex behavior exhibited at high temperatures by cast superalloys. The classical separation of plasticity and creep employed in most of the finite element codes in use today is known to be deficient in modeling elevated temperature time dependent phenomena. Rate dependent, unified constitutive theories can overcome many of these difficulties. A new unified constitutive theory was developed to model the high temperature, time dependent behavior of Rene' 80 which is a cast turbine blade and vane nickel base superalloy. Considerations in model development included the cyclic softening behavior of Rene' 80, rate independence at lower temperatures and the development of a new model for static recovery.

Ramaswamy, V. G.↗

Experimental flutter boundaries with unsteady pressure distributions for the NACA 0012 Benchmark Model

The Structural Dynamics Division at NASA Langley Research Center has started a wind tunnel activity referred to as the Benchmark Models Program. The primary objective of the program is to acquire test data that will be useful for developing and evaluating aeroelastic type CFD codes currently in use of under development. This paper describes the progress achieved in testing the first model in the Benchmark Models Program. Experimental flutter boundaries are presented for a rigid semispan model (NACA 0012 airfoil section) mounted on a flexible mount system. Also, steady and unsteady pressure measurements taken at the flutter condition are presented. The pressure data were acquired over the entire model chord located at the 60 percent span station.

Rivera, Jose A.↗

Hypersonic Laminar Viscous Flow Past Spinning Cones at Angle of Attack

Computational results are presented for hypersonic viscous flow past spinning sharp and blunt cones of angle of attack, obtained with a parabolic Navier-Stokes marching code. The code takes into account the asymmetries in the flowfield resulting from spinning motion and computes the asymmetric shock shape, cross-flow and streamwise shear, heat transfer, cross-flow separation, and vortex structure. The Magnus force and moments are also computed. Comparisons are made with other theoretical analyses based on boundary-layer and boundary-region equations, and an anomaly is discovered in the displacement thickness contribution to the Magnus force when compared with boundary-layer results.

Agarwal, Ramesh↗

Systematic analysis of coding and noncoding DNA sequences using methods of statistical linguistics

We compare the statistical properties of coding and noncoding regions in eukaryotic and viral DNA sequences by adapting two tests developed for the analysis of natural languages and symbolic sequences. The data set comprises all 30 sequences of length above 50 000 base pairs in GenBank Release No. 81.0, as well as the recently published sequences of C. elegans chromosome III (2.2 Mbp) and yeast chromosome XI (661 Kbp). We find that for the three chromosomes we studied the statistical properties of noncoding regions appear to be closer to those observed in natural languages than those of coding regions. In particular, (i) a n-tuple Zipf analysis of noncoding regions reveals a regime close to power-law behavior while the coding regions show logarithmic behavior over a wide interval, while (ii) an n-gram entropy measurement shows that the noncoding regions have a lower n-gram entropy (and hence a larger "n-gram redundancy") than the coding regions. In contrast to the three chromosomes, we find that for vertebrates such as primates and rodents and for viral DNA, the difference between the statistical properties of coding and noncoding regions is not pronounced and therefore the results of the analyses of the investigated sequences are less conclusive. After noting the intrinsic limitations of the n-gram redundancy analysis, we also briefly discuss the failure of the zeroth- and first-order Markovian models or simple nucleotide repeats to account fully for these "linguistic" features of DNA. Finally, we emphasize that our results by no means prove the existence of a "language" in noncoding DNA.

NASA Discipline Number 14-10↗

Reducing costs of managing and accessing navigation and ancillary data by relying on the extensive capabilities of NASA's spice system

The SPICE system of navigation and ancillary data possesses a number of traits that make its use in modern space missions of all types highly cost efficient. The core of the system is a software library providing API interfaces for storing and retrieving such data as trajectories, orientations, time conversions, and instrument geometry parameters. Applications used at any stage of a mission life cycle can call SPICE APIs to access this data and compute geometric quantities required for observation planning, engineering assessment and science data analysis. SPICE is implemented in three different languages, supported on 20+ computer environments, and distributed with complete source code and documentation. It includes capabilities that are extensively tested by everyday use in many active projects and are applicable to all types of space missions - flyby, orbiters, observatories, landers and rovers. While a customer's initial SPICE adaptation for the first mission or experiment requires a modest effort, this initial effort pays off because adaptation for subsequent missions/experiments is just a small fraction of the initial investment, with the majority of tools based on SPICE requiring no or very minor changes.

ancillary data↗

Reply to "Comment on 'A Self-Consistent Model of the Interacting Ring Current Ions and Electromagnetic Ion Cyclotron Waves, Initial Results: Waves and Precipitation Fluxes' and 'Self-Consistent Model of the Magnetospheric Ring Current and Propagating Electromagnetic Ion Cyclotron Waves: Waves in Multi-Ion Magnetosphere' by Khazanov et al. et al."

It is well-known that the effects of electromagnetic ion cyclotron (EMIC) waves on ring current (RC) ion and radiation belt (RB) electron dynamics strongly depend on such particle/wave characteristics as the phase-space distribution function, frequency, wavenormal angle, wave energy, and the form of wave spectral energy density. The consequence is that accurate modeling of EMIC waves and RC particles requires robust inclusion of the interdependent dynamics of wave growth/damping, wave propagation, and[ particles. Such a self-consistent model is being progressively developed by Khazanov et al. [2002, 2006, 2007]. This model is based on a system of coupled kinetic equations for the RC and EMIC wave power spectral density along with the ray tracing equations. Thome and Home [2007] (hereafter referred to as TH2007) call the Khazanov et al. [2002, 2006] results into question in their Comment. The points in contention can be summarized as follows. TH2007 claim that: (1) "the important damping of waves by thermal heavy ions is completely ignored", and Landau damping during resonant interaction with thermal electrons is not included in our model; (2) EMIC wave damping due to RC O + is not included in our simulation; (3) non-linear processes limiting EMIC wave amplitude are not included in our model; (4) growth of the background fluctuations to a physically significantamplitude"must occur during a single transit of the unstable region" with subsequent damping below bi-ion latitudes,and consequently"the bounce averaged wave kinetic equation employed in the code contains a physically erroneous 'assumption". Our reply will address each of these points as well as other criticisms mentioned in the Comment. TH2007 are focused on two of our papers that are separated by four years. Significant progress in the self-consistent treatment of the RC-EMIC wave system has been achieved during those years. The paper by Khazanov et al. [2006] presents the latest version of our model, and in this Reply we refer mostly to this paper.

Khazanov, G. V.↗

A Requirements-Driven Optimization Method for Acoustic Liners Using Analytic Derivatives

More than ever, there is flexibility and freedom in acoustic liner design. Subject to practical considerations, liner design variables may be manipulated to achieve a target attenuation spectrum. But characteristics of the ideal attenuation spectrum can be difficult to know. Many multidisciplinary system effects govern how engine noise sources contribute to community noise. Given a hardwall fan noise source to be suppressed, and using an analytical certification noise model to compute a community noise measure of merit, the optimal attenuation spectrum can be derived using multidisciplinary systems analysis methods. In a previous paper on this subject, a method deriving the ideal target attenuation spectrum that minimizes noise perceived by observers on the ground was described. A simple code-wrapping approach was used to evaluate a community noise objective function for an external optimizer. Gradients were evaluated using a finite difference formula. The subject of this paper is an application of analytic derivatives that supply precise gradients to an optimization process. Analytic derivatives improve the efficiency and accuracy of gradient-based optimization methods and allow consideration of more design variables. In addition, the benefit of variable impedance liners is explored using a multi-objective optimization.

Aircraft Noise↗

A 2D Kaleidoscope of Electron Heat Fluxes Driven by Auroral Electron Precipitation

Electron heat flux is an important value for ionospheric space weather modeling networks. Utilizing the 2D array of Time History of Events and Macroscale Interactions during Substorms all-sky-imager (ASI) observations, Gabrielse et al. (2021, https://doi.org/10.3389/fphy.2021.744298) described a new method that estimates the auroral scale sizes of intense precipitating electron energy fluxes and their mean energies during two substorms on 16 February 2010. These parameters in combination with SuperThermal Electron Transport code were used to develop a new methodology to calculate electron thermal fluxes from data inputs in 2D during one of the substorms at 09:40:00 UT across Canada and Alaska. To test the effect of various precipitation lifetimes on electron heat flux values, boxcar averages ranging from 0 to 900 s were applied to the ASI data. These data are then combined with the newly developed kinetic simulation to determine the thermal fluxes associated with the observed diffuse and discrete precipitation.

Electron precipitation in aurora↗

Building an Integrated Ecosystem of Computational and Observational Facilities to Accelerate Scientific Discovery

Future scientific discoveries will rely on flexible ecosystems that incorporate modern scientific instruments, high performance computing resources, parallel distributed data storage, and performant networks across multiple, independent facilities. In addition to connecting physical resources, such an ecosystem presents many challenges in logistics and accessibility, especially in orchestrating computations and experiments that span across leadership computing systems and experimental instruments. Past efforts have typically been application-specific or limited to interfaces for computing resources. This paper proposes a general framework for integrating computation resources and instrument operations, addressing challenges in code development/execution, data staging and collection, software stack, control mechanisms, resource authorization and governance, and hardware integration. We also describe a demonstration use case wherein a Bayesian optimization algorithm running on an edge computing resource guides a scanning probe microscope to autonomously and intelligently characterize a material sample. This science edge ecosystem framework will provide a blueprint for federating multi-institutional, disparate resources and orchestrating scientific workflows across them to enable next-generation discoveries.

Somnath, Suhas↗

Erratum to: Measurements of higher-order cumulants of multiplicity and net-electric charge distributions in inelastic proton-proton interactions by NA61/SHINE

This Erratum replaces, due to a discovery of coding mistakes, the following quantities: κ 3 /κ 1 of the h + − h − distribution presented in Fig. 6 and Table 6, κ 4 listed in Table 4, and $\hat{C}$ 4 presented in Fig. 7 and Table 5. All mentioned figures and tables were updated.

High-Energy Particle Collision Data Analysis↗

Numerical solution of Space Shuttle Orbiter flow field

The supersonic, viscous laminar flow around the Space Shuttle Orbiter forebody has been computed with a parabolized Navier-Stokes code using a generalized coordinate transformation. The initial solution for the nose part of the Orbiter geometry was obtained with a three-dimensional time-dependent Navier-Stokes solver. It was necessary to employ a wind axis oriented coordinate system to obtain the initial solution with the time-dependent code. The generalized PNS technique used in this study allows the solution to be marched from the given initial data surface to any desired surface downstream. A grid point clustering scheme was employed to accurately describe the body shape by clustering points at the wing tip and at the wing body juncture. The computed heat transfer coefficients, pressure coefficients, and shock shapes are compared with the available experimental data for 0 degrees and 30 degrees angle of attack.

Venkatapathy, E.↗

Effects of cosmic rays on single event upsets

Assistance was provided to the Brookhaven Single Event Upset (SEU) Test Facility. Computer codes were developed for fragmentation and secondary radiation affecting Very Large Scale Integration (VLSI) in space. A computer controlled CV (HP4192) test was developed for Terman analysis. Also developed were high speed parametric tests which are independent of operator judgment and a charge pumping technique for measurement of D(sub it) (E). The X-ray secondary effects, and parametric degradation as a function of dose rate were simulated. The SPICE simulation of static RAMs with various resistor filters was tested.

Venable, D. D.↗

Probabilistic analysis of bladed turbine disks and the effect of mistuning

Probabilistic assessment of the maximum blade response on a mistuned rotor disk is performed using the computer code NESSUS. The uncertainties in natural frequency, excitation frequency, amplitude of excitation and damping are included to obtain the cumulative distribution function (CDF) of blade responses. Advanced mean value first order analysis is used to compute CDF. The sensitivities of different random variables are identified. Effect of the number of blades on a rotor on mistuning is evaluated. It is shown that the uncertainties associated with the forcing function parameters have significant effect on the response distribution of the bladed rotor.

Shah, A. R.↗