SEARCH · Engineering Papers

Results for “surface 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,801 records · Page 100

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Automated bubble analysis of high-speed subcooled flow boiling images using U-net transfer learning and global optical flow

Capturing and analyzing the bubble dynamics is crucial to improving the understanding of boiling heat transfer mechanisms and predicting boiling heat transfer coefficient and boiling crisis. High speed video (HSV) imaging has been used for decades towards this end. Still, there is no universal approach to quantitatively analyze bubble dynamics from HSV images. In this study, we propose a data-driven post-processing approach to segment, track, and identify wall-attached vapor bubbles from HSV images of the boiling process in subcooled flow conditions. Firstly, we employ a transfer learning framework with a U-Net-based convolution neural network (CNN) architecture to detect and segment bubbles in HSV images of diverse contrast and surface texture using very little data (e.g., 10 images) for training. Then, we evaluate the trained CNN model with 100 ground-truth images, and the validation results show that the model accuracy and precision in detecting the optical footprint of bubbles are higher than 90%. Finally, we suggest a criterion to identify a condensing bubble based on the divergence of the bubble displacement, which is calculated from sequential segmented bubble images using a global optical flow code. Using this combination of machine learning and optical flow, we can identify nucleation sites and track the growth of bubbles nucleating at each site to quantify nucleation site density, nucleation frequency, and other fundamental boiling parameters. The proposed system is validated using results obtained on a special heater, which enables both infrared (IR) thermometry and HSV imaging on a metallic surface. We compare the fundamental boiling parameters obtained by the two different diagnostics. The results show good agreement. In conclusion, the difference between the measurements of nucleation site density, averaged nucleation frequency, and averaged growth time performed with the two techniques is always within ± 20% and mostly ± 10% of the values measured with IR thermometry.

42 ENGINEERING↗

NASA NTRS1982

Additional extensions to the NASCAP computer code, volume 2

Particular attention is given to comparison of the actural response of the SCATHA (Spacecraft Charging AT High Altitudes) P78-2 satellite with theoretical (NASCAP) predictions. Extensive comparisons for a variety of environmental conditions confirm the validity of the NASCAP model. A summary of the capabilities and range of validity of NASCAP is presented, with extensive reference to previously published applications. It is shown that NASCAP is capable of providing quantitatively accurate results when the object and environment are adequately represented and fall within the range of conditions for which NASCAP was intended. Three dimensional electric field affects play an important role in determining the potential of dielectric surfaces and electrically isolated conducting surfaces, particularly in the presence of artificially imposed high voltages. A theory for such phenomena is presented and applied to the active control experiments carried out in SCATHA, as well as other space and laboratory experiments. Finally, some preliminary work toward modeling large spacecraft in polar Earth orbit is presented. An initial physical model is presented including charge emission. A simple code based upon the model is described along with code test results.

Stannard, P. R.↗

NASA NTRS2007

Simulation of a Synthetic Jet in Quiescent Air Using TLNS3D Flow Code

Although the actuator geometry is highly three-dimensional, the outer flowfield is nominally two-dimensional because of the high aspect ratio of the rectangular slot. For the present study, this configuration is modeled as a two-dimensional problem. A multi-block structured grid available at the CFDVAL2004 website is used as a baseline grid. The periodic motion of the diaphragm is simulated by specifying a sinusoidal velocity at the diaphragm surface with a frequency of 450 Hz, corresponding to the experimental setup. The amplitude is chosen so that the maximum Mach number at the jet exit is approximately 0.1, to replicate the experimental conditions. At the solid walls zero slip, zero injection, adiabatic temperature and zero pressure gradient conditions are imposed. In the external region, symmetry conditions are imposed on the side (vertical) boundaries and far-field conditions are imposed on the top boundary. A nominal free-stream Mach number of 0.001 is imposed in the free stream to simulate incompressible flow conditions in the TLNS3D code, which solves compressible flow equations. The code was run in unsteady (URANS) mode until the periodicity was established. The time-mean quantities were obtained by running the code for at least another 15 periods and averaging the flow quantities over these periods. The phase-locked average of flow quantities were assumed to be coincident with their values during the last full time period.

Vatsa, Veer N.↗

NASA NTRS2003

Surface Segregation Energies of BCC Binaries from Ab Initio and Quantum Approximate Calculations

We compare dilute-limit segregation energies for selected BCC transition metal binaries computed using ab initio and quantum approximate energy method. Ab initio calculations are carried out using the CASTEP plane-wave pseudopotential computer code, while quantum approximate results are computed using the Bozzolo-Ferrante-Smith (BFS) method with the most recent parameterization. Quantum approximate segregation energies are computed with and without atomistic relaxation. The ab initio calculations are performed without relaxation for the most part, but predicted relaxations from quantum approximate calculations are used in selected cases to compute approximate relaxed ab initio segregation energies. Results are discussed within the context of segregation models driven by strain and bond-breaking effects. We compare our results with other quantum approximate and ab initio theoretical work, and available experimental results.

Good, Brian S.↗

NASA NTRS2010

Performance Analysis of Apollo Navigational Starter Routine

The focus of this project is to recreate and analyze the effectiveness of the original Apollo Starter Routine (ASR) which was used to generate the state vector of the Apollo spacecraft based on a series of radiometric observations. The original Apollo navigation software is unavailable in a modern programming language and the original coding has not been preserved. This necessitates its recreation using the original software documentation. Space Shuttle navigation software does not typically use the ASR or an algorithm like it since the Shuttle s state vector is easily deduced from GPS information or other sources. However, this tactic will be ineffective when trying to determine the state vector of a craft approaching, departing or in orbit around the Moon since the GPS network faces the surface of the Earth, not outer space. The recreation of the ASR from the original documentation is therefore vital as a simulation baseline for the navigation software under development for the Constellation program. The algorithms that make up the ASR will be extracted from the original documentation and adapted for and then implemented in a modern programming language; the majority of it will be coded in Matlab. The ASR s effectiveness will then be tested using simulated tracking data. The ability of the ASR to handle realistically noisy data and the accuracy with which it generates state vectors were analyzed. The ASR proved to be robust enough to process data with range and angle noise as large as 10,000 meters and 10(exp -6) radians together and 300,000 meters and 5x10(exp -4) radians separately at Lunar distances. The ASR was able to handle marginally more noise at distances closer to the Earth where the angle noise was less significant. The ASR is capable of effectively processing 40-80 data points gathered at a rate of one per 20 seconds at close Earth orbit and up to 28-40 data points gathered at a rate of one per minute at distant Earth orbit and Lunar orbit.

Ivanov, Stoyan I.↗

NASA NTRSDate not supplied

Lunar Nuclear Reactor Neutron Fluence and Gamma Dose Estimates

NASA aims to deploy a fission reactor on the lunar surface by 2030 to generate 100 kWe of power for at least 10 years without maintenance or refueling. To assess system and subsystem survivability according to those requirements, a model of the reactor pallet was created in the radiation transport code MCNP to determine the neutron fluences and gamma doses that the materials and systems would be exposed to. The model was based on the 2025 government design. The model indicated that the control drum motors, Brayton cycle engines, and pipes and valves would be exposed to neutron fluences of 1×10 15 – 1×10 16 n/cm 2 per year and gamma doses ranging from 10 to 500 Mrad per year. The electronics box housing the instrumentation and control elements is expected to receive 1×10 14 – 1×10 15 n/cm 2 per year and gamma doses in the range of 1–5 Mrad per year. These doses necessitate that the control drum motors, Brayton cycle engines, valve components, control electronics, and other components be radiation-hardened (>10 16 n/cm 2 and >10 Mrad gamma dose) or additionally shielded to survive the 10-year lunar mission.

Radiation-Hardened Materials↗

DOE OSTI2020

Soil microbial function and fungal community data along an ectomycorrhizal basal area gradient, Vermont, May 2019

Mycorrhizal fungi can affect soil organic matter cycling through several mechanisms including priming, nutrient competition, and direct enzyme production. Differences in nutrient foraging strategies between ectomycorrhizal (EcM) and arbuscular mycorrhizal (AM) fungi produce divergent belowground dynamics: where EcM can take up organic nitrogen and directly break down SOM by producing enzymes, AM fungi are limited to scavenging mineral N. EcM-associated tree species also have leaf litter with relatively higher ratios of carbon to nitrogen (C:N), and belowground saprotrophic communities more dominated by fungi. Consequently, free-living microbes in EcM-dominated soils should experience nitrogen limitation, with subsequent increases in enzyme production and decreased carbon use efficiency (CUE). However, the relative importance of the effects of substrate quality and fungal community composition on enzyme production and CUE are unclear. To assess this distinction, we sampled the organic horizon and 10 cm of the mineral horizon in northern temperate forest soils along a gradient of EcM dominance (see SCxTS_DBH.csv for full tree species information). We characterized fungal community composition by measuring EcM relative abundances from extracted fungal DNA (To run raw ASV data through R code, use SCxTS_ASV_taxonomy.csv, and R code ITS_processing.Rmd; to see processed data with EcM relative abundance, use file SCxTS_sum_afitch_Feb18.20.csv) and the fungal to bacterial ratios (see data file SCxTS_plfa.csv) from phospholipid fatty acid analysis. Soil microbial functions were measured as potential activities of five hydrolytic and two oxidative enzymes ( To run raw data through R code us SCxTS_metadata.csv, hydrolytic_oxidative_raw.csv, MUBraw.csv, Enzyme_soilweight.csv, and incubationtime.csv, and R code SCxTS_enzymes.Rmd; to use Enzyme data ready to run through statistics code, use file Enzyme.data.csv) and microbial carbon use efficiency (see file SCxTS_CUE.csv and for full calculations file 18O_CUEcalcs_SCxTS_2019.xlsx). We assessed soil substrate quality as the soil carbon:nitrogen ratio (See file Organic_CN.csv). Statistical analysis can be run using R code from file SCxTS_analysis.Rmd.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Slender-body approach for computing second-order wave loads in the frequency domain

This work presents a slender-body approach to evaluate the second-order wave loads acting on a floating structure in the frequency domain. The approach is in the same spirit as the common use of Morison’s equation to approximate the wave loads without solving the radiation/diffraction problem. To do so, we employ Rainey’s equation, which can be seen as an extension of the inertial part of Morison’s equation to include nonlinear effects. We introduce modifications to Rainey’s formulation in order to evaluate wave kinematics at the mean body position instead of the original approach of considering instantaneous displacements. We also propose a simple approximation to partially account for wave scattering effects on the second-order loads based on the analytical solution of a surface-piercing bottom-mounted vertical circular cylinder. Though limited to structures composed of cylinders, this slender-body approach is orders of magnitude faster than computing second-order wave coefficients with a radiation/diffraction code. We implemented this approach for difference-frequency (slow drift) loads in an open-source frequency-domain floating wind turbine model. We present comparisons against results obtained with radiation/diffraction theory for three reference floating wind turbine designs: the OC3-Hywind spar, the OC4-DeepCwind semisubmersible, and the VolturnUS-S semisubmersible. In general, the results show that the proposed slender-body approach with the correction to approximate wave scattering effects provides useful estimations of the difference-frequency wave loads and the resulting motions of the floater.

17 WIND ENERGY↗

DOE OSTI2023

Real-space solution to the electronic structure problem for nearly a million electrons

We report a Kohn–Sham density functional theory calculation of a system with more than 200 000 atoms and 800 000 electrons using a real-space high-order finite-difference method to investigate the electronic structure of large spherical silicon nanoclusters. Our system of choice was a 20 nm large spherical nanocluster with 202 617 silicon atoms and 13 836 hydrogen atoms used to passivate the dangling surface bonds. To speed up the convergence of the eigenspace, we utilized Chebyshev-filtered subspace iteration, and for sparse matrix–vector multiplications, we used blockwise Hilbert space-filling curves, implemented in the PARSEC code. For this calculation, we also replaced our orthonormalization + Rayleigh–Ritz step with a generalized eigenvalue problem step. We utilized all of the 8192 nodes (458 752 processors) on the Frontera machine at the Texas Advanced Computing Center. We achieved two Chebyshev-filtered subspace iterations, yielding a good approximation of the electronic density of states. Our work pushes the limits on the capabilities of the current electronic structure solvers to nearly 106 electrons and demonstrates the potential of the real-space approach to efficiently parallelize large calculations on modern high-performance computing platforms.

Chemistry↗

DOE OSTI2021

R Coronae Borealis Star Evolution: Simulating 3D Merger Events to 1D Stellar Evolution Including Large-scale Nucleosynthesis

R Coronae Borealis (RCB) stars are rare hydrogen-deficient carbon-rich variable supergiants thought to be the result of dynamically unstable white dwarf mergers. In this article, we attempt to model RCB stars through all the relevant timescales by simulating a merger event in Octo-tiger, a 3D adaptive mesh refinement (AMR) hydrodynamics code, and mapping the post-merger object into MESA, a 1D stellar evolution code. We then post-process the nucleosynthesis on a much larger nuclear reaction network to study the enhancement of s-process elements. We present models that match observations or previous studies in most surface abundances, isotopic ratios, early evolution, and lifetimes. We also observe similar mixing behavior to previous modeling attempts that result in the partial He-burning products visible on the surface in observations. However, we do note that our subsolar models lack any enhancement in s-process elements, which we attribute to a lack of hydrogen in the envelope. We also find that the 16 O/ 18 O isotopic ratio is very sensitive to initial hydrogen abundance and increases outside of the acceptable range with a hydrogen mass fraction greater than 10 -4 .

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1975

Calculation of three-dimensional compressible laminar and turbulent boundary flows. Three-dimensional compressible boundary layers of reacting gases over realistic configurations

A three-dimensional boundary-layer code was developed for particular application to realistic hypersonic aircraft. It is very general and can be applied to a wide variety of boundary-layer flows. Laminar, transitional, and fully turbulent flows of compressible, reacting gases are efficiently calculated by use of the code. A body-oriented orthogonal coordinate system is used for the calculation and the user has complete freedom in specifying the coordinate system within the restrictions that one coordinate must be normal to the surface and the three coordinates must be mutually orthogonal.

Kendall, R. M.↗

NASA NTRS1980

Finite difference grid generation by multivariate blending function interpolation

The General Interpolants Method (GIM) code which solves the multidimensional Navier-Stokes equations for arbitrary geometric domains is described. The geometry module in the GIM code generates two and three dimensional grids over specified flow regimes, establishes boundary condition information and computes finite difference analogs for use in the GIM code numerical solution module. The technique can be classified as an algebraic equation approach. The geometry package uses multivariate blending function interpolation of vector-values functions which define the shapes of the edges and surfaces bounding the flow domain. By employing blending functions which conform to the cardinality conditions the flow domain may be mapped onto a unit square (2-D) or unit cube (3-D), thus producing an intrinsic coordinate system for the region of interest. The intrinsic coordinate system facilitates grid spacing control to allow for optimum distribution of nodes in the flow domain.

Anderson, P. G.↗

NASA NTRS1986

Computation of transonic flow about helicopter rotor blades

An inviscid, nonconservative, three-dimensional full-potential flow code, ROT22, has been developed for computing the quasi-steady flow about a lifting rotor blade. The code is valid throughout the subsonic and transonic regime. Calculations from the code are compared with detailed laser velocimeter measurements made in the tip region of a nonlifting rotor at a tip Mach number of 0.95 and zero advance ratio. In addition, comparisons are made with chordwise surface pressure measurements obtained in a wind tunnel for a nonlifting rotor blade at transonic tip speeds at advance ratios from 0.40 to 0.50. The overall agreement between theoretical calculations and experiment is very good. A typical run on a CRAY X-MP computer requires about 30 CPU seconds for one rotor position at transonic tip speed.

Arieli, R.↗

NASA NTRS1987

Use of gamma ray radiation to parallel the plates of a Fabry-Perot interferometer

The use of gamma radiation to parallel the plates of a Fabry-Perot etalon is examined. The method for determining the etalon parallelism, and the procedure for irradiating the posts are described. Changes in effective gap for the etalon over the surface are utilized to measure the parallelism of the Fabry-Perot etalon. An example in which this technique is applied to an etalon of fused silica plates, which are 132 mm in diameter and coded with zinc sulfide and cryolite, with Zerodur spaces 2 cm in length. The effect of the irradiation of the posts on the thermal performance of the etalon is investigated.

Skinner, Wilbert R.↗

NASA NTRS1990

Numerical simulations of supersonic flow through oscillating cascade sections

A finite difference code was developed for modeling inviscid, unsteady supersonic flow by solution of the compressible Euler equations. The code uses a deforming grid technique to capture the motion of the airfoils and can model oscillating cascades with any arbitrary interblade phase angle. A flat plate cascade is analyzed, and results are compared with results from a small perturbation theory. The results show very good agreement for both the unsteady pressure distributions and the integrated force predictions. The reason for using the numerical Euler code over a small perturbation theory is the ability to model real airfoils that have thickness and camber. Sample predictions are presented for a cascade of loaded airfoils and show appreciable differences in the unsteady surface pressure distributions when compared with the flat plate results.

Huff, Dennis L.↗

NASA NTRS1993

Mathmatical modeling for diffractive optics

We consider a 'diffractive optic' to be a biperiodic surface separating two half-spaces, each having constant constitutive parameters; within a unit cell of the periodic surface and across the transition zone between the two half-spaces, the constitutive parameters can be a continuous, complex-valued function. Mathematical models for diffractive optics have been developed, and implemented as numerical codes, both for the 'direct' problem and for the 'inverse' problem. In problems of the 'direct' class, the diffractive optic is specified, and the full set of Maxwell's equations is cast in a variational form and solved numerically by a finite element approach. This approach is well-posed in the sense that existence and uniqueness of the solution can be proved and specific convergence conditions can be derived. An example of a metallic grating at a Wood anomaly is presented as a case where other approaches are known to have convergence problems. In problems of the 'inverse' class, some information about the diffracted field (e.g., the far-field intensity) is given, and the problem is to find the periodic structure in some optimal sense. Two approaches are described: phase reconstruction in the far-field approximation; and relaxed optimal design based on the Helmholtz equation. Practical examples are discussed for each approach to the inverse problem, including array generators in the far-field case and antireflective structures for the relaxed optimal design.

Dobson, David↗

NASA NTRS1995

Minimization of the Vibration Energy of Thin-Plate Structures and the Application to the Reduction of Gearbox Vibration

While the vibration analysis of gear systems has been developed, a systematic approach to the reduction of gearbox vibration has been lacking. The technique of reducing vibration by shifting natural frequencies is proposed here for gearboxes and other thin-plate structures using the theories of finite elements, modal analysis, and optimization. A triangular shell element with 18 degrees of freedom is developed for structural and dynamic analysis. To optimize, the overall vibration energy is adopted as the objective function to be minimized at the excitation frequency by varying the design variable (element thickness) under the constraint of overall constant weight. Modal analysis is used to determine the sensitivity of the vibration energy as a function of the eigenvalues and eigenvectors. The optimum design is found by the gradient projection method and a unidimensional search procedure. By applying the computer code to design problems for beams and plates, it was verified that the proposed method is effective in reducing vibration energy. The computer code is also applied to redesign the NASA Lewis gear noise rig test gearbox housing. As one example, only the shape of the top plate is varied, and the vibration energy levels of all the surfaces are reduced, yielding an overall reduction of 1/5 compared to the initial design. As a second example, the shapes of the top and two side plates are varied to yield an overall reduction in vibration energy of 1/30.

Inoue, Katsumi↗

NASA NTRS2003

Flutter and Divergence Analysis using the Generalized Aeroelastic Analysis Method

The Generalized Aeroelastic Analysis Method (GAAM) is applied to the analysis of three well-studied checkcases: restrained and unrestrained airfoil models, and a wing model. An eigenvalue iteration procedure is used for converging upon roots of the complex stability matrix. For the airfoil models, exact root loci are given which clearly illustrate the nature of the flutter and divergence instabilities. The singularities involved are enumerated, including an additional pole at the origin for the unrestrained airfoil case and the emergence of an additional pole on the positive real axis at the divergence speed for the restrained airfoil case. Inconsistencies and differences among published aeroelastic root loci and the new, exact results are discussed and resolved. The generalization of a Doublet Lattice Method computer code is described and the code is applied to the calculation of root loci for the wing model for incompressible and for subsonic flow conditions. The error introduced in the reduction of the singular integral equation underlying the unsteady lifting surface theory to a linear algebraic equation is discussed. Acknowledging this inherent error, the solutions of the algebraic equation by GAAM are termed 'exact.' The singularities of the problem are discussed and exponential series approximations used in the evaluation of the kernel function shown to introduce a dense collection of poles and zeroes on the negative real axis. Again, inconsistencies and differences among published aeroelastic root loci and the new 'exact' results are discussed and resolved. In all cases, aeroelastic flutter and divergence speeds and frequencies are in good agreement with published results. The GAAM solution procedure allows complete control over Mach number, velocity, density, and complex frequency. Thus all points on the computed root loci can be matched-point, consistent solutions without recourse to complex mode tracking logic or dataset interpolation, as in the k and p-k solution methods.

Edwards, John W.↗

Refine your search to explore more results.