Search NASA⌕ Search

SEARCH · Search NASA

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Numerical analysis of HgCdTe dual-band infrared detector

In this paper we describe recent progress in the HgCdTe multilayer heterostructures grown by metaloorganic chemical vapour deposition on GaAs substrates for photodetectors operated above 200 K. We present example of design and characterization of the new classical N{sup +}–n–P{sup +}–p–N{sup +} back to back HgCdTe dual-band photodiode structure, that operates within the mid-wavelength infrared band in sequential mode. A numerical modelling was used for investigation of the device design on the current responsivity and dark currents. The program based on the solution of the system of the carrier transport equations including the whole spectrum of various generation and recombination mechanisms consisting of Shockley–Read–Hall, Auger and radiative generation-recombination terms. Additionally, tunneling effects such as band-to-band and trap-assisted tunnelling models are included in the continuity equations by incorporating them as additional generation–recombination processes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1971

Numerical analysis of shells of revolution.

Summary of some of the major computer programs in existence for the analysis of shells of revolution based on numerical integration and finite-difference procedures. Selected results obtained with several of these programs are presented. A number of significant points is made on the basis of these results.

Fulton, R. E.↗
DOE OSTI2016

Numerical analysis of electrostatic ion cyclotron instability in an electron-positron-ion plasma

This paper presents a theoretical study of the effects of positron density on the electrostatic ion cyclotron instability in an electron-positron-ion plasma using the kinetic theory approach. It is supposed that positrons and electrons can drift parallel to the magnetic field either in the same or the opposite directions. The dispersion relation for the electrostatic ion cyclotron waves in an electron-positron-ion plasma is derived, and the numerical results are investigated. It is found that an increase in positron concentration increases the critical drift velocity for excitation of the instability in both configurations. It is also found that as the positron concentration increases the growth rate of instability decreases. In addition, it is shown that at low velocities the maximum instability growth rate for the unidirectional case is higher than the counter-streaming case; however, after a certain velocity, the maximum growth rate in the counter-streaming case dominates that of the unidirectional case.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Numerical analysis of a thermal contact problem with adhesion

In this paper, we study a class of dynamic thermal problems involving a frictional normal compliance adhesive contact model and a non-clamped condition for visco-elastic materials. The variational formulation of the problem leads to a general system defined by a second-order quasi-variational evolution inequality on the displacement field coupled with the two first-order evolution equations describing the evolution of the temperature and adhesion. We establish an existence and uniqueness result of the solution on displacement, temperature and adhesion. Then we provide a fully discrete numerical scheme of approximation and derive an error estimate. Finally, various numerical computations are included.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Numerical Analysis of Blast Effects and Mitigation in the Far-Field from Small Explosions

Requirements for explosive safety are often given in terms of a “K-Factor”, correlating incident blast effects with the distance and TNT equivalent weight of an explosion. Traditionally, this is conducted using empirical correlations to experimental measurements (e.g., the Kingery–Bulmash equations). In the far-field, empirical verification of incident overpressure and impulse magnitudes can be difficult; extrapolations from data give expected values at reasonable standoff distances that sometimes are too small to measure on available equipment but are larger than some regulations require. The present paper describes the results of numerical hydrocode analysis to verify the expected incident overpressure and impulse from small hemispherical ground charges of TNT at these relatively large distances. Furthermore, the dynamic effect of incident blast waves on lightweight, modular mitigation barriers is studied to gauge their effectiveness at providing safety standard compliance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Shape effect on drying behavior of cement-based materials: Mechanisms and numerical analysis

The shape of the specimen or structure has a direct impact on the behavior of drying concrete. This effect is not taken into account in most of shrinkage models and codes. Only the current B4 model and its previous version B3 model incorporate a shape factor. The determination of this coefficient is based on the nonlinear moisture diffusion analysis, which can explain the shape effect on the shrinkage rate. However, the model integrates this coefficient in the ultimate shrinkage values as well. An experimental study has been performed to confirm and understand this behavior. The mass loss and drying shrinkage of cylinders showed higher rate and ultimate values than the prisms, for the same effective cross section thickness. To explain this result, the drying, the differential shrinkage, the internal stress and the damaged depth were determined for both shapes through numerical modeling. The main input data of the model were the measured mass loss and shrinkage. The calculated desorption isotherm curve of the cylinder was below that of the prism, which is consistent with its higher ultimate mass-loss. The determination of stress profiles indicates slower stress mitigation for the cylinder. The cylinder develops a larger relative damaged area, defined as the part of the specimen where cracking is likely to occur due to tensile stresses. The results actually show that differential shrinkage between the edge and the core is more pronounced for the cylinder and decreases with time.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Numerical Analysis of Flow Instability in Once-through Steam Generator

Two-phase flow instability is a critical area to consider in the design of steam generators. Oscillation of the flow rate and the system pressure caused by heating may cause several problems and is unfavorable in an extreme case because it causes a burnout of the heat transfer surface. A large amplitude variation in the heat transfer tube wall temperature generally occurs under density wave oscillation (DWO) conditions, and a continuous cycling tube wall temperature may lead to damage to the tube, causing thermal fatigue problems. Therefore, instability of the flow should be avoided, and a two-phase flow system should be designed to have a sufficient margin to avoid this. A computer program was developed for predicting a two-phase flow instability in a steam generator heated by liquid sodium, and the results of numerical studies for several parameters are presented in this paper.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1994

Numerical analysis of tangential slot blowing on a generic chined forebody

A numerical study is performed to investigate the effects of tangential slot blowing on a generic chined forebody. The Reynolds-averaged, thin-layer, Navier-Stokes equations are solved to obtain the high-angle-of-attack viscous flow field about a generic chined forebody. Tangential slot blowing is investigated as a means of forebody flow control to generate side force and yawing moment on the forebody. The effects of jet mass flow ratios, angle of attack, and blowing slot location in the axial and circumferential directions are studied. The computed results are compared with available wind tunnel experimental data. The solutions with and without blowing are also analyzed using helicity density contours, surface flow patterns, and off-surface instantaneous streamlines. The results of this analysis provide details of the flow field about the generic chined forebody, as well as show that tangential slot blowing can be used as a means of forebody flow control to generate side force and yawing moment.

Agosta, Roxana M.↗
NASA NTRS2018

Nosetip Bluntness Effects on Transition at Hypersonic Speeds: Experimental and Numerical Analysis Under NATO STO AVT-240

The existing database of transition measurements in hypersonic ground facilities has established that the onset of boundary layer transition over a circular cone at zero angle of attack shifts downstream as the nosetip bluntness is increased with respect to a sharp cone. However, this trend is reversed at sufficiently large values of the nosetip Reynolds number, so that the transition onset location eventually moves upstream with a further increase in nosetip bluntness. This transition reversal phenomenon, which cannot be explained on the basis of linear stability theory, was the focus of a collaborative investigation under the NATO STO group AVT-240 on Hypersonic Boundary-Layer Transition Prediction. The current paper provides an overview of that effort, which included wind tunnel measurements in three different facilities and theoretical analysis related to modal and nonmodal amplification of boundary layer disturbances. Because neither first and second-mode waves nor entropy-layer instabilities are found to be substantially amplified to initiate transition at large bluntness values, transient (i.e., nonmodal) disturbance growth has been investigated as the potential basis for a physics-based model for the transition reversal phenomenon. Results of the transient growth analysis indicate that disturbances that are initiated within the nosetip or in the vicinity of the juncture between the nosetip and the frustum can undergo relatively significant nonmodal amplification and that the maximum energy gain increases nonlinearly with the nose radius of the cone. This finding does not provide a definitive link between transient growth and the onset of transition, but it is qualitatively consistent with the experimental observations that frustum transition during the reversal regime was highly sensitive to wall roughness, and furthermore, was dominated by disturbances that originated near the nosetip.

Numerical analysis↗
NASA NTRS1992

Numerical analysis of a thermal deicer

An algorithm has been developed to numerically model the concurrent phenomena of two-dimensional transient heat transfer, ice accretion and ice shedding which arise from the use of an electrothermal pad. The Alternating Direction Implicit method is used to simultaneously solve the heat transfer and accretion equations occurring in a multilayered body covered with ice. In order to model the phase change between ice and water, a technique was used which assumes a phase for each node. This allows the equations to be linearized such that a direct solution is possible. This technique requires an iterative procedure to find the correct phase at each node. The computer program developed to find this solution has been integrated with the NASA/Lewis flow/trajectory code LEWICE.

Wright, W. B.↗
NASA NTRS1992

Numerical analysis of hypersonic turbulent film cooling flows

As a building block, numerical capabilities for predicting heat flux and turbulent flowfields of hypersonic vehicles require extensive model validations. Computational procedures for calculating turbulent flows and heat fluxes for supersonic film cooling with parallel slot injections are described in this study. Two injectant mass flow rates with matched and unmatched pressure conditions using the database of Holden et al. (1990) are considered. To avoid uncertainties associated with the boundary conditions in testing turbulence models, detailed three-dimensional flowfields of the injection nozzle were calculated. Two computational fluid dynamics codes, GASP and FDNS, with the algebraic Baldwin-Lomax and k-epsilon models with compressibility corrections were used. It was found that the B-L model which resolves near-wall viscous sublayer is very sensitive to the inlet boundary conditions at the nozzle exit face. The k-epsilon models with improved wall functions are less sensitive to the inlet boundary conditions. The testings show that compressibility corrections are necessary for the k-epsilon model to realistically predict the heat fluxes of the hypersonic film cooling problems.

Chen, Y. S.↗
DOE OSTI2018

Numerical Analysis of Electron Runaway in the Presence of Enhanced Field in the Vicinity of a Microtip

Transition of field-emission electrons to the runaway regime in the region of enhanced electric field determined by the configuration of a microtip on a cathode is studied at several pressures of gas medium. The problem is solved using simulation of electron motion in the presence of nonuniform electric field with the aid of the Monte Carlo procedure in the 2D configuration. Nitrogen is used as a working gas. Passage through a relatively small region of the enhanced field in the vicinity of the microtip may substantially facilitate electron escape to the runaway regime, especially, at pressures of greater than 10 atm. In our opinion, the resulting runaway electrons may provide preionization of gas medium and formation of the initial stage of a 3D discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Experimental and Numerical Analysis of Liquid Metal Embrittlement Crack Location

Advanced high-strength steels used in automotive structural components are commonly protected using zinc coatings. However, the steel/zinc system creates the potential for liquid metal embrittlement during welding. Although liquid metal embrittlement cracks are known to form, the current literature does not include crack location when assessing crack severity. In this work, TRIP1100 joints showed that LME cracks decreased strength from 7.8 to 42.2%, depending on location, between the coated (cracked) and uncoated (non-cracked) condition. Liquid metal embrittlement cracks in critical locations were observed to propagate until fracture from lap shear testing. However, cracks in non-critically loaded areas were not part of the fracture path and did not result in a significant loss in strength. This shows LME crack location can be controlled to improve joint performance and vehicle safety. In addition, a model of lap shear testing in a cracked sample showed how the presence of a crack can affect the internal stress field depending on its location.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Effects of anisotropic interaction-induced properties of hydrogen-rare gas compounds on rototranslational Raman scattering spectra: Comprehensive theoretical and numerical analysis

A comprehensive study is presented of many aspects of the depolarized anisotropic collision induced (CI) component of light scattered by weakly bound compounds composed of a dihydrogen molecule and a rare gas (Rg) atom, H{sub 2}–Rg. The work continues a series of earlier projects marking the revival of interest in linear light scattering following the development of new highly advanced tools of quantum chemistry and other theoretical, computational, and experimental means of spectral analyses. Sophisticated ab initio computing procedures are applied in order to obtain the anisotropic polarizability component’s dependence on the H{sub 2}–Rg geometry. These data are then used to evaluate the CI spectral lines for all types of Rg atoms ranging from He to Xe (Rn excluded). Evolution of the properties of CI spectra with growing polarizability/masses of the complexes studied is observed. Special attention is given to the heaviest, Kr and Xe based, scatterers. The influence of specific factors shaping the spectral lines (e.g., bound and metastable contribution, potential anisotropy) is discussed. Also the share of pressure broadened allowed rotational transitions in the overall spectral profile is taken into account and the extent to which it is separable from the pure CI contribution is discussed. We finish with a brief comparison between the obtained results and available experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Design of novel dual-port tapered waveguide plasma apparatus by numerical analysis

Microwave plasma apparatus is often of particular interest due to their superiority of low cost, electrode contamination free, and suitability for industrial production. However, there exist problems of unstable plasma and low electron density in conventional waveguide apparatus based on single port, due to low strength and non-uniformity of microwave field. This study proposes a novel dual-port tapered waveguide plasma apparatus based on power-combining technique, to improve the strength and uniformity of microwave field for the applications of plasma. A 3D model of microwave-induced plasma (field frequency 2.45 GHz) in argon at atmospheric pressure is presented. On the condition that the total input power is 500 W, simulations indicate that coherent power-combining will maximize the electric-field strength to 3.32 × 10{sup 5 }V/m and improve the uniformity of distributed microwave field, which raised 36.7% and 47.2%, respectively, compared to conventional waveguide apparatus of single port. To study the optimum conditions for industrial application, a 2D argon fluid model based on above structure is presented. It demonstrates that relatively uniform and high-density plasma is obtained at an argon flow rate of 200 ml/min. The contrastive result of electric-field distribution, electron density, and gas temperature is also valid and clearly proves the superiority of coherent power-combining to conventional technique in flow field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Numerical Analysis of Molten Pool Behavior and Spatter Formation with Evaporation During Selective Laser Melting of 316L Stainless Steel

During the selective laser-melting process, material evaporation and resultant spatter are common phenomena that bring about many defects. However, the underlying physical phenomena such as molten pool behavior and spatter formation, as evaporation occurred, are sparsely understood and difficult to observe during the process. Thus, a three–dimensional powder-scale model was established to investigate the thermal and flow behavior of the molten pool, the morphology evolution of the molten pool and keyhole, and the spatter formation with evaporation in the selective laser-melting processing of 316L stainless steel. Phase transitions and variations in interfacial force were taken into account in this model. The modified phase-field method was applied to trace the melt–gas interface. The results showed that keyhole formed in molten pool under recoil pressure, and that there were some differences in the temperature distribution and flow behavior inside and outside the keyhole. In addition, the dimension and surface morphology of the molten pool and keyhole depth altered and gradually stabilized during the process. Moreover, droplet spatter formation at rear of molten pool proceeded roughly as follows: the backward ejection of melt vapor would induce a depression in molten pool and a bump on the molten pool surface, the bump moved backward and subsequently a liquid column formed, then the liquid column formed a neck and gradually pinched off, resulting in the droplet spatter. Furthermore, the characteristic morphology of the scan track and dimension of the molten pool were obtained through experiments, which showed good agreement with the simulation results.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Numerical analysis of temperature and stress fields in hybrid indium antimonide arrays detector with laser irradiation

Thermal-stress effects are the major cause of failure in infrared focal-plane arrays detector, during laser irradiation. Based on the established three-dimensional structural model of hybrid indium antimonide infrared focal-plane arrays, the temperature and stress fields of hybrid indium antimonide detectors irradiated by 1.064-μm Gauss pulsed laser are studied, considering the temperature-dependent material parameters. The results indicate that the temperature increase and thermal stress effects are different, in each layer. Especially in indium antimonide chip, which is the uppermost layer of the indium antimonide detector that directly absorbs laser energy, the temperature shows not a smooth decrease as laser intensity from centre to outside, but a concentric-ringed ripple decrease with discontinuous high temperature extremum areas; which induces its own unique stress distributions: during the laser irradiation region, the thermal stress in areas above indium bumps is much lower than that in areas above underfill, but the stress distribution outside the laser irradiation region is exactly the opposite. The temperature and stress distribution in other materials are also different from that in indium antimonide chip. The main reason is that each material has different thermal properties, especially indium bump and underfill, which are alternative distribution in middle layers, have quite different thermal parameters. All these make the temperature and stress distribution in each layer having own characteristics. In addition, the change trend of the maximum temperature and stress in each layer in accordance with the laser pulse number are also studied. The maximum stress as well as temperature always occurs in the indium antimonide chip.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Numerical Analysis of Fatigue Behavior of ITER-Like Monoblock Divertor Interlayer Under Coupled Heat Loads

High-confinement mode is a very prominent operation style for future fusion device due to its unique advantages. However, the conjuncted edge localized modes (ELMs) are very difficult to control so that divertor plates are very prone to suffer both stationary high heat flux (HHF) loads of long-pulse operating mode and transient shock loads of ELMs. Most previous researches focus on degradation of plasma facing material (PFM), however, as a layer joining PFM and cooling tube, the soft copper interlayer suffers concentrated thermal stress loads due to mismatched thermal expansion of PFM and cooling tube. Its thermal fatigue behavior under such coupled loads is also of great significance to structural safety of divertor component. With such a motivation, the reduction effects on fatigue life time of a typical interlayer of monoblock divertor under series of coupled HHF and ELMs shock loading conditions are investigated. It is found that: (1) The transient shock feature of ELMs loading is propagated into interlayer with less sharp pattern. The increase of damage induced by coupled ELMs loading is limited in single cycle, while the accumulated damage of multiple consecutive coupled loading cycles is increased nonlinearly. (2) Under the coupled HHF and ELMs loading, the fatigue life time of interlayer is generally decreasing. The magnitude of decrease is increasing nonlinearly with the magnitude of ELMs peak and averaged heat flux. (3) For three characteristic parameters of ELMs shock loading such as frequency, duration and peak heat flux, the peak heat flux and frequency are two parameters more sensitive to determine coupled reduction effects on fatigue lifetime of the interlayer, while for high frequency case, time averaged heat flux takes the lead.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗