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 1,459 records · Page 81
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Use of the TM tasseled cap transform for interpretation of spectral contrasts in an urban scene

Investigations are being conducted with the objective to develop automated numerical image analysis procedures. In this context, an examination is performed of physically-based multispectral data transforms as a means to incorporate a priori knowledge of land radiance properties in the analysis process. A physically-based transform of TM observations was developed. This transform extends the Landsat MSS Tasseled Cap transform reported by Kauth and Thomas (1976) to TM data observations. The present study has the aim to examine the utility of the TM Tasseled Cap transform as applied to TM data from an urban landscape. The analysis conducted is based on 512 x 512 subset of the Washington, DC November 2, 1982 TM scene, centered on Springfield, VA. It appears that the TM tasseled cap transformation provides a good means to explain land physical attributes of the Washington scene. This result provides a suggestion regarding a direction by which a priori knowledge of landscape spectral patterns may be incorporated into numerical image analysis.

Goward, S. N.↗
DOE OSTI2025

Investigation of residual stress distribution in wire-arc directed energy deposited refractory molybdenum alloy utilizing numerical thermo-mechanical analysis and neutron diffraction method

Directed energy deposition (DED), a metal additive manufacturing (AM) technique, offers higher deposition rates and energy efficiency, making it suitable for fabricating components from refractory molybdenum alloys, such as molybdenum-titanium-zirconium (TZM). However, large thermal gradients and non-equilibrium thermal cycles in DED could generate high residual stress in the component, potentially deteriorating quality and performance. Thus, this study aims to investigate residual stress generation and its distribution in wire-arc DED of TZM thin-wall, utilizing thermo-mechanical analysis and high-fidelity neutron diffraction (ND) method. Two interpass temperatures (50°C and 200°C) have been considered to investigate their impact on residual stress formation. During experiments, in-situ thermal data has been recorded using thermocouples, which have been utilized for calibrating the thermal model. Thermocouple data shows a good agreement with the simulation results, having a difference of less than 10 %. Post-deposition part deformation has been observed, which is measured using a coordinate measuring machine, showing maximum values of 0.93 mm and 0.78 mm for interpass temperatures of 50°C and 200°C, respectively. Numerical predictions of distortion deviated by less than 15% from the experimental results. ND measurement and simulation results indicate that residual stress magnitude and evolution vary across the TZM deposits, revealing microstructural anisotropy in both conditions. Notably, lower interpass temperatures resulted in higher residual stresses, confirmed by experimental and simulation data. Further, this study demonstrated that an integrated experimental and thermo-mechanical analysis can potentially reveal the temperature history, part deformation, and residual stress formation in wire-arc DED TZM alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Magnetoacoustic tomography with magnetic induction for high-resolution bioimepedance imaging through vector source reconstruction under the static field of MRI magnet

Purpose: Magnetoacoustic tomography with magnetic induction (MAT-MI) is an imaging modality to reconstruct the electrical conductivity of biological tissue based on the acoustic measurements of Lorentz force induced tissue vibration. This study presents the feasibility of the authors' new MAT-MI system and vector source imaging algorithm to perform a complete reconstruction of the conductivity distribution of real biological tissues with ultrasound spatial resolution. Methods: In the present study, using ultrasound beamformation, imaging point spread functions are designed to reconstruct the induced vector source in the object which is used to estimate the object conductivity distribution. Both numerical studies and phantom experiments are performed to demonstrate the merits of the proposed method. Also, through the numerical simulations, the full width half maximum of the imaging point spread function is calculated to estimate of the spatial resolution. The tissue phantom experiments are performed with a MAT-MI imaging system in the static field of a 9.4 T magnetic resonance imaging magnet. Results: The image reconstruction through vector beamformation in the numerical and experimental studies gives a reliable estimate of the conductivity distribution in the object with a ∼1.5 mm spatial resolution corresponding to the imaging system frequency of 500 kHz ultrasound. In addition, the experiment results suggest that MAT-MI under high static magnetic field environment is able to reconstruct images of tissue-mimicking gel phantoms and real tissue samples with reliable conductivity contrast. Conclusions: The results demonstrate that MAT-MI is able to image the electrical conductivity properties of biological tissues with better than 2 mm spatial resolution at 500 kHz, and the imaging with MAT-MI under a high static magnetic field environment is able to provide improved imaging contrast for biological tissue conductivity reconstruction.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
NASA NTRS1989

Diagnostic analysis of observed and numerically simulated cloud and precipitation structures during pre-ERICA

A preliminary 140-km simulation of a rapid cyclogenesis, using the Limited Area Mesoscale Prediction System (LAMPS) is discussed. A three-dimensional version of LAMPS used in the simulation is described, along with a trajectory model utilizing the three-hourly output history files from the LAMPS model. Preliminary integrations of the LAMPS code as well as passive microwave digital data from SSM/I aboard the DMSP satellite are used for exploring the evolving baroclinic structure and moisture field of the cyclone event. A trajectory analysis of the model histories reveals a significant role for preconditioning of the low-level upstream air, which is then drawn into the ascending warm conveyor belt ahead of the storm.

Perkey, Donald J.↗
DOE OSTI2020

Numerical modeling and performance analysis of a membrane-based air dehumidifier using ionic liquid desiccant

Membrane-based Liquid Desiccant Air Dehumidifier (MLDAD) has the potential to overcome the shortcomings of the conventional open-tower liquid desiccant air dehumidifier. The MLDAD uses membranes with high water vapor permeability to separate air stream from the liquid desiccant and thus eliminates carry-over and the resulting corrosion issues. Recently, the ionic liquid desiccant becomes a promising alternative to conventional liquid desiccant. The ionic liquid desiccant has a large potential of dehumidification, which is also non-corrosive to metals and non-crystallizable. A two-dimensional numerical heat and mass transfer model of the MLDAD using a recently identified ionic liquid desiccant, [EMIM]OAc, is presented in this paper. We report this model can simulate the performance of the MLDAD with various designs, including different selections of liquid desiccant and membrane materials, dimensions of the MLDAD, and flow patterns. For porous membranes, this model accounts for several micro-scale mass transfer mechanisms of vapor transportation across the membrane. It can also model the mass transfer performance of dehumidifiers using non-porous membranes based on the experimentally measured membrane permeance. The numerical model was validated against the performance data available from literature and experimental tests. The maximum discrepancy of the latent effectiveness between the measured and the model-predicted results was about 6% when porous membranes were used. A parametric study was conducted with the numerical model. The results indicated that the membrane permeability, air path height, and solution to air mass ratio were the most critical parameters, determining the dehumidification performance of the MLDAD.

42 ENGINEERING↗
DOE OSTI2020

Numerical modeling and performance analysis of a membrane-based air dehumidifier using ionic liquid desiccant

Membrane-based Liquid Desiccant Air Dehumidifier (MLDAD) has the potential to overcome the shortcomings of the conventional open-tower liquid desiccant air dehumidifier. The MLDAD uses membranes with high water vapor permeability to separate air stream from the liquid desiccant and thus eliminates carry-over and the resulting corrosion issues. Recently, the ionic liquid desiccant becomes a promising alternative to conventional liquid desiccant. The ionic liquid desiccant has a large potential of dehumidification, which is also non-corrosive to metals and non-crystallizable. A two-dimensional numerical heat and mass transfer model of the MLDAD using a recently identified ionic liquid desiccant, [EMIM]OAc, is presented here. This model can simulate the performance of the MLDAD with various designs, including different selections of liquid desiccant and membrane materials, dimensions of the MLDAD, and flow patterns. For porous membranes, this model accounts for several micro-scale mass transfer mechanisms of vapor transportation across the membrane. It can also model the mass transfer performance of dehumidifiers using non-porous membranes based on the experimentally measured membrane permeance. The numerical model was validated against the performance data available from literature and experimental tests. The maximum discrepancy of the latent effectiveness between the measured and the model-predicted results was about 6% when porous membranes were used. A parametric study was conducted with the numerical model. Finally, the results indicated that the membrane permeability, air path height, and solution to air mass ratio were the most critical parameters, determining the dehumidification performance of the MLDAD.

42 ENGINEERING↗
DOE OSTI2015

Non-monotonic resonance in a spatially forced Lengyel-Epstein model

We study resonant spatially periodic solutions of the Lengyel-Epstein model modified to describe the chlorine dioxide-iodine-malonic acid reaction under spatially periodic illumination. Using multiple-scale analysis and numerical simulations, we obtain the stability ranges of 2:1 resonant solutions, i.e., solutions with wavenumbers that are exactly half of the forcing wavenumber. We show that the width of resonant wavenumber response is a non-monotonic function of the forcing strength, and diminishes to zero at sufficiently strong forcing. We further show that strong forcing may result in a π/2 phase shift of the resonant solutions, and argue that the nonequilibrium Ising-Bloch front bifurcation can be reversed. We attribute these behaviors to an inherent property of forcing by periodic illumination, namely, the increase of the mean spatial illumination as the forcing amplitude is increased.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Finite spatial-grid effects in energy-conserving particle-in-cell algorithms

Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. We also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length. Finite-grid (or aliasing) instabilities are pervasive in particle-in-cell (PIC) plasma simulation algorithms, and force the modeler to resolve the smallest (Debye) length scale in the problem regardless of dynamical relevance. These instabilities originate in the aliasing of interpolation errors between mesh quantities and particles (which live in the space–time continuum). Recently, strictly energy-conserving PIC (EC-PIC) algorithms have been developed that promise enhanced robustness against aliasing instabilities. In this study, we confirm by analysis that EC-PIC is stable against aliasing instabilities for stationary plasmas. For drifting plasmas, we demonstrate by analysis and numerical experiments that, while EC-PIC algorithms are not free from these instabilities in principle, they feature a benign stability threshold for finite-temperature plasmas that make them usable in practice for a large class of problems (featuring ambipolarity and realistic ion-electron mass ratios) without the need to consider the size of the Debye length. Finally, we also demonstrate that this threshold is absent for the popular momentum-conserving PIC algorithms, which are therefore unstable for both drifting and stationary plasmas beyond a threshold in cell size compared to Debye length.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

Development of a unified numerical procedure for free vibration analysis of structures

This paper presents the details of a unified numerical algorithm and the associated computer program developed for the efficient determination of natural frequencies and modes of free vibration of structures. Both spinning and nonspinning structures with or without viscous and/or structural damping may be solved by the current procedure; in addition, the program is capable of solving static problems with multiple load cases as well as the quadratic matrix eigenproblem associated with a finite dynamic element structural discretization.

Gupta, K. K.↗
NASA NTRS1986

A numerical method for the design and analysis of counter-rotating propellers

A numerical method has been developed using the techniques of Lock and Theodorsen as described by Davidson to design and analyze counter-rotating propellers. The design method develops the optimum propeller geometry by calculating the planform and twist distribution for each propeller disk through the use of specific inputs of engine shaft horsepower, diameter, and disk spacing. The analysis method calculates the performance of a given counter-rotating propeller system at any flight condition. Using the NACA four-digit airfoil family, the performance of a counter-rotating propeller design for a given flight condition was investigated in the design and analysis mode.

Playle, S. C.↗
DOE OSTI2016

A high-throughput technique for determining grain boundary character non-destructively in microstructures with through-thickness grains

Grain boundaries (GBs) govern many properties of polycrystalline materials. However, because of their structural variability, our knowledge of GB constitutive relations is still very limited. We present a novel method to characterise the complete crystallography of individual GBs non-destructively, with high-throughput, and using commercially available tools. This method combines electron diffraction, optical reflectance and numerical image analysis to determine all five crystallographic parameters of numerous GBs in samples with through-thickness grains. We demonstrate the technique by measuring the crystallographic character of about 1,000 individual GBs in aluminum in a single run. Our method enables cost- and time-effective assembly of crystallography–property databases for thousands of individual GBs. Furthermore, such databases are essential for identifying GB constitutive relations and for predicting GB-related behaviours of polycrystalline solids.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Analytical Solutions for Seismic Fluid-Structure Interaction of Head-Supported Cylindrical Tanks

Fluid-structure interaction (FSI) of earthquake shaking is a design consideration in the civil, nuclear, chemical, and mechanical industries. Preliminary sizing of structures and components in these industries will often use analytical solutions because there is insufficient information to warrant the analysis of numerical models. Here, this study extends prior analytical work on seismic FSI of base-supported tanks in the 1980s to accommodate head-supported tanks with geometries similar to those proposed for new-build, liquid (fluid)-filled, advanced reactors. Analytical solutions are derived for flexible, head-supported, cylindrical tanks subjected to small-amplitude, unidirectional, horizontal seismic motion, and address frequencies of lateral motions of the tank, hydrodynamic pressures, frequencies of waves, wave heights, and reactions at the head support. Solutions are presented for a range of tank dimensions and could be used for a) preliminary design of head-supported cylindrical tanks and their supporting structures with idealized seismic input, and b) verification of numerical models to be used for final seismic analysis, design, and probabilistic risk assessment.

42 ENGINEERING↗
DOE OSTI2020

FPDetect: Efficient Reasoning About Stencil Programs Using Selective Direct Evaluation

We present FPDetect, a low-overhead approach for detecting logical errors and soft errors affecting stencil computations without generating false positives. We develop an offline analysis that tightly estimates the number of floating-point bits preserved across stencil applications. This estimate rigorously bounds the values expected in the data space of the computation. Violations of this bound can be attributed with certainty to errors. FPDetect helps synthesize error detectors customized for user-specified levels of accuracy and coverage. FPDetect also enables overhead reduction techniques based on deploying these detectors coarsely in space and time. Experimental evaluations demonstrate the practicality of our approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Advancing cloud lifecycle representation in numerical models using innovative analysis methods that bridge arm observations over a breadth of scales

From its location on the subtropics-midlatitude boundary, the Azores is influenced by both the subtropical high pressure and the midlatitude baroclinic storm regimes, and therefore experiences a wide range of cloud structures, from fair-weather scenes to stratocumulus sheets to deep convective systems. This project combined three types of data sets to study cloud variability in the Azores: a satellite analysis of cloud regimes, a reanalysis characterization of storminess, and a 19-month field campaign that occurred on Graciosa Island. Combined analysis of the three data sets provides a detailed picture of cloud variability and the respective dynamic influences, with emphasis on low clouds that constitute a major uncertainty source in climate model simulations. The satellite cloud regime analysis shows that the Azores cloud distribution is similar to the mean global distribution and can therefore be used to evaluate cloud simulation in global models. Regime analysis of low clouds shows that stratocumulus decks occur under the influence of the Azores high-pressure system, while shallow cumulus clouds are sustained by cold-air outbreaks, as revealed by their preference for post-frontal environments and northwesterly flows. An evaluation of CMIP5 climate model cloud regimes over the Azores shows that all models severely underpredict shallow cumulus clouds, while most models also underpredict the occurrence of stratocumulus cloud decks. It is demonstrated that carefully selected case studies can be related through regime analysis to climatological cloud distributions, and a methodology is suggested utilizing process-resolving model simulations of individual cases to better understand cloud-dynamics interactions and attempt to explain and correct climate model cloud deficiencies.

58 GEOSCIENCES↗
DOE OSTI2018

Numerical modeling and sensitivity analysis of seawater intrusion in a dual-permeability coastal karst aquifer with conduit networks

Long-distance seawater intrusion has been widely observed through the subsurface conduit system in coastal karst aquifers as a source of groundwater contaminant. In this study, seawater intrusion in a dual-permeability karst aquifer with conduit networks is studied by the two-dimensional density-dependent flow and transport SEAWAT model. Local and global sensitivity analyses are used to evaluate the impacts of boundary conditions and hydrological characteristics on modeling seawater intrusion in a karst aquifer, including hydraulic conductivity, effective porosity, specific storage, and dispersivity of the conduit network and of the porous medium. The local sensitivity analysis evaluates the parameters' sensitivities for modeling seawater intrusion, specifically in the Woodville Karst Plain (WKP). A more comprehensive interpretation of parameter sensitivities, including the nonlinear relationship between simulations and parameters, and/or parameter interactions, is addressed in the global sensitivity analysis. The conduit parameters and boundary conditions are important to the simulations in the porous medium because of the dynamical exchanges between the two systems. The sensitivity study indicates that salinity and head simulations in the karst features, such as the conduit system and submarine springs, are critical for understanding seawater intrusion in a coastal karst aquifer. The evaluation of hydraulic conductivity sensitivity in the continuum SEAWAT model may be biased since the conduit flow velocity is not accurately calculated by Darcy's equation as a function of head difference and hydraulic conductivity. In addition, dispersivity is no longer an important parameter in an advection-dominated karst aquifer with a conduit system, compared to the sensitivity results in a porous medium aquifer. In the end, the extents of seawater intrusion are quantitatively evaluated and measured under different scenarios with the variabilities of important parameters identified from sensitivity results, including salinity at the submarine spring with rainfall recharge, sea level rise, and a longer simulation time under an extended low rainfall period.

58 GEOSCIENCES↗
NASA NTRS1976

Inlet spillage drag tests and numerical flow-field analysis at subsonic and transonic speeds of a 1/8-scale, two-dimensional, external-compression, variable-geometry, supersonic inlet configuration

Accurate spillage drag and pressure data are presented for a realistic supersonic inlet configuration. Results are compared with predictions from a finite-differencing, inviscid analysis computer procedure. The analytical technique shows good promise for the evaluation of inlet drag, but necessary refinements were identified. A detailed description of the analytical procedure is contained in the Appendix.

Hawkins, J. E.↗