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 OSTI2018

Numerical Analysis of the Effect of Infrared Radiation on Cryogenic Inertial Confinement Fusion Targets

The present work applies the finite element method to calculate the maximum allowable time that cryogenic inertial confinement fusion (ICF) targets can be exposed to infrared radiation (IR). Hence, a 3-D numerical model integrated with discrete coordinate radiation model was developed to investigate the influence of transmittance of the laser entrance holes (LEHs) and boundary conditions on the temperature field distribution and the maximum DT layer deterioration time for CH, Be, and diamond capsules. Our study shows that introducing such a radiation model can accurately obtain more detailed spatial and temporal distribution information in the ICF targets. The simulation results demonstrate that the Be and diamond capsules provided much better temperature field homogenization than the CH capsule under equivalent boundary conditions, but the CH capsule was heated more by IR radiation than the Be and diamond. In addition, the maximum DT layer deterioration time was significantly increased to 3 s when decreasing the transmittance of the LEH from 0.2 to 0.01. However, either reducing the capsule IR absorption or increasing the inner hohlraum IR absorption demonstrated no conclusive increase in the maximum DT layer deterioration time. These results are expected to provide useful parameters in the design of cryogenic targets and shroud systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Numerical analysis of amplification of picosecond pulses in a THL-100 laser system with an increase in the pump energy of the XeF(C – A) amplifier

Energy characteristics of laser radiation with a pulse width of 50 ps at an elevated pump energy of the XeF(C – A) amplifier of a hybrid THL-100 laser system are analysed numerically. The dynamics of the change in the energy and maximum intensity of laser radiation with an increase in the pump energy of the XeF(C – A) amplifier from 270 to 400 J is investigated. The results of studying the influence of the input beam divergence on the energy characteristics of the output beam are presented. It is shown that, for the existing system of mirrors, an increase in the pump energy to 400 J leads to an increase in the output energy from 3.2 to 5.5 J at a maximum radiation intensity of 57 GW cm{sup −2}. A system of amplifier mirrors with 27 laser beam passes and enlarged divergence angle of the amplified beam is considered. Theoretically, the proposed system of mirrors allows one to increase the laser pulse energy to 7.5 J at a maximum intensity of no more than 14.8 GW cm{sup −2}. The calculated efficiency of the conversion of the pump energy absorbed in the amplifier gas chamber into the lasing energy exceeds 3 % in this regime. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2015

Numerical analysis of sandwich beam with corrugated core under three-point bending

The strength problem of sandwich beam with corrugated core under three-point bending is presented.The beam are made of steel and formed by three mutually orthogonal corrugated layers. The finite element analysis (FEA) of the sandwich beam is performed with the use of the FEM system - ABAQUS. The relationship between the applied load and deflection in three-point bending is considered.

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 OSTI2014

Numerical analysis of micro-/nanoscale gas-film lubrication of sliding surface with complicated structure

It has been reported that the friction between a partially polished diamond-coated surface and a metal surface was drastically reduced to zero when they are slid at a few m/s. Since the sliding was noiseless, it seems that the diamond-coated surface was levitated over the counter surface and the sliding mechanism was the gas film lubrication. Recently, the mechanism of levitation of a slider with a micro/nanoscale surface structure on a rotating disk was theoretically clarified [S. Yonemura et al., Tribol. Lett., (2014), doi:10.1007/s11249-014-0368-2]. Probably, the partially polished diamond-coated surface may be levitated by high gas pressure generated by the micro/nanoscale surface structure on it. In this study, in order to verify our deduction, we performed numerical simulations of sliding of partially polished diamond-coated surface by reproducing its complicated surface structure using the data measured by an atomic force microscope (AFM). As a result, we obtained the lift force which is large enough to levitate the slider used in the experiment.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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.↗
DOE OSTI2014

Numerical analysis of flows of rarefied gases in long channels with octagonal cross section shapes

Isothermal, pressure driven rarefied gas flows through long channels with octagonal cross section shapes are analyzed computationally. The capillary is between inlet and outlet reservoirs. The cross section is constant along the axial direction. The boundary condition at the solid-gas interface is assumed to be diffuse reflection. Since the channel is long, the gaseous velocity is small compared to the average molecular speed. Consequently, a linearized description can be used. The flow is described by the linearized Bhatnagar-Gross-Krook kinetic model. The solution of the problem is divided into two stages. First, the local flow field is determined by assuming the local pressure gradient. Secondly, the global flow behavior is deduced by the consideration of the conservation of the mass along the axis of the capillary. The kinetic equation is solved by the discrete velocity method on the cross section. Both spatial and velocity spaces are discretized. A body fitted rectangular grid is used for the spatial space. Near the boundary, first-order, while in the interior part of the flow domain, second-order finite-differences are applied to approximate the spatial derivatives. This combination results into an efficient and straightforward numerical treatment. The velocity space is represented by a Gauss-Legendre quadrature. The kinetic equation is solved in an iterative manner. The local dimensionless flow rate is calculated and tabulated for a wide range of the gaseous rarefaction for octagonal cross sections with various geometrical parameters. It exhibits the Knudsen minimum phenomenon. The flow rates in the octagonal channel are compared to those through capillaries with circular and square cross sections. Typical velocity profiles are also shown. The mass flow rate and the distribution of the pressure are determined and presented for global pressure driven flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗