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

Strain relaxation of GaAs/Ge crystals on patterned Si substrates

We report on the mask-less integration of GaAs crystals several microns in size on patterned Si substrates by metal organic vapor phase epitaxy. The lattice parameter mismatch is bridged by first growing 2-μm-tall intermediate Ge mesas on 8-μm-tall Si pillars by low-energy plasma enhanced chemical vapor deposition. We investigate the morphological evolution of the GaAs crystals towards full pyramids exhibiting energetically stable (111) facets with decreasing Si pillar size. The release of the strain induced by the mismatch of thermal expansion coefficients in the GaAs crystals has been studied by X-ray diffraction and photoluminescence measurements. The strain release mechanism is discussed within the framework of linear elasticity theory by Finite Element Method simulations, based on realistic geometries extracted from scanning electron microscopy images.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Finite element approximation of the radiative transport equation in a medium with piece-wise constant refractive index

The radiative transport equation can be used as a light transport model in a medium with scattering particles, such as biological tissues. In the radiative transport equation, the refractive index is assumed to be constant within the medium. However, in biomedical media, changes in the refractive index can occur between different tissue types. In this work, light propagation in a medium with piece-wise constant refractive index is considered. Light propagation in each sub-domain with a constant refractive index is modeled using the radiative transport equation and the equations are coupled using boundary conditions describing Fresnel reflection and refraction phenomena on the interfaces between the sub-domains. The resulting coupled system of radiative transport equations is numerically solved using a finite element method. The approach is tested with simulations. The results show that this coupled system describes light propagation accurately through comparison with the Monte Carlo method. It is also shown that neglecting the internal changes of the refractive index can lead to erroneous boundary measurements of scattered light.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Derivation of a continuum model and the energy law for moving contact lines with insoluble surfactants

A continuous model is derived for the dynamics of two immiscible fluids with moving contact lines and insoluble surfactants based on thermodynamic principles. The continuum model consists of the Navier-Stokes equations for the dynamics of the two fluids and a convection-diffusion equation for the evolution of the surfactant on the fluid interface. The interface condition, the boundary condition for the slip velocity, and the condition for the dynamic contact angle are derived from the consideration of energy dissipations. Different types of energy dissipations, including the viscous dissipation, the dissipations on the solid wall and at the contact line, as well as the dissipation due to the diffusion of surfactant, are identified from the analysis. A finite element method is developed for the continuum model. Numerical experiments are performed to demonstrate the influence of surfactant on the contact line dynamics. The different types of energy dissipations are compared numerically.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

GaAs/Ge crystals grown on Si substrates patterned down to the micron scale

Monolithic integration of III-V compounds into high density Si integrated circuits is a key technological challenge for the next generation of optoelectronic devices. In this work, we report on the metal organic vapor phase epitaxy growth of strain-free GaAs crystals on Si substrates patterned down to the micron scale. The differences in thermal expansion coefficient and lattice parameter are adapted by a 2-μm-thick intermediate Ge layer grown by low-energy plasma enhanced chemical vapor deposition. The GaAs crystals evolve during growth towards a pyramidal shape, with lateral facets composed of (111) planes and an apex formed by (137) and (001) surfaces. The influence of the anisotropic GaAs growth kinetics on the final morphology is highlighted by means of scanning and transmission electron microscopy measurements. The effect of the Si pattern geometry, substrate orientation, and crystal aspect ratio on the GaAs structural properties was investigated by means of high resolution X-ray diffraction. The thermal strain relaxation process of GaAs crystals with different aspect ratio is discussed within the framework of linear elasticity theory by Finite Element Method simulations based on realistic geometries extracted from cross-sectional scanning electron microscopy images.

36 MATERIALS SCIENCE↗
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↗