Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite element method”

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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Prediction method of basic domain structure in Fe3%Si(110) single crystal with grooved surface

This paper proposes the method to accurately predict the 180° basic domain width (D{sub w}) in demagnetized states of the grooved Fe3%Si(110) single crystal with the tilt angle of [001] out of the sheet surface (β). The evaluation of D{sub w} enables the estimation of the anomalous eddy current losses. In this paper, D{sub w} is optimized to minimize the magnetic Gibbs free energy represented by vector potentials using the finite element method and the conjugate gradient method. The μ*-method is adopted to approximate the magnetization relaxation. The stray field energy generated by the magnetic charges occurring on both grooves cross section and sheet surfaces is considered in our proposed method. The validity of the proposed method was confirmed by comparison with the observed D{sub w}. As a result, we could reveal the β dependence of D{sub w} against the groove depth. Moreover, the theoretical threshold of the domain refinement due to the grooves has been suggested.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

An FEM-MLMC algorithm for a moving shutter diffraction in time stochastic model

We consider a moving shutter and non-deterministic generalization of the diffraction in time model introduced by Moshinsky several decades ago to study a class of quantum transients. We first develop a moving-mesh finite element method (FEM) to simulate the determisitic version of the model. We then apply the FEM and multilevel Monte Carlo (MLMC) algorithm to the stochastic moving-domain model for simulation of approximate statistical moments of the density profile of the stochastic transients.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

Hierarchical Multiscale Process Modeling of A Textile Composite Y-Joint for the Aurora D8 Aircraft

Woven polymer matrix composites (PMCs) are widely used in aerospace applications for their favorable specific properties. However, the behavior of the inter- and intra-tow matrix during manufacturing of woven PMCs is difficult to isolate experimentally due to geometric and processing complexities. Classical closed-form homogenization approaches allow for composite thermos-chemo-mechanical properties to be characterized at any fiber volume fraction. In this work, virtual curing is applied to the inter- and intra-tow matrix materials of a plain weave repeating unit cell (RUC). Thermo-chemo-mechanical property evolution within the tows is calculated using the Concentric Cylinder Model (CCM) and the Rule of Mixtures (ROM). Curing is implemented through user-written subroutines within the commercial finite element method software Abaqus. Virtual curing results are compared with those from a linear-elastic thermal cool down simulation in Abaqus. Results indicate that tow property evolution during curing has a non-negligible effect on the final stress state.

curing↗
DOE OSTI2016

Numerical and experimental study of Lamb wave propagation in a two-dimensional acoustic black hole

The propagation of laser-generated Lamb waves in a two-dimensional acoustic black-hole structure was studied numerically and experimentally. The geometrical acoustic theory has been applied to calculate the beam trajectories in the region of the acoustic black hole. The finite element method was also used to study the time evolution of propagating waves. An optical system based on the laser-Doppler vibration method was assembled. The effect of the focusing wave and the reduction in wave speed of the acoustic black hole has been validated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

The Young's modulus of high-aspect-ratio carbon/carbon nanotube composite microcantilevers by experimental and modeling validation

This paper reports the Young's modulus of a carbon nanotube (CNT)-reinforced carbon/CNT (C/CNT) composite microcantilevers measured by laser Doppler vibrometer and validated by finite element method. Also, the microfabrication process of the high-aspect-ratio C/CNT microcantilever arrays based on silicon micromolding and pyrolysis is presented in detail. With the in-plane natural resonant frequencies of the microcantilevers measured by a laser Doppler vibrometer, a single degree of freedom (SDoF) model based on Euler-Bernoulli (E-B) beam theory is used to calculate the Young's modulus of this composite. To figure out whether this SDoF model can be applied to these composite microcantilevers, the finite element (FE) simulation of these microcantilevers was performed. The Young's modulus of C/CNT composite microcantilevers fabricated by the pyrolysis process at 600 °C is 9391 MPa, and a good agreement between the results from experiments and FE simulation is obtained.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2015

Alignment of graphene oxide nanostructures between microgap electrodes via dielectrophoresis for hydrogen gas sensing applications

Graphene oxide (GO) nanostructures have been aligned between conducting electrodes via dielectrophoresis (DEP) with different electrical configurations. The arrangement of ground with respect to peak-to-peak voltage (V{sub pp}) plays a crucial role in manipulating the GO nanostructures. Grounds on both sides of the V{sub pp} electrode give an excellent linking of GO nanostructures which is explained by scanning electron microscopy and current-voltage characteristics. A finite element method simulation explains the electric field and voltage variation profile during DEP process. The optimized aligned GO nanostructures are used as hydrogen gas sensor with a sensitivity of 6.0% for 800 ppm hydrogen gas.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Crash simulation of UNS electric vehicle under frontal front impact

Sebelas Maret University has been developing an Electric Vehicle namely SmarT-EV UNS. The main structure of the car are chasis and body. The chasis is made from steel and the body is made from fiberglass composite. To ensure the safety of the car, both static and dynamic tests were carried out to these structures, including their materials, like: tensile test, bending test, and impact test. Another test needed by this vehicle is crashworthiness test. To perform the test, it is needed complex equipments and it is quite expensive. Another way to obtain vehicle crashworthiness behaviour is by simulate it. The purpose of this study was to simulate the response of the Smart-EV UNS electric vehicle main structure when crashing rigid barrier from the front. The crash simulation was done in according to the NHTSA (National Highway Traffic Safety Administration) within the speed of the vehicle of 35 mph. The UNS Electric Vehicle was modelled using SolidWorks software, and the simulation process was done by finite element method using ANSYS software. The simulation result showed that the most internal impact energy was absorbed by chassis part. It absorbed 76.2% of impact energy, then the base absorbed 11.3 %, while the front body absorbed 2.5 %, and the rest was absorbed by fender, hood, and other parts.

42 ENGINEERING↗
DOE OSTI2014

TH-A-BRF-08: Deformable Registration of MRI and CT Images for MRI-Guided Radiation Therapy

Purpose: To evaluate the quality of a commercially available MRI-CT image registration algorithm and then develop a method to improve the performance of this algorithm for MRI-guided prostate radiotherapy. Methods: Prostate contours were delineated on ten pairs of MRI and CT images using Eclipse. Each pair of MRI and CT images was registered with an intensity-based B-spline algorithm implemented in Velocity. A rectangular prism that contains the prostate volume was partitioned into a tetrahedral mesh which was aligned to the CT image. A finite element method (FEM) was developed on the mesh with the boundary constraints assigned from the Velocity generated displacement vector field (DVF). The resultant FEM displacements were used to adjust the Velocity DVF within the prism. Point correspondences between the CT and MR images identified within the prism could be used as additional boundary constraints to enforce the model deformation. The FEM deformation field is smooth in the interior of the prism, and equal to the Velocity displacements at the boundary of the prism. To evaluate the Velocity and FEM registration results, three criteria were used: prostate volume conservation and center consistence under contour mapping, and unbalanced energy of their deformation maps. Results: With the DVFs generated by the Velocity and FEM simulations, the prostate contours were warped from MRI to CT images. With the Velocity DVFs, the prostate volumes changed 10.2% on average, in contrast to 1.8% induced by the FEM DVFs. The average of the center deviations was 0.36 and 0.27 cm, and the unbalance energy was 2.65 and 0.38 mJ/cc3 for the Velocity and FEM registrations, respectively. Conclusion: The adaptive FEM method developed can be used to reduce the error of the MIbased registration algorithm implemented in Velocity in the prostate region, and consequently may help improve the quality of MRI-guided radiation therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Effects of Nucleus Density and Dendritic Growth Influenced by Ti and Nb on Solidification Structure of Fe-18 Pct Cr Ferritic Stainless Steel

The effects of nucleus density and dendritic growth influenced by Ti and Nb on the solidification structure of Fe-18 pct Cr ferritic stainless steel (FSS) ingots were investigated. The equiaxed zone ratios (EZR) were 8.5, 28.8, 67.5, and 57.3 pct in the blank, Nb-bearing, Ti-bearing, and Ti-Nb-bearing FSS ingots, respectively. The grain size was the smallest in the Ti-Nb-bearing ingot. Silicon oxide was present in the blank sample and Nb-bearing ingot, but the complex nuclei that were composed of TiO{sub x}-Ti(C, N) and TiO{sub x}-(Ti, Nb){sub x}(C, N){sub y} existed in the Ti-bearing and Ti-Nb-bearing ingots, respectively. The bulk density of the complex nucleus decreased by 70 pct when the Ti content decreased from 0.16 to 0.14 mass pct. Thermodynamic calculations indicate that the phase transformation from Ti{sub 3}O{sub 5} to Ti{sub 2}O{sub 3} was delayed in the ingot with a lower Ti content and that the Ti{sub 2}O{sub 3} shows more potential for TiN nucleation compared with Ti{sub 3}O{sub 5} according to the disregistry theory. The increased EZR mechanism was illustrated by a cellular-automaton–finite-element method with actual distribution statistics of the complex nucleus and dendritic-growth-kinetic analysis. The heterogeneous nucleation effect of the complex nucleus by Ti addition and the variations in dendritic-growth kinetics by Ti and Nb contribute to the EZR increase of FSS.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Formation of the domain structure in CLN under the pyroelectric field induced by pulse infrared laser heating

The evolution of the self-assembled quasi-regular micro- and nanodomain structures after pulse infrared laser irradiation in congruent lithium niobate crystal was studied by in situ optical observation. Several scenarios of domain kinetics represented covering of the irradiated zone by nets of the separated domain chains and rays have been revealed. The time dependence of the total domain length was analyzed in terms of modified Kolmogorov-Avrami theory. The domain structure evolution was attributed to the action of pyroelectric field appeared during cooling. The time dependence of the spatial distribution of the pyroelectric field during pulse laser heating and subsequent cooling was calculated by finite element method. The results of computer simulation allowed us to explain the experimental results and can be used for creation of tailored domain structures thus opening the new abilities of the submicron-scale domain engineering in ferroelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Numerical Research on Magnetic Field, Temperature Field and Flow Field During Melting and Directionally Solidifying TiAl Alloys by Electromagnetic Cold Crucible

The electromagnetic cold crucible (EMCC) technique is an effective method to melt and directionally solidify reactive and high-temperature materials without contamination. The temperature field and fluid flow induced by the electromagnetic field are very important for melting and controlling the microstructure. In this article, a 3D EMCC model for calculating the magnetic field in the charges (TiAl alloys) using the T-Ω finite element method was established and verified. Magnetic fields in the charge under different electrical parameters, positions and dimensions of the charge were calculated and analyzed. The calculated results show that the magnetic field concentrates in the skin layer, and the magnetic flux density (B) increases with increasing of the frequency, charge diameter and current. The maximum B in the charge is affected by the position of the charge in EMCC (h{sub 1}) and the charge height (h{sub 2}), which emerges at the middle of coils (h{sub c}) when the relationship of h{sub c} < h{sub 1} + h{sub 2} < h{sub c} + δ is satisfied. Lower frequency and smaller charge diameter can improve the uniformity of the magnetic field in the charge. Consequently, the induced uniform electromagnetic stirring weakens the turbulence and improves temperature uniformity in the vicinity of the solid/liquid (S/L) interface, which is beneficial to forming a planar S/L interface during directional solidification. Based on the above conclusions, the TiAlNb alloy was successfully melted with lower power consumption and directionally solidified by the square EMCC.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A Fully Computable A Posteriori Error Estimate for the Stokes Equations on Polytopal Meshes

Here in this paper, we present a simple a posteriori error estimate for the weak Galerkin finite element method for the Stokes equation. This residual type estimator can be applied to general meshes such as polytopal mesh or meshes with hanging nodes. The reliability and efficiency of the estimator are proved in this paper. Five numerical tests demonstrate the effectiveness and flexibility of the adaptive mesh refinement guided by the designed error estimator.

finite element methods↗
DOE OSTI2019

An Integrated Numerical Method to Predict the Strain Sensing Behavior of Flexible Conductive Fibers under Electric–Thermal Conditions

Flexible conductive textiles have attracted an increasing amount of attention due to being flexible and stretchable with multi-functions, in addition to the strain sensing function. In this paper, an integrated numerical method is established to predict the strain sensing behavior of PPy-coated flexible conductive fibers. Based on our previous investigation on the strain sensing behavior of the PPy-coated elastic conductive fibers, it has been known that the main strain sensing mechanisms of the conductive fibers are micro-cracks opening under stretching and closing under unloading. Therefore, in the numerical approach, extended finite element method is used to simulate the multi-crack initiation and propagation. Then, by establishing a representative element model including a single crack, the resistance variation with the increase in crack width/depth can be obtained using the electric–thermal package in ABAQUS software. Combining the relationship between crack width and strain, the relationship between fractional increment in resistance (ΔR/R{sub 0}) and the strain can be obtained. The prediction of the fractional increment in resistance (ΔR/R{sub 0}) versus strain is consistent with the experimental results of PPy-coated Lycra fibers. The new integrated method is further verified by experimental results of PPy-coated XLA fibers, showing that the integrated numerical approach can capture the complicated strain sensing mechanisms of PPy-coated elastic fibers and predict the strain sensing behavior of flexible conductive fibers. This approach will make it possible to shorten the period, and reduce the cost in R&D of the PPy-coated electrically conductive fibers, comparing with the corresponding experimental work.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Theoretical realization of three-dimensional nanolattice structure fabrication based on high-order waveguide-mode interference and sample rotation

A theoretical method of using high-order waveguide-mode interference combined with sample rotation is proposed to fabricate different kinds of three-dimensional nanolattice structures. The lithography sample is assisted by an asymmetric metal-cladding dielectric waveguide structure. High-order waveguide modes can be excited in the photoresist layer by both transverse magnetic and transverse electric polarized light. By utilizing multiple exposures of high-order waveguide-mode interference combined with sample-rotation, various three-dimensional nanolattice structures can be obtained. The resulting optical field distributions are simulated using the finite element method in this study, and sample rotation is expressed by coordinate matrix transformation. As examples, the three-dimensional optical field distributions resulting from fifth-order waveguide-mode interference were simulated with 90° and 60° sample rotations and multiple exposures. The results show that a quasi-cuboid structure with a simple tetragonal arrangement and a quasi-hexagonal structure with a hexagonal close-packed lattice can be obtained. Moreover, the numerical simulation results revealed that the shapes, sizes, arrangements, and periods of the structures can be controlled by the rotation method, photoresist thickness, waveguide modes used for exposure, and so on.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2014

A 2D Particle in Cell model for ion extraction and focusing in electrostatic accelerators

Negative ions are fundamental to produce intense and high energy neutral beams used to heat the plasma in fusion devices. The processes regulating the ion extraction involve the formation of a sheath on a scale comparable to the Debye length of the plasma. On the other hand, the ion acceleration as a beam is obtained on distances greater than λ{sub D}. The paper presents a model for both the phases of ion extraction and acceleration of the ions and its implementation in a numerical code. The space charge of particles is deposited following usual Particle in Cell codes technique, while the field is solved with finite element methods. Some hypotheses on the beam plasma transition are described, allowing to model both regions at the same time. The code was tested with the geometry of the NIO1 negative ions source, and the results are compared with existing ray tracing codes and discussed.

43 PARTICLE ACCELERATORS↗
NASA NTRS1996

Computing Radiation From Axisymmetric Waveguide Feed Horns

Hybrid finite-element method developed for computing radiation patterns and reflection coefficients of axisymmetric waveguide feed horns that include dielectric bodies like rods and radomes. Dielectric inhomogeneities introduced by such bodies and even anisotropies of some dielectric materials pose no special problem in method because requisite dielectric properties readily incorporated into finite elements used to model bodies.

Chinn, Gilbert C.↗
DOE OSTI2016

Standalone BISON Through VERA: User's Manual

Activities to incorporate fuel performance capabilities into the Virtual Environment for Reactor Applications (VERA) are receiving increasing attention. The multiphysics emphasis is expanding as the neutronics (MPACT) and thermal-hydraulics (CTF) packages are becoming more mature. Capturing the finer details of fuel phenomena (swelling, densification, relocation, gap closure, etc.) is the natural next step in the VERA Core Simulator (VERA-CS) development process since these phenomena are currently not directly taken into account. While several codes could be used to accomplish this, the BISON fuel performance code being developed by the Idaho National Laboratory (INL) is the focus of ongoing work in the Consortium for Advanced Simulation of Light Water Reactors (CASL). Built on INL’s MOOSE framework, BISON uses the finite element method for geometric representation and a Jacobian-free Newton-Krylov (JFNK) scheme to solve systems of partial differential equations for various fuel characteristic relationships. There are several modes of operation in BISON, but, for this work, it uses a 2D azimuthally symmetric (R-Z) smeared-pellet model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1990

Finite-element models of continental extension

Numerical models of the initial deformation of extending continental lithosphere, computed to investigate the control of preexisting thermal and mechanical heterogeneities on the style of deformation, are presented. The finite element method is used to calculate deformation with a viscoelastic-plastic model for the lithosphere. Comparisons of the results of analytic models and finite-element models using this method show that good results may be obtained by the numerical technique, even with elements containing both brittle and viscoelastic sampling points. It is shown that the gross style of initial extensional deformation is controlled by the depth and width of the initial heterogeneity which localizes deformation.

Lynch, H. David↗