Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite element methods”

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

Finite element solution of MHD power-law fluid with slip velocity effect and non-uniform heat source/sink

The laminar steady boundary layer flow of an incompressible MHD Power-law fluid past a continuously moving surface is investigated numerically. The study involves the influence of surface slip and non-uniform heat source/sink on the flow and heat transfer. The governing boundary layer flow equations are transformed into non-dimensional, non-linear coupled ordinary differential equations with the help of suitable similarity transformations. The Galerkin finite element method is implemented to crack the resulting system. The impact of different involved physical parameters is exhibited on the dimensionless velocity profile, temperature distributions and rate of heat transfer in graphical and tabular forms for pseudoplastic and dilatant fluids. The local Nusselt number is found to be the decreasing function of slip parameter, temperature and space-dependent heat sink parameter whereas it increases with increasing values of temperature and space dependent heat source parameter. The problem has important application in attaining the sustainable heat transfer rate for the cooling of fluids, especially in heat exchangers used frequently in chemical industry, in order to increase the trustworthiness of a system, as it removes high heat loads from these systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Smoothed boundary method for simulating incompressible flow in complex geometries

Simulating flow through porous media with explicit considerations of complex microstructures is very challenging using conventional sharp-interface methods because of the difficulties in generating meshes conformal to complex geometries. In this work, a diffuse interface embedded boundary method known as the Smoothed Boundary Method (SBM) is utilized to facilitate simulations of fluid dynamics involving complex geometries. In diffuse-interface methods, the geometry is described by a domain parameter. The SBM allows the straightforward reformulation of the time-dependent Navier-Stokes equations in terms of this domain parameter, using only algebraic identities. Thus, enforcing the appropriate boundary conditions at the irregular embedded boundary is greatly simplified. Adaptive mesh refinement is used to increase the accuracy of the diffuse interface method by allowing thinner interfacial thickness to be used in the domain parameter. Here, the SBM-formulated Navier-Stokes equations are solved with the Finite Difference Method on refined mesh systems. Sharp-interface Finite Element Method simulations using the commercial software COMSOL on body-conforming meshes are also provided for comparison. Favorable agreement between the two methods is observed. Since it is no longer necessary for the mesh to conform to the complex geometry, the grid system for the SBM simulations can be generated rapidly and without additional manual interventions, making the entire simulation process more expedient.

36 MATERIALS SCIENCE↗
DOE OSTI2016

High Order Finite Elements S{sub N} Transport in X-Y Geometry on Meshes with Curved Surfaces

Radiation hydrodynamics is useful to model astrophysics, ICF, and glass manufacturing. A large amount of research has been conducted to understand how to solve hydrodynamics equations using high-order finite elements. Researchers at LLNL have developed a general high-order FEM library, Modular Finite Element Methods (MFEM), to enable this research. The BLAST code uses MFEM to evolve the Lagrangian hydrodynamics equations, resulting in meshes with higher order curved surfaces. Our research focuses on the requirements for high order finite element solutions of the radiation transport equation with the eventual goal of studying the complete high-order radiation hydrodynamics solution methodology. For this research we used MFEM to develop a proof-of-concept code that solves the high order finite element transport equations in Cartesian coordinates on meshes with curved surfaces. We have solved several test problems to assess the behavior of our method. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Bandgaps and directional properties of two-dimensional square beam-like zigzag lattices

In this paper we propose four kinds of two-dimensional square beam-like zigzag lattice structures and study their bandgaps and directional propagation of elastic waves. The band structures are calculated by using the finite element method. Both the in-plane and out-of-plane waves are investigated simultaneously via the three-dimensional Euler beam elements. The mechanism of the bandgap generation is analyzed by studying the vibration modes at the bandgap edges. The effects of the geometry parameters of the xy- and z-zigzag lattices on the bandgaps are investigated and discussed. Multiple complete bandgaps are found owing to the separation of the degeneracy by introducing bending arms. The bandgaps are sensitive to the geometry parameters of the periodic systems. The deformed displacement fields of the harmonic responses of a finite lattice structure subjected to harmonic loads at different positions are illustrated to show the directional wave propagation. An extension of the proposed concept to the hexagonal lattices is also presented. The research work in this paper is relevant to the practical design of cellular structures with enhanced vibro-acoustics performance.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Geometry optimization for micro-pressure sensor considering dynamic interference

Presented is the geometry optimization for piezoresistive absolute micro-pressure sensor. A figure of merit called the performance factor (PF) is defined as a quantitative index to describe the comprehensive performances of a sensor including sensitivity, resonant frequency, and acceleration interference. Three geometries are proposed through introducing islands and sensitive beams into typical flat diaphragm. The stress distributions of sensitive elements are analyzed by finite element method. Multivariate fittings based on ANSYS simulation results are performed to establish the equations about surface stress, deflection, and resonant frequency. Optimization by MATLAB is carried out to determine the dimensions of the geometries. Convex corner undercutting is evaluated. Each PF of the three geometries with the determined dimensions is calculated and compared. Silicon bulk micromachining is utilized to fabricate the prototypes of the sensors. The outputs of the sensors under both static and dynamic conditions are tested. Experimental results demonstrate the rationality of the defined performance factor and reveal that the geometry with quad islands presents the highest PF of 210.947 Hz{sup 1/4}. The favorable overall performances enable the sensor more suitable for altimetry.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2016

Thermo-mechanical simulations of CO{sub 2} laser–fused silica interactions

CO{sub 2} laser heating of silica glass is used in many scientific and industrial applications. Particularly, localized CO{sub 2} laser heating of silica glass has demonstrated its ability to mitigate surface damage on optics used for high power laser applications. To develop such applications, the control of temperature, heat affected area, and resulting mechanical stresses are critical. Therefore, it is necessary to understand the silica transformation, the material ejection, and the thermo-mechanical stresses induced by the laser heating and subsequent cooling. In this paper, we detail the development of comprehensive thermo-mechanical numerical simulations of these physical processes, based on finite-element method. The approach is developed for 2D or 3D cases to tackle the case of a moving beam at the surface of the sample, and we particularly discuss the choice of the different parameters based on bibliographic inputs. The thermal and mechanical numerical results have been compared to different dedicated experimental studies: infrared thermography measurements at the surface of the irradiated area, optical profilometry measurements of the laser-processed sites, and photo-elastic measurements. Very consistent results are obtained between numerical and experimental results for the description of the temperature gradients, the material ejection, and the residual stresses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Opposing flow in square porous annulus: Influence of Dufour effect

Heat and mass transfer in porous medium is very important area of research which is also termed as double diffusive convection or thermo-solutal convection. The buoyancy ratio which is the ratio of thermal to concentration buoyancy can have negative values thus leading to opposing flow. This article is aimed to study the influence of Dufour effect on the opposing flow in a square porous annulus. The partial differential equations that govern the heat and mass transfer behavior inside porous medium are solved using finite element method. A three node triangular element is used to divide the porous domain into smaller elements. Results are presented with respect to geometric and physical parameters such as duct diameter ratio, Rayleigh number, radiation parameter etc. It is found that the heat transfer increase with increase in Rayleigh number and radiation parameter. It is observed that Dufour coefficient has more influence on velocity profile.

42 ENGINEERING↗
DOE OSTI2018

Fast Algorithm for Evaluating Critical Current of High-Temperature Superconducting Pancake Coil

Second-generation (2G) high-temperature superconducting (HTS) tapes are now capable of carrying very high transport current and promising for a wide range of applications. The critical current of HTS coils is important for applications, such as superconducting electric machines, superconducting magnetic energy storage, and superconducting magnets. Therefore, precisely and quickly calculating critical current of HTS coils is very important for designing HTS devices. This paper provides a fast algorithm for evaluating critical current of HTS pancake coil. The fast algorithm is realized through a stationary model, which is based on finite element method (FEM) software. The stationary model means that the model is solved by stationary study instead of time-dependent study. To validate this method, a pancake HTS coil was wound and its critical current was measured. Meanwhile, an axial symmetric stationary model was built according to the geometry of the measured HTS coil. By comparing measured and calculated results, the effectiveness of the stationary model was demonstrated. Moreover, the stationary model is compared with H formulation model. The calculated results by the two models are nearly the same. However, by using stationary calculation, the stationary model can remarkably speed up the computational process. Due to the advantage of calculating speed, the stationary model can be used to characterize and design large-scale HTS applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Research on Superconducting Magnet in a Superconducting Synchronous Generator

Design and test of superconducting magnet (SM) for a superconducting synchronous generator (SSG) are studied in this paper. First, a suitable superconducting material is selected to make SM by analyzing the characteristics and application background of the five practical superconducting materials, and the applicability and feasibility of the selected superconducting materials are demonstrated. Second, the initial structural parameters of the SM are obtained by finite element method (FEM). In addition, the influence of structural parameters of SM on the performance of SSG is analyzed by simulation. Then, an optimal algorithm combining response surface methodology (RSM) and particle swarm optimization (PSO) is proposed to optimize the structure of SM in order to give full play to the advantages of SM. The simulated results verify that the proposed optimal algorithm is feasible and effective. Finally, the prototype of SM is made by the optimized size, and the SM is tested at a low temperature to obtain the critical temperature and critical current of the SM; the experimental results verify that the designed SM can meet the requirements of SSG.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Tunable and angle-insensitive plasmon resonances in graphene ribbon arrays with multispectral diffraction response

Plasmon resonances in graphene ribbon arrays are investigated numerically by means of the Finite Element Method. Numerical analysis shows that a series of multipolar resonances take place when graphene ribbon arrays are illuminated by a TM polarized electromagnetic wave. Moreover, these resonances are angle-independent, and can be tuned greatly by the width and the doping level of the graphene ribbons. Specifically, we demonstrate that for graphene arrays with several sets of graphene ribbons, which have different widths or doping levels, each of these multipolar resonances will be split into several ones. In addition, as plasmon resonances can confine electromagnetic field at the ribbon edges, graphene ribbons with different widths or doping levels offer intriguing application for electrically tunable spectral imaging.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Effect of annealing on the laser induced damage of polished and CO{sub 2} laser-processed fused silica surfaces

We investigate the effect of different heat treatments on the laser-induced damage probabilities of fused silica samples. Isothermal annealing in a furnace is applied, with different temperatures in the range 700–1100 °C and 12 h annealing time, to super-polished fused silica samples. The surface flatness and laser damage probabilities at 3 ns, 351 nm are measured before and after the different annealing procedures. We have found a significant improvement of the initial laser damage probabilities of the silica surface after annealing at 1050 °C for 12 h. A similar study has been conducted on CO{sub 2} laser-processed sites on the surface of the samples. Before and after annealing, we have studied the morphology of the sites, the evolution of residual stress, and the laser-induced damage threshold measured at 351 nm, 3 ns. In this case, we observe that the laser damage resistance of the laser created craters can reach the damage level of the bare fused silica surface after the annealing process, with a complete stress relieve. The obtained results are then compared to the case of local annealing process by CO{sub 2} laser irradiation during 1 s, and we found similar improvements in both cases. The different results obtained in the study are compared to numerical simulations made with a thermo-mechanical model based on finite-element method that allows the simulation of the isothermal or the local annealing process, the evolution of stress and fictive temperature. The simulation results were found to be very consistent with experimental observations for the stresses evolution after annealing and estimation of the heat affected area during laser-processing based on the density dependence with fictive temperature. Following this work, the temperature for local annealing should reach 1330–1470 °C for an optimized reduction of damage probability and be below the threshold for material removal, whereas furnace annealing should be kept below the annealing point to avoid sample deformation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Compact spherical neutron polarimeter using high-T{sub c} YBCO films

We describe a simple, compact device for spherical neutron polarimetry measurements at small neutron scattering angles. The device consists of a sample chamber with very low (<0.01 G) magnetic field flanked by regions within which the neutron polarization can be manipulated in a controlled manner. This allows any selected initial and final polarization direction of the neutrons to be obtained. We have constructed a prototype device using high-T{sub c} superconducting films and mu-metal to isolate regions with different magnetic fields and tested device performance in transmission geometry. Finite-element methods were used to simulate the device’s field profile and these have been verified by experiment using a small solenoid as a test sample. Measurements are reported using both monochromatic and polychromatic neutron sources. The results show that the device is capable of extracting sample information and distinguishing small angular variations of the sample magnetic field. As a more realistic test, we present results on the characterization of a 10 μm thick Permalloy film in zero magnetic field, as well as its response to an external magnetic field.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2017

A New Fokker–Planck Approach for the Relaxation-driven Evolution of Galactic Nuclei

We present an approach for simulating the collisional evolution of spherical isotropic stellar systems based on the one-dimensional Fokker–Planck equation. A novel aspect is that we use the phase volume as the argument of the distribution function instead of the traditionally used energy, which facilitates the solution. The publicly available code PhaseFlow implements a high-accuracy finite-element method for the Fokker–Planck equation, and can handle multiple-component systems, optionally with the central black hole and taking into account loss-cone effects and star formation. We discuss the energy balance in the general setting, and in application to the Bahcall–Wolf cusp around a central black hole, for which we derive a perturbative solution. We stress that the cusp is not a steady-state structure, but rather evolves in amplitude while retaining an approximately ρ∝r{sup −7/4} density profile. Finally, we apply the method to the nuclear star cluster of the milky Way, and illustrate a possible evolutionary scenario in which a two-component system of lighter main-sequence stars and stellar-mass black holes develops a Bahcall–Wolf cusp in the heavier component and a weaker ρ∝r{sup −3/2} cusp in the lighter, visible component, over the period of several Gyr. The present-day density profile is consistent with the recently detected mild cusp inside the central parsec, and is weakly sensitive to initial conditions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Modelling and simulation of a porous core photonic crystal fibre for terahertz wave propagation

A porous core photonic crystal fibre based on conventional hexagonal lattice cladding is proposed for propagating terahertz radiation. The structure is designed and theoretically investigated using full vectorial finite element method. Simulation results show that at 300 µm core diameter, with a high porosity of 85%, and an operating frequency of 1.3 THz, the proposed fibre reduces the bulk absorption loss of cyclic olefin copolymer (TOPAS) by about 81%, which corresponds to an ultra-low effective material loss value of 0.039 cm{sup −1}. Furthermore, the proposed fibre shows near zero dispersion coefficient of 0.47 ps/THz/cm with an extremely small variation of 0.05 over a broad 1.3 THz bandwidth; with confinement and bending losses investigated and found to be negligibly low. It is anticipated that the proposed waveguide can potentially be used for short range transmission of terahertz radiation in the communication window.

42 ENGINEERING↗
DOE OSTI2014

Arbitrary-level hanging nodes for adaptive hphp-FEM approximations in 3D

In this paper we discuss constrained approximation with arbitrary-level hanging nodes in adaptive higher-order finite element methods (hphp-FEM) for three-dimensional problems. This technique enables using highly irregular meshes, and it greatly simplifies the design of adaptive algorithms as it prevents refinements from propagating recursively through the finite element mesh. The technique makes it possible to design efficient adaptive algorithms for purely hexahedral meshes. We present a detailed mathematical description of the method and illustrate it with numerical examples.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2015

Structure-dependent localized surface plasmon resonance characteristics and surface enhanced Raman scattering performances of quasi-periodic nanoarrays: Measurements and analysis

A set of periodic and quasi-periodic Au nanoarrays with different morphologies have been fabricated by using electron beam lithography technique, and their optical properties have been examined experimentally and analyzed theoretically by scanning near-field optical microscope and finite element method, respectively. Results present that the localized surface plasmon resonance of the as-prepared Au nanoarrays exhibit the structure-depended characteristics. Comparing with the periodic nanoarrays, the quasi-periodic ones demonstrate stronger electric field enhancement, especially for Thue-Morse nanoarray. Meanwhile, the surface enhanced Raman scattering (SERS) spectra of 4-mercaptobenzoic acid molecular labeled nanoarrays show that the quasi-periodic nanoarrays exhibit distinct SERS enhancement, for example, a higher enhancement factor of ∼10{sup 7} is obtained for the Thue-Morse nanoarray consisted of square pillars of 100 nm size. Therefore, it is significant to optimally design and fabricate the chip-scale quasi-periodic nanoarrays with high localized electric field enhancement for SERS applications in biosensing field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Magnetic and Thermal Design of HTS Quadrupole Magnet for Newly Developed Superconducting Proton Cyclotron Beam Line

A magnetic and thermal design of high-temperature superconducting (HTS) quadrupole magnet for newly developed superconducting proton cyclotron was presented in this research. With superconducting technology, the design can reduce the magnet size and remove the heat loads more efficiently. Calculations are conducted with finite element method (FEM) to study the quadrupole magnet’s magnetic and thermal properties. For the magnet cold-mass and cryostat system, the heat load is mainly generated by conduction and radiation. Multilayer thermal insulation and G10 supports are used to restrict them. The magnetic field distribution in the HTS coils is calculated to consider the critical current degradation of the YBCO tapes. The effects on the field gradient quality in relation to the pole tip profile and end chamfer are analyzed. The sixth multipole component in the integral field can be improved with end chamfer. Finally, a 20 T/m HTS quadrupole magnet with integral field uniformity less than 0.05% out to 75% of the inscribed radius is proposed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

A MATLAB-based finite-element visualization of quantum reactive scattering. I. Collinear atom-diatom reactions

We restate the application of the finite element method to collinear triatomic reactive scattering dynamics with a novel treatment of the scattering boundary conditions. The method provides directly the reactive scattering wave function and, subsequently, the probability current density field. Visualizing these quantities provides additional insight into the quantum dynamics of simple chemical reactions beyond simplistic one-dimensional models. Application is made here to a symmetric reaction (H+H{sub 2}), a heavy-light-light reaction (F+H{sub 2}), and a heavy-light-heavy reaction (F+HCl). To accompany this article, we have written a MATLAB code which is fast, simple enough to be accessible to a wide audience, as well as generally applicable to any problem that can be mapped onto a collinear atom-diatom reaction. The code and user's manual are available for download from http://www2.chem.umd.edu/groups/alexander/FEM.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗