Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical analysis”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Hurricane Debby - An illustration of the complementary nature of VAS soundings and cloud and water vapor motion winds

The utility of VISSR Atmospheric Sounder (VAS) temperature and moisture soundings and cloud and water vapor motion winds in defining a storm and its surroundings at subsynoptic scales has been examined using a numerical analysis and prognosis system. It is shown that the VAS temperature and moisture data, which specify temperature and moisture well in cloud-free areas, are complemented by cloud and water vapor motion data generated in the cloudy areas. The cloud and water vapor 'winds' provide thermal gradient information for interpolating the soundings across cloudy regions. The loss of analysis integrity due to the reduction of VAS sounding density in the cloudy regions associated with synoptic activity is ameliorated by using cloud and water vapor motion winds. The improvement in numerical forecasts resulting from the addition of these data to the numerical analysis is recorded.

Le Marshall, J. F.↗
DOE OSTI2023

Conjugate Heat Transfer Analysis in Pressurized Helium Gas Cooling Channels

Los Alamos National Laboratory (LANL) is working with NorthStar Medical Radioisotopes (NMR) on their efforts to produce Mo-99 from the irradiation of Mo-100 targets. The NMR target consists of an Inconel window that allows the electron beam to penetrate and irradiate a stack of Mo-100 disks. The irradiation process generates large amounts of heat and is cooled using pressurized helium gas flowing at 400 g/s. LANL provides both numerical analysis support and experimental support. Computational fluid dynamics (CFD) and conjugate heat transfer (CHT) analysis are performed on the NMR target to analyze the helium flow and volumetric heating through the NMR system. In this iteration of numerical analysis, we update the Mo-100 disk thickness to 0.74 mm and the gap between the Inconel window and the first disk to 0.69 mm. Estimates of the heat flux from the electron beams from Monte Carlo N-Particle code (MCNP) analysis were combined with CFD to perform the CHT analysis. The pressure drop through the targets, force applied on the first disk, beam window and target disk temperatures are some of the quantities estimated using the simulations. The pressure drop through the target was estimated as 0.198 MPa, while the force due to venturi effects across the first disk was estimated to be 100 N. The results show maximum window temperatures over 600 °C and peak disk temperatures over 1600 °C. Past work has shown that beyond 650-700 °C, the Inconel window becomes susceptible to yield and rupture. The increased temperature in these results reduces the margin of error possible on beam spot size. The disk temperatures showed a sinusoidal distribution with the 6 th – 8 th disk from each window exhibiting highest disk temperatures. The peak temperature of the stainless steel housing was close to 316 °C and that of the laminations and spacers was approximately 582 °C.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

Design, analysis and numerical experiments for the virtual element p and hp approximations of elliptic eigenvalue problems

We discuss the p- and the hp-versions of the virtual element method for the approximation of the eigenvalues of elliptic operators with a potential term on polygonal meshes. An application of this model is provided by the Schr¨odinger equation with a pseudo-potential term. We present in details the analysis of the p-version of the method, proving the exponential convergence in the case of analytic eigenfunctions. The theoretical results are supplied with a wide set of experiments. We also show numerically that, in the case of eigenfunctions with finite Sobolev regularity, an exponential approximation of the eigenvalues in terms of the cubic root of the number of degrees of freedom can be obtained by employing hp-refinements. Importantly, the geometric flexibility of polygonal meshes is exploited in the construction of the hp-spaces.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Final Scientific/Technical Report, Multi-Sensor Inspection and Robotic Systems for Dry Storage Casks

Once removed from service, fuel rod assemblies are cooled in a pool and then placed in dry storage until they can be transported to a repository for disposal. The extended duration that used nuclear fuel spends in dry storage makes it prudent to be able to nondestructively inspect the dry storage casks. The degradation mechanisms of most concern are chloride-induced stress corrosion cracking of the stainless steel canister that confines the spent fuel and distributed cracking of the concrete in the overpack structure. Multi-sensor inspection and robotic systems were researched, developed, built, and tested. One system operates completely outside the cask to inspect the concrete. The other system, known as the PRINSE system, for Proactive Robotic Inspection of Nuclear Storage Enclosures, enters the cask through an outlet vent and navigates around the outer surface of the stainless steel canister. The PRINSE system delivers sensors for: (i) real-time characterization of the canister surface in order to assess whether chlorides are present, (ii) detecting the presence of cracks using shear-horizontal ultrasonic guided waves, (iii) measuring the temperature of the canister surface, and (iv) a gamma radiation counter. The four sensing modalities are applicable to any cask containing a stainless steel canister, however the robotic delivery system is designed specifically for Holtec’s HI-STORM family of casks. Horizontal axis casks have exposed concrete surfaces. A fully contactless air-coupled scanning ultrasonic test platform and a frequency-wavenumber domain data analysis method are used for these concrete structures. Analysis, numerical simulation and experimental results demonstrate that distributed cracks in concrete can be detected and cracked regions localized within a scanned region of concrete. Vertical axis casks (like the HI-STORM family) have steel-clad concrete overpacks. In a separate effort, a semi-coupled test platform and a guided wave attenuation-based data analysis method are used for disbond detection. Analysis, numerical simulation and experimental results show that the bonding conditions of steel-concrete interfaces can be quantitatively characterized using the guided wave signal analysis scheme. The PRINSE system is designed for vertical axis dry storage casks, specifically the HI-STORM family. The double pulse laser-induced breakdown spectroscopy (LIBS) system takes advantage of an optical fiber to deliver intense light to the canister surface that enables surface composition characterization in-situ in real time. Chlorine concentrations as low as 10 mg/m2 were detected on stainless steel. Additionally, a simpler laser setup was created to detect sodium as a surrogate for chlorine, and sodium was detectable when the chlorine concentration was as low as 0.5 mg/m2. The electromagnetic acoustic transducers (EMATs) enable detection of stress corrosion cracks. The EMATs generate shear-horizontal waves that propagate around the canister and are reflected by cracks enabling 100% of the canister weld lines to be inspected. The EMAT-based crack detection system was tested at the Pacific Northwest National Laboratory on a canister containing blind defects and was shown to detect the crack sizes it was designed to detect. The PRINSE system is a novel marsupial robot. The command center and winch outside the cask control the delivery arm and train of cars that traverse the outside of the canister from the cask ventilation system. The sensors are housed in the cargo bays of three cars that comprise the train, which is lowered to the bottom of the canister by the winch. Encoders on the winch enable the sensors to be delivered to a prescribed location, e.g., for re-inspection of a problem area. The PRINSE system does not mar the canister surface, can be extracted manually in the case of a power failure, and is designed for an elevated temperature and gamma radiation environment, although not all components of the prototype are made from the final materials. The PRINSE system was tested on a full canister that was placed inside a HI-STORM cask at the Holtec Manufacturing Division. It was also demonstrated on a full-height mockup at Penn State to a live audience.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1992

Numerical sensitivity analysis of passive EHF and SMMW channels to tropospheric water vapor, clouds, and precipitation

Potential uses of specific extremely High Frequency (EHF) and Sub-Millimeter-Wave (SMMW) channels at 90, 166, 183, 220, 325, 340, and 410 GHz for passive spaceborne remote sensing of the troposphere and lower stratosphere are investigated using an iterative numerical radiative transfer model. Collectively, these channels offer potential for high spatial resolution imaging using diffraction-limited apertures of practical size, along with the ability to profile water vapor, map precipitation beneath optically opaque cloud cover, and to measure nonprecipitating cloud (e.g., cirrus) parameters. A widely-spaced set of EHF and SMMW channels can yield observable degrees of freedom related to clouds and precipitation not available by exclusively using the more thoroughly studied microwave channels below 183 GHz. A new passive airborne imaging instrument for tropospheric meteorological sensing is described.

Gasiewski, A. J.↗
NASA NTRS1999

Gas Evolution Dynamics in Godunov-Type Schemes and Analysis of Numerical Shock Instability

In this paper we are going to study the gas evolution dynamics of the exact and approximate Riemann solvers, e.g., the Flux Vector Splitting (FVS) and the Flux Difference Splitting (FDS) schemes. Since the FVS scheme and the Kinetic Flux Vector Splitting (KFVS) scheme have the same physical mechanism and similar flux function, based on the analysis of the discretized KFVS scheme the weakness and advantage of the FVS scheme are closely observed. The subtle dissipative mechanism of the Godunov method in the 2D case is also analyzed, and the physical reason for shock instability, i.e., carbuncle phenomena and odd-even decoupling, is presented.

Xu, Kun↗
DOE OSTI2015

Numerical—experimental analysis of two floor platform designs for auxiliary construction truck

In an auxiliary construction truck, the permanent deformations created by impact of rocks against the floor platform during the daily work generate the premature failure and need for replacement of the platform floor. This article shows the process followed in the analysis of two floor platforms for auxiliary construction trucks. This analysis has been carried out by means of the methodology developed in collaboration with truck manufacturer (ZAMARBU S.L) and the research group “New Technologies Applied on Vehicles and Road Safety” (VEHIVIAL) of the University of Zaragoza. The methodology applied in the floor platform development has been based on the application of numerical techniques by means of the finite element method (FEM) with explicit integration of dynamic equilibrium equation and the validation of the numerical results by means of experimental test. Applied load cases correspond to the impact of a sphere of 50, 150, and 250 kg. These load cases reproduce the impact caused by stones falling on the floor platform, which are the most restrictive for the floor design, because the vehicle’s rigidity and resistance against torsion and bending are provided for the chassis. The results of the analysis served to demonstrate the efficiency of the methodology developed for the optimization of the floor platform for auxiliary construction truck, getting an innovative floor platform design that improves the current models.

42 ENGINEERING↗
DOE OSTI2014

Numerical stability analysis of the pseudo-spectral analytical time-domain PIC algorithm

The pseudo-spectral analytical time-domain (PSATD) particle-in-cell (PIC) algorithm solves the vacuum Maxwell's equations exactly, has no Courant time-step limit (as conventionally defined), and offers substantial flexibility in plasma and particle beam simulations. It is, however, not free of the usual numerical instabilities, including the numerical Cherenkov instability, when applied to relativistic beam simulations. This paper derives and solves the numerical dispersion relation for the PSATD algorithm and compares the results with corresponding behavior of the more conventional pseudo-spectral time-domain (PSTD) and finite difference time-domain (FDTD) algorithms. In general, PSATD offers superior stability properties over a reasonable range of time steps. More importantly, one version of the PSATD algorithm, when combined with digital filtering, is almost completely free of the numerical Cherenkov instability for time steps (scaled to the speed of light) comparable to or smaller than the axial cell size.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

A new technique for aerodynamic noise calculation

A novel method for the numerical analysis of aerodynamic noise generation is presented. The method involves first solving for the time-dependent incompressible flow for the given geometry. This fully nonlinear method that is tailored to extract the relevant acoustic fluctuations seems to be an efficient approach to the numerical analysis of aerodynamic noise generation.

Hardin, J. C.↗
DOE OSTI2023

Fracture of PC and PMMA cylinders in the Taylor tests

The present paper analyzes failure observed in polymer rods collided with a steel anvil under the Taylor cylinder impact experiments. The experimental and numerical analysis demonstrates that association of the failure patterns with the loading mode depends on the rod material. Two materials are considered in the paper, specifically, polymethylmethacrylate (PMMA) and polycarbonate (PC). Damage in PMMA rods at the Taylor tests is shown to be dominated by the shear load resulting in plastic deformation. In contrast, the major mechanism of damage in PC rods is governed by omni-directional (pressure) mode of load. Analysis of the damage in sections of the rods associated with the glass transition is also addressed in the present theoretical study. Failure in these polymers follows the traditional fracture wave mechanism for brittle glass-type materials due to load in the corresponding mode. Numerical analysis employs a two-phase model with improved inter-phase thermal balance accounting for the inter-phase mass, momentum, and heat exchanges during the glass phase transition, while damage being accumulated in time. In conclusion, the simulation results demonstrate a good correlation with the experimental results.

36 MATERIALS SCIENCE↗
NASA NTRS1988

Nonequilibrium flow computations. 1: An analysis of numerical formulations of conservation laws

Modern numerical techniques employing properties of flux Jacobian matrices are extended to general, nonequilibrium flows. Generalizations of the Beam-Warming scheme, Steger-Warming and van Leer Flux-vector splittings, and Roe's approximate Riemann solver are presented for 3-D, time-varying grids. The analysis is based on a thermodynamic model that includes the most general thermal and chemical nonequilibrium flow of an arbitrary gas. Various special cases are also discussed.

Liu, Yen↗
NASA NTRS1989

Nonequilibrium flow computations. I - An analysis of numerical formulations of conservation laws

Modern numerical techniques employing properties of flux Jacobian matrices are extended to general, nonequilibrium flows. Generalizations of the Beam-Warming scheme, Steger-Warming and van Leer Flux-vector splittings, and Roe's approximate Riemann solver are presented for 3-D, time-varying grids. The analysis is based on a thermodynamic model that includes the most general thermal and chemical nonequilibrium flow of an arbitrary gas. Various special cases are also discussed.

Liu, Yen↗
NASA NTRS2005

The Effects of Instrumental Elliptical Polarization on Stellar Point Spread Function Fine Structure

We present procedures and preliminary results from a study on the effects of instrumental polarization on the fine structure of the stellar point spread function (PSF). These effects are important to understand because the the aberration caused by instrumental polarization on an otherwise diffraction-limited will likely have have severe consequences for extreme high contrast imaging systems such as NASA's planned Terrestrial Planet Finder (TPF) mission and the proposed NASA Eclipse mission. The report here, describing our efforts to examine these effects, includes two parts: 1) a numerical analysis of the effect of metallic reflection, with some polarization-specific retardation, on a spherical wavefront; 2) an experimental approach for observing this effect, along with some preliminary laboratory results. While the experimental phase of this study requires more fine-tuning to produce meaningful results, the numerical analysis indicates that the inclusion of polarization-specific phase effects (retardation) results in a point spread function (PSF) aberration more severe than the amplitude (reflectivity) effects previously recorded in the literature.

high angular resolution↗
DOE OSTI2017

Deep Vadose Zone Contaminant Flux Evaluation at the Hanford BY-Cribs Site using Forward and Imposed Concentration Modeling Approaches

For sites with vadose zone contamination, the nature and extent of groundwater contaminant plumes are related to the flux into groundwater. In thick vadose zones at arid sites, transport is likely relative slow, complicating direct measurement of contaminant flux. A predictive numerical analysis was applied to estimate contaminant flux to groundwater at the Hanford BY Cribs site. The site received large amounts of waste water in the 1950’s, containing contaminants such as technicium-99 and nitrate. Recently, large increases in groundwater concentrations have been observed in site groundwater wells and, for potential remediation purposes, the assumption that the changes are the result of increased vadose zone fluxes into groundwater needed to be evaluated. The analysis approach uses forward simulations, tracking contaminants during and after disposal, and imposed concentration (IC) simulations, where concentrations, derived from electrical conductivity data obtained in 2007, are imposed on the flow field. A total of three IC methods, relating electrical conductivity to nitrate concentrations, are developed and tested. The numerical analysis shows that the forward simulations predict high fluxes during and directly after disposal, with decreasing trends starting not later than 1970. Although the IC simulations show differences in the nitrate vadose zone mass imposed in 2007 compared to forward simulations, the IC simulations also show declining nitrate fluxes over time. The predicted fluxes of the forward and IC simulations after 2007 are in reasonable agreement. Given the declining nature of the contaminant fluxes, it can be concluded that the observed concentration increases in the groundwater below the BY Cribs site cannot be explained as resulting from an increase in contaminant flux from the vadose zone. An alternative explanation, based on the observed decrease in aquifer thickness below the site, is presented.

Oostrom, Martinus↗
NASA NTRS1991

Three-dimensional flows in a transonic compressor rotor

This study involves an experimental and numerical investigation of the three-dimensional flows in a transonic compressor rotor. A variety of data which could be used, in a complementary fashion, to validate/calibrate the computational fluid dynamics turbomachinery code and improve understanding of the flow physics, were acquired. Detailed radial survey data which consisted of total pressure, total temperature, static pressure and flow angle were obtained at stations upstream and downstream of the rotor blade. Detailed velocity and turbulence profiles were obtained upstream of the rotor and used as the upstream boundary conditions for the numerical analysis. Calibrated flush-mounted hot film probes were used to measure wall shear stress on the hub and casing walls upstream of the rotor. The blade-to-blade shear-stress angle distributions were obtained at two axial locations on the rotor casing, using flush-mounted hot film probes. A numerical analysis conducted using a three-dimensional Navier-Stokes code was compared with the experimental results.

Reid, Lonnie↗
DOE OSTI2023

Learning continuous models for continuous physics

Abstract Dynamical systems that evolve continuously over time are ubiquitous throughout science and engineering. Machine learning (ML) provides data-driven approaches to model and predict the dynamics of such systems. A core issue with this approach is that ML models are typically trained on discrete data, using ML methodologies that are not aware of underlying continuity properties. This results in models that often do not capture any underlying continuous dynamics—either of the system of interest, or indeed of any related system. To address this challenge, we develop a convergence test based on numerical analysis theory. Our test verifies whether a model has learned a function that accurately approximates an underlying continuous dynamics. Models that fail this test fail to capture relevant dynamics, rendering them of limited utility for many scientific prediction tasks; while models that pass this test enable both better interpolation and better extrapolation in multiple ways. Our results illustrate how principled numerical analysis methods can be coupled with existing ML training/testing methodologies to validate models for science and engineering applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Sensitivity analysis of numerical modeling input parameters on floating offshore wind turbine loads

Abstract. Floating wind turbines must withstand a unique and challenging set of loads from the wind and ocean environment. To de-risk development, accurate predictions of these loads are necessary. Uncertainty in modeling predictions leads to larger required safety factors, increasing project costs and the levelized cost of energy. Complex aero-hydro-elastic modeling tools use many input parameters to represent the wind, waves, current, aerodynamic loads, hydrodynamic loads, and structural properties. It is helpful to understand which of these parameters ultimately drives a design. In this work, an ultimate and fatigue-proxy load sensitivity analysis was performed with 35 different input parameters, using an elementary effects approach to identify the most influential parameters for a case study involving the National Renewable Energy Laboratory (NREL) 5 MW baseline wind turbine atop the OC4-DeepCwind semisubmersible during normal operation. The importance of each parameter was evaluated using 14 response quantities of interest across three operational wind speed conditions. The study concludes that turbulent wind velocity standard deviation is the parameter with the strongest sensitivity; this value is important not just for turbine loads, but also for the global system response. The system center of mass in the wind direction is found to have the highest impact on the system rotation and tower loads. The current velocity is found to be the most dominating parameter for the system global motion and consequently the mooring loads. All tested wind turbulence parameters in addition to the standard deviation are also found to be influential. Wave characteristics are influential for some fatigue-proxy loading but do not significantly impact the extreme ultimate loads in these operational load cases. The required number of random seeds for stochastic environmental conditions is considered to ensure that the sensitivities are due to the input parameters and not due to the seed. The required number of analysis points in the parameter space is identified so that the conclusions represent a global sensitivity. The results are specific to the platform, turbine, and choice of parameter ranges, but the demonstrated approach can be applied widely to guide focus in parameter uncertainty.

17 WIND ENERGY↗