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,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Vibrations of Elastic Shells of Revolution Partially Filled with Ideal Liquid

We propose an algorithm for the numerical analyses of vibrations of elastic shells of revolution partially filled with ideal incompressible liquids. In the solution of this problem, the wave motions of liquid on its free surface are taken into account. The solution of the problem of hydroelasticity is based on the application of the method of decomposition of the domain of integration of equations of the theory of shells with the use of the variational statement of the problem and on the approximate construction of the operator inverse to the operator of the hydrodynamic part of the problem. We construct a generalized functional of displacements of the shell for which the role of conditions of matching for the solutions in different subdomains is played by the natural boundary conditions. The obtained numerical results are compared with the existing exact solutions of the analyzed problem for a shell in the form of a straight circular cylinder.

42 ENGINEERING↗
DOE OSTI2018

Numerical Method for Fractional Advection-Diffusion Equation with Heredity

We propose a method of construction of difference schemes for fractional partial differential equations with delay in time. For the fractional equation with two-sided diffusion, fractional transfer in time, and a functional aftereffect, we construct an implicit difference scheme. We use the shifted Grünwald–Letnikov formulas for the approximation of fractional derivatives with respect to spatial variables and the L1-algorithm for the approximation of fractional derivatives in time. Also we use piecewise constant interpolation and extrapolation by extending the discrete prehistory of the model in time. The algorithm is a fractional analog of a purely implicit method; on each time step, it is reduced to the solution of linear algebraic systems. We prove the stability of the method and find its order of convergence.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Specific Properties of One-Dimensional Pseudorepresentations of Groups

We obtain assertions concerning general properties of one-dimensional (not necessarily bounded) pseudorepresentations of groups. In particular, we obtain a quantitative condition on the numerical defect of a given pseudorepresentation which is sufficient for the pseudorepresentation to be pure, i.e., for the restriction of the given pseudorepresentation to every amenable subgroup be an ordinary character of this subgroup.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Modified Steklov Functions and Numerical Differentiation Formulas

We consider an approximation method based on Steklov functions of the first and second order. We obtain estimates for the norms in the space C of continuous periodic functions and clarify how they connect with numerical differentiation formulas.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Effect of Yagi–Uda nano-antenna element shape on the directivity and radiation efficiency

In this paper, different optimal designs of Yagi–Uda nano-antennas (NAs) are introduced and analyzed based on different element shapes to maximize the directivity and radiation efficiency. The studied element shapes include triangular, hexagonal, cube, square, rotated square, rectangular, ellipsoid, and elliptical cylinder shapes. The NAs parameters for each design are optimized using particle swarm optimization algorithm to achieve high directivity at a wavelength of 500 nm for wireless point-to-point applications. The designs are numerically analysed using 3-D finite difference time domain method (3-D FDTD) via Lumerical software package. To illustrate the performance of the proposed antennas, different radiation parameters such as radiation pattern, directivity, and radiation efficiency have been studied. Moreover, a tolerance of ± 5% is calculated to consider the imperfections that may appear in the fabrication process at a nanoscale. The optimized rectangular prism NA achieves the highest directivity of 24.21 with an efficiency of 50.74%. However, the optimized triangular prism NA has the highest efficiency of 73% with a directivity of 20.89. This enhancement is mainly attributed to the efficient coupling between array elements along with low side lobe level. However, a trade-off between directivity and radiation efficiency is obtained by the hexagonal prism NA with directivity of 21.9 and radiation efficiency of 69.1%.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2019

Refractive index from negative to positive with frequencies at the Dirac-like cone in a photonic crystal

In this paper, we numerically study the threefold accidental degeneracy conical dispersion (Driac-like cone) at the Brillouin zone center of the two-dimensional photonic crystals, which are composed of silicon pillars arranging in a triangular lattice. The effective permittivity and permeability near the Dirac-like point evolve from negative to positive by using the method of eigen-field averaging. Also, the isotropic behaviour of the Dirac-like cone is revealed by analysing the isofrequency contours. Moreover, we carry out numerical simulations including the reverse Snell’s law effect, negative Goos–Hänchen shifts and superfocusing lens to verify the negative refractive index characteristics of the designed structure. The proposed structure might find significant applications in the on-chip photonic interconnect and the photonic integrated circuit techniques.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Self-focusing of cosh-Gaussian laser beam and its effect on the excitation of ion-acoustic wave and stimulated Brillouin backscattering in collisionless plasma

An analytical and numerical study has been carried out for self-focusing of an intense cosh-Gaussian laser beam in collisionless plasma and its impact on the excitation of ion-acoustic wave and stimulated Brillouin backscattering process. The analytical model has been developed under Wentzel–Kramers–Brillouin and paraxial ray approximations. The nonlinearities of ponderomotive force on electron and the relativistic oscillation of the electron mass have been used in this study. The nonlinear differential equations have been set up for the beam width parameters of the main beam, ion-acoustic wave, backscattered wave and back reflectivity of stimulated Brillouin scattering (SBS). These equations have been solved numerically for different values of decentred parameter (b), relative plasma density (ω{sub p0}/ω{sub 0}) and incident laser intensity (a). The results have been compared with only relativistic nonlinearity and Gaussian profile of laser beam. The focusing of laser beam, ion-acoustic wave and scattered wave are found to be strong under relativistic-ponderomotive regime compared to only relativistic regime. Further, it is observed that focusing/intensity of main laser beam, ion acoustic wave and SBS back reflectivity increases with increasing the values of b and ω{sub p0}/ω{sub 0}. Itis also found that back reflectivity of SBS process gets suppressed with the increase in the value of a. This study may be useful in laser induced fusion scheme where back scattering of SBS plays very important role.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Thermal Modeling of the Injection of Standard and Thermally Insulated Cored Wire

Cored wire injection is a widespread method used to perform alloying additions during ferrous and non-ferrous liquid metal treatment. The wire consists of a metal casing that is tightly wrapped around a core of material; the casing delays the release of the material as the wire is immersed into the melt. This method of addition presents advantages such as higher repeatability and yield of cored material with respect to bulk additions. Experimental and numerical work has been performed by several authors on the subject of alloy additions, spherical and cylindrical geometries being mainly considered. Surprisingly this has not been the case for cored wire, where the reported experimental or numerical studies are scarce. This work presents a 1-D finite volume numerical model aimed for the simulation of the thermal phenomena which occurs when the wire is injected into a liquid metal bath. It is currently being used as a design tool for the conception of new types of cored wire. A parametric study on the effect of injection velocity and steel casing thickness for an Al cored wire immersed into a steel melt at 1863 K (1590 °C) is presented. The standard single casing wire is further compared against a wire with multiple casings. Numerical results show that over a certain range of injection velocities, the core contents’ release is delayed in the multiple casing when compared to a single casing wire.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Magnesium Loss in Nd:YAG Pulsed Laser Welding of Aluminum Alloys

In the current investigation, the loss of magnesium (Mg) in Nd:YAG pulsed laser welding of a number of aluminum (Al) alloys was studied both experimentally and numerically. The experimental results showed that at a fixed average laser power, with the increasing pulse frequency, Mg concentration of the weld metals decreased. A numerical model was developed for the calculation of Mg evaporation from the weld pool. The model predicted a trend in the changes of Mg concentration of the weld pool, working best for alloys containing the higher levels of Mg. The calculations showed that with the increasing pulse frequency from 15 to 40 Hz, due to a decrease in the average temperature of the weld pool, Mg loss in terms of mass per unit area of the weld spot actually decreased from 230 × 10{sup −5} to 120 × 10{sup −5} kg/m{sup 2}. However, due to a decrease in the weld penetration/volume, the effective Mg loss of the weld metal in terms of concentration increased from 110 × 10{sup −1} to 260 × 10{sup −1} kg/m{sup 3}. It was shown that due to interactions between Mg and other main alloying elements, the activity coefficient of Mg had to be modified for the model to accurately predict Mg evaporation of the weld pool.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Skin and Proximity Effects in Electrodes and Furnace Shells

A review of two-dimensional (2D) analytical models of skin and proximity effects in large industrial furnaces with three electrodes arranged in an equilateral triangle is given. The models cover three different cases: one electrode only, three electrodes where two are approximated by line currents, and induced shell currents where all electrodes are approximated by line currents. The first two models show how the skin and proximity effects depend on electrode material properties and size, and the distance between the electrodes. The third model shows how the strength of the induced shell currents will depend on electrode position and furnace size. These models are compared to numerical studies including distributed electrodes and shell currents. The analytical models are accurate when induced shell currents can be disregarded. However, strong shell currents may have a significant impact on the current distribution within the electrodes. This electrode-shell proximity effect competes with the electrode-electrode proximity effect. Finally, the 2D models have been compared with three-dimensional (3D) case studies of large industrial furnaces. In 3D, the shell currents are significantly smaller than what are predicted by the 2D models, but they are sufficiently strong to cause a significant correction of the electrode current density.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Experimental and Numerical Analyses of the Consolidation Process of AA 7075–2 wt.% ZrO{sub 2} Powders by Equal Channel Angular Pressing

AA7075 matrix composite was obtained by mechanical alloying using ZrO{sub 2} particles as reinforcement. After that, the powders were cold-compacted into an aluminum tube and consolidated by equal channel angular pressing at 220 °C. The temperature influence was analyzed, and it was found by scanning electron microscopy that this route gave the minimum pore percentage. Besides, it was noticeable that the hardness was wide superior in comparison with the AA 7075 aged bulk alloy. In addition, there was no evidence of any chemical reaction between ZrO{sub 2} and alloying elements at the temperature employed. On the other hand, it was numerically and experimentally demonstrated that the more strained regions exhibited a higher relative density values.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Experimental Investigation and Numerical Study on the Elimination of Porosity in Aluminum Alloy Laser Welding and Laser–GMA Welding

Keyhole-induced porosity can strongly weaken the joint strength during aluminum alloy laser welding. Two factors are proposed to be responsible for the porosity formation: bubble formation caused by the instability of the keyhole and bubble captured by the solidification interface. Based on numerical and experimental results, the porosity is suppressed with preheating, due to the capillary stability of the keyhole and easy escape of bubble. N{sub 2} shielding gas has a minor effect on keyhole stability, but it reduces the porosity number as N{sub 2} bubble can dissolve into the liquid aluminum alloy melt pool quickly. The laser–gas metal arc (laser–GMA) welding process can stabilize the metallic plume and enlarge the weld pool. As a result, few spatters and bubbles are formed, and bubbles are easily escaped, both of which contribute to the suppression of porosity.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Modeling of Keyhole-Induced Pore Formation in Laser-Arc Hybrid Welding of Aluminum Alloy with a Horizontal Fillet Joint

A three-dimensional transient model is proposed to investigate the weld pool dynamic behavior in laser + metal inert gas (MIG) hybrid fillet welding of aluminum alloy in the horizontal position, which allows for the joint configuration and coupling of the keyhole, droplet and weld pool as well as the heat and mass exchange between gas and liquid phases and is able to simulate the temperature distribution, fluid flow and formation process of a keyhole-induced pore in hybrid welding with a horizontal fillet joint. The weld porosity is also measured using x-ray nondestructive testing technology. Keyhole behavior and the formation mechanism of keyhole-induced porosity were analyzed. The calculated results are in generally good agreement with the experimental ones. A clockwise vortex always exists at the middle part of the weld pool. The formation and growth of the molten metal bulge on the keyhole wall are responsible for the occurrence of a gas bubble, which has a variation in size and shape and can be split during welding. The keyhole collapses easily at its middle or upper part in horizontal fillet welding, and the capture of the bubble by the upper molten pool boundary enhances the possibility of porosity formation to some degree. The keyhole-induced pore is mainly formed at the regions near the keyhole bottom and the upper fusion line of the weld pool.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Effect of natural fractures on determining closure pressure

Abstract In Utah FORGE reservoir, eight pressure transient tests (microhydraulic fracturing and DFIT™) show natural fracture/pressure-dependent leakoff. This behavior may lead to misinterpretation of the closure pressure (proxy for minimum principal stress). The closure pressures obtained from DFIT™ (or microhydraulic fracturing) testing may reflect shear failure along natural fractures or discontinuities rather than tensile failure and lead to inaccurate estimates of the minimum principal stress. In pressure tests conducted at the Utah FORGE site showed that reactivation or opening of natural fractures intersecting the wellbore and were suggested by multiple closure events in DFIT™ tests. In addition, comparisons between the pre- and post-well-tests FMI logs show that there are populations of induced fracture and reactivated natural fractures which are mostly vertical and sub-vertical. In this study, DFIT™ (or microhydraulic fracturing) test analysis and numerical simulations were used to suggest that hydraulic shearing of critically stressed natural fractures can contribute to multiple closure signatures and possibly lead to incorrect determination of the minimum principal stress. As other authors have previously advocated, the DFIT™ test analyses and numerical simulations suggest that better estimations of minimum principal stress may be derived by injecting at relatively high rate and pressure and insuring that tensile breakdown is reached. In addition, extended shut-in period is required to determine accurate reservoir characteristics and fluid flow regime.

Nadimi, S.↗
DOE OSTI2018

Large-eddy simulations for hill terrains: validation with wind-tunnel and field measurements

Development and validation of large-eddy simulation (LES) framework to study atmospheric boundary layer flows over complex terrains is reported here. The LES model uses the fourth-order time-accurate Runge–Kutta scheme and a fractional step method. The inflow boundary conditions are generated by using the so-called recycling (or mapping) technique. Evaluation of potential in-land wind park locations is the main application of the validated model. In this work, LES are carried out for turbulent boundary-layer flows over both simple and complex hill geometries. The prediction of the flow separation and reattachment length for a steeper wind-tunnel hill was closer to the measurements than previous numerical studies reported by other authors for the same hill geometry. For the complex hill case, the LES model showed good agreement with full-scale measurements.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Multiobjective optimization to a TB-HIV/AIDS coinfection optimal control problem

We consider a recent coinfection model for Tuberculosis (TB), Human Immunodeficiency Virus (HIV) infection, and Acquired Immunodeficiency Syndrome (AIDS) proposed in Silva and Torres (Discr Contin Dyn Syst 35(9):4639–4663, 2015). We introduce and analyze a multiobjective formulation of an optimal control problem, where the two conflicting objectives are minimization of the number of HIV-infected individuals with AIDS clinical symptoms and coinfected with AIDS and active TB; and costs related to prevention and treatment of HIV and/or TB measures. The proposed approach eliminates some limitations of previous works. The results of the numerical study provide comprehensive insights about the optimal treatment policies and the population dynamics resulting from their implementation. Some nonintuitive conclusions are drawn. Overall, the simulation results demonstrate the usefulness and validity of the proposed approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Numerical study of a smoothing algorithm for the complementarity system over the second-order cone

This paper deals with the complementarity system over the second-order cone (denoted by CSSOC) which contains a wide class of problems. We extend a class of regularized Chen–Harker–Kanzow–Smale smoothing functions studied by Huang and Sun (Appl Math Optim 52:237–262, 2005) for the linear complementarity problem to the CSSOC. Based on this class of functions, we propose a smoothing algorithm for solving the CSSOC. Under weak assumptions, we prove that the proposed algorithm has global and local quadratic convergence. The proposed algorithm is different from existing smoothing algorithms for solving the CSSOC because it adopts a new nonmonotone line search rule. In addition, our algorithm solves a new equation reformulation of the CSSOC. Numerical experiments indicate that the proposed algorithm is quite effective.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Oblique water wave damping by two submerged thin vertical porous plates of different heights

Oblique water wave damping by two fully submerged vertical parallel porous plates of different heights is investigated within the framework of linear water wave theory. A channel of infinite horizontal extent but of finite depth is divided into three regions in each of which three different velocity potentials are considered and subsequently the corresponding boundary value problems are formulated. By means of eigenfunction expansion and least square method, complete analytical solution in each of the regions is obtained. With the help of the matching conditions along the vertical boundaries between any two successive regions, the reflection and transmission coefficients, and hence the energy loss (in %), are obtained numerically by applying a matrix method. Numerical study is carried out for various relevant parameters and it is found that with an appropriate choice of the parameters and positioning of the porous plates, both of the reflection and transmission coefficients become significantly low and consequently the porous plates can be utilized as an effecting wave absorber. Comparison of this work with earlier available work is made, and the good agreement in the results suggests the applicability of this method in further engineering purposes. A particular case of two free surface piercing plates is also investigated on similar line and the results are presented in order to observe the difference with the earlier case.

97 MATHEMATICS AND COMPUTING↗