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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
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↗
DOE OSTI2026

Performing Numerical Analysis of Cybersecurity Options Using Dynamic Risk Analysis Tool EMRALD

Cyberattacks can have many different attack paths, durations, and goals. There are also many different mitigation options involving hardware, software, and/or humans. Considering a cyber threat should involve defense-in-depth methods and a quantitative or numerical evaluation of overall effectiveness against dynamic, time-dependent attacks to make cost and risk-informed decisions. Typical cyberattack modeling methods only provide a qualitative evaluation. The main areas of cybersecurity are confidentiality, integrity, and availability. For companies with cyber-physical systems such as advanced nuclear reactors, cyber-related safety is a requirement set by North American Electric Reliability and the U.S. Nuclear Regulatory Commission. They are also concerned about availability or reliability as a business case. As cyber threats are evolving to a business-for-hire structure, more attacks may focus on disrupting business success and reliability, causing financial and economic stability risk. Companies want to know business reliability and recovery from those threats, and that requires modeling physical behavior of the targets. Dynamic-state-based and Markov-based modeling provides a method for better cyber scenario modeling with different tools having issues such as state-base explosion. Dynamic modeling enables time and conditional features not found in other numerical evaluation methods. EMRALD (Event Modeling Risk Assessment using Lined Diagrams) is a dynamic risk analysis modeling and simulation tool and has features that reduce modeling issues. It has been used to model different time-dependent events including plant behavior and operator procedures. As a general modeling tool, EMRALD can also be used to model cyberattack scenarios with varying mitigation options and quantify effectiveness, producing numerical data for risk-informed decisions. This paper uses EMRALD to demonstrate that dynamic numerical risk analysis can be used for cyber threat modeling to provide insights for design decision-making and optimize defense strategies. Keywords: cyber modeling; cyber-physical systems; numerical cyber modeling

97 - MATHEMATICS AND COMPUTING↗
NASA NTRS1986

Numerical analysis of some supersonic viscous flows related to inlet and nozzle systems

A numerical method originally developed for three-dimensional supersonic inlet flow calculations is extended and applied to the study of two-dimensional and three-dimensional flows associated with arbitrary propulsion systems. The method is based on the forward spatial marching solution of a reduced form of the three-dimensional steady Navier-Stokes equations in which streamwise pressure gradients are retained in both the subsonic and supersonic regions. The present paper briefly describes the analysis and then shows three applications. In the first application, a wall transpiration study has been performed for the two-dimensional shock wave/turbulent boundary layer interaction flow field with application to an inlet configuration. The second application treats a rectangular high speed inlet with a swept sideplate including the effects of sideplate spillage. Finally, the method is utilized to analyze the interaction of an under-expanded supersonic jet with an ambient flow. Computed results are examined and compared with available experimental measurements. It is demonstrated that the present numerical method is capable of numerically simulating complex two- and three-dimensional flows relevant to hypersonic propulsion systems in a manner which both shows good agreement with data when such data is available, and which shows the complex flow features in the absence of data.

Kim, Y.-N.↗
NASA NTRS2008

Numerical Analysis of Rocket Exhaust Cratering

Supersonic jet exhaust impinging onto a flat surface is a fundamental flow encountered in space or with a missile launch vehicle system. The flow is important because it can endanger launch operations. The purpose of this study is to evaluate the effect of a landing rocket s exhaust on soils. From numerical simulations and analysis, we developed characteristic expressions and curves, which we can use, along with rocket nozzle performance, to predict cratering effects during a soft-soil landing. We conducted a series of multiphase flow simulations with two phases: exhaust gas and sand particles. The main objective of the simulation was to obtain the numerical results as close to the experimental results as possible. After several simulating test runs, the results showed that packing limit and the angle of internal friction are the two critical and dominant factors in the simulations.

Source record↗
DOE OSTI2017

Numerical Analysis of Experiments on Thermally Induced Dissociation of Methane Hydrates in Porous Media

Numerical simulation is essential for the prediction and evaluation of hydrocarbon reservoir performance. Numerical simulators developed for the description of the behavior of hydrates under production and the corresponding flow of fluids and heat accounting for all known processes are powerful, but they need validation through comparison to field or experimental data in order to instill confidence in their predictions. In this study, we analyze by means of numerical simulation the results of an experiment of methane hydrate dissociation by thermal stimulation in unconsolidated porous media heated through the vessel walls. The physics captured by the model include multicomponent heat and mass transfer, multiphase flow through porous media, and the phase behavior of the CH 4 + H 2 O system involved in methane hydrate formation and dissociation. The set of governing equations consists of the mass and energy conservation equations coupled with constitutive relationships, i.e., the dissolution of gas in H 2 O, relative permeability and capillary pressure models, composite thermal conductivity models, and methane hydrate phase equilibria. The model geometry describes accurately the hydrate reactor used in a recent experimental study investigating methane hydrate dissociation behavior [Chong et. al. Appl. Energy 2016, 177, 409-421]. The cumulative gas production is estimated and validated against three tests of experimental data involving different boundary temperatures, showing a good agreement between observations and numerical predictions. The predicted evolution of the spatial distributions of different phases over time shows that hydrate dissociation progresses inward from the reactor boundary to the center, methane gas accumulates to the top of the reactor because of buoyancy, and water migrates down to the bottom of the reactor because of gravity. A sharp hydrate dissociation front is predicted, and the estimated location of hydrate dissociation front suggests a linear relationship with the square root of time. A sensitivity analysis on the thermal conductivity of sand under fully saturated conditions is conducted to elucidate its effect on the gas production behavior. In addition, the energy efficiency ratio computed from the simulation of this boundary-wall heating technique varies from 14.0 to 16.2. Deviations between observations and predictions of the evolution of the temperature profile are attributed to initial heterogeneous distribution of the hydrate phase in the hydrate reactor.

42 ENGINEERING↗
DOE OSTI2015

Numerical Analysis of Coolant Flow and Heat Transfer in ITER Diagnostic First Wall

We performed numerical simulations of the ITER Diagnostic First Wall (DFW) using ANSYS workbench. During operation DFW will include solid main body as well as liquid coolant. Thus thermal and hydraulic analysis of the DFW was performed using conjugated heat transfer approach, in which heat transfer was resolved in both solid and liquid parts, and simultaneously fluid dynamics analysis was performed only in the liquid part. This approach includes interface between solid and liquid part of the systemAnalysis was performed using ANSYS CFX software. CFX software allows solution of heat transfer equations in solid and liquid part, and solution of the flow equations in the liquid part. Coolant flow in the DFW was assumed turbulent and was resolved using Reynolds averaged Navier-Stokes equations with Shear Stress Transport turbulence model. Meshing was performed using CFX method available within ANSYS. The data cloud for thermal loading consisting of volumetric heating and surface heating was imported into CFX Volumetric heating source was generated using Attila software. Surface heating was obtained using radiation heat transfer analysis. Our results allowed us to identify areas of excessive heating. Proposals for cooling channel relocation were made. Additional suggestions were made to improve hydraulic performance of the cooling system.

42 ENGINEERING↗
NASA NTRS1984

Manufacturing in space: Fluid dynamics numerical analysis

Thermocapillary convection in an air-jet Marangoni flow suppressing technique was numerically simulated for flow under various conditions of fluid properties, gravity, and temperature gradient. Results presented include computer generated plots of streamlines, velocities, and temperatures throughout the contained fluid flow, and plots of Marangoni flow velocities over the free surface.

Robertson, S. J.↗
DOE OSTI2015

Numerical analysis for finite-range multitype stochastic contact financial market dynamic systems

In an attempt to reproduce and study the dynamics of financial markets, a random agent-based financial price model is developed and investigated by the finite-range multitype contact dynamic system, in which the interaction and dispersal of different types of investment attitudes in a stock market are imitated by viruses spreading. With different parameters of birth rates and finite-range, the normalized return series are simulated by Monte Carlo simulation method and numerical studied by power-law distribution analysis and autocorrelation analysis. To better understand the nonlinear dynamics of the return series, a q-order autocorrelation function and a multi-autocorrelation function are also defined in this work. The comparisons of statistical behaviors of return series from the agent-based model and the daily historical market returns of Shanghai Composite Index and Shenzhen Component Index indicate that the proposed model is a reasonable qualitative explanation for the price formation process of stock market systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1977

Numerical analysis of free vibrations of damped rotating structures

This paper is concerned with the efficient numerical solution of damped and undamped free vibration problems of rotating structures. While structural discretization is achieved by the finite element method, the associated eigenproblem solution is effected by a combined Sturm sequence and inverse iteration technique that enables the computation of a few required roots only without having to compute any other. For structures of complex configurations, a modal synthesis technique is also presented, which is based on appropriate combinations of eigenproblem solution of various structural components. Such numerical procedures are general in nature, which fully exploit matrix sparsity inherent in finite element discretizations, and prove to be most efficient for the vibration analysis of any damped rotating structure, such as rotating machineries, helicopter and turbine blades, spinning space stations, among others.

Gupta, K. K.↗
NASA NTRS2001

A Hybrid Numerical Analysis Method for Structural Health Monitoring

A new hybrid surface-integral-finite-element numerical scheme has been developed to model a three-dimensional crack propagating through a thin, multi-layered coating. The finite element method was used to model the physical state of the coating (far field), and the surface integral method was used to model the fatigue crack growth. The two formulations are coupled through the need to satisfy boundary conditions on the crack surface and the external boundary. The coupling is sufficiently weak that the surface integral mesh of the crack surface and the finite element mesh of the uncracked volume can be set up independently. Thus when modeling crack growth, the finite element mesh can remain fixed for the duration of the simulation as the crack mesh is advanced. This method was implemented to evaluate the feasibility of fabricating a structural health monitoring system for real-time detection of surface cracks propagating in engine components. In this work, the authors formulate the hybrid surface-integral-finite-element method and discuss the mechanical issues of implementing a structural health monitoring system in an aircraft engine environment.

Forth, Scott C.↗
DOE OSTI2015

Numerical analysis of modified Central Solenoid insert design

The United States ITER Project Office (USIPO) is responsible for fabrication of the Central Solenoid (CS) for ITER project. The ITER machine is currently under construction by seven parties in Cadarache, France. The CS Insert (CSI) project should provide a verification of the conductor performance in relevant conditions of temperature, field, currents and mechanical strain. The US IPO designed the CSI that will be tested at the Central Solenoid Model Coil (CSMC) Test Facility at JAEA, Naka. To validate the modified design we performed three-dimensional numerical simulations using coupled solver for simultaneous structural, thermal and electromagnetic analysis. Thermal and electromagnetic simulations supported structural calculations providing necessary loads and strains. According to current analysis design of the modified coil satisfies ITER magnet structural design criteria for the following conditions: (1) room temperature, no current, (2) temperature 4K, no current, (3) temperature 4K, current 60 kA direct charge, and (4) temperature 4K, current 60 kA reverse charge. Fatigue life assessment analysis is performed for the alternating conditions of: temperature 4K, no current, and temperature 4K, current 45 kA direct charge. Results of fatigue analysis show that parts of the coil assembly can be qualified for up to 1 million cycles. Distributions of the Current Sharing Temperature (TCS) in the superconductor were obtained from numerical results using parameterization of the critical surface in the form similar to that proposed for ITER. Lastly, special ADPL scripts were developed for ANSYS allowing one-dimensional representation of TCS along the cable, as well as three-dimensional fields of TCS in superconductor material. Published by Elsevier B.V.

42 ENGINEERING↗
NASA NTRS1982

Numerical analysis of Weyl's method for integrating boundary layer equations

A fast method for accurate numerical integration of Blasius equation is proposed. It is based on the limit interchange in Weyl's fixed point method formulated as an iterated limit process. Each inner limit represents convergence to a discrete solution. It is shown that the error in a discrete solution admits asymptotic expansion in even powers of step size. An extrapolation process is set up to operate on a sequence of discrete solutions to reach the outer limit. Finally, this method is extended to related boundary layer equations.

Najfeld, I.↗
DOE OSTI2018

Numerical analysis of experimental studies of methane hydrate formation in a sandy porous medium

Here we analyse numerically an earlier experimental study that involved the formation of methane hydrates by the excess water method in a small reactor filled with a sandy porous medium, and seek to address questions about the type of the hydration reaction and the phase heterogeneity in the resulting hydrate-bearing sand. Using a fine discretization describing the reactor assembly, the experimental process is faithfully replicated numerically. The multi-stage process of hydrate formation is subdivided in 7 steps. The experimental data from the continuously-monitored pressure and temperature during each step are used for comparison against the numerical predictions, the identification of the dominant processes and the determination of the associated parameters through a history-matching process that minimizes deviations between observations and simulation results. The results of this first-ever study on this subject demonstrate unequivocally that the hydration reaction is a kinetic (as opposed to an equilibrium) process, and that the spatial distributions of the various phases (aqueous, gas and hydrate) at the end of the formation process are strongly heterogeneous. This has serious implications in simulation studies of hydrate dissociation that assume uniform initial phase saturation distributions. Finally, the history-matching process indicates that (a) the system behaviour is sensitive to some flow parameters (porosity and irreducible water saturation) only during the first water injection, (b) it is insensitive to the sand intrinsic permeability during all steps of the study, and (c) thermal processes dominate after the first water injection, yielding estimates of the thermal properties of the sand and of time-variable key parameters of the kinetic reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1980

Numerical analysis of natural convection in two-dimensional square and circular containers in low gravity

A numerical study of natural convection in circular cylinder and square enclosures shows that the analytic low Rayleigh number theory of previous investigators is valid for Rayleigh numbers up to 1000. For a Rayleigh number of 5000, steady state values of maximum fluid velocity differ by 20 percent. This deviation between analytic theory and numerical results increases for higher Rayleigh numbers. In addition, the low Rayleigh number theory is shown to be valid for higher Rayleigh numbers for a portion of the transient phase before significant deviation becomes apparent. It is also shown that square shaped experimental configurations may be analytically approximated with good accuracy by circular cylinders of equal cross sectional area for the prediction of convection velocities and flow patterns at low Rayleigh number.

Robertson, S. J.↗